Basin‐Scale Transmissivity and Storativity Estimation Using Hydraulic Tomography

While tomographic inversion has been successfully applied to laboratory- and field-scale tests, here we address the new issue of scale that arises when extending the method to a basin. Specifically, we apply the hydraulic tomography (HT) concept to jointly interpret four multiwell aquifer tests in a synthetic basin to illustrate the superiority of this approach to a more traditional Theis analysis of the same tests. Transmissivity and storativity are estimated for each element of a regional numerical model using the geostatistically based sequential successive linear estimator (SSLE) inverse solution method. We find that HT inversion is an effective strategy for incorporating data from potentially disparate aquifer tests into a basin-wide aquifer property estimate. The robustness of the SSLE algorithm is investigated by considering the effects of noisy observations, changing the variance of the true aquifer parameters, and supplying incorrect initial and boundary conditions to the inverse model. Ground water flow velocities and total confined storage are used as metrics to compare true and estimated parameter fields; they quantify the effectiveness of HT and SSLE compared to a Theis solution methodology. We discuss alternative software that can be used for implementing tomography inversion.

[1]  S. P. Neuman,et al.  Three‐dimensional numerical inversion of pneumatic cross‐hole tests in unsaturated fractured tuff: 2. Equivalent parameters, high‐resolution stochastic imaging and scale effects , 2001 .

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

[3]  John D Bredehoeft From models to performance assessment: the conceptualization problem. , 2003, Ground water.

[4]  P. Hansma,et al.  The scanning ion-conductance microscope. , 1989, Science.

[5]  John D. Bredehoeft,et al.  Ground-water models cannot be validated , 1992 .

[6]  R. W. Andrews,et al.  Sensitivity Analysis for Steady State Groundwater Flow Using Adjoint Operators , 1985 .

[7]  J. Xiang Stochastic Estimation of Hydraulic Spatial Properties , 2007 .

[8]  William W.-G. Yeh,et al.  A stochastic inverse solution for transient groundwater flow: Parameter identification and reliability analysis , 1992 .

[9]  T.-C. Jim Yeh,et al.  Stochastic modelling of groundwater flow and solute transport in aquifers , 1992 .

[10]  Yonghong Hao,et al.  Hydraulic Tomography for Detecting Fracture Zone Connectivity , 2008, Ground water.

[11]  T. Yeh,et al.  Analysis of hydraulic tomography using temporal moments of drawdown recovery data , 2006 .

[12]  Shuyun Liu,et al.  Effectiveness of hydraulic tomography: Sandbox experiments , 2001 .

[13]  V. Pauwels,et al.  Improvement of the PEST parameter estimation algorithm through Extended Kalman Filtering , 2007 .

[14]  T.-C. Jim Yeh,et al.  A Structured Approach for Calibrating Steady‐State Ground‐Water Flow Models , 1996 .

[15]  Keiji Kojima,et al.  Hydropulse Tomography for Identifying 3-D Permeability Distribution , 1993 .

[16]  Minghui Jin,et al.  AN ITERATIVE STOCHASTIC INVERSE METHOD: CONDITIONAL EFFECTIVE TRANSMISSIVITY AND HYDRAULIC HEAD FIELDS , 1995 .

[17]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[18]  T. Harter,et al.  A Numerical Model for Water Flow and Chemical Transport in Variably Saturated Porous Media , 1993 .

[19]  T. Yeh,et al.  Hydraulic tomography: Development of a new aquifer test method , 2000 .

[20]  Cheng-Haw Lee,et al.  Time to Change the Way We Collect and Analyze Data for Aquifer Characterization , 2007, Ground water.

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

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

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

[24]  T.-C. Jim Yeh,et al.  An inverse model for three‐dimensional flow in variably saturated porous media , 2000 .

[25]  Walter A. Illman,et al.  Steady-state hydraulic tomography in a laboratory aquifer with deterministic heterogeneity: Multi-method and multiscale validation of hydraulic conductivity tomograms , 2007 .

[26]  Mary C. Hill,et al.  UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API , 2006 .

[27]  Yan Chen,et al.  Data assimilation for transient flow in geologic formations via ensemble Kalman filter , 2006 .

[28]  T.-C. Jim Yeh,et al.  An iterative geostatistical inverse method for steady flow in the vadose zone , 1996 .

[29]  Peter Dietrich,et al.  A travel time based hydraulic tomographic approach , 2003 .

[30]  N Oreskes,et al.  Verification, Validation, and Confirmation of Numerical Models in the Earth Sciences , 1994, Science.

[31]  Peter Dietrich,et al.  Identification of the permeability distribution in soil by hydraulic tomography , 1995 .

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

[33]  Heidi Christiansen Barlebo,et al.  Investigating the Macrodispersion Experiment (MADE) site in Columbus, Mississippi, using a three‐dimensional inverse flow and transport model , 2004 .

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

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

[36]  Michael Edward Hohn,et al.  An Introduction to Applied Geostatistics: by Edward H. Isaaks and R. Mohan Srivastava, 1989, Oxford University Press, New York, 561 p., ISBN 0-19-505012-6, ISBN 0-19-505013-4 (paperback), $55.00 cloth, $35.00 paper (US) , 1991 .

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

[38]  M. Marietta,et al.  Pilot Point Methodology for Automated Calibration of an Ensemble of conditionally Simulated Transmissivity Fields: 1. Theory and Computational Experiments , 1995 .

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