A robust approach for iterative contaminant source location and release history recovery.

Contamination source identification is a crucial step in environmental remediation. The exact contaminant source locations and release histories are often unknown due to lack of records and therefore must be identified through inversion. Coupled source location and release history identification is a complex nonlinear optimization problem. Existing strategies for contaminant source identification have important practical limitations. In many studies, analytical solutions for point sources are used; the problem is often formulated and solved via nonlinear optimization; and model uncertainty is seldom considered. In practice, model uncertainty can be significant because of the uncertainty in model structure and parameters, and the error in numerical solutions. An inaccurate model can lead to erroneous inversion of contaminant sources. In this work, a constrained robust least squares (CRLS) estimator is combined with a branch-and-bound global optimization solver for iteratively identifying source release histories and source locations. CRLS is used for source release history recovery and the global optimization solver is used for location search. CRLS is a robust estimator that was developed to incorporate directly a modeler's prior knowledge of model uncertainty and measurement error. The robustness of CRLS is essential for systems that are ill-conditioned. Because of this decoupling, the total solution time can be reduced significantly. Our numerical experiments show that the combination of CRLS with the global optimization solver achieved better performance than the combination of a non-robust estimator, i.e., the nonnegative least squares (NNLS) method, with the same solver.

[1]  C. Zheng,et al.  A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems , 1997 .

[2]  J. Bear Hydraulics of Groundwater , 1979 .

[3]  L. J. Campbell,et al.  Groundwater Flow Systems , 2000 .

[4]  Michael J. Todd,et al.  Self-Scaled Barriers and Interior-Point Methods for Convex Programming , 1997, Math. Oper. Res..

[5]  M. G. Anderson Encyclopedia of hydrological sciences. , 2005 .

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

[7]  S. Gorelick,et al.  Identifying sources of groundwater pollution: An optimization approach , 1983 .

[8]  Bruce E. Hajek,et al.  Cooling Schedules for Optimal Annealing , 1988, Math. Oper. Res..

[9]  János D. Pintér,et al.  Global optimization in action , 1995 .

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

[11]  Dimitri P. Solomatine,et al.  Groundwater Remediation Strategy Using Global Optimization Algorithms , 2002 .

[12]  G. Golub,et al.  Efficient algorithms for least squares type problems with bounded uncertainties , 1997 .

[13]  Jos F. Sturm,et al.  A Matlab toolbox for optimization over symmetric cones , 1999 .

[14]  A. Mayer,et al.  Pump‐and‐treat optimization using well locations and pumping rates as decision variables , 1997 .

[15]  Tao Wang Global Optimization for Constrained Nonlinear Programming , 2000 .

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

[17]  A. Neumaier Acta Numerica 2004: Complete search in continuous global optimization and constraint satisfaction , 2004 .

[18]  Bithin Datta,et al.  Optimal Identification of Ground-Water Pollution Sources and Parameter Estimation , 2001 .

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

[20]  A. Bagtzoglou,et al.  State of the Art Report on Mathematical Methods for Groundwater Pollution Source Identification , 2001 .

[21]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

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

[23]  J. Guan,et al.  Optimal remediation with well locations and pumping rates selected as continuous decision variables , 1999 .

[24]  Federico Thomas,et al.  An ellipsoidal calculus based on propagation and fusion , 2002, IEEE Trans. Syst. Man Cybern. Part B.

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

[26]  Laurent El Ghaoui,et al.  Robust Solutions to Least-Squares Problems with Uncertain Data , 1997, SIAM J. Matrix Anal. Appl..

[27]  Masakazu Kojima,et al.  Exploiting sparsity in primal-dual interior-point methods for semidefinite programming , 1997, Math. Program..

[28]  N. Sun Mathematical Modeling of Groundwater Pollution , 1995 .

[29]  Arkadi Nemirovski,et al.  Robust Truss Topology Design via Semidefinite Programming , 1997, SIAM J. Optim..

[30]  Donald Goldfarb,et al.  Second-order cone programming , 2003, Math. Program..

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

[32]  Roseanna M. Neupauer,et al.  Adjoint method for obtaining backward‐in‐time location and travel time probabilities of a conservative groundwater contaminant , 1999 .

[33]  Arnold Neumaier,et al.  SNOBFIT -- Stable Noisy Optimization by Branch and Fit , 2008, TOMS.

[34]  N. Sun Inverse problems in groundwater modeling , 1994 .

[35]  Per Christian Hansen,et al.  REGULARIZATION TOOLS: A Matlab package for analysis and solution of discrete ill-posed problems , 1994, Numerical Algorithms.

[36]  You‐Kuan Zhang Stochastic Methods for Flow in Porous Media: Coping with Uncertainties , 2001 .

[37]  A. Kurzhanski,et al.  Ellipsoidal Calculus for Estimation and Control , 1996 .

[38]  Gene H. Golub,et al.  Matrix computations , 1983 .

[39]  G. Mahinthakumar,et al.  Hybrid Genetic Algorithm—Local Search Methods for Solving Groundwater Source Identification Inverse Problems , 2005 .

[40]  Melvyn Sim,et al.  The Price of Robustness , 2004, Oper. Res..

[41]  Mustafa M. Aral,et al.  Identification of Contaminant Source Location and Release History in Aquifers , 2001 .

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

[43]  Alexander Y. Sun,et al.  Inverse Methods for Parameter Estimations , 2006 .

[44]  A. Neumaier Complete search in continuous global optimization and constraint satisfaction , 2004, Acta Numerica.

[45]  T. Skaggs,et al.  Limitations in recovering the history of a groundwater contaminant plume , 1998 .

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

[47]  R. Neupauer,et al.  Adjoint‐derived location and travel time probabilities for a multidimensional groundwater system , 2001 .

[48]  S. P. Neuman,et al.  Combined Estimation of Hydrogeologic Conceptual Model and Parameter Uncertainty , 2004 .

[49]  Paul Sas,et al.  ROBUST DESIGN AND ROBUST STABILITY ANALYSIS OF ACTIVE NOISE CONTROL SYSTEMS , 2001 .

[50]  Stephen P. Boyd,et al.  Applications of second-order cone programming , 1998 .