一种基于岩相约束的AVO反演油气检测方法技术

本发明专利技术提供了一种基于岩相约束的AVO反演油气检测方法，属于流体检测领域。所述方法包括以下步骤：(1)弹性参数分布特征统计；(2)蒙特卡罗仿真随机模拟；(3)对叠前地震道集进行叠前道集预处理得到预处理后的叠前道集；(4)入射角道集形成；(5)先验信息约束的最优AVO反演实现。本发明专利技术方法比常规AVO方程具有更高的精度和更大的适用范围，更适用于大入射角及强反射界面情况下的流体性质检测，并提高了流体预测精度。

【技术实现步骤摘要】
—种基于岩相约束的AVO反演油气检测方法
本专利技术属于地震资料油气检测，具体涉及一种基于岩相约束的AVO(振幅随偏移距的变化)反演油气检测方法，可应用于复杂储层，特别是适用于强反射界面情况的流体检测应用。
技术介绍
自Ostrander (1984)提出AVO方法可以帮助描绘储层的流体含量以来，奠定了以AVO属性为基础的油气预测技术基础。地球物理学家分别从不同角度对Zoeppritz方程进行简化以建立不同形式的AVO方程，实现AVO属性反演及流体检测。Smi th和Gidlow(1987)组合Ak1-Richards方程使用泥岩线识别有碳氢化合物储层指示的地震数据的异常，形成了流体异常识别的截距和梯度交汇图法；Goodway等(1997)提出了一种流体异常识别的LMR法，能够获得拉梅参数和密度；Hedllin (2000)提出的孔隙模量方法涉及到P波和S波阻抗的组合；Russell等(2003)根据Biot-Gassmann理论提取了一种新的流体识别因子，基本上总结了前人的结果。 利用AVO属性分析及反演技术开展流体检测，国内学者也提出多种方法开展AVO方程近似研究及流体识别因子建立以提高流体检测精度。宁中华等人提出的基于地震资料的高灵敏度流体识别因子；李景叶等人提出的不受孔隙度变化影响的新流体识别因子；贺振华等人提出的扩展流体识别因子等，这些流体识别因子在流体识别特别是气藏识别方面取得了良好的效果。由于直接求解AVO属性赖以建立的Zc^ppritz方程极为复杂，并且难以给出清楚的物理概念，AVO近似方程大多是建立在一系列假设基础之上的，且与实际地质情况关联性较差，因此影响AVO属性进行油气检测的精度。随着油气勘探开发的不断深入，油气储集体越来越复杂，勘探重点由原来的构造油气藏向岩性油气藏、隐蔽油气藏转移。一方面，由于勘探目标逐渐变深、变小和变得复杂，勘探难度加大，对储层预测的精度要求也越来越高。另一方面，储层预测的最终目的除了判别储层的存在与否，预测储层发育程度的好与坏，同时也需要判定储层中所含流体的性质，流体识别已成为复杂储层预测中遇到的新的挑战和瓶颈问题。例如中石化西部塔河油田碳酸盐岩储层流体预测不同于碎屑岩储层的流体预测，其埋藏深、非均质强、储层与围岩阻抗差异大，造成流体难以区分、油气预测难度大。
技术实现思路
本专利技术的目的在于解决上述现有技术中存在的难题，提供一种基于岩相约束的AVO反演油气检测方法，解决中下奥陶统碳酸盐岩缝洞型储层与围岩阻抗差异大、非均质性强及流体检测难度大等问题。本专利技术是通过以下技术方案实现的:所述方法包括以下步骤:(I)弹性参数分布特征统计:依据目标区的岩芯、测井信息开展岩石物理分析，建立目标层的不同岩相组合弹性参数分布特征；(2)蒙特卡罗仿真随机模拟:针对步骤(1)得到的所述不同岩相组合弹性参数分布特征，基于蒙特卡罗仿真随机模拟，建立一个与求解有关的概率模型，并使随机抽样结果的解与该概率模型的统计量相一致，然后通过对随机抽样结果进行不同岩相正演模拟求得基函数f ( Θ );(3)对叠前地震道集进行叠前道集预处理得到预处理后的叠前道集；预处理包括道编辑、带通滤波、真振幅恢复、静校正、速度分析、剩余静校正、地表振幅一致性补偿、叠前反褶积及动校正；(4)入射角道集形成:基于步骤(3)得到的预处理后的叠前道集及叠加速度体或层速度体，实现偏移距道集到角度道集的转换，获得叠前角度道集；(5)先验信息约束的最优AVO反演实现:基于最优AVO方程反演目标函数如下:本文档来自技高网...

【技术保护点】
一种基于岩相约束的AVO反演油气检测方法，其特征在于：所述方法包括以下步骤：(1)弹性参数分布特征统计：依据目标区的岩芯、测井信息开展岩石物理分析，建立目标层的不同岩相组合弹性参数分布特征；(2)蒙特卡罗仿真随机模拟：针对步骤(1)得到的所述不同岩相组合弹性参数分布特征，基于蒙特卡罗仿真随机模拟，建立一个与求解有关的概率模型，并使随机抽样结果的解与该概率模型的统计量相一致，然后通过对随机抽样结果进行不同岩相正演模拟求得基函数f(θ)；(3)对叠前地震道集进行叠前道集预处理得到预处理后的叠前道集；(4)入射角道集形成：基于步骤(3)得到的预处理后的叠前道集及叠加速度体或层速度体，实现偏移距道集到角度道集的转换，获得叠前角度道集；(5)先验信息约束的最优AVO反演实现：基于最优AVO方程反演目标函数如下： | | R r ( θ ) - c r f ( θ ) | | 2 ⇒ min 其中，其中，Rr(θ)为叠前道集反射系数，f(θ)为基函数，cr为待反演属性；在基函数f(θ))的约束下，由叠前角度道集开展AVO属性反演，实现弹性参数提取及流体性质检测。...

【技术特征摘要】
1.一种基于岩相约束的AVO反演油气检测方法，其特征在于:所述方法包括以下步骤: (1)弹性参数分布特征统计:依据目标区的岩芯、测井信息开展岩石物理分析，建立目标层的不同岩相组合弹性参数分布特征； (2)蒙特卡罗仿真随机模拟:针对步骤(1)得到的所述不同岩相组合弹性参数分布特征，基于蒙特卡罗仿真随机模拟，建立一个与求解有关的概率模型，并使随机抽样结果的解与该概率模型的统计量相一致， 然后通过对随机抽样结果进行不同岩相正演模拟求得基函数 f(Θ ); (3)对叠前地震道集进行叠前道集预处理得到预处理后的叠前道集； (4)入射角道集形成:基于步骤(3)得到的预处理后的叠前道集及叠加速度体或层速度体，实现偏移距道集到角度道集的转换，获得叠前角度道集； (5)先验信息约束的最优AVO反演实现:基于最优AVO方程反演目标函数如下: 2.根据权利要求1所述的基于岩相约束的AVO反演油气检测方法，其特征在于:所述步骤(1)中的所述岩芯、测井信息包括声波、密度、中子、深侧向电阻率、浅侧向电阻率、自然伽玛和井径； 所述步骤(1)中的所述开展岩石物理分析包括以下步骤: (Al)，依据测井资料与不同矿物成份、流体物性响应关系，获得地层中骨架矿物成份的体积含量、孔隙度及含流体饱和度，具体如下: 假定有N-1条测井曲线，利用下面的线性超定方程组计算得到所述体积含量、孔隙度及含流体饱和度: 3.根据权利要求2所述的基于岩相约束的AVO反演油气检测方法，其特征在于: 所述步骤(1)中的所述建立目标层的不同岩相组合弹性参数分布特征是这样实现的:分别以纵波速度、横波速度及密度为横坐标，以相应样点数为纵坐标而形成的直方图和概率分布曲线。4.根据权利要求3所述的基于岩相约束的AVO反演油气检测方法，其特征在于: 所述步骤(2)中的所述针对步骤(1)得到的所述不同岩相组合弹性参数分布特征，基于蒙特卡罗仿真随机模拟，建立一个与求解有关的概率模型具体，使随机抽样结果的解与该概率模型的统计量相一致包括以下步骤: (BI)，对不同岩相的纵波速度、横波速度及密度进行分析； (B2)，依据不同岩相的纵波速度、横波速度及密度...

【专利技术属性】
技术研发人员：唐金良，王世星，周单，朱克，
申请(专利权)人：中国石油化工股份有限公司， 中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11