In-well time-of-travel approach to evaluate optimal purge duration during low-flow sampling of monitoring wells

A common assumption with groundwater sampling is that low (<0.5 L/min) pumping rates during well purging and sampling captures primarily lateral flow from the formation through the well-screened interval at a depth coincident with the pump intake. However, if the intake is adjacent to a low hydraulic conductivity part of the screened formation, this scenario will induce vertical groundwater flow to the pump intake from parts of the screened interval with high hydraulic conductivity. Because less formation water will initially be captured during pumping, a substantial volume of water already in the well (preexisting screen water or screen storage) will be captured during this initial time until inflow from the high hydraulic conductivity part of the screened formation can travel vertically in the well to the pump intake. Therefore, the length of the time needed for adequate purging prior to sample collection (called optimal purge duration) is controlled by the in-well, vertical travel times. A preliminary, simple analytical model was used to provide information on the relation between purge duration and capture of formation water for different gross levels of heterogeneity (contrast between low and high hydraulic conductivity layers). The model was then used to compare these time–volume relations to purge data (pumping rates and drawdown) collected at several representative monitoring wells from multiple sites. Results showed that computation of time-dependent capture of formation water (as opposed to capture of preexisting screen water), which were based on vertical travel times in the well, compares favorably with the time required to achieve field parameter stabilization. If field parameter stabilization is an indicator of arrival time of formation water, which has been postulated, then in-well, vertical flow may be an important factor at wells where low-flow sampling is the sample method of choice.     

Detection of free vertical convection and double-diffusion in groundwater monitoring wells with geophysical borehole measurements

Two algorithms for in-situ detection and identification of vertical free convective and double-diffusive flows in groundwater monitoring wells or boreholes are proposed. With one algorithm the causes (driving forces) and with the other one the effects (convection or double-diffusion) of vertical transport processes can be detected based on geophysical borehole measurements in the water column. Five density-driven flow processes are identified: thermal, solutal, and thermosolutal convection leading to an equalization, as well as saltfingers and diffusive layering leading to an intensification of a vertical density gradient. The occurrence of density-driven transport processes could be proven in many groundwater monitoring wells and boreholes; especially shallow sections of boreholes or groundwater monitoring wells are affected dramatically by such vertical flows. Deep sections are also impaired as the critical threshold for the onset of a density-driven flow is considerably low. In monitoring wells or boreholes, several sections with different types of density-driven vertical flows may exist at the same time. Results from experimental investigations in a medium-scale testing facility with high aspect ratio (height/radius = 19) and from numerical modeling of a water column agree well with paramters of in-situ detected convection cells.     

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.     

Site characterization and monitoring of natural attenuation indicator parameters in a fuel contaminated coastal aquifer: Karaduvar (Mersin,SE Turkey)

This paper presents results from a site characterization and monitoring study at Karaduvar area (Mersin, SE Turkey), where high concentrations of refined petroleum products have been detected in domestic and irrigation water wells. The saturated and unsaturated zones in the deltaic aquifer are contaminated by large quantities of gasoline and diesel range fuel hydrocarbons (GRHs and DRHs) released from diverse sources that include accidental spills, storage tank fires, pipeline breaks, deliberate discharge of waste petroleum products from slop tanks and illegal tanker truck washing facilities. At the site, due to the complex nature of the pollution sources, overlapping contaminant plumes exist and cover an area of about 0.5 km². In both polluted and unpolluted parts of the aquifer, monitoring of groundwater physicochemical parameters in a total of 55 sampling points was carried out between 2006 and 2007. The results show that the terminal electron acceptors (e.g. dissolved oxygen, nitrate, Mn(IV), Fe(III), sulfate) were reduced near the source area(s) indicating presence of actively operating biodegradation processes at the site. Close to the contaminant source area(s), conditions in the plume are highly anoxic and reducing; where high amounts of transformation products (e.g. bicarbonate, dissolved iron, and manganese) are present in solution. Additionally, at the site, excessive pumping, careless land use, and deliberate wastewater discharges significantly deteriorated the quality and quantity of groundwater. Excessive groundwater pumping for industrial and agricultural uses has resulted in substantial water level declines (2–3 m) near the coastal part where seawater intrusion threatens the groundwater resources.     

Iron-hydroxide clogging of public supply wells receiving artificial recharge: near-well and in-well hydrological and hydrochemical observations

Clogging of water wells by iron-hydroxide incrustations due to mixing of anoxic and oxic groundwater is a common well-ageing problem. The relation between well operation (on and off), the spatial and temporal variations in hydrochemistry outside and inside a supply well, and the distribution of clogging iron-hydroxides were studied in an artificial recharge well field in the Netherlands. Camera inspection, high-resolution multi-level water sampling outside the well and detailed in-well pH/EC/O₂ profiles revealed remarkable patterns. During pumping, the top of the upper well screen abstracted oxic filtrate, although the larger part of the in-well water column was anoxic. The column rapidly turned oxic after shutdown due to a downward short-circuiting of oxic water via the well. Within 15 d it became anoxic due to the slow advance of anoxic lake filtrate created by local changes in flow direction as the neighboring wells continued to pump. Severe clogging occurred where the oxic filtrate entered the well, while half-clogging of the upper well screen occurred due to less inflow of oxic filtrate on the lake side. Transport of iron flocs and bacterial slimes after shutdown seemed to clog the lower part of the well screen. Frequent on/off switching should be avoided in iron-clogged wells.     

Dynamic process analyses of saline intrusion with over-exploitation by coupled flow and dispersive modelling

Understanding of saline intrusion into coastal aquifer is an important issue in management and protection of groundwater resource, which can be well achieved by groundwater modelling. To explain some phenomena of correlation between groundwater level and salinity in observation wells in coastal area, the authors compare the migration velocity of saline particles and transmission velocity of water pressure and derive analytical equations of these two velocities for plane and radial flows. The driving force and resistance of saline intrusion were analysed based on the analytical modelling. The destruction and reconstruction of the equilibrium between fresh water seepage towards sea and saline dispersion to inland were considered as an essence of whole intrusion process. The dynamic process of seepage and dispersion at different stages of saline intrusion were analysed under groundwater over-exploitation. The basic equations of saline intrusion were derived and the mechanism of transitional zone movement was discussed. These constitute coupled seepage–dispersion theory of saline intrusion, which becomes an important supplement to existing theory of saline intrusion. As a case study on saline intrusion in Guangrao County in Eastern China, this theory was applied to predict the development of saline intrusion in the study area.     

0-1整数规划在水源地开采井最优布局中的应用研究

新建水源地如何最优布井,对水源地开采运行具有重大影响。本文介绍了0-1整数规划及其求解方法,并结合一个实际问题进行说明。作者应用0-1整数规划方法建立了北京密怀顺平原区拟建水源地最优布井的地下水管理模型,采用响应矩阵法将地下水模拟模型和管理模型相耦合,运用分支-定界法求出最优井位,并将优化井位与原拟定井位的计算结果进行了对比。结果表明,0-1整数规划能很好地解决最优布井问题,对水源地规划、基坑降水、地下水人工回灌等方面都具有一定的指导意义。     

青岛李哥庄地下咸水体恢复方案的数值优化

海水入侵导致地下水变咸对青岛李哥庄地区经济发展产生了不利影响,急需开展治理该地区咸水体的工作。根据研究区已有的降雨、蒸发及水文地质等资料,结合现场调查和监测,利用地下水模拟软件Visual-MODFLOW中的SEAWAT模块建立了该地区的数值模型,并利用实测资料识别和验证了该模型参数,表明所建模型能反映李哥庄地区的实际水文地质条件。为了模拟研究区地下咸水体的恢复治理情况,提出了连续抽水、间歇式抽水和抽注水结合三种咸水恢复方案,并利用该模型对三种恢复方案进行了模拟优化。结果表明,抽取地下咸水只是咸水恢复的一个因素,周边淡水的驱替作用则是咸水恢复另一个的因素,抽注水结合恢复方案为三种方案中最优的方案。     

Recovery of groundwater in the Sanriku region contaminated by the tsunami inundation from the 2011 Tohoku earthquake

Many wells in the Sanriku region used as sources for water supply systems were heavily contaminated by the tsunami of the 2011 great Tohoku earthquake on March 11 in 2011. To better understand the nature of the groundwater contamination by the tsunami inundation and to clarify the recovery process of contaminated groundwater at the study wells, groundwater monitoring has been conducted once or twice yearly since early summer in 2011. High and abnormal values of electric conductivity (EC), chloride ion concentration (CIC), Na⁺, Ca⁺, heavy metal ions, and heavier isotopes of the contaminated groundwater were also obtained in April and June 2011. The chemical elements have rapidly and exponentially decreased as a result of effective pumping of the contaminated groundwater from the study wells and because of abundant rainfall in 2011. In April 2015 (about 4 years after the tsunami inundation), the CIC and EC of the contaminated groundwater of two study wells in Minamisanriku town had reached pre-inundation values. The estimated residence times of groundwater of the two study wells were 105–118 days in the full-day pumping stage and 910–1000 days in the daytime-only pumping stage.     

Fluids’ dynamics in transient air sparging of a heterogeneous unconfined aquifer

Water table dynamics, dissolved oxygen (DO) content, electrical resistivity (ER) in monitoring wells and air pressure in the vadose zone are monitored in air sparging (AS) accompanied by soil vapor extraction (SVE) at a hydrocarbon-contaminated groundwater site in Oman, where a diesel spillover affected a heterogeneous unconfined aquifer. The formation of a groundwater mound at the early stage of air injection and potential lateral migration of contaminants from the mound apex called for an additional hydrodynamic barrier constructed as a pair of pump-and-treat (P&T) wells whose recirculation zone encompassed the AS and SVE wells. In all monitored piezometers the phreatic surface showed a rapid and distinct peak, which is attributed to the time of air breakthrough from the injection point to the vadose zone and a relatively mild recession limb interpreted as a decay of the mound. Tracer tests showed a layer of a relatively low hydraulic conductivity at an intermediate depth of the screened interval of the wells. Increased levels of DO and borehole air pressure that have been observed (as far as 50 m away) are likely mitigated by SVE and P&T. Radius of influence can be indirectly inferred from ER and DO changes in the AS operation zone. Salt tracer tests have shown that groundwater velocity within the AS zone decreases with the increase of air injection rate.     

Chemical and mechanical clogging of groundwater abstraction wells at well field Heel,the Netherlands

Well field Heel, in the south east of the Netherlands, consists of a row of wells drilled in an anoxic pyrite-containing aquifer alongside a former gravel pit, which now serves as a recharge basin, where water is actively aerated. All wells are seriously affected by chemical (screen slot) and/or mechanical (well bore) clogging. The objective of this study is to explain this combined occurrence. A combination of chemical, hydraulic and well-maintenance data indicate three groundwater quality types: (1) oxic basin water, (2) anoxic iron-containing basin water after oxidation of the traversed aquifer, and (3) deeply anoxic native groundwater. Wells abstracting a mixture of oxic basin water and anoxic basin water and/or native groundwater experience chemical well clogging, whereas wells abstracting (only or partly) native groundwater are vulnerable to mechanical well clogging. In the end, after oxic basin water has completely oxidized the traversed the aquifer, only two groundwater quality types will be present. Wells abstracting only oxic basin water will show no clogging, and wells abstracting a mixture of native groundwater and oxic basin water will experience chemical and possibly also mechanical well clogging. In this reasoning, the sequence in abstracted groundwater quality types coincides with a sequence in well clogging: from mechanical to chemical to no clogging. As well field Heel is situated in sloping terrain, the interplay between regional hydraulic gradient and different water qualities results in one-sided chemical clogging in the upper part of the well screen during abstraction, and in the lower part during the resting phase.     

Effects of slow recovery rates on water column geochemistry in aquitard wells

Monitoring wells are often installed in aquitards to verify effectiveness for preventing migration of surface contaminants to underlying aquifers. However, water sampling of aquitard wells presents a challenge due to the slow recovery times for water recharging the wells, which can take as long as weeks, months or years to recharge depending on the sample volume needed. In this study, downhole profiling and sampling of aquitard wells was used to assess geochemical changes that occur in aquitard wells during water level recovery. Wells were sampled on three occasions spanning 11 years, 1 year and 1 week after they were purged and casing water showed substantial water chemistry variations. Temperature decreased with depth, whereas pH and specific conductance increased with depth in the water column after 11 years of water level recovery. Less stable parameters such as dissolved O₂ (DO) and Eh showed strong zonation in the well column, with DO stratification occurring as the groundwater slowly entered the well. Oxidation of reduced till groundwater along with degassing of CO₂ from till pore water affects mineral solubility and dissolved solid concentrations. Recommendations for sampling slowly recovering aquitard wells include identifying the zone of DO and Eh stratification in the well column and collecting water samples from below the boundary to better measure unstable geochemical parameters.     

基于TOUGH2数值模拟软件的延津地区热储采灌研究

本文采用TOUGH2数值模拟软件,根据研究区对井在预设情景下运转30年的模拟结果,对河南省延津县下新近系热储层压力场、温度场进行模拟,在模型验证的基础上,对不同回灌率及不同采灌井间距工况进行模拟研究,结果显示:回灌井的回灌过程会使地下水水温降低,回灌时间越长水温降低幅度越大,影响范围也越大,根据模拟抽水井和回灌井不同间距条件下(55 m、150 m、250 m、300 m)温度场变化情况,建议区域下新近系地层地热回灌项目抽水井和回灌井间距应大于等于250m,以保证抽水井出水温度;而不同回灌率工况下分析化学场和压力场的模拟数据可以看出,回灌率越大,地热水水质、压力影响越大,基本呈线性关系。为区域地热资源回灌-开采提供示范和指导。     

An analytical study on groundwater flow in drainage basins with horizontal wells

Analytical studies on release/capture zones are often limited to a uniform background groundwater flow. In fact, for basin-scale problems, the undulating water table would lead to the development of hierarchically nested flow systems, which are more complex than a uniform flow. Under the premise that the water table is a replica of undulating topography and hardly influenced by wells, an analytical solution of hydraulic head is derived for a two-dimensional cross section of a drainage basin with horizontal injection/pumping wells. Based on the analytical solution, distributions of hydraulic head, stagnation points and flow systems (including release/capture zones) are explored. The superposition of injection/pumping wells onto the background flow field leads to the development of new internal stagnation points and new flow systems (including release/capture zones). Generally speaking, the existence of n injection/pumping wells would result in up to n new internal stagnation points and up to 2n new flow systems (including release/capture zones). The analytical study presented, which integrates traditional well hydraulics with the theory of regional groundwater flow, is useful in understanding basin-scale groundwater flow influenced by human activities.     

Groundwater resource sustainability in the Wadi Watir delta, Gulf of Aqaba, Sinai, Egypt

The Wadi Watir delta, in the arid Sinai Peninsula, Egypt, contains an alluvial aquifer underlain by impermeable Precambrian basement rock. The scarcity of rainfall during the last decade, combined with high pumping rates, resulted in degradation of water quality in the main supply wells along the mountain front, which has resulted in reduced groundwater pumping. Additionally, seawater intrusion along the coast has increased salinity in some wells. A three-dimensional (3D) groundwater flow model (MODFLOW) was calibrated using groundwater-level changes and pumping rates from 1982 to 2009; the groundwater recharge rate was estimated to be 1.58?×?10⁶ m³/year. A variable-density flow model (SEAWAT) was used to evaluate seawater intrusion for different pumping rates and well-field locations. Water chemistry and stable isotope data were used to calculate seawater mixing with groundwater along the coast. Geochemical modeling (NETPATH) determined the sources and mixing of different groundwaters from the mountainous recharge areas and within the delta aquifers; results showed that the groundwater salinity is controlled by dissolution of minerals and salts in the aquifers along flow paths and mixing of chemically different waters, including upwelling of saline groundwater and seawater intrusion. Future groundwater pumping must be closely monitored to limit these effects.     

Application of multivariate statistical procedures on major ions and trace elements in a multilayered coastal aquifer: the case of the south Rhodope coastal aquifer

Multivariate statistical techniques including cluster analysis and principal components analysis were applied on 22 variables consisted of 3 physicochemical parameters, 8 major ions and 11 trace elements. Samples were collected from the south Rhodope multilayered coastal aquifer in north Greece which is facing saltwater intrusion and anthropogenic contamination over the last 35 years. Cluster analysis grouped the variables into five main groups while principal components analysis revealed four distinct hydrochemical processes in the aquifer system, explaining 84.5 % of the total variance between the variables. The identified processes correspond to, saltwater intrusion and subsequent reverse cation exchange, the presence of deep connate groundwater masses, application of fertilizers in shallow wells and anthropogenic contamination with heavy metals nearby an improperly constructed landfill. The wells categorized with the above techniques were grouped and five constituent ratios Na/Cl, (Mg + Ca)/Cl, Ca/(HCO₃ + SO₄), Ca/SO₄ and Ca/Mg were utilized to identify the ones which enable the more accurate distinction between the group cases. The results of stepwise discriminant analysis showed that the calculated classification function can distinguish almost 80 % of groundwater samples with the Na/Cl ratio being the most statistically significant grouping variable. All the aforementioned statistical models managed to successfully identify numerous hydrochemical processes in a complex multilayered aquifer system and to explicitly attribute them for every investigated well, allowing a deeper insight into groundwater chemical characteristics with the use of an optimized smaller number of variables.     

Application of Numerical Modeling for Groundwater Flow System Analysis in the Akaki Catchment,Central Ethiopia

A three dimensional steady-state finite difference groundwater flow model is used to quantify the groundwater fluxes and analyze the subsurface hydrodynamics in the Akaki catchment by giving particular emphasis to the well field that supplies water to the city of Addis Ababa. The area is characterized by Tertiary volcanics covered with thick residual and alluvial soils. The model is calibrated using head observations from 131 wells. The simulation is made in a two layer unconfined aquifer with spatially variable recharge and hydraulic conductivities under well-defined boundary conditions. The calibrated model is used to forecast groundwater flow pattern, the interaction of groundwater and surface water, and the effect of pumping on the well field under different scenarios. The result indicates that the groundwater flows regionally to the south converging to the major well field. Reservoirs and rivers play an important role in recharging the aquifer. Simulations made under different pumping rate indicate that an increase in pumping rate results in substantial regional groundwater level decline, which will lead to the drying of springs and shallow hand dug wells. Also, it has implications of reversal of flow from contaminated rivers into productive aquifers close to main river courses. The scenario analysis shows that the groundwater potential is not enough to sustain the ever-growing water demand of the city of Addis Ababa. The sensitivity and scenario analysis provided important information on the data gaps and the specific sites to be selected for monitoring, and may be of great help for transient model development. This study has laid the foundation for developing detailed predictive groundwater model, which can be readily used for groundwater management practices.     

Enhanced pumping test using physicochemical tracers to determine surface-water/groundwater interactions in an alluvial island aquifer,river Rhône,France

Groundwater extracted from the Barthelasse Island aquifer, surrounded by the river Rhône (southeastern France), contributes to the drinking water supplies of 180,000 inhabitants. Owing to its location close to the river and the presence of two backwaters (oxbow lakes), the pumped groundwater is highly vulnerable to river pollution. A pumping test was conducted over 24 h to analyse and quantify the water exchange processes between the river, backwaters and groundwater. During the pumping test, isotopic (δ¹⁸O, δ²H and ²²²Rn), hydrochemical and hydrophysical monitoring of the groundwater was undertaken. Hydraulic heads were measured in pumping wells and at a piezometer located between the wells. Discrete water samples were collected at several observation points in the field, including the backwater and river. The results show mixing between three end-members, as defined by the deuterium excess and silica concentration, led by river Rhône water which had been affected by water–rock interactions over time and mixing with surface evaporated waters. The pumped water resulted from mixing between three end-members, all of which depended on the river Rhône but differed in terms of residence time in the system. Although the groundwater pumping wells are close to each other (<70 m) and have similar depths, the changes in the contributions from end-member waters at each well were different during the pumping test. Comparing isotopic tracers and geochemistry made it possible to quantify the different hydrological compartments that contribute to the groundwater pumped from the boreholes, which is critical in constructing a conceptual flow model.

     

Temporal trends in concentrations of DBCP and nitrate in groundwater in the eastern San Joaquin Valley,California, USA

Temporal monitoring of the pesticide 1,2-dibromo-3-chloropropane (DBCP) and nitrate and indicators of mean groundwater age were used to evaluate the transport and fate of agricultural chemicals in groundwater and to predict the long-term effects in the regional aquifer system in the eastern San Joaquin Valley, California. Twenty monitoring wells were installed on a transect along an approximate groundwater flow path. Concentrations of DBCP and nitrate in the wells were compared to concentrations in regional areal monitoring networks. DBCP persists at concentrations above the US Environmental Protection Agency’s maximum contaminant level (MCL) at depths of nearly 40 m below the water table, more than 25 years after it was banned. Nitrate concentrations above the MCL reached depths of more than 20 m below the water table. Because of the intensive pumping and irrigation recharge, vertical flow paths are dominant. High concentrations (above MCLs) in the shallow part of the regional aquifer system will likely move deeper in the system, affecting both domestic and public-supply wells. The large fraction of old water (unaffected by agricultural chemicals) in deep monitoring wells suggests that it could take decades for concentrations to reach MCLs in deep, long-screened public-supply wells, however.     

Effectiveness of time-series analysis for thermal plume propagation assessment in an open-loop groundwater heat pump plant

Thermal perturbation produced in the subsurface by open-loop groundwater heat pumps (GWHPs) represents a complex transport phenomenon that is affected by several factors, including intrinsic characteristics of the exploited aquifer, abstraction and reinjection well features, and the temporal dynamics of the accessed groundwater. Post-GWHP water may have become warmed or cooled before being reinjected into the aquifer, thereby creating a thermal plume, known as the thermal affected zone (TAZ), which can alter aquifer temperature. The TAZ is propagated mainly by advection, after which the plume tends to degrade via conductive heat transport and convection within moving water. Groundwater monitoring and multiparametric probes are used to check the dynamics of plume propagation and whether a system’s thermal plumes are generating unsuitable interference with wells, subsurface infrastructure, or land use. Analyses of time-series groundwater monitoring data can be used to monitor TAZ movement. In this paper, the thermal plume velocity was calculated by both an analytical solution and cross-correlation. Cross-correlation calculated between temperature measured in the reinjection well and control downstream piezometers can reveal plume dynamics and demonstrate the importance of advective transport in aquifer heat transfer.