胡華鋒,胡 起,林正良

(中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

異常地層壓力分布是含油氣沉積盆地普遍存在的現(xiàn)象,有效預(yù)測地層超壓對于了解油氣富集、運移狀況和保證鉆井安全具有重要意義[1]。導(dǎo)致地層異常高壓的原因有很多,可以總結(jié)為三大類:第一類是欠壓實作用引起地層高壓,即當(dāng)?shù)貙映练e速度過快時,地層中的孔隙流體不能及時排出,導(dǎo)致孔隙流體壓力增加。第二類是孔隙流體體積膨脹引起超壓,即在水受熱、油氣生成、油裂解成氣、蒙脫石向伊利石轉(zhuǎn)化等因素影響下,孔隙流體體積發(fā)生膨脹,從而引起地層壓力升高。此外,構(gòu)造作用、地層抬升和地層剝蝕等因素也會引起地層壓力的變化[2]。

從20世紀(jì)60年代開始,就有學(xué)者嘗試對異常地層壓力進行預(yù)測[3-6]。有效應(yīng)力原理[7]是地層壓力預(yù)測的理論基礎(chǔ),即上覆地層壓力與巖石骨架上的垂直有效應(yīng)力和孔隙流體壓力之和形成平衡。由于地層欠壓實可以使得有效應(yīng)力不增加或者增加緩慢,但不會引起有效應(yīng)力的降低,而孔隙流體體積膨脹可以引起有效應(yīng)力的降低,進而導(dǎo)致速度曲線的倒轉(zhuǎn),因此,針對欠壓實成因提出的地層壓力預(yù)測方法通常假設(shè)淺層泥巖段為正常壓實,通過擬合淺層泥巖段速度趨勢構(gòu)建常壓趨勢線,結(jié)合Eaton方程、Fillippone方程等,計算地層壓力。在頁巖氣儲層中,泥頁巖段由于生烴增壓的影響通常不是常壓,因此常規(guī)地層壓力預(yù)測方法并不適合流體體積膨脹成因的地層壓力預(yù)測。

隨著四川盆地五峰-龍馬溪組頁巖氣藏勘探開發(fā)的不斷深入,人們逐漸認(rèn)識到地層超壓對于頁巖氣富集具有重要意義。吳波等[8]通過綜合應(yīng)用主頻較高的疊后阻抗體和地震層速度體,得到高低頻閉合速度體,最終利用Fillippone壓力預(yù)測模型計算出頁巖氣儲層的地層孔隙壓力。李金磊[9]利用高角度裂縫密度參數(shù)修正復(fù)雜構(gòu)造區(qū)地層壓力系數(shù),并預(yù)測了涪陵焦石壩頁巖地層壓力。但是,有機質(zhì)生烴是頁巖氣形成的重要地質(zhì)過程,也是頁巖氣儲層形成超壓的重要原因,常規(guī)基于欠壓實成因的地層壓力預(yù)測方法并不適合直接應(yīng)用于頁巖氣儲層的地層壓力預(yù)測。本文以Eaton方程為基礎(chǔ),考慮干酪根生烴的影響,基于巖石物理建模的思路,提出一種新的正常壓實趨勢線(normal compaction trend,NCT)構(gòu)建方法,形成新的頁巖氣儲層地層壓力預(yù)測模型,并將其與地震地層壓力預(yù)測相結(jié)合,有效提高頁巖氣儲層的鉆前地層壓力預(yù)測精度。

1 方法原理

1.1 基于測井資料的頁巖氣儲層地層壓力預(yù)測模型構(gòu)建方法

根據(jù)Terzaghi有效應(yīng)力原理,EATON[4]給出了地層壓力與地震層速度的關(guān)系:

(1)

式中:Pov為上覆地層壓力;Phy為靜水壓力;v為地震層速度;vnct為正常壓實條件下地震層速度;n為Eaton常數(shù),由實際數(shù)據(jù)擬合得到。Pov和Phy可以通過巖石密度和地層水密度在深度上積分獲得,v可以通過聲波測井或者地震資料處理獲得。該方法在實際應(yīng)用中的難點是vnct的獲取,通常的做法是假設(shè)淺層泥巖段為正常壓實,拾取淺層泥巖段速度,擬合其與深度的關(guān)系,得到一條平滑的正常壓實速度趨勢線[10],如圖1所示。但是很多情況下正常壓實趨勢線難以建立,比如深水環(huán)境中常常在泥線以下就進入高壓狀態(tài);又比如在四川盆地,淺層為灰?guī)r出露,深層泥頁巖普遍超壓。

圖1 常規(guī)正常壓實速度趨勢線

HUFFMAN等[11]分析了不同巖性的壓力加載和卸載曲線,指出不同巖性的有效應(yīng)力對速度的敏感程度不同。本文研究認(rèn)為,常規(guī)基于擬合方法得到的平滑正常壓實速度曲線是對真實正常壓實速度的一種高度近似,而正常壓實速度應(yīng)該與巖石的礦物組分、礦物磨圓度、膠結(jié)類型等有關(guān)。近年來,PERVUKHINA等[12-14]提出一種名為Clay-Plus-Silt(CPS)的泥頁巖巖石物理模型,并通過實際應(yīng)用證明了該模型可用于求取正常壓實情況下的地層速度,相對基于數(shù)據(jù)擬合得到正常壓實速度趨勢線具有更高的壓力預(yù)測精度[15]。但CPS模型是針對普通泥頁巖建立的,并未考慮生烴物質(zhì)(有機質(zhì))的影響,有機質(zhì)的存在不僅會大大降低巖石的硬度,而且其定向排列的空間分布形態(tài)往往會進一步增強巖石的各向異性,要想利用CPS模型構(gòu)建的正常壓實速度趨勢線進行富含有機質(zhì)的頁巖氣儲層地層壓力預(yù)測,需對該模型加以改進。本文通過在CPS模型中加入有機質(zhì),在一定程度上考慮有機質(zhì)生烴對異常壓力的影響,使之能更好地適用于富有機質(zhì)頁巖氣儲層的地層壓力預(yù)測。

CPS模型能計算正常聲波時差的理論基礎(chǔ)在于模型假設(shè)巖石孔隙均與粘土相關(guān),孔隙流體與粘土顆粒構(gòu)成“濕粘土”混合物。實驗數(shù)據(jù)表明,濕粘土混合物的彈性張量與濕粘土孔隙度(濕粘土中流體的體積含量)呈線性負(fù)相關(guān)關(guān)系,而與粘土的礦物組分無關(guān)。因此,濕粘土的彈性張量可由濕粘土孔隙度單一確定,整個泥頁巖的彈性張量可由濕粘土孔隙度和砂質(zhì)混合物(除粘土以外的其它礦物組分)的體積含量共同確定。與此同時,地層壓力的變化只會影響巖石中軟孔隙的開啟和關(guān)閉,而軟孔隙對總孔隙度的影響可忽略不計,故影響泥頁巖彈性張量的兩個因素均不受異常壓力的影響,從而可計算出正常壓實情況下巖石的速度或聲波時差。

本文在改進CPS模型中依然假設(shè)粘土與孔隙流體構(gòu)成濕粘土混合物,而石英、長石、方解石、黃鐵礦等硬性礦物組成砂質(zhì)混合物;與CPS模型不同的是,模型的構(gòu)成組分由原來的濕粘土-砂質(zhì)混合物兩相變?yōu)闈裾惩?砂質(zhì)混合物-有機質(zhì)三相。相比原模型中采用微分等效介質(zhì)理論(DEM模型)來求取兩相混合物的等效彈性張量,改進模型中選用Backus平均公式[16]來求取濕粘土-砂質(zhì)混合物-有機質(zhì)三相所構(gòu)成的等效介質(zhì)。由于Backus平均公式具有顯式表達(dá)形式,比需要迭代求解的DEM模型具有更高的計算效率,因此被很多學(xué)者應(yīng)用于含有機質(zhì)頁巖的巖石物理建模流程中[17-20]。

本文基于改進CPS模型的頁巖氣儲層單井地層壓力預(yù)測流程如圖2所示,具體實現(xiàn)步驟如下。

圖2 基于改進CPS模型的頁巖氣儲層單井地層壓力預(yù)測流程

1) 用經(jīng)驗公式計算濕粘土彈性張量:

(2)

(3)

2) 用Voight-Reuss-Hill模型計算砂質(zhì)混合物彈性張量:

(4)

(5)

3) 用Backus平均公式計算濕粘土-砂質(zhì)混合物-有機質(zhì)構(gòu)成的等效頁巖彈性張量:

(6a)

(6b)

C12=C11-〈c11〉+〈c12〉

(6c)

(6d)

(6e)

C66=〈c66〉

(6f)

式中:Cij表示等效頁巖的彈性剛度分量,cij為每一相的彈性剛度分量,尖括號〈〉表示對其內(nèi)屬性按體積比進行加權(quán)平均。

4) 將等效頁巖的張量元素轉(zhuǎn)化為正常壓實下的速度:

(7)

將公式(7)帶入公式(1),計算地層壓力。

1.2 頁巖氣儲層地震地層壓力預(yù)測

由于正常壓實趨勢線通常難以獲得,因此,在利用地震資料進行地層壓力預(yù)測時多采用Fillippone法[5]。該方法通過獲取高精度的速度場來盡可能提高基于地震資料的鉆前地層壓力預(yù)測精度,但地震資料處理過程中獲得的速度精度有限,無法滿足頁巖氣儲層甜點預(yù)測的精度要求。石萬忠等[21]指出,充分利用地震資料解釋過程中獲得的信息(包括地震屬性及疊前AVO反演獲得的速度資料),可以有效提高地震地層壓力預(yù)測的精度。

RASOLOFOSAON等[22]基于有效應(yīng)力原理,提出了基于波阻抗的地層壓力預(yù)測方法。該方法的優(yōu)勢如下:①地層壓力與縱波阻抗具有良好的相關(guān)性;②基于地震資料獲取縱波阻抗的技術(shù)成熟、手段豐富、分辨率高。本文基于上述方法,結(jié)合1.1節(jié)介紹的頁巖氣儲層地層壓力預(yù)測模型,建立了以疊后波阻抗反演為基礎(chǔ)的地震地層壓力預(yù)測流程,如圖3所示。具體實現(xiàn)步驟如下:

1) 基于常規(guī)波阻抗反演獲得高精度波阻抗數(shù)據(jù)Ip。

2) 基于1.1節(jié)單井計算獲得的正常壓實速度vnct,結(jié)合地質(zhì)層位信息,利用克里金插值方法構(gòu)建正常壓實速度趨勢體。

3) 計算3D密度體ρ,并計算上覆地層壓力Pov和靜水壓力Phy。

4) 將Ip,vnct,ρ,Pov,Phy代入Eaton壓力預(yù)測模型,計算地層壓力及壓力系數(shù)。

(8)

圖3 頁巖氣儲層地震地層壓力預(yù)測流程

2 實際應(yīng)用

四川盆地五峰—龍馬溪組頁巖層系普遍存在超壓現(xiàn)象,某頁巖氣區(qū)塊A井五峰-龍馬溪組實測壓力系數(shù)區(qū)間為1.20～1.55[23]。以A井實測壓力系數(shù)范圍為約束,采用不同的方法構(gòu)建了正常壓實速度趨勢線,如圖4a所示。其中紅色曲線(NCT-model)為基于本文改進的CPS模型,通過巖石物理模型正演計算得到的正常壓實速度曲線;綠色實線(NCT-fit)為通過常規(guī)數(shù)據(jù)擬合方法得到的正常壓實速度曲線;藍(lán)色曲線為泥巖實際聲波速度。對比可見,通過巖石物理模型構(gòu)建的正常壓實速度曲線能夠更好地反映巖性細(xì)微變化對巖石骨架應(yīng)力的影響,這種巖性差異會導(dǎo)致正常壓實速度的細(xì)微變化。圖4b為A井靜水壓力(黑色實線)、上覆地層壓力(品紅色實線)及地層壓力計算結(jié)果,其中綠色曲線為基于NCT-fit曲線結(jié)合Eaton方程計算得到的地層壓力,紅色曲線為基于NCT-model曲線結(jié)合Eaton方程計算得到的地層壓力。圖4c為兩種不同方法計算的壓力系數(shù)(壓力系數(shù)=地層壓力/靜水壓力)?？梢钥闯?即使在最佳擬合的情況下,采用NCT-fit曲線計算出的地層壓力曲線依然有部分?jǐn)?shù)據(jù)點偏離了實測壓力系數(shù)區(qū)間,而基于NCT-model曲線計算的地層壓力曲線數(shù)據(jù)點變化范圍與實測壓力區(qū)間非常匹配。

圖4 A井地層壓力預(yù)測結(jié)果a 改進的CPS模型計算的NCT與擬合方法計算的NCT對比; b 基于兩種方法求取的地層壓力; c 基于兩種方法求取的地層壓力系數(shù)

以A井為模型校驗井,確定頁巖氣儲層地層壓力預(yù)測模型相關(guān)參數(shù),并將該模型應(yīng)用到區(qū)塊內(nèi)B井地層壓力預(yù)測中。圖5a為B井目的層段靜水壓力(藍(lán)色實線)、上覆地層壓力(紅色實線)和地層壓力(品紅色實線)計算結(jié)果。圖5b對比了地層壓力系數(shù)與鉆井液泥漿比重,可以看到壓力系數(shù)與泥漿比重具有較好的相關(guān)性,在2320m處地層壓力系數(shù)增大,對應(yīng)的泥漿比重也增大,說明通過A井構(gòu)建的地層壓力預(yù)測模型能夠在全區(qū)塊很好地推廣應(yīng)用。

W1井區(qū)是四川盆地另一頁巖氣新探區(qū),由于實測壓力數(shù)據(jù)缺乏,利用傳統(tǒng)方法建立壓實趨勢線具有較大的不確定性,也無法利用實測壓力數(shù)據(jù)來回歸其它經(jīng)驗公式法中的經(jīng)驗系數(shù),給地層壓力預(yù)測帶來了很大困難。采用本文提出的改進CPS模型對該區(qū)進行了壓力預(yù)測,圖6為W1井的壓力預(yù)測結(jié)果,圖6a中藍(lán)線、綠線和紅線分別表示計算出的靜水壓力、地層壓力和上覆地層壓力,圖6b給出了預(yù)測壓力系數(shù)?？梢钥闯?目的層優(yōu)質(zhì)頁巖段預(yù)測地層壓力最高可達(dá)66.97MPa,對應(yīng)壓力系數(shù)為1.92,高于上覆普通泥頁巖段和下伏碳酸鹽巖段。該井目的層實測壓力為68.69MPa,對應(yīng)實測壓力系數(shù)為1.94,說明本文給出的單井地層壓力預(yù)測方法具有較高的預(yù)測精度。

圖5 B井地層壓力預(yù)測結(jié)果a 地層壓力; b 壓力系數(shù)與泥漿比重

利用W1井構(gòu)建的地層壓力預(yù)測模型,對該井區(qū)進行鉆前地震地層壓力預(yù)測,圖7為W1井區(qū)地震地層壓力系數(shù)剖面,可見五峰-龍馬溪組底部(黑色實線)地層壓力系數(shù)明顯偏高,介于1.70～1.96之間,屬于超壓頁巖氣藏,且深凹區(qū)壓力系數(shù)更大,保存條件相對隆起區(qū)要好。對該區(qū)W1,W2,W3井目的層地層壓力系數(shù)預(yù)測結(jié)果進行了統(tǒng)計(表1),誤差小于6%,地震地層壓力系數(shù)預(yù)測結(jié)果與實測地層壓力系數(shù)吻合較好。

圖6 W1井地層壓力預(yù)測結(jié)果a 地層壓力; b 壓力系數(shù)

圖7 W1井區(qū)地震地層壓力系數(shù)剖面

表1 W1井區(qū)3口井五峰-龍馬溪組頁巖氣儲層地層壓力預(yù)測結(jié)果統(tǒng)計

3 結(jié)束語

本文基于頁巖氣儲層異常高壓成因分析和CPS模型,考慮干酪根生烴的影響,提出一種新的正常壓實速度趨勢線構(gòu)建方法,結(jié)合Eaton方程,形成新的頁巖氣儲層地層壓力預(yù)測模型。以該模型為基礎(chǔ),結(jié)合有效應(yīng)力原理,建立了基于波阻抗的頁巖氣儲層地層壓力地震預(yù)測技術(shù)流程。在四川盆地多個頁巖氣區(qū)塊的實際應(yīng)用表明:①該流程以測井解釋的礦物組分、孔隙度等參數(shù)為基礎(chǔ),基于巖石物理建模構(gòu)建頁巖氣儲層壓力預(yù)測模型,通常只需要知道目的層段地層壓力系數(shù),就能構(gòu)建出比較合理的模型,因此能更好地適應(yīng)頁巖氣儲層實測壓力數(shù)據(jù)較少的情況。②通過巖石物理建模構(gòu)建頁巖氣儲層壓力預(yù)測模型,能更有效地反映巖石骨架有效應(yīng)力隨巖性的細(xì)微變化,結(jié)合高精度波阻抗反演結(jié)果,能有效提高頁巖氣儲層鉆前地層壓力預(yù)測精度。

參 考 文 獻(xiàn)

[1] 張冰,符力耘,魏偉,等.異常地層壓力聲波衰減分析及地層壓力地震數(shù)據(jù)外推反演[J].地球物理學(xué)報,2014,57(10):3373-3388

ZHANG B,FU L Y,WEI W,et al.Acoustic wave attenuation analysis of the abnormal geo-pressure and seismic extrapolation inversion of the formation pressure[J].Chinese Journal of Geophysics,2014,57(10):3373-3388

[2] 孫武亮,孫開峰.地震地層壓力預(yù)測綜述[J].勘探地球物理進展,2007,30(6):428-432

SUN W L,SUN K F.Summarization of geopressure prediction using seismic data[J].Progress in Exploration Geophysics,2007,30(6):428-432

[3] PENNEBAKER E S J R.Seismic data indicate depth,magnitude of abnormal pressure[J].World Oil,1968,166:73-78

[4] EATON B A.Graphical method predicts geopressure worldwide[J].World Oil,1972,182:51-56

[5] FILLIPPONE W R.Estimation of formation parameters and the prediction of overpressure from seismic data[J].Expanded Abstracts of 52ndAnnual Internat SEG Mtg,1982:502-503

[6] DUTTA N C.Geopressure prediction using seismic data:current status and the road ahead[J]:Geophysics,2002,67(6):2012-2041

[7] TERZAGHI K.Theoretical soil mechanics[M].New York:John Wiley and Sons,1943:528-530

[8] 吳波,王薦,潘樹林,等.基于高低頻速度閉合技術(shù)的地層壓力預(yù)測[J].石油物探,2017,56(4):575-580

WU B,WANG J,PAN S L,et al.Formation pressure prediction based on a closed velocity body by merging the high frequency velocity with the low frequency velocity[J].Geophysical Prospecting for Petroleum,2017,56(4):575-580

[9] 李金磊.涪陵焦石壩頁巖氣層壓力預(yù)測技術(shù)研究[J].石油物探,2017,56(4):567-574

LI J L.Pressure prediction of Jiaoshiba shale gas reservoir in the Fuling Shale Gasfield[J].Geophysical Prospecting for Petroleum,2017,56(4):567-574

[10] SUWANNASRI K,PROMRAK W,UTITSAN S,et al.Reducing the variation of Eaton’s exponent for overpressure prediction in a basin affected by multiple overpressure mechanisms[J].Expanded Abstracts of 84thAnnual Internat SEG Mtg,2014:SB57-SB68

[11] HUFFMAN A R,BOWERS G L.The future of pressure prediction using geophysical methods[J].AAPG Memoir,2002,76:217-233

[12] PERVUKHINA M,DEWHURST D,GUREVICH B,et al.Stress-dependent elastic properties of shales:measurement and modeling[J].The Leading Edge,2008,27(6):772-779

[13] PERVUKHINA M,GOLODONIUC P,GUREVICH B,et al.An estimation of sonic velocities in shale using clay and silt fractions from the elemental capture spectroscopy log[J].Extended Abstract of 74thEAGE Conference & Exhibition,2012:4-7

[14] PERVUKHINA M,CLAUDIO D P,DEWHURST D N,et al.An estimation of pore pressure in shales from sonic velocities[J].Expanded Abstracts of 83rdAnnual Internat SEG Mtg,2013:2653-2657

[15] HAN T,PERVUKHINA M,CLENNELL M B.Overpressure prediction in shales[J].Extended Abstract of 76thEAGE Conference & Exhibition,2014:16-19

[16] BACKUS G E.Long-wave elastic anisotropy produced by horizontal layering[J].Journal of Geophysical Research,1962,67(11):4427-4440

[17] SAYERS C M.The effect of kerogen on the elastic anisotropy of organic-rich shales[J].Geophysics,2013,78(2):D65-D74

[18] GUO Z Q,CHAPMAN M,LI X Y.A shale rock physics model and its application in the prediction of brittleness index,mineralogy,and porosity of the Barnett Shale[J].Expanded Abstracts of 82ndAnnual Internat SEG Mtg,2012:1-5

[19] VERNIK L,NUR A.Ultrasonic velocity and anisotropy of hydrocarbon source rocks[J].Geophysics,1992,57(5):727-735

[20] VERNIK L,LANDIS C.Elastic anisotropy of source rocks,implications for hydrocarbon generation and primary migration[J].AAPG Bulletin,1996,80(4):531-544

[21] 石萬忠,何生,陳紅漢.多地震屬性聯(lián)合反演在地層壓力預(yù)測中的應(yīng)用[J].石油物探,2006,45(6):580-585

SHI W Z,HE S,CHEN H H.Application of joint inverse modeling of multi-seismic-attributes in predicting formation pressure[J].Geophysical Prospecting for Petroleum,2006,45(6):580-585

[22] RASOLOFOSAON P,TONELLOT T.Method for quantitative evaluation of fluid pressures and detection of overpressures in an underground medium:US 20090093963A1[P].2008-10-02

[23] 郭旭升.涪陵頁巖氣田礁石垻區(qū)塊富集機理與勘探技術(shù)[M].北京:科學(xué)出版社,2014:267-271

GUO X S.Accumulation mechanism and exploration technology for the Jiaoshiba block of Fuling shale gas field [M].Beijing:Science Press,2014:267-271