一種泥漿脈沖信號識別方法

李傳偉，楊亮，祝環芬，王芝峰

（1.西北工業大學自動化學院，陜西 西安710072；2.中國石油集團測井有限公司，陜西 西安710077）

0 引 言

隨鉆測井與傳統電纜測井及存儲式測井的主要區別是獲取地層參數的時效性，隨鉆測井能夠在剛鉆開地層時獲得反映原狀地層評價參數，有助于及時綜合評價地層，實時指導鉆進過程，提高儲層鉆遇率。隨鉆測井在水平井和大斜度井勘探開發中得到廣泛應用。隨著隨鉆測井方法不斷完善和儀器的研制，井下測量數據越來越多，實時傳輸到地面成為制約隨鉆測井技術應用的瓶頸［1］。本文針對廣泛使用的泥漿脈沖傳輸方式中脈沖信號識別方法開展研究，提出先使用脈沖有效上升沿法判斷脈沖大體位置，再用滑動窗口最大值法精確調整和校準脈沖位置的方法，在隨鉆測井過程中有效地識別出泥漿壓力波形中的脈沖信號。

1 信號產生及特征

隨鉆測井過程中井下探測的信號經過采集轉化為測量數據后，通過數據編碼策略驅動脈沖發生器閥體運動改變泥漿循環通道的截面積，泥漿壓力發生變化，在地面探測的泥漿壓力波形將會產生相應的脈沖［如圖1（a）所示］。這一過程將泥漿壓力變化通過測井數據傳輸到地面，在地面通過檢測壓力波形中的脈沖，根據位置關系解碼獲得相應的測量參數［2］。

泥漿脈沖信號的傳輸過程是一種能量轉換過程，傳輸特性受泥漿介質的影響，存在明顯的衰減，在地面采集檢測的壓力是一個非平穩的信號，表現在壓力變化波峰（波谷）是1條斜線，包含各種噪聲和干擾，其中噪聲表現為寬帶白噪聲，干擾表現為與泵噪聲特性相關的周期性脈沖。研究人員采取多種方法對信號進行濾波去噪處理，如中值濾波、小波濾波等去除信號中的各種噪聲。作者曾使用小波變換方法對泥漿壓力信號進行濾波處理［3－4］，效果見圖1（b）所示。圖1（b）中紅色為采集的原始信號，黑色為經過小波濾波處理后的信號，藍色為實際脈沖位置對應的信號。

從濾波處理結果看出，濾波后的壓力波形有明顯的起伏變化，從中能夠識別出脈沖信號、獲取位置關系，實現井下數據的解碼恢復。然而在環境噪聲較大、泵壓不穩時，處理后的壓力波形波幅不等、寬窄不均，使得脈沖檢測十分困難，影響井下測量數據解碼。

圖1 泥漿脈沖信號與處理

2 算法的設計與實現

傳統的脈沖位置檢測方法是設定門檻值，高于門檻的波形就認為是1個脈沖。該方法適用于波形變化較小、基值平穩的信號，但對傳輸過程中泵壓變化導致脈沖幅度變化較大、存在基線漂移時，脈沖檢測效果較差，解碼數值準確性不高［5］。

本文介紹基于濾波去噪后的壓力波形進行脈沖識別的方法。首先求取信號幅度的一次差分和差分結果的二次差分，在信號中尋找出極值點（波峰、波谷）；其次根據脈沖信號寬度特征及約束條件從極值點中篩選出波峰，并利用閾值判斷出有效上升沿；最后采取滑動窗口內面積最大原則確定脈沖所在的位置，轉換成方波序列，得到與井下編碼規則一致的脈沖信號。

2.1 波峰、波谷查找算法

脈沖識別的前提是先在波形中查找出所有局部極大值和極小值，記錄極值的位置和幅值，即確定波形點的x和y坐標。連續函數求極值的算法一般采用一階導數駐點法，求連續函數一階導數為0的點；對于采集的泥漿壓力信號的離散值函數，則采用鄰域差值［6］，依次計算相鄰2個波形間的波幅差值diff＝y2－y1，計算結果有3種情況，即diff＞0、diff＜0、diff＝0。

根據結果判斷極值點，如果diff＝0并且diff正負特性發生改變的波形點即為極值點，將所有的極值點存入極值點數組中。某組測試數據處理結果見圖2（a）。圖2（a）中綠色點為極大值和極小值點。

2.2 查找波形有效沿算法

波形的上升沿和下降沿是任意相鄰的2個極值點之間的連接線，這條連接線的2個端點為V1（x1，y1）和V2（x2，y2），都是波形上的極值點，其中x1＜x2，如果y2＞y1則為上升沿，否則為下降沿。有效沿是指y2－y1的絕對值大于一定閾值的沿，分為有效上升沿和有效下降沿；利用有效沿可以排除識別到的部分虛假波峰，剔除干擾尖峰。

依次遍歷極值點數組中的每個元素，將相鄰2個點連接作為1個波形沿，并計算波形沿的波幅差值，如果大于閾值（測試中取經驗值400）則為有效沿，記入有效沿數組。試驗中發現下降沿對于判定方波位置的作用較小，而且容易受到噪聲影響，所以在算法中均以有效上升沿為判定波峰位置的主要依據，而拋棄了有效下降沿數據。某組測試數據處理結果見圖2（b），其中綠色線表示有效上升沿。

圖2 波形信號識別

2.3 滑動窗口面積法確定方波

利用有效上升沿法確定的類正弦波波峰中，采用波峰的頂點作為中心位置確定出對應的脈沖（見圖3綠色實線脈沖）。圖3表明由于波形畸變，該方法確定的方波不能準確反映波形代表的脈沖位置，以波峰頂點作為中心點確定的脈沖及位置效果不理想。

本文在用波峰頂點確定脈沖后，采用滑動窗口面積法進一步校正脈沖位置?；瑒哟翱诿娣e法就是以波峰頂點左一個脈寬位置為起始x坐標，用固定寬度值為滑動窗口，沿x軸遞增改變x坐標，計算窗口左右2個邊界及下邊界和波形所圍成的面積，確定掃過一個波峰中面積最大時脈沖的位置。為了簡化計算，設窗口的下邊界為0，計算面積值也就轉化為計算窗口內所有點波幅度的和。求出這一組面積值的最大值，其窗口所在位置即為方波所在位置（窗口寬度等于方波脈沖寬度）（見圖3黑色虛線脈沖所示）。

圖3 波峰轉換方波算法示意圖

3 測試應用效果

采取文中介紹的方法識別壓力波形中的脈沖信號，對比分析解碼測量數據的應用分為地面模擬信號測試和現場測井采集信號試驗2個部分［7］。

3.1 模擬信號測試

測試方法是連接地面數據采集系統與井下儀器，利用井下脈沖發生器測試盒產生標準的壓力脈沖信號，地面系統采集的波形和處理結果如圖4所示。圖4中綠色為模擬產生的泥漿信號，藍色為濾波處理后的波形信號，紅色為通過脈沖識別方法識別形成的方波序列。圖4中，模擬信號沒有經過泥漿通道傳輸，噪聲很小，是1組波形規則的方波信號，經過濾波處理形成類正弦波，通過本文算法識別出的方波序列能準確反映出模擬信號中方波位置。

圖4 地面模擬脈沖識別結果

根據圖4中的紅色方波序列和編碼規則解碼得到的數據與井下存儲數據對比結果如表1所示，數值表明解碼的工具面、井斜、方位和溫度值與井下儀器中存儲的結果一致。

表1 解碼結果與儀器內存儲數據對比

3.2 現場測井信號試驗

通過在長慶油田慶陽區塊長－××井的上井試驗，地面數據采集系統采集的立管壓力傳感器探測的壓力波形和處理結果如圖5所示。圖5中綠色為原始采樣信號，藍色為小波濾波去噪處理后的信號，紅色方波是對濾波后的泥漿壓力信號進行脈沖識別結果，形成方波序列。

圖5 泥漿壓力信號與識別脈沖的對應關系

從圖5中看出，現場采集的信號噪聲較大，經過濾波處理后的波形基本反映出數據脈沖產生的變化，本文算法識別出的方波能夠反映出波形對應的脈沖位置。

根據圖5中的紅色方波序列和編碼規則解碼得到的數據與井下存儲數據對比結果如表2所示，數值表明解碼得到的工具面、井斜、方位和溫度值與井下儀器中存儲的結果一致。

表2 解碼結果與儀器內存儲數據對比

4 結 論

（1）泥漿壓力波形傳輸受鉆井環境和井深的影響較大，直接關系到隨鉆測井施工效果和測量數值的準確性，對信號處理方法提出了較高要求。

（2）通過對模擬信號和大量現場記錄數據的處理分析結果表明，采用本文的脈沖識別方法能夠從波形中識別出脈沖，為數據解碼提供準確的脈沖位置，解碼得到的結果與井下編碼數據準確率達到90%以上。

（3）有效提高數據解碼準確率，為隨鉆測井系統在各種工況環境應用奠定了基礎，已經在現場施工中進行了試驗和應用。

［1］ 賈朋.鉆井液連續波發生器設計和信號傳輸特性試驗研究［D］.北京：中國石油大學，2010.

［2］ 劉修善，蘇義腦.鉆井液脈沖信號傳輸特性研究［J］.石油鉆采工藝，2000（4）：8－10

［3］ 張偉，師奕兵，盧濤.小波神經網絡在隨鉆測量系統在泥漿信號監測中的應用研究［J］.電子測量與儀器學報，2008（6）：43－46

［4］ 李傳偉，李安宗.一種小波變換信號處理方法［J］.西安電子科技大學學報，2009，36（4）：751－755

［5］ 遠方，鞏憲鋒.隨鉆泥漿脈沖器信號處理技術研究［J］.煤礦機械，2011，32（10）：72－74

［6］ 同濟大學數學系.高等數學［M］.6版.北京：高等教育出版社，2007.

［7］ 中國石油集團測井有限公司.FELWD隨鉆測井系統操作手冊［Z］.2011.