基于Metropolis抽样的弹性阻抗随机反演.pdfVIP

第39卷第 1期 物 探 与 化 探 Vo1．39．No．1 2015年 2月 GEOPHYSICAL ＆GEOCHEMICALEXPLORATION Feb．，2015 doi110．11720／wtyht．2015．1-33 孙瑞莹，印兴耀，王保丽，等．基于Metropolis抽样的弹性阻抗随机反演 [J]．物探与化探，2015，39(1)：203—210．http：／／doi．org／10．11720／wtyht． 2015．1．33 SunRY，YinXY，WangBL，eta1．StochasticinversionofelasticimpedancebasedonMetropolissamplingalgorithm[J]．GeophysicalandGeochemical Exploration，2015，39(1)：203—210．http：／／doi．org／lO．11720／wtyht．2015．1．33 基于 Metropolis抽样的弹性阻抗随机反演 孙瑞莹，印兴耀 ，王保丽，张广智 (中国石油大学(华东)地球科学与技术学院，山东 青岛 266580) 摘 要：笔者提出了基于Metropolis抽样的弹性阻抗随机反演方法 ，该方法是一种基于蒙特卡洛的非线性反演方法， 能够有效地融合测井资料中的高频信息，提高反演结果的分辨率。应用贝叶斯理论框架，首先通过快速傅里叶滑 动平均模拟算法 (fastFouriertransform．movingaverage，FFT—MA)和逐渐变形算法 (gradualdeformationmethod，GDM) 得到基于地质统计学的先验信息 ，然后用Metropolis抽样算法对后验概率密度进行抽样，得到反演 问题的解。其中 FFr—MA模拟作为一种高效的频率域模拟方法 ，结合GDM后，在保持模拟空间结构不变的前提下，可以连续修改储 层模型，直至满足实际观测地震记录。数值试验表明：Ffvr—MA模拟有效地提高了计算效率，融人 GDM更新算法， 可以保证反演结果有效地收敛，并且反演结果与理论模型吻合较好 ，具有较高的分辨率，该方法采用两步法反演弹 性参数，运算速度较快。 关键词：快速傅里叶滑动平均模拟算法；逐渐变形算法 ；弹性阻抗；高分辨率；Metropolis抽样 中图分类号：P631．4 文献标识码 ：A 文章编号：1000—8918(2015)叭一0203—08 地震勘探 已经从简单的构造识别转 向复杂构 原因，使用范围比较小_7]。但是近些年，由于随机 造、薄储层以及老油田剩余油的研究。勘探难度的 反演能够在空间相关性和井约束下，模拟 出地震频 加大，也对能得到更高信噪比和分辨率的地震资料 带以外的信息，分辨率高于常规的确定性反演方法。 提出了挑战。Chopra等发现实际的地震子波主频 随机反演又逐渐被人们所关注。而弹性阻抗反演采 大约在 30Hz。很难区分 25m 以下的薄层顶底_】]。 用叠前角度部分叠加道集，反演过程中考虑了子波 但是 ，工业 目标地质体的厚度一般在 10111甚至 10 随偏移距的变化，相 比叠后波阻抗反演更能反映地 m以下，所以提高分辨率以反演薄层信息，在储层识 层岩性和流体特征，比叠前 AVO反演方法的抗噪能 别和油藏描述中起到了重要的作用 ．并且随着地球 力更强，因而在地震勘探领域具有广泛应用。1999 物理研究的深人及反演精度 的提高，非线性反演方 年，Connolly[8发表弹性阻抗的文章后。在国内外地 法更容易得到真值。 球物理学界掀起 了弹性阻抗反演研究的热潮 ，提出 随机反演在储层识别和油藏描述中具有重要作 了许多扩展的弹性阻抗概念 ，并对阻抗域的储层预 用。许多学者对确定性反演和随机反演进行了研究 测和流体识别等做 了大量的分析研究 于是为 对 比_2j，结果表