地形对不同反演的影响-地质与勘探.pdfVIP

第５０卷　第 １期 地质与勘探 Ｖｏｌ．５０　Ｎｏ．１ ２０１４年 １月 　　　　　　　　　　　　　ＧＥＯＬＯＧＹＡＮＤＥＸＰＬＯＲＡＴＩＯＮ　　　　　　　　　　　　　 Ｊａｎｕａｒｙ，２０１４ 地形对不同反演的影响 朱崇利 （成都理工大学地球探测与信息技术教育部重点实验室，四川成都　６１００５９） ［摘　要］根据观测数据和模型参数之间的关系，通常将求误差函数的最优问题的方法通常分为线 性反演和非线性反演；线性反演要求有相当准确的背景介质模型的条件，这和地球物理问题中的实际 情况相矛盾。特殊情况下，线性反演甚至出现发散的情况。而非线性反演不用求目标函数的偏导数，全 局随机搜索模型空间的最优解的办法，避免了问题的线性化，在山区复杂地形条件下，通过对比论证，非 线性反演能够获取更合理的真实地电模型，从而说明由于非线性反演适应性较强、精度较高已被广泛地 应用于地球物理的各个领域。 ［关键词］　地形　非线性共轭梯度　反演　线性 ［中图分类号］Ｐ６３１　　　［文献标识码］Ａ　［文章编号］０４９５－５３３１（２０１４）０１－０１６２－０５ ＺｈｕＣｈｏｎｇｌｉ．Ｉｎｆｌｕｅｎｃｅｏｆｔｅｒｒａｉｎｏｎｄｉｆｆｅｒｅｎｔｉｎｖｅｒｓｉｏｎｓ［Ｊ］．ＧｅｏｌｏｇｙａｎｄＥｘｐｌｏｒａｔｉｏｎ，２０１４，５０ （１）：０１６２－０１６６． ０　引言 １　非线性共轭梯度反演 目前 ＭＴ法已经应用于地质探矿各个领域，而 在二维地电模型中，给定边界条件后，一般通过 反演解释才是 ＭＴ法的终极目的。反演问题最终形 有限元法或者有限差分法进行正演模拟，反演问题 Ｔ 式都是求目标函数的最优化问题，即最小值问题，在 可表示为：ｄ＝Ｆ（ｍ）＋ｅ。式中：ｄ＝［ｄ，ｄ，…，ｄ］ １ ２ Ｎ 转化的过程中所用的手段不同，导致形成的算法也 为观测数据向量；Ｆ（ｍ）为模型正演响应函数，ｍ＝ Ｔ 不同。如当今在 ＭＴ反演中占据主导地位的牛顿 ［ｍ，ｍ，…，ｍ ］为模型参数向量；ｅ为残差向量。 １ ２ Ｎ 法、拟牛顿法、 观测数据与理论模型响应的拟合差可表示为 Ｏｃｃａｍ等均采用新点、旧点之间按照 φ ＝ ｄ 某个方向不断的更新迭代，直到收敛为止，即达到最 ２ Ｔ －１ Ｗｅ 。定义：（ｍ）＝（ｄ－Ｆ（ｍ））－Ｕ （ｄ－Ｆ ‖ ‖ φ ｄ 优目的。但是它们在根据目标函数梯度信息的线性 Ｔ Ｔ