Jointly Mapping Hydraulic Conductivity and Porosity by Assimilating Concentration Data via Ensemble Kalman Filter

Summary Real-time data from on-line sensors offer the possibility to update environmental simulation models in real-time. Information from on-line sensors concerning contaminant concentrations in groundwater allow for the real-time characterization and control of a contaminant plume. In this paper it is proposed to use the CPU-efficient Ensemble Kalman Filter (EnKF) method, a data assimilation algorithm, for jointly updating the flow and transport parameters (hydraulic conductivity and porosity) and state variables (piezometric head and concentration) of a groundwater flow and contaminant transport problem. A synthetic experiment is used to demonstrate the capability of the EnKF to estimate hydraulic conductivity and porosity by assimilating dynamic head and multiple concentration data in a transient flow and transport model. In this work the worth of hydraulic conductivity, porosity, piezometric head, and concentration data is analyzed in the context of aquifer characterization and prediction uncertainty reduction. The results indicate that the characterization of the hydraulic conductivity and porosity fields is continuously improved as more data are assimilated. Also, groundwater flow and mass transport predictions are improved as more and different types of data are assimilated. The beneficial impact of accounting for multiple concentration data is patent.

[1]  Andrés Sahuquillo,et al.  Joint simulation of transmissivity and storativity fields conditional to steady-state and transient hydraulic head data , 1999 .

[2]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

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

[4]  J. Jaime Gómez-Hernández,et al.  Groundwater flow inverse modeling in non-MultiGaussian media: Performance assessment of the normal-score Ensemble Kalman Filter , 2011 .

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

[6]  Harrie-Jan Hendricks Franssen,et al.  Ensemble Kalman filtering versus sequential self-calibration for inverse modelling of dynamic groundwater flow systems , 2009 .

[7]  J. Jaimegomezhernandez,et al.  ISIM3D: An ANSI-C three-dimensional multiple indicator conditional simulation program , 1990 .

[9]  Dongxiao Zhang,et al.  Investigation of flow and transport processes at the MADE site using ensemble Kalman filter , 2008 .

[10]  Arlen W. Harbaugh,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process , 2000 .

[11]  Bill X. Hu,et al.  Using data assimilation method to calibrate a heterogeneous conductivity field and improve solute transport prediction with an unknown contamination source , 2009 .

[12]  Liangping Li,et al.  An approach to handling non-Gaussianity of parameters and state variables in ensemble Kalman filtering , 2011 .

[13]  G. Evensen,et al.  Analysis Scheme in the Ensemble Kalman Filter , 1998 .

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

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

[16]  Sebastien Strebelle,et al.  Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics , 2002 .

[17]  J. Jaime Gómez-Hernández,et al.  ISIM3D: and ANSI-C three-dimensional multiple indicator conditional simulation program , 1990 .

[18]  Qi Han,et al.  On Integrating Groundwater Transport Models with Wireless Sensor Networks , 2010, Ground water.

[19]  Wolfgang Nowak,et al.  Parameter Estimation by Ensemble Kalman Filters with Transformed Data , 2010 .

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

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

[22]  Dean S. Oliver,et al.  The Ensemble Kalman Filter for Continuous Updating of Reservoir Simulation Models , 2006 .

[23]  Harrie-Jan Hendricks Franssen,et al.  Modeling transient groundwater flow by coupling ensemble Kalman filtering and upscaling , 2012 .

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

[25]  J. J. Gómez-Hernández,et al.  Probabilistic assessment of travel times in groundwater modeling , 1994 .

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

[27]  Jeffrey P. Walker,et al.  Extended versus Ensemble Kalman Filtering for Land Data Assimilation , 2002 .

[28]  A. Journel.  Geostatistics for Conditional Simulation of Ore Bodies , 1974 .

[29]  Carlos Llopis-Albert,et al.  Gradual conditioning of non-Gaussian transmissivity fields to flow and mass transport data: 1. Theory , 2009 .

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

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

[32]  R. L. Naff,et al.  High‐resolution Monte Carlo simulation of flow and conservative transport in heterogeneous porous media: 2. Transport results , 1998 .

[33]  Liangping Li,et al.  Transport upscaling using multi-rate mass transfer in three-dimensional highly heterogeneous porous media , 2011 .

[34]  J. Jaime Gómez-Hernández,et al.  Three-dimensional hydraulic conductivity upscaling in groundwater modeling , 2010, Comput. Geosci..

[35]  J. Gómez-Hernández,et al.  Simulation of Non-Gaussian Transmissivity Fields Honoring Piezometric Data and Integrating Soft and Secondary Information , 1999 .

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

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

[38]  Alexander Y. Sun,et al.  Sequential updating of multimodal hydrogeologic parameter fields using localization and clustering techniques , 2009 .

[39]  S. Ge,et al.  Simultaneous rejuvenation and aging of groundwater in basins due to depth‐decaying hydraulic conductivity and porosity , 2010 .

[40]  J. Gómez-Hernández,et al.  Joint Sequential Simulation of MultiGaussian Fields , 1993 .

[41]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

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

[43]  Xian-Huan Wen,et al.  Stochastic inverse mapping of hydraulic conductivity and sorption partitioning coefficient fields conditioning on nonreactive and reactive tracer test data , 2004 .

[44]  Liangping Li,et al.  A comparative study of three-dimensional hydraulic conductivity upscaling at the macro-dispersion experiment (MADE) site, Columbus Air Force Base, Mississippi (USA) , 2011 .

[45]  Wolfgang Nowak,et al.  Best unbiased ensemble linearization and the quasi‐linear Kalman ensemble generator , 2009 .

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

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

[48]  Alberto Guadagnini,et al.  Relative importance of geostatistical and transport models in describing heavily tailed breakthrough curves at the Lauswiesen site. , 2008, Journal of contaminant hydrology.

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

[50]  Jesús Carrera,et al.  Coupled estimation of flow and solute transport parameters , 1996 .

[51]  Xian-Huan Wen,et al.  Real-Time Reservoir Model Updating Using Ensemble Kalman Filter , 2005 .

[52]  Bill X. Hu,et al.  Examining the influence of heterogeneous porosity fields on conservative solute transport. , 2009, Journal of contaminant hydrology.

[53]  Harrie-Jan Hendricks Franssen,et al.  Operational real‐time modeling with ensemble Kalman filter of variably saturated subsurface flow including stream‐aquifer interaction and parameter updating , 2011 .

[54]  A. Hassan.  Water flow and solute mass flux in heterogeneous porous formations with spatially random porosity , 2001 .

[55]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[56]  R. L. Naff,et al.  High‐resolution Monte Carlo simulation of flow and conservative transport in heterogeneous porous media: 1. Methodology and flow results , 1998 .

[57]  W. Yeh.  Review of Parameter Identification Procedures in Groundwater Hydrology: The Inverse Problem , 1986 .

[58]  Geir Nævdal,et al.  Reservoir Monitoring and Continuous Model Updating Using Ensemble Kalman Filter , 2005 .

[59]  W. Kinzelbach,et al.  Real‐time groundwater flow modeling with the Ensemble Kalman Filter: Joint estimation of states and parameters and the filter inbreeding problem , 2008 .

[60]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

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

[62]  C. Llopis-Albert,et al.  Gradual conditioning of non-Gaussian transmissivity fields to flow and mass transport data: 2. Demonstration on a synthetic aquifer , 2009 .

[63]  Olaf A. Cirpka,et al.  Use of steady-state concentration measurements in geostatistical inversion , 2009 .

[64]  Jeffrey L. Anderson.  An Ensemble Adjustment Kalman Filter for Data Assimilation , 2001 .

[65]  R. L. Naff,et al.  Nonreactive and reactive solute transport in three-dimensional heterogeneous porous media: Mean displacement, plume spreading, and uncertainty , 1994 .

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