[发明专利]基于贝叶斯反演框架的反射系数估计方法和装置

摘要

本发明公开了基于贝叶斯反演框架的反射系数估计方法和装置。反褶积是一种用来压缩地震子波、提高地震资料分辨率的手段，通常的反褶积方法往往要对地震子波和反射系数进行一系列假设从而获得反射系数的估计，而这些假设在很多情况下并不合理，不合理的假设阻碍了反褶积的应用效果。本发明中，基于贝叶斯反演框架进行稀疏盲反褶积，通过地质统计学，可以减少对子波和反射系数的不合理约束，从而获得更为精准的反射系数估计，提高地震记录的分辨率。

主权项

一种基于贝叶斯反演框架的反射系数估计方法，包括：步骤101，基于下式求解子波序列w＝[w0,…,wp]T：其中，R表示反射系数序列r的褶积矩阵，r＝[r1,r2,...,rM]T,r0为初始反射系数序列，M是反射系数序列r中的元素个数，Qw是由[μJr(r)+αJI(r)]T构成的对角阵，ρr(ri)表示反射系数序列r中元素ri的先验概率，C为积分算子矩阵，表示为ξ＝Cr，μ和α均为给定的约束因子，d＝[d1,d2,...,dN]T是观测到的一道地震数据,d＝Gr+n,G是N×M维子波褶积矩阵,n＝[n1,n2,...,nN]T表示观测噪声,进入步骤102；步骤102，基于当前得到的w重新计算反射系数序列r，进入步骤103；步骤103，基于当前得到的w和r判断盲反褶积目标函数J(w,r)是否收敛，如果收敛则确定当前得到的r即为估计的反射系数序列；如果不收敛，则将r0更新为当前得到的r，返回步骤101，进行下一轮迭代，直至盲反褶积目标函数J(w,r)收敛，其中，盲反褶积目标函数J(w,r)为：J(w,r)＝J0(w,r)+μJr(r)+αJI(r)，其中，