Research Article  Open Access
A StraightLine Method for Analyzing Residual Drawdowns at an Observation Well
Abstract
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 straightline through the field data and it helps to calculate the parameters quickly without any need for longterm 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 . The method can be applied reliably to extensive and homogeneous confined aquifers resulting in different storage coefficients during the pumping and recovery periods . 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 latetime 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 straightline through the residual drawdowns () versus plot ( indicates the total time since pump start while is the time since pump shut). This plot on the semilogarithmic graph paper passes below the origin (, ) giving the value of 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 timerecovery drawdown graph to find the time of zero recovery. The method depends on the extension of the timedrawdown pumping test data, which can also be applied to the timerecovery graph. Later, USDI [7] gave an alternative method for determining the storage coefficient () as follows: where is the transmissivity, is pumping period drawdown projected to time at any radial distance, , is residual drawdown at time , is recovery at time , and is slope of the timerecovery graph. In this formulation corresponds to the antilogarithm as .
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 largediameter 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 largediameter 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 wellfunction solution. These type curves depend on the dimensionless duration of pumping. Ballukraya and Sharma [15] proposed an approach derived from the CooperJacob 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 (), 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 straightline method based on the CooperJacob 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 . Singh [19] presented an optimization method based on nonlinear leastsquares 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 drawdownderivative 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 LevenbergMarquardt 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 by using a modified version of Banton and Bangoy [15] method. This method is basically fitting a straightline to a plot of residual drawdown versus square of radial distance at the same time.
In this study, another straightline 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 , 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, .
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:where is the residual drawdown at any distance and at any recovery time , is the constant rate (discharge) towards the pumping well during the pumping and recovery periods, is the transmissivity, is the dimensionless time factor for the pumping period, is another dimensionless time factor for the recovery period, and are the storage coefficients of aquifer during the pumping and recovery periods, , and is the time of pumping. This expression can be considered after the first three terms of the exponential function series as follows:
The error involved in adopting (3b) instead of (2) is less than 1% for . Theis [1] proposed the first two terms of the series in (2) by considering that in order to estimate the aquifer transmissivity only. Theis approach is valid for . Hence, (3b) considers more recovery data than Theis method. Equation (3b) can be rewritten to estimate the aquifer parameters asEquation (4) presents a straightline between and . 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 straightline 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, , should be determined especially for the first data values. If is greater than 0.2, then the straightline should be rearranged. Furthermore, the ratio of may be easily investigated with various straightlines.
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 m^{3}/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 m^{2}/min with the Theis recovery method (for and ) and the storage coefficient as 0.07 according to (1) (for m and m). Figure 1(a) explicitly shows a difficulty at fitting a straightline to the data on a semilogarithmic graph plot between residual drawdown and . By the application of (5) and (7) for to the observed recovery data except for the data at , 540, and 600 min and after fitting a straightline, a relationship similar to (4) is obtained (Figure 1(b)). From the straightline parameters, the transmissivity and the storage coefficients are calculated as 2.8596 m^{2}/min and 0.0661, respectively. For , the recovery time, , 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 m^{2}/min and 0.0666 (for and ), respectively. These parameters yield 0.568 m as a close value to 0.567 m, which is the last drawdown at 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 ( m^{2}/min and ) which uses a straightline 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 in the rightside of (3b). Figure 2 shows the measured and simulated residual drawdowns versus time.
(a)
(b)
The second set of recovery data is produced synthetically for m^{3}/min, m, min, m, m^{2}/min, , and . Table 1 shows the residual drawdowns for this data. Figure 3(a) presents a nonlinear relation between and for , while Figure 3(b) shows a linear relation for . Figure 3(a) implies that the rate of has a big effect on the latetime residual drawdowns during the recovery period. From the straightline parameters at Figure 3(b), the transmissivity and the storage coefficients during the recovery and pumping periods are calculated as 1.0 m^{2}/min from (6), 0.0048 from (8), and 0.0053 from , respectively. For , the recovery time, , is determined as 14.99 min from these estimations ( and ). For this reason, the data after 14.99 min is reconsidered, and the transmissivity and storage coefficients are found as 1.0 m^{2}/min, 0.00494, and 0.00543 (for and ), 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).

(a)
(b)
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 . This approach considers a lot of residual drawdown data than the classical Theis recovery method. The procedure depends on a straightline 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 straightline 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 ().
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
References
 C. V. Theis, “The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage,” Transactions, American Geophysical Union, vol. 16, no. 2, pp. 519–524, 1935. View at: Publisher Site  Google Scholar
 Z. Şen, Applied Hydrogeology for Scientists and Engineers, CRC Press, New York, NY, USA, 1995.
 V. Batu, Aquifer Hydraulics: A Comprehensive Guide to Hydrogeologic Data Analysis, John Wiley & Sons, New York, NY, USA, 1998.
 H. H. Cooper and C. E. Jacob, “A generalized graphical method for evaluating formation constants and summarizing wellfield history,” Transactions, American Geophysical Union, vol. 27, no. 4, pp. 526–534, 1946. View at: Publisher Site  Google Scholar
 C. E. Jacob, “The recovery method for determining the coefficient of transmissibility,” in Methods of Determining Permeability, Transmissibility and Drawdown, USGS Water Supply Paper 1536I, pp. 283–292, 1963. View at: Google Scholar
 J. Bruin and H. E. Hudson, Selected Methods for Pumping Test Analysis, Report of Investigation no. 25, Illinois State Water Survey, Champaign, Ill, USA, 1955.
 USDI, Ground Water Manual, US Department of Interior, Bureau of Reclamation, Washington, DC, USA, 1981.
 C. M. Case, W. W. Pidcoe, and P. R. Fenske, “Theis equation analysis of residual drawdown data,” Water Resources Research, vol. 10, no. 6, pp. 1253–1256, 1974. View at: Publisher Site  Google Scholar
 R. G. Agarwal, “A new method to account for producing time effects when drawdown type curves are used to analyze pressure buildup and other test data,” in Proceedings of the 55th Annual Fall Technical Conference and Exhibition of the Society of Petroleum Engineers, SPE Paper 9289, 1980. View at: Google Scholar
 H. J. Ramey, “A drawdown and buildup type curve for interference testing,” in Proceedings of the 3rd Invitational WellTesting Symposium, T. W. Doc and W. J. Schwarz, Eds., pp. 130–134, Berkeley, Calif, USA, 1980. View at: Google Scholar
 G. C. Mishra and A. G. Chachadi, “Analysis of flow to a largediameter well during the recovery period,” Ground Water, vol. 23, no. 5, pp. 646–651, 1985. View at: Publisher Site  Google Scholar
 Z. Şen, “Drawdown distribution during recovery around a large diameter well,” Nordic Hydrology, vol. 22, no. 4, pp. 253–264, 1991. View at: Google Scholar
 H.D. Yeh and C.T. Wang, “A semianalytical solution for residual drawdown at a finite diameter well in a confined aquifer,” Journal of the American Water Resources Association, vol. 49, no. 4, pp. 966–972, 2013. View at: Publisher Site  Google Scholar
 D. J. Goode, “Composite recovery type curves in normalized time from Theis' exact solution,” Ground Water, vol. 35, no. 4, pp. 672–678, 1997. View at: Publisher Site  Google Scholar
 P. N. Ballukraya and K. K. Sharma, “Estimation of storativity from recovery data,” Ground Water, vol. 29, no. 4, pp. 495–498, 1991. View at: Publisher Site  Google Scholar
 O. Banton and L. M. Bangoy, “A new method to determine storage coefficient from pumping test recovery data,” Ground Water, vol. 34, no. 5, pp. 772–777, 1996. View at: Publisher Site  Google Scholar
 S. K. Singh, “Storage coefficient and transmissivity from residual drawdowns,” Journal of Hydraulic Engineering, vol. 129, no. 8, pp. 637–644, 2003. View at: Publisher Site  Google Scholar
 L. Zheng, J.Q. Guo, and Y. Lei, “An improved straightline fitting method for analyzing pumping test recovery data,” Ground Water, vol. 43, no. 6, pp. 939–942, 2005. View at: Publisher Site  Google Scholar
 S. K. Singh, “Identification of aquifer parameters from residual drawdowns: an optimization approach,” Hydrological Sciences Journal, vol. 51, no. 6, pp. 1139–1148, 2006. View at: Publisher Site  Google Scholar
 N. Samani, M. Pasandi, and D. A. Barry, “Characterizing a heterogeneous aquifer by derivative analysis of pumping and recovery test data,” Journal of Geological Society of Iran, vol. 1, pp. 29–41, 2006. View at: Google Scholar
 B. V. N. P. Kambhammettu and J. P. King, “Estimation of aquifer parameters from residual drawdowns,” Proceedings of the Institution of Civil Engineers: Water Management, vol. 163, no. 7, pp. 361–365, 2010. View at: Publisher Site  Google Scholar
 J. Ashjari, “Determination of storage coefficients during pumping and recovery,” GroundWater, vol. 51, no. 1, pp. 122–127, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Mesut Çimen. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.