弹性波阻抗反演方法、装置、电子设备及存储介质.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN113933896A(43)申请公布日2022.01.14(21)申请号202010609676.6(22)申请日2020.06.29(71)申请人中国石油化工股份有限公司地址100728北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油物探技术研究院(72)发明人张渝悦纪永祯韩波王彬权(74)专利代理机构北京聿宏知识产权代理有限公司11372代理人吴大建金淼(51)Int.Cl.G01V1/30(2006.01)权利要求书3页说明书16页附图3页(54)发明名称弹性波阻抗反演方法、装置、电子设备及存储介质(57)摘要本公开提供的一种弹性波阻抗反演方法、装置、电子设备及存储介质，所述方法包括根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式；其中，所述自适应先验分布中独立参数的数量与所述目标地震道的采样点的数量相同；通过自适应的方法，得到各个入射角的弹性阻抗反射系数反演结果；根据各个入射角的弹性阻抗反射系数反演结果，获得各个入射角的弹性阻抗反演结果；根据所有入射角的弹性阻抗反演结果，获得所述目标地震道的弹性参数反演结果。该方法使得弹性波阻抗反演结果具有更好的准确性和分辨率。CN113933896ACN113933896A权利要求书1/3页1.一种弹性波阻抗反演方法，其特征在于，所述方法包括：获取目标地震道中不同入射角的叠前地震数据；根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式；其中，所述自适应先验分布中独立参数的数量与所述目标地震道的采样点的数量相同；通过自适应的方法，求解所述自适应先验分布中的所有独立参数，并将其代入各个入射角的弹性阻抗反射系数反演结果的计算式中，以得到各个入射角的弹性阻抗反射系数反演结果；根据各个入射角的弹性阻抗反射系数反演结果，获得各个入射角的弹性阻抗反演结果；根据所有入射角的弹性阻抗反演结果，获得所述目标地震道的弹性参数反演结果。2.根据权利要求1所述的方法，其特征在于，所述自适应的方法为边缘似然函数最大化的方法。3.根据权利要求1所述的方法，其特征在于，根据各个入射角的叠前地震数据，基于贝叶斯框架并以自适应先验分布为约束，构建各个入射角的后验概率分布，并根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式，包括以下步骤：根据各个入射角的叠前地震数据，基于贝叶斯框架，构建各个入射角的似然函数；根据各个入射角的似然函数，并以自适应先验分布为约束，构建各个入射角的后验概率分布；根据各个入射角的后验概率分布，获得各个入射角的弹性阻抗反射系数反演结果的计算式。4.根据权利要求3所述的方法，其特征在于，各个入射角的叠前地震数据的噪声向量的均值为0；根据各个入射角的叠前地震数据，基于贝叶斯框架，构建各个入射角的似然函数，包括以下步骤：根据各个入射角的叠前地震数据，基于贝叶斯的框架，构建各个入射角的似然函数为：22T2p(d|m，σ，θ)＝(2πσ)-Nexp[-(d-Gm)(d-Gm)/(2σ)]；其中，p为似然概率，d为入射角θ的叠前地震数据，G是入射角θ的子波褶积矩阵，m为入射角θ的弹性阻抗反射系数向量，σ2为入射角θ的叠前地震数据的方差，N为所述目标地震道的采样点的数量。5.根据权利要求4所述的方法，其特征在于，所述自适应先验分布为：T其中，h为[h1，h2，…，hN]向量的转置矩阵，h＝[h1，h2，…，hN]，hn为所述自适应先验分布中的第n个独立参数，所述自适应先验分布中独立参数的数量与所述目标地震道的采样点的数量一致，mn为第n个采样点处的弹性阻抗反射系数；根据各个入射角的似然函数，并以自适应先验分布为约束，构建各个入射角的后验概2CN113933896A权利要求书2/3页率分布，包括以下步骤：根据各个入射角的似然函数，并以所述自适应先验分布为约束，构建各个入射角的后验概率分布为：p(m|d，h，σ2，θ)＝C|∑|-1/2exp[-1/2(m-R)T∑-1(m-R)]；其中，C是常数，∑＝(H+σ-2GTG)-1，R是入射角θ的弹性阻抗反射系数反演结果，弹性阻抗-2T反射系数反演结果的计算式为R＝σ∑Gd，H为[h1，h2，…，hN]向量的对角矩阵，H＝diag(h1，h2，…，hN)。6.根据权利要求1所述的方法，其特征在于，根据各个入射角的弹性阻抗反射系数反演结果，获得各个入射角的弹性阻抗反演结果，包括以下步骤：根据各