An aquifer consisting of a skin zone and a formation zone is considered as a two-zone aquifer. Existing solutions for the problem of constant-flux pumping in a two-zone confined aquifer involve laborious calculation. This study develops a new approximate solution for the problem based on a mathematical model describing steady-state radial and vertical flows in a two-zone aquifer. Hydraulic parameters in these two zones can be different but are assumed homogeneous in each zone. A partially penetrating well may be treated as the Neumann condition with a known flux along the screened part and zero flux along the unscreened part. The aquifer domain is finite with an outer circle boundary treated as the Dirichlet condition. The steady-state drawdown solution of the model is derived by the finite Fourier cosine transform. Then, an approximate transient solution is developed by replacing the radius of the aquifer domain in the steady-state solution with an analytical expression for a dimensionless time-dependent radius of influence. The approximate solution is capable of predicting good temporal drawdown distributions over the whole pumping period except at the early stage. A quantitative criterion for the validity of neglecting the vertical flow due to a partially penetrating well is also provided. Conventional models considering radial flow without the vertical component for the constant-flux pumping have good accuracy if satisfying the criterion.
The constant-flux pumping (CFP) test is a widely used well test for characterizing the aquifer properties such as transmissivity and storage coefficient. The test is performed with a constant pumping rate at a fully or partially penetrating well in either a confined or unconfined aquifer. Existing analytical solutions for the CFP in a homogenous confined aquifer are briefly reviewed herein. Theis (1935) was the first article in the groundwater literature to present an analytical solution for aquifer drawdown due to pumping in a fully penetrating well with an infinitesimal radius. Carslaw and Jaeger (1959) presented analytical solutions for the three kinds of heat conduction problems which can be analogous to the CFP problems including the aquifers of the infinite domain with a finite-radius well, finite domain with a finite-radius well, and finite domain with an infinitesimal-radius well. Hantush (1962) developed an analytical solution of drawdown induced by a partially penetrating well for the CFP. Papadopulos and Cooper (1967) obtained an analytical solution of drawdown considering the effects of well radius and wellbore storage. They provided a quantitative criterion of time for neglecting the effects. The criterion will be stated in the next section. Chen (1984) derived an analytical solution for drawdown in a circular aquifer with the Dirichlet boundary condition of zero drawdown and provided a quantitative criterion describing the beginning time of the boundary effect on the drawdown. Yang et al. (2006) developed an analytical solution describing aquifer drawdown due to a partially penetrating well with a finite radius. The effect of partial penetration on temporal drawdown distributions was discussed. Wang and Yeh (2008) provided a quantitative criterion for the beginning time of the boundary effect on drawdown induced by the CFP and constant-head pumping. Yeh and Chang (2013) provided a comprehensive review on analytical solutions for the CFP in unconfined and multilayered aquifer systems.
Drilling an aquifer to install a well may decrease or increase the permeability of the formation around the wellbore. The perturbed formation, called skin zone, extends from a few millimeters to several meters. A positive skin zone means that its permeability is lower than the original formation. On the other hand, a negative skin zone is of a higher permeability than the original formation. Existing solutions accounting for the CFP in a two-zone confined aquifer consisting of the skin zone and formation zone are reviewed. Novakowski (1989) developed a semi-analytical solution of drawdown with the wellbore storage effect and investigated the effect of an infinitesimally thin skin on temporal drawdown curves. Hemker (1999) proposed an analytical–numerical solution describing pumping drawdown in a multilayered aquifer system where the radial flow was analytically treated and the vertical one was handled by a finite difference method. The flux along the well screen was non-uniform through an infinitesimal thin skin, and the flow was subject to the wellbore storage effect. Kabala and El-Sayegh (2002) presented a semi-analytical solution for the transient flowmeter test in a multilayered aquifer system where the radial flow was considered in each layer assuming no vertical flow component and uniform flux along the well screen. Predictions from the solution were compared with those from a numerical solution which relaxes those two assumptions. Yeh et al. (2003) obtained an analytical solution for pumping drawdown induced by a finite-radius well in a two-zone confined aquifer and discussed the error caused by neglecting the well radius. Chen and Chang (2006) developed a semi-analytical solution for the CFP on the basis of the Gram–Schmidt method to deal with the non-uniform skin effect represented by an arbitrary piecewise function of elevation. They indicated that flow near a pumping well is three dimensional due to the effect and away from the well is radial. Perina and Lee (2006) proposed a general well function for transient flow toward a partially penetrating well considering the wellbore storage effect and non-uniform flux between the screen and skin zone in a confined, unconfined, or leaky aquifer. Chiu et al. (2007) developed a semi-analytical solution for the CFP at a partially penetrating well in a two-zone confined aquifer. They indicated that the influence of the partial penetration on drawdown is more significant for a negative skin zone than a positive one. C. T. Wang et al. (2012) provided an analytical solution of drawdown for the CFP in a two-zone confined aquifer of finite extent with an outer boundary under the Dirichlet condition of zero drawdown. They also derived a large-time drawdown solution which reduces to the Thiem solution in the absence of the skin zone. X. Wang et al. (2012) presented a finite layer method (FLM) based on Galerkin's technique for simulating radial and vertical flows toward a partially penetrating well in a multilayered aquifer system. The FLM was verified by an analytical solution and finite difference solution.
It is informative to classify the above solutions into two groups, i.e., homogeneous aquifer and two-zone aquifer systems in Table 1. The solutions in each group are categorized according to the well penetration, well radius, and wellbore storage.
At present, a time-domain analytical solution of drawdown for flow induced by the CFP at a finite-radius partially penetrating well in a two-zone confined aquifer has not been developed. The Laplace-domain result of the above-mentioned problem was presented by Chiu et al. (2007) resorting to a numerical inversion scheme called the Crump method. The application of their solution may therefore be inconvenient for those who are not familiar with numerical approaches. The purpose of this note is to develop a new approximate transient solution for the problem in a way similar to the work of Yang et al. (2014). A mathematical model for steady-state flow due to a partially penetrating well in a finite-extent two-zone confined aquifer is built. The flow equations describing spatial drawdowns in the skin and formation zones are employed. The outer boundary of the aquifer is specified as the Dirichlet condition of zero drawdown. The well is treated as the Neumann condition with a constant flux for the screened part and zero flux for the unscreened part. The steady-state solution of the model for drawdown is derived by the method of finite Fourier cosine transform. The approximate transient solution of drawdown is then obtained on the basis of the steady-state solution and a time-dependent radius of influence. The transient solution is in terms of simple series with advantages of fast convergence, simplicity, and good accuracy from a practical viewpoint. It can be used as a convenient tool to estimate temporal and spatial drawdown distributions for the constant-flux pumping and to explore physical insights into the flow behavior affected by hydrogeological properties and aquifer configuration. The accuracy of the solution is investigated in comparison with the Chiu et al. (2007) solution. In addition, the condition of neglecting the effect of the vertical flow on temporal drawdown distributions is investigated.
Categorization of the solutions for the constant-flux pumping in confined aquifers.
A schematic diagram of the constant-flux pumping at a partially penetrating well in a cylinder two-zone confined aquifer with the Dirichlet boundary. (The symbols of the variables are defined in Table 2.)
This section introduces a new mathematical model for steady-state flow due
to the CFP at a finite-radius partially penetrating well in a radial
two-zone confined aquifer. The symbols representing variables and parameters
for the model are listed in Table A1. The hydraulic parameters in the two
zones are different but in each zone are assumed homogeneous. The outer
boundary is considered to be under the Dirichlet condition of
The governing equations describing steady-state dimensionless drawdown
distributions in the skin and formation zones are expressed, respectively, as
A new solution derived by the application of the finite Fourier cosine
transform to the model can be written as
The inverse Laplace transform to Chiu et al. (2007) semi-analytical solution
of drawdown leads to a time-domain result for the CFP in a two-zone aquifer
system; however, the resultant solution involves laborious calculations. We
therefore develop an approximate transient solution of drawdown for the CFP
problem. The idea originated from the concept of a time-dependent diffusion
layer for the solution of the diffusion equation in the field of
electrochemistry (Fang et al., 2009). The approximate transient solution is
obtained by replacing the
Predicted drawdowns by the Chiu et al. (2007) solution and
the approximate solution, Eqs. (16) and (17), with
When
When
On the basis of the comparison of predictions from the approximate solution
and Chiu et al. (2007) Laplace-domain solution, we have concluded that the
accuracy of the present solution depends only on dimensionless time
The vertical flow induced by well partial penetration is strongly dependent
on both dimensionless lumped parameters
Temporal drawdown distributions predicted by the
approximate solution, Eq. (17), with
This study presents an approximate drawdown solution, Eqs. (16) and (17), in
terms of a simple series for the CFP at a partially penetrating well in a
radial two-zone confined aquifer. The solution is developed on the basis of
the steady-state drawdown solution with an outer boundary represented by the
time-dependent radius of influence. The comparison with the Chiu et al. (2007)
solution reveals that the approximate solution gives accurate
temporal drawdown distributions in an observation well over the entire
pumping period except at an early time when the dynamic radius of influence
reaches the well (i.e.,
The finite Fourier cosine transform is defined, in our notation, as
Summary of symbols used in the text and their definitions.
Research leading to this paper has been partially supported by the grants from the Taiwan Ministry of Science and Technology under the contract numbers NSC 102-2221-E-009-072-MY2 and MOST 103-2221-E-009-156. Edited by: S. Attinger