3-D ground‐penetrating radar simulation and plane‐wave theory in anisotropic media

Modeling ground‐penetrating radar (GPR) waves requires simulation of the 3-D full wavefield and the correct description of the electromagnetic (EM) properties. Magnetic and dielectric relaxations are described by relaxation functions associated with each principal component of the respective tensorial property. Anisotropy is modeled up to orthorhombic symmetry, i.e., the principal coordinate systems of the three EM material tensors coincide, and each property is described by three different principal components. The algorithm uses the pseudospectral method for computing the spatial derivatives and a second‐order finite difference in time. A complete plane‐wave analysis, including energy balance, gives the expressions of measurable quantities such as the EM-wave velocity and the quality factor as a function of frequency and propagation direction. The algorithm reproduces the wavefront shape and attenuation predicted by the plane‐wave analysis. In addition, the results are in excellent agreement with an ana...

[1]  Edip Baysal,et al.  Forward modeling by a Fourier method , 1982 .

[2]  George A. McMechan,et al.  Acquisition and processing of wide-aperture ground-penetrating radar data , 1992 .

[3]  P. Saraf,et al.  Anisotropy in Geoelectromagnetism , 1990 .

[4]  Gary R. Olhoeft,et al.  Petrophysical causes of electromagnetic dispersion , 1994 .

[5]  A. Taflove.  The Finite-Difference Time-Domain Method , 1995 .

[6]  José M. Carcione,et al.  Ground‐penetrating radar: Wave theory and numerical simulation in lossy anisotropic media , 1996 .

[7]  J. Daniels,et al.  Modeling near-field GPR in three dimensions using the FDTD method , 1997 .

[8]  George A. McMechan,et al.  GPR attenuation and its numerical simulation in 2.5 dimensions , 1997 .

[9]  M. Carcione,et al.  Wave propagation simulation in a linear viscoacoustic medium , 1997 .

[10]  Peter G. Petropoulos,et al.  The wave hierarchy for propagation in relaxing dielectrics , 1995 .

[11]  Imaging Pipelines In 3-D By Ground-Penetrating Radar , 1992 .

[12]  S. Tillard.  Radar experiments in isotropic and anisotropic geological formations (granite and schists)1 , 1994 .

[13]  George A. McMechan,et al.  Comparison of ray and Fourier methods for modeling monostatic ground‐penetrating radar profiles , 1995 .

[14]  J. Carcione.  Radiation patterns for 2-D GPR forward modeling , 1998 .

[15]  Géza Seriani,et al.  A spectral scheme for wave propagation simulation in 3-D elastic-anisotropic media , 1992 .

[16]  Gilles Grandjean,et al.  GPR data processing for 3D fracture mapping in a marble quarry (Thassos, Greece) , 1996 .

[17]  J. Carcione,et al.  Generalized Mechanical Model Analogies of Linear Viscoelastic Behaviour , 1992 .

[18]  David A. Casper,et al.  Simulation of ground-penetrating radar waves in a 2-D soil model , 1996 .

[19]  W. Chew.  Waves and Fields in Inhomogeneous Media , 1990 .

[20]  Moshe Reshef,et al.  Three-dimensional elastic modeling by the Fourier method , 1988 .

[21]  J. Carcione.  Ground radar simulation for archaeological applications1 , 1996 .

[22]  Mark Grasmueck,et al.  3-D ground‐penetrating radar applied to fracture imaging in gneiss , 1996 .