張殿偉 孫 煒 李雙建 郝運(yùn)輕 劉 玲

中國(guó)石化石油勘探開(kāi)發(fā)研究院

0 引言

四川盆地涪陵頁(yè)巖氣田的頁(yè)巖氣勘探開(kāi)發(fā)實(shí)踐表明，保存條件的好壞是影響頁(yè)巖氣資源富集程度的關(guān)鍵因素[1-5]。針對(duì)頁(yè)巖氣儲(chǔ)層的保存條件的研究和評(píng)價(jià)工作，前人從地質(zhì)及構(gòu)造分析的角度做了大量的研究。魏祥峰等[6]通過(guò)研究建立了四川盆地及周緣復(fù)雜構(gòu)造區(qū)下古生界海相頁(yè)巖氣5大類、28項(xiàng)參數(shù)的保存條件綜合評(píng)價(jià)指標(biāo)體系，并指出壓力系數(shù)可在一定程度上指示保存條件的優(yōu)劣；董清源等[7]綜合成藏條件、構(gòu)造特征、蓋層條件、地層水特征及頁(yè)巖氣組分等因素對(duì)上奧陶統(tǒng)五峰組—下志留統(tǒng)龍馬溪組頁(yè)巖氣保存條件進(jìn)行了探討，并建立了保存條件評(píng)價(jià)指標(biāo)，優(yōu)選了保存條件有利區(qū)；何順等[8]從構(gòu)造因素、頂?shù)装鍡l件、埋深、地層壓力系數(shù)、水文地質(zhì)及氣體組分條件等非構(gòu)造因素出發(fā)，分析了川東南丁山地區(qū)頁(yè)巖氣的保存條件。

從前人的研究成果來(lái)看，目前關(guān)于保存條件的研究主要是從地質(zhì)因素考慮建立保存條件評(píng)價(jià)標(biāo)準(zhǔn)，并優(yōu)選保存條件有利區(qū)，但是利用地球物理方法預(yù)測(cè)保存條件特征的研究則相對(duì)較少。為此，筆者從頁(yè)巖氣保存條件的兩項(xiàng)指示參數(shù)——裂縫和地層壓力系數(shù)出發(fā)，利用地球物理方法對(duì)頁(yè)巖氣層的保存條件進(jìn)行評(píng)價(jià)，形成了一套能夠有效預(yù)測(cè)頁(yè)巖氣保存條件分布特征的新方法，并基于涪陵頁(yè)巖氣田焦石壩區(qū)塊五峰組—龍馬溪組頁(yè)巖氣井的生產(chǎn)數(shù)據(jù)，對(duì)比探討了新方法的適用性。

1 方法原理

1.1 基于方位各向異性的裂縫預(yù)測(cè)方法

泥頁(yè)巖中的裂縫依其產(chǎn)狀可分為水平縫（或低角度縫）與高角度裂縫。水平縫一般沿著裂縫的順層薄弱面發(fā)育，有利于頁(yè)巖氣的保存；高角度裂縫的密集發(fā)育，容易造成頁(yè)巖氣的散失，是保存條件被破壞的重要特征。從目前的地球物理預(yù)測(cè)技術(shù)現(xiàn)狀來(lái)看，水平縫的預(yù)測(cè)技術(shù)仍不成熟，相比之下，基于方位各向異性的高角度裂縫預(yù)測(cè)技術(shù)應(yīng)用范圍較廣。因此，筆者利用基于方位各向異性的裂縫預(yù)測(cè)技術(shù)描述高角度縫的分布特征，作為保存條件是否被破壞的依據(jù)。

基于具有水平對(duì)稱軸的橫向各向同性（Transversely Isotropy with a Horizontal Axis of Symmetry，簡(jiǎn)稱HTI）介質(zhì)的方位各向異性裂縫預(yù)測(cè)技術(shù)，目前應(yīng)用最廣泛的方法是Rüger[9-10]基于弱各向異性的概念，結(jié)合Thomsen的各向異性系數(shù)，在小入射角的前提下得到各向異性介質(zhì)中反射系數(shù)隨方位角和入射角變化的公式，即

式中θ表示入射角，(e)；表示方位角（圖1），(e)；RP(θ,)表示與入射角θ和方位角相關(guān)的縱波反射系數(shù)；Z=ρα表示縱波波阻抗，g/cm3gm/s；ρ表示介質(zhì)密度，g/cm3；α表示縱波速度，m/s；G=ρβ2表示剪切模量，Pa；β表示橫波速度，m/s；γ、δ、ε表示Thomsen各向異性系數(shù)；、分別表示上、下界面物理量的差值、均值，其中分別表示α、β、γ、δ、ε、G、Z。

圖1 縱波在HTI介質(zhì)模型中的入射示意圖

Rüger在小入射角的前提下，對(duì)式（1）進(jìn)一步簡(jiǎn)化，并引入了各向同性項(xiàng)Biso和各向異性項(xiàng)Bani，將地震反射系數(shù)變?yōu)锽iso、Bani以及裂縫方位角的關(guān)系式[10-11]，即

基于式（2），可以對(duì)方位各向異性強(qiáng)度進(jìn)行求取。具體的實(shí)現(xiàn)方法一般包括兩種：①利用3個(gè)方位或3個(gè)方位以上的地震振幅數(shù)據(jù)，對(duì)3個(gè)未知量進(jìn)行超定方程組的求解[12]，將求出的各向同性項(xiàng)Biso和各向異性項(xiàng)Bani進(jìn)行比值運(yùn)算，即可得到方位各向異性強(qiáng)度，同時(shí)求解出的方位角為裂縫的對(duì)稱軸方向；②對(duì)原始方位角道集進(jìn)行分方位角疊加、偏移處理，然后在不同方位角的地震數(shù)據(jù)上進(jìn)行地震屬性求取，最后在方位角—地震屬性的坐標(biāo)系中進(jìn)行方位各向異性橢圓的擬合，橢圓的長(zhǎng)軸或短軸方向代表著裂縫的走向，且裂縫密度越大、各向異性越強(qiáng)，橢圓的扁率越小，橢圓越扁[12-17]。

利用上述兩種方法，均可實(shí)現(xiàn)對(duì)高角度裂縫的裂縫發(fā)育程度以及裂縫方向的預(yù)測(cè)，筆者采用方法②以實(shí)現(xiàn)對(duì)各向異性強(qiáng)度及方向的求取。

1.2 基于疊前縱波波阻抗反演的地層壓力預(yù)測(cè)方法

地層壓力即孔隙流體壓力，是由巖石孔隙中地層水、石油、天然氣等地質(zhì)流體作用產(chǎn)生的壓力。在任何地質(zhì)背景下，正常地層壓力與從地表到目的層的靜水柱壓力相等，偏離正常壓力趨勢(shì)線的壓力被認(rèn)為是異常地層壓力。

在對(duì)焦石壩區(qū)塊龍馬溪組頁(yè)巖氣儲(chǔ)層的研究中發(fā)現(xiàn)，異常地層壓力是頁(yè)巖氣富集高產(chǎn)的重要標(biāo)志[1-5]。地層壓力出現(xiàn)異常高壓現(xiàn)象的根本原因是天然氣在封閉空間（如有機(jī)孔等）內(nèi)的大量聚集，異常高壓說(shuō)明頁(yè)巖儲(chǔ)層的保存條件較好。

目前，在勘探階段針對(duì)地層壓力評(píng)價(jià)的方法主要是綜合地震數(shù)據(jù)和地層壓力計(jì)算公式進(jìn)行區(qū)域地層壓力計(jì)算，并進(jìn)一步結(jié)合靜水壓力實(shí)現(xiàn)地層壓力系數(shù)的預(yù)測(cè)。地層壓力計(jì)算公式有等效深度法、Eaton公式法及其改進(jìn)方法和Fillippone公式法及其改進(jìn)方法等。等效深度法和Eaton公式法中的正常壓實(shí)趨勢(shì)線往往難于確定，而Fillippone公式法雖然避免了建立正常壓實(shí)趨勢(shì)線，但是，其地層壓力預(yù)測(cè)精度受控于求取的層速度精度的影響，分辨率較低。

筆者提出了基于疊前縱波波阻抗反演的新的地層壓力系數(shù)預(yù)測(cè)方法（以下簡(jiǎn)稱波阻抗反演法），如式（3）：

式中ap表示地層壓力系數(shù)，無(wú)量綱；p地層表示地層壓力，MPa；p靜水表示地層靜水柱壓力，MPa；p視巖表示基于縱波阻抗構(gòu)建的地層壓力，MPa；p靜巖表示地層靜巖壓力，MPa；A表示壓力校正系數(shù)，無(wú)量綱。

對(duì)于有實(shí)測(cè)地層壓力數(shù)據(jù)的井段，可以根據(jù)測(cè)井曲線計(jì)算出靜水柱壓力值、靜巖壓力值和視巖壓力值，結(jié)合實(shí)測(cè)地層壓力可得到對(duì)應(yīng)的校正系數(shù)A；再分別計(jì)算出靜水柱壓力數(shù)據(jù)體、靜巖壓力數(shù)據(jù)體和視巖壓力數(shù)據(jù)體，結(jié)合式（3）可計(jì)算目的層段的地層壓力系數(shù)數(shù)據(jù)體，從而描述目的層地層壓力的分布特征。筆者研究中將新公式與傳統(tǒng)Fillippone法進(jìn)行了單井預(yù)測(cè)精度的對(duì)比，結(jié)果表明新公式的精度更高，從而證實(shí)了該方法的有效性（圖2）。

圖2 JY1井地層壓力系數(shù)預(yù)測(cè)與實(shí)測(cè)結(jié)果對(duì)比圖

2 應(yīng)用實(shí)例

2.1 研究區(qū)地質(zhì)背景

涪陵頁(yè)巖氣田焦石壩構(gòu)造屬于四川盆地東部川東褶皺帶，構(gòu)造上位于萬(wàn)縣復(fù)向斜的南部與方斗山背斜帶西側(cè)的交匯處，邊界主要受北東向大耳山西斷層、石門斷層、吊水巖斷層、天臺(tái)場(chǎng)斷層等斷層及北西向?yàn)踅瓟鄬涌刂芠18-19]，主體構(gòu)造特征總體表現(xiàn)為南寬北窄、中部寬緩的特點(diǎn)，北東走向，主體為平緩的箱狀背斜，背斜兩側(cè)發(fā)育多條逆斷層（圖3）。從應(yīng)力特征來(lái)看，位于正向構(gòu)造焦石壩斷背斜的構(gòu)造高部位以張應(yīng)力為主[18-19]。筆者研究的地震數(shù)據(jù)面積為464 km2，目的層為五峰組—龍馬溪組的一套優(yōu)質(zhì)泥頁(yè)巖（五峰組與龍馬溪組巖性相同，且五峰組地層厚度小于10 m[1-2]，地震剖面上兩套地層無(wú)法區(qū)分）。

頁(yè)巖氣儲(chǔ)層有多種孔隙類型，但都是以微米—納米級(jí)孔隙為主。Yang等[20-21]針對(duì)泥頁(yè)巖特低孔、特低滲的特點(diǎn)，利用流體示蹤實(shí)驗(yàn)分析了泥頁(yè)巖地層的孔隙結(jié)構(gòu)特征，認(rèn)為：①頁(yè)巖內(nèi)部連通孔隙非常少，裂縫是頁(yè)巖氣擴(kuò)散的主要通道，裂縫是影響頁(yè)巖氣保存的關(guān)鍵因素；②大量的不連通微米—納米孔隙存在于泥頁(yè)巖烴源巖中，這些不連通孔隙內(nèi)氣體的不斷聚集可能是頁(yè)巖氣地層形成超壓的重要成因。

圖3 焦石壩構(gòu)造五峰組底界構(gòu)造圖

2.2 裂縫發(fā)育特征

利用各向異性來(lái)預(yù)測(cè)高角度裂縫，就是利用不同方位角的地震屬性差異來(lái)擬合方位各向異性橢圓，并進(jìn)一步預(yù)測(cè)裂縫的發(fā)育程度和方向。賀振華等[17]認(rèn)為裂縫特別是高角度裂縫的地震波動(dòng)力學(xué)屬性（頻率、衰減等）比其運(yùn)動(dòng)學(xué)屬性（速度、旅行時(shí)間）更為敏感。此外，在動(dòng)力學(xué)屬性中，衰減相關(guān)屬性與裂縫的關(guān)系更為密切[13]。筆者以收集到的成像測(cè)井?dāng)?shù)據(jù)為基礎(chǔ)進(jìn)行屬性優(yōu)選，最終確定使用85%衰減頻率屬性。該屬性是指在有效頻帶范圍內(nèi)，能量衰減到85%時(shí)對(duì)應(yīng)的頻率值（fm1），計(jì)算公式為：

式中fh和f1分別表示有效頻帶的最高頻率和最低頻率，Hz；p(f)表示頻譜函數(shù)。

利用頻率衰減屬性進(jìn)行各向異性橢圓的擬合得到研究區(qū)各向異性強(qiáng)度的數(shù)據(jù)體，可用于描述高角度裂縫的分布特征。首先從剖面上驗(yàn)證預(yù)測(cè)結(jié)果的準(zhǔn)確性，將預(yù)測(cè)結(jié)果的剖面圖在井點(diǎn)處放大，并與對(duì)應(yīng)深度段的成像測(cè)井進(jìn)行對(duì)比（圖4-a～d）。從對(duì)比情況來(lái)看，JY4井的無(wú)裂縫段在預(yù)測(cè)結(jié)果上表現(xiàn)為局部低值，JY1井、JY2和JY3井的裂縫段均對(duì)應(yīng)局部高值區(qū)，預(yù)測(cè)結(jié)果與成像測(cè)井的裂縫發(fā)育情況較為一致，表明預(yù)測(cè)結(jié)果的準(zhǔn)確性。

圖4 預(yù)測(cè)各向異性強(qiáng)度與單井成像測(cè)井對(duì)比圖

對(duì)得到的高角度裂縫數(shù)據(jù)體進(jìn)行平面成圖，得到焦石壩區(qū)塊五峰組—龍馬溪組高角度裂縫的平面分布特征（圖5）。圖5中平面圖的色標(biāo)表示各向異性強(qiáng)度，用于表征高角度縫的發(fā)育情況，暖色調(diào)代表高角度縫發(fā)育區(qū)；玫瑰圖的暖色調(diào)表示預(yù)測(cè)裂縫的優(yōu)勢(shì)方向。從圖5可以看出，裂縫方向與最大水平主應(yīng)力方向一致，以近東西向?yàn)橹?，北東東—南西西向?yàn)檩o；高角度縫主要發(fā)育在斷層附近，箱狀背斜區(qū)高角度裂縫不發(fā)育。此外，基于4口井的成像測(cè)井資料，對(duì)預(yù)測(cè)的裂縫方向進(jìn)行進(jìn)一步對(duì)比分析發(fā)現(xiàn)，單井預(yù)測(cè)裂縫方向與成像測(cè)井的最大水平主應(yīng)力方向較為一致，預(yù)測(cè)的裂縫方向準(zhǔn)確（表1）。

圖5 焦石壩構(gòu)造五峰組—龍馬溪組各向異性強(qiáng)度平面圖

2.3 地層壓力系數(shù)分布特征

筆者收集了8口的實(shí)測(cè)地層壓力數(shù)據(jù)（表2），以實(shí)測(cè)地層壓力數(shù)據(jù)為基礎(chǔ)，分別利用傳統(tǒng)的Fillippone公式法以及波阻抗反演法來(lái)預(yù)測(cè)研究區(qū)龍一1亞段的地層壓力系數(shù)分布特征（圖6-a、b），并基于實(shí)測(cè)地層壓力數(shù)據(jù)，進(jìn)行兩種方法預(yù)測(cè)結(jié)果的誤差對(duì)比（表2），驗(yàn)證兩種方法預(yù)測(cè)結(jié)果的準(zhǔn)確性。

結(jié)合圖6和表2的數(shù)據(jù)分析認(rèn)為：①?gòu)牡貙訅毫ο禂?shù)分布圖的分辨率來(lái)看，F(xiàn)illippone公式法預(yù)測(cè)的超壓發(fā)育帶呈大片連續(xù)分布，分辨率較低，波阻抗反演法的預(yù)測(cè)結(jié)果分辨率更高；②從地層壓力系數(shù)的分布特征來(lái)看，兩種方法預(yù)測(cè)的地層壓力系數(shù)整體趨勢(shì)較為一致，預(yù)測(cè)的超壓均主要分布在研究區(qū)的中部和東北部地區(qū)；③從預(yù)測(cè)精度的誤差分析來(lái)看，F(xiàn)illippone公式法預(yù)測(cè)的結(jié)果誤差較大；波阻抗反演法預(yù)測(cè)的平均相對(duì)誤差絕對(duì)值為12.86%，單井最大相對(duì)誤差絕對(duì)值為18.18%，預(yù)測(cè)精度相對(duì)更高。

2.4 保存系數(shù)

采用波阻抗反演法預(yù)測(cè)的地層壓力系數(shù)與開(kāi)發(fā)井的無(wú)阻流量關(guān)系如圖7所示。利用波阻抗反演法得出的地層壓力系數(shù)與單井無(wú)阻流量線性相關(guān)關(guān)系總體較好，但個(gè)別井、特別是低產(chǎn)井（JY3-3HF、JY7-2HF及JY1-3HF）的地層壓力系數(shù)與單井無(wú)阻流量一致性并不好，這說(shuō)明頁(yè)巖氣開(kāi)發(fā)井產(chǎn)量除了受地層壓力的影響，還與其他因素有關(guān)。

表1 疊前方位各向異性預(yù)測(cè)裂縫方向與單井裂縫方向?qū)Ρ缺?/p>

表2 兩種方法的地層壓力系數(shù)預(yù)測(cè)結(jié)果分析表

圖6 兩種方法預(yù)測(cè)的焦石壩構(gòu)造龍一1亞段地層壓力系數(shù)平面圖

圖7 焦石壩區(qū)塊部分頁(yè)巖氣井無(wú)阻流量與波阻抗反演法地層壓力系數(shù)對(duì)比圖

考慮到裂縫發(fā)育可能會(huì)對(duì)保存條件產(chǎn)生影響，因此，結(jié)合井上無(wú)阻流量數(shù)據(jù)，將高角度裂縫預(yù)測(cè)結(jié)果與地層壓力系數(shù)的預(yù)測(cè)結(jié)果進(jìn)行綜合對(duì)比分析，并通過(guò)各向異性強(qiáng)度與地層壓力系數(shù)構(gòu)建新的參數(shù)——保存系數(shù)（Preservation Indicator）來(lái)進(jìn)行保存條件與頁(yè)巖氣富集高產(chǎn)的關(guān)系分析，即

式中PIimp表示保存系數(shù)；Fden表示各向異性強(qiáng)度；pc-imp表示波阻抗反演法預(yù)測(cè)的地層壓力系數(shù)。

將前述研究中得到的各向異性強(qiáng)度數(shù)據(jù)體以及地層壓力系數(shù)數(shù)據(jù)體代入到式（5），即可得到保存系數(shù)的分布特征（圖8），將計(jì)算所得的保存系數(shù)與氣井無(wú)阻流量進(jìn)行進(jìn)一步對(duì)比（圖9）可以看出，保存系數(shù)與頁(yè)巖氣無(wú)阻流量對(duì)應(yīng)性較好，并且彌補(bǔ)了低產(chǎn)井與無(wú)阻流量之間的不吻合關(guān)系，從而驗(yàn)證了新的頁(yè)巖氣保存條件預(yù)測(cè)參數(shù)的有效性。

圖8 焦石壩構(gòu)造龍一1亞段保存系數(shù)平面圖

圖9 焦石壩區(qū)塊部分頁(yè)巖氣井保存系數(shù)與無(wú)阻流量對(duì)比圖

3 結(jié)論

1）頁(yè)巖氣藏的保存條件與產(chǎn)量的關(guān)系密切，綜合裂縫發(fā)育強(qiáng)度和地層壓力系數(shù)提出的新參數(shù)——保存系數(shù)，能夠較好的評(píng)價(jià)頁(yè)巖氣的保存條件。

2）疊前方位各向異性裂縫預(yù)測(cè)方法和疊前縱波阻抗反演方法兩種地球物理技術(shù)方法對(duì)疊前地震數(shù)據(jù)有一定要求，需要疊前數(shù)據(jù)有足夠?qū)挼姆轿唤羌拜^大的偏移距范圍，以保證預(yù)測(cè)結(jié)果的可靠性。

3）保存條件評(píng)價(jià)對(duì)于頁(yè)巖氣勘探開(kāi)發(fā)具有重要的參考意義，應(yīng)進(jìn)一步將地質(zhì)與地球物理相結(jié)合，形成更為完善的頁(yè)巖氣保存條件評(píng)價(jià)體系，從而推動(dòng)我國(guó)頁(yè)巖氣勘探開(kāi)發(fā)的進(jìn)程。