Limit Drainage of Separated Tunnel Based on Ecological Groundwater Table: A Case Study

After the excavation of the tunnel in water-rich 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 area-well 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 Wt within the drainage influence range decreases, and the total groundwater discharge quantity Qt also decreases gradually. Since the decreasing rate of Wt is greater than that of Qt, 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 q0 of the separated double tunnel is 0.48m/(m·d) considering ecological groundwater table.


Introduction
With the development of economy and the improvement of transportation facilities in various regions, the construction of large-scale 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][2][3][4][5][6][7][8][9]. During tunnel excavation and operation in water-rich 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]. e 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]. e construction of Zhongliangshan tunnel in Chongqing caused 29 ground surface collapses and dried up 48 springs [13]. e 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. e 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 waterproofing-drainage measures in the past for tunnels passing through the water-rich regions were mainly free drainage to reduce the influence of water pressure on the stability of lining structure [18][19][20]. With the improvement of environmental protection awareness, the guiding ideology of pressure-regulated drainage has gradually been applied in engineering practice [21][22][23].
e 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 cross-sea 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]. ere 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 area-well method in groundwater dynamics, the relationship between the groundwater drawdown and the time, and drainage quantity of the tunnel are obtained. en through the balance analysis of the total groundwater discharge Q t and the total groundwater supply W t mainly from rainfall infiltration supply, the method for calculating the drainage to maintain the balance of groundwater is proposed. Finally, taking a separated double-hole tunnel under construction in Fuzhou as an example, its limit drainage is determined.  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

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 area-well method in groundwater dynamics, the groundwater seepage model of a single-hole 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. en the rainfall supplement quantity W t within the influence range is obtained by using the method of rainfall infiltration coefficient, and then W t is compared to the total groundwater discharge quantity Q t . rough adjusting groundwater discharge q 0 to make W t equal to Q t , q 0 is the maximum allowable groundwater discharge in the tunnel site. 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. erefore, 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.

Relationship between Tunnel Drainage and Groundwater
e area-well 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 area-well 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 horseshoe-shaped and needs to be converted into a circle by area equivalence principle. (2) e water inflow on the tunnel wall is uniform, stable, and equal everywhere. e 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. e 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.
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 by where 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). erefore, water inflow per unit area can be expressed by where ε is the water inflow per unit area (m/d). e coordinate system is taken as shown in Figure 4(b). e x ′ oy ′ coordinate system represents the position of the microelement, and the xoy 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. eis well function [38] is the most widely used formula involving the unsteady flow of groundwater, as expressed by

Advances in Civil Engineering
Infinite pumping wells

Advances in Civil Engineering
And the dimensionless parameter u is defined by e 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). e mathematical substitution of eis well function is carried out. Setting en we get the differential of x as expressed by us, equation (5) can be expressed by en, 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 (x ′ , y ′ ) as the center and getting microelement dx ′ dy ′ on the seepage surface, the water flow of the microelement is εdx ′ dy ′ . Point M (x, y) under the seepage action of differential surface dx ′ dy ′ , its differential drawdown ds can be expressed by the following equation based on equation (9): e drawdown of point M (x, y) under the action of seepage face can be obtained by integrating in the whole face, as expressed by e mathematical substitution of equation (11) is carried out. Setting en we get the differential of η as expressed by erefore, the following equation can be obtained: where In the same way, equation (15) can be obtained as follows: Substituting equations (14) and (15) into equation (11) and noticing that T � aμ, equation (16) can be obtained as follows: e mathematical substitution of equation (16) is carried out. Setting τ � ξ/t, then we get the differential dτ � dξ/t, and equation (16) is converted into Advances in Civil Engineering Define the following function: en, 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 S * (α, β) value in any of the scopes of α and β is developed based on MATLAB. For the center point of a single-hole 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.

e Groundwater Drawdown under the Drainage of the Separated Double-Hole Tunnel.
For the separated doublehole tunnel (Figure 5(a)), the area-well method is also adopted, and the two tunnels are equivalent to two area wells, as shown in Figure 5(b).
For a single-hole 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 . e diameter, length, and width of the left tunnel and right tunnel are the same, and they are r, L, and B, respectively. e 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 y-coordinate direction. erefore, 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. erefore, 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).

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. e 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. e 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. e 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). en, 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

Advances in Civil Engineering
When s h � 0, the b is obtained, which is the ellipse short radius in the influence range of the tunnel.

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. erefore, 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.

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:

Determination of Rainfall Infiltration Recharge.
For tunnel engineering, rainfall infiltration is usually the main recharge source of groundwater. e 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 W r to the corresponding precipitation W; that is, p � W r /W [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 W r is the daily average rainfall infiltration supply (m 3 / d); p is rainfall infiltration coefficient; W is annual rainfall; F is the area replenished by rainfall infiltration (m 3 ).

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 e 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 W t are obtained. (1) If Q t � W t , the groundwater is just in balance, and the groundwater table will no longer continue to decline. (2) If Q t > W t , the drainage is too large. In this case, the groundwater will be drained and the ecological environment will be affected. erefore, water inflow should be reduced to balance the discharge and supply. (3) If Q t < W t , the discharge is appropriate, and no action is required. e 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]. e 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. e lithology of this area is relatively broken, so the rainfall infiltration coefficient is taken as p � 0.2. e allowable drainage calculation parameters are summarized in Table 3. 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.

Ecological Groundwater
At present, there is no research on ecological groundwater table for these trees. erefore, when comprehensively considering ecological protection and construction cost, the ecological groundwater table is 5 m referring to  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.

Time t for the Groundwater
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. e 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 single-hole 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). e 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.

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 three-dimensional diagram of drawdown was obtained by MATLAB, as shown in Figure 10.
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.

Influence Range Caused by Tunnel Drainage.
According to equation (32), when t � 9 d, set s h ≈ 0. 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 s z ≈ 0. 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.

Determination of Groundwater Total Discharge.
e 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 .

Determination of Rainfall Infiltration Recharge.
First, the area of precipitation supply area is obtained by equation (36), so we can get F � πlb � 588106 m 3 .
Second, the daily precipitation supply in the influence area can be obtained by equation (37), so we can get W r � ρWF/365 � 437.71 m 3 .
ird, the total precipitation supply in 9 d can be obtained by equation (39), so we can get W t � W r t � 3939 m 3 .
Fourth, it can be found that Q t is greater than W t , which shows that long-term discharge will cause groundwater drying up, so appropriate measures should be taken to reduce groundwater discharge, so as to achieve groundwater balance.

Determination of Balanced Drainage for Maintaining
Groundwater. According to Section 5.4, when the discharge of a single-hole tunnel is q � 1.24 m 3 /(m·d), Q t is greater than W t . erefore, the value of q needs to be reduced to make Q t � W 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, W t − Q t (m 3          When q � 0.49, W t − 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. e 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. e 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.

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 area-well 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). e 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 area-well 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 W t 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 W t 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 two-dimensional 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 . e total rainfall recharge W t 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
e authors declare that there are no conflicts of interest.