摘要：近年来,探地雷达（GPR）凭借其快速、高效、无破损等特点,已经广泛应用于浅地层目标探测中.数值模拟是研究探地雷达电磁波在地下结构中传播规律的有效手段.辛算法是一种保持Hamilton系统总能量不变的时域数值计算方法.本文提出了基于一阶显式辛分块龙格库塔方法的探地雷达数值模拟方法.通过对比本文算法与时域有限差分方法计算结果可知,在同等计算精度下,本文算法可以节省25%的计算时间.并基于本文算法对两个复杂GPR模型进行正演模拟,得到模拟GPR探测wiggle图,这有助于更好的理解和分析实测雷达数据.

关键词：探地雷达 一阶辛方法 透射边界 正演模拟 

单位：大连理工大学建设工程学部 大连116024 大连理工大学海岸和近海工程国家重点实验室 大连116024 郑州大学水利与环境学院 郑州450001