基于upml 吸收边界条件的gpr 有限元数值模拟 - 中国有色金属学报.pdf

第 23 卷第 7 期 中国有色金属学报  2013 年 7 月  Vol.23 No.7  The Chinese Journal of Nonferrous Metals  July2013  文章编号：1004­0609(2013)07­2003­09  基于 UPML 吸收边界条件的 GPR 有限元数值模拟 1  1, 2  王洪华  ，戴前伟  (1.  中南大学 地球科学与信息物理学院，长沙 410083；  2.  中南大学 地球科学与信息物理学院 有色金属成矿预测教育部重点实验室，长沙  410083)  摘 要：为了进一步减小截断边界对探地雷达(GPR)数值模拟精度的影响、提高复杂 GPR地电模型数值模拟的精 度， 采用具有良好宽频带吸收特性的单轴各向异性完全匹配层(UPML)吸收边界条件与有限单元法(FEM)相结合的 算法来模拟 GPR波在复杂地电模型中的传播。从 GPR满足的波动方程出发， 利用伽辽金法导出二维 GPR有限元 方程。然后，对 UPML 区域满足的两个频率域旋度方程采用傅里叶变换，推导 GPR 时间域有限元波动方程及其 求解方法，并建立相应的 GPR数值模拟算法。应用以上算法编制的程序对两个复杂 GPR地电模型进行了正演模 拟，得到了相应的雷达剖面图。模拟结果表明：UPML 吸收边界条件能够对截断边界处的强烈反射波进行充分的 吸收，大大减弱了截断边界处的强反射；FEM 能够对复杂的 GPR地电模型进行高精度的模拟。基于 UPML 吸收 边界条件的 FEM 算法能够对复杂的 GPR地电模型进行快速、高效的模拟。 关键词：UPML；吸收边界；探地雷达；数值模拟；有限单元法 中图分类号：P631  文献标志码：A  Finite elementnumerical simulation for GPR based on  UPML boundary condition  1  1,2  WANG Hong­hua , DAI Qian­wei  (1. School of Geosciences and Info­Physics, Central South University, Changsha410083,China;  2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education,  School of Geosciences and Info­Physics, Central South University, Changsha 410083, China)  Abstract: To decrease the influence of truncation boundaries on the GPR numerical simulation precision availably and  improve the precision of complex GPR geo­electric model numerical simulation, the FEM, co