Many so-called sandstone aquifers are actually multiple-aquifer systems consisting of discontinuous sand bodies distributed complexly in a matrix of lower-permeability silts and clays. The arrangement and Interconnectedness of these various lithofacies strongly influence spatial patterns of hydraulic conductivity (K) and, in turn, groundwater flow and mass transport. A promising technique of estimating such patterns of K involves careful analysis of both subsurface geologic and subsurface hydrologic data. In this study the three-dimensional distribution of K was estimated for a numerical flow model of part of the Wilcox aquifer system in Texas, using K data from core samples and pumping tests and more than 100 geophysical logs. The aquifer system, which is up to 320 m thick, consists of multiple, elongate sand bodies and silts and clays deposited in a fluvial environment and is similar to many other systems found in the Gulf Coast and other sedimentary basins. The resulting deterministic-conceptual flow model demonstrates the importance and methods of incorporating geologic information in groundwater models. Flow in the aquifer is shown to be controlled not so much by K of the sands as by their continuity and Interconnectedness. Much of the aquifer system consists of large zones in which the fluvial channel-fill sands are sparse and apparently disconnected, resulting in groundwater flow rates lower by a factor of 101 to 103 than in adjacent, well-interconnected belts of fluvial channel-fill sand belts. Modeling results also raise serious doubts regarding our ability to predict regional scale flow and mass transport in complex aquifers such as the Wilcox, using current technology. Though sand body Interconnectedness is critically important, it is also very difficult to estimate. One or two well-connected sands among a system of otherwise disconnected sands can completely alter a velocity field. This is particularly true if the sands are connected vertically and nonzero vertical hydraulic gradients exist. Because the model is three-dimensional, sensitivity of hydraulic head to heterogeneity or Interconnectedness is much less than normally observed in two-dimensional models, and therefore heads computed by the model give little to no indication of the location of well-interconnected zones. Thus such zones can easily go undetected, even in carefully calibrated models which yield reasonably accurate hydraulic heads. This is a significant point for modeling of solute transport.
[1]
William F. Guyton.
GROUND-WATER CONDITIONS IN ANDERSON, CHEROKEE, FREESTONE, AND HENDERSON COUNTIES, TEXAS
,
1972
.
[2]
W. E. Galloway,et al.
Catahoula Formation of the Texas Coastal Plain: Depositional Systems, Composition, Structural Development, Ground-Water Flow History, and Uranium Distribution
,
1977
.
[3]
N. T.,et al.
NUMERICAL MODEL FOR SATURATED-UNSATURATED FLOW IN * DEFORMABLE POROUS MEDIA, PART I: THEORY
,
2007
.
[4]
H. Reineck,et al.
Depositional sedimentary environments
,
1973
.
[5]
G. Dagan.
Stochastic modeling of groundwater flow by unconditional and conditional probabilities: 1. Conditional simulation and the direct problem
,
1982
.
[6]
Franklin W. Schwartz,et al.
Mass transport: 1. A stochastic analysis of macroscopic dispersion
,
1980
.
[7]
Allan L. Gutjahr,et al.
Stochastic analysis of spatial variability in subsurface flows: 1. Comparison of one‐ and three‐dimensional flows
,
1978
.
[8]
W. Pryor.
Permeability-Porosity Patterns and Variations in Some Holocene Sand Bodies
,
1973
.
[9]
R. Freeze.
A stochastic‐conceptual analysis of one‐dimensional groundwater flow in nonuniform homogeneous media
,
1975
.
[10]
O. Weres,et al.
DOCUMENTATION FOR PROGRAM OGRE
,
1978
.
[11]
Franklin W. Schwartz,et al.
Mass transport: 3. Role of hydraulic conductivity data in prediction
,
1981
.
[12]
W. R. Kaiser,et al.
Sand-Body Geometry and the Occurrence of Lignite in the Eocene of Texas
,
1978
.
[13]
R. E. Eddy,et al.
Predicting Reservoir Rock Geometry and Continuity in Pennsylvanian Reservoirs, Elk City Field, Oklahoma
,
1977
.
[14]
H. S. Price,et al.
Flow in Heterogeneous Porous Media
,
1961
.
[15]
Franklin W. Schwartz,et al.
Macroscopic dispersion in porous media: The controlling factors
,
1977
.
[16]
A. J. Scott,et al.
Reservoir characterization for numerical simulation of Mesaverde meanderbelt sandstone, northwestern Colorado
,
1984
.
[17]
Paul A. Witherspoon,et al.
Numerical model for saturated‐unsaturated flow in deformable porous media: 3. Applications
,
1978
.
[18]
D. Bottjer,et al.
Sedimentary Environments and Facies
,
1978
.
[19]
C. Axness,et al.
Three‐dimensional stochastic analysis of macrodispersion in aquifers
,
1983
.