[发明专利]叠前弹性参数反演方法、装置、计算机设备和存储介质

权利要求书

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横波阻抗的趋势，λ₁,λ₂分别为趋势约束的权重参数，F为离散傅里叶变换矩阵，D为差分矩阵，IP和IS分别表示拟反演的纵波阻抗和横波阻抗参数的对数，

A(θ)＝diag[A(t₁,θ) A(t₂,θ)…A(t_K,θ)]，

B(θ)＝diag[B(t₁,θ) B(t₂,θ)…B(t_K,θ)]，diag[…]为对角矩阵。

4.根据权利要求3所述的叠前弹性参数反演方法，其特征在于，还包括：将约束后的所述正演矩阵表达式简化成如下形式：

d_(M*N+2*K)＝G_{(M*N+2*K)×(2K+2)}m_2K+2+n

其中d代表G包含傅里叶变换矩阵、角度相关矩阵、阻抗趋势提取矩阵和差分矩阵，m代表待估计的纵横波阻抗参数组成的向量，n代表噪声向量。

5.根据权利要求1所述的叠前弹性参数反演方法，其特征在于，所述获取弹性参数的反射系数的步骤，包括：

构建纵横波阻抗的弹性参数模型；

根据所述弹性参数模型，得到所述反射系数。

6.根据权利要求1所述的叠前弹性参数反演方法，其特征在于，所述目标函数的表达式为：

p(m|d,h,σ²,θ)＝C|Σ|^-1/2exp[-1/2(m-μ)^TΣ^-1(m-μ)]

其中，Σ＝(H+σ^-2G^TG)^-1，μ＝σ^-2ΣG^Td，C为一个常数，μ为待估计纵横波阻抗弹性参数的后验概率均值，H＝diag(h₁,h₂,…,h_2K)为一个对角矩阵。

7.根据权利要求1所述的叠前弹性参数反演方法，其特征在于，所述先验信息的表达式为：

其中，h＝[h₁,h₂,…,h_2K]^T包含2K个独立的参数。