Assessment of uncertainty associated with the estimation of well catchments by moment equations

Abstract Non-local stochastic moment equations are used successfully to analyze groundwater flow in randomly heterogeneous media. Here we present a moment equations-based approach to quantify the uncertainty associated with the estimation of well catchments. Our approach is based on the development of a complete second order formalism which allows obtaining the first statistical moments of the trajectories of conservative solute particles advected in a generally non-uniform groundwater flow. Approximate equations of moments of particles’ trajectories are then derived on the basis of a second order expansion in terms of the standard deviation of the aquifer log hydraulic conductivity. Analytical expressions are then obtained for the predictors of locations of mean stagnation points, together with their associated uncertainties. We implement our approach on heterogeneous media in bounded two-dimensional domains, with and without including the effect of conditioning on hydraulic conductivity information. The impact of domain size, boundary conditions, heterogeneity and non-stationarity of hydraulic conductivity on the prediction of a well catchment is explored. The results are compared against Monte Carlo simulations and semi-analytical solutions available in the literature. The methodology is applicable to both infinite and bounded domains and is free of distributional assumptions (and so applies to both Gaussian and non-Gaussian log hydraulic conductivity fields) and formally includes the effect of conditioning on available information.

[1]  Matthijs van Leeuwen,et al.  Stochastic determination of well capture zones , 1998 .

[2]  Dongxiao Zhang,et al.  Solute spreading in nonstationary flows in bounded, heterogeneous unsaturated‐saturated media , 2003 .

[3]  Matthijs van Leeuwen,et al.  Stochastic determination of well capture zones conditioned on regular grids of transmissivity measurements , 2000 .

[4]  Alberto Guadagnini,et al.  Travel time and trajectory moments of conservative solutes in three dimensional heterogeneous porous media under mean uniform flow , 2005 .

[5]  Dongxiao Zhang,et al.  On stochastic study of well capture zones in bounded, randomly heterogeneous media , 2003 .

[6]  Alberto Guadagnini,et al.  Time‐Related Capture Zones for Contaminants in Randomly Heterogeneous Formations , 1999 .

[7]  Alberto Guadagnini,et al.  Time-related capture zones for radial flow in two dimensional randomly heterogeneous media , 1999 .

[8]  W. Kinzelbach,et al.  Determination of a well head protection zone by stochastic inverse modelling , 1998 .

[9]  Alberto Guadagnini,et al.  Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains: 1. Theory and computational approach , 1999 .

[10]  Gour-Tsyh Yeh,et al.  On the computation of Darcian velocity and mass balance in the finite element modeling of groundwater flow , 1981 .

[11]  V. Cvetkovic,et al.  Modeling solute advection coupled with sorption kinetics in heterogeneous formations , 1992 .

[12]  Alberto Guadagnini,et al.  Probabilistic estimation of well catchments in heterogeneous aquifers , 1996 .

[13]  E. Scott Bair,et al.  A Monte Carlo-Based Approach for Determining Traveltime-Related Capture Zones of Wells Using Convex Hulls as Confidence Regions , 1991 .

[14]  Stephen E. Silliman,et al.  Capture Zones for Passive Wells in Heterogeneous Unconfined Aquifers , 1997 .

[15]  G. Dagan.  Solute transport in heterogeneous porous formations , 1984, Journal of Fluid Mechanics.

[16]  Alberto Guadagnini,et al.  Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion , 2003 .

[17]  Alberto Guadagnini,et al.  Mean travel time of conservative solutes in randomly heterogeneous unbounded domains under mean uniform flow , 2003 .

[18]  Wolfgang Kinzelbach,et al.  Determination of groundwater catchment areas in two and three spatial dimensions , 1992 .

[19]  Xian-Huan Wen,et al.  Analysis of nonlinear effects on tracer migration in heterogeneous aquifers using Lagrangian travel time statistics , 1996 .

[20]  G. Dagan,et al.  A solute flux approach to transport in heterogeneous formations: 1. The general framework , 1992 .

[21]  Wolfgang Kinzelbach,et al.  Semianalytical uncertainty estimation of well catchments: Conditioning by head and transmissivity data , 2004 .

[22]  Gedeon Dagan,et al.  Theory of Solute Transport by Groundwater , 1987 .

[23]  Gedeon Dagan,et al.  Solute transport in divergent radial flow through heterogeneous porous media , 1999, Journal of Fluid Mechanics.

[24]  Keith Beven,et al.  Stochastic capture zone delineation within the generalized likelihood uncertainty estimation methodology: Conditioning on head observations , 2001 .

[25]  Alberto Guadagnini,et al.  Convergence assessment of numerical Monte Carlo simulations in groundwater hydrology , 2004 .

[26]  Alberto Guadagnini,et al.  Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains: 2. Computational examples , 1999 .

[27]  Y. Rubin.  Stochastic modeling of macrodispersion in heterogeneous porous media , 1990 .

[28]  Daniel M. Tartakovsky,et al.  Nonlocal and localized analyses of conditional mean transient flow in bounded, randomly heterogeneous porous media , 2004 .

[29]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[30]  Bill X. Hu,et al.  A three-dimensional numerical method of moments for groundwater flow and solute transport in a nonstationary conductivity field , 2003 .

[31]  Alexander Y. Sun,et al.  Solute flux approach to transport through spatially nonstationary flow in porous media , 2000 .

[32]  Uncertainty estimation of well catchments in heterogeneous aquifers , 2002 .