Finite difference time domain method forward simulation of complex geoelectricity ground penetrating radar model

Abstract

The ground penetrating radar (GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.

Keywords

ground penetrating radar / finite difference time domain method / forward simulation / ideal frequency dispersion relationship

References

[1] LambotS, SlobE C, VanB, et al.. Modeling of ground-penetrating radar for accurate characterization of subsurface electric properties [J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(11): 2555-2568

[2] GuangyouF. FDTD and optimization approach to time-domain inversing problem for underground multiple objects [J]. Microwave and Optical Technology Letters, 2001, 31(5): 384-387

[3] MeinckeP. Linear GPR inversion for lossy soil and a planar air-soil interface[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(12): 2713-2721

[4] KowalskyM B, DietrichP, TeutschG, et al.. Forward modeling of ground-penetrating radar data using digitized outcrop images and multiple scenarios of water saturation[J]. Water Resources Research, 2001, 37(6): 1615-1625

[5] EisenO, WilhelmsF, NixdorfU, et al.. Revealing the nature of radar reflections in ice: DEP-based FDTD forward modeling[J]. Geophysical Research Letters, 2003, 30(5): 1218-1224

[6] HubbardS, GroteK. Mapping the volumetric soil water content of a California vineyard using high-frequency GPR ground wave data[J]. Society of Exploration Geophysicists, 2002, 21(6): 556-559

[7] CullJ, MassieD, RobertsJ. Complex impedance mapping using GPR survey methods[A]. International Geoscience and Remote Sensing Symposium[C], 2003, Toulouse, IEEE: 2909-2911

[8] HasegawaY, YokoeK, KawaiY, et al.. GPR-based adaptive sensing GPR manipulation according to terrain configurations [A]. International Conference on Intelligent Robots and Systems [C], 2004, Sendai, IEEE: 3021-3026

[9] 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

[10] CarcioneJ M. Radiation patterns for 2-D GPR forward modeling[J]. Geophysics, 1998, 63(2): 424-430

[11] BernardiP, CavagnaroM, AtanasioP, et al.. FDTD, multiple-region/FDTD, ray-tracing/FDTD: A comparison on their applicability for human exposure evaluation [J]. International Journal of Numerical Modelling, 2002, 15(5–6): 579-593

[12] RejibaF, CamerlynckC, MechlerP. FDTD-SUPMLADE simulation for ground-penetrating radar modeling [J]. Radio Science, 2003, 38(1): 511-513

[13] YarovoyA G, VazourasC N, FikiorisJ G, et al.. Numerical simulations of the scattered field near a statistically rough air-ground interface [J]. IEEE Transactions on Antennas and Propagation, 2004, 52(3): 780-789

[14] MeiK K, FangJ. Superabsorption—a method to improve absorbing boundary conditions [J]. IEEE Transactions on Antennas and Propagation, 1992, 40(9): 1001-1010

[15] CuiT J, ChewW C, AydinerA A, et al.. Fast-for-ward solvers for the low-frequency detection of buried dielectric objects[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(9): 2026-2036

[16] JoaquimF G. A novel 3-D subsurface radar imaging technique[J]. IEEE Transactions on Geoscience and Remote Sensing, 2002, 40(2): 443-452