Estimation of co-conditional moments of transmissivity, hydraulic head, and velocity fields

An iterative co-conditional Monte Carlo simulation (IMCS) approach is developed. This approach derives co-conditional means and variances of transmissivity (T), head (f), and Darcy’s velocity ( q), based on sparse measurements of T and f in heterogeneous, confined aquifers under steady-state conditions. It employs the classical co-conditional Monte Carlo simulation technique (MCS) and a successive linear estimator that takes advantage of our prior knowledge of the covariances of T and f and their cross-covariance. In each co-conditional simulation, a linear estimate of T is improved by solving the governing steady-state flow equation, and by updating residual covariance functions iteratively. These residual covariance functions consist of the covariance of T and f and the cross-covariance function between T and f .A s a result, the non-linear relationship between T and f is incorporated in the co-conditional realizations of T and f. Once the T and f fields are generated, a corresponding velocity field is also calculated. The average of the co-conditioned realizations of T, f, and q yields the co-conditional mean fields. In turn, the co-conditional variances of T, f, and q, which measure the reduction in uncertainty due to measurements of T and f, are derived. Results of our numerical experiments show that the co-conditional means from IMCS for T and f fields have smaller mean square errors (MSE) than those from a non-iterative Monte Carlo simulation (NIMCS). Finally, the co-conditional mean fields from IMCS are compared with the co-conditional effective fields from a direct approach developed by Yehet al .[ Water Resources Research, 32(1), 85‐92, 1996]. q 1998 Elsevier Science Limited.

[1]  Allan L. Gutjahr,et al.  Co-kriging for stochastic flow models , 1989 .

[2]  G. Dagan,et al.  Stochastic identification of transmissivity and effective recharge in steady groundwater flow: 2. Case study , 1987 .

[3]  P. Kitanidis,et al.  An Application of the Geostatistical Approach to the Inverse Problem in Two-Dimensional Groundwater Modeling , 1984 .

[4]  Allan L. Gutjahr,et al.  An Iterative Cokriging‐Like Technique for Ground‐Water Flow Modeling , 1995 .

[5]  G. Dagan Stochastic Modeling of Groundwater Flow by Unconditional and Conditional Probabilities: The Inverse Problem , 1985 .

[6]  S. P. Neuman A statistical approach to the inverse problem of aquifer hydrology: 3. Improved solution method and added perspective , 1980 .

[7]  S. P. Neuman,et al.  Estimation of aquifer parameters under transient and steady-state conditions: 2 , 1986 .

[8]  Minghui Jin,et al.  AN ITERATIVE STOCHASTIC INVERSE METHOD: CONDITIONAL EFFECTIVE TRANSMISSIVITY AND HYDRAULIC HEAD FIELDS , 1995 .

[9]  Michael D. Dettinger,et al.  First order analysis of uncertainty in numerical models of groundwater flow part: 1. Mathematical development , 1981 .

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

[11]  William W.-G. Yeh,et al.  A stochastic inverse solution for transient groundwater flow: Parameter identification and reliability analysis , 1992 .

[12]  Samuel Kamel. Hanna An iterative Monte Carlo technique for estimating conditional means and variances of transmissivity and hydraulic head fields. , 1995 .

[13]  T. Harter,et al.  Conditional stochastic analysis of solute transport in heterogeneous, variably saturated soils , 1996 .

[14]  S. P. Neuman,et al.  A statistical approach to the inverse problem of aquifer hydrology: 1. Theory , 1979 .

[15]  Allan L. Gutjahr,et al.  Stochastic analysis of spatial variability in two‐dimensional steady groundwater flow assuming stationary and nonstationary heads , 1982 .

[16]  G. Dagan Stochastic modeling of groundwater flow by unconditional and conditional probabilities: 1. Conditional simulation and the direct problem , 1982 .

[17]  P. Kitanidis Quasi‐Linear Geostatistical Theory for Inversing , 1995 .

[18]  S. P. Neuman,et al.  Effects of kriging and inverse modeling on conditional simulation of the Avra Valley Aquifer in southern Arizona , 1982 .

[19]  R. Willis,et al.  Groundwater Systems Planning and Management , 1987 .

[20]  J. Carrera,et al.  On geostatistical formulations of the groundwater flow inverse problem , 1991 .

[21]  Richard L. Cooley,et al.  Incorporation of prior information on parameters into nonlinear regression groundwater flow models: 1. Theory , 1982 .

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

[23]  Peter K. Kitanidis,et al.  Prediction of transmissivities, heads, and seepage velocities using mathematical modeling and geostatistics , 1989 .

[24]  Sean Hatch,et al.  Joint conditional simulations and the spectral approach for flow modeling , 1994 .

[25]  C. R. Dietrich,et al.  A stability analysis of the geostatistical approach to aquifer transmissivity identification , 1989 .

[26]  E. G. Vomvoris,et al.  A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulations , 1983 .

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

[28]  R. Clapp,et al.  Calibration of groundwater flow models using Monte Carlo simulations and geostatistics. , 1990 .

[29]  P. C. Shah,et al.  Reservoir History Matching by Bayesian Estimation , 1976 .

[30]  T. Harter,et al.  A Numerical Model for Water Flow and Chemical Transport in Variably Saturated Porous Media , 1993 .