承压含水层井流-盖层弯曲效应的解析理论

指出承压含水层盖层的弯曲变形与开采井周围的径向地下水运动存在相互作用, 而这一效应在传统的井流理论中没有被认识到.通过引入弹性薄板理论, 建立了无越流的承压含水层井流-顶板弯曲效应的解析模型, 同时考虑了含水层和水的压缩性, 结果表明Theis井流方程给出的抽水降深偏小.在此基础上推导了有越流承压含水层井流-盖层弯曲效应的偏微分方程, 求出了解析解, 并与传统理论的结果进行了对比, 表明Hantush-Jacob公式计算的降深也是偏小的.在抽水井附近和抽水初期, 传统理论可能导致显著的相对误差.      

水资源和环境工程中水平井研究简介

简单介绍了水平井在水资源和环境工程中的近期研究进展, 重点介绍了在不同含水层(如承压含水层、潜水含水层、越流含水层及河下含水层) 中, 水平井抽水条件下降深的半解析解, 同时分析了降深随抽水时间变化的标准曲线和微分标准曲线.这些半解析解可用于分析小流量水平井在中长时段的降深特性.分析了用于排水和供水的大流量水平井的水力学特征, 并介绍了求解渗流-管流耦合井流系统的一种新方法.同时介绍了在非饱和含水层中水平井抽取气体的动力学特征, 并分析了地面覆盖和未覆盖2种情况下的气体清除效率, 报告了在不同饱和含水层条件下水平井的捕获区和捕获时间的计算.      

Hydrogeologic and hydraulic characterization of aquifer and nonaquifer layers in a lateritic terrain (West Bengal,India)

The hydrogeologic and hydraulic characteristics of a lateritic terrain in West Bengal, India, were investigated. Test drilling was conducted at ten sites and grain-size distribution curves (GSDCs) were prepared for 275 geologic samples. Performance evaluation of eight grain-size-analysis (GSA) methods was carried out to estimate the hydraulic conductivity (K) of subsurface formations. Finally, the GSA results were validated against pumping-test data. The GSDCs indicated that shallow aquifer layers are coarser than the deeper aquifer layers (uniformity coefficient 0.19–11.4). Stratigraphy analysis revealed that both shallow and deep aquifers of varying thickness exist at depths 9–40 and 40–79 m, respectively. The mean K estimates by the GSA methods are 3.62–292.86 m/day for shallow aquifer layers and 0.97–209.93 m/day for the deeper aquifer layers, suggesting significant aquifer heterogeneity. Pumping-test data indicated that the deeper aquifers are leaky confined with transmissivity 122.69–693.79 m²/day, storage coefficient 1.01?×?10^?7–2.13?×?10^?4 and leakance 2.01?×?10^?7–34.56?×?10^?2 day^?1. Although the K values yielded by the GSA methods are generally larger than those obtained from the pumping tests, the Slichter, Harleman and US Bureau Reclamation (USBR) GSA methods yielded reasonable values at most of the sites (1–3 times higher than K estimates by the pumping-test method). In conclusion, more reliable aquifers exist at deeper depths that can be tapped for dependable water supply. GSA methods such as Slichter, Harleman and USBR can be used for the preliminary assessment of K in lateritic terrains in the absence of reliable field methods.     

Automatic estimation of aquifer parameters using long-term water supply pumping and injection records

Analyses are presented of long-term hydrographs perturbed by variable pumping/injection events in a confined aquifer at a municipal water-supply well field in the Region of Waterloo, Ontario (Canada). Such records are typically not considered for aquifer test analysis. Here, the water-level variations are fingerprinted to pumping/injection rate changes using the Theis model implemented in the WELLS code coupled with PEST. Analyses of these records yield a set of transmissivity (T) and storativity (S) estimates between each monitoring and production borehole. These individual estimates are found to poorly predict water-level variations at nearby monitoring boreholes not used in the calibration effort. On the other hand, the geometric means of the individual T and S estimates are similar to those obtained from previous pumping tests conducted at the same site and adequately predict water-level variations in other boreholes. The analyses reveal that long-term municipal water-level records are amenable to analyses using a simple analytical solution to estimate aquifer parameters. However, uniform parameters estimated with analytical solutions should be considered as first rough estimates. More accurate hydraulic parameters should be obtained by calibrating a three-dimensional numerical model that rigorously captures the complexities of the site with these data.     

Analytical solutions for constant-rate test in bounded confined aquifers with non-Darcian effect

This paper proposes a simplified analytical solution considering non-Darcian and wellbore storage effect to investigate the pumping flow in a confined aquifer with barrier and recharge boundaries. The mathematical modelling for the pumping-induced flow in aquifers with different boundaries is developed by employing image-well theory with the superposition principle, of which the non-Darcian effect is characterized by Izbash’s equation. The solutions are derived by Boltzmann and dimensionless transformations. Then, the non-Darcian effect and wellbore storage are especially investigated according to the proposed solution. The results show that the aquifer boundaries have non-negligible effects on pumping, and ignoring the wellbore storage can lead to an over-estimation of the drawdown in the first 10 minutes of pumping. The higher the degree of non-Darcian, the smaller the drawdown.     

Multi-aquifer pumping test to determine cutoff wall length for groundwater flow control during site excavation in Tokyo, Japan

Multi-aquifer pumping tests, using a multi-screen pumping well and multi-level piezometers, were carried out for groundwater flow control in a large-scale excavation site in Tokyo, Japan. The site was underlain by multi-layered confined aquifers. In the tests, pumping was carried out using a multi-aquifer pumping well in which a screen depth was chosen arbitrarily. Changes in groundwater pressure heads in each aquifer were measured at each screen position of the multi-aquifer pumping well. Hydraulic conductivity (K) and specific storage (S _s) of not only aquifers, but also for low permeability layers between the aquifers, were estimated using the Cooper-Jacob method, and calibrated by a finite element method (FEM) groundwater model. Four different cutoff wall lengths were assumed for final excavation depth, and correlations among wall length, pumping discharge and drawdown at the back of the cutoff wall were obtained from simulations using the K and S _s parameters in the FEM model. Then, the most suitable wall length was selected based on the simulated correlations considering environmental condition, construction period and cost of the cutoff wall.     

A generalized solution for groundwater head fluctuation in a tidal leaky aquifer system

A new analytical solution is developed for describing groundwater level fluctuations in a coupled leaky confined aquifer system which consists of an unconfined aquifer, confined aquifer, and an aquitard in between. The aquifer system has a tidal boundary at the seashore, a no flow boundary at remote inland side, and a confined aquifer extending under the sea and terminated with an outlet-capping. This new solution has shown to be a generalisation of most existing analytical solutions for a tidal aquifer system which includes single confined and leaky confined aquifers. In addition, the solution is used to explore the influences of the dimensionless leakance of the outlet-capping, the dimensionless hydraulic diffusivities, and the leakages of the inland and offshore aquitards on the head responses in the leaky confined aquifer.     

Geochemistry of heavily exploited aquifers in the Emilia-Romagna region (Po Valley, northern Italy)

 More than 5 800 chemical analyses on water samples collected during 1987–1995 from 528 monitoring wells located in the southernmost part of the Po Valley (Emilia-Romagna region, northern Italy), one of the most urbanized, industrialized and agriculturally developed areas of Italy, have been processed. The analysis of data showed that: (1) waters are discharging from both confined and unconfined aquifers; (2) the water in the unconfined aquifer(s) is Ca(Mg)-HCO₃ in composition while confined ones are Na-Cl and/or Na-(HCO₃); (3) both confined and unconfined aquifer samples have δ¹⁸O and δD isotopic values of meteoric signature; (4) waters from both the aquifers are at least 40 years old; (5) the pumping rate has caused subsidence, particularly where the aquifer(s) is (are) unconfined; (6) the unconfined aquifer(s) is exposed to the risk of NO₃ pollution; (7) considering the present "pressure" (i.e. pumping rate) on this natural environment by human activity, care must be taken in the future to preserve this "strategic" resource. Received: 27 October 1997 · Accepted: 12 March 1998     

Identification of aquifer and well parameters from step-drawdown tests

A new analysis technique has been proposed for interpreting transient step-drawdown test data. The proposed method is based on taking the derivative of the drawdown with respect to time for the entire pumping test period to eliminate the time-independent well-loss terms. The derivative function is subsequently integrated to obtain the time-dependent aquifer drawdown as a continuous function. The well-loss parameters are then obtained with higher accuracy once the aquifer behavior is identified. The proposed method is applicable for analyzing data obtained not only from ideal confined aquifers but also from other aquifer types (i.e. unconfined) and non-ideal aquifers (i.e. heterogeneous). The technique was tested for synthetically generated and field data; the proposed approach was noted to provide accurate aquifer and well-loss parameter estimates. The results of the proposed method were compared with those of some of the existing methods for analyzing step-drawdown test data and were found to be more reliable and robust.     

Analysis of groundwater-level fluctuation in a coastal confined aquifer induced by sea-level variation

Groundwater responses to tide fluctuations in different hydrogeological situations have been investigated for many years. Various solutions have been derived using the assumption that tides are composed of sinusoidal components, neglecting non-periodic variables. This approach has resulted in the introduction of some inaccuracy in predictions of groundwater responses to sea-level variation. To resolve this problem, this study used the Fourier sine transform method to derive an analytical solution based on the measured sea-level boundary. Compared with an analytical solution based on a sinusoidal assumption and a numerical solution generated by MODFLOW, this solution provided better performance in groundwater-level prediction in a coastal confined aquifer in Zhuhai City, China. The hydrogeological parameters estimated by the three aforementioned methods fitted well with those estimated by field surveys and pumping tests. The introduced analytical solution not only reflects the physical mechanisms of tide-induced groundwater-level fluctuation, but also reveals the non-periodic fluctuation of groundwater level caused by sea-level variation. This solution may also be used to evaluate groundwater responses to random mechanical stresses.     

利用不对称的地下水位潮汐波动确定滨海含水层参数

在滨海地区,地下水水位受潮汐波动影响较大,使得传统的抽水试验、水位回复试验等方法确定含水层参数存在困难且花费较大。通过对广西北海市滨海含水层地下水位动态资料进行分析,发现其上升段和下降段是不对称的。基于海岸带承压含水层正弦潮汐波的传播理论,提出了确定含水层参数的分段法,并与振幅衰减法和滞后时间法进行对比,各种方法求出的储水系数与导水系数之比(S/T)很接近,说明分段法是有效的。对于北海市滨海含水层,上升段求出的S/T值比下降段要大,其成因机理还有待进一步分析。     

Regional hydrostratigraphy and groundwater flow modeling in the arsenic-affected areas of the western Bengal basin, West Bengal, India

The first documented interpretation of the regional-scale hydrostratigraphy and groundwater flow is presented for a ～21,000-km² area of the arsenic-affected districts of West Bengal [Murshidabad, Nadia, North 24 Parganas and South 24 Parganas (including Calcutta)], India. A hydrostratigraphic model demonstrates the presence of a continuous, semi-confined sand aquifer underlain by a thick clay aquitard. The aquifer thickens toward the east and south. In the south, discontinuous clay layers locally divide the near-surface aquifer into several deeper, laterally connected, confined aquifers. Eight 22-layer model scenarios of regional groundwater flow were developed based on the observed topography, seasonal conditions, and inferred hydrostratigraphy. The models suggest the existence of seasonally variable, regional, north–south flow across the basin prior to the onset of extensive pumping in the 1970s. Pumping has severely distorted the flow pattern, inducing high vertical hydraulic gradients across wide cones of depression. Pumping has also increased total recharge (including irrigational return flow), inflow from rivers, and sea water intrusion. Consequently, downward flow of arsenic contaminated shallow groundwater appears to have resulted in contamination of previously safe aquifers by a combination of mechanical mixing and changes in chemical equilibrium.     

Characterization of the physical and hydraulic properties of the sediment in karst aquifers of the Springfield Plateau,Central Missouri,USA

Physical and hydraulic properties of sediment from two karst aquifers were measured to determine (1) the similarity of sediment between karst aquifer systems and (2) the importance of sediment in modeling flow through karst aquifers. The sediment from the two systems was similar in size and composition. Within both aquifers, the silt-sized sediment was composed primarily of quartz, with minor amounts of plagioclase and clays. Hydraulic conductivity of the sediment measured directly (falling-head test) ranged from 1.61×10⁻⁷ to 1.33×10⁻⁶ m s^–1 and estimated using the Campbell equation ranged from 8.30×10⁻⁸ to 8.98×10⁻⁷ m s^–1. These values of hydraulic conductivity fall within the span of hydraulic conductivities for carbonate rocks, indicating that the sediment and carbonate matrix could be represented as one mathematical unit in modeling flow through karst aquifers. Statistical agreement in the hydraulic conductivity values generated by the two methods indicates that the estimation technique could be used to calculate hydraulic conductivities; thus allowing karst scientist to collect bulk sediment samples instead of having to collect cores from within karst aquifers. Electronic Publication     

A comparative study of hydraulic conductivity estimations using geostatistics

Three approaches for estimating the hydraulic conductivity (K) of the Trifa aquifer, Morocco were investigated: (1) kriging of the K values obtained from pumping tests, (2) cokriging of the pumping test data with electrical resistivity data as a secondary variable, and (3) cokriging of the pumping test data with the slope of the water table. Gauss-transformed values of the variables are used because they provide more robust variograms and transformed values of the primary and secondary variables show correlations higher than the raw values, which is beneficial in cokriging. In cokriging with electrical resistivity, two zones are considered since the geological deposits are different from the north to the south of the aquifer, which is reflected in different correlations between the variables. Comparison of the three approaches is based mainly on the estimation errors, and to a lesser degree on the cross-validations of the corresponding variogram models and general considerations, like the measurements’ reliability and aquifer make-up. The best-estimated K is given by cokriging with the slope of the water table and is therefore preferred for further use in groundwater flow modeling. Thus, electrical resistivity or the slope of the water table can both be used as secondary variables to estimate K, especially in heterogeneous aquifers with lateral variations in lithology, as is the case of the Trifa aquifer.     

Vertical movement of water in a High Plains Aquifer induced by a pumping well

Field observation and numerical simulations were carried out to evaluate the hydraulic relationship between the shallow and deep aquifer of a High Plains Aquifer system, in which shallow and deep aquifers are separated by an aquitard. Pumping from the lower aquifer resulted in a small drawdown in the upper aquifer and a larger drawdown in the aquitard; pumping from the shallow aquifer caused a small drawdown in the aquitard and the deep aquifer. Analysis of pumping test data gives the values of the hydraulic conductivity of the aquitard and the deep aquifer. Long-term observation of groundwater levels in the shallow and deep aquifers showed that a strong downward hydraulic gradient was maintained during an irrigation season. Numerical simulations were used to calculate the induced leakage of water from the shallow to the deep aquifer. Water budget analyses suggested that after pumping continues for a couple of days, the leakage from the overlying layers begins to supply the majority of the withdrawal from the deep aquifer. However, the induced leakage from the upper shallow aquifer can travel only a few meters into the aquitard, and it can not reach the lower aquifer during a 90 day pumping period. The major portion of the induced leakage occurred during the pumping period, but a small leakage can continue as a residual effect after the pumping period. The vertical hydraulic conductivity of the aquitard plays a major role in partitioning the ratio of the induced leakage for the pumping and after-pumping periods.     

Determination of porosity and storage capacity of a calcareous aquifer (France) by correlation and spectral analyses of time series

Time-series analyses were used to investigate the relationships between barometric pressure changes, earth tides, and water-level fluctuations in a confined aquifer. The method was applied to data from the fractured aquifer at the Hydrogeological Experimental Site in Poitiers (France) and used to yield estimates of the aquifer’s storage capacity, porosity and barometric efficiency. The aim is to address relevance of these analyses for an aquifer showing both fracture draining and confined karstic flow in thin strata. Cross-correlation and spectral analysis are used to compare water-level head and atmospheric-pressure fluctuations. Porosity and storage capacity are calculated using this method and compared to results from petrophysical measurements and hydraulic pumping tests, respectively. The storage coefficients calculated by the time-series analyses are in agreement with those obtained by interpretation of the interference pumping tests. Conversely, porosity values calculated by time-series analyses are underestimated compared to those obtained by other methods.     

Monitoring artificially stimulated fluid movement in the Cretaceous Dakota aquifer,western Kansas

Aquifer properties can be evaluated by monitoring artificially stimulated fluid movements between wells, if the fluid is heated. Changes in the temperature profile recorded in observation wells indicate the flow path of the heated fluid, which in effect acts as a tracer. A fluid-flow experiment in the Cretaceous Dakota Formation at the Hodgeman County site, west-central Kansas, demonstrated the advantage of using the distributed optical-fiber temperature sensing method for monitoring transient temperature conditions in this hydrological application. The fluid flow in the aquifer was increased by producing water from a pumping well and injecting heated water in an injection well 13 m (43 ft) distant from the pumping well. The time-temperature series data obtained and compared with results from previous pumping tests point to interwell heterogeneity of the aquifer and to a zone in the sandstone aquifer of high hydraulic conductivity. However, the experiment would have allowed further clarification of aquifer heterogeneity and thermal properties if at least one observation well had been present between the injection and production wells. Electronic Publication     

Numerical groundwater-flow modeling to evaluate potential effects of pumping and recharge: implications for sustainable groundwater management in the Mahanadi delta region,India

Process-based groundwater models are useful to understand complex aquifer systems and make predictions about their response to hydrological changes. A conceptual model for evaluating responses to environmental changes is presented, considering the hydrogeologic framework, flow processes, aquifer hydraulic properties, boundary conditions, and sources and sinks of the groundwater system. Based on this conceptual model, a quasi-three-dimensional transient groundwater flow model was designed using MODFLOW to simulate the groundwater system of Mahanadi River delta, eastern India. The model was constructed in the context of an upper unconfined aquifer and lower confined aquifer, separated by an aquitard. Hydraulic heads of 13 shallow wells and 11 deep wells were used to calibrate transient groundwater conditions during 1997–2006, followed by validation (2007–2011). The aquifer and aquitard hydraulic properties were obtained by pumping tests and were calibrated along with the rainfall recharge. The statistical and graphical performance indicators suggested a reasonably good simulation of groundwater flow over the study area. Sensitivity analysis revealed that groundwater level is most sensitive to the hydraulic conductivities of both the aquifers, followed by vertical hydraulic conductivity of the confining layer. The calibrated model was then employed to explore groundwater-flow dynamics in response to changes in pumping and recharge conditions. The simulation results indicate that pumping has a substantial effect on the confined aquifer flow regime as compared to the unconfined aquifer. The results and insights from this study have important implications for other regional groundwater modeling studies, especially in multi-layered aquifer systems.     

Review: “Jacob’s Zoo”—how using Jacob’s method for aquifer testing leads to more intuitive understanding of aquifer characteristics

The interpretation of aquifer responses to pumping tests is an important tool for assessing aquifer geometry and properties, which are critical in the assessment of water resources or in environmental remediation. However, the responses of aquifers, measured by time-drawdown relationships in monitoring wells, are nonunique solutions that are affected by many factors. Jacob’s Zoo is a collection of graphical interpretations that allows students and practitioners to develop an intuitive feel for how natural hydrogeological systems work, and develop a set of skills that provide a better understanding of aquifer properties far beyond interpretation of pumping tests. Jacob’s Zoo, based on the work of Jacob (1950), fosters a deeper understanding, although few practitioners realize the full utility of the method. Jacob CE (1950) Flow of groundwater, In: Rouse H (ed) Engineering Hydraulics, Wiley, New York. P 321–386.

     

Field experiments and numerical simulations of whirlpool foundation pit dewatering

A whirlpool foundation pit is a small-diameter, deep circular pit. Because of its depth and small diameter, a large drawdown is required, and a limited number of wells can be installed inside the pit. During excavation, partially penetrating wells inside and outside the foundation pit have to be installed to lower the water level when the aquifer is too thick. However, partially penetrating wells near partially penetrating curtains cannot be treated by analytical methods. Therefore, it is necessary to use numerical methods to predict dewatering during excavation. Field experiments were performed on whirlpool foundation pit 1880 of Baosteel Group, Shanghai, China, to obtain pumping rates and drawdown, pumping with a single well and two wells in the confined aquifer. The results indicate that the drawdown inside the pit induced by pumping wells outside the foundation pit was small, whereas it was large for pumping wells inside the pit. The pumping wells inside and outside the pit had to be combined to lower the water level. A three-dimensional numerical model was developed to simulate the dewatering process. The hydraulic conductivities of the confined aquifers were inversed by using the pumping tests. Operation schedules were simulated with the corrected model for different combinations of wells inside and outside the pit. The results suggest that different schedules and operation conditions affect drawdown. The monitored results during dewatering indicate that the simulation and field measurements were in agreement. The results can be applied to similar situations.