Research Article  Open Access
Jianyou Yu, Shensheng Ge, Lingchao Hou, Yiteng Zhang, Yuzhu Liu, "Limit Drainage of Separated Tunnel Based on Ecological Groundwater Table: A Case Study", Advances in Civil Engineering, vol. 2021, Article ID 2669924, 16 pages, 2021. https://doi.org/10.1155/2021/2669924
Limit Drainage of Separated Tunnel Based on Ecological Groundwater Table: A Case Study
Abstract
After the excavation of the tunnel in waterrich regions, the groundwater is continuously discharged into the tunnel. Excessive discharge will cause the groundwater table to drop, which will destroy the normal growth of vegetation and ecological balance. In order to protect the ecological environment, the key is to develop an effective method to determine the tunnel drainage. In this paper, a method for calculating the limit drainage of the tunnel is proposed based on areawell method and the concept of ecological groundwater table in ecology and agronomy. Some conclusions are drawn as follows. (1) When the ecological groundwater table is constant, with the increase of groundwater discharge, the drainage influence range of the tunnel decreases, the rainfall supplement quantity within the drainage influence range decreases, and the total groundwater discharge quantity Q_{t} also decreases gradually. Since the decreasing rate of is greater than that of Q_{t}, there is a unit drainage that makes the total groundwater discharge equal to rainfall supplement. (2) A separated tunnel under construction in Fuzhou City is taken as an example; the limit drainage q_{0} of the separated double tunnel is 0.48 m^{3}/(m·d) considering ecological groundwater table.
1. Introduction
With the development of economy and the improvement of transportation facilities in various regions, the construction of largescale tunnels is increasing year by year, which not only brings new opportunities to the development of transportation and economic growth, but also brings challenges to the local ecological, hydrological, and geological environment [1–9]. During tunnel excavation and operation in waterrich regions, groundwater leaks through the tunnel, so the groundwater table drops [8, 10, 11] as shown in Figure 1, leading to the reduction of surface water storage and even springs drying up. After the tunnel construction in Seoul, the annual discharge of groundwater was 63 million m^{3}, and the groundwater table decreased from 16.85 m in 1996 to 20.40 m in 2003 [12]. The construction of a railway tunnel through the Hallandsås ridge in the southwest of Sweden resulted in drying up private wells and springs close to the tunnel axis due to groundwater drawdown [6]. The construction of Zhongliangshan tunnel in Chongqing caused 29 ground surface collapses and dried up 48 springs [13]. The influence area of surface water was about 7.2 km^{2}, which seriously endangered the living and production water of local residents [14]. At the same time, the decrease of surface water and groundwater caused the shortage of water supply in the local ecosystem, which lead to the deterioration of plant growth. The width of tree rings in the influence area of Zhongliangshan tunnel in Chongqing continued to decrease [15], and the growth rate of dominant plants (Pinus massoniana, Cunninghamia lanceolata, etc.) declined in varying degrees due to water deficiency [16]. Tunnels excavation also decreased the soil water contents and then changed plants water uptake patterns in karst areas [17].
Such adverse effects of tunnel drainage occurred, because the waterproofingdrainage measures in the past for tunnels passing through the waterrich regions were mainly free drainage to reduce the influence of water pressure on the stability of lining structure [18–20]. With the improvement of environmental protection awareness, the guiding ideology of pressureregulated drainage has gradually been applied in engineering practice [21–23]. The principle of pressureregulated drainage is that a portion of groundwater is allowed to discharge out of the tunnel and the level of the groundwater is regulated either by a regulating value or by a raised outlet opening.
For tunnels, there have been some research on the calculation methods and standards of the groundwater drainage. For instance, the allowable drainage of Oslofjord crosssea tunnel in Norway was determined according to the drainage capacity and economic benefits, and the maximum allowable drainage of the land section under residential and recreational areas was 0.288 m^{3}/(m·d) [24]. There existed ancient buildings in the upper part of Storhaug Tunnel in Norway; the water drainage range was set to 0.0432–0.144 m^{3}/(m·d) to protect them [25]. Due to the sensitivity of urban recreational area to groundwater drainage, the maximum allowable drainage of Baneheia Tunnel in Norway was determined to be 0.03024 m^{3}/(m·d) [25]. To protect the water environment and vegetation, the maximum allowable drainage of Geleshan Tunnel in China was 1.0 m^{3}/(m·d), and the actual drainage in operation was 0.95 m^{3}/(m·d), which met the design allowable drainage [26]. In order to reduce the loss of surface water around Qiyueshan Tunnel, the allowable maximum drainage was 3.0 m^{3}/(m·d), which was determined based on the engineering geological conditions, hydrological conditions, and the drainage capacity of equipment [22]. Based on the ecological water demand of vegetation, [27] studied the drainage funnel after tunnel drainage and obtained the calculation method for tunnel drainage to maintain the ecological balance of vegetation and groundwater.
At present, the study of groundwater in coal mining is more mature, proposing water inrush solution measures [28] and giving the protection of aquifers [29]. However, there is no relatively unified calculation method of tunnel drainage and the limit standard with wider applicability; the previous provisions of tunnel drainage have their own considerations and specific provisions in terms of the use function and geological environment. In this research, from the perspective of vegetation demand and drawing on the areawell method in groundwater dynamics, the relationship between the groundwater drawdown and the time, and drainage quantity of the tunnel are obtained. Then through the balance analysis of the total groundwater discharge Q_{t} and the total groundwater supply mainly from rainfall infiltration supply, the method for calculating the drainage to maintain the balance of groundwater is proposed. Finally, taking a separated doublehole tunnel under construction in Fuzhou as an example, its limit drainage is determined.
2. Relationship between Groundwater Table and Vegetation
Groundwater provides plants with water, salt, organic nutrients, and heat. When the groundwater table drops enough, the root system cannot absorb enough water from the capillary zone, leading to vegetation degradation. If water is in short supply for a long time, vegetation will disappear. When the groundwater table is enough shallow, it may cause soil salinization and threaten the growth of vegetation. According to this phenomenon, the concept of ecological groundwater table is put forward, which means that the groundwater table meets the normal growth of plants and the problems of vegetation degradation, land desertification, and soil salinization will not occur under no irrigation [30].
Yang et al. [31] used the Yushen coal mine district in Mu Us Desert to study quantitative analysis of the relationship between vegetation and groundwater table in different ecological and geological (ecogeological) environments. Liu et al. [32] used the Yushenfu mining area as a case study; the ecosafety groundwater table was determined to be 1.5∼5.0 m based on the curve describing the relationship between groundwater table and NDVI (the Normalized Difference Vegetation Index, as obtained by optical satellites, is one of the most popular indicators for the status of vegetation). Horton et al. [33] analyzed the physiological response of three different plants to the change of groundwater table and proposed the range of groundwater table to maintain the normal growth of two plants. Newman et al. [34] summarized the research results of many scholars and proposed the relationship between groundwater and vegetation growth from the perspective of ecology hydrology. Shi et al. [35] studied the relationship between groundwater and vegetation based on vegetation coverage and proposed the groundwater table suitable for different vegetation growth in different regions. Cheng et al. [36] analyzed the relationship between different plant populations, plant coverage, and groundwater table, pointed out that the groundwater table and soil moisture content in aeration zone were important factors affecting vegetation growth, and proposed the groundwater table suitable for different vegetation.
The above analysis shows that vegetation is closely related to groundwater table. Therefore, when determining the tunnel drainage, it is necessary to introduce the ecological groundwater table. Groundwater drawdown caused by tunnel drainage must be ensured not exceeding the ecological groundwater table to protect the vegetation of the tunnel site.
3. Theoretic Methods for Calculating Tunnel Drainage
Based on the ecological groundwater table, the specific progress of the theoretical method for calculating tunnel drainage is shown in Figure 2. First, based on areawell method in groundwater dynamics, the groundwater seepage model of a singlehole tunnel and separated tunnel is established, which reveals this relationship between the water inflow quantity and groundwater table drawdown. Second, according to the proposed groundwater seepage model, the influence range of groundwater drainage funnel is derived under the premise of controlling the maximum groundwater drawdown within the range required by normal growth of vegetation. Then the rainfall supplement quantity within the influence range is obtained by using the method of rainfall infiltration coefficient, and then is compared to the total groundwater discharge quantity Q_{t}. Through adjusting groundwater discharge q_{0} to make equal to Q_{t}, q_{0} is the maximum allowable groundwater discharge in the tunnel site.
3.1. Relationship between Tunnel Drainage and Groundwater Drawdown
3.1.1. The Groundwater Drawdown under the Drainage of the SingleHole Tunnel
Due to the similarity between tunnel seepage and well pumping, the current calculation method for groundwater drawdown after tunnel drainage is generally based on well flow theory in groundwater dynamics. In the well flow theory, the assumption is that (1) the well diameter is infinitely small; (2) the well is vertical; and (3) the drawdown occurs around the wellhead. While the tunnel is equivalent to the horizontal well, the drawdown funnel will appear at the top of the tunnel. Therefore, it can be assumed that the tunnel wall is composed of infinite pumping wells, and the superposition of infinite pumping wells is considered as shown in Figure 3 [14].
For the superposition of infinite pumping wells, the areawell method is adopted to obtain the final groundwater drawdown. The areawell method is that when the flow distribution of a well hole is relatively uniform on a certain geometric shape (such as rectangle and circle), the well group can be regarded as a whole, and the pumping from each well is regarded as the uniform “evaporation” on the distribution area of the well group [37].
If the areawell method is adopted for determining tunnel drainage, the following assumptions should be made: (1) the tunnel cross section is assumed to be standard circular, and the tunnel as a whole is cylindrical. However, the shape of the tunnel cross section is usually horseshoeshaped and needs to be converted into a circle by area equivalence principle. (2) The water inflow on the tunnel wall is uniform, stable, and equal everywhere.
The length of the tunnel is L and the radius is R_{0}, as shown in Figure 4(a). When the tunnel wall surface is expanded into a plane, the plane size is L = 2l_{x} and the width is 2πR_{0} = 2l_{y}, as shown in Figure 4(b) [27]. Assuming that the plane is composed of numerous tiny seepage surfaces, the water inflow on the plane is equal everywhere. The action of microelement on any point M can be regarded as the action of a single well, so the action of the whole face to the point M can be obtained by integrating the action of microelement.
(a)
(b)
Assuming that the daily water inflow per linear meter of a tunnel is q, the total water inflow per unit time on the seepage surface can be expressed bywhere Q is the daily tunnel water inflow (m^{3}/d); q is the daily water inflow per linear meter of a tunnel (m^{3}·d^{−1}m^{−1}); L is the length of the tunnel (m).
Therefore, water inflow per unit area can be expressed by where is the water inflow per unit area (m/d).
The coordinate system is taken as shown in Figure 4(b). The coordinate system represents the position of the microelement, and the coordinate system represents the position of calculation point for drawdown.
Because it takes a long time for the groundwater around the tunnel to form a stable flow state, the calculation method of unsteady groundwater flow fits the actual situation better. Theis well function [38] is the most widely used formula involving the unsteady flow of groundwater, as expressed by
And the dimensionless parameter is defined by
The well function is defined by where s is the groundwater drawdown (m); Q is the daily tunnel water inflow (m^{3}/d); T is the transmissibility coefficient (m^{2}/d); a is the coefficient of pressure conductivity, where a = T/μ, and μ is the specific yield; t is time from the start of leakage (d); r is the distance between the groundwater drawdown point and the well (m).
The mathematical substitution of Theis well function is carried out. Setting
Then we get the differential of x as expressed by
Thus, equation (5) can be expressed by
Then, equation (8) (the well function after substitution) is brought into equation (3), and the drawdown at any point can be expressed by
Assuming that there is any point M (x, y) in the tunnel area, the drawdown under the action of seepage surface with discharge Q is solved by the following method.
Taking any point as the center and getting microelement on the seepage surface, the water flow of the microelement is . Point M (x, y) under the seepage action of differential surface , its differential drawdown ds can be expressed by the following equation based on equation (9):
The drawdown of point M (x, y) under the action of seepage face can be obtained by integrating in the whole face, as expressed by
The mathematical substitution of equation (11) is carried out. Setting
Then we get the differential of as expressed by
Therefore, the following equation can be obtained:where is called Gaussian error function.
In the same way, equation (15) can be obtained as follows:
Substituting equations (14) and (15) into equation (11) and noticing that , equation (16) can be obtained as follows:
The mathematical substitution of equation (16) is carried out. Setting , then we get the differential , and equation (16) is converted into
Define the following function:
Then, the drawdown at any point M (x, y) under the action of rectangular area well can be expressed by
Here, A_{r} is expressed by where A_{r} can be called the well function of rectangular area well. In this paper, a program for calculating value in any of the scopes of and is developed based on MATLAB. For the center point of a singlehole tunnel, setting x = 0 and y = 0 in equation (19), the groundwater drawdown of the tunnel center point can be expressed by where s_{c} is the groundwater drawdown of the tunnel center point.
Obviously, the maximum groundwater drawdown is at the center of the tunnel, so equation (21) is the calculation formula of the maximum groundwater table drawdown caused by tunnel drainage. s_{c} is a function of t, which can be calculated when s_{c} and tunnel water inflow q are given.
3.1.2. The Groundwater Drawdown under the Drainage of the Separated DoubleHole Tunnel
For the separated doublehole tunnel (Figure 5(a)), the areawell method is also adopted, and the two tunnels are equivalent to two area wells, as shown in Figure 5(b).
(a)
(b)
For a singlehole tunnel, the drawdown equation of any point M (x, y) under the action of rectangular area wells is expressed by equations (1), (2), (19), and (20). For separated tunnels, the drawdown of point M is the sum of the drawdown caused by the drainage of left tunnel and right tunnel, as shown in Figure 5(b). Setting the drainage quantity of the left tunnel as q_{1}, the drainage quantity of the right tunnel is q_{2}. The diameter, length, and width of the left tunnel and right tunnel are the same, and they are r, L, and B, respectively. The distance between the left and right tunnel is D, and the coordinate origin is set in the center of the left tunnel, as shown in Figure 5(b). According to equations (1), (2), (19), and (20), the drawdown s_{1} due to drainage in left tunnel to point M is expressed by
For the left tunnel, the right tunnel is equivalent to a right translation (D + B) along the ycoordinate direction. Therefore, the drawdown s_{2} due to seepage in right tunnel to point M is expressed by
For the drawdown s of M point, s = s_{1} + s_{2} is formed by superposition of the drawdown of the left tunnel and right tunnel. Except that x and y are independent variables and other parameters in s function are constants, the maximum value of s appears in the center of the left tunnel and the right tunnel through mathematical means. Therefore, the center drawdown of any tunnel is selected as the groundwater control drawdown. Setting that the drawdown at the center of the left tunnel is the control drawdown, when x = 0 and y = 0, the drawdown at the center of the left tunnel can be expressed by where s_{l} is the groundwater drawdown at the center of the left tunnel due to left tunnel seepage.
Similarly, the drawdown at the center of the left tunnel due to right tunnel seepage can be expressed by where s_{r} is the groundwater drawdown at the center of the left tunnel due to right tunnel seepage.
If the maximum allowable drawdown s is given, the time t to reach the maximum drawdown can be obtained by equations (28) and (29).
3.2. Drainage Influence Range around Tunnel
3.2.1. Horizontal Influence Range
For the tunnel sites without obvious recharge boundary and watershed, a place with groundwater drawdown of approximately 0 m is the boundary of the influence range. Pumping over time creates a pumping funnel around the well. Depending on the flow characteristics of the well, a descending funnel will form above the tunnel when the discharge is too high. However, unlike a well, the length of the tunnel is much greater than the width; therefore, the shape of the pumping funnel is an inverted elliptical cone. The elliptical scope of the drawdown funnel is the influence scope of the tunnel and the rainfall infiltration recharge area [27], as shown in Figure 6.
The drawdown of each point is the superposition of double tunnels; then, the influence range in the middle of the tunnel is the largest in transverse direction. The coordinate of the horizontal groundwater drawdown of approximately 0 m is taken as the ellipse short radius b in the influence range.
Taking the center point of the left tunnel as the origin, the coordinate system is shown in Figure 5(b). Then, the expression of drawdown at any point is the sum of the left tunnel’s drawdown s_{h1} and the right tunnel’s drawdown s_{h2}. Setting x = 0 and y = −b in equation (24), groundwater drawdown s_{h1} at transverse −b in the middle of the tunnel caused by the left tunnel’s seepage can be expressed by
Setting x = 0 and y = −b in equation (27), groundwater drawdown s_{h2} at transverse −b in the middle of the tunnel caused by the right tunnel’s seepage can be expressed by
After considering the superposition of two tunnels, the groundwater drawdown at y = −b in the middle of the tunnel can be expressed by
When s_{h} = 0, the b is obtained, which is the ellipse short radius in the influence range of the tunnel.
3.2.2. Longitudinal Influence Range
In the longitudinal direction, considering the symmetry of the left tunnel and the right tunnel, the groundwater drawdown on the center line of the left tunnel and the right tunnel is the same. Therefore, a place with drawdown of approximately 0 m in the direction of x axis is taken as the ellipse long radius l in the influence range.
Setting x = l and y = 0 in equation (24), groundwater drawdown s_{z1} at longitudinal l caused by the left tunnel’s seepage can be expressed by
Setting x = l and y = 0 in equation (27), groundwater drawdown s_{z1} at longitudinal l caused by the right tunnel’s seepage can be expressed by
After considering the superposition of double tunnels, the groundwater drawdown at x = l can be expressed by
When s_{z} = 0, the l is obtained, which is the ellipse long radius in the influence range.
3.2.3. Influence Range of Tunnel Seepage Funnel
After the short radius and long radius of the ellipse are respectively determined by equations (32) and (35), the influence range of the tunnel drainage funnel can be determined, as shown in Figure 6. Accordingly, the ellipse area F can be obtained, as expressed by equation (36). F is also the rainfall infiltration recharge area:
3.3. Determination of Rainfall Infiltration Recharge
For tunnel engineering, rainfall infiltration is usually the main recharge source of groundwater. The general process of rainfall infiltration is that rainwater first replenishes water in the vadose zone above the groundwater surface, and when the water content in the vadose zone reaches the field water capacity, rainwater will be replenishment to the groundwater. In order to accurately understand the recharge of rainfall to groundwater, rainfall infiltration coefficient is an important evaluation parameter. Rainfall infiltration coefficient is the ratio of rainfall infiltration supplement to the corresponding precipitation ; that is, [39].
In plain, hill, and mountain areas, when there is no test data, the precipitation infiltration coefficient of aquifer can be approximately determined according to Table 1 [40].

Rainfall can be obtained according to local meteorological statistics. According to the area, annual rainfall and infiltration coefficient are obtained in the recharge regions, and the amount of daily recharge within the influence scope can be obtained, as expressed by [39]where is the daily average rainfall infiltration supply (m^{3}/d); p is rainfall infiltration coefficient; is annual rainfall; F is the area replenished by rainfall infiltration (m^{3}).
3.4. Groundwater Balance Analysis
In the process of groundwater drawdown calculation, the water inflow per unit time is q; then, the total water inflow in time t can be expressed by
The total supply of precipitation in time t can be expressed by
After reaching allowable drawdown, the total amount of drainage Q_{t} and the amount of recharge are obtained. (1) If Q_{t} = , the groundwater is just in balance, and the groundwater table will no longer continue to decline. (2) If Q_{t} > , the drainage is too large. In this case, the groundwater will be drained and the ecological environment will be affected. Therefore, water inflow should be reduced to balance the discharge and supply. (3) If Q_{t} < , the discharge is appropriate, and no action is required.
4. Case Study
4.1. Tunnel Description
The proposed tunnel is a separated doublehole tunnel with an average distance of about 20 m. The total length of the left and right lines is 1176 and 1140 m, respectively, and about 1000 m of them pass through the aquifer. The tunnel has a width of 20.02 m, a height of 13.53 m, and an equivalent circle radius of 10.8 m. The main surrounding rock of the tunnel site is unstable completely weathered granite and moderately weathered granite. The average thickness of the underground aquifer is H = 30 m. The permeability coefficient of the surrounding rock is taken from the geotechnical engineering investigation report of the detailed investigation stage of the extension line project of North Industrial Road (section I) in Fuzhou City [41]. The comprehensive value of the permeability coefficient in this paper is 0.3 m/d, the specific yield is μ = 0.02 according to the empirical value of common rock specific yield, and the water conductivity coefficient of aquifer is T = kH = 9 m^{2}/d. According to Table 2, the annual average precipitation of Fuzhou is 1359 mm. The lithology of this area is relatively broken, so the rainfall infiltration coefficient is taken as p = 0.2. The allowable drainage calculation parameters are summarized in Table 3.


4.2. Ecological Groundwater Table
According to the field survey, the main plants above the tunnel are camphor tree, mango tree, longan tree, and banyan tree, among which the diameter of camphor tree is 0.3 m–0.9 m, mango tree is 0.3–0.8 m, longan tree is 0.2–0.4 m, and banyan tree is 0.7–1.5 m.
At present, there is no research on ecological groundwater table for these trees. Therefore, when comprehensively considering ecological protection and construction cost, the ecological groundwater table is 5 m referring to Table 4 [42].

5. Results and Discussion
5.1. Time t for the Groundwater Table to Reach the Ecological Balance of Groundwater
According to Section 4.2, the maximum allowable groundwater drawdown in this area is s = 5 m.
According to the geotechnical engineering investigation report of the detailed investigation stage of the extension line project of North Industrial Road (section I) in Fuzhou City [41], the normal water inflow of the left tunnel is 543.77 m^{3}/d and the maximum water inflow is 1882.91 m^{3}/d. The normal water inflow of the right tunnel is 506.20 m^{3}/d and the maximum water inflow is 2011.31 m^{3}/d. Take the average water inflow of a singlehole tunnel q = 1.24 m^{3}/(m·d).
Given that q = 1.24 m^{3}/(m·d) and allowable drawdown s = s_{l} + s_{r} = 5 m, the time t of drawdown at the center of the left tunnel can be calculated by equations (28) and (29). The relationship between s and time t at the center of the tunnel is obtained through MATLAB programming, as shown in Figure 7. It is obtained that the allowable drawdown reaches 5 m when t = 9 d (integer). In other words, when there is no supply, the groundwater table in the center of the left tunnel drops by 5 m when t = 9 d.
5.2. Configuration and Expansion Form of Groundwater Falling Funnel
For the drawdown s of any point, s = s_{1} + s_{2} is formed by superposition of the drawdown of the left tunnel and right tunnel. Setting x = 0, the drawdown of groundwater table in the transverse section of the tunnel center can be obtained, as shown in Figure 8.
Setting y = 0, the drawdown of groundwater table in the longitudinal section of the tunnel center can be obtained, as shown in Figure 9.
In order to more directly observe the trend of groundwater table changing with time, a threedimensional diagram of drawdown was obtained by MATLAB, as shown in Figure 10.
(a)
(b)
(c)
It can be seen from Figures 8–10 that with the increase of time, the drawdown of groundwater table increases continuously, and the drawdown at the center of the two tunnels is the largest.
5.3. Influence Range Caused by Tunnel Drainage
According to equation (32), when t = 9 d, set . When s_{h} < 0.01 m, b = 260 m is obtained, that is, the ellipse short radius of tunnel influence area.
According to equation (35), when t = 9 d, set . When s_{z} < 0.01 m, l = 720 m is obtained, that is, the ellipse long radius of tunnel influence area.
Here the precision of l and b is 1e − 2, which meet the engineering calculation accuracy.
5.4. Analysis of Groundwater Balance
5.4.1. Determination of Groundwater Total Discharge
The total discharge of groundwater is determined by equation (38), so we can get Q_{t} = (Q_{1} + Q_{2}) = (q_{1}Lt + q_{2}Lt) = 22320 m^{3}.
5.4.2. Determination of Rainfall Infiltration Recharge
First, the area of precipitation supply area is obtained by equation (36), so we can get m^{3}.
Second, the daily precipitation supply in the influence area can be obtained by equation (37), so we can get m^{3}.
Third, the total precipitation supply in 9 d can be obtained by equation (39), so we can get m^{3}.
Fourth, it can be found that Q_{t} is greater than , which shows that longterm discharge will cause groundwater drying up, so appropriate measures should be taken to reduce groundwater discharge, so as to achieve groundwater balance.
5.5. Determination of Balanced Drainage for Maintaining Groundwater
According to Section 5.4, when the discharge of a singlehole tunnel is q = 1.24 m^{3}/(m·d), Q_{t} is greater than . Therefore, the value of q needs to be reduced to make Q_{t} = . Because the range of q cannot be obtained before calculation, and the calculation efficiency will be serious if the range is too large, so the Bisection method is adopted.
Specific steps: (1) first determine the value interval of q (0, 1.24) and take the median value of 0.62 in the interval for trial calculation. In this case, − Q_{t} (m^{3}) < 0. (2) Determine the value interval of q (0, 0.62) again and take the median value within the interval 0.31 for trial calculation. In this case, − Q_{t} (m^{3}) > 0. (3) Determine the value interval of q (0.31, 0.62) again and take the median value of 0.465 in the interval for trial calculation. In this case, − Q_{t} (m^{3}) > 0. (4) Determine the value range of q (0.465, 0.62) again. Since q = 0.465 m^{3}/(m·d), is close to Q_{t}, so q = 0.48, q = 0.49, and q = 0.50 are selected for trial calculation. When q = 0.48, − Q_{t} = 1283 m^{3} > 0. When q = 0.49, − Q_{t} = −237 m^{3} < 0. Table 5 and Figures 11–13 can be obtained by selecting different q values for trial calculation.

From Figures 11–13, it can be seen that with the increase of unit water inflow q_{0}, the time t required to reach a certain drawdown gradually decreases, indicating that the larger the groundwater discharge is, the faster the groundwater table drops. The influence range gradually decreases because the influence range is a function of time t, and the decrease of t leads to the decrease of the influence range of tunnel. The total recharge and discharge decrease with the increase of unit water inflow when reaching the established drawdown. It can be seen from Figure 13 that when q_{0} is about 0.48 m^{3}/(m·d), the total discharge is equal to the total recharge, and when q is less than 0.48 m^{3}/(m·d), the recharge is greater than the discharge, which is beneficial to the groundwater balance, and when q is greater than 0.48 m^{3}/(m·d), the total discharge is greater than the recharge, which means that the recharge cannot reach the balance with the discharge, so the unit water inflow q_{0} = 0.48 m^{3}/(m·d) is more suitable for discharge.
After one year’s operation of the tunnel, the vegetation above the tunnel is growing well (Figure 14), which proves that the calculation method for this limit drainage is reliable.
6. Conclusion
By introducing the concept of ecological groundwater table in ecology, based on the method of groundwater dynamics and considering rainfall recharge, this paper analyzes the calculating method of tunnel drainage to maintain groundwater balance and draws the following conclusions:(1)From the perspective of protecting ecological vegetation from destruction, the relationship between groundwater drawdown and drainage in a singlehole tunnel is deduced based on the areawell method, as shown in equation (19). On that basis, the relationship between groundwater drawdown and drainage of the separated double tunnel is deduced by considering the superposition effect, as shown in equations (24) and (27). The purpose is to control the groundwater drawdown within the range required for the normal growth of vegetation and to maintain the ecological balance of vegetation in the tunnel site.(2)Based on areawell method, the relationship between the water inflow quantity and groundwater drawdown is obtained. It is also used to calculate the influence scope of the tunnel with a drawdown of approximately 0 m; thus, the area replenished by rainfall infiltration is determined. Based on the rainfall infiltration coefficient, the rainfall infiltration supplement is obtained, which is compared with the groundwater discharge Q_{t}. By continuously adjusting the unit water inflow q_{0}, the total groundwater discharge Q_{t} and the rainfall infiltration supplement are equal. When they are equal, the unit water inflow q_{0} can maintain the groundwater balance.(3)Based on the actual example of a separated double tunnel under construction, the paper puts forward the unit water inflow of the tunnel to keep the groundwater balance, which is 0.48 m^{3}/(m·d). In the process of analysis, the twodimensional and threedimensional diagram of groundwater drawdown with time change are described by MATLAB programming, which gives an intuitive understanding of the shape of groundwater drawdown funnel.(4)When the groundwater drawdown is constant, the time t and influence range gradually decrease with the increase of unit water inflow q_{0}. The total rainfall recharge and tunnel discharge Q_{t} also decrease with the increase of unit water inflow q_{0}.
Data Availability
Data are available upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
This work was financially supported by a grant from the Hebei Provincial Expressway YanqingChongli Construction Office (grant no. YC201903). In addition, the authors would like to thank the Fuzhou Municipal Construction Development Co., Ltd. (grant no. 20191061), which provided the construction site and geological exploration data.
References
 C. Butscher, P. Huggenberger, and E. Zechner, “Impact of tunneling on regional groundwater flow and implications for swelling of clay–sulfate rocks,” Engineering Geology, vol. 117, no. 34, pp. 198–206, 2011. View at: Publisher Site  Google Scholar
 D. Ma, H. Duan, X. Li, Z. Li, Z. Zhou, and T. Li, “Effects of seepageinduced erosion on nonlinear hydraulic properties of broken red sandstones,” Tunnelling and Underground Space Technology, vol. 91, 2019a. View at: Publisher Site  Google Scholar
 F. Mossmark, L. O. Ericsson, M. Norin, and L.O. Dahlström, “Hydrochemical changes caused by underground constructions  a case study of the Kattleberg rail tunnel,” Engineering Geology, vol. 191, pp. 86–98, 2015. View at: Publisher Site  Google Scholar
 E. Pujades, E. VázquezSuñé, L. Culí, J. Carrera, A. Ledesma, and A. Jurado, “Hydrogeological impact assessment by tunnelling at sites of high sensitivity,” Engineering Geology, vol. 193, pp. 421–434, 2015. View at: Publisher Site  Google Scholar
 J. R. Raposo, J. Molinero, and J. Dafonte, “Quantitative evaluation of hydrogeological impact produced by tunnel construction using water balance models,” Engineering Geology, vol. 116, no. 34, pp. 323–332, 2010. View at: Publisher Site  Google Scholar
 A. SjölanderLindqvist, “Conflicting perspectives on water in a Swedish railway tunnel project,” Environmental Values, vol. 14, no. 2, pp. 221–239, 2005. View at: Publisher Site  Google Scholar
 T. Strozzi, R. Caduff, U. Wegmüller, H. Raetzo, and M. Hauser, “Widespread surface subsidence induced in Alpine hard rocks by the construction of the 57 kilometreslong Gotthard Base Tunnel (Switzerland) observed with satellite SAR interferometry,” Procedia Technology, vol. 16, pp. 69–73, 2014. View at: Publisher Site  Google Scholar
 V. Vincenzi, A. Gargini, and N. Goldscheider, “Using tracer tests and hydrological observations to evaluate effects of tunnel drainage on groundwater and surface waters in the Northern Apennines (Italy),” Hydrogeology Journal, vol. 17, no. 1, pp. 135–150, 2009. View at: Publisher Site  Google Scholar
 Q. Xie, S. Dan, and W. Yang, “Evaluation of groundwater environmental disturbance by mountain highway tunnel construction,” Geological Journal of China Universities, vol. 25, no. 3, pp. 437–443, 2019. View at: Google Scholar
 K. Su, Y. Zhou, H. Wu, C. Shi, and L. Zhou, “An analytical method for groundwater inflow into a drained circular tunnel,” Groundwater, vol. 55, no. 5, pp. 712–721, 2017. View at: Publisher Site  Google Scholar
 C. Yoo, “Interaction between tunneling and groundwaternumerical investigation using three dimensional stresspore pressure coupled analysis,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, no. 2, pp. 240–250, 2005. View at: Publisher Site  Google Scholar
 G. T. Chae, S. T. Yun, B. Y. Choi et al., “Hydrochemistry of urban groundwater, Seoul, Korea: the impact of subway tunnels on groundwater quality,” Journal of Contaminant Hydrology, vol. 101, no. 14, pp. 42–52, 2008. View at: Publisher Site  Google Scholar
 J. Luo and W. Fu, “Destruction caused by tunneling works to the circulatory system of underground water in karst areas  a case study of the Zhongliang mountains in Chongqing,” Journal of Southwest Agricultural University, vol. 27, no. 4, pp. 432–435, 2005. View at: Google Scholar
 P. Cheng, L. Zhao, Z. Luo et al., “Analytical solution for the limiting drainage of a mountain tunnel based on areawell theory,” Tunnelling and Underground Space Technology, vol. 84, pp. 22–30, 2019. View at: Publisher Site  Google Scholar
 W. Zheng, X. Wang, Y. Tang, H. Liu, M. Wang, and L. Zhang, “Use of tree rings as indicator for groundwater level drawdown caused by tunnel excavation in Zhongliang Mountains, Chongqing, Southwest China,” Environmental Earth Sciences, vol. 76, no. 15, p. 522, 2017. View at: Publisher Site  Google Scholar
 F. Wang, W. Zheng, H. Xu, and J. Chen, “Influence of groundwater loss induced by karst tunnel construction on growth of surrounding vegetation and its countermeasures,” Tunnel Construction, vol. 38, no. 6, pp. 915–923, 2018. View at: Google Scholar
 J. Liu, L. Shen, Z. Wang et al., “Response of plants water uptake patterns to tunnels excavation based on stable isotopes in a karst trough valley,” Journal Of Hydrology, vol. 571, pp. 485–493, 2019. View at: Publisher Site  Google Scholar
 P. Arjnoi, J.H. Jeong, C.Y. Kim, and K.H. Park, “Effect of drainage conditions on porewater pressure distributions and lining stresses in drained tunnels,” Tunnelling and Underground Space Technology, vol. 24, no. 4, pp. 376–389, 2009. View at: Publisher Site  Google Scholar
 X. Wang, z. Tan, M. Wang, M. Zhang, and H. Ming, “Theoretical and experimental study of external water pressure on tunnel lining in controlled drainage under high water level,” Tunnelling and Underground Space Technology, vol. 23, no. 5, pp. 552–560, 2008. View at: Publisher Site  Google Scholar
 G. Xu, C. He, C. Qi, Z. Zhang, and C. Dai, “Study on Water Pressure Distribution and Inner Force of Drainaged Segment Lining,” in Proceedings of the ICPTT 2012: Better Pipeline Infrastructure for a Better Life, pp. 1502–1511, Wuhan, China, October 2013. View at: Publisher Site  Google Scholar
 C. Zhang, D. Zhang, M. Wang, and Y. Xiang, “Study on appropriate parameters of grouting circle for tunnels with limiting discharge lining in high water pressure and waterenriched region Chinese,” Journal of Rock Mechanics and Engineering, vol. 26, no. 11, pp. 2270–2276, 2007. View at: Google Scholar
 M.Q. Zhang and H.J. Huang, “Technology for blocking water by grouting for highpressure cranny water of qiyueshan tunnel,” Journal of Railway Engineering Society, vol. 136, no. 1, pp. 68–73, 2010. View at: Google Scholar
 Z. Zou, “Application of waterproof technology by curtain grouting to Geleshan Tunnel construction,” Modern Tunnelling Technology, vol. 40, no. 1, pp. 44–49, 2003. View at: Google Scholar
 O. Blindheim and E. Ovstedal, “Design Principles and Construction Methods for Water Control in Subsea Road Tunnels in rock,” Water Control in Norwegian Tunnelling, vol. 12, 2002. View at: Google Scholar
 E. Grov, “Water Control in Norwegian Tunnelling,” in Proceedings of the AITESITA Downunder 2002: 28th ITA General Assembly and World Tunnel Congress, p. 656, Institution of Engineers, Sydney, Australia, March 2002, Modern TunnelsChallenges and Solutions. View at: Google Scholar
 N. Li, Study of the Impact on the Water Environment during the Construction of the Geleshan Railway Tunnel, Southwest Jiaotong University, Chengdu, China, 2004.
 P. Cheng, L.h. Zhao, L. Li, J.f. Zou, and W. Luo, “Limiting drainage criterion for groundwater of mountain tunnel,” Journal of Central South University, vol. 21, no. 12, pp. 4660–4668, 2014. View at: Publisher Site  Google Scholar
 D. Ma, H. Duan, J. Liu, X. Li, and Z. Zhou, “The role of gangue on the mitigation of mininginduced hazards and environmental pollution: an experimental investigation,” The Science of the Total Environment, vol. 664, pp. 436–448, 2019. View at: Publisher Site  Google Scholar
 D. Ma, J. Zhang, H. Duan et al., “Reutilization of gangue wastes in underground backfilling mining: overburden aquifer protection,” Chemosphere, vol. 264, no. 1, Article ID 128400, 2021. View at: Publisher Site  Google Scholar
 H. Zhang, “The ecological equilibrium embedment depth of groundwater in arid region,” Site Investigation Science and Technology, vol. 6, no. 20, pp. 9–13, 1992. View at: Google Scholar
 Z. Yang, W. Li, X. Li, and J. He, “Quantitative analysis of the relationship between vegetation and groundwater buried depth: a case study of a coal mine district in Western China,” Ecological Indicators, vol. 102, pp. 770–782, 2019. View at: Publisher Site  Google Scholar
 S. Liu, W. Li, W. Qiao, Q. Wang, Y. Hu, and Z. Wang, “Effect of natural conditions and mining activities on vegetation variations in arid and semiarid mining regions,” Ecological Indicators, vol. 103, pp. 331–345, 2019. View at: Publisher Site  Google Scholar
 J. L. Horton, T. E. Kolb, and S. C. Hart, “Physiological response to groundwater depth varies among species and with river flow regulation,” Ecological Applications, vol. 11, no. 4, pp. 1046–1059, 2001. View at: Publisher Site  Google Scholar
 B. D. Newman, B. P. Wilcox, S. R. Archer et al., “Ecohydrology of water‐limited environments: a scientific vision,” Water Resources Research, vol. 42, no. 6, 2006. View at: Publisher Site  Google Scholar
 R. Shi, X. Li, X. Dong, P.n. Yang, and L. Liu, “Research on the relationship between natural vegetation growth and groundwater in Yanqi basin,” Journal of Natural Resources, vol. 24, no. 12, pp. 2096–2103, 2009. View at: Google Scholar
 D.H. Cheng, W.K. Wang, G.C. Hou, H.B. Yang, Y. Li, and E.Y. Zhang, “Relationship between vegetation and groundwater in Mu Us desert,” Journal of Jilin University(Earth Science Edition), vol. 42, no. 1, pp. 184–189, 2012. View at: Google Scholar
 C. Chen and M. Lin, Groundwater Dynamics, China University of Geosciences Press, Wuhan, China, 1998.
 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
 R. Zhang, X. Liang, M. Jin, L. Wan, and Q. Yu, Fundamentals of Hydrogeology, Geological Publishing House, Beijing, China, 2011.
 D. Zhu and Q. Li, “Predication of outflow in tunnel based on experimental permeating coefficient of precipitation,” Journal of Railway Engineering Society, vol. 45, no. 1, pp. 100–102, 1995. View at: Google Scholar
 Z. Xiao and J. Zhang, The Geotechnical Engineering Investigation Report of the Detailed Investigation Stage of the Extension Line Project of North Industrial Road (Section I) in Fuzhou City, Fuzhou Investigation and Surveying Institute, Fuzhou, China, 2018.
 S. Wang, Q. Huang, L. Fan, and W. Wang, Coal Mining and Ecological Water Level protection in Ecologically Fragile Areas, Science China Press, Beijing, China, 2010.
Copyright
Copyright © 2021 Jianyou Yu et al. 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.