探地雷达信号的高阶时间域有限差分模拟精选.pdfVIP

第 3 1卷 　第 4 期 物探化探计算技术 　　 2009年 7月 　 　 文章编号 : 100 1—1749 (2009) 04 —0308—06 　 　 探地雷达信号的高阶时间域有限差分模拟 吴丰收 , 曾昭发 , 黄 　玲 , 张玲玲 (吉林大学 　地球探测科学与技术学院 ,吉林 长春 　130026) ( ( 摘 　要 : 探地雷达信号时域有限差分法模拟 , 多采用二阶精度的中心差分法近似 FD TD 2 , 2) ) ,虽然计算简单 ,但数值色散误差较大 ,影响了模拟精度 。在解决较复杂介质分布的探地雷 达信号时 ,不能很好地反映信号的精细变化 。而高阶时域有限差分模拟能减少数值色散带来的 ( ) 误差 ,提高了模拟的精度 。采用各向异性完全匹配层 U PML 作为吸收边界条件 ,实现了高阶计 算 ,有效地吸收边界电磁波反射 ,而且提高了计算效率 。通过模拟的结果分析可知 ,高阶时域有 限差分法能很好地提高模拟精度 。 关键词 : 探地雷达 ;时域有限差分 ;各向异性完全匹配层 ;正演模拟 中图分类号 : P 63 1. 4　　文献标识码 : A ( ( ) ) 具有二阶精度的中心差分法近似 FD TD 2 , 2 。 0　前言 虽然形式简单 ,但数值色散误差较大 [ 9～11 ] 。数值 色散是影响 FD TD 法模拟电磁现象精度的重要因 探地 雷 达 ( Ground Penetrating R adar, 简 称 素 ,直接关系到计算所需要时间及对具体对象模拟 ) 的可行性 。目前解决数值色散影响 ,是通过减小离 GPR 是根据高频脉冲电磁波在地下介质中传播的 特性 ,来探测地下结构和特性的方法 。因其具有快 散过程所取时间和空间步长 ,但这意味着要增加计 捷 、分辨率高 、非破坏性等特点 , 已广泛应用于工程 算网格的空间数 目,相应地增加了对计算机存储空 勘察 [ 1 ] ,环境监测 ,水文地质调查 [ 2、3 ] 及考古 [ 4 ] 等 间和计算时间的要求 [ 10、11 ] ,使该方法在应用方面 领域 ,成为浅层勘探的重要工具 [ 5 ] 。探地雷达的 受到 了一定 限制 。而高阶 FD TD 法 ( FD TD ( 2 , 正演模拟是探地雷达资料解释的重要内容 ,通过对 2M ) )具有良好的数值色散特性 ,提高了电磁场的 电磁波模拟 ,可以加深对探地雷达数据和电磁波传 时域模拟精度 [ 12 ] 。 播规律的认识 ,提高解释精度 [ 6 ] 。由于计算技术 作者以二维为例 ,在时间上采用二阶差分近 (