佘鈺蔚，朱仕軍，朱鵬宇，黨洲濤，唐緒磊，楊維磊，胡桂林

(1.西南石油大學(xué) 地球科學(xué)與技術(shù)學(xué)院，四川 成都610500;2.中國(guó)石油川慶鉆探工程公司地球物理勘探公司，四川 成都610213;3.中國(guó)石油長(zhǎng)慶油田公司，陜西 西安710018)

0 引 言

AVO(amplitude versus offset)技術(shù)依據(jù)振幅隨炮間距的變化規(guī)律所包含的儲(chǔ)層巖性及其孔隙流體的性質(zhì)的信息，能直接進(jìn)行油氣預(yù)測(cè)［1］。1984年Ostander［2］提出AVO 技術(shù)后，作為一種重要的流體檢測(cè)工具，得到了迅速的發(fā)展和應(yīng)用［3-4］。其中，Zoeppritz 方程［5］是傳統(tǒng)AVO 的理論基礎(chǔ)，Shuey 兩項(xiàng)近似式［6］的出現(xiàn)，促使AVO 技術(shù)有革命性的發(fā)展。在近角度情況下截距和梯度能較好的刻畫AVO 特征，是巖石骨架、孔隙度及孔隙流體的綜合響應(yīng)，其交會(huì)圖分析技術(shù)能直觀簡(jiǎn)單的刻畫含氣性，因此在儲(chǔ)層預(yù)測(cè)中被廣泛使用［7-8］。

中國(guó)油氣田多發(fā)育薄層或薄互層儲(chǔ)層，隨著油氣勘探和開發(fā)的不斷深入，砂泥巖薄互層儲(chǔ)層無疑會(huì)成為勘探重點(diǎn)。近年來，Bakke［9］，Yinbin Liu［10］，Jingye Li［11］，姚 陳［12］、周 麗［13］，Wenyong Pan［14］等學(xué)者在單一薄層AVO 響應(yīng)特征方面取得了許多重要成果，這些成果為地層中單一薄層流體檢測(cè)定性解釋奠定了良好的基礎(chǔ)，但作為地層中沉積更為普遍更具工業(yè)開采價(jià)值的地質(zhì)結(jié)構(gòu)——砂泥巖薄互層，因其組合關(guān)系多，內(nèi)部結(jié)構(gòu)多樣，AVO 響應(yīng)復(fù)雜，導(dǎo)致相關(guān)的研究相對(duì)較少［15］，因此對(duì)薄互層的AVO 響應(yīng)特征進(jìn)行研究具有重要的意義。

薄互層的頂?shù)捉缑媸菦Q定反射特征的主要因素，利用已知的等厚薄互層巖性、油、氣和水的模型，分別在聲學(xué)介質(zhì)和彈性介質(zhì)條件下進(jìn)行數(shù)值模擬，優(yōu)選出更能精確刻畫厚薄互層頂界面AVO響應(yīng)特征的數(shù)值模擬方法。從而采用該波動(dòng)方程進(jìn)一步研究薄層數(shù)和單層厚度對(duì)薄互層截距、梯度及其交會(huì)屬性的影響規(guī)律，有助于從實(shí)際疊前地震數(shù)據(jù)中識(shí)別含氣薄互層有利區(qū)，提高儲(chǔ)層預(yù)測(cè)的精度。

1 正演模擬

設(shè)計(jì)如圖1 所示的地質(zhì)模型，該模型由均勻泥巖背景中包含不同互層數(shù)的砂泥巖等厚薄互層地質(zhì)體組成，砂泥巖的互層數(shù)分別設(shè)定為:1，2，4，6，8，10，借以考察不同互層數(shù)對(duì)薄互層的AVO 響應(yīng)的影響。其中砂巖骨架Vp=6 000 m/s，Vs=3 000 m/s，ρ=2.65 g/cm3. 泥巖阻抗為儲(chǔ)層骨架阻抗的80%.地震子波采用零相位雷克子波，主頻為30 Hz.采用交錯(cuò)有限差分法正演模擬獲取單炮記錄;采取中間放炮、兩邊接收的的采集方式，地震道采樣點(diǎn)數(shù)為4 000，采樣間隔0.001 s，道間距20 m，共400 道，單層厚度從1 以2 m 為間隔變化到35 m.對(duì)孔隙度12%，含氣飽和度80%的儲(chǔ)層模擬了這6 種薄互層模型得到108 套地震數(shù)據(jù)體。

2 聲學(xué)介質(zhì)與彈性介質(zhì)

Yinbin Liu 和Schmitt［7］已研究單一薄層在聲學(xué)介質(zhì)下的AVO 響應(yīng)。Wenyong Pan 和Kristopher A［11］在前者的基礎(chǔ)上，進(jìn)行單一薄層在彈性介質(zhì)下的AVO/AVF 特征研究，表明聲學(xué)介質(zhì)不適合單一薄層AVO 研究。類似的，針對(duì)薄互層的AVO 響應(yīng)，利用Grassmann 方程進(jìn)行流體替代，分別在聲學(xué)介質(zhì)和彈性介質(zhì)下，設(shè)計(jì)等厚四層5 m 地質(zhì)模型(參數(shù)表(1))進(jìn)行數(shù)值模擬。在0°入射的情況下，聲學(xué)介質(zhì)和彈性介質(zhì)下砂巖頂界面的振幅值基本相同。然而，當(dāng)入射角大于20°兩種介質(zhì)情況下AVO 差異越來越明顯。含水砂巖薄互層的情況下，聲學(xué)介質(zhì)中振幅隨入射角增大出現(xiàn)了砂巖底界面假亮點(diǎn)的現(xiàn)象，而彈性介質(zhì)中出現(xiàn)第二類AVO 響應(yīng)。將含水砂巖替換成含氣砂巖后，聲學(xué)介質(zhì)中在50°出現(xiàn)極性反轉(zhuǎn)的現(xiàn)象，而彈性介質(zhì)則無此現(xiàn)象。因此，在預(yù)測(cè)流體性質(zhì)時(shí)上述解釋陷阱的出現(xiàn)表明在進(jìn)行薄互層的AVO 分析時(shí)采用彈性介質(zhì)更加精確。

3 數(shù)據(jù)分析與認(rèn)識(shí)

3.1 薄互層AVO 特征與互層數(shù)之間的關(guān)系

從蓋層為儲(chǔ)層骨架阻抗的80%，儲(chǔ)層孔隙度為12%，含水飽和度為20%的等厚薄互層彈性波動(dòng)方程正演結(jié)果圖3 分析，薄互層仍然能產(chǎn)生清晰的AVO 異常，且具有一定的規(guī)律。因此能夠利用薄互層AVO 響應(yīng)特征進(jìn)行烴類檢測(cè)。為刻畫薄互層振幅隨偏移距變化的特征，采用Shuey 兩項(xiàng)式的思想，把sin2(θ)視為變量，對(duì)0° ～30°的薄互層AVO 響應(yīng)進(jìn)行線性擬合。從表2 知擬合相關(guān)系數(shù)均大于0.973，由此說明截距和梯度仍然可以描述近角度入射情況下薄互層AVO 的特征。

單層砂體厚度/m 一層砂 二層砂 四層砂 六層砂 八層砂 十層砂0.992 0.997 0.997 0.998 0.993 0.998 3 0.993 0.997 0.997 0.996 0.994 0.999 5 0.997 0.996 0.997 0.995 0.997 0.998 7 0.993 0.996 0.997 0.997 0.997 0.998 9 0.994 0.998 0.982 0.992 0.996 0.996 11 0.996 0.999 0.996 0.997 0.997 0.997 13 0.998 0.996 0.992 0.993 0.994 0.993 15 0.998 0.996 0.998 0.999 0.991 0.999 17 0.998 0.999 0.999 0.999 0.999 0.996 19 0.998 0.997 0.997 0.999 0.999 0.999 21 0.999 0.984 0.993 0.992 0.985 0.994 23 0.998 0.964 0.989 0.988 0.987 0.984 25 0.998 0.973 0.995 0.994 0.993 0.991 27 0.998 0.984 0.994 0.994 0.995 0.995 29 0.994 0.982 0.993 0.990 0.994 0.987 31 0.998 0.997 0.997 0.997 0.991 0.994 33 0.998 0.997 0.996 0.982 0.998 0.994 1 35 0.998 0.995 0.997 0.999 0.999 0.999

3.1.1 截距

單一薄層法線入射時(shí)AVO 的截距和梯度的絕對(duì)值與砂體厚度密切相關(guān)。當(dāng)砂體厚度小于1/64波長(zhǎng)時(shí)，沒有AVO 響應(yīng);小于1/4 波長(zhǎng)范圍內(nèi)，單一砂體法線入射時(shí)AVO 的截距和梯度的絕對(duì)值隨著砂體厚度的增加單調(diào)遞增(圖4)。

等厚薄互層頂界面AVO 的截距絕對(duì)值也與互層數(shù)密切相關(guān)。當(dāng)單層厚度小于1/64 波長(zhǎng)時(shí)，單一薄砂層沒有AVO 響應(yīng)，但單層厚度小于1/64 波長(zhǎng)的多個(gè)薄互層也可形成較強(qiáng)的AVO 響應(yīng)(圖3(a)、圖4)。當(dāng)單層厚度大于1/64 波長(zhǎng)且小于1/8 波長(zhǎng)時(shí)，單一砂體厚度下，截距絕對(duì)值隨互層數(shù)的增加，呈調(diào)諧效應(yīng)，先增加后減小，調(diào)諧振幅大于單一薄層時(shí)的振幅值(圖3(b)、圖4)。其中1/32 到1/27 波長(zhǎng)之間4 層薄互層處在調(diào)諧厚度;1/27 到1/10 波長(zhǎng)之間2 層薄互層處在調(diào)諧厚度。因此，薄互層隨著單層砂體厚度增加，出現(xiàn)調(diào)諧截距的層數(shù)在遞減。當(dāng)單層厚度大于1/8 波長(zhǎng)且小于等于1/4 波長(zhǎng)時(shí)，薄互層截距幾乎不受砂體層數(shù)的影響，且薄互層截距一般小于單一薄砂巖截距值(圖3(e)、圖4)。

3.1.2 梯度

等厚薄互層頂界面AVO 的梯度絕對(duì)值也與互層數(shù)密切相關(guān)。當(dāng)單層厚度小于1/64 波長(zhǎng)時(shí)，10層截距和梯度值最大(圖4，圖5)。當(dāng)單層厚度大于1/64 波長(zhǎng)且小于1/8 波長(zhǎng)時(shí)，固定單層砂體厚度時(shí)，薄互層梯度絕對(duì)值隨互層數(shù)的增加，呈調(diào)諧效應(yīng)，先增加后減小，調(diào)諧梯度大于單一薄層時(shí)的梯度值(圖3(b)、圖5)。其中1/64 到1/32 波長(zhǎng)8層薄互層產(chǎn)生調(diào)諧梯度;1/32 到1/16 波長(zhǎng)之間4層薄互層處在調(diào)諧厚度，出現(xiàn)梯度最大值;而1/16到1/8 波長(zhǎng)之間2 層薄互層產(chǎn)生調(diào)諧梯度。因此，薄互層隨著單層砂體厚度增加，出現(xiàn)調(diào)諧梯度的層數(shù)在遞減。當(dāng)單層厚度大于1/8 波長(zhǎng)且小于1/4 波長(zhǎng)時(shí)，梯度幾乎不受互層數(shù)的影響，且薄互層梯度小于單一薄砂巖梯度值(圖3(d)、圖5)。

薄互層梯度與截距變化大體趨勢(shì)相似，但又有區(qū)別。隨著薄互層單層砂體厚度增加，出現(xiàn)調(diào)諧梯度和梯度的層數(shù)均在遞減，且截距變小更快。單層砂巖厚度大于1/8 波長(zhǎng)后，梯度不受薄互層層數(shù)的影響，而截距出現(xiàn)在1/10 波長(zhǎng)處。因此，截距對(duì)薄互層單砂巖厚度更加敏感。與單砂體進(jìn)行對(duì)比，受薄互層單層厚度和層數(shù)變化的調(diào)諧影響，薄互層單層厚度小于1/8 波長(zhǎng)時(shí)，薄互層調(diào)諧截距和梯度均大于等于單砂巖的值，對(duì)AVO 效應(yīng)起到加強(qiáng)的作用，更易識(shí)別;但當(dāng)薄互層單層厚度大于1/8 波長(zhǎng)小于1/4 波長(zhǎng)時(shí)，AVO 效應(yīng)較單層砂巖相同厚度時(shí)減弱。

3.2 薄互層P-G 交會(huì)屬性特征

經(jīng)上述薄互層截距和梯度隨單層厚度變化特征分析知，可按照截距的變化規(guī)律，將小于1/4 波長(zhǎng)單層厚度分為5 個(gè)區(qū)間，即(0，1/64)(1/64，1/32)(1/32，1/16)(1/16，1/8)(1/8，1/4)，且每個(gè)區(qū)間內(nèi)薄互層截距和梯度隨互層數(shù)的變化規(guī)律基本一致。因此，可按照單層砂體厚度將小于1/4 波長(zhǎng)薄互層按區(qū)間劃分為5 種類型。尋找該5 種類型薄互層的規(guī)律，便能指導(dǎo)劃分出有利含氣區(qū)域。因此，分別進(jìn)行單層厚度1，5，11，21 和33 m 的薄互層進(jìn)行P-G 交會(huì)屬性研究。

將5 種含氣薄互層與含氣厚層和含水厚層進(jìn)行對(duì)比(圖6)，發(fā)現(xiàn)薄互層在P -G 交會(huì)圖上分布具有一定的規(guī)律。當(dāng)薄互層單砂體厚度落入(0，1/64)波長(zhǎng)，含氣薄互層P -G 分布在厚層的右上方黃色區(qū)域(圖6(a))。當(dāng)薄互層單砂體厚度為(1/64，1/32)波長(zhǎng)時(shí)，含氣薄互層P -G 分布在厚層的右上方黃色區(qū)域，非常接近含水厚層位置(圖6(b))。而單砂體厚度落入(1/32，1/16)波長(zhǎng)范圍內(nèi)，含氣薄互層P-G 分布在厚層的右下方黃色區(qū)域，接近含水厚層位置(圖6(c))。當(dāng)單砂體厚度為(1/16，1/8))波長(zhǎng)時(shí)，含氣薄互層P -G 分布以厚層為中心的黃色區(qū)域(圖6(d))，開始遠(yuǎn)離含水厚層區(qū)域。當(dāng)單砂體厚度落入(1/8，1/4)波長(zhǎng)時(shí)，含氣薄互層P-G 分布在厚層左上方黃色區(qū)域(圖6(e))，遠(yuǎn)離含水厚層。最后，當(dāng)薄互層單砂體厚度等于1/4 波長(zhǎng)時(shí)，含氣薄互層P-G 略小于單層含氣砂巖，且遠(yuǎn)離含氣厚層和含水厚層，出現(xiàn)典型的Ⅲ類AVO 現(xiàn)象，有利區(qū)分布在厚層左下方黃色區(qū)域(圖6(f))。由此說明，隨著薄互層單層厚度的增加，P-G 由右上區(qū)域變化到左下區(qū)域。

4 實(shí)際資料的分析

選擇蘇里格地區(qū)某氣田高產(chǎn)氣井13 ～51 井，進(jìn)行AVO 正演模擬研究。該井儲(chǔ)蓋條件、含氣飽和度和孔隙度與設(shè)計(jì)模型大致相當(dāng)，13 ～51 井具有2 層約21 m 的砂巖儲(chǔ)層。其正演模擬及實(shí)際井旁的地震偏移距道集和P -G 交會(huì)圖分別如圖7，圖8 所示?？梢钥闯?3 -51 井正演模擬及實(shí)際井旁道集的P-G 屬性均分布在左上區(qū)域，與設(shè)計(jì)模型單層厚度為21 m 時(shí)P -G 屬性(圖6(e))分布規(guī)律一致。

5 結(jié) 論

1)進(jìn)行薄互層的AVO 數(shù)值模擬分析時(shí)聲學(xué)介質(zhì)條件下會(huì)出現(xiàn)假亮點(diǎn)，在彈性介質(zhì)中進(jìn)行數(shù)值模擬更加合理;

2)薄互層單層厚度小于1/8 波長(zhǎng)時(shí)，P -G 均大于等于單層砂相同厚度的P-G 值，薄互層對(duì)頂面AVO 效應(yīng)起到加強(qiáng)的作用。當(dāng)薄互層單層厚度大于1/8 波長(zhǎng)小于1/4 波長(zhǎng)時(shí)，AVO 效應(yīng)較薄層情況減弱;

3)隨著薄互層單層厚度的增加，P -G 由右上區(qū)域變化到左下區(qū)域，有一定的分布規(guī)律。因此可以根據(jù)工區(qū)薄互層單層厚度和層數(shù)確定出薄互層發(fā)育的P-G 屬性的分布范圍，從而指導(dǎo)儲(chǔ)層預(yù)測(cè)，為勘探開發(fā)提供依據(jù);

4)實(shí)際井資料正演模擬及實(shí)際井旁道集的P-G 屬性均分布在左上區(qū)域，分布規(guī)律一致。

References

［1］ 孫鵬遠(yuǎn).AVO 技術(shù)新進(jìn)展［J］. 勘探地球物理進(jìn)展，2005，28(6):432 -438.SUN Peng-yuan.Advances in AVO techniques［J］.Progress in Exploration Geophysics，2005，28(6):432 -438.

［2］ Ostrander W.Plane-wave reflection coefficients for gas sands at nonormal angles of incidence［J］. Geophysics，1984，49(6):1 637 -1 648.

［3］ Rutherford Williams.Amplitudes-versus-offset in gas sands［J］.Geophysics，1989(26):138 -157.

［4］ 王寶江，李 美，于 強(qiáng)，等. 復(fù)雜氣藏AVO 屬性交會(huì)含氣性檢測(cè)應(yīng)用研究［J］.地球物理學(xué)進(jìn)展，2013，28(1):310 -319.WANG Bao-jiang，LI Mei，YU Qiang，et al. Predicting ans detecting of the gas sand by AVO attributes crossplot in complex reservoir［J］.Progress in Geoghys，2013，28(1):310 -319.

［5］ Zoeppritz K，Uber Erdbebenwellen.On the reflection and propagation of seismic waves［J］. Gottingen Nachricten der Konigl，1919(2):66 -84.

［6］ Shuey R T. A simplification of the Zoeppritz equations［J］.Geophysics，1985，50(4):539 -546.

［7］ 潘仁芳，陳思路，張利萍，等. 含氣砂巖AVO 正演的半定量分析［J］. 石油地球物理勘探，2013，48(1):103 -108.PAN Ren-fang，CHEN Si-lu，ZHANG Li-ping，et al.Semi-quantitative analysis of gas bearing sand by AVO forward modeling［J］.Oil Geophysical Prospecting，2013，48(1):103-108.

［8］ 張德英，陳全紅.疊前彈性波阻抗反演技術(shù)在赤道幾內(nèi)亞X 區(qū)塊儲(chǔ)層預(yù)測(cè)中的應(yīng)用［J］.西安科技大學(xué)學(xué)報(bào)，2013，33(2):173 -177.ZHANG De-ying，CHEN Quan-hong.Application of prestack elastic wave impedance inversion technique in reservoir prediction of X block in equatorial guinea［J］.Journal of Xi’An University of Science and Technology，2013，33(2):173 -177.

［9］ Bakke N E，Ursin B.Thin-bed AVO effects［J］.Geophysical Prospecting，1998，46(6):571 -587.

［10］LIU Y，Schmitt D R.Amplitude and AVO responses of a single thin bed［J］. Geophysics，2003，68(4):1 161 -1 168.

［11］LI Jing-ye，CHEN Xiao-hong.Seismic AVO modeling and analysis for thin inter-layer reservoir［C］//SEG Houston 2009 International Exposition and Annual Meeting，2009:316 -320.

［12］姚 陳，蔡明剛，王 贄.各向同性薄層反射理論地震圖［J］.地球物理學(xué)報(bào)，2010，53(1):164 -170.YAO Chen，CAI Ming-gang，WANG Zhi. Synthetic seismograms of reflection from isotropic thin layer［J］.Chinese Journal of Geophysics，2010，53(1):164 -170.

［13］周 麗，顧漢明，馬靈偉，等.基于波動(dòng)方程正演模擬分析薄砂層含不同流體的AVO 特征［J］. 地質(zhì)科技情報(bào)，2013，32(2):169 -173.ZHOU Li，GU Han-ming，MA Ling-wei，et al.Analysis of AVO characteristics on thin san reservoir with various fluid based seismic modeling［J］. Geological Science and Technology Information，2013，32(2):169 -173.

［14］ PAN Wen-yong，Kristopher A. AVO/AVF analysis of thin-bed in elastic media［J］.SEG Houston 2013 Annual Meeting，2013:373 -377.

［15］趙 偉，陳小宏，李景葉.薄互層調(diào)諧效應(yīng)對(duì)AVO 的影響［J］.石油物探，2006，45(6):570 -573.ZHAO Wei，CHEN Xiao-hong，LI Jing-ye. Analysis of impact of thin effect on AVO［J］.Geophysical Prospecting for Petroleum，2006，45(6):570 -573.