Spreadsheet Modeling of Slug Tests Using the van der Kamp Method

Slug testing is frequently employed to calculate aquifer transmissivity and hydraulic conductivity. The van der Kamp technique for interpreting slug test data which experience force-free water level oscillations is not routinely employed because it requires adjusting equations to match the observed well response data. This adjustment can be rapid and convenient when a commercial spreadsheet is employed.     

Predicting Collector Well Yields with MODFLOW

Groundwater flow models are commonly used to design new wells and wellfields. As the spatial scale of the problem is large and much local‐scale detail is not needed, modelers often utilize two‐dimensional (2D) or quasi three‐dimensional models based on the Dupuit‐Forchheimer assumption. Dupuit models offer a robust set of tools for simulating regional groundwater flow including interactions with surface waters, the potential for well interference, and varying aquifer properties and recharge rates. However, given an assumed operating water level or drawdown at a well screen, Dupuit models systematically overpredict well yields. For design purposes, this discrepancy is unacceptable, and a method for predicting accurate well yields is needed. While published methods exist for vertical wells, little guidance is available for predicting yields in horizontal screens or collector wells. In plan view, a horizontal screen has a linear geometry, and will likely extend over several neighboring cells that may not align with rows or columns in a numerical model. Furthermore, the model must account for the effects of converging three‐dimensional (3D) flow to the well screens and hydraulic interference among the well screens; these all depend on the design of a specific well. This paper presents a new method for simulating the yield of angled or horizontal well screens in numerical groundwater flow models, specifically using the USGS code MODFLOW. The new method is compared to a detailed, 3D analytic element model of a collector well in a field of uniform flow.     

Modeling multiaquifer wells with MODFLOW

Multiaquifer wells, i.e., wells that are open across more than one aquifer, can have a profound effect on the hydraulics of a ground water system. These wells change the physical system by establishing direct hydraulic links between otherwise isolated strata. Several methods are available to simulate multiaquifer wells in the context of comprehensive ground water flow simulators. In this paper, we review four methods to represent multiaquifer wells with the widely used code MODFLOW. These methods include a specialized code developed, but never formally released, by the U.S. Geological Survey (USGS), the Multi-Aquifer Well (MAW1) Package. An expanded implementation of the techniques in the MAW1 Package has been incorporated in the Multi-Node Well Package released recently by the USGS (Halford and Hanson 2002). We examine the performance of the methods in the context of a benchmarking study against the analytical solutions of Papadopulos (1966) and Sokol (1963). Our results demonstrate that results obtained with the MAW1 Package closely match exact solutions for pumping and nonpumping conditions, using both coarse and refined grids.     

Modeling Surface Water-Groundwater Interaction with MODFLOW: Some Considerations

The accuracy with which MODFLOW simulates surface water-groundwater interaction is examined for connected and disconnected losing streams. We compare the effect of different vertical and horizontal discretization within MODFLOW and also compare MODFLOW simulations with those produced by HydroGeoSphere. HydroGeoSphere is able to simulate both saturated and unsaturated flow, as well as surface water, groundwater and the full coupling between them in a physical way, and so is used as a reference code to quantify the influence of some of the simplifying assumptions of MODFLOW. In particular, we show that (1) the inability to simulate negative pressures beneath disconnected streams in MODFLOW results in an underestimation of the infiltration flux; (2) a river in MODFLOW is either fully connected or fully disconnected, while in reality transitional stages between the two flow regimes exist; (3) limitations in the horizontal discretization of the river can cause a mismatch between river width and cell width, resulting in an error in the water table position under the river; and (4) because coarse vertical discretization of the aquifer is often used to avoid the drying out of cells, this may result in an error in simulating the height of the groundwater mound. Conditions under which these errors are significant are investigated.     

Slug Tests in the Presence of Background Head Trends

We extend Bouwer and Rice (1976) slug test theory to incorporate background head trends that may be important in incompressible material of low permeability k. The extension, which features a convolution integral of the background head, is closed form for linear trends. A sensitivity study suggests that a rising background head can diminish the head changes associated with a slug-out test and underestimate k if it is ignored, as does falling background trend with a slug-in test. A falling background head can reinforce slug-in test head change and, if ignored, can overestimate k, as does a rising background head with a slug-out test. The simple extension is verified by field tests in glacial till and stratified drift deposits in eastern Massachusetts.     

Parallelization of MODFLOW Using a GPU Library

A new method based on a graphics processing unit (GPU) library is proposed in the paper to parallelize MODFLOW. Two programs, GetAb_CG and CG_GPU, have been developed to reorganize the equations in MODFLOW and solve them with the GPU library. Experimental tests using the NVIDIA Tesla C1060 show that a 1.6‐ to 10.6‐fold speedup can be achieved for models with more than 10⁵ cells. The efficiency can be further improved by using up‐to‐date GPU devices.     

Measurement of LNAPL Flux Using Single‐Well Intermittent Mixing Tracer Dilution Tests

The stability of subsurface Light Nonaqueous Phase Liquids (LNAPLs) is a key factor driving expectations for remedial measures at LNAPL sites. The conventional approach to resolving LNAPL stability has been to apply Darcy's Equation. This paper explores an alternative approach wherein single‐well tracer dilution tests with intermittent mixing are used to resolve LNAPL stability. As a first step, an implicit solution for single‐well intermittent mixing tracer dilution tests is derived. This includes key assumptions and limits on the allowable time between intermittent mixing events. Second, single‐well tracer dilution tests with intermittent mixing are conducted under conditions of known LNAPL flux. This includes a laboratory sand tank study and two field tests at active LNAPL recovery wells. Results from the sand tank studies indicate that LNAPL fluxes in wells can be transformed into formation fluxes using corrections for (1) LNAPL thicknesses in the well and formation and (2) convergence of flow to the well. Using the apparent convergence factor from the sand tank experiment, the average error between the known and measured LNAPL fluxes is 4%. Results from the field studies show nearly identical known and measured LNAPL fluxes at one well. At the second well the measured fluxes appear to exceed the known value by a factor of two. Agreement between the known and measured LNAPL fluxes, within a factor of two, indicates that single‐well tracer dilution tests with intermittent mixing can be a viable means of resolving LNAPL stability.     

Study of Water Relative Permeability in Fractures Using Well Tests and Radon: Gas Bubbles Effect

Few published data are available for two-phase flow in fractures from field studies. All measurements of relative permeability reported in the literature were done in laboratory-scale. The in situ water saturations are normally not known for multiphase flow in natural fractures; therefore, the direct measurements of relative permeability are difficult in field-scale. With the help of a case study before and after the 2008 M_w 5.4 Antung earthquake, groundwater radon was used as a tracer to determine the gas and water saturations in a small naturally fractured aquifer. Well tests were also conducted to estimate aquifer transmissivity before and after the 2008 Antung earthquake. Anomalous declines in both groundwater radon concentration and transmissivity were observed precursory to the 2008 Antung earthquake. Both declines are two precursory phenomena having a common effect of gas bubbles. Using the data from well tests and radon tracer, one data point of water relative permeability can be obtained for in situ fractures. This data point reveals strong phase interference between water and gas bubbles for multiphase flow in natural fractures. Both the data of well tests and radon tracer are essential to gain an improved understanding of mass transfer behavior of groundwater-dissolved gases between water and gas phases.