基于COMSOL的起伏地形三维大地电磁数值模拟北大核心CSCD

COMSOL是一个多物理场仿真软件,笔者尝试将其应用于大地电磁法的数值模拟,研究起伏地形对大地电磁响应的影响及解决方案。本文由麦克斯韦方程组出发,施加第一类边界条件,将地形起伏数字化为参数化曲面,导入COMSOL进行计算。在对COMSOL的精度进行了验证的基础上,分析不同频点所需内存以及自由度,然后选择合适的频点进行研究。将三维起伏地形的模拟结果与无地形起伏时的结果进行对比,总结出山谷、山峰地形对三维大地电磁响应的影响,发现在XY模式下,地形对视电阻率的影响大于相位,在视电阻率等值线图中能够清晰地看出起伏地形的轮廓,而在YX模式下,地形对相位的影响大于对视电阻率的影响。最后,利用“比值法”进行地形校正,校正前由于起伏地形影响,在异常体附近会出现部分视电阻率数值跳跃现象,通过校正有效抑制了起伏地形对异常体附近视电阻率的影响,校正后的视电阻率曲线整体上移,并且异常体周围的视电阻率值也趋于背景场值,获得了较好的效果。