[发明专利]弹性波阻抗反演方法、装置、电子设备及存储介质

权利要求书

1.一种弹性波阻抗反演方法，其特征在于，所述方法包括：

获取目标地震道中不同入射角的叠前地震数据；

根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式；其中，所述自适应先验分布中独立参数的数量与所述目标地震道的采样点的数量相同；

通过自适应的方法，求解所述自适应先验分布中的所有独立参数，并将其代入各个入射角的弹性阻抗反射系数反演结果的计算式中，以得到各个入射角的弹性阻抗反射系数反演结果；

根据各个入射角的弹性阻抗反射系数反演结果，获得各个入射角的弹性阻抗反演结果；

根据所有入射角的弹性阻抗反演结果，获得所述目标地震道的弹性参数反演结果。

2.根据权利要求1所述的方法，其特征在于，所述自适应的方法为边缘似然函数最大化的方法。

3.根据权利要求1所述的方法，其特征在于，根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式，包括以下步骤：

根据各个入射角的叠前地震数据，基于贝叶斯框架，构建各个入射角的似然函数；

根据各个入射角的似然函数，并以自适应先验分布为约束，构建各个入射角的后验概率分布；

根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式。

4.根据权利要求3所述的方法，其特征在于，各个入射角的叠前地震数据的噪声向量的均值为0；根据各个入射角的叠前地震数据，基于贝叶斯框架，构建各个入射角的似然函数，包括以下步骤：

根据各个入射角的叠前地震数据，基于贝叶斯的框架，构建各个入射角的似然函数为：

p(d|m，σ²，θ)＝(2πσ²)-N_exp[-(d-Gm)^T(d-Gm)/(2σ²)]；

其中，p为似然概率，d为入射角θ的叠前地震数据，G是入射角θ的子波褶积矩阵，m为入射角θ的弹性阻抗反射系数向量，σ²为入射角θ的叠前地震数据的方差，N为所述目标地震道的采样点的数量。

5.根据权利要求4所述的方法，其特征在于，所述自适应先验分布为：

其中，h为[h₁，h₂，…，h_N]向量的转置矩阵，h＝[h₁，h₂，…，h_N]^T，h_n为所述自适应先验分布中的第n个独立参数，所述自适应先验分布中独立参数的数量与所述目标地震道的采样点的数量一致，m_n为第n个采样点处的弹性阻抗反射系数；

根据各个入射角的似然函数，并以自适应先验分布为约束，构建各个入射角的后验概率分布，包括以下步骤：

根据各个入射角的似然函数，并以所述自适应先验分布为约束，构建各个入射角的后验概率分布为：

p(m|d，h，σ²，θ)＝C|∑|^-1/2exp[-1/2(m-R)^T∑^-1(m-R)]；

其中，C是常数，∑＝(H+σ^-2G^TG)^-1，R是入射角θ的弹性阻抗反射系数反演结果，弹性阻抗反射系数反演结果的计算式为R＝σ^-2∑G^Td，H为[h₁，h₂，…，h_N]向量的对角矩阵，H＝diag(h₁，h₂，…，h_N)。