測井信號的模極大值小波去噪與交替投影重構

高丙坤 周泳含 孫靈川 張 莉

(1.東北石油大學電氣信息工程學院，黑龍江 大慶 163318；2.北京郵電大學信息與通信工程學院，北京 100876)

核磁共振測井可以有效地探測、評估剩余油分布范圍、采油率及儲油效率等，該技術所獲取的原始數據是自旋回波串，包含滲透率、孔隙度、可動流體及不可動流體等眾多的物理信息[1]。由于測井中熱力和重力的作用，加之井下環境惡劣，干擾因素多，致使核磁共振測井得到的自回波信號較微弱，信噪比低。核磁共振測井回波信號為非線性非平穩性信號，故筆者通過小波變換對信號進行去噪處理，以提高分辨率，對測井信號做出準準的解釋。

核磁共振測井儀器探測到的自旋回波為指數衰減信號，實測的1 000組樣本形成的回波信號如圖1所示，有用信號被大量噪聲湮沒，難以表征信號的物理信息。

圖1 核磁共振測井儀器實測信號

信號的奇異性需要利用Lipschitz指數進行衡量[2]。函數f(t)在某一點的Lipschitz指數α越大，該點就越光滑；否則，該點的奇異性越大。核磁共振測井自旋回波信號是連續變換的，沒有間斷點出現，其有效信號的Lipschitz指數α>0，噪聲信號的Lipschitz指數α<0。基于模極大值的小波變換去噪方法就是根據這一特性來檢測信號的局部奇異性[3]，具體方法為：

a. 求各尺度下的模極大值，其與信噪比和重構誤差均有關系。信號重構時，誤差與尺度j成正比關系，故具體的選擇需根據信號的信噪比確定。

b. 確定尺度j上的模極大值點。對信號做二進制小波分解后，選定一個閾值A和最大分解尺度j，并將j上噪聲對應的模極大值點與閾值A作比較，若其對應的幅值比閾值A小，就將其刪除，否則將其保留。

c. 尋找尺度j-1上的傳播點，這個傳播點要和尺度j上小波變換模極大值點相對應，從而保留有效信號的極值點，將其他(如噪聲引起)極值點去除。

2 測井自旋回波信號的交替投影重構

在尺度1到尺度j下，利用模極大值進行信號重構。為減小誤差，采用交替投影算法[4]恢復小波系數。

尋求某一個函數h(x)∈L2(R)，使得該函數在每一個尺度j下，其二進小波變換Wjh(x)的模極大值與Wjf(x)具有相同的模極大值。h(x)應滿足以下兩個約束條件：

a. 在各個尺度j下Wjf(x)的任意模極大值

b. {xnj}是Wjf(x)局部極大值的全體集合。

從這兩個條件可以看出，第二個條件比第一個條件實現起來要困難，因此不要求信號的局部極值點全部位于{xnj}內，而要求|Wjh(x)|的平均值達到最小。另一方面，模極大值的點數取決于Wjh(x)。因此，為了保證在第二個條件改變的情況下Wjh(x)的模極大值點數的準確性，使Wjh(x)在集合{xnj}以外的點數盡可能地少，且能量最小，引入以下范數：

(1)

對于第一個條件，設K為所有的函數序列{gj(x)}j∈Z組成的空間，則通過式(1)可得：

(2)

一個由所有函數序列{gj(x)}j∈Z組成的集合Γ，Γ?K，則有：

(3)

(4)

令PΓ和PV分別為凸集Γ和空間V上的兩個投影算子。這兩個算子分別對范數‖*‖是正

交的。令P=PV°PΓ是兩個空間V和Γ的交替投影算子。迭代n次后P(n)=P°P°…°P,可以證明，對任意空間K中的序列Y={gj(x)}j∈Z，有：

(5)

即對空間V和Γ的交替投影收斂于Hilbert空間V的正交投影。因此，經過若干次迭代后，就會收斂于所要求的函數。

3 測井自旋回波信號去噪與重構效果

對前述核磁共振測井自旋回波信號進行分解，最大層數為4，交替投影算法的迭代次數為6。經過小波模極大值處理后，得到的結果如圖2～4所示。

圖2 回波信號小波分解近似系數和細節系數

圖3 模極大值的尺度傳播特性

圖4 模極大值重構的回波信號

由圖2可以看出，信號在經過模極大值三層分解后，近似系數和原始信號較為相似，能夠反映出原始信號的特征信息。細節系數分布較緊密，表明信號噪聲較集中。

由圖3可以看出，尺度越大，小波變換的模極大值就越稀少，但還是會很好地吻合原始信號的突變點，較好地反映突變點的大小和位置，由這些保留下來的模極大值點重構的信號能夠很好地反應信號的特性，有助于提高回波信號測量的精確度。

由圖4可以看出，模極大值去噪處理的效果很明顯，探針信號中的絕大部分噪聲都被濾除。重構信號光滑、平穩，去噪后的信號與原始信號很相似。

4 結論

4.1小波模極大值去噪方法有效地保留了核磁共振測井自旋回波原始信號的奇異點信息，且該算法對噪聲的依賴性相對較低。

4.2交替投影重構信號光滑、平穩，能很好地顯示去噪后的信號與原始信號的相近性。

[1] 肖立志,謝然紅,丁娛嬌,等．核磁共振測井儀器的最新進展與未來發展方向[J].測井技術,2003，27(4)：265～269．

[2] 華春紅,任章,張敏虎.基于自適應閾值估計的模極大值去噪方法[J].航天控制,2011,29(1)：37～40,47.

[3] 段煒.基于小波變換的探地雷達信號去噪方法研究[D].長沙：中南大學,2008．

[4] 孫豐榮,翟廣濤,李艷玲,等.由小波變換模極大值重構信號的快速算法[J].小型微型計算機系統,2005,26(12):2147～2149.