薛国强,常江浩,雷康信,陈 康

(1.中国科学院地质与地球物理研究所 中国科学院矿产资源研究重点实验室，北京 100029；2.中国科学院大学 地球与行星科学学院，北京 100049；3.中国科学院地球科学研究院，北京 100029；4.青海省第三地质勘查院，青海 西宁 810000；5.河北地质大学 河北省战略性关键矿产资源重点实验室，河北 石家庄 050031；6.中国水利水电科学研究院 流域水循环模拟与调控国家重点实验室，北京 100038)

0 引 言

瞬变电磁法(Transient Electromagnetic Method，TEM)具有纯二次场观测、施工效率高等特点，在金属和非金属矿、工程地质和煤田水文地质等领域得到了广泛应用。该方法发射信号频带宽、频谱信息丰富，一次激发可覆盖探测所需的频段[1]。由于探测环境越来越复杂，发展瞬变电磁场精细模拟方法就变得更加重要[2]。

自20世纪70年代开始，国内外学者就开始研究瞬变电磁法响应的一维正演模拟。一维正演模拟采用的思路一般是先在频率域采用分离变量法解出层状介质条件下的频域电磁响应，然后将频域电磁响应转变到时间域。朴化荣通过快速Hankel变换求解贝塞尔函数得到一维模型在频率域的解，并采用G-S逆拉氏变换法将电磁场从拉普拉斯域转换到时间域，实现了瞬变电磁响应一维正演[3]。一维层状模型只适用于简单的地层结构；对于局部异常和复杂目标体的分析，则需要利用高维数值模拟的手段，根据不同的地质条件建立模型，以便获得更加符合实际情况的电磁响应[4-7]。

Kuo等采用二维有限元法对实际矿体模型的瞬变电磁响应进行了研究[8]；……