[发明专利]基于VTI介质的页岩储层AVO反演表征方法及系统

说明书

本发明提供一种基于VTI介质的页岩储层AVO反演表征方法及系统，属于页岩气勘探技术领域，确定地震数据的振幅缩放因子；构建VTI介质新属性参数的反射系数方程；基于测井数据，统计新属性参数的先验信息，结合反射系数方程，构建新属性参数的协方差矩阵；基于测井数据，结合振幅缩放因子、反射系数方程以及协方差矩阵，建立新属性参数模型；根据新属性参数模型，评估页岩储层位置及识别因子分布特征。本发明推导出了基于整合流体参数的VTI介质反射系数方程，获得含4个新属性参数的线性反射系数方程；实现了新属性参数的准确预测；通过建立新属性参数间的交会关系，识别了页岩储层位置及强衰减储层段，实现流体因子、衰减因子和各向异性特征的准确识别。

技术领域

本发明涉及页岩气勘探技术领域，具体涉及一种基于VTI介质的页岩储层AVO反演表征方法及系统。

背景技术

地震反演是获得地下介质内部图像、对储层进行精细描述的有效方法，也是高分辨率地震勘探的最终表现形式，地震数据反演很大程度上提高了储层表征的价值。

随着地震研究的重点由勘探逐渐向开发转移，以及由常规油气向非常规油气尤其是页岩气的转移，通过地震反演等手段来揭示地下油气藏的精细分布特征，对油气藏储层进行精雕细刻的技术也受到越来越多的关注。

页岩储层的水平薄互层构造可以等效为VTI介质。从目前的研究来看，地震反演研究主要涉及基于Zoeppritz方程的AVO反演方法和全波形反演方法。基于精确Zoeppritz方程和全波形反演方法能利用丰富的叠前信息预测模型参数，但其计算量巨大，在反演尺度和计算效率上不能满足实际油藏储层精细表征要求，尤其对于实际三维大偏移距地震数据。

针对VTI介质，传统的基于Zoeppritz方程近似反射系数是一个包含新属性参数和弱各向异性参数的表达式，虽然通过等价改写能重新表示为5个模型参数的线性表达式，但是5个未知参数造成的病态反演问题限制了方法的应用。并且，常规反演获得的新属性参数难以直接对储层识别因子进行刻画，流体因子、衰减因子等页岩储层识别因子的间接计算会引入累积误差。

发明内容

本发明的目的在于提供一种基于VTI介质的页岩储层AVO反演表征方法及系统，以解决上述背景技术中存在的至少一项技术问题。

为了实现上述目的，本发明采取了如下技术方案：

第一方面，本发明提供一种基于VTI介质的页岩储层AVO反演表征方法，包括如下步骤：

确定地震数据的振幅缩放因子；

构建VTI介质新属性参数的反射系数方程；

基于测井数据，统计新属性参数的先验信息，结合反射系数方程，构建新属性参数的协方差矩阵；

基于测井数据，结合振幅缩放因子、反射系数方程以及协方差矩阵，建立新属性参数模型；

根据新属性参数模型，评估页岩储层位置及识别因子分布特征；所述识别因子包括流体因子和衰减因子。

优选的，基于地震数据提取角度依赖的子波，结合实际井旁地震数据确定振幅缩放因子，包括：

基于地震数据提取角度依赖的地震子波；基于测井数据和VTI介质反射系数方程正演模拟角道集，并结合实际井旁地震数据确定振幅缩放因子；

其中，当地震数据信噪比高于阈值时，为角道集的每一道使用统一的振幅缩放因子；当信噪比低于阈值时，分别计算近端偏移距、中间偏移距、远端偏移距对应的振幅缩放因子。

优选的，构建VTI介质新属性参数的反射系数方程包括：

利用弹性波逆散射理论、波恩近似及稳相法，推导出基于纵横波模量的VTI介质地震波反射系数方程，对地震波反射系数方程进行整合化简，获得VTI介质新属性参数的线性反射系数方程。