二维大地电磁数据的聚焦反演算法探讨.pdfVIP

　2007 年 6 月 　　　　　　　　　　　　　　石油地球物理勘探　　　　　　　　　　　　 第 42 卷 　第 3 期 　 ·非地震 · 二维大地电磁数据的聚焦反演算法探讨 刘小军 　王家林 　陈　冰 　于 　鹏 ( 同济大学海洋地质国家重点实验室) 刘小军 ,王家林 ,陈冰 ,于鹏. 二维大地电磁数据的聚焦反演算法探讨. 石油地球物理勘探 ,2007 ,42 (3) :338～342 摘要 　文中简要分析了当前地球物理反演中几种常用的模型 目标函数 ,并根据聚焦反演成像方法的基本理论 , 在大地电磁反演 目标函数中引入了一个新的稳定泛函 , 以代替传统最小构造或最大平滑准则下的模型 目标函 数 。通过对 自适应正则化共轭梯度法解最优化问题过程的研究 , 以及对几个典型模型的试算 ,并将其结果与传 统反演算法结果进行对比分析 ,表明了本文方法的正确性和较高的计算精度 。最后对雷 —琼拗陷一条实测的测线 的 M T 资料进行聚焦反演成像处理 ,并与同一测线的地震资料进行对比表明,聚焦反演结果具有很高的准确性 。 关键词 　大地电磁 　目标函数 　稳定泛函　聚焦反演成像 　共轭梯度法 际资料处理中取得了较好的效果 。上述方法在选取 1 　引言 模型目标函数时都是基于最小构造或最大平滑准 则 ,这种思路在迭代次数 、收敛性等方面仍存在一定 大地 电磁数据反演是整个数据处理的核心部 不足 ,主要在于对地下电性出现突变情况下 ,所得到 分 ,将直接影响后期地质构造解释的效果 。传统的 的地电断面图不能很准确地分辨出地质体的分界 反演方法大多是建立在正则化理论基础上 ,其基本 面 ,给后期的地质解释带来了一定的困难 。 [ 7 ] 思想是在反演 目标函数中加入一个约束泛函来减小 聚焦反演成像 是 Port niaguine 等提出的一种 解的奇异性 。现在常用的约束稳定泛函有模型参数 对不同地质体界面进行成像的方法 ,它的基本思想 的范数 、最大平滑稳定泛函和最小模型泛函等 。其 是将最小支撑泛函的稳定器与惩罚泛函相结合 , 引 中比较著名的有 Co n st able[ 1 ] 开发的 OCCA M 法 , 入一个新的稳定泛函来描述模型参数变化较大和不 它是在极小可能构造下寻找符合数据的光滑模型 。 连续的区域 ,这样在反演中模型 目标泛函会集中到 在解 目标函数的方法上 , 国内外很多学者从提高反 最小的面积 ,从而显示出物性差别较大的界面 。这 演速度 、求优化问题的方法等方面提出了很多改进 一新方法已在非地震资料处理中得到了较好的应 [2 ] ( ) 用 ,在重力反演中对局部密度异常反应较灵敏[ 8 ] ,用 方法 ,如 Smit h 等 的快速松弛法 RR I ,其基本思 想是假设电阻率主要以垂向变化为主 ,这样大大加 于三维井间电磁模型试算精度较高[9 ] 。本文根据聚 [3 ] 焦反演的特性 ,将其应用到 M T 数据反演的目标函 快了求取灵敏矩阵的速度 ; Ro di 等 提出了非线性 共轭梯度法 ,抛弃了线性反演的框架 ,不需要直接求