基于UPML边界条件的LOD-FDTD方法GPR波的数值模拟

西安科技大学 通信与信息工程学院,陕西 西安 710054

GPR; UPML; LOD-FDTD方法; 波场快照; 雷达剖面图

GPR wave numerical simulation based on the UPML boundary conditions of LOD-FDTD method
HAN Xiao-bing,WEI Hai-liang,LI Qi

(College of Communication and Information Engineering,Xi'an University of Science and Technology,Xi'an 710054,China)

GPR; UPML; LOD-FDTD algorithm; wave field snap; radar profile

DOI: 10.13800/j.cnki.xakjdxxb.2016.0516

备注

提出局部一维时域有限差分方法(LOD-FDTD)和单轴各向异性理想匹配层(UPML)边界条件相结合的LOD-FDTD探地雷达(GPR)电磁波传播数值模拟算法。通过将Maxwell方程离散化,导出三维GPR波LOD-FDTD算法迭代差分形式,引入UPML边界条件、给出公式推导。在此基础上,建立UPML边界条件应用模型、LOD-FDTD算法性能分析模型,验证UPML边界条件在数值模拟中对GPR波具有良好的吸收特性和算法在时间步不受CFL条件的限制下计算的高效性。建立LOD-FDTD方法的GPR波的正演模型,得到波场快照及其雷达剖面图,为了解电磁波在复杂介质中的传播特性和雷达数据的解释提供了依据。

For simulating the transmission of Ground Penetrating Radar's electromagnetic wave in complex medium,a LOD-FDTD algorithm of the combination of LOD-FDTD and UPML boundary conditions is put forward.Through the discretization of the three-dimentional Maxwell's equation,sub-time step iterative finite differential form of LOD-FDTD algorithm is derived.Introducing UPML boundary conditions and its formula derivation,then the UPML boundary condition application model and LOD-FDTD algorithm precision analysis model are established to prove the absorbing validation of UPML boundary condition in numerical simulation of GPR wave and the efficiency of the algorithm under the time step being not restricted by CFL condition.Based on establishing the forward model of GPR wave by LOD-FDTD method,wave field snap and shot radar profile are captured,which can be used to better understand the transmission characteristics of electromagnetic wave in complex medium and the radar data interpretation.