1887
  1. Home
  2. A-Z Publications
  3. Near Surface Geophysics
  4. Volume 8, Issue 6
  5. Article
Volume 8, Issue 6
  • ISSN: 1569-4445
  • E-ISSN: 1873-0604

Abstract

ABSTRACT

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.

© European Association of Geoscientists and Engineers © 2010 European Association of Geoscientists and Engineers
Loading

Article metrics loading...

/content/journals/10.3997/1873-0604.2010054
2010-08-01
2024-04-26

  • From This Site

    /content/journals/10.3997/1873-0604.2010054
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. ArchieG.E.1942. The electrical resistivity log as an aid in determining some reservoir characteristics. Transactions of the American Institute of Mining and Metallurgical Engineers/Petroleum Division146, 54–62.
     [Google Scholar]
  2. BarrashW. and ClemoT.2002. Hierarchical geostatistics and multifacies systems: Boise Hydrogeophysical Research Site, Boise, Idaho. Water Resources Research38, 1196.
     [Google Scholar]
  3. BinleyA., CassianiG., MiddletonR. and WinshipP.2002a. Vadose zone flow model parameterisation using cross‐borehole radar and resistivity imaging. Journal of Hydrology267, 147–159.
     [Google Scholar]
  4. BinleyA., WinshipP., WestL.J., PokarM. and MiddletonR.2002b. Seasonal variation of moisture content in unsaturated sandstone inferred from borehole radar and resistivity profiles. Journal of Hydrology267, 160–172.
     [Google Scholar]
  5. BleisteinN.1986. 2‐1/2 dimensional in‐plane wave‐propagation. Geophysical Prospecting34, 686–703.
     [Google Scholar]
  6. CirpkaO.A., FienenM.N., HoferM., HoehnE., TessariniA., KipferR. and KitanidisP.K.2007. Analyzing bank filtration by deconvolving time series of electric conductivity. Ground Water45, 318–328.
     [Google Scholar]
  7. DanielsJ.J., AllredB., BinleyA., LabrecqueD. and AlumbaughD.2005. Hydrogeophysical case studies in the vadose zone. In: Hydrogeophysics (eds Y.Rubin and S.S.Hubbard ), pp. 413–440. Springer.
     [Google Scholar]
  8. DiemS., VogtT. and HoehnE.2010. Räumliche Charakterisierung der hydraulischen Leitfähigkeit in alluvialen Schotter‐Grundwasserleitern: Ein Methodenvergleich. Grundwasser (in press).
     [Google Scholar]
  9. DoetschJ.A., CosciaI., GreenhalghS., LindeN., GreenA. and GüntherT.2010a. The borehole‐fluid effect in electrical resistivity imaging. Geophysics75, F107–F114.
     [Google Scholar]
  10. DoetschJ., LindeN., CosciaI., GreenhalghS. and GreenA.2010b. Zonation for 3D aquifer characterization based on joint inversion of multi‐method crosshole geophysical data. Geophysics75 (in press).
     [Google Scholar]
  11. DoetschJ., LindeN. and GreenA.G.2009. Joint inversion improves zonation for aquifer characterization. EGU meeting, 19–24 April, Vienna, Austria, Expanded Abstracts, EGU 2009–2303.
     [Google Scholar]
  12. ErnstJ.R.2007. 2‐D finite‐difference time‐domain full‐waveform inversion of crosshole georadar data. PhD thesis, Swiss Federal Institute of Technology Zurich.
     [Google Scholar]
  13. ErnstJ.R., MaurerH., GreenA.G. and HolligerK.2007a. Full‐waveform inversion of crosshole radar data based on 2‐D finite‐difference time‐domain solutions of Maxwells equations. IEEE Transactions on Geoscience and Remote Sensing45, 2807–2828.
     [Google Scholar]
  14. ErnstJ.R., MaurerH., GreenA.G. and HolligerK.2007b. Application of a new 2D time‐domain full‐waveform inversion scheme to crosshole radar data. Geophysics72, J53–J64.
     [Google Scholar]
  15. FuchtbauerH.1988. Sedimente und Sedimentgesteine: Sandsteine, 4th edn. Schweizerbart Stuttgart.
     [Google Scholar]
  16. GaramboisS., SenechalP. and PerroudH.2002. On the use of combined geophysical methods to assess water content and water conductivity of near‐surface formations. Journal of Hydrology259, 32–48.
     [Google Scholar]
  17. HagreyS.A. and MuellerC.2000. GPR study of pore water content and salinity in sand. Geophysical Prospecting48, 63–85.
     [Google Scholar]
  18. HolligerK., MusilM. and MaurerH.2001. Ray‐based amplitude tomography for crosshole georadar data: A numerical assessment. Journal of Applied Geophysics47, 285–298.
     [Google Scholar]
  19. KnightR.J. and NurA.1987. The dielectric constant of sandstones, 60 kHz to 4 MHz. Geophysics52, 644–654.
     [Google Scholar]
  20. KurodaS., TakeuchiM. and KimH.J.2007. Full‐waveform inversion algorithm for interpreting crosshole radar data: A theoretical approach. Korean Geosciences Journal11, 211–217.
     [Google Scholar]
  21. LindeN., BinleyA., TryggvasonA., PedersenL.B. and RevilA.2006. Improved hydrogeophysical characterization using joint inversion of cross‐hole electrical resistance and ground‐penetrating radar travel‐time data. Water Resources Research42, W12404.
     [Google Scholar]
  22. LoomsM.C., JensenK.H., BinleyA. and NielsenL.2008. Monitoring unsaturated flow and transport using cross‐borehole geophysical methods. Vadose Zone Journal7, 227–237.
     [Google Scholar]
  23. MaurerH. and MusilM.2004. Effects and removal of systematic errors in crosshole georadar attenuation tomography. Journal of Applied Geophysics55, 261–270.
     [Google Scholar]
  24. MelesA.G., van der KrukJ., GreenhalghS.A., ErnstJ.R., MaurerH. and GreenA.G.2010. A new vector waveform inversion algorithm for simultaneous updating of conductivity and permittivity parameters from combination crosshole (borehole‐to‐surface GPR data. IEEE Transactions Geoscience and Remote Sensing48, 3391–3407.
     [Google Scholar]
  25. MoraP.1987. Nonlinear two‐dimensional elastic inversion of multi‐offset seismic data. Geophysics52, 1211–1228.
     [Google Scholar]
  26. PrattR.G.1999. Seismic waveform inversion in the frequency domain. Part 1: Theory and verification in a physical scale model. Geophysics64, 888–901.
     [Google Scholar]
  27. PrattR.G. and ShippR.M.1999. Seismic waveform inversion in the frequency domain, part 2: Fault delineation in sediments using crosshole data. Geophysics64, 902–914.
     [Google Scholar]
  28. StreichR. and van der KrukJ.2007a. Accurate imaging of multicomponent GPR data based on exact radiation patterns. IEEE Transactions on Geoscience and Remote Sensing45, 93–103.
     [Google Scholar]
  29. StreichR. and van der KrukJ.2007b. Characterization a GPR‐antenna system by near‐field electrical field measurements. Geophysics72, A51–A55.
     [Google Scholar]
  30. StreichR., van der KrukJ. and GreenA.G.2007. Vector‐migration of standard copolarized 3D GPR data. Geophysics72, 65–75.
     [Google Scholar]
  31. TarantolaA.1984a. Inversion of seismic‐reflection data in the acoustic approximation. Geophysics49, 1259–1266.
     [Google Scholar]
  32. TarantolaA.1984b. Linearized inversion of seismic reflection data. Geophysical Prospecting32, 998–1015.
     [Google Scholar]
  33. TarantolaA.1986. A strategy for non‐linear elastic inversion of seismic‐reflection data. Geophysics51, 1893–1903.
     [Google Scholar]
  34. ToppG.C., DavisJ.L. and AnnanA.P.1980. Electromagnetic determination of soil water content: Measurements in coaxial transmission lines. Water Resources Research16, 574–582.
     [Google Scholar]
  35. TronickeJ., DietrichP., WahligU. and AppelE.2002. Integrating surface georadar and crosshole radar tomography: A validation experiment in braided stream deposits. Geophysics67, 1516–1523.
     [Google Scholar]
  36. TronickeJ. and HolligerK.2004. Effects of gas‐ and water‐filled boreholes on the amplitudes of crosshole georadar data as inferred from experimental evidences. Geophysics69, 1255–1260.
     [Google Scholar]
  37. TuressonA.2006. Water content and porosity estimated from ground‐penetrating radar and resistivity. Journal of Applied Geophysics58, 99–111.
     [Google Scholar]
  38. WatanabeT.K., NiheiK.T., NakagawaS. and MyerL.R.2004. Viscoacoustic waveform inversion of transmission data for velocity and attenuation. Journal of the Acoustical Society of America115, 3059–3067.
     [Google Scholar]
  39. WaxmanM.H. and SmitsL.J.M.1968. Electrical conductivities in oil bearing sands. Journal of the society of Petroleum Engineering8,107–122.
     [Google Scholar]
  40. WinshipP., BinleyA. and GomezD.2006. Flow and transport in the unsaturated Sherwood sandstone: Characterization using cross‐borehole geophysical methods. Geological Society, London, Special Publications263, 219–231. doi:10.1144/GSL.SP.2006.263.01.12
     [Google Scholar]
  41. WorthingtonP.F.1993. The uses and abuses of the Archie equations: 1 the formation factor–porosity relationship. Journal of Applied Geophysics30,215–228.
     [Google Scholar]
  42. ZhouB. and GreenhalghS.A.1998a. Crosshole acoustic velocity imaging with the full‐waveform spectral data: 2.5‐D numerical simulations. Exploration Geophysics29, 680–684.
     [Google Scholar]
  43. ZhouB. and GreenhalghS.A.1998b. A damping method for 2.5‐D Green’s function for arbitrary acoustic media. Geophysical Journal International144, 111–120.
     [Google Scholar]
  44. ZhouB. and GreenhalghS.A.2003. Crosshole seismic inversion with normalized full‐waveform amplitude data. Geophysics68, 1320–1330.
     [Google Scholar]
http://instance.metastore.ingenta.com/content/journals/10.3997/1873-0604.2010054
Loading
Full‐waveform inversion of cross‐hole ground‐penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland
Near Surface Geophysics 8, 635 (2010); https://doi.org/10.3997/1873-0604.2010054
/content/journals/10.3997/1873-0604.2010054
/content/journals/10.3997/1873-0604.2010054
Loading

Data & Media loading...

  • Article Type: Research Article

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error