斯通利波在滲透地層和裂縫帶的反射及透射:解析和數值方法模擬對比

夏飛月， 蘇遠大,3， 唐曉明,3*

1 中國石油大學(華東)地球科學與技術學院， 山東青島 266580 2 中國石油大學(華東)深層油氣重點實驗室， 山東青島 266580 3 青島海洋科學與技術國家實驗室海洋礦產資源評價與探測技術功能實驗室， 山東青島 266071

0 引言

探測地下滲透結構的滲透性及分布位置在油氣勘探開發中具有重要意義，水力壓裂縫帶滲透性的評估一直是這方面的一項重要研究課題.人們通常利用VSP(Vertical Seismic Profile，垂直地震剖面)(Minato et al., 2017)和聲波測井(Tang and Cheng, 1993)中的斯通利波進行裂縫帶和滲透地層的識別和評價.斯通利波是沿著井壁傳播的導波，在井徑不規則處(Tezuka et al., 1997)、地層彈性性質或滲透性能改變處均會發生反射及透射.Tang和Cheng(1993)將簡化的Biot-Rosenbaum(BR)理論(Tang et al., 1991a)與井筒波導的一維傳播理論結合得到的等效波數法應用于滲透地層的斯通利波反射及透射問題.一維等效波數法的中心思想和巧妙之處在于：無論裂縫帶的情況有多么復雜，只要確定聲波在裂縫地層傳播的等效波數，就可以有效地模擬裂縫地層對波傳播的影響(即波的反射和透射).這一思想在Kostek等(1998a,b)對單個和多個平板裂縫模型的模擬結果中得到了很好的驗證.簡化的BR理論能夠快速高效地模擬低頻井孔斯通利波對滲透地層的響應，在地層(Tang and Cheng, 1996; 莊春喜等, 2019)和裂縫滲透率反演(Tang et al., 1991b)有較好應用.然而，基于該理論的一維等效波數法主要考慮了地層孔隙流體流動而對骨架的影響考慮不足(唐曉明和鄭傳漢, 2004)，再加上斯通利波一維路徑傳播的近似處理，該方法的有效性需要進行評估.裂縫地層對測井聲波影響復雜，研究和分析裂縫帶的測井聲波響應對于測井數據解釋有重要意義(閻守國等, 2015).對于較為復雜的徑向和軸向都有變化的滲透結構，如低滲透地層中徑向有限長延伸的裂縫帶(Hornby et al., 1992)，一維等效波數法難以進行模擬計算.

由于所涉問題和方法的重要性，有學者利用實驗(Zhu et al., 1992)和現場數據(Tang et al., 2011)對一維等效波數法進行過驗證，但基于基礎理論和嚴格數值模擬的對比和驗證工作還不多見.Bakulin等(2005)和Alexandrov等(2007)利用一維等效波數法與基于孔彈理論的有限差分法來研究斯通利波在裸眼井和套管井滲透地層模型的反射及透射，兩種方法計算的斯通利波反射系數在0～1 kHz吻合較好.本文的目的是將前人的工作進一步深入，利用孔彈有限差分方法(Guan and Hu, 2011; Ou and Wang, 2019)和一維等效波數法(Tang and Cheng, 1993)研究不同地層厚度、滲透率和孔隙度的滲透地層對斯通利波反射及透射的影響.兩種方法計算結果驗證了有限差分法模擬結果的正確性和一維等效波數法在0～2 kHz頻段(Tang et al., 2011)內的可靠性.不同滲透率、不同結構裂縫帶的有限差分數值模擬結果有助于分析和解釋斯通利波在高滲透裂縫帶的反射及透射現象，為聲波探測方法的實際應用和分析方法適用性提供了理論依據.

1 滲透地層測井斯通利波傳播的理論和方法

1.1 Biot理論的有限差分近似方法

Biot(1956a，b，1962)對飽含流體孔隙介質建立了彈性波動理論.將井外地層作為飽含流體的孔隙介質，可以用Biot理論來求解滲透地層的聲波測井問題.Guan和Hu(2011)給出了頻域均勻各向同性孔隙介質彈性波動方程組，

(1)

(2)

(3)

(4)

其中w是流體和固體的相對位移，u是固體的位移，τ是應力張量，I是克羅內克張量，p是孔隙流體壓力.ω是角頻率，ρf和η分別是孔隙流體密度和流體黏度，ρ是地層有效密度，滿足ρ=ρs(1-φ)+ρfφ，φ是地層孔隙度.C，M和H是孔隙介質參數的組合(Biot, 1962)，

(5)

C=(1-Kb/Ks)M,

(6)

H=(1-Kb/Ks)2M+Kb+4G/3,

(7)

其中Ks是固體顆粒的體積模量，Kf是孔隙流體體積模量，Kb和G分別是巖石干燥體積模量和剪切模量.

Johnson等(1987)通過理論分析給出了描述Biot介質孔隙流體流動特征的動態滲透率κ(ω)如下：

(8)

其中ωc=φη/(α∞ρfκ0)為Biot特征頻率，κ0為靜態達西滲透率.在特征頻率以下的頻段(ω?ωc)，κ(ω)→κ0，孔隙流體的運動是靜態壓力驅動下的擴散型黏滯流動，即達西滲流；在特征頻率以上的頻段(ω?ωc)，κ(ω)→iηφ/(α∞ρfω)，流體運動是一種傳播型波動.其中α∞是孔隙彎曲率，m=φΛ2/(α∞κ0)是一個無綱量組合,其中的Λ是孔隙體面比.對于裂縫滲透地層，α∞=1，m=12；對于孔隙滲透地層，α∞=3，m=8(Tang et al., 1991a; Tang and Cheng, 1993).

由于動態滲透率中含有根號項，不便于用有限差分法將方程(1)求解.為解決這個問題，Masson等(2006)將κ(ω)用以下形式替換

(9)

將(8)式中的根號取泰勒展開的一階項便可得到(9)式，該近似滿足條件|4ω/(mωc)|<1.一般情況下，式(8)和(9)在由低到高的整個頻段內都十分符合(Masson et al., 2006; Guan et al., 2009)，因此可以將(9)式看作是為實現Biot理論有限差分方法所做的一個有效近似.由此可以將頻域中的Biot方程組變換到時域的速度-應力形式，整理得到

(10)

(11)

(12)

(13)

以上方程組的系數是線性的(非時間或頻率的函數)，很適合用有限差分方法求解.對于聲波測井問題，可以將公式(10)—(13)在柱坐標系表示(Randall et al., 1991; He et al., 2012)，將介質的模型參數在空間域內網格化，并利用交錯網格和中心差分格式對以上方程組進行離散化求解，時間和空間的差分均采用二階近似.前人對以上孔彈方程組的有限差分求解過程進行了詳細闡述(Guan and Hu, 2011; Ou and Wang, 2019)，因此本文將不再贅述.模型的吸收邊界采用不分裂的完美匹配層(Nonsplitting Perfectly Matched Layer, NPML),具體公式推導可參見Wang和Tang(2003)和Song等(2005)；對于模型中流體、彈性固體和孔隙固體共存的介質，通過設置參數仍然可以用Biot方程組((10)—(13))來表示(Guan and Hu, 2011; 閻守國等, 2015; Ou and Wang, 2019).在流體中，ρ=ρf，φ=1，G=0，C=M=H=Kf，C2=ρf，C1=0，τ=-pI，vw=0,方程(11)和(12)簡化為流體中表示的聲壓波動方程組.在彈性固體中，ρf=ρ=ρs，φ=0，G=Gs，C=Ks，H=Ks+4Gs/3，M=∞，C1=∞，C2=∞，p=0，vw=0，方程(11)和(12)退化為彈性固體中的應力速度表示的波動方程組，其中ρs，Ks和Gs分別表示彈性固體的密度、體積模量和剪切模量.因此，波動方程(10)—(13)可統一描述三種介質的聲場，不同介質銜接處可采用物質參數平均技術(Guan and Hu, 2011).

1.2 一維等效波數法

井中斯通利波具有良好的導波性質，在低頻范圍內井孔聲場是相對穩定的.因此，可以利用導波一維傳播理論(White, 1983; Hornby et al., 1989; Tang and Cheng, 1993; Tezuka et al., 1997; Bakku et al., 2013)來描述斯通利波的傳播特征，這包括了兩個重要的假定，一是斯通利波的能量主要集中在井內；二是斯通利波在滲透地層中的傳播主要受井孔與地層的流體交換影響，受地層彈性性質影響較小.將一維傳播理論用于非均勻波導的反射和透射時，下述的等效波數理論使問題的求解變得十分簡便.對于單個滲透夾層，Tang和Cheng(1993)推導得到斯通利波的反射系數R及透射系數T如下：

(14)

(15)

其中L表示夾層厚度，k1表示夾層上下無限大地層的波數，k2表示中間夾層的波數，a1和a2是井孔截面積.對于多個夾層的情況，可以利用一維傳播矩陣計算(Tezuka et al., 1997),每一夾層中的波傳播由該層的波數控制.對于彈性地層，采用無限大地層斯通利波波數(White, 1983)；對于滲透地層，采用簡化BR理論推導的斯通利波波數(Tang et al., 1991a)，

(16)

R是井孔半徑，a是儀器半徑(本文未考慮儀器的情況，所以取a=0)，ke是無限大彈性地層的斯通利波波數(White，1983)，K1和K0分別是一階和零階變型貝塞爾函數，D是孔隙流體的擴散率(唐曉明和鄭傳漢, 2004).

將一維等效波數的解析方法和有限差分法數值方法用來研究斯通利波在滲透地層的反射及透射有以下意義：等效波數法有明確的物理意義和簡潔的數學描述，可以快捷地得到模擬結果，但該方法主要模擬地層滲透性的影響而對地層彈性的影響考慮不足(唐曉明和鄭傳漢, 2004)，再加上所做的簡化和近似，其結果的精度和有效性需要評估.孔彈有限差分方法對地層的彈性和滲透性的影響均包括在內，可對分層界面處波在井內外的能量分布和傳遞進行模擬，但相對于前一方法而言計算費工耗時.通過兩種方法計算結果的比較，可以加深對斯通利波在滲透層的反射和透射問題的物理機制的理解并相互驗證兩種方法的結果，計算結果對照的差異可以區分地層彈性性質和滲透率的影響.對于等效波數法，在驗證其有效性和適用范圍后，可以利用該方法的快速計算優點來反演模型參數(Tang et al., 1991b)；而孔彈有限差分方法更有用的功能是做復雜的模型計算，如對于前一方法不能模擬的軸向和徑向分層變化的復雜情況.

1.3 有限差分法對均勻滲透地層的驗證

數值算法的準確性是衡量時域有限差分方法(Finite-Difference Time-Domain method,FDTD)的重要標準，對于裸眼井均勻滲透地層模型，可用實軸積分法(Real-Axis Integral,RAI)(Schmitt,1988; Chen et al., 2014)計算結果對數值方法進行驗證.表1給出不同地層的流體和彈性參數，滲透地層的干燥體積模量和剪切模量可以用縱橫波速度計算(Norris, 1989; Tang et al., 1991a).為了驗證有限差分在中高滲透率范圍的可行性，計算的裸眼井地層模型分別為孔隙滲透地層和裂縫滲透地層.對于孔隙滲透地層，孔隙度及滲透率分別為φ=0.25和κ0≈10-12m2，一般砂巖地層孔隙曲率α∞=3，無量綱常數m=8，由此可得κ(ω)中的Biot特征頻率ωc=50 kHz.對于裂縫滲透地層，裂縫孔隙度和滲透率與地層所含裂縫開度總和相關(Tang and Cheng, 1993)，定義式如下所示：

表1 彈性和滲透地層參數表

(17)

該式表示厚度為L的地層含有n條開度為L0的裂縫，一般情況下裂縫開度范圍為0.001 mm到1 mm(Hornby et al., 1992).裂縫滲透地層孔隙曲率α∞=1，表示流體在直的裂縫通道內流動；由條件Λ=L0得到無量綱常數m=12.Tang等(1991a)理論證明通過以上參數設置，裂縫導通率(Tang and Cheng, 1989)可以用動態達西滲透率(式(8))來表示，因此該參數下孔隙滲透介質可以表征多條裂縫的斯通利波響應(Tang and Cheng, 1993; Kostek et al., 1998b; Minato et al., 2017).對于人工壓裂區的高滲透裂縫帶，取平均裂縫寬度L0=0.049 mm，裂縫孔隙度φ=0.25，計算得到平均滲透率κ0≈5×10-11m2，此時Biot特征頻率ωc=5 kHz.

(18)

f0是聲源主頻，Tc代表聲源脈沖寬度.為激發井孔中的低頻斯通利波，取聲源主頻f0=1 kHz，聲源脈寬Tc=2 ms.

圖1a和1b分別是孔隙和裂縫滲透地層陣列聲壓波形圖，紅色虛線為有限差分(FDTD)解，黑色實線為實軸積分(RAI)解，兩種方法計算的陣列波形以第一道波幅最大值做歸一化.聲源到第一個接收器的距離為0.5 m，接收器間距為0.5 m(圖1a)和0.2 m(圖1b)，裂縫滲透地層具有較高的滲透率，波幅衰減很大，因此縮短接收間距以顯示波形特征.對于孔隙和裂縫滲透地層模型，有限差分數值解和解析解計算的波形相位和幅度符合十分好.盡管為實現有限差分解做了式(8)的近似，中、高滲透率地層的斯通利波傳播模擬仍具有較高的計算精度.因此，有限差分方法提供了較為精確的Biot理論數值解，可以作為驗證一維等效波數法的理論參照.

圖1 兩種方法(紅色虛線：有限差分；黑色實線：實軸積分)計算的滲透地層中沿井壁傳播的斯通利波陣列聲壓波形圖

2 對復雜滲透地層模型的應用

由于在實際測井數據處理中主要提取0～2 kHz的斯通利波數據進行滲透結構的識別和評價(Tang et al., 2011)，本文主要研究0～2 kHz測井斯通利波在復雜滲透地層的反射及透射.首先利用有限差分方法和一維等效波數法模擬斯通利波在單個和多個滲透夾層的反射及透射，考察等效波數法在各條件下的適用性.對于等效波數法難以模擬的復雜滲透模型，如井孔周圍徑向和軸向分布的裂縫帶，可采用有限差分方法進行模擬.一維等效波數法通過公式(14)和(15)直接給出頻域中單個滲透夾層的斯通利波反射和透射系數；有限差分得到的是時域波形，需要用反射、透射和直達(即入射)斯通利波傅里葉變換后的振幅譜來計算這些系數，如(19)式：

(19)

式中的RWV、TWV和DWV分別為反射、透射和直達斯通利波的振幅譜.

2.1 斯通利波在孔隙滲透地層的反射及透射

首先考慮Tang和Cheng(1993)文章中的單個滲透夾層模型：彈性地層中含有1 m厚孔隙滲透地層(表1)，井孔半徑為0.13 m.改變孔隙滲透地層的孔隙度和滲透率，由低孔低滲(φ=0.08，κ≈10-15m2)到高孔高滲(φ=0.35，κ≈2×10-12m2)變化.圖2為兩種方法計算的斯通利波反射及透射系數，有限差分法(實線)和等效波數法(虛線)計算結果在0～2 kHz都十分吻合，但在大于1 kHz的較高頻段出現一些差異.這些差異隨著孔隙度和滲透率升高逐漸變小，其變化趨勢可以在斯通利波透射系數圖(圖2b)明顯看到.當地層孔隙度和滲透率增加，斯通利波反射系數在整個頻段內增加，低頻(<0.8 kHz)增加較為明顯；透射系數在整個頻段內降低，高頻(2 kHz)降低較為明顯.圖2所示的低頻反射增強和高頻透射減少現象反映的是斯通利波在滲透地層的波速降低(主要發生在低頻)和波幅衰減(隨頻率增加)的物理本質，在兩種方法的模擬結果中都得到很好的體現.反射系數隨頻率變化的峰、谷現象是夾層上、下界面反射波相互干涉造成的，稍后有進一步討論.

考慮圖3所示的多個滲透夾層的情況，在無限大彈性地層中，厚度為2.4 m的低孔滲地層(φ1=0.08，κ1≈10-15m2)含有一定厚度的高孔滲夾層(φ2=0.35，κ2≈2×10-12m2)，井孔半徑為0.13 m.圖4為高孔滲夾層在不同厚度(0.4 m, 0.8 m, 1.0 m)下兩種方法計算得到的斯通利波反射及透射系數，有限差分法(實線)和等效波數法(虛線)的結果吻合良好.值得注意的是等效波數法僅考慮井壁上(忽略滲透夾層之間)流體的交換(Tang and Cheng, 1993)，但計算結果仍與考慮了夾層間流體交換的有限差分結果十分一致.隨著高孔滲夾層厚度增加，斯通利波透射系數在整個頻段內單調降低，反射系數的峰值向低頻移動，這是由于高孔滲夾層厚度變化引起的上、下反射波相互干涉造成的.

接下來討論圖3模型中高孔滲夾層的滲透性變化對斯通利波反射及透射的影響，保持高孔滲夾層厚度為1 m不變，改變地層的孔滲參數(φ=0.15和κ≈2×10-13m2,φ=0.25和κ≈10-12m2，φ=0.35和κ≈2×10-12m2)，其他參數與圖4所用參數一致.圖5是高滲透夾層在3種不同孔滲參數下兩種方法計算的斯通利波反射(a)及透射(b)系數，有限差分法(實線)和等效波數法(虛線)的結果吻合較好，但2 kHz附近存在較大差異.隨著中間高孔滲夾層的滲透性增強，該差異有所減小，此規律與圖2相同.并且與圖2單個滲透性夾層對比，圖5中高孔滲夾層在相同參數(φ=0.35和κ≈2×10-12m2)下的斯通利波反射及透射系數曲線形態十分相似，說明在低孔滲背景下，斯通利波的反射及透射由高孔滲地層控制.

圖2 4種孔滲參數下兩種方法計算的單個孔隙滲透夾層(Tang and Cheng, 1993)的斯通利波反射(a)及透射(b)系數圖實線是有限差分計算結果；虛線是一維等效波數法計算結果.

圖3 多個滲透夾層模型.在無限大彈性地層中，2.4 m厚的低孔滲地層(φ1=0.08，κ1≈10-15m2)含有一定厚度的高孔滲(φ2=0.35，κ2≈2×10-12m2)夾層

圖4 圖3模型中高孔滲夾層在3種厚度(0.4 m, 0.8 m, 1.0 m)下的斯通利波反射(a)及透射(b)系數圖

圖5 圖3模型中高孔滲夾層在3種不同孔滲參數下的斯通利波反射(a)及透射(b)系數圖

2.2 斯通利波在高滲透裂縫帶的反射及透射

斯通利波分析在裂縫及裂縫帶的識別和評價有著重要應用(Tang and Cheng, 1993; Tang et al., 2011).天然和人造裂縫往往會在井孔附近形成軸向和徑向延伸的裂縫帶(Hornby et al., 1992; Medlin and Schmitt, 1994),裂縫帶的存在大大增加了井壁附近的滲透性，使得斯通利波在裂縫帶產生明顯的衰減和較強的反射.對于軸向和徑向都有變化的裂縫帶，我們用有限差分數值方法來研究斯通利波的傳播特征.地層模型如圖6所示，位于無限大彈性地層中厚度為2.4 m的低孔滲地層含有徑向和軸向有限延伸的裂縫帶.彈性地層和低孔滲地層骨架及流體參數與圖3一致，裂縫帶骨架參數為表1滲透地層參數，裂縫孔隙度φ=0.25，孔隙曲率α∞=1，m=12.

圖6 高滲透裂縫帶模型.位于無限大彈性地層中厚度為2.4 m的低孔滲地層(φ1=0.08，κ1≈10-15m2)含有軸向延伸長度為L，徑向延伸長度為w的高孔滲裂縫帶

圖7是裂縫帶在不同滲透率下有限差分計算的斯通利波反射(a)及透射(b)系數.裂縫帶徑向延伸長度w=0.2 m，軸向延伸長度L=1 m，滲透率分別取8×10-12m2、2×10-11m2和5×10-11m2，對應的Biot特征頻率分別為31.25 kHz、12.5 kHz和5 kHz.當裂縫帶滲透率由8×10-12m2變化到2×10-11m2時，斯通利波反射(透射)系數在高頻范圍(大于0.6 kHz)增加(降低)顯著.當滲透率由2×10-11m2變化到5×10-11m2時，斯通利波反射系數在整個頻段處于高值且0～1 kHz和1～2 kHz的差異變小，波峰與波谷較為清晰；透射系數曲線變化較為復雜，在低頻段降低后隨頻率增加稍有升高.由于高滲透率條件下裂縫帶流體滲流接近Biot特征頻率，由黏滯流動向傳播波動過渡(Johnson et al., 1987)，因此波的透射系數隨頻率的變化并不是單調降低.

圖7 圖6模型中裂縫帶在不同滲透率(8×10-12m2，2×10-11m2，5×10-11m2)下的斯通利波的反射(a)及透射(b)系數圖

考察裂縫帶徑向延伸長度的影響.保持圖6模型參數不變，裂縫帶滲透率κ≈2×10-11m2和軸向延伸長度L=1 m，計算裂縫帶徑向延伸長度w在不同取值(0.05 m, 0.1 m, 0.3 m, 0.5 m, 無限長)下的斯通利波響應.圖8為裂縫帶的斯通利波反射及透射系數計算結果.當0.05 m≤w≤0.3 m時，斯通利波反射系數先在整個頻段內明顯升高，0.1 m≤w≤0.3 m時僅在低頻(<1 kHz)升高較為明顯，透射系數在整個頻段內大幅降低.當w>0.3 m時，斯通利波反射系數在低頻較小幅度上升，透射系數在低頻降低高頻稍有升高，最后反射和透射系數趨近于無限長延伸裂縫帶的曲線.模擬結果表明斯通利波激發的孔隙流體滲流的穿透深度有限，是一種與頻率有關的趨膚效應(Zhao et al., 1993)，裂縫帶徑向長度在0.1 m內增加時，斯通利波反射系數在0～2 kHz頻段內均升高，反射系數僅低頻部分(<1 kHz)對徑向長度大于0.1 m的裂縫帶有響應；斯通利波透射系數則在全頻段(0～2 kHz)對0.3 m內徑向長度變化的裂縫帶響應靈敏.

最后，考察圖6模型中裂縫帶軸向延伸長度變化對斯通利波的影響.固定裂縫帶徑向長度w=0.2 m，其他參數與圖8一致，計算軸向長度L不同取值(0.4 m, 0.8 m, 1.0 m, 1.5 m)下斯通利波響應，結果如圖9所示.裂縫帶軸向變化引起的斯通利波反射及透射系數變化規律與滲透地層類似(圖4)；斯通利波反射系數峰/谷值及其變化與裂縫帶L值變化引起的上、下層界面反射波干涉有關；透射系數隨L增加在整個頻段單調降低.這是因為無論裂縫帶的徑向長度如何，其產生的波幅衰減隨波在衰減帶的傳播距離增加而增加.

圖8 圖6模型中裂縫帶在不同徑向延伸長度(0.05 m，0.1 m，0.3 m，0.5 m，無限長)下的斯通利波的反射(a)及透射(b)系數圖

圖9 圖6模型中裂縫帶在不同軸向延伸長度(0.4 m，0.8 m，1.0 m，1.5 m)下的斯通利波的反射(a)及透射(b)系數圖

3 結論

本文利用Biot理論研究了測井斯通利波通過滲透性地層時的傳播特征.通過有限差分和一維等效波數法對滲透地層的數值模擬結果對比，分析和驗證了一維等效波數法在滲透地層模擬的有效性.對于等效波數法難以模擬的復雜裂縫帶模型，有限差分方法給出了不同滲透率、軸向和徑向延伸長度下裂縫帶的斯通利波響應.主要結論如下：

(1) Biot理論的有限差分數值方法對滲透地層的測井斯通利波模擬具有較高精度；通過與一維等效波數法的計算結果對比，驗證了后者在0～2 kHz范圍內模擬結果的可靠性，為等效波數法的實際應用提供理論依據.

(2) 斯通利波的反射及透射由高滲透地層和裂縫帶控制.波在滲透地層的波速降低(主要發生在低頻)和波幅衰減(隨頻率增加)使得斯通利波反射系數在低頻段升高，透射系數隨頻率單調降低.斯通利波反射系數峰/谷值隨地層上、下層界面距離增加向低頻移動，透射系數則在整個頻段內降低.

(3) 斯通利波在滲透地層激發的流體滲流是一種穿透深度有限的趨膚效應，要達到0.1 m以上的探測深度，需要使用較低的激發頻率(0～2 kHz)，此頻段內的斯通利波反射和透射系數對井壁附近裂縫帶的滲透性能均有響應.隨著滲透率升高，Biot特征頻率降低，流體滲流具有波動特征，斯通利波反射系數，特別是透射系數呈較復雜變化趨勢.

本文的方法和結果為利用測井斯通利波數據分析和識別裂縫型滲透地層以及反演滲透性能提供了理論基礎，有助于水力壓裂的高滲透地層和裂縫帶的識別和評價.