A field proof‐of‐concept of aquifer imaging using 3‐D transient hydraulic tomography with modular, temporarily‐emplaced equipment

[1] Hydraulic tomography is a field scale aquifer characterization method capable of estimating 3-D heterogeneous parameter distributions, and is directly sensitive to hydraulic conductivity (K), thus providing a useful data source for improving flow and transport models. We present results from a proof-of-concept field and modeling study in which we apply 3-D transient hydraulic tomography (3DTHT) to the relatively high-K and moderately heterogeneous unconfined aquifer at the Boise Hydrogeophysical Research Site. Short-duration (20 min) partially penetrating pumping tests, for which observed responses do not reach steady state, are used as the aquifer stimulation. To collect field data, we utilize a system of temporarily emplaced packer equipment to isolate multiple discrete intervals in boreholes. To analyze the data, we utilize MODFLOW combined with geostatistical inversion code based on the quasilinear approach of Kitanidis (1995). This combination of practical software allows inversion of large datasets (>250 drawdown curves, and almost 1000 individual data points) and estimation of K at >100,000 locations; reasonable runtimes are obtained using a single multicore computer with 12 GB of RAM. The K heterogeneity results from 3DTHT are cross-validated against K characterization from a large set of partially penetrating slug tests, and found to be quite consistent. The use of portable, modular equipment for field implementation means that 3DTHT data collection can be performed (including mobilization/demobilization) within a matter of days. Likewise, use of a practical, efficient and scalable numerical modeling and inversion strategy means that computational effort is drastically reduced, such that 3-D aquifer property distributions can be estimated quickly.

[1]  Michael A. Celia,et al.  Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 3. Hydraulic c , 1992 .

[2]  Andrew Binley,et al.  Structural joint inversion of time‐lapse crosshole ERT and GPR traveltime data , 2010 .

[3]  Peter Dietrich,et al.  A Rapid Method for Hydraulic Profiling in Unconsolidated Formations , 2008, Ground water.

[4]  Robert P Chapuis,et al.  "Characterizing hydraulic conductivity with the direct-push permeameter," by James J. Butler Jr., Peter Dietrich, Volker Wittig, and Tom Christy, July–August 2007, v. 45, no. 4:409-419. , 2010, Ground water.

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

[6]  W. W. Wood,et al.  Large-Scale Natural Gradient Tracer Test in Sand and Gravel, , 1991 .

[7]  Thomas Ptak,et al.  Dual-tracer transport experiments in a physically and chemically heterogeneous porous aquifer: effective transport parameters and spatial variability , 1996 .

[8]  Walter A. Illman,et al.  Three‐dimensional transient hydraulic tomography in a highly heterogeneous glaciofluvial aquifer‐aquitard system , 2011 .

[9]  R. Knight,et al.  Improving crosshole radar velocity tomograms: A new approach to incorporating high-angle traveltime data , 2007 .

[10]  Michael D. Knoll,et al.  Capacitive conductivity logging and electrical stratigraphy in a high-resistivity aquifer, Boise Hydrogeophysical Research Site , 2009 .

[11]  W. Li,et al.  Two‐dimensional characterization of hydraulic heterogeneity by multiple pumping tests , 2007 .

[12]  L. Toran,et al.  Stream Bottom Resistivity Tomography to Map Ground Water Discharge , 2008, Ground water.

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

[14]  John M Healey,et al.  Hydraulic Tests with Direct‐Push Equipment , 2002, Ground water.

[15]  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 .

[16]  P. Gilbert Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

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

[18]  Junfeng Zhu,et al.  Traditional analysis of aquifer tests: Comparing apples to oranges? , 2005 .

[19]  James Irving,et al.  Inversion of Multiple Intersecting High-Resolution Crosshole GPR Profiles for Hydrological Characterization at the Boise Hydrogeophysical Research Site , 2011 .

[20]  Denis R. LeBlanc,et al.  A Statistical Evaluation of Formation Disturbance Produced by Well‐Casing Installation Methods , 1988 .

[21]  Phil Long,et al.  Geophysical monitoring of hydrological and biogeochemical transformations associated with Cr(VI) bioremediation. , 2008, Environmental science & technology.

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

[23]  Michael D. Knoll,et al.  VSP traveltime inversion: Near‐surface issues , 2004 .

[24]  D. Genereux,et al.  A borehole flowmeter investigation of small‐scale hydraulic conductivity variation in the Biscayne Aquifer, Florida , 2001 .

[25]  E. Sudicky A natural gradient experiment on solute transport in a sand aquifer: Spatial variability of hydraulic conductivity and its role in the dispersion process , 1986 .

[26]  Tom Clemo,et al.  Field, laboratory, and modeling investigation of the skin effect at wells with slotted casing, Boise Hydrogeophysical Research Site , 2006 .

[27]  Tian-Chyi J. Yeh,et al.  Characterization of aquifer heterogeneity using transient hydraulic tomography , 2004 .

[28]  Thomas Aigner,et al.  Towards realistic aquifer models: Three-dimensional georadar surveys of Quaternary gravel deltas (Si , 1999 .

[29]  David L. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 2. Spatial moments and the advection and dispersion of nonreactive tracers , 1986 .

[30]  K. Karasaki,et al.  Estimation of reservoir properties using transient pressure data: An asymptotic approach , 2000 .

[31]  Linda M. Abriola,et al.  The influence of field-scale heterogeneity on the infiltration and entrapment of dense nonaqueous phase liquids in saturated formations , 2000 .

[32]  James J Butler,et al.  A Simple Correction for Slug Tests in Small‐Diameter Wells , 2002, Ground water.

[33]  Martin Sauter,et al.  Cross-well slug interference tests: An effective characterization method for resolving aquifer heterogeneity , 2010 .

[34]  D. W. Vasco,et al.  Zeroth-order inversion of transient pressure observations , 2008 .

[35]  Michael D. Annable,et al.  Estimating nonaqueous phase liquid spatial variability using partitioning tracer higher temporal moments , 2003 .

[36]  Michael D. Knoll,et al.  Multivariate analysis of cross‐hole georadar velocity and attenuation tomograms for aquifer zonation , 2004 .

[37]  D. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 1. Approach and overview of plume movement , 1986 .

[38]  Z. J. Kabala,et al.  The dipole flow test: A new single‐borehole test for aquifer characterization , 1993 .

[39]  James J. Butler,et al.  Steady shape analysis of tomographic pumping tests for characterization of aquifer heterogeneities , 2002 .

[40]  John A. Cherry,et al.  A New Multilevel Ground Water Monitoring System Using Multichannel Tubing , 2002 .

[41]  G. Fogg,et al.  Dispersion of groundwater age in an alluvial aquifer system , 2002 .

[42]  Frederick L. Paillet,et al.  Flow modeling and permeability estimation using borehole flow logs in heterogeneous fractured formations , 1998 .

[43]  H. Bouwer,et al.  A slug test for determining hydraulic conductivity of unconfined aquifers with completely or partially penetrating wells , 1976 .

[44]  Richard L. Smith,et al.  Transport of microspheres and indigenous bacteria through a sandy aquifer: Results of natural- and forced-gradient tracer experiments , 1989 .

[45]  Jian Yang,et al.  Natural Gradient , 2009, Encyclopedia of Biometrics.

[46]  H. Vereecken,et al.  Imaging and characterisation of subsurface solute transport using electrical resistivity tomography (ERT) and equivalent transport models , 2002 .

[47]  Bwalya Malama,et al.  Information content of slug tests for estimating hydraulic properties in realistic, high-conductivity aquifer scenarios , 2011 .

[48]  Cheng-Haw Lee,et al.  River stage tomography: A new approach for characterizing groundwater basins , 2009 .

[49]  Brian R. Zurbuchen,et al.  Dipole Probe: Design and Field Applications of a Single‐Borehole Device for Measurements of Vertical Variations of Hydraulic Conductivity , 1998 .

[50]  S. Hubbard,et al.  Joint inversion of crosshole radar and seismic traveltimes , 2008 .

[51]  Alan G. Green,et al.  Mapping the architecture of glaciofluvial sediments with three-dimensional georadar , 1995 .

[52]  Vitaly A. Zlotnik,et al.  Multi-level slug tests in highly permeable formations : 1. Modification of the Springer-Gelhar (SG) model , 1998 .

[53]  D. W. Vasco,et al.  Interpretation and inversion of low‐frequency head observations , 2006 .

[54]  Allen M. Shapiro,et al.  Crosswell seismic investigation of hydraulically conductive, fractured bedrock near Mirror Lake, New Hampshire , 2002 .

[55]  Geoffrey C. Bohling,et al.  A field assessment of the value of steady shape hydraulic tomography for characterization of aquifer heterogeneities , 2007 .

[56]  A. W. Harbaugh MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process , 2005 .

[57]  Philippe Renard,et al.  Three-dimensional high resolution fluvio-glacial aquifer analog: Part 1: Field study , 2011 .

[58]  Michael D. Knoll,et al.  Investigating the stratigraphy of an alluvial aquifer using crosswell seismic traveltime tomography , 2006 .

[59]  Don W. Vasco,et al.  Inversion of pressure observations: an integral formulation , 2001 .

[60]  Warren Barrash,et al.  Hydrological parameter estimations from a conservative tracer test with variable‐density effects at the Boise Hydrogeophysical Research Site , 2011 .

[61]  Hiromitsu Saegusa,et al.  Hydraulic tomography in fractured granite: Mizunami Underground Research site, Japan , 2009 .

[62]  J. J. Butler,et al.  The Design, Performance, and Analysis of Slug Tests , 1997 .

[63]  Paul V. Roberts,et al.  A natural gradient experiment on solute transport in a sand aquifer: 4. Sorption of organic solutes and its influence on mobility , 1986 .

[64]  Kenichi Soga,et al.  The significance of heterogeneity on mass flux from DNAPL source zones: an experimental investigation. , 2007, Journal of contaminant hydrology.

[65]  Peter Dietrich,et al.  Characterizing Hydraulic Conductivity with the Direct‐Push Permeameter , 2007, Ground water.

[66]  P. Kitanidis Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[67]  Brady Johnson,et al.  Analytical Modeling of Saturated Zone Head Response to Evapotranspiration and River-Stage Fluctuations , 2010 .

[68]  Tian-Chyi J. Yeh,et al.  Robustness of joint interpretation of sequential pumping tests: Numerical and field experiments , 2011 .

[69]  P. Kitanidis,et al.  An Application of Bayesian Inverse Methods to Vertical Deconvolution of Hydraulic Conductivity in a Heterogeneous Aquifer at Oak Ridge National Laboratory , 2004 .

[70]  Charles F. Harvey,et al.  What does a slug test measure: An investigation of instrument response and the effects of heterogeneity , 2002 .

[71]  L. Slater Near Surface Electrical Characterization of Hydraulic Conductivity: From Petrophysical Properties to Aquifer Geometries—A Review , 2007 .

[72]  Michael D. Knoll,et al.  Reflectivity Modeling of a Ground-Penetrating-Radar Profile of a Saturated Fluvial Formation , 2006 .

[73]  Craig T. Simmons,et al.  Groundwater ages in fractured rock aquifers , 2005 .

[74]  M. Riva,et al.  Joint Inversion of Steady-State Hydrologic and Self-Potential Data for 3D Hydraulic Conductivity Distribution at the Boise Hydrogeophysical Research Site , 2011 .

[75]  Tom Clemo,et al.  MODFLOW-2005 Ground-Water Model - User Guide to the Adjoint State Based Sensitivity Process (ADJ) , 2007 .

[76]  Junfeng Zhu,et al.  Laboratory sandbox validation of transient hydraulic tomography , 2007 .

[77]  M. Celia,et al.  Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts, 2, Analysis of spatial moments for a nonreactive tracer , 1991 .

[78]  Andrew Binley,et al.  Electrical resistivity imaging of the architecture of substream sediments , 2008 .

[79]  Peter Dietrich,et al.  Direct Push-Technologies , 2006 .

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

[81]  S. Gorelick,et al.  Estimating lithologic and transport properties in three dimensions using seismic and tracer data , 1996 .

[82]  Peter Dietrich,et al.  A field assessment of high‐resolution aquifer characterization based on hydraulic travel time and hydraulic attenuation tomography , 2011 .

[83]  Andrew Binley,et al.  Geophysical characterization of riverbed hydrostratigraphy using electrical resistance tomography , 2010 .

[84]  Alfred E. Hess,et al.  Identifying hydraulically conductive fractures with a slow-velocity borehole flowmeter , 1986 .

[85]  John E. Peterson,et al.  Crosswell seismic imaging in a contaminated basalt aquifer , 2001 .

[86]  Linda M. Abriola,et al.  The influence of field-scale heterogeneity on the surfactant-enhanced remediation of entrapped nonaqueous phase liquids , 2000 .

[87]  David E. Schaad,et al.  The dipole-flow test with a tracer: a new single-borehole tracer test for aquifer characterization , 2000 .

[88]  W. P. Clement,et al.  Crosshole Radar Tomography in a Fluvial Aquifer Near Boise, Idaho , 2006 .

[89]  A Revil,et al.  A Potential‐Based Inversion of Unconfined Steady‐State Hydraulic Tomography , 2009, Ground water.

[90]  Walter A Illman,et al.  Field Study of Hydrogeologic Characterization Methods in a Heterogeneous Aquifer , 2011, Ground water.

[91]  R. Glass,et al.  Bench-scale visualization of DNAPL remediation processes in analog heterogeneous aquifers: surfactant floods and in situ oxidation using permanganate. , 2002, Journal of contaminant hydrology.

[92]  S. Gorelick,et al.  Saline tracer visualized with three‐dimensional electrical resistivity tomography: Field‐scale spatial moment analysis , 2005 .

[93]  Andrew Binley,et al.  An Investigation of the Ability of Induced Polarization to Resolve Aquifer Heterogeneity in an Unconsolidated Sedimentary Aquifer , 2011 .

[94]  Michael D. Knoll,et al.  Traveltime inversion of vertical radar profiles , 2006 .

[95]  Frederick L. Paillet,et al.  Using flowmeter pulse tests to define hydraulic connections in the subsurface: a fractured shale example , 2002 .

[96]  P. Jardine,et al.  Plume and lithologic profiling with surface resistivity and seismic tomography , 2005, Ground water.

[97]  John A. Cherry,et al.  Tritium and helium: 3 as groundwater age tracers in the Borden Aquifer , 1992 .

[98]  S. C. Young,et al.  Borehole flowmeters: field application and data analysis , 1994 .

[99]  M. Cardiff,et al.  3‐D transient hydraulic tomography in unconfined aquifers with fast drainage response , 2011 .

[100]  Mark N. Goltz,et al.  A natural gradient experiment on solute transport in a sand aquifer: 3. Retardation estimates and mass balances for organic solutes , 1986 .

[101]  Geoffrey C. Bohling,et al.  Sensitivity and resolution of tomographic pumping tests in an alluvial aquifer , 2007 .

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

[103]  A. Binley,et al.  Cross-hole electrical imaging of a controlled saline tracer injection , 2000 .

[104]  John T. Wilson,et al.  Measuring vertical profiles of hydraulic conductivity with in situ direct-push methods , 2000 .

[105]  W. Barrash,et al.  Significance of porosity for stratigraphy and textural composition in subsurface, coarse fluvial deposits: Boise Hydrogeophysical Research Site , 2004 .

[106]  Junfeng Zhu,et al.  Sequential aquifer tests at a well field, Montalto Uffugo Scalo, Italy , 2007 .

[107]  John H. Bradford,et al.  Estimating porosity with ground‐penetrating radar reflection tomography: A controlled 3‐D experiment at the Boise Hydrogeophysical Research Site , 2009 .