A Semi-Analytical Model and Parameter Analysis of a Collaborative Drainage Scheme for a Deeply Buried Tunnel and Parallel Adit in Water-Rich Ground

: For a railway or highway tunnel under high water pressure during operation, various factors such as the design of the drainage system, material aging, and pipeline blockage must be considered for the tunnels to work with the parallel adit to drain and control the external water pressure on the tunnel lining. A simpliﬁed steady-state seepage model in a semi-inﬁnite multi-connected domain for the tunnel and parallel adit was established and was solved iteratively using the complex variable method and the Schwartz alternating method. After verifying the numerical simulation, parametric analysis, orthogonal tests, and multivariate nonlinear regression were also carried out. Results show that the simpliﬁed theoretical model and its semi-analytical algorithm have a fast convergence speed, and the obtained regression formula is simple, which is suitable for calculation and parameter analysis. A scheme that primarily relies on the parallel adit for drainage would make the external water pressure of the lining facing the parallel adit side less than that of the opposite side. Therefore, to reduce pressure uniformly and meet the requirements of surrounding rock stability, the horizontal net distance between the parallel adit and the tunnel should be no less than the tunnel diameter. Drainage volume of the parallel adit is linearly negatively correlated with tunnel water pressure on the lining and has the most signiﬁcant effect on pressure reduction. The inﬂuence of the vertical distance between the parallel adit and the tunnel on water pressure is small. ends. Theoretical calculations show that the convergence rate is related to the relative positions of the parallel adit and tunnel, burial depth, and drainage volume. ﬁtting ﬁtting The ﬁtting results are compared with the semi-analytical The ﬁtting results are in good agreement with the theoretical calculation results, and their deviations are within 0.8%, showing that the ﬁtting Formula (25) is suitable for the calculation and parameter analysis of the drainage scheme for tunnel and parallel adit.


Introduction
With the ongoing development of China's economy and society, railway and highway network systems have been continuously improved and extended into the western region. However, the mountainous terrain in China's western regions is undulating and has complex geologic conditions. Railroad and highway lines through these areas commonly require the construction of deeply buried tunnels. If groundwater is abundant in these areas, the construction and operation of deep tunnels also face the threat of high external water pressure [1]. High water pressure not only causes water inflow in the tunnel but also leads to significant damage to the lining during tunnel operation.
For example, after the completion of the Meihuashan Tunnel (Liupanshui, China) along the Guikun Railway, the high water pressure caused cracks in the lining, and a large amount of water passed through the cracks [2]. The drainage system of the Xiushan Tunnel (Chongqing, China) gradually aged during operation, resulting in a reduced drainage capacity, causing water pressure surges that led to tunnel lining cracking [3] The Sanxia Jijiapo Tunnel (Yichang, China), Yuanliangshan Tunnel (Chongqing, China), and Qiyueshan Tunnel (Enshi, China) have all undergone damage from water pressure [3]. Therefore, reducing damage caused by high water pressure during tunnel operation has become a necessity.
Water pressure on a lining will increase during tunnel operation due to an inadequate drainage system design or aging and clogging of the drainage system, resulting in lining damage and cracking. To address such problems, designing tunnels and parallel adits using a collaborative drainage scheme is effective [4].
A parallel adit is generally parallel to the main tunnel and can not only assist in tunnel construction but also affect the seepage field around the tunnel through drainage during operation. Water pressure on the tunnel lining is also released. In recent years, scholars have analyzed drainage schemes for tunnels and parallel adits. In the construction of the No. 2 Branch Tunnel of the Northwest Water Supply Project, Sun Shuang et al. [5] added a parallel adit to the side of the main tunnel for drainage, solving the problem of water inflow and collapse. During the construction of the Zhengpantai Tunnel, Wang Lei et al. [6] designed auxiliary pits, such as parallel adits and inclined shafts, in order to reduce water pressure in response to a water surge problem. All of the above studies have empirically investigated the drainage scheme for tunnels and parallel adits but lack quantitative analysis. Sara Zingg. et al. [7,8] studied the effectiveness of various advanced drainage schemes (including drainage through an external parallel adit) on the stabilization of tunnel faces. Several other scholars have also focused on advanced drainage's effect on tunnel face stability [9,10]. These studies are mainly concerned with the influence of parallel adit drainage during the construction phase. Zhao Jinpeng. et al. [11] took the Gongbei Tunnel (Zhuhai, China) as the research object and discussed the waterproofing and drainage system designs. In addition, they analyzed the relationship between the lining external water pressure and drainage volume and the distribution of water pressure around the tunnel with an indoor simulation test and field test. Wang Chunmei et al. [12] used the numerical software FLAC-3D to quantify the effects of water head, envelope rock level, and relative position of the parallel adit and tunnel cross-section on tunnel water pressure considering parallel adit drainage. However, they did not analyze key factors such as the drainage volume of the parallel adit and tunnel.
An analytical model for a tunnel and parallel adit drainage scheme in tunnel operation has yet to be reported. Analytical and semi-analytical methods not only verify the results of numerical simulations but also avoid the complex process of numerical modeling. The water pressure at the tunnel boundary can be easily obtained by entering the relevant parameters. The analytical method can be used to analyze parameters as well as optimize design schemes. Therefore, an analytical model of the drainage scheme for a tunnel and parallel adit has a certain application value.
Tunnel and parallel adit drainage schemes involve the analytical study of seepage fields for parallel-twin tunnels. At present, the research on tunnel seepage fields primarily focuses on single tunnel problems [13][14][15][16][17][18][19]. For seepage in twin tunnels, Zhu Chengwei et al. [20,21] solved the seepage field of a twin tunnel using conformal mapping and the superposition method considering the tunnel lining as isotropic porous material. In addition, Zhang Bingqiang [22] used the mirror image method to obtain the seepage field for twin tunnels. Guo Yufeng et al. [23] used conformal mapping and the Schwartz iteration method to obtain the seepage field of an underwater shallowburied double parallel tunnel under the condition of equal water pressure at the tunnel boundary. The above studies are all based on a given fixed water pressure at the tunnel boundary. In the collaborative drainage scheme for a tunnel and parallel adit, the parallel adit actively controls the drainage volume through the drain borehole, reducing the external water pressure on the tunnel lining. Therefore, drainage volume is the appropriate boundary condition for a parallel adit. In addition, the influence of parallel adit drainage on different parts of the lining varies significantly, resulting in uneven water pressure at the tunnel boundary. Using equal water pressure or an equal total head boundary around the parallel adit and tunnel is not consistent with field conditions. Through the above introduction, previous studies on parallel adit drainage mainly focused on the influence of drainage on tunnel inflow and the stabilization of tunnel faces during construction. The analytical model of the drainage scheme for tunnel and parallel adit in the operation period has not been reported. In this paper, an analytical model of a collaborative drainage scheme is established for the actual drainage conditions of parallel adits and tunnels during operation. A semi-analytical algorithm is obtained using the complex variable method and the Schwartz method. The semi-analytical algorithm has a fast convergence speed. The boundary conditions can be accurately satisfied after a few iterations. In addition, the semi-analytical results are compared with numerical simulation results based on a case study, demonstrating the convergence and accuracy of the algorithm. This study also analyzes the relevant parameters of the drainage scheme. The tunnel crown pressure head fitting formula is also obtained using the orthogonal test and multivariate nonlinear regression method. This study provides a reference for the quantitative design of the operation drainage scheme of deeply buried tunnel in the presence of groundwater.

Simplified Theoretical Model
The collaborative drainage system consists of a tunnel blind pipe network and parallel adit drain borehole. The tunnel is drained through the blind pipe network, and the parallel adit is drained using the drain boreholes ( Figure 1). boundary. Using equal water pressure or an equal total head boundary aroun parallel adit and tunnel is not consistent with field conditions. Through the above introduction, previous studies on parallel adit drainage m focused on the influence of drainage on tunnel inflow and the stabilization of tunne during construction. The analytical model of the drainage scheme for tunnel and p adit in the operation period has not been reported. In this paper, an analytical mod collaborative drainage scheme is established for the actual drainage conditions of p adits and tunnels during operation. A semi-analytical algorithm is obtained usi complex variable method and the Schwartz method. The semi-analytical algorithm fast convergence speed. The boundary conditions can be accurately satisfied after iterations. In addition, the semi-analytical results are compared with num simulation results based on a case study, demonstrating the convergence and accur the algorithm. This study also analyzes the relevant parameters of the drainage sc The tunnel crown pressure head fitting formula is also obtained using the orthogon and multivariate nonlinear regression method. This study provides a reference f quantitative design of the operation drainage scheme of deeply buried tunnel presence of groundwater.

Simplified Theoretical Model
The collaborative drainage system consists of a tunnel blind pipe networ parallel adit drain borehole. The tunnel is drained through the blind pipe networ the parallel adit is drained using the drain boreholes ( Figure 1). According to the equal permeable area principle, the sections of tunnel and p adit are equivalent to circular sections, and the crown position of the tunnel re unchanged ( Figure 2). In addition, the following assumptions are proposed in this (1) Only the external water pressure of the lining (without considering the thickness lining) is considered, and the permeability coefficient of surrounding rock is isotrop constant; (2) The fluid is incompressible and has a state of steady flow; (3) The wate is horizontal and remains constant; (4) The radial seepage velocities are equal at bo tunnel and parallel adit boundaries.
A global Cartesian coordinate system is established with the water table line x-axis and the vertical symmetry axis of the parallel adit as the y-axis. The half-p below the water table contains the half-plane Zp and half-plane Zt, which corresp the coordinate system xpopyp above the parallel adit and xtotyt above the t respectively; Qp and Qt are the drainage volume of the parallel adit and t respectively; rp and rt are the radius of the equivalent circle of the parallel adit and t respectively; dp, dt, and dw are the burial depths of the parallel adit, tunnel, and water respectively; Ht is the pressure head at the tunnel boundary when the parallel ad tunnel are drained; H 0 t is the initial hydrostatic pressure head at the tunnel bound  According to the equal permeable area principle, the sections of tunnel and parallel adit are equivalent to circular sections, and the crown position of the tunnel remains unchanged ( Figure 2). In addition, the following assumptions are proposed in this paper.
(1) Only the external water pressure of the lining (without considering the thickness of the lining) is considered, and the permeability coefficient of surrounding rock is isotropic and constant; (2) The fluid is incompressible and has a state of steady flow; (3) The water table is horizontal and remains constant; (4) The radial seepage velocities are equal at both the tunnel and parallel adit boundaries.
A global Cartesian coordinate system is established with the water table line as the x-axis and the vertical symmetry axis of the parallel adit as the y-axis. The half-plane Z below the water table contains the half-plane Z p and half-plane Z t , which correspond to the coordinate system x p o p y p above the parallel adit and x t o t y t above the tunnel, respectively; Q p and Q t are the drainage volume of the parallel adit and tunnel, respectively; r p and r t are the radius of the equivalent circle of the parallel adit and tunnel, respectively; d p , d t , and d w are the burial depths of the parallel adit, tunnel, and water table, respectively; H t is the pressure head at the tunnel boundary when the parallel adit and tunnel are drained; H 0 t is the initial hydrostatic pressure head at the tunnel boundary; S h is the horizontal net distance between the tunnel and the parallel adit; and S v is the height difference between the parallel adit and the tunnel. The points A, B, C, and D represent the position of tunnel crown, left wall, inverted arch, and right wall respectively. is the horizontal net distance between the tunnel and the parallel adit; and Sv is th difference between the parallel adit and the tunnel. The points A, B, C, and D re the position of tunnel crown, left wall, inverted arch, and right wall respectively. According to Darcy's law, mass conservation law, and the above assumpti total head satisfies the Laplace equation: where h is the total head, given by the sum of the elevation head and pressu (kinetic head is neglected because groundwater flow is assumed to be slow): h = u u is the pore water pressure, γw is the unit weight of water, u/γw is the pressure he y is the elevation head. Using the basic assumptions, the boundary conditions were determined as fo At the water table: At the parallel adit circumference:  where cp = dp + rp − dw, ct = dt + rt − dw, l = rt + rp + Sh; np and nt are the radial directio circumference of the parallel adit and the tunnel, respectively (the outer no positive); K is rock permeability coefficient; and ip and it are the radial hydraulic g of the parallel adit and tunnel circumference, respectively.

Semi-Analytical Solution
The Z plane includes both a tunnel and a parallel adit, which is a steady-state problem in a multi-connected domain and difficult to solve directly. To solv According to Darcy's law, mass conservation law, and the above assumptions, the total head satisfies the Laplace equation: where h is the total head, given by the sum of the elevation head and pressure head (kinetic head is neglected because groundwater flow is assumed to be slow): h = u/γ w + y, u is the pore water pressure, γ w is the unit weight of water, u/γ w is the pressure head, and y is the elevation head.
Using the basic assumptions, the boundary conditions were determined as follows: At the water table: At the parallel adit circumference: At the tunnel circumference: where c p = d p + r p − d w , c t = d t + r t − d w , l = r t + r p + S h ; n p and n t are the radial direction of the circumference of the parallel adit and the tunnel, respectively (the outer normal is positive); K is rock permeability coefficient; and i p and i t are the radial hydraulic gradient of the parallel adit and tunnel circumference, respectively.

Semi-Analytical Solution
The Z plane includes both a tunnel and a parallel adit, which is a steady-state seepage problem in a multi-connected domain and difficult to solve directly. To solve such problems, the Schwartz alternative method can be used to analyze the interaction between the tunnel and parallel adit drainage [23]. The Schwartz alternating method transforms the multi-connected domain problem into multiple single-connected domain problems [24]. The interaction between these single connected domain problems is then iteratively considered. Using this method, the first step involves solving the seepage field for single tunnel drainage in a semi-infinite plane.

General Solution of Single Tunnel Drainage
The water table and tunnel (or parallel adit) circumferences in the Z j plane are conformally mapped to two circles with radius 1 and α j in the ζ j plane by mapping functions (where j = p, t; p represents the parallel adit, and t represents the tunnel. For convenience, follow-up discussions also follow this agreement) ( Figure 3) [14]. The mapping function is: is the point of the mapping plane, and i is the imaginary unit.
Appl. Sci. 2022, 12, x FOR PEER REVIEW 5 of 18 the tunnel and parallel adit drainage [23]. The Schwartz alternating method transforms the multi-connected domain problem into multiple single-connected domain problems [24]. The interaction between these single connected domain problems is then iteratively considered. Using this method, the first step involves solving the seepage field for single tunnel drainage in a semi-infinite plane.

General Solution of Single Tunnel Drainage
The water table and tunnel (or parallel adit) circumferences in the Zj plane are conformally mapped to two circles with radius 1 and αj in the ζj plane by mapping functions (where j = p, t; p represents the parallel adit, and t represents the tunnel. For convenience, follow-up discussions also follow this agreement) ( Figure 3) [14]. The mapping function is: the point of the mapping plane, and i is the imaginary unit. From the mapping Function (5), the radius of the parallel adit and tunnel in the ζj plane is: Using conformal mapping, Equation (1) can be expressed in the ζj plane in polar coordinates as: where hj,k is the total head, and ρj and φj are the polar radius and polar angle in ζj plane, respectively. Subscripts j and k represent the serial number of the tunnels and the order of iterations, respectively. In the ζj plane, Equation (6) (the Laplace equation in polar coordinates) can be solved by using the method of separating variables [25]. Its general solution is shown in Equation (7). When the solution domain is an annulus, the coefficients in the general solution are easily obtained by the boundary conditions. From the mapping Function (5), the radius of the parallel adit and tunnel in the ζ j plane is: Using conformal mapping, Equation (1) can be expressed in the ζ j plane in polar coordinates as: where h j,k is the total head, and ρ j and ϕ j are the polar radius and polar angle in ζ j plane, respectively. Subscripts j and k represent the serial number of the tunnels and the order of iterations, respectively. In the ζ j plane, Equation (6) (the Laplace equation in polar coordinates) can be solved by using the method of separating variables [25]. Its general solution is shown in Equation (7). When the solution domain is an annulus, the coefficients in the general solution are easily obtained by the boundary conditions.
where E (j,k) , F (j,k) , A (j,k)m , C (j,k)m , B (j,k)m , D (j,k)m are coefficients, which can be determined by boundary conditions. The meanings of subscripts j and k are as described above, and the subscript m denotes the order of the Fourier series.
According to mapping Function (5), ρ j and ϕ j in the ζ j plane have the following relationship with x and y in the Z plane: In addition, the radial hydraulic gradient along the tunnel (or parallel adit) boundary needs to be transformed from the Z j plane to the ζ j plane in subsequent calculations, and the relationship between them is: | is the modulus of the derivative of the mapping function, and n j is the radial direction on the circumference of the tunnel (or parallel adit) in the Z j plane.

Iterative Methods
The flow chart of the Schwartz iterative solution is shown in Figure 4.
where E(j,k), F(j,k), A(j,k)m, C(j,k)m, B(j,k)m, D(j,k)m are coefficients, which can be determined by boundary conditions. The meanings of subscripts j and k are as described above, and the subscript m denotes the order of the Fourier series.
According to mapping Function (5), ρj and φj in the ζj plane have the following relationship with x and y in the Z plane: In addition, the radial hydraulic gradient along the tunnel (or parallel adit) boundary needs to be transformed from the Zj plane to the ζj plane in subsequent calculations, and the relationship between them is: is the modulus of the derivative of the mapping function, and nj i the radial direction on the circumference of the tunnel (or parallel adit) in the Zj plane.

Iterative Methods
The flow chart of the Schwartz iterative solution is shown in Figure 4.

The first iteration
Input Parameters:dw, d t , r t , r p , S h , Sv. Start Only the parallel adit is contained in the half-plane Z p and the radial hydraulic gradient at the boundary of the parallel adit is i p . The total head h p,1 and the additional radial hydraulic gradient i t,1 on the tunnel boundary are obtained.
Only the tunnel is contained in the half-plane Z t and the radial hydraulic gradient at the boundary of the tunnel is i t -i t,1 . The total head h t,1 and the additional radial hydraulic gradient i p,1 on the parallel adit boundary are obtained End Only the parallel adit is contained in the half-plane Z p and the radial hydraulic gradient at the boundary of the parallel adit is -i p,k-1 . The total head h p,k and the additional radial hydraulic gradient i t,k on the tunnel boundary are obtained.
Only the tunnel is contained in the half-plane Z t and the radial hydraulic gradient at the boundary of the tunnel is -i t,k . The total head h t,k and the additional radial hydraulic gradient i p,k on the parallel adit boundary are obtained Step 1 of the first iteration: Considering that the hydraulic gradient at the boundary of the parallel adit is i p (i p is the initial radial hydraulic gradient of the parallel adit boundary) and only the parallel adit is contained in the Z p plane, the total head h p,1 can be obtained using the method described in Section 3.1. The additional radial hydraulic gradient i t,1 at the tunnel boundary (virtual boundary) can then be determined using h p,1 .
Step 2 of the first iteration: Considering that the hydraulic gradient at the boundary of the tunnel is i t -i t,1 (i t is the initial radial hydraulic gradient of tunnel boundary) and only the tunnel is contained in the Z t plane, the total head h t,1 can also be obtained using the method described in Section 3.1. The additional radial hydraulic gradient i p,1 at the parallel adit boundary (virtual boundary) can then be determined using h t,1 .
After the first iteration, the seepage field obtained by summing the total head h p,1 and h t,1 fully satisfies the boundary conditions of the tunnel, but an additional radial hydraulic gradient i p,1 at the parallel adit boundary exists. If |i p,1 | max is less than the required value, the iteration ends; otherwise, the next iteration begins. The iteration termination condition is: Step 1 of the kth iteration: Considering that the hydraulic gradient at the boundary of the parallel adit is −i p,k−1 and only the parallel adit is contained in the Z p plane, the total head h p,k can be obtained. The additional radial hydraulic gradient i t,k at the tunnel boundary can then be determined using h p,k .
Step 2 of the kth iteration: Considering that the hydraulic gradient at the boundary of the tunnel is −i t,k , and only the tunnel is contained in the Z t plane, the total head h t,k can be obtained. The additional radial hydraulic gradient i p,k at the parallel adit boundary can then be determined using h t,k . The boundary conditions for the iteration steps are shown in Figure 5. After each iteration, Equation (11) is used to determine if the iteration continues.
of the parallel adit is ip (ip is the initial radial hydraulic gradient of the parallel adit boundary) and only the parallel adit is contained in the Zp plane, the total head hp,1 can be obtained using the method described in Section 3.1. The additional radial hydraulic gradient it,1 at the tunnel boundary (virtual boundary) can then be determined using hp,1.
Step 2 of the first iteration: Considering that the hydraulic gradient at the boundary of the tunnel is it-it,1 (it is the initial radial hydraulic gradient of tunnel boundary) and only the tunnel is contained in the Zt plane, the total head ht,1 can also be obtained using the method described in Section 3.1. The additional radial hydraulic gradient ip,1 at the parallel adit boundary (virtual boundary) can then be determined using ht,1.
After the first iteration, the seepage field obtained by summing the total head hp,1 and ht,1 fully satisfies the boundary conditions of the tunnel, but an additional radial hydraulic gradient ip,1 at the parallel adit boundary exists. If |ip,1|max is less than the required value, the iteration ends; otherwise, the next iteration begins. The iteration termination condition is: Step 1 of the kth iteration: Considering that the hydraulic gradient at the boundary of the parallel adit is −ip,k−1 and only the parallel adit is contained in the Zp plane, the total head hp,k can be obtained. The additional radial hydraulic gradient it,k at the tunnel boundary can then be determined using hp,k.
Step 2 of the kth iteration: Considering that the hydraulic gradient at the boundary of the tunnel is −it,k, and only the tunnel is contained in the Zt plane, the total head ht,k can be obtained. The additional radial hydraulic gradient ip,k at the parallel adit boundary can then be determined using ht,k. The boundary conditions for the iteration steps are shown in Figure 5. After each iteration, Equation (11) is used to determine if the iteration continues. As the iterative process continues, the additional hydraulic gradient around the parallel adit will decay, eventually reaching 0. After the iteration is terminated, the total water head obtained by each iteration is superimposed to obtain the final water head after n iterations: Step 1 of the first iteration Step 2 of the first iteration Step 1 of the kth iteration Step 2 of the kth iteration As the iterative process continues, the additional hydraulic gradient around the parallel adit will decay, eventually reaching 0. After the iteration is terminated, the total water head obtained by each iteration is superimposed to obtain the final water head after n iterations:

Calculation Process
Boundary conditions for step 1 of the first iteration (in the ζ p plane): At the water table: At the tunnel circumference: The boundary condition (14) can be expanded using the Fourier series: where a (p, After substituting Equation (7) into Equations (13) and (14), the coefficients can be obtained: Thus, h p,1 (x,y) can be obtained by substituting Equations (8) and (9) into h p,1 (ρ j ,ϕ j ). The additional radial hydraulic gradient i t,1 at the tunnel boundary can then be determined using h p,1 (x,y).
Boundary conditions for step 2 of the first iteration (in the ζ t plane): At the water table: At the tunnel circumference: The boundary condition (17) can be expanded by the Fourier series: where After substituting Equation (7) into Equations (16) and (17), the coefficients can be obtained: Substitute i p,1 into the iteration termination condition (11). If the termination condition is met, the iteration ends; otherwise, the next iteration begins.
Boundary conditions for step 1 of the kth (k = 2,3,4 . . . ) iteration (in the ζ p plane): At the water table: At the parallel adit circumference: Boundary conditions for step 2 of the kth (k = 2,3,4 . . . ) iteration (in the ζ t plane): At the water table: At the tunnel circumference: The coefficients E (j,k) , F (j,k) , A (j,k)m , C (j,k)m , B (j,k)m , and D (j,k)m can be obtained using the boundary condition Equations (19)- (22). After iterations are completed, all the water head h j,k in the iteration process are superposed to obtain the final total water head h(x,y). Finally, the pressure head around the tunnel is:

Comparison with Numerical Results
The Jiajiaoshan Tunnel (Chongqing, China) proposed for the China Chengdawan high-speed railway passes through the deep slow flow zone of soluble rock ( Figure 6). The permeability coefficient of the surrounding rock in the water-rich high-pressure section is 1.0 × 10 −6 -9.0 × 10 −6 m/s. The predicted maximum water pressure is 3 MPa and the maximum water inflow is 32.5 × 10 4 m 3 /d. A parallel adit is proposed to assist the construction of the tunnel and drain during operation to reduce the water pressure on the tunnel lining. The cross-sectional dimensions of the tunnel and parallel adit are shown in Figure 7. According to the equal permeable area principle (equal circumference of the outer profile of the section), the Jiajiaoshan Tunnel and the proposed parallel adit are equivalent to two circular sections with radii of 6.9 m and 3.9 m, respectively. The crown position of the tunnel remains unchanged in the equivalent process. Buried depth of tunnel and water table and position parameters of parallel adit relative to the tunnel are shown in Table 1.        The following calculations are based on the parameters of the Jiajiaoshan Tunnel (Table 1). In order to verify the accuracy of the semi-analytical results, a finite element model is established using ABAQUS (2020)  simulate a semi-infinite domain with a deeply buried tunnel and a parallel adit (Figure 8).      Figure 9 shows the variation of the maximum additional hydraulic gradient |i j,k | max at the parallel adit and tunnel boundary during the iteration. Horizontal coordinates represent the iterative steps, i.e., 1-1 represents the first step of the first iteration, 1-2 represents the second step of the first iteration. The additional radial hydraulic gradient decreases quickly, indicating that the analytical algorithm has a fast convergence speed. After two iterations, |i p,1 | max /i p = 0.00796%. Termination condition (11) is met, and the iteration ends. Theoretical calculations show that the convergence rate is related to the relative positions of the parallel adit and tunnel, burial depth, and drainage volume.
at the parallel adit and tunnel boundary during the iteration. Horizontal coordi represent the iterative steps, i.e., 1-1 represents the first step of the first iteration represents the second step of the first iteration. The additional radial hydraulic gra decreases quickly, indicating that the analytical algorithm has a fast convergence sp After two iterations, |ip,1|max/ip = 0.00796%. Termination condition (11) is met, and iteration ends. Theoretical calculations show that the convergence rate is related t relative positions of the parallel adit and tunnel, burial depth, and drainage volume     represents the second step of the first iteration. The additional radial hydraulic gr decreases quickly, indicating that the analytical algorithm has a fast convergence After two iterations, |ip,1|max/ip = 0.00796%. Termination condition (11) is met, a iteration ends. Theoretical calculations show that the convergence rate is related relative positions of the parallel adit and tunnel, burial depth, and drainage volum

Parametric Analysis
Parameters should be non-dimensionalized prior to performing parametric analysis. The pressure head at the tunnel boundary is non-dimensionalized by the pressure head that the tunnel and the parallel adit have when fully waterproofed. The values of S h and S v are non-dimensionalized by tunnel equivalent circle diameter D t . The drainage volume of the parallel adit and tunnel are non-dimensionalized by the maximum drainage volume Q j,max that they drainage alone. The maximum drainage Q j,max is obtained when the tunnel or the parallel adit is all fully drained, that is, the water pressure at their boundary is zero. The calculation method is as follows [14]: According to the parameters in Table 1, Q p,max = 157.9 m 3 ·(d·m) −1 , and Q t,max = 176.1 m 3 ·(d·m) −1 .
In order to more intuitively reflect the influence of parallel adit drainage on the water pressure at tunnel boundary, a parameter analysis is carried out for the undrained tunnel (Q t = 0) in Sections 5.1-5.3. Figure 11 illustrates the relationship between the pressure head of the tunnel crown and horizontal distance S h with different parallel adit discharges. The ordinate represents the relative pressure head at the crown of the tunnel (H tc is the pressure head of the tunnel crown when the collaborative drainage system works, and H 0 tc is the pressure head at the crown when the tunnel and parallel adit are all fully waterproofed). The horizontal distance S h is linearly positively correlated with the tunnel crown pressure head. For the three different parallel adit drainage conditions (Q p = 0.95 Q p,max , Q p = 0.75 Q p,max , Q p = 0.55 Q p,max ), S h changes from 1.9 D t to 0.6 D t , and the relative pressure head decreases from 0.47, 0.59, 0.70 to 0.38, 0.51, 0.64, respectively. The pressure head on the lining can be reduced to a certain extent by reducing S h ; however, considering the construction conditions and surrounding rock stability, the distance between the parallel adit and tunnel should be no less than the tunnel diameter.

Effect of Horizontal Distance between the Parallel Adit and Tunnel
that the tunnel and the parallel adit have when fully waterproofed. The values of Sv are non-dimensionalized by tunnel equivalent circle diameter Dt. The drainage v of the parallel adit and tunnel are non-dimensionalized by the maximum drainage v Qj,max that they drainage alone. The maximum drainage Qj,max is obtained when the or the parallel adit is all fully drained, that is, the water pressure at their boundary i The calculation method is as follows [14]: According to the parameters in Table 1, Qp,max = 157.9 m 3 ·(d·m) −1 , and Qt,max = m 3 ·(d·m) −1 .
In order to more intuitively reflect the influence of parallel adit drainage on the pressure at tunnel boundary, a parameter analysis is carried out for the undrained (Qt = 0) in Sections 5.1-5.3. Figure 11 illustrates the relationship between the pressure head of the tunnel and horizontal distance Sh with different parallel adit discharges. The ordinate repr the relative pressure head at the crown of the tunnel (Htc is the pressure head of the crown when the collaborative drainage system works, and H 0 tc is the pressure head crown when the tunnel and parallel adit are all fully waterproofed). The hori distance Sh is linearly positively correlated with the tunnel crown pressure head. F three different parallel adit drainage conditions (Qp = 0.95Qp,max, Qp = 0.75Qp,max 0.55Qp,max), Sh changes from 1.9Dt to 0.6Dt, and the relative pressure head decrease 0.47, 0.59, 0.70 to 0.38, 0.51, 0.64, respectively. The pressure head on the lining reduced to a certain extent by reducing Sh; however, considering the constr conditions and surrounding rock stability, the distance between the parallel ad tunnel should be no less than the tunnel diameter. Figure 11. Relative pressure head on the tunnel crown versus horizontal distance betwe parallel adit and tunnel. Figure 12 shows the influence of the relative height difference between the p adit and tunnel on the pressure head of the tunnel crown. The relative height diff Sv changes within the conventional range, and the relative pressure head of the crown change slightly, indicating that the variation of Sv in a certain range has little on the water pressure on the lining. Therefore, Sv is not the primary factor in the  Figure 11. Relative pressure head on the tunnel crown versus horizontal distance between the parallel adit and tunnel. Figure 12 shows the influence of the relative height difference between the parallel adit and tunnel on the pressure head of the tunnel crown. The relative height difference S v changes within the conventional range, and the relative pressure head of the tunnel crown change slightly, indicating that the variation of S v in a certain range has little effect on the water pressure on the lining. Therefore, S v is not the primary factor in the design of the drainage scheme. However, when the tunnel encounters large-scale water inflow, water needs to be introduced into the parallel adit through the horizontal channel between the tunnel and the parallel adit. Therefore, the bottom of the parallel adit should be lower than the bottom of the tunnel drainage ditch.

Effect of the Height Difference between the Parallel Adit and Tunnel
of the drainage scheme. However, when the tunnel encounters large-scale w water needs to be introduced into the parallel adit through the horizontal chan the tunnel and the parallel adit. Therefore, the bottom of the parallel adit sho than the bottom of the tunnel drainage ditch.  Figure 13 shows the pressure head contour map when drainage prima the parallel adit. There is a clear pressure drop funnel around the parallel influence radiates to the surrounding area, which releases the water pressur tunnel. However, the water pressure around the tunnel is unevenly released accurately evaluate this inhomogeneity, Figure 14 shows the relative pr around the tunnel. The water pressure on the side of the tunnel facing the pa significantly lower than that on the other side. The greater the discharge of adit, the more significant the heterogeneity. When the parallel adit drainage the pressure heads of the four points A, B, C, and D at different parts of the t are reduced to 0.42, 0.28, 0.40, and 0.49 times of their initial water heads, Point B has the largest pressure release, which is 1.41 times that of point smallest pressure release. In addition, the closer to the parallel adit, the greate of the contour lines and the greater the pressure gradient, indicating that between the parallel adit and the tunnel decreases, and the inhomogeneity pressure around the tunnel increases ( Figure 13).   Figure 13 shows the pressure head contour map when drainage primarily relies on the parallel adit. There is a clear pressure drop funnel around the parallel adit, and its influence radiates to the surrounding area, which releases the water pressure around the tunnel. However, the water pressure around the tunnel is unevenly released. In order to accurately evaluate this inhomogeneity, Figure 14 shows the relative pressure head around the tunnel. The water pressure on the side of the tunnel facing the parallel adit is significantly lower than that on the other side. The greater the discharge of the parallel adit, the more significant the heterogeneity. When the parallel adit drainage Q p = 0.95 Q p , the pressure heads of the four points A, B, C, and D at different parts of the tunnel lining are reduced to 0.42, 0.28, 0.40, and 0.49 times of their initial water heads, respectively. Point B has the largest pressure release, which is 1.41 times that of point D, with the smallest pressure release. In addition, the closer to the parallel adit, the greater the density of the contour lines and the greater the pressure gradient, indicating that the distance between the parallel adit and the tunnel decreases, and the inhomogeneity of the water pressure around the tunnel increases ( Figure 13).  Figure 13 shows the pressure head contour map when drainage pr the parallel adit. There is a clear pressure drop funnel around the para influence radiates to the surrounding area, which releases the water pre tunnel. However, the water pressure around the tunnel is unevenly rele accurately evaluate this inhomogeneity, Figure 14 shows the relativ around the tunnel. The water pressure on the side of the tunnel facing th significantly lower than that on the other side. The greater the discharg adit, the more significant the heterogeneity. When the parallel adit drain the pressure heads of the four points A, B, C, and D at different parts of are reduced to 0.42, 0.28, 0.40, and 0.49 times of their initial water hea Point B has the largest pressure release, which is 1.41 times that of p smallest pressure release. In addition, the closer to the parallel adit, the gr of the contour lines and the greater the pressure gradient, indicating between the parallel adit and the tunnel decreases, and the inhomogen pressure around the tunnel increases ( Figure 13).

Effect of Drainage Volume
When the tunnel drainage volume is close to Qt,max, th also close to 0. At this time, the drainage of the parallel capacity of the tunnel drainage system design is insuffici drainage volume is less than Qt,max, and parallel adit drain Therefore, in this section, considering the tunnel drainage 0.75Qt,max, the influence of the parallel adit drainage on the Keeping the tunnel drainage volume unchanged, the adit is linearly negatively correlated with the tunnel correlation coefficient (Figure 15). Compared with positio drainage volume on water pressure is more significant. Wh = 0.25Qt,max, the relative pressure head of the tunnel crown 0.42 and 0.16, respectively, when parallel adit drainage water pressure on the tunnel crown decreases by approxim only relying on tunnel drainage, which is a significant re tunnel drainage is Qt = 0.50Qt,max, Qt = 0.75Qt,max, the rate of same as the above two conditions. With increasing dra pressure head continues to decrease until it reaches 0.

Effect of Drainage Volume
When the tunnel drainage volume is close to Q t,max , the pressure head on the lining is also close to 0. At this time, the drainage of the parallel adit is unnecessary. When the capacity of the tunnel drainage system design is insufficient, old, or blocked, the tunnel drainage volume is less than Q t,max , and parallel adit drainage will have a positive effect. Therefore, in this section, considering the tunnel drainage for Q t = 0, 0.25 Q t,max , 0.50 Q t,max , 0.75 Q t,max , the influence of the parallel adit drainage on the water pressure is analyzed.
Keeping the tunnel drainage volume unchanged, the drainage volume of the parallel adit is linearly negatively correlated with the tunnel pressure head with a strong correlation coefficient ( Figure 15). Compared with position parameters, the influence of drainage volume on water pressure is more significant. When tunnel drainage is Q t = 0, Q t = 0.25 Q t,max , the relative pressure head of the tunnel crown decreases from 1.0 and 0.74 to 0.42 and 0.16, respectively, when parallel adit drainage volume reaches the limit. The water pressure on the tunnel crown decreases by approximately 58% compared with that only relying on tunnel drainage, which is a significant release in water pressure. When tunnel drainage is Q t = 0.50 Q t,max , Q t = 0.75 Q t,max , the rate of water pressure reduction is the same as the above two conditions. With increasing drainage of the parallel adit, the pressure head continues to decrease until it reaches 0.

Effect of Drainage Volume
When the tunnel drainage volume is close to Qt,max, the pressure head on also close to 0. At this time, the drainage of the parallel adit is unnecessar capacity of the tunnel drainage system design is insufficient, old, or blocked drainage volume is less than Qt,max, and parallel adit drainage will have a po Therefore, in this section, considering the tunnel drainage for Qt = 0, 0.25Qt,m 0.75Qt,max, the influence of the parallel adit drainage on the water pressure is Keeping the tunnel drainage volume unchanged, the drainage volume o adit is linearly negatively correlated with the tunnel pressure head w correlation coefficient ( Figure 15). Compared with position parameters, the drainage volume on water pressure is more significant. When tunnel drainag = 0.25Qt,max, the relative pressure head of the tunnel crown decreases from 1. 0.42 and 0.16, respectively, when parallel adit drainage volume reaches th water pressure on the tunnel crown decreases by approximately 58% compar only relying on tunnel drainage, which is a significant release in water pre tunnel drainage is Qt = 0.50Qt,max, Qt = 0.75Qt,max, the rate of water pressure red same as the above two conditions. With increasing drainage of the para pressure head continues to decrease until it reaches 0.   Figure 15. Relative pressure head of the tunnel crown versus drainage volume of the parallel adit.

Pressure Head Fitting Formula
Under the given drainage boundary conditions, the pressure head at the crown of the tunnel is linearly related to the drainage volume of the parallel adit and tunnel. Therefore, the empirical formula for the calculation of the water pressure around the tunnel can be fitted using the semi-analytical results, so that the calculation and parameter analysis of the drainage scheme can be achieved without the iterative process. The fitting formula can be expressed as (which still takes the crown of the tunnel as an example): where δ t and δ p are the influence coefficient of the tunnel and parallel adit drainage on tunnel crown pressure head, respectively, and the dimension is 1. They are determined using the variables r p , r t , d tw , S h , and S v , where d tw = d t − d w is the underwater depth of tunnel crown; other symbols are defined as before. The expressions of the influence coefficients δ t and δ p are obtained using multivariate nonlinear regression, that is, the relationship between a single variable and δ t and δ p is analyzed one by one (other parameters remain unchanged) to obtain the optimal one-variable fitting function [26]. The one-variable fitting functions are linearly superposed to obtain the multivariate fitting function. According to the above method, the optimal fitting functions of δ t and r p , r t , d tw , S h , and S v are a linear function, logarithmic function, logarithmic function, quadratic polynomial function, and quadratic polynomial function, respectively. The optimal univariate fitting functions of δ p and r p , r t , d tw , S h , and S v are quadratic polynomial function, linear function, logarithmic function, linear function, and linear function, respectively. Set the fitting function of δ t , δ p as (26), (27). The undetermined coefficients of Formulas (26) and (27) are obtained using the orthogonal test to reduce the amount of calculation. The level of each factor is shown in Table 2, and the design scheme and results of the orthogonal test are shown in Table 3. δ t = a 0 + a 1 r p + a 2 ln(r t ) + a 3 ln(d tw ) + a 4 S h 2 + a 5 S h + a 6 S v 2 + a 7 S v (26) where a 0 , a 1 ... a 7 , b 0 , b 1 ... b 6 are undetermined coefficients. Based on the data in Table 3, SPSS software was used for multivariate nonlinear regression, and the fitting formulas of δ t and δ p were obtained, such as (28) and (29). Their R 2 was 0.99 and 0.98, respectively, indicating that the fitting results were satisfactory. δ t = 1.3689 × 10 −1 + 1.9872 × 10 −3 r p − 1.5918 × 10 −1 ln(r t ) + 1.5280 × 10 −1 ln(d tw ) +2.9999 × 10 −6 S h 2 + 9 × 10 −5 S h − 1.1520 × 10 −3 S v 2 + 4.8384 × 10 −3 S v (28) δ p = −1.9872 × 10 −1 + 4.0954 × 10 −3 r p 2 − 3.3273 × 10 −2 r p − 1.1362 × 10 −3 r t +1.5552 In order to test the accuracy of the fitting formula, a set of parameters were randomly selected in the fitting range (r p = 3.0 m, r t = 7.2 m, d tw = 300 m, S h = 18 m, S v = 1.0 m, K = 4 × 10 −6 m/s) and calculated. The fitting results are compared with the semi-analytical results ( Figure 16). The fitting results are in good agreement with the theoretical calculation results, and their deviations are within 0.8%, showing that the fitting Formula (25) is suitable for the calculation and parameter analysis of the drainage scheme for tunnel and parallel adit. Based on the data in Table 3, SPSS software was used for multivariate regression, and the fitting formulas of δt and δp were obtained, such as (28) and ( R 2 was 0.99 and 0.98, respectively, indicating that the fitting results were satisfa

Conclusions
1. The theoretical model and semi-analytical algorithm of the collaborative scheme for a tunnel and parallel adit proposed here can accurately meet drainage boundary conditions in several iterations. The calculation re consistent with the numerical simulation results, which are suitable calculation and analysis of drainage schemes. 2. The drainage volume of the parallel adit negatively correlates with t pressure at the tunnel boundary with a strong correlation coefficient, influence on water pressure is the most significant. The horizontal distance parallel adit and tunnel has a linear positive correlation with water press tunnel boundary; however, its influence is less than drainage volume. T
The theoretical model and semi-analytical algorithm of the collaborative drainage scheme for a tunnel and parallel adit proposed here can accurately meet the preset drainage boundary conditions in several iterations. The calculation results are consistent with the numerical simulation results, which are suitable for the calculation and analysis of drainage schemes.

2.
The drainage volume of the parallel adit negatively correlates with the water pressure at the tunnel boundary with a strong correlation coefficient, and the influence on water pressure is the most significant. The horizontal distance between parallel adit and tunnel has a linear positive correlation with water pressure at the tunnel boundary; however, its influence is less than drainage volume. The water pressure of the lining facing the parallel adit side is significantly less than that on the other side when the parallel adit is collaboratively drained. The smaller the horizontal distance between parallel adit and tunnel, the larger the drainage volume of the parallel adit, and the more significant the inhomogeneity. Therefore, in order to reduce water pressure uniformly and meet the requirements of surrounding rock stability, the horizontal distance between the parallel adit and tunnel should be no less than the tunnel diameter. The height difference between the parallel adit and the tunnel has little effect on the water pressure on the lining; however, it only needs to meet other construction and operation requirements.

3.
Based on the semi-analytical results, the pressure head formula of tunnel crown fitted by using nonlinear regression theory has a quality fit, which can be used to calculate and analyze the parameters of the collaborative drainage scheme for a tunnel and parallel adit. The formula can also intuitively express the contribution of water pressure release to each parameter.