基于Metropolis抽样的弹性阻抗随机反演-物探与化探.PDFVIP

　 第３９卷第 １期 物　 探　 与　 化　 探 Ｖｏｌ．３９，Ｎｏ．１　 　 ２０１５年２月 ＧＥＯＰＨＹＳＩＣＡＬ ＆ ＧＥＯＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ Ｆｅｂ．，２０１５　 ｄｏｉ：１０．１１７２０／ ｗｔｙｈｔ．２０１５．１．３３ － 孙瑞莹，印兴耀，王保丽，等．基于Ｍｅｔｒｏｐｏｌｉｓ抽样的弹性阻抗随机反演［Ｊ］．物探与化探，２０１５，３９（１）：２０３ ２１０．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ． ２０１５．１．３３ Ｓｕｎ Ｒ Ｙ，Ｙｉｎ Ｘ Ｙ，Ｗａｎｇ Ｂ Ｌ，ｅｔ ａｌ．Ｓｔｏｃｈａｓｔｉｃ ｉｎｖｅｒｓｉｏｎ ｏｆ ｅｌａｓｔｉｃ ｉｍｐｅｄａｎｃｅ ｂａｓｅｄ ｏｎ Ｍｅｔｒｏｐｏｌｉｓ ｓａｍｐｌｉｎｇ ａｌｇｏｒｉｔｈｍ［Ｊ］．Ｇｅｏｐｈｙｓｉｃａｌ ａｎｄ Ｇｅｏｃｈｅｍｉｃａｌ － Ｅｘｐｌｏｒａｔｉｏｎ，２０１５，３９（１）：２０３ ２１０．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１５．１．３３ 基于Ｍｅｔｒｏｐｏｌｉｓ抽样的弹性阻抗随机反演 孙瑞莹，印兴耀，王保丽，张广智 （中国石油大学（华东） 地球科学与技术学院，山东 青岛　 ２６６５８０） 摘 要：笔者提出了基于Ｍｅｔｒｏｐｏｌｉｓ抽样的弹性阻抗随机反演方法，该方法是一种基于蒙特卡洛的非线性反演方法， 能够有效地融合测井资料中的高频信息，提高反演结果的分辨率。 应用贝叶斯理论框架，首先通过快速傅里叶滑 动平均模拟算法（ｆａｓｔ Ｆｏｕｒｉｅｒ ｔｒａｎｓｆｏｒｍ⁃ｍｏｖｉｎｇ ａｖｅｒａｇｅ，ＦＦＴ⁃ＭＡ）和逐渐变形算法（ｇｒａｄｕａｌ ｄｅｆｏｒｍａｔｉｏｎ ｍｅｔｈｏｄ，ＧＤＭ） 得到基于地质统计学的先验信息，然后用Ｍｅｔｒｏｐｏｌｉｓ抽样算法对后验概率密度进行抽样，得到反演问题的解。 其中 ＦＦＴ⁃ＭＡ模拟作为一种高效的频率域模拟方法，结合ＧＤＭ后，在保持模拟空间结构不变的前提下，可以连续修改储 层模型，直至满足实际观测地震记录。 数值试验表明：ＦＦＴ⁃ＭＡ模拟有效地提高了计算效率，融入 ＧＤＭ更新算法， 可以保证反演结果有效地收敛，并且反演结果与理论模型吻合较好，具有较高的分辨率，该方法采用两步法反演弹 性参数，运算速度较快。 关键词：快速傅里叶滑动平均模拟算法；逐渐变形算法；弹性阻抗；高分辨率；Ｍｅｔｒｏｐｏｌｉｓ抽样 － － － 中图分类号：Ｐ６３１．４　 　 　 文献标识码：Ａ　 　 　 文章编号：１０００ ８９１８（２０１５）０１ ０２０３ ０８ ［７］ 　 　 地震勘探已经从简单的构造识别转向复杂构 原因，使用范围比较小 。 但是近些年，由于随机 造、薄储层以及老油田剩余油的研究。 勘探难度的 反演能够在空间相关性和井约束下，模拟出地震频 加大，也对能得到更高信噪比和分辨率的地震资料 带以外的信息，分辨率高于常规的确定性反演方法， 提出了挑战。 Ｃｈｏｐｒａ 等发现实际的地震子波主频 随机反演又逐渐被人们所关注。 而弹性阻抗反演采 ［１］ 大约在３０ Ｈｚ，很难区分２５ ｍ 以下的薄层顶底 。 用叠前角度部分叠加道集，反演过程中考虑了子波 但是，工业目标地质体的厚度一般在 １０ ｍ甚至 １０ 随偏移距的变化，相比叠后波阻抗反演更能反映地 ｍ 以下，