三角网格剖分时域有限元法的探地雷达正演模拟

摘要

在探地雷达正演方法中,射线追踪法只能获取了波的运动学特征,时域有限差分法虽然能得到波的动力学特征,但是由于时空离散采用规则的结构化Yee单元,不能模拟复杂地质构造下波的传播规律.为了构建更复杂的地电模型和获取雷达波在不同地下构造中传播的波场特征,本文采用三角网格剖分的时域有限元法对探地雷达做正演模拟.时域有限元法可采用非结构化的拓扑网格对求解区域进行网格划分,选择三角形网格剖分物性区域可进一步解决矩形网格剖分对物性参数分布复杂以及几何特征不规则模型适应性差的问题.时域有限元法的求解首先以电磁波的麦克斯韦方程为基础,根据电磁场的本构关系,导出雷达波的波动方程,将波动方程作为微分控制方程并定义时域电磁场中的边值问题,对求解区域进行三角网格剖分后,采用线性插值基函数,应用Galerkin加权余量法对微分方程求解形成线性矩阵方程组,并在解线性方程组时采用LU分解的双共轭梯度法,加快了收敛速度,减小了迭代次数.在时间差分近似时采用Newmark-β法消除由时间步长过大带来的数值色散,实现无条件稳定.在求解区域边界外加完全匹配层吸收外向传播的电磁波,压制来自截断边界处的反射波.通过对构建的地电模型做正演模拟,其研究结果表明:使用三角网格剖分的时域有限元法进行复杂地质构造下的GPR正演数值模拟,构建的正演模型更加接近实际地质构造,具有更高的模拟精度,能进一步解译雷达数据剖面.