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.     

Conservative identification of nodal transmissivities: theory,computational experiments and application to the Rharb coastal aquifer (Morocco)

Abstract

In order to calculate the transmissivity from the inverse problem corresponding to the groundwater flow in an isotropic horizontal aquifer, a numerical conservative approach is tested. The method deals with triangulation of the domain and applies the conservation of mass to elements of the mesh using the harmonic mean for internodal transmissivities. An optimal sweeping algorithm is used to evaluate nodal transmissivities from one element to another with a minimal relative error accumulation. The practical importance of the method is demonstrated through two synthetic examples representing those experienced in the field, then through application to a Moroccan aquifer. The computed hydraulic head is well fitted to the reference one, which confirms the validity of the identified transmissivity model.     

Hydraulic diffusivity in a coastal aquifer: spectral analysis of groundwater level in responses to marine system

The hydraulic diffusivity gives a measure of diffusion speed of pressure disturbances in groundwater system; large values of hydraulic diffusivity lead to fast propagation of signals in aquifer. This research provides a novel design and derives spectral representation to determine hydraulic diffusivity using spectral analysis of groundwater levels coupled with time-dependent boundary adjacent to marine system and no flow boundary in aquifer system. To validate the proposed method, water levels of fluctuated boundary and groundwater well in a sandy confined aquifer were collected. The hydraulic diffusivity is then obtained by an inverse process in the non-linear complex form of spectral relationship. The method essentially is constructed on the conceptual design of natural forcing transmitted in large aquifer. It is unlike the conventional field pumping test which is only used to determine hydraulic properties of groundwater in small range around the well. Hydraulic diffusivity of the confined aquifer is determined using real observation and then checked by comparing to the published range. It suggests that without local aquifer test to estimate hydraulic diffusivity in a coastal aquifer using spectral representation with its relevant flow system and boundary has become feasible.     

Aquifer Imaging with Oscillatory Hydraulic Tomography: Application at the Field Scale

Modeling and laboratory experiments have demonstrated the ability of oscillatory hydraulic tomography (OHT) to characterize heterogeneity in aquifer hydraulic properties. In OHT, a location is stressed via periodic pumping/injection at a set frequency, and the resulting head signal is measured at a number of monitoring locations. The source of oscillations is repeatedly moved, allowing tomographic imaging of aquifer properties. Changing the period of oscillation also results in observations with additional information. In theory, OHT is comparable to other hydraulic tomography methods in that distributed pressure change measurements provide characterization information. In practice, OHT has several benefits including: (1) little to no water injected into or extracted from the aquifer; and (2) an observational signal at a set period that can be easily extracted in the presence of noise. We report the first field application of OHT, carried out at the Boise Hydrogeophysical Research Site (BHRS) using an oscillating signal generator with a very small cycling volume of <2 L, and a period range of 5 to 70 s. For these tests, signals were detected at distances of over 15 m. After processing to extract periodic signal properties, we perform tomography using a frequency-domain numerical model for groundwater flow. In comparing results against prior characterization results from the BHRS, we find moderate to strong positive correlations between K profiles estimated via different methods at multiple wells, with moderate overall correlation between estimated three-dimensional (3D) K volumes.     

B.—UNESCO

Summary

Methods are outlined for determining the hydraulic properties of an aquifer, namely, the transmissivity and storativity from a well located near one or more impermeable barriers. The procedures are based on the theory of the nonsteady radial flow toward a well steadily discharging from an aquifer of semiinfinite areal extent. Two different flow systems have been considered, namely a well near an impermeable barrier and a well located between two impermeable barriers.     

Determination of confined phreatic aquifer anisotropy

The determination of anisotropy in a phreatic aquifer by means of a pumping test can be carried out, under given conditions, with common analysis methods. This paper develops and applies an amended method based on a technique proposed by Boulton (1970). This new method enables a quick, simple determination of vertical hydraulic conductivity, given drawdown data from a pumping well and at least one observation well, that is applicable to confined phreatic aquifers.     

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

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

Modeling of Plume Capture by Continuous, Low-Permeability Barriers

Modeling was performed to simulate ground water flow through reactive barriers of lower hydraulic conductivity than the surrounding aquifer to determine the plume capture widths. As a plume approaches such a barrier, it spreads laterally. Therefore, to intercept an entire plume, the barrier must be wider than the upgradient width of the undisturbed plume. The results indicate that, for practical values of barrier thickness and plume width, hydraulic conductivities ten-fold less than that of the aquifer can be accommodated by making the width of the barrier approximately 20% greater than the upgradient width of the plume. Barrier hydraulic conductivities one-hundred-fold less than that of the aquifer may require barrier widths up to twice the width of the upgradient plume for plumes 100 feet wide (33 m) and as little as 1.1 times for plumes 1000 feet wide (325 m). The results presented here lend support to the view that novel emplacement methods that create zones of slightly lower hydraulic conductivity than the native aquifer may be viable alternatives to the excavation-and-backfill approaches which have thus far been used for installing permeable reactive barriers.     

Estimating hydraulic parameters when poroelastic effects are significant

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

Flow‐Through Stream Modeling with MODFLOW and MT3D: Certainties and Limitations

This paper aims to assess MODFLOW and MT3D capabilities for simulating the spread of contaminants from a river exhibiting an unusual relationship with an alluvial aquifer, with the groundwater head higher than the river head on one side and lower on the other (flow‐through stream). A series of simulation tests is conducted using a simple hypothetical model so as to characterize and quantify these limitations. Simulation results show that the expected contaminant spread could be achieved with a specific configuration composed of two sets of parameters: (1) modeled object parameters (hydraulic groundwater gradient, hydraulic conductivity values of aquifer and streambed), and (2) modeling parameters (vertical discretization of aquifer, horizontal refinement of stream modeled with River [RIV] package). The influence of these various parameters on simulation results is investigated, and potential complications and errors are identified. Contaminant spread from stream to aquifer is not always reproduced by MT3D due to the RIV package's inability to simulate lateral exchange fluxes between stream and aquifer. This paper identifies the need for a MODFLOW streamflow package allowing lateral stream‐aquifer interactions and streamflow routine calculations. Such developments could be of particular interest for modeling contaminated flow‐through streams.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

一种测求水井含水层导水系数的新方法

本文给出了一种测求水井含水层导水系数的新方法。利用Cooper理论和振动理论,通过简单的试验,可以测求出水井含水层的导水系数。用这种方法计算出珍珠泉井含水层的导水系数为2439m~2/d,用抽水试验法测得该系数为2618m~2/d,两者符合得较好。     

Implementing a Method of Screening One-Well Hydraulic Barrier Design Alternatives

This article provides details of applying the method developed by the authors ( Rubin et al. 2008b ) for screening one-well hydraulic barrier design alternatives. The present article with its supporting information (manual and electronic spreadsheets with a case history example) provides the reader complete details and examples of solving the set of nonlinear equations developed by Rubin et al. (2008b) . It allows proper use of the analytical solutions and also depicting the various charts given by Rubin et al. (2008b) . The final outputs of the calculations are the required position and the discharge of the pumping well. If the contaminant source is nonaqueous phase liquid (NAPL) entrapped within the aquifer, then the method provides an estimate of the aquifer remediation progress (which is a by-product) due to operating the hydraulic barrier.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.     

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

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

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

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

Steady flow to a partially penetrating blind‐wall well in a confined aquifer

Wells in aquifers of loose collapsible sediment are cased so that they have a blind wall and gain water only from the bottom. The hydraulic gradient established at the bottom of these wells during pumping brings the aquifer materials in a quicksand state, which may cause abrasion of pipes and pumps and even the destruction of well structure. To examine the quicksand occurrence, an analytical solution for the steady flow to a partially penetrating blind‐wall well in a confined aquifer is developed. The validity of the proposed solution is evaluated numerically. The sensitivity of maximum vertical gradient along the well bottom in response to aquifer and well parameters is examined. The solution is presented in the form of dimensionless‐type curves and equations that can be easily used to design the safe pumping rate and optimum well geometry to protect the well against sand production. The solution incorporates the anisotropy of aquifer materials and can also be used to determine the hydraulic conductivity of the aquifer. Copyright © 2012 John Wiley & Sons, Ltd.     

Tensor Hydraulic Conductivity Estimation for Heterogeneous Aquifers under Unknown Boundary Conditions

A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements.     

Application of a Bayesian approach to stochastic delineation of capture zones

This paper presents a Bayesian Monte Carlo method for evaluating the uncertainty in the delineation of well capture zones and its application to a wellfield in a heterogeneous, multiaquifer system. In the method presented, Bayes' rule is used to update prior distributions for the unknown parameters of the stochastic model for the hydraulic conductivity, and to calculate probability-based weights for parameter realizations using head residuals. These weights are then assigned to the corresponding capture zones obtained using forward particle tracking. Statistical analysis of the set of weighted protection zones results in a probability distribution for the capture zones. The suitability of the Bayesian stochastic method for a multilayered system is investigated, using the wellfield Het Rot at Nieuwrode, Belgium, located in a three-layered aquifer system, as an example. The hydraulic conductivity of the production aquifer is modeled as a spatially correlated random function with uncertain parameters. The aquitard and overlying unconfined aquifer are assigned random, homogeneous conductivities. The stochastic results are compared with deterministic capture zones obtained with a calibrated model for the area. The predictions of the stochastic approach are more conservative and indicate that parameter uncertainty should be taken into account in the delineation of well capture zones.