一種隨鑽泥漿脈衝信號的處理方法

摘要：無線隨鑽測量系統中的泥漿脈衝信號受到各種噪聲的幹擾，需要對採集到的信號進行處理還原，以實時監測井底狀況。研究了泥漿脈衝信號特徵，設計了對其基於最大似然估計閾值去噪、平滑及去除基線漂移的信號處理方法。利用該方法進行信號處理，能較好的恢復信號的特徵。
關鍵詞：泥漿脈衝信號；小波變換；最大似然估計閩值；平滑；基線漂移

利用鑽井液脈衝進行信息傳輸，可靠性較高，傳輸距離遠，更符合鑽井的實際情況，是國內通用的傳輸方式。在泥漿信道MWD中，由於現場測量條件的影響，安裝在立管上檢測泥漿壓力波動的壓力傳感器的輸出，不僅包含井下傳來的有用信號，而且含有由於泥漿泵壓縮泥漿而引起的大幅度周期性壓力脈動，其他各種機械作用所引起的壓力波動以及隨機噪聲。因此，在井口處所採集到的信號中。有用信號完全淹沒在各種噪聲中，所以，對泥漿脈衝信號進行去噪處理就成了一個重要的研究課題。當前信號去噪處理的一般方法是通過小波變化進行去噪。小波變換突破了傅立葉變換在時域沒有任何分辨力的限制，可以對指定頻帶和時間段內的信號成份進行分析。在時域和頻域同時具有良好的局部化性質，並且由於對頻率成份採用逐漸精細的時域或頻域取樣步長，從而可以聚焦到信號的任何細節。而實際採集到的信號存在基線漂移，去噪後的信號又需要經過平滑處理，因此，選取合一種合適的信號處理方法是一個關鍵問題。

1 小波變換
1．1 概念介紹

在尺度上的伸縮和時域上的平移來分析信號，具有很好的時域和頻域局部特性。
1．2 小波變換去噪
基於小波變換在頻域和時域的特性，小波變換在信號去噪方面有著很好的應用。常用的去噪方法是閾值去噪，基於小波變換的閾值去噪方法步驟如下：
1)選擇合適的小波，對給定的信號進行小波變換，得到小波變換係數W：
2)計算閾值，選擇合適的閾值方法(硬閾值或軟閾值)對小波係數進行取捨，得到新的小波係數Wδ；
3)對得到的係數進行逆變換，得到去噪後數據。
其中，硬閾值的函數表達式為η(ω)=ωI(|ω|>T)，如圖1所示。軟閾值的函數表達式η(ω)=(ω-sgn(ω)T)I(|ω|>T)，如圖2所示。

本文引用地址：http://www.eepw.com.cn/article/201610/309098.htm

硬閾值的缺點是在某些點會產生間斷，軟閾值使得重建信號比較平滑，但可能會造成邊緣模糊等失真現象。因此，提出來基於最大似然函數的閾值估計方法。
1．3 最大似然估計閾值法去噪
稀疏性的概念來自資訊理論，其中採用熵測度來量化分布的稀疏性，與之相對應的是均勻性。熵值越小，說明分布越稀疏，均勻性則越差。稀疏分布並沒有一個定量的定義，通常是指概率密度函數在零點處存在一個尖峰的分布。在井口採集到的信號是泥漿正脈衝信號與噪聲的加性混合信號。脈衝信號是一種常見的信號，其概率密度函數為典型的稀疏分布。對符合稀疏分布形式的信號，Hyvarinen根據最大似然原則得到以下閾值準則：

下圖給出了最大似然估計原則閾值與硬、軟閾值處理結果比較，通過比較發現，當小波係數小於閾值時，其處理結果都是將小波係數置於零，當小波係數大於閾值時，其處理結果在軟硬閾值之間。

1．4 模擬仿真
根據最大似然估計閾值去噪原理，對MATLAB中自帶信號源noisbump分別用最大似然估計閾值、軟閾值、硬閾值進行去噪處理結果如圖4所示。

定義信噪比如下：，其中s(i)是原始信號，sn(i)是最大似然估計後去噪信號，同理ss(i)是軟閾值去噪後信號，sh(i)是硬閾值去噪後信號。
定義均方差公式如下：，由於原始信號是噪聲信號，因此，去噪後信號與原始信號的均方差越大，說明去噪效果越好。

通過對比SNR和MSE可以發現，最大似然估計閾值去燥處理後的結果明顯優於軟閾值和硬閾值處理後的結果。

2 泥漿脈衝信號平滑
為了消弱幹擾信號的影響，提高曲線的光滑度，需要對去噪後的信號進行數據平滑處理。常用的信號平滑方法有直線滑動平滑法。直線滑動平均法就是利用最小二乘原理對離散數據進行線性平滑的方法，該方法主要根據某點臨近的採樣點的波幅來對該點進行波幅修正，從而達到對波形光滑去噪的目的。一般取5個近鄰點的數據點進行計算，其計算公式如下式所示：


3 泥漿脈衝信號的基線漂移校正
3．1 去基線目的
由於鑽井現場各種電機、磁場的影響，加之在傳輸過程中受到泵壓的幹擾，地面通過壓力傳感器接收到的鑽井液脈衝信號完全被各種噪聲淹沒，因此鑽井液脈衝信號中的成分非常複雜，信號不穩定，存在基線漂移，不能直接用到鑽井現場中解碼得到井下的各種參數，在解碼的過程前需要對信號進行矯正處理。
3．2 去基線方法
文中採用移動窗口中值濾波法實現基線漂移校正。中值濾波是基於排序統計理論的一種能有效抑制噪聲的非線性信號處理技術，中值濾波的基本原理是把數字序列中一點的值用該點的一個鄰域中各點值的中值代替，讓周圍的像素值接近的真實值，從而消除孤立的噪聲點。
假設待處理原始泥漿脈衝信號為S1，信號長度為L，處理過程具體描述如下：
1)選擇合適的窗寬M，為一般情況下K取奇數，如M=101；
2)為了防止出現邊緣效應，對原始信號S1的兩端進行延拓，得到S2信號，處理過程如下式所示：

延拓後信號S2的長度為L+M-1；
3)對S2信號加窗，對窗口內信號進行中值濾波，即對窗口內信號進行排序，然後用中值來取代窗口中心點的值，移動該窗口，遍歷S2信號，擬合出漂移了的基線BL，BL與S2之間的關係如下式所示：BL(i)=median[S2(i)：S2(i+2*M)]其中，median()函數表示取中值操作，擬合得到的基線BL的長度與原始信號S1的長度相同；
4)從原始信號S1中減去BL，得到消除基線漂移後的信號S3。

4 處理結果
根據以上數據處理流程，在MATLAB編寫程序，實現相應算法，根據文獻已有的結果，選取haar小波基進行去噪處理。圖4是在中國石油某鑽井施工過程中採集到的一段完整的原始數據，包括開泵信號，數據脈衝信號及關泵信號。圖5是對原始數據截取其中一段去噪後的信號，圖6是經過平滑後得到的信號，圖7是經過去滑動平滑後得到信號，圖8是去基線後得到的信號波形。

5 結論
針對泥漿脈衝信號特徵，採用上述信號處理流程方法對基於最大似然估計小波閾值的方法對其進行去噪、平滑及去基線，可以較好的恢復出原始脈衝信號的面貌特徵，為下一步的解碼處理打下基礎，實時監測井下鑽進情況。