A multipopulation genetic algorithm to solve the inverse problem in hydrogeology

The inverse problem of groundwater flow is treated with an automatic method that can produce several alternative solutions at once. During their joint optimization, these solutions can exchange information in order to maintain some diversity and thus avoid a systematic premature convergence toward a single local minimum. Although genetic algorithms are capable of doing this, they have not often been used in groundwater inverse optimization. First, a specific multipopulation genetic algorithm is developed. It is then tested on two synthetic cases of steady state flow with transmissivity values extending over 4 orders of magnitude. The first test is nonparametric and optimizes as many parameters as those used to define the reference case. The second test uses a sort of “pilot point” parametrization. The optimization is carried out on a limited number of perturbations that are interpolated and superimposed on an initial transmissivity field. In view of the good quality of the results, these initial attempts provide incentives to further develop genetic algorithms in groundwater inverse problems.

[1]  Heinz Mühlenbein,et al.  The parallel genetic algorithm as function optimizer , 1991, Parallel Comput..

[2]  David W. Pollock,et al.  A Controlled Experiment in Ground Water Flow Model Calibration , 1998 .

[3]  Jens Lienig,et al.  A parallel genetic algorithm for performance-driven VLSI routing , 1997, IEEE Trans. Evol. Comput..

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

[5]  S. P. Neuman,et al.  Estimation of aquifer parameters under transient and steady-state conditions: 2 , 1986 .

[6]  A. Simpson,et al.  An Improved Genetic Algorithm for Pipe Network Optimization , 1996 .

[7]  Alexandru Agapie,et al.  Modelling Genetic Algorithms: From Markov Chains to Dependence with Complete Connections , 1998, PPSN.

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

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

[10]  Mary C. Hill,et al.  Hydrogeologic Evaluation and Numerical Simulation of the Death Valley Regional Ground-Water Flow System, Nevada and California , 1998 .

[11]  Alexandru Agapie,et al.  Genetic Algorithms: Minimal Conditions for Convergence , 1997, Artificial Evolution.

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

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

[14]  R. P. Oliveira,et al.  Operating rules for multireservoir systems , 1997 .

[15]  Frederick Delay,et al.  Transmissivity and morphological features in a chalk aquifer: a geostatistical approach of their relationship , 1997 .

[16]  Ralf Salomon,et al.  Evolutionary algorithms and gradient search: similarities and differences , 1998, IEEE Trans. Evol. Comput..

[17]  Richard L. Cooley,et al.  Nonlinear‐regression groundwater flow modeling of a deep regional aquifer system , 1986 .

[18]  Thomas C. Peachey,et al.  The Nature of Mutation in Genetic Algorithms , 1995, ICGA.

[19]  Alex S. Mayer,et al.  Development and application of a coupled-process parameter inversion model based on the maximum likelihood estimation method , 1999 .

[20]  Lothar Thiele,et al.  A Mathematical Analysis of Tournament Selection , 1995, ICGA.

[21]  Jagath J. Kaluarachchi,et al.  Application of artificial neural network and genetic algorithm in flow and transport simulations , 1998 .

[22]  P. C. Shah,et al.  Reservoir History Matching by Bayesian Estimation , 1976 .

[23]  P. Khanna,et al.  Genetic algorithm for optimization of water distribution systems , 1999, Environ. Model. Softw..

[24]  L. Townley,et al.  Computationally Efficient Algorithms for Parameter Estimation and Uncertainty Propagation in Numerical Models of Groundwater Flow , 1985 .

[25]  K. Katsifarakis,et al.  Minimization of pumping cost in zoned aquifers by means of genetic algorithms , 1998 .

[26]  J. Kaluarachchi,et al.  Parameter estimation using artificial neural network and genetic algorithm for free‐product migration and recovery , 1998 .

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

[28]  G. Newsam,et al.  Sufficient conditions for identifying transmissivity in a confined aquifer , 1990 .

[29]  Catherine Certes,et al.  Application of the pilot point method to the identification of aquifer transmissivities , 1991 .

[30]  P. Kitanidis,et al.  Maximum likelihood parameter estimation of hydrologic spatial processes by the Gauss-Newton method , 1985 .

[31]  Shumeet Baluja,et al.  Structure and Performance of Fine-Grain Parallelism in Genetic Search , 1993, ICGA.

[32]  David Levine,et al.  Application of a hybrid genetic algorithm to airline crew scheduling , 1996, Comput. Oper. Res..

[33]  L. L. Rogers,et al.  Optimization of groundwater remediation using artificial neural networks with parallel solute transport modeling , 1994 .

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

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

[36]  Brian J. Wagner,et al.  Sampling Design Methods For Groundwater Modeling Under Uncertainty , 1995 .

[37]  G. Chavent,et al.  History Matching by Use of Optimal Theory , 1975 .

[38]  Dana S. Richards,et al.  Distributed genetic algorithms for the floorplan design problem , 1991, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[39]  E. G. Vomvoris,et al.  A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulations , 1983 .

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

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

[42]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 3. Application to Synthetic and Field Data , 1986 .

[43]  Zbigniew Michalewicz,et al.  Evolutionary Algorithms in Engineering Applications , 1997, Springer Berlin Heidelberg.

[44]  Mary C. Hill,et al.  Hydrogeologic evaluation and numerical simulation of the Death Valley regional ground-water flow system, Nevada and California , 1997 .

[45]  José Carlos Príncipe,et al.  A Markov Chain Framework for the Simple Genetic Algorithm , 1993, Evolutionary Computation.

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