Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform cross-hole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4–10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

[1]  M. H. Waxman,et al.  Electrical Conductivities in Oil-Bearing Shaly Sands , 1968 .

[2]  A. P. Annan,et al.  Electromagnetic determination of soil water content: Measurements in coaxial transmission lines , 1980 .

[3]  N. Bleistein Two-and-One-Half Dimensional In-Plane Wave Propagation. , 1984 .

[4]  A. Tarantola Inversion of seismic reflection data in the acoustic approximation , 1984 .

[5]  A. Tarantola LINEARIZED INVERSION OF SEISMIC REFLECTION DATA , 1984 .

[6]  A. Tarantola A strategy for nonlinear elastic inversion of seismic reflection data , 1986 .

[7]  P. Mora Nonlinear two-dimensional elastic inversion of multioffset seismic data , 1987 .

[8]  R. Knight,et al.  The dielectric constant of sandstones, 60 kHz to 4 MHz , 1987 .

[9]  Paul F. Worthington,et al.  The uses and abuses of the Archie equations, 1: The formation factor-porosity relationship , 1993 .

[10]  Z. Bing,et al.  Crosshole acoustic velocity imaging with full-waveform spectral data: 2.5-D numerical simulations , 1998 .

[11]  R. Shipp,et al.  Seismic waveform inversion in the frequency domain, Part 2: Fault delineation in sediments using crosshole data , 1999 .

[12]  R. Pratt Seismic waveform inversion in the frequency domain; Part 1, Theory and verification in a physical scale model , 1999 .

[13]  Müller,et al.  GPR study of pore water content and salinity in sand , 2000 .

[14]  H. Maurer,et al.  Ray-based amplitude tomography for crosshole georadar data: a numerical assessment , 2001 .

[15]  A. Binley,et al.  Vadose zone flow model parameterisation using cross-borehole radar and resistivity imaging , 2001 .

[16]  Herve Perroud,et al.  On the use of combined geophysical methods to assess water content and water conductivity of near-surface formations , 2002 .

[17]  P. Dietrich,et al.  Integrating surface georadar and crosshole radar tomography: A validation experiment in braided stream deposits , 2002 .

[18]  W. Barrash,et al.  Hierarchical geostatistics and multifacies systems: Boise Hydrogeophysical Research Site, Boise, Idaho , 2002 .

[19]  A. Binley,et al.  Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles. , 2002 .

[20]  S. Greenhalgh,et al.  Crosshole seismic inversion with normalized full‐waveform amplitude data , 2003 .

[21]  Martin Musil,et al.  Effects and removal of systematic errors in crosshole georadar attenuation tomography , 2004 .

[22]  K. Nihei,et al.  Viscoacoustic wave form inversion of transmission data for velocity and attenuation , 2004 .

[23]  Jens Tronicke,et al.  Effects of gas- and water-filled boreholes on the amplitudes of crosshole georadar data as inferred from experimental evidence , 2004 .

[24]  Barry J. Allred,et al.  Hydrogeophysical Case Studies in the Vadose Zone , 2005 .

[25]  Anita Turesson,et al.  Water content and porosity estimated from ground-penetrating radar and resistivity , 2006 .

[26]  D. Gomez,et al.  Flow and transport in the unsaturated Sherwood Sandstone: characterization using cross-borehole geophysical methods , 2006, Geological Society, London, Special Publications.

[27]  A. Binley,et al.  Improved hydrogeophysical characterization using joint inversion of cross‐hole electrical resistance and ground‐penetrating radar traveltime data , 2006 .

[28]  Jan Van der Kruk,et al.  Accurate Imaging of Multicomponent GPR Data Based on Exact Radiation Patterns , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[29]  Jacques R. Ernst,et al.  Application of a new 2D time-domain full-waveform inversion scheme to crosshole radar data , 2007 .

[30]  Seiichiro Kuroda,et al.  Full-waveform inversion algorithm for interpreting crosshole radar data: a theoretical approach , 2007 .

[31]  Jacques R. Ernst,et al.  Full-Waveform Inversion of Crosshole Radar Data Based on 2-D Finite-Difference Time-Domain Solutions of Maxwell's Equations , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[32]  Markus Hofer,et al.  Analyzing Bank Filtration by Deconvoluting Time Series of Electric Conductivity , 2007, Ground water.

[33]  Jan van der Kruk,et al.  Characterizing a GPR antenna system by near-field electric field measurements , 2007 .

[34]  Alan G. Green,et al.  Vector-migration of standard copolarized 3D GPR data , 2007 .

[35]  Andrew Binley,et al.  Monitoring Unsaturated Flow and Transport Using Cross‐Borehole Geophysical Methods , 2008 .

[36]  Niklas Linde,et al.  The borehole-fluid effect in electrical resistivity imaging , 2010 .

[37]  Samuel Diem,et al.  Räumliche Charakterisierung der hydraulischen Leitfähigkeit in alluvialen Schotter-Grundwasserleitern: Ein Methodenvergleich , 2010 .

[38]  S. Greenhalgh,et al.  Zonation for 3D aquifer characterization based on joint inversions of multimethod crosshole geophysical data , 2010 .

[39]  Jacques R. Ernst,et al.  A New Vector Waveform Inversion Algorithm for Simultaneous Updating of Conductivity and Permittivity Parameters From Combination Crosshole/Borehole-to-Surface GPR Data , 2010, IEEE Transactions on Geoscience and Remote Sensing.