Contaminant release history identification in 2-D heterogeneous aquifers through a minimum relative entropy approach

The minimum relative entropy (MRE) method has been applied in a wide variety of fields since it was first introduced. Woodbury and Ulrych (Water Resour Res 29(8): 2847–2860, 1993, Water Resour Res 32(9): 2671–2681, 1996) adopted and improved this method to solve linear inverse problems in aquifers. In this work, the MRE method was improved to detect the source release history in 2-D aquifer characterized by a non-uniform flow-field. The approach was tested on two cases: a 2-D homogeneous conductivity field and a heterogeneous one (the hydraulic conductivity presents three orders of magnitude in terms of variability). In the latter case the transfer function cannot be described with an analytical formulation, thus, the transfer functions were estimated by means of a numerical procedure. In order to analyze the method performance in different conditions, two datasets have been used: observations collected at the same time at 20 different monitoring points, and observations collected at 2 monitoring points at several times. The observed data have been processed with and without a random error and the Boxcar and Gaussian probability distribution functions were considered as a priori information. The agreement between the true and the estimated data has been evaluated through the calculation of the normalized Root Mean Square error. The approach was able to recover the release history even in the most difficult case.

[1]  Kurt Roth,et al.  Transfer Functions and Solute Movement Through Soil: Theory and Applications , 1991 .

[2]  T. Skaggs,et al.  Recovering the release history of a groundwater contaminant , 1994 .

[3]  A. Zanini,et al.  Contaminant source and release history identification in groundwater: a multi-step approach. , 2014, Journal of contaminant hydrology.

[4]  T. Ulrych,et al.  Minimum relative entropy: Forward probabilistic modeling , 1993 .

[5]  I. Butera,et al.  Simultaneous identification of the pollutant release history and the source location in groundwater by means of a geostatistical approach , 2013, Stochastic Environmental Research and Risk Assessment.

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

[7]  P. Kitanidis,et al.  A geostatistical approach to contaminant source identification , 1997 .

[8]  A. Sun,et al.  A constrained robust least squares approach for contaminant release history identification , 2006 .

[9]  T. Ulrych,et al.  Extensions to minimum relative entropy inversion for noisy data. , 2003, Journal of contaminant hydrology.

[10]  T. Ulrych,et al.  Minimum Relative Entropy Inversion: Theory and Application to Recovering the Release History of a Groundwater Contaminant , 1996 .

[11]  P. Kitanidis,et al.  Geostatistical inversing for large-contrast transmissivity fields , 2009 .

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

[13]  P. Kitanidis,et al.  Aquifer heterogeneity characterization with oscillatory pumping: Sensitivity analysis and imaging potential , 2013 .

[14]  Solomon Kullback,et al.  Information Theory and Statistics , 1960 .

[15]  A. Zanini,et al.  Laboratory sandbox validation of pollutant source location methods , 2014, Stochastic Environmental Research and Risk Assessment.

[16]  T. Skaggs,et al.  Comment on “Minimum relative entropy inversion: Theory and application to recovering the release history of a groundwater contaminant” by Allan D. Woodbury and Tadeusz J. Ulrych , 1998 .

[17]  Allan D. Woodbury,et al.  Three-dimensional plume source reconstruction using minimum relative entropy inversion , 1998 .

[18]  Amvrossios C. Bagtzoglou,et al.  Pollution source identification in heterogeneous porous media , 2001 .

[19]  Amvrossios C. Bagtzoglou,et al.  Probabilistic Simulation for Reliable Solute Source Identification in Heterogeneous Porous Media , 1991 .

[20]  P. Kitanidis,et al.  Estimation of historical groundwater contaminant distribution using the adjoint state method applied to geostatistical inverse modeling , 2004 .

[21]  Allan D. Woodbury,et al.  Minimum relative entropy, Bayes and Kapur , 2011 .

[22]  Brian J. Wagner,et al.  Simultaneous parameter estimation and contaminant source characterization for coupled groundwater flow and contaminant transport modelling , 1992 .

[23]  Brian Borchers,et al.  Comparison of inverse methods for reconstructing the release history of a groundwater contamination source , 2000 .

[24]  I. Butera,et al.  Use of numerical modelling to identify the transfer function and application to the geostatistical procedure in the solution of inverse problems in groundwater , 2006 .

[25]  Stephan L. Mintz,et al.  Fundamental Theories in Physics , 2012, Studies in the Natural Sciences.

[26]  J. Bear,et al.  Modeling groundwater flow and pollution , 1987 .

[27]  Rodney W. Johnson Relative-Entropy Minimization with Uncertain Constraints: Theory and Application to Spectrum Analysis , 1987 .

[28]  I. Butera,et al.  A geostatistical approach to recover the release history of groundwater pollutants , 2003 .

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

[30]  A. Bagtzoglou,et al.  Application of particle methods to reliable identification of groundwater pollution sources , 1992 .

[31]  Brian Borchers,et al.  A MATLAB implementation of the minimum relative entropy method for linear inverse problems , 2001 .