Multi-objective optimization in variably saturated fluid flow

A general purpose multi-objective inverse modeling strategy was developed and implemented to quantify fluid flow parameters in variably saturated porous materials. The strategy combines a robust and mass-conservative numerical simulator of the flow equation with an optimization algorithm and an experimental data set to estimate the parameters. The numerical simulator of the direct problem shows excellent agreement with a reference solution and conserves global mass with near perfection. An adaptive method was proposed in which the sensitivity matrix was calculated by one-sided finite difference approximation at the early stages of the optimization and the more accurate two-sided differentiation as the search approaches the minimum. A combined termination criterium was developed to stop the inverse code at the solution. The results of the multi-objective optimization were compared with those of single-objective minimization. While single-objective optimization generates reasonable results for either the fluid pressure head profile or the fluid content data, the proposed multi-objective optimization shows excellent agreement with both profiles simultaneously.

[1]  Kouroush Sadegh Zadeh.  MULTI-SCALE INVERSE MODELING IN BIOLOGICAL MASS TRANSPORT PROCESSES , 2006 .

[2]  R. Carsel,et al.  Developing joint probability distributions of soil water retention characteristics , 1988 .

[3]  Y. Mualem.  A New Model for Predicting the Hydraulic Conductivity , 1976 .

[4]  N Oreskes,et al.  Verification, Validation, and Confirmation of Numerical Models in the Earth Sciences , 1994, Science.

[5]  Cass T. Miller,et al.  Accurate and economical solution of the pressure-head form of Richards' equation by the method of lines , 1997 .

[6]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[7]  T. Brubaker,et al.  Nonlinear Parameter Estimation , 1979 .

[8]  Karsten Pruess,et al.  Robust numerical methods for saturated-unsaturated flow with dry initial conditions in heterogeneous media , 1995 .

[9]  P. Milly,et al.  A mass-conservative procedure for time-stepping in models of unsaturated flow , 1985 .

[10]  P. Sabatier,et al.  Past and future of inverse problems , 2000 .

[11]  Benny Malengier,et al.  Parameter identification in stationary groundwater flow problems in drainage basins , 2004 .

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

[13]  R. McCuen.  Modeling Hydrologic Change: Statistical Methods , 2002 .

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

[15]  J. Seinfeld,et al.  Estimation of absolute and relative permeabilities in petroleum reservoirs , 1987 .

[16]  Feike J. Leij,et al.  The RETC code for quantifying the hydraulic functions of unsaturated soils , 1992 .

[17]  L. A. Richards.  Capillary conduction of liquids through porous mediums , 1931 .

[18]  Aimin Yan,et al.  On the recovery of transport parameters in groundwater modelling , 2004 .

[19]  Jesús Carrera,et al.  Maximum-likelihood adjoint-state finite-element estimation of groundwater parameters under steady- and nonsteady-state conditions , 1985 .

[20]  A. Yan,et al.  On the recovery of multiple flow parameters from transient head data , 2004 .

[21]  Richard J. Cleary,et al.  Statistical Methods for Engineers , 1999 .

[22]  John H. Seinfeld,et al.  Estimation of two-phase petroleum reservoir properties by regularization , 1987 .

[23]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[24]  G. W. Snedecor.  Statistical Methods , 1964 .

[25]  Adel Shirmohammadi,et al.  Evaluation of infiltration models in contaminated landscape , 2007, Journal of environmental science and health. Part A, Toxic/hazardous substances & environmental engineering.

[26]  M. Celia,et al.  A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation , 1990 .

[27]  Adel Shirmohammadi,et al.  Identification of biomolecule mass transport and binding rate parameters in living cells by inverse modeling , 2006, Theoretical Biology and Medical Modelling.

[28]  M. R. Osborne.  Nonlinear least squares — the Levenberg algorithm revisited , 1976, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[29]  Adel Shirmohammadi,et al.  A finite element model for protein transport in vivo , 2007, Biomedical engineering online.

[30]  M. Th. van Genuchten,et al.  Parameter estimation for unsaturated flow and transport models — A review , 1987 .

[31]  T. Reilly,et al.  Flow and Storage in Groundwater Systems , 2002, Science.