田 佳 張 勇 胡佳男

(1. 西安石油大學 地球科學與工程學院/陜西省油氣成藏地質(zhì)學重點實驗室， 西安 710065；2. 吉林油田公司英臺采油廠， 吉林 松原 138000；3. 吉林油田公司扶余采油廠， 吉林 松原 138000)

0 前 言

松遼盆地屬于東部陸相含油氣盆地，油氣資源類型多、分布廣、儲量豐富、產(chǎn)能較高。經(jīng)過長期的勘探與開發(fā)，其厚層三角洲支河道砂體富集的含油氣層系逐漸被摸清，目前有待進一步挖掘致密層系的開發(fā)潛力[1]。新民油田位于松遼盆地南部扶新隆起帶北部斜坡帶，以低孔 — 低滲型儲層為主，儲層物性較差，非均質(zhì)性嚴重，需要采取水力壓裂等增產(chǎn)措施以獲取商業(yè)油氣流[2]。

近年來，微地震監(jiān)測技術(shù)被廣泛應(yīng)用于致密油氣藏的水力壓裂監(jiān)測、礦山巖爆監(jiān)測、二氧化碳地質(zhì)封存監(jiān)測及常規(guī)油氣田的注采誘發(fā)地震監(jiān)測[3-5]。尤其是在致密油氣藏開發(fā)中，微地震水力壓裂監(jiān)測是唯一能夠?qū)毫芽p幾何產(chǎn)狀實現(xiàn)成像的遠場技術(shù)[6-7]。目前，可以采用地面或井下實時監(jiān)測的微地震數(shù)據(jù)與背景噪聲壓制、有效事件自動拾取、震源定位與震源機制反演等數(shù)據(jù)處理技術(shù)，估算水力壓裂縫的幾何產(chǎn)狀、幾何規(guī)模、增產(chǎn)體積、地應(yīng)力變化和網(wǎng)格連通性等參數(shù)，并據(jù)此確定壓裂井的最優(yōu)射孔區(qū)間，調(diào)整壓裂施工方案[8-9]。

本次研究優(yōu)選了新民油田3口新鉆井，針對扶余油層進行水力壓裂，采用井下微地震監(jiān)測壓裂過程中壓裂縫的幾何產(chǎn)狀，擬明確目標區(qū)域的最大水平主應(yīng)力方位，分析壓裂施工參數(shù)對裂縫擴展的影響因素，以加深對裂縫擴展規(guī)律和裂縫延伸形態(tài)的認識。

1 研究區(qū)概況

研究區(qū)塊位于新民油田東壘南部，為由南東向北西傾沒的單斜構(gòu)造。區(qū)塊內(nèi)斷層發(fā)育較好，砂體以粉砂巖為主，平均鉆遇砂巖厚度約38.0 m，有效厚度約4.4 m。區(qū)塊內(nèi)開發(fā)的目的層主要為扶余油層，其儲層平均滲透率為1.5×10-3μm2，物性較差。平面上油層連續(xù)性較差，但主力油層分布相對穩(wěn)定。

優(yōu)選研究區(qū)內(nèi)3口新鉆井(民X-1井、民X-2井和民X-3井)進行水力壓裂作業(yè)，每口井設(shè)置2個壓裂層，其壓裂層均屬于扶楊層系扶余泉四段油層。在鄰井民M井布置了10級井下三分量檢波器，用于逐次監(jiān)測各井水力壓裂過程中所產(chǎn)生的微地震信號。圖1所示為這3口井的平面井位分布與井下微地震觀測系統(tǒng)設(shè)置。觀測系統(tǒng)中以監(jiān)測井的井口坐標為坐標原點，檢波器陣列中點與民X-1井、民X-2井和民X-3井最遠射孔位置之間的距離均小于300 m，以滿足井下微地震監(jiān)測條件。

圖1 研究區(qū)3口井的平面井位分布與井下微地震觀測系統(tǒng)設(shè)置

2 微地震數(shù)據(jù)處理與解釋

針對每個單井，采用圖2所示流程進行井下微地震數(shù)據(jù)處理。根據(jù)震源定位結(jié)果，分析各壓裂層段壓裂縫的長、寬、高和方位等幾何特征與裂縫規(guī)模，定量評價壓裂效果[10-11]。

圖2 井下微地震監(jiān)測數(shù)據(jù)處理流程

2.1 頻譜分析與噪聲壓制

運用常規(guī)的微地震頻譜分析法，即通過傅里葉變換將時域原始信號轉(zhuǎn)換到頻域，再根據(jù)頻域響應(yīng)特征選取有效信號的頻帶。該方法忽略了微地震瞬態(tài)信號在時域的分布特征，且易受到單頻噪聲的干擾，無法準確獲取有效信號的頻帶。為解決這一問題，在此通過S變換得到原始微地震信號的時頻譜，然后根據(jù)時頻譜中高能量區(qū)域的分布特征來確定有效信號的頻帶范圍[12]。據(jù)圖3所示地震原始信號及其對應(yīng)的時頻譜，確定了本次微地震有效信號的頻帶范圍為50～250 Hz，因此，可采用此頻帶范圍的帶通濾波對原始波形記錄進行噪聲壓制。

2.2 微地震有效事件自動拾取

目前，常采用長短時窗能量比法(STA/LTA)或者自回歸 Akaike 信息準則(AIC)法，基于信號與噪聲在時域分布或統(tǒng)計特性的差異進行自動拾取[13]。這類方法僅利用了信號在時域分布的特征，而忽略了其在頻域分布的特性，不利于對弱事件的拾取[14]。因此，采用基于疊加S變換的時頻譜拾取法實現(xiàn)對有效事件的自動識別與拾取[15]。首先，對原始微地震記錄進行S變換，得到時頻譜；其次，將時頻譜在有效信號頻帶范圍內(nèi)按照每個頻率成分進行疊加，將其疊加譜作為特征函數(shù)，以此構(gòu)造長短時窗，從而實現(xiàn)對有效事件的自動拾取。若拾取到有效事件，則采用AIC算法進一步精細提取縱波和橫波初至。

基于上述微地震原始信號和時頻譜完成了有效事件自動拾取，結(jié)果如圖4所示。其中圖4a所示為圖3中微地震原始信號的STA/LTA比值曲線及設(shè)定閾值，從中僅能分辨出S波初至；圖4b所示為圖3中時頻譜的STA/LTA比值曲線及設(shè)定閾值，從中可以拾取到清楚的P波和S波初至；圖4c所示為采用AIC算法進一步精細提取的P波和S波初至?？梢钥闯?，基于S變換疊加譜的自動拾取算法能夠識別弱事件的起跳信息。

圖3 某道微地震原始信號及其時頻譜

圖4 某道微地震有效事件自動拾取

2.3 速度模型建立及校正

基于民M井的聲波測井曲線，實現(xiàn)對壓裂目的層段的速度分層，并建立了一維P波層狀速度模型(見圖5a)。隨著油田開采的持續(xù)，儲層的壓力和空隙內(nèi)流體的性質(zhì)會逐漸發(fā)生改變，使得基于聲波測井信息建立的速度模型無法滿足微地震震源定位對精度的需求。為此，采用Pei等人提出的速度模型優(yōu)化算法[16]，基于射孔信號的P波初至建立了優(yōu)化速度模型的目標函數(shù)，如式(1)所得:

(1)

式中：E(v) —— 目標函數(shù)；

v—— 一維層狀速度向量；

Nr—— 接收器總數(shù)；

由于射孔位置已知，因此可調(diào)整各層速度以使目標函數(shù)值最小，此時的速度向量v為最優(yōu)速度。圖5b所示為基于射孔信號校正的速度模型。由于缺失橫波測井信息，S波速度模型由校正后的P波速度模型與研究區(qū)域的縱橫波比經(jīng)驗值估算而得。

圖5 一維P波層狀速度模型

2.4 偏振分析與震源定位

由于單井監(jiān)測在水平面上缺乏對震源方位的約束，導致潛在震源最終分布在以檢波器水平坐標為圓心、以有效事件縱橫波初至差估算的震源距離為半徑的圓弧上。在此，采用矢端曲線法將x分量與y分量波形分別投影到x-y平面，通過最佳線性擬合質(zhì)點位移的方位來估算P波的偏振方向，從而確定微地震有效事件的震源傳播方位。

獲取縱橫波初至與P波偏振方向后，采用概率密度坍塌網(wǎng)格搜索法進行震源定位[17]。將研究區(qū)域的速度模型離散化，通過拾取的P波、S波初至來構(gòu)建反演震源位置的非線性目標函數(shù)：

(2)

式中：L(x) —— 非線性目標函數(shù)；

N—— 接收器數(shù)量；

σP、σS—— P波、S波拾取初至的標準偏差。

坍塌網(wǎng)格搜索的策略是：首先，搜索較大網(wǎng)格的目標函數(shù)最小值；然后，再將該網(wǎng)格不斷二分，直到其滿足算法終止標準；最后，通過搜索拾取初至與理論初至殘差的全局最小值來確定震源位置。

對3口壓裂井共6個壓裂層段的原始微地震監(jiān)測數(shù)據(jù)作了定位處理。以民X-1井第1壓裂層段的處理結(jié)果為例，該壓裂層段共拾取了56個有效事件，其震源定位如圖6所示。首先，通過事件位置的空間分布，估算出壓裂縫的縫長為159.02 m，縫高為36.98 m，縫寬為50.63 m，裂縫總體延伸方位為NE63.06°。同時，由此判斷該區(qū)域的最大水平主應(yīng)力方向為北東向，所形成的有效面積約8 045 m2。此外，觀測到2組沿西北 — 東南方向延伸的分支裂縫形成的復雜共軛裂縫結(jié)構(gòu)。

圖6 民X-1井第1層水力壓裂微地震震源定位

圖7所示為民X-1井第1層壓裂縫動態(tài)擴展過程，壓裂前期事件點數(shù)較少，僅形成了一條主裂縫。裂縫整體沿北東向進行延伸，壓裂過程中事件點數(shù)不斷增多，裂縫規(guī)模不斷擴大，并構(gòu)成2條與主應(yīng)力方向垂直的共軛裂縫。為了定量評價壓裂縫網(wǎng)的復雜程度，基于震源在平面上的空間分布進行直線擬合，將擬合系數(shù)作為評價復雜程度的依據(jù)。擬合系數(shù)越大，表示縫網(wǎng)復雜指數(shù)越低，線性擬合關(guān)系如式(3)所示：

圖7 民X-1井第1層壓裂縫動態(tài)擴展過程

IFC=1-R2

(3)

式中：IFC—— 裂縫復雜指數(shù)；

R2—— 擬合系數(shù)。

經(jīng)擬合，該壓裂層的裂縫復雜指數(shù)為0.63。

3 壓裂縫幾何形態(tài)對比

將研究區(qū)域3口壓裂井產(chǎn)生的微地震震源位置在空間上進行疊放對比，對縫長、縫寬與縫高進行歸一化(見圖8)。經(jīng)統(tǒng)計，可得研究區(qū)域內(nèi)3口壓裂井的裂縫形態(tài)特征：

圖8 研究區(qū)3口壓裂井微地震裂縫規(guī)模對比

(1) 縫長。民X-3井的壓裂縫最短，民X-2井次之，民X-1井的壓裂縫最長。

(2) 縫寬。民X-3井的壓裂縫最寬，民X-2井次之，民X-1井的壓裂縫最窄。

(3) 縫高。民X-3井和民X-2井高度相似，均小于民X-1井。

(4) 裂縫復雜指數(shù)。民X-3井裂縫復雜指數(shù)最高，民X-2井次之，民X-1井的裂縫復雜指數(shù)最低。

總體上，沿著西北 — 東南方向，3口井的壓裂縫規(guī)模遞增，裂縫復雜指數(shù)遞減。其中，民X-3井的壓裂縫規(guī)模較小，裂縫復雜指數(shù)較高，形成粗短的復雜縫網(wǎng)結(jié)構(gòu)；民X-1井的壓裂縫規(guī)模較大，裂縫復雜指數(shù)較低，形成細長的條帶狀縫網(wǎng)結(jié)構(gòu)，沿著北 — 東方向不斷延伸；民X-2井的壓裂縫規(guī)模與裂縫復雜程度則處在民X-3井與民X-1井之間。

4 壓裂縫幾何形態(tài)影響因素分析

上述3口壓裂井總計6個壓裂層段均屬于扶余油層的泉4段，巖性基本相同，主要由泥巖與粉砂巖、細砂巖、泥質(zhì)粉砂巖、粉砂質(zhì)泥巖互層組成。為進一步探索影響壓裂縫幾何形態(tài)的因素，對這3口壓裂井的水力壓裂施工參數(shù)(加砂總量與累計液量)與裂縫形態(tài)特征作了定量對比分析，并總結(jié)壓裂施工參數(shù)對壓裂縫延伸形態(tài)的影響規(guī)律。

根據(jù)第1、2層水力壓裂期間的加砂總量、累計液量、最大壓力統(tǒng)計對比結(jié)果(見圖9)，以及水力壓裂期間的加砂總量、累計液量與裂縫規(guī)模、裂縫復雜指數(shù)的關(guān)系(見圖10)，得到以下認識：

圖9 研究區(qū)3口壓裂井加砂總量、累計液量、最大壓力統(tǒng)計對比

圖10 研究區(qū)3口壓裂井加砂總量、累計液量與裂縫規(guī)模、裂縫復雜指數(shù)的關(guān)系

(1) 壓裂縫的規(guī)模與累計液量、加砂總量之間具有明顯的正相關(guān)性。

(2) 隨著累計液量、加砂總量的增加，壓裂縫長度也明顯增大，而壓裂縫高度增幅不明顯。這表明該區(qū)域的垂向應(yīng)力大于水平應(yīng)力，而水力壓裂主要提供的是剪切方向上的應(yīng)力，因此加大排量和加砂總量無法使裂縫高度明顯增大。

(3) 壓裂縫的規(guī)模與最大壓力之間無明顯的相關(guān)性，因為最大壓力只起到破裂地層的作用，當?shù)貙悠屏押笾瘟芽p的開啟主要是受排量和加砂總量的影響。

由于加砂總量、累計液量與縫長、縫高具有較明顯的正相關(guān)性，因此可通過對壓裂施工參數(shù)與壓裂縫規(guī)模的最小二乘擬合推導出壓裂縫長和縫高的預(yù)測公式：

(4)

式中：lf—— 預(yù)測的裂縫長度；

hf—— 預(yù)測的裂縫高度；

Vf—— 壓裂期間的累計液量；

Vs—— 加砂總量。

由式(4)可知，裂縫長度、裂縫高度總體上與累計液量、加砂總量具有較強的線性正相關(guān)性。由于裂縫寬度、裂縫復雜指數(shù)與累計液量、加砂總量的相關(guān)性不顯著，且研究區(qū)壓裂數(shù)據(jù)規(guī)模有限，因此無法得到明確的擬合關(guān)系。

5 結(jié) 語

通過井下微地震監(jiān)測，對松遼盆地新民油田南部的3口壓裂井的裂縫幾何形態(tài)特征作了定量對比，對壓裂縫形態(tài)、復雜指數(shù)與壓裂施工參數(shù)之間的關(guān)系作了定性分析，并推導出研究區(qū)內(nèi)基于壓裂施工參數(shù)的壓裂縫長、縫高的預(yù)測關(guān)系式。經(jīng)過分析，得到以下認識：

(1) 研究區(qū)內(nèi)的水力壓裂縫主體沿著NE66.3° — NE77°方向延伸，最大水平主應(yīng)力方向為北東向。

(2) 研究區(qū)內(nèi)的水力壓裂縫基本沿最大主應(yīng)力方向擴展，除沿最小主應(yīng)力延伸的分支裂縫外，未見裂縫主體方位發(fā)生大的轉(zhuǎn)變。這意味著在壓裂地區(qū)，壓裂縫受天然裂縫的影響較小。

(3) 壓裂縫整體呈近東西向的條帶狀裂縫，但其內(nèi)部分布著沿北西 — 南東方向延伸的分支裂縫和沿北東 — 南西方向延伸的分支裂縫，形成網(wǎng)狀交叉裂縫。

(4) 民X-1井、民X-2的第1層壓裂段微地震事件在深度上有重疊，存在重復壓裂現(xiàn)象，為此需要進一步調(diào)整優(yōu)化壓裂施工參數(shù)。

(5) 研究區(qū)內(nèi)加砂總量、累計液量的增加都可使裂縫的長度、高度增大，但是對于裂縫的寬度并無影響。