基于zondmt2d大地电磁反演算法对比研究参考.pdfVIP

78 西部探矿工程 2017年第3期 基于ZondMT2D大地电磁反演算法对比研究 刘 祥* （东华理工大学地球物理与测控技术学院，江西南昌330013） 摘 要：在大地电磁反演中，如何选择合适的反演算法是反演过程的关键，将直接影响着反演效果的 好坏。针对当前大地电磁反演软件反演算法种类繁多的现状，基于大地电磁二维数据处理软件 ZondMT2D平台，开展了平滑约束（Smoothness）算法、奥克姆（Occam）算法、聚焦（Focused）算法、马 奎特（Marquardt）算法等4种经典大地电磁反演算法的对比研究，分析4种反演算法在不同模型下的 应用效果，总结各算法的特点及缺陷，为实际资料处理提供参考。 关键词：大地电磁；反演算法；对比；分析 中图分类号：P631 文献标识码：A 文章编号：1004-5716(2017)03-0078-05 1 概述 m——模型参数向量； 大地电磁测深（MT）具有较大的勘探深度，并且不 μ——正则化参数； 受高阻层屏蔽影响，对低阻层也有较高的分辨能力，已 Δf——实测数据与预测数据之差。 经是深部地质构造研究和矿产资源勘查的基本探测方 2.2 奥克姆算法 法。大地电磁二维反演算法研究一直是大地电磁测深 奥克姆算法是加入平滑算子的最小二乘法和对比 方法研究的一个重要部分，选择合适的反演算法是反演 度最小化的一种反演方法。该方法由Constable等人[2] 过程的关键，将直接影响着其应用效果。针对当前大地 在1987年提出，它为了寻找在有极小可能构造意义下 电磁反演软件反演算法种类繁多的现状，本文基于大地 拟合数据的模型，模型与原始测量曲线最大拟合且模 电磁二维数据处理软件ZondMT2D平台，从平滑约束 型数据最平滑，需引入模型粗糙度压制非数据的模型 （Smoothness）、奥克姆（Occam）、聚焦（Focused）、马奎 构造求得模型的最光滑解。它引入Lagrange乘子构造 [3] 特（Marquardt）等4种经典大地电磁反演算法的基本原 一个不受约束的目标函数U(m) ： 理出发，通过对比其在几种典型的地电模型的收敛速 ( ) 2 2 -1 ( ) 2 2 U m = ∂ m + ∂ m +μ Wd-WF m -X      y Z { *} 度、误差下降情况及反演结果等方面，以此总结各算法 (2) 的特点及缺陷，为实际资料处理提供参考。 式中：W——数据标准差； 2 大地电磁二维反演算法 d——观测数据； 2.1 平滑约束算法 F(m)——模型响应； 平滑约束算法又称为最小偏置算法，是通过加入