A hybrid global optimization method for inverse estimation of hydraulic parameters: Annealing‐Simplex Method

Inverse estimation of unsaturated hydraulic parameters is often a highly nonlinear optimization problem with multiple parameters. The objective functions involved are often topographically complex and may contain many local minima. Because of these reasons, the inverse solutions are commonly very sensitive to the initial guess of the parameters when conventional optimizers are used. This paper presents an annealing-simplex method that incorporates simulated annealing strategies into a classical downhill simplex method. An upward infiltration experiment was used as an example of inverse estimation to test the method. Numerical experiments of both minimizing an algebraic function and inversion of upward infiltration data showed that the new method successfully converged to the global minimum in all cases, irrespective of the initial hydraulic parameter estimates, while the classical downhill method often converged to unfavorable local minima. The CPU times needed for the annealing-simplex method to estimate 5 and 7 hydraulic parameters simultaneously are about a half hour and 1 hour on a PC, respectively. Additionally, no special requirements are imposed on the objective function, and the method is independent of the details of the simulation submodel. Therefore the proposed method should be applicable to other optimization problems in water resources when it is important to have a robust global search capability.

[1]  G. Pope,et al.  Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing , 1995 .

[2]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[3]  Janusz S. Kowalik,et al.  Iterative methods for nonlinear optimization problems , 1972 .

[4]  K. Abbaspour,et al.  A sequential uncertainty domain inverse procedure for estimating subsurface flow and transport parameters , 1997 .

[5]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[6]  Antoine Mensch,et al.  An inverse technique for developing models for fluid flow in fracture systems using simulated annealing , 1993 .

[7]  Uri Shani,et al.  Analyses of infiltration events in relation to determining soil hydraulic properties by inverse problem methodology , 1991 .

[8]  Pengcheng Liu,et al.  An Improved Simulated Annealing-Downhill Simplex Hybrid Global Inverse Algorithm. , 1995 .

[9]  S. Hartzell,et al.  Calculation of earthquake rupture histories using a hybrid global search algorithm: application to the 1992 Landers, California, earthquake , 1996 .

[10]  Lehua Pan,et al.  A Transformed Pressure Head‐Based Approach to Solve Richards' Equation for Variably Saturated Soils , 1995 .

[11]  J. Parker,et al.  Analysis of the inverse problem for transient unsaturated flow , 1988 .

[12]  Robert A. Marryott,et al.  Optimal Groundwater Management: 1. Simulated Annealing , 1991 .

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

[14]  R. Hills,et al.  Parameter estimation of hydraulic properties from one‐step outflow data , 1992 .

[15]  Van Genuchten,et al.  A closed-form equation for predicting the hydraulic conductivity of unsaturated soils , 1980 .

[16]  M. T. van Genuchten,et al.  Estimating Unsaturated Soil Hydraulic Properties from Tension Disc Infiltrometer Data by Numerical Inversion , 1996 .

[17]  G. Parkin,et al.  Unsaturated Hydraulic Conductivity Measured by Time Domain Reflectometry Under a Rainfall Simulator , 1995 .

[18]  Jack C. Parker,et al.  Determining Soil Hydraulic Properties from One-step Outflow Experiments by Parameter Estimation: I. Theory and Numerical Studies1 , 1985 .

[19]  R. G. Hills,et al.  Unsaturated Hydraulic Properties from Upward Flow into Soil Cores , 1996 .

[20]  R. G. Hills,et al.  Modeling one‐dimensional infiltration into very dry soils: 1. Model development and evaluation , 1989 .