Effect of Conditioning Randomly Heterogeneous Transmissivity on Temporal Hydraulic Head Measurements in Transient Two-Dimensional Aquifer Flow

Abstract We investigated the effect of conditioning transient, two-dimensional groundwater flow simulations, where the transmissivity was a spatial random field, on time dependent head data. The random fields, representing perturbations in log transmissivity, were generated using a known covariance function and then conditioned to match head data by iteratively cokriging and solving the flow model numerically. A new approximation to the cross-covariance of log transmissivity perturbations with time dependent head data and head data at different times, that greatly increased the computational efficiency, was introduced. The most noticeable effect of head data on the estimation of head and log transmissivity perturbations occurred from conditioning only on spatially distributed head measurements during steady flow. The additional improvement in the estimation of the log transmissivity and head perturbations obtained by conditioning on time dependent head data was fairly small. On the other hand, conditioning on temporal head data had a significant effect on particle tracks and reduced the lateral spreading around the center of the paths.

[1]  A. Mantoglou,et al.  The Turning Bands Method for simulation of random fields using line generation by a spectral method , 1982 .

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

[3]  John L. Wilson,et al.  Ground‐Water Quality Modelinga , 1974 .

[4]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

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

[6]  Allan L. Gutjahr,et al.  Stochastic models of subsurface flow: Log linearized Gaussian models are “exact” for covariances , 1984 .

[7]  Peter K. Kitanidis,et al.  Analysis of the Spatial Structure of Properties of Selected Aquifers , 1985 .

[8]  John L. Wilson,et al.  Efficient and accurate front tracking for two‐dimensional groundwater flow models , 1991 .

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

[10]  R. Freeze.  A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media , 1975 .

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

[12]  John L. Wilson,et al.  Illustration and Verification of Adjoint Sensitivity Theory for Steady State Groundwater Flow , 1985 .

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

[14]  T.-C. Jim Yeh,et al.  An iterative geostatistical inverse method for steady flow in the vadose zone , 1996 .

[15]  P. Kitanidis.  On the geostatistical approach to the inverse problem , 1996 .

[16]  A. Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data—I. Theory , 1997 .

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

[18]  Gedeon Dagan,et al.  A Note on Higher‐Order Corrections of the Head Covariances in Steady Aquifer Flow , 1985 .

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

[20]  M. Mariño,et al.  Continuous Time Stochastic Analysis of Groundwater Flow in Heterogeneous Aquifers , 1995 .