Documentation of a computer program to simulate transient leakage from confining units using the modular finite-difference, ground-water flow model

This report presents a new method of simulating transient leakage from confining units using the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model. Transient leakage into or out of a compressible fine-grained confining unit results from ground-water storage changes within the unit. The importance of fine-grained units in analyses of transient ground-water flow is not always recognized. The new method of simulating transient leakage in the Modular Finite-Difference Ground-Water Flow Model is referred to as the Transient-Leakage Package, version 1, or the TLK1 Package. The TLK1 Package solves integrodifferential equations that describe the flow components across the upper and lower boundaries of confining units. The exact equations are approximated to allow efficient solution for the flow components. The flow components are incorporated into the finite-difference equations for model cells that are adjacent to confining units. Vertical hydraulic conductivity, thickness, and specific storage are specified in input arrays for each confining unit. Confining-unit properties can differ from cell to cell and a confining unit need not be present at all locations; however, a confining unit must be bounded above and below by model layers in which head is calculated or specified. The TLK1 Package was used for an example problem to simulate drawdown around a pumped well in a system with two aquifers separated by a confining unit. The pumping was limited to one of the two aquifers. The exact drawdown values determined from an analytical solution were compared with the simulated drawdown. The solution using the TLK1 Package closely matched the exact solution for drawdown values in excess of 1 centimeter in the pumped and unpumped aquifers. The problem also was simulated without the TLK1 Package by using a separate model layer to represent the confining unit. That simulation was further refined by using two model layers to represent the confining unit. The simulation made using the TLK1 Package was faster and more accurate than either of the simulations using model layers to represent the confining unit.

[1]  J. C. Jaeger,et al.  Conduction of Heat in Solids , 1952 .

[2]  S. A. Leake.  Simulation of vertical compaction in models of regional ground-water flow , 1991 .

[3]  R. L. Cooley.  Numerical simulation of flow in an aquifer overlain by a water table aquitard , 1972 .

[4]  S. P. Neuman,et al.  Theory of Flow in a Confined Two Aquifer System , 1969 .

[5]  Arlen W. Harbaugh,et al.  A modular three-dimensional finite-difference ground-water flow model , 1984 .

[6]  Richard L. Cooley,et al.  Derivation of finite-element equations and comparisons with analytical solutions , 1992 .

[7]  C. E. Jacob.  On the flow of water in an elastic artesian aquifer , 1940 .

[8]  Ismael Herrera,et al.  Integrodifferential equations for systems of leaky aquifers and applications: 1. The nature of approximate theories , 1973 .

[9]  David E. Prudic,et al.  Documentation of a computer program to simulate aquifer-system compaction using the modular finite d , 1988 .

[10]  Extension of the integrodifferential approach to inhomogeneous multiaquifer systems , 1981 .

[11]  A Conjugate Gradient Finite Element Model of Flow for Large Multiaquifer Systems , 1986 .

[12]  S. Leake.  Interbed storage changes and compaction in models of regional groundwater flow , 1990 .

[13]  J. D. Bredehoeft,et al.  Digital Analysis of Areal Flow in Multiaquifer Groundwater Systems: A Quasi Three‐Dimensional Model , 1970 .

[14]  M. S. Hantush.  Modification of the theory of leaky aquifers , 1960 .

[15]  Emmanuel Ledoux,et al.  Modelling of large multilayered aquifer systems: Theory and applications , 1978 .

[16]  I. Herrera,et al.  A Correspondence Principle for the Theory of Leaky Aquifers , 1969 .

[17]  I. Herrera.  Theory of Multiple Leaky Aquifers , 1970 .

[18]  J. Premchitt.  A technique in using integrodifferential equations for model simulation of multiaquifer systems , 1981 .

[19]  R. Yates,et al.  Integrodifferential equations for systems of leaky aquifers and applications 3. A numerical method of unlimited applicability , 1977 .

[20]  Mary C. Hill,et al.  PRECONDITIONED CONJUGATE-GRADIENT 2 (PCG2), a computer program for solving ground-water flow equations , 1990 .

[21]  George F. Pinder,et al.  Finite difference model for aquifer simulation in two dimensions with results of numerical experiments , 1976 .

[22]  J. D. Bredehoeft,et al.  On the Maintenance of Anomalous Fluid Pressures: II. Source Layer at Depth , 1968 .