A Bayesian approach to stochastic capture zone delineation incorporating tracer arrival times, conductivity measurements, and hydraulic head observations

This paper presents a methodology to invert tracer arrival times and to incorporate travel time data in the delineation of well capture zones. Within a Bayesian framework the observed arrival times are used to obtain probability-based weights for each realization of the hydraulic conductivity field. Realizations that closely reproduce the observed arrival times are more likely to represent the "true" or real conductivity field than realizations yielding worse predictions, and are consequently characterized by a higher conditional probability. In the prediction of the capture zones the realizations are weighted by their respective probability. We combine the arrival time data with conductivity measurements and head observations to delineate stochastic capture zones. The conductivity measurements update the prior distributions specified for the unknown structural parameters of the conductivity field, and are used as conditioning data in the generation of conditional conductivity fields. The parameter values used to generate the conductivity realizations are sampled by Monte Carlo from the updated parameter distributions. The head and travel time observations are subsequently used to assign probability-based weights to the conductivity realizations by applying Bayes' theorem. The hyperparameters of the error model are assumed unknown, and their effect is accounted for by marginalization. A synthetic flow setup consisting of a single extraction well in a fully confined aquifer under a regional gradient is used to illustrate the method. We evaluate the relative worth of the different data types by introducing them sequentially in the inverse methodology. Results indicate that the different data types are complementary and that the travel time data are essential to improve the predictions of the capture zones. This is confirmed by the results for the case where uncertainty in the homogeneous porosity is accounted for.

[1]  Steven M. Gorelick,et al.  Inferring the relation between seismic slowness and hydraulic conductivity in heterogeneous aquifers , 2000 .

[2]  Y. Rubin,et al.  Bayesian Method for hydrogeological site characterization using borehole and geophysical survey data: Theory and application to the Lawrence Livermore National Laboratory Superfund Site , 1999 .

[3]  Geostatistical Formulations of Groundwater Coupled Inverse Problems , 1993 .

[4]  Keith Beven,et al.  Bayesian methodology for stochastic capture zone delineation incorporating transmissivity measurements and hydraulic head observations , 2003 .

[5]  Yoram Rubin,et al.  Mapping permeability in heterogeneous aquifers using hydrologic and seismic data , 1992 .

[6]  A. S. Cullick,et al.  Construction of geostatistical aquifer models integrating dynamic flow and tracer data using inverse technique , 2002 .

[7]  Peter J. Diggle,et al.  An Introduction to Model-Based Geostatistics , 2003 .

[8]  Y. Rubin,et al.  The moments of the breakthrough curves of instantaneously and kinetically sorbing solutes in heterogeneous geologic media: Prediction and parameter inference from field measurements , 1997 .

[9]  Y. Rubin,et al.  A Full‐Bayesian Approach to parameter inference from tracer travel time moments and investigation of scale effects at the Cape Cod Experimental Site , 2000 .

[10]  Y. Rubin,et al.  Spatial correlation structure estimation using geophysical and hydrogeological data , 1999 .

[11]  J. C. Ramírez,et al.  Estimation of aquifer parameters under transient and steady-state conditions , 1984 .

[12]  J. H. Cushman,et al.  Inverse methods for subsurface flow: A critical review of stochastic techniques , 1990 .

[13]  D. McLaughlin,et al.  A Reassessment of the Groundwater Inverse Problem , 1996 .

[14]  G. Dagan,et al.  Conditional estimation of solute travel time in heterogeneous formations: Impact of transmissivity measurements , 1992 .

[15]  Alberto Guadagnini,et al.  Time‐Related Capture Zones for Contaminants in Randomly Heterogeneous Formations , 1999 .

[16]  J. Bear,et al.  On the movement of water bodies injected into aquifers , 1965 .

[17]  L. Hu.  Gradual Deformation and Iterative Calibration of Gaussian-Related Stochastic Models , 2000 .

[18]  Matthijs van Leeuwen,et al.  Stochastic determination of well capture zones , 1998 .

[19]  Angelos N. Findikakis,et al.  Quantitative Estimates of the Uncertainty in the Evaluation of Ground Water Remediation Schemes , 2000 .

[20]  D. Oliver,et al.  Markov chain Monte Carlo methods for conditioning a permeability field to pressure data , 1997 .

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

[22]  J. H. Cushman,et al.  A FAST FOURIER TRANSFORM STOCHASTIC ANALYSIS OF THE CONTAMINANT TRANSPORT PROBLEM , 1993 .

[23]  P. Diggle,et al.  Bayesian Inference in Gaussian Model-based Geostatistics , 2002 .

[24]  M. G. Marietta,et al.  Pilot Point Methodology for Automated Calibration of an Ensemble of Conditionally Simulated Transmissivity Fields: 2. Application , 1995 .

[25]  A. Dassargues,et al.  Delineation of groundwater protection zones based on tracer tests and transport modeling in alluvial sediments , 1998 .

[26]  Allan D. Woodbury,et al.  Inversion of the Borden Tracer Experiment data: Investigation of stochastic moment models , 1992 .

[27]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[28]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

[29]  J. Harris,et al.  Coupled seismic and tracer test inversion for aquifer property characterization , 1993 .

[30]  Charles F. Harvey,et al.  Mapping Hydraulic Conductivity: Sequential Conditioning with Measurements of Solute Arrival Time, Hydraulic Head, and Local Conductivity , 1995 .

[31]  B. Ursin,et al.  Prediction of Reservoir Variables Based on Seismic Data and Well Observations , 2002 .

[32]  Y. Rubin,et al.  Analysis of the Cape Cod tracer data , 1997 .

[33]  J. W. Barker,et al.  Quantifying Uncertainty in Production Forecasts: Another Look at the PUNQ-S3 Problem , 2001 .

[34]  Matthijs van Leeuwen,et al.  Stochastic determination of well capture zones conditioned on regular grids of transmissivity measurements , 2000 .

[35]  Mike West,et al.  Markov Random Field Models for High-Dimensional Parameters in Simulations of Fluid Flow in Porous Media , 2002, Technometrics.

[36]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

[37]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

[38]  Keith Beven,et al.  Stochastic capture zone delineation within the generalized likelihood uncertainty estimation methodology: Conditioning on head observations , 2001 .

[39]  J. M. Shafer,et al.  Assessment of Uncertainty in Time‐Related Capture Zones Using Conditional Simulation of Hydraulic Conductivity , 1991 .

[40]  Peter J. Diggle,et al.  Bayesian methodology to stochastic capture zone determination: Conditioning on transmissivity measurements , 2002 .

[41]  G. Dagan,et al.  A solute flux approach to transport in heterogeneous formations: 1. The general framework , 1992 .