Determination of the hydraulic parameters (transmissivity and storage coefficients) of a confined aquifer is important for effective groundwater resources. For this purpose, the residual drawdowns have been in use to estimate the aquifer parameters by the classical Theis recovery method. The proposed method of this paper depends on a straight-line through the field data and it helps to calculate the parameters quickly without any need for long-term pumping data. It is based on the expansion series of the Theis well function by consideration of three terms, and this approach is valid for the dimensionless time factor u′=S′r2/4Tt′≤0.2. The method can be applied reliably to extensive and homogeneous confined aquifers resulting in different storage coefficients during the pumping and recovery periods S≠S′. It presents a strength methodology for the parameters decision making from the residual data in the groundwater field of civil engineering.
1. Introduction
One of the practical ways to estimate the aquifer parameter is to measure the water level rise by time in the production or observation wells after the pumping test stoppage. This is referred to as the recovery test which starts just after the pump shut. The recovery method serves as a check and alternative to the pumping test. The parameters’ estimations from both tests are practically equal to each other if the Theis [1] assumptions are satisfied.
The residual drawdown measurement at any time during the recovery period is the difference between the observed water level and the prepumping static water level. The recovery drawdown is known as the difference between the total drawdown at the end of pumping and the residual drawdown [2, 3]. With the Theis recovery method, the transmissivity can be estimated easily using pumping well recovery data, but the storage coefficient cannot be calculated due to wellbore storage effects, unknown effective radius, and difficulty in finding the time of zero recovery as needed for the application of Cooper and Jacob [4] method. However, in cases of measurements from the observation wells there will not be such restrictive effects. The Theis recovery method considers the late-time residual drawdowns with the Cooper and Jacob formulations and it estimates the transmissivity and the ratio of storativity values during pumping and recovery periods. Theis [1] observed that a straight-line through the residual drawdowns (s′) versus t/t′ plot (t indicates the total time since pump start while t′ is the time since pump shut). This plot on the semilogarithmic graph paper passes below the origin (t/t′=1, s′=0) giving the value of t/t′>1 for a zero residual drawdown and that is the reason why different storage coefficient estimations are valid for the pumping and recovery periods. On the other hand, Jacob [5] observed that the storage coefficient estimation is generally greater during the pumping period than the recovery period.
Bruin and Hudson [6] proposed the time-recovery drawdown graph to find the time of zero recovery. The method depends on the extension of the time-drawdown pumping test data, which can also be applied to the time-recovery graph. Later, USDI [7] gave an alternative method for determining the storage coefficient (S) as follows: (1)S=2.25Tt′/r2log-1sp-s′/Δsp-s′,where T is the transmissivity, sp is pumping period drawdown projected to time t′ at any radial distance, r, s′ is residual drawdown at time t′, sp-s′ is recovery at time t′, and Δsp-s′ is slope of the time-recovery graph. In this formulation log-1x corresponds to the antilogarithm as 10x.
The recovery analysis is also investigated by many researchers. For instance, Case et al. [8] developed convenient equations in the forms of a series based on the Theis recovery equation using the residual drawdown data. Agarwal [9] developed a method for recovery test analysis which is widely used by petroleum engineers. Through a simple transformation of the data from a recovery test, Agarwal method allows one to apply the same diagnostic principles and type curves used for drawdown analysis in the interpretation of recovery data. Ramey [10] presented the type curves for drawdown during the pumping and recovery periods, where the recovery times are plotted as large times. Mishra and Chachadi [11] obtained the recovery type curves for large-diameter pumping wells by discrete kernel approach while Şen [12] presented an analytical solution and a set of type curves for the drawdown distribution in a large-diameter well recovery. Later, Yeh and Wang [13] developed a mathematical model for describing the residual drawdown by taking into consideration the pumping drawdown distribution in addition to the effects of well radius and wellbore storage. They obtained the Laplace domain solution for the residual drawdown. Goode [14] proposed a set of graphical recovery type curves based on Theis’ [1] exact well-function solution. These type curves depend on the dimensionless duration of pumping. Ballukraya and Sharma [15] proposed an approach derived from the Cooper-Jacob equation for estimating storativity by using residual drawdown measurements. Banton and Bangoy [16] presented a graphical method with the first three terms of the Theis series approximation. The method involves three separate plots with the equality of the storage coefficients in pumping and recovery periods (S=S′), but this approach requires at least two observation wells. Singh [17] proposed a numerical method by considering the derivative of the Theis recovery equation. Zheng et al. [18] suggested a straight-line method based on the Cooper-Jacob approximation for the extended pumping period and the first three terms in the expansion are from the well function for the recovery period. The method considers that S=S′. Singh [19] presented an optimization method based on nonlinear least-squares for the identification of the transmissivity and the storage coefficients in the pumping and recovery periods. Samani et al. [20] used a derivative analysis of pumping and recovery test data to estimate the hydraulic parameters in a heterogeneous aquifer. They showed that the drawdown-derivative analysis improves estimation of aquifer parameters and identification of different forms of heterogeneity. Kambhammettu and King [21] estimated the transmissivity and storage coefficient using a generalized MATLAB code with the conventional Levenberg-Marquardt algorithm. They considered the residual drawdowns measurements from a single observation well. Ashjari [22] determined the transmissivity and storage coefficients from residual data in case of S≠S′ by using a modified version of Banton and Bangoy [15] method. This method is basically fitting a straight-line to a plot of residual drawdown versus square of radial distance at the same time.
In this study, another straight-line method is proposed using the first three terms from the expansion of the well function for the pumping and recovery periods. The method offers the use of spreadsheet for the inequality S≠S′, which implies different storage coefficients during pumping and recovery periods. The procedure involves a linear regression line and its coefficients’ estimations based on a set of recovery data from a single observation well. It is valid for the dimensionless time factor, u′≤0.2.
2. Proposed Method
In a homogeneous isotropic confined aquifer with infinite domain without the well storage, Theis [1] gave the residual drawdown expression for an observation well as follows:(2)sr,t′=Q4πT∫u∞e-xxdx-∫u′∞e-xxdx,where s(r,t′) is the residual drawdown at any distance r and at any recovery time t′, Q is the constant rate (discharge) towards the pumping well during the pumping and recovery periods, T is the transmissivity, u=r2S/4Tt is the dimensionless time factor for the pumping period, u′=r2S′/4Tt′ is another dimensionless time factor for the recovery period, S and S′ are the storage coefficients of aquifer during the pumping and recovery periods, t=tp+t′, and tp is the time of pumping. This expression can be considered after the first three terms of the exponential function series as follows:(3a)sr,t′=Q4πT-0.5772-lnu+u-0.5772-lnu′+u′,(3b)sr,t′=Q4πTlnS′Stt′+r2S′4TSS′1t-1t′.
The error involved in adopting (3b) instead of (2) is less than 1% for u′≤0.2. Theis [1] proposed the first two terms of the series in (2) by considering that S/S′=1 in order to estimate the aquifer transmissivity only. Theis approach is valid for u′≤0.01. Hence, (3b) considers more recovery data than Theis method. Equation (3b) can be rewritten to estimate the aquifer parameters as(4)W=aτ-b,(5)W=sr,t′1/t′-S/S′1/t,(6)a=Q4πT,(7)τ=lnS′/St/t′1/t′-S/S′1/t,(8)b=ar2S′4T.Equation (4) presents a straight-line between τ and W. The first approximation of the aquifer parameter estimations can be obtained from the time residual drawdown data on a spreadsheet with the exception of a few early time instances, and the aquifer parameters can be estimated after fitting a straight-line to the field data and the coefficients of the regression line result from (6) and (8). By considering the calculated parameters, the dimensionless time factors, u′, should be determined especially for the first data values. If u′ is greater than 0.2, then the straight-line should be rearranged. Furthermore, the ratio of S/S′ may be easily investigated with various straight-lines.
3. Application and Discussion
Two data sets are used to illustrate the application of the proposed method. The first set of data is taken from the USDI [7]. The data is recorded in an observation well located at 30.48 m from the pumping well. The well is pumped during 800 min with a constant discharge rate of 4.613 m3/min, and the recovery period is also recorded as 800 min after the pump is turned off. The last record of pumping data is 0.567 m at 800 min. USDI [7] estimated the transmissivity as 2.982 m2/min with the Theis recovery method (for S=S′ and u′≤0.01) and the storage coefficient as 0.07 according to (1) (for sp-s′=0.533 m and Δsp-s′=0.302 m). Figure 1(a) explicitly shows a difficulty at fitting a straight-line to the data on a semilogarithmic graph plot between residual drawdown and t/t′. By the application of (5) and (7) for S/S′=1 to the observed recovery data except for the data at t′=0, 540, and 600 min and after fitting a straight-line, a relationship similar to (4) is obtained (Figure 1(b)). From the straight-line parameters, the transmissivity and the storage coefficients are calculated as 2.8596 m2/min and 0.0661, respectively. For u′=0.2, the recovery time, t′, is determined as 26.8 min from these estimations. For this reason, the data after 26.8 min is reconsidered, and the transmissivity and the storage coefficients are recalculated as 2.8594 m2/min and 0.0666 (for a=0.12838 and b=0.694922), respectively. These parameters yield 0.568 m as a close value to 0.567 m, which is the last drawdown at t=800 min during the pumping period. According to USDI’s [7] parameter values, the last drawdown is 0.544 m, which is far away from 0.567 m. Zheng et al. [18] method (T=2.855 m2/min and S=0.0696) which uses a straight-line similar to the proposed methodology produces 0.563 m. The reason of the difference between the values of Zheng et al. [18] and the methodology of this paper may be due to the lack of S/S′t in the right-side of (3b). Figure 2 shows the measured and simulated residual drawdowns versus time.
Straight-line fit to recovery data: (a) Theis recovery method and (b) proposed method.
Observed and simulated residual drawdowns.
The second set of recovery data is produced synthetically for Q=3 m3/min, r=50 m, tp=70 min, sr,tp=0.396 m, T=1 m2/min, S=0.0055, and S′=0.005. Table 1 shows the residual drawdowns for this data. Figure 3(a) presents a nonlinear relation between τ and W for S=S′, while Figure 3(b) shows a linear relation for S/S′=1.1. Figure 3(a) implies that the rate of S/S′ has a big effect on the late-time residual drawdowns during the recovery period. From the straight-line parameters at Figure 3(b), the transmissivity and the storage coefficients during the recovery and pumping periods are calculated as 1.0 m2/min from (6), 0.0048 from (8), and 0.0053 from S/S′=1.1, respectively. For u′=0.2, the recovery time, t′, is determined as 14.99 min from these estimations (T and S′). For this reason, the data after 14.99 min is reconsidered, and the transmissivity and storage coefficients are found as 1.0 m2/min, 0.00494, and 0.00543 (for a=0.238856 and b=0.737662), respectively. The errors in the obtained storage coefficients with respect to the other storage coefficients are due to the rounded recovery values, which are calculated from (2).
Synthetic residual drawdowns.
t′ (min)
s(r,t′) (m)
0
0.593
5
0.506
10
0.415
15
0.354
20
0.309
25
0.276
30
0.249
40
0.208
50
0.178
70
0.138
90
0.111
120
0.085
150
0.068
210
0.045
270
0.032
330
0.023
Plots of W versus τ of the synthetic data: (a) for S=S′ and (b) straight-line fit to recovery data for S/S′=1.1 and fitting parameters.
4. Conclusion
An effective method has been proposed for decision making to the transmissivity and storage coefficients estimations from the residual drawdowns during the recovery period. The methodology is valid for the confined aquifers and considers the expanding series of Theis well function with the first three terms and the maximum dimensionless time as u′≤0.2. This approach considers a lot of residual drawdown data than the classical Theis recovery method. The procedure depends on a straight-line through the field data and calculates rather easily the aquifer parameters without the pumping data (if it is unavailable). Validity of the procedure is presented by considering actual field data, while the Theis recovery method has some difficulties at fitting a straight-line to a given field data. This method of this paper is very effective for the aquifer parameters estimation and it can be reliably applied to the residual data at an observation well in the extensive and homogeneous confined aquifers with different storage coefficients during the pumping and recovery periods (S≠S′).
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
TheisC. V.The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground-water storage193516251952410.1029/tr016i002p00519ŞenZ.1995New York, NY, USACRC PressBatuV.1998New York, NY, USAJohn Wiley & SonsCooperH. H.JacobC. E.A generalized graphical method for evaluating formation constants and summarizing well-field history194627452653410.1029/tr027i004p00526JacobC. E.The recovery method for determining the coefficient of transmissibility1963283292USGS Water Supply Paper 1536-IBruinJ.HudsonH. E.1955Champaign, Ill, USAIllinois State Water SurveyReport of Investigation no. 25USDI1981Washington, DC, USAUS Department of Interior, Bureau of ReclamationCaseC. M.PidcoeW. W.FenskeP. R.Theis equation analysis of residual drawdown data19741061253125610.1029/wr010i006p012532-s2.0-0016326694AgarwalR. G.A new method to account for producing time effects when drawdown type curves are used to analyze pressure buildup and other test dataProceedings of the 55th Annual Fall Technical Conference and Exhibition of the Society of Petroleum Engineers1980SPE Paper 9289RameyH. J.DocT. W.SchwarzW. J.A drawdown and build-up type curve for interference testingProceedings of the 3rd Invitational Well-Testing Symposium1980Berkeley, Calif, USA130134MishraG. C.ChachadiA. G.Analysis of flow to a large-diameter well during the recovery period198523564665110.1111/j.1745-6584.1985.tb01513.x2-s2.0-0021825850ŞenZ.Drawdown distribution during recovery around a large diameter well19912242532642-s2.0-0026359356YehH.-D.WangC.-T.A semianalytical solution for residual drawdown at a finite diameter well in a confined aquifer201349496697210.1111/jawr.120472-s2.0-84881109441GoodeD. J.Composite recovery type curves in normalized time from Theis' exact solution199735467267810.1111/j.1745-6584.1997.tb00133.x2-s2.0-0031193668BallukrayaP. N.SharmaK. K.Estimation of storativity from recovery data199129449549810.1111/j.1745-6584.1991.tb00540.x2-s2.0-0026189456BantonO.BangoyL. M.A new method to determine storage coefficient from pumping test recovery data199634577277710.1111/j.1745-6584.1996.tb02069.x2-s2.0-0030243028SinghS. K.Storage coefficient and transmissivity from residual drawdowns2003129863764410.1061/(ASCE)0733-9429(2003)129:8(637)2-s2.0-0042322285ZhengL.GuoJ.-Q.LeiY.An improved straight-line fitting method for analyzing pumping test recovery data200543693994210.1111/j.1745-6584.2005.00094.x2-s2.0-27944468177SinghS. K.Identification of aquifer parameters from residual drawdowns: an optimization approach20065161139114810.1623/hysj.51.6.11392-s2.0-33845316910SamaniN.PasandiM.BarryD. A.Characterizing a heterogeneous aquifer by derivative analysis of pumping and recovery test data200612941KambhammettuB. V. N. P.KingJ. P.Estimation of aquifer parameters from residual drawdowns2010163736136510.1680/wama.2010.163.7.3612-s2.0-79958214916AshjariJ.Determination of storage coefficients during pumping and recovery201351112212710.1111/j.1745-6584.2012.00917.x2-s2.0-84871918413