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

详细信息    查看全文

• 作者：James Irving and Rosemary Knight
• 关键词：Ground-penetrating radar ; GPR ; Numerical modeling ; Electromagnetics ; Finite-difference time domain ; FDTD
• 刊名：Computers and Geosciences
• 出版年：2006
• 出版时间：November 2006
• 年：2006
• 卷：32
• 期：9
• 页码：1247-1258
• 全文大小：546 K

文摘

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.