基于单轴各向异性完全匹配层边界条件的adifdtd-中南大学学报.pdf

第 42卷第 8期 中南大学学报(自然科学版)  Vol.42  No.8  2011 年 8 月  Journal of Central South University (Science and Technology)  Aug. 2011  基于单轴各向异性完全匹配层边界条件的  ADI­FDTD 三维 GPR 全波场正演 冯德山，谢源 (中南大学 地球科学与信息物理学院，湖南 长沙，410083) 摘要：基于交替方向隐式差分(ADI­FDTD)法突破了 CFL 稳定性条件的约束，具有无条件稳定性的特点，而单轴 各向异性完全匹配层(UPML)边界条件具有宽频带吸收特性、不需要对电场和磁场进行分裂、迭代公式简单、便 于编程的特点，提出基于 UPML 边界条件的 ADI­FDTD 三维探地雷达数值模拟算法；通过对三维 Maxwell 方程 进行离散化，推导 GPR波的 ADI­FDTD及其 UPML边界条件 2个子时间步的迭代差分公式，以及计算步骤。开 发 ADI­FDTD 三维模拟程序，并对三维复杂 GPR 模型进行正演模拟，分析该模型的剖面图、切片图及不同时刻  6 个分量的全波场快照。研究结果表明：由全波场快照可以了解雷达波形在空间中的传播过程及变化规律，有助 于对雷达资料进行更可靠、更准确的解释；基于 UPML边界条件的 ADI­FDTD算法消除了截断边界处的强反射， 能对三维复杂的 GPR模型进行快速、高效模拟。 关键词：探地雷达；交替方向隐式有限差分法；单轴各向异性完全匹配层；正演 中图分类号：P631  文献标志码：A  文章编号：1672−7207(2011)08−2363−09  Three dimensional GPR numerical simulation of full wave field based on  UPML boundary condition of ADI­FDTD  FENG De­shan, XIEYuan  (School of Geosciences and Info­Physics, Central South University, Changsha410083, China)  Abstract: Based on the fact that ADI­FDTD can break through the restriction of the CFL stability condition and its  unlimited stable, however, UPML boundary condition has the characteristic of absorbing wide frequency and without the  necessity for separating electric field from magnetic field, and its simple iterative formula is convenient for compiling  program.  The  GPR numerical  simulation  arithmetic  of  ADI­FDTD  was put  forward based  upon  UPML  boundary  condition.  Through  the  discretization  of  the  three­dimension  Maxwell  equations,  two  sub­time  step  iterative  finite