張國棟，莊春喜，黑　創(chuàng)

(1.中海石油(中國)有限公司上海分公司，上海 200335；2.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島 266580)

?

東海西湖凹陷探井儲層壓后縫高評價新方法

張國棟1，莊春喜2，黑創(chuàng)2

(1.中海石油(中國)有限公司上海分公司，上海 200335；2.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島 266580)

針對東海西湖凹陷多口探井壓裂后，陣列聲波交叉偶極子各向異性方法無法準(zhǔn)確評價裂縫縱向高度的技術(shù)難題，研究了評價探井儲層壓裂后縫高的新方法。通過研究壓后裂縫徑向變化特點，利用縱波走時層析成像技術(shù)重建壓裂前后井壁附近地層的二維速度剖面，如二維速度剖面壓裂后較壓裂前有明顯差異，說明儲層被壓開；同時，壓裂后會在主裂縫周圍形成一系列微裂縫，通過偶極聲波遠(yuǎn)探測成像技術(shù)對其進(jìn)行成像，可以評價裂縫的縱向高度和橫向深度。研究得知，與壓裂前相比，壓裂后形成的裂縫使井壁巖石聲速降低20%，井周裂縫形成的強(qiáng)反射和散射區(qū)域表明井眼周圍至少20 m的范圍內(nèi)形成了明顯的“壓裂體積改造”，通過對比可以準(zhǔn)確判斷裂縫的縱向高度，識別精度可達(dá)0.152 4 m。研究結(jié)果表明，縱波徑向?qū)游龀上窈团紭O聲波遠(yuǎn)探測成像2項技術(shù)可以解決無法準(zhǔn)確評價探井裂縫高度的難題，提升壓裂后的綜合評價效果。

探井；縱波層析成像；偶極橫波遠(yuǎn)探測成像；壓裂；裂縫高度；西湖凹陷

西湖凹陷是東海陸架盆地天然氣勘探的重點區(qū)域，近年在花港組發(fā)現(xiàn)了多套巨厚低孔低滲天然氣儲層[1-2]，但是DST測試產(chǎn)能低，達(dá)不到商業(yè)產(chǎn)能開采下限。為了提高產(chǎn)能，對多口探井儲層進(jìn)行了壓裂改造，但多數(shù)探井的壓裂效果不太理想，有些儲層物性和含氣性較好，但壓裂后天然氣產(chǎn)能低或直接出水。對于此類儲層，認(rèn)為有2方面的原因：1)壓裂參數(shù)設(shè)計和壓裂液選擇不合適；2)壓裂縫高控制不當(dāng)，導(dǎo)致裂縫縱向延伸過大，壓穿了上下水層導(dǎo)致出水。因此，研究裂縫縱向有效發(fā)育高度評價方法就顯得異常重要。

目前，評價井筒壓裂裂縫高度常用的測井方法主要有井溫測井、同位素測井、注硼中子測井、補(bǔ)償中子測井、偶極聲波測井等，其思路是根據(jù)井壁存在壓裂裂縫與井壁不存在壓裂裂縫時在測井資料上的不同響應(yīng)特征，來獲得壓裂裂縫高度信息，評價水力壓裂效果[3]。其中：井溫測井精度低，人為因素影響大，一般為輔助手段；同位素測井和注硼中子測井作業(yè)操作流程和評價原理基本相同，但兩者都有一定的放射性污染，目前應(yīng)用較少；補(bǔ)償中子測井可以采用非放射性示蹤陶粒，但受壓裂工藝(如脫砂的影響)、壓裂前后儀器類型及刻度的影響，測量結(jié)果具有一定的不確定性[4]；偶極聲波測井主要通過對比壓裂前后時差的各向異性來評價裂縫高度，該方法技術(shù)成熟，實施便捷，無污染，判別時直觀快速，是目前評價壓裂裂縫高度最常用的測井方法[5-8]。

2013年以來，西湖凹陷一直使用橫波各向異性來判斷壓裂裂縫高度，但應(yīng)用效果差別較大，當(dāng)壓裂后形成的是單一方向裂縫時，壓裂前后時差的各向異性差別明顯，容易識別裂縫高度；但當(dāng)壓裂后形成網(wǎng)狀縫時，儲層壓裂前后聲波的各向異性無明顯差異，就難以判斷裂縫高度。為此，筆者于2015年初利用縱波走時層析徑向成像技術(shù)重建了壓裂前后井壁附近地層的二維速度剖面，提出了一種利用壓裂前后縱波層析成像的差異來識別壓裂裂縫高度的方法，并結(jié)合偶極聲波遠(yuǎn)探測成像技術(shù)，準(zhǔn)確評價壓裂后各種裂縫的高度，而且還能在一定程度上評價裂縫的展布深度(井外20～30 m)，從而大大提升儲層壓裂評價效果。

1　縫高評價方法及原理

1.1縱波層析成像法

儲層被壓裂后，井壁周圍巖石發(fā)生破碎，會產(chǎn)生大量裂縫，使巖石擴(kuò)容，越靠近井壁破碎程度越大，而井壁巖石的波速降低程度則取決于巖石中裂縫的密度[9-10]。因此，從井壁到地層深處，地層波速度呈現(xiàn)出由低到高的徑向變化，據(jù)此可以得到井壁附近裂縫的密度及其影響范圍。若巖石被壓開，縱波徑向速度變化率與壓前相比則會增大；若巖石未被壓開，縱波徑向速度變化率與壓前相比則不會有明顯變化，因此該方法可用于評價壓裂裂縫的高度。

陣列聲波測井儀具有多個接收器(如圖1所示)，激發(fā)的聲波在波速徑向增加的地層中傳播；不同接收器接收聲波的徑向探測深度有所不同(見圖1中初至波的射線路徑)，源距越短探測深度越淺，源距越長探測的徑向深度越深，所以波的走時就含有地層變化的信息，從而可以通過確定地層縱波速度的變化，求出縱波沿井筒附近至原狀地層的速度變化率，利用縱波速度的徑向變化就可以表征儲層在徑向上的破損程度。

圖1　速度徑向變化地層中的聲波射線路徑示意Fig.1　Pathways of acoustic rays in formations with velocity radial changes

筆者采用B.E.Hornby提出的射線追蹤法[11]建立縱波徑向速度剖面，用以描述波速沿井眼軸向z和徑向r的變化?？v波徑向?qū)游龀上裉幚淼闹饕鞒倘缦拢?/p>

1) 基于貝葉斯準(zhǔn)則和自回歸模型，利用高精度時差提取方法獲得實測陣列波形的走時Tn，其中n表示某個接收器，n=1，2，3，…。走時的準(zhǔn)確提取是徑向?qū)游龀上竦年P(guān)鍵，具體可參閱文獻(xiàn)[12]。

2) 采用文獻(xiàn)[11]給出的射線追蹤法，選取二維速度剖面初始迭代模型(共有m×s個網(wǎng)格點)，計算聲波走時tn：

(1)

式中：Sn(vn)為依賴于地層速度模型vn(r,z)的聲波傳播的最短路徑。

3) 利用最小二乘法，求取滿足實測走時‖Tn-tn‖min的速度模型Vn(r,z)，其中Vn(r,z)?vn(r,z)。

4) 將步驟3)得到的模型設(shè)定為新的初始迭代模型，重復(fù)步驟3)，直至遍歷每個網(wǎng)格點，即可得到徑向聲速層析剖面。

圖2為XX-3井H4層壓裂前后縱波徑向?qū)游龀上駥Ρ葓D。其中，第4道為壓裂前的縱波層析成像圖，第5道為壓裂后的縱波層析成像圖，圖中的變化是徑向上聲速與原狀地層聲速相對差別的百分比，由深紅色到藍(lán)色，對應(yīng)的速度變化率為20%～0，探測范圍為井壁至井眼周圍1 m。通過壓裂前后縱波層析圖可以明顯看到，壓裂后縱波徑向速度發(fā)生了明顯變化，這表明此區(qū)間的巖石發(fā)生了破碎，裂縫在縱向上的延伸范圍為4 032.00～4 124.00 m，裂縫縱向延伸較長，在壓裂液體積有限的情況下(海上壓裂規(guī)模較小，一般壓裂液體積為250～400 m3)，裂縫徑向延伸會較淺，無法充分改善儲層的滲流條件，導(dǎo)致壓裂后產(chǎn)氣效果不佳，該井壓裂后僅見到微量氣。第3道是壓裂前后時差各向異性的對比，可以看出壓裂前后無明顯差別，難以判斷縫高。而壓裂前后縱波徑向?qū)游龀上駝t有顯著的差異性，能夠比較準(zhǔn)確地判斷裂縫高度，所以該方法明顯優(yōu)于時差各向異性對比法。

圖2　XX-3井H4層壓后聲波時差各向異性和縱波層析成像評價組合Fig.2　Acoustic time anisotropy and radial compression-wave velocity variation imaging for fractures in H4 layer of Well XX-3

1.2偶極聲波遠(yuǎn)探測成像法

X.M.Tang[13]于2004年提出了偶極橫波遠(yuǎn)探測聲波成像方法，該方法采用偶極子聲源在井中向井外發(fā)射橫波，并在井中接收由井外地質(zhì)體反射回來的橫波。圖3為單井偶極橫波遠(yuǎn)探測的原理示意圖。從圖3可以看出，偶極子聲源向井外輻射的彎曲波，可分解為平行于反射體的水平偏振的橫波(SH波)和垂直于反射體的垂向偏振的橫波(SV波)，SH波和SV波遇到反射體后會形成反射波而被接收器接收到，通過對接收到的反射波信息進(jìn)行提取，就可以得到反射體的位置及走向[13]。偶極子聲源的頻率約為2～5 kHz，橫向探測深度較深(可達(dá)20～30 m)。此外，由于偶極子聲源的指向具有方向性，采用多分量(實際中為四分量)的偶極發(fā)射和接收，可以確定反射體界面的方位[14]。但該方法不能確定該反射體是在井的右側(cè)還是井的左側(cè)，換言之，該方法只能確定反射體的走向，不能確定其傾向(具有180°不確定性)。

圖3　單井偶極橫波遠(yuǎn)探測的基本原理Fig.3　Fundamental principles of dipole acoustic imaging

儲層被壓裂后，一般會有單向縫和網(wǎng)狀縫2種裂縫形態(tài)，不管形成哪種類型的裂縫，在主裂縫的周圍都會產(chǎn)生大量走向不過井眼的微裂縫[15]，裂縫切面如圖4所示。這些裂縫就像在井眼周圍發(fā)育的一個個反射界面，而在偶極聲波遠(yuǎn)探測成像圖上，儲層被壓開處會看到明顯的反射界面，沒有壓開的地方則與壓裂前的偶極聲波遠(yuǎn)探測成像圖基本保持一致，通過對比壓裂前后偶極聲波遠(yuǎn)探測成像有助于評價裂縫井眼軸向上的發(fā)育高度。

圖4　壓裂后裂縫發(fā)育示意(切面圖)Fig.4　Fracture development after fracturing operations (cross section)

2　現(xiàn)場試驗

縱波徑向?qū)游龀上窈团紭O聲波遠(yuǎn)探測成像技術(shù)在東海西湖凹陷2口探井儲層壓后裂縫高度評價中進(jìn)行了現(xiàn)場試驗，取得了較好的效果。2口井均進(jìn)行了壓裂后陣列聲波測井，下面以XX-8井為例介紹其具體試驗情況。

對XX-8井P4砂層進(jìn)行了水力壓裂改造，壓裂后初期平均產(chǎn)氣量約4.0×104m3/d，壓裂改造較為成功，壓裂后進(jìn)行了偶極陣列聲波測井。該井壓裂前后聲波時差各向異性分析結(jié)果如圖5第3道所示，可以看出壓裂前后時差的各向異性并沒有明顯差別，無法準(zhǔn)確判斷壓裂裂縫的高度。為此，對該井進(jìn)行了縱波徑向?qū)游龀上窈团紭O聲波遠(yuǎn)探測技術(shù)處理，結(jié)果見圖5。其中，第4道和第5道為壓裂前后縱波徑向?qū)游龀上駥Ρ?，可見?4 195.00～4 260.00 m井段壓裂前后有明顯差異，顯示該井段井壁附近的巖石被壓開，所以壓裂裂縫的高度為65.00 m。

圖5　XX-8井P4層壓后縱波徑向?qū)游龀上窈团紭O聲波遠(yuǎn)探測成像組合Fig.5　Combination of radial compression-wave velocity variation image and dipole acoustic image after fracturing in Layer P4 of Well XX-8

為了得到該井壓裂后裂縫更準(zhǔn)確的發(fā)育高度，對陣列聲波數(shù)據(jù)進(jìn)行了偶極橫波遠(yuǎn)探測成像處理，其主要處理流程(見圖6)為：1)數(shù)據(jù)加載，該井?dāng)?shù)據(jù)為四分量正交偶極數(shù)據(jù)，數(shù)據(jù)質(zhì)量較好；2)數(shù)據(jù)預(yù)處理，主要是針對不同振型的波形，進(jìn)行濾波處理，去除一些低頻、高頻和噪聲，并對周期較長的彎曲波進(jìn)行壓陣處理；3)成像處理，主要是波場分離和偏移成像；4)成像后處理，通過數(shù)值圖像處理技術(shù)，壓制干擾噪聲，突出有效反射信號[16]。由于該井陣列聲波數(shù)據(jù)質(zhì)量較好，遠(yuǎn)探測成像質(zhì)量較好，如圖5第7道和第8道所示，4 200.00～4 240.00 m井段壓裂后有明顯的強(qiáng)反射界面，顯示該井段巖石被壓開，裂縫高度約40.00 m。

圖6　偶極橫波遠(yuǎn)探測成像處理流程Fig.6　Processing of dipole acoustic images

從分析效果看，縱波層析成像和偶極橫波遠(yuǎn)探測成像兩者的評價結(jié)果存在一定差別，前者顯示裂縫發(fā)育井段為4 195.00～4 260.00 m，而后者則顯示發(fā)育井段為4 200.00～4 240.00 m。分析認(rèn)為，其原因是縱波層析成像主要利用直達(dá)縱波，反映的是井眼周圍約1 m范圍內(nèi)的地層速度變化率；而偶極聲波遠(yuǎn)探測成像利用的是反射橫波，反映的是遠(yuǎn)井筒處的反射面，兩者的分析對象有所不同，所以需要將兩項技術(shù)結(jié)合起來進(jìn)行綜合評價。

綜合分析認(rèn)為，4 195.00～4 260.00井段都發(fā)生了破裂，但是橫向延伸主裂縫主要在4 200.00～4 240.00 m井段發(fā)育，裂縫橫向延伸遠(yuǎn)，至少在井眼周圍20 m范圍內(nèi)看到了明顯的微裂縫發(fā)育。但是4 200.00～4 240.00 m井段儲層未被充分壓開，改造程度不夠，僅在井壁附近發(fā)生了破裂，這在一定程度上影響了測試產(chǎn)能。

通過兩口井的實際應(yīng)用可以看出，壓裂后形成的裂縫使原狀地層巖石聲速降低20%左右，井周裂縫形成的強(qiáng)反射和散射區(qū)域表明井眼周圍至少20 m范圍內(nèi)形成了明顯的“壓裂體積改造”，通過對比可以準(zhǔn)確判斷裂縫的縱向高度，識別精度可達(dá)0.152 4 m。需要說明的是，由于儀器的采集記錄時長有限，目前的常規(guī)陣列聲波儀器的偶極聲波遠(yuǎn)探測成像只能探測到井外20～30 m，此深度之外裂縫的橫向延伸情況目前無法評價。

3　結(jié)論及認(rèn)識

1) 陣列聲波測井是評價東海西湖凹陷探井壓裂裂縫高度的主要手段，實踐表明，壓裂后形成單一方向的裂縫時，通過對比壓裂前后時差的各向異性評價裂縫高度效果較好。如果壓裂后形成的是網(wǎng)狀縫，則不同方向時差各向異性無明顯差別，借助于時差各向異性難以評價裂縫高度。

2) 縱波層析成像和偶極聲波遠(yuǎn)探測成像2種陣列聲波測井處理方法，能夠彌補(bǔ)聲波時差各向異性的不足，不僅可以有效判別壓裂裂縫高度，還可以在一定范圍內(nèi)評價壓裂裂縫的徑向長度，進(jìn)一步提升了壓后評價效果。目前該技術(shù)已經(jīng)在西湖凹陷2口探井進(jìn)行了試驗，取得了較好的效果。

3) 采用陣列聲波測井方法評價壓裂裂縫高度時，建議綜合使用聲波時差各向異性、縱波層析成像和偶極聲波遠(yuǎn)探測成像等處理方法，以提高評價精度。

[1]張國華.西湖凹陷高壓形成機(jī)制及其對油氣成藏的影響[J].中國海上油氣，2013，25(2)：1-8.

ZHANG Guohua.Origin mechanism of high formation pressure and its influence on hydrocarbon accumulation in Xihu Sag[J].China Offshore Oil and Gas,2013,25(2):1-8.

[2]施榮富.西湖凹陷低孔低滲儲層壓裂改造技術(shù)體系探索與實踐[J].中國海上油氣，2013，25(2)：79-82.

SHI Rongfu.Exploration and application of fracturing reconstruction technique for low porosity and low permeability reservoirs in Xihu Sag[J].China Offshore Oil and Gas,2013,25(2):79-82.

[3]易新民，唐雪萍，梁濤，等.利用測井資料預(yù)測判斷水力壓裂裂縫高度[J].西南石油大學(xué)學(xué)報(自然科學(xué)版)，2009，31(5)：21-24.

YI Xinmin,TANG Xueping,LIANG Tao,et al.Prediction and assessment of fracture height of hydraulic fracturing with logging data[J].Journal of Southwest Petroleum University(Science & Technology Edition),2009,31(5):21-24.

[4]溫偉.補(bǔ)償中子測井在水力壓裂縫高檢測中的應(yīng)用[J].遼寧化工，2013，42(7)：828-830.

WEN Wei.Application of compensated neutron log in fracture height detection after hydraulic fracturing[J].Liaoning Chemical Industry,2013,42(7):828-830.

[5]竇偉坦，侯雨庭.利用偶極聲波測井進(jìn)行儲層壓裂效果評價[J].中國石油勘探，2007，12(3)：58-63.

DOU Weitan,HOU Yuting.Evaluation on fracturing effect of reservoirs based on dipole sonic logging[J].China Petroleum Exploration,2007,12(3):58-63.

[6]TANG X M,PATTERSON D,HINDS M.Evaluating hydraulic fracturing in cased holes with cross-dipole acoustic technology[J].SPE Reservoir Evaluation & Engineering,2001,4(4):281-288.

[7]陳颙，黃庭芳，劉恩儒.巖石物理學(xué)[M].合肥：中國科學(xué)技術(shù)大學(xué)出版社，2009：118-119.

CHEN Yong,HUANG Tingfang,LIU Enru.Rock physics[M].Hefei:China University of Science and Technology Press,2009:118-119.

[8]WINKLER W K.Borehole damage indicator from stress-induced velocity variations[J].Geophysics,2005,70(1):F11-F16.

[9]THOMSEN L.Biot-consistent elastic moduli of porous rocks:low-frequency limit[J].Geophysics,1985,50(12):2797-2807.

[10]TANG X M,CHEN X L,XU X K.A cracked porous medium elastic wave theory and its application to interpreting acoustic data from tight formations[J].Geophysics,2012,77(6):D245-D252.

[11]HORNBY B E.Tomographic reconstruction of near-borehole slowness using refracted borehole sonic arrivals[J].Geophysics,1993,58(12):1726-1738.

[12]郎曉政，蘇遠(yuǎn)大，莊春喜，等.高精度時差提取技術(shù)在薄互層中的應(yīng)用[J].應(yīng)用聲學(xué)，2014，33(3)：252-257.

LANG Xiaozheng,SU Yuanda,ZHUANG Chunxi,et al.Application of high-precision slowness technic in the process of thin-bed sand-shale sequence model[J].Journal of Applied Acoustics,2014,33(3):252-257.

[13]TANG X M.Imaging near borehole structure using directional acoustic wave measurement[J].Geophysics,2004,69(6):1378-1386.

[14]TANG X M,PATTERSON D.Single-well S-wave imaging using multi-compoent dipole acoustic log data[J].Geophysics,2009,74(6):211-223.

[15]唐曉明，魏周拓.利用井中偶極聲源遠(yuǎn)場輻射特性的遠(yuǎn)探測測井[J].地球物理學(xué)報，2012，55(8)：2798-2807.

TANG Xiaoming,WEI Zhouta.Single-well acoustic reflection imaging using far-field radiation characteristics of a borehole dipole source[J].Chinese Journal of Geophysics,2012,55(8):2798-2807.

[16]莊春喜，燕菲，孫志峰，等.偶極橫波遠(yuǎn)探測測井?dāng)?shù)據(jù)處理及應(yīng)用[J].測井技術(shù)，2014，38(3)：330-336.

ZHUANG Chunxi,YAN Fei,SUN Zhifeng,et al.Data processing and applications of dipole shear-wave imaging logging[J].Well Logging Technology,2014,38(3):330-336.

[編輯令文學(xué)]

New Techniques for Fracture Height Determination in Exploration Wells Drilled in the Xihu Sag,East China Sea

ZHANG Guodong1,ZHUANG Chunxi2,HEI Chuang2

(1.ShanghaiBranchofCNOOCLtd.,Shanghai,200335,China; 2.SchoolofGeosciences&Technology,ChinaUniversityofPetroleum(Huadong),Qingdao,Shandong,266580,China)

Anisotropy analyses of conventional array acoustic cross dipole data failed to determine the vertical heights of fractures formed in reservoir formations in exploration wells drilled in the Xihu Sag,East China Sea.To eliminate such technical challenges,innovative techniques have been developed to determine the fracture heights in reservoir formations after fracturing in such exploration wells.To determine the features in radial directions after fracturing operations,2D velocity profiles for formations in vicinity of sidewalls before and after fracturing were created by using the longitudinal wave travel time tomography technique.It was found that such fracturing operations might make a huge difference such velocity profiles if the formation were fractured successfully.Meanwhile,the hydraulic fracturing operations might also create a series of micro-fractures around the major fractures.Capable of detecting those micro-fractures,the dipole acoustic far-field images could be used to determine heights and widths of the resulting fractures.Research results showed that the acoustic speed reduced by 20% after the fracture were created near the wellbore,whereas the strong acoustic reflection and scattering confirm the presence of volumetric fractures were created within 20 m around the well.The dipole acoustic far-field images could identify the fracture heights with high resolution up to 0.152 4 m.Combination of the longitudinal wave travel time tomography technique and the dipole acoustic far-field images could effectively enhance accuracy in the determination of fracture heights in exploration wells.Eventually,the comprehensive performances of fracturing operations can be assessed more accurately.

exploratory well; radial compression-wave velocity variation imaging; dipole acoustic imaging; fracturing; fracture height; Xihu Sag

2016-01-18；改回日期：2016-08-01。

張國棟(1984—)，男，山東成武人，2006年畢業(yè)于中國石油大學(xué)(華東)資源勘查專業(yè)，工程師，從事測井研究與管理工作。E-mail:zhanggd4@cnooc.com.cn。

中國海洋石油總公司重大科研項目“東海低孔低滲氣藏勘探開發(fā)關(guān)鍵技術(shù)研究與實踐”(編號:CNOOC-KJ 125 ZDXM 07 LTD 04 SH 2011)、山東省自然科學(xué)基金“非常規(guī)儲層壓裂裂縫延伸長度的過套管聲波測井評價應(yīng)用基礎(chǔ)研究”(編號:ZR2014DL009)資助。

?測井錄井?doi:10.11911/syztjs.201605021

P631.8+14；TE357.1+1

A

1001-0890(2016)05-0122-05