摘要
近年来,探地雷达(GPR)凭借其快速、高效、无破损等特点,已经广泛应用于浅地层目标探测中.数值模拟是研究探地雷达电磁波在地下结构中传播规律的有效手段.辛算法是一种保持Hamilton系统总能量不变的时域数值计算方法.本文提出了基于一阶显式辛分块龙格库塔方法的探地雷达数值模拟方法.通过对比本文算法与时域有限差分方法计算结果可知,在同等计算精度下,本文算法可以节省25%的计算时间.并基于本文算法对两个复杂GPR模型进行正演模拟,得到模拟GPR探测wiggle图,这有助于更好的理解和分析实测雷达数据.
In recent years, the ground penetrating radar （GPR）, which is characterized by its speediness, high efficiency, non-destruction, etc. , has been widely used for detection in shallow subsurface. Numerical simulation is an effective measure for the research on GPR wave propagation in underground structure. The symplectic method is a kind of time-domain numerical scheme designed to conserve the total energy of the Hamiltonian system. This paper presented a numerical method for GPR simulation based on the 1st-order symplectic partitioned Runge-Kutta （SPRK） method. Comparison of the numerical results is made between SPRK method and FDTD method. It can be found that the proposed method can reach the same level of accuracy as FDTD scheme. But the CPU time is reduced by 25 % of that by FDTD scheme. The wiggle map of two complicated GPR model is simulated by the proposed algorithm, which has a better understanding of the real GPR profiles.
出处
《地球物理学报》
SCIE
EI
CAS
CSCD
北大核心
2013年第2期653-659,共7页
Chinese Journal of Geophysics
基金
中德科学基金项目(GZ566)
国家自然科学基金重点项目(51138001)资助
关键词
探地雷达
一阶辛方法
透射边界
正演模拟
Ground penetrating radar, lst-order symplectic method, Transmitting boundary,Forward simulation
