Simulation of ground penetrating radar on dispersive media by a finite element time domain algorithm

Hai Liu, Bangan Xing, Honghua Wang, Jie Cui, B F Spencer

Research output: Contribution to journalArticle

Abstract

This paper presents a two-dimensional finite element time domain (FETD) algorithm for simulation of ground penetrating radar (GPR) on dispersive and lossy media. The medium of linear dispersion is modeled by a multi-pole Debye model, which is incorporated into the FETD scheme through the Fourier transform. We apply the finite element Galerkin method to discretize the computation domain using triangular elements. A perfectly matched layer is extended to match the dispersive media and is used as an absorbing boundary condition to truncate the computation domain. The proposed FETD algorithm is proven to have comparable accuracy with the finite difference time domain (FDTD) method by comparing the simulated electromagnetic waves propagating in a homogeneous dispersive medium with the corresponding analytical solution. At an expense of computational costs, the FETD approach can yield more accurate results when simulating a complex geometry than the FDTD approach, which suffers from the staircase error. Moreover, we analyze the influence of material dispersion on GPR data acquired by different survey geometries. The results demonstrate that material dispersion will bring a great challenge for quantitative interpretation of GPR data. Numerical simulation of GPR data with dispersion does help in accurate interpretation of GPR data.

Original languageEnglish (US)
Article number103821
JournalJournal of Applied Geophysics
Volume170
DOIs
StatePublished - Nov 2019

Fingerprint

ground penetrating radar
radar data
simulation
absorbing boundary
geometry
Galerkin method
homogeneous medium
lossy media
electromagnetic wave
perfectly matched layers
stairways
Fourier transform
finite difference time domain method
boundary condition
electromagnetic radiation
poles
boundary conditions
costs
cost

Keywords

  • Debye model
  • Finite element time domain method (FETD)
  • Ground penetrating radar (GPR)
  • Medium dispersion
  • Perfect matched layer (PML)

ASJC Scopus subject areas

  • Geophysics

Cite this

Simulation of ground penetrating radar on dispersive media by a finite element time domain algorithm. / Liu, Hai; Xing, Bangan; Wang, Honghua; Cui, Jie; Spencer, B F.

In: Journal of Applied Geophysics, Vol. 170, 103821, 11.2019.

Research output: Contribution to journalArticle

@article{4ad36d74dd7543a285fd5b9f46729381,
title = "Simulation of ground penetrating radar on dispersive media by a finite element time domain algorithm",
abstract = "This paper presents a two-dimensional finite element time domain (FETD) algorithm for simulation of ground penetrating radar (GPR) on dispersive and lossy media. The medium of linear dispersion is modeled by a multi-pole Debye model, which is incorporated into the FETD scheme through the Fourier transform. We apply the finite element Galerkin method to discretize the computation domain using triangular elements. A perfectly matched layer is extended to match the dispersive media and is used as an absorbing boundary condition to truncate the computation domain. The proposed FETD algorithm is proven to have comparable accuracy with the finite difference time domain (FDTD) method by comparing the simulated electromagnetic waves propagating in a homogeneous dispersive medium with the corresponding analytical solution. At an expense of computational costs, the FETD approach can yield more accurate results when simulating a complex geometry than the FDTD approach, which suffers from the staircase error. Moreover, we analyze the influence of material dispersion on GPR data acquired by different survey geometries. The results demonstrate that material dispersion will bring a great challenge for quantitative interpretation of GPR data. Numerical simulation of GPR data with dispersion does help in accurate interpretation of GPR data.",
keywords = "Debye model, Finite element time domain method (FETD), Ground penetrating radar (GPR), Medium dispersion, Perfect matched layer (PML)",
author = "Hai Liu and Bangan Xing and Honghua Wang and Jie Cui and Spencer, {B F}",
year = "2019",
month = "11",
doi = "10.1016/j.jappgeo.2019.103821",
language = "English (US)",
volume = "170",
journal = "Journal of Applied Geophysics",
issn = "0926-9851",
publisher = "Elsevier",

}

TY - JOUR

T1 - Simulation of ground penetrating radar on dispersive media by a finite element time domain algorithm

AU - Liu, Hai

AU - Xing, Bangan

AU - Wang, Honghua

AU - Cui, Jie

AU - Spencer, B F

PY - 2019/11

Y1 - 2019/11

N2 - This paper presents a two-dimensional finite element time domain (FETD) algorithm for simulation of ground penetrating radar (GPR) on dispersive and lossy media. The medium of linear dispersion is modeled by a multi-pole Debye model, which is incorporated into the FETD scheme through the Fourier transform. We apply the finite element Galerkin method to discretize the computation domain using triangular elements. A perfectly matched layer is extended to match the dispersive media and is used as an absorbing boundary condition to truncate the computation domain. The proposed FETD algorithm is proven to have comparable accuracy with the finite difference time domain (FDTD) method by comparing the simulated electromagnetic waves propagating in a homogeneous dispersive medium with the corresponding analytical solution. At an expense of computational costs, the FETD approach can yield more accurate results when simulating a complex geometry than the FDTD approach, which suffers from the staircase error. Moreover, we analyze the influence of material dispersion on GPR data acquired by different survey geometries. The results demonstrate that material dispersion will bring a great challenge for quantitative interpretation of GPR data. Numerical simulation of GPR data with dispersion does help in accurate interpretation of GPR data.

AB - This paper presents a two-dimensional finite element time domain (FETD) algorithm for simulation of ground penetrating radar (GPR) on dispersive and lossy media. The medium of linear dispersion is modeled by a multi-pole Debye model, which is incorporated into the FETD scheme through the Fourier transform. We apply the finite element Galerkin method to discretize the computation domain using triangular elements. A perfectly matched layer is extended to match the dispersive media and is used as an absorbing boundary condition to truncate the computation domain. The proposed FETD algorithm is proven to have comparable accuracy with the finite difference time domain (FDTD) method by comparing the simulated electromagnetic waves propagating in a homogeneous dispersive medium with the corresponding analytical solution. At an expense of computational costs, the FETD approach can yield more accurate results when simulating a complex geometry than the FDTD approach, which suffers from the staircase error. Moreover, we analyze the influence of material dispersion on GPR data acquired by different survey geometries. The results demonstrate that material dispersion will bring a great challenge for quantitative interpretation of GPR data. Numerical simulation of GPR data with dispersion does help in accurate interpretation of GPR data.

KW - Debye model

KW - Finite element time domain method (FETD)

KW - Ground penetrating radar (GPR)

KW - Medium dispersion

KW - Perfect matched layer (PML)

UR - http://www.scopus.com/inward/record.url?scp=85070543874&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85070543874&partnerID=8YFLogxK

U2 - 10.1016/j.jappgeo.2019.103821

DO - 10.1016/j.jappgeo.2019.103821

M3 - Article

AN - SCOPUS:85070543874

VL - 170

JO - Journal of Applied Geophysics

JF - Journal of Applied Geophysics

SN - 0926-9851

M1 - 103821

ER -