Confined aquifer flow modeling whithout temporary discretization
DOI:
https://doi.org/10.21701/bolgeomin.119.2.006Keywords:
confined aquifer, flow model, groundwater, Laplace transform, StehfestAbstract
A interactive program for flow modeling in confined homogeneous or heterogeneous aquifer is presented (MGAC). The equation from which the program has been elaborated is the one obtained by applying the Laplace transform to the differential flow equation on the two directions of the plane. The models spatially discretize the aquifer in square cells and pose the system of transformed equations in the form of finite differences, the time variable not being present. The SOR method is used to solve the system. The solution obtained is numerically inverted to obtain the values of the drawdown, si, instead of their transform, s̄i, and to recover the time variable instead of the inversion parameter p. The inversion method used is that of Stehfest (1970). The application of time is performed once and in its totality, so that the models do not need time discretization. The use of the Stehfest method entails an error for small values of Tt/r2S. The validity of the method for homogeneous confined aquifer is approximated by comparing the values obtained with the inversion of the transformed analytic Theis equation and the own Theis expression. The analytical results are also compared with those obtained for a generical model. The truncation errors and the limitations originated by the lack of precision in the initial potentials or in the fixed real potentials in transient flow are reviewed. The required memory and the computer times for different examples are shown.
Downloads
References
Ahlfeld D.P. and Mulligan A.E. 2000. Optimal Management of Flow in Ground water Systems. Academic Press. San Diego.
Andrés M.F. 2003. Modelos de flujo y de gestión hidraúlica de acuíferos. Simulación cuasitridimensional sin iteraciones temporales. Tesis Doctoral. Universidad de Salamanca. 625 p.
Cheng A.H-D. and Ou K. 1989. An efficient Laplace transform solution for multiaquifer systems. Water Resources Research, 25(4), 742-748. https://doi.org/10.1029/WR025i004p00742
Cheng A.H-D. and Morohunfola O. 1990. Boundary element formulation for multiaquifer systems. In: Computational Engineering with Boundary Elements, vol. I. Fluid and Potential Problems, edited by S. Grilli, C. A. Brebbia, and A. H-D. Cheng, p. 145-156, Computational Mechanics Publications, Billerica, Mass.
Cheng A.H-D. and Morohunfola O.K. 1993a. Multilayered leaky aquifer systems. I. Pumping Well Solution. Water Resources Research, 29(8), 2787-2800. https://doi.org/10.1029/93WR00768
Cheng A.H-D. and Morohunfola O.K. 1993b. Multilayered leaky aquifer systems, II, Boundary element solutions, Water Resources Research, 29(8), 2801-2811. https://doi.org/10.1029/93WR00769
Cheng A.H.D., Sidauruk P. and Abousleiman Y. 1994. Aproximate Inversion of the Laplace Transform. The Mathematica Journal, 4(2), 76-82.
De Marsily G., Delay F., Gonçalves J., Renard P., Teles V. and Violette S. 2005. Dealing with spacial heterogeneity. Hydrogeology Journal, 13(1), 161- 183. https://doi.org/10.1007/s10040-004-0432-3
Desbarats A.J. 1992. Spatial averaging of the transmissivity in heterogeneous fields with flow toward a well. Water Resources Research, 27(5), 667-698.
Gómez-Hernández J.J, Sovero H. and Sahuquillo A. 1995. Some issues on the analysis of pumping test in heterogeneous aquifers. Hydrogeologie, 3, 13-18.
Gorelick S.M. 1983. A review of distributed parameter groundwater management modeling methods. Water Resources Research, 19(2), 305-319. https://doi.org/10.1029/WR019i002p00305
Herrera I. and Figueroa G.E. 1969. A correspondence principe for the theory of leaky aquifers. Water Resources Research, 5(4), 900-904. https://doi.org/10.1029/WR005i004p00900
Herrera I. 1970. Theory of multiple leaky aquifers. Water Resources Research, 6(1), 185-193. https://doi.org/10.1029/WR006i001p00185
Lachassagne P., Ledoux E. and de Marsily G. 1989. Evaluation of hydrogeological parameters in heterogeneous porous media. Groundwater Management, Quantity and Quality. IAHS Publ. nº 188. 3-18.
Matheron G. 1967. Composition des perméabilités en milieu poreux hétérogène e de Schwydler et règles de pondération. Révue de l'Institute Français du Pètrole, 3, XXII, 443-446.
Morohunfola O.K. 1992. Analytical and numerical solutions for multiple leaky aquifer systems. Ph. D. dissertation. Univ. of Delaware, Newark.
Peng-Hsiang T. and Tien-Chang L. 1998. Numerical evaluation of exponential integral: Theis well function approximation. Journal of Hydrology, 205(1-2), 38-51. https://doi.org/10.1016/S0022-1694(97)00134-0
Piessens R. 1975. A bibliography on numerical inversion of the Lapace transform and applications. Journal of Computational and Applied Mathematics, 1, 115-126. https://doi.org/10.1016/0771-050X(75)90029-7
Piessens R. and Dang N.D.P. 1976. A bibliography on numerical inversion of the Laplace transform and applications, A supplement, Journal of Computational and Applied Mathematics, 2, 225-228. https://doi.org/10.1016/0771-050X(76)90009-7
Prickett T.A. 1967. Designing pumped well characteristics into electric analog model. Ground Water, 5(4), 38-46. https://doi.org/10.1111/j.1745-6584.1967.tb01625.x
Sahuquillo A. 1983. An Eigenvalue Numerical Technique for Solving Unsteady Linear Groundwater Models Continuously in Time. Water Resources Research, 19(1), 87-93. https://doi.org/10.1029/WR019i001p00087
Sánchez-Vila X., Carrera J., Jorge P. and Girardi J.P. 1996. Scale effects in transmissivity. Journal of Hydrology, 183(1-2), 1-22. https://doi.org/10.1016/S0022-1694(96)80031-X
Stehfest H. 1970. Numerical inversion of Laplace transforms. Communications of the ACM., 13 (1), 47-49. https://doi.org/10.1145/361953.361969
Theis C.V. 1935. The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage. Transactions American Geophysical Union Annual Meeting, 16th, 519-524. https://doi.org/10.1029/TR016i002p00519
Thiery D., Schwartz J., Berge J., Fotoohi F., Konstantopedos K. and Lambert M. 1995. Un système d'aide à la gestion des ressources en eaux souterraines. Aplication au site de Bordeaux. Hydrogeologie, 1, 129-139.
Wang H.F. and Anderson M.P. 1982. Introduction to Groundwater Modeling: Finite Difference and Finite Element Methods. W.H. Freeman and Co., N. York.
Wolfram S. 1991. Mathematica. A System for doing Mathematics by Computer. Addison-Wesley Pub.Co. 2nd ed. N York.
Zheng, C. 1997. ModGA, Documentation and User's Guide, Technical Report to DuPont Company, Hydrogeology Group, Univ. of Alabama.
Downloads
Published
How to Cite
Issue
Section
License
Copyright (c) 2026 Consejo Superior de Investigaciones Científicas (CSIC)

This work is licensed under a Creative Commons Attribution 4.0 International License.
© CSIC. Manuscripts published in both the print and online versions of this journal are the property of the Consejo Superior de Investigaciones Científicas, and quoting this source is a requirement for any partial or full reproduction.
All contents of this electronic edition, except where otherwise noted, are distributed under a Creative Commons Attribution 4.0 International (CC BY 4.0) licence. You may read the basic information and the legal text of the licence. The indication of the CC BY 4.0 licence must be expressly stated in this way when necessary.
Self-archiving in repositories, personal webpages or similar, of any version other than the final version of the work produced by the publisher, is not allowed.






