王祝文，徐方慧，劉菁華，寧琴琴，于 洋

吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院，長春 130026

0 引言

聲波測井的全波列資料可以用來評價(jià)地層性質(zhì)[1]。隨著人們對孔隙地層中油氣的不斷開發(fā)，非常規(guī)油氣藏的地位開始變得越來越重要。非常規(guī)儲層的主要特征是低孔隙度、低滲透率，并且裂隙比較發(fā)育，從而這類地層中的體波和模式波會有頻散和衰減現(xiàn)象[2]。因此，研究彈性波對孔隙、裂隙介質(zhì)及油氣的響應(yīng)特征顯得十分重要。唐曉明等[3-4]利用含孔隙、裂隙介質(zhì)彈性波動的統(tǒng)一理論研究了孔隙、裂隙地層的聲波測井響應(yīng)特征，結(jié)果表明隨著裂隙密度的增加，縱、橫波速度出現(xiàn)頻散現(xiàn)象，斯通利波、偽瑞利波以及彎曲波的能量衰減和速度頻散受裂隙密度的影響，裂隙對聲場影響的頻率區(qū)間由裂隙縱橫比控制。Hornby等[5]用平板狀裂縫模型研究了斯通利波通過單一裂縫時(shí)的傳播特性，結(jié)果表明裂縫的存在會產(chǎn)生反射斯通利波。隨后Tang等[6]用滲透帶理論模型研究了斯通利波的傳播規(guī)律，指出在滲透帶地層同樣會產(chǎn)生反射斯通利波。由于斯通利波對滲透性地層或者開口裂縫十分敏感，所以斯通利波的相關(guān)信息能較好地確定裂縫的位置，評價(jià)裂縫帶的滲透性，而這種滲透性地層往往是大量油氣的儲存地點(diǎn)，因此斯通利波對地層的響應(yīng)特征在很大程度上可以指示油氣。潘保芝等[7]利用常規(guī)測井資料和電成像測井?dāng)?shù)據(jù)建立了一種識別火成巖儲層有效裂縫的方法，可快速識別儲層裂縫。黃文新等[8]進(jìn)一步證明了反射斯通利波和反射率能夠用于評估儲層滲透率，但未曾使用反射系數(shù)等參數(shù)對實(shí)際井段的陣列聲波測井?dāng)?shù)據(jù)進(jìn)行系統(tǒng)評價(jià)。

我們對信號進(jìn)行傅里葉變換可以得到聲波信號在不同性質(zhì)地層上的頻率分布特征，傅里葉變換適用于具有周期性或平穩(wěn)性的信號，所以在處理具有時(shí)變特性的陣列聲波測井信號時(shí)不免會有一定的局限性。對于非平穩(wěn)的陣列聲波測井信號，我們不僅需要從總體上了解它的頻率成分，還需要了解每一時(shí)刻信號中所包含的頻率成分。希爾伯特黃變換(Hilbert-Huang transform, HHT)可以將信號分解為有限個固有模態(tài)函數(shù)(intrinsic mode function, IMF)，然后對各模態(tài)分量進(jìn)行變換從而得到信號能量在時(shí)間尺度上的分布規(guī)律，實(shí)現(xiàn)信號動力特性的量化提取[9-12]。熊學(xué)軍等[13]利用經(jīng)驗(yàn)?zāi)B(tài)分解(empirical mode decomposition, EMD)方法對海浪觀測資料進(jìn)行了處理，結(jié)果表明，各模態(tài)在HHT譜中的分布趨勢和傅里葉譜中譜線的變化趨勢是一致的，第一模態(tài)的中心頻率與傅里葉頻率譜的譜峰頻率對應(yīng)，EMD方法是一種對非線性、非平穩(wěn)數(shù)據(jù)處理的有效方法。王祝文等[14-15]基于HHT對陣列聲波測井信號進(jìn)行了時(shí)頻分析，從HHT幅值譜中得到了邊際譜、能量譜和瞬時(shí)能量譜，認(rèn)為HHT譜表現(xiàn)出的相關(guān)特征與實(shí)際波列中所包含的縱波、橫波、斯通利波和偽瑞利波有關(guān)，不同HHT譜中的峰值與儲集層巖石結(jié)構(gòu)構(gòu)造可能存在內(nèi)在聯(lián)系；但是還沒有對它們之間的具體關(guān)系進(jìn)行更深入地分析。

前人大多是對不同結(jié)構(gòu)構(gòu)造的地層進(jìn)行理論上的數(shù)值模擬，得到的結(jié)論和規(guī)律都是理想化的，在實(shí)際工作中由于地層性質(zhì)的復(fù)雜多變往往會導(dǎo)致最終的結(jié)果與理論分析存在一定的偏差。為了加強(qiáng)聲波測井理論與實(shí)際測井?dāng)?shù)據(jù)之間的聯(lián)系，從理論分析得到對實(shí)際問題實(shí)質(zhì)的理解，進(jìn)一步弄清楚斯通利波在含氣裂縫性地層的響應(yīng)特征和HHT與地層性質(zhì)的關(guān)系，筆者在前人研究結(jié)果的基礎(chǔ)上，將中值濾波、傅里葉變換和HHT應(yīng)用在遼河盆地東部凹陷含氣孔、裂隙地層的實(shí)測陣列聲波測井?dāng)?shù)據(jù)上，研究了裂隙和油氣對斯通利波及其反射波和高、低頻波的影響，同時(shí)分析了HHT處理波列數(shù)據(jù)的優(yōu)點(diǎn)和孔、裂隙地層陣列聲波信號的時(shí)頻分布特征。

1 斯通利波波場分離

斯通利波是沿井身在流體中傳播的壓力波。在全波列記錄中，斯通利波主要存在于5 kHz以下的低頻部分。我們可以用頻率域低通濾波器得到所需的低頻斯通利波波形；再經(jīng)過中值濾波，把斯通利波的直達(dá)波與反射波分離，就可以利用反射波和直達(dá)波進(jìn)行反射系數(shù)、反射能量等參數(shù)的求取。

求取反射系數(shù)的公式[16]為

(1)

式中：D(k)是直達(dá)斯通利波的頻譜；R(k)是反射斯通利波的頻譜。如果D(k)接近0，那么r(k)便會出現(xiàn)極值；為了避免這種情況，可以用下面的公式來代替：

(2)

式中：D*(k)為D(k)的共軛頻譜；E為D(k)D*(k)的峰值。

對直達(dá)斯通利波和反射斯通利波進(jìn)行傅里葉變換，得到頻率譜，把每個深度各頻率幅值的平方相加，就求得了反射斯通利波和直達(dá)斯通利波的能量曲線。

2 希爾伯特黃變換原理

2.1 EMD原理

EMD方法的思路是用波動上、下包絡(luò)的平均值去確定“瞬時(shí)平衡位置”，進(jìn)而提取內(nèi)在模函數(shù)IMF。將原始信號x(t)分解成如下形式[14-15]：

(3)

式中：cj(t)為第j個固有模態(tài)函數(shù)；rn(t)為殘量，為一個單調(diào)函數(shù)或小于某個預(yù)定值的數(shù)。

2.2 HHT和HHT譜計(jì)算

對于任一固有模態(tài)函數(shù)c(t)，其HHT結(jié)果g(t)定義為

(4)

式中：P為柯西主值；t為時(shí)間變量。

c(t)的解析信號z(t)可定義為

z(t)=c(t)+ig(t)=a(t)eiθ(t)。

(5)

式中，a(t)和θ(t)分別為陣列聲波信號x(t)的瞬時(shí)振幅和瞬時(shí)相位，按下式計(jì)算：

(6)

(7)

由瞬時(shí)相位可得到陣列聲波信號的瞬時(shí)頻率ω(t)：

(8)

可見，經(jīng)過HHT得到的振幅和頻率都是時(shí)間的單值函數(shù)，它突破了傅里葉變換固定幅值、固定頻率的限制和整個信號長度全局量的限制。如果把振幅顯示在頻率-時(shí)間平面上，就可以得到HHT幅值譜H(ω,t)。

(9)

3 含油氣孔、裂隙地層斯通利波響應(yīng)

本文采用的是遼河盆地東部凹陷HXXX井單極子陣列聲波數(shù)據(jù)。由實(shí)測資料可知3 800～3 870 m為不含油氣地層，3 370～3 390 m為工業(yè)油氣流井段。圖1原始波列中3 812～3 822 m層段的波列幅值明顯高于3 372～3 386 m層段和3 860～3 868 m層段，3 860～3 868 m層段有“V”字型反射斯通利波。對原始波列低通濾波后得到了斯通利波，從圖1中可以看出：3 812～3 822 m斯通利波幅值穩(wěn)定，基本無反射斯通利波，為致密地層；3 860～3 868 m斯通利波的幅值明顯降低，反射斯通利波“V”字型明顯，常規(guī)測井井徑曲線顯示非擴(kuò)徑，說明可能有裂隙發(fā)育。圖1中的電成像測井圖(圖1c)證明了該層為裂隙地層。

3 372～3 386 m層段斯通利波幅值降低，特別是3 372～3 380 m；但是3 372～3 380 m沒有明顯的反射斯通利波，3 380～3 386 m出現(xiàn)了反射斯通利波，而斯通利波幅值降低不明顯。為了解釋這一現(xiàn)象，我們給出了圖1右側(cè)的3 372～3 380 m層段和3 380～ 3 386 m層段中部分電成像測井圖。

英寸(in)為非法定計(jì)量單位，1 in=2.54 cm，下同。圖1 斯通利波波場分離Fig.1 Wavefield separation of Stoneley wave

1)電成像圖1a(3 377～3 380 m)中裂隙非常發(fā)育，屬于網(wǎng)狀裂隙，該段地層被解釋為含油氣的孔、裂隙地層。該層段沒有出現(xiàn)“V”字型反射斯通利波，究其原因可能有2個：1)裂隙的高度發(fā)育致使斯通利波能量衰減非常嚴(yán)重，從而使反射斯通利波的顯示不明顯，甚至完全消失；2)網(wǎng)狀裂縫中不同裂縫(波阻抗)界面產(chǎn)生的反射斯通利波難以有效地區(qū)分，并且非同相位反射斯通利波的疊加會使信號強(qiáng)度減弱[17]。

2)電成像圖1b(3 382～3 385 m)中裂隙不明顯，而薄互層地層、巖性界面較明顯。薄互層地層、巖性界面也會產(chǎn)生反射斯通利波，但對斯通利波能量的影響不大。由于該層段為滲透性孔隙地層，斯通利波能量也會發(fā)生衰減并產(chǎn)生相應(yīng)的反射斯通利波，所以該層段可見到“V”字型反射波并且斯通利波幅值下降。這也說明在實(shí)際應(yīng)用中應(yīng)注意裂縫與薄互層地層、巖性界面的區(qū)別。

圖2中3 860～3 868 m的裂隙地層相較于3 812～3 822 m的致密地層，3 372～3 380 m的含油氣裂隙地層相較于3 380～3 386 m的含油氣孔隙地層，它們的斯通利波能量衰減明顯。這是由于斯通利波是沿井身在流體中傳播的壓力波，當(dāng)井眼穿過裂隙時(shí)井內(nèi)流體進(jìn)入裂隙，流體相對于骨架振動，在這些地層中產(chǎn)生黏滯擴(kuò)散，使斯通利波的能量明顯衰減。

如果以2 kHz頻率為分界點(diǎn)，我們可以得到“直達(dá)的”高頻斯通利波(2～4 kHz)和低頻斯通利波(0～2 kHz)的能量曲線，其數(shù)值代表不同地層斯通利波能量的相對大小(圖2)。在3 812～3 822 m致密地層區(qū)間，斯通利波的高頻能量明顯大于低頻能量；3 860～3 868 m裂隙地層高頻能量發(fā)生了明顯衰減，特別是裂隙密集段3 860～3 863 m，但低頻斯通利波能量幾乎沒有變化。此現(xiàn)象說明裂隙地層中斯通利波能量的衰減主要發(fā)生在高頻段，故會出現(xiàn)主頻下降的現(xiàn)象。圖2中的能量比是高頻斯通利波能量與低頻斯通利波能量的比值。從能量比中不僅可以看出高頻、低頻斯通利波能量的相對大小，還可以得到高、低頻斯通利波能量衰減幅度的信息：3 812～3 822 m致密地層高頻斯通利波能量大約是低頻斯通利波能量的5倍，3 860～3 868 m裂隙地層能量比的數(shù)值明顯降低且?guī)缀跣∮?；說明高頻斯通利波的能量衰減嚴(yán)重且已低于低頻，可見能量比更直觀地反映出了裂隙對于高、低頻斯通利波能量的影響。

英尺(ft)為非法定計(jì)量單位，1 ft=0.304 8 m，下同。圖2 斯通利波綜合圖Fig.2 Comprehensive achievement of Stoneley wave

同樣，裂隙對地層性質(zhì)的改變會使斯通利波產(chǎn)生相對時(shí)滯。這是因?yàn)榫兴雇ɡㄓ陕暶}沖與井壁接觸時(shí)產(chǎn)生，其同時(shí)受到彈性分界面兩側(cè)介質(zhì)性質(zhì)的影響，地層裂隙的出現(xiàn)改變了地層的性質(zhì)，降低了地層的體積模量、剪切模量等參數(shù)，致使斯通利波的速度下降，從而使波至?xí)r間推遲。因此，斯通利波在滲透性地層(裂縫)的能量衰減、主頻下降和到時(shí)滯后三者之間具有相關(guān)性，這種現(xiàn)象在3 372～3 380 m含油氣裂隙地層尤其明顯。

通過上文的討論我們知道，3 812～3 822 m為無油氣較致密地層，3 380～3 386 m為含油氣孔隙儲層，基本可以認(rèn)為油氣是這兩個層段的唯一差異。圖2含油氣儲層(3 380～3 386 m)斯通利波的能量相較于不含氣地層(3 812～3 822 m)來說是減小的，證明了油氣確實(shí)會導(dǎo)致斯通利波能量的衰減。同理，我們也可以得到油氣會使斯通利波主頻下降、到時(shí)滯后的結(jié)論。這是由于油氣的殘留會使孔隙流體的流通性和可壓縮性變大，斯通利波的衰減和走時(shí)滯后現(xiàn)象增強(qiáng)，所以油氣和裂隙對斯通利波的影響是相似的，并且裂隙對斯通利波的影響比油氣更大一些。從能量比的角度思考，相對于3 812～3 822 m致密地層，3 372～3 380 m含油氣裂隙地層的能量比幾乎為零，說明高頻斯通利波的能量衰減相當(dāng)嚴(yán)重(裂隙占主要貢獻(xiàn)，也有油氣的小部分貢獻(xiàn))；3 380～3 386 m含油氣孔隙地層的能量比雖然較小，但是其值仍然大于1，說明在同樣有油氣貢獻(xiàn)的情況下，滲透性孔隙地層高頻斯通利波能量的衰減程度遠(yuǎn)小于裂隙地層。即能量比接近0時(shí)，判斷地層為滲透性裂隙地層；能量比較小但有一定數(shù)值(略大于1)時(shí)，判斷地層為滲透性孔隙地層。能量比能反映出高頻、低頻斯通利波能量各自的衰減程度，故可從另一個側(cè)面衡量滲透性地層的性質(zhì)。

通過對斯通利波進(jìn)行波場分離，我們可以得到直達(dá)波、反射波的一些特征。由式(2)可知，反射系數(shù)為反射斯通利波和直達(dá)斯通利波主頻峰值的比值，可選取直達(dá)斯通利波主頻峰值位置進(jìn)行計(jì)算，其值與反射斯通利波能量成正比，與直達(dá)斯通利波能量成反比。圖3中：3 812～3 822 m致密地層的直達(dá)斯通利波能量比較穩(wěn)定，反射能量很小，故反射系數(shù)幾乎為零；3 860～3 868 m裂隙地層反射斯通利波能量增強(qiáng)(低頻部分增強(qiáng)明顯)，直達(dá)斯通利波能量降低，由式(2)可知該層反射系數(shù)變大并出現(xiàn)了峰值，因此反射系數(shù)表現(xiàn)為尖銳的峰值一般指示裂縫； 3 372～3 380 m層段是典型的含油氣孔、裂隙地層，由前文知3 372～3 380 m層段沒有出現(xiàn)“V”字型反射斯通利波，原因是高度發(fā)育的網(wǎng)狀裂縫使非同相位的反射斯通利波相互疊加或者斯通利波能量的嚴(yán)重衰減致使反射波的能量較低，式(2)中分母的衰減程度(直達(dá)波的能量)比分子(反射波的能量)要大得多，所以盡管沒有“V”字型反射斯通利波的出現(xiàn)，該層段的反射系數(shù)還是很大；3 380～3 386 m為含油氣孔隙地層，油氣層會產(chǎn)生一定的反射波(包括薄互層地層和巖性界面產(chǎn)生的)，最重要的是孔隙(滲透性)地層會使斯通利波能量衰減，所以一般含油氣孔隙地層也會產(chǎn)生相應(yīng)的反射系數(shù)，但其數(shù)值理論上不會很大。延伸到巖性界面和薄互層較多的致密地層，由于致密地層中的斯通利波能量幾乎不會發(fā)生衰減，即便有反射斯通利波的產(chǎn)生，反射斯通利波能量與直達(dá)斯通利波能量量級的差別過大同樣會導(dǎo)致反射系數(shù)趨近于零；說明只有裂縫產(chǎn)生的反射斯通利波才會形成高值反射系數(shù)。因此，反射系數(shù)在判斷地層的滲透性時(shí)有較高的可靠性。

4 經(jīng)驗(yàn)?zāi)B(tài)分解與希爾伯特黃變換時(shí)頻分析

4.1 波列信號的EMD和HHT

圖4是HXXX井中致密地層陣列聲波測井原始波列信號的EMD圖和HHT幅值譜。圖4中每個IMF分量都對應(yīng)著不同的振幅和頻率，從上到下逐一減小，并且每個頻率不是定值，而是圍繞一個中心頻率波動。頻率高的IMF分量對應(yīng)原始信號中先到達(dá)的縱波和橫波，頻率低的IMF分量則表示后到達(dá)的斯通利波，故IMF分量不僅包含原始信號固有的物理性質(zhì)，又具有時(shí)間上的局域化特征。將每個IMF按照式(4)—(9)進(jìn)行HHT可以得到波列信號的瞬時(shí)振幅和瞬時(shí)頻率，把各內(nèi)在模函數(shù)的振幅表現(xiàn)在頻率-時(shí)間二維圖上就可以得到三維HHT幅值譜H(ω,t)。致密地層HHT幅值譜能量主要集中在0～5 kHz，低頻能量較大，高頻能量在整個時(shí)間尺度內(nèi)均有分布，主要集中在500～2 000 μs。

HHT幅值譜直觀地顯示了原始信號的時(shí)頻分布特征。如果將HHT幅值譜對時(shí)間積分，就可以得到HHT邊際譜，邊際譜代表每個頻率的振幅在整個時(shí)間長度的累積和；將振幅的平方對時(shí)間積分可得到HHT能量譜；將振幅的平方對頻率積分可得到HHT瞬時(shí)能量譜，瞬時(shí)能量譜與圖1中原始波列給出的信息類似，它提供了陣列聲波信號能量隨時(shí)間變化的規(guī)律。

圖3 直達(dá)斯通利波和反射斯通利波對比Fig.3 Comparison of direct and reflected Stoneley wave

res.信號本身微弱的趨勢或儀器的漂零。圖4 陣列聲波信號的IMF分量(a)及HHT幅值譜(b)Fig.4 IMF components (a) and Hilbert-Huang amplitude spectrum(b) of array acoustic logging signal

由圖4中HHT幅值譜得到圖5所示的HHT邊際譜。HHT和傅里葉變換都能給出原始信號的頻率分布。由圖5可知，聲波信號的能量主要集中在0～10 kHz，并且0～5 kHz的低頻能量大于5～10 kHz的高頻能量。這是由于致密地層中各種成分的波都無明顯衰減，并且高頻縱、橫波的能量量級遠(yuǎn)小于低頻的斯通利波，高頻的偽瑞利波能量較大，所以致使低頻能量比高頻能量大，但仍保持在同一個量級上。HHT邊際譜把原始信號的頻率劃分得非常細(xì)致，得到的頻率分辨率明顯優(yōu)于傅里葉變換譜，并且HHT邊際譜對高、低頻能量的分級相對合理，即致密地層斯通利波能量較其他波動能量強(qiáng)得多；而在傅里葉變換頻率譜中得到的高頻能量較大。陣列聲波信號的一大特點(diǎn)就是無周期性，邊際譜只包含由原波動分解出的具有瞬時(shí)頻率的各個子波，從而改善了傅里葉變換的假頻現(xiàn)象，不僅提高了信號表示的效率，而且能夠表示可變的頻率，故能更好地反映出原始信號中固有的物理性質(zhì)。

圖5 信號的HHT邊際譜(a)和傅里葉頻率譜(b)Fig.5 HHT marginal spectrum (a) and Fourier transform frequency spectrum (b) of signal

圖6給出了致密地層的HHT瞬時(shí)能量譜與原始波列圖。HHT瞬時(shí)能量譜給出的信息與原始波列圖類似，并且突出了能量的集中時(shí)間段。由原始波列圖可知，斯通利波的到時(shí)在2.5 ms左右，與HHT瞬時(shí)能量譜2.5 ms能量驟增為峰值完全對應(yīng)；粗略地看出偽瑞利波在原始波列圖波至?xí)r間為1.8 ms，對比HHT瞬時(shí)能量譜發(fā)現(xiàn)與1.8～2.5 ms偽瑞利波能量正好吻合，從而驗(yàn)證了HHT瞬時(shí)能量譜的正確性；在0～1.8 ms瞬時(shí)能量譜的幅值幾乎為零，這是由于縱橫波與斯通利波能量量級相差太大而造成的。HHT瞬時(shí)能量譜包含了原始波列的一切信息并著重給出了不同時(shí)間段的能量信息，還突出了斯通利波、偽瑞利波和縱波、橫波的能量分級。由于斯通利波對地層性質(zhì)的變化非常敏感，故HHT瞬時(shí)能量譜能更好地揭示井中聲波的傳播規(guī)律，放大儲集層結(jié)構(gòu)的響應(yīng)特征，特別是孔、裂隙滲透性地層。從上面的敘述中可以看出，HHT譜從不同的方面突出了陣列聲波的動力特性，它們與斯通利波和縱波、橫波的特征相關(guān)聯(lián)，從而我們可以根據(jù)HHT譜的特征對地層屬性進(jìn)行進(jìn)一步評價(jià)。

4.2 陣列聲波測井信號的HHT時(shí)頻特征

圖7為HXXX井陣列聲波測井信號的HHT瞬時(shí)能量譜、HHT邊際譜和HHT能量譜。3 372～3 380 m和3 860～3 868 m裂隙地層原始波列的幅值比3 810～3 822 m致密地層低。由于原始波列的時(shí)間尺度較大，在整個井段不能清晰地看出幅值衰減的具體時(shí)間，這個問題在HHT瞬時(shí)能量譜中得到了較好的解決。HHT瞬時(shí)能量譜中能量具有明顯的時(shí)間集中特征，集中在1.8～3.0 ms，主要是斯通利波和偽瑞利波的能量。由于縱、橫波與斯通利波的能量數(shù)量級相差很大，所以在圖中顯示不明顯。在3 860～3 868 m裂隙地層中，整個能量段發(fā)生嚴(yán)重衰減。在3 372～3 380 m含油氣地層中，網(wǎng)狀裂縫的存在導(dǎo)致能量段衰減程度很大。3 380～3 386 m含油氣的孔隙地層造成的斯通利波能量衰減在HHT瞬時(shí)能量譜中也顯示得比較清楚。原始波列中能量衰減不明顯的問題在HHT瞬時(shí)能量譜中得到了較好的解決，體現(xiàn)了HHT瞬時(shí)能量譜的優(yōu)點(diǎn)。

在HHT邊際譜和HHT能量譜中，致密地層的波動能量主要集中在0～5和8～12 kHz這兩個部分，其中0～5 kHz主要由斯通利波能量組成，8～12 kHz主要代表縱、橫波和高頻偽瑞利波。3 860～3 868和3 372～3 380 m可明顯觀測到由油氣或裂隙所造成的低頻能量衰減(主頻降低)，其中高頻偽瑞利波似乎也發(fā)生了嚴(yán)重衰減，而偽瑞利波與地層性質(zhì)的關(guān)系還有待研究。HHT放大了陣列聲波信號對地層性質(zhì)的響應(yīng)特征，給出的地層性質(zhì)信息非常豐富。當(dāng)?shù)貙有再|(zhì)(滲透率)發(fā)生改變時(shí)，陣列聲波信號的能量和頻率變化在HHT譜上表現(xiàn)得相當(dāng)明顯，這種直觀程度明顯高于圖2綜合多種譜線的分析。

圖6 信號的聲波信號全波列(a)和HHT瞬時(shí)能量譜(b)Fig.6 Full waveform (a) and HHT instantaneous energy spectrum (b) of signal

圖7 原始數(shù)據(jù)HHT譜Fig.7 HHT spectrum of raw data

由前文可知，斯通利波對滲透性地層的敏感度非常大，為了進(jìn)一步說明HHT更能突出斯通利波對含油氣裂隙地層的響應(yīng)特征，我們對斯通利波進(jìn)行經(jīng)驗(yàn)?zāi)B(tài)分解和希爾伯特黃變換。圖8是斯通利波HHT譜。由于濾除了偽瑞利波和縱波、橫波，斯通利波HHT譜的能量整體上小于原始波列的HHT譜中的能量。

HHT瞬時(shí)能量譜能量集中在2.5～3.0 ms，2.5 ms對應(yīng)著斯通利波的波至?xí)r間。在3 812～3 822 m致密地層中斯通利波能量很大，在3 860～3 868 m裂隙地層中斯通利波能量發(fā)生嚴(yán)重衰減，到時(shí)滯后明顯。從圖8中的原始波列圖能夠看出能量的衰減，但不能看出波至?xí)r間的變化，并且能量的衰減特征沒有在HHT瞬時(shí)能量譜中表現(xiàn)得明顯。在3 372～3 380 m含油氣裂隙地層中，網(wǎng)狀裂縫的存在導(dǎo)致斯通利波衰減程度非常大；3 380～3 386 m含油氣孔隙地層造成的斯通利波能量衰減也顯示得非常清楚。但是，若直接從圖8中的斯通利波波列圖或圖2中的斯通利波能量曲線分析，能量的衰減雖然可見卻非常有限，故HHT瞬時(shí)能量譜對于粗略地尋找孔、裂隙滲透性地層的效果可能更好。

圖8 斯通利波HHT譜圖Fig.8 HHT spectrum of Stoneley wave

在HHT邊際譜和HHT能量譜中，致密地層的斯通利波能量主要集中在0～5 kHz，這與傅里葉變換得到的頻率譜表現(xiàn)一致。在3 860～3 868 m和 3 372～3 380 m可明顯觀測到由裂隙或油氣所造成的斯通利波能量衰減，主頻下降非常明顯。3 380～3 386 m滲透性孔隙地層造成的斯通利波能量衰減和主頻下降雖不及裂隙地層程度大，但是相對圖2的斯通利波頻率譜明顯了很多。從HHT邊際譜和HHT能量譜中不僅可以得到斯通利波能量和頻率的信息，還可以清晰地看出地層性質(zhì)變化導(dǎo)致斯通利波能量、頻率的變化規(guī)律，波的能量、主頻和到時(shí)的變化特征僅僅在HHT譜上就能全部表示出來，非常簡單直觀。

5 結(jié)論

本文應(yīng)用希爾伯特黃變換對含氣孔、裂隙地層的斯通利波進(jìn)行了處理，以研究斯通利波時(shí)頻分布特征，可得到以下的結(jié)論。

1)油氣和裂隙會引起斯通利波能量衰減、主頻下降和到時(shí)滯后，三者之間具有相關(guān)性。斯通利波的能量衰減主要發(fā)生在高頻部分。油氣對斯通利波的影響程度小于裂隙對斯通利波的影響程度。

2)能量比(高頻與低頻斯通利波的比值)接近零時(shí)判斷地層為滲透性裂隙地層，能量比較小但有一定數(shù)值時(shí)判斷地層為滲透性孔隙地層。

3)裂隙會使反射斯通利波能量增強(qiáng)、反射系數(shù)變大，計(jì)算的高值反射系數(shù)一般只會出現(xiàn)在裂縫性地層，故利用反射系數(shù)評價(jià)地層性質(zhì)的可靠性更高。

4)HHT幅值譜直觀地顯示了陣列聲波測井信號的時(shí)頻分布。HHT邊際譜和HHT能量譜在裂縫性地層清晰地顯示出主頻能量的大幅度下降；HHT瞬時(shí)能量譜表現(xiàn)出裂縫性儲層聲波能量的衰減規(guī)律，同時(shí)又可觀測到到時(shí)滯后。

5)HHT方法是處理時(shí)變信號的一種比較適用的方法，它能夠從原始信號中分解出若干個具有不同瞬時(shí)頻率的IMF分量，既保證了原始信號固有的物理性質(zhì)，又具有時(shí)間上的局域化特征。