Three‐dimensional Bayesian geostatistical aquifer characterization at the Hanford 300 Area using tracer test data

[1] Tracer tests performed under natural or forced gradient flow conditions can provide useful information for characterizing subsurface properties, through monitoring, modeling, and interpretation of the tracer plume migration in an aquifer. Nonreactive tracer experiments were conducted at the Hanford 300 Area, along with constant-rate injection tests and electromagnetic borehole flowmeter tests. A Bayesian data assimilation technique, the method of anchored distributions (MAD) (Rubin et al., 2010), was applied to assimilate the experimental tracer test data with the other types of data and to infer the three-dimensional heterogeneous structure of the hydraulic conductivity in the saturated zone of the Hanford formation.In this study, the Bayesian prior information on the underlying random hydraulic conductivity field was obtained from previous field characterization efforts using constant-rate injection and borehole flowmeter test data. The posterior distribution of the conductivity field was obtained by further conditioning the field on the temporal moments of tracer breakthrough curves at various observation wells. MAD was implemented with the massively parallel three-dimensional flow and transport code PFLOTRAN to cope with the highly transient flow boundary conditions at the site and to meet the computational demands of MAD. A synthetic study proved that the proposed method could effectively invert tracer test data to capture the essential spatial heterogeneity of the three-dimensional hydraulic conductivity field. Application of MAD to actual field tracer data at the Hanford 300 Area demonstrates that inverting for spatial heterogeneity of hydraulic conductivity under transient flow conditions is challenging and more work is needed.

[1]  Walter A. Illman,et al.  Heterogeneity in hydraulic conductivity and its role on the macroscale transport of a solute plume: From measurements to a practical application of stochastic flow and transport theory , 2010 .

[2]  Paul D. Thorne,et al.  Borehole Completion and Conceptual Hydrogeologic Model for the IFRC Well Field, 300 Area, Hanford Site , 2009 .

[3]  T D Scheibe,et al.  Use of Quantitative Models to Design Microbial Transport Experiments in a Sandy Aquifer , 2001, Ground water.

[4]  A. Bellin,et al.  On the use of peak concentration arrival times for the inference of hydrogeological parameters , 2004 .

[5]  H. H. Franssen,et al.  A comparison of seven methods for the inverse modelling of groundwater flow. Application to the characterisation of well catchments , 2009 .

[6]  Danielle L. Saunders,et al.  Soil Water Balance and Recharge Monitoring at the Hanford Site – FY 2010 Status Report , 2010 .

[7]  J. Zachara,et al.  River‐Induced Flow Dynamics in Long‐Screen Wells and Impact on Aqueous Samples , 2011, Ground water.

[8]  R. A. Behrens,et al.  Incorporating seismic attribute maps in 3‐D reservoir models , 1996 .

[9]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[10]  Keith Beven,et al.  Towards integrated environmental models of everywhere: uncertainty, data and modelling as a learning process , 2007 .

[11]  S. Attinger,et al.  Effective dispersion in heterogeneous media under random transient flow conditions , 2003 .

[12]  Charles F. Harvey,et al.  Temporal Moment‐Generating Equations: Modeling Transport and Mass Transfer in Heterogeneous Aquifers , 1995 .

[13]  Keith Beven,et al.  The future of distributed models: model calibration and uncertainty prediction. , 1992 .

[14]  P. Mantovan,et al.  Hydrological forecasting uncertainty assessment: Incoherence of the GLUE methodology , 2006 .

[15]  Peter K. Kitanidis,et al.  An interactive Bayesian geostatistical inverse protocol for hydraulic tomography , 2006 .

[16]  Peter K. Kitanidis,et al.  Sensitivity of temporal moments calculated by the adjoint-state method and joint inversing of head and tracer data , 2000 .

[17]  Yoram Rubin,et al.  The concept of comparative information yield curves and its application to risk‐based site characterization , 2009 .

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

[19]  G. Evensen.  Sequential data assimilation with a nonlinear quasi‐geostrophic model using Monte Carlo methods to forecast error statistics , 1994 .

[20]  R. A. Behrens,et al.  Incorporating Seismic Attribute Maps in 3D Reservoir Models , 1998 .

[21]  J. Gómez-Hernández,et al.  Uncertainty assessment and data worth in groundwater flow and mass transport modeling using a blocking Markov chain Monte Carlo method. , 2009 .

[22]  P. Reed,et al.  Addressing model bias and uncertainty in three dimensional groundwater transport forecasts for a physical aquifer experiment , 2008 .

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

[24]  Michael N. Fienen,et al.  Obtaining parsimonious hydraulic conductivity fields using head and transport observations: A Bayesian geostatistical parameter estimation approach: COUPLED FLOW AND PARAMETER ESTIMATION , 2009 .

[25]  P. Kitanidis,et al.  An Application of the Geostatistical Approach to the Inverse Problem in Two-Dimensional Groundwater Modeling , 1984 .

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

[27]  Y. Rubin,et al.  Bayesian approach for three-dimensional aquifer characterization at the Hanford 300 Area , 2010 .

[28]  Charles F. Harvey,et al.  When good statistical models of aquifer heterogeneity go bad: A comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields , 2003 .

[29]  Yoram Rubin,et al.  Characterization of aquifer heterogeneity using indicator variables for solute concentrations , 2002 .

[30]  Y. Rubin,et al.  A Bayesian approach for inverse modeling, data assimilation, and conditional simulation of spatial random fields , 2010 .

[31]  Yousu Chen,et al.  Three-Dimensional Groundwater Models of the 300 Area at the Hanford Site, Washington State , 2008 .

[32]  A. Bellin,et al.  The concept of block-effective macrodispersivity and a unified approach for grid-scale- and plume-scale-dependent transport , 1999, Journal of Fluid Mechanics.

[33]  J. Vrugt,et al.  Inverse Modeling of Subsurface Flow and Transport Properties: A Review with New Developments , 2008 .

[34]  Jesús Carrera,et al.  Effective solute transport in temporally fluctuating flow through heterogeneous media , 2005 .

[35]  Jeffrey S. Racine,et al.  Nonparametric Econometrics: The np Package , 2008 .

[36]  P. Kitanidis,et al.  Characterization of mixing and dilution in heterogeneous aquifers by means of local temporal moments , 2000 .

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

[38]  Wolfgang Nowak,et al.  Efficient Computation of Linearized Cross-Covariance and Auto-Covariance Matrices of Interdependent Quantities , 2003 .

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

[40]  Paul D. Thorne,et al.  Uranium Contamination in the Subsurface Beneath the 300 Area, Hanford Site, Washington , 2008 .

[41]  Kuolin Hsu,et al.  Uncertainty assessment of hydrologic model states and parameters: Sequential data assimilation using the particle filter , 2005 .

[42]  Glenn E. Hammond,et al.  Field‐scale model for the natural attenuation of uranium at the Hanford 300 Area using high‐performance computing , 2010 .

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

[44]  Daniel J. Goode,et al.  Apparent dispersion in transient groundwater flow , 1990 .

[45]  C. Zheng,et al.  Multi-Scale Mass Transfer Processes Controlling Natural Attenuation and Engineered Remediation: An IFRC Focused on Hanford’s 300 Area Uranium Plume , 2007 .

[46]  J. Carrera,et al.  Geostatistical inversion of coupled problems: dealing with computational burden and different types of data , 2003 .

[47]  A. Bellin,et al.  On the use of block‐effective macrodispersion for numerical simulations of transport in heterogeneous formations , 2003 .

[48]  A. Michalak,et al.  A geostatistical framework for incorporating transport information in estimating the distribution of a groundwater contaminant plume , 2005 .

[49]  Yoram Rubin,et al.  Modelling of block-scale macrodispersion as a random function , 2011, Journal of Fluid Mechanics.

[50]  G. Dagan.  Stochastic Modeling of Groundwater Flow by Unconditional and Conditional Probabilities: The Inverse Problem , 1985 .

[51]  Nando de Freitas,et al.  An Introduction to MCMC for Machine Learning , 2004, Machine Learning.

[52]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[53]  Andres Alcolea,et al.  Pilot points method incorporating prior information for solving the groundwater flow inverse problem , 2006 .

[54]  Y. Rubin,et al.  Bayesian geostatistical design: Task‐driven optimal site investigation when the geostatistical model is uncertain , 2010 .

[55]  W. Nowak,et al.  Geostatistical inference of hydraulic conductivity and dispersivities from hydraulic heads and tracer data , 2006 .

[56]  Ghislain de Marsily,et al.  Three‐dimensional interference test interpretation in a fractured aquifer using the Pilot Point Inverse Method , 2001 .

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

[58]  Patrick M. Reed,et al.  Many‐objective groundwater monitoring network design using bias‐aware ensemble Kalman filtering, evolutionary optimization, and visual analytics , 2011 .

[59]  Peter K. Kitanidis,et al.  How Observations and Structure Affect the Geostatistical Solution to the Steady‐State Inverse Problem , 1998 .

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

[61]  John Doherty,et al.  Non-uniqueness of inverse transmissivity field calibration and predictive transport modeling , 2003 .

[62]  Keith Beven,et al.  Comment on "Hydrological forecasting uncertainty assessment: incoherence of the GLUE methodology" by Pietro Mantovan and Ezio Todini , 2007 .

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

[64]  Wei Li,et al.  Three‐Dimensional Geostatistical Inversion of Flowmeter and Pumping Test Data , 2008, Ground water.

[65]  A. Ullah,et al.  Nonparametric Econometrics , 1999 .

[66]  Yoram Rubin,et al.  A Geostatistical Approach to the Conditional Estimation of Spatially Distributed Solute Concentration and Notes on the Use of Tracer Data in the Inverse Problem , 1996 .

[67]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[68]  N. T. Burdine.  Relative Permeability Calculations From Pore Size Distribution Data , 1953 .

[69]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data 2. Demonstration on a synthetic aquifer , 1997 .

[70]  G. Dagan,et al.  Stochastic identification of transmissivity and effective recharge in steady groundwater flow: 2. Case study , 1987 .

[71]  Olaf A. Cirpka,et al.  Fully coupled hydrogeophysical inversion of synthetic salt tracer experiments , 2010 .

[72]  Vincent R. Vermeul,et al.  Vertical Wellbore Flow Monitoring for Assessing Spatial and Temporal Flow Relationships with a Dynamic River Boundary , 2010 .

[73]  D. Cox,et al.  An Analysis of Transformations , 1964 .

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

[75]  Stephan R. Sain,et al.  Multi-dimensional Density Estimation , 2004 .