童孝忠，何 婷，赵理芳，李爱勇

(1.中南大学 地球科学与信息物理学院，湖南 长沙 410083；2.中南大学 有色资源与地质灾害探查湖南省重点实验室，湖南 长沙 410083；3.江苏省有色金属华东地质勘查局 八一四队，江苏 镇江 212005)

1 引 言

大地电磁测深(Magnetotelluric，简称MT)是以天然电磁场为场源来研究地球内部电性结构的一种重要的地球物理手段，在深部地球动力学研究、矿产与油气资源勘探、地热资源勘探等领域得到广泛应用[1]。但大地电磁勘探过程中，需要将电磁信号转换为勘探目标的电性变化，进而达到查明断裂、构造层以及岩体的空间分布特征。这一过程离不开大地电磁勘探数据的处理与定性定量解释，其中包括数据预处理、正演模拟与反演成像等。大地电磁正反演解释需要以勘探对象的实际电性结构变化为基础，可根据实际地质资料针对性地开展一维、二维和三维正反演解释。

大地电磁测深的正演问题求解与其他地球物理勘探方法的正演问题类似，不外乎有物理模拟和数值模拟两种方法，但物理模拟对制作者要求较高，并且很难实现复杂模型的模拟，而伴随着计算机的快速发展，大数据量的存储和大型线性方程组的求解已成为可能，所以数值模拟对于复杂模型更容易实现且模拟精度较高。大地电磁测深正演模拟的数值方法主要有3种：有限差分法[2-5]、有限单元法[6-8]和积分方程法[9，10]。前两种方法常常被应用于二维地电模型的数值模拟，而第3种方法主要被应用于三维地电模型的数值计算。……