Interpretation of transmissivity estimates from single-well pumping aquifer tests

Interpretation of single-well tests with the Cooper-Jacob method remains more reasonable than most alternatives. Drawdowns from 628 simulated single-well tests where transmissivity was specified were interpreted with the Cooper-Jacob straight-line method to estimate transmissivity. Error and bias as a function of vertical anisotropy, partial penetration, specific yield, and interpretive technique were investigated for transmissivities that ranged from 10 to 10,000 m(2)/d. Cooper-Jacob transmissivity estimates in confined aquifers were affected minimally by partial penetration, vertical anisotropy, or analyst. Cooper-Jacob transmissivity estimates of simulated unconfined aquifers averaged twice the known values. Transmissivity estimates of unconfined aquifers were not improved by interpreting results with an unconfined aquifer solution. Judicious interpretation of late-time data consistently improved estimates where transmissivity exceeded 250 m(2)/d in unconfined aquifers.     

Applications of the time-naught term in the Cooper and Jacob (1946) equation

The ability to manipulate analytical expressions for aquifer drawdown can provide insights into groundwater flow processes and assist with assessing strengths and weaknesses of aquifer parameter estimation methods. In the Cooper and Jacob (1946) parameter estimation method, the antilog of the horizontal-axis intercept in a plot of drawdown vs. log(time) is referred to as time naught (t(0)), which is used for estimating storativity. This article briefly reviews traditional uses of the time-naught concept and then spends time introducing new insights and applications involving (1) time-naught/distance relationships, including ways to compensate for certain missing data; (2) use of time naught in a simple method providing a quick visual check of which data in a Cooper-Jacob plot are suitable for use in linear regression; (3) application of time naught, as determined for one well, in estimating the later minimum time for which data from a distant well can be used in the Cooper-Jacob method; (4) development of relationships between drawdown and time naught; (5) use of time naught in a simple algebraic equation to estimate drawdown at smaller times than feasible using the Cooper-Jacob method; and (6) employment of time naught and a vertical-axis intercept on a plot of drawdown vs. log(time) for evaluating storativity. This information may be useful to new hydrogeologists or others interested in further developing their analytical well hydraulics skills.     

Determination of confined phreatic aquifer anisotropy

The determination of anisotropy in a phreatic aquifer by means of a pumping test can be carried out, under given conditions, with common analysis methods. This paper develops and applies an amended method based on a technique proposed by Boulton (1970). This new method enables a quick, simple determination of vertical hydraulic conductivity, given drawdown data from a pumping well and at least one observation well, that is applicable to confined phreatic aquifers.     

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.     

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.     

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.     

Methods to determine storativity of infinite confined aquifers from a recovery test

Starting from the equations of Theis and Cooper-Jacob, two new mathematical methods are proposed for interpreting the residual drawdown data for an infinite confined aquifer. Under Theis' assumptions and using the Cooper-Jacob approximation, the principal aquifer characteristics of transmissivity, pumping storativity, and recovery storativity are expressed without any correction or additional assumption. An actual case is used for illustration and confirms the validity of proposed equations and methods.     

Aquifer parameter estimation using an incremental area method

Theoretical well functions have been derived over the years to predict ground water level behaviour in aquifer systems under stress owing to groundwater extraction. The drawdown data collected during pump tests are typically analysed using graphical curve‐matching procedures to estimate aquifer parameters based on these well functions. Difficulty in aquifer characteristic identification and parameter estimation may arise when the field data do not perfectly match the drawdown curves obtained from the well functions. The present study provides a new method for the interpretation of aquifer pump tests which supplements the existing curve‐matching procedures in case ideal conditions do not exist; the proposed method provides a greater degree of flexibility in the data analysis for diagnostic tool purposes. The method, referred to as the Incremental Area Method (IAM) is based on integrating the logarithmic‐based drawdown curves within a discrete time and matching the results with a corresponding time integral of the Theis ( 1935 ) Well Function which governs ideal confined aquifers. The application of the proposed method to synthetically generated data and field data showed that IAM represents a viable method which yields information on potential non‐idealness of the aquifer and provides aquifer parameter estimates thus potentially overcoming drawdown data curve‐matching difficulties. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite layer method for flow in layered radial two-zone aquifer systems

A new finite layer method (FLM) is presented in this paper for transient flow analysis in layered radial two-zone aquifer systems. A radial two-zone system is an aquifer configuration in which a circular aquifer with finite radius is surrounded by a matrix possessing different permeability and storage properties. The aquifers can be pumped from fully or partially penetrating wells of infinitesimal radius. The trial function for drawdown is obtained through the use of piecewise linear correction functions in the present method. The trial function can satisfy the continuity conditions of flow and possess an appropriate continuity of C(0) at the two-zone interface. On the basis of Galerkin's method and the continuity condition of flow, the finite layer formulation is derived. The proposed method can cope with the anisotropy and layered heterogeneity in radial two-zone aquifer systems. Several numerical examples are presented to verify the validity of the present method through comparison with the analytical solution and the numerical results based on the finite difference method, in which a test of three-dimensional (3D) flow to a partially penetrating well in anisotropic two-zone aquifers is included. Furthermore, an additional application in simulating the two-zone flow in aquitard-aquifer systems is presented to demonstrate the applicability of FLM in modeling flow in more complex aquifer systems.     

稳定井流抽水含水层参数计算及富水性评价

针对目前单孔稳定流求参存在的问题,本文在分析新生界松散含水层条件及三次降深抽水过程基础上,利用其抽水试验恢复阶段的数据,分别求得各含水层多个参数,其值真实反映了含水层的实际情况。利用多元回归方法,求得降深与流量关系,通过其系数值大小分析,间接得出各含水层的富水性程度,为地下水的勘探与评价提供一定借鉴。     

Obtaining the steady‐state drawdown solutions of constant‐head and constant‐flux tests

The solutions of constant‐head and constant‐flux tests are commonly used to predict the temporal or spatial drawdown distribution or to determine aquifer parameters. Theis and Thiem equations, for instance, are well‐known transient and steady‐state drawdown solutions, respectively, of the constant‐flux test. It is known that the Theis equation is not applicable to the case where the aquifer has a finite boundary or the pumping time tends to infinity. On the other hand, the Thiem equation does not apply to the case where the aquifer boundary is infinite. However, the issue of obtaining the Thiem equation from the transient drawdown solution has not previously been addressed. In this paper, the drawdown solutions for constant‐head and constant‐flux tests conducted in finite or infinite confined aquifers with or without consideration of the effect of the well radius are examined comprehensively. Mathematical verification and physical interpretation of the solutions to these two tests converging or not converging to the Thiem equation are presented. The result shows that there are some finite‐domain solutions for these two tests that can converge to the Thiem equation when the time becomes infinitely large. In addition, the time criteria to give a good approximation to the finite‐domain solution by the infinite‐domain solution and the Thiem equation are investigated and presented. Copyright © 2008 John Wiley & Sons, Ltd.     

Assessing radial transmissivity variation in heterogeneous aquifers using analytical techniques

Pumping tests are one of the most commonly used in situ testing techniques for assessing aquifer hydraulic properties. Numerous researches have been conducted to predict the effects of aquifer heterogeneity on the groundwater levels during pumping tests. The objectives of the present work were as follows: (1) to predict drawdown conditions and to estimate aquifer properties during pumping tests undertaken in radially symmetric heterogeneous aquifers, and (2) to identify a method for assessing the transmissivity field along the radial coordinate in radially symmetric and fully heterogeneous transmissivity fields. The first objective was achieved by expanding an existing analytical drawdown formulation that was valid for a radially symmetric confined aquifer with two concentric zones around the pumping well to an N concentric zone confined aquifer having a constant transmissivity value within each zone. The formulation was evaluated for aquifers with three and four concentric zones to assess the effects of the transmissivity field on the drawdown conditions. The specific conditions under which aquifer properties could be identified using traditional methods of analysis were also evaluated. The second objective was achieved by implementing the inverse solution algorithm (ISA), which was developed for petroleum reservoirs to groundwater aquifer settings. The results showed that the drawdown values are influenced by a volumetric integral of a weighting function and the transmissivity field within the cone of depression. The weighting function migrates in tandem with the expanding cone of depression. The ability of the ISA to predict radially symmetric and log‐normally distributed transmissivity fields was assessed against analytical and numerical benchmarks. The results of this investigation indicated that the ISA method is a viable technique for evaluating the radial transmissivity variations of heterogeneous aquifer settings. Copyright © 2013 John Wiley & Sons, Ltd.     

Monitoring of drawdown pattern during pumping in an unconfined physical aquifer model

In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.     

Determination of Horizontal Aquifer Anisotropy with Three Wells

Existing methods for the determination of horizontal aquifer anisotropy by means of pumping tests require a minimum of four wells, one for water withdrawal and three for drawdown observations. This paper shows how the same methods can be used to determine anisotropy with as few as three wells, if at least two of them can be pumped in sequence. A field example is included. A method of analyzing data from more wells than the above minimum, by least squares, is also described.     

A new method for aquifer system identification and parameter estimation

The standard practice for assessing aquifer parameters is to match groundwater drawdown data obtained during pumping tests against theoretical well function curves specific to the aquifer system being tested. The shape of the curve derived from the logarithmic time derivative of the drawdown data is also very frequently used as a diagnostic tool to identify the aquifer system in which the pumping test is being conducted. The present study investigates the incremental area method (IAM) to serve as an alternative diagnostic tool for the aquifer system identification as well as a supplement to the aquifer parameter estimation procedure. The IAM based diagnostic curves for ideal confined, leaky, bounded and unconfined aquifers have been derived as part of this study, and individual features of the plots have been identified. These features were noted to be unique to each aquifer setting, which could be used for rapid evaluation of the aquifer system. The effectiveness of the IAM methodology was investigated by analyzing field data for various aquifer settings including leaky, unconfined, bounded and heterogeneous conditions. The results showed that the proposed approach is a viable method for use as a diagnostic tool to identify the aquifer system characteristics as well as to support the estimation of the hydraulic parameters obtained from standard curve matching procedures. Copyright © 2012 John Wiley & Sons, Ltd.     

An Analytical Method to Determine Ground Water Supply Well Network Designs

An analytical method is provided where the ground water practitioner can quickly determine the size (number of wells) and spacing of a well network capable of meeting a known ground water demand. In order to apply the method, two new parameters are derived that relate theoretical drawdown to the maximum drawdown that is achievable without mining the aquifer. The size of a well network is shown to be proportional to the ground water demand and inversely proportional to the transmissivity and available head. The spacing between wells in a supply well network is shown to be most sensitive to a derived parameter r _{HA/ 3}, which is related to the available head and the propagation of drawdown away from a theoretical well if the total ground water demand was applied to that well. The method can be used to quickly determine the required spacing between wells in well networks of various sizes that are completed in confined aquifers with no leakance.     

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.     

RATIO METHOD FOR DETERMINING CHARACTERISTICS OF IDEAL,LEAKY AND BOUNDED AQUIFERS

ABSTRACT

A Ratio Method is presented, which dispenses with the classical curve-matching procedure for computing characteristics of Ideal, Leaky and Bounded aquifers. This simple method consists of selecting any three values of observed drawdown, (t₀, s₀), (t₁, s₁) and (t₂, S₂), such that t₁ = 2 t₀ and t₂ = 100 t₀. Then, using the ratios, s_{0/s 1} and s₀/s₂, the values of u₀, W_(u0) as well as r/B (in the case of leaky aquifers) and r_i (in the case of aquifers with one boundary) are extrapolated from the charts presented. The extrapolated values are substituted in the standard equations to determine the aquifer characteristics. The Ratio Method may also have useful application in other branches of science where the curve-matching procedure is employed for calculation of parameters. (Key Words: Ratio Method; Aquifer characteristics; Ideal aquifer; Leaky aquifer; Bounded aquifer)     

The unusual and large drawdown response of buried-valley aquifers to pumping

The buried-valley aquifers that are common in the glacial deposits of the northern hemisphere are a typical case of the strip aquifers that occur in many parts of the world. Pumping from a narrow strip aquifer leads to much greater drawdown and much more distant drawdown effects then would occur in a sheet aquifer with a similar transmissivity and storage coefficient. Widely used theories for radial flow to wells, such as the Theis equation, are not appropriate for narrow strip aquifers. Previously published theory for flow to wells in semiconfined strip aquifers is reviewed and a practical format of the type curves for pumping-test analysis is described. The drawdown response of strip aquifers to pumping tests is distinctive, especially for observation wells near the pumped well. A case study is presented, based on extensive pumping test experience for the Estevan Valley Aquifer in southern Saskatchewan, Canada. Evaluation of groundwater resources in such buried-valley aquifers needs to take into account the unusually large drawdowns in response to pumping.     

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.