Abstract: This paper describes how to apply the boundary element method to solve 2-D magneto-telluric sounding problem and then couple the BEM with finite element method to deal with inhomogeneous geological structure and its relation to regional stress regime. The BEM is able to calculate precisely the derivative of field variables. Such a new approach can more flexibly than others set up the upper boundary condition of the air region. Two practical examples are given. One is the study of topographical effect on the MT fields using boundary element method. There is an inverse relationship between the change of apparent resistivity and the topographical variation. The other is a coupled model of a basin in Nei Monngol Autonomous Region. A lower resistivity zone must exist in the crust and should be thicker below the central portion of the basin. This seems to agree with the regional stress regime changed from extension in Mesozoic to compression in Cenozoic.

Key words: Magnetotelluric sounding, Boundary element, Finite element, Two dimensional coupling model

摘要： 本文介绍了在大地电磁测深问题中边界单元和有限单元耦合模拟的具体过程。通过与解析解、标准大地电磁测深的模型曲线对比,证明了这个耦合方法能够有效地模拟大地电磁场的变化规律。它能够直接计算场的导数,提高辅助场值的计算精度,对保证实测视电阻率曲线解释结果的可靠性是具有重要意义的。此外,它既能灵活地处理地下非均质的电性结构,又能研究半无限空气域内的电场变化,这为解决E极化数值模型的边界选择提供了有效途径。

关键词: 大地电磁测深, 边界元, 有限元, 二维耦合模拟