Inverse procedure for calibrating parameters that control land subsidence caused by subsurface fluid withdrawal: 1. Methods

[1] Land subsidence prediction depends on a selected theoretical model and the calibration of its parameter values. An inverse model is developed that calibrates five parameters of a compacting confined aquifer system: K′v (vertical hydraulic conductivity of the aquitards), S′skv (aquitard nonrecoverable skeletal specific storage), S′ske (aquitard recoverable skeletal specific storage), Ssk (skeletal specific storage of the aquifer), and p′max 0 (initial preconsolidation stress within the aquitards). This inverse model combines a procedure for finding an initial set of parameter values and an inverse adjustment procedure. Initial values are estimated by using a new graphic method for the analysis of field data. The follow-up inverse procedure is a computer algorithm that combines the Newton-Raphson method with Helm's one-dimensional finite difference subsidence model (COMPAC). COMPAC offers two options: a constant-parameter simulation and a stress-dependent parameter simulation. The classical least squares criterion is employed as the objective function of this inverse procedure. Three characteristics of the objective function are considered. This inverse model is applied to calculations of an idealized compacting confined aquifer system using Helm's model. This investigation shows that the final values for the five parameters lead to simulated compaction error of less than 1%. This contrasts to the initial error of more than 60%. These results indicate that the proposed inverse procedure (model) can be applied to actual field data.

[1]  I. Khan.  Inverse Problem in Ground Water: Model Application , 1986 .

[2]  Donald C. Helm,et al.  One‐dimensional simulation of aquifer system compaction near Pixley, California: 1. Constant parameters , 1975 .

[3]  David E. Prudic,et al.  Documentation of a computer program to simulate aquifer-system compaction using the modular finite d , 1988 .

[4]  D. C. Helm.  Field Verification of a One-Dimensional Mathematical Model for Transient Compaction and Expansion of a Confined Aquifer System , 1978 .

[5]  William W.-G. Yeh,et al.  Coupled inverse problems in groundwater modeling: 2. Identifiability and experimental design , 1990 .

[6]  Devin L. Galloway,et al.  Land subsidence in the United States , 1999 .

[7]  Giuseppe Gambolati,et al.  Mathematical simulation of the subsidence of Venice: 1 , 1973 .

[8]  D. C. Helm.  Field-based computational techniques for predicting subsidence due to fluid withdrawal , 1984 .

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

[10]  Giuseppe Gambolati,et al.  Mathematical simulation of the subsidence of Venice: 2. Results , 1974 .

[11]  E. Poeter,et al.  Documentation of UCODE; a computer code for universal inverse modeling , 1998 .

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

[13]  Donald C. Helm,et al.  One‐dimensional simulation of aquifer system compaction near Pixley, California: 2. Stress‐Dependent Parameters , 1976 .

[14]  Irfan A. Khan,et al.  Inverse Problem in Ground Water: Model Development , 1986 .

[15]  P. Witherspoon,et al.  The role of aquitards in multiple‐aquifer systems , 1972 .

[16]  William W.-G. Yeh,et al.  Aquifer parameter identification with optimum dimension in parameterization , 1981 .

[17]  Jackie P. Hallberg,et al.  Parameter Estimation Tools for Hydrologic and Hydraulic Simulations , 2006 .

[18]  Miguel A. Mariño,et al.  Prediction of optimal safe ground water yield and land subsidence in the Los Banos-Kettleman City area, California, using a calibrated numerical simulation model , 2001 .

[19]  Kenneth W. Hudnut,et al.  Detection of aquifer system compaction and land subsidence using interferometric synthetic aperture radar, Antelope Valley, Mojave Desert, California , 1998 .

[20]  Andres Alcolea,et al.  Inverse problem in hydrogeology , 2005 .

[21]  William W.-G. Yeh,et al.  Coupled inverse problems in groundwater modeling - 1. Sensitivity analysis and parameter identification. , 1990 .

[22]  Yi Liu,et al.  Inverse procedure for calibrating parameters that control land subsidence caused by subsurface fluid withdrawal: 2. Field application , 2008 .

[23]  T. Burbey.  Stress‐Strain Analyses for Aquifer‐System Characterization , 2001 .

[24]  Howard A. Zebker,et al.  Inverse modeling of interbed storage parameters using land subsidence observations, Antelope Valley, California , 2003 .

[25]  W. T. Pecora,et al.  Physical and hydrologic properties of water-bearing deposits in subsiding areas in central California , 1966 .

[26]  Paul J. Ossenbruggen.  Systems analysis for civil engineers , 1984 .

[27]  S. A. Leake,et al.  MODFLOW-2000 Ground-Water Model?User Guide to the Subsidence and Aquifer-System Compaction (SUB) Package , 2003 .