A single recovery type curve from Theis'' exact solution

The Theis type curve matching method and the Cooper-Jacob semilog method are commonly used for estimation of transmissivity and storage coefficient of infinite, homogeneous, isotropic, confined aquifers from drawdown data of a constant rate pumping test. Although these methods are based on drawdown data, they are often applied indiscriminately to analyze both drawdown and recovery data. Moreover, the limitations of drawdown type curve to analyze recovery data collected after short pumping times are not well understood by the practicing engineers. This often may result in an erroneous interpretation of such recovery data. In this paper, a novel but simple method is proposed to determine the storage coefficient as well as transmissivity from recovery data measured after the pumping period of an aquifer test. The method eliminates the dependence on pumping time effects and has the advantage of employing only one single recovery type curve. The method based on the conversion of residual drawdown to recovered drawdown (buildup) data plotted versus a new equivalent time (delta(t) x t(p)/t(p) + delta(t)). The method uses the recovery data in one observation point only, and does not need the initial water level h0, which may be unknown. The accuracy of the method is checked with three sets of field data. This method appears to be complementary to the Cooper-Jacob and Theis methods, as it provides values of both storage coefficient and transmissivity from recovery data, regardless of pumping duration.     

A Potential-Based Inversion of Unconfined Steady-State Hydraulic Tomography

The importance of estimating spatially variable aquifer parameters such as transmissivity is widely recognized for studies in resource evaluation and contaminant transport. A useful approach for mapping such parameters is inverse modeling of data from series of pumping tests, that is, via hydraulic tomography. This inversion of field hydraulic tomographic data requires development of numerical forward models that can accurately represent test conditions while maintaining computational efficiency. One issue this presents is specification of boundary and initial conditions, whose location, type, and value may be poorly constrained. To circumvent this issue when modeling unconfined steady-state pumping tests, we present a strategy that analyzes field data using a potential difference method and that uses dipole pumping tests as the aquifer stimulation. By using our potential difference approach, which is similar to modeling drawdown in confined settings, we remove the need for specifying poorly known boundary condition values and natural source/sink terms within the problem domain. Dipole pumping tests are complementary to this strategy in that they can be more realistically modeled than single-well tests due to their conservative nature, quick achievement of steady state, and the insensitivity of near-field response to far-field boundary conditions. After developing the mathematical theory, our approach is first validated through a synthetic example. We then apply our method to the inversion of data from a field campaign at the Boise Hydrogeophysical Research Site. Results from inversion of nine pumping tests show expected geologic features, and uncertainty bounds indicate that hydraulic conductivity is well constrained within the central site area.     

Determination of Storage Coefficients During Pumping and Recovery

An aquifer test is used mostly to determine the storage coefficient and transmissivity. Although residual drawdown data are widely used in estimating the transmissivity of aquifers, the estimation of storage coefficients with recovery data is controversial. Some researchers have proposed methods to estimate storage coefficients with recovery data by assuming equality of storage coefficients for the recovery and pumping periods (S = S′). The aim of this study is to determine storage coefficients without such an assumption, that is, S≠S′. The method is a modified version of Banton‐Bangoy's method without considering drawdown data due to pumping. Drawdown is plotted vs. the logarithmic ratio (t′/t) or time since pumping stopped to the duration of pumping and the ratio of storage coefficient during recovery to the storage coefficient from the pumping period (S′/S). The method is verified with one case study and two synthetic examples. Thus, it is possible to determine storage coefficient of pumping period accurately without any data from pumping period by recovery data.     

A Method for Estimating Field-Scale Mass Transfer Rate Parameters and Assessing Aquifer Cleanup Times

A general method for estimating ground-water solute mass transfer rate parameters from field test data is presented. The method entails matching solute concentration and hydraulic head data collected during the recovery phase of a pumping test through application of a simulation-regression technique. Estimation of hydraulic conductivity and mass transfer rate parameter values is performed by fitting model simulations to the data. Parameter estimates are utilized to assess cleanup times for pump-and-treat aquifer remediation scenarios. Uncertainty in the cleanup time estimate is evaluated using statistical information obtained with the parameter estimation technique. Application of the method is demonstrated using a hypothetical ground-water flow and solute transport system. Simulations of field testing, parameter estimation, and remedial time frames are performed to evaluate the usefulness of the method. Sets of random noise that signify potential field and laboratory measurement errors are combined with the hypothetical data to provide rigorous testing of the method. Field tests are simulated using ranges of values for data noise, the mass transfer rate parameters, the test pumping rates, and the duration of recovery monitoring to evaluate their respective influence on parameter and cleanup time estimates. The demonstration indicates the method is capable of yielding accurate estimates of the solute mass transfer rate parameters. When the parameter values for the hypothetical system are well estimated, cleanup time predictions are shown to be more accurate than when calculated using the local equilibrium assumption.     

Analytical and Semi‐Analytical Tools for the Design of Oscillatory Pumping Tests

Oscillatory pumping tests—in which flow is varied in a periodic fashion—provide a method for understanding aquifer heterogeneity that is complementary to strategies such as slug testing and constant‐rate pumping tests. During oscillatory testing, pressure data collected at non‐pumping wells can be processed to extract metrics, such as signal amplitude and phase lag, from a time series. These metrics are robust against common sensor problems (including drift and noise) and have been shown to provide information about aquifer heterogeneity. Field implementations of oscillatory pumping tests for characterization, however, are not common and thus there are few guidelines for their design and implementation. Here, we use available analytical solutions from the literature to develop design guidelines for oscillatory pumping tests, while considering practical field constraints. We present two key analytical results for design and analysis of oscillatory pumping tests. First, we provide methods for choosing testing frequencies and flow rates which maximize the signal amplitude that can be expected at a distance from an oscillating pumping well, given design constraints such as maximum/minimum oscillator frequency and maximum volume cycled. Preliminary data from field testing helps to validate the methodology. Second, we develop a semi‐analytical method for computing the sensitivity of oscillatory signals to spatially distributed aquifer flow parameters. This method can be quickly applied to understand the “sensed” extent of an aquifer at a given testing frequency. Both results can be applied given only bulk aquifer parameter estimates, and can help to optimize design of oscillatory pumping test campaigns.     

A GENERAL PURPOSE MICROCOMPUTER AQUIFER TEST EVALUATION TECHNIQUE

Abstract. Although determination of aquifer characteristics from pumping test data is generally carried out using type curves or other graphical techniques, a number of computer methods have been developed recently for this purpose. Based on the principle of least squares, these methods of nonlinear regression analysis can be applied to any flow system for which analytical expressions of the drawdown distribution are known. In view of the growing general interest in the application of microcomputers in ground-water hydrology, a BASIC routine has been developed for estimating any number of aquifer parameters. The least squares solution is calculated by Marquardt's algorithm, using the singular-value decomposition of the Jacobian matrix. The robust computing method obtained can be applied to all kinds of pumping tests. Aquifer characteristics as well as their standard deviations are computed with optimal speed and accuracy. The technique is demonstrated by a simple application to steady flow in a leaky aquifer and an example is provided. Other applications are easily implemented and programs for unsteady-state aquifer tests, recovery tests and multiple aquifer tests are available.     

A review of thermal response test analysis using pumping test concepts

The design of ground-coupled heat pump systems requires knowledge of the thermal properties of the subsurface and boreholes. These properties can be measured with in situ thermal response tests (TRT), where a heat transfer fluid flowing in a ground heat exchanger is heated with an electric element and the resulting temperature perturbation is monitored. These tests are analogous to standard pumping tests conducted in hydrogeology, because a system that is initially assumed at equilibrium is perturbed and the response is monitored in time, to assess the system's properties with inverse modeling. Although pumping test analysis is a mature topic in hydrogeology, the current analysis of temperature measurements in the context of TRTs is comparatively a new topic and it could benefit from the application of concepts related to pumping tests. The purpose of this work is to review the methodology of TRTs and improve their analysis using pumping test concepts, such as the well function, the superposition principle, and the radius of influence. The improvements are demonstrated with three TRTs. The first test was conducted in unsaturated waste rock at an active mine and the other two tests aimed at evaluating the performance of thermally enhanced pipe installed in a fully saturated sedimentary rock formation. The concepts borrowed from pumping tests allowed the planning of the duration of the TRTs and the analysis of variable heat injection rate tests accounting for external heat transfer and temperature recovery, which reduces the uncertainty in the estimation of thermal properties.     

A Simple, Low-Cost Method to Monitor Duration of Ground Water Pumping

Monitoring ground water withdrawals for agriculture is a difficult task, while agricultural development leads frequently to overexploitation of the aquifers. To fix the problem, sustainable management is required based on the knowledge of water uses. This paper introduces a simple and inexpensive direct method to determine the duration of pumping of a well by measuring the temperature of its water outlet pipe. A pumping phase is characterized by a steady temperature value close to ground water temperature. The method involves recording the temperature of the outlet pipe and identifying the different stages of pumping. It is based on the use of the low-cost and small-size Thermochron iButton temperature logger and can be applied to any well, provided that a water outlet pipe is accessible. The temperature time series are analyzed to determine the duration of pumping through manual and automatic posttreatments. The method was tested and applied in South India for irrigation wells using electricity-powered pumps. The duration of pumping obtained by the iButton method is fully consistent with the duration of power supply (1.5% difference).     

天津地区含水层类型判别与参数分析

根据含水层数学模型的降深-时间双对数和半对数曲线特征,可以确定含水层类型;依据选定的含水层模型,调整其中水文地质参数以拟合抽水试验中观测井降深-时间实测数据,可以反演水文地质参数。将以上方法应用于天津地区6组工程抽水试验,获得试验含水层类型及其水文地质参数。该方法可用于其他相关工程,天津地区和工程建设有关的50m埋深范围内含水层属性及其水文地质参数取值范围,可参考借鉴本文计算结果。     

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.     

Identification of Pumping Influences in Long‐Term Water Level Fluctuations

Identification of the pumping influences at monitoring wells caused by spatially and temporally variable water supply pumping can be a challenging, yet an important hydrogeological task. The information that can be obtained can be critical for conceptualization of the hydrogeological conditions and indications of the zone of influence of the individual pumping wells. However, the pumping influences are often intermittent and small in magnitude with variable production rates from multiple pumping wells. While these difficulties may support an inclination to abandon the existing dataset and conduct a dedicated cross‐hole pumping test, that option can be challenging and expensive to coordinate and execute. This paper presents a method that utilizes a simple analytical modeling approach for analysis of a long‐term water level record utilizing an inverse modeling approach. The methodology allows the identification of pumping wells influencing the water level fluctuations. Thus, the analysis provides an efficient and cost‐effective alternative to designed and coordinated cross‐hole pumping tests. We apply this method on a dataset from the Los Alamos National Laboratory site. Our analysis also provides (1) an evaluation of the information content of the transient water level data; (2) indications of potential structures of the aquifer heterogeneity inhibiting or promoting pressure propagation; and (3) guidance for the development of more complicated models requiring detailed specification of the aquifer heterogeneity.     

Hydraulic parameters of a multilayered aquifer system in Kuwait City

The Kuwait Group consists mainly of clastic sediments overlying unconformably the Dammam Formation of Tertiary age. The Kuwait Group is generally divided into three main hydrostratigraphic units: the upper and lower aquifers separated by an aquitard. The upper aquifer is further divided into the water table aquifer, an aquitard and a semiconfined aquifer. This semiconfined unit was pumped and the drawdowns were observed in piezometers screened in various subunits of the Kuwait Group. Some pumping tests of short duration were carried out in the top water table aquifer as well. These tests showed that the subunits of the Kuwait Group are hydraulically interconnected to a varying degree.

The pumping test data were analysed using conventional analytical solutions. The semiconfined pumping test was also simulated by a quasi-three-dimensional model using a leaky multiaquifer modelling technique. The initial hydraulic parameters were improved manually in the model till best fit drawdowns were obtained.

The final parameters obtained by simulation of the pumping tests were used in designing a pilot drainage system for the control of a rising groundwater table in parts of Kuwait City.     

A Comparative Study of Specific Yield Determinations for a Shallow Sand Aquifer

Using the type-curve methods of Boulton (1963) and Neuman (1972), and comparisons, at various times, of the cumulative volume of water pumped to the volume of the water-table drawdown cone (volume-balance method), values of specific yield were obtained from pumping test data from numerous piezometers in an unconfined sand aquifer. The long-term value of specific yield for the aquifer was determined from measurements of the laboratory drainage curve of the aquifer material. The volume-balance method gave specific yield values of 0.02, 0.05, 0.12, 0.20, 0.23, and 0.25 at times of 0.25, 0.66, 10, 26, 45, and 65 hours, respectively, indicating a gradual increase in specific yield and an asymptotic approach to the long-term value of 0.30 determined from the laboratory method. The type-curve methods provided values of 0.07 and 0.08, which correspond to the volume-balance values at early times, but which are less than one-third of the value obtained from the laboratory method and from the volume-balance method applied at the end of the pumping test (2.7 days). The type-curve procedures therefore provide unrealistically low values of specific yield for application to problems concerning the long-term yield characteristics of the aquifer. The observed trend towards increasing values of specific yield with increasing duration of pumping, and the vertical hydraulic head profiles that were measured during the pumping test indicate that both delayed drainage from above the water table and downward hydraulic gradients in the saturated zone can be important hydraulic effects contributing to the delayed-drawdown segment that is characteristic of time-drawdown graphs for unconfined aquifers.     

Estimating Flow Using Tracers and Hydraulics in Synthetic Heterogeneous Aquifers

Regional ground water flow is most usually estimated using Darcy's law, with hydraulic conductivities estimated from pumping tests, but can also be estimated using ground water residence times derived from radioactive tracers. The two methods agree reasonably well in relatively homogeneous aquifers but it is not clear which is likely to produce more reliable estimates of ground water flow rates in heterogeneous systems. The aim of this paper is to compare bias and uncertainty of tracer and hydraulic approaches to assess ground water flow in heterogeneous aquifers. Synthetic two-dimensional aquifers with different levels of heterogeneity (correlation lengths, variances) are used to simulate ground water flow, pumping tests, and transport of radioactive tracers. Results show that bias and uncertainty of flow rates increase with the variance of the hydraulic conductivity for both methods. The bias resulting from the nonlinearity of the concentration–time relationship can be reduced by choosing a tracer with a decay rate similar to the mean ground water residence time. The bias on flow rates estimated from pumping tests is reduced when performing long duration tests. The uncertainty on ground water flow is minimized when the sampling volume is large compared to the correlation length. For tracers, the uncertainty is related to the ratio of correlation length to the distance between sampling wells. For pumping tests, it is related to the ratio of correlation length to the pumping test's radius of influence. In regional systems, it may be easier to minimize this ratio for tracers than for pumping tests.     

Geophysical Assessment of Groundwater Potential: A Case Study from Mian Channu Area,Pakistan

An integrated study using geophysical method in combination with pumping tests and geochemical method was carried out to delineate groundwater potential zones in Mian Channu area of Pakistan. Vertical electrical soundings (VES) using Schlumberger configuration with maximum current electrode spacing (AB/2 = 200 m) were conducted at 50 stations and 10 pumping tests at borehole sites were performed in close proximity to 10 of the VES stations. The aim of this study is to establish a correlation between the hydraulic parameters obtained from geophysical method and pumping tests so that the aquifer potential can be estimated from the geoelectrical surface measurements where no pumping tests exist. The aquifer parameters, namely, transmissivity and hydraulic conductivity were estimated from Dar Zarrouyk parameters by interpreting the layer parameters such as true resistivities and thicknesses. Geoelectrical succession of five‐layer strata (i.e., topsoil, clay, clay sand, sand, and sand gravel) with sand as a dominant lithology was found in the study area. Physicochemical parameters interpreted by World Health Organization and Food and Agriculture Organization were well correlated with the aquifer parameters obtained by geoelectrical method and pumping tests. The aquifer potential zones identified by modeled resistivity, Dar Zarrouk parameters, pumped aquifer parameters, and physicochemical parameters reveal that sand and gravel sand with high values of transmissivity and hydraulic conductivity are highly promising water bearing layers in northwest of the study area. Strong correlation between estimated and pumped aquifer parameters suggest that, in case of sparse well data, geophysical technique is useful to estimate the hydraulic potential of the aquifer with varying lithology.     

Comparison of Two Computational Methods for Estimating Transmissivity Based on the Picking Equation

A comparison is presented of two computational methods, PICKINGmodel and PPC-Recovery, to estimate transmissivities based on the Picking equation using water-level recovery data from brief pumping tests of relatively low-yielding domestic wells. The tests were performed by the United States Geological Survey (USGS) in 50 domestic bedrock wells in south-central New York State, and USGS staff performed the analysis using PICKINGmodel based on the Picking equation. The results indicated that the estimated transmissivities ranged from 0.86 to 2900 ft²/d (0.080 to 270 m²/d) with a median of 41 ft²/d (3.8 m²/d). The same data were later analyzed using PPC-Recovery also based on the Picking equation. The two sets of estimated transmissivities were compared and statistically had the same median value at a probability of 95%. In another analysis, the PPC-Recovery method was applied to the same data that had been truncated at the point when the slope of the recovery data curve began to deviate from a straight line aligned with the middle portion of the recovery data. Comparing these resulting estimates of transmissivity with values originally obtained using the PICKINGmodel, the two had statistically the same median value for transmissivity at a probability of 95%. It was concluded that using PPC-Recovery in this manner to estimate transmissivity in low-yielding domestic wells will yield transmissivity values sufficiently close to the results had PICKINGmodel been used, and with less time and effort.     

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.     

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.     

Identification of aquifer parameters from residual drawdowns: an optimization approach

Abstract

An optimization method based on nonlinear least squares is proposed for the identification of aquifer transmissivity and two storage coefficients, one for the pumping phase and the other for the recovery phase, using only residual drawdowns at an observation well. The method is applicable irrespective of whether pumping and recovery storage coefficients are the same or are different and the argument of well function during recovery is less than or greater than 0.01. The method can check for the changed value of storage coefficient during the recovery phase and can quantify this change. Use of the method is illustrated through examples, which show that the parameters are estimated reliably.     

On Simulation and Analysis of Variable‐Rate Pumping Tests

Analytical solutions for constant‐rate pumping tests are widely used to infer aquifer properties. In this note, we implement a methodology that approximates the time‐varying pumping record as a series of segments with linearly varying pumping rates. We validate our approach using an analytical solution for a sinusoidally varying pumping test. We also apply our methodology to analyze synthetic test data and compare the results with those from a commonly used method where rate variations are represented by a series of constant‐rate steps.