Bayesian prestack seismic inversion with a self-adaptive Huber-Markov random-field edge protection scheme  

Huber-Markov随机场自适应边缘保护贝叶斯地震资料叠前反演(英文)

作　　者：田玉昆 周辉[1] 陈汉明 邹雅铭 关守军 

机构地区：[1]中国石油大学(北京)油气资源与探测国家重点实验室,CNPC物探重点实验室,北京102249

出　　处：《Applied Geophysics》2013年第4期453-460,512,共9页应用地球物理（英文版）

基　　金：supported by the National Basic Research Program of China(973 Program)(No.2013CB228603);National Science and Technology major projects(No.2011ZX05024 and 2011ZX05010);the National Natural Science Foundation of China(No.41174119)

摘　　要：Seismic inversion is a highly ill-posed problem, due to many factors such as the limited seismic frequency bandwidth and inappropriate forward modeling. To obtain a unique solution, some smoothing constraints, e.g., the Tikhonov regularization are usually applied. The Tikhonov method can maintain a global smooth solution, but cause a fuzzy structure edge. In this paper we use Huber-Markov random-field edge protection method in the procedure of inverting three parameters, P-velocity, S-velocity and density. The method can avoid blurring the structure edge and resist noise. For the parameter to be inverted, the Huber- Markov random-field constructs a neighborhood system, which further acts as the vertical and lateral constraints. We use a quadratic Huber edge penalty function within the layer to suppress noise and a linear one on the edges to avoid a fuzzy result. The effectiveness of our method is proved by inverting the synthetic data without and with noises. The relationship between the adopted constraints and the inversion results is analyzed as well.由于地震数据频带范围的有限性以及噪声、正演近似等因素的影响,地震反演是一个高度不适定的反问题。为了得到稳定唯一的解,通常在反演目标函数中加入光滑化约束,常用的是Tikhonov正则化约束。这种正则化对解是全局光滑的,在构造边缘位置会产生模糊。为了解决边缘模糊的问题,同时保证抗噪的作用,本文采用Huber-Markov随机场约束边缘保护的方法求解叠前纵横波速度和密度三参数反演问题。对于待反演的参数,通过Markov邻域系统建立起纵向和横向上的约束。通过Huber边缘惩罚函数阈值的合理设置,达到在层内区域采用二次型函数抑制噪声得到平滑的结果;在边缘位置则采用线性函数,避免边缘的模糊化。通过无噪和有噪合成地震记录的反演,验证了方法的正确性,分析了正则项与反演效果的关系。

关 键 词：Huber edge punishment function markov random-field bayesian framework prestack inversion 

分 类 号：P631.44[天文地球—地质矿产勘探]