Ambient well-bore mixing, aquifer cross-contamination, pumping stress, and water quality from long-screened wells: What is sampled and what is not?

Water quality tests were performed on two long-screened alluvial aquifer wells (15–30 m of screen) that had been completed in a heterogeneous aquifer that exhibits extreme temporal water quality variability. When stressed, the total dissolved solids (TDS) in one well decreased from 10,600 to 3,500 mg/L and in another well the TDS increased from 136 to 2,255 mg/L. Nested short-screened monitoring wells were constructed in chemically distinct horizons affecting each well. Water level measurements and solute and isotopic samples were obtained from the production wells and the monitoring wells during a water quality test. Results of a time drawdown tests demonstrate transmissivity differences between horizons. Ambient water quality in the production wells and aquifer cross-contamination are controlled by well-bore mixing due to head differences of as little as 0.01 m between chemically distinct horizons which are linked by the production well screen. During non-stress periods, the ambient well-bore chemistry is controlled by the horizon with the greatest hydraulic head, whereas during stressed conditions, horizon transmissivity controls the well-bore chemistry. In one well, aquifer cross-contamination, driven by an ambient head differential of 1.2 m, persisted until about 1,600 well-bore volumes were purged.     

Characterization of fracture connectivity in a siliciclastic bedrock aquifer near a public supply well (Wisconsin, USA)

In order to protect public supply wells from a wide range of contaminants, it is imperative to understand physical flow and transport mechanisms in the aquifer system. Although flow through fractures has typically been associated with either crystalline or carbonate rocks, there is growing evidence that it can be an important component of flow in relatively permeable sandstone formations. The objective of this work is to determine the role that fractures serve in the transport of near-surface contaminants such as wastewater from leaking sewers, to public supply wells in a deep bedrock aquifer. A part of the Cambrian aquifer system in Madison, Wisconsin (USA), was studied using a combination of geophysical, geochemical, and hydraulic testing in a borehole adjacent to a public supply well. Data suggest that bedrock fractures are important transport pathways from the surface to the deep aquifer. These fractured intervals have transmissivity values several orders of magnitude higher than non-fractured intervals. With respect to rapid transport of contaminants, high transmissivity values of individual fractures make them the most likely preferential flow pathways. Results suggest that in a siliciclastic aquifer near a public supply well, fractures may have an important role in the transport of sewer-derived wastewater contaminants.     

Hydrogeological studies on the Nubian sandstone aquifer in El-Bahariya Oasis, Western Desert, Egypt

The current study aimed to evaluate hydrogeologically the Nubian sandstone aquifer in El-Bahariya Oasis. It represents the main water-bearing horizon in the study area and consists of continental elastic sediments, mainly sandstone alternating with shale and clays. The general flow lines are directed from SW to NE direction, as detected from the constructed potentiometric head contour map. The piezometric surface reaches 149 m in El-Heiz area at the southern part, while it reaches 90 m at the northern, reflecting higher pressure head of the aquifer in the southern part. The map also illustrates that the southern part is considered as the most promising location for development. The structural elements play an important role in the deposition and distribution of the sedimentary succession of the Nubian sandstone sediments. Consequently, this sedimentary pattern affects the occurrences and movements of the groundwater within the aquifer system. Along the structurally high areas, in the study area, the piezometric head increases, while the reverse is recorded along the structurally low areas. The step-drawdown tests data were carried out by calculating the aquifer loss coefficient (B) and the well loss constant (C). The B values are smaller compared with C values, indicating that the aquifer under pressure has a behavior of leaky aquifer; therefore, it shows hydraulic connection with surrounding formation. The values of well efficiency range from 78.50% to 87.76%. Analysis of 12 pumping test data (constant discharge tests) was carried out in order to calculate the Nubian aquifer hydraulic parameters (transmissivity, hydraulic conductivity, and storage coefficient). The transmissivity values decrease from 3,045 m²/day in the southern part (El-Heiz area) to 236 m²/day in the northeastern part (El-Harra area). Accordingly, the aquifer classified as a high to moderate potentiality. Transmissivity contour map observes gradual increase of transmissivity values from the southern to northeastern direction. This may be due to the increase of shale or clay content in the concerned aquifer in that direction. The storage coefficient values range between 1.04 × 10^?4 and 5.22 × 10^?3, as obtained from the results of pumping test analysis, which ensure that the Nubian sandstone aquifer is classified as semi-confined to confined aquifer type. The S values show a decrease from southwest to northeast direction as detected from S-map. The hydraulic conductivity values vary from to 0.46 m/day in the northern part to 10.88 m/day in the southern part with an average of 5.67 m/day. According to the classification based on K values, the aquifer is mainly composed of coarse sand.     

Modeling approaches and strategies for data-scarce aquifers: example of the Dar es Salaam aquifer in Tanzania

Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.     

自流井有限差分模拟的校正模型

自流井是井孔-含水层系统的一种典型水文地质现象, 但在地下水资源评价和管理模型中研究较少.目前地下水数值模拟对自流井的3种处理方法均不合理, 也缺乏验证.常规地下水流有限差分法既不能准确地求出抽水井的井中水位, 也不能准确地求出自流井的流量, 必须进行校正.考虑自流井周围的径向流, 通过与井点附近地下水流的解析解相结合, 提出了矩形网格有限差分模拟的自流井校正模型, 其中考虑了自流井周围含水层存贮量的释放.算例表明, 即使用千米级大尺度网格, 本校正模型也能够获得相当精确的自流井流量, 绝大多数情况下相对误差小于5%.      

Artesian Flowing Wells Field of Phu Tok Aquifer

Phu tok Aquifer is the most productive water bearing rock in the Northeastern region of Thailand. Generally, well drilled in this aquifer yields approximately 20–50 m3/hr with some wells yield over 150 m3/hr. This aquifer characterized both unconfined and confined aquifer. Data collected from core sampling and TV bore hole show water bearing fractures in the fine grained massive sandstone intercalated with thin shale layers. These intermittently bedding plane fracture zones were found vertically at 60-65 and 95-110 m below ground surface. Three directions of vertically fractures were also found at 103-104 m in the test well with underlain densely sandstone until 120 m of depth. Wells drilled to 60-70 m fractures at the area where surface elevation between 150 and 170 m above mean sea level were artesian well with rising water about 1-3 m. Natural discharge rate from 2-5 inches of diameter casing is 5-10 m3/hr. Fractures at 90-110 m yield 150 m3/hr through 6 inch casing well with 5-6 m of risen water above the ground. It is found from pumping test of Phu Tok Aquifer that hydraulic conductivity of unconfined aquifer at shallower than 50 m is 0.005-17 m/d. Transmissivity and storage coefficient are 0.05-20 m2/d and 7×10-3-0.725 respectively. The confined aquifer at depth not exceed 90 m has hydraulic conductivity value of 0.2-10 m/d while transmissivity and storage coefficient are 3.19-150 m2/d and 1×10-10-1.6×10-2. Another confined aquifer at 90-120 m of depth have hydraulic conductivity value of 0.08-15 m/d and transmissivity and storage coefficient values of 1.7-178 m2/d and 4×10-7-4.5×10-3 respectively.     

论“单位涌水量就是导水系数”成立的特殊条件

针对截取侧向径流的抽水井,兰太权提出一个解释单位涌水量的公式,认为“单位涌水量就是导水系数”(q=T),并以砂箱试验的结果作为证据。然而,这个等式的推导过程缺少严谨的水动力学基础。砂箱试验的边界条件等价于含有1个或2个定水头边界的矩形承压含水层,单井抽水的稳定流场可以用镜像井点叠加法求解,从理论上表明单位涌水量与导水系数成正比,而比例系数取决于含水层的形状和抽水井的半径。在特定条件下,含水层和抽水井的尺寸可以造成单位涌水量与导水系数相等的情况,兰太权发表于文献中的2组砂箱试验恰好符合这种条件,但并不能证明q=T具有普适性。新的砂箱试验避免了几何尺寸的这种特殊组合,得到了显著偏离q=T的结果。     

Aquifer system response to intensive pumping in urban areas of the Gangetic plains,India: the case study of Patna

A significant component of domestic demand for water of urban areas located in the Gangetic plains is met by heavy pumping of groundwater. The present study is focused on the Patna municipal area, inhabited by 17 million people and spanning over 134 km², where entire urban water demand is catered from pumping by wells of various capacities and designs. The present study examines the nature of the aquifer system within the urban area, the temporal changes in the water/piezometric level and the recharge mechanism of the deeper aquifers. The aquifer system is made up of medium-to-coarse unconsolidated sand, lying under a ~40-m-thick predominantly argillaceous unit holding 8- to 13-m-thick localised sand layers and continues up to 220 m below ground. Groundwater occurs under semi-confined condition, with transmissivity of aquifers in 5,500–9,200 m² day^?1 range. Hydraulic head of the deeper aquifer remains in 9–19 m range below ground, in contrast to 1–9 m range of that of the upper aquitard zone. The estimated annual groundwater extraction from the deeper aquifer is ~212.0 million m³, which has created a decline of 3.9 m in the piezometric level of the deeper aquifer during the past 30 years. Unregulated construction of deep tube wells with mushrooming of apartment culture may further exacerbate the problem. The sand layers within the aquitard zone are experiencing an annual extraction of 14.5 million m³ and have exhibited stable water level trend for past one and half decades. This unit is recharged from monsoon rainfall, besides contribution from water supply pipe line leakage and seepage from unlined storm water drains.     

A conceptual and numerical model for groundwater management: a case study on a coastal aquifer in southern Tuscany, Italy

Ongoing hydrogeological research aims to develop a correct management model for the Plio-Pleistocene multi-aquifer system of the Albegna River coastal plain (southern Tuscany, Italy); overexploitation of this aquifer for irrigation and tourism has caused seawater intrusion. The conceptual model is based on field and laboratory data collected during the 1995–2003 period. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer. Natural outflow occurs through a diffuse flow into the sea and river; artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage, closing of aquitard) and artificially (multiscreened wells). The aquifer was represented by a three-dimensional finite element model using the FEFLOW numerical code. The model was calibrated for steady-state and transient conditions by matching computed and measured piezometric levels (February 1995–February 1996). The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it occurs above all in the Osa-Albegna sector, as well as along the river that at times feeds the aquifer. The effects of hypothetical aquifer exploitation were assessed in terms of water budget and hydraulic head evolution.     

A Simple Approach to Account for Radial Flow and Boundary Conditions When Kriging Hydraulic Head Fields for Confined Aquifers

The estimation and mapping of realistic hydraulic head fields, hence of flow paths, is a major goal of many hydrogeological studies. The most widely used method to obtain reliable head fields is the inverse approach. This approach relies on the numerical approximation of the flow equation and requires specifying boundary conditions and the transmissivity of each grid element. Boundary conditions are often unknown or poorly known, yet they impose a strong signature on the head fields obtained by inverse analysis. A simpler alternative to the inverse approach is the direct kriging of the head field using the measurements obtained at observation wells. The kriging must be modified to incorporate the available information. Use of the dual kriging formalism enables simultaneously estimating the head field, the aquifer mean transmissivity, and the regional hydraulic gradient from head data in steady or transient state conditions. In transient state conditions, an estimate of the storage coefficient can be obtained. We test the approach on simple analytical cases, on synthetic cases with solutions obtained numerically using a finite element flow simulator, and on a real aquifer. For homogeneous aquifers, infinite or bounded, the kriging estimate retrieves the exact solution of the head field, the exact hydrogeological parameters and the flow net. With heterogeneous aquifers, kriging accurately estimates the head field with prediction errors of the same magnitude as typical head measurement errors. The transmissivities are also accurately estimated by kriging. Moreover, if inversion is required, the kriged head along boundaries can be used as realistic boundary conditions for flow simulation.     

Groundwater modeling applied to production well management in a contaminated aquifer

A study was conducted to evaluate production strategies for a well field system near a source of groundwater contamination. Numerical modeling of groundwater flow was employed to generate hydraulic head configurations for different production scenarios. For a given scenario, an evaluation of contamination susceptibility was made by comparing head distributions in two aquifer units to the positions of the contaminant source and discharging water supply wells. The results of this study suggest that groundwater flow modeling can be a useful technique for planning the production of water supply wells in aquifers at risk of contamination from anthropogenic pollution sources.     

抽水井周围饱水介质场三维水位动态观测及地下水流运动特征

利用设计建造的抽水试验装置,采取多项潜水及承压水的抽水试验,获得大量的抽水井周围饱水介质场三维动态观测数据。试验发现：无论是潜水还是承压水状态,抽水井附近观测井的水位变化,均具有距抽水井愈远愈高、愈近愈低的漏斗状分布,抽水井附件的地下水位降落漏斗是显著的;同时发现每个观测井不同高度上的观测点的水位,均呈现上高下低的特征,说明地下水具有由上到下的运动分量,地下水由周边的供水边界向抽水井聚集运移时,运动方向是下斜向的,潜水如此,承压水也如此。分析认为,这是水头压力与水体自身重力同时起作用的结果。由此推断,在抽水时,靠近承压水含水层顶板存在强势水流,且越靠近补给源,这种强势水流越强。     

Strategic groundwater resources in the Tagliamento River basin (northern Italy): hydrogeological investigation integrated with geophysical exploration

The western sector of the Tagliamento River basin (Friuli Venezia-Giulia Region; northern Italy) is characterized by important water resources, both superficial and underground. In particular, in the Quaternary deposits of the plain, up to a depth of 500 m, six artesian aquifer systems exist. A large amount of lithological, geomorphological, and hydrogeological data is presented, allowing for definition of (a) the principal aquifer system of the area to a depth of about 500 m; (b) geometrical characteristics of the aquifers (thickness, lateral extension, etc); (c) the hydraulic parameters (hydraulic conductivity, transmissivity); (d) chemico-physical characteristics of the water; and (e) vulnerability to pollution of the aquifer systems. In a test area, where many boreholes were drilled for fresh water supply, the conceptual hydrogeological model was integrated by a detailed and repeated three-dimensional (3D) resistivity survey using combined electrical resistivity tomography (ERT) and time-domain electromagnetic (TDEM) soundings. ERT investigation was mainly used to get detailed information about geometry and porosity of the overburden and to calibrate the shallowest TDEM information.     

Simultaneous parameter identification of a heterogeneous aquifer system using artificial neural networks

An artificial neural network (ANN) model is proposed for the simultaneous determination of transmissivity and storativity distributions of a heterogeneous aquifer system. ANNs may be useful tools for parameter identification problems due to their ability to solve complex nonlinear problems. As an extension of previous study—Karahan H, Ayvaz MT (2006) Forecasting aquifer parameters using artificial neural networks, J Porous Media 9(5):429–444—the performance of the proposed ANN model is tested on a two-dimensional hypothetical aquifer system for transient flow conditions. In the proposed ANN model, Cartesian coordinates of observation wells, associated piezometric heads and observation time are used as inputs while corresponding transmissivity and storativity values are used as outputs. The training, validation and testing processes of the ANN model are performed under two scenarios. In scenario 1, all the sampled data are used through the simulation time. However, in the scenario 2, there are data gaps due to irregular observations. By using the determined synaptic network weights, transmissivity and storativity distributions are predicted. In addition, the performance of the proposed ANN is tested for different noise data conditions. Results showed that the developed ANN model may be used in simultaneous aquifer parameter estimation problems.     

地下水位变化对透—阻型岩溶塌陷影响的分析

文章从一维地下水运动和渗透力学的角度,分析比较潜水位上升与承压水位下降对岩溶地区透—阻型盖层中阻水层渗透稳定性的影响,重点讨论了承压水位下降速度(降速)与下降深度（降深）对阻水层中渗透坡降的影响,结果表明:地下水位变化（潜水位上升或承压水位下降）产生的非稳定渗流不利于岩溶洞穴开口上方阻水层的稳定,承压水位的下降对岩溶开口附近处阻水层中渗透力的影响远大于潜水位的变化;在承压层水位最大降深确定的条件下,承压水位下降速度愈快,岩溶开口附近处阻水层中向下渗透力愈大。因此,在覆盖型岩溶地区抽取地下水时,为了减缓或避免覆盖型塌陷的发生,应同时控制好最大降深和最大开采速度。      

Pattern Recognition in a Bimodal Aquifer Using the Normal-Score Ensemble Kalman Filter

The ensemble Kalman filter (EnKF) is now widely used in diverse disciplines to estimate model parameters and update model states by integrating observed data. The EnKF is known to perform optimally only for multi-Gaussian distributed states and parameters. A new approach, the normal-score EnKF (NS-EnKF), has been recently proposed to handle complex aquifers with non-Gaussian distributed parameters. In this work, we aim at investigating the capacity of the NS-EnKF to identify patterns in the spatial distribution of the model parameters (hydraulic conductivities) by assimilating dynamic observations in the absence of direct measurements of the parameters themselves. In some situations, hydraulic conductivity measurements (hard data) may not be available, which requires the estimation of conductivities from indirect observations, such as piezometric heads. We show how the NS-EnKF is capable of retrieving the bimodal nature of a synthetic aquifer solely from piezometric head data. By comparison with a more standard implementation of the EnKF, the NS-EnKF gives better results with regard to histogram preservation, uncertainty assessment, and transport predictions.     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.     

A numerical investigation to illustrate the consequences of hydraulic connections between granular and fractured-rock aquifers

Natural or artificial hydraulic connections between a granular aquifer in contact with a fractured-rock aquifer can have significant physical and chemical impacts at both a local and a regional scale. In this study, numerical simulations are conducted in order to illustrate the hydrogeological consequences of such hydraulic relationships. The numerical investigation, based on a conceptual model, focuses on the effects of the hydraulic connections when conducting a pumping test in a well that is set into a granular confined aquifer overlying a fractured-rock aquifer which presents a few fractures directly connected to the granular aquifer. It is illustrated that when interpreting the pumping test with the conventional methods consisting of plotting the drawdown versus time, a bias is introduced on the estimation of the transmissivity of the granular aquifer due to groundwater flowing from the fractured-rock aquifer via connecting fractures. However, it is underlined that plotting drawdown log-derivative versus time helps to diagnose the existence of these hydraulic relationships and therefore avoids committing a bias on the transmissivity estimation of the granular aquifer. Numerical simulation results also illustrate that hydraulic connections between the two aquifers can have significant impacts on the hydrogeochemical signature of the granular aquifer under investigation.     

Numerical simulation of a managed aquifer recharge system designed to supply drinking water to the city of Amsterdam,The Netherlands

Managed aquifer recharge (MAR) is increasingly used to secure drinking water supply worldwide. The city of Amsterdam (The Netherlands) depends largely on the MAR in coastal dunes for water supply. A new MAR scheme is proposed for the production of 10 × 10⁶ m³/year, as required in the next decade. The designed MAR system consists of 10 infiltration ponds in an artificially created sandbank, and 25 recovery wells placed beneath the ponds in a productive aquifer. Several criteria were met for the design, such as a minimum residence time of 60 days and maximum drawdown of 5 cm. Steady-state and transient flow models were calibrated. The flow model computed the infiltration capacity of the ponds and drawdowns caused by the MAR. A hypothetical tracer transport model was used to compute the travel times from the ponds to the wells and recovery efficiency of the wells. The results demonstrated that 98% of the infiltrated water was captured by the recovery wells which accounted for 65.3% of the total abstraction. Other sources include recharge from precipitation (6.7%), leakages from surface water (13.1%), and natural groundwater reserve (14.9%). Sensitivity analysis indicated that the pond conductance and hydraulic conductivity of the sand aquifer in between the ponds and wells are important for the infiltration capacity. The temperature simulation showed that the recovered water in the wells has a stable temperature of 9.8–12.5 °C which is beneficial for post-treatment processes. The numerical modelling approach is useful and helps to gain insights for implementation of the MAR.