專利名稱：一種地層回波信號提取方法及裝置的制作方法
技術領域：
本發(fā)明涉及ー種井下核磁共振測井儀的地層回波信號的提取技木，尤其涉及ー種地層回波信號提取方法及裝置。
背景技術：
井下核磁共振儀器測井過程中，由于地層回波信號非常微弱，使得回波信號的提取十分困難。現(xiàn)有核磁共振測井儀中普遍使用數字相敏檢波器(DPSD)的檢波算法或者小波提取算法，對輸入信號進行幅度與相位的提取。但是目前的提取方法誤差較大，精度較低，且極易受到外界噪聲的干擾影響回波信號的提取質量。

發(fā)明內容
本發(fā)明所要解決的技術問題是克服目前提取核磁共振的地層回波信號誤差較大、精度較低的缺陷。 為了解決上述技術問題，本發(fā)明提供了一種地層回波信號提取方法，包括如下步驟采集核磁共振的回波原始信號；將所述回波原始信號與同頻率的正弦信號和余弦信號分別進行正交乘法運算，獲得所述回波原始信號的正交分量和同相分量；對所述正交分量進行小波分解和小波重構，獲得正交輸出信號；對所述同相分量進行小波分解和小波重構，獲得同相輸出信號；采用數字相敏檢波技術對所述正交輸出信號和同相輸出信號進行合成，獲得核磁共振回波信號的初始相位和包絡信號幅值。優(yōu)選地，對所述正交分量進行小波分解和小波重構的步驟，包括對所述正交分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構；對所述同相分量進行小波分解和小波重構的步驟，包括對所述同相分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構。優(yōu)選地，對所述正交分量或同相分量以db5小波基進行七層小波分解的步驟中，每ー層小波分解均采用如下表達式進行d) = d\' + a(s\1 + )；s) ニ S1' + j3yd) +；= d) + /(S'l + (I)；= S11 + SlCi^z + §2d;n ；St=i.S;-;以及dt=d;/^；
其中，第一層小波分解中，<和ぜ分別為對所述正交分量或同相分量進行奇偶分裂操作得到的奇地址數據和偶地址數據；第二至第七層小波分解中，<和^*分別為對上ー層小波分解得到的dt進行奇偶分裂操作得到的奇地址數據和偶地址數據；<+1表示<向左移動ー個時間單位后得到的數據；表示式向右平移一個時間単位后得到的數據；表示ゼ向右平移一個時間単位后得到的數據；尤2表示ゴム向左延遲ー個時間周期后得到的數據；ゴム表示<向左平移一個時間周期后得到的數據；
a , ^ ^ 2> y , 6 6 2 以及 2,為系數。優(yōu)選地，a= -I. 5861343420693648 ； ^ : = I. 0796367753628087 ； ^ 2=-0. 0529801185718856 ； y = -0. 8829110755411875 ^1 = O. 4435068520511142 ； S2 =
I.576237461483639 ;以及 I = -I. 1496043988602418。本發(fā)明提供了一種地層回波信號提取裝置，包括采集模塊，用于采集核磁共振的回波原始信號；運算模塊，用于將所述回波原始信號與同頻率的正弦信號和余弦信號分別進行正交乘法運算，獲得所述回波原始信號的正交分量和同相分量；分解重構模塊，用于對所述正交分量進行小波分解和小波重構，獲得正交輸出信號；用于對所述同相分量進行小波分解和小波重構，獲得同相輸出信號；合成模塊，用于采用數字相敏檢波技術對所述正交輸出信號和同相輸出信號進行合成，獲得核磁共振回波信號的初始相位和包絡信號幅值。優(yōu)選地，所述分解重構模塊用于對所述正交分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構；所述分解重構模塊用于對所述同相分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構。優(yōu)選地，所述分解重構模塊用于采用如下表達式進行七層小波分解中的每ー層小波分解d) = 4° + び(< + <+1)；S1t = sf + pyd\ + fi2d)_y ；=d) + y{s]+s]^)\sf = sj + + S2d;+l ；ち=#.42;以及dt=df ;其中，第一層小波分解中，<和sf分別為對所述正交分量或同相分量進行奇偶分裂操作得到的奇地址數據和偶地址數據；第二至第七層小波分解中，<和 < 分別為對上ー層小波分解得到的dt進行奇偶分裂操作得到的奇地址數據和偶地址數據；<+1表示^>向左移動ー個時間單位后得到的數據；
表示式向右平移一個時間単位后得到的數據；表示べ向右平移一個時間単位后得到的數據；<+2表示向左延遲一個時間周期后得到的數據；表示<向左平移一個時間周期后得到的數據；a , ^ ^ 2> y , 6 6 2 以及(為系數。優(yōu)選地，a= -I. 5861343420693648 ； ^ : = I. 0796367753628087 ； ^ 2=-0. 0529801185718856 ； y = -0. 8829110755411875 J1 = O. 4435068520511142 ； S2 =I. 576237461483639 ;以及 I = -I. 1496043988602418。與現(xiàn)有技術相比，本發(fā)明的實施例將提升小波算法與DPSD相敏檢波算法兩者結合在一起應用于核磁共振儀器微弱信號的提取設計，兼容了兩種算法各自的優(yōu)點。通過將 兩種算法結合設計來實現(xiàn)對核磁微弱信號的提取，即具有小波算法的優(yōu)點又同時兼?zhèn)銬PSD算法的優(yōu)點，可以高效地濾除信號中的噪聲干擾，比傳統(tǒng)核磁共振測井儀所使用的信號提取算法效率更高，提取有用信號的信噪比更大，大大提高了核磁共振測井中回波信號的提取質量。本發(fā)明的其它特征和優(yōu)點將在隨后的說明書中闡述，并且，部分地從說明書中變得顯而易見，或者通過實施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點可通過在說明書、權利要求書以及附圖中所特別指出的結構來實現(xiàn)和獲得。


附圖用來提供對本發(fā)明技術方案的進ー步理解，并且構成說明書的一部分，與本發(fā)明的實施例一起用于解釋本發(fā)明的技術方案，并不構成對本發(fā)明技術方案的限制。圖I為本發(fā)明實施例的地層回波信號提取方法的流程示意圖。圖2a為本發(fā)明實施例中對回波信號的正交分量和同相分量進行小波分解和小波重構的流程示意圖。圖2b為本發(fā)明實施例中對回波信號進行小波分解和小波重構的原理示意圖。圖3為本發(fā)明實施例中對回波信號的正交分量和同相分量進行小波分解的分解模塊的原理示意圖。圖4為本發(fā)明實施例中進行小波重構的重構模塊的原理示意圖。圖5為本發(fā)明實施例的地層回波信號提取裝置的結構示意圖。
具體實施例方式以下將結合附圖及實施例來詳細說明本發(fā)明的實施方式，借此對本發(fā)明如何應用技術手段來解決技術問題，并達成技術效果的實現(xiàn)過程能充分理解并據以實施。本發(fā)明實施例以及實施例中的各個特征在不相沖突前提下的相互結合，均在本發(fā)明的保護范圍之內。如圖I所示，本發(fā)明實施例的地層回波信號提取方法主要包括如下步驟。步驟S110，采集獲得核磁共振的回波原始信號x(n)。步驟S120，將回波原始信號x(n)與同頻率的正弦信號Sin(CociIi)和余弦信號cos(co0n)分別進行正交乘法運算，獲得回波原始信號的正交分量Q(n)和同相分量I(n)。
步驟S130，對正交分量Q(n)進行小波分解和小波重構，獲得正交輸出信號Q_out(n);對同相分量I(n)進行小波分解和小波重構，獲得同相輸出信號I_out(n)。步驟S140，采用數字相敏檢波技術對正交輸出信號Q_out(n)和同相輸出信號1_out (n)進行相位合成和幅值合成的處理，獲得核磁共振回波信號的初始相位デ(《)和包絡信號幅值y(n)。本發(fā)明的實施例中，步驟S130中，對回波信號的正交分量Q(n)和同相分量I (n)進行小波分解和小波重構的過程可以采用如圖2a所示的流程來進行。步驟S131，將回波信號的正交分量Q(n)和同相分量I (n)分別以db5小波基進行七層小波分解，各自獲得七組細節(jié)系數和ー組最底層近似系數(第七層分解后的近似系數)。步驟S132，將七組細節(jié)系數全部置為零。
步驟S133，利用該最底層近似系數以及置為零的七組細節(jié)系數進行小波重構，得到兩路原始信號しout (n)與Q_out(n)以完成對其的噪聲濾波處理。如圖2b所示，將輸入信號按照圖2b所示進行抽取后，得到的奇偶兩組序列與行列式/"7(一 f進行乘法運算，最后得出的值即為分解后的奇地址分解數據dt和偶地址分解數據
StO圖2b 中
r n …、「 I Op (A+ル)]「I 0]「1 +0],、./(J=4 ハ,^n - 式
_a(l + z ) I 0 I _ ブ(1 + z) I 0I0 i
(A)P(z~')r是對表達式A中行列式P(Z)進行轉置求倒后得出的，保證POT1)7'* P (め=I。其中，a、運” @2、Y、S2以及(等都為系數。z表示時域信號作頻域Z變換的基本變量。圖2b中的戶(ff與P(Z)的意義表示輸入的時域信號進行頻域z變換后與或者p(z)作行列式頻域乘法運算，得出的值在轉換到相對應的時域上，即得到圖3與圖4所示信號分解與重構的時域操作流程圖。小波分解過程即是完成將輸入信號的奇偶兩組序列與戶(廠1/進行行列式相乘得出兩組分解后的系數。進ー步將其輸入信號的奇偶兩組序列與POT1,行列式相乘的運算進行數學上的分步拆分化簡，最后得出了如圖3所示小波分解的步驟與相關計算公式，如圖4所示的重構過程也類似與分解過程。將分解后所得的最底層近似系數與置為零的細節(jié)系數乘以表達式P(Z)進行重構，將該步行列式乘法運算進行數學的分解化簡后，得出如圖4所示結構與相關計算公式。上述步驟S131中以db5小波基對正交分量Q (n)和同相分量I (n)進行七層小波分解中，每ー層小波分解都可以采用如圖3所示的分解模塊來進行。如圖3所示，分解模塊主要包括信號分離單元310，第一、第二、...、第八延時單元321、322、 ,328,以及第一、第二、 、第六預測更新單元 331、332、 、336。在第一層小波分解中，信號分離單兀(SPLIT) 310用于對輸入的正交分量Q(n)或同相分量I (n)進行奇偶分裂操作；按照對正交分量Q (n)或同相分景T (n)的采樣信號x(t)的地址編號的奇偶性，將采樣信號分為偶地址數據 < 和奇地址數據J
其中=x(/)',式(I)4。=沖):,ト1 式⑵表達式(I)和表達式(2)中的下腳標t為大于等于0的整數。x(t)2t表示偶地址編號的采樣信號x(t)的集合，x(t)2t+1表示奇地址編號的采樣信號x(t)的集合。偶地址數據<輸入到第一延時單元321進行第一延時處理，經過第一延時處理的偶地址數據^)輸入到第三延時單元323進行第三延時處理；奇地址數據c/ 入到第二延時単元322進行第二延時處理。經過第一延時處理的偶地址數據^>、經過第三延時處理的偶地址數據4以及經過第二延時處理的奇地址數據ゴf一起輸入到第一預測更新単元331。第一預測更新単元331，分別完成預測、更新及歸ー化操作。通過如下表達式進行計算，得到第一奇地址預測數據劣。 rf; =<+び(<+A)式⑶表達式(3)中，系數a在本發(fā)明的實施例中取值為-I. 5861343420693648 C1表示偶地址數據<向左移動ー個時間單位(超前預測)后得到的數據，比如^序列的當前值是< 序列下ー個時刻值。第一預測更新単元331的輸出信號分別輸入到第四延時單元324和第二預測更新單元332。第二預測更新單元332接入第三延時單元323的輸出、第一預測更新單元331的輸出以及第四延時單元324的輸出，根據如下表達式進行計算，得到第一偶地址預測數據= S0t + 賦 + ル/ム式(4)表達式(4)中，參數P1在本發(fā)明的實施例中取值為1.0796367753628087 ;參數 在本發(fā)明的實施例中取值為-0. 0529801185718856 表示劣向右平移ー個時間單位
后得到的數據，比如 < 序列當前值是<4序列下一個時刻的值。第二預測更新単元332的輸出信號分別輸入到第五延時單元325和第三預測更新単元333。第五延時單元325的輸出信號輸入到第六延時單元326進行第六延時處理。第三預測更新単元333接入第四延時單元324的輸出、第二預測更新単元332的輸出以及第五延時單元325的輸出，根據如下表達式進行計算，得到第二奇地址預測數據42。42 = 41 + YisI + 4—1 試(5)表達式(5)中，參數Y在本發(fā)明的實施例中取值為-0. 8829110755411875 表示4向右平移一個時間単位后得到的數據，比如く i序列對應的值是べ序列對應上一個時間時刻的值。第三預測更新単元333的輸出信號分別輸入到第七延時單元327和第四預測更新単元334。第七延時單元327的輸出信號輸入到第八延時單元328進行第八延時處理。第四預測更新単元334接入第六延時單元326的輸出、第三預測更新単元333的輸出以及第七延時單元327的輸出，根據如下表達式進行計算，得到第二偶地址預測數據e2
°sf = s) + S'dl2 + S2d~t ェ式(6)表達式(6)中，參數S1在本發(fā)明的實施例中取值為o. 4435068520511142 ;參數 S2在本發(fā)明的實施例中取值為I. 576237461483639。ゴム表示Jw向左延遲一個時間周期后得到的數據，當前4序列對應的值是Ai1序列對應下ー時刻的值；式^表示或2向左平移ー個時間周期后得到的數據，當前<+1序列對應的值是 < 序列對應下ー時刻的值。第四預測更新単元334將第二偶地址預測數據ろ2輸入到第五預測更新単元335，第五預測更新単元335根據如下表達式(7)進行計算，輸出一次小波分解后的偶地址分解數據st。第八延時單元328的輸出信號輸入到第六預測更新単元336，第六預測更新単元336根據如下表達式(8)進行計算，輸出一次小波分解后的奇地址分解數據dt。ろ=彳式(7)4=42 K式⑶表達式(7)和表達式⑶中，參數I在本發(fā)明的實施例中取值為-I. 1496043988602418。將上述的奇地址分解數據dt再次進行信號采樣形成采樣信號X (t)，輸入到圖3所示的分解模塊繼續(xù)進行分解操作，后續(xù)每ー層小波分解的輸入信號，均是上ー層小波分解得到的奇地址分解數據。直到對正交分量Q(n)和同相分量I(n)進行七層小波分解，最后得到最底層近似系數。本發(fā)明的實施例中，各個預測更新單元采用位移的方式來實現(xiàn)乘法運算，如此可以高效地降低資源量，提高運行速度，方便硬件實現(xiàn)。步驟S133的小波重構過程是步驟S131的小波分解過程的逆過程，可以完成對信號的重構操作。上述步驟S133中，采用如圖4所示的重構模塊來進行小波重構，得到兩路原始信號 I_out (n)與 Q_out(n)。經過七層分解后，得到的奇地址分解數據dt和偶地址分解數據St分別輸入到如圖4所示的第一預測更新単元411和第二預測更新単元412進行小波重構運算。此處的小波重構運算是前述小波分解運算的逆運算。在這個逆運算的處理過程中，使用前述小波分解所得到的最底層近似系數，以及置為零的各細節(jié)系數。第一預測更新単元411的輸出信號<輸入到第一延時單元421以及第三預測更新単元413，第一延時單元421輸出的信號也輸入到第三預測更新単元413。第二預測更新單元412的輸出信號#輸入到第三預測更新単元413。第一延時單元421的輸出信號輸入到第三延時單元423，第三延時單元423的輸出信號輸入到第四預測更新単元414。第三預測更新単元413產生數據^分別發(fā)送給第四延時單元424和第四預測更新単元414，第四延時單元424的輸出信號也輸入到第四預測更新單元414。第四預測更新単元414產生數據J分別發(fā)送給第五延時單元425和第五預測更新單兀415。第四延時單兀424的輸出信號以及第五延時單兀425的輸出信號也輸入到第五預測更新単元415。第五延時單元425的輸出信號還輸入到第六延時單元426。
第五預測更新単元415產生數據<分別發(fā)送給第七延時單元427和第六預測更新単元416，第四延時單元424的輸出信號也輸入到第六預測更新単元416，第七延時單元427的輸出信號在輸入到第八延時單元428的同時也輸入到第六預測更新単元416。第六預測更新単元416產生數據<以及第八延時單元428輸出的數據^一起輸入到信號合成単元430中產生合成信號y (t)，據此可以獲得對應的原始信號I_out(n)與Q_
out (nノ。本發(fā)明的實施例通過上述小波分解與重構計算處理后的正交輸出信號Q_out(n)和同相輸出信號I_out(n)，從回波原始信號中濾除了混雜在其中的噪聲信號，大大提高了輸出信號的信噪比。如圖5所示，本發(fā)明實施例的地層回波信號提取裝置，主要包括采集模塊510、運算模塊520、分解重構模塊530以及合成模塊540，其中采集模塊510，用于采集核磁共振的回波原始信號；運算模塊520，與采集模塊510相連，用于將所述回波原始信號與同頻率的正弦信號和余弦信號分別進行正交乘法運算，獲得所述回波原始信號的正交分量和同相分量；分解重構模塊530，與運算模塊520相連，用于對所述正交分量進行小波分解和小波重構，獲得正交輸出信號；用于對所述同相分量進行小波分解和小波重構，獲得同相輸出信號;合成模塊540，與分解重構模塊530相連，用于采用數字相敏檢波技術對所述正交輸出信號和同相輸出信號進行合成，獲得核磁共振回波信號的初始相位和包絡信號幅值。所述分解重構模塊530用于對所述正交分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構。所述分解重構模塊530用于對所述同相分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構。所述分解重構模塊530用于采用如下表達式進行七層小波分解中的每ー層小波分解d)=《+ a{ sf + C1)；s) = sf + p'dt + JS2Jll ；d; = d) + y(s] + )；sマ=sj + SlCi;: + S2d~+l ;St = K ;以及dt=(i； I其中，第一層小波分解中，<和^>分別為對所述正交分量或同相分量進行奇偶分裂操作得到的奇地址數據和偶地址數據；第二至第七層小波分解中，<和<分別為對上ー層小波分解得到的dt進行奇偶分裂操作得到的奇地址數據和偶地址數據；sf+1表示<向左移動ー個時間單位后得到的數據；足1表示式向右平移一個時間單位后得到的數據；
表示べ向右平移一個時間単位后得到的數據；表示向左延遲ー個時間周期后得到的數據；式2+1表示<向左平移一個時間周期后得到的數據；a、P1. P 2、Y、S2以及る為系數，本發(fā)明實施例中的取值請參考前述本發(fā)明實施例地層回波信號提取的描述。本發(fā)明的實施例結合提升式小波算法與DPSD相敏檢波算法的優(yōu)點，高效地濾除了回波中含有的噪聲信號，有效地提取其中的回波信息，更好地反映地層流體信息，將噪聲的影響減小到最低程度。本發(fā)明的實施例能夠大大提高回波信號的提取效率與精確度，為后續(xù)更好地進行地層流體分析奠定了良好的基礎。本領域的技術人員應該明白，上述的本發(fā)明實施例所提供的裝置和/或系統(tǒng)的各組成部分，以及方法中的各步驟，可以用通用的計算裝置來實現(xiàn)，它們可以集中在單個的計 算裝置上，或者分布在多個計算裝置所組成的網絡上，可選地，它們可以用計算裝置可執(zhí)行的程序代碼來實現(xiàn)，從而，可以將它們存儲在存儲裝置中由計算裝置來執(zhí)行，或者將它們分別制作成各個集成電路模塊，或者將它們中的多個模塊或步驟制作成單個集成電路模塊來實現(xiàn)。這樣，本發(fā)明不限制于任何特定的硬件和軟件結合。雖然本發(fā)明所揭露的實施方式如上，但所述的內容只是為了便于理解本發(fā)明而采用的實施方式，并非用以限定本發(fā)明。任何本發(fā)明所屬領域內的技術人員，在不脫離本發(fā)明所揭露的精神和范圍的前提下，可以在實施的形式上及細節(jié)上作任何的修改與變化，但本發(fā)明的專利保護范圍，仍須以所附的權利要求書所界定的范圍為準。
權利要求
1.一種地層回波信號提取方法，包括如下步驟 采集核磁共振的回波原始信號； 將所述回波原始信號與同頻率的正弦信號和余弦信號分別進行正交乘法運算，獲得所述回波原始信號的正交分量和同相分量； 對所述正交分量進行小波分解和小波重構，獲得正交輸出信號； 對所述同相分量進行小波分解和小波重構，獲得同相輸出信號； 采用數字相敏檢波技術對所述正交輸出信號和同相輸出信號進行合成，獲得核磁共振回波信號的初始相位和包絡信號幅值。
2.根據權利要求I所述的方法，其中 對所述正交分量進行小波分解和小波重構的步驟，包括對所述正交分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構； 對所述同相分量進行小波分解和小波重構的步驟，包括對所述同相分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構。
3.根據權利要求2所述的方法，其中 對所述正交分量或同相分量以db5小波基進行七層小波分解的步驟中，每ー層小波分解均采用如下表達式進行d)=《+ a(s]' + C1);^1 = 5/0+M1+M1-I ； =d)+y{s]+s)_,)； s； = sj + SjdfI2 + 6zdf+l ； A ろ2;以及 4=42,#; 其中，第一層小波分解中，< 和4分別為對所述正交分量或同相分量進行奇偶分裂操作得到的奇地址數據和偶地址數據；第二至第七層小波分解中，4°和ぜ分別為對上ー層小波分解得到的dt進行奇偶分裂操作得到的奇地址數據和偶地址數據； 表示<向左移動ー個時間單位后得到的數據； 表示41向右平移一個時間単位后得到的數據； 表示ろ1向右平移一個時間単位后得到的數據； 式I2表示或2+1向左延遲一個時間周期后得到的數據； 表示42向左平移一個時間周期后得到的數據； a > ^ ^ 2> y > 8 6 2 以及 2,為系數。
4.根據權利要求3所述的方法，其中 a = -I. 5861343420693648 ；^1 = L 0796367753628087 ；^2 = -O. 0529801185718856 ；Y = -0. 8829110755411875 ；S1 = O. 4435068520511142 ； S2=L 576237461483639 ；以及 I = -I. 1496043988602418。
5.一種地層回波信號提取裝置，包括 采集模塊，用于采集核磁共振的回波原始信號； 運算模塊，用于將所述回波原始信號與同頻率的正弦信號和余弦信號分別進行正交乘法運算，獲得所述回波原始信號的正交分量和同相分量； 分解重構模塊，用于對所述正交分量進行小波分解和小波重構，獲得正交輸出信號；用于對所述同相分量進行小波分解和小波重構，獲得同相輸出信號； 合成模塊，用于采用數字相敏檢波技術對所述正交輸出信號和同相輸出信號進行合成，獲得核磁共振回波信號的初始相位和包絡信號幅值。
6.根據權利要求5所述的裝置，其中 所述分解重構模塊用于對所述正交分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構； 所述分解重構模塊用于對所述同相分量以db5小波基進行七層小波分解，獲得七組細節(jié)系數和ー組最底層近似系數；將七組細節(jié)系數置為零，利用所述最底層近似系數以及置為零的七組細節(jié)系數進行小波重構。
7.根據權利要求6所述的裝置，其中 所述分解重構模塊用于采用如下表達式進行七層小波分解中的每ー層小波分解
8.根據權利要求7所述的裝置，其中 a = -I. 5861343420693648 ；^1 = L 0796367753628087 ；β 2 = -0. 0529801185718856 ； γ = -0.8829110755411875 ；δ1: = 0.4435068520511142 ；δ2 = 1. 576237461483639 ；以及 ξ = -1. 1496043988602418。
全文摘要
本發(fā)明公開了一種地層回波信號提取方法及裝置，克服目前提取核磁共振的地層回波信號誤差較大、精度較低的缺陷，該方法包括采集核磁共振的回波原始信號；將所述回波原始信號與同頻率的正弦信號和余弦信號分別進行正交乘法運算，獲得所述回波原始信號的正交分量和同相分量；對所述正交分量進行小波分解和小波重構，獲得正交輸出信號；對所述同相分量進行小波分解和小波重構，獲得同相輸出信號；采用數字相敏檢波技術對所述正交輸出信號和同相輸出信號進行合成，獲得核磁共振回波信號的初始相位和包絡信號幅值。本發(fā)明可以高效地濾除信號中的噪聲干擾。
文檔編號G01V3/38GK102830435SQ20121029234
公開日2012年12月19日 申請日期2012年8月16日 優(yōu)先權日2012年8月16日
發(fā)明者張嘉偉, 宋公仆, 薛志波 申請人:中國海洋石油總公司, 中海油田服務股份有限公司