各向异性介质GPR非结构化网格有限元正演

为了更准确地认识电磁波在各向异性介质中的传播特征,从麦克斯韦方程组出发,推导电磁场在单斜各向异性介质中的波动方程。采用Delaunay非结构化网格有限元法进行空间域离散,采用中心差分公式进行时间域离散,导出单斜各向异性介质中GPR时空域有限元方程。在此基础上,建立各向异性介质GPR非结构化网格有限元正演算法。利用以上算法编制的程序对均匀各向异性介质模型进行正演计算,并与解析解对比,验证该算法的正确性和有效性;对均匀各向异性介质模型和圆形异常体模型进行计算,并与背景介质为各向同性介质的计算结果进行对比。研究结果表明:与FDTD计算结果对比,非结构化网格有限元法的拟合误差更小,精度更高;电磁波传播受介质各向异性的影响,波前面呈椭圆状向外扩散传播,其能量衰减较慢;与各向同性介质中圆形异常体反射波相比,各向异性介质的圆形异常体反射波的曲率、能量和双程走时,受不同方向上电磁波速度影响会发生显著变化。