Next Article in Journal
Microphysical Characteristics of Winter Precipitation in Eastern China from 2014 to 2019
Next Article in Special Issue
Water Flow, Solute and Heat Transfer in Groundwater
Previous Article in Journal
SPH Simulation of Interior and Exterior Flow Field Characteristics of Porous Media
Previous Article in Special Issue
The Impact of Capillary Trapping of Air on Satiated Hydraulic Conductivity of Sands Interpreted by X-ray Microtomography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seepage Characteristics of a Single Ascending Relief Well Dewatering an Overlying Aquifer

1
Henan Province Key Laboratory of Rock and Soil Mechanics and Structural Engineering, North China University of Water Resources and Electric Power, 136 Jinshui East Rd, Zhengzhou 450045, China
2
Energy Geosciences Division, Earth and Environmental Sciences Area, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA 94720, USA
3
Hydrology and Water Resources Bureau of Henan Province, Zhengzhou 450003, China
*
Authors to whom correspondence should be addressed.
Water 2020, 12(3), 919; https://doi.org/10.3390/w12030919
Submission received: 10 February 2020 / Revised: 16 March 2020 / Accepted: 17 March 2020 / Published: 24 March 2020
(This article belongs to the Special Issue Water Flow, Solute and Heat Transfer in Groundwater)

Abstract

:
The application of groundwater relief, i.e., dewatering, ascending wells, drilled upward from the mining tunnel into the overlying aquifer, is common in underground mining engineering. In this study, the seepage characteristics of single ascending partially and fully penetrating relief wells are investigated using a series of laboratory sand-tank experiments and numerical simulations. The seepage characteristics of ascending wells dewatering an overlying aquifer are different from those of conventional pumping wells descending from the ground surface into the underlying aquifer, because of the pronounced influence of the seepage face boundary condition along the seepage boundary of the ascending dewatering well. The seepage face of the ascending well is formed as the well casing remains open and water is discharged under the action of gravity through the well casing. The results of laboratory sand-tank experiments and modeling show that when the degree of penetration of an ascending relief well does not exceed a critical value, the effect of the seepage face cannot be ignored. In particular, the seepage flux increases as the degree of penetration increases following an exponential function, and the relationship between the seepage flux and the well radius can be described using a power law function. The results of numerical simulations are used to develop a series of type curves to evaluate the effects of the critical degree of penetration for different well radii and different aquifer water levels. Modified versions of the Dupuit and Dupuit–Thiem formulae for a single ascending partially well for the degree of penetration less than the critical one for the unconfined, confined, and confined-unconfined aquifers are developed.

1. Introduction

Construction of many underground facilities and tunnels, deep foundation pits, and underground coal mines, which may be affected by underground water intrusion, require dewatering of an overlying aquifers [1,2,3,4,5,6,7,8,9,10,11]. A common method to reduce the groundwater level to prevent inundation of underground excavations is the use of water pumping wells, which are drilled from the land surface into the aquifer [5,7,8]. Another method applied in the mining industry for dewatering and depressurization of overlying aquifers is to drill dewatering wells from the underground mines’ tunnels upward into the overlying aquifer [1,11,12,13,14,15,16,17,18,19].
In this paper, the dewatering wells drilled upward into the overlying aquifer are called “ascending relief wells” (ARWs), and they are used to partially or completely dewater an aquifer above an active or abandoned underground mine. The ARWs can also be used in combination with pumping wells. For example, a comprehensive drainage from pumping wells and the panels was used in the 8th mining area of Taiping Coal Mine, China [17]. The ascending wells were drilled below the mining face to test hydrological parameters in many underground coal mines [11,15]. The drainage water discharged under gravity from the overlying aquifer through ARWs is collected in the underground tunnel and then pumped out through the drainage system. Schematic diagrams of the partially penetrating ARWs in an unconfined aquifer, a confined aquifer, and a confined-unconfined aquifer are shown in Figure 1.
Contrary to the conventional pumping wells, the casing of the ARW remains open to the ambient air with no standing water. The seepage characteristics of the ARWs are different from those of conventional pumping wells, because of downward water flow by the action of gravity into the underground tunnel. The physics of seepage through the walls of the ARWs is generally similar to that of seepage through the walls of the open tunnels. In particular, because of the formation of the seepage face at the tunnel walls [20,21,22,23,24], boundary conditions for numerical simulations of seepage into the tunnels are based on assigning the atmospheric pressure at the tunnel walls and a constant hydraulic head boundary condition along the tunnel perimeter [25,26,27,28]. The literature review shows that there is no generalized theory to assess the performance of ARWs [11,29,30], and determination of seepage characteristics of ARWs is based on using methods developed for conventional pumping wells [4,31,32,33,34,35,36,37,38,39].
The overall objectives of this paper are (1) to study the evolution of the seepage characteristics of single ascending partially and fully penetrating relief wells, and (2) to determine the critical degree of penetration of a partially penetrating ARW needed to obtain a maximum outflow under different far-field water level boundary conditions, (3) to assess the applicability and to modify the Dupuit and Dupuit–Thiem formulae for a single ascending partially well, in an unconfined, homogeneous and isotropic aquifer. The Dupuit and Dupuit–Thiem formulae were considered as a basis to calculate seepage fluxes of the ARW, and laboratory sand-tank experiments and numerical simulations were also conducted.
The following simplifying assumptions were taken into consideration in this study: (1) the aquifer is homogeneous, isotropic, and with no regional flow; (2) there is no water leakage through the top and bottom of the aquifer, i.e., through the underlying and overlying aquicludes of the confined aquifer; (3) the ARW is vertical with a screen along its entire length of the ARW penetration into the overlying aquifer, and no water column inside the ARW, allowing for free drainage (under the gravity), through the open well bottom, into the underlying tunnel, and there is no enforced pumping; (4) the atmospheric pressure boundary condition at the well-aquifer interface; and (5) groundwater flow in the aquifer toward the well is axisymmetric, and can be described by the Darcy law.

2. Methods of Analytical, Experimental, and Modeling Investigations

2.1. Dupuit and Dupuit–Thiem Formulae

Analysis of pumping wells drilled from the surface, in general, relies on the Dupuit assumption of a constant hydraulic head along the vertical profile of the aquifer. This assumption makes possible to use a depth-integrated equation for the flux evaluation into the well. Based on the original analytical solution for unconfined flow by Dupuit [40], Thiem [41] was likely the first to estimate aquifer parameter from pumping tests in a confined aquifer. Boulton [42] developed the first transient well test solution to analyze the unconfined aquifer. Streltsova [43] proposed the solution to consider unsteady radial flow in an unconfined aquifer. Neuman [44,45] developed solutions considering both confined storage and delayed yield from the unconfined aquifer. A variety of new analytical approaches were applied to deal with well hydraulics models involving mixed and complicated boundaries, which gained more accurate solutions for particular conditions [7,36,39,46,47,48,49,50,51]. Mishra and Kuhlman [52] have recently summarized the application of the Dupuit theory for an unconfined aquifer from Dupuit to the present. Charnyi [53] showed that the Dupuit formula for an unconfined aquifer is valid not only for a hydraulic approximation, i.e., assuming that the groundwater velocity is independent of the height above the aquitard, but also for rigorous hydrodynamic calculations. Although, the Dupuit formula is not accurate for calculations of the phreatic line for r < H, it provides accurate calculations of the seepage flux [53]. Shercliff [54] conducted one—and two-dimensional modeling of steady seepage flow in unconfined aquifers, and provided a proof of Charnyi’s result that one—and two-dimensional theory yield the same value for the flow rate in a horizontal aquifer or porous bed between vertical ends, and showed the extent to which it can be generalized to non-uniform or anisotropic media.
The forms of the Dupuit and Dupuit–Thiem formulae are summarized in Table 1, and were used in this paper as a basis to estimate seepage fluxes into ARWs, and were then modified based on the results of laboratory sand-tank experiments and modeling (see Section 3).

2.2. Laboratory Experiments

2.2.1. Seepage Sand-Tank and Boundary Conditions

Laboratory sand-tank axisymmetric radial water flow experiments were conducted to investigate the seepage to the ARW in an unconfined aquifer. The outer boundary head was at the elevation of 60 cm and kept constant for the entire test. The water entering the ARW was discharged freely under gravity, and the piezometric head inside the ARW was zero. The boundary conditions are shown Figure 2. The sand-tank model is generally a downscaled model of the aquifer [55].
The tank was fabricated from steel plates and acrylic sheets, representing a 60° sector, i.e., a 1/6 portion of a circular flow system for saving experimental materials, with a radius of 210 cm and height of 70 cm. A design schematic and a photograph of the sand-tank assembly are shown in Figure 3. The outer boundary water head was maintained at the elevation of 60 cm above the bottom of the sand-tank, and water was supplied into the sand-tank from a plastic cistern. The outflow from the sand-tank was monitored by an electronic scale. The overflowing water and the well discharge were collected in separate water collection vessels for cyclic utilization. To measure the hydraulic pressure, 96 piezometers with water-level monitoring tubes were installed in six layers at different elevations and distances from the modeled ARW well. The inside diameter of the piezometers was 4 mm, shown in Figure A1 (see Appendix A). The piezometer storage effect is irrelevant because steady state water levels were measured. Water levels in the piezometer tubes were recorded using a camera, which is shown in Figure 3a.

2.2.2. Ascending Relief Well (ARW) Design

Laboratory tests were conducted to investigate the performance of 24 types of ARWs with different radii—4.5 cm, 9.5 cm, 14.5 cm and 19.5 cm, and lengths—2 cm, 5 cm, 10 cm, 20 cm, 30 cm and 60 cm, corresponding to the degrees of penetration 1/30, 1/12, 1/6, 1/3, 1/2 and 1, as shown in Figure A2. To prevent migration of sand particles into the testing wells, 1-mm radius holes with 4 mm pitches were made in the well screens and tops, and the well screens were wrapped up by stainless steel filters with 0.25 mm diameter openings and 0.18 mm diameter steel wire.

2.2.3. Laboratory Testing Protocol

Loading the Sand-Tank

The tank was filled up with preliminary washed river coarse sand, the particle size distribution curve of a sand sample used for the sand-tank experiment is shown in Figure A3. To preclude possible seepage along the sand-tank wall, a thin layer of vaseline was first besmeared along the inner surface of the tank wall. Six 10 cm thick sand layers were loaded into the tank layer by layer (in total 1.38 m3), and 16 piezometers were installed between each two layers–in all 96 piezometers.

Sand Saturation

To saturate the sand, water was injected through two bottom inlet ports, while the water supply level was set at the 60 cm elevation above the tank bottom. During the tank saturation, the rising water level in the tank allowed for air escaping upward from the sand, which reduced the volume of entrapped air [56]. When the water table reached the elevation of 60 cm, water injection was continued for about two hours to allow for expelling the remaining entrapped air from the sand. Then, the bottom water inlets were shut off and kept closed during the dewatering process, and the water was diverted from the plastic cistern to the water pressure-stabilizing tank to maintain the stable upper boundary water level. The estimated saturated hydraulic conductivity of sand was about 0.6 cm/s. During dewatering, to assess the flux, water samples were collected and weighed, simultaneously with camera recording and measurements of piezometric heads in monitoring tubes.

2.3. Numerical Simulations

Numerical simulations of flow in the aquifer and the seepage from the ARW were carried out using a software package, MIDAS GTS NX. This software package is based on finite-element simulations for geotechnical design applications, such as deep foundations, excavations, complex tunnel systems, seepage, consolidation analysis, embankment design, dynamic and slope stability analysis—see http://en.midasuser.com/product/gtsnx_overview.asp. Two types of 3-D numerical simulations were conducted: one was to replicate the sand-tank laboratory experiment (Figure A4a), and the other one was to simulate the extended flow field domain (Figure A4b).
The 3-D extended modeling domain was 1200 m by 1200 m in the plan and 120 m in the vertical direction, with a 60 m homogeneous aquifer and a 60 m aquiclude overlying the aquifer. The ARW was located at the center of the domain, shown in Figure A4b. Simulations were performed to assess water flow to a single partially ARW and a fully penetrating ARW with different radii (0.1 m, 0.2 m, 0.5 m, 1 m, 2 m, 4 m, 6 m and 10 m), the well length (from 1 to 60 m with the 1 interval), and the outer boundary heads (60 m, 70 m, 90 m, 120 m, 150 m, 180 m, 210 m, and 240 m). To simulate the discharge from the ARW under gravity, the piezometric head inside the ARW was set to zero. The aquifer saturated hydraulic conductivity was 0.4 m/d. The fixed water pressure head (Dirichlet-type) boundary conditions were assigned at the outer boundary of the flow domain and the ARW well screen.

3. Results of Laboratory Sand-Tank Tests

To present the results of investigations, all distances and hydraulic heads are given as dimensionless values being normalized by the aquifer thickness M for the confined aquifer or the original water level, H, above the unconfined aquifer bottom. For example, for the confined aquifer: R nor = R x M , H nor = H M , p nor = p M , l nor = l w M , s nor = s M and h nor = h M , where Rx is the radial distance, H is the outer boundary water level, p is the piezometric head, lw is the ARW penetrating length into the aquifer, s is the drawdown, and h is the height above the aquifer bottom. Here, lnor is also the ARW degree of penetration into the aquifer. The well radius rw was normalized as r nor = r w r 1 , using r1 = 1 cm for the lab experiment, and 1 m for numerical modeling. The seepage fluxes Qr1 for the fully penetrating well of radius r1 = 1 cm (for lab experiment) and 1 m (numerical modeling), which were calculated for the case of complete drawdown, were used to determine the dimensionless seepage fluxes given by Q nor = Q Q r 1 . The results are given below using dimensionless parameters.

3.1. Piezometric Head Distribution

Figure 4 demonstrates the piezometric head distribution vs. the distance from the well in the sand tank experiments for the case of the large diameter fully penetrating well, with rnor = 19.5 and lnor = 1. For large values of the well radius and length, the piezometric heads were zero at different heights above the ARW center with a maximum seepage flux. The seepage characteristics in the vicinity and at the ARW are practically the same as those in the case of a fully penetrating pumping well. Figure 5 shows the piezometric head distribution vs. distance from a partially penetrating ARW with a small degree of penetration of lnor = 1/30 and a radius of rnor = 4.5. The figure shows the zero piezometric head inside and at the top of the ARW, and the piezometric heads above the well top being greater than zero, indicating the downward flow. Figure 6 depicts the piezometric head distribution along the height of the sand-tank at the well center Rnor = 0 and at the radial distance Rnor = 1/6. In particular, the piezometric heads at the heights of hnor = 1/6 and hnor = 1/3 above the ARW center were nearly the same.
Figure 7 shows the 2-D plots of the piezometric head distribution (for the case of rnor = 4.5) for different values of the ARW degree of penetration. The experimental results show that the piezometric head inside the ARW was zero, and a piezometric head dropped in the vicinity of the partially penetrating ARW. The water table is determined at the elevation where the piezometric heads is zero. Figure 5 and Figure 6 show that the piezometric heads in the vicinity and above the partially penetrating ARW were greater in the middle and smaller at the top and bottom of the well. With the increase of the radial distance from the ARW, this effect gradually disappeared, and the piezometric heads at different depths increased with depth (see Figure 4 and Figure 5).
The effect of the well length on the water head distribution can be described as follows. For lnor ≥ 1/2, the piezometric head above the well is zero. For lnor ≤ 1/3, the piezometric head above the well is greater than zero. Because experiments were conducted for only six values of ARW length for rnor = 4.5, it can be deduced that the critical length of the ARW penetration lnor ranges from 1/3 to 1/2.
The critical degree of ARW penetration of a confined aquifer ζ c = l c / M is defined as the ratio of the well critical length to the aquifer thickness M. For an unconfined aquifer, ζ c = l c / H , where H is the original water level above the unconfined aquifer bottom.
The influence of the well radius, length and boundary aquifer head on the critical degree of penetration of an ARW will be evaluated using the results of numerical modeling summarized below in Section “Numerical simulation results.”

3.2. Seepage Flux

Figure 8 demonstrates that the seepage flux (Qnor) increases asymptotically with the increase of the degree of penetration lnor, which can be described by an exponential relationship (with R2 > 0.96) given by:
Q nor = Q nor 0 + A 0 × exp ( B 0 × l nor ) ,   for   l nor < ζ c Q nor = Q nor 1 + A 1 × exp ( B 1 × ζ c ) ,   for   l nor ζ c
where Q nor 0 , A0 and B0, Q nor 1 , A1 and B1 are fitting parameters, and ζ c is the critical degree of penetration.
Figure 8 demonstrates that for the small well radii, the seepage flux increased initially rapidly, and then practically stabilized for ζ c > 0.35. When the well radius is large, the seepage flux reaches the maximum value at a relatively small degree of penetration. The critical degree of penetration ζ c decreases with the increase in the well radius, such as ζ c = 0.1 for rnor = 14.5, ζ c = 0.15 for rnor = 9.5, and ζ c = 0.35 for rnor = 4.5. (The dependence of the critical degree of penetration ζ c on the well radius and the boundary head, based on the results of numerical simulations, is described below in Section “Numerical simulation results”). Figure 9 demonstrates a relationship between the seepage flux Qnor and the well radius rnor given by
Q nor = Q nor 2 + A 2 × r nor B 2
with R2 > 0.96 and B2 < 1, where Q nor 2 , A2 and B2 are fitting parameters, rnor is the well radius.

4. Numerical Simulation Results

4.1. Relationships between Seepage Flux, Well Length and Radius

The results of numerical simulations of the relationships between the seepage flux, well length and radius for different values of the boundary hydraulic head are shown in Figure 10 and Figure 11. The relationship between the seepage flux and the degree of penetration can be described using an exponential function given by Equation (1) with R2 > 0.99. The seepage flux increases with the increase in the well radius according to a power law relationship given by Equation (2) with R2 = 1. The results of numerical simulations correspond to those from the laboratory experiments. Figure 10 and Figure 11 also show that the seepage flux increases with the increase in the degree of penetration, and asymptotically reaches the maximum value for l nor ζ c . When the degree of penetration is greater than the critical value of ζ c , the seepage flux reaches the maximum value and remains practically constant as the degree of penetration continues to increase. The critical degree of penetration ζ c decreases with the increase in the well radius, and increases with the increase of the boundary head.

4.2. Seepage Characteristics

The water piezometric head pnor distributions above the well for the case Hnor = 1 and rnor = 1.0 in an unconfined aquifer are shown in Figure 12. For the partially penetrating ARW (for lnor from 1/60 to ½), shown in the subfigures (1)–(7) of Figure 12, the lines of the piezometric head of zero (i.e., indicating the water level) are located above the well screen (at x = 0), which correspond to the degree of penetration being less than the critical value of ζ c . For the partially penetrating wells (for lnor = 2/3 and 5/6), shown in subfigures (8) and (9) and a fully penetrating well, shown in subfigure (10) of Figure 12, the lines of the zero piezometric head are below the top of the well, which correspond to conditions exceeding the critical value of ζ c .
When the well length penetration exceeds the critical length lc, the seepage flux of a partially penetrating ARW is the same as that of a fully penetrating well generating a maximum seepage flux. A similar pattern of the seepage flux was determined for a descending partially penetrating pumping well [57]. The piezometric head distribution characteristics determined using numerical modeling are consistent with the experimental results in the sand-tank described above in Section “Results of Laboratory Sand-Tank Tests.”
Figure 13 depicts the flow net composed of streamlines and equipotential lines for the boundary head of Hnor = 1 and rnor = 1.0. The flow net indicates that the equipotential lines above the ARW are concaved in the vicinity above the ARW. The curvature of the arched shape equipotential lines gradually increases with the well length increase, and when the equipotential line intersects the water table line, the equipotential line bifurcates into two lines, symmetrically distributed around the ARW. When the ARW’s degree of penetration of exceeds the critical value of ζ c , the equipotential lines are symmetrically distributed around the ARW. The streamlines shown in Figure 13 indicate that a single ARW can drain the entire aquifer regardless of the degree of penetration, and even for the degree of the penetration less than ζ c , groundwater flow from the entire aquifer thickness is directed into the ARW.

4.3. ARW’s Critical Degree of Penetration

A summary of the calculations from the results of numerical simulation scenarios is presented in Table A2 and Figure 14. Figure 14 shows that ζ c increases with Hnor for different rnor, according to the relationship given by:
ζ c = ( 0.0488 × r nor + 0.98 ) + ( 0.0273 r nor + 0.1118 ) × ln ( H nor + 0.0569 r nor 1.0214 )
where rnor is the ARW normalized well radius, Hnor = H/M for the case of a confined aquifer, and Hnor = H/H = 1 for an unconfined aquifer.
For large Hnor and a small well radius, for example, when Hnor > 2.0 and rnor < 0.2, even for a fully penetrating well, the piezometric head above the well still exceeds zero, and the maximum seepage flux is not reached, as shown in Figure 14. Figure 15 presents a contour map of ζ c for different Hnor and rnor based on the results of numerical simulations. These results can be used for designing optimal ARWs to gain maximum seepage flux with minimum drilling lengths.

4.4. Model Validation

Validation of the numerical model was conducted based on a comparison of the simulation results, performed using the code MIDAS GTS NX and calculations using the Dupuit formulae for a fully penetrating well in an unconfined aquifer, to the results of laboratory experiments, which are shown in Figure 16. The results of calculations of errors are listed in Table A1. Based on the results of calculations, the absolute relatively errors of the numerical simulations and Dupuit formulae calculations compared to the results of laboratory experiments are <7.65% and <4.49%, respectively.

5. Modified Dupuit and Dupuit–Thiem Formulae for a Single ARW

5.1. Unconfined Aquifer

Table A3 presents the results of simulations of the seepage flux from numerical simulations (Qn) and calculations using the Dupuit formula (Qs) for a single partially penetrating ARW in an unconfined aquifer, shown in Figure 1a, for rnor = 1.0 and the boundary head Hnor = 1. One can see that the seepage flux calculated using the Dupuit equation is smaller than that determined from numerical simulations. For lw < lc, the effective water level drawdown above the ARW can be determined as the value seff = s+lw, and a modified version of the Dupuit equation for the seepage flux Qr can be given by:
Q r = 1.366 K ( 2 H s eff ) s eff lg R r w = 1.366 K ( 2 H s eff ) s eff lg 10 s eff K r w
with parameters identified on Figure 1a.
The seepage flux Q nor r calculated using the modified Dupuit Equation (4) is also given in Table A3, and it is practically the same as that from numerical simulations, as shown in Figure 17a. Note that for lwlc, the ARW seepage flux can be calculated from the Dupuit equation, as for a fully penetrating pumping well.

5.2. Confined Aquifer

For the case of a partially penetrating ARW in a confined aquifer, shown in Figure 1b, the seepage flux can be calculated using a modified version of the Dupuit–Thiem formula given by:
Q r = 1.366 K M s eff lg R r w = 1.366 K M s eff lg 10 s eff K r w
where the value of seff = s+lw is the effective piezometric head drawdown.
As an example, the results of calculations of the seepage flux for rnor = 1.0 and the boundary head Hnor = 2 for different values of the well penetration are summarized in Table A4 and Figure 17b. Table A4 includes the seepage flux Q nor n from numerical simulations, Q nor s calculated using the Dupuit–Thiem formula, and Q nor r from Equation (5). Table A4 also shows that the relative error of calculations Q nor r using Equation (5) compared to the Q nor n is from −24.88% to 1.1%, while the relative error of using the Dupuit–Thiem formula is from −33.17% to −51.87%.

5.3. Confined-Unconfined Aquifer

For the case of a partially penetrating ARW in a confined-unconfined aquifer, shown in Figure 1c, the results of simulations of the seepage flux can be described using a modified version of the Dupuit–Thiem formula given by:
Q r = 1.366 K ( 2 H M M 2 h w 2 ) lg R r w = 1.366 K ( 2 H M M 2 h w 2 ) lg 10 s eff K r w
As an example, the results of numerical simulations, calculations using the Dupuit–Thiem formula, and formula (7), for rnor = 1.0 and the boundary head Hnor = 7/6, are summarized in Table A5 and shown in Figure 17c. The relative error of calculations of Qr based on Equation (6), compared to numerical simulations, ranges from 2.37% to 8.87%.

5.4. Comparison of Modified Dupuit–Thiem Formulae with Results of Laboratory Sand-Tank Experiments

A comparison of calculations of the seepage flux based on the original and modified Dupuit formulae and sand-tank experiments are summarized in Figure 18 and Table A6. Figure 18 indicates that the seepage flux calculated using the modified Dupuit formula (Equation 6) is closer to the results of laboratory experiments. The seepage flux calculated from Dupuit formula is less than that from a modified Dupuit formula, with larger values being from laboratory experiments. However, the relative errors of Q nor s compared to Q nor r depend on the well length and the water table drawdown. As the well degree of penetration increases and the water table drawdown decreases, the relative error of Q nor s compared to Q nor r is increasing, as shown in Table A6. The data given in Table A4, Table A5 and Table A6 confirm the validity of the application of a modified Dupuit–Thiem formula for calculations of the seepage flux for a single partially penetrating ARW.

6. Concluding Remarks

Laboratory sand-tank experiments and numerical simulations were conducted to assess the seepage characteristics of single partially and fully penetrating ARWs in homogeneous, isotropic unconfined, confined, and unconfined-confined aquifers. The seepage characteristics of the ARWs are affected by the seepage face boundary condition along an outflow boundary of the ARW, because the interior of the well casing remains open to the ambient air and is affected by the atmospheric pressure, so that water is discharging under gravity through the open bottom of the well casing. The experimental and modeling results indicate that (1) the relationships between the seepage flux and the degree of penetration, and seepage flux and the well radius, can be described using the exponential function (Equation 1) and the power function (Equation 2), respectively; (2) for the ARW’s degree of penetration less than the critical value ζ c , the piezometric head above the well is greater than zero with a piezometric head descent zone formed near the well.
The modeling results were used to develop (a) a series of type curves to evaluate the effect of the critical degree of penetration for different well radii and different aquifer water levels; and (b) the modified Dupuit and Dupuit–Thiem formulae by replacing the drawdown s by an effective drawdown seff = s+lw. It is shown that for the ARWs in an unconfined aquifer with the degree of penetration exceeding a critical value, the seepage flux can be calculated based on the original Dupuit formula, and for confined or unconfined-confined aquifers—based on the Dupuit–Thiem formulae for a fully penetrating well.

Author Contributions

Conceptualization, W.W., T.J., J.D., and B.F.; Data curation, W.W.; Formal analysis, W.W., B.F. and J.D.; Funding acquisition, W.W.; Experimental investigations, W.W. and Y.L.; Methodology, W.W., B.F. and T.J.; Project administration, W.W.; Resources, W.W. and Y.L.; Software, W.W. and T.J.; Supervision, B.F.; Validation, W.W., B.F., J.D. and Y.L.; Writing—original draft, W.W. and J.D.; Writing—review and editing, W.W., B.F., T.J. and J.D. All authors have read and agreed to the published version of the manuscript.

Funding

Wenxue Wang: 41602298; Wenxue Wang: 2020HYTP014; Boris Faybishenko: DE-AC02-05CH11231.

Acknowledgments

Wang, W.X. would like to acknowledge the support of the National Natural Science Foundation of China under Grant No. 41.602.298 and the Young Talent Promotion Project of Henan Province (2020HYTP014). B. Faybishenko was supported by the U.S. Department of Energy, Office of Science, Office of Advanced Scientific Computing under Contract No. DE-AC02-05CH11231.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Results of calculations of the errors of the seepage flux Qnor (for lnor = 1 and rnor = 4.5, 9.5, 14.5 and 19.5) from numerical simulations and Dupuit formula in comparison to the laboratory sand-tank experiments.
Table A1. Results of calculations of the errors of the seepage flux Qnor (for lnor = 1 and rnor = 4.5, 9.5, 14.5 and 19.5) from numerical simulations and Dupuit formula in comparison to the laboratory sand-tank experiments.
Well Radius/rnorSeepage Flux QnorErrors (%)
Lab ExperimentsDupuit FormulaNumerical Simulations
Q nor L Q nor s Q nor n Errors   of   Q nor n Errors   of   Q nor s
4.51.2941.161.195−7.653.02
9.51.5171.4251.467−3.302.95
14.51.6361.6361.6963.673.67
19.51.7921.8261.9086.474.49
Table A2. For different well radii and the boundary head.
Table A2. For different well radii and the boundary head.
Hnor11.52.02.53.03.54.0
rnor
0.10.650.96
0.20.630.93
0.50.580.91
10.550.840.970.98
20.490.760.930.970.98
40.420.660.850.940.960.980.99
60.350.570.780.900.930.960.98
100.290.480.670.800.90.930.95
Table A3. Seepage flux Qnor calculations for an unconfined aquifer with rnor = 1.0 and Hnor = 1.
Table A3. Seepage flux Qnor calculations for an unconfined aquifer with rnor = 1.0 and Hnor = 1.
lnorsnorsnor+lnorSeepage Flux QnorErrors (%)
Q nor n Q nor r Q nor s Error   of   Q nor r Error   of   Q nor s
1/20.3620.8620.9370.9360.657−0.10−29.88
1/30.2820.6150.8460.8600.5631.60−33.41
1/60.1860.3530.6350.6470.4271.93−32.73
1/100.1430.2430.4970.5110.3562.82−28.45
1/150.1110.1780.4090.4140.2991.22−27.04
1/300.0800.1130.3020.3030.2370.39−21.37
1/600.0610.0780.2350.2330.197−0.79−16.04
Notes: Q nor s represents the results of calculations using the Dupuit formula, Q nor n represents the results of numerical simulations, and Q nor r represents the results of calculations using the modified Dupuit formula.
Table A4. Parameters used for calculations of the seepage flux for a confined aquifer with rnor = 1.0 and Hnor = 2.
Table A4. Parameters used for calculations of the seepage flux for a confined aquifer with rnor = 1.0 and Hnor = 2.
lnorsnorsnor+lnorQnorErrors (%)
Q nor n Q nor r Q nor s Error   of   Q nor r Error   of   Q nor s
10.6521.6521.0501.1530.5239.83−50.18
5/60.6281.4611.0311.0420.5071.10−50.81
2/30.5751.2420.9720.9130.472−6.10−51.46
1/20.4990.9990.8740.7670.421−12.20−51.87
1/30.4020.7360.7230.6050.354−16.35−51.11
1/60.2690.4360.5040.3900.257−22.72−49.02
1/100.2030.3030.3840.2890.206−24.75−46.30
1/150.1700.2370.3120.2410.181−22.87−42.19
1/300.1200.1530.2270.1720.139−24.47−38.88
1/600.0960.1120.1760.1320.118−24.88−33.17
Table A5. Parameters used for calculations of the seepage flux for a confined-unconfined aquifer with rnor = 1.0 and Hnor = 7/6.
Table A5. Parameters used for calculations of the seepage flux for a confined-unconfined aquifer with rnor = 1.0 and Hnor = 7/6.
lnorsnor h nor w snor+lnorQnorErrors (%)
Q nor n Q nor r Q nor s Error   of   Q nor r Error   of   Q nor s
10.3910.0000.7790.9651.0001.2193.6426.28
5/60.4110.0000.7890.9771.0001.2072.3723.50
2/30.4050.0950.7860.9771.0071.2023.1623.09
1/20.3500.3170.7580.9280.9761.1535.1624.25
1/30.2790.5550.7230.8050.8601.0066.8424.88
1/60.1920.8080.6790.5810.6330.7258.8724.73
1/100.1470.9200.6570.4510.4900.5548.8222.90
Table A6. Comparison of seepage flux calculations using data from the laboratory experiments, original and modified Dupuit formulae. Errors of original and modified Dupuit formulae are calculated in comparison to the results of laboratory sand-tank experiments.
Table A6. Comparison of seepage flux calculations using data from the laboratory experiments, original and modified Dupuit formulae. Errors of original and modified Dupuit formulae are calculated in comparison to the results of laboratory sand-tank experiments.
rnor(lnor).snorQnorErrors (%)
Lab ExperimentsModified Dupuit FormulaDupuit Formula
Q nor L Q nor r Q nor s Errors   of   Q nor r Errors   of   Q nor s
4.5(1/30)0.200.9880.8530.801−13.66−18.93
4.5(1/12)0.231.0611.0110.875−4.71−17.53
4.5(1/6)0.331.2171.2641.0663.86−12.41
9.5(1/30)0.371.3131.2671.220−3.50−7.08
9.5(1/12)0.401.4111.3661.267−3.19−10.21
Figure A1. Photograph of the monitoring piezometer tubes.
Figure A1. Photograph of the monitoring piezometer tubes.
Water 12 00919 g0a1
Figure A2. Photographs of seepage well screens (a) and well covers (b) used in the sand-tank experiments.
Figure A2. Photographs of seepage well screens (a) and well covers (b) used in the sand-tank experiments.
Water 12 00919 g0a2
Figure A3. Particle size distribution curve of a sand used for the sand-tank experiments.
Figure A3. Particle size distribution curve of a sand used for the sand-tank experiments.
Water 12 00919 g0a3
Figure A4. 3-D numerical discretization of the flow domains for simulations of: (a) the sand-tank experiments, and (b) the confined aquifer.
Figure A4. 3-D numerical discretization of the flow domains for simulations of: (a) the sand-tank experiments, and (b) the confined aquifer.
Water 12 00919 g0a4

References

  1. Schubert, J.P. Reducing Water Leakage into Underground Coal Mines by Aquifer Dewatering; No. ANL/LRP-CP-9; CONF-7809134-1; Argonne National Lab.: Argonne, IL, USA, 1978.
  2. Sgambat, J.P.; LaBella, E.A.; Roebuck, S. Effects of Coal Mining on Ground-Water in the Eastern United States; EPA 60017-80-120; U.S. Environmental Protection Agency: Washington, DC, USA, 1980; p. 183.
  3. Niskovskiy, Y.; Vasianovich, A. Investigation of possibility to apply untraditional and ecologically good methods of coal mining under sea bed. In Proceedings of the Sixth International Offshore and Polar Engineering Conference, Los Angeles, CA, USA, 26–31 May 1996. [Google Scholar]
  4. Wu, Q.; Zhou, W.F. Prediction of inflow from overlying aquifers into coalmines: A case study in Jinggezhuang Coalmine, Kailuan, China. Environ. Geol. 2008, 55, 775–780. [Google Scholar] [CrossRef]
  5. Zhou, N.Q.; Vermeer, P.A.; Lou, R.; Tang, Y.; Jiang, S. Numerical simulation of deep foundation pit dewatering and optimization of controlling land subsidence. Eng. Geol. 2010, 114, 251–260. [Google Scholar] [CrossRef]
  6. Struzina, M.; Müller, M.; Drebenstedt, C.; Mansel, H.; Jolas, P. Dewatering of Multi-aquifer Unconsolidated Rock Opencast Mines: Alternative Solutions with Horizontal Wells. Mine Water Environ. 2011, 30, 90–104. [Google Scholar] [CrossRef]
  7. Pujades, E.; López, A.; Carrera, J.; Vázquez-Suñé, E.; Jurado, A. Barrier effect of underground structures on aquifers. Eng. Geol. 2012, 145, 41–49. [Google Scholar] [CrossRef]
  8. Wu, Y.X.; Shen, S.L.; Yuan, D.J. Characteristics of dewatering induced drawdown curve under blocking effect of retaining wall in aquifer. J. Hydrol. 2016, 539, 554–566. [Google Scholar] [CrossRef]
  9. Brownlow, J.W.; James, S.C.; Yelderman, J.C. Influence of Hydraulic Fracturing on Overlying Aquifers in the Presence of Leaky Abandoned Wells. Groundwater 2016, 54, 781–792. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, W.X.; Sui, W.H.; Faybishenko, B.; Stringfellow, W.T. Permeability variations within mining-induced fractured rock mass and its influence on groundwater inrush. Environ. Earth Sci. 2016, 75, 1–15. [Google Scholar] [CrossRef]
  11. Xu, J.P.; Chen, M.Y.; Hong, H.; Xu, X.Q.; Liu, B. Differences of hydraulic conductivities from pumping and draining tests with intensity effects. J. Hydrol. 2019, 579, 124165. [Google Scholar] [CrossRef]
  12. Cartwright, K.; Hunt, C.S. Hydrogeologic Aspects of Coal Mining in Illinois: An Overview. Champaign, III; Illinois State Geological Survey: Champaign, IL, USA, 1981; p. 19.
  13. Hill, J.G.; Price, D.R. The impact of deep mining on an overlying aquifer in western Pennsylvania. Groundw. Monit. Remediat. 1983, 3, 138–143. [Google Scholar] [CrossRef]
  14. Younger, P.L.; Wolkersdorfer, C. Mining impacts on the fresh water environment: Technical and managerial guidelines for catchment scale management. Mine Water Environ. 2004, 23, 2–80. [Google Scholar]
  15. Zhang, J.C. Investigations of water inrushes from aquifers under coal seams. Int. J. Rock Mech. Min. Sci. 2005, 42, 350–360. [Google Scholar] [CrossRef]
  16. Booth, C.J. Groundwater as an environmental constraint of longwall coal mining. Environ. Geol. 2006, 49, 796–803. [Google Scholar] [CrossRef]
  17. Hang, Y.; Zhang, G.L.; Yang, G.Y. Numerical simulation of dewatering thick unconsolidated aquifers for safety of underground coal mining. Min. Sci. Technol. 2009, 19, 312–316. [Google Scholar] [CrossRef]
  18. Barnes, M.R.; Vermeulen, P.D. Guide to groundwater monitoring for the coal industry. Water SA 2012, 38, 831–836. [Google Scholar] [CrossRef] [Green Version]
  19. Newman, C.; Agioutantis, Z.; Leon, G.B.J. Assessment of potential impacts to surface and subsurface water bodies due to longwall mining. Int. J. Min. Sci. Technol. 2017, 27, 57–64. [Google Scholar] [CrossRef]
  20. Meiri, D. Unconfined groundwater flow calculation into a tunnel. J. Hydrol. 1985, 82, 69–75. [Google Scholar] [CrossRef]
  21. Liu, X.X.; Shen, S.; Xu, Y.S. Analytical approach for time-dependent groundwater inflow into shield tunnel face in confined aquifer. Int. J. Numer. Anal. Methods Geomech. 2018, 42, 655–773. [Google Scholar] [CrossRef]
  22. Hwang, J.H.; Lu, C.C. A semi-analytical method for analyzing the tunnel water inflow. Tunn. Undergr. Space Technol. 2007, 22, 39–46. [Google Scholar] [CrossRef]
  23. Pulido-Velazquez, D.; Llopis-Albert, C.; Peña-Haro, S.; Pulido-Velazquez, M. Efficient conceptual model for simulating the effect of aquifer heterogeneity on natural groundwater discharge to rivers. Adv. Water Res. 2011, 34, 1377–1389. [Google Scholar] [CrossRef]
  24. Mahdi Rasouli, M. Groundwater Seepage Rate (GSR); a new method for prediction of groundwater inflow into jointed rock tunnels. Tunn. Undergr. Space Technol. 2018, 71, 505–517. [Google Scholar] [CrossRef]
  25. Goodman, R. Groundwater inflows during tunnel driving. Eng. Geol. 1965, 2, 39–56. [Google Scholar]
  26. Shizhong, L. An Analytical Solution for Steady Flow into a Tunnel. Groundwater 1999, 37, 23–26. [Google Scholar]
  27. El Tani, M. Circular tunnel in a semi-infinite aquifer. Tunn. Undergr. Space Technol. 2003, 18, 49–55. [Google Scholar] [CrossRef]
  28. Kolymbas, D.; Wagner, P. Groundwater ingress to tunnels–The exact analytical solution. Tunn. Undergr. Space Technol. 2007, 22, 23–27. [Google Scholar] [CrossRef]
  29. Corps of Engineers US Army (USACE). Design, Construction, and Maintenance of Relief Wells; Engineer Manual: Washington, DC, USA, 1994.
  30. Tammetta, P. Estimation of the Change in Storage Capacity above Mined Longwall Panels. Groundwater 2016, 54, 646–655. [Google Scholar] [CrossRef] [PubMed]
  31. Shamsai, A.; Narasimhan, T.N. A numerical investigation of free surface-seepage face relationship under steady state flow conditions. Water Resour. Res. 1991, 27, 409–421. [Google Scholar] [CrossRef]
  32. Butler, J.J.; Healey, J.M. Relationship between Pumping-Test and Slug-Test Parameters: Scale Effect or Artifact? Groundwater 1998, 36, 305–312. [Google Scholar] [CrossRef]
  33. Tsoflias, G.P.; Halihan, T.; Sharp, J.M. Monitoring pumping test response in a fractured aquifer using ground-penetrating radar. Water Resour. Res. 2001, 37, 1221–1229. [Google Scholar] [CrossRef] [Green Version]
  34. Simpson, M.J.; Clement, T.P.; Gallop, T.A. Laboratory and Numerical Investigation of Flow and Transport near a Seepage-Face Boundary. Groundwater 2003, 41, 690–700. [Google Scholar] [CrossRef]
  35. Li, W.; Englert, A.; Cirpka, O.A.; Vereecken, H. Three-dimensional geostatistical inversion of flowmeter and pumping test data. Groundwater 2008, 46, 193–201. [Google Scholar] [CrossRef]
  36. Sethi, R. A dual-well step drawdown method for the estimation of linear and nonlinear flow parameters and wellbore skin factor in confined aquifer systems. J. Hydrol. 2011, 400, 187–194. [Google Scholar] [CrossRef]
  37. Ji, S.H.; Koh, Y.K. Nonlinear groundwater flow during a slug test in fractured rock. J. Hydrol. 2015, 520, 30–36. [Google Scholar] [CrossRef]
  38. Sun, H. A semi-analytical solution for slug tests in an unconfined aquifer considering unsaturated flow. J. Hydrol. 2016, 532, 29–36. [Google Scholar] [CrossRef]
  39. Wu, Y.X.; Shen, J.S.; Chen, W.C.; Hino, T. Semi-analytical solution to pumping test data with barrier, wellbore storage, and partial penetration effects. Eng. Geol. 2017, 226, 44–51. [Google Scholar] [CrossRef]
  40. Dupuit, J.É.J. Études Théoriques et Pratiques sur le Mouvement des Eaux Dans les Canaux Découverts et à Travers les Terrains Perméables: Avec des Considérations Relatives au Régime des Grandes Eaux, au Débouché à leur Donner, et à la Marche des Alluvions dans les Rivières à Fond Mobile; Dunod: Paris, France, 1863. [Google Scholar]
  41. Thiem, G. Hydrologische Methoden: Dissertation zur Erlangung der Wurde eines; JM Gebhardt: Leipzig, Germany, 1906. [Google Scholar]
  42. Boulton, N.S. The drawdown of the water-table under non-steady conditions near a pumped well in an unconfined formation. Proc. Inst. Civ. Eng. 1954, 3, 564–579. [Google Scholar] [CrossRef]
  43. Streltsova, T. Unsteady radial flow in an unconfined aquifer. Water Resour. Res. 1972, 8, 1059–1066. [Google Scholar] [CrossRef]
  44. Neuman, S.P. Effect of partial penetration on flow in unconfined aquifers considering delayed gravity response. Water Resour. Res. 1974, 10, 303–312. [Google Scholar] [CrossRef]
  45. Neuman, S.P. Theory of flow in unconfined aquifers considering delayed response of the water table. Water Resour. Res. 1972, 8, 1031–1045. [Google Scholar] [CrossRef]
  46. Ruud, N.C.; Kabala, Z.J. Response of a partially penetrating well in a heterogeneous aquifer: Integrated well-face flux vs. uniform well-face flux boundary conditions. J. Hydrol. 1997, 194, 76–94. [Google Scholar] [CrossRef]
  47. Cassiani, G.; Kabala, Z.J. Hydraulics of a partially penetrating well: Solution to a mixed-type boundary value problem via dual integral equations. J. Hydrol. 1998, 211, 100–111. [Google Scholar] [CrossRef]
  48. Hemker, C.J. Transient well flow in vertically heterogeneous aquifers. J. Hydrol. 1999, 225, 1–18. [Google Scholar] [CrossRef]
  49. Chang, C.C.; Chen, C.S. A flowing partially penetrating well in a finite-thickness aquifer: A mixed-type initial boundary value problem. J. Hydrol. 2003, 271, 101–118. [Google Scholar] [CrossRef]
  50. Perina, T.; Lee, T.C. General well function for pumping from a confined, leaky, or unconfined aquifer. J. Hydrol. 2006, 317, 239–260. [Google Scholar] [CrossRef]
  51. Feng, Q.G.; Zhan, H.B. On the aquitard–aquifer interface flow and the drawdown sensitivity with a partially penetrating pumping well in an anisotropic leaky confined aquifer. J. Hydrol. 2015, 521, 74–83. [Google Scholar] [CrossRef]
  52. Mishra, P.K.; Kuhlman, K.L. Unconfined aquifer flow theory: From Dupuit to present. In Advances in Hydrogeology; Springer: New York, NY, USA, 2013; pp. 185–202. [Google Scholar]
  53. Charnyi, I.A. A rigorous derivation of Dupuit’s formula for unconfined seepage with seepage surface. Dokl. Akad. Nauk Mosc. 1951, 79, 937–940. (In Russian) [Google Scholar]
  54. Shercliff, J.A. Seepage flow in unconfined aquifers. J. Fluid Mech. 1975, 71, 181–192. [Google Scholar] [CrossRef]
  55. Bhunya, P.K.; Anjaneyulu, B. Study of flow analysis in semiconfined aquifer using a sand tank model. ISH J. Hydraul. Eng. 2001, 7, 13–21. [Google Scholar] [CrossRef]
  56. Faybishenko, B.A. Hydraulic behavior of quasi-saturated soils in the presence of entrapped air: Laboratory experiments. Water Resour. Res. 1995, 31, 2421–2435. [Google Scholar] [CrossRef]
  57. Faybishenko, B.A.; Javandel, I.; Witherspoon, P.A. Hydrodynamics of the Capture Zone for Contaminant Transport with Partially Penetrating Well in a Confined Aquifer. Water Resour. Res. 1995, 31, 859–866. [Google Scholar] [CrossRef]
Figure 1. Schematic diagrams of a partially penetrating ascending relief well Ascending Relief Wells (ARW) in (a) an unconfined aquifer, (b) a confined aquifer, and (c) a confined-unconfined aquifer. On the figures, s is the drawdown of water table, hw is the height of the water table above the well, lw is well screen length, rw is well radius, H is the initial piezometric head, M is the aquifer thickness, and K is hydraulic conductivity.
Figure 1. Schematic diagrams of a partially penetrating ascending relief well Ascending Relief Wells (ARW) in (a) an unconfined aquifer, (b) a confined aquifer, and (c) a confined-unconfined aquifer. On the figures, s is the drawdown of water table, hw is the height of the water table above the well, lw is well screen length, rw is well radius, H is the initial piezometric head, M is the aquifer thickness, and K is hydraulic conductivity.
Water 12 00919 g001
Figure 2. Schematic depicting initial and boundary conditions of the radial flow to a partially penetrating well in the sand-tank tests.
Figure 2. Schematic depicting initial and boundary conditions of the radial flow to a partially penetrating well in the sand-tank tests.
Water 12 00919 g002
Figure 3. Schematic design diagram (a) and a photograph (b) of the seepage sand-tank used for laboratory experiments.
Figure 3. Schematic design diagram (a) and a photograph (b) of the seepage sand-tank used for laboratory experiments.
Water 12 00919 g003
Figure 4. Piezometric head pnor distribution for the case of rnor = 19.5 and lnor = 1, i.e., a fully penetrating well in the laboratory sand-tank experiments.
Figure 4. Piezometric head pnor distribution for the case of rnor = 19.5 and lnor = 1, i.e., a fully penetrating well in the laboratory sand-tank experiments.
Water 12 00919 g004
Figure 5. Piezometric head pnor distributions for the case of rnor = 4.5 and lnor = 1/30 in the laboratory sand-tank experiments.
Figure 5. Piezometric head pnor distributions for the case of rnor = 4.5 and lnor = 1/30 in the laboratory sand-tank experiments.
Water 12 00919 g005
Figure 6. Piezometric head pnor distribution at radial distances Rnor of 0 and 1/6 for the case of rnor = 4.5 and lnor = 1/30 in the laboratory experiment.
Figure 6. Piezometric head pnor distribution at radial distances Rnor of 0 and 1/6 for the case of rnor = 4.5 and lnor = 1/30 in the laboratory experiment.
Water 12 00919 g006
Figure 7. Piezometric head pnor distribution for the case of rnor = 4.5 for different penetrating lengths lnor. For partially penetrating wells: (a) lnor = 1/30, (b) lnor = 1/12, (c) lnor = 1/6, (d) lnor = 1/3, (e) lnor = 1/2, and for a fully penetrating well: (f) lnor = 1. The dashed red lines indicate the length of the well screened interval.
Figure 7. Piezometric head pnor distribution for the case of rnor = 4.5 for different penetrating lengths lnor. For partially penetrating wells: (a) lnor = 1/30, (b) lnor = 1/12, (c) lnor = 1/6, (d) lnor = 1/3, (e) lnor = 1/2, and for a fully penetrating well: (f) lnor = 1. The dashed red lines indicate the length of the well screened interval.
Water 12 00919 g007
Figure 8. Relationship between the seepage flux Qnor and the degree of penetration lnor for different well radii.
Figure 8. Relationship between the seepage flux Qnor and the degree of penetration lnor for different well radii.
Water 12 00919 g008
Figure 9. Relationship between seepage flux Qnor and well radius rnor for different degree of penetration lnor, based on the results of the sand-tank experiments.
Figure 9. Relationship between seepage flux Qnor and well radius rnor for different degree of penetration lnor, based on the results of the sand-tank experiments.
Water 12 00919 g009
Figure 10. Relationships between the seepage flux Qnor and the degree of well penetration lnor for different boundary conditions: (a) Hnor = 1, (b) Hnor = 2, (c) Hnor = 3, and (d) Hnor = 4.
Figure 10. Relationships between the seepage flux Qnor and the degree of well penetration lnor for different boundary conditions: (a) Hnor = 1, (b) Hnor = 2, (c) Hnor = 3, and (d) Hnor = 4.
Water 12 00919 g010aWater 12 00919 g010b
Figure 11. Relationships between the seepage flux Qnor and the well radius rnor for different boundary conditions: (a) Hnor = 1, (b) Hnor = 2, (c) Hnor = 3, (d) Hnor = 4.
Figure 11. Relationships between the seepage flux Qnor and the well radius rnor for different boundary conditions: (a) Hnor = 1, (b) Hnor = 2, (c) Hnor = 3, (d) Hnor = 4.
Water 12 00919 g011
Figure 12. Results of simulations of the water piezometric head pnor distribution around the ARW in the unconfined aquifer for different well length lnor: (1) lnor = 1/60, (2) lnor = 1/30, (3) lnor = 1/15, (4) lnor = 1/10, (5) lnor = 1/6, (6) lnor = 1/3, (7) lnor = 1/2, (8) lnor = 2/3, (9) lnor = 5/6, (10) lnor = 1 (for the case of Hnor = 1 and rnor = 1.0). The x-axis is the normalized radial distance from the well.
Figure 12. Results of simulations of the water piezometric head pnor distribution around the ARW in the unconfined aquifer for different well length lnor: (1) lnor = 1/60, (2) lnor = 1/30, (3) lnor = 1/15, (4) lnor = 1/10, (5) lnor = 1/6, (6) lnor = 1/3, (7) lnor = 1/2, (8) lnor = 2/3, (9) lnor = 5/6, (10) lnor = 1 (for the case of Hnor = 1 and rnor = 1.0). The x-axis is the normalized radial distance from the well.
Water 12 00919 g012
Figure 13. The flow net around the ARW in the unconfined aquifer for different well length lnor: (1) lnor = 1/60, (2) lnor = 1/30, (3) lnor = 1/15, (4) lnor = 1/10, (5) lnor = 1/6, (6) lnor = 1/3, (7) lnor = 1/2, (8) lnor = 2/3, (9) lnor = 5/6, (10) lnor = 1 (for the case Hnor = 1 and r = 1.0). The x-axis is the radial distance from the well.
Figure 13. The flow net around the ARW in the unconfined aquifer for different well length lnor: (1) lnor = 1/60, (2) lnor = 1/30, (3) lnor = 1/15, (4) lnor = 1/10, (5) lnor = 1/6, (6) lnor = 1/3, (7) lnor = 1/2, (8) lnor = 2/3, (9) lnor = 5/6, (10) lnor = 1 (for the case Hnor = 1 and r = 1.0). The x-axis is the radial distance from the well.
Water 12 00919 g013
Figure 14. Relationship between ζ c and Hnor for different well radii.
Figure 14. Relationship between ζ c and Hnor for different well radii.
Water 12 00919 g014
Figure 15. Type curves of ζ c for as a function of Hnor and rnor.
Figure 15. Type curves of ζ c for as a function of Hnor and rnor.
Water 12 00919 g015
Figure 16. Comparison of the relationships between the seepage flux Qnor with lnor = 1 and rnor = 4.5, 9.5, 14.5 and 19.5 from numerical simulations, Dupuit formula and laboratory sand-tank experiments.
Figure 16. Comparison of the relationships between the seepage flux Qnor with lnor = 1 and rnor = 4.5, 9.5, 14.5 and 19.5 from numerical simulations, Dupuit formula and laboratory sand-tank experiments.
Water 12 00919 g016
Figure 17. Seepage flux calculations using the results of numerical modeling (Qn), original Dupuit (Qs) and modified Dupuit formula (Qr) for the cases of: (a) an unconfined aquifer, (b) confined aquifer, and (c) confined-unconfined aquifer.
Figure 17. Seepage flux calculations using the results of numerical modeling (Qn), original Dupuit (Qs) and modified Dupuit formula (Qr) for the cases of: (a) an unconfined aquifer, (b) confined aquifer, and (c) confined-unconfined aquifer.
Water 12 00919 g017
Figure 18. Comparison of the results of calculations of the seepage flux from the laboratory sand-tank experiments, original and modified Dupuit formulae.
Figure 18. Comparison of the results of calculations of the seepage flux from the laboratory sand-tank experiments, original and modified Dupuit formulae.
Water 12 00919 g018
Table 1. Formulae for steady-state seepage flux in unconfined, confined, and confined-unconfined aquifers.
Table 1. Formulae for steady-state seepage flux in unconfined, confined, and confined-unconfined aquifers.
Types of AquifersSeepage Flux (Q)
Unconfined
(Dupuit)
Q = 1.366 K ( 2 H s ) s lg R r w
Confined
(Dupuit–Thiem)
Q = 2.73 K M s lg R r w
Confined-Unconfined
(Dupuit–Thiem)
Q = 1.366 K ( 2 H M M 2 h w 2 ) lg R r w
Notes: Q is the well seepage flux, K is hydraulic conductivity, H is the initial piezometric head, s is the drawdown (i.e., dewatering depth), rw is the well radius, M is the thickness of the aquifer, R is the radius of influence of the pumping well, hw is the water level in the pumping well.

Share and Cite

MDPI and ACS Style

Wang, W.; Faybishenko, B.; Jiang, T.; Dong, J.; Li, Y. Seepage Characteristics of a Single Ascending Relief Well Dewatering an Overlying Aquifer. Water 2020, 12, 919. https://doi.org/10.3390/w12030919

AMA Style

Wang W, Faybishenko B, Jiang T, Dong J, Li Y. Seepage Characteristics of a Single Ascending Relief Well Dewatering an Overlying Aquifer. Water. 2020; 12(3):919. https://doi.org/10.3390/w12030919

Chicago/Turabian Style

Wang, Wenxue, Boris Faybishenko, Tong Jiang, Jinyu Dong, and Yang Li. 2020. "Seepage Characteristics of a Single Ascending Relief Well Dewatering an Overlying Aquifer" Water 12, no. 3: 919. https://doi.org/10.3390/w12030919

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop