粟　琪　戴世坤　赵东东

(①中南大学地球科学与信息物理学院，湖南长沙 410083；②中南大学有色金属成矿预测与地质环境监测教育部重点实验室，湖南长沙 410083)

1　引言

随着频率域可控源电磁法和计算机技术的不断发展，二维和三维电磁场的数值模拟及反演得到了很大的重视，而2.5维表现为地球物理模型的二维性和场源的三维性，相比于二维场源模拟更加接近于实际场源勘探的结果。早期关于2.5维可控源电磁法的研究大多基于各向同性介质，Coggon[1]首先把有限单元法应用于2.5维电磁法的模拟中； Stoyer等[2]利用有限差分法进行频率域2.5维磁偶极子源电磁响应的正演模拟； Doherty[3]运用积分方程法求解2.5维磁偶极子源的电磁响应； Unsworth等[4]基于有限元详细研究了2.5维电性源激发产生的电磁场响应特征； Torres等[5]提出了一种2.5维数值模拟的快速算法； Mitsuhata[6]实现了基于总场的可控源电磁法2.5维正演模拟，利用网格单元场函数插值公式与等参单元有效模拟了地形；Li等[7-9]等基于非结构化三角网格、采用预制条件的共轭梯度法和准最小残差法求解方程，实现了2.5维海洋电磁自适应有限元正演模拟； 柳建新等[10]重点研究了基于双二次插值的频率域可控源电磁法有限元正演模拟，并对波数选取进行研究，给出了一组兼顾计算速度和模拟精度的波数； Unsworth等[11]在有限元正演的基础上提出了一种子空间海洋可控源电磁的二维反演算法，引入二维OCCAM算法并得以成功应用；Abubakar等[12]研究了2.5维低频电磁法正反演，实现了一种快速可靠的反演算法；……