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

说明书

本公开提供的一种弹性波阻抗反演方法、装置、电子设备及存储介质，所述方法包括根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式；其中，所述自适应先验分布中独立参数的数量与所述目标地震道的采样点的数量相同；通过自适应的方法，得到各个入射角的弹性阻抗反射系数反演结果；根据各个入射角的弹性阻抗反射系数反演结果，获得各个入射角的弹性阻抗反演结果；根据所有入射角的弹性阻抗反演结果，获得所述目标地震道的弹性参数反演结果。该方法使得弹性波阻抗反演结果具有更好的准确性和分辨率。

技术领域

本公开涉及地震勘探技术领域，特别地涉及一种弹性波阻抗反演方法、装置、电子设备及存储介质。

背景技术

弹性波阻抗反演方法是叠前反演方法之一，是获得地下岩石弹性参数的最重要手段。相较于叠后反演，叠前反演方法由于利用了不同入射角数据的信息，可以获得更多样的弹性参数估计，在储层描述中起着至关重要的作用。弹性波阻抗的概念可以理解为叠后波阻抗概念在叠前的拓展，引入了不同入射角的概念，叠后波阻抗可以看作是入射角为0度，即地震波垂直入射时的弹性阻抗。基于弹性阻抗的叠前反演方法的实施步骤可以分为两步：第一步是基于不同入射角的叠前地震数据获取不同入射角的弹性阻抗；第二步是基于不同入射角的弹性阻抗获取弹性参数的估计结果。在第一步中，常见的约束方法是在贝叶斯框架下，假设弹性阻抗反射系数具有某种分布特征，对弹性波阻抗反演进行约束，获得弹性波阻抗反演结果。常用的分布特征如高斯分布、柯西分布和修正柯西分布等，这些分布的表征大体是通过两个参数进行，如高斯分布受均值和方差的控制，柯西分布受位置参数和尺度参数的控制。

进行弹性波阻抗反演时，先利用基于测井资料计算获得的弹性阻抗反射系数的分布拟合出这两个参数，进而在数据的反演中进行约束，一旦拟合出两个参数后相应分布即固定，并运用于所有地震道的弹性波阻抗反演中进行约束。拟合出的参数本质上与正则化参数内涵一致。但由于弹性阻抗反射系数受到不同地区(即不同地震道)的地质因素影响，因此采用某一种固定的分布来约束反演，会产生误差；并且，弹性阻抗反射系数并不完全符合这种假设的分布，也会产生误差。

发明内容

针对上述问题，本公开提供一种弹性波阻抗反演方法、装置、电子设备及存储介质，解决了现有技术中采用固定参数进行约束导致的反演结果误差较大的问题。

第一方面，本公开提供了一种弹性波阻抗反演方法，所述方法包括：

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

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

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

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

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

根据本公开的实施例，可选的，上述弹性波阻抗反演方法中，所述自适应的方法为边缘似然函数最大化的方法。

根据本公开的实施例，可选的，上述弹性波阻抗反演方法中，根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式，包括以下步骤：

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