Slug Tests in Wells Screened Across the Water Table: Some Additional Considerations

The majority of slug tests done at sites of shallow groundwater contamination are performed in wells screened across the water table and are affected by mechanisms beyond those considered in the standard slug‐test models. These additional mechanisms give rise to a number of practical issues that are yet to be fully resolved; four of these are addressed here. The wells in which slug tests are performed were rarely installed for that purpose, so the well design can result in problematic (small signal to noise ratio) test data. The suitability of a particular well design should thus always be assessed prior to field testing. In slug tests of short duration, it can be difficult to identify which portion of the test represents filter‐pack drainage and which represents formation response; application of a mass balance can help confirm that test phases have been correctly identified. A key parameter required for all slug test models is the casing radius. However, in this setting, the effective casing radius (borehole radius corrected for filter‐pack porosity), not the nominal well radius, is required; this effective radius is best estimated directly from test data. Finally, although conventional slug‐test models do not consider filter‐pack drainage, these models will yield reasonable hydraulic conductivity estimates when applied to the formation‐response phase of a test from an appropriately developed well.     

Assessing the Hydraulic Conductivity Anisotropy and Skin-Effect Based on Data of Slug Tests in Partially Penetrating Wells

Water Resources - A semianalytic and approximate analytical solution is given to the problem of slug test in a partially penetrating well in a confined or unconfined anisotropic aquifer. An...     

Hydraulic Tomography: Continuity and Discontinuity of High‐K and Low‐K Zones

Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model.     

Numerical Modeling of the Near‐Field Hydraulics of Water Wells

Numerical flow models can be a useful tool for dimensioning water wells and to investigate the hydraulics in their near‐field. Fully laminar flow can be assumed for all models calculated up to the screen. Therefore models can be used to predict—at least qualitatively, neglecting turbulent losses inside the well—the spatial distribution of inflow into the well and the overall hydraulic performance of different combinations of aquifer parameters and technical installations. Models for both horizontal (plan view) and vertical flow (cross section) to wells were calculated for a variety of setups. For the latter, this included variations of hydraulic conductivity of the screen, pump position, and aquifer heterogeneity. Models of suction flow control devices showed that they indeed can homogenize inflow, albeit at the cost of elevated entrance losses.     

A Simple Field Method for Assessing Near‐Surface Saturated Hydraulic Conductivity

This technical note describes an effective and inexpensive field technique for measuring the saturated hydraulic conductivity of both undisturbed cores and repacked soil samples. The method requires no specialized equipment; everything that is required can be obtained in a hardware store. The method is a straightforward field implementation of the widely used falling‐head laboratory analysis directly derived from Darcy's law. As such, it sidesteps the need for empirical assumptions about soil texture and the relationship between saturated and unsaturated flow components which many permeameter‐based methods rely upon. The method is shown to produce results that are consistent with K values obtained elsewhere in the same homogeneous sand formation. Furthermore, the proposed method is useful for measuring hydraulic conductivity in drill cuttings obtained from direct push or auguring drill techniques, which cannot be done with any other field method. The range of hydraulic conductivity values that this test is appropriate for is on the order of 1E ? 7 m/s to 1E ? 3 m/s.     

Understanding the Geometry of Connected Fracture Flow with Multiperiod Oscillatory Hydraulic Tests

An understanding of the spatial and hydraulic properties of fast preferential flow pathways in the subsurface is necessary in applications ranging from contaminant fate and transport modeling to design of energy extraction systems. One method for the characterization of fracture properties over interwellbore scales is Multiperiod Oscillatory Hydraulic (MOH) testing, in which the aquifer response to oscillatory pressure stimulations is observed. MOH tests were conducted on isolated intervals of wells in siliciclastic and carbonate aquifers in southern Wisconsin. The goal was to characterize the spatial properties of discrete fractures over interwellbore scales. MOH tests were conducted on two discrete fractured intervals intersecting two boreholes at one field site, and a nest of three piezometers at another field site. Fracture diffusivity estimates were obtained using analytical solutions that relate diffusivity to observed phase lag and amplitude decay. In addition, MOH tests were used to investigate the spatial extent of flow using different conceptual models of fracture geometry. Results indicated that fracture geometry at both field sites can be approximated by permeable two‐dimensional fracture planes, oriented near‐horizontally at one site, and near‐vertically at the other. The technique used on MOH field data to characterize fracture geometry shows promise in revealing fracture network characteristics important to groundwater flow and transport.     

Determination of Hydraulic Conductivity from Grain‐Size Distribution for Different Depositional Environments

Over 400 unlithified sediment samples were collected from four different depositional environments in global locations and the grain‐size distribution, porosity, and hydraulic conductivity were measured using standard methods. The measured hydraulic conductivity values were then compared to values calculated using 20 different empirical equations (e.g., Hazen, Carman‐Kozeny) commonly used to estimate hydraulic conductivity from grain‐size distribution. It was found that most of the hydraulic conductivity values estimated from the empirical equations correlated very poorly to the measured hydraulic conductivity values with errors ranging to over 500%. To improve the empirical estimation methodology, the samples were grouped by depositional environment and subdivided into subgroups based on lithology and mud percentage. The empirical methods were then analyzed to assess which methods best estimated the measured values. Modifications of the empirical equations, including changes to special coefficients and addition of offsets, were made to produce modified equations that considerably improve the hydraulic conductivity estimates from grain size data for beach, dune, offshore marine, and river sediments. Estimated hydraulic conductivity errors were reduced to 6 to 7.1 m/day for the beach subgroups, 3.4 to 7.1 m/day for dune subgroups, and 2.2 to 11 m/day for offshore sediments subgroups. Improvements were made for river environments, but still produced high errors between 13 and 23 m/day.     

MODFLOW‐Based Coupled Surface Water Routing and Groundwater‐Flow Simulation

In this paper, we present a flexible approach for simulating one‐ and two‐dimensional routing of surface water using a numerical surface water routing (SWR) code implicitly coupled to the groundwater‐flow process in MODFLOW. Surface water routing in SWR can be simulated using a diffusive‐wave approximation of the Saint‐Venant equations and/or a simplified level‐pool approach. SWR can account for surface water flow controlled by backwater conditions caused by small water‐surface gradients or surface water control structures. A number of typical surface water control structures, such as culverts, weirs, and gates, can be represented, and it is possible to implement operational rules to manage surface water stages and streamflow. The nonlinear system of surface water flow equations formulated in SWR is solved by using Newton methods and direct or iterative solvers. SWR was tested by simulating the (1) Lal axisymmetric overland flow, (2) V‐catchment, and (3) modified Pinder‐Sauer problems. Simulated results for these problems compare well with other published results and indicate that SWR provides accurate results for surface water‐only and coupled surface water/groundwater problems. Results for an application of SWR and MODFLOW to the Snapper Creek area of Miami‐Dade County, Florida, USA are also presented and demonstrate the value of coupled surface water and groundwater simulation in managed, low‐relief coastal settings.     

Small Purge Method to Sample Vapor from Groundwater Monitoring Wells Screened Across the Water Table

Groundwater monitoring wells are present at most hydrocarbon release sites that are being assessed for cleanup. If screened across the vadose zone, these wells provide an opportunity to collect vapor samples that can be used in the evaluation of vapor movement and biodegradation processes occurring at such sites. This paper presents a low purge volume method (modified after that developed by the U.S. EPA) for sampling vapor from monitoring wells that is easy to implement and can provide an assessment of the soil gas total petroleum hydrocarbon (TPH) and O₂ concentrations at the base of the vadose zone. As a result, the small purge method allows for sampling of vapor from monitoring wells to support petroleum vapor intrusion (PVI) risk assessment. The small purge volume method was field tested at the Hal's service station site in Green River, Utah. This site is well‐known for numerous soil gas measurements containing high O₂ and high TPH vapor concentrations in the same samples which is inconsistent with well‐accepted biodegradation models for the vapor pathway. Using the low purge volume method, monitoring wells were sampled over, upgradient, and downgradient of the light nonaqueous phase liquid (LNAPL) footprint. Results from our testing at Hal's show that vapor from monitoring wells over LNAPL contained very low O₂ and high TPH concentrations. In contrast, vapor from monitoring wells not over LNAPL contained high O₂ and low TPH concentrations. The results of this study show that a low purge volume method is consistent with biodegradation models especially for sampling at sites where low permeability soils exist in and around a LNAPL source zone.     

Design for a Packer/Vacuum Slug Test System for Estimating Hydraulic Conductivity in Wells with LNAPLs, Casing Leaks, or Water Tables Intersecting the Screens

"Results indicate that the packer/vacuum system is an adequate method for obtaining values of hydraulic conductivity."     

An Integrated Laboratory Method to Measure and Verify Directional Hydraulic Conductivity in Fine‐to‐Medium Sandy Sediments

The constant‐head permeameter test (CHPT) is widely used in sandy samples as a standard method in the laboratory to investigate hydraulic conductivity (K). However, it neither can be used to consistently determine directional hydraulic conductivity (DHC) nor guarantee the comparability of measured K values of samples with different sizes. Therefore, this paper proposes an integrated laboratory method, called modified CHPT (MCHPT), for the efficient determination and verification of consistent DHC values in fine‐to‐medium sandy sediments, based on a new methodological framework. A precise and standardized procedure for preparing the experimental setup of MCHPT was conducted, based on the integrated experimental setup of CHPT and tracer tests. Moreover, a formula was yielded for the time‐optimized sample saturation control. In comparison with grain size‐based methods, the validity of consistent K_h and K_v values determined by MCHPT was convincing.     

Evaluating the Use of In‐Well Heat Tracer Tests to Measure Borehole Flow Rates

Recent research has demonstrated the use of in‐well heat tracer tests monitored by a fiber optic distributed temperature sensing (DTS) system to characterize borehole flow conditions in open bedrock boreholes. However, the accuracy of borehole flow rates determined from in‐well heat tracer tests has not been evaluated. The purpose of the research presented here is to determine whether borehole flow rates obtained using DTS‐monitored in‐well heat tracer tests are reasonable, and to evaluate the range of flow rates measureable with this method. To accomplish this, borehole flow rates measured using in‐well heat tracer tests are compared to borehole flow rates measured in the same boreholes using an impeller or heat pulse flowmeter. A comparison of flow rates measured using in‐well heat tracer tests to flow rates measured with an impeller flowmeter under the same conditions showed good agreement. A comparison of in‐well heat tracer test flow rate measurements to previously‐collected heat pulse flowmeter measurements indicates that the heat tracer test results produced borehole flow rates and flow profiles similar to those measured with the heat pulse flowmeter. The results of this study indicate that borehole flow rates determined from DTS‐monitored in‐well heat tracer tests are reasonable estimates of actual borehole flow rates. In addition, the range of borehole flow rates measurable by in‐well heat tracer tests spans from less than 10^?1 m/min to approximately 10¹ m/min, overlapping the ranges typically measurable with an impeller flowmeter or a heat pulse flowmeter, making in‐well heat tracer testing a versatile borehole flow logging tool.     

Evaluation of Two Low‐Cost–High‐Performance Adsorbent Materials in the Waste‐to‐Product Approach for the Removal of Pesticides from Drinking Water

This study evaluates the performance of two low cost and high performance adsorption materials, i.e., activated carbon produced from two natural waste products: Bamboo and coconut shell, in the removal of three pesticides from drinking water sources. Due to the fact that bamboo and coconut shell are abundant and inexpensive materials in many parts of the world, they respond to the “low‐cost” aspect. The adsorption capacities of two local adsorbents have been compared with commercial activated carbon to explore their potential to respond to the “high quality” aspect. Two pesticides were selected, namely dieldrin and chlorpyrifos, because they are commonly used in agriculture activities, and may remain in high concentrations in surface water used as drinking water sources. The results indicate that the adsorption of pesticides on activated carbons is influenced by physico‐chemical properties of the activated carbon and the pesticides such as the presence of an aromatic ring, and their molar mass. The activated carbon produced from bamboo can be employed as low‐cost and high performance adsorbent, alternative to commercial activated carbon for the removal of pesticides during drinking water production. The performance of activated carbon from bamboo was better due to its relatively large macroporosity and planar surface. The effect of adsorbent and pesticide characteristics on the performance was derived from batch experiments in which the adsorption behavior was studied on the basis of Freundlich isotherms.     

A Semi‐Structured MODFLOW‐USG Model to Evaluate Local Water Sources to Wells for Decision Support

In order to better represent the configuration of the stream network and simulate local groundwater‐surface water interactions, a version of MODFLOW with refined spacing in the topmost layer was applied to a Lake Michigan Basin (LMB) regional groundwater‐flow model developed by the U.S. Geological. Regional MODFLOW models commonly use coarse grids over large areas; this coarse spacing precludes model application to local management issues (e.g., surface‐water depletion by wells) without recourse to labor‐intensive inset models. Implementation of an unstructured formulation within the MODFLOW framework (MODFLOW‐USG) allows application of regional models to address local problems. A “semi‐structured” approach (uniform lateral spacing within layers, different lateral spacing among layers) was tested using the LMB regional model. The parent 20‐layer model with uniform 5000‐foot (1524‐m) lateral spacing was converted to 4 layers with 500‐foot (152‐m) spacing in the top glacial (Quaternary) layer, where surface water features are located, overlying coarser resolution layers representing deeper deposits. This semi‐structured version of the LMB model reproduces regional flow conditions, whereas the finer resolution in the top layer improves the accuracy of the simulated response of surface water to shallow wells. One application of the semi‐structured LMB model is to provide statistical measures of the correlation between modeled inputs and the simulated amount of water that wells derive from local surface water. The relations identified in this paper serve as the basis for metamodels to predict (with uncertainty) surface‐water depletion in response to shallow pumping within and potentially beyond the modeled area, see Fienen et al. (2015a).