一种基于对比函数的块状声波阻抗反演方法、装置及设备制造方法及图纸

本文涉及油气资源勘探地震参数反演技术领域，尤其涉及一种基于对比函数的块状声波阻抗反演方法、装置及设备。其方法包括，根据地震数据和测井数据提取地震子波；根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。通过本文实施例，减小了反演的计算量，提高了反演的精度。提高了反演的精度。提高了反演的精度。

【技术实现步骤摘要】
一种基于对比函数的块状声波阻抗反演方法、装置及设备


[0001]本文涉及油气资源勘探地震参数反演
，尤其涉及一种基于对比函数的块状声波阻抗反演方法、装置及设备。

技术介绍

[0002]声波阻抗是一种非常重要的地震属性，精确的声波阻抗可以用来储层表征和岩相分类，因此声波阻抗反演成为了油气勘探中的一项常规工作流程，由于较高的计算效率，常用确定性反演方法求解声波阻抗，如何获得块状分布且具有较高频带的声波阻抗是本领域密切关注的问题。
[0003]现有技术中常用声波阻抗反演方法分为两类：1、稀疏脉冲反演方法，需在稀疏约束下精确地反演高分辨率的反射系数。反射系数只有用零范数约束才能得到非常精确的反演结果，因此，通常用柯西分布和拉普拉斯分布等近似稀疏的分布来约束反射系数的求解过程，反演得到的反射系数的统计分布会趋近于指定的某个近似稀疏的分布，这就使得反射系数的稀疏性变弱；或用现有的一范数和二范数的组合来逼近零范数约束，然后采用复杂的混合算法求解这种混合范数目标函数，这使得求解反射系数的过程变得非常复杂。2、基于模型的反演方法，通常需要一个光滑的先验模型约束来增强反问题求解的稳定性，减弱反问题的病态性；同时，还需要加入边缘保护正则化项来保持声波阻抗的块状特性。这个过程同样会导致在目标函数中引入多个正则化参数，从而使得反演过程变得复杂。
[0004]现在亟需一种高分辨率的块状声波阻抗反演方法，解决现有技术中块状声波阻抗反演精度低且计算量大的问题。

技术实现思路

[0005]为解决现有技术中块状声波阻抗反演精度低且计算量大的问题，本文实施例提供了一种基于对比函数的块状声波阻抗反演方法、装置及设备，能够得到较高精度的反演结果、增强反演的稳定性并减小计算量。
[0006]本文提供了一种基于对比函数的块状声波阻抗反演方法，包括，
[0007]根据地震数据和测井数据提取地震子波；
[0008]根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
[0009]根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；
[0010]根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；
[0011]根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；
[0012]根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。
[0013]本文实施例还提供了一种基于对比函数的块状声波阻抗反演装置，包括，
[0014]地震子波提取单元，根据地震数据和测井数据提取地震子波；
[0015]初始声波阻抗模型建立单元，根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；
[0016]相对地震记录计算单元，根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；
[0017]正演算子构造单元，根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；
[0018]声波阻抗对比函数计算单元，根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；
[0019]声波阻抗计算单元，根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。
[0020]本文实施例还提供了一种计算机设备，包括存储器、处理器、以及存储在所述存储器上的计算机程序，所述处理器执行所述计算机程序时实现上述的方法。
[0021]本文实施例还提供了一种计算机存储介质，其上存储有计算机程序，所述计算机程序被计算机设备的处理器运行时，执行上述的方法。
[0022]利用本文实施例，根据地震子波和初始声波阻抗模型计算相对地震记录，提高了块状声波阻抗反演的分辨率，根据地震数据的最大时间采样点数以及地震子波构造正演算子，根据相对地震记录和正演算子反演的声波阻抗对比函数与真实声波阻抗相比具有更强的稀疏性，因此声波阻抗对比函数的分布更加接近于稀疏分布，不需要给定很强的稀疏约束就能够得到精确的反演结果，提高了反演的精度。此外，本文实施例反演的声波阻抗对比函数与声波阻抗之间具有简单的数学关系，通过由测井数据以及地震构造解释资料建立初始声波阻抗模型和反演的声波阻抗对比函数即可求得声波阻抗，与传统的稀疏脉冲反演方法相比，不需要在反演得到高精度的反射系数序列之后再通过反演求得具有块状特征的声波阻抗，减小了反演的计算量。
附图说明
[0023]为了更清楚地说明本文实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本文的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0024]图1所示为本文实施例一种基于对比函数的块状声波阻抗反演方法的流程图；
[0025]图2所示为本文实施例一种基于对比函数的块状声波阻抗反演装置的结构示意图；
[0026]图3所示为本文实施例一种基于对比函数的块状声波阻抗反演装置的详细结构图；
[0027]图4所示为本文实施例基于对比函数的块状声波阻抗反演方法的逻辑示意图；
[0028]图5所示为本文实施例计算机设备的结构示意图；
[0029]图6所示为本文实施例声波阻抗、声波阻抗对比函数和反射系数对应的频数统计图、核密度估计以及标准高斯分布和标准拉普拉斯分布拟合结果。
[0030]【附图标记说明】：
[0031]201、地震子波提取单元；
[0032]202、初始声波阻抗模型建立单元；
[0033]203、相对地震记录计算单元；
[0034]204、正演算子构造单元；
[0035]205、声波阻抗对比函数计算单元；
[0036]206、声波阻抗计算单元；
[0037]301、地震子波提取单元；
[0038]3011、地震子波提取模块；
[0039]3012、振幅调整模块；
[0040]302、初始声波阻抗模型建立单元；
[0041]3021、地质层位模型构建模块；
[0042]3022、声波阻抗值计算模块；
[0043]3023、初始声波阻抗模型建立模块；
[0044]303、相对地震记录计算单元；
[0045]3031、初始地震记录计算模块；
[0046]3032、相对地震记录计算模块；
[0047]304、正演算子构造单元；
[0048]3041、差分算子计算模块；
[0049]3042、子波矩阵构建模块；
[0050]3043、正演算子构造模块；
[0051]305、声波阻抗对比函数计算单元；
[0052]3051、反演目标函数构建模块；
[0053]3052、声波阻抗对比函数计算模块；
[本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于对比函数的块状声波阻抗反演方法，其特征在于，所述方法包括，根据地震数据和测井数据提取地震子波；根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型；根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录；根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子，所述正演算子用于联系所述相对地震记录与声波阻抗对比函数；根据所述相对地震记录以及所述正演算子反演所述声波阻抗对比函数；根据所述初始声波阻抗模型和所述声波阻抗对比函数计算声波阻抗。2.根据权利要求1所述的基于对比函数的块状声波阻抗反演方法，其特征在于，根据地震数据和测井数据提取地震子波之后，还包括，根据所述测井数据通过褶积模型进行正演模拟，得到模拟地震记录；将所述模拟地震记录与所述地震数据相对比，计算振幅缩放因子；根据所述振幅缩放因子，调整所述地震子波的振幅与所述地震数据中实际的地震振幅相一致。3.根据权利要求1所述的基于对比函数的块状声波阻抗反演方法，其特征在于，根据所述测井数据以及地震构造解释资料建立初始声波阻抗模型进一步包括，利用散点差值方法对所述地震构造解释资料进行差值，得到地质层位模型；利用散点差值方法对所述测井数据进行横向差值，得到地下每个点的声波阻抗值；根据所述地质层位模型和所述地下每个点的声波阻抗值建立所述初始声波阻抗模型。4.根据权利要求1所述的基于对比函数的块状声波阻抗反演方法，其特征在于，根据所述地震子波、所述初始声波阻抗模型以及所述地震数据计算相对地震记录进一步包括，根据所述地震子波以及所述初始声波阻抗模型，通过褶积模型计算初始地震记录；根据所述初始地震记录和所述地震数据通过时间域的散射理论计算相对地震记录。5.根据权利要求1所述的基于对比函数的块状声波阻抗反演方法，其特征在于，根据所述地震数据的最大时间采样点数以及所述地震子波构造正演算子的公式为，s
rela
＝G
χ
(χ
z
)G
χ
(χ
z
)＝0.5WD ln(1+χ
z
)其中，s
rela
表示所述相对地震记录，G
χ
表示所述正演算子，χ
z
表示声波阻抗对比函数，W表示由所...

【专利技术属性】
技术研发人员：李景叶，郭康康，陈小宏，马继涛，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：