改進的雙變量閾值函數(shù)在動液面信號去噪中的應(yīng)用

岳武峰 高丙坤

(東北石油大學(xué)電氣信息工程學(xué)院，黑龍江 大慶 163318)

改進的雙變量閾值函數(shù)在動液面信號去噪中的應(yīng)用

岳武峰 高丙坤

(東北石油大學(xué)電氣信息工程學(xué)院，黑龍江 大慶 163318)

以油井動液面聲波反射信號為研究對象，引入一種改進的雙變量閾值函數(shù)對小波變換的閾值選擇方法進行改造。該函數(shù)克服了硬閾值函數(shù)不連續(xù)和軟閾值函數(shù)存在恒定偏差的缺點，并在小波系數(shù)的絕對值小于閾值的區(qū)間上進行壓縮處理。對比去噪的結(jié)果表明：改進的方法在MSE、SNR比較中均優(yōu)于傳統(tǒng)的軟、硬閾值函數(shù)。

動液面信號 去噪 閾值函數(shù)

隨著大多數(shù)油田步入開發(fā)的后期，原油開采由于受到地層壓力下降等因素的影響，需要通過人工舉升的方法才能將原油抽取到地面。在這個過程中，油井動液面深度值作為判斷油井工作狀態(tài)的重要參數(shù)之一，可以為油田生產(chǎn)計劃的規(guī)劃及油井開發(fā)壽命的預(yù)測等提供重要的參考價值，并對提高油井產(chǎn)量，降低能耗和生產(chǎn)成本及減少油井事故發(fā)生率等有著重要的意義[1]。

傳統(tǒng)測量動液面深度值的方法有浮標測量法、井下壓力測量法等[2]，其中浮標測量法利用在井底放置浮標來測量，這種方法操作簡單，但是由于受到浮標體積、重量及油井壓力等因素的影響，能夠使用的環(huán)境十分有限。井下壓力測量法通過在井底安置壓力計，利用獲得的壓力值來計算動液面值，這種方法相比浮標法應(yīng)用更加廣泛且測量精度更高，但設(shè)備長時間置于井底會發(fā)生腐蝕、老化等現(xiàn)象，需要經(jīng)常更換，這導(dǎo)致成本相對較高，而且測量時需要停止油井工作，會耽誤正常開采，因此該方法也無法廣泛應(yīng)用于油田的實際生產(chǎn)中。隨著信號處理技術(shù)和聲學(xué)傳感器設(shè)備的不斷發(fā)展，現(xiàn)階段油田廣泛利用聲波反射法采集動液面數(shù)據(jù)。然而在實際應(yīng)用中，由于油井正常工作時井下環(huán)境復(fù)雜，采集到的聲波信號受到多種噪聲的干擾[3]，往往無法成功取得清晰準確的回波信號，增加了動液面深度值測量的難度。因此對動液面信號去噪方法的研究成為油田生產(chǎn)活動中待解決的重要問題之一。

1 動液面背景噪聲分析

油井動液面深度從地面到井底數(shù)千米不等，且油井向下的軌跡會出現(xiàn)彎曲或轉(zhuǎn)折，因此聲波向下傳播的時候會產(chǎn)生折射及漫反射等；而且在油井正常工作狀態(tài)下，井筒環(huán)空中套管氣會產(chǎn)生聲響，還包含大量的機械噪聲、自然噪聲及電磁噪聲等[4]；再隨著生產(chǎn)的進行，油井內(nèi)壓力和溫度逐漸變化使得井中各種流體、固體發(fā)生形狀和狀態(tài)的變化，容易產(chǎn)生死油帽子、泡沫油層及結(jié)蠟等現(xiàn)象。這都對動液面反射波的接收產(chǎn)生嚴重的影響，造成真實的反射信號被上述噪聲覆蓋，使得液面反射波和接箍波不易辨別，增加了信號處理的難度。

上述動液面反射信號的背景噪聲大部分具有高頻特性，具體表現(xiàn)為信號幅度隨著時間增長快速變化，比如在油井套管中流動的氣體產(chǎn)生的聲音、電機轉(zhuǎn)動產(chǎn)生的聲音及機械器件摩擦產(chǎn)生的聲音等，這些噪聲會產(chǎn)生尖銳復(fù)雜的振蕩波形夾雜在所要測量的反射聲波波形中，對反射聲波的接收產(chǎn)生很大影響，對液面反射波和接箍波的干擾極其嚴重。傳統(tǒng)的傅里葉變換方法是將整個時間軸的信號進行分析，信號在某一時間位置處任意小的變化都會使頻譜發(fā)生較大變化[5]，這導(dǎo)致傅里葉變換缺乏時頻局部化能力，不適合處理這類非平穩(wěn)信號，而近些年來發(fā)展起來的小波變換在處理瞬時信號方面顯示出明顯優(yōu)勢。

2 小波閾值去噪方法

小波變換是在短時傅里葉變換基礎(chǔ)上發(fā)展起來的一種信號時頻分析方法，具有良好的時頻局部化處理能力。目前，小波變換方法主要可分為3類：基于小波變換模極大值小波去噪、基于小波變換相關(guān)去噪和閾值去噪法[6]。其中小波閾值去噪法由于重構(gòu)信號原理簡單、計算量小且去噪效果較好，因此得到廣泛的應(yīng)用。

小波閾值去噪主要分為如下幾個步驟：選取一個小波基對含噪信號做小波變換得到低頻小波系數(shù)和高頻小波系數(shù)；選擇閾值函數(shù)和合適的閾值對各層分解的高頻小波系數(shù)進行閾值處理；對分解后的低頻系數(shù)和各層分解后經(jīng)過閾值處理的高頻小波系數(shù)做小波重構(gòu)，得到去噪信號。小波閾值去噪流程如圖1所示。

圖1 小波閾值去噪流程

小波閾值去噪法主要體現(xiàn)在閾值處理部分，閾值處理包括閾值的選擇和閾值函數(shù)的構(gòu)造，其中閾值函數(shù)的構(gòu)造是影響信號去噪效果的重點和難點。常用的閾值函數(shù)有兩種：硬閾值函數(shù)和軟閾值函數(shù)。

硬閾值函數(shù)表達式如下：

(1)

由式(1)可以看出硬閾值函數(shù)在λ和-λ處不連續(xù)，這會導(dǎo)致重構(gòu)后的信號出現(xiàn)振蕩現(xiàn)象[7]。

軟閾值函數(shù)表達式如下：

(2)

(3)

其中，0≤α≤1，μ取自然數(shù)，可以根據(jù)信號實際情況自由調(diào)節(jié)。

圖2 3種閾值函數(shù)示意圖

3 改進閾值函數(shù)去噪實驗與結(jié)果分析

為了證明改進雙變量閾值函數(shù)的有效性，對比硬閾值函數(shù)、軟閾值函數(shù)進行去噪仿真實驗，然后比較實驗結(jié)果。為了定量評價各閾值函數(shù)對在小波閾值去噪中的有效性和優(yōu)越性，采用信噪比SNR和均方誤差MSE作為對比指標。信噪比SNR和均方誤差MSE的定義如下：

式中f(t)——原始信號；

N——信號長度。

圖3 3種閾值函數(shù)對實驗室采集信號的去噪對比

從圖3中可以看出，硬閾值函數(shù)處理后的波形出現(xiàn)Pseudo-Gibbs振蕩現(xiàn)象，波形失真嚴重；軟閾值函數(shù)處理后的波形過于平滑，損失了很多原始信號特征值；改進閾值函數(shù)處理得到的波形連續(xù)性好并且很好地保留了原始信號的特征值，便于下一步信號的特征提取。表1是3種閾值函數(shù)對實驗室采集信號去噪所得SNR和MSE的對比，從表1可以看出改進閾值函數(shù)的SNR最大且MSE最小。結(jié)合圖3得出的結(jié)果，驗證了筆者改進的雙變量閾值函數(shù)對實驗室采集動液面信號的去噪效果。

表1 3種閾值函數(shù)對實驗室采集信號去噪所得SNR和MSE對比

從圖4中可以看出，采用改進閾值函數(shù)去噪得到的波形在視覺效果上要好于軟、硬閾值函數(shù)，它既保留了動液面反射信號的特征點，又能有效濾除噪聲。表2是3種閾值函數(shù)對動液面反射信號去噪所得的SNR和MSE對比，從表2對比可以看出，改進閾值函數(shù)去噪效果表現(xiàn)優(yōu)越。因此得出結(jié)論，筆者改進的方法能有效去除油井動液面反射信號的噪聲。

圖4 3種閾值函數(shù)對動液面反射信號的去噪對比

表2 3種閾值函數(shù)對動液面反射信號去噪所得SNR和MSE對比

4 結(jié)束語

研究了油井動液面聲波反射信號的主要特性，并根據(jù)其特點提出一種改進的雙變量小波閾值函數(shù)。該函數(shù)克服了硬閾值函數(shù)不連續(xù)和軟閾值函數(shù)存在固定偏差的缺點，并在小波系數(shù)的絕對值小于閾值的區(qū)間上進行壓縮處理。通過對比軟、硬閾值函數(shù)對油井動液面聲波反射信號進行去噪實驗，得到的結(jié)果表明，筆者提出的閾值函數(shù)的去噪效果較傳統(tǒng)軟、硬閾值函數(shù)更加明顯，驗證了該閾值函數(shù)對解決油井動液面反射信號去噪問題的有效性。

[1] 金立旸. 基于盲源分離的油井動液面連續(xù)監(jiān)測系統(tǒng)[D].大慶：東北石油大學(xué)，2014.

[2] 呂思平.油井動液面測量系統(tǒng)的研制[D].青島：中國石油大學(xué)，2011.

[3] 劉迎新,楊亦春,韓寶坤，等.低頻聲波油井液面檢測方法研究[J].應(yīng)用聲學(xué),2015,34(1):24～30.

[4] 王海文,林立星.基于小波變換的聲波法測油井動液面信號去噪[J].工業(yè)儀表與自動化裝備,2011,(6):56～58.

[5] 闞玲玲,高丙坤,梁洪衛(wèi)，等.小波去噪在油井動液面檢測中的應(yīng)用[J].化工自動化及儀表,2014,41(9):1009～1011.

[6] 姚建紅,林娜,付強.基于多目標函數(shù)的粒子群算法優(yōu)化小波閾值的去噪方法研究[J].化工自動化及儀表,2013,40(2):154～157.

[7] 周云龍,張慧冬,李洪偉，等.改進的小波消噪閾值方法在油氣水三相流圖像信號中的應(yīng)用[J].化工自動化及儀表,2009,36(4):49～53.

[8] Sanam T F, Shahnaz C.Noisy Speech Enhancement Based on an Adaptive Threshold and a Modified Hard Threshold Function in Wavelet Packet Domain[J].Digital Signal Processing,2013,23(3):941～951.

[9] 鄧玉娟.基于小波變換的語音閾值去噪算法研究[D].重慶:重慶大學(xué),2009.

(Continued on Page 312)

ApplicationofImprovedDouble-variableThresholdFunctioninDenoisingWorkingLevelSignals

YUE Wu-feng, GAO Bing-kun

(SchoolofElectricalEngineering&Information,NortheastPetroleumUniversity,Daqing163318,China)

Taking acoustic wave’s reflected signal in oil well’s working level as the research object, an improved double-variable threshold function was introduced to modify threshold selection method of the wavelet transform. This function overcomes the hard threshold function’s discontinuousness and soft threshold func-

TH816

A

1000-3932(2016)03-0268-05

2016-01-14(修改稿)

東北石油大學(xué)研究生創(chuàng)新科研項目(YJSCX2014-030NEPU)