高密度电阻率法是近些年广泛应用于浅层物探的方法，隶属于直流电祖率法的范 畴，它基于圆滑约束最小二乘法的计算机反演计算程序, 使用了准牛顿最优化非线性 最小二乘新算法, 使得大数据量下的计算速度较常规最小二乘法快十倍以上, 且占 用内存较少。而大地电磁测深(MT)技术, 在寻找含油气构造、划分沉积岩相带和研究 沉积盆地基底构造中应用广泛，其资料反演须计算雅可比偏导数矩阵，对于二维模型 其反演计算量大，误差难于控制。求解大地电磁测深（MT）二维反演问题中雅克比 偏导数矩阵的逆运算量很大,并且大型反演方程组解的非唯一性非常严重，使得反演 模型常发生畸变，带来解释上的错误。而高密度电阻率法基于圆滑约束最小二乘法的 计算机反演计算程序和准牛顿最优化非线性最小二乘新算法,不仅求解唯一，而且计 算速度要快得多。 研究使用高密度电阻率反演方法对 MT 资料进行处理解释，主要解决三个问题： 第一、将非阵列数据进行插值转化为阵列数据，研究其可行性；第二、将有方向性参 数的电阻率数据用处理非方向性数据的方法来反演，研究其方法的可实现性；第三、 确定两种方法在探测深度上的转化关系。最后将两种反演模型进行对比，区分其对地 质构造反应的优劣，重点比照水平层状地层的形状和深度以及断裂带的位置。在经过 了理论的研究及实践验证之后，本文的主要成果包括：①从理论上部分证实了将大地 电磁测深资料进行高密度电阻率法反演的正确性和可行性；②从实践上验证了相对于 常规二维 RRI 反演模型来说，高密度电阻率法反演模型在确定水平地质构造的产状 及深度上有一定的效果，对局部小断层的识别灵敏度更高。③在处理尝试过程中，用 VB 语言自编了数据转化程序软件。 本文通过将ZB区块的五条MT测线进行高密度电阻率法反演尝试，并与常规大地电 磁二维反演结果进行对比，发现这种应用不仅可行，而且运算量大大降低、误差可以 控制，并且在确定高阻基底的深度和对地质构造电性特征反映上有一定的效果，在确 定目标体横向电阻率的差异和断层位置上有参考效果，可以作为大地电磁资料处理和 解释的有益补充。 关键词：高密度电阻率法 MT 二维反演 效果比较Research of the feasibility of high-density resistivity method Inversion used in Magnetotelluric（MT ）sounding data process Introduction of the author: Wang Huai-kun, male, was born in December, 1982 whose tutor was Professor Lei Wan. He graduated from Chengdu University of Technology in Earth exploration and information technology and was granted the Master Degree in June, 2009. Abstract high-density resistivity method is a part of direct-current Geophysical methods which is widely used in shallow geophysical In recent years， It is a Computer inversion and computational procedure which is based on the smooth restraint least squares method ，Because it used the accurate Newton optimization misalignment to be youngest two rides the new algorithm ，s o under the great data quantity the computation speed is ten times faster than conventional least squares method ，and takes fewer memory.Magnetotelluric sounding is widely used In the search for oil and gas structures, division of sedimentary facies sedimentary basin and basement structure that its data inversion need to be calculated Jacobian partial derivative matrix,Inversion of two-dimensional model for calculating the volume is very large, and difficult to control error. The research of high-density resistivity method Inversion used in Magnetotelluric（ MT ）sounding data process,mainly to solve the three problems:First, carries on the interpolation the non-array data to transform as the array data,researching the feasibility;Second, will have the directive parameter electronic resistivity data with to process the non-directive data the method to come the inversion ,studies the realizability of its method;Third, determines two methods in relations of the detection depth transformed.finally, comparing the two inversion model,differentiates it to the geologic structure response fit and unfit quality,focus on the shape of layered strata and depth, as well as the location of fault. This article through survey line carries on ZB the sub-area five MT the high density resistivity method inversion attempt ,and compared with the conventional two-dimensional magnetotelluric inversion results,found that this application is not only feasible,but also greatly reduce the computation and erro can easily control.It has a better effect in the determination high-resistance basis's depth,and is clearer to the geologic structure level layered characteristic reflection,it has the reference effect in the determination goal body crosswise electronic resistivity's difference and the fault position,May as the Magnetotelluric data processing explanation beneficial supplement. Keywords: high-density resistivity method MT Two-dimensional inversion The effect compares目 录 摘要 ····························································································································Ⅰ ABSTRACT···················································································································· Ⅱ 引言 ································································································································1 第 1 章 高密度电阻率法反演理论 ·············································································2 1.1 高密度电阻率法的发展历程·················································································2 1.2 基本方法原理·········································································································2 1.2.1 采集方式·········································································································2 1.2.2 视电阻率值及最佳有效探测深度··································································3 1.2.3 圆滑约束最小二乘法的基本原理··································································3 第 2 章 大地电磁测深基本理论·················································································5 2.1 测量原理·················································································································5 2.2 MT 资料常规处理和反演······················································································6 2.2.1 资料常规处理解释流程·················································································6 2.2.2 资料的预处理·································································································6 2.2.2.1 编辑与转轴·······························································································7 2.2.2.2 原始曲线 TE、TM 极化模式的确定··························································7 2.2.2.3 曲线的圆滑·······························································································7 2.2.2.4 曲线畸变分析···························································································7 2.3 资料的定性和半定量反演解释 ··············································································9 2.3.1 总纵向电导 S 参数的计算 ·············································································9 2.3.2 半定量解释方法·····························································································10 2.4 定量反演解释方法··································································································11 2.4.1 一维定量反演解释方法·················································································11 2.4.2 二维定量反演解释方法·················································································12 第 3 章 大地电磁测深进行高密度电阻率法反演问题探讨········································15 3.1 视电阻率的相似性··································································································15 3.1.1 S 等值性··········································································································15 3.1.2 T 等值性 ·········································································································16 3.2 数据类型的相似性及处理转化 ·············································································173.2.1 MT 数据的阵列插值······················································································18 3.2.2 大地电磁的方向性·························································································20 3.2.3 三种模式的共同点·························································································23 第 4 章 ZB 区块 MT 资料高密度反演尝试 ·································································25 4.1 测区地球物理概况··································································································25 4.1.1 测线长度及合格率·························································································25 4.1.2 自然地理概况·································································································25 4.1.3 区域构造概况·································································································26 4.2 MT 资料进行圆滑约束二乘反演的具体实施·························································29 4.2.1 处理流程·········································································································29 4.2.2 数据编辑·········································································································30 4.2.3 MT01 实测资料的处理与解释 ······································································30 4.2.4 电磁数据的圆滑约束二乘反演 ·····································································35 4.3 反演解释及推断······································································································37 4.4 小结 ·························································································································41 结 论 ····························································································································42 致 谢 ····························································································································44