A Simple Approach to Account for Radial Flow and Boundary Conditions When Kriging Hydraulic Head Fields for Confined Aquifers

The estimation and mapping of realistic hydraulic head fields, hence of flow paths, is a major goal of many hydrogeological studies. The most widely used method to obtain reliable head fields is the inverse approach. This approach relies on the numerical approximation of the flow equation and requires specifying boundary conditions and the transmissivity of each grid element. Boundary conditions are often unknown or poorly known, yet they impose a strong signature on the head fields obtained by inverse analysis. A simpler alternative to the inverse approach is the direct kriging of the head field using the measurements obtained at observation wells. The kriging must be modified to incorporate the available information. Use of the dual kriging formalism enables simultaneously estimating the head field, the aquifer mean transmissivity, and the regional hydraulic gradient from head data in steady or transient state conditions. In transient state conditions, an estimate of the storage coefficient can be obtained. We test the approach on simple analytical cases, on synthetic cases with solutions obtained numerically using a finite element flow simulator, and on a real aquifer. For homogeneous aquifers, infinite or bounded, the kriging estimate retrieves the exact solution of the head field, the exact hydrogeological parameters and the flow net. With heterogeneous aquifers, kriging accurately estimates the head field with prediction errors of the same magnitude as typical head measurement errors. The transmissivities are also accurately estimated by kriging. Moreover, if inversion is required, the kriged head along boundaries can be used as realistic boundary conditions for flow simulation.

[1]  C. E. Jacob,et al.  A generalized graphical method for evaluating formation constants and summarizing well‐field history , 1946 .

[2]  G. Marsily.  Quantitative Hydrogeology: Groundwater Hydrology for Engineers , 1986 .

[3]  J. H. Cushman,et al.  Role of boundary conditions in convergence and nonlocality of solutions to stochastic flow problems in bounded domains , 2000 .

[4]  Harald Osnes,et al.  Stochastic analysis of head spatial variability in bounded rectangular heterogeneous aquifers , 1995 .

[5]  Dongxiao Zhang.  Numerical solutions to statistical moment equations of groundwater flow in nonstationary, bounded, heterogeneous media , 1998 .

[6]  G. Dagan,et al.  Correlation structure of flow variables for steady flow toward a well with application to highly anisotropic heterogeneous formations , 1998 .

[7]  J. J. Royer,et al.  Dual Formalism of Kriging , 1984 .

[8]  J. Delhomme.  Kriging in the hydrosciences , 1978 .

[9]  J. Chilès,et al.  Geostatistics: Modeling Spatial Uncertainty , 1999 .

[10]  Gravity Data Transformation by Kriging , 1993 .

[11]  D. A. Zimmerman,et al.  A comparison of seven geostatistically based inverse approaches to estimate transmissivities for modeling advective transport by groundwater flow , 1998 .

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

[13]  S. P. Neuman,et al.  Role of Geostatistics in Subsurface Hydrology , 1984 .

[14]  X. Sanchez‐Vila.  Radially convergent flow in heterogeneous porous media , 1997 .

[15]  Frederick Delay,et al.  Empirical Orthogonal Functions Analysis Applied to the Inverse Problem in Hydrogeology: Evaluation of Uncertainty and Simulation of New Solutions , 2001 .

[16]  S. P. Neuman,et al.  Radial Flow in a Bounded Randomly Heterogeneous Aquifer , 2001 .

[17]  Geostatistics Applied to Hydrogeology’s Inverse Problem: Taking Boundary Conditions Into Account , 1997 .

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

[19]  M. W. Davis,et al.  Contouring Very Large Datasets Using Kriging , 1984 .

[20]  A. Maréchal.  Kriging Seismic Data in Presence of Faults , 1984 .

[21]  J. Bear,et al.  On the movement of water bodies injected into aquifers , 1965 .

[22]  Matthew J Tonkin,et al.  Kriging Water Levels with a Regional‐Linear and Point‐Logarithmic Drift , 2002, Ground water.