GPR modelling by the Fourier method: improvement of the algorithm

We improve two aspects of the modelling scheme for the simulation of electromagnetic radio waves, based on the Fourier pseudospectral method. When there are large contrasts in the material properties, use of the standard algorithm (regular grid) causes a series of artefacts, as, for instance, ringing and acausal events. These problems, due to the non-locality of the differential operator, are solved by using the staggered Fourier method (staggered grid). Realistic radiation patterns can be obtained from simple combinations of magnetic and electric sources. If the directivity pattern of the antenna is known, from either a finite-difference simulation or an analytic evaluation or an experimental characterization, it can then be simulated by a composite-source concept. This effective source is implemented in the modelling algorithm by means of a perturbation technique, which first computes the intensity and directional spectra of the single electromagnetic sources. Their location is optimized to obtain the best fit with a minimum number of sources. The approach is, in principle, valid for the far-field radiation pattern of the antenna.

[1]  Waymond R. Scott,et al.  Accurate computation of the radiation from simple antennas using the finite-difference time-domain method , 1989 .

[2]  Charles Elachi,et al.  Radiation patterns of interfacial dipole antennas , 1982 .

[3]  E. L. Harder,et al.  The Institute of Electrical and Electronics Engineers, Inc. , 2019, 2019 IEEE International Conference on Software Architecture Companion (ICSA-C).

[4]  Roger F. Harrington,et al.  Field computation by moment methods , 1968 .

[5]  Greg Turner.  Modelling antenna-ground interactions , 1994 .

[6]  A. P. Annan,et al.  Radio Interferometry Depth Sounding: Part I—THEORETICAL Discussion , 1973 .

[7]  Steven A. Arcone,et al.  Numerical studies of the radiation patterns of resistively loaded dipoles , 1995 .

[8]  J. Bourgeois.  A Fully Three-Dimensional Simulation of a Ground-Penetrating Radar: FDT with Experiment Theory Compared , 1996 .

[9]  Glenn S. Smith,et al.  A fully three-dimensional simulation of a ground-penetrating radar: FDTD theory compared with experiment , 1996, IEEE Trans. Geosci. Remote. Sens..

[10]  B. Fornberg.  High-order finite differences and the pseudospectral method on staggered grids , 1990 .

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

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

[13]  José M. Carcione.  Ground‐penetrating radar: Wave theory and numerical simulation in lossy anisotropic media , 1995 .

[14]  Salvatore Piro,et al.  Field experiments for characterization of GPR antenna and pulse propagation , 1995 .

[15]  J. Robertsson,et al.  Finite‐difference modeling of electromagnetic wave propagation in dispersive and attenuating media , 1998 .

[16]  J. Deen,et al.  MEASURED UNDERWATER NEAR‐FIELD E‐PATTERNS OF A PULSED, HORIZONTAL DIPOLE ANTENNA IN AIR: COMPARISON WITH THE THEORY OF THE CONTINUOUS WAVE, INFINITESIMAL ELECTRIC DIPOLE1 , 1990 .

[17]  G. McMechan,et al.  Causes and reduction of numerical artefacts in pseudo-spectral wavefield extrapolation , 1996 .

[18]  M. Landrø,et al.  Implementing measured source signatures in a coarse-grid, finite-difference modeling scheme , 1993 .

[19]  George A. McMechan,et al.  Algorithms for staggered‐grid computations for poroelastic, elastic, acoustic, and scalar wave equations , 1997 .

[20]  Glenn S. Smith.  Directive properties of antennas for transmission into a material half-space , 1984 .

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