曹鴻飛， 章成廣*， 蔡明， 陳義群， 蔡德洋

(1．長(zhǎng)江大學(xué)， 油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室， 武漢 430100; 2．中石油塔里木油田分公司勘探開發(fā)研究院， 庫(kù)爾勒 841000)

隨著中外石油天然氣工業(yè)的迅猛發(fā)展和油氣勘探的逐漸深入，常規(guī)油氣產(chǎn)量日漸減少，目前，非常規(guī)儲(chǔ)層具有越來越廣闊的勘探開發(fā)前景。中外有關(guān)致密儲(chǔ)層的研究表明，由于致密砂巖儲(chǔ)層的脆性較強(qiáng)，裂縫在致密砂巖儲(chǔ)層中普遍存在且往往比較發(fā)育，在裂縫性致密砂巖儲(chǔ)層中，裂縫不僅作為流通油氣的通道，而且有儲(chǔ)存和滲透油氣的作用，裂縫的研究、評(píng)價(jià)和預(yù)測(cè)成為了研究人員的關(guān)注重點(diǎn)[1-11]。利用現(xiàn)有的常規(guī)測(cè)井?dāng)?shù)據(jù)很難精確地評(píng)價(jià)裂縫性致密儲(chǔ)層，特殊測(cè)井資料評(píng)價(jià)裂縫的方法也各有缺陷，例如，聲、電成像測(cè)井測(cè)量范圍較窄，不能準(zhǔn)確地反映出井外裂縫的延伸。陣列聲波測(cè)井技術(shù)雖可以克服成像測(cè)井評(píng)價(jià)方法的缺陷，但由于裂縫性致密儲(chǔ)層的裂縫的復(fù)雜性和裂縫位置的不確定性，使裂縫地層對(duì)陣列聲波波形的干擾增大[12]。由于地層中的裂縫常具有高角度和粗糙性等復(fù)雜特征，利用現(xiàn)有的陣列聲波測(cè)井技術(shù)和成像測(cè)井技術(shù)對(duì)地層中的裂縫進(jìn)行定量地精細(xì)評(píng)價(jià)具有重要的實(shí)際意義[9,13]。

斯通利波是一種比井內(nèi)流體聲速稍低的低頻散導(dǎo)波，具有低頻率、高能量和高振幅的特征。當(dāng)?shù)貙又辛芽p與井孔相交時(shí)，由于地層中的泥漿會(huì)和地層中的液體互相流通，從而使斯通利波能量流失，并在非連續(xù)界面上產(chǎn)生反射和衰減[14]。在陣列聲波測(cè)井中，由于斯通利波反射系數(shù)對(duì)裂縫性地層中的裂縫寬度具有較高的靈敏度[15]，因此可以用斯通利波對(duì)裂縫進(jìn)行識(shí)別和分析。許多學(xué)者利用多種方法研究了不同裂縫因素對(duì)斯通利波傳播影響的研究。Ou等[16]通過有限差分的方法研究水平裂縫的斯通利波反射，在多個(gè)模型中模擬了斯通利波的傳播，說明了利用斯通利波檢測(cè)裂縫的可行性。楊博等[17]通過水槽實(shí)驗(yàn)得到裂縫張開度和斯通利波幅度衰減系數(shù)之間的關(guān)系，結(jié)合由陣列聲波測(cè)井資料得到的斯通利波衰減系數(shù)，計(jì)算裂縫張開度，并反演得到裂縫滲透率。劉黎等[18]利用傾斜裂縫地層模型和巖心實(shí)驗(yàn)研究了斯通利波穿過傾斜光滑裂縫時(shí)的響應(yīng)情況。Li等[19]研究了裂縫寬度、傾角和裂縫充填物等因素對(duì)斯通利波傳播的影響。

上述學(xué)者所研究的裂縫都為光滑裂縫，由于裂縫壁面的粗糙性會(huì)使測(cè)井裂縫評(píng)價(jià)結(jié)果和實(shí)際產(chǎn)生偏差，為了減小因裂縫的粗糙性造成的誤差，有必要對(duì)粗糙裂縫進(jìn)行研究。粗糙裂縫的粗糙性相對(duì)難以量化，常用的二維和三維表面粗糙度表征方法，分為統(tǒng)計(jì)，分形和方向特征[20]。中外學(xué)者用多種方法對(duì)粗糙裂縫進(jìn)行構(gòu)建和描述，Brown等[21]通過分形維數(shù)、均方根粗糙度和描述兩個(gè)裂縫面之間不匹配程度的長(zhǎng)度來描述粗糙裂縫，建立了一種巖石粗糙裂縫的數(shù)學(xué)模型。曲冠政等[22]利用三維掃描原理，計(jì)算裂縫的迂曲度、傾角、粗糙度等特征參數(shù)，模擬了頁(yè)巖粗糙裂縫中流體的滲流。李大奇等[23]通過調(diào)整分形理論產(chǎn)生粗糙裂縫面和裂縫基準(zhǔn)面之間的位置來描述粗糙裂縫的縫寬分布，從而構(gòu)建粗糙裂縫。彭嬌等[24]對(duì)裂縫結(jié)構(gòu)進(jìn)行分析，通過采用滲流速度描述裂縫結(jié)構(gòu)中流體滲流的迂曲度計(jì)算公式，用于描述裂縫的粗糙性。都凱等[25]運(yùn)用節(jié)理粗糙度系數(shù)(joint roughness coefficient，JRC)來計(jì)算粗糙裂縫滲透率。Ban等[26]用表征巖石表面局部和整體特征的兩個(gè)粗糙度參數(shù)來替代裂縫粗糙度。上述有關(guān)粗糙裂縫的研究大多集中在粗糙裂縫的描述和滲流方面，在測(cè)井裂縫評(píng)價(jià)結(jié)果中的應(yīng)用較少，由于地層中存在粗糙裂縫，研究粗糙裂縫對(duì)陣列聲波的影響有利于提高測(cè)井裂縫評(píng)價(jià)結(jié)果的精度。

1 實(shí)驗(yàn)原理

理想全波列波形包括縱波、橫波和斯通利波。首先是縱波先到達(dá)，其波速最快。其次是橫波，橫波波速低于縱波而大于斯通利波，其振幅受地層影響較大。斯通利波是一種井筒波，波速比縱、橫波速度低，頻率比縱波頻率小，井眼環(huán)境對(duì)斯通利波的振幅影響比較大：在規(guī)則的井眼環(huán)境中，斯通利波振幅遠(yuǎn)大于縱波和橫波；在擴(kuò)徑的井眼環(huán)境中，斯通利波振幅會(huì)根據(jù)擴(kuò)徑的程度的改變而變小。全波列波形經(jīng)過裂縫后的波形發(fā)生變化，當(dāng)聲波經(jīng)過裂縫時(shí)，全波列的振幅衰減非常明顯。本文實(shí)驗(yàn)主要研究斯通利波和橫波的衰減。實(shí)驗(yàn)中斯通利波是由單極源低頻激發(fā)產(chǎn)生，由單級(jí)接收器接收記錄，其頻率范圍在3～200 kHz。

在數(shù)據(jù)分析時(shí)利用衰減系數(shù)對(duì)全波列波形經(jīng)過裂縫后波形發(fā)生的變化進(jìn)行描述。通過提取橫波、斯通利波的幅度，用來計(jì)算衰減系數(shù)α[17,27-28]，計(jì)算公式為

(1)

式(1)中：A0為沒有裂縫時(shí)讀取的斯通利波幅度，mV；A為裂縫存在時(shí)讀取的斯通利波幅度，mV；l為源距，m。

2 具有不同粗糙度裂縫模型實(shí)驗(yàn)

2.1 模型井裂縫粗糙度實(shí)驗(yàn)

模型井是根據(jù)真實(shí)地層的特征加工而成，用于模擬實(shí)際儲(chǔ)集層中復(fù)雜的地質(zhì)特征[29]，能系統(tǒng)地研究井眼周圍的裂縫和地層孔隙度等地質(zhì)特征對(duì)測(cè)井結(jié)果的影響。

在實(shí)驗(yàn)室研究中，通常根據(jù)相似理論，將井孔和換能器的尺寸縮小一定倍數(shù)(5倍)，而聲源頻率提高5倍，來模擬實(shí)際井中彈性波的傳播。與聲系尺寸相比，實(shí)際井的深度可以看作是無限的，井壁地層厚度也可視為無限。模型井的井壁厚度設(shè)計(jì)主要取決于聲波的探測(cè)深度的大小。在單極子聲波測(cè)井中，縱波的探測(cè)深度與聲系源距有關(guān)，當(dāng)源距為3～4 in(1 in=25.4 mm)時(shí)，聲速探測(cè)深度為3～4 in，即10 cm左右。因此，可近似地認(rèn)為，在井壁層厚度至少為15 cm的模擬井中，在短源距的條件下，可以測(cè)得斯通利波的傳播速度。

根據(jù)以上分析，實(shí)驗(yàn)設(shè)定模型井長(zhǎng)、寬、高都為40 cm，井孔直徑為4 cm，如圖1所示。

圖2為花崗巖平板狀裂縫模型，花崗巖的骨架密度為2. 710 5 g/cm3，縱波速度為5 400 m/s，橫波速度為2 905 m/s。實(shí)驗(yàn)用到的設(shè)備有水槽和起吊架、尺子和游標(biāo)卡尺、手搖式千斤頂、一對(duì)60 kHz超聲探頭、CTS-8077PR超聲儀和DS1102E型數(shù)字示波器、TST3206動(dòng)態(tài)測(cè)試分析儀。

用2塊砂巖模型，制作凹凸鋸齒狀槽，如圖3所示，單邊刻槽，寬度為2 mm，深度為1 mm。在井孔周圍6個(gè)波長(zhǎng)范圍內(nèi)刻槽，井孔直徑為4 cm。測(cè)量源距為25 cm。

2.2 模型井實(shí)驗(yàn)數(shù)據(jù)分析

實(shí)測(cè)數(shù)據(jù)包括刻槽深度為1 mm，寬度為2 mm，源距為25 cm，頻率探頭為60 kHz。對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行分析，得到裂縫寬度與斯通利波衰減關(guān)系，包括刻槽前后的衰減曲線，如圖4所示。

圖1 模型井結(jié)構(gòu)Fig.1 Model well structure

圖2 平板狀裂縫模型Fig.2 Flat crack model

圖3 凹凸鋸齒狀槽模型Fig.3 Concave-convex serrated groove model

將刻槽前和刻槽后的斯通利波衰減少系數(shù)進(jìn)行比較，如圖5所示?？梢钥吹剑Ｐ途畬?shí)驗(yàn)的刻槽前和刻槽后的斯通利波斯通利波衰減系數(shù)在形態(tài)和趨勢(shì)上很接近。不同的是刻槽后的斯通利波衰減系數(shù)相比刻槽前的斯通利波衰減系數(shù)有小幅度的增加。這說明了裂縫粗糙度會(huì)造成斯通利波幅度衰減，但斯通利波衰減隨裂縫寬度的變化關(guān)系與平板光滑裂縫相比基本一致。

圖4 光滑裂縫和粗糙裂縫的衰減曲線Fig.4 Attenuation curves of smooth cracks and rough cracks with groove

圖5 刻槽前、后的斯通利波衰減擬合曲線Fig.5 Fitting curve of stoneley wave attenuation before and after grooving

3 模擬結(jié)果驗(yàn)證

3.1 數(shù)值模型與算法的實(shí)現(xiàn)

圖6為平板狀裂縫地層井眼模型示意圖。模型由井外地層、井眼和水平裂縫組成，裂縫傾角為0°，井眼和裂縫內(nèi)的填充物為流體，流體波速V=1 500 m/s，密度ρ=1 000 kg/m3。井外地層部分為均勻彈性固體地層，地層縱波速度VP=5 400 m/s，地層橫波速度VS=2 905 m/s，地層密度ρ1=2 710.5 kg/m3。地層的外邊界采用PML(perfectly matched layer)降低聲波反射。井眼半徑為0.101 6 m，聲源與接收點(diǎn)之間的距離為6英寸。激發(fā)聲源類型為變形瑞克子波，表達(dá)式為

式(2)中:f(t)為幅值，mV；t為時(shí)間，s；ω為最大能量頻率，rad/s。為了使接收到的斯通利波受縱橫波的影響最小，將聲源頻率設(shè)置為2.5 kHz。

為了研究斯通利波穿過裂縫的響應(yīng)情況，在研究模型的裂縫中采用非均勻的網(wǎng)格對(duì)裂縫進(jìn)行剖分，采用非均勻網(wǎng)格可以在基本不降低計(jì)算精度的情況下減小計(jì)算量[18]。變網(wǎng)格差分算子在空間分布上具有四階精度，圖7為變網(wǎng)格差分示意圖。圖7中節(jié)點(diǎn)i處的場(chǎng)量ψ對(duì)空間的導(dǎo)數(shù)?ψ/?x可以用4個(gè)節(jié)點(diǎn)處場(chǎng)量的組合表達(dá)式來表示[11]，即

圖6 平板狀裂縫地層井眼模型示意圖Fig.6 Schematic diagram of borehole model in flat fracture formation

圖7 變網(wǎng)格差分示意圖Fig.7 Variable grid difference schematic

(3)

式(3)為四階空間精度差分算子在變網(wǎng)格時(shí)的表現(xiàn)形式。在網(wǎng)格為均勻時(shí)，當(dāng)s1=r1=r2=Δr時(shí)，變網(wǎng)格差分算子改變?yōu)榫哂兴碾A空間精度均勻網(wǎng)格差分算子，即

(4)

3.2 模擬數(shù)據(jù)處理

3.2.1 粗糙裂縫模型

在進(jìn)行物理實(shí)驗(yàn)時(shí)，在裂縫表面刻凹凸鋸齒狀方槽，而天然裂縫表面并非是規(guī)則凹凸鋸齒狀方槽，所以在平板狀裂縫地層井眼模型的裂縫上下表面制作正弦形狀的槽，裂縫上下表面對(duì)稱，將三維裂縫模型簡(jiǎn)化為二維截面，如圖8所示。

迂曲度即流體實(shí)際滲流路徑與其視滲流路徑的比值，是描述粗糙裂縫結(jié)構(gòu)滲流的一個(gè)重要參數(shù)[30]，其計(jì)算公式為

(5)

式(5)中：τ為迂曲度；L1為流體實(shí)際滲流距離，m，在簡(jiǎn)化的二維模型中L1為彎曲曲線的長(zhǎng)度；L0為流體視滲流距離，m；在簡(jiǎn)化的二維模型中L0為裂縫的長(zhǎng)度，如圖8所示。對(duì)于槽寬L=5 mm，槽深A(yù)=2 mm的裂縫，經(jīng)過計(jì)算得到其迂曲度為1.092 4。

3.2.2 粗糙裂縫模型計(jì)算

對(duì)粗糙裂縫進(jìn)行數(shù)值模擬，圖9為槽寬L=5 mm，槽深A(yù)=2 mm的粗糙裂縫隨寬度變化的波形圖，圖10為槽寬L=5 mm，槽深A(yù)=2 mm的粗糙裂縫的不同裂縫寬度斯通利波幅度衰減與平板光滑裂縫的斯通利波衰減隨寬度變化的對(duì)比圖?？梢钥闯?，在裂縫寬度較小時(shí)，粗糙裂縫斯通利波的幅度衰減比平板光滑裂縫的斯通利波的幅度衰減明顯要小，這主要是因?yàn)榇植诘牧芽p表面在裂縫寬度較小時(shí)明顯地降低了裂縫的最小的滲流寬度。

A為槽的深度，L為槽寬，w為粗糙裂縫表面的最大裂縫寬度， w0為粗糙裂縫表面的最小裂縫寬度 圖8 粗糙裂縫二維截面Fig.8 Two-dimensional section of rough crack

圖9 正弦粗糙的裂縫隨寬度變化的波形圖Fig.9 Waveform diagram of sinusoidal rough crack changing with width

圖10 數(shù)值模擬平板光滑裂縫和粗糙裂縫的 斯通利波衰減對(duì)比Fig.10 Comparison of stoneley wave attenuation between smooth cracks and rough cracks in numerical simulation plate

3.2.3 迂曲度的影響

對(duì)于裂縫表面形如正弦的粗糙裂縫，迂曲度由槽寬和槽深決定。若固定槽深A(yù)，改變槽寬L相當(dāng)于粗糙裂縫的最大寬度和最小寬度不變，改變了最大寬度與最小寬度之間的填充程度。若固定槽寬L，改變槽深A(yù)，相當(dāng)于粗糙裂縫的最大寬度不變，改變了裂縫的最小寬度，最大寬度與最小寬度之間的填充程度的改變不大。這兩種情況都會(huì)造成迂曲度的改變，下面分別對(duì)這兩種情況進(jìn)行分析。

為了研究槽寬L對(duì)斯通利波的影響，取裂縫寬度為6、10、16、22 mm，取槽寬L分別為0.5、1、3、5、6、8、10 mm及槽寬為無窮大(相當(dāng)于光滑裂縫)，槽深為2 mm的4種不同裂縫寬度的粗糙裂縫進(jìn)行模擬，槽寬越小，迂曲度越大，計(jì)算結(jié)果用迂曲度表示，如圖11所示?？梢钥吹剑瑢?duì)于迂曲度的對(duì)數(shù)lgτ小于0.091 4的裂縫，隨著迂曲度的增加，斯通利波的衰減系數(shù)減小的速度很快，在迂曲度的對(duì)數(shù)lgτ大于0.091 4后，斯通利波的衰減系數(shù)減小的速度變慢。

當(dāng)槽寬為無窮小時(shí)，lgτ趨于無窮大，相當(dāng)于裂縫寬度由最大寬度w變?yōu)樽钚挾葁0。對(duì)于裂縫寬度為6、10、16、22 mm的粗糙裂縫，當(dāng)lgτ趨于無窮大時(shí)，通過圖11的擬合曲線可以計(jì)算得到斯通利波衰減系數(shù)和斯通利波等效寬度，并在圖10中標(biāo)出，分別對(duì)應(yīng)的點(diǎn)是y1、y2、y3、y4，在圖11中可以看到這些點(diǎn)都在對(duì)應(yīng)裂縫的最大、最小寬度之間，可以證明模擬的結(jié)果的正確性。

為了研究槽深A(yù)對(duì)斯通利波的影響，取裂縫寬度為22 mm，取槽深A(yù)分別為1、2、3、4、6、8、10 mm，槽寬為5 mm的粗糙裂縫進(jìn)行模擬，結(jié)果如圖12所示。結(jié)果顯示為隨著迂曲度的增大，斯通利波的衰減減小，造成這一結(jié)果的原因是因?yàn)橛厍鹊脑龃?，?dǎo)致裂縫之間的最小的滲流寬度減小，斯通利波的衰減減小。

圖11 寬為6、10、16、22 mm的裂縫隨迂曲度變化的 斯通利波衰減圖Fig.11 Stoneley wave attenuation plots of cracks with widths of 6, 10, 16, and 22 mm as a function of tortuosity

圖12 寬為22、26、28、30 mm的裂縫隨迂 曲度變化的斯通利波衰減圖Fig.12 Stoneley wave attenuation plots of fractures with widths of 22, 26, 28, and 30 mm as a function of tortuosity

選取槽寬為5 mm的54組不同的裂縫寬度和不同迂曲度的粗糙裂縫模型進(jìn)行計(jì)算，通過斯通利波衰減系數(shù)計(jì)算粗糙裂縫的斯通利波等效寬度(表1)。其中，等效裂縫寬度通過斯通利波衰減曲線換算，刻槽后測(cè)量的斯通利波衰減值在未刻槽測(cè)量的斯通利波衰減曲線上查到對(duì)應(yīng)的橫坐標(biāo)值，這個(gè)橫坐標(biāo)值就是對(duì)應(yīng)的等效裂縫寬度。

通過二元擬合建立粗糙裂縫的斯通利波等效寬度與裂縫寬度、迂曲度之間的關(guān)系，得到的粗糙裂縫等效寬度計(jì)算公式為

w′=6.492-7.261τ+0.749 9w+1.467τ2-0.227 2τw+0.016 63w2，R2=0.999 3

(6)

3.3 與物理模型試驗(yàn)結(jié)果對(duì)比

選取相同粗糙度的裂縫進(jìn)行模擬，對(duì)于模型井實(shí)驗(yàn)中槽寬為2 mm，槽深為1 mm的裂縫，源距選擇為25 cm，通過計(jì)算得到迂曲度為1.508。在數(shù)值模擬時(shí)，選取槽寬為5 mm，槽深為5.3 mm的裂縫進(jìn)行數(shù)值模擬，通過計(jì)算此時(shí)的迂曲度為1.508。將實(shí)驗(yàn)數(shù)據(jù)歸一化處理，圖13為迂曲度為1.508時(shí)的實(shí)驗(yàn)和數(shù)值模擬的斯通利波衰減擬合曲線。

表1 斯通利波等效寬度和裂縫寬度、迂曲度數(shù)據(jù)表Table 1 Data table of equivalent width, crack width and tortuosity of stoneley wave

圖13 實(shí)驗(yàn)和數(shù)值模擬的斯通利波衰減擬合曲線Fig.13 Stoneley wave attenuation fitting curve of experiment and numerical simulation

可以看到對(duì)于相同迂曲度的粗糙裂縫，實(shí)驗(yàn)和模擬的斯通利波衰減系數(shù)在形態(tài)和趨勢(shì)上基本一致。

4 結(jié)論

裂縫的粗糙性會(huì)對(duì)聲波的響應(yīng)產(chǎn)生影響，使得計(jì)算結(jié)果和光滑裂縫的計(jì)算結(jié)果產(chǎn)生偏差。通過物理實(shí)驗(yàn)和數(shù)值模擬分別研究了粗糙裂縫對(duì)斯通利波響應(yīng)的影響，得到粗糙裂縫寬度和斯通利波衰減的關(guān)系。分別對(duì)影響迂曲度的兩個(gè)參數(shù)(槽寬和槽深)進(jìn)行分析，研究了粗糙裂縫迂曲度、寬度對(duì)斯通利波衰減系數(shù)的影響，建立斯通利波衰減系數(shù)和迂曲度、裂縫寬度之間的關(guān)系。通過對(duì)實(shí)驗(yàn)和數(shù)值模擬結(jié)果的分析得到以下的結(jié)論。

(1)模型井粗糙裂縫實(shí)驗(yàn)的結(jié)果表明，在探頭頻率為60 kHz時(shí)，光滑裂縫和粗糙裂縫的斯通利波衰減系數(shù)和裂縫寬度的關(guān)系都為對(duì)數(shù)關(guān)系，粗糙裂縫的斯通利波幅度隨裂縫寬度變大而減小，斯通利波衰減系數(shù)隨裂縫寬度的增大而增大。

(2)對(duì)正弦結(jié)構(gòu)粗糙裂縫模型進(jìn)行數(shù)值模擬計(jì)算，將得到的裂縫寬度與斯通利波的衰減圖和光滑裂縫的結(jié)果進(jìn)行對(duì)比，結(jié)果顯示粗糙裂縫和光滑裂縫的斯通利波衰減系數(shù)和裂縫寬度的關(guān)系都為對(duì)數(shù)衰減，不同的是，裂縫的寬度越小，粗糙裂縫和光滑裂縫的斯通利波衰減系數(shù)之間的差值越大，造成這一現(xiàn)象的原因是因?yàn)椴凵畈蛔兊那闆r下隨著裂縫寬度的減小，槽深相對(duì)于裂縫寬度增大，減小了粗糙裂縫的滲流寬度。

(3)對(duì)正弦結(jié)構(gòu)粗糙裂縫模型的數(shù)值模擬結(jié)果顯示，當(dāng)槽深一定時(shí)，槽寬越小，迂曲度越大，隨著迂曲度的增加，斯通利波的衰減系數(shù)減小的速度由快變慢。當(dāng)槽寬一定時(shí)，槽深越大，迂曲度越大，隨著迂曲度的增大，斯通利波的衰減系數(shù)減小。