An Approximation of Inner Boundary Conditions for Wells Intersecting Highly Conductive Structures

Inner boundary conditions describe the interaction of groundwater wells with the surrounding aquifer during pumping and are associated with well-skin damage that limits water production and water derived from wellbore storage. Pumping test evaluations of wells during immediate and early time flow require assignment of inner boundary conditions. Originally, these concepts were developed for vertical well screens, and later transferred to wellbores intersecting highly conductive structures, such as preferential flow zones in fractured and karstic systems. Conceptual models for pumping test analysis in complex bedrock geology are often simplified. Classic analytical solutions generally lump or ignore conditions that limit or enhance well productivity along the well screen at the onset of pumping. Numerical solutions can represent well drawdowns in complex geological settings, such as karst systems, more precisely than many analytical solutions by accounting for additional physical processes and avoiding assumptions and simplifications. Suitable numerical tools for flow simulations in karst are discrete pipe-continuum models that account for various physical processes such as the transient hydraulics of wellbores intersecting highly conductive structures during pumping.     

Effect of wellbore storage and finite thickness skin on flow to a partially penetrating well in a phreatic aquifer

An analytical model is presented for the analysis of constant flux tests conducted in a phreatic aquifer having a partially penetrating well with a finite thickness skin. The solution is derived in the Laplace transform domain for the drawdown in the pumping well, skin and formation regions. The time-domain solution in terms of the aquifer drawdown is then obtained from the numerical inversion of the Laplace transform and presented as dimensionless drawdown–time curves. The derived solution is used to investigate the effects of the hydraulic conductivity contrast between the skin and formation, in addition to wellbore storage, skin thickness, delayed yield, partial penetration and distance to the observation well. The results of the developed solution were compared with those from an existing solution for the case of an infinitesimally thin skin. The latter solution can never approximate that for the developed finite skin. Dimensionless drawdown–time curves were compared with the other published results for a confined aquifer. Positive skin effects are reflected in the early time and disappear in the intermediate and late time aquifer responses. But in the case of negative skin this is reversed and the negative skin also tends to disguise the wellbore storage effect. A thick negative skin lowers the overall drawdown in the aquifer and leads to more persistent delayed drainage. Partial penetration increases the drawdown in the case of a positive skin; however its effect is masked by the negative skin. The influence of a negative skin is pronounced over a broad range of radial distances. At distant observation points the influence of a positive skin is too small to be reflected in early and intermediate time pumping test data and consequently the type curve takes its asymptotic form.     

Evaluation of Pumping Induced Flow in Observation Wells During Aquifer Testing

The vertical variation of drawdown around pumping wells generates an induced flow in the observation wells. A set of governing equations is presented to couple the drawdown variation and the vertical flux distribution in observation wells. A numerical example is performed to justify the governing equations and to verify the solution methods used by the simulation software WT. The example analyzes the effect of skin loss, wellbore storage, and vertical segmentation on the drawdown and induced flow in observation well during pumping. The evaluation of the Fairborn pumping test involves a vertically homogeneous and anisotropic water table aquifer, uniform well‐face drawdown conditions in the pumping well and simulation of the drawdown evolution in the observation well with and without the effect of induced flow. The computer calibrations resulted in small differences between the measured and simulated drawdown curves.     

Introduction to hydromechanical well tests in fractured rock aquifers

This article introduces hydromechanical well tests as a viable field method for characterizing fractured rock aquifers. These tests involve measuring and analyzing small displacements along with pressure transients. Recent developments in equipment and analyses have simplified hydromechanical well tests, and this article describes initial field results and interpretations during slug and constant-rate pumping tests conducted at a site underlain by fractured biotite gneiss in South Carolina. The field data are characterized by displacements of 0.3 μm to more than 10 μm during head changes up to 10 m. Displacements are a hysteretic function of hydraulic head in the wellbore, with displacements late in a well test always exceeding those at similar wellbore pressures early in the test. Displacement measurements show that hydraulic aperture changes during well tests, and both scaling analyses and field data suggest that T changed by a few percent per meter of drawdown during slug and pumping tests at our field site. Preliminary analyses suggest that displacement data can be used to improve estimates of storativity and to reduce nonuniqueness during hydraulic well tests involving single wells.     

A new analytical solution solved by triple series equations method for constant-head tests in confined aquifers

The constant-head pumping tests are usually employed to determine the aquifer parameters and they can be performed in fully or partially penetrating wells. Generally, the Dirichlet condition is prescribed along the well screen and the Neumann type no-flow condition is specified over the unscreened part of the test well. The mathematical model describing the aquifer response to a constant-head test performed in a fully penetrating well can be easily solved by the conventional integral transform technique under the uniform Dirichlet-type condition along the rim of wellbore. However, the boundary condition for a test well with partial penetration should be considered as a mixed-type condition. This mixed boundary value problem in a confined aquifer system of infinite radial extent and finite vertical extent is solved by the Laplace and finite Fourier transforms in conjunction with the triple series equations method. This approach provides analytical results for the drawdown in a partially penetrating well for arbitrary location of the well screen in a finite thickness aquifer. The semi-analytical solutions are particularly useful for the practical applications from the computational point of view.     

Measuring Flow Rate in Crystalline Bedrock Wells Using the Dissolved Oxygen Alteration Method

Determination of vertical flow rates in a fractured bedrock well can aid in planning and implementing hydraulic tests, water quality sampling, and improving interpretations of water quality data. Although flowmeters are highly accurate in flow rate measurement, the high cost and logistics may be limiting. In this study the dissolved oxygen alteration method (DOAM) is expanded upon as a low‐cost alternative to determine vertical flow rates in crystalline bedrock wells. The method entails altering the dissolved oxygen content in the wellbore through bubbler aeration, and monitoring the vertical advective movement of the dissolved oxygen over time. Measurements were taken for upward and downward flows, and under ambient and pumping conditions. Vertical flow rates from 0.06 to 2.30 Lpm were measured. To validate the method, flow rates determined with the DOAM were compared to pump discharge rates and found to be in agreement within 2.5%.     

Analytical Solutions for Determination of Non-Steady- State and Steady-State Capture Zones

This paper presents analytical solutions for determining non-steady-state capture zones produced by a single recovery well and steady-state capture zones produced by multiple recovery wells. Analysis of non-steady-slate capture zones is based on the lime-dependent location of caplure zone stagnation points and the geometric similarity between steady-slate and non-steady-state capture zones. The analytical solution of steady-state capture zones is obtained from spatial variations of discharge potential across the capture zone boundary. Both capture zone analyses are based on the assumptions of uniform flow field with a constant hydraulic conductivity, the Dupuit assumption of insignificant vertical flow, a negligible delayed yield, and a fully penetrating well with a constant pumping rate. For a ground water pump-and-trcat remediation program, the pumping rate and well location design variables can be adjusted to ensure containment of the ground water contaminant plume.     

Three-dimensional semi-analytical solution to groundwater flow in confined and unconfined wedge-shaped aquifers

The Laplace domain solutions have been obtained for three-dimensional groundwater flow to a well in confined and unconfined wedge-shaped aquifers. The solutions take into account partial penetration effects, instantaneous drainage or delayed yield, vertical anisotropy and the water table boundary condition. As a basis, the Laplace domain solutions for drawdown created by a point source in uniform, anisotropic confined and unconfined wedge-shaped aquifers are first derived. Then, by the principle of superposition the point source solutions are extended to the cases of partially and fully penetrating wells. Unlike the previous solution for the confined aquifer that contains improper integrals arising from the Hankel transform [Yeh HD, Chang YC. New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions. Adv Water Resour 2006;26:471–80], numerical evaluation of our solution is relatively easy using well known numerical Laplace inversion methods. The effects of wedge angle, pumping well location and observation point location on drawdown and the effects of partial penetration, screen location and delay index on the wedge boundary hydraulic gradient in unconfined aquifers have also been investigated. The results are presented in the form of dimensionless drawdown-time and boundary gradient-time type curves. The curves are useful for parameter identification, calculation of stream depletion rates and the assessment of water budgets in river basins.     

Leaky Aquifer Models to Simulate the Well Flow in Fractured Aquifers with Linear Interporosity Flow

Following Hemker and Maas (1987) the models of two or three leaky aquifers are applied to simulate the flow to vertical wells operating in the fractured or dual porosity aquifers. The software WellTest (WT) (Székely 2015) is used for calculating the drawdown and discharge rate variation. The comparative analysis with the independent analytical solutions by Boulton and Streltsova-Adams (1978), Warren and Root (1963), Kazemi et al. (1969) concluded with acceptable agreement between the WT simulation and the alternate calculation methods. The selected field tests have been conducted in fractured limestone aquifers. The pumping test west of Copenhagen shows an example of fractured aquifer with considerable negative skin effect at the well face. The flowing well Wafra W1 in Kuwait operates in the two-zone aquifer exhibiting sufficient vertical recharge via leakage beyond a circular domain of estimated radius of 2460 m.     

A Type-Curve Method for the Analysis of Pumping Tests with Piecewise-Linear Pumping Rates

Aquifer hydraulic parameters are commonly inferred from constant-rate pumping tests, while variable pumping rates are frequently encountered in actual field conditions. In this study, we propose a generally applicable dimensionless form of the analytical solution for variable-rate pumping tests in confined aquifers. In particular, we adopt a piecewise-linear fitting of variable pumping rates and propose a new type-curve method for estimating the hydraulic conductivity (K ) and specific storage (S_s ) of the investigated confined aquifer. For each test, a series of type curves, which depend on the variable pumping rates, the location of observation wells and the introduced first dimensionless inflection time, need to be provided for matching the observed drawdown data on a log-log graph. We first demonstrate the applicability and robustness of this method through a synthetic pumping test. Subsequently, we apply this method to analyze drawdown data from four pumping tests conducted within a multilayered aquifer/aquitard system in Wuxi city, Jiangsu Province, China. The parameter estimates are then compared with those reported by PEST. The K and S_s values estimated by the new type-curve method are found to be quite close to PEST-based estimates. Parameter estimation results demonstrate the difference in K and S_s values between observation wells. The difference could be attributed to the spatial heterogeneity in K and S_s . A future research topic may focus on the characterization of K and S_s heterogeneity with the currently available drawdown data from variable-rate pumping tests.     

Saltwater Upconing Due to Cyclic Pumping by Horizontal Wells in Freshwater Lenses

This article deals with the quantification of saltwater upconing below horizontal wells in freshwater lenses using analytical solutions as a computationally fast alternative to numerical simulations. Comparisons between analytical calculations and numerical simulations are presented regarding three aspects: (1) cyclic pumping; (2) dispersion; and (3) finite horizontal wells in a finite domain (a freshwater lens). Various hydrogeological conditions and pumping regimes within a dry half year are considered. The results show that the influence of elastic and phreatic storage (which are not taken into account in the analytical solutions) on the upconing of the interface is minimal. Furthermore, the analytical calculations based on the interface approach compare well with numerical simulations as long as the dimensionless interface upconing is below 1/3, which is in line with previous studies on steady pumping. Superimposing an analytical solution for mixing by dispersion below the well over an analytical solution based on the interface approach is appropriate in case the vertical flow velocity around the interface is nearly constant but should not be used for estimating the salinity of the pumped groundwater. The analytical calculations of interface upconing below a finite horizontal well compare well with the numerical simulations in case the distance between the horizontal well and the initial interface does not vary significantly along the well and in case the natural fluctuation of the freshwater lens is small. In order to maintain a low level of salinity in the well during a dry half year, the dimensionless analytically calculated interface upconing should stay below 0.25.     

DNAPL accumulation in wells and DNAPL recovery from wells: Model development and application to a laboratory study

Dense nonaqueous phase liquid (DNAPL) accumulation and recovery from wells cannot be accurately modeled through typical pressure or flux boundary conditions due to gravity segregation of water and DNAPL in the wellbore, the effects of wellbore storage, and variations of wellbore inflow and outflow rates with depth, particularly in heterogeneous formations. A discrete wellbore formulation is presented for numerical modeling of DNAPL accumulation in observation wells and DNAPL removal from recovery wells. The formulation includes fluid segregation, changing water and DNAPL levels in the well and the corresponding changes in fluid storage in the wellbore. The method was added to a three-dimensional finite difference model (CompSim) for three phase (water, gas, DNAPL) flow. The model predictions are compared to three-dimensional pilot scale experiments of DNAPL (benzyl alcohol) infiltration, redistribution, recovery, and water flushing. Model predictions match experimental results well, indicating the appropriateness of the model formulation. Characterization of mixing in the extraction well is important for predicting removal of highly soluble organic compounds like benzyl alcohol. A sensitivity analysis shows that the incorporation of hysteresis is critical for accurate prediction. Among the multiphase flow and transport parameters required for modeling, results are most sensitive to soil intrinsic permeability.     

A spreadsheet program for two-well tracer test data analysis

Two-well tracer tests are often conducted to investigate subsurface solute transport in the field. Analyzing breakthrough curves in extraction and monitoring wells using numerical methods is nontrivial due to highly nonuniform flow conditions. We extended approximate analytical solutions for the advection-dispersion equation for an injection-extraction well doublet in a homogeneous confined aquifer under steady-state flow conditions for equal injection and extraction rates with no transverse dispersion and negligible ambient flow, and implemented the solutions in Microsoft Excel using Visual Basic for Application (VBA). Functions were implemented to calculate concentrations in extraction and monitoring wells at any location due to a step or pulse injection. Type curves for a step injection were compared with those calculated by numerically integrating the solution for a pulse injection. The results from the two approaches are similar when the dispersivity is small. As the dispersivity increases, the latter was found to be more accurate but requires more computing time. The code was verified by comparing the results with published-type curves and applied to analyze data from the literature. The method can be used as a first approximation for two-well tracer test design and data analysis, and to check accuracy of numerical solutions. The code and example files are publicly available.     

Boundary depletion rate and drawdown in leaky wedge‐shaped aquifers

This study presents analytical solutions of the three‐dimensional groundwater flow to a well in leaky confined and leaky water table wedge‐shaped aquifers. Leaky wedge‐shaped aquifers with and without storage in the aquitard are considered, and both transient and steady‐state drawdown solutions are derived. Unlike the previous solutions of the wedge‐shaped aquifers, the leakages from aquitard are considered in these solutions and unlike similar previous work for leaky aquifers, leakage from aquitards and from the water table are treated as the lower and upper boundary conditions. A special form of finite Fourier transforms is used to transform the z‐coordinate in deriving the solutions. The leakage induced by a partially penetrating pumping well in a wedge‐shaped aquifer depends on aquitard hydraulic parameters, the wedge‐shaped aquifer parameters, as well as the pumping well parameters. We calculate lateral boundary dimensionless flux at a representative line and investigate its sensitivity to the aquitard hydraulic parameters. We also investigate the effects of wedge angle, partial penetration, screen location and piezometer location on the steady‐state dimensionless drawdown for different leakage parameters. Results of our study are presented in the form of dimensionless flux‐dimensionless time and dimensionless drawdown‐leakage parameter type curves. The results are useful for evaluating the relative role of lateral wedge boundaries and leakage source on flow in wedge‐shaped aquifers. This is very useful for water management problems and for assessing groundwater pollution. The presented analytical solutions can also be used in parameter identification and in calculating stream depletion rate and volume. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

A COMPUTERIZED TECHNIQUE FOR ESTIMATING THE HYDRAULIC CONDUCTIVITY OF AQUIFERS FROM SPECIFIC CAPACITY DATA

Abstract. Specific capacity data obtained from well construction reports can provide useful estimates of hydraulic conductivity (K). A simple computer program has been developed which can correct specific capacity data for partial penetration and well loss and, using an iterative technique, provide rapid estimates of K at hundreds of data points. The program allows easy data handling and is easily linked with existing statistical programs or contour mapping routines. The method was tested at two field sites in Wisconsin, one underlain by a sandy outwash aquifer, the other by fractured dolomite. In both areas, estimates of K from corrected specific capacity data agree reasonably well with data from pumping tests.     

A method for evaluating horizontal well pumping tests

Predicting the future performance of horizontal wells under varying pumping conditions requires estimates of basic aquifer parameters, notably transmissivity and storativity. For vertical wells, there are well-established methods for estimating these parameters, typically based on either the recovery from induced head changes in a well or from the head response in observation wells to pumping in a test well. Comparable aquifer parameter estimation methods for horizontal wells have not been presented in the ground water literature. Formation parameter estimation methods based on measurements of pressure in horizontal wells have been presented in the petroleum industry literature, but these methods have limited applicability for ground water evaluation and are based on pressure measurements in only the horizontal well borehole, rather than in observation wells. This paper presents a simple and versatile method by which pumping test procedures developed for vertical wells can be applied to horizontal well pumping tests. The method presented here uses the principle of superposition to represent the horizontal well as a series of partially penetrating vertical wells. This concept is used to estimate a distance from an observation well at which a vertical well that has the same total pumping rate as the horizontal well will produce the same drawdown as the horizontal well. This equivalent distance may then be associated with an observation well for use in pumping test algorithms and type curves developed for vertical wells. The method is shown to produce good results for confined aquifers and unconfined aquifers in the absence of delayed yield response. For unconfined aquifers, the presence of delayed yield response increases the method error.     

SUPRPUMP: An Interactive Program for Well Test Analysis and Design

Abstract. We have developed a program which aids in the design and analysis of pumping tests and slug tests. In design mode, the program emphasizes calculation and plotting of the sensitivities of drawdown (or head) to well function parameters. In analysis mode, the program can analyze a given set of experimental data. For pumping tests, the program allows multiple observation wells and multiple variable-rate pumping wells. The program is written in a modular fashion, allowing easy addition of well functions to the currently existing library. An example based on a hypothetical pumping test illustrates the utility of sensitivity analysis for well test design.     

Well storage effect during pumping tests in an aquifer of low permeability

Abstract

Well storage effect is generally considered while interpreting pumping test data from large diameter wells. However, in an aquifer of low permeability, the well storage is found to be significant during pumping tests conducted on bore wells. The interpretation of such data gives ambiguous results unless well storage effect is taken into account. A field example is presented to illustrate the difficulty in interpretation of the pumping test data. In order to take into account the well storage effect, a finite difference approach of interpreting pumping test data is suggested.