陣列聲波測井海上油氣田水力壓裂效果評價技術(shù)

張明，劉峰，張博，張璋，祁曉，李東，張聰慧

1.中海石油(中國)有限公司天津分公司，天津 300459

2.中海油田服務(wù)股份有限公司油田技術(shù)事業(yè)部，天津 300459

經(jīng)過60余年的發(fā)展與完善，壓裂技術(shù)經(jīng)歷了探索、研究與開發(fā)的階段，逐步發(fā)展成熟[1-2]。該技術(shù)由早期單純以增產(chǎn)為目的的單項技術(shù)，逐漸發(fā)展為與勘探、開發(fā)、采油等專業(yè)相結(jié)合的綜合技術(shù)體系，其配套發(fā)展的壓裂效果評價技術(shù)更是有力推動了壓裂工藝、技術(shù)、設(shè)備的不斷進(jìn)步。從油氣田勘探開發(fā)的角度來看，深入了解和認(rèn)識壓裂裂縫的幾何形態(tài)、延伸情況、發(fā)育程度、分布及方位等信息，對于后續(xù)油氣勘探開發(fā)至關(guān)重要，充分說明了壓裂效果評價的重要性，國內(nèi)外各大油氣田都開展了大量的研究工作。然而，由于壓裂系統(tǒng)本身所具有的多目標(biāo)性、多層次性、動態(tài)性、信息不完備等特點，使壓裂效果評價成為一個復(fù)雜而又困難的問題[3]。

壓裂效果監(jiān)測常采用井溫法、放射性測量法等常規(guī)壓裂監(jiān)測技術(shù)，但在實際應(yīng)用過程中上述方法均表現(xiàn)出壓裂層段定位不準(zhǔn)、環(huán)境污染等局限性[4]。目前非常規(guī)油氣勘探主要應(yīng)用微地震監(jiān)測技術(shù)進(jìn)行壓裂效果監(jiān)測，該技術(shù)主要通過在鄰井中(或地面)布設(shè)檢波器，通過監(jiān)測壓裂井在壓裂過程中誘發(fā)的微地震波來描述壓裂過程中裂縫生長的幾何形狀和空間展布，能夠?qū)崟r提供壓裂施工過程中產(chǎn)生的裂隙高度、長度和方位角等信息[5]。但是由于海上作業(yè)成本問題以及作業(yè)難度的限制，微地震技術(shù)目前無法應(yīng)用于海上壓裂效果評價[6]。利用正交偶極各向異性反演評價壓裂縫高度也是較為傳統(tǒng)的壓裂效果評價技術(shù)[7]，在世界范圍內(nèi)獲得了較為廣泛的應(yīng)用，該方法直觀準(zhǔn)確，但也存在很多局限性，會受壓裂縫網(wǎng)形態(tài)[8]、地層相對于井眼的角度等影響，且手段單一，無法精細(xì)刻畫壓裂對井壁附近地層的改造作用，不能滿足海上儲層壓裂效果評價的需求。

為了滿足渤海油田儲層壓裂效果評價需求，筆者提出了一種基于陣列聲波測井資料的水力壓裂效果評價技術(shù)，該技術(shù)主要采用在壓裂前、后分別進(jìn)行陣列聲波測井，再利用壓裂前、后的聲波時差，橫波幅度，徑向速度剖面的差異評價壓裂縫高度，結(jié)合偶極聲波遠(yuǎn)探測成像技術(shù)，在一定程度上評價壓裂縫高度，實現(xiàn)了對儲層壓裂縫高度的定性評價。該評價技術(shù)在渤海油田開發(fā)井中取得了良好的應(yīng)用效果。

1 水力壓裂效果評價技術(shù)

1.1 聲波時差、橫波幅度響應(yīng)特征分析

為研究地層在壓裂作業(yè)前后，井孔聲場中地層時差、橫波幅度的變化情況，本文建立了壓裂地層模型，采用有限差分?jǐn)?shù)值模擬壓裂后地層時差、橫波幅度的測井響應(yīng)特征。實際地層中介質(zhì)通常是非均質(zhì)的，地層在壓裂作業(yè)之后，由于水力壓裂區(qū)域的地層通常會呈現(xiàn)出各向異性和隨機特性，因此本次研究利用隨機介質(zhì)模型來描述巖石中的非均質(zhì)起伏變化[9-11]。二維井孔水力壓裂模型如圖1所示，設(shè)置垂直于井眼方向X的長度為20 m，沿井眼方向Z的長度為30 m，建立坐標(biāo)軸，井眼中心軸所在位置為X=0，經(jīng)過模擬水力壓裂區(qū)域的橢圓隨機介質(zhì)中心點坐標(biāo)為(0，15)，X方向軸長為14 m，Z方向軸長為10 m，模型具體參數(shù)見表1。

表1 模型計算參數(shù)

圖1 二維井孔水力壓裂模型

將位移運動方程離散化得到二維彈性波動方程的有限差分形式[12]，如式(1)所示：

(1)

模擬例子中假設(shè)點源和接收陣列均處于井眼中心軸的流體中(X=0)。首先分別模擬壓裂前的均質(zhì)地層與壓裂后的非均質(zhì)地層中，單極子、偶極子兩種聲源下的單極聲波測井陣列波形、偶極聲波測井陣列波形，計算所用聲源類型為雷克子波，聲源中心頻率分別為10、3 kHz，聲波儀器源距為3 m，間距為0.152 4 m。數(shù)值模擬結(jié)果如圖2所示。圖2(a)為壓裂前、后的單極聲波測井理論模擬陣列波形。與壓裂前的地層縱波(紅色曲線)相比，壓裂后地層中的地層縱波(黑色曲線)到時明顯滯后，說明地層受水力壓裂作用，地層破碎帶的地層縱波速度明顯降低，且壓裂后的縱波初至斜率明顯小于壓裂前的，即地層縱波時差明顯增大；圖2(b)為壓裂前、后的偶極聲波測井理論模擬陣列波形，壓裂前、后的偶極橫波初至與圖2(a)中單極縱波初至呈現(xiàn)相同特征，說明地層破碎帶處地層橫波時差同樣明顯增大。因此，地層受水力壓裂作用之后，壓裂前、后的地層縱波、橫波時差的差異可以反映地層的壓裂情況。通過對比壓裂前、后的地層縱、橫波時差變化，時差差異明顯井段即為壓裂縫存在的位置，該井段長度即壓裂縫沿井壁延伸的高度。利用聲波時差評價地層的壓裂情況存在一定的局限性，根據(jù)本多夫定律[13]，聲波在射線出射處沿井軸方向的視速度等于聲波在其最大穿透深度處的地層速度，該視速度即為波在陣列中傳播的平均速度[14]，該視速度的倒數(shù)即為聲波時差。當(dāng)?shù)貙訅毫岩?guī)模不大時，壓裂前、后的地層縱、橫波時差均為原狀地層縱、橫波時差，沒有明顯差異，此時通過對比壓裂前、后的地層縱、橫波時差變化不能評價地層的壓裂情況。

圖2 壓裂前、后的單極、偶極聲波測井陣列波形理論模擬

為了進(jìn)一步研究偶極橫波幅度在壓裂前后的響應(yīng)特征，本文對模擬所得壓裂前、后的單道偶極聲波測井波形進(jìn)行了分析，如圖3所示。其中紅色波形為模擬所得壓裂前均質(zhì)地層的單道偶極聲波波形，黑色波形為模擬所得壓裂后非均質(zhì)地層的單道偶極聲波波形。與壓裂前相比，壓裂后的偶極橫波幅度明顯降低，說明壓裂使井壁周圍巖石破碎，井壁附近地層非均質(zhì)性增加，聲波從聲源到接收器傳播時，受地層非均質(zhì)性影響，偶極橫波幅度明顯降低，因此壓裂前、后的橫波波幅變化可反映地層的壓裂情況。

圖3 壓裂前、后的單道偶極子聲波測井理論模擬波形

基于上述數(shù)值模擬結(jié)果，利用壓裂前、后的聲波時差、橫波波幅差異進(jìn)行儲層壓裂效果評價。本次研究采用均方根法計算壓裂前、后陣列聲波的橫波幅度，公式如下：

(2)

由于壓裂前、后的井況存在差異，所以僅利用橫波幅度絕對值難以對壓裂程度進(jìn)行準(zhǔn)確描述，對此，引入橫波幅度相對值。在壓裂未波及且井況良好的井段，橫波幅度在壓裂前、后可視為無變化，可將該井段壓裂前的橫波幅度曲線設(shè)置為基準(zhǔn)線，通過線性移動該井段壓裂后的橫波幅度曲線，使之與基準(zhǔn)線重合，即可獲得橫波幅度校正量，然后對全井段壓裂后的橫波幅度曲線采用該校正量進(jìn)行校正，校正完畢后與壓裂前的橫波幅度曲線進(jìn)行對比。

當(dāng)儲層未被壓裂時，壓裂前、后的橫波幅度無明顯差異；當(dāng)儲層被壓裂后產(chǎn)生明顯壓裂縫時，壓裂后的橫波幅度明顯下降，通過對比壓裂前、后的橫波幅度變化，橫波幅度差異明顯井段即為壓裂縫存在的位置，該井段長度即壓裂縫沿井壁延伸的高度。

利用橫波幅度評價地層的壓裂情況存在一定的局限性，當(dāng)壓裂前、后測井作業(yè)選用的聲波測井儀器不同時，壓裂前、后的橫波幅度受儀器類型影響存在明顯，此時通過對比壓裂前、后的橫波幅度變化已不能真實地反映地層的壓裂情況。

1.2 聲波各向異性反演技術(shù)

在裸眼井情況下，采用聲波各向異性反演技術(shù)可進(jìn)行壓裂效果評價，主要利用偶極橫波在各向異性地層中產(chǎn)生的橫波分裂現(xiàn)象，通過對比壓裂前、后的地層各向異性結(jié)果能夠?qū)毫芽p沿井壁縱向發(fā)育范圍(通常稱為“壓裂縫高度”)進(jìn)行準(zhǔn)確評價，但利用該方法進(jìn)行壓裂效果評價仍存在問題。當(dāng)壓裂作業(yè)所造成的縫網(wǎng)系統(tǒng)不是主縫，而是沒有優(yōu)勢方向的網(wǎng)狀裂縫時，偶極橫波不發(fā)生分裂，導(dǎo)致壓裂前后各向異性基本無變化，此時采用聲波各向異性反演技術(shù)難以對壓裂效果進(jìn)行有效評價。

現(xiàn)階段大多數(shù)壓裂作業(yè)均在套管井中進(jìn)行，為了進(jìn)行壓裂效果評價，需要在套管井中進(jìn)行陣列聲波測井，這時需要考慮套管對陣列聲波測井資料的影響。在套管井情況下，現(xiàn)已證明透過套管測量地層橫波各向異性的可行性[15]，如果套管與水泥層、地層之間固井質(zhì)量良好，則采用偶極聲波測井可以測量地層的橫波各向異性，但如果套管與水泥、地層之間固井質(zhì)量不好，套管中激發(fā)產(chǎn)生的管波與偶極子波互相干擾，尤其當(dāng)管波速度與偶極子波速度相當(dāng)時，會嚴(yán)重影響偶極聲波測量結(jié)果。此時，利用聲波各向異性反演技術(shù)獲得的地層各向異性結(jié)果中包含了套管的影響，不能真實地反映地層壓裂效果，為該項技術(shù)的應(yīng)用帶來了一定困擾。

1.3 縱波走時層析成像技術(shù)

縱波走時層析成像技術(shù)能夠很好地解決上述問題，該技術(shù)主要采用唐曉明等[14]、HORNBY[16]提出的射線追蹤法建立縱波徑向速度剖面，對于均質(zhì)地層，縱波以沿著井壁的滑行縱波傳播，其到達(dá)每一個接收器的波至?xí)r間(即波的走時)應(yīng)為等差特征，包含地層的真實速度信息。對于徑向速度發(fā)生變化的地層，射線由淺到深折射進(jìn)入地層后，以滑行縱波傳播，再由深到淺折射回接收器，具體傳播路徑如圖4所示[17]。由縱波的傳播路徑可知，縱波走時中包含了儀器探測范圍內(nèi)地層速度變化情況，因此采用縱波走時層析成像技術(shù)可以重建井壁附近地層的速度變化剖面。

注：Rdf表示儀器到深部地層的距離，Rmf表示儀器到中部地層的距離，Rsf表示儀器到淺部地層的距離，Rbh表示儀器到井壁的距離。

壓裂能夠造成井壁周圍巖石破碎并產(chǎn)生大量壓裂縫，進(jìn)而形成巖石孔隙體積的增加(一般稱為“巖石的擴容”)，而巖石擴容會造成巖石彈性波速下降已被大量實驗數(shù)據(jù)證實[18]。隨著徑向距離的增加，巖石破碎擴容程度也逐漸減弱，而縱波走時層析成像技術(shù)恰好能夠?qū)υ摤F(xiàn)象進(jìn)行良好表征。當(dāng)壓裂效果差或井壁不存在壓裂縫時，由于徑向速度未發(fā)生明顯變化，因此聲波徑向速度剖面成像結(jié)果無異常顯示。當(dāng)壓裂效果較好、近井壁破碎程度較高、破碎程度隨徑向距離增大而減弱時，聲波徑向速度剖面成像結(jié)果會有明顯表現(xiàn)，通過分析壓裂前、后聲波徑向速度剖面的差異，主要包括徑向深度是否明顯加深、徑向速度變化量是否變大，從而確定壓裂縫位置，徑向速度發(fā)生變化的井段長度即為壓裂縫沿井壁延伸的高度，同時該技術(shù)不受儀器類型、地層壓裂規(guī)模、壓裂縫形態(tài)等因素的影響，能夠直觀顯示近井壁附近地層的壓裂情況。

1.4 偶極聲波遠(yuǎn)探測成像技術(shù)

TANG[19]于2004年提出了偶極聲波遠(yuǎn)探測技術(shù)，首次把偶極子用于單井反射聲波測井中。該技術(shù)主要采用偶極子聲源從井內(nèi)向井外發(fā)射橫波，將井中接收到的由井外地質(zhì)體反射回來的橫波進(jìn)行處理，對井旁地質(zhì)構(gòu)造進(jìn)行成像。利用該技術(shù)徑向探測深度遠(yuǎn)這一優(yōu)勢，在壓裂效果評價中引入該技術(shù)，通過對比壓裂前、后各方位上，聲波遠(yuǎn)探測成像差異，對遠(yuǎn)井壁壓裂改造程度進(jìn)行定性評價。

地層被壓裂之后，一般會在地層中產(chǎn)生壓裂縫，在主壓裂縫周圍都會產(chǎn)生大量的井旁微裂縫(見圖5)。為了進(jìn)一步分析這些壓裂縫在聲波遠(yuǎn)探測成像中的成像特征，通過正演模擬了網(wǎng)狀縫形態(tài)下的反射波聲場，然后采用偶極橫波遠(yuǎn)探測技術(shù)對所測聲波測井資料進(jìn)行處理，獲得相應(yīng)的遠(yuǎn)探測成像圖，結(jié)果如圖6所示。圖6(a)為正演模擬建立的網(wǎng)狀縫模型，利用Born近似計算得到該模型的偶極散射波場響應(yīng)特征[20]，圖6(b)為偶極聲波遠(yuǎn)探測成像圖。從數(shù)值模擬結(jié)果可知，網(wǎng)狀縫的形態(tài)主要由近井壁附近形成與井軸垂直的主壓裂縫、遠(yuǎn)井壁處主壓裂縫周圍產(chǎn)生大量的井旁微裂縫組成。由于主壓裂縫與井軸垂直，井內(nèi)接收不到其反射回來的橫波信息，在遠(yuǎn)探測成像圖中近井壁附近沒有主壓裂縫的反射信息，在遠(yuǎn)井壁處可以看到與地質(zhì)模型尺寸大小相當(dāng)?shù)摹伴L弧狀”反射信息，這是沿主裂縫上下擴展形成的井旁微裂縫界面的反射信息。

圖5 壓裂后裂縫發(fā)育切面圖

圖6 網(wǎng)狀縫數(shù)值模擬及成像結(jié)果

基于上述研究，地層在受到水力壓裂作用后，在近井壁附近，壓裂前、后的地層縱、橫波時差的差異、地層橫波幅度差異、聲波各向異性、徑向速度剖面差異均可以反映地層在近井壁處的壓裂情況，尤其是壓裂縫沿井壁延伸的高度。但前3種參數(shù)在應(yīng)用中具有一定的局限性，徑向速度變化具有更好的適用性，因此，當(dāng)不同方法給出的壓裂縫高度不同時，以徑向速度剖面差異的評價結(jié)果為準(zhǔn)；在遠(yuǎn)井壁處，采用偶極聲波遠(yuǎn)探測成像技術(shù)能夠評價壓裂縫在遠(yuǎn)處的展布情況。鑒于上述分析，本文提出了一種水力壓裂效果綜合評價技術(shù)，主要利用壓裂前、后的聲波時差、橫波幅度、徑向速度剖面的差異對近井壁處壓裂縫高度進(jìn)行評價，同時結(jié)合偶極聲波遠(yuǎn)探測成像技術(shù)，在一定程度上評價遠(yuǎn)井壁處的壓裂縫高度，通過“遠(yuǎn)近結(jié)合”的方式，實現(xiàn)了對儲層壓裂縫高度及壓裂程度的定性評價。

2 應(yīng)用實例

水力壓裂效果綜合評價技術(shù)在渤海油田開發(fā)井儲層壓裂效果評價中進(jìn)行了現(xiàn)場應(yīng)用，取得了較好的應(yīng)用效果。渤海油田某開發(fā)區(qū)塊目前共計完鉆13口探井，在新近系明化鎮(zhèn)組明下段、館陶組以及古近系的東營組和沙河街組均鉆遇了良好的油氣顯示，且該油田在開發(fā)階段獲得了良好的油氣產(chǎn)出。為了完善注采井網(wǎng)、提高油田開發(fā)效果，設(shè)計了一口開發(fā)井，該井主力目的層為沙河街組，其儲層孔隙度范圍為12.9%～19.7%，平均孔隙度為16.5%，滲透率范圍為2.5～39.9 mD，平均滲透率為19.7 mD，屬于中低孔滲儲層。僅依靠其自身的產(chǎn)能難以形成工業(yè)油氣，為了釋放產(chǎn)能獲得更好的油氣產(chǎn)出，對該井沙河街組沙Ⅰ段中3 941.0～3 943.0 m井段，采用負(fù)壓返涌方式進(jìn)行了射孔壓裂作業(yè)。

該井于壓裂作業(yè)前、后的套管井中分別進(jìn)行了陣列聲波測井，采用水力壓裂效果綜合評價技術(shù)對壓裂前、后的陣列聲波資料進(jìn)行了處理。圖7是該井壓裂前、后的水力壓裂效果綜合評價技術(shù)處理成果圖。該井壓裂前全井各向異性較穩(wěn)定，各向異性度為4%～8%，壓裂后各向異性度變化不大，無法判斷壓裂縫高度，需要進(jìn)一步對壓裂前、后的聲波時差、橫波幅度、徑向速度剖面、聲波遠(yuǎn)探測結(jié)果進(jìn)行綜合分析。

在3 932.2～3 946.0 m井段，壓裂后的地層縱波、橫波時差明顯增大，地層橫波幅度明顯降低，其中縱波時差平均增量為6 μs/ft，橫波時差平均增量為16 μs/ft，橫波波幅平均降低47%；壓裂前、后的地層徑向速度剖面差異明顯，對壓裂前、后速度剖面進(jìn)行做差處理，如圖7中的第7道所示，該井段徑向速度剖面變化最為顯著，表現(xiàn)為縱波速度明顯小于原狀地層縱波速度。綜合上述測井響應(yīng)，說明3 932.2～3 946.0 m井段近井壁附近的地層被壓開，形成了一定規(guī)模的壓裂裂縫，為壓裂主要作用段，近井壁壓裂縫高度為13.8 m。為了進(jìn)一步得到該井壓裂后遠(yuǎn)井壁處的壓裂縫發(fā)育情況，對壓裂前、后的陣列聲波數(shù)據(jù)進(jìn)行了偶極橫波遠(yuǎn)探測成像處理，結(jié)果如圖7中的第8道和第9道所示。在3 929.0～3 947.0 m井段壓裂后井旁存在明顯的強反射界面，顯示該井段地層被壓開，近井壁壓裂縫高度約為18 m。

從分析結(jié)果來看，在3 932.2～3 946.0 m井段地層在近井壁和遠(yuǎn)井范圍內(nèi)都發(fā)生了破裂，且地層橫向延伸主裂縫發(fā)育在該井段內(nèi)，隨著主裂縫在地層中延伸，主壓裂縫周圍產(chǎn)生了大量的井旁微裂縫；遠(yuǎn)探測成像圖中顯示，至少在10 m范圍內(nèi)看到了明顯的“長弧狀”反射信息，結(jié)合正演模擬結(jié)果，該反射信息為井旁微裂縫反射界面，說明該井至少在10 m范圍內(nèi)形成了明顯的“壓裂體積改造”。

在3 903.0～3 911.0 m井段，壓裂前、后的徑向速度剖面差異明顯，然而該井段單極全波中首波為能量較強的套管波，說明本井段固井質(zhì)量較差，認(rèn)為該差異可能受到套管波的影響不能完全反映地層壓裂效果。進(jìn)一步結(jié)合聲波時差、橫波幅度曲線、聲波遠(yuǎn)探測結(jié)果，地層縱波、橫波時差有所增大，地層橫波幅度明顯降低，遠(yuǎn)探測成像沒有明顯差異，認(rèn)為該井受壓裂作業(yè)影響，水泥環(huán)中出現(xiàn)流體流動通道，壓裂液沿著流動通道進(jìn)入3 903.0～3 911.0 m儲層中，僅對該儲層近井壁附近造成了小規(guī)模改造，近井壁壓裂縫高度為8 m，遠(yuǎn)井壁處未受到壓裂作業(yè)影響。

總體來看，該井在近井壁附近產(chǎn)生壓裂縫高度為21.8 m，遠(yuǎn)井壁至少在10 m范圍內(nèi)形成了明顯的“壓裂體積改造”。該井壓裂前沒有自然產(chǎn)能，壓裂改造后產(chǎn)油120 m3/d，產(chǎn)氣110 m3/d，其中壓裂液返排率高達(dá)90%，壓裂效果較好，不論是壓裂前、后產(chǎn)能的變化還是壓裂液返排率，均驗證了水力壓裂效果綜合評價技術(shù)在儲層壓裂效果評價中的正確性及有效性。由于陣列聲波測井儀器采集記錄時間有限，結(jié)合該井聲波時差信息，該井聲波遠(yuǎn)探測成像只能探測到井外10～15 m，因此，無法對該深度范圍之外區(qū)域進(jìn)行評價。

3 結(jié)論

1)陣列聲波測井是海上儲層壓裂效果評價的重要手段，應(yīng)用結(jié)果表明，當(dāng)?shù)貙舆M(jìn)行射孔、壓裂等儲層改造作業(yè)時，若形成單一方向的裂縫，通過對比壓裂前后的時差各向異性即可較好地評價裂縫高度，若形成的是網(wǎng)狀縫，那么利用各向異性技術(shù)就難以進(jìn)行壓裂效果評價。

2)為了滿足海上儲層壓裂效果評價需求，提出了一種基于陣列聲波測井資料的水力壓裂效果綜合評價技術(shù)。通過數(shù)值模擬，壓裂后的聲波時差有所增大、橫波幅度明顯降低。因此，可利用壓裂前、后聲波時差、橫波幅度的差異、聲波各向異性及縱波走時層析成像技術(shù)，來定性評價近井壁處儲層壓裂縫高度，同時結(jié)合偶極聲波遠(yuǎn)探測成像技術(shù)，在一定程度上評價遠(yuǎn)井壁處壓裂縫展布深度，通過“遠(yuǎn)近結(jié)合”的方式實現(xiàn)對儲層壓裂縫高度及壓裂程度的綜合評價。

3)通過采用基于陣列聲波測井資料的水力壓裂效果綜合評價技術(shù)處理渤海油田某開發(fā)井資料，表明在各向異性反演等方法失效的情況下，該技術(shù)在固井質(zhì)量、擴徑等因素干擾的影響下，依然能夠評價儲層壓裂縫高度，最終結(jié)合現(xiàn)場壓裂產(chǎn)能及壓裂液返排信息，驗證了該技術(shù)的正確性及有效性。