Numerical modeling of ground-penetrating radar in 2-D using MATLAB

被引:162
作者
Irving, James [1 ]
Knight, Rosemary [1 ]
机构
[1] Stanford Univ, Dept Geophys, Stanford, CA 94305 USA
基金
美国国家科学基金会;
关键词
ground-penetrating radar; GPR; numerical modeling; electromagnetics; finite-difference time domain; FDTD;
D O I
10.1016/j.cageo.2005.11.006
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We present MATLAB codes for finite-difference time-domain (FDTD) modeling of ground-penetrating radar (GPR) in two dimensions. Surface-based reflection GPR is modeled using a transverse magnetic (TM-) mode formulation. Crosshole and vertical radar profiling (VRP) geometries are modeled using a transverse electric (TE-) mode formulation. Matrix notation is used in the codes wherever possible to optimize them for speed in the MATLAB environment. To absorb waves at the edges of the modeling grid, we implement perfectly matched layer (PML) absorbing boundaries. Although our codes are two-dimensional and do not incorporate features such as dispersion in electrical properties, they capture many of the important elements of GPR surveying and run at a fraction of the computational cost of more elaborate algorithms. In addition, the codes are well commented, relatively easy to understand, and can be easily modified for the user's specific purpose. (c) 2005 Elsevier Ltd. All rights reserved.
引用
收藏
页码:1247 / 1258
页数:12
相关论文
共 30 条
[1]   A PERFECTLY MATCHED LAYER FOR THE ABSORPTION OF ELECTROMAGNETIC-WAVES [J].
BERENGER, JP .
JOURNAL OF COMPUTATIONAL PHYSICS, 1994, 114 (02) :185-200
[2]   Finite-difference modeling of electromagnetic wave propagation in dispersive and attenuating media [J].
Bergmann, T ;
Robertsson, JOA ;
Holliger, K .
GEOPHYSICS, 1998, 63 (03) :856-867
[3]   Numerical properties of staggered finite-difference solutions of Maxwell's equations for ground-penetrating radar modeling [J].
Bergmann, T ;
Robertsson, JOA ;
Holliger, K .
GEOPHYSICAL RESEARCH LETTERS, 1996, 23 (01) :45-48
[4]   A fully three-dimensional simulation of a ground-penetrating radar: FDTD theory compared with experiment [J].
Bourgeois, JM ;
Smith, GS .
IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, 1996, 34 (01) :36-44
[5]   RAY-BASED SYNTHESIS OF BISTATIC GROUND-PENETRATING RADAR PROFILES [J].
CAI, J ;
MCMECHAN, GA .
GEOPHYSICS, 1995, 60 (01) :87-96
[6]   Ground-penetrating radar: Wave theory and numerical simulation in lossy anisotropic media [J].
Carcione, JM .
GEOPHYSICS, 1996, 61 (06) :1664-1677
[7]   Simulation of ground-penetrating radar waves in a 2-D soil model [J].
Casper, DA ;
Kung, KJS .
GEOPHYSICS, 1996, 61 (04) :1034-1049
[8]   Application of perfectly matched layers to the transient modeling of subsurface EM problems [J].
Chen, YH ;
Chew, WC ;
Oristaglio, ML .
GEOPHYSICS, 1997, 62 (06) :1730-1736
[9]   A 3D PERFECTLY MATCHED MEDIUM FROM MODIFIED MAXWELLS EQUATIONS WITH STRETCHED COORDINATES [J].
CHEW, WC ;
WEEDON, WH .
MICROWAVE AND OPTICAL TECHNOLOGY LETTERS, 1994, 7 (13) :599-604
[10]  
Deutsch C. V., 1992, GSLIB GEOSTATISTICAL