基于Debye模型的时间域航空电磁激电效应正演模拟

时间域航空电磁法是一种高效的地球物理勘查方法。由于激电效应的存在，晚期道电磁数据经常会出现符号反转现象。此时，传统的不考虑激电效应的成像和反演方法无法恢复地下真实电性结构，因此需要开发含激电效应的电磁正反演模拟技术。本文提出了一种基于非结构有限元的时间域航空电磁激电效应三维正演算法。为描述介质的激电效应，引入Debye模型，利用频时转换将其变成时间域导数形式，进而采用无条件稳定的二阶后推欧拉方法进行时间离散。之后，采用基于非结构四面体网格的矢量有限元法，对麦克斯韦方程进行空间离散，实现含激电效应的时间域航空电磁正演问题求解。通过与均匀半空间模型的一维半解析解进行对比，验证了该算法的准确性。最后，计算了半空间和三维地质体具有不同极化条件时的电磁响应，发现可以根据航空电磁响应曲线的衰减特征，有效识别围岩及异常体的极化特征。此外，激电参数对时间域航空电磁系统的影响特征研究表明，地质体电阻率越高或者充电率越强，极化效应影响越严重。