doi:10.1016/S0022-1694(98)00223-6
Copyright © 1998 Elsevier Science B.V. All rights reserved.
Hydraulics of a partially penetrating well: solution to a mixed-type boundary value problem via dual integral equations
G. Cassiani1 and Z. J. Kabala*
Department of Civil and Environmental Engineering, Duke University, Durham NC 27708, USA
Received 24 April 1997;
accepted 24 August 1998.
Available online 18 January 1999.
References and further reading may be available for this article. To view references and further reading you must
purchase this article.
Abstract
New semi-analytic solutions are obtained for well response to the pumping test and slug test performed on a partially penetrating well. The solutions account not only for the wellbore storage, infinitesimal skin, and aquifer anisotropy, but also for the mixed-type boundary condition at the well face, which is novel. The solutions are obtained via the method of dual integral equations (DE). The new solutions are computationally robust and efficient, about one to two orders of magnitude faster than the corresponding finite difference solutions. Existing approximate solutions obtained with flux–flux discontinuous boundary conditions are compared to our DE solutions. The accuracy of the approximate solutions appears to be adequate for slender well screens. Our DE solution is computationally more efficient than the approximate solutions. In the range where the approximate solutions are less accurate the DE solution is about an order of magnitude faster. More important, the new solutions provide the correct distribution of the point flux (local velocity) along the well screen, unlike all existing solutions. This feature is essential in cases where vertical variations of hydraulic conductivity are sought (e.g. in flowmeter tests), and for tracer tests.
Author Keywords: Dual integral equations; Mixed-type boundary value problem; Partially penetrating well
Fig. 1.A partially penetrating well in a confined aquifer of semi-infinite vertical extent.
Fig. 2.The point flux for a pumped well with no wellbore storage. α is the dimensionless time. The figure shows an example for dimensionless screen length ξ=50. The uniform point flux used in classical approximations is shown for comparison.
Fig. 3.The dimensionless drawdown using the DE approach for pumped wells in confined aquifers without wellbore storage. The solution for fully penetrating wells with finite well radius, but no wellbore storage (Hantush, 1964), and the (Theis, 1935) solution are shown for comparison.
Fig. 4.The dimensionless drawdown using the DE approach with and without wellbore storage. The solution for fully penetrating wells in confined aquifers with finite well radius, but no wellbore storage (Hantush, 1964), and the Papadopulos and Cooper, 1967 solution (f.p.) for fully penetrating wells with wellbore storage, are shown for comparison (f.p.=fully penetrating st.=storage).
Fig. 5.The dimensionless drawdown using the DE approach for slug tests. The solution for fully penetrating wells by Cooper et al., 1967 is shown for reference.
Fig. 6.The dimensionless drawdown using the DE approach, compared with the approximate solution for constant point flux at the well face, and infinitesimal well radius (Hantush, 1961a). The solution for fully penetrating wells with finite well radius, but no wellbore storage (Hantush, 1964), and the Theis, 1935 solution are shown for comparison.
Fig. 7.The comparison between the DE approach and the Dougherty and Babu, 1984 (D and B) approach for pumped tests. The agreement is good. The solution for fully penetrating wells with finite well radius, but no wellbore storage (Hantush, 1964), and the Theis, 1935 solution are shown for comparison.
Fig. 8.The comparison between the DE approach and the Dougherty and Babu, 1984 (D and B) approach for slug tests. The agreement is excellent. The solution for fully penetrating wells by Cooper et al., 1967 is shown for reference. For all curves the wellbore storage coefficient is αc=10-5.