Joint estimation of transmissivities and recharges—application: stochastic characterization of well capture zones

This paper is concerned about aquifer models that are affected by the uncertainty on hydraulic conductivity and recharge. In particular, we study the impact of these sources of uncertainty on the stochastic characterization of drinking water well capture zones. The estimation of a drinking water well capture zone is important for the design of a protection zone around the well. The estimation is associated with uncertainty due to imperfect knowledge of the input parameters of the groundwater flow and mass transport equations. A first synthetic study shows that the uncertainty of the mean recharge has a larger impact on the well catchment uncertainty than the uncertainty of the spatial variability of recharge. Therefore, it is concluded that in practice the recharge can be estimated for one zone or a limited number of zones, without the need to estimate the recharge on the grid cell scale. The sequential self-calibrated method for the Monte Carlo type inverse modelling of groundwater flow and mass transport was extended for the estimation of recharge rates. Recharge rates can be estimated together with hydraulic conductivities, and possibly other parameters, using both head and concentration data. In a second synthetic study the joint estimation of transmissivities and recharge rates is illustrated. It is found that the transmissivity, hydraulic head and recharge fields are better characterised after the inverse conditioning. In order to stabilize the inverse estimation, it is essential that a regularization term is included in the objective function in case parameters are estimated zone-wise.

[1]  A. Soares,et al.  Geostatistics Tróia '92 , 1993 .

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

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

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

[5]  Stochastic delineation of well capture zones , 2004 .

[6]  Andrés Sahuquillo,et al.  Joint simulation of transmissivity and storativity fields conditional to steady-state and transient hydraulic head data , 1999 .

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

[8]  Andrés Sahuquillo,et al.  Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data , 2003 .

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

[10]  J. C. Ramírez,et al.  Estimation of aquifer parameters under transient and steady-state conditions , 1984 .

[11]  Andrés Sahuquillo,et al.  Stochastic simulation of transmissivity fields conditional to both transmissivity and piezometric data 2. Demonstration on a synthetic aquifer , 1997 .

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

[13]  W. R. R. Blain,et al.  Hydraulic Engineering Software IV , 1990 .