大地电磁测深数据一维自动反演.docVIP

大地电磁测深数据一维自动反演 　　摘要：在学科实践中，反演是一种相当普遍的手段。地球物理反演则是通过对地面、地下、空间，甚至海洋上的观测（地震仪、重力仪、地电及地磁仪）资料进行分析计算，来推断地球内部介质的地震波速度、密度、电导率等参数的分布，从而得到地球内部介质分布的二维或三维结构图像。本文根据阮百尧教授在其“电阻率激发极化法测深数据的一维最优化反演方法”一文中所述，导出了在大地电磁测深中与电阻率测深一维最优反演相似的反演方法，通过自动迭代反演出地下模型参数。计算表明，这是一种快速和实用的自动反演方法。 　　Abstract: In practice, inversion is a method used widely. Geophysics inversion deduces the distribution of seismic-wave speed, density, conductivity and so on in earth interior, then obtains the 2D or 3D structure images for the distribution of earth interior medium, by analyzing the observed date on surface, in subsurface and in space, even on the sea. Referring to the correlate articles by Profession Ruan B.Y, Derived in the magnetotelluric sounding and resistivity sounding similar to the one-dimensional inversion method optimal inversion, Inverted through the automatic iterative model parameters. Calculations show that this is a fast and practical method for the automatic inversion. 　　关键词：大地电磁测深；水平层状；一维正反演；快速反演 　　Key words: magnetotelluric sounding；horizontal layered；1D forward and inversion；rapid inversion 　　 中图分类号：P631.3文献标识码：A 文章编号：1006-4311（2011）25-0285-02 　　0 引言 　　对于简单的地电模型（一维）一般能够计算其解析解。然而对于复杂的地电模型（二维或三维）则需要计算其数值解。二维或三维电阻率正演问题包括各种不同的数值计算方法，如积分方程法、有限单元法、边界单元法等。1971年，J.H.Coggon[1]首先将有限单元法引入到电法勘探正演计算中来，对于点源二维地电断面的情况，他导出了等价的变分方程，并实现了有限单元法计算。1977年，L.R.Rijo[2]改进了点源二维有限单元法，引入了一个通用性网格，使计算速度和精度大为提高，使有限单元法解决电法正演问题达到实用的阶段。Hohmann[3]在1975年首先提出激发极化和电磁场三维数值模拟方法。1978年，Pridmore[4]在博士论文中第一次应用有限单元法计算了复杂三维地电断面上的视电阻率。Dey和Morrision[5]在1979年又采用有限差分方法实现了三维数值模拟。Poirmeur和 Vasseur[6]在1988年以及Lesur[7]在1999年用边界单元法进行了三维电法数值模拟。 　　近几年来，随着野外资料采集仪器不断改进，观测仪器分辨率和抗干扰能力不断增强，电阻率测深法也从大剖面勘探发展到面积勘探，从概查应用到详查。由于计算手段的不断进步，其处理方法也由一维反演发展到目前的三维反演。为了不断改进资料处理的效果，在反演环节上追求更高的精度，要进行二维和三维反演方法的研究，同时也不能忽视对一维反演方法的更深入研究。在大地电磁测深中，Bostick基于测深曲线低频渐近线的性质提出了一种直接反演方法[8]。 　　本文根据阮百尧教授在其“电阻率激发极化法测深数据的一维最优化反演方法”[9]一文中所述，导出了在大地电磁测深中与电阻率测深一维最优反演相似的反演方法，通过自动迭代反演出地下模型参数。计算表明，这是一种快速和实用的自动反演方法。 　　1 正反演方法介绍 　　1.1 正演方法介绍 大地电磁测深的理论曲线是指在给定地下介质电阻率分布的情况下，通过计算得出的视电阻率和信号周期之