Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

中国大陆及邻区上地幔P波各向异性结构

利用分布在中国大陆及邻近地区的213个地震台站记录到的远震P波走时数据和弱各向异性条件下P波速度扰动调和分析方法,研究了中国大陆上地幔P波各向异性结构.研究结果表明中国大陆西部上地幔变形主要受印度大陆俯冲的影响.印度大陆的P波快波方向总体为NNE方向,与绝对板块运动方向一致,这表明印度大陆上地幔流动方向与板块运动方向一致.青藏高原内部、东天山的P波快波方向与主压应力方向接近,而在青藏高原南缘、北缘及东北缘等块体边界地区P波快波方向与主压应力方向垂直.中国大陆东部上地幔变形主要受菲律宾板块和太平洋板块俯冲的影响.在扬子板块内部P波快波方向为SE方向,这与绝对板块的运动方向一致.华北地区的各向异性结构较为复杂,可能与华北克拉通裂解有关.中国大陆东北的东部平均方向为SE,而在兴安岭一侧为SSW方向,即平行于构造线方向.根据各向异性的倾角,中国大陆及邻区上地幔各向异性结构大体可分为三块:1)青藏新疆地区的各向异性倾角接近水平,推测该区形变力源主要为上地幔物质水平流动.2)南北带地区的各向异性倾角较大,特别是在青藏东缘地区的倾角约为40°,这可能是由于青藏向东挤出过程中受华南地块和鄂尔多斯地块的阻挡,在板块边界地区产生了垂直变形.3)中国东部地区各向异性结构较为复杂,在中国大陆东北部各向异性倾角接近水平,这可能是该区上地幔变形主要受太平洋板块俯冲的影响,而在太行山、大别-苏鲁地区各向异性倾角较大,这表明该区上地幔以垂直变形为主.     

High-permeability layers for remediation of ground water; go wide, not deep

A nitrate-reactive porous media layer comprising wood particles with very high hydraulic conductivity (K approximately 1 cm/s) was used to successfully treat nitrate in a shallow sand-and-gravel aquifer in southern Ontario. Nitrate concentrations of 1.3 to 14 mg/L as N in the aquifer were attenuated to <0.5 mg/L as N in the reactive layer. Borehole dilution testing indicated that ground water velocities in the reactive layer, although variable, averaged five times higher than in the surrounding aquifer, suggesting that the layer was capturing ground water flow from deeper in the aquifer. The use of high-K reactive media opens up the possibility of installing permeable reactive barriers as horizontal layers in the shallow water table zone that do not necessarily have to penetrate the full depth of a contaminant plume to be effective. Model simulations show that the depth of capture of a high-K layer increases as the layer width in the direction of flow increases. Shallower emplacement could decrease barrier costs at some sites.     

Estimation of anisotropic hydraulic conductivity using geophysical data in a coastal aquifer of Karnataka,India

Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models.     

Anisotropic seismic ground responses identified from the Hualien vertical array

This paper presents the results of ground response analyses using the free-field motions of six earthquakes recorded at the site of the Hualien Large Scale Seismic Test. From the phase spectra of transfer functions between the ground motions recorded at the surface and various depths, the predominant frequencies of soil layers are identified and used to estimate the shear wave velocities of soil layers. Results show that the soils at the Hualien LSST site are anisotropic, such that the earthquake responses exhibit different characteristics in two horizontal principal directions.     

Three-dimensional deformation and strain induced by municipal pumping, Part 2: Numerical analysis

Field measurements consisting of water levels from a municipal well and three-dimensional surface deformations and strains from high-precision GPS measurements at various radial distances from the well were collected as part of a 62-day controlled aquifer test at Mesquite, NV. These measurements were used as observations in several numerical models and a parameter estimation code to characterize and constrain hydraulic and mechanical properties of a 400 m thick basin-fill aquifer. A parsimonious approach was used in conceptualizing the aquifer system. Nonetheless, results from the calibrated deformation and flow models accurately reproduced the observed head and deformations during the first 20 days of pumping, the time at which a new equilibrium was achieved. Surface deformations were shown to reflect hydraulic anisotropy and direction of principal conductivity. In addition, the radius of influence and cone of depression from pumping was approximated in spite of the fact that no monitoring well data existed at the site. Sensitivity analysis indicates that cyclical head values are most sensitive to changes in horizontal hydraulic conductivity, while time-dependent vertical deformations are most sensitive to changes in skeletal specific storage. This investigation shows that GPS monitoring can be used in place of costly monitoring wells to characterize aquifers for water-management purposes where skeletal deformation tends to be elastic.     

DARCY'S COEFFICIENT OF PERMEABILITY AS SYMMETRIC TENSOR OF SECOND RANK

Synopsis

The dynamic equation of motion that governs the laminar flow of water through soils is the empirical equation of Darcy. According to Darcy's equation the velocity of the flowing water is proportional to the hydraulic gradient under which the water is flowing, with the constant of proportionality being the coefficient of permeability. The interesting question arising is whether or not the coefficient of permeability is a scalar quantity (having only a magnitude) or a vector (having both magnitude and direction). It is proved, in the present paper, that the permeability coefficient is neither a scalar nor a vector but a symmetric tensor of second rank. The fact that the permeability tensor is symmetric gives rise to great simplifications and permits a simple graphical construction of the tensor ellipsoid. Having the tensor ellipsoid, the determination of the direction at which the water will flow under a known imposed hydraulic gradient can be found graphically. In case of isotropic soils (the permeability coefficient has the same value along any direction) the ellipsoid reduces to a sphere and the tensor becomes a scalar. In the general case of anisotropic soils the permeability tensor is an entity with nine elements, six of which are independent representing pure extension or contraction along the three principal coordinate axes, thus transforming the permeability sphere into an ellipsoid and vice versa. It should be noted that in anisotropic soils the only directions along which the flow takes place in the direction of the hydraulic gradient are those of the principal axes of the tensor ellipsoid.

Permeability tests were conducted on anisotropic sandstone samples taken at different directions with respect to rectangular coordinates. The permeability coefficient values plotted on a two-dimensional polar coordinate graph paper give rise to an ellipse substantiating therefore the tensor concept of the permeability coefficient. The graphical construction of the tensor ellipse and the use of it in order to obtain the direction of flow by knowing the direction of the hydraulic gradient is also shown.     

Analysis of Transient Flow to an Extended Fully Penetrating Well at Constant-Head Pumping

Vertical wells with radial extension at the well bottom can improve the rate of water production. No study has yet investigated the effects of the transient state and anisotropy in directional hydraulic conductivities on the wellbore flux rate for this type of well. This study derives a semianalytical transient drawdown solution for constant-head pumping at a fully penetrating well radially extended at the bottom of a confined, anisotropic aquifer by applying Laplace transform and separation of variables as well as conducting a Fourier analysis. The results of this new solution indicate that transient and steady-state wellbore flux rates can be increased by a factor of two for greater radial extension of the well. Compared with an isotropic aquifer (a ratio of vertical and horizontal hydraulic conductivities equal to one), an anisotropic aquifer with the ratio less than one may produce a higher transient wellbore flux rate and lower steady-state wellbore flux rate. Moreover, the time required to achieve the steady-state wellbore flux rate can be substantially affected by anisotropy of the aquifer.     

Three-dimensional flow in the storative semiconfining layers of a leaky aquifer

An analytical solution for three-dimensional (3D) flow in the storative semiconfining layers of a leaky aquifer fully penetrated by a production well is developed in this article to provide a method from which accurate hydraulic parameters in the semiconfining layers can be derived from aquifer test data. The analysis of synthetic aquifer test data with the 3D analytical solution in the semiconfining layers provided more accurate optimal hydraulic parameters than those derived using the available quasi-two-dimensional (2D) solution. Differences between the 3D and 2D flow solutions in the semiconfining layers become larger when a no flow boundary condition is imposed at either at the top of the upper semiconfining layer or at the bottom of the lower semiconfining layer or when the hydraulic conductivity ratio of the semiconfining layer to the aquifer is larger than 0.001. In addition, differences between the 3D and 2D flow solutions in the semiconfining layers are illustrated when the thickness ratio of the semiconfining layer to the aquifer is changed. Analysis of water level data from two hypothetical and one real aquifer test showed that the 3D solution in the semiconfining layers provides lower correlation coefficients among hydraulic parameters than the 2D solution.     

A Finite Layer Formulation for Groundwater Flow to Horizontal Wells

A finite layer approach for the general problem of three‐dimensional (3D) flow to horizontal wells in multilayered aquifer systems is presented, in which the unconfined flow can be taken into account. The flow is approximated by an integration of the standard finite element method in vertical direction and the analytical techniques in the other spatial directions. Because only the vertical discretization is involved, the horizontal wells can be completely contained in one specific nodal plane without discretization. Moreover, due to the analytical eigenfunctions introduced in the formulation, the weighted residual equations can be decoupled, and the formulas for the global matrices and flow vector corresponding to horizontal wells can be obtained explicitly. Consequently, the bandwidth of the global matrices and computational cost rising from 3D analysis can be significantly reduced. Two comparisons to the existing solutions are made to verify the validity of the formulation, including transient flow to horizontal wells in confined and unconfined aquifers. Furthermore, an additional numerical application to horizontal wells in three‐layered systems is presented to demonstrate the applicability of the present method in modeling flow in more complex aquifer systems.     

Aquifer response to stream-stage and recharge variations. II. Convolution method and applications

In this second of two papers, analytical step-response functions, developed in the companion paper for several cases of transient hydraulic interaction between a fully penetrating stream and a confined, leaky, or water-table aquifer, are used in the convolution integral to calculate aquifer heads, streambank seepage rates, and bank storage that occur in response to stream-stage fluctuations and basinwide recharge or evapotranspiration. Two computer programs developed on the basis of these step-response functions and the convolution integral are applied to the analysis of hydraulic interaction of two alluvial stream–aquifer systems in the northeastern and central United States. These applications demonstrate the utility of the analytical functions and computer programs for estimating aquifer and streambank hydraulic properties, recharge rates, streambank seepage rates, and bank storage. Analysis of the water-table aquifer adjacent to the Blackstone River in Massachusetts suggests that the very shallow depth of water table and associated thin unsaturated zone at the site cause the aquifer to behave like a confined aquifer (negligible specific yield). This finding is consistent with previous studies that have shown that the effective specific yield of an unconfined aquifer approaches zero when the capillary fringe, where sediment pores are saturated by tension, extends to land surface. Under this condition, the aquifer's response is determined by elastic storage only. Estimates of horizontal and vertical hydraulic conductivity, specific yield, specific storage, and recharge for a water-table aquifer adjacent to the Cedar River in eastern Iowa, determined by the use of analytical methods, are in close agreement with those estimated by use of a more complex, multilayer numerical model of the aquifer. Streambank leakance of the semipervious streambank materials also was estimated for the site. The streambank-leakance parameter may be considered to be a general (or lumped) parameter that accounts not only for the resistance of flow at the river–aquifer boundary, but also for the effects of partial penetration of the river and other near-stream flow phenomena not included in the theoretical development of the step-response functions.     

Determining anisotropic transmissivity using a simplified Papadopulos method

The straight-line method presented by Papadopulos requires a minimum of three observation wells for determining the transmissivity tensor of a homogeneous and anisotropic aquifer. A simplification of this method was developed for fractured aquifers where the principal directions of the transmissivity tensor are known prior to implementation, such as when fracture patterns on outcropping portions of the aquifer may be used to infer the principal directions. This new method assumes that observation wells are drilled along the two principal directions from the pumped well, thus reducing the required number of observation wells to two. This method was applied for an aquifer test in the fractured Navajo Sandstone of southwestern Utah and yielded minimum and maximum principal transmissivity values of 70 and 1800 m(2)/d, respectively, indicating an anisotropy ratio of approximately 24 to 1.     

Measurement of 3-D hydraulic conductivity in aquifer cores at in situ effective stresses

An innovative and nondestructive method to measure the hydraulic conductivity of drill core samples in horizontal and vertical directions within a triaxial cell has been developed. This has been applied to characterizing anisotropy and heterogeneity of a confined consolidated limestone aquifer. Most of the cores tested were isotropic, but hydraulic conductivity varied considerably and the core samples with lowest values were also the most anisotropic. Hydraulic conductivity decreased with increasing effective stress due to closure of microfractures caused by sampling for all core samples. This demonstrates the importance of replicating in situ effective stresses when measuring hydraulic conductivity of cores of deep aquifers in the laboratory.     

Performance of the steady-state dipole flow test in layered aquifers

Steady-state numerical simulations of the dipole flow test in layered aquifers demonstrate that the test produces a good estimate of the equivalent hydraulic conductivity anisotropy ratio for the part of the aquifer spanned by the well chambers. The effects of chamber size, different conductivity of layers and layer location on the estimated anisotropy ratios are presented. The steady-state dipole flow test, when performed at different levels in the well, can yield estimates of the down-hole anisotropy ratio distribution. Numerical simulations also illustrate that the skin effect can significantly distort the anisotropy estimates produced by the dipole flow test. © 1997 John Wiley & Sons, Ltd.     

A Method to Estimate the Hydraulic Conductivity of the Ground by TRT Analysis

The knowledge of hydraulic properties of aquifers is important in many engineering applications. Careful design of ground‐coupled heat exchangers requires that the hydraulic characteristics and thermal properties of the aquifer must be well understood. Knowledge of groundwater flow rate and aquifer thermal properties is the basis for proper design of such plants. Different methods have been developed in order to estimate hydraulic conductivity by evaluating the transport of various tracers (chemical, heat etc.); thermal response testing (TRT) is a specific type of heat tracer that allows including the hydraulic properties in an effective thermal conductivity value. Starting from these considerations, an expeditious, graphical method was proposed to estimate the hydraulic conductivity of the aquifer, using TRT data and plausible assumption. Suggested method, which is not yet verified or proven to be reliable, should be encouraging further studies and development in this direction.     

Effects of Improper Characterization of Aquifer Thickness on Estimates of Hydraulic Conductivity from Pumping Tests

Pumping test data for surficial aquifers are commonly analyzed under the assumption that the base of the aquifer corresponds to the bottom of the test wells (i.e., the aquifer is truncated). This practice can lead to inaccurate hydraulic conductivity estimates, resulting from the use of low saturated thickness values with transmissivity estimates, and not accounting for the effects of partially penetrating wells. Theoretical time-drawdown data were generated at an observation well in a hypothetical unconfined aquifer for various values of saturated thickness and were analyzed by standard curve-matching techniques. The base of the aquifer was assumed to be the bottom of the pumping and observation wells. The overestimation of horizontal hydraulic conductivity was found to be directly proportional to the error in assumed saturated thickness, and to the (actual) ratio of vertical to horizontal hydraulic conductivity (K_v/K_h). Inaccurately high estimates of hydraulic conductivity obtained by aquifer truncation can lead to overestimates of ground water velocity and contaminant plume spreading, narrow capture zone configuration estimates, and overestimates of available ground water resources.     

General Steady-State Shape Factor for a Partially Penetrating Well

We present the closed form of a general steady-state shape factor for a partially penetrating well in a uniform anisotropic aquifer. Our simple analytical expression for the shape factor has a uniform representation for full range of parameters and meets or exceeds the accuracy of previous results obtained through semiempirical methods (e.g., Bouwer and Rice [1976] equations). This general shape factor pertains to the flow of fluids (water or air) in subsurface formations when the upper formation boundary has constant potential and the lower boundary is impermeable. The results of our investigation are directly applicable to analyses of (1) slug tests with falling or rising head and (2) injection/extraction tests with constant head, essential techniques for the characterization of hydraulic conductivity of aquifers, streambeds, or lakebeds as well as the design of aquifer and soil remediation systems.     

Estimating hydraulic parameters when poroelastic effects are significant

For almost 80 years, deformation-induced head changes caused by poroelastic effects have been observed during pumping tests in multilayered aquifer-aquitard systems. As water in the aquifer is released from compressive storage during pumping, the aquifer is deformed both in the horizontal and vertical directions. This deformation in the pumped aquifer causes deformation in the adjacent layers, resulting in changes in pore pressure that may produce drawdown curves that differ significantly from those predicted by traditional groundwater theory. Although these deformation-induced head changes have been analyzed in several studies by poroelasticity theory, there are at present no practical guidelines for the interpretation of pumping test data influenced by these effects. To investigate the impact that poroelastic effects during pumping tests have on the estimation of hydraulic parameters, we generate synthetic data for three different aquifer-aquitard settings using a poroelasticity model, and then analyze the synthetic data using type curves and parameter estimation techniques, both of which are based on traditional groundwater theory and do not account for poroelastic effects. Results show that even when poroelastic effects result in significant deformation-induced head changes, it is possible to obtain reasonable estimates of hydraulic parameters using methods based on traditional groundwater theory, as long as pumping is sufficiently long so that deformation-induced effects have largely dissipated.