Abstract
Flow to wells in nonuniform geologic formations is of central interest to hydrogeologists and petroleum engineers. There are, however, very few mathematical analyses of such flow. We present analytical expressions for leading statistical moments of vertically averaged hydraulic head and flux under steady-state flow to a well that pumps water from a bounded, randomly heterogeneous aquifer. Like in the widely used Thiem equation, we prescribe a constant pumping rate deterministically at the well and a constant head at a circular outer boundary of radius L. We model the natural logarithm Y = lnT of aquifer transmissivity T as a statistically homogeneous random field with a Gaussian spatial correlation function. Our solution is based on exact nonlocal moment equations for multidimensional steady state flow in bounded, randomly heterogeneous porous media. Perturbation of these nonlocal equations leads to a system of local recursive moment equations that we solve analytically to second order in the standard deviation of Y. In contrast to most stochastic analyses of flow, which require that log transmissivity be multivariate Gaussian, our solution is free of any distributional assumptions. It yields expected values of head and flux, and the variance–covariance of these quantities, as functions of distance from the well. It also yields an apparent transmissivity, T a, defined as the negative ratio between expected flux and head gradient at any radial distance. The solution is supported by numerical Monte Carlo simulations, which demonstrate that it is applicable to strongly heterogeneous aquifers, characterized by large values of log transmissivity variance. The two-dimensional nature of our solution renders it useful for relatively thin aquifers in which vertical heterogeneity tends to be of minor concern relative to that in the horizontal plane. It also applies to thicker aquifers when information about their vertical heterogeneity is lacking, as is commonly the case when measurements of head and flow rate are done in wells that penetrate much of the aquifer thickness. Potential uses include the analysis of pumping tests and tracer test conducted in such wells, the statistical delineation of their respective capture zones, and the analysis of contaminant transport toward fully penetrating wells.
Similar content being viewed by others
References
Ababou, R. and Wood, E. F.: 1990, Comment on ‘Effective groundwater model parameter values: Influence of spatial variability of hydraulic conductivity, leakance and recharge’ by J. J. Gomez-Hernandez and S. M. Gorelick, Water Resour. Res. 26, 1643-1846.
Abramowitz, M. and Stegun, I. A.: 1972, Handbook of Mathematical Functions, Dover.
Butler, J. J.: 1991, A stochastic analysis of pumping tests in laterally nonuniform media. Water Resour. Res. 27, 2401-2414.
Dagan, G.: 1982, Stochastic modeling of groundwater flow by unconditional and conditional probabilities, 1. Conditional simulation and the direct problem. Water Resour. Res. 18, 813-833.
Dagan, G.: 1989, Flow and Transport in Porous Formations, Springer.
Dagan, G. and Neuman, S. P. (eds): 1997, Subsurface Flow and Transport: A Stochastic Approach, Cambridge University Press.
Desbarats, A. J.: 1992, Spatial averaging of transmissivity in heterogeneous fields with flow toward a well, Water Resour. Res. 28, 757-767.
Desbarats, A. J.: 1994, Spatial averaging of hydraulic conductivities under radial flow conditions. Math. Geol. 26, 1-21.
Deutsch, C. V. and Journel, A. G.: 1998, GSLIB Geostatistical Software Library and User's Guide, Oxford University Press.
Fiori, A., Indelman, P. and Dagan, G.: 1998, Correlation structure of flow variables for steady flow toward a well with application to highly anisotropic heterogeneous formations, Water Resour. Res. 34, 699-708.
Franzetti, S., Guadagnini, A. and Orsi, E.: 1996, Monte Carlo simulation and effective conductivity in confined radial flow fields, In: Kovar, K. and Van der Heijde, P. (eds.), Calibration and Reliability in Groundwater Modelling–ModelCARE 96, IAHS Publication No. 237, pp. 463-471.
Franzetti, S., Guadagnini, A., and Riva, M.: 1997, Conduttivit`a idraulica di un campo confinato eterogeneo con flusso mediamente radiale (in Italian), In: Proc. IX Nat. Congress of Geology, Vol. 1, Flash Print SAS, pp. 405-409.
Gomez-Hernandez, J. J. and Gorelick, S. M.: 1989, Effective groundwater model parameter values: Influence of spatial variability of hydraulic conductivity, leakance and recharge, Water Resour. Res. 25, 405-419.
Gradshteyn, I. S. and Ryzhik, I. M.: 1994, Table of Integrals, Series, and Products, Academic.
Guadagnini, A. and Neuman, S. P.: 1997, Nonlocal and localized finite element solution of conditional mean flow in randomly heterogeneous media, Tech. Rep. HWR-97-100, Department of Hydrology and Water Resources, The University of Arizona, Tucson.
Guadagnini, A. and Neuman, S. P.: 1999a, Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains, 1, theory and computational approach, Water Resour. Res. 35, 2999-3018.
Guadagnini, A. and Neuman, S. P.: 1999b, Nonlocal and localized analyses of conditional mean steady state flow in bounded, randomly nonuniform domains, 2, computational examples, Water Resour. Res. 35, 3019-3039.
Guadagnini, A. and Neuman, S. P.: 2001, Recursive conditional moment equations for advective transport in randomly heterogeneous velocity fields, Transport in Porous Media 42(1/2), 37-67.
Indelman, P. and Dagan, G.: 1999, Solute transport in divergent radial flow through heterogeneous porous media. J. Fluid Mech. 384, 159-182.
Indelman, P. and Abramovich, B.: 1994, Nonlocal properties of nonuniform averaged flows in heterogeneous media, Water Resour. Res. 30, 3385-3393.
Indelman, P., Fiori, A. and Dagan, G.: 1996, Steady flow toward wells in heterogeneous formations: mean head and equivalent conductivity. Water Resour. Res. 32, 1975-1983.
de Marsily, G.: 1986, Quantitative Hydrogeology, Academic Press, New York, 440 pp.
Matheron, G.: 1967, Elements Pour Une Theorie des Milieux Poreux, Masson et Cie.
McDonald, M. G. and Harbaugh, A. W.: 1988, A modular three-dimensional finite-difference groundwater flow model, Manual 83-875. U.S. Geological Survey.
Naff, R. L.: 1991, Radial flow in heterogeneous porous media: An analysis of specific discharge, Water Resour. Res. 27, 307-316.
Neuman, S. P.: 1993, Eulerian-Lagrangian theory of transport in space-time nonstationary velocity fields: exact nonlocal formalism by conditional moments and weak approximation, Water Resour. Res. 29, 633-645.
Neuman, S. P. and Orr, S.: 1993, Prediction of Steady State Flow in Nonuniform Geologic Media by Conditional Moments: Exact Nonlocal Formalism, Effective Conductivities, and Weak Approximation, Water Resour. Res. 29, 341-364.
Neuman, S. P., Tartakovsky, D. M., Wallstrom, T. C. and Winter, C. L.: 1996, Correction to ‘Prediction of steady state flow in nonuniform geologic media by conditional moments: Exact nonlocal formalism, effective conductivities, and weak approximation’, Water Resour. Res. 32, 1479-1480.
Oliver, D. S.: 1993, The influence of nonuniform transmissivity and storativity on drawdown, Water Resour. Res. 29, 169-178.
Osnes, H.: 1995, Stochastic analysis of head spatial variability in bounded rectangular heterogeneous aquifers. Water Resour. Res. 31, 2981-2990.
Osnes, H.: 1998, Stochastic analysis of velocity spatial variability in bounded rectangular heterogeneous aquifers, Adv Water Res. 21, 203-215.
Riva, M.: 1999, Modellazione stocastica di flussi radiali bi-e tri-dimensionali in formazioni porose eterogenee (in italian), Ph.D. Thesis. Politecnico di Milano.
Riva, M., Guadagnini, A., Tartakovky, D. M. and Neuman, S. P.: 1999, Three-dimensional steady state flow to a well in a randomly heterogeneous aquifer, In: Stauffer, F., Kinzedlbach, W., Kovar, K. and Hoehn, E. (eds), Calibration and Reliability in Groundwater Modelling, ModelCARE99, IAHS Publication No. 265, pp. 131-136.
Sánchez-Vila, X.: 1997, Radially convergent flow in heterogeneous porous media,Water Resour. Res. 33, 1633-1641.
Shvidler, M. I.: 1964, Filtration Flows in Heterogeneous Media (A Statistical Approach), Authorized Translation from the Russian, Consultants Bureau Enterprises.
Tartakovsky, D. M. and Neuman, S. P.: 1998, Transient flow in bounded randomly heterogeneous domains: 1. exact conditional moment equations and recursive approximations, Water Resour. Res. 34, 1-12.
Tartakovsky, D. M. and Neuman, S. P.: 1999, Extension of ‘Transient flow in bounded randomly heterogeneous domains, 1. Exact conditional moment equations and recursive approximations’, Water Resour. Res. 35, 1921-1925.
Weinberger, H. F.: 1965, A First Course in Partial Differential Equations: with Complex Parameter and Transform Methods, Wiley.
Zhang, D. and Neuman, S. P.: 1996, Effect of local dispersion on solute transport in randomly heterogeneous media, Water Resour. Res. 32, 2715-2723.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Riva, M., Guadagnini, A., Neuman, S.P. et al. Radial Flow in a Bounded Randomly Heterogeneous Aquifer. Transport in Porous Media 45, 139–193 (2001). https://doi.org/10.1023/A:1011880602668
Issue Date:
DOI: https://doi.org/10.1023/A:1011880602668