其中＝ f 被称为灵敏度矩阵元素为讨论方 - 地质与勘探.pdfVIP

维普资讯 第42卷 第5期 地 质 与 勘 探 Vo1．42 No．5 2007年9月 GE0L0GYAND PR0SPECTING September，2007 Occam反演及其在瞬变 电磁测深 中_的应用 翁爱华 (吉林大学 应用地球物理系，长春 130026) [摘 要]Occam反演技术被应用中心回线瞬变电磁测深数据的反演中。理论模型和实际数据的反 演结果表明，利用Occam反演技术可以较为准确地获得地电断面的电阻率分布。一般5次迭代 目标函 数就可以收敛到5×10～。 [关键词】瞬变电磁测深 中心回线 Occam反演 GeoElectro软件系统 [中图分类号]P631．3；P628 [文献标识码】A [文章编号]0495—5331(2007)05—0074—03 {d，i=1，…，Ⅳ}，它可表示为相应的模型理论响应 0 引 言 与观测误差的和，即 ． 瞬变电磁测深在石油勘探、矿产勘查、甚至水文 dj= Q，n)+ ，i=1，…Ⅳ． (1) 调查工程勘察的广泛应用得益于一系列优点…。 其中 为方差为or的0均值数据独立随机噪音 ，n 但由于瞬变电磁技术理论模拟的复杂性，对其观测 为长度为 的未知的层状电阻率分布矢量。反演 数据 目前仍然主要进行定性分析和一维定量半定量 时，要求模型电阻率使下面的目标函数达到最小，即 解释 。在对模型进行一维定量反演时 ，通常从 x2(nest，d)：∑N( 1 (2) 模型的某个猜测开始，采用迭代技术 ，对模型参数进 l 、 (J ， 行逐步的修正，最终得到反演结果_5j。 由于感应电动势是电阻率的非线性函数，对于上面 尽管一般的迭代过程经过修正，能减轻反演对 的极小化问题可从一个初始猜测开始，通过迭代加 初始模型的依赖性，但在实际工作中，摆脱对初始模 以解决 。假设初始模型为n。接近真解，在该猜测 型的依赖 仍然是非常有意义 的问题。Constable 的一个小的邻域 内，将模型响应展开，并忽略高次 等 和DeGroot—Hedlin等将传统的非线性反演 项，有 迭代过程作了改进，提出了著名的Occam反演方 法，实现了反演过程的自动化，并将其应用到大地电 di= Qi，no)+∑ △ ，i=I，…N，(3) 磁测深、电阻率测深的一维、二维模型反演中。由于 一 J=1 Occam反演的反演结果独立于初始模型猜测，国内 其中=f。被称为灵敏度矩阵元素。为讨论方 也将其应用到大地电磁测深资料的反演中 j。 便，假设数据方差为 1。将上式代入(1)式并对模型 本文将Occam反演方法应用到瞬变 电磁测深