齐彦福， 智庆全,3， 李貅*， 景旭, 戚志鹏，孙乃泉， 周建美， 刘文韬

1 长安大学地质工程与测绘学院， 西安 7100542 长安大学地球物理场多参数综合模拟实验室， 西安 7100543 中国地质科学院地球物理地球化学勘查研究所， 河北 廊坊 065000

0 引言

随着经济的快速发展，我国对矿产资源的需求不断增加，然而经过数十年的不间断开采，地形相对平坦、地质条件简单的浅部矿藏已经开采殆尽，找矿工作正在向地形起伏严重的山区进军.地面瞬变电磁法作为一种低成本、高效率的地球物理勘探方法，已经在矿产资源勘查领域获得广泛应用(薛国强等，2007；底青云等，2019).然而，在过去的数十年中受到计算机条件和数据解释技术的制约，地面瞬变电磁数据解释工作主要以视电阻率成像和一维反演方法为主(薛国强等，2008；石显新等,2009；武军杰等，2013)，此类方法在地形平坦、地下介质成层性明显的地区可以获得很好的应用效果，然而在地形起伏严重、三维地电结构发育的山区，利用一维处理解释技术很难获取准确的地下三维电性分布信息.开发适用于起伏地形的瞬变电磁三维反演算法，已成为瞬变电磁数据精细化解释的关键.

作为反演的基础，前人在瞬变电磁三维正演算法方面开展了大量工作，开发了一系列三维正演方法(汤井田等，2007)，例如频时转换法(殷长春等，2013)、时间域有限差分法(Wang and Hohmann，1993；Commer and Newman，2004)、时间域有限体积法(Haber et al.，2007；Oldenburg et al.，2013)和时间域有限元法(Yin et al.，2016；Cai et al.，2017).频时转换法首先计算三维频率域电磁响应，然后采用频……