深度域地震反演的低波数初始模型建立方法、装置及介质与流程

本发明涉及一种深度域地震反演的低波数初始模型建立方法、装置及介质，属于油气地震勘探领域。

背景技术：

1、由于深度偏移处理技术可以对地下地质构造进行更为精确的成像，对实现更加可靠的储层预测与流体识别具有独特的意义。因此，人们对用于油藏描述的、直接的深度域反演方法(如阻抗反演、弹性参数反演等)的需求越来越多。

2、然而，深度域地震数据是带限数据，缺少高、低波数成分。其中，低波数成分包含来自地下深部的反射信息，能够反映地下结构的变化趋势，这对基于模型驱动的深度域地震弹性参数反演方法和地震解释显得非常重要。

3、在深度域地震弹性参数反演中，低波数初始模型的精度和正确性会严重影响反演结果的正确性和可靠性。低波数初始阻抗模型的构建方法主要是利用测井信息在层位约束下进行外拓和插值，相应的常规方法有距离倒数加权法、克里金插值法、协同克里金插值法等。这三种常规的方法是基于地质统计学的插值方法，该类插值方法需要丰富的测井资料作为已知样本，少井或井集中都会导致对远井区低波数初始模型的约束力不强。此外，初始模型在构建时依赖于层位约束，层位解释越精细，初模型构建越精细，越有利于反演。在没有精细层位且构造复杂区域，层位的插值结果可能会引起层位抖动，进而导致初始模型精度不高，增加了反演结果的不确定性。

技术实现思路

1、针对上述技术问题，本发明提供一种深度域地震反演的低波数初始模型建立方法、装置及介质，该方法在建模的过程中利用直接在深度域提取强相关的地震属性数据与测井数据间的相关关系，使得该方法在能够快速地适用于深度域多种低波数初始模型建模的同时，具有较高的精度和较强的稳定性。

2、为实现上述目的，本发明采取以下技术方案：

3、一种深度域地震反演的低波数初始模型建立方法，包括：

4、提取目标层段的测井数据，剔除测井纵、横波速度数据和测井密度数据中的空值和异常值，得到有效测井纵波速度数据和有效测井密度数据，计算待建低波数初始模型的测井属性；

5、根据井震标定后的测井深度和地震记录深度关系，提取井旁测井段地震数据，并按照测井深度采样间距进行上采样，得到测井深度上的井旁段地震数据；

6、计算井旁段地震数据的振幅属性；

7、以井旁段地震数据的振幅属性为自变量，待建低波数初始模型的测井属性为函数值，进行线性拟合，得到拟合结果；

8、从拟合结果中选出与最优线性相关的振幅属性，得到最优斜率参数和最优截距参数；

9、根据与最优线性相关的振幅属性，计算相应的深度域地震数据体最小振幅属性；

10、根据最优斜率参数和最优截距参数，将深度域地震数据体振幅属性转化成待建低波数初始模型属性；

11、对待建低波数初始模型属性进行高斯滤波，得到最终的低波数初始模型。

12、所述的低波数初始模型建立方法，优选地，待建低波数初始模型的测井属性为测井属性纵波阻抗ip，ip＝vp·ρ，其中，vp为有效测井纵波速度；ρ为有效测井密度数据。

13、所述的低波数初始模型建立方法，优选地，计算井旁段地震数据的振幅属性，计算公式如下：

14、

15、式中，为井旁段地震数据的振幅属性；i为索引号，取值范围为1～nlog，nlog为测井段地震数据长度；fun为地震振幅属性的计算方法；为-90°相位旋转后的井旁段地震数据；l1为地震数据属性计算窗长。

16、所述的低波数初始模型建立方法，优选地，以井旁段地震数据的振幅属性为自变量，待建低波数初始模型的测井属性为函数值，进行线性拟合，拟合过程如下：

17、

18、式中，iplog为测井纵波阻抗属性；a和b分别为线性拟合的斜率参数和截距参数。

19、所述的低波数初始模型建立方法，优选地，根据与最优线性相关的振幅属性，计算相应的深度域地震数据体最小振幅属性，具体步骤如下：

20、

21、式中，ampmin(j,k)表示第j道，第k个采样点的最小振幅属性；l2为地震数据属性计算窗长；funmin为地震最小振幅属性的计算方法；为相位旋转-90°后的第j道地震记录。

22、所述的低波数初始模型建立方法，优选地，根据最优斜率参数和最优截距参数，将深度域地震数据体振幅属性转化成待建低波数初始模型属性，具体步骤如下：

23、ip′(j,k)＝a′·ampmin(j,k)+b′

24、式中，ip′为利用最小地震振幅属性转换成的纵波阻抗属性；a′为最优斜率参数；b′为最优截距参数。

25、所述的低波数初始模型建立方法，优选地，对待建低波数初始模型属性进行高斯滤波，得到最终的低波数初始模型，具体步骤如下：

26、

27、式中，g为二维高斯窗函数；m为窗长；σ为高斯分布的标准差；u,v为索引号；取值范围均为1～m；ip′为待建低波数纵波阻抗初始模型属性；ip″为低波数纵波阻抗初始模型。

28、本发明由于采取以上技术方案，其具有以下优点：

29、1、本发明方法在充分利用测井和地震属性数据的强相关性的同时，避免了测井数据因重采样而造成的数据信息丢失，能够在测井尺度上极大地利用测井信息；避免了因少井或井集中而导致在远井区对低波数初始模型约束力不够；避免了因层位解释不细致而导致初始模型构建不准确。

30、2、本发明方法在建立起测井数据属性和井旁相应段地震数据属性的线性拟合关系后，能够根据地震数据属性和待建模的数据属性之间的强线性关系快速地建立起较为准确的低波数初始模型，从而提高在复杂构造条件下初始模型的精确度，提高反演结果的确定性。

技术特征：

1.一种深度域地震反演的低波数初始模型建立方法，其特征在于，包括：

2.根据权利要求1所述的低波数初始模型建立方法，其特征在于，待建低波数初始模型的测井属性为测井属性纵波阻抗ip，ip＝vp·ρ，其中，vp为有效测井纵波速度；ρ为有效测井密度数据。

3.根据权利要求2所述的低波数初始模型建立方法，其特征在于，计算井旁段地震数据的振幅属性，计算公式如下：

4.根据权利要求3所述的低波数初始模型建立方法，其特征在于，以井旁段地震数据的振幅属性为自变量，待建低波数初始模型的测井属性为函数值，进行线性拟合，拟合过程如下：

5.根据权利要求4所述的低波数初始模型建立方法，其特征在于，根据与最优线性相关的振幅属性，计算相应的深度域地震数据体最小振幅属性，具体步骤如下：

6.根据权利要求5所述的低波数初始模型建立方法，其特征在于，根据最优斜率参数和最优截距参数，将深度域地震数据体振幅属性转化成待建低波数初始模型属性，具体步骤如下：

7.根据权利要求6所述的低波数初始模型建立方法，其特征在于，对待建低波数初始模型属性进行高斯滤波，得到最终的低波数初始模型，具体步骤如下：

8.一种深度域地震反演的低波数初始模型建立装置，其特征在于，包括：

9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-8任意一项所述深度域地震反演的低波数初始模型建立方法的步骤。

10.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1-8任意一项所述深度域地震反演的低波数初始模型建立方法的步骤。

技术总结
本发明涉及一种深度域地震反演的低波数初始模型建立方法、装置及介质，方法包括：提取目标层段的测井数据，计算待建低波数初始模型的测井属性；根据井震标定结果，提取井旁测井段地震数据，并按照测井深度采样间距进行上采样，得到测井深度上的井旁段地震数据；计算井旁段地震数据的振幅属性；以井旁段地震数据的振幅属性为自变量，待建低波数初始模型的测井属性为函数值，进行线性拟合；从拟合结果中选出与最优线性相关的振幅属性，得到最优斜率参数和最优截距参数；计算相应的深度域地震数据体最小振幅属性；将深度域地震数据体振幅属性转化成待建低波数初始模型属性；对待建低波数初始模型属性进行高斯滤波，得到最终的低波数初始模型。

技术研发人员：王建花,牛聪,黄饶,凌云,陈剑军,王艳冬,张洪亮,王迪,郑颖,张玉华,李超
受保护的技术使用者：中海石油（中国）有限公司
技术研发日：
技术公布日：2024/1/15