[发明专利]一种叠前地震反演方法及装置

说明书

本申请实施例公开了一种叠前地震反演方法及装置。所述方法包括：获取目的工区的叠前地震数据和测井数据；基于所述叠前地震数据和测井数据，确定所述目的工区的初始正演地震数据；基于所述叠前地震数据和所述初始正演地震数据，确定误差函数项；基于所述误差函数项，以及与所述目的工区的弹性参数的改变量相关联的第一特征项、第二特征项和第三特征项，构建所述目的工区的反演目标函数；基于所述反演目标函数，确定所述目的工区的弹性参数的目标值。可以提高叠前地震反演所得到的反演结果的准确度。

技术领域

本申请涉及地震勘探技术领域，特别涉及一种叠前地震反演方法及装置。

背景技术

将地表接收到的来自地下同一反射点的地震波数据经过动校正后进行叠加，叠加后的地震波数据可以称为叠后地震数据。叠加前的地震波数据可以称为叠前地震数据。相对于叠后地震数据，叠前地震数据保留了更多的地下地质信息。例如，与储层中某一地层位置处的纵波速度、横波速度或密度等相关的信息。正因如此，基于叠前地震数据的叠前AVO(Amplitude Versus Offset，振幅随偏移距的变化)反演可以提供储层中目的层位的纵波速度、横波速度或密度，其反演精度对储层预测和油气藏描述有至关重要的作用。

佐普利兹(Zoeppritz)方程是叠前AVO反演的理论基础。但由于其复杂性，通常将Zoeppritz方程简化后得到用于表示地震波反射系数与弹性参数之间的映射关系的Zoeppritz近似公式，基于Zoeppritz近似公式进行叠前AVO反演。弹性参数可以包括：纵波速度、横波速度和密度。

现有技术中通常采用的叠前AVO反演方法的主要过程是：获取目的工区的叠前地震数据；在叠前地震数据的纵波入射角为30度以内的假设条件下，基于Zoeppritz近似公式，可以构建目的工区中目的层位的地震波反射系数与弹性参数之间的映射关系；基于叠前地震数据和映射关系，进行AVO反演处理，从而确定目的工区中目的层位的弹性参数的目标值，即纵波速度、横波速度和密度的目标值。

发明人发现现有技术中至少存在如下问题：针对纵波入射角大于30度的叠前地震数据，基于Zoeppritz近似公式得到的地震波反射系数的准确度较低，可能导致最终得到的反演结果存在较大偏差，即得到目的工区中目的层位的纵波速度、横波速度和密度的准确度较低。

发明内容

本申请实施例的目的是提供一种叠前地震反演方法及装置，以提高叠前AVO反演所得到的反演结果的准确度。

为解决上述技术问题，本申请实施例提供一种叠前地震反演方法及装置是这样实现的：

一种叠前地震反演方法，包括：

获取目的工区的叠前地震数据和测井数据；基于所述叠前地震数据和测井数据，确定所述目的工区的初始正演地震数据；

基于所述叠前地震数据和所述初始正演地震数据，确定误差函数项；

基于所述误差函数项，以及与所述目的工区的弹性参数的改变量相关联的第一特征项、第二特征项和第三特征项，构建所述目的工区的反演目标函数；第一特征项用于调整所述反演目标函数的矩阵条件数；第二特征项用于控制所述弹性参数的改变量的范围；第三特征项用于对所述弹性参数的值进行平滑处理；

基于所述反演目标函数，确定所述目的工区的弹性参数的目标值。

优选方案中，所述基于叠前地震数据和测井数据，确定所述目的工区的初始正演地震数据，包括：

基于所述叠前地震数据，确定所述目的工区的地震子波数据；

基于所述测井数据，确定所述目的工区中各个层位的弹性参数的初始值；

根据预设纵波入射角和所述各个层位的弹性参数的初始值，确定所述目的工区中各个层位的地震波反射系数的初始值；