Reduced order parameter estimation using quasilinearization and quadratic programming

[1] The ability of a particular model to accurately predict how a system responds to forcing is predicated on various model parameters that must be appropriately identified. There are many algorithms whose purpose is to solve this inverse problem, which is often computationally intensive. In this study, we propose a new algorithm that significantly reduces the computational burden associated with parameter identification. The algorithm is an extension of the quasilinearization approach where the governing system of differential equations is linearized with respect to the parameters. The resulting inverse problem therefore becomes a linear regression or quadratic programming problem (QP) for minimizing the sum of squared residuals; the solution becomes an update on the parameter set. This process of linearization and regression is repeated until convergence takes place. This algorithm has not received much attention, as the QPs can become quite large, often infeasible for real-world systems. To alleviate this drawback, proper orthogonal decomposition is applied to reduce the size of the linearized model, thereby reducing the computational burden of solving each QP. In fact, this study shows that the snapshots need only be calculated once at the very beginning of the algorithm, after which no further calculations of the reduced-model subspace are required. The proposed algorithm therefore only requires one linearized full-model run per parameter at the first iteration followed by a series of reduced-order QPs. The method is applied to a groundwater model with about 30,000 computation nodes where as many as 15 zones of hydraulic conductivity are estimated.

[1]  D. Oliver,et al.  Recent progress on reservoir history matching: a review , 2011 .

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

[3]  W. Yeh,et al.  Quasilinearization and the Identification of Aquifer Parameters , 1971 .

[4]  W. Yeh,et al.  Optimal Identification of Aquifer Diffusivity Using Quasilinearization , 1971 .

[5]  Richard L. Cooley,et al.  A Comparison of Several Methods of Solving Nonlinear Regression Groundwater Flow Problems , 1985 .

[6]  William W.-G. Yeh,et al.  Groundwater Management Using Model Reduction via Empirical Orthogonal Functions , 2008 .

[7]  R. Willis.  A planning model for the management of groundwater quality , 1979 .

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

[9]  Richard Bellman,et al.  Introduction to Matrix Analysis , 1972 .

[10]  J. Peraire,et al.  Balanced Model Reduction via the Proper Orthogonal Decomposition , 2002 .

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

[12]  H. Park,et al.  Estimation of parameters in flow reactors using the Karhunen-Loève decomposition , 1998 .

[13]  Arthur Veldman,et al.  Proper orthogonal decomposition and low-dimensional models for driven cavity flows , 1998 .

[14]  The influence of a confining layer on saltwater intrusion under surface recharge and groundwater extraction conditions , 2002 .

[15]  William W.-G. Yeh.  Optimal identification of parameters in an inhomogeneous medium with quadratic programming , 1975 .

[16]  Arnold Heemink,et al.  Reduced models for linear groundwater flow models using empirical orthogonal functions , 2004 .

[17]  Arnold W. Heemink,et al.  Inverse modeling of groundwater flow using model reduction , 2005 .

[18]  Adam J. Siade,et al.  Snapshot selection for groundwater model reduction using proper orthogonal decomposition , 2010 .

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

[20]  Mary C. Hill,et al.  UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API , 2006 .