Design Screening Tools for Passive Funnel and Gate Systems

The funnel and gate remediation concept (Star and Cherry 1993) represents a promising, yet relatively under-developed, technology for the passive control and in situ remediation of contaminated ground water. Effective design and implementation of such a system may, however, prove difficult under conditions of large or unpredictable variations in contaminant migration or ground water flow.
Numerical modeling of two-dimensional ground water flow has been used to predict the hydraulic performance of passive, straight, or winged funnel and gate configurations over a range of hydrogeologic and ambient ground water flow conditions. The results of these analyses were used to construct generic correlation diagrams relating upstream capture zone or gale through put to the barrier, gale, and aquifer characteristics. These diagrams serve as useful screening tools to (1) quantitatively estimate the capture zone of pre-determined funnel and gale configurations, or (2) develop preliminary funnel and gale designs that will yield a desired capture zone, independent of aquifer characteristics.     

Reactive barriers: hydraulic performance and design enhancements

The remediation of contaminated ground water is a multibillion-dollar global industry. Permeable reactive barriers (PRBs) are one of the leading technologies being developed in the search for alternatives to the pump-and-treat method. Improving the hydraulic performance of these PRBs is an important part of maximizing their potential to the industry. Optimization of the hydraulic performance of a PRB can be defined in terms of finding the balance between capture, residence time, and PRB longevity that produces a minimum-cost acceptable design. Three-dimensional particle tracking was used to estimate capture zone and residence time distributions. Volumetric flow analysis was used for estimation of flow distribution across a PRB and in the identification of flow regimes that may affect the permeability or reactivity of portions of the PRB over time. Capture zone measurements extended below the base of partially penetrating PRBs and were measured upgradient from the portion of aquifer influenced by PRB emplacement. Hydraulic performance analysis of standard PRB designs confirmed previously presented research that identified the potential for significant variation in residence time and capture zone. These variations can result in the need to oversize the PRB to ensure that downgradient contaminant concentrations do not exceed imposed standards. The most useful PRB design enhancements for controlling residence time and capture variation were found to be customized downgradient gate faces, velocity equalization walls, deeper emplacement of the funnel than the gate, and careful manipulation of the hydraulic conductivity ratio between the gate and the aquifer.     

Analytical Solutions for Flow Fields near Drain-and-Gate Reactive Barriers

Permeable reactive barriers (PRBs) are a popular technology for passive contaminant remediation in aquifers through installation of reactive materials in the pathway of a plume. Of fundamental importance are the degree of remediation inside the reactor (residence time) and the portion of groundwater intercepted by a PRB (capture width). Based on a two-dimensional conformal mapping approach (previously used in related work), the latter is studied in the present work for drain-and-gate (DG) PRBs, which may possess a collector and a distributor drain (“full” configuration) or a collector drain only (“simple” configuration). Inherent assumptions are a homogeneous unbounded aquifer with a uniform far field, in which highly permeable drains establish constant head boundaries. Solutions for aquifer flow fields in terms of the complex potential are derived, illustrated, and analyzed for doubly symmetric DG configurations and arbitrary reactor hydraulic resistance as well as ambient groundwater flow direction. A series of practitioner-friendly charts for capture width is given to assist in PRB design and optimization without requiring complex mathematics. DG PRBs are identified as more susceptible to flow divergence around the reactor than configurations using impermeable side structures (e.g., funnel-and-gate), and deployment of impermeable walls on drains is seen to mitigate this problem under certain circumstances.     

Combining pump-and-treat and physical barriers for contaminant plume control

A detailed analysis is presented of the hydraulic efficiency of plume management alternatives that combine a conventional pump-and-treat system with vertical, physical hydraulic barriers such as slurry walls or sheet piles. Various design settings are examined for their potential to reduce the pumping rate needed to obtain a complete capture of a given contaminated area. Using established modeling techniques for flow and transport, those barrier configurations (specified by location, shape, and length) that yield a maximum reduction of the pumping rate are identified assuming homogeneous aquifer conditions. Selected configurations are further analyzed concerning their hydraulic performance under heterogeneous aquifer conditions by means of a stochastic approach (Monte Carlo simulations) with aquifer transmissivity as a random space function. The results show that physical barriers are an appropriate means to decrease expected (mean) pumping rates, as well as the variance of the corresponding pumping rate distribution at any given degree of heterogeneity. The methodology presented can be transferred easily to other aquifer scenarios, provided some basic premises are fulfilled, and may serve as a basis for reducing the pumping rate in existing pump-and-treat systems.     

Efficiency verification of a horizontal flow barrier via flowmeter tests and multilevel sampling

The remediation strategy for an industrial site located in a coastal area involves a pump and treat system and a horizontal flow barrier (HFB) penetrating the main aquifer. To validate the groundwater flow conceptual model and to verify the efficiency of the remediation systems, we carried out piezometric measurements, slug tests, pumping tests, flowmeter tests and multilevel sampling. Flowmeter tests are used to infer vertical groundwater flow directions, and base exchange index is used to infer horizontal flow directions at a metric scale. The selected wells are located both upstream and downstream of the HFB. The installation of the HFB produced constraints to the groundwater flow. A stagnant zone of contaminated freshwater floating over the salt wedge in the upper portion of the aquifer is detected downstream of the HFB. This study confirms that the adopted remediation system is efficiently working in the area upstream of the HFB and even downstream in the bottom part of the aquifer. At the same time, it has also confirmed that hot spots are still present in stagnant zones located downstream of the HFB in the upper part of the aquifer, requiring a different approach to accomplish remediation targets. The integrated approach for flow quantification used in this study allows to discriminate the direction and the magnitude of groundwater fluxes near an HFB in a coastal aquifer. Copyright © 2012 John Wiley & Sons, Ltd.     

Capture and release zones of permeable reactive barriers under the influence of advective–dispersive transport in the aquifer

The problem of permeable reactive barrier (PRB) capture and release behavior is investigated by means of an approximate analytical approach exploring the invariance of steady-state solutions of the advection–dispersion equation to conformal mapping. PRB configurations considered are doubly-symmetric funnel-and-gate as well as less frequent drain-and-gate systems. The effect of aquifer heterogeneity on contaminant plume spreading is hereby incorporated through an effective transverse macro-dispersion coefficient, which has to be known. Results are normalized and graphically represented in terms of a relative capture efficiency M of contaminant mass or groundwater passing a control plane (transect) at a sufficient distance up-stream of a PRB as to comply with underlying assumptions. Factors of safety FS are given as the ratios of required capture width under advective–dispersive and purely advective transport for achieving equal capture efficiency M. It is found that M also applies to the release behavior down-stream of a PRB, i.e., it describes the spreading and dilution of PRB treated groundwater possibly containing incompletely remediated contamination and/or remediation reaction products. Hypothetical examples are given to demonstrate results.     

Time-related capture zones for radial flow in two dimensional randomly heterogeneous media

We consider the effect of randomly heterogeneous hydraulic conductivity on the spatial location of time-related capture zones (isochrones) for a non-reactive tracer in the steady-state radial flow field due to a pumping well in a confined aquifer. A Monte Carlo (MC) procedure is used in conjunction with FFT-based spectral methods. The log hydraulic conductivity field is assumed to be Gaussian and stationary, with isotropic exponential correlation. Various degrees of domain heterogeneity are considered and stability and accuracy of the MC procedure is examined. The location of an isochrone becomes uncertain due to heterogeneity, and it is strongly influenced by hydraulic conductivity variance. The probability that a particle released at a point in the aquifer is pumped by the well within a given time is identified. We propose a new expression for the probabilistic spatial distribution of isochrones, which is formally similar to the analytical solution for a uniform medium and takes into account the effects of heterogeneity.     

Conventional and combined pump-and-treat systems under nonuniform background flow

We investigate the performance of vertical hydraulic barriers in combination with extraction wells for the partial hydraulic isolation of contaminated aquifer areas. The potential advantage of such combinations compared to a conventional pump-and-treat system has already been demonstrated in a previous study. Here we extend the scope of the performance analysis to the impact of uncertainty in the regional flow direction as well as to highly heterogeneous aquifer transmissivity distributions. In addition, two new well-barrier scenarios are proposed and analyzed. The hydraulic efficiency of the scenarios is rated based on the expected (mean) reduction of the pumping rate that is required to achieve downgradient contaminant capture. The uncertain spatial distribution of aquifer transmissivity is considered by means of unconditioned Monte Carlo simulations. The significance of uncertain background flow conditions is incorporated by computing minimized pumping rates for deviations of the regional flow direction up to 30 degrees from a normative base case. The results give an answer on how pumping rates have to be changed for each barrier-well combination in order to achieve robust systems. It is exposed that in comparison to installing exclusively wells, the barrier-supported approach generally yields savings in the (average) pumping rate. The particular efficiency is shown to be highly dependent on the interaction of variance and integral scale of transmissivity distribution, well and barrier position, as well as direction of background flow.     

Cost Comparisons of Aquifer Heterogeneity Characterization Methods

The characterization of heterogeneity in hydraulic conductivity (K) is a major challenge for subsurface remediation projects. There are a number of field studies that compare the K estimates obtained using various techniques, but to our knowledge, no field‐based studies exists that compare the performance of estimated K heterogeneity fields or the associated characterization costs. In this paper, we compare the costs of characterizing the three‐dimensional K heterogeneity and its uncertainty estimates of a glaciofluvial aquifer‐aquitard sequence at a 15 m × 15 m × 18 m field site situated on the University of Waterloo campus. We compare geostatistical analysis of high resolution permeameter K data obtained from repacked core samples in five boreholes and hydraulic tomography analysis of four pumping tests consisting of up to 41 monitoring points per test. Aside from the comparison of costs, we also assess the performance of each method by predicting several pumping tests. Our analysis reveals that hydraulic tomography is somewhat more costly than the geostatistical analysis of high resolution permeameter K data due to the higher capital costs associated with the method. However, the equipment may be reused at other sites; hence these costs may be recovered over the life of the equipment. More significantly, hydraulic tomography is able to capture the most important features of the aquifer‐aquitard sequence leading to more accurate predictions of independent pumping tests. This suggests that more robust remediation systems may be designed if site characterization is performed with hydraulic tomography.     

The Horizontal Reactive Media Treatment Well (HRX Well®) for Passive In‐Situ Remediation

A new in‐situ remediation concept termed a Horizontal Reactive Media Treatment Well (HRX Well®) is presented that utilizes horizontal wells filled with reactive media to passively treat contaminated groundwater in‐situ. The approach involves the use of large‐diameter directionally drilled horizontal wells filled with granular reactive media generally installed parallel to the direction of groundwater flow. The design leverages natural “flow‐focusing” behavior induced by the high in‐well hydraulic conductivity of the reactive media relative to the aquifer hydraulic conductivity to passively capture and treat proportionally large volumes of groundwater within the well. Clean groundwater then exits the horizontal well along its downgradient sections. Many different types of solid granular reactive media are already available (e.g., zero valent iron, activated carbon, ion exchange resins, zeolite, apatite, chitin); therefore, this concept could be used to address a wide range of contaminants. Three‐dimensional flow and transport simulations were completed to assess the general hydraulic performance, capture zones, residence times, effects of aquifer heterogeneity, and treatment effectiveness of the concept. The results demonstrate that capture and treatment widths of up to tens of feet can be achieved for many aquifer settings, and that reductions in downgradient concentrations and contaminant mass flux are nearly immediate. For a representative example, the predicted treatment zone width for the HRX Well is approximately 27 to 44 feet, and contaminant concentrations immediately downgradient of the HRX Well decreased an order of magnitude within 10 days. A series of laboratory‐scale physical tests (i.e., tank tests) were completed that further demonstrate the concept and confirm model prediction performance. For example, the breakthrough time, peak concentration and total mass recovery of methylene blue (reactive tracer) was about 2, 35, and 20 times (respectively) less than chloride (conservative tracer) at the outlet of the tank‐scale HRX Well.     

Feasibility Study of In-Well Vapor Stripping Using Airlift Pumping

This study tests the feasibility of an aquifer remediation concept proposed by Gvirtzman and Gorelick (1992) that involves the removal of volatile organic compounds (VOCs) dissolved in ground water. The principal is 10 inject air into a well, creating airlift pumping, which is used as a means of in-well vapor stripping. The partially treated water is diverted away from the well and infiltrates back to the water table, thus allowing remediation of a larger aquifer volume.
A remediation well prototype, constructed in a laboratory aquifer model, was used to demonstrate the processes involved. The removal rates of trichloroethylene, toluene, and chloroform were monitored using eight triple-level observation wells. The continuous decrease of VOC concentrations during the short-term experiment has yielded macroscopic evidence that the process offers some promise. It was found that the flow field in the saturated zone. involving the continuous water circulation between the pumping well and the recharging area, caused temporal and spatial variation in remediation efficiency.     

An analytical approach to capture zone delineation for a well near a stream with a leaky layer

Abstract

An analytical solution is developed to delineate the capture zone of a pumping well in an aquifer with a regional flow perpendicular to a stream, assuming a leaky layer between the stream and the aquifer. Three different scenarios are considered for different pumping rates. At low pumping rates, the capture zone boundary will be completely contained in the aquifer. At medium pumping rates, the tip of the capture zone boundary will intrude into the leaky layer. Under these two scenarios, all the pumped water is supplied from the regional groundwater flow in the aquifer. At high pumping rates, however, the capture zone boundary intersects the stream and pumped water is supplied from both the aquifer and the stream. The two critical pumping rates which separate these three scenarios, as well as the proportion of pumped water from the stream and the aquifer, are determined for different hydraulic settings.

Editor D. Koutsoyiannis; Associate editor A. Koussis

Citation Asadi-Aghbolaghi, M., Rakhshandehroo, G.R., and Kompani-Zare, M., 2013. An analytical approach to capture zone delineation for a well near a stream with a leaky layer. Hydrological Sciences Journal, 58 (8), 1813–1823.     

Assessing radial transmissivity variation in heterogeneous aquifers using analytical techniques

Pumping tests are one of the most commonly used in situ testing techniques for assessing aquifer hydraulic properties. Numerous researches have been conducted to predict the effects of aquifer heterogeneity on the groundwater levels during pumping tests. The objectives of the present work were as follows: (1) to predict drawdown conditions and to estimate aquifer properties during pumping tests undertaken in radially symmetric heterogeneous aquifers, and (2) to identify a method for assessing the transmissivity field along the radial coordinate in radially symmetric and fully heterogeneous transmissivity fields. The first objective was achieved by expanding an existing analytical drawdown formulation that was valid for a radially symmetric confined aquifer with two concentric zones around the pumping well to an N concentric zone confined aquifer having a constant transmissivity value within each zone. The formulation was evaluated for aquifers with three and four concentric zones to assess the effects of the transmissivity field on the drawdown conditions. The specific conditions under which aquifer properties could be identified using traditional methods of analysis were also evaluated. The second objective was achieved by implementing the inverse solution algorithm (ISA), which was developed for petroleum reservoirs to groundwater aquifer settings. The results showed that the drawdown values are influenced by a volumetric integral of a weighting function and the transmissivity field within the cone of depression. The weighting function migrates in tandem with the expanding cone of depression. The ability of the ISA to predict radially symmetric and log‐normally distributed transmissivity fields was assessed against analytical and numerical benchmarks. The results of this investigation indicated that the ISA method is a viable technique for evaluating the radial transmissivity variations of heterogeneous aquifer settings. Copyright © 2013 John Wiley & Sons, Ltd.     

Modeling of Contaminant Movement Near Pumping Wells: Saturated-Unsaturated Flow with Particle Tracking

A transient axisymmetric saturated-unsaturated numerical flow model was coupled with a particle tracking model to investigate the movement of contaminants when a shallow unconfined aquifer is pumped at a constant rate. The particle tracking model keeps track of locations and masses of solutes in the aquifer, and the time of capture by the well. At the end of each time-step the flow model solves the Richard's equation for the hydraulic head distribution from which elemental velocities are calculated. Solutes are then displaced for a period equivalent to the time-step using both the magnitude and direction of the elemental velocities. Numerical experiments were performed to investigate effluent concentrations in wells with screens of different length and in different positions relative to zones of stratified contamination. At early times of pumping the effluent concentrations were similar to the concentrations adjacent to the well screen, but at late times, the concentrations approached the vertically averaged concentration in the aquifer. Time to attain the vertically averaged concentration was determined by the well geometry, initial location of the contaminant plume in relation to the well screen, and hydraulic properties of the aquifer. The results are consistent with the hydraulics of flow to a pumping well and of particular importance, they demonstrate that short-term pump tests could give erroneous design concentrations for pump-and-treat systems. The model provides a means of quantifying arrival times and mixing ratios. It could therefore provide a useful means of designing production wells in aquifers with stratified contamination and more efficient recovery systems for aquifer remediation.     

Hydrogeology of the Brunswick (Passaic) Formation and Implications for Ground Water Monitoring Practice

Fractured shales of the Brunswick Formation provide a major aquifer in the most industrialized region of New Jersey. Numerous cases of ground water contamination have been documented in this formation. However, effectiveness of monitoring and remediation efforts is often hampered by the use of inappropriate concepts regarding ground water flow controls in this complex aquifer system. One such concept presumes that near-vertical fractures parallel to the strike of beds provide principal passages for the flow and produce an anisotropic response to pumping stress. Field evidence presented in this paper confirms that the Brunswick Formation hosts a gently dipping, multiunit, leaky aquifer system that consists of thin water-bearing units and thick intervening aquitards. The water-bearing units are associated with major bedding partings and/or intensely fractured seams. Layered heterogeneity of such a dipping multiunit aquifer system produces an anisotropic flow pattern with preferential flow along the strike of beds. Within the weathered zone, the permeability of the water-bearing units can be greatly reduced. The commonly used hydrogeologic model of the Brunswick as a one-aquifer system, sometimes with vaguely defined "shallow" and "deep" zones, often leads to the development of inadvertent cross-flows within monitoring wells. If undetected, cross-flows may promote contaminant spread into deeper units and impair the quality of hydrogeologic data. Hydrogeologic characterization of the Brunswick shales at any given site should be aimed primarily at identification of the major water-bearing and aquitard units. Recommended techniques for this characterization include fluid logging and other in-well tests.     

Steady flow to a well near a stream with a leaky bed

We present an explicit analytic solution for steady, two-dimensional ground water flow to a well near a leaky streambed that penetrates the aquifer partially. Leakage from the stream is approximated as occurring along the centerline of the stream. The problem domain is infinite and pumping on one side of the stream induces flow on the other side. The solution includes the effects of uniform flow in the far field and a sloping hydraulic head in the stream. We use the solution to investigate the interaction between ground water and surface water in the stream, the effects of pumping on the opposite side of the stream, and the effects of the leaky streambed on the capture zone envelope of the well. We develop a relationship between parameters such that the pumping well will not capture water from the stream, or from the opposite side of the stream. When the discharge of the well is large enough to capture water from the stream, the shape of the capture zone envelope depends on flow conditions on the side of the stream opposite the well.     

A post audit of a model-designed ground water extraction system

Model post audits test the predictive capabilities of ground water models and shed light on their practical limitations. In the work presented here, ground water model predictions were used to design an extraction/treatment/injection system at a military ammunition facility and then were re-evaluated using site-specific water-level data collected approximately one year after system startup. The water-level data indicated that performance specifications for the design, i.e., containment, had been achieved over the required area, but that predicted water-level changes were greater than observed, particularly in the deeper zones of the aquifer. Probable model error was investigated by determining the changes that were required to obtain an improved match to observed water-level changes. This analysis suggests that the originally estimated hydraulic properties were in error by a factor of two to five. These errors may have resulted from attributing less importance to data from deeper zones of the aquifer and from applying pumping test results to a volume of material that was larger than the volume affected by the pumping test. To determine the importance of these errors to the predictions of interest, the models were used to simulate the capture zones resulting from the originally estimated and updated parameter values. The study suggests that, despite the model error, the ground water model contributed positively to the design of the remediation system.     

Modeling and Remediation of Ground Water Contamination at the Engelse Werk Wellfield

Contaminants have been threatening the Engelse Werk wellfield located between the town of Zwolle and the IJssel River in the Netherlands. Chemical analysis of water samples taken in production wells, both at the IJssel River and near the Zwolle railway station, indicated elevated concentrations of mainly organic contaminants including benzene, bentazon, acenaftene, trichloroethane, and bromacil. Immediate contaminant prevention and remediation measures are needed to safeguard the production wells. Ground water flow and transport models were developed to assist in the design of remediation strategies. Ground water flow models indicated that the IJssel River and a waste disposal ditch at the railway station are within the capture zone of the wellfield. A chloride transport model simulated minimum travel times in the order of four to 13 years for contaminants in the IJssel River to reach the production wells of the wellfield. A transport model for benzene was set up to advise on the remediation measures to be taken at the waste disposal ditch to clean up the contamination in the upper aquifer between this site and the Engelse Werk wellfield. The designed remediation system consists of 12 pumping wells with a combined capacity of 1650 m³/day. The system is capable of reducing the benzene levels at the threatened production wells at the Engelse Werk wellfield to a permissible level below 0.1 μg/L within a period of 5 years.     

The capture efficiency map: the capture zone under time-varying flow

The capture zone or contributing area of a ground water extraction well can be defined as that portion of the aquifer from which the well draws its water. Accurate delineation of capture zones is important in many ground water remediation applications and in the definition of wellhead protection areas. Their mathematical delineation is often simplified by using quasi-steady-state models based on time-weighted average pumping rates and background hydraulic gradients. We present a new semianalytic approach for the definition of capture zones under transient-flow conditions. We then use this approach to evaluate the effects of time variations in the direction of the background hydraulic gradient on capture. Results are presented in the form of capture efficiency maps (CEMs). Although the area contributing to a given well is found to generally expand relative to the steady-state average capture zone when the gradient direction varies, the zone of 100% capture may expand or contract depending on site-specific conditions. We illustrate our CEM approach by applying it to the design of a plume containment system.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.