# Analytical Solution of Seepage Field in Karst Tunnel

Authors: Chong Jiang; Han-song Xie; Jia-li He; Wen-yan Wu; Zhi-chao Zhang
Journal: Advances in Civil Engineering
 (2018)

---

## Abstract

An analytical solution for the seepage field in water-filled karst tunnel is derived based on the inversion of complex function and groundwater hydraulics theory. The solution considers the distance between the tunnel and the cavern, the size of the cavern, and the properties of the lining structure, such as the permeability coefficient as well as the radius of the grouting ring. This paper also performed numerical simulations for two cases: the application of gravity and the absence of gravity. The numerical solution was obtained to verify the analytical solution, and a good agreement was found. Then, the effect of parameters is discussed in detail, including the distance between the tunnel and the cavern, the radius of the cavern, the grouting ring, and the initial support. The results show that when the radius of the cavern is constant, the pressure head and seepage flow decrease as the distance between the tunnel and the cavern increases. When the distance is constant, the pressure head and seepage flow increase with the increase of the radius of the cavern. In addition, the pressure head and the seepage flow decrease with the increase of the thickness of the grouting ring and decrease with the decrease of the permeability coefficient. As the thickness of the initial support increases, the pressure head gradually increases and the percolation decreases. Furthermore, due to the great influence of the grouting ring and initial support on the pressure head and seepage flow, the thickness and permeability coefficient of the grouting ring and initial support should be taken into account carefully during construction.


---

## Body

## 1. Introduction

In recent years, with the development of the Belt and Road strategy, the construction of tunnels has developed rapidly. And it is possible to build a tunnel in the karst zone in southwestern China. The karst water has a great influence on the construction of the tunnel. Its seepage causes karst tunnel water gushing, and it is easy to cause tunnel deformation and other disasters when the water volume is too large. Therefore, it is important to study the distribution of seepage field in karst tunnels.The analysis of seepage field in karst tunnels generally uses analytical method and numerical simulation. Based on well mapping theory and superposition principle, Zhang [1] deduced the potential function and distribution of water head in the karst tunnel, but did not study the influence of karst water on tunnel lining structure. Zhang [2] derived the analytical solution of the seepage field of a semiinfinite parallel double-hole tunnel through the seepage mechanics and image method. Yang et al. [3] studied the effect of seismic loading on the double-sided slope of high water content sandy clay. Wu [4] deduced the analytic solution of the complex variable function for the distribution of high water pressure seepage field according to the conformal mapping theory. Park et al. [5] compared two existing analytical solutions for the steady-state groundwater inflow into a drained circular tunnel in a semiinfinite aquifer. Some researchers used the theory of hydraulics and complex functions to solve the analytical solution of seepage field of underwater tunnel composed of surrounding rocks and lining structures [6–9]. Based on the theory of the growth of pore pressure caused by linear waves, Xiong et al. [10] developed an analytical solution of the seepage field of a tunnel affected by the wave in the semiinfinite aquifer. Ma et al. [11, 12] carried out some experimental investigations on hydraulic characteristics and seepage characteristics. Lei [13] acquired an analytical solution for steady flow into a tunnel based on the image method. Huangfu et al. [14] studied analytical solutions for steady groundwater flowing into a horizontal tunnel in the semiinfinite aquifer.However, most of the aforementioned literature have studied high hydraulic pressure tunnels or high hydraulic parallel tunnels in a semiinfinite plane. There is less analytical formula for the distribution of seepage field in karst tunnels. Therefore, based on the inversion of complex function and groundwater hydraulics theory, this paper deduces the distribution of seepage field in karst tunnels and calculates the amount of seepage. The formula related to pressure head and tunnel seepage is deduced from the relationship between the grouting ring, initial support, and secondary lining. The results are compared with the numerical solutions of the two cases of gravity and no gravity. The change of lining structure characteristics such as grouting circle and initial support is discussed in detail. It is helpful for us to have a deeper understanding of underground water seepage in karst tunnels, and to take effective measures to reduce the harmful effects of karst water on tunnels, so as to provide references for such problems in the future.


## 2. Calculation Model

The schematic diagram of the karst tunnel in an infinite plane is shown in Figure1. Here, q is the water pressure in the karst cavern, a is the tunnel radius, and d is the distance between the center of the tunnel and the karst cavern. Additionally, the basic assumptions are as follows:Figure 1
Schematic diagram of karst tunnel in infinite plane. A: tunnel; B: surrounding rocks; C: cavern.(1)
The surrounding rock of the tunnel is homogeneous with isotropic permeability [5]
(2)
The aquifer and karst water are incompressible. The flow is in a steady state and is governed by Darcy’s law [13]
(3)
Due to the long distance from the tunnel to the surface, it is assumed that the plane where the karst tunnel is located is an infinite plane, and the seepage field is only affected by karst water
(4)
When studying the lining structure, the pore pressure is constant on the same ring, and the water-filling hole is equal to the water pressure head [7]


## 3. Analytical Solution of Seepage Field

### 3.1. Analytical Method

Considering the characteristics of karst tunnels, the inversion of complex functions and groundwater hydraulics is used to analyze the seepage field of karst tunnels. The solution steps are as follows:(a)
According to the complex potential function of the point source and the inverse transformation relation, the mirror image pointsa2/z0¯ of z0 in the tunnel are obtained, and the complex potential function at any place in the seepage plane of the karst tunnel is solved(b)
Through the resetting potential at any place in the seepage plane of the karst tunnel, the potential function, flow function, and complex velocity of the seepage field in the karst tunnel are solved(c)
The seepage flow and pressure head are obtained by the boundary conditions, and the head pressure and the seepage flow in each lining structure were solved according to the relationship of the lining structures


### 3.2. Analysis of Complex Potential Function

The water-filled cavern existing outside the tunnel is treated as a constant point source, and the inversion of the cavern and tunnel is shown in Figure2. In the infinite plane where the tunnel radius is a, there is a cavern with the water pressure q at z0=x0+iy0 outside the tunnel. The cavern is considered as a point source. After the inversion transformation, a cavern with the same intensity is generated at a2/z0¯, which is calculated as a point sink. The complex potential function of the point z=x+iy in the infinite plane of the point source z0=x0+iy0 is(1)fz=q2πlnz=q2πlnaeiθ,where the potential function is ϕ=q/2πlna and the flow function is φ=q/2πθ.Figure 2
Schematic diagram of inversion of caverns and karst tunnels.Based on the complex potential function of the point source in the infinite plane and the inversion transformation relationship, the conjugate complex potential function of the point sink is obtained:(2)f¯a2z=−q2πlnz−a2z0¯.The point source in the unbounded fluid has a reset potential offz, and the conjugate reset potential of the point converted according to the inversion is f¯a2/z. All the points are in the region of z>a outside the tunnel. In the flow field behind the tunnel z=a, the tunnel wall can be a resetting potential of a streamline. According to the superposition principle of resetting potential in seepage mechanics, the complex potential function of the point seeping in the infinite plane of the tunnel is obtained:(3)Fz=fz+f¯a2z.


### 3.3. Analysis of Seepage Field

According to the complex potential formula (3) of the point in the infinite plane of the tunnel, formulas (1) and (2) are plugged into (3), as shown in the following:(4)Fz=q2πlnz−z0−q2πlnz−a2z0¯.For simplicity, without loss of generality, the point source can be selected to the right of the horizontal direction of the tunnel; that is, thex-axis passes through the point source. Then, there is z0=z0¯=x0≡d, as shown in Figure 3. Therefore, Formula (4) can be written as(5)Fz=q2πlnz−d−q2πlnz−a2d,=q2πlnx−d2+y2−iq2πarctgyx−d,−q2πlnx−a2d2+y2+iq2πarctgyx−a2/d.Figure 3
Schematic diagram of the relationship between each lining structure.Complex velocity:(6)Wz=dFdz=q/2πlnz−d−q/2πlnz−a2/ddz,=q2π1z−d−1z−a2/d.The potential function of the point source in the infinite plane of the tunnel:(7)ϕ=q4πlnx−d2+y2x−a2/d2+y2.The flow function of the point source in the infinite plane of the tunnel:(8)φ=−q2πarctgyx−d+q2πarctgyx−a2/d.Boundary conditions:Cavern boundary:(9)Φw=−q2πlnrwl.Tunnel boundary:(10)Φs=−q2πlnda,where l is the distance between the mirror point a2/z0¯ and z0, rw is the radius of the cavern, Φw is the cavern potential function, and Φs is the tunnel potential function.According to the boundary conditions and the steady-state percolation of the single-phase liquid [15], it can be obtained that(11)q=2πΦw−Φs−lnarw/ld=2πΦw−Φslnd2−a2/arw,Hw=pwρg,Hs=psρg,Q=2πKρgHw−Hslnd2−a2/arw,where K is the permeability, Pw and Ps are the head pressures of the cavern and tunnel, respectively, Q is the unit permeate flow, Hw is the head pressure of the cavern, and Hs is the head pressure of the tunnel boundary, that is, the head pressure of the surface of the grouting ring.


### 3.4. Analysis of Seepage Field of Lining Structure

The seepage field in the lining and grouting circle cannot be solved by the complex function, but it can be solved by using the groundwater seepage mechanics and the shaft theory [6]. The relationship between the initial support, the secondary lining, and the grouting ring pressure head is shown in Figure 3, and the formulas are as follows:(12)Hs−H1=Q2πk1lnar1,H1−H2=Q2πk2lnr1r2,H2−H3=Q2πk3lnr2r3,where Hs is the head pressure of the surface of the grouting ring; H1, H2,and H3 are the head pressures of the joints of the grouting ring and the initial support, the initial support and the secondary lining joint, and the secondary lining inner surface, respectively; r1, r2, and r3 are the radius of the grouting ring, initial support, and internal of secondary lining, respectively; k1, k2, and k3 are the permeability coefficients of the grouting circle, initial support, and secondary lining, respectively.Because of the same seepage, the above equations can be integrated to obtain the tunnel seepage when the grouting circle, the initial support, and the two lining are completed:(13)Q=2πHw−H31/k2lnr1/r2+1/k1lna/r1+1/k3lnr2/r3+1/klnd2−a2/arw.Head pressure of initial support and secondary lining joint:(14)H2=H3+Hw−H3k3Alnr2r3.Head pressure of grouting circle and initial support connection:(15)H1=H3+Hw−H3k3Alnr2r3+Hw−H3k2Alnr1r2.In the formula,(16)A=1k2lnr1r2+1k1lnar1+1k3lnr2r3+1klnd2−a2arw,where k is the permeability coefficients of the surrounding rock.According to different stages of construction such as tunnel holes and grouting, the different water pressure distribution can be obtained. When tunnel grouting and initial support are completed, the tunnel seepage is(17)Q=2πHw−H21/k2lnr1/r2+1/k1lnr/r1+1/klnd2−a2/arw.Head pressure of grouting circle and initial support connection:(18)H1=H2+Hw−H2lnr1/r2lnr1/r2+k2/k1lnr/r1+k2/klnd2−a2/arw.


## 3.1. Analytical Method

Considering the characteristics of karst tunnels, the inversion of complex functions and groundwater hydraulics is used to analyze the seepage field of karst tunnels. The solution steps are as follows:(a)
According to the complex potential function of the point source and the inverse transformation relation, the mirror image pointsa2/z0¯ of z0 in the tunnel are obtained, and the complex potential function at any place in the seepage plane of the karst tunnel is solved(b)
Through the resetting potential at any place in the seepage plane of the karst tunnel, the potential function, flow function, and complex velocity of the seepage field in the karst tunnel are solved(c)
The seepage flow and pressure head are obtained by the boundary conditions, and the head pressure and the seepage flow in each lining structure were solved according to the relationship of the lining structures


## 3.2. Analysis of Complex Potential Function

The water-filled cavern existing outside the tunnel is treated as a constant point source, and the inversion of the cavern and tunnel is shown in Figure2. In the infinite plane where the tunnel radius is a, there is a cavern with the water pressure q at z0=x0+iy0 outside the tunnel. The cavern is considered as a point source. After the inversion transformation, a cavern with the same intensity is generated at a2/z0¯, which is calculated as a point sink. The complex potential function of the point z=x+iy in the infinite plane of the point source z0=x0+iy0 is(1)fz=q2πlnz=q2πlnaeiθ,where the potential function is ϕ=q/2πlna and the flow function is φ=q/2πθ.Figure 2
Schematic diagram of inversion of caverns and karst tunnels.Based on the complex potential function of the point source in the infinite plane and the inversion transformation relationship, the conjugate complex potential function of the point sink is obtained:(2)f¯a2z=−q2πlnz−a2z0¯.The point source in the unbounded fluid has a reset potential offz, and the conjugate reset potential of the point converted according to the inversion is f¯a2/z. All the points are in the region of z>a outside the tunnel. In the flow field behind the tunnel z=a, the tunnel wall can be a resetting potential of a streamline. According to the superposition principle of resetting potential in seepage mechanics, the complex potential function of the point seeping in the infinite plane of the tunnel is obtained:(3)Fz=fz+f¯a2z.


## 3.3. Analysis of Seepage Field

According to the complex potential formula (3) of the point in the infinite plane of the tunnel, formulas (1) and (2) are plugged into (3), as shown in the following:(4)Fz=q2πlnz−z0−q2πlnz−a2z0¯.For simplicity, without loss of generality, the point source can be selected to the right of the horizontal direction of the tunnel; that is, thex-axis passes through the point source. Then, there is z0=z0¯=x0≡d, as shown in Figure 3. Therefore, Formula (4) can be written as(5)Fz=q2πlnz−d−q2πlnz−a2d,=q2πlnx−d2+y2−iq2πarctgyx−d,−q2πlnx−a2d2+y2+iq2πarctgyx−a2/d.Figure 3
Schematic diagram of the relationship between each lining structure.Complex velocity:(6)Wz=dFdz=q/2πlnz−d−q/2πlnz−a2/ddz,=q2π1z−d−1z−a2/d.The potential function of the point source in the infinite plane of the tunnel:(7)ϕ=q4πlnx−d2+y2x−a2/d2+y2.The flow function of the point source in the infinite plane of the tunnel:(8)φ=−q2πarctgyx−d+q2πarctgyx−a2/d.Boundary conditions:Cavern boundary:(9)Φw=−q2πlnrwl.Tunnel boundary:(10)Φs=−q2πlnda,where l is the distance between the mirror point a2/z0¯ and z0, rw is the radius of the cavern, Φw is the cavern potential function, and Φs is the tunnel potential function.According to the boundary conditions and the steady-state percolation of the single-phase liquid [15], it can be obtained that(11)q=2πΦw−Φs−lnarw/ld=2πΦw−Φslnd2−a2/arw,Hw=pwρg,Hs=psρg,Q=2πKρgHw−Hslnd2−a2/arw,where K is the permeability, Pw and Ps are the head pressures of the cavern and tunnel, respectively, Q is the unit permeate flow, Hw is the head pressure of the cavern, and Hs is the head pressure of the tunnel boundary, that is, the head pressure of the surface of the grouting ring.


## 3.4. Analysis of Seepage Field of Lining Structure

The seepage field in the lining and grouting circle cannot be solved by the complex function, but it can be solved by using the groundwater seepage mechanics and the shaft theory [6]. The relationship between the initial support, the secondary lining, and the grouting ring pressure head is shown in Figure 3, and the formulas are as follows:(12)Hs−H1=Q2πk1lnar1,H1−H2=Q2πk2lnr1r2,H2−H3=Q2πk3lnr2r3,where Hs is the head pressure of the surface of the grouting ring; H1, H2,and H3 are the head pressures of the joints of the grouting ring and the initial support, the initial support and the secondary lining joint, and the secondary lining inner surface, respectively; r1, r2, and r3 are the radius of the grouting ring, initial support, and internal of secondary lining, respectively; k1, k2, and k3 are the permeability coefficients of the grouting circle, initial support, and secondary lining, respectively.Because of the same seepage, the above equations can be integrated to obtain the tunnel seepage when the grouting circle, the initial support, and the two lining are completed:(13)Q=2πHw−H31/k2lnr1/r2+1/k1lna/r1+1/k3lnr2/r3+1/klnd2−a2/arw.Head pressure of initial support and secondary lining joint:(14)H2=H3+Hw−H3k3Alnr2r3.Head pressure of grouting circle and initial support connection:(15)H1=H3+Hw−H3k3Alnr2r3+Hw−H3k2Alnr1r2.In the formula,(16)A=1k2lnr1r2+1k1lnar1+1k3lnr2r3+1klnd2−a2arw,where k is the permeability coefficients of the surrounding rock.According to different stages of construction such as tunnel holes and grouting, the different water pressure distribution can be obtained. When tunnel grouting and initial support are completed, the tunnel seepage is(17)Q=2πHw−H21/k2lnr1/r2+1/k1lnr/r1+1/klnd2−a2/arw.Head pressure of grouting circle and initial support connection:(18)H1=H2+Hw−H2lnr1/r2lnr1/r2+k2/k1lnr/r1+k2/klnd2−a2/arw.


## 4. Verification and Discussion

### 4.1. Numerical Verification

Based on the complex function and groundwater hydraulics theory, the distribution of seepage field suitable for homogeneous and isotropic karst tunnels and lining structures were deduced. In order to ensure the accuracy of the analytical solution of the karst tunnel, the numerical solution simulated by the finite difference software FLAC3D is compared with the theoretical solution obtained by the method used in this paper.The calculation conditions are listed as follows: the head pressure of the cavern is 54 m; the radius of the cavern is 4 m; the radius of the tunnel is 7.25 m; the radius of the grouting circle is 2.25 m; and the initial support radius is 0.25 m; the permeability coefficient of surrounding rocks is 1.5 × 10−6 cm·s−1; the permeability coefficient of grouting ring is 10−7 cm·s−1; the initial support permeability coefficient of initial support is 10−8 cm·s−1.Six groups of 100 m × 100 m models of caverns on the right side of the tunnel were established and calculated with FLAC3D. The pore water pressure of the tunnel with gravity and without gravity was obtained, respectively. The six groups of models are d=15m, d=20m, d=25m, d=30m, d=35m, and d=40m. When d=20m, the karst tunnel model consists of 4236 elements, as shown in Figure 4. Figures 5 and 6 are contours of zone pore pressure with gravity and without gravity, respectively.Figure 4
Finite element model with circular karst tunnel.Figure 5
Contour of pore water pressure with gravity ford=20, a=7.25, and rw=4.Figure 6
Contour of pore water pressure without gravity ford=20, a=7.25, and rw=4.The pore water pressures extracted from the six groups of numerical simulation data in Figure5 with gravity and Figure 6 without gravity are compared with the theoretical solutions. The results are shown in Figure 7. The law obtained by numerical and theoretical solutions is that when the grouting ring and the initial support are completed and the secondary lining is not yet completed, H1 decreases with the increase of d. Moreover, the curves have a good consistency. When d=15m, the analytical solution for H1 is 25.78 m, the numerical solution with gravity is 25.67 m, and the numerical solution without gravity is 25.77 m. Therefore, the analytical method in this paper is suitable for solving the distribution of seepage field in karst tunnels.Figure 7
Comparison of numerical solution and theoretical solution fora=7.25 and rw=4.


### 4.2. The Effect of Distance and Cavern Radius

Figure8 shows the relationship between H1 and d. Figure 9 shows the relationship between Q and rw. In the infinite plane, when Hs is constant, the main factors affecting H1 and Q in the tunnel are rw and d. When rw is constant, H1 and Q decrease with the increase of d. When d is constant, H1 and Q gradually increase as rw increases.Figure 8
Relationship between head pressure and distance.Figure 9
Relationship between seepage flow and cavern radius.


### 4.3. The Effect of Grouting Ring

When the diameter of the grouting ring and the initial support radius is constant, the tunnel radius and the thickness of the grouting ring increase. Figures10–13 show that H1 and Q gradually decrease with the increase of the thickness of the grouting ring when the ratio between the distance and the permeability coefficient is constant. When the thickness of the grouting ring was 0.25 m and a reached 6 m, Q and H1 decreased to 41.94% and 41.86%, respectively. When a is constant, H1 and Q gradually decrease as the permeability coefficient decreases. When a increases to 7 m, the rate of H1 and Q decrease gradually. From Figures 12 and 13, when a is constant, H1 becomes smaller as the seepage coefficient of the grouting ring decreases. When kr=k/k1 increases to 80, Q decreases to 22.4% of Q when kr=5, and the Q reduction rate gradually decreases. In summary, increasing the thickness of the grouting ring and reducing the permeability coefficient of the grouting ring can reduce Q and H1. However, excessively increasing the thickness of the grouting ring will not only increase the cost but also reduce the impact on the pressure head and seepage flow. Therefore, a reasonable grouting material and construction process should be selected to ensure the thickness of the grouting ring when the permeability coefficient of the grouting ring is less than the specified value.Figure 10
Effect ofd on the relationship between pressure head and tunnel radius.Figure 11
Relationship between seepage flow and tunnel radius.Figure 12
Effect ofkr on the relationship between pressure head and tunnel radius.Figure 13
Relationship between seepage flow and permeability coefficient.


### 4.4. The Effect of Initial Support

When the inner diameterr1 of the grouting ring is constant, increasing r2 is to reduce the thickness of the initial support. Figures 14 and 15 show that as the thickness of the initial support decreases, H1 decreases and Q increases gradually. When r2=2.25m, Q and H1 are 13.3% of Q and 207% of H1 when r2=4.75m, respectively. It shows that the initial support thickness has a great influence on the pressure head. During construction, we should pay attention to the safety range of the initial support thickness.Figure 14
Relationship between pressure head and initial support radius.Figure 15
Relationship between seepage flow and initial support radius.If the initial diameterr2 and a are constant, increasing r1 means increasing the thickness of the initial support and reducing the thickness of the grouting ring. In Figures 16 and 17, with the increase of r1, H1 gradually increases and Q gradually decreases.Figure 16
Relationship between pressure head and radius of grouting ring.Figure 17
Relationship between seepage flow and radius of grouting ring.As shown in Figures18 and 19, the permeability coefficient of grouting ring is constant, and H1 gradually decreases and Q gradually increases with the increase of k2. The change rate of Q and H1 is larger when k2<10−8cm·s−1. When k2=2×10−8 cm·s−1, both Q and H1 tend to be stable. When the other conditions are constant, the change of k2 on the magnitude of 10−8 has little effect on Q and H1, but the influence on the magnitude of 10−9 is greater.Figure 18
Relationship between pressure head and permeability coefficient of grouting ring.Figure 19
Relationship between seepage flow and permeability coefficient of grouting ring.


## 4.1. Numerical Verification

Based on the complex function and groundwater hydraulics theory, the distribution of seepage field suitable for homogeneous and isotropic karst tunnels and lining structures were deduced. In order to ensure the accuracy of the analytical solution of the karst tunnel, the numerical solution simulated by the finite difference software FLAC3D is compared with the theoretical solution obtained by the method used in this paper.The calculation conditions are listed as follows: the head pressure of the cavern is 54 m; the radius of the cavern is 4 m; the radius of the tunnel is 7.25 m; the radius of the grouting circle is 2.25 m; and the initial support radius is 0.25 m; the permeability coefficient of surrounding rocks is 1.5 × 10−6 cm·s−1; the permeability coefficient of grouting ring is 10−7 cm·s−1; the initial support permeability coefficient of initial support is 10−8 cm·s−1.Six groups of 100 m × 100 m models of caverns on the right side of the tunnel were established and calculated with FLAC3D. The pore water pressure of the tunnel with gravity and without gravity was obtained, respectively. The six groups of models are d=15m, d=20m, d=25m, d=30m, d=35m, and d=40m. When d=20m, the karst tunnel model consists of 4236 elements, as shown in Figure 4. Figures 5 and 6 are contours of zone pore pressure with gravity and without gravity, respectively.Figure 4
Finite element model with circular karst tunnel.Figure 5
Contour of pore water pressure with gravity ford=20, a=7.25, and rw=4.Figure 6
Contour of pore water pressure without gravity ford=20, a=7.25, and rw=4.The pore water pressures extracted from the six groups of numerical simulation data in Figure5 with gravity and Figure 6 without gravity are compared with the theoretical solutions. The results are shown in Figure 7. The law obtained by numerical and theoretical solutions is that when the grouting ring and the initial support are completed and the secondary lining is not yet completed, H1 decreases with the increase of d. Moreover, the curves have a good consistency. When d=15m, the analytical solution for H1 is 25.78 m, the numerical solution with gravity is 25.67 m, and the numerical solution without gravity is 25.77 m. Therefore, the analytical method in this paper is suitable for solving the distribution of seepage field in karst tunnels.Figure 7
Comparison of numerical solution and theoretical solution fora=7.25 and rw=4.


## 4.2. The Effect of Distance and Cavern Radius

Figure8 shows the relationship between H1 and d. Figure 9 shows the relationship between Q and rw. In the infinite plane, when Hs is constant, the main factors affecting H1 and Q in the tunnel are rw and d. When rw is constant, H1 and Q decrease with the increase of d. When d is constant, H1 and Q gradually increase as rw increases.Figure 8
Relationship between head pressure and distance.Figure 9
Relationship between seepage flow and cavern radius.


## 4.3. The Effect of Grouting Ring

When the diameter of the grouting ring and the initial support radius is constant, the tunnel radius and the thickness of the grouting ring increase. Figures10–13 show that H1 and Q gradually decrease with the increase of the thickness of the grouting ring when the ratio between the distance and the permeability coefficient is constant. When the thickness of the grouting ring was 0.25 m and a reached 6 m, Q and H1 decreased to 41.94% and 41.86%, respectively. When a is constant, H1 and Q gradually decrease as the permeability coefficient decreases. When a increases to 7 m, the rate of H1 and Q decrease gradually. From Figures 12 and 13, when a is constant, H1 becomes smaller as the seepage coefficient of the grouting ring decreases. When kr=k/k1 increases to 80, Q decreases to 22.4% of Q when kr=5, and the Q reduction rate gradually decreases. In summary, increasing the thickness of the grouting ring and reducing the permeability coefficient of the grouting ring can reduce Q and H1. However, excessively increasing the thickness of the grouting ring will not only increase the cost but also reduce the impact on the pressure head and seepage flow. Therefore, a reasonable grouting material and construction process should be selected to ensure the thickness of the grouting ring when the permeability coefficient of the grouting ring is less than the specified value.Figure 10
Effect ofd on the relationship between pressure head and tunnel radius.Figure 11
Relationship between seepage flow and tunnel radius.Figure 12
Effect ofkr on the relationship between pressure head and tunnel radius.Figure 13
Relationship between seepage flow and permeability coefficient.


## 4.4. The Effect of Initial Support

When the inner diameterr1 of the grouting ring is constant, increasing r2 is to reduce the thickness of the initial support. Figures 14 and 15 show that as the thickness of the initial support decreases, H1 decreases and Q increases gradually. When r2=2.25m, Q and H1 are 13.3% of Q and 207% of H1 when r2=4.75m, respectively. It shows that the initial support thickness has a great influence on the pressure head. During construction, we should pay attention to the safety range of the initial support thickness.Figure 14
Relationship between pressure head and initial support radius.Figure 15
Relationship between seepage flow and initial support radius.If the initial diameterr2 and a are constant, increasing r1 means increasing the thickness of the initial support and reducing the thickness of the grouting ring. In Figures 16 and 17, with the increase of r1, H1 gradually increases and Q gradually decreases.Figure 16
Relationship between pressure head and radius of grouting ring.Figure 17
Relationship between seepage flow and radius of grouting ring.As shown in Figures18 and 19, the permeability coefficient of grouting ring is constant, and H1 gradually decreases and Q gradually increases with the increase of k2. The change rate of Q and H1 is larger when k2<10−8cm·s−1. When k2=2×10−8 cm·s−1, both Q and H1 tend to be stable. When the other conditions are constant, the change of k2 on the magnitude of 10−8 has little effect on Q and H1, but the influence on the magnitude of 10−9 is greater.Figure 18
Relationship between pressure head and permeability coefficient of grouting ring.Figure 19
Relationship between seepage flow and permeability coefficient of grouting ring.


## 5. Conclusion

Based on the inversion of complex function and groundwater hydraulics theory, this paper derives the analytical solution for the karst seepage field in the karst tunnel in an infinite plane. The potential function, the flow function, and the complex velocity at any place in the infinite plane are solved. According to the distribution of the seepage field and the boundary conditions, the seepage flow of the karst tunnel is obtained. The seepage field of the lining structure is deduced from the relationship between the grouting ring, the initial support, and the secondary lining. The new solution was verified with numerical simulation via the software FLAC3D. The conclusions obtained from this study are summarized as follows:(1)
WhenHs is constant for karst tunnels, rw and d are the main factors affecting the pressure head and seepage flow of the tunnel. When rw is constant, H1 and Q decrease with the increase of d. When d is constant, H1 and Q increase with the increase of rw.(2)
When the initial support is completed and the secondary lining has not been constructed yet,H1 and Q will gradually decrease with the increase of the thickness and the decrease of the permeability coefficient of the grouting ring. When the permeability coefficient of grouting ring is less than the specified value, reasonable grouting material and construction technology should be selected to ensure the thickness of grouting ring.(3)
H1 gradually decreases with the decrease of the thickness and the increase of k2. However, Q gradually increases with the decrease of the thickness and the increase of k2. During construction, we should pay attention to the safety range of the initial support thickness.



---

*Source: 9215472-2018-09-20.xml*