基于修正柯西分布的弹性参数反演方法.pdfVIP

第 52卷第 6期 石 油 物 探 V01．52，No．6 2013年 11月 GE()PHYSICALPROSPECTING FORPETR0LEUM NOV．，2013 文章编号：1000—1441(2013)06—0609一O8 基于修正柯西分布的弹性参数反演方法 桂金咏 ，印兴耀 ，高建虎 ，雍学善 ，李胜军 (1．中国石油勘探开发研究院西北分院，甘肃兰州730020~2．中国石油大学(华东)地球科学与技 术学院，山东青岛266580) 摘要：岩石弹性参数能够直接反映储层流体的变化 ，因而在含油气储层流体识别中应用十分广泛。传统的弹性 参数反演方法通常是基于不同角度的弹性阻抗体进行线性求解 ，对噪声较为敏感，稳健性较差。其主要原因在 于传统方法所采用的最小平方线性回归的应用前提是噪声扰动服从高斯分布，而实际上 ，在地球物理勘探领域 ， 观测资料中噪声的分布特征更接近于修正柯西分布。在传统的弹性参数反演方法原理的基础上，从最大似然估 计理论出发，分析了最小平方线性回归的应用局限，并基于修正柯西概率密度分布函数，推导出了一种高稳健性 的线性回归方法，将其应用到弹性参数反演中，取得了较好的应用效果。 关键词：弹性参数 ；反演方法；稳健性；高斯分布；修正柯西分布 DOI：10．3969／j．issi~．1000—1441．2013．06．008 中图分类号：P631．4 文献标识码：A 不同岩石以及同一种岩石含有不同流体时的 抗反演得到不同角度的弹性阻抗体，然后再进行线 弹性性质不同，具体表现为弹性参数值的变化不 性回归，才能得到最终的弹性参数反演结果。前人 同I1≈]。弹性参数如纵、横波速度、拉梅参数与泊松 根据不同的弹性参数反演需要，推导出了多种表达 n一 比等是储层预测中较为重要的参数，其反演精度一 形式的弹性阻抗方程，可统一表述为l7 直是地球物理工作者比较关心的问题。实际应用 IEc一 ()()()“c1 中考虑到抗噪性以及反演效率等要求，通常利用弹 性阻抗理论进行弹性参数的反演：基于3个及 以上 其中，A，B，C为不同的弹性参数；A。，B。，C。为 目 不同角度的弹性阻抗体，求解线性矩阵，即可反演 的层段相应的弹性参数平均值； 为相应的尺度 出纵、横波速度及密度等多种弹性参数。这也是 目 标准化因子；()，6()，f()为随入射角 变化的 前应用最为广泛的弹性参数反演方法之一[4]。但 函数。 由于地层的复杂性、噪声的干扰、截断误差的累积 由于系数项n()，6()，c()只与角度有关，因 等多种原因，常会造成弹性阻抗体存在无法预见的 此在角度相同的情况下，各采样点t处所对应的系 异常，进而对弹性参数反演的结果造成较大的影 数项也相同。对于 个采样点，(1)式两边取对数： 响，这就要求反演方法具有较高的稳健性l8]。 m (￡1，6}) m B(￡1) ， C(1) — 一 In — 常规利用弹性阻抗理论进行弹性参数反演的 L加 方法涉及到线性回归，最为常用的是最小平方线性 ( ，) B(2) C(2) m — 一 m百 m — 一 回归。最小平方线性回归的假设前提是噪声扰动 服从高斯分布，然而在地球物理勘探领域，资料中 的异常扰动并不都严格满足高斯分布的假设条件， ， ( ，) B(