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 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 . 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 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 () 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 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 () 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 time-recovery 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 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 (), 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 . 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 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 , 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 straight-line 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 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, , should be determined especially for the first data values. If is greater than 0.2, then the straight-line should be rearranged. Furthermore, the ratio of 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 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 straight-line 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 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 , 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 m2/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 ( m2/min and ) 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 in the right-side 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 m3/min, m, min, m, m2/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 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 , 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 m2/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 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 ().
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.