Towards 3D full-waveform inversion of crosshole GPR data

A. Mozaffari, A. Klotzsche, Guowei He, H. Vereecken, J. Van Der Kruk, C. Warren, A. Giannopoulos

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

2D crosshole ground penetrating (GPR) full-1 waveform inversion (FWI) has shown superior resolution 1 compared to ray-based inversion tomograms for synthetic and experimental data. To invert measured 3D data with a 2D j model that has a reduced geometrical spreading and assumes infinite source and receiver dimensions perpendicular to the 2D inversion plane, the Bleistein transformation can be used to convert the measured 3D data to 2D. This far-field conversion consists of a phase shift and amplitude correction that is based on the first arrival travel time of each trace. In the case of late arrival and high amplitude events that 1 can occur due to low-velocity waveguides andfast propagating refracted waves, this transformation can introduce errors especially in the amplitude such that the inverted conductivity is less accurate. To overcome these problems, we have replaced the 2D finite difference time domain FDTD forward model in the FWI scheme with the well-known gprMax3D FDTD modeling program. In this way, we do not need to use the Bleistein 3D-to-2D filter with its far-field approximation, and we can deal with the correct geometrical spreading and approximate better realistic point source and receivers. Currently, the 2D FWI algorithm has been extended to 2.5D by replacing the 2D FDTD with the gprMax 3D FDTD modeling program. The first test of the 2.5D FWI is based on inversion results of a data set acquired at the Widen site in Switzerland that contained a low-velocity waveguide. The new 2.5D FWI performed well in reconstructing the models. The new 2.5D FWI enables a more reliable re-construction of the subsurface image, especially the electrical conductivity tomograms and to the full use of all the modeling possibilities of the gprMax modeling tool.

Original languageEnglish (US)
Title of host publicationProceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016
PublisherInstitute of Electrical and Electronics Engineers Inc.
ISBN (Electronic)9781509051816
DOIs
StatePublished - Sep 20 2016
Event16th International Conference of Ground Penetrating Radar, GPR 2016 - Hong Kong, Hong Kong
Duration: Jun 13 2016Jun 16 2016

Other

Other16th International Conference of Ground Penetrating Radar, GPR 2016
CountryHong Kong
CityHong Kong
Period6/13/166/16/16

Fingerprint

waveforms
inversions
Waveguides
finite difference time domain method
Travel time
Phase shift
low speed
arrivals
far fields
receivers
refracted waves
waveguides
Switzerland
point sources
travel
rays
phase shift
filters
conductivity
electrical resistivity

Keywords

  • 3D to 2D transformation
  • Cross-hole GPR
  • EM waves
  • Full-waveform inversion
  • Wave propagation

ASJC Scopus subject areas

  • Signal Processing
  • Instrumentation

Cite this

Mozaffari, A., Klotzsche, A., He, G., Vereecken, H., Van Der Kruk, J., Warren, C., & Giannopoulos, A. (2016). Towards 3D full-waveform inversion of crosshole GPR data. In Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016 [7572687] Institute of Electrical and Electronics Engineers Inc.. https://doi.org/10.1109/ICGPR.2016.7572687

Towards 3D full-waveform inversion of crosshole GPR data. / Mozaffari, A.; Klotzsche, A.; He, Guowei; Vereecken, H.; Van Der Kruk, J.; Warren, C.; Giannopoulos, A.

Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016. Institute of Electrical and Electronics Engineers Inc., 2016. 7572687.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Mozaffari, A, Klotzsche, A, He, G, Vereecken, H, Van Der Kruk, J, Warren, C & Giannopoulos, A 2016, Towards 3D full-waveform inversion of crosshole GPR data. in Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016., 7572687, Institute of Electrical and Electronics Engineers Inc., 16th International Conference of Ground Penetrating Radar, GPR 2016, Hong Kong, Hong Kong, 6/13/16. https://doi.org/10.1109/ICGPR.2016.7572687
Mozaffari A, Klotzsche A, He G, Vereecken H, Van Der Kruk J, Warren C et al. Towards 3D full-waveform inversion of crosshole GPR data. In Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016. Institute of Electrical and Electronics Engineers Inc. 2016. 7572687 https://doi.org/10.1109/ICGPR.2016.7572687
Mozaffari, A. ; Klotzsche, A. ; He, Guowei ; Vereecken, H. ; Van Der Kruk, J. ; Warren, C. ; Giannopoulos, A. / Towards 3D full-waveform inversion of crosshole GPR data. Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016. Institute of Electrical and Electronics Engineers Inc., 2016.
@inproceedings{16354f07eeab4ad1a417713979e3d34f,
title = "Towards 3D full-waveform inversion of crosshole GPR data",
abstract = "2D crosshole ground penetrating (GPR) full-1 waveform inversion (FWI) has shown superior resolution 1 compared to ray-based inversion tomograms for synthetic and experimental data. To invert measured 3D data with a 2D j model that has a reduced geometrical spreading and assumes infinite source and receiver dimensions perpendicular to the 2D inversion plane, the Bleistein transformation can be used to convert the measured 3D data to 2D. This far-field conversion consists of a phase shift and amplitude correction that is based on the first arrival travel time of each trace. In the case of late arrival and high amplitude events that 1 can occur due to low-velocity waveguides andfast propagating refracted waves, this transformation can introduce errors especially in the amplitude such that the inverted conductivity is less accurate. To overcome these problems, we have replaced the 2D finite difference time domain FDTD forward model in the FWI scheme with the well-known gprMax3D FDTD modeling program. In this way, we do not need to use the Bleistein 3D-to-2D filter with its far-field approximation, and we can deal with the correct geometrical spreading and approximate better realistic point source and receivers. Currently, the 2D FWI algorithm has been extended to 2.5D by replacing the 2D FDTD with the gprMax 3D FDTD modeling program. The first test of the 2.5D FWI is based on inversion results of a data set acquired at the Widen site in Switzerland that contained a low-velocity waveguide. The new 2.5D FWI performed well in reconstructing the models. The new 2.5D FWI enables a more reliable re-construction of the subsurface image, especially the electrical conductivity tomograms and to the full use of all the modeling possibilities of the gprMax modeling tool.",
keywords = "3D to 2D transformation, Cross-hole GPR, EM waves, Full-waveform inversion, Wave propagation",
author = "A. Mozaffari and A. Klotzsche and Guowei He and H. Vereecken and {Van Der Kruk}, J. and C. Warren and A. Giannopoulos",
year = "2016",
month = "9",
day = "20",
doi = "10.1109/ICGPR.2016.7572687",
language = "English (US)",
booktitle = "Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016",
publisher = "Institute of Electrical and Electronics Engineers Inc.",

}

TY - GEN

T1 - Towards 3D full-waveform inversion of crosshole GPR data

AU - Mozaffari, A.

AU - Klotzsche, A.

AU - He, Guowei

AU - Vereecken, H.

AU - Van Der Kruk, J.

AU - Warren, C.

AU - Giannopoulos, A.

PY - 2016/9/20

Y1 - 2016/9/20

N2 - 2D crosshole ground penetrating (GPR) full-1 waveform inversion (FWI) has shown superior resolution 1 compared to ray-based inversion tomograms for synthetic and experimental data. To invert measured 3D data with a 2D j model that has a reduced geometrical spreading and assumes infinite source and receiver dimensions perpendicular to the 2D inversion plane, the Bleistein transformation can be used to convert the measured 3D data to 2D. This far-field conversion consists of a phase shift and amplitude correction that is based on the first arrival travel time of each trace. In the case of late arrival and high amplitude events that 1 can occur due to low-velocity waveguides andfast propagating refracted waves, this transformation can introduce errors especially in the amplitude such that the inverted conductivity is less accurate. To overcome these problems, we have replaced the 2D finite difference time domain FDTD forward model in the FWI scheme with the well-known gprMax3D FDTD modeling program. In this way, we do not need to use the Bleistein 3D-to-2D filter with its far-field approximation, and we can deal with the correct geometrical spreading and approximate better realistic point source and receivers. Currently, the 2D FWI algorithm has been extended to 2.5D by replacing the 2D FDTD with the gprMax 3D FDTD modeling program. The first test of the 2.5D FWI is based on inversion results of a data set acquired at the Widen site in Switzerland that contained a low-velocity waveguide. The new 2.5D FWI performed well in reconstructing the models. The new 2.5D FWI enables a more reliable re-construction of the subsurface image, especially the electrical conductivity tomograms and to the full use of all the modeling possibilities of the gprMax modeling tool.

AB - 2D crosshole ground penetrating (GPR) full-1 waveform inversion (FWI) has shown superior resolution 1 compared to ray-based inversion tomograms for synthetic and experimental data. To invert measured 3D data with a 2D j model that has a reduced geometrical spreading and assumes infinite source and receiver dimensions perpendicular to the 2D inversion plane, the Bleistein transformation can be used to convert the measured 3D data to 2D. This far-field conversion consists of a phase shift and amplitude correction that is based on the first arrival travel time of each trace. In the case of late arrival and high amplitude events that 1 can occur due to low-velocity waveguides andfast propagating refracted waves, this transformation can introduce errors especially in the amplitude such that the inverted conductivity is less accurate. To overcome these problems, we have replaced the 2D finite difference time domain FDTD forward model in the FWI scheme with the well-known gprMax3D FDTD modeling program. In this way, we do not need to use the Bleistein 3D-to-2D filter with its far-field approximation, and we can deal with the correct geometrical spreading and approximate better realistic point source and receivers. Currently, the 2D FWI algorithm has been extended to 2.5D by replacing the 2D FDTD with the gprMax 3D FDTD modeling program. The first test of the 2.5D FWI is based on inversion results of a data set acquired at the Widen site in Switzerland that contained a low-velocity waveguide. The new 2.5D FWI performed well in reconstructing the models. The new 2.5D FWI enables a more reliable re-construction of the subsurface image, especially the electrical conductivity tomograms and to the full use of all the modeling possibilities of the gprMax modeling tool.

KW - 3D to 2D transformation

KW - Cross-hole GPR

KW - EM waves

KW - Full-waveform inversion

KW - Wave propagation

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

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

U2 - 10.1109/ICGPR.2016.7572687

DO - 10.1109/ICGPR.2016.7572687

M3 - Conference contribution

BT - Proceedings of 2016 16th International Conference of Ground Penetrating Radar, GPR 2016

PB - Institute of Electrical and Electronics Engineers Inc.

ER -