Borehole-GPR numerical simulation of full wave field based on convolutional perfect matched layer boundary

Abstract

The absorbing boundary is the key in numerical simulation of borehole radar. Perfect match layer (PML) was chosen as the absorbing boundary in numerical simulation of GPR. But CPML (convolutional perfect match layer) approach that we have chosen has the advantage of being media independent. Beginning with the Maxwell equations in a two-dimensional structure, numerical formulas of finite-difference time-domain (FDTD) method with CPML boundary condition for transverse electric (TE) or transverse magnetic (TM) wave are presented in details. Also, there are three models for borehole-GPR simulation. By analyzing the simulation results, the features of targets in GPR are obtained, which can provide a better interpretation of real radar data. The results show that CPML is well suited for the simulation of borehole-GPR.

Keywords

borehole-GPR / numerical simulation / convolutional perfect match layer / finite-difference time-domain method

References

[1] HarryM JGround Penetrating Radar: Theory and Applications [M], 2009AmsterdamElsevier Science1-25

[2] MichaelA MMonitoring dielectrical property changes and groundwater velocity changes in a contaminated aquifer using borehole ground penetrating radar and multi-level point velocity probes [D], 2004ScotlandUniversity of St. Andrews19-321

[3] CaiJ, McmechanG A. Ray-based synthesis of bistatic ground-penetrating radar profiles [J]. Geophysics, 1995, 60(1): 87-96

[4] ThomasS, HansruediM, HeikeW. Characterization of an unstable rock mass based on borehole logs and diverse borehole radar data [J]. Journal of Applied Geophysics, 2007, 61(1): 16-38

[5] EllefsenK J. Effects of layered sediments on the guided wave in cross-well radar data [J]. Geophysics, 1999, 64(6): 1698-1707

[6] CarcioneJ M. Ground-penetrating radar: Wave theory and numerical simulation in lossy anisotropic media [J]. Geophysics, 1996, 61(6): 1664-1677

[7] BorchertO, BehaimanotK, GlasmachersA. Directional borehole radar calibration [J]. Journal of Applied Geophysics, 2009, 67(4): 352-360

[8] BergmannT, RobertssonJ O A, HolligerK. Numerical properties of staggered finite-difference solutions of Maxwell’s equations for ground penetrating radar modeling [J]. Geophysical Research Letters, 1996, 23(1): 45-48

[9] BergmannT, RobertssonJ O A, HolligerK. Finite difference modeling of electromagnetic wave propagation in dispersive and attenuating media [J]. Geophysics, 1998, 63(3): 856-867

[10] JungH K, SamG P, MyeongJ Y. Borehole radar investigations for locating ice ring formed by cryogenic condition in an underground cavern [J]. Journal of Applied Geophysics, 2007, 62(3): 204-214

[11] ChenY H, ChewW C. Application of perfectly matched layers to the transient modeling of subsurface EM problems [J]. Geophysics, 1997, 62(6): 1730-1736

[12] FangJ, WuZ. Generalized perfectly matched layer for the absorption of propagating and evanescent waves in lossless and lossy media [J]. IEEE Transactions on Microwave Theory and Techniques, 1996, 44(12): 2216-2222

[13] GiannopoulosA. Modelling ground penetrating radar by GprMax [J]. Construction and Building Materials, 2005, 19(10): 755-762

[14] CassidyN J, MillintonT M. The application of finite-difference time-domain modeling for the assessment of GPR in magnetically lossy materials [J]. Journal of Applied Geophysics, 2009, 67(4): 296-308

[15] TafloveA. Review of the formulation and applications of the finite-difference time-domain method for numerical modeling of electromagnetic wave interactions with arbitrary structures [J]. Wave Motion, 1988, 10(6): 547-582

[16] GiannopoulosA. Modeling ground penetrating radar by GPR-Max [J]. Construction and Building Materials, 2005, 19(10): 755-762

[17] LiuS-x, ZengS-f, XuBo. Subsurface fractures detection by borehole radar [J]. Journal of Jilin University: Earth Science Edition, 2005, 35(S1): 128-131

[18] WangJ, ChenW-m, ZhangPeng. Application of borehole radar to site characterization of high-level radioactive waste repository [J]. Uranium Geology, 2005, 21(6): 360-363

[19] ZhongSheng. Key issues of dynamic exploration survey based on the borehole radar and digital imaging [D]. Institute of Rock and soil Mechanics, Chinese Academy of Sciences, 200825-56

[20] ZhouJ-feng. Numerical simulation of borehole radar for metal mine detection [D]. College of Geo-exploration Science and Technology, 20091-15

[21] JamesI, RosemaryK. Numerical modeling of ground-penetrating radar in 2-D using MATLAB [J]. Computers & Geosciences, 2006, 32(9): 1247-1258

[22] YeeK. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media [J]. IEEE Transactions on Antennas and propagation, 1966, 14(3): 302-307

[23] LuebbersR, HunsbergerF. FDTD for n-th-order dispersive media [J]. IEEE Transactions on Antennas and Propagation, 1992, 40(11): 1297-1301

[24] GeorgakopoulosS V, BirtcherC R, BalanisC A. Higher-order finite-difference schemes for electromagnetic radiation, scattering, and penetration, Part 1: Theory [J]. IEEE Antennas and Propagation Magazine, 2002, 44(2): 134-142

[25] ChenJ-s, ChenC-xin. The review of borehole radar technology [J]. Progress in Geophysics, 2008, 23(5): 1634-1640