基于空变阻抗约束的贝叶斯稀疏脉冲反演-CT理论与应用研究.PDFVIP

第24卷 第 6 期 CT理论与应用研究 Vol.24, No.6 2015年11月（809-817） CT Theory and Applications Nov., 2015 肖开攀, 张达, 张繁昌. 基于空变阻抗约束的贝叶斯稀疏脉冲反演[J]. CT 理论与应用研究, 2015, 24(6): 809-817. doi:10.15953/j.1004-4140.2015.24.06.06. Xiao KP, Zhang D, Zhang FC. The Bayes sparse spike inversion based on spatial variable impedance coefficient[J]. CT Theory and Applications, 2015, 24(6): 809-817. (in Chinese). doi:10.15953/j.1004-4140.2015.24.06.06. 基于空变阻抗约束的贝叶斯稀疏脉冲反演  1 2 1 肖开攀 ，张达 ，张繁昌 1.中国石油大学（华东）地球科学与技术学院，山东 青岛 266580 2.中国石化胜利物探研究院，山东 东营 257022 摘要：常规的基于贝叶斯理论的稀疏脉冲反演中，各约束项的拉格朗日算子均采用的是恒定的 常系数。反演实际资料发现，波阻抗反演剖面与钻井资料的油气显示并不能很好地对应。考虑 到反演不同地震道数据时，低频趋势模型所起的约束作用应当不同。本文在常规反演的基础上 做出了改进，假定阻抗约束系数是一个空间变量，由各地震道的实际地震数据与合成记录之间 的振幅残差来确定该道的阻抗约束系数。实际资料应用表明，改进后的反演结果更稳定，能更 准确地反映地下阻抗信息。 关键词：贝叶斯理论；稀疏脉冲反演；空变系数 doi:10.15953/j.1004-4140.2015.24.06.06 中图分类号：P631 文献标志码：A 地震反演可以实现地震剖面到岩性剖面的转换，获得地下介质的波阻抗、纵横波速度 及密度等信息，是进行储层预测和油藏描述的主要手段[1－3]。而波阻抗作为一种重要的岩石 物性指标，不仅能反映地层单元空间展布及上、下地层之间的差异性，也可直接与钻井资 料作对比[4－5]。狭义的地震反演通常就是指的波阻抗反演。 反演问题通常不适定，具有不稳定性和多解性，可用正则化和概率化方法求解[6－8]。贝 叶斯概率化方法不仅考虑了观测数据中存在的不确定问题，而且也考虑了待求解参数的先 验信息，具有明显的优越性[8－10]。本文波阻抗反演方法即为基于贝叶斯理论的概率化方法， 其实质是在地质规律、钻井和测井资料的约束下，根据野外地震记录求解地下岩层波阻抗 信息，以研究地下地层物理性质和空间结构。考虑到实际资料通常含噪及算法本身的收敛 性问题等，反演实际是在误差允许下求解出真实情况的一个逼近。当前勘探环境对反演稳 定性、反演精度和分辨率都提出了更高的要求，一些地球物理工作者近年来提出了相应的 方法和理论，如文献[9]的地层构造约束反演，文献[11]的边界保持块约束反演，文献[12] 和文献[13]的全变差正则化约束反演等，都取得了较好的应用效果。 常规反演中模型约束项的拉格朗日算子通常取恒定常数，实际数据应用发现不能准确 反映阻抗信息。本文在前人研究基础上，通过约束模型与地震资料之间的修正作用，将阻 抗约束项的权系数由常数改进为空间变量。实际资料应用表明，改善效果明显。 收稿日期：2015-10-08。 基金项目：国家重点基础研究发展计划（2013CB228604）；山东省中青年科学家奖励基金（2014BSE28009）。 810 CT理论与应用研究 24卷 1 贝叶斯框架下的反演目标函数 1.1 贝叶斯理论 贝叶斯理论的基本思路[14－16]是：已知观测样本d 、待估计参数r 的先