Stochastic Solution to Inverse Problem in Ground Water
Publication: Journal of Hydraulic Engineering
Volume 123, Issue 12
Abstract
A solution to the inverse problem in groundwater is presented with a geostatistical framework, using Kalman filtering and a nonlinear gradient-based search technique. The Kalman filtering recursions are based on a newly developed and linear state-space equation that relates aquifer head perturbations to stochastic perturbations of log-aquifer properties and effective recharge. The Davidon-Fletcher-Powell (DFP) search algorithm is used to identify the mean and the variance of log-aquifer transmissivity and storativity, by minimizing the joint negative log-likelihood function of the innovations (prediction errors). Application to a numerical experiment indicates that the methodology performs well for log-transmissivity integral scale smaller than aquifer dimensions. The results underline the need for conditioning on point measurements of transmissivity and storativity, if the objective is to estimate the variance parameters. While head measurements are instrumental for estimating the geometric means (large-scale parameters), they, however, may not be sufficient for inferring the variance of log-aquifer properties (small-scale parameters).
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Aboufirassi, M., and Mariño, M. A.(1984). “Cokriging of aquifer transmissivity from field measurements of transmissivity and specific capacity.”Mathematical Geology, 16(1), 19–35.
2.
Andricevic, R.(1990). “Cost-effective network design for groundwater flow monitoring.”Stochastic Hydrol. Hydr., 4(1), 27–41.
3.
Bakr, A. A., Gelhar, L. W., Gutjahr, A. L., and McMillan, J. R.(1978). “Stochastic analysis of spatial variability in subsurface flow, 1, Comparison of one- and three-dimensional flows.”Water Resour. Res., 14(2), 263–272.
4.
Bras, R. L., and Rodríguez-Iturbe, I. (1985). Random functions and hydrology. Addison-Wesley Publishing Company, Inc., Reading, Mass.
5.
Carrera, J., and Glorioso, L.(1991). “On geostatistical formulation of the groundwater flow inverse problem.”Advanced Water Resour., 14(5), 273–280.
6.
Carrera, J., and Neuman, S. P.(1986). “Estimation of aquifer parameters under transient and steady state conditions, 1, Maximum likelihood method incorporating prior information.”Water Resour. Res., 22(2), 199–210.
7.
Clifton, P. M., and Neuman, S. P.(1982). “Effects of kriging and inverse modeling on conditional simulation of the Avra valley aquifer in southern Arizona.”Water Resour. Res., 18(4), 1215–1234.
8.
Cooley, R. L.(1977). “A method of estimating parameters and assessing reliability for models of steady state groundwater flow, 1, Theory and numerical properties.”Water Resour. Res., 13(2), 318–324.
9.
Dagan, G.(1985). “Stochastic modeling of groundwater flow by unconditional and conditional probabilities: The inverse problem.”Water Resour., 21(1), 65–72.
10.
Dagan, G. (1986). “Statistical theory of groundwater flow and transport: pore to laboratory, laboratory to formation and formation to regional scale.”Water Resour. Res., 22(9), 120S–135S.
11.
Dagan, G. and Rubin(1988). “Stochastic identification of recharge, transmissivity, and storativity in aquifer transient flow: a quasi-steady approach.”Water Resour. Res., 24(10), 1698–1710.
12.
Ginn, T. R., and Cushman, J. H.(1990). “Inverse methods for subsurface flow: a critical review of stochastic techniques.”Stochastic Hydrol. Hydr., 4, 1–26.
13.
Gorelick, S. M., Voss, C. I., Gill, P. E., Murray, W., Saunders, M. A., and Wright, M. H.(1984). “Aquifer reclamation design: the use of contaminant transport simulation combined with nonlinear programming.”Water Resour. Res., 20(4), 415–427.
14.
Gupta, N. K., and Mehra, R. K. (1974). “Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations.”IEEE Trans. Automatic Control, AC-19(6), 774–783.
15.
Hantush, M. M., and Mariño, M. A.(1994). “Two-dimensional stochastic analysis and optimal estimation in aquifers: random recharge.”Water Resour. Res., 30(2), 559–569.
16.
Hantush, M. M., and Mariño, M. A.(1995). “Continuous-time stochastic analysis of groundwater flow in heterogeneous aquifers.”Water Resour. Res., 31(3), 565–575.
17.
Hantush, M. M., and Mariño, M. A.(1997). “Estimation of spatially-variable aquifer hydraulic properties using Kalman filtering.”J. Hydr. Engrg., ASCE, 123(11), 1027–1035.
18.
Hoeksema, R. J. and Kitanidis(1984). “An application of the geostatistical approach to the inverse problem in two-dimensional groundwater modeling.”Water Resour. Res., 20(7), 1009–1020.
19.
Hoeksema, R. J., and Kitanidis, P. K.(1985). “Comparison of Gaussian conditional mean and kriging estimation in the geostatistical solution of the inverse problem.”Water Resour. Res., 21(6), 825–836.
20.
IMSL. (1984). International Mathematics and Statistics Libraries. IMSL, Houston, Tex.
21.
Jazwinski, A. H. (1970). Stochastic processes and filtering theory. Academic Press, Inc., New York, N.Y.
22.
Kitanidis, P. K., and Lane, R. W.(1985). “Maximum likelihood parameter estimation of hydrologic spatial processes by the Gauss-Newton method.”J. Hydrol., 79(1), 53–71.
23.
Kitanidis, P. K., and Vomvoris, E. G.(1983). “A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one dimensional simulations.”Water Resour. Res., 19(3), 677–690.
24.
Kuiper, L. K.(1973). “Analytical solution of spatially discretized groundwater flow equations.”Water Resour. Res., 9(4), 1094–1097.
25.
Lee, K. K., and Cushman, J. H.(1993). “Multiscale adaptive estimation of conductivity field using pump test data.”Stochastic Hydrol. Hydr, 7(3), 241–253.
26.
Loaiciga, H. A., and Mariño, M. A.(1987). “The inverse problem for confined aquifer flow: identification and estimation with extensions.”Water Resour. Res., 23(1), 92–104.
27.
Rubin, Y., and Dagan, G.(1987). “Stochastic identification of transmissivity and effective recharge in steady groundwater flow: 1. Theory.”Water Resour. Res., 23(7), 1185–1192.
28.
Sun, N.-Z., and Yeh, W. W.-G.(1992). “A stochastic inverse solution for transient groundwater flow: parameter identification and reliability analysis.”Water Resour. Res., 28(12), 3269–3280.
29.
Townley, L. R., and Wilson, J. L.(1985). “Computationally efficient algorithms for parameter estimation and uncertainty propagation in numerical models of groundwater flow.”Water Resour. Res., 21(12), 1851–1860.
30.
Vecchia, A. D., and Cooley, R. L.(1987). “Simultaneous confidence and prediction intervals for nonlinear regression models with application to a groundwater flow model.”Water Resour. Res., 23(7), 1237–1250.
31.
Yeh, W. W.-G.(1986). “Review of parameter identification procedures in groundwater hydrology: the inverse problem.”Water Resour. Res., 22(1), 95–108.
32.
Yeh, W. W.-G., and Yoon, Y. S.(1981). “Aquifer parameter identification with optimum dimension in parameterization.”Water Resour. Res., 17(3), 664–672.
33.
Yeh, W. W.-G., Yoon, Y. S., and Lee, K. S.(1983). “Aquifer parameter identification with kriging and optimum parameterization.”Water Resour. Res., 19(1), 225–233.
Information & Authors
Information
Published In
Journal of Hydraulic Engineering
Volume 123 • Issue 12 • December 1997
Pages: 1139 - 1146
Copyright
Copyright © 1997 American Society of Civil Engineers.
History
Published online: Dec 1, 1997
Published in print: Dec 1997
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.