[发明专利]基于叠前叠后联合反演的井段气层识别方法及系统

说明书

本发明公开了一种基于叠前叠后联合反演的井段气层识别方法及系统，包括获取井段的叠前地震数据，所述叠前地震数据包括二维地震剖面；获取二维地震剖面中井段的角度道集；从所述角度道集中选择角度子波；对各角度子波进行叠前地震反演，得到对应的横波速度；将所述横波速度与测井数据进行吻合，根据吻合结果选择最优角度子波；通过Zoeppritz方程计算井段的反射系数；对所述反射系数、最优角度子波进行褶积，得到直接反演横波速度；根据所述直接反演横波速度和约束公式计算得到纵横波速度比；根据所述纵横波速度比和二维地震剖面中单个测井的含气饱和度计算二维地震剖面的含气饱和度；通过所述二维地震剖面的含气饱和度对气层进行识别。

技术领域

本发明涉及油气地球物理勘探技术领域，具体涉及一种基于叠前叠后联合反演的井段气层识别方法及系统。

背景技术

叠前地震反演是一种基于褶积模型的反演，它基于振幅随炮检距变化理论，可以同步获得纵波速度、横波速度和密度等多种弹性参数，可以有效识别储层流体类型和构造特征。因此，叠前地震反演是目前应用最广泛、发展最成熟的地震反演技术之一，在油气资源的勘探和开发过程中发挥重要作用。

叠后地震反演是油气储层预测的核心环节，反演过程中综合利用了地质、地震和测井资料，使波阻抗反演结果实现了横向连续性与纵向分辨率的统一，为研究有利储层的空间变化及分布特征提供重要依据。叠后地震反演方法主要有道积分反演、递推反演、约束稀疏脉冲反演和地质统计学反演等。道积分反演和递推反演直接利用地震数据，从上至下依次做积分运算或递推运算，得到地层的相对波阻抗，约束稀疏脉冲反演以地震数据为基础，利用测井数据和初始模型进行约束，通过道合并得到地层的绝对波阻抗。

目前，在储层区域评价中多参数同步反演是一种必不可少的特殊处理手段，叠前地震资料的多参数同步反演是行之有效的方法之一。但是，由于多参数之间的量级的差异、敏感性的差异，致使同步获得多参数的精度难以保证，尤其是对于一些复杂储层和品质一般的叠前地震资料。为此，可以尝试叠后和叠前联合多参数反演。以实际资料的气层识别为契机，开展叠后和叠前联合多参数反演，通过叠后反演获得纵波速度和密度参数，叠前反演获得横波速度、纵横波速度比、泊松比等参数。

目前，大部分井钻会遇到气层，但部分井上也会遇到气水同层、含气水层及纯水层。目前，研究的目的是找到有效气层的分布范围，而气水同层、含气水层虽然含气，但商业价值较低，是后续打开发井需要规避的风险。为此，需要通过新方法的研究，计算含气饱和度，以便有效的区分气层与气水同层。但是，直接通过地震资料定量计算含气饱和度的难度非常大，其精度无法满足要求，只能识别气层(或气水同层)与水层。

发明内容

针对现有技术的不足，本发明的目的在于提供一种基于叠前叠后联合反演的井段气层识别方法及系统，以解决现有技术中存在的计算含气饱和度精度不高的问题。

为解决上述技术问题，本发明采用的技术方案为：

一种井段气层识别方法，包括如下步骤：

获取井段的叠前地震数据，所述叠前地震数据包括二维地震剖面；

获取二维地震剖面中井段的角度道集；

从所述角度道集中选择角度子波；

对各角度子波进行叠前地震反演，得到对应的横波速度；

将所述横波速度与测井数据进行吻合，根据吻合结果选择最优角度子波；

通过Zoeppritz方程计算井段的反射系数；

对所述反射系数、最优角度子波进行褶积，得到直接反演横波速度；

根据所述直接反演横波速度和约束公式计算得到纵横波速度比；

根据所述纵横波速度比和二维地震剖面中单个测井的含气饱和度计算二维地震剖面的含气饱和度；