3-D ground-penetrating radar simulation: Design of the modeling algorithm
José M. Carcione and Michael Schoenberg
Event name: 5th International Congress of the Brazilian Geophysical Society
Session: Use of Ground Penetrating Radar - GPR
Publication date: 07 November 1997
Info: Extended abstract, PDF ( 1.71Mb )
Modeling ground-penetrating radar (GPR) waves requires the simulation of the 3-D full wavefield and the correct description of the different electromagnetic properties. Magnetic and dielectric relaxation and deviations from Ohm's lay are described by viscoelastic-type relaxation functions associated to each principal component of the respective tensorial property. Anisotropy is modeled up to orthorhombic symmetry, i.e., the principal coordinate systems of the three electromagnetic material tensors coincide, and each property is described by three different principal components. The algorithm scheme is based on the Fourier pseudospectral method for computing the spatial derivatives and a staggered second-order time stepping algorithm. A complete plane wave analysis, including energy balance, gives the expressions of measurable quantities such as the electromagentic velocity and the quality factor as a function of frequency and propagation direction. The results are preliminary and verify the modeling results by comparing snapshots with plane wave theory.