大地电磁三维畸变特征及校正新技术的应用研究
|
摘要
大地电磁研究和应用中,面临的实际地质结构多数是三维的,但是当前三维反演技术还不成熟,仍需要利用二维模型进行近似反演解释。在进行二维反演解释时,如何分析认识三维局部畸变的性质,最大限度地克服三维局部畸变的影响,获得尽可能可信的地质解释结果,如何认识三维模型的二维近似程度,如何利用观测的TE模式和TM模式的视电阻率和相位资料等问题,是受到普遍关注的科学问题和实际难题,也成为当前影响大地电磁向前发展的关键问题之一。
本论文正是针对上面存在的问题,采用三维数值模拟方法模拟大地电磁局部畸变,在此基础上,研究局部畸变的影响特征,提出新的畸变校正方法和新的电性结构维性判定参数,并对三维数据的二维近似反演从数据差异性和反演结果合理性的角度做了详细的研究,最后通过对实测资料的应用检验本文的研究成果。论文主要研究内容和取得的成果概述如下:
1.大地电磁局部畸变基础理论及研究状况归纳总结
从大地电磁局部畸变的电磁场控制方程出发,推导了异常电场和异常磁场的积分表达式。论述了局部畸变模型成立的物理假设,并详细推导了阻抗张量畸变和倾子矢量畸变的计算表达式。介绍了局部畸变的传统分类中电流型畸变和感应型畸变的概念,它们对二维TE和TM模式视电阻率曲线的影响,最后简单分析了近地表三维局部畸变影响的一些特点。介绍了大地电磁构造维性分析中常用的阻抗张量基本旋转不变量、WAL旋转不变量以及基于这些不变量的构造维性判定方法、B-Q构造维性判定参数,指出应用这些不变量分析时应注意的问题。详细介绍了目前已有的常用畸变校正方法,对与本论文关系比较密切的算法如Swift、Bahr、GB方法做了详细的推导说明,并指出其优点和不足之处,为本文提出新的畸变校正方法提供了研究的出发点。另外,对电、磁场全畸变公式中的错误也作了修正。
2.大地电磁局部畸变特征的三维模拟分析
在第一部分理论研究的基础上,设计典型的局部畸变模型,采用改进的交错网格有限差分法进行了三维局部畸变的数值模拟,对计算结果进行了定性的分析,获得一些新的认识。研究发现TM模式视电阻率曲线在穿越低阻异常体界面时,曲线会先上移然后再下移,而在穿越高阻异常体界面时,曲线会先下移再上移。这与电性分界面处积累面电荷产生的二次电场有关。当异常体为高阻时,在其表面积累的面电荷产生的二次电场,在异常体上方其与一次电场方向相同,则观测总场是增大的。在异常体的两端边界附近,二次电场与一次电场的方向是相反的,因此导致观测的总场是减小的;当异常体为良导时,在异常体上方二次电场与一次电场的方向相反,则观测总场是减小的。在异常体的两端边界附近,二次电场与一次电场是同向的,观测总场是增大的。
区域二维构造走向为南北方向时,在YX模式视电阻率平面图中三维低阻异常体在南北方向被拉伸,在东西边界处视电阻率分别存在一个增大的区域。而XY模式视电阻率平面图中三维低阻体异常体在东西方向被拉伸,在南北边界处存在视电阻率增大区域。在XY模式相位平面图中,三维低阻异常体在东西方向被拉长,在南北边界处存在相位值减小区。在YX模式相位平面图中,三维低阻异常体在南北方向被拉长,在东西边界处存在相位值减小区。该现象与二维TM模式曲线在穿越电性分界面时的规律是一致的,当测线沿东西方向展布时,此时的三维响应的YX模式对应二维情况下的TM模式。而当测线沿南北方向展布时,此时三维响应XY模式对应二维情况下的TM模式。
低阻小异常体对二维区域响应的畸变比高阻小异常体要严重的多,高阻体对区域响应TM模式的畸变要稍大于TE模式,而低阻体则是对TE模式的畸变大于TM模式。同一极化模式,相位曲线相比视电阻率曲线而言,受畸变的影响要小。
位于小异常体中心上方测点的三维畸变响应与对应二维区域响应的差异比较大,但可以等效于其他的二维模型响应,这种特殊畸变在资料解释时要特别注意。
3.大地电磁局部畸变校正新方法研究
介绍了国外研究者新近提出的得到广泛应用的相位张量,并在国内首次将该方法应用于理论数据和实测资料的解释。提出了一种新的单点单频畸变校正方法(SSSF)和单点多频校正方法(SSMF),推导了其理论计算公式,在此研究的基础上提出新的构造维性判定参数,编制了相应的局部畸变分析软件。首先采用合成理论数据验证了新方法的正确性,然后采用Swift、Bahr、GB、相位张量、WAL、SSSF和SSMF方法对三维局部畸变模拟的结果进行了定量的分析,最后将新方法应用于实测资料的解释。研究结果表明SSSF和SSMF法的计算结果与其他张量分解方法基本一致或更好,其提供的构造维性参数和畸变因子为大地电磁资料的分析提供了有力的补充。另外根据实际应用中存在的问题,修正了WAL方法维性判定准则。
4.三维畸变数据的二维反演中极化模式的选择
针对三维和二维模型响应数据的差异、如何认识二维模式的近似程度,如何利用观测的TE模式和TM模式的视电阻率和相位资料等问题,对三维模型数据的二维近似反演问题的研究,主要从三维模型和二维模型响应差异性和二维反演结果合理性两个方面进行了研究。研究结果表明:在三维模型条件下,利用二维模型进行反演时,TE模式对模型的二维的近似程度要求远远高于TM模式;当三维结构影响较明显时,仅利用TM模式数据进行二维反演比仅利用TE模式或利用TE+TM模式联合反演都更合理,反演结果中的虚假结构明显减少;对于TM模式,相位受三维畸变影响较小,视电阻率较大,所以反演中可适当加大相位的权;对于实测数据的二维反演,应优先考虑采用TM模式数据进行二维反演,其次是TM+TE模式,一般不要单独采用TE模式。本章对上述情况提出了物理分析的依据。
5.论文研究成果在实测资料中的综合应用
将论文研究成果综合应用于青藏高原东边缘带的美姑-绥江剖面大地电磁实测资料的解释中,对区域构造维性采用相位张量、SSSF方法和WAL方法进行综合判定,采用SSSF方法确定区域主轴方位角,局部畸变校正选用了论文新提出的SSMF方法。二维反演采用TM模式数据进行反演,在反演中还采用了多项新技术,其一是利用测点中心网格的自动生成技术构建带地形的二维反演初始模型,其二是正则化因子τ的分段取值。采用复杂模型大地电磁有效探测深度技术计算了二维反演结果的有效趋肤深度,同时参考相位张量不变量参数的结果,综合评价了反演结果的可靠性。取得了比较满意的反演结果,为研究区地电模型的建立和地质解释奠定了坚实的基础。
In MT surveys, the actual structure of the Earth is often three-dimensional (3D). However the current state-of-the-art for MT data interpretation is two-dimensional (2D) inversion and 3D inversion will become available in the near future. When interpretation of 3D data with 2-D techniques, some problems are still exist, such as how to analyze the features of 3D galvanic distortion, how to overcome the influences of 3D local distortion as much as possible and get more reliable geological results, how to understand 2D approximate interpretation of 3D models, and how to make best use of apparent resistivity and phase of TE mode and TM mode. These are the focus of concern scientific problems and the key factors which delay the developing of MT. This thesis focuses on the problems in the above paragraph. MT responses of local distortion are computed using a 3D finite-difference on a staggered grid code. The effect of local distortion has been studied. New correction algorithm and new parameters for determination of geoelectric dimensionality are brought out. The study about 2D approximate inversion versus 3D data is done. The research is divided into two parts: one is the difference between 3D and 2D responses; another is validity of 2D inversion results. At last the research results of this thesis are applied to the interpretation of filed data. The conclusions of this thesis are summarized below.
1. Basic theory and summary of research situation of galvanic distortion in MT
Based on the governing equation for the EM field of local distortion, the integral expression of anomalous electric and magnetic field are deduced. The physics hypothesis for local distortion model is discussed and the expression of impedance distortion tensor and vertical magnetic distortion vector are derived from integral equation in detail. Galvanic distortion and induction distortion in classic classification of local distortion and their effect on apparent resistivity and phase curves of TE mode and TM mode are illustrated. Some features of local distortion caused by small-scale surface inhomogeneous are briefly discussed. Basic invariants of impedance tensor and WAL rotational invariants used in geoelectric dimensionality analysis are introduced. Some methods and B-Q parameters for determination of geoelectric dimensionality based on these invariants are presented, and some issues needing attention using these invariants are indicated. Several approaches close with this thesis such as Swift, Bahr, GB methods are reviewed in detail. The merits and shortcoming of these methods are point out. These work provide a research beginning for developing novel distortion correction algorithm in this thesis. The bug in formula describing distortion of electric and magnetic fields is fixed.
2. Three-dimensional modeling and analysis of local distortion in MT
Based on theory in the first part, several 3D local distortion models are designed and model responses are computed using a improved finite-difference on a staggered grid code. Qualitative analysis for model responses is done and some new insights are obtained. When apparent resistivity curve of TM mode across small conductive surface body, the curve will upward shift at first and then downward shift. While across resistive body, the curve will downward shift at first and then upward shift.. This phenomenon is associated with secondary electric field caused by the accumulative surface electric charges at the boundaries of the body. When small body is resistive or conductive, the polarity of the charges results in a secondary field additive to the primary field or opposing the primary fields. Upon vectorial addition of the primary field and secondary field, one can estimate the sense of distortion in the resulting total electric fields. In the resistive body, the total field is increased directly over the body, it is decreased off the ends, and an increase occurs along the sides of the body. With a conductive body, the total field is reduced directly over the body, it is enhanced off the ends of the body, and it is diminished along the sides.
When the strike direction of 2D regional structure is NS direction, in the plane map of apparent resistivity of YX mode, 3D conductive body is elongated in NS direction, two zones with the values increasing are located respectively at east and west boundary. While in the plane map of apparent resistivity of XY mode, 3D conductive body is elongated in EW direction, two zones with the values increasing are located respectively at north and south boundary. For phase of XY mode, 3D conductive body is elongated in EW direction; two zones with the values decreasing are located respectively at north and south boundary. For phase of YX mode, 3D conductive body is elongated in NS direction; two zones with the values decreasing are located respectively at east and west boundary. This phenomenon is consistent with the case that curves of 2D TM mode across the electric boundaries. When the profile runs along EW direction, YX mode response in 3D corresponds to TM mode response in 2D. If the profile runs along NS direction, XY mode response in 3D corresponds to TM mode response in 2D.
The distortion on 2D regional response caused by small conductive inhomogeneous are more severe than that of resistive inhomogeneous. The distortion on TM mode response is slightly larger than that on TE mode caused by resistive body. While for conductive body, the distortion on TE mode response is larger than that of TM mode. In the same mode the distortion on phase is less than that on apparent resistivity.
At the position located directly over the small inhomogeneous, the difference between 3D distortion response and 2D corresponding regional response is much large. But this 3D distortion response can be equivalent to other 2D response. It is more important to pay more attention to this special distortion in MT interpretation.
3. Research of new correction method for local distortion in MT
Phase tensor approach put forward by overseas geophysicist is introduced and it is the first time applied to the interpretation of synthetic data and field data at home. New single site, single frequency correction algorithm (SSSF) and single site, multi-frequency algorithm (SSMF) are brought out. The formula for them are present in detail. Based on these research, new parameters for determination of geoelectric dimensionality are defined and a program to perform MT distortion analysis is developed. Firstly, synthetic data is used to check the validity of new methods. Secondly, quantitative analysis on 3D local distortion response is done using Swift, Bahr, GB, phase tensor, WAL, SSSF and SSMF approach. Finally, new methods are applied to field data. The research results manifest that the application of new methods are comparable to or even superior to other decomposition methods. The parameters for determination of geoelectric dimensionality and distortion factors provided by new methods are important complement to the analysis of MT data. According to problem occurred in application for WAL method, WAL invariants criteria is corrected.
4. The choice of polarization mode in 2D inversion versus 3D MT distortion data
In order to understand the difference between 3D response and 2D response, answer the questions how to understand 2D approximate interpretation of 3D models, and how to make best use of apparent resistivity and phase of TE mode and TM mode. the study about 2D approximate inversion versus 3D data is done. The research is divided into two parts: one is the difference between 3D response and 2D response, another is validity of 2D inversion results.
The results manifest that in the 3D scenario the demand for approximate 2D structure in TE mode is much stricter than in TM mode. When 3D effect is relatively severe, by applying 2D inversion to 3D data, the results show that the inversion model using only TM mode data is more validity, and the phantom structures are less than using TE mode or TE+TM mode. In TM mode, the distortion on phase is less than apparent resistivity, so the weight of phase in inversion can be increased. For 2D inversion of field data, TM mode is first choice, next is TE+TM, TE mode is not proposed. The physics analysis basis is provided in this chapter.
5. Application of research results in filed data
The research results of this thesis are applied to the interpretation of Megu-Suijiang field data in eastern margin of Tibetan plateau. The dimensionality analysis is carried out by phase tensor, SSSF method and WAL invariants criteria. The strike direction of regional structure is determined by SSSF approach. SSMF method is used in local distortion correction. TM mode data is used in 2D inversion. Some new approaches are adopted in inversion, one is Site-Center Mesh (SCM) which used in constructing initial model, another is the regularization parametersτwhich determined through segmental variation. The validity of inversion model is checked by using effective skin depth of complex structure method and invariants of phase tensor. A good inversion results is obtained and it has laid a solid foundation for the construction of electric structure and geological interpretation.
本论文正是针对上面存在的问题,采用三维数值模拟方法模拟大地电磁局部畸变,在此基础上,研究局部畸变的影响特征,提出新的畸变校正方法和新的电性结构维性判定参数,并对三维数据的二维近似反演从数据差异性和反演结果合理性的角度做了详细的研究,最后通过对实测资料的应用检验本文的研究成果。论文主要研究内容和取得的成果概述如下:
1.大地电磁局部畸变基础理论及研究状况归纳总结
从大地电磁局部畸变的电磁场控制方程出发,推导了异常电场和异常磁场的积分表达式。论述了局部畸变模型成立的物理假设,并详细推导了阻抗张量畸变和倾子矢量畸变的计算表达式。介绍了局部畸变的传统分类中电流型畸变和感应型畸变的概念,它们对二维TE和TM模式视电阻率曲线的影响,最后简单分析了近地表三维局部畸变影响的一些特点。介绍了大地电磁构造维性分析中常用的阻抗张量基本旋转不变量、WAL旋转不变量以及基于这些不变量的构造维性判定方法、B-Q构造维性判定参数,指出应用这些不变量分析时应注意的问题。详细介绍了目前已有的常用畸变校正方法,对与本论文关系比较密切的算法如Swift、Bahr、GB方法做了详细的推导说明,并指出其优点和不足之处,为本文提出新的畸变校正方法提供了研究的出发点。另外,对电、磁场全畸变公式中的错误也作了修正。
2.大地电磁局部畸变特征的三维模拟分析
在第一部分理论研究的基础上,设计典型的局部畸变模型,采用改进的交错网格有限差分法进行了三维局部畸变的数值模拟,对计算结果进行了定性的分析,获得一些新的认识。研究发现TM模式视电阻率曲线在穿越低阻异常体界面时,曲线会先上移然后再下移,而在穿越高阻异常体界面时,曲线会先下移再上移。这与电性分界面处积累面电荷产生的二次电场有关。当异常体为高阻时,在其表面积累的面电荷产生的二次电场,在异常体上方其与一次电场方向相同,则观测总场是增大的。在异常体的两端边界附近,二次电场与一次电场的方向是相反的,因此导致观测的总场是减小的;当异常体为良导时,在异常体上方二次电场与一次电场的方向相反,则观测总场是减小的。在异常体的两端边界附近,二次电场与一次电场是同向的,观测总场是增大的。
区域二维构造走向为南北方向时,在YX模式视电阻率平面图中三维低阻异常体在南北方向被拉伸,在东西边界处视电阻率分别存在一个增大的区域。而XY模式视电阻率平面图中三维低阻体异常体在东西方向被拉伸,在南北边界处存在视电阻率增大区域。在XY模式相位平面图中,三维低阻异常体在东西方向被拉长,在南北边界处存在相位值减小区。在YX模式相位平面图中,三维低阻异常体在南北方向被拉长,在东西边界处存在相位值减小区。该现象与二维TM模式曲线在穿越电性分界面时的规律是一致的,当测线沿东西方向展布时,此时的三维响应的YX模式对应二维情况下的TM模式。而当测线沿南北方向展布时,此时三维响应XY模式对应二维情况下的TM模式。
低阻小异常体对二维区域响应的畸变比高阻小异常体要严重的多,高阻体对区域响应TM模式的畸变要稍大于TE模式,而低阻体则是对TE模式的畸变大于TM模式。同一极化模式,相位曲线相比视电阻率曲线而言,受畸变的影响要小。
位于小异常体中心上方测点的三维畸变响应与对应二维区域响应的差异比较大,但可以等效于其他的二维模型响应,这种特殊畸变在资料解释时要特别注意。
3.大地电磁局部畸变校正新方法研究
介绍了国外研究者新近提出的得到广泛应用的相位张量,并在国内首次将该方法应用于理论数据和实测资料的解释。提出了一种新的单点单频畸变校正方法(SSSF)和单点多频校正方法(SSMF),推导了其理论计算公式,在此研究的基础上提出新的构造维性判定参数,编制了相应的局部畸变分析软件。首先采用合成理论数据验证了新方法的正确性,然后采用Swift、Bahr、GB、相位张量、WAL、SSSF和SSMF方法对三维局部畸变模拟的结果进行了定量的分析,最后将新方法应用于实测资料的解释。研究结果表明SSSF和SSMF法的计算结果与其他张量分解方法基本一致或更好,其提供的构造维性参数和畸变因子为大地电磁资料的分析提供了有力的补充。另外根据实际应用中存在的问题,修正了WAL方法维性判定准则。
4.三维畸变数据的二维反演中极化模式的选择
针对三维和二维模型响应数据的差异、如何认识二维模式的近似程度,如何利用观测的TE模式和TM模式的视电阻率和相位资料等问题,对三维模型数据的二维近似反演问题的研究,主要从三维模型和二维模型响应差异性和二维反演结果合理性两个方面进行了研究。研究结果表明:在三维模型条件下,利用二维模型进行反演时,TE模式对模型的二维的近似程度要求远远高于TM模式;当三维结构影响较明显时,仅利用TM模式数据进行二维反演比仅利用TE模式或利用TE+TM模式联合反演都更合理,反演结果中的虚假结构明显减少;对于TM模式,相位受三维畸变影响较小,视电阻率较大,所以反演中可适当加大相位的权;对于实测数据的二维反演,应优先考虑采用TM模式数据进行二维反演,其次是TM+TE模式,一般不要单独采用TE模式。本章对上述情况提出了物理分析的依据。
5.论文研究成果在实测资料中的综合应用
将论文研究成果综合应用于青藏高原东边缘带的美姑-绥江剖面大地电磁实测资料的解释中,对区域构造维性采用相位张量、SSSF方法和WAL方法进行综合判定,采用SSSF方法确定区域主轴方位角,局部畸变校正选用了论文新提出的SSMF方法。二维反演采用TM模式数据进行反演,在反演中还采用了多项新技术,其一是利用测点中心网格的自动生成技术构建带地形的二维反演初始模型,其二是正则化因子τ的分段取值。采用复杂模型大地电磁有效探测深度技术计算了二维反演结果的有效趋肤深度,同时参考相位张量不变量参数的结果,综合评价了反演结果的可靠性。取得了比较满意的反演结果,为研究区地电模型的建立和地质解释奠定了坚实的基础。
In MT surveys, the actual structure of the Earth is often three-dimensional (3D). However the current state-of-the-art for MT data interpretation is two-dimensional (2D) inversion and 3D inversion will become available in the near future. When interpretation of 3D data with 2-D techniques, some problems are still exist, such as how to analyze the features of 3D galvanic distortion, how to overcome the influences of 3D local distortion as much as possible and get more reliable geological results, how to understand 2D approximate interpretation of 3D models, and how to make best use of apparent resistivity and phase of TE mode and TM mode. These are the focus of concern scientific problems and the key factors which delay the developing of MT. This thesis focuses on the problems in the above paragraph. MT responses of local distortion are computed using a 3D finite-difference on a staggered grid code. The effect of local distortion has been studied. New correction algorithm and new parameters for determination of geoelectric dimensionality are brought out. The study about 2D approximate inversion versus 3D data is done. The research is divided into two parts: one is the difference between 3D and 2D responses; another is validity of 2D inversion results. At last the research results of this thesis are applied to the interpretation of filed data. The conclusions of this thesis are summarized below.
1. Basic theory and summary of research situation of galvanic distortion in MT
Based on the governing equation for the EM field of local distortion, the integral expression of anomalous electric and magnetic field are deduced. The physics hypothesis for local distortion model is discussed and the expression of impedance distortion tensor and vertical magnetic distortion vector are derived from integral equation in detail. Galvanic distortion and induction distortion in classic classification of local distortion and their effect on apparent resistivity and phase curves of TE mode and TM mode are illustrated. Some features of local distortion caused by small-scale surface inhomogeneous are briefly discussed. Basic invariants of impedance tensor and WAL rotational invariants used in geoelectric dimensionality analysis are introduced. Some methods and B-Q parameters for determination of geoelectric dimensionality based on these invariants are presented, and some issues needing attention using these invariants are indicated. Several approaches close with this thesis such as Swift, Bahr, GB methods are reviewed in detail. The merits and shortcoming of these methods are point out. These work provide a research beginning for developing novel distortion correction algorithm in this thesis. The bug in formula describing distortion of electric and magnetic fields is fixed.
2. Three-dimensional modeling and analysis of local distortion in MT
Based on theory in the first part, several 3D local distortion models are designed and model responses are computed using a improved finite-difference on a staggered grid code. Qualitative analysis for model responses is done and some new insights are obtained. When apparent resistivity curve of TM mode across small conductive surface body, the curve will upward shift at first and then downward shift. While across resistive body, the curve will downward shift at first and then upward shift.. This phenomenon is associated with secondary electric field caused by the accumulative surface electric charges at the boundaries of the body. When small body is resistive or conductive, the polarity of the charges results in a secondary field additive to the primary field or opposing the primary fields. Upon vectorial addition of the primary field and secondary field, one can estimate the sense of distortion in the resulting total electric fields. In the resistive body, the total field is increased directly over the body, it is decreased off the ends, and an increase occurs along the sides of the body. With a conductive body, the total field is reduced directly over the body, it is enhanced off the ends of the body, and it is diminished along the sides.
When the strike direction of 2D regional structure is NS direction, in the plane map of apparent resistivity of YX mode, 3D conductive body is elongated in NS direction, two zones with the values increasing are located respectively at east and west boundary. While in the plane map of apparent resistivity of XY mode, 3D conductive body is elongated in EW direction, two zones with the values increasing are located respectively at north and south boundary. For phase of XY mode, 3D conductive body is elongated in EW direction; two zones with the values decreasing are located respectively at north and south boundary. For phase of YX mode, 3D conductive body is elongated in NS direction; two zones with the values decreasing are located respectively at east and west boundary. This phenomenon is consistent with the case that curves of 2D TM mode across the electric boundaries. When the profile runs along EW direction, YX mode response in 3D corresponds to TM mode response in 2D. If the profile runs along NS direction, XY mode response in 3D corresponds to TM mode response in 2D.
The distortion on 2D regional response caused by small conductive inhomogeneous are more severe than that of resistive inhomogeneous. The distortion on TM mode response is slightly larger than that on TE mode caused by resistive body. While for conductive body, the distortion on TE mode response is larger than that of TM mode. In the same mode the distortion on phase is less than that on apparent resistivity.
At the position located directly over the small inhomogeneous, the difference between 3D distortion response and 2D corresponding regional response is much large. But this 3D distortion response can be equivalent to other 2D response. It is more important to pay more attention to this special distortion in MT interpretation.
3. Research of new correction method for local distortion in MT
Phase tensor approach put forward by overseas geophysicist is introduced and it is the first time applied to the interpretation of synthetic data and field data at home. New single site, single frequency correction algorithm (SSSF) and single site, multi-frequency algorithm (SSMF) are brought out. The formula for them are present in detail. Based on these research, new parameters for determination of geoelectric dimensionality are defined and a program to perform MT distortion analysis is developed. Firstly, synthetic data is used to check the validity of new methods. Secondly, quantitative analysis on 3D local distortion response is done using Swift, Bahr, GB, phase tensor, WAL, SSSF and SSMF approach. Finally, new methods are applied to field data. The research results manifest that the application of new methods are comparable to or even superior to other decomposition methods. The parameters for determination of geoelectric dimensionality and distortion factors provided by new methods are important complement to the analysis of MT data. According to problem occurred in application for WAL method, WAL invariants criteria is corrected.
4. The choice of polarization mode in 2D inversion versus 3D MT distortion data
In order to understand the difference between 3D response and 2D response, answer the questions how to understand 2D approximate interpretation of 3D models, and how to make best use of apparent resistivity and phase of TE mode and TM mode. the study about 2D approximate inversion versus 3D data is done. The research is divided into two parts: one is the difference between 3D response and 2D response, another is validity of 2D inversion results.
The results manifest that in the 3D scenario the demand for approximate 2D structure in TE mode is much stricter than in TM mode. When 3D effect is relatively severe, by applying 2D inversion to 3D data, the results show that the inversion model using only TM mode data is more validity, and the phantom structures are less than using TE mode or TE+TM mode. In TM mode, the distortion on phase is less than apparent resistivity, so the weight of phase in inversion can be increased. For 2D inversion of field data, TM mode is first choice, next is TE+TM, TE mode is not proposed. The physics analysis basis is provided in this chapter.
5. Application of research results in filed data
The research results of this thesis are applied to the interpretation of Megu-Suijiang field data in eastern margin of Tibetan plateau. The dimensionality analysis is carried out by phase tensor, SSSF method and WAL invariants criteria. The strike direction of regional structure is determined by SSSF approach. SSMF method is used in local distortion correction. TM mode data is used in 2D inversion. Some new approaches are adopted in inversion, one is Site-Center Mesh (SCM) which used in constructing initial model, another is the regularization parametersτwhich determined through segmental variation. The validity of inversion model is checked by using effective skin depth of complex structure method and invariants of phase tensor. A good inversion results is obtained and it has laid a solid foundation for the construction of electric structure and geological interpretation.
引文
[1] Tikhonov A N. On the determination of electric characteristics of deep layers of the earth’s crust[J]. Dok. Akad. Nauk. 1950, 73(2): 295-297.
[2] Cagniard L. Basic theory of the magneto-telluric method of geophysicial prospecting[J]. Geophysics. 1953, 18(3): 605-635.
[3]刘国栋,陈乐寿.大地电磁测深研究[M].北京:地震出版社, 1984.
[4]陈乐寿,王光鄂.大地电磁测深法[M].北京:地质出版社, 1990.
[5] Madden T, Nelson P. A defense of Cagniard's magnetotelluric method[R]. Massachusetts Institute Of Technology, 1964.
[6] Berdichevsky M N, Dmitriev V I. Basic principles of interpretation of magnetotelluric sounding curves[C]. Geoelectric and Geothermal Studies, KAPG Geophysical Monograph., Akademiai Kiado, Budapest, 1976.
[7] Morrison H F, Phillips R J, Brien D P. Quantitative interpretation of transient electromagnetic fields over a layered half-space[J]. Geophysical Prospecting. 1969, 17(1): 82-101.
[8] Cantwell T, Madden T R. Preliminary report on crustal magnetotelluric measurements[J]. Journal of Geophysical Research. 1960, 65: 4202-4205.
[9] Bostick F X, Smith H W. Investigation of large-scale inhomogenities in the earth by the magnetotelluric method[J]. Proceeding of the IRE(IEEE). 1962, 50: 2339-2346.
[10] Swift C M. A magnetotelluric investigation of an electrical conductivity anomaly in the southwestern United States[D]. Massachusetts Institute of Technology, 1967.
[11] Sims W E, Bostick F X, Smith H W. The estimation of magnetotelluric impedance tensor elements from measured data[J]. Geophysics. 1971, 36(5): 938-942.
[12] Gamble T D. Remote reference magnetotellurics with SQUIDS[D]. Univ. of California, Berkeley, 1978.
[13] Gamble T D, Goubau W M, Clarke J. Error analysis for remote reference magnetotellurics[J]. Geophysics. 1979a, 44(5): 959-968.
[14] Gamble T D, Goubau W M, Clarke J. Magnetotellurics with a remote magnetic reference[J]. Geophysics. 1979b, 44(1): 53-68.
[15] Sutarno D, Vozoff K. Phase-smoothed robust M-estimation of magnetotelluric impedance functions[J]. Geophysics. 1991, 56(12): 1999-2007.
[16] Egbert G D, Livelybrooks D W. Single station magnetotelluric impedance estimation: Coherence weighting and the regression M-estimate[J]. Geophysics. 1996, 61(4): 964-970.
[17] Larsen J C, Mackie R L, Manzella A, et al. Robust smooth magneto tell uric transfer functions[J]. Geophysical Journal International. 1996, 124(3): 801-819.
[18] Oettinger G, Haak V, Larsen J C. Noise reduction in magnetotelluric time-series with a new signal-noise separation method and its application to a field experiment in the Saxonian Granulite Massif[J]. Geophysical Journal International. 2001, 146: 659-669.
[19] Smirnov M Y. Magnetotelluric data processing with a robust statistical procedure having a high breakdown point[J]. Geophysical Journal International. 2003, 152(1): 1-7.
[20] Chave A D, Thomson D J. Robust, controlled leverage estimation of magnetotelluric responsefunctions[C]. 11th Workshop on EM Induction in the Earth, Wellington, New Zealand, 1992.
[21] Chave A D, Thomson D J. Some comments on magnetotelluric response function estimation[J]. Journal of Geophysical research. 1989, 94(B10): 215-225.
[22] Chave A D, Thomson D J, Ander M E. On the Robust Estimation of Power Spectra, Coherences, and Transfer Functions[J]. J. Geophys. Res. 1987, 92(B1): 633-648.
[23] Egbert G D, Booker J R. Robust estimation of geomagnetic transfer functions[J]. Geophysical Journal International. 1986, 87(1): 173-194.
[24] Chant I J, Hastie L M. Time-frequency analysis of magnetotelluric data[J]. Geophysical Journal International. 1992, 111(2): 399-413.
[25] Raiche A P. An integral equation approach to 3-D modeling[J]. Geoph. J. R. Astr. Soc. 1974, 55: 507-511.
[26] Hohmann G W. Three-dimensional induced polarization and electromagnetic modeling[J]. Geophysics. 1975, 40(2): 309-324.
[27] Ting S C, Hohmann G W. Integral equation modeling of three-dimensional magnetotelluric response[J]. Geophysics. 1981, 46(2): 182-197.
[28] Wannamaker P E. Advances in three-dimensional magnetotelluric modeling using integral equations[J]. Geophysics. 1991, 56(11): 1716-1728.
[29] Wannamaker P E, Hohmann G W, Ward S H. Magnetotelluric responses of three-dimensional bodies in layered earths[J]. Geophysics. 1984, 49(9): 1517-1533.
[30] Zhdanov M S, Fang S. Quasi-linear approximation in 3-D electromagnetic modeling[J]. Geophysics. 1996, 61(3): 646-665.
[31] Xiong Z. Electromagnetic modeling of 3-D structures by the method of system iteration using integral equations[J]. Geophysics. 1992, 57(12): 1556-1561.
[32] Xiong Z, Tripp A C. Electromagnetic scattering of large structures in layered earths using integral equations[J]. Radio Sci. 1995, 30(4): 921-929.
[33] Jones F W, Vozoff K. The calculation of magnetotelluric quantities for three-dimensional conductivity inhomogeneities[J]. Geophysics. 1978, 43(6): 1167-1175.
[34] Lines L R, Jones F W. The perturbation of alternating geomagnetic fields by three-dimensional island structure[J]. Geophys. J. Roy. Astr. Soc. 1973, 32: 133-154.
[35] Pu X, Agarwal A K, Weaver J T. Magnetic field solutions of E-polarization induction problems[J]. J. Geomagn. Geoelectr. 1993, 45: 859-872.
[36]陈小斌,汤吉,赵国泽. MT三维数值模拟中辅助场的求解[C].中国地球物理学会年会,中国昆明2001.
[37]汤吉.长白山天池火山区地壳深部电性结构和大地电磁三维数值模拟研究[D].中国地震局地质研究所, 2001.
[38]陈乐寿.有限元法在大地电磁场正演计算中的应用及改进[J].地球科学-中国地质大学学报. 1981.
[39]陈小斌.有限元直接迭代算法及其在电磁模拟中的应用研究[D].江汉石油学院, 2000a.
[40]李予国,徐世浙,刘斌,等.电导率分块连续变化的二维MT有限元模拟(II)[J].高校地质学报. 1996, 2(4): 448-452.
[41]史明娟,徐世浙,刘斌.大地电磁二次函数插值的有限元法正演模拟[J].地球物理学报. 1997, 40(3): 421-430.
[42]徐世浙,刘斌.电导率分层连续变化的水平层的大地电磁正演[J].地球物理学报. 1995a, 38(2): 262-268.
[43]徐世浙,于涛,李予国,等.电导率分块连续变化的二维MT有限元模拟(I)[J].高校地质学报. 1995b, 1(2): 65-73.
[44]赵生凯,徐世浙.有限单元法计算良导嵌入体的大地电磁测深曲线[J].物探化探计算技术. 1983, 5(1): 47-55.
[45]胡建德,蔡纲.用三角形二次插值法计算二维大地电磁测深曲线[J].石油地球物理勘探. 1984, 19(4): 358-367.
[46]罗焕炎,张淑霞.二维大地电磁测深的边界元与有限元耦合模拟[J].地震地质. 1987, 9(2): 1-9.
[47]王宝钧,高维安.用有限元正演京津塘地区二维构造的大地电磁测深资料[J].地震地质. 1984, 6(1): 69-80.
[48]陈乐寿,孙必俊.有限元法在大地电磁测深中正演计算的改进[J].石油地球物理勘探. 1982, 17(3): 69-72.
[49] Wannamaker P E, Stodt J A, Rijo L. Two-dimensional topographic responses in magnetotellurics modeled using finite elements[J]. Geophysics. 1986, 51(11): 2131-2144.
[50] Silvester P, Hastam C R S. Magnetotelluirc modeling by the finite element method[J]. Geophysical Prospecting. 1972, 20: 872-891.
[51] Rodi W L. A Technique for Improving the Accuracy of Finite Element Solutions for Magnetotelluric Data[J]. Geophysical Journal International. 1976, 44(2): 483-506.
[52] Reddy I K, Rankin D. Magnetotelluric response of a two-dimensional sloping contact by the finite element method [J]. Pure and Applied Geophysics. 1973, 105(1): 847-857.
[53] Philip E Wannamaker J A S L. A stable finite element solution for two-dimensional magnetotelluric modelling[J]. Geophysical Journal International. 1987, 88(1): 277-296.
[54] Coggon J H. Electromagnetic and electrical modeling by the finite element method[J]. Geophysics. 1971, 36(1): 132-155.
[55]谭捍东,余钦范,魏文博.大地电磁法三维交错采样有限差分数值模拟[J].地球物理学报. 2003a, 46(5): 705-711.
[56] Mackie R L, Madden T R, Wannamaker P E. Three-dimensional magnetotelluric modeling using difference equations---Theory and comparisons to integral equation solutions[J]. Geophysics. 1993a, 58(2): 215-226.
[57] Mackie R L, Smith J T, Madden T R. Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example[J]. Radio Sci. 1994, 29(4): 923-935.
[58] Mackie R L, Madden T R, Park S K. A three-dimensional magnetotelluric investigation of the California Basin and Range[J]. J. Geophys. Res. 1996, 101(B7): 221-239.
[59] Zhang J, Mackie R L, Madden T R. 3-D resistivity forward modeling and inversion using conjugate gradients[J]. Geophysics. 1995, 60(5): 1313-1325.
[60] Smith J T. Conservative modeling of 3-D electromagnetic fields, Part I: Properties and error analysis[J]. Geophysics. 1996a, 61(5): 1308-1318.
[61] Smith J T. Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator[J]. Geophysics. 1996b, 61(5): 1319-1324.
[62] Mitsuhata Y, Uchida T. 3D magnetotelluric modeling using the T-Omega finite-element method[J]. Geophysics. 2004, 69(1): 108-119.
[63] Shi X, Utada H, Wang J, et al. Three dimensional magnetotelluric forward modeling using vector finite element method combined with divergence corrections (VFE++)[Z]. IAGA WG 1.2 on Electromagnetic Induction in the Earth,Proceedings of the 17th Workshop, Hyderabad, India,2004.
[64] Smith J T, Booker J R. Rapid Inversion of Two- and Three-Dimensional Magnetotelluric Data[J]. J. Geophys. Res. 1991, 96(B3): 3905–3922.
[65] Degroot-Hedlin C, Constable S. Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data[J]. Geophysics. 1990, 55(12): 1613-1624.
[66] Constable S C, Parker R L, Constable C G. Occam's inversion: A practical algorithm for generating smooth models from electromagnetic sounding data[J]. Geophysics. 1987, 52(3): 289-300.
[67] Siripunvaraporn W, Egbert G. An efficient data-subspace inversion method for 2-D magnetotelluric data[J]. Geophysics. 2000, 65(3): 791-803.
[68] Uchida T. Smooth 2-D inversion for magnetotelluric data based on statistical criterion ABIC[J]. Journal of Geomagnetism and Geoelectricity. 1993a, 45: 841-858.
[69] Uchida T. Inversion of COPROD2 magnetotelluric data by use of ABIC minimization method[J]. Journal of Geomagnetism and Geoelectricity. 1993b, 45: 1063-1071.
[70] Rodi W, Mackie R L. Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J]. Geophysics. 2001, 66(1): 174-187.
[71] Mackie R L, Madden T R. Three-dimensional magnetotelluric inversion using conjugate gradients[J]. Geophysical Journal International. 1993b, 115(1): 215-229.
[72] Newman G A, Alumbaugh D L. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J]. Geophysical Journal International. 2000, 140(2): 410-424.
[73] Zhdanov M S, Fang S, Hursan G. Electromagnetic inversion using quasi-linear approximation[J]. Geophysics. 2000, 65(5): 1501-1513.
[74]胡祖志.大地电磁拟三维反演研究[D].中国地质大学(武汉), 2005.
[75] Avdeeva A, Avdeev D. Three-dimensional magnetotelluric inversion using quasi-Newton minimization[Z]. 4th International Symposium on Three-Dimensional Electromagnetics, Freiberg, Germany, 2007.
[76]谭捍东,余钦范,魏文博.大地电磁法三维快速松弛反演[J].地球物理学报. 2003b, 46(6): 850-855.
[77] Spichak V, Menvielle M, Roussignol M. Three-Dimensional Inversion of MT Fields Using Bayesian Statistics[C]. Three-Dimensional Electromagnetics, Tulsa, USA, 1999.
[78] Spichak V, Popova I. Artificial neural network inversion of magnetotelluric data in terms of three-dimensional earth macroparameters[J]. Geophysical Journal International. 2000, 142(1): 15-26.
[79] Siripunvaraporn W, Egbert G, Lenbury Y, et al. Three-dimensional magnetotelluric inversion: data-space method[J]. Physics of The Earth and Planetary Interiors Electromagnetic Induction in the Earth. 2005a, 150(1-3): 3-14.
[80]刘羽,王家映,孟永良.基于PC机群的大地电磁Occam反演并行计算研究[J].石油物探. 2006, 45(3): 311-315.
[81]谭捍东.电法勘探三维正反演问题研究[D].中国科学院地质与地球物理研究所, 2006.
[82]赵国泽,刘国栋.一种二维大地电磁反演方法[J].地震地质. 1985, 7(1): 56-58.
[83]吴小平,吴广耀,胡建德.二维最平模型的大地电磁快速反演[J].地球科学. 1994, 19(6): 821-830.
[84]张大海,徐世浙.带相位信息的一维大地电磁曲线对比反演法[J].地震地质. 2001a, 23(2):227-231.
[85]张大海,徐世浙.二维M丁快速曲线对比反演方法的可行性研究[J].地震地质. 2001b, 23(2): 232-237.
[86] Meju M A. Joint inversion of TEM and distorted MT soundings: Some effective practical considerations[J]. Geophysics. 1996, 61(1): 56-65.
[87] Moorkamp M, Jones A G, Eaton D W. Joint inversion of teleseismic receiver functions and magnetotelluric data using a genetic algorithm: Are seismic velocities and electrical conductivities compatible?[J]. Geophysical Research Letters. 2007, 34(L16311): doi:10.1029/2007GL030519.
[88] Sasaki Y. Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data[J]. Geophysics. 1989, 54(2): 254-262.
[89]过仲阳,王家林,吴健生.应用遗传算法联合反演地震一大地电磁测深数据[J].石油物探. 1999, 38(1): 102-108.
[90]杨辉,王家林,吴健生,等.大地电磁与地震资料仿真退火约束联合反演[J].地球物理学报. 2002, 45(5): 723-734.
[91]于鹏,张罗磊,王家林,等.重磁电震资料联合反演黔中隆起物性结构[J].同济大学学报(自然科学版). 2008, 36(3): 406-412.
[92]周辉,刘财,渠广学.地震一大地电磁数据非线性联合反演方法研究[J].长春地质学院学报. 1995, 25(2): 205-211.
[93]周辉,杨宝俊,刘财.地震一大地电磁测深资料综合反演[J].地球物理学报. 1994, 37(增刊I): 477-485.
[94]王家映, Oldenburg Douglas, Levy Shlomo.大地电磁的拟地震解释法[J].石油地球物理勘探. 1985, 20(1): 66-79.
[95]于鹏,王家林,吴健生.大地电磁全频段偏移成像结合反演技术[J].同济大学学报(自然科学版). 2005, 33(4): 540-544.
[96]于鹏,王家林,吴健生.有限差分法大地电磁多参数偏移成像[J].地球物理学报. 2001, 44(4): 552-562.
[97]王绪本,毛立峰,高永才.电磁导数场概率成像方法研究[J].成都理工大学学报(自然科学版). 2004, 31(6): 679-684.
[98]魏胜,王家映.二维大地电磁资料的偏移[J].地球物理学报. 1993, 36(2): 256-262.
[99]师学明.二维大地电磁拟波动方程反演法研究[D].中国地质大学, 2000.
[100] Kaufman A A. Reduction of the geological noise in magnetotelluric soundings[J]. Geoexploration. 1988, 25(2): 145-161.
[101] Chave A D, Smith J T. On electric and magnetic galvanic distortion tensor decompositions[J]. J. Geophys. Res. 1994, 99(B3): 4669-4682.
[102] Bahr K. Interpretation of the magnetotelluric impedance tensor: regional induction and local telluric distortion[J]. J. Geophys. 1988, 62: 119-127.
[103] Jiracek G R. Near-surface and topographic distortions in electromagnetic induction[J]. Surveys in Geophysics. 1990, 11(2): 163-203.
[104] Berdichevskiy M N, Dmitriev V I. Distortion of magnetic and electrical fields by near-surface lateral inhomogeneities[J]. Acta Geod., Geophys. Montanist. Acad. Sci. Hung. 1976, 11(3-4): 447-483.
[105] Larsen J C. Low Frequency (0.1-6.0 cpd) Electromagnetic Study of Deep Mantle Electrical Conductivity Beneath the Hawaiian Islands[J]. Geophysical Journal International. 1975, 43(1): 17-46.
[106] Zhang P. Magnetotelluric strike rules[J]. Geophysics. 1987, 52(3): 267-278.
[107] Zhang P, Pedersen L B, Mareschal M, et al. Channelling contribution to tipper vectors: a magnetic equivalent to electrical distortion[J]. Geophysical Journal International. 1993, 113(3): 693-700.
[108] Bahr K. Geological noise in magnetotelluric data: a classification of distortion types[J]. Physics of the Earth and Planetary Interiors. 1991, 66: 24-38.
[109] Pracser E, Szarka L. A correction ro Bahr's "phase deviation" method for tensor decomposition[J]. Earth Planets Space. 1999, 51: 1019-1022.
[110] Groom R W, Bailey R C. Decomposition of Magnetotelluric Impedance Tensors in the Presence of Local Three-Dimensional Galvanic Distortion[J]. J. Geophys. Res. 1989, 94(B2): 1913-1925.
[111] Mcneice G W, Jones A G. Multisite, multifrequency tensor decomposition of magnetotelluric data[J]. Geophysics. 2001, 66(1): 158-173.
[112] Smith J T. Understanding telluric distortion matrices[J]. Geophysical Journal International. 1995, 122(1): 219-226.
[113] Smith J T. Estimating galvanic-distortion magnetic fields in magnetotellurics[J]. Geophysical Journal International. 1997, 130(1): 65-72.
[114] Lilley F E M. Magnetotelluric analysis using Mohr circles[J]. Geophysics. 1993, 58(10): 1498-1506.
[115] Lilley F E M. Magnetotelluric tensor decomposition: Part I, Theory for a basic procedure[J]. Geophysics. 1998a, 63(6): 1885-1897.
[116] Lilley F E M. Magnetotelluric tensor decomposition: Part II, Examples of a basic procedure[J]. Geophysics. 1998b, 63(6): 1898-1907.
[117] Caldwell T G, Bibby H M, Brown C. The magnetotelluric phase tensor[J]. Geophysical Journal International. 2004, 158(2): 457-469.
[118] Bibby H M, Caldwell T G, Brown C. Determinable and non-determinable parameters of galvanic distortion in magnetotellurics[J]. Geophysical Journal International. 2005, 163(3): 915-930.
[119] Utada H, Munekane H. On galvanic distortion of regional three-dimensional magnetotelluric impedances[J]. Geophysical Journal International. 2000, 140(2): 385-398.
[120]赵国泽,汤吉,刘铁胜,等.山西阳高~河北容城剖面大地电磁资料的初步解释-阻抗张量分解技术及其应用[J].地震地质. 1996, 18(1): 66-74.
[121]晋光文,江钊,孙洁.区域三维大地电磁阻抗张量畸变效应的Mohr圆分析[J].地震地质. 1996, 18(4): 382-390.
[122]晋光文,孙洁,江钊.大地电磁阻抗张量不变量及其Mohr圆分析[J].地震地质. 1995, 17(4): 439-445.
[123]孙洁,晋光文,白登海,等. Mohr圆分析方法在川西一藏东MT剖面资料解释中的应用[J].地震地质. 2003, 25(3): 385-393.
[124]王书明.表面局部三维大地电磁曲线畸变校正—MT畸变校正阻抗张量分解方法[J].西北地震学报. 1998, 20(4): 1-11.
[125]杨长福,林长佑,徐世浙.大地电磁GB张量分解法的改进[J].地球物理学报. 2002, 45(增刊): 356-364.
[126]王立凤,晋光文,孙洁,等.一种简单的大地电磁阻抗张量畸变分解方法[J].西北地震学报. 2001, 23(2): 172-180.
[127]高红伟,张胜业.阻抗张量分解进行大地电磁静校正的研究[J].地质科技情报. 1998, 17(1): 91-96.
[128]胡玉平,鲍光淑,贾继标.基于畸变3个因子的MT畸变曲线改正方法与应用[J].地质与勘探. 2005, 41(6): 75-79.
[129] Park S K, Livelybrooks D W. Quantitative interpretation of rotationally invariant parameters in magnetotellurics[J]. Geophysics. 1989, 54(11): 1483-1490.
[130] Szarka L, Menvielle M. Analysis of rotational invariants of the MT impedance tensor[J]. Geophysical Journal International. 1997, 129(1): 133-142.
[131] Szarka L,ádám A, Menvielle M. A field test of imaging properties of rotational invariants of the magnetotelluric impedance tensor[J]. Geophysical Prospecting. 2005, 53(3): 325-334.
[132] Weaver J T, Agarwal A K, Lilley F E M. Characterization of the magnetotelluric tensor in terms of its invariants[J]. Geophysical Journal International. 2000, 141(2): 321-336.
[133] Weaver J T, Agarwal A K, Lilley F E M. The relationship between the magnetotelluric tensor invariants and the phase tensor of Caldwell, Bibby, and Brown[J]. Exploration Geophysics. 2006, 37(3): 261-267.
[134] Yaghjian A D. Electric dyadic Green's functions in the source region[J]. Proceedings of the IEEE. 1980, 68(2): 248- 263.
[135] Habashy T M, Groom R W, Spies B R. Beyond the Born and Rytov Approximations: A Nonlinear Approach to Electromagnetic Scattering[J]. J. Geophys. Res. 1993, 98(B2): 1759-1775.
[136] Marquis G, Jones A G, Hyndman R D. Coincident conductive and reflective middle and lower crust in southern British Columbia[J]. Geophysical Journal International. 1995, 120(1): 111-131.
[137] Park S K, Mackie R L. Resistive (dry?) lower crust in an active orogen, Nanga Parbat, northern Pakistan[J]. Tectonophysics. 2000, 316(3-4): 359-380.
[138] Wannamaker P E. Affordable magnetotellurics: interpretation in natural environments[C]. Three-dimensional Electromagnetics, Tulsa, USA, 1999.
[139] Ledo J, Jones A G. Upper mantle temperature determined from combining mineral composition, electrical conductivity laboratory studies and magnetotelluric field observations: Application to the intermontane belt, Northern Canadian Cordillera[J]. 2005, 236(1-2): 258-268.
[140] Ledo J, Queralt P, Mart A, et al. Two-dimensional interpretation of three-dimensional magnetotelluric data: an example of limitations and resolution[J]. Geophysical Journal International. 2002, 150(1): 127-139.
[141] Fischer G, Masero W. Rotational properties of the magnetotelluric impedance tensor: the example of the Araguainha impact crater, Brazil[J]. Geophysical Journal International. 1994, 119(2): 548-560.
[142] Marti A, Queralt P, Jones A G, et al. Improving Bahr’s invariant parameters using theWAL approach[J]. Geophys. J. Int. 2005, 163: 38-41.
[143] Romo J M, Gomez-Trevino E, Esparza F J. An invariant representation for the magnetic transfer function in magnetotellurics[J]. Geophysics. 1999, 64(5): 1418-1428.
[144] Vozoff K. The magnetotelluric method[C]. Soc. Explor. Geophys, 1991.
[145] MartíA. A Magnetotelluric Investigation of Geoelectrical Dimensionality and Study of the Central Betic Crustal Structure[D]. Universitat de Barcelona, 2006.
[146] Mart A, Queralt P, Roca E. Geoelectric dimensionality in complex geological areas: application to the Spanish Betic Chain[J]. Geophysical Journal International. 2004, 157(3): 961-974.
[147] Kellett R L, Mareschal M, Kurtz R D. A model of lower crustal electrical anisotropy for the Pontiac Subprovince of the Canadian Shield[J]. Geophysical Journal International. 1992, 111(1): 141-150.
[148] Ritter P, Banks R J. Separation of local and regional information in distorted GDS response functions by hypothetical event analysis[J]. Geophysical Journal International. 1998, 135(3): 923-942.
[149] Lezaeta P. Distortion Analysis and 3-D Modeling of Magnetotelluric Data in the Southern Central Andes[D]. Universit¨at Berlin, 2001.
[150] Cerv V, Pek J. Magnetotelluric decomposition via stochastic sampling[Z]. IAGA WG 1.2 on Electromagnetic Induction in the Earth Available at Proceedings of the 17th Workshop, Hyderabad, India, 2004: _volume, S.8-P7.
[151] Becken M, Burkhardt H. An ellipticity criterion in magnetotelluric tensor analysis[J]. Geophysical Journal International. 2004, 159(1): 69-82.
[152] Groom R W, Bailey R C. Analytic investigations of the effects of near-surface three-dimensional galvanic scatterers on MT tensor decompositions[J]. Geophysics. 1991, 56(4): 496-518.
[153] Jones A G, Groom R W. Strike-angle determination from the magnetotelluric impedance tensor in the presence of noise and local distortion: rotate at your peril![J]. Geophysical Journal International. 1993, 113(2): 524-534.
[154] Chave A D, Luther D S, Meinen C S. Correction of Motional Electric Field Measurements for Galvanic Distortion[J]. Journal of Atmospheric and Oceanic Technology. 2004, 21(2): 317-330.
[155] Agarwal A K, Weaver J T. Magnetic distortion of the magnetotelluric tensor: a numerical study[J]. Earth Planets Space. 2000, 52(5): 347-353.
[156]陈小斌.有限元直接迭代算法及其在电磁模拟中的应用研究[D].江汉石油学院, 2000b.
[157] Wannamker P E, Hohmann G W. Electromagnetic induction studies[J]. Reviews of Geophysics. 1991, Supl: 405-415.
[158] Kaufman A A. Tutorial: Distribution of alternating electrical charges in a conducting medium[J]. Geophysical Prosp. 1985, 33: 171-184.
[159] Heise W, Caldwell T G, Bibby H M, et al. Three-dimensional modelling of magnetotelluric data from the Rotokawa geothermal field, Taupo Volcanic Zone, New Zealand[J]. Geophysical Journal International. 2008, 173(2): 740-750.
[160] Thiel S, Heinson G, Gray D R, et al. Ophiolite emplacement in NE Oman: constraints from magnetotelluric sounding[J]. Geophysical Journal International. 2009, 176(3): 753-766.
[161] Stagpoole V M, Bennie S L, Bibby H M, et al. Deep structure of a major subduction back thrust: Magneto-telluric investigations of the Taranaki Fault, New Zealand[J]. Tectonophysics. 2009, 463(1-4): 77-85.
[162] Berdichevsky M N. Marginal Notes On Magnetotellurics[J]. Surveys in Geophysics. 1999, 20(3~4): 341-375.
[163] Groom R W, Bahr K. Corrections for near surface effects: Decomposition of the magnetotelluric impedance tensor and scaling corrections for regional resistivities: A tutorial[J]. Surveys in Geophysics. 1992, 13(4): 341-379.
[164] Jones A G, Ledo J, Ferguson I J. Electromagnetic images of the Trans-Hudson orogen: The North American Central Plains(NACP) anomaly revealed[J]. Canadian Journal of Earth Sciences. 2005, 42: 457-478.
[165] Siripunvaraporn W, Egbert G, Uyeshima M. Interpretation of two-dimensional magnetotelluric profile data with three-dimensional inversion: synthetic examples[J]. Geophysical Journal International. 2005b, 160(3): 804-814.
[166]陈小斌.大地电磁正反演新算法研究及资料处理与解释的可视化集成系统开发[D].中国地震局地质研究所, 2003.
[167]陈小斌,赵国泽,马霄.大地电磁三维模型的一维二维反演近似问题研究[J].工程地球物理学报. 2006, 3(1): 9-15.
[168]陈小斌,赵国泽.大地电磁二维反演的测点“漂移”问题及其解决方法[C].第8届中国国际地球电磁学讨论会,中国荆州2007a.
[169]陈小斌,赵国泽,詹艳.复杂模型下大地电磁有效探测深度的研究初步[C].中国地球物理学会年会,中国成都2007b.
[2] Cagniard L. Basic theory of the magneto-telluric method of geophysicial prospecting[J]. Geophysics. 1953, 18(3): 605-635.
[3]刘国栋,陈乐寿.大地电磁测深研究[M].北京:地震出版社, 1984.
[4]陈乐寿,王光鄂.大地电磁测深法[M].北京:地质出版社, 1990.
[5] Madden T, Nelson P. A defense of Cagniard's magnetotelluric method[R]. Massachusetts Institute Of Technology, 1964.
[6] Berdichevsky M N, Dmitriev V I. Basic principles of interpretation of magnetotelluric sounding curves[C]. Geoelectric and Geothermal Studies, KAPG Geophysical Monograph., Akademiai Kiado, Budapest, 1976.
[7] Morrison H F, Phillips R J, Brien D P. Quantitative interpretation of transient electromagnetic fields over a layered half-space[J]. Geophysical Prospecting. 1969, 17(1): 82-101.
[8] Cantwell T, Madden T R. Preliminary report on crustal magnetotelluric measurements[J]. Journal of Geophysical Research. 1960, 65: 4202-4205.
[9] Bostick F X, Smith H W. Investigation of large-scale inhomogenities in the earth by the magnetotelluric method[J]. Proceeding of the IRE(IEEE). 1962, 50: 2339-2346.
[10] Swift C M. A magnetotelluric investigation of an electrical conductivity anomaly in the southwestern United States[D]. Massachusetts Institute of Technology, 1967.
[11] Sims W E, Bostick F X, Smith H W. The estimation of magnetotelluric impedance tensor elements from measured data[J]. Geophysics. 1971, 36(5): 938-942.
[12] Gamble T D. Remote reference magnetotellurics with SQUIDS[D]. Univ. of California, Berkeley, 1978.
[13] Gamble T D, Goubau W M, Clarke J. Error analysis for remote reference magnetotellurics[J]. Geophysics. 1979a, 44(5): 959-968.
[14] Gamble T D, Goubau W M, Clarke J. Magnetotellurics with a remote magnetic reference[J]. Geophysics. 1979b, 44(1): 53-68.
[15] Sutarno D, Vozoff K. Phase-smoothed robust M-estimation of magnetotelluric impedance functions[J]. Geophysics. 1991, 56(12): 1999-2007.
[16] Egbert G D, Livelybrooks D W. Single station magnetotelluric impedance estimation: Coherence weighting and the regression M-estimate[J]. Geophysics. 1996, 61(4): 964-970.
[17] Larsen J C, Mackie R L, Manzella A, et al. Robust smooth magneto tell uric transfer functions[J]. Geophysical Journal International. 1996, 124(3): 801-819.
[18] Oettinger G, Haak V, Larsen J C. Noise reduction in magnetotelluric time-series with a new signal-noise separation method and its application to a field experiment in the Saxonian Granulite Massif[J]. Geophysical Journal International. 2001, 146: 659-669.
[19] Smirnov M Y. Magnetotelluric data processing with a robust statistical procedure having a high breakdown point[J]. Geophysical Journal International. 2003, 152(1): 1-7.
[20] Chave A D, Thomson D J. Robust, controlled leverage estimation of magnetotelluric responsefunctions[C]. 11th Workshop on EM Induction in the Earth, Wellington, New Zealand, 1992.
[21] Chave A D, Thomson D J. Some comments on magnetotelluric response function estimation[J]. Journal of Geophysical research. 1989, 94(B10): 215-225.
[22] Chave A D, Thomson D J, Ander M E. On the Robust Estimation of Power Spectra, Coherences, and Transfer Functions[J]. J. Geophys. Res. 1987, 92(B1): 633-648.
[23] Egbert G D, Booker J R. Robust estimation of geomagnetic transfer functions[J]. Geophysical Journal International. 1986, 87(1): 173-194.
[24] Chant I J, Hastie L M. Time-frequency analysis of magnetotelluric data[J]. Geophysical Journal International. 1992, 111(2): 399-413.
[25] Raiche A P. An integral equation approach to 3-D modeling[J]. Geoph. J. R. Astr. Soc. 1974, 55: 507-511.
[26] Hohmann G W. Three-dimensional induced polarization and electromagnetic modeling[J]. Geophysics. 1975, 40(2): 309-324.
[27] Ting S C, Hohmann G W. Integral equation modeling of three-dimensional magnetotelluric response[J]. Geophysics. 1981, 46(2): 182-197.
[28] Wannamaker P E. Advances in three-dimensional magnetotelluric modeling using integral equations[J]. Geophysics. 1991, 56(11): 1716-1728.
[29] Wannamaker P E, Hohmann G W, Ward S H. Magnetotelluric responses of three-dimensional bodies in layered earths[J]. Geophysics. 1984, 49(9): 1517-1533.
[30] Zhdanov M S, Fang S. Quasi-linear approximation in 3-D electromagnetic modeling[J]. Geophysics. 1996, 61(3): 646-665.
[31] Xiong Z. Electromagnetic modeling of 3-D structures by the method of system iteration using integral equations[J]. Geophysics. 1992, 57(12): 1556-1561.
[32] Xiong Z, Tripp A C. Electromagnetic scattering of large structures in layered earths using integral equations[J]. Radio Sci. 1995, 30(4): 921-929.
[33] Jones F W, Vozoff K. The calculation of magnetotelluric quantities for three-dimensional conductivity inhomogeneities[J]. Geophysics. 1978, 43(6): 1167-1175.
[34] Lines L R, Jones F W. The perturbation of alternating geomagnetic fields by three-dimensional island structure[J]. Geophys. J. Roy. Astr. Soc. 1973, 32: 133-154.
[35] Pu X, Agarwal A K, Weaver J T. Magnetic field solutions of E-polarization induction problems[J]. J. Geomagn. Geoelectr. 1993, 45: 859-872.
[36]陈小斌,汤吉,赵国泽. MT三维数值模拟中辅助场的求解[C].中国地球物理学会年会,中国昆明2001.
[37]汤吉.长白山天池火山区地壳深部电性结构和大地电磁三维数值模拟研究[D].中国地震局地质研究所, 2001.
[38]陈乐寿.有限元法在大地电磁场正演计算中的应用及改进[J].地球科学-中国地质大学学报. 1981.
[39]陈小斌.有限元直接迭代算法及其在电磁模拟中的应用研究[D].江汉石油学院, 2000a.
[40]李予国,徐世浙,刘斌,等.电导率分块连续变化的二维MT有限元模拟(II)[J].高校地质学报. 1996, 2(4): 448-452.
[41]史明娟,徐世浙,刘斌.大地电磁二次函数插值的有限元法正演模拟[J].地球物理学报. 1997, 40(3): 421-430.
[42]徐世浙,刘斌.电导率分层连续变化的水平层的大地电磁正演[J].地球物理学报. 1995a, 38(2): 262-268.
[43]徐世浙,于涛,李予国,等.电导率分块连续变化的二维MT有限元模拟(I)[J].高校地质学报. 1995b, 1(2): 65-73.
[44]赵生凯,徐世浙.有限单元法计算良导嵌入体的大地电磁测深曲线[J].物探化探计算技术. 1983, 5(1): 47-55.
[45]胡建德,蔡纲.用三角形二次插值法计算二维大地电磁测深曲线[J].石油地球物理勘探. 1984, 19(4): 358-367.
[46]罗焕炎,张淑霞.二维大地电磁测深的边界元与有限元耦合模拟[J].地震地质. 1987, 9(2): 1-9.
[47]王宝钧,高维安.用有限元正演京津塘地区二维构造的大地电磁测深资料[J].地震地质. 1984, 6(1): 69-80.
[48]陈乐寿,孙必俊.有限元法在大地电磁测深中正演计算的改进[J].石油地球物理勘探. 1982, 17(3): 69-72.
[49] Wannamaker P E, Stodt J A, Rijo L. Two-dimensional topographic responses in magnetotellurics modeled using finite elements[J]. Geophysics. 1986, 51(11): 2131-2144.
[50] Silvester P, Hastam C R S. Magnetotelluirc modeling by the finite element method[J]. Geophysical Prospecting. 1972, 20: 872-891.
[51] Rodi W L. A Technique for Improving the Accuracy of Finite Element Solutions for Magnetotelluric Data[J]. Geophysical Journal International. 1976, 44(2): 483-506.
[52] Reddy I K, Rankin D. Magnetotelluric response of a two-dimensional sloping contact by the finite element method [J]. Pure and Applied Geophysics. 1973, 105(1): 847-857.
[53] Philip E Wannamaker J A S L. A stable finite element solution for two-dimensional magnetotelluric modelling[J]. Geophysical Journal International. 1987, 88(1): 277-296.
[54] Coggon J H. Electromagnetic and electrical modeling by the finite element method[J]. Geophysics. 1971, 36(1): 132-155.
[55]谭捍东,余钦范,魏文博.大地电磁法三维交错采样有限差分数值模拟[J].地球物理学报. 2003a, 46(5): 705-711.
[56] Mackie R L, Madden T R, Wannamaker P E. Three-dimensional magnetotelluric modeling using difference equations---Theory and comparisons to integral equation solutions[J]. Geophysics. 1993a, 58(2): 215-226.
[57] Mackie R L, Smith J T, Madden T R. Three-dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example[J]. Radio Sci. 1994, 29(4): 923-935.
[58] Mackie R L, Madden T R, Park S K. A three-dimensional magnetotelluric investigation of the California Basin and Range[J]. J. Geophys. Res. 1996, 101(B7): 221-239.
[59] Zhang J, Mackie R L, Madden T R. 3-D resistivity forward modeling and inversion using conjugate gradients[J]. Geophysics. 1995, 60(5): 1313-1325.
[60] Smith J T. Conservative modeling of 3-D electromagnetic fields, Part I: Properties and error analysis[J]. Geophysics. 1996a, 61(5): 1308-1318.
[61] Smith J T. Conservative modeling of 3-D electromagnetic fields, Part II: Biconjugate gradient solution and an accelerator[J]. Geophysics. 1996b, 61(5): 1319-1324.
[62] Mitsuhata Y, Uchida T. 3D magnetotelluric modeling using the T-Omega finite-element method[J]. Geophysics. 2004, 69(1): 108-119.
[63] Shi X, Utada H, Wang J, et al. Three dimensional magnetotelluric forward modeling using vector finite element method combined with divergence corrections (VFE++)[Z]. IAGA WG 1.2 on Electromagnetic Induction in the Earth,Proceedings of the 17th Workshop, Hyderabad, India,2004.
[64] Smith J T, Booker J R. Rapid Inversion of Two- and Three-Dimensional Magnetotelluric Data[J]. J. Geophys. Res. 1991, 96(B3): 3905–3922.
[65] Degroot-Hedlin C, Constable S. Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data[J]. Geophysics. 1990, 55(12): 1613-1624.
[66] Constable S C, Parker R L, Constable C G. Occam's inversion: A practical algorithm for generating smooth models from electromagnetic sounding data[J]. Geophysics. 1987, 52(3): 289-300.
[67] Siripunvaraporn W, Egbert G. An efficient data-subspace inversion method for 2-D magnetotelluric data[J]. Geophysics. 2000, 65(3): 791-803.
[68] Uchida T. Smooth 2-D inversion for magnetotelluric data based on statistical criterion ABIC[J]. Journal of Geomagnetism and Geoelectricity. 1993a, 45: 841-858.
[69] Uchida T. Inversion of COPROD2 magnetotelluric data by use of ABIC minimization method[J]. Journal of Geomagnetism and Geoelectricity. 1993b, 45: 1063-1071.
[70] Rodi W, Mackie R L. Nonlinear conjugate gradients algorithm for 2-D magnetotelluric inversion[J]. Geophysics. 2001, 66(1): 174-187.
[71] Mackie R L, Madden T R. Three-dimensional magnetotelluric inversion using conjugate gradients[J]. Geophysical Journal International. 1993b, 115(1): 215-229.
[72] Newman G A, Alumbaugh D L. Three-dimensional magnetotelluric inversion using non-linear conjugate gradients[J]. Geophysical Journal International. 2000, 140(2): 410-424.
[73] Zhdanov M S, Fang S, Hursan G. Electromagnetic inversion using quasi-linear approximation[J]. Geophysics. 2000, 65(5): 1501-1513.
[74]胡祖志.大地电磁拟三维反演研究[D].中国地质大学(武汉), 2005.
[75] Avdeeva A, Avdeev D. Three-dimensional magnetotelluric inversion using quasi-Newton minimization[Z]. 4th International Symposium on Three-Dimensional Electromagnetics, Freiberg, Germany, 2007.
[76]谭捍东,余钦范,魏文博.大地电磁法三维快速松弛反演[J].地球物理学报. 2003b, 46(6): 850-855.
[77] Spichak V, Menvielle M, Roussignol M. Three-Dimensional Inversion of MT Fields Using Bayesian Statistics[C]. Three-Dimensional Electromagnetics, Tulsa, USA, 1999.
[78] Spichak V, Popova I. Artificial neural network inversion of magnetotelluric data in terms of three-dimensional earth macroparameters[J]. Geophysical Journal International. 2000, 142(1): 15-26.
[79] Siripunvaraporn W, Egbert G, Lenbury Y, et al. Three-dimensional magnetotelluric inversion: data-space method[J]. Physics of The Earth and Planetary Interiors Electromagnetic Induction in the Earth. 2005a, 150(1-3): 3-14.
[80]刘羽,王家映,孟永良.基于PC机群的大地电磁Occam反演并行计算研究[J].石油物探. 2006, 45(3): 311-315.
[81]谭捍东.电法勘探三维正反演问题研究[D].中国科学院地质与地球物理研究所, 2006.
[82]赵国泽,刘国栋.一种二维大地电磁反演方法[J].地震地质. 1985, 7(1): 56-58.
[83]吴小平,吴广耀,胡建德.二维最平模型的大地电磁快速反演[J].地球科学. 1994, 19(6): 821-830.
[84]张大海,徐世浙.带相位信息的一维大地电磁曲线对比反演法[J].地震地质. 2001a, 23(2):227-231.
[85]张大海,徐世浙.二维M丁快速曲线对比反演方法的可行性研究[J].地震地质. 2001b, 23(2): 232-237.
[86] Meju M A. Joint inversion of TEM and distorted MT soundings: Some effective practical considerations[J]. Geophysics. 1996, 61(1): 56-65.
[87] Moorkamp M, Jones A G, Eaton D W. Joint inversion of teleseismic receiver functions and magnetotelluric data using a genetic algorithm: Are seismic velocities and electrical conductivities compatible?[J]. Geophysical Research Letters. 2007, 34(L16311): doi:10.1029/2007GL030519.
[88] Sasaki Y. Two-dimensional joint inversion of magnetotelluric and dipole-dipole resistivity data[J]. Geophysics. 1989, 54(2): 254-262.
[89]过仲阳,王家林,吴健生.应用遗传算法联合反演地震一大地电磁测深数据[J].石油物探. 1999, 38(1): 102-108.
[90]杨辉,王家林,吴健生,等.大地电磁与地震资料仿真退火约束联合反演[J].地球物理学报. 2002, 45(5): 723-734.
[91]于鹏,张罗磊,王家林,等.重磁电震资料联合反演黔中隆起物性结构[J].同济大学学报(自然科学版). 2008, 36(3): 406-412.
[92]周辉,刘财,渠广学.地震一大地电磁数据非线性联合反演方法研究[J].长春地质学院学报. 1995, 25(2): 205-211.
[93]周辉,杨宝俊,刘财.地震一大地电磁测深资料综合反演[J].地球物理学报. 1994, 37(增刊I): 477-485.
[94]王家映, Oldenburg Douglas, Levy Shlomo.大地电磁的拟地震解释法[J].石油地球物理勘探. 1985, 20(1): 66-79.
[95]于鹏,王家林,吴健生.大地电磁全频段偏移成像结合反演技术[J].同济大学学报(自然科学版). 2005, 33(4): 540-544.
[96]于鹏,王家林,吴健生.有限差分法大地电磁多参数偏移成像[J].地球物理学报. 2001, 44(4): 552-562.
[97]王绪本,毛立峰,高永才.电磁导数场概率成像方法研究[J].成都理工大学学报(自然科学版). 2004, 31(6): 679-684.
[98]魏胜,王家映.二维大地电磁资料的偏移[J].地球物理学报. 1993, 36(2): 256-262.
[99]师学明.二维大地电磁拟波动方程反演法研究[D].中国地质大学, 2000.
[100] Kaufman A A. Reduction of the geological noise in magnetotelluric soundings[J]. Geoexploration. 1988, 25(2): 145-161.
[101] Chave A D, Smith J T. On electric and magnetic galvanic distortion tensor decompositions[J]. J. Geophys. Res. 1994, 99(B3): 4669-4682.
[102] Bahr K. Interpretation of the magnetotelluric impedance tensor: regional induction and local telluric distortion[J]. J. Geophys. 1988, 62: 119-127.
[103] Jiracek G R. Near-surface and topographic distortions in electromagnetic induction[J]. Surveys in Geophysics. 1990, 11(2): 163-203.
[104] Berdichevskiy M N, Dmitriev V I. Distortion of magnetic and electrical fields by near-surface lateral inhomogeneities[J]. Acta Geod., Geophys. Montanist. Acad. Sci. Hung. 1976, 11(3-4): 447-483.
[105] Larsen J C. Low Frequency (0.1-6.0 cpd) Electromagnetic Study of Deep Mantle Electrical Conductivity Beneath the Hawaiian Islands[J]. Geophysical Journal International. 1975, 43(1): 17-46.
[106] Zhang P. Magnetotelluric strike rules[J]. Geophysics. 1987, 52(3): 267-278.
[107] Zhang P, Pedersen L B, Mareschal M, et al. Channelling contribution to tipper vectors: a magnetic equivalent to electrical distortion[J]. Geophysical Journal International. 1993, 113(3): 693-700.
[108] Bahr K. Geological noise in magnetotelluric data: a classification of distortion types[J]. Physics of the Earth and Planetary Interiors. 1991, 66: 24-38.
[109] Pracser E, Szarka L. A correction ro Bahr's "phase deviation" method for tensor decomposition[J]. Earth Planets Space. 1999, 51: 1019-1022.
[110] Groom R W, Bailey R C. Decomposition of Magnetotelluric Impedance Tensors in the Presence of Local Three-Dimensional Galvanic Distortion[J]. J. Geophys. Res. 1989, 94(B2): 1913-1925.
[111] Mcneice G W, Jones A G. Multisite, multifrequency tensor decomposition of magnetotelluric data[J]. Geophysics. 2001, 66(1): 158-173.
[112] Smith J T. Understanding telluric distortion matrices[J]. Geophysical Journal International. 1995, 122(1): 219-226.
[113] Smith J T. Estimating galvanic-distortion magnetic fields in magnetotellurics[J]. Geophysical Journal International. 1997, 130(1): 65-72.
[114] Lilley F E M. Magnetotelluric analysis using Mohr circles[J]. Geophysics. 1993, 58(10): 1498-1506.
[115] Lilley F E M. Magnetotelluric tensor decomposition: Part I, Theory for a basic procedure[J]. Geophysics. 1998a, 63(6): 1885-1897.
[116] Lilley F E M. Magnetotelluric tensor decomposition: Part II, Examples of a basic procedure[J]. Geophysics. 1998b, 63(6): 1898-1907.
[117] Caldwell T G, Bibby H M, Brown C. The magnetotelluric phase tensor[J]. Geophysical Journal International. 2004, 158(2): 457-469.
[118] Bibby H M, Caldwell T G, Brown C. Determinable and non-determinable parameters of galvanic distortion in magnetotellurics[J]. Geophysical Journal International. 2005, 163(3): 915-930.
[119] Utada H, Munekane H. On galvanic distortion of regional three-dimensional magnetotelluric impedances[J]. Geophysical Journal International. 2000, 140(2): 385-398.
[120]赵国泽,汤吉,刘铁胜,等.山西阳高~河北容城剖面大地电磁资料的初步解释-阻抗张量分解技术及其应用[J].地震地质. 1996, 18(1): 66-74.
[121]晋光文,江钊,孙洁.区域三维大地电磁阻抗张量畸变效应的Mohr圆分析[J].地震地质. 1996, 18(4): 382-390.
[122]晋光文,孙洁,江钊.大地电磁阻抗张量不变量及其Mohr圆分析[J].地震地质. 1995, 17(4): 439-445.
[123]孙洁,晋光文,白登海,等. Mohr圆分析方法在川西一藏东MT剖面资料解释中的应用[J].地震地质. 2003, 25(3): 385-393.
[124]王书明.表面局部三维大地电磁曲线畸变校正—MT畸变校正阻抗张量分解方法[J].西北地震学报. 1998, 20(4): 1-11.
[125]杨长福,林长佑,徐世浙.大地电磁GB张量分解法的改进[J].地球物理学报. 2002, 45(增刊): 356-364.
[126]王立凤,晋光文,孙洁,等.一种简单的大地电磁阻抗张量畸变分解方法[J].西北地震学报. 2001, 23(2): 172-180.
[127]高红伟,张胜业.阻抗张量分解进行大地电磁静校正的研究[J].地质科技情报. 1998, 17(1): 91-96.
[128]胡玉平,鲍光淑,贾继标.基于畸变3个因子的MT畸变曲线改正方法与应用[J].地质与勘探. 2005, 41(6): 75-79.
[129] Park S K, Livelybrooks D W. Quantitative interpretation of rotationally invariant parameters in magnetotellurics[J]. Geophysics. 1989, 54(11): 1483-1490.
[130] Szarka L, Menvielle M. Analysis of rotational invariants of the MT impedance tensor[J]. Geophysical Journal International. 1997, 129(1): 133-142.
[131] Szarka L,ádám A, Menvielle M. A field test of imaging properties of rotational invariants of the magnetotelluric impedance tensor[J]. Geophysical Prospecting. 2005, 53(3): 325-334.
[132] Weaver J T, Agarwal A K, Lilley F E M. Characterization of the magnetotelluric tensor in terms of its invariants[J]. Geophysical Journal International. 2000, 141(2): 321-336.
[133] Weaver J T, Agarwal A K, Lilley F E M. The relationship between the magnetotelluric tensor invariants and the phase tensor of Caldwell, Bibby, and Brown[J]. Exploration Geophysics. 2006, 37(3): 261-267.
[134] Yaghjian A D. Electric dyadic Green's functions in the source region[J]. Proceedings of the IEEE. 1980, 68(2): 248- 263.
[135] Habashy T M, Groom R W, Spies B R. Beyond the Born and Rytov Approximations: A Nonlinear Approach to Electromagnetic Scattering[J]. J. Geophys. Res. 1993, 98(B2): 1759-1775.
[136] Marquis G, Jones A G, Hyndman R D. Coincident conductive and reflective middle and lower crust in southern British Columbia[J]. Geophysical Journal International. 1995, 120(1): 111-131.
[137] Park S K, Mackie R L. Resistive (dry?) lower crust in an active orogen, Nanga Parbat, northern Pakistan[J]. Tectonophysics. 2000, 316(3-4): 359-380.
[138] Wannamaker P E. Affordable magnetotellurics: interpretation in natural environments[C]. Three-dimensional Electromagnetics, Tulsa, USA, 1999.
[139] Ledo J, Jones A G. Upper mantle temperature determined from combining mineral composition, electrical conductivity laboratory studies and magnetotelluric field observations: Application to the intermontane belt, Northern Canadian Cordillera[J]. 2005, 236(1-2): 258-268.
[140] Ledo J, Queralt P, Mart A, et al. Two-dimensional interpretation of three-dimensional magnetotelluric data: an example of limitations and resolution[J]. Geophysical Journal International. 2002, 150(1): 127-139.
[141] Fischer G, Masero W. Rotational properties of the magnetotelluric impedance tensor: the example of the Araguainha impact crater, Brazil[J]. Geophysical Journal International. 1994, 119(2): 548-560.
[142] Marti A, Queralt P, Jones A G, et al. Improving Bahr’s invariant parameters using theWAL approach[J]. Geophys. J. Int. 2005, 163: 38-41.
[143] Romo J M, Gomez-Trevino E, Esparza F J. An invariant representation for the magnetic transfer function in magnetotellurics[J]. Geophysics. 1999, 64(5): 1418-1428.
[144] Vozoff K. The magnetotelluric method[C]. Soc. Explor. Geophys, 1991.
[145] MartíA. A Magnetotelluric Investigation of Geoelectrical Dimensionality and Study of the Central Betic Crustal Structure[D]. Universitat de Barcelona, 2006.
[146] Mart A, Queralt P, Roca E. Geoelectric dimensionality in complex geological areas: application to the Spanish Betic Chain[J]. Geophysical Journal International. 2004, 157(3): 961-974.
[147] Kellett R L, Mareschal M, Kurtz R D. A model of lower crustal electrical anisotropy for the Pontiac Subprovince of the Canadian Shield[J]. Geophysical Journal International. 1992, 111(1): 141-150.
[148] Ritter P, Banks R J. Separation of local and regional information in distorted GDS response functions by hypothetical event analysis[J]. Geophysical Journal International. 1998, 135(3): 923-942.
[149] Lezaeta P. Distortion Analysis and 3-D Modeling of Magnetotelluric Data in the Southern Central Andes[D]. Universit¨at Berlin, 2001.
[150] Cerv V, Pek J. Magnetotelluric decomposition via stochastic sampling[Z]. IAGA WG 1.2 on Electromagnetic Induction in the Earth Available at Proceedings of the 17th Workshop, Hyderabad, India, 2004: _volume, S.8-P7.
[151] Becken M, Burkhardt H. An ellipticity criterion in magnetotelluric tensor analysis[J]. Geophysical Journal International. 2004, 159(1): 69-82.
[152] Groom R W, Bailey R C. Analytic investigations of the effects of near-surface three-dimensional galvanic scatterers on MT tensor decompositions[J]. Geophysics. 1991, 56(4): 496-518.
[153] Jones A G, Groom R W. Strike-angle determination from the magnetotelluric impedance tensor in the presence of noise and local distortion: rotate at your peril![J]. Geophysical Journal International. 1993, 113(2): 524-534.
[154] Chave A D, Luther D S, Meinen C S. Correction of Motional Electric Field Measurements for Galvanic Distortion[J]. Journal of Atmospheric and Oceanic Technology. 2004, 21(2): 317-330.
[155] Agarwal A K, Weaver J T. Magnetic distortion of the magnetotelluric tensor: a numerical study[J]. Earth Planets Space. 2000, 52(5): 347-353.
[156]陈小斌.有限元直接迭代算法及其在电磁模拟中的应用研究[D].江汉石油学院, 2000b.
[157] Wannamker P E, Hohmann G W. Electromagnetic induction studies[J]. Reviews of Geophysics. 1991, Supl: 405-415.
[158] Kaufman A A. Tutorial: Distribution of alternating electrical charges in a conducting medium[J]. Geophysical Prosp. 1985, 33: 171-184.
[159] Heise W, Caldwell T G, Bibby H M, et al. Three-dimensional modelling of magnetotelluric data from the Rotokawa geothermal field, Taupo Volcanic Zone, New Zealand[J]. Geophysical Journal International. 2008, 173(2): 740-750.
[160] Thiel S, Heinson G, Gray D R, et al. Ophiolite emplacement in NE Oman: constraints from magnetotelluric sounding[J]. Geophysical Journal International. 2009, 176(3): 753-766.
[161] Stagpoole V M, Bennie S L, Bibby H M, et al. Deep structure of a major subduction back thrust: Magneto-telluric investigations of the Taranaki Fault, New Zealand[J]. Tectonophysics. 2009, 463(1-4): 77-85.
[162] Berdichevsky M N. Marginal Notes On Magnetotellurics[J]. Surveys in Geophysics. 1999, 20(3~4): 341-375.
[163] Groom R W, Bahr K. Corrections for near surface effects: Decomposition of the magnetotelluric impedance tensor and scaling corrections for regional resistivities: A tutorial[J]. Surveys in Geophysics. 1992, 13(4): 341-379.
[164] Jones A G, Ledo J, Ferguson I J. Electromagnetic images of the Trans-Hudson orogen: The North American Central Plains(NACP) anomaly revealed[J]. Canadian Journal of Earth Sciences. 2005, 42: 457-478.
[165] Siripunvaraporn W, Egbert G, Uyeshima M. Interpretation of two-dimensional magnetotelluric profile data with three-dimensional inversion: synthetic examples[J]. Geophysical Journal International. 2005b, 160(3): 804-814.
[166]陈小斌.大地电磁正反演新算法研究及资料处理与解释的可视化集成系统开发[D].中国地震局地质研究所, 2003.
[167]陈小斌,赵国泽,马霄.大地电磁三维模型的一维二维反演近似问题研究[J].工程地球物理学报. 2006, 3(1): 9-15.
[168]陈小斌,赵国泽.大地电磁二维反演的测点“漂移”问题及其解决方法[C].第8届中国国际地球电磁学讨论会,中国荆州2007a.
[169]陈小斌,赵国泽,詹艳.复杂模型下大地电磁有效探测深度的研究初步[C].中国地球物理学会年会,中国成都2007b.
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |