叠前弹性参数反演方法、装置、计算机设备和存储介质制造方法及图纸

本发明专利技术提供一种叠前弹性参数反演方法、装置、计算机设备和存储介质，该方法包括获取叠前地震数据；获取弹性参数的反射系数；构建叠前地震数据的正演矩阵；根据正演矩阵，得到地震数据的似然函数；根据似然函数及先验信息；得到弹性参数反演的目标函数；基于最大似然估计法，将似然函数转化成最大化边缘似然函数；对最大化边缘似然函数进行求解，得到估计因子；根据估计因子及目标函数，反演得到弹性参数，通过引入自动相关判别先验信息，将弹性参数的先验信息融合到反演中，不需要对弹性参数反射系数的总体分布作出假设，使得弹性参数的反射系数保持稀疏、更好的体现岩性边界特征，从而提高叠前弹性参数反演结果分辨率和准确度。度。度。

【技术实现步骤摘要】
叠前弹性参数反演方法、装置、计算机设备和存储介质


[0001]本专利技术涉及地震勘探
，特别涉及一种叠前弹性参数反演方法、装置、计算机设备和存储介质。

技术介绍

[0002]地震勘探是利用地震信息进行石油勘探的一种方法，由于地震信息可以反映出储层参数的变化趋势，因此，可以使用该方法来预测储层参数。地震数据分为叠前和叠后两种。叠前AVA(Amplitude Versus Angle，振幅随入射度变化)反演技术，是指利用叠前地震数据的振幅入随角度变化特征，通过反演的手段获取岩石弹性参数的过程。相较于叠后反演，叠前反演由于利用了不同入射角数据的信息，可以获得更多样的弹性参数估计，在储层描述中起着至关重要的作用。因此近年来，基于叠前地震数据的反演一直是地震勘探领域中的热门话题。但由于利用的叠前地震数据与弹性参数的映射关系更加复杂，信噪比更低，叠前反演方法需要应对更多难题。其中一个难题就是降低解的非唯一性，通常的解决办法是在目标函数中加入正则化项或在贝叶斯框架下加入待估计参数的先验信息约束，达到在许多可能的解中选择一个最优或最符合假设的解的目的。<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种叠前弹性参数反演方法，其特征在于，包括：获取叠前地震数据；获取弹性参数的反射系数；基于Fatti近似方程，根据所述叠前地震数据及所述反射系数，构建叠前地震数据的正演矩阵；根据所述正演矩阵，得到地震数据的似然函数；获取自动相关判别的先验信息；基于贝叶斯理论，根据所述似然函数及所述先验信息，得到所述弹性参数反演的目标函数；基于最大似然估计法，将似然函数转化成最大化边缘似然函数；对所述最大化边缘似然函数进行求解，得到估计因子；根据所述估计因子及所述目标函数，反演得到弹性参数。2.根据权利要求1所述的叠前弹性参数反演方法，其特征在于，在所述根据所述正演矩阵，得到地震数据的似然函数的步骤之前，还包括：将所述正演矩阵中加入趋势约束，得到约束后的所述正演矩阵；所述根据所述正演矩阵，得到地震数据的似然函数的步骤，包括：根据约束后的所述正演矩阵，得到所述似然函数。3.根据权利要求2所述的叠前弹性参数反演方法，其特征在于，约束后的所述正演矩阵表达式为：其中，Ψ＝F
‑1ΛFC∈R
K
×
K
,Λ＝diag[H]，Λ为以汉宁窗函数向量作为对角元素的对角矩阵，H为汉宁窗函数向量，C为积分矩阵，LowIP＝F
‑1ΛF(log(IP))∈R
K
×1为纵波阻抗的趋势，LowIS＝F
‑1ΛF(log(IS))∈R
K
×1横波阻抗的趋势，λ1,λ2分别为趋势约束的权重参数，F为离散傅里叶变换矩阵，D为差分矩阵，IP和IS分别表示拟反演的纵波阻抗和横波阻抗参数的对数，A(θ)＝diag[A(t1,θ) A(t2,θ)
…
A(t
K
,θ)]，B(θ)＝diag[B(t1,θ) B(t2,θ)
…
B(t
K
,θ)]，diag[
…
]为对角矩阵。4.根据权利要求3所述的叠前弹性参数反演方法，其特征在于，还包括：将约束后的所述正演矩阵表达式简化成如下形式：d
(M*N+2*K)
＝G
(M*N+2*K)
×
(2K+2)
m
2K+2
+n
其中d代表G包含傅里叶变换矩阵、角度相关矩阵、阻抗趋势提取矩阵和差...

【专利技术属性】
技术研发人员：纪永祯，张渝悦，胡华锋，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：