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) pr operties. Magnetic and dielectric relaxations are described by relaxation f unctions associated with each principal component of the respective tensori al property. Anisotropy is modeled up to orthorhombic symmetry, i.e., the p rincipal 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 deri vatives and a second-order finite difference in time. A complete plane-wave analysis, including energy balance, gives the expressions of measurable qu antities such as the EM-wave velocity and the quality factor as a function of frequency and propagation direction. The algorithm reproduces the wavefr ont shape and attenuation predicted by the plane-wave analysis. In addition , the results are in excellent agreement with an analytical 3-D transient s olution. The modeling is applied to evaluation of the EM response of two pi pes buried in an anisotropic sand overlying a clay layer. The results demon strate that anisotropy and loss mechanisms cause significant amplitude and phase distortions.