Application of a nonlinear slug test model

Knowledge of the hydraulic conductivity distribution is of utmost importance in understanding the dynamics of an aquifer and in planning the consequences of any action taken upon that aquifer. Slug tests have been used extensively to measure hydraulic conductivity in the last 50 years since Hvorslev's (1951) work. A general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the wellbore, and a Hvorslev model for the aquifer has been implemented in this work. The nonlinear model has three parameters: beta, which is related primarily to radius changes in the water column; A, which is related to the nonlinear head losses; and K, the hydraulic conductivity. An additional parameter has been added representing the initial velocity of the water column at slug initiation and is incorporated into an analytical solution to generate the first time step before a sequential numerical solution generates the remainder of the time solution. Corrections are made to the model output for acceleration before it is compared to the experimental data. Sensitivity analysis and least squares fitting are used to estimate the aquifer parameters and produce some diagnostic results, which indicate the accuracy of the fit. Finally, an example of field data has been presented to illustrate the application of the model to data sets that exhibit nonlinear behavior. Multiple slug tests should be taken at a given location to test for nonlinear effects and to determine repeatability.     

Influence of ancient thrust faults on the hydrogeology of the Blue Ridge Province

The Blue Ridge Province contains ubiquitous northeast-southwest-trending thrust faults or smaller thrust "slivers" that greatly impact the nature and character of ground water flow in this region. Detailed investigations at a field site in Floyd County, Virginia, indicate that high-permeability zones occur in the brittle crystalline rocks above these thrust faults. Surface and borehole geophysics, aquifer tests, and chlorofluorocarbon and geochemical data reveal that the shallow saprolite aquifer is separated from the deeper fault-zone aquifer by a low-fracture permeability bedrock confining unit, the hydraulic conductivity of which has been estimated to be six orders of magnitude less than the conductivity of the fault zones at the test site. Within the Blue Ridge Province, these fault zones can occur at depths of 300 m or more, can contain a significant amount of storage, and yield significant quantities of water to wells. Furthermore, it is expected that these faults may compartmentalize the deep aquifer system. Recharge to and discharge from the deep aquifer occurs by slow leakage through the confining unit or through localized breach zones that occur where quartz accumulated in high concentrations during metamorphism and later became extensively fractured during episodes of deformation. The results of this investigation stress the importance of thrust faults in this region and suggest that hydrogeologic models for the Blue Ridge Province include these ancient structural features. Faults in crystalline-rock environments may not only influence the hydrology, they may dominate the flow characteristics of a region.     

A simple approach using Bouwer and Rice's method for slug test data analysis

Slug test data obtained from tests performed in an unconfined aquifer are commonly analyzed by graphical or numerical approaches to determine the aquifer parameters. This paper derives three fourth-degree polynomials to represent the relationship between Bouwer and Rice's coefficients and the ratio of the screen length to the radius of the gravel envelope. A numerical approach using the nonlinear least squares and Newton's method is used to determine hydraulic conductivity from the best fit of the slug test data. The method of nonlinear least squares minimizes the sum of the squares of the differences between the predicted and observed water levels inside the well. With the polynomials, the hydraulic conductivity can be obtained by simply solving the nonlinear least squares equation by Newton's method. A computer code, SLUGBR, was developed from the derived polynomials using the proposed numerical approach. The results of analyzing two slug test datasets show that SLUGBR can determine hydraulic conductivity with very good accuracy.     

Delineating alluvial aquifer heterogeneity using resistivity and GPR data

Conceptual geological models based on geophysical data can elucidate aquifer architecture and heterogeneity at meter and smaller scales, which can lead to better predictions of preferential flow pathways. The macrodispersion experiment (MADE) site, with >2000 measurements of hydraulic conductivity obtained and three tracer tests conducted, serves as an ideal natural laboratory for examining relationships between subsurface flow characteristics and geophysical attributes in fluvial aquifers. The spatial variation of hydraulic conductivity measurements indicates a large degree of site heterogeneity. To evaluate the usefulness of geophysical methods for better delineating fluvial aquifer heterogeneities and distribution of preferential flow paths, a surface grid of two-dimensional ground penetrating radar (GPR) and direct current (DC) resistivity data were collected. A geological model was developed from these data that delineate four stratigraphic units with distinct electrical and radar properties including (from top to bottom) (1) a meandering fluvial system (MFS); (2) a braided fluvial system (BFS); (3) fine-grained sands; and (4) a clay-rich interval. A paleochannel, inferred by other authors to affect flow, was mapped in the MFS with both DC resistivity and GPR data. The channel is 2 to 4 m deep and, based on resistivity values, is predominantly filled with clay and silt. Comparing previously collected hydraulic conductivity measurements and tracer-plume migration patterns to the geological model indicates that flow primarily occurs in the BFS and that the channel mapped in the MFS has no influence on plume migration patterns.     

Delineation of regional arid karstic aquifers: an integrative data approach

This research integrates data procedures for the delineation of regional ground water flow systems in arid karstic basins with sparse hydrogeologic data using surface topography data, geologic mapping, permeability data, chloride concentrations of ground water and precipitation, and measured discharge data. This integrative data analysis framework can be applied to evaluate arid karstic aquifer systems globally. The accurate delineation of ground water recharge areas in developing aquifer systems with sparse hydrogeologic data is essential for their effective long-term development and management. We illustrate the use of this approach in the Cuatrociénegas Basin (CCB) of Mexico. Aquifers are characterized using geographic information systems for ground water catchment delineation, an analytical model for interbasin flow evaluation, a chloride balance approach for recharge estimation, and a water budget for mapping contributing catchments over a large region. The test study area includes the CCB of Coahuila, Mexico, a UNESCO World Biosphere Reserve containing more than 500 springs that support ground water-dependent ecosystems with more than 70 endemic organisms and irrigated agriculture. We define recharge areas that contribute local and regional ground water discharge to springs and the regional flow system. Results show that the regional aquifer system follows a topographic gradient that during past pluvial periods may have linked the Río Nazas and the Río Aguanaval of the Sierra Madre Occidental to the Río Grande via the CCB and other large, currently dry, upgradient lakes.     

Fracture control of ground water flow and water chemistry in a rock aquitard

There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/S(s)) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies.     

Numerical Modeling of Slug Tests with MODFLOW Using Equivalent Well Blocks

Slug tests are a widely used technique to estimate aquifer hydraulic parameters and the test data are generally interpreted with analytical solutions under various assumptions. However, these solutions are not convenient when slug tests are required to be analyzed in a three‐dimensional model for complex aquifer‐aquitard systems. In this study, equivalent well blocks (EWB) are proposed in numerical modeling of slug test data with MODFLOW. Multi‐well slug tests in partially penetrating wells with skin zones can be simulated. Accuracy of the numerical method is demonstrated by benchmarking with analytical solutions. The EWB method is applied in a case study on slug tests in aquitards in the Pearl River Delta, China.     

Application of Electromagnetic Logging to Contamination Investigations in Glacial Sand-and-Gravel Aquifers

Electromagnetic (EM) logging provides an efficient method for high-resolution, vertical delineation of electrically conductive contamination in glacial sand-and-gravel aquifers. LM. gamma, and lithologic logs and specific conductance data from sand-and-gravel aquifers at five sites in the northeastern United States were analyzed to define the relation of KM conductivity to aquifer lithology and water quality. Municipal waste disposal, septic waste discharge, or highway deicing salt application at these sites has caused contaminant plumes in which the dissolved solids concentration and specific conductance of ground water exceed background levels by as much as 10 to 20 limes.
The major hydrogeologic factors that affected KM log response at the five sites were the dissolved solids concentration of the ground water and the silt and clay content in the aquifer. KM conductivity of sand and gravel with uncontaminated water ranged from less than 5 to about 10 millisiemens per meter (mS/m); that of silt and clay zones ranged from about 15 to 45 mS/m: and that of the more highly contaminated zones in sand and gravel ranged from about 10 to more than 80 mS/m. Specific conductance of water samples from screened intervals in sand and gravel at selected monitoring well installations was significantly correlated with KM conductivity.
CM logging can be used in glacial sand-and-gravel aquifer investigations to (1) determine optimum depths for the placement of monitoring well screens: (2) provide a nearly continuous vertical profile of specific conductance to complement depth-specific water quality samples; and (3) identify temporal changes in water quality through sequential logging. Detailed lithologic or gamma logs, preferably both, need to be collected along with the F.M logs to define zones in which elevated EM conductivity is caused by the presence of sill and clay beds rather than contamination.     

Characterization of High-K Pathways by Borehole Flowmeter and Tracer Tests

This paper reviews different borehole flowmeter analysis methods and evaluates their applicability to a test site composed of fluvial deposits. Results from tracer and aquifer tests indicate that the aquifer is highly heterogeneous and that low-K skin effects exist at the wells. Borehole flowmeter tests were performed at 37 wells. An appropriate method for calculating borehole flowmeter K values was developed based on results from multiwell pumping tests, single-well pumping tests, and slug tests. The flowmeter data produced 881 K values. The trends and the magnitude of the K values are consistent with results from geologic investigations, recirculating tracer tests, and large-scale multiwell pumping tests. The field tests illustrate that high-K deposits can significantly affect ground-water flows in some heterogeneous fluvial aquifers.     

A simple correction for slug tests in small-diameter wells

A simple procedure is presented for correcting hydraulic conductivity (K) estimates obtained from slug tests performed in small-diameter installations screened in highly permeable aquifers. Previously reported discrepancies between results from slug tests in small-diameter installations and those from tests in nearby larger-diameter wells are primarily a product of frictional losses within the small-diameter pipe. These frictional losses are readily incorporated into existing models for slug tests in high-K aquifers, which then serve as the basis of a straightforward procedure for correcting previously obtained K estimates. A demonstration of the proposed procedure using data from a series of slug tests performed in a controlled field setting confirms the validity of the approach. The results of this demonstration also reveal the detailed view of spatial variations in K that can be obtained using slug tests in small-diameter installations.     

On using simple time-of-travel capture zone delineation methods

As part of its Wellhead Protection Program, the U.S. EPA mandates the delineation of "time-of-travel capture zones" as the basis for the definition of wellhead protection zones surrounding drinking water production wells. Depending on circumstances the capture zones may be determined using methods that range from simply drawing a circle around the well to sophisticated ground water flow and transport modeling. The simpler methods are attractive when faced with the delineation of hundreds or thousands of capture zones for small public drinking water supply wells. On the other hand, a circular capture zone may not be adequate in the presence of an ambient ground water flow regime. A dimensionless time-of-travel parameter T is used to determine when calculated fixed-radius capture zones can be used for drinking water production wells. The parameter incorporates aquifer properties, the magnitude of the ambient ground water flow field, and the travel time criterion for the time-of-travel capture zone. In the absence of interfering flow features, three different simple capture zones can be used depending on the value of T . A modified calculated fixed-radius capture zone proves protective when T < 0.1, while a more elongated capture zone must be used when T > 1. For values of T between 0.1 and 1, a circular capture zone can be used that is eccentric with respect to the well. Finally, calculating T allows for a quick assessment of the validity of circular capture zones without redoing the delineation with a computer model.     

Natural Gradient Tracer Test to Evaluate Natural Attenuation of MTBE Under Anaerobic Conditions

A natural gradient tracer test using perdeuterated MTBE was conducted in an anaerobic aquifer to determine the relative importance of dispersion and degradation in reducing MTBE concentrations in ground water. Preliminary ground water chemistry and hydraulic conductivity data were used to place the tracer within an existing dissolved MTBE plume at Port Hueneme, California. Following one year of transport, the tracer plume was characterized in detail.
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m²/day and 0.08 m²/day, respectively). A mass-balance analysis, carried out after one year of transport, accounted for 110% of the injected mass and indicated that no significant mass loss occurred. The plume structure created by zones of higher and lower hydraulic conductivity at the site was complex, consisting of several localized areas of high tracer concentration in a lower concentration plume. This is important because the aquifer has generally been characterized as exhibiting fairly minor heterogeneity. In addition, the tracer plume followed a curved flowpath that deviated from the more macroscopic direction of ground water flow inferred from local ground water elevation measurements and the behavior of the existing plume. Understanding the mass balance, plume structure, curvature of the tracer plume, and consequently natural attenuation behavior required the detailed sampling approach employed in this study. These data imply that a detailed understanding of site hydrogeology and an extensive sampling network may be critical for the correct interpretation of monitored natural attenuation of MTBE.     

Multiscale characterization of a heterogeneous aquifer using an ASR operation

Heterogeneity in the physical properties of an aquifer can significantly affect the viability of aquifer storage and recovery (ASR) by reducing the recoverable proportion of low-salinity water where the ambient ground water is brackish or saline. This study investigated the relationship between knowledge of heterogeneity and predictions of solute transport and recovery efficiency by combining permeability and ASR-based tracer testing with modeling. Multiscale permeability testing of a sandy limestone aquifer at an ASR trial site showed that small-scale core data give lower-bound estimates of aquifer hydraulic conductivity (K), intermediate-scale downhole flowmeter data offer valuable information on variations in K with depth, and large-scale pumping test data provide an integrated measure of the effective K that is useful to constrain ground water models. Chloride breakthrough and thermal profiling data measured during two cycles of ASR showed that the movement of injected water is predominantly within two stratigraphic layers identified from the flowmeter data. The behavior of the injectant was reasonably well simulated with a four-layer numerical model that required minimal calibration. Verification in the second cycle achieved acceptable results given the model's simplicity. Without accounting for the aquifer's layered structure, high precision could be achieved on either piezometer breakthrough or recovered water quality, but not both. This study demonstrates the merit of an integrated approach to characterizing aquifers targeted for ASR.     

Evaluation of Combined Direct-Push Methods Used for Aquifer Model Generation

Most established methods to characterize aquifer structure and hydraulic conductivities of hydrostratigraphical units are not capable of delivering sufficient information in the spatial resolution that is desired for sophisticated numerical contaminant transport modeling and adapted remediation design. With hydraulic investigation methods based on the direct-push (DP) technology such as DP slug tests, DP injection logging, and the hydraulic profiling tool, it is possible to rapidly delineate hydrogeological structures and estimate their hydraulic conductivity in shallow unconsolidated aquifers without the need for wells. A combined application of these tools was used for the investigation of a contaminated German refinery site and for the setup of hydraulic aquifer models. The quality of DP investigation and the models was evaluated by comparisons of tracer transport simulations using these models and measured breakthroughs of two natural gradient tracer tests. Model scenarios considering the information of all tools together showed good reproduction of the measured breakthroughs, indicating the suitability of the approach and a minor impact of potential technical limitations. Using the DP slug tests alone yielded significantly higher deviations for the determined hydraulic conductivities compared to considering two or three of the tools. Realistic aquifer models developed on basis of such combined DP investigation approaches can help optimize remediation concepts or identify flow regimes for aquifers with a complex structure.     

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.     

Semi‐analytical solution for a slug test in partially penetrating wells including the effect of finite‐thickness skin

This paper presents a new semi‐analytical solution for a slug test in a well partially penetrating a confined aquifer, accounting for the skin effect. This solution is developed based on the solution for a constant‐flux pumping test and a formula given by Peres and co‐workers in 1989. The solution agrees with that of Cooper and co‐workers and the KGS model when the well is fully penetrating. The present solution can be applied to simulate the temporal and spatial head distributions in both the skin and formation zones. It can also be used to demonstrate the influences of skin type or skin thickness on the well water level and to estimate the hydraulic parameters of the skin and formation zones using a least‐squares approach. The results of this study indicate that the determination of hydraulic conductivity using a conventional slug‐test data analysis that neglects the presence of a skin zone will give an incorrect result if the aquifer has a skin zone. Copyright © 2008 John Wiley & Sons, Ltd.     

Effects of water use on arsenic release to well water in a confined aquifer

Field-based experiments were designed to investigate the release of naturally occurring, low to moderate (< 50 microg/L) arsenic concentrations to well water in a confined sandstone aquifer in northeastern Wisconsin. Geologic, geochemical, and hydrogeologic data collected from a 115 m2 site demonstrate that arsenic concentrations in ground water are heterogeneous at the scale of the field site, and that the distribution of arsenic in ground water correlates to solid-phase arsenic in aquifer materials. Arsenic concentrations in a test well varied from 1.8 to 22 microg/L during experiments conducted under no, low, and high pumping rates. The quality of ground water consumed from wells under typical domestic water use patterns differs from that of ground water in the aquifer because of reactions that occur within the well. Redox conditions in the well can change rapidly in response to ground water withdrawals. The well borehole is an environment conducive to microbiological growth, and biogeochemical reactions also affect borehole chemistry. While oxidation of sulfide minerals appears to release arsenic to ground water in zones within the aquifer, reduction of arsenic-bearing iron (hydr)oxides is a likely mechanism of arsenic release to water having a long residence time in the well borehole.     

An Evaluation of Interpolation Methodologies for Generating Three-Dimensional Hydraulic Property Distributions from Measured Data

The process of attempting to model ground-water systems requires a good understanding of the spatial variation of aquifer hydraulic properties. The capabilities of the more recent innovative flowmeters such as the electromagnetic and heat pulse flowmeters provide the sensitivity to measure ambient flows and pump-induced flows. These flowmeters provide the measurements of pump-induced vertical flows which are analyzed to obtain vertical variations in horizontal hydraulic conductivity, K(z). With discrete areal K-values, K(x, y), and vertical profiles of K, provided by multiwell testing, the essential elements are present to produce a three-dimensional hydraulic conductivity field. The advent of these new flow measuring devices has contributed much to the motivation behind this paper. This paper presents the results of applying deterministic and stochastic methodology to the three-dimensional interpolation of hydraulic properties, specifically, hydraulic conductivity, K. Three of the approaches applied in this paper are deterministic in nature, inverse-distance weighting, inverse-distance-squared weighting, and ordinary kriging, while the fourth is a stochastic approach based on self-affine fractals. All of the methods are applied to measured data collected from 14 wells at a site in the United States near Mobile, Alabama. The three-dimensional K-distributions generated by each of the methods are used as inputs to an advective based transport model with the resulting model output compared to a two-well tracer study run previously at the same site.     

Site characterization in densely fractured dolomite: comparison of methods

One of the challenges in characterizing fractured-rock aquifers is determining whether the equivalent porous medium approximation is valid at the problem scale. Detailed hydrogeologic characterization completed at a small study site in a densely fractured dolomite has yielded an extensive data set that was used to evaluate the utility of the continuum and discrete-fracture approaches to aquifer characterization. There are two near-vertical sets of fractures at the site; near-horizontal bedding-plane partings constitute a third fracture set. Eighteen boreholes, including five coreholes, were drilled to a depth of approximately 10.6 m. Borehole geophysical logs revealed several laterally extensive horizontal fractures and dissolution zones. Flowmeter and short-interval packer testing identified which of these features were hydraulically important. A monitoring system, consisting of short-interval piezometers and multilevel samplers, was designed to monitor four horizontal fractures and two dissolution zones. The resulting network consisted of >70 sampling points and allowed detailed monitoring of head distributions in three dimensions. Comparison of distributions of hydraulic head and hydraulic conductivity determined by these two approaches suggests that even in a densely fractured-carbonate aquifer, a characterization approach using traditional long-interval monitoring wells is inadequate to characterize ground water movement for the purposes of regulatory monitoring or site remediation. In addition, traditional multiwell pumping tests yield an average or bulk hydraulic conductivity that is not adequate for predicting rapid ground water travel times through the fracture network, and the pumping test response does not appear to be an adequate tool for assessing whether the porous medium approximation is valid.     

Slug test in confined aquifers, the over-damped case: quasi-steady flow analysis

In the great majority of slug tests performed in wells fully penetrating confined geologic formations, and for over-damped conditions, the response data are evaluated with the transient-flow model of Cooper et al. (1967) when the radial hydraulic conductivity K(r) and the coefficient of specific storage S(s) are to be estimated. That particular analytical solution, however, is computationally involved and awkward to use. Thus, groundwater professionals often use a few pre-prepared type-curves to fit the data by a rough matching procedure, visually or computationally. On the other hand, the method of Hvorslev (1951), which assumes the flow to be quasi-steady, is much simpler but yields only K(r) estimates. In this work, we develop a complete quasi-steady flow model that includes a storage balance inside the aquifer and allows estimating K(r) and S(s) simultaneously, through matching of the well response data to a type-curve. The new model approximates the model of Cooper et al. closely and has the practical advantage that its solution type-curves are generated easily using an electronic spreadsheet, so that the optimal fit of data by a type-curve can be readily automated.