A field comparison of Fresnel zone and ray-based GPR attenuation-difference tomography for time-lapse imaging of electrically anomalous tracer or contaminant plumes

Ground-penetrating radar GPR attenuation-difference tomographyisausefultoolforimagingthemigrationofelectrically anomalous tracer or contaminant plumes. Attenuation-difference tomography uses the difference in the trace amplitudes of tomographic data sets collected at different times to image the distribution of bulk-conductivity changes within the medium. Themostcommonapproachforcomputingthetomographicsensitivities uses ray theory, which is well understood and leads to efficient computations. However, ray theory requires the assumption that waves propagate at infinite frequency, and thus sensitivities are distributed along a line between the source and receiver. The infinite-frequency assumption in ray theory leads to a significant loss of resolution both spatially and in terms of amplitude of the recovered image. We use scattering theory to approximate the sensitivity of electromagnetic EM wave amplitude to changes in bulk conductivity within the medium. These sensitivities occupy the first Fresnel zone, account for the finite frequency nature of propagating EM waves, and are valid whenvelocityvariationswithinthemediumdonotcausesignificant ray bending. We evaluate the scattering theory sensitivities by imaging a bromide tracer plume as it migrates through a coarse alluvial aquifer over two successive days. The scattering theory tomograms display a significant improvement in resolution over the ray-based counterparts, as shown by a direct comparison of the tomograms and also by a comparison of the verticalfluidconductivitydistributionmeasuredinamonitoringwell, located within the tomographic plane. By improving resolution, the scattering theory sensitivities increase the utility of GPR attenuation-difference tomography for monitoring the movement of electrically anomalous plumes. In addition, the improved accuracy of information gathered through attenuation-difference tomography using scattering theory is a positive step toward futuredevelopmentsinusingGPRdatatohelpcharacterizethedistributionofhydrogeologicproperties.

[1]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[2]  Frederick D. Day-Lewis,et al.  Combined interpretation of radar, hydraulic, and tracer data from a fractured-rock aquifer near Mirror Lake, New Hampshire, USA , 2006 .

[3]  R. Lytle,et al.  Investigation of fracture flow paths using alterant geophysical tomography , 1986 .

[4]  Frederick D. Day-Lewis,et al.  Time‐lapse inversion of crosswell radar data , 2002 .

[5]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[6]  O. Olsson,et al.  BOREHOLE RADAR APPLIED TO THE CHARACTERIZATION OF HYDRAULICALLY CONDUCTIVE FRACTURE ZONES IN CRYSTALLINE ROCK1 , 1992 .

[7]  Frederick D. Day-Lewis,et al.  Time‐lapse imaging of saline‐tracer transport in fractured rock using difference‐attenuation radar tomography , 2003 .

[8]  W. Barrash,et al.  Water Sampling and Analysis for the Tracer/Time-Lapse Radar Imaging Test at the Boise Hydrogeophysical Research Site , 2002 .

[9]  D. Hyndman,et al.  Tracer/Time-Lapse Radar Imaging Test; Design, Operation, and Preliminary Results , 2002 .

[10]  Per Christian ANALYSIS OF DISCRETE ILL-POSED PROBLEMS BY MEANS OF THE L-CURVE* , 1992 .

[11]  P. Routh,et al.  Fresnel volume georadar attenuation‐difference tomography , 2005 .

[12]  Andrew Binley,et al.  Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution‐dependent limitations , 2005 .

[13]  K. Holliger,et al.  Numerical modeling of borehole georadar data , 2002 .

[14]  E. L. Majer,et al.  Wavepath traveltime tomography , 1993 .

[15]  Tom Clemo,et al.  Boise Hydrogeophysical Research Site (Bhrs): Objectives, Design, Initial Geostatistical Results , 1999 .

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

[17]  Vlastislav Cerveny,et al.  Fresnel volume ray tracing , 1992 .