资源描述：
《AVO地震反演在油气预测中的应用.pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、第9卷第4期特种油气藏Vol19No142002年8月SpecialOilandGasReservoirsAug12002文章编号:100626535(2002)0420022203AVO地震反演在油气预测中的应用斌1211方 贝,柯善珍,陈家敏,刘 敬(11中油辽河油田分公司勘探开发研究院,辽宁 盘锦 124010;21中石化中原油田分公司信息中心,河南 濮阳 473000)摘要:利用AVO分析和地震反演相结合的三维多道全局优化地震反演技术和方法进行叠前反演。结果得到纵横波波阻抗、纵横波速度比、泊松
2、比等相关数据;结合岩性贡献概率量板,最终得到岩性数据体。将其用于岩性预测和地下含油气信息检测后效果良好。由此认为:本方法能处理各种道集数据,是一种新的流体预测和岩性预测方法。关键词:AVO反演;全局优化;角度叠加;截距梯度;低频模型;门槛值中图分类号:P631 文献标识码:A引 言AVO(AmplitudeVersusOffset)分析,是提取多种岩性参数的新技术,可用于寻找“亮点型”的气藏。目前,在利用不同岩性的地质模型来研究其AVO特性方面,前人已做了大量研究工作,而利用AVO分析和地震反演相结
3、合来研究储层的物理特性才刚刚开始。AVO分析和地震反演相结合可以提供一种新的流体预测和岩性预测的方法。AVO地震数据的物性特征反演可直接对岩性进行解释,预测地下岩性和流体的变化。它所用的反演方法是三维多道全局优化地震反演技术,其子波的频率和能量是空变的,可采用偏移距叠加、角道集、梯度和截距资料,使用叠加和偏移速度,实现对巨大数据体的准确的AVO反演。采用AVO地震资料的实际反演增强了资料的实用性,从而使获得的的物性剖图1 入射P波与其产生的反射波和透射波的关系面更加可靠。本文给出了AVO三维多道全局优
4、化Fig11RelationbetweenincidentPwaveandits地震反演技术的基本原理、软件实现的基本方法及backandtransmittedwaves处理的流程和研究实例。波系数)、ATS(透射横波系数);再假定入射波的能[1]1 基本原理量为1,且各波振幅大小满足物理条件:斯奈尔定律;位移的连续性;应力的连续性,那么,可以推如图1所示,假定有两层水平各向同性介质导出佐普里兹(Zoeppritz)方程。该方程能精确地反v1、ρ1、w1、σ1和v2、ρ2、w2、σ2,当纵波从介质1到
5、映各波振幅与入射角的关系,但由于所含变量较达介质2表面时将产生4种波,其振幅分别为ARP多,不便于分析求解,为此R1T·Shuey导出了一个(反射纵波系数)、ARS(反射横波系数)、ATP(透射纵[2]简化公式: 收稿日期:2001211219;改回日期:2002204209 作者简介:方斌贝(19652),男,工程师,1989年毕业于中国地质大学,现从事科技期刊编辑工作。电话:042727823579。斌第4期 方 贝等:AVO地震反演在油气预测中的应用232RP(θ)≈
6、I+Gsinθ(1)可以看出,通过综合应用纵波波阻抗和泊松比,可式中:RP(θ)相当于CDP道集内各道的振幅;I为以将砂泥岩区分开来,进而可以将含气砂岩和含水垂直入射的反射系数;G为振幅随入射角(炮检距)砂岩区分开来。的变化率。变化率G具有特殊的物理意义:(1)当正反射(I>0)时,如果G>0,则振幅随炮检距的增加而增大;如果G<0,振幅随炮检距的增加而减小。(2)当负反射(I<0)时,如果G>0,则振幅随炮检距的增加而衰减;如果G<0,则振幅随炮检距的增加而增大。2 实现方法用于AVO反演的数据可以
7、是截距梯度数据、角度叠加数据或偏移叠加数据。主要是基于Shuey的近似方程,采用伽德纳常量和从测井曲线数据估算的低频VP/VS模型:图2 纵波阻抗与泊松比的经验关系(图中画出了2RP(θ)≈I+Gsinθ(2)砂岩和泥岩的经验曲线,纯流体的一般变化趋势和约束带用箭头表示)对于截距梯度数据,第一步是计算纵横波反射Fig12EmpiricalrelationofP-waveimpedanceandPoisson’s系数剖面(方程3和4),然后估算纵横波的子波,ratio(empiricalcurveofs
8、andstoneandmudstoneisshown,normal主要用于反演。反演结果是纵横波波阻抗,可用来trendofpurefluidandconstrainingbeltareshownbyarrows)计算用于指示岩性的数据体,如泊松比、纵横波比、图3显示出截距资料(RP)、计算的横波资料拉梅常数等。(RS)及没有加低频信息的纵横波反演结果。其中RP=I(3)横波资料是用截距和梯度资料采用方程(4)计算出RS=[I/(1+α)-G]/[2(2
