劉仕友 李洋森 沈利霞 閆安菊

(中海石油(中國(guó))有限公司湛江分公司 廣東湛江 524000)

鶯歌海盆地的勘探儲(chǔ)量發(fā)現(xiàn)集中于盆地中央凹陷帶，尤其在東方區(qū)中深層獲得了重大突破，先后發(fā)現(xiàn)了東方13-1、東方13-2兩個(gè)整裝大氣田，徹底消除了中深層高溫高壓條件下天然氣不能大規(guī)模成藏的顧慮[1-2]。而樂(lè)東區(qū)中深層與東方區(qū)具有類似的成藏條件，主要存在中央峽谷-水道和三角洲前緣-海底扇兩大儲(chǔ)集體系[3]，但其后在同一水道上部署的LD1X-1-1、LD1X-1-2等多口深井效果并不理想，且產(chǎn)能差異極大，目的層段砂巖整體表現(xiàn)為低孔—特低滲特征[4]，優(yōu)質(zhì)甜點(diǎn)儲(chǔ)層預(yù)測(cè)是當(dāng)前樂(lè)東區(qū)高溫高壓中深層勘探開(kāi)發(fā)面臨的主要難題。

前人針對(duì)低孔低滲甜點(diǎn)地球物理評(píng)價(jià)技術(shù)作了大量研究。在儲(chǔ)層預(yù)測(cè)方面，開(kāi)展了疊前反演方法研究，挖掘可以有效識(shí)別儲(chǔ)層物性的橫波信息，形成了以相控約束疊前同時(shí)反演[5]、擴(kuò)展彈性阻抗反演[6]以及基于馬爾可夫鏈地質(zhì)統(tǒng)計(jì)學(xué)反演[7]等為代表的技術(shù)方法。但對(duì)于儲(chǔ)層孔隙度預(yù)測(cè)，常通過(guò)建立波阻抗與孔隙度的線性關(guān)系進(jìn)行轉(zhuǎn)換，很難得到反映儲(chǔ)層縱橫向空間物性變化的準(zhǔn)確巖石物理模型，導(dǎo)致由地震波阻抗轉(zhuǎn)換得到的孔隙度預(yù)測(cè)精度不高，甜點(diǎn)儲(chǔ)層的預(yù)測(cè)精度也較低。在AVO烴類檢測(cè)方面，Ostrander[8]發(fā)現(xiàn)了含氣砂巖反射振幅隨偏移距增加而增大，含水砂巖反射振幅隨偏移距增加而減小的現(xiàn)象，首先提出利用反射系數(shù)隨入射角的變化識(shí)別亮點(diǎn)型含氣砂巖。但受區(qū)域地質(zhì)條件、埋深等因素限制，含氣儲(chǔ)層在深層低孔低滲環(huán)境中往往會(huì)表現(xiàn)為暗點(diǎn)特征，而并非亮點(diǎn)特征，直接利用AVO分析識(shí)別儲(chǔ)層的含氣性往往存在一定的不確定性，為此針對(duì)暗點(diǎn)型油氣藏需要建立更為可靠的流體檢測(cè)方法，有效降低鉆前AVO烴類檢測(cè)多解性的風(fēng)險(xiǎn)[9]。而暗點(diǎn)在地震剖面上以弱振幅形式隱含在干擾背景中，在識(shí)別技術(shù)上存在一定難度。這也是鶯歌海盆地樂(lè)東區(qū)中深層未獲重大突破的重要原因之一。

目前樂(lè)東區(qū)高溫高壓中深層黃流組和梅山組是天然氣增儲(chǔ)上產(chǎn)的主要層段，其優(yōu)質(zhì)儲(chǔ)層預(yù)測(cè)一直是本區(qū)勘探探索研究的關(guān)鍵，尋找有利潛在目標(biāo)是該區(qū)勘探開(kāi)發(fā)的首要目標(biāo)[10]。本文從區(qū)域巖石物理特征出發(fā)，揭示了樂(lè)東區(qū)暗點(diǎn)型油氣藏的形成原因，明確了暗點(diǎn)型油氣藏AVO的影響因素，并創(chuàng)新建立了II類AVO地球物理評(píng)價(jià)技術(shù)，在實(shí)際應(yīng)用中取得了良好效果。

1 暗點(diǎn)型氣藏形成原因

鶯歌海盆地樂(lè)東區(qū)中深層為高溫、異常高壓地層，具有獨(dú)特的地震速度特征(圖1)。該區(qū)已鉆井速度、密度測(cè)井曲線統(tǒng)計(jì)表明，鶯歌海組就開(kāi)始出現(xiàn)欠壓實(shí)低速泥巖，并且不同的沉降速率存在不同的欠壓實(shí)趨勢(shì)線[10-11]，黃流組一段低速泥巖中砂質(zhì)含量的增加導(dǎo)致低速泥巖速度向正常壓實(shí)趨勢(shì)線靠近，黃流組二段低速泥巖速度隨著砂質(zhì)含量的減少而又迅速降低，低速泥巖速度多為2 500～3 500 m/s，凹陷內(nèi)部低至1 900 m/s。而樂(lè)東區(qū)水道主力目的層位于黃流組二段，砂巖速度分布較為穩(wěn)定,為4 000～4 500 m/s，由于砂巖相對(duì)泥巖表現(xiàn)為低密度特征，因此樂(lè)東區(qū)深層黃流組二段含氣砂巖多表現(xiàn)為阻抗疊置或高阻抗特征，具備形成暗點(diǎn)氣藏的基本巖石物理?xiàng)l件。

圖1 鶯歌海盆地樂(lè)東區(qū)已鉆井巖石物理特征分析Fig.1 Characteristics analysis of rock physics for drilled wells in Ledong area of Yinggehai basin

2 暗點(diǎn)型氣藏AVO影響因素

由于樂(lè)東區(qū)低速泥巖速度變化范圍大，在利用AVO技術(shù)對(duì)暗點(diǎn)進(jìn)行分析時(shí)，必須明確不同泥巖速度對(duì)AVO類型的影響。選擇樂(lè)東區(qū)已鉆井進(jìn)行不同孔隙度、不同泥巖背景條件下含氣砂巖的AVO正演模擬，具體參數(shù)如表1所示。從AVO正演模擬結(jié)果上看(圖2)，不同蓋層泥巖速度對(duì)AVO特征影響較大，當(dāng)上覆高速泥巖或正常泥巖時(shí)，含氣儲(chǔ)層若表現(xiàn)為II類AVO時(shí)，極性反轉(zhuǎn)時(shí)孔隙度約為10%左右；而當(dāng)上覆欠壓實(shí)低速泥巖時(shí)，含氣儲(chǔ)層若表現(xiàn)為II類AVO時(shí)，極性反轉(zhuǎn)時(shí)孔隙度約為14%左右。

表1 不同泥巖背景條件下AVO正演模擬參數(shù)Table 1 AVO modeling parameters under different mudstone background conditions

圖2 不同泥巖背景條件下AVO正演模擬結(jié)果Fig.2 Simulation results of AVO modeling under different mudstone background conditions

為了進(jìn)一步定量分析AVO變化規(guī)律，建立了樂(lè)東區(qū)不同泥巖速度條件下的含氣砂巖AVO-孔隙度量版(圖3)。可以看出，不同上覆泥巖速度下，同一類AVO特征代表的儲(chǔ)層孔隙度不同，上覆泥巖速度越低，同一類AVO代表的儲(chǔ)層孔隙度越高；當(dāng)上覆泥巖速度不變時(shí)，隨著孔隙度的增加，含氣砂巖AVO類型逐漸從I類向II類變化。因此，相同AVO類型孔隙度大小受上覆泥巖背景的影響，通過(guò)分析高溫高壓含氣儲(chǔ)層上覆泥巖速度的大小可以較為準(zhǔn)確地應(yīng)用AVO類型判斷含氣砂巖孔隙度的相對(duì)大小。所以，在低速泥巖背景下，尋找II類AVO特征的潛在目標(biāo)可以幫助該區(qū)實(shí)現(xiàn)勘探突破。

圖3 樂(lè)東區(qū)AVO-孔隙度量版Fig.3 AVO-porosity template in Ledong area

3 暗點(diǎn)型氣藏預(yù)測(cè)技術(shù)

開(kāi)展樂(lè)東區(qū)暗點(diǎn)型氣藏預(yù)測(cè)具有重要的現(xiàn)實(shí)意義，其難點(diǎn)主要體現(xiàn)在2個(gè)方面，一是在全疊加剖面上多為弱振幅、弱連續(xù)的反射特征，無(wú)明顯的地質(zhì)體外觀形態(tài)，難以識(shí)別和精細(xì)落實(shí)；二是由于截距較小，用指示亮點(diǎn)的AVO乘積顯然已不適應(yīng)弱振幅II類AVO異常分析。因此，需要提出不同的方法來(lái)刻畫II類AVO暗點(diǎn)響應(yīng)特征。

3.1 技術(shù)原理

Rutherford和Williams[12]將暗點(diǎn)歸類為聲阻抗差趨于零的砂體，其AVO類型細(xì)分為2種類型：一種存在相位反轉(zhuǎn)(IIa類AVO)，即對(duì)于較小的入射角，其反射波振幅是減小的，隨著角度的增大，相位發(fā)生反轉(zhuǎn)，反射波振幅相應(yīng)增加；另一種沒(méi)有相位反轉(zhuǎn)(IIb類AVO)，具有較小的負(fù)法向入射反射系數(shù)，其振幅隨著入射角的增加而增加。

首先，對(duì)于存在相位反轉(zhuǎn)的IIa類AVO，為了有效刻畫樂(lè)東區(qū)甜點(diǎn)儲(chǔ)層的II類AVO特征，本文采用Fstack項(xiàng)來(lái)描述II類AVO屬性(圖4)：利用遠(yuǎn)近炮檢距的振幅差異來(lái)描述，用遠(yuǎn)炮檢距疊加振幅減去近炮檢距疊加振幅，在I類、III類、IV類AVO疊加響應(yīng)被減弱的同時(shí)，II類AVO含氣響應(yīng)則被加強(qiáng)，即

(1)

式(1)中：an和af分別為近道、遠(yuǎn)道振幅的疊加；θn、θf(wàn)和θmax分別為限定的近角孔徑、遠(yuǎn)角孔徑以及最大有效角度孔徑，其界限值主要取決于該區(qū)儲(chǔ)層響應(yīng)特征；c1為由響應(yīng)類別所確定的常數(shù)(當(dāng)為IIa類AVO含氣砂巖，則c1=1;當(dāng)為IIb類AVO含氣砂巖，由于近道疊加振幅較低，常?？珊雎?，則c1=0)。

圖4 II類AVO曲線示意圖Fig.4 Diagram of II class AVO curve

其次，對(duì)于沒(méi)有相位反轉(zhuǎn)的IIb類AVO含氣砂巖，F(xiàn)stack屬性雖然也具有指示意義，但較第IIa類AVO含氣砂巖的Fstack屬性異常稍弱。由于截距接近零值，傳統(tǒng)的AVO處理方法(如P×G屬性)顯然得不到油氣異常顯示，即使出現(xiàn)所謂的第II類異常也不是油氣的真實(shí)反映。在P-G反射坐標(biāo)系中，亮點(diǎn)型低阻含氣砂巖P×G屬性為強(qiáng)正值,位于第三象限內(nèi)；亮點(diǎn)型高阻含氣砂巖P×G屬性為強(qiáng)負(fù)值，位于第一象限內(nèi)；而暗點(diǎn)則位于第二、三象限內(nèi)，即有相位反轉(zhuǎn)和沒(méi)有相位反轉(zhuǎn)的II類AVO暗點(diǎn)型含氣砂巖都基本在G軸附近。為了能夠用P×G屬性刻畫II類AVO暗點(diǎn)型油氣藏，可以利用坐標(biāo)旋轉(zhuǎn)(圖5)將II類AVO的暗點(diǎn)型油氣藏旋轉(zhuǎn)至第三象限內(nèi)，即由坐標(biāo)系統(tǒng)P-G轉(zhuǎn)換至坐標(biāo)系統(tǒng)P′-G′，這時(shí)就可以借助常規(guī)亮點(diǎn)型P′×G′屬性來(lái)解釋暗點(diǎn)型油氣藏。

圖5 AVO分析坐標(biāo)旋轉(zhuǎn)示意圖Fig.5 Diagram of coordinate rotation in AVO analysis

根據(jù)Zoeppritz方程簡(jiǎn)化出截距和梯度項(xiàng)的反射系數(shù)公式，即Shuey公式為

Rpp=P+Gsin2θ

(2)

式(2)中：P為截距；G為梯度；θ為入射角。

坐標(biāo)旋轉(zhuǎn)后變換公式可以表示為

P′=Pcosα+Gsinα

(3)

G′=Gcosα-Psinα

(4)

式(3)、(4)中：α為旋轉(zhuǎn)角度。

利用樂(lè)東區(qū)已鉆井的含氣砂巖AVO關(guān)系可求得該旋轉(zhuǎn)角度，當(dāng)旋轉(zhuǎn)角度為35°～40°時(shí)，可將坐標(biāo)系統(tǒng)P-G系統(tǒng)下II類AVO的響應(yīng)特征旋轉(zhuǎn)至坐標(biāo)系統(tǒng)P′-G′下的III類AVO響應(yīng)特征，利用旋轉(zhuǎn)后P′×G′屬性可以有效評(píng)價(jià)暗點(diǎn)型油氣藏是否真實(shí)可靠，有效剔除假的Fstack屬性異常。

為進(jìn)一步驗(yàn)證Fstack屬性、旋轉(zhuǎn)P′×G′屬性分析的可靠性，分別設(shè)計(jì)I、IIa、IIb、III、IV類AVO含氣砂巖模型(表2)，提取不同屬性進(jìn)行對(duì)比分析(圖6)。在旋轉(zhuǎn)前P×G屬性可以較好地識(shí)別I、III、IV類AVO含氣砂巖，對(duì)IIa、IIb類AVO含氣砂巖，由于截距值較小，其P×G屬性也基本接近于零，對(duì)暗點(diǎn)型氣藏識(shí)別基本無(wú)效；而旋轉(zhuǎn)后的P′×G′屬性、Fstack屬性均對(duì)IIa、IIb類AVO含氣砂巖較為敏感，可以明顯地反映出暗點(diǎn)型II類AVO含氣砂巖的特征，為暗點(diǎn)型氣藏評(píng)價(jià)提供新的技術(shù)評(píng)價(jià)途徑。

表2 不同AVO類型模型參數(shù)及AVO屬性Table 2 Model parameters and attributes of different AVO types

圖6 不同AVO類型正演模擬及其屬性對(duì)比Fig.6 Forward simulation of different AVO types and comparison of their attributes

3.2 應(yīng)用效果

在實(shí)際勘探應(yīng)用過(guò)程中，往往需要借助屬性切片技術(shù)通過(guò)對(duì)部分疊加體提取Fstack屬性沿層切片進(jìn)行分析，在樂(lè)東區(qū)黃流組二段發(fā)現(xiàn)了一套具有Fstack屬性異常的暗點(diǎn)型潛在有利目標(biāo)(圖7)。過(guò)該屬性異常剖面發(fā)現(xiàn)，該砂體具有與NW—SE向軸向海底扇水道體系相似特征(圖8a)，在目的層段砂體頂面表現(xiàn)為明顯的II類AVO特征，在遠(yuǎn)偏移距出現(xiàn)明顯的極性反轉(zhuǎn)現(xiàn)象(圖8b)，進(jìn)一步說(shuō)明衍生Fstack屬性能夠指示II類AVO屬性異常的分布。另外,疊前反演結(jié)果表明該暗點(diǎn)型目標(biāo)具有明顯的低Vp/Vs異常(圖8c)，砂體尖滅清晰，且處于“泥包砂”的背景，具有較好的生、儲(chǔ)、蓋條件。在烴類檢測(cè)方面，通過(guò)截距-梯度坐標(biāo)旋轉(zhuǎn)后，P′×G′屬性由II類AVO暗點(diǎn)特征變換為III類AVO亮點(diǎn)特征(圖9)，可以較好地反映該套有利目標(biāo)砂體的含氣概率，為樂(lè)東區(qū)甜點(diǎn)儲(chǔ)層預(yù)測(cè)提供了有效技術(shù)支撐。2017年11月在該區(qū)實(shí)施LD1X-2-1井鉆探，在構(gòu)造邊緣部位鉆遇到氣層8 m，孔隙度為9.3%，新增探明儲(chǔ)量34.6×108m3。該套氣層的發(fā)現(xiàn)，進(jìn)一步證實(shí)了本文預(yù)測(cè)技術(shù)的有效性，展現(xiàn)了樂(lè)東區(qū)深層暗點(diǎn)型油氣藏的巨大勘探潛力。

圖7 鶯歌海盆地樂(lè)東區(qū)黃流組二段Fstack屬性切片F(xiàn)ig.7 Fstack attribute slice of the 2nd Member of Huangliu Formation in Ledong area of Yinggehai basin

圖8 過(guò)LD1X-2-1井純波地震剖面、Fstack屬性剖面及反演剖面對(duì)比(剖面位置見(jiàn)圖7)Fig.8 Comparison of raw seismic profile,Fstack attribute profile and inversion result of Well LD1X-2-1(see Fig.7 for location)

圖9 鶯歌海盆地樂(lè)東區(qū)黃流組二段坐標(biāo)旋轉(zhuǎn)前后P×G 屬性平面圖Fig.9 Plane diagram of P×G attribute before and after coordinate rotation of the 2nd Member of Huangliu Formation in Ledong area of Yinggehai baasin

4 結(jié)論

1) 巖石物理特征分析表明，鶯歌海盆地樂(lè)東區(qū)深層低速泥巖快速沉降，導(dǎo)致泥巖速度特別低，砂泥巖阻抗疊置或砂巖呈高阻抗特征，具備形成暗點(diǎn)型油氣藏的條件。

2) 不同地層條件下砂巖AVO-孔隙度量版顯示，受上覆泥巖背景速度影響，泥巖速度越低，II類AVO含氣砂巖代表的孔隙度越高，物性越好，因此含氣優(yōu)質(zhì)儲(chǔ)層應(yīng)表現(xiàn)為II類AVO異常。

3) 借助近遠(yuǎn)道疊加道集關(guān)系建立了II類AVO暗點(diǎn)型氣藏識(shí)別技術(shù)，通過(guò)提取Fstack屬性可以有效識(shí)別樂(lè)東區(qū)中深層有利潛在目標(biāo)。鉆井結(jié)果證實(shí)本文預(yù)測(cè)技術(shù)的有效性和可靠性，進(jìn)一步拓展了II類AVO流體指示因子在烴類檢測(cè)中的應(yīng)用。