MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2019/1269064 1269064 Research Article Influence Zone Division and Risk Assessment of Underwater Tunnel Adjacent Constructions http://orcid.org/0000-0002-5776-4621 Zhou Zhong 1 http://orcid.org/0000-0003-2380-5859 Gao Wenyuan 1 Liu Zhuangzhuang 1 Zhang Chengcheng 1 Weber Gerhard-Wilhelm School of Civil Engineering Central South University Changsha 410075 China csu.edu.cn 2019 1512019 2019 25 08 2018 20 11 2018 24 12 2018 1512019 2019 Copyright © 2019 Zhong Zhou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Compared with ordinary tunnels, the influence analysis of underwater tunnel adjacent constructions is more complicated. At present, the empirical method used to divide the influence zone of the tunnel adjacent constructions has great uncertainty. So it is of great significance for actual construction and design to determine the influence zone accurately according to theoretical calculations. In this paper, based on the Hoek–Brown nonlinear failure criterion of rock mass and taking seepage factor into account, the stress state of rock mass around the underwater tunnel adjacent constructions can be deduced by elastoplastic theory. Then combined with the concept of “loose zone-bearing zone”, the influence zone division method of underwater tunnel adjacent constructions is proposed, and it is applied to the analysis of engineering examples. Through the deduced theoretical formulas, the influence zone of underwater tunnel adjacent construction can be divided into extensively strong, strong, relatively strong, weak, and noninfluence zones. Corresponding the influence zones with the risk levels in the code, different control measures are adopted for different risk levels, which can provide certain guidance for the design and construction of tunnel in practical engineering.

National Natural Science Foundation of China 50908234
1. Introduction

Presently, the development of urban subway traffic construction is rapid. As the utilization ratio of land resources and building space in the city is gradually improved, tunnel construction inevitably will be close to existing structures, which makes tunnel adjacent construction become the key link in current tunnel design and construction, and underwater tunnel is no exception. However, due to its more difficult construction, more complicated conditions, and more stringent requirements, the underwater tunnel adjacent construction is different from ordinary tunnels. How to evaluate the impact of underwater tunnel construction on adjacent structures is of great guiding significance for tunnel line selection and site construction.

Summarizing previous studies, it is found that the current research on the influence zone of tunnel adjacent construction is rarely about the underwater tunnel. Using the Japanese tunnel construction guidelines as examples, all of them are comprehensively evaluated on the basis of distance, engineering geology, and construction design. When applied to underwater tunnels, it is insufficient due to the lack of consideration of water weakening and permeation to surrounding rock. In addition, most scholars use the Mohr–Coulomb failure criterion for the theoretical analysis. Based on the development of constitutive models of rock mass in recent years, the Hoek–Brown failure criterion is considered as a more universal model. The Hoek–Brown nonlinear theory is extensively used in various kinds of problems due to its relatively accurate solutions. In this study, based on the Hoek-Brown yield criterion, considering the influence of water seepage, the elastoplastic analysis of the underwater tunnel construction is carried out, and the method of zoning of the underwater tunnel adjacent construction is further proposed. Corresponding the influence zones with the risk levels in the code, then the appropriate control measures are selected on the basis of the risk level to guide the site construction.

2. Elastoplastic Analysis of Surrounding Rock under Seepage Conditions 2.1. Hoek–Brown Failure Criterion and Basic Assumptions

Based on the Griffith brittle fracture theory, the Hoek–Brown failure criterion believed that the cause of rock failure is the deformation and expansion of the existing cracks, which has no significant relationship with the complete rock body yield strength. It is assumed that the rock mass cracks are irregular and the whole rock is isotropic. After years of research and development, the Hoek–Brown failure criterion has been amended several times. Finally, the relationship between rock parameters m, s, α and GSI is established, and its mathematical expression is as follows:(1)σ1=σ3+σcmbσ3σc+sα.

In the formula, σ1 and σ3 are the maximum and minimum principal stresses when rock mass is destroyed (the compressive stress is positive), σc is the uniaxial compressive strength of complete rock mass, and m, s and α are the dimensionless parameter, which have the following equivalent expressions:(2)mb=mi·expGSI-10028-14D.(3)s=expGSI-1009-3D.(4)α=12+e-GSI/15-e-20/36.

In the formula, D is the disturbance factor that reflects the type of rock mass, and GSI reflects the integrity of rock mass.

Based on the revised Hoek–Brown failure criterion, the elastic–plastic stress state of underwater tunnel excavation is analyzed, and the calculation model is shown in Figure 1. In the model, the tunnel section shape is circular, the radius of the tunnel is a, the inner water head is ha, the surrounding rock plastic zone radius is Rp, the water head at the edge of the plastic zone is hp, the normal stress at the interface between the elastic zone and the plastic zone is Pp, and the elastic zone radius is Re, the water head far enough is h0, the original rock stress of the surrounding rock is P0, and the lining support force is Pa. In order to simplify the calculation model and the solution process, the following basic assumptions are made during the solution of the circular tunnel: (1) assuming that the surrounding rocks around the tunnel are homogeneous rocks, and ignoring the self–weight of the calculation unit, the lateral pressure coefficient is 1.0; (2) with the same permeability coefficient of rocks, the seepage direction is radial to form a stable seepage field; (3) the tunnel lining support force is evenly distributed along the radial direction; (4) the length of the tunnel is long enough to handle the problem as axially symmetric plane strains during calculation; (5) finally, the elastoplastic analysis of the surrounding rock follows the Hoek–Brown failure criterion.

Calculation model of surrounding rock.

2.2. Seepage Field Status Calculation

According to Darcy’s law and the continuity equation of seepage flow in the underwater tunnel, the seepage continuous equilibrium differential equation (5) can be obtained. Equations (6) and (7) are the seepage flow boundary conditions which are expressed as follows:(5)2Hr2+1rHr=0,(6)Hrr=a=ha,(7)Hrr=Re=h0.

Combined with above equations, the water head H (r) can be obtained:(8)Hr=1lnRe/ah0lnra+halnRer.

2.3. Elastic–Plastic Zone Stress State Calculation

When the stress state of the elastic–plastic zone of the surrounding rock is calculated, the seepage water pressure acts on the unit in the form of volume force, and the expression is as follows:(9)fr=-γwdξHdr=γwξha-h0rlnk.

Considering the effect of seepage volume force, according to elastoplastic calculation principle, the balance differential equation of elastic zone unit is expressed as follows:(10)dσrdr+σr-σθr+γwξha-h0rlnk=0.

In the formula, σr is radial stress; σθ is tangential stress (tensile stress is positive and compressive stress is negative). γw is the weight of water, ξ is the equivalent pore water pressure coefficient, and k=Re/a.

After the unit is deformed, according to the displacement–strain geometric equation (11) and the physical stress–strain physical equation (12), then combined with the equilibrium differential equation, (13) can be obtained:(11)εθ=ur,εr=dudr,γrθ=0,(12)σr=1-μE1+μ1-2μεr+μ1+μεθσθ=1-μE1+μ1-2μεθ+μ1+μεr.(13)u=C3r+C4r+1+μ1-2μ21-μE×ξγwh0-halnkrlnr.

After introducing the elastic boundary condition equation (14), the expression of stress state in the elastic zone of surrounding rock considering the seepage condition can be solved as (15) and (16):(14)σrr=Rp=Pp,σrr=Rε=P0.(15)σr=P0+ξγwh0-halnr/Re21-μln1/k+Rp2Rp2-Re2×Re2r2-1P0-Pp+ξγwh0-halnRp/Re21-μln1/k.(16)σθ=P0+ξγwh0-halnr/Re+2μ-121-μln1/k-Rp2Rp2-Re2×Re2r2+1×P0-Pp+ξγwh0-halnRp/Re21-μln1/k.

When the stress state of the plastic zone is calculated, the surrounding rock follows the Hoek–Brown failure criterion in the plastic zone. According to the force analysis when λ=1, it can be known the tangential pressure σθ is the maximum principal stress σ1 and the radial stress σr is the minimum principal stress σ3, the Hoek–Brown failure criterion is expressed as follows:(17)σθ=σr+σcmbσrσc+sα.

The equilibrium equation of the unit can be changed into (18) as follows:(18)dσrdr-σcmbσr/σc+sαr+γwξh0-harlnk=0.

After integrating (18), (19) and (20) can be obtained. According to (20), the amount of change of the radial stress along the radial direction in the plastic zone can be determined. Furthermore, according to (17), the amount of change of the tangential stress along the radial direction in the plastic zone can be determined:(19)Paσrdσrσcmbσr/σc+sα-γwξh0-ha/lnk=ardrr=lnra.(20)r=a×expPaσrdσrσcmbσr/σc+sα-γwξh0-ha/lnk.

2.4. Calculation of the Radius of Loose Zone–Bearing Zone

According to the definition of the “loose area”, it is assumed that the tangential stress on the boundary of the loose area is the original rock stress, that is σθ=σp, then it can be obtained according to the Hoek-Brown failure criterion:(21)σθ=σl+σcmbσlσc+sα=P0.

Substituting (21) into (20), the radius of the loose zone can be obtained as follows:(22)RL=a·expPaσldσrP0-σr-γwξh0-ha/lnk=aP0-Pa-γwξh0-ha/lnkP0-σl-γwξh0-ha/lnk.

It can be seen from (22) that the mechanical parameters of rock mass, groundwater pressure, and seepage effect have a great influence on the radius of the loose zone. Meanwhile, the larger the section size of the excavation tunnel, the larger the radius, whereas the larger the tunnel support force Pa, the smaller the radius.

The radius of the bearing zone can be determined by the following:(23)RpRb0σθe-P0dr=12RpReσθe-P0dr.

By substituting (12) into (23), the radius of the bearing zone Rb0 can be obtained.(24)RpRb0ξγwh0-halnr/Re+2μ-121-μln1/k+Rp2Rp2-Re2×Re2r2+1×P0-Pp+ξγwh0-halnRp/Re21-μln1/k=12RpReξγwh0-halnr/Re+2μ-121-μln1/k+Rp2Rp2-Re2×Re2r2+1×P0-Pp+ξγwh0-halnRp/Re21-μln1/k.

3. Underwater Tunnel Adjacent Construction Influence Zones Division

After the excavation of the underwater tunnel, the initial stress equilibrium state of the surrounding rock is damaged and the stress is redistributed. The rock stress around the tunnel exceeds the yield stress of the rock mass and the rock becomes damaged or enters the plastic state. Subsequently, the stress of the plastic zone around the tunnel is partly transferred to the deep rock mass due to the stress redistribution, while the other part is relieved and eliminated due to the deformation. As the distance from the cave wall increases, the minimum radial principal stress increases and the bearing capacity of the rock mass also increases, which makes the stress state of the surrounding rock transition from plastic state to elastic state in space.

Compared with the initial state of stress, the plastic zone of surrounding rock can be divided into the bearing zone and the loose zone after excavation. The area in the deep plastic zone where the stress is higher than the initial stress and the area with the higher stress in the elastic area of the surrounding rock are combined as the “bearing zone”, while the area whose stress is lower than the initial stress is the “loose zone”. In the loose zone, the stress and crack propagation of the surrounding rock increase, and the plastic slip is obvious. When the surrounding rock reaches the plastic state, all mechanical parameters are deteriorated. At the same time, the “weak” surrounding rock structure is further weakened by the action of groundwater through pore hydrostatic pressure and hydrodynamic pressure, which makes it easier to damage.

According to the above analysis of the stress state of the surrounding rock after the underwater tunnel excavation and the stress graph of the surrounding rock, the influence zone of the adjacent construction is divided. The loose and bearing zones are considered to be the strong influence zone. In the elastic zone, under the elastic stress state, the range of tangential stress at the hole edge greater than or equal to 1.01 times of the initial stress is defined as the elastic stress concentration zone. Therefore, the stress concentration zone in the elastic zone can be regarded as the weak influence zone and the outer zone is an initial stress zone. The influence zones division is shown in Figure 2.

Elastic–plastic stress status of the surrounding rock in underwater tunnel. In the figure, 1 and 2 are the plastic, and 3 and 4 are the elastic zones; 1 is the loose zone, 2 and 3 are the bearing zones, and 4 is the elastic stress concentration zone.

Based on the Hoek-Brown failure criterion, the loose zone can be divided into three regions, the strong influence, weak influence, and no influence zones, and the strong influence zone can be divided into extensively strong influence, strong influence, and fairly strong influence zones.

4. Engineering Example

The research on the influence zone division of the adjacent construction has a strong guiding role in the design and construction of underwater tunnels. For example, when the adjacent construction is in a strong influence zone, reinforcement and support are required. When it is in a weak influence zone, monitoring measures are needed. When it is in a no influence zone, no measures are needed. Figure 3 is an example of underwater shield tunnel engineering.

Location of the engineering example.

The underwater tunnel adopts Φ6250 earth pressure balance shield machine construction. Under the river, the left tunnel is parallel to the existing tunnel main line; the minimum spacing between the outer contour of the tunnel and the outer contour of the existing tunnel is 14.700m. Then below the river bank, new tunnel underpasses the existing tunnel; the vertical distance from the outer contour of the new tunnel to the main floor of the existing tunnel is 5.442 m. The riverbed strata are mainly highly weathered conglomerate, and moderately weathered conglomerate exists locally. Two typical sections are selected for calculation and analysis. According to geological data, the water depth of the river in section A is 6 m, the depth of the shield tunnel is 10.906 m, and the whole tunnel is in strong weathered conglomerate. The water depth of river in section B is 1 m, the depth of the shield tunnel is 26.73 m, and the tunnel is also in strong weathered conglomerate.

The relevant parameters at section A are as follows: the tunnel radius is a = 3.14 m, the inner water head is ha = 0 m, the water head far enough is h0 = 20 m, the original rock stress of surrounding rock is P0 = 0.4 MPa, the equivalent pore water pressure coefficient is ξ = 1.0, the surrounding rock elastic modulus is E = 45 MPa, Poisson’s ratio is μ = 0.3, the cohesion is c = 50 kPa, and the internal friction angle is φ = 37°.

By substituting r = Rp in Equation (19), the radial stress Pp of the surrounding rock at the elastic–plastic zone can be obtained. Then the plastic zone radius Re can be obtained. Combined with the relevant parameters of surrounding rock and tunnel, the radius of the plastic zone after the excavation of tunnel is 5.76m.

The radii of the loose zone and the bearing zone are further solved according to (22). When substituting the surrounding rock and tunnel parameters into the equation, the following can be obtained: the loose and bearing zones radii are 4.58 m and 9.53 m. Therefore, the range of extensively strong influence zone is 3–4.58 m, the range of strong influence zone is 4.58–5.76 m, and the range of fairly strong influence zone is 5.76–9.53 m. After solving the radius of the loose and bearing zones, according to (15) and (16), then combining the calculation result of the plastic zone, the following can be obtained: the radius of stress concentration in the elastic zone is 18.44 m, so the radius of the weak influence zone is 18.44 m. The influence zones of section A are shown in Figure 4.

Influence zones of section A.

The relevant parameters at section B are as follows: the tunnel radius is a = 3.14 m, the inner water head is ha = 0 m, the water head far enough is h0 = 10 m, the original rock stress of surrounding rock is P0 = 0.5 MPa, and the equivalent pore water pressure coefficient is ξ = 1.0; and the surrounding rock elastic modulus is E = 50 MPa, the Poisson’s ratio is μ = 0.2, the cohesion is c = 55 kPa, and the internal friction angle is φ = 40°.

Based on the process of calculating the radius of the influence zones in section A, the radius of the influence zones in section B can be calculated as follows.

It can be obtained from (19) that the radius of the tunnel excavation plasticity influence zone is 6.17m, then the radius of the loose zone is 4.72m and the radius of the bearing zone is 8.28m. Therefore, the range of extensively strong influence zone is 3–4.72 m, the range of strong influence zone is 4.72–6.17 m, and the range of fairly strong influence zone is 6.17–8.28 m. Similarly, the radius of the stress concentration zone in the elastic zone is 15.22m, so the radius of the weak influence zone is 15.22m.The influence zones of section B are shown in Figure 5.

Influence zones of section B.

5. Adjacent Construction Influence Zones and Risk Classification

On the basis of the adjacent construction influence zone division, different influence zones of the adjacent construction are associated with different risk levels in the risk assessment, so that the risk level can be determined according to different influence zones during construction. Then the corresponding measures can be taken. Different countermeasures are adopted for the adjacent construction with different risk levels, which has guiding significance for the actual project.

Referring to GB50652–2011 underground rail transit urban construction risk management practices , GB / T 50839–2013 urban rail transit safety control technical specifications , and CJJ / T 202–2013 urban rail transit structural safety protection technical specifications , the influence zones obtained above and risk levels are corresponded. Different risk levels should be taken different risk disposal guidelines and control plans, as shown in Table 1.

Adjacent construction influence zone and risk classification.

No. Adjacent construction influence zone Risk level Accepted guidelines Disposal guidelines Control plans
1 Extensively influence zone Level I Unacceptable Risk control measures must be implemented to reduce the risk, at least to a level that is either acceptable or unwilling Risk warning and emergency response programs, or program amendments or adjustments should be prepared

2 Strong influence zone Level II Unwilling to accept Risk management should be implemented to reduce risk, and the cost of risk reduction should not be higher than the risk after the risk of loss Risk prevention and monitoring should be implemented and risk management measures should be formulated

3 Fairly Strong influence zone Level III Acceptable Risk management should be implemented and risk management measures can be taken The daily management and monitoring frequency should be increased

4 Weak influence zone Level IV Ignorable Can implement risk management Examination can be conducted daily

5 No influence zone None No No need to consider the risk General management should be checked
6. Control Measures

For the adjacent construction, three aspects can be controlled: one is to reinforce the existing structure, the other is to strengthen the new structure, and the third is to take measures against the soil between the existing and the new structures. Therefore, the following measures can be taken for the adjacent construction.

( 1 ) Measures to be taken for the existing tunnels or projects: the types of countermeasures include basic measures, strengthening measures, and repairing measures. The basic measures include backfilling grouting and prevent tunnel lining off the block. Strengthening measures include the use of arches, anchorage, and crossbar to strengthen tunnel linings. Repairing measures include stripping of floating blocks that may fall off, cleaning the surface of the tunnel lining, repairing drainage ditch, and preventing water leakage.

( 2 ) Measures to be taken for the new tunnels or projects: when construction has to be conducted in a strong-influence zone, in order to reduce the impact on the existing tunnel, a more perfect construction scheme or new countermeasures should be developed. The concrete measures include adjusting the driving parameters and changing the supporting structure of lining.

( 3 ) Measures to be taken for the surrounding rock between the existing and the new construction or projects: when the adjacent construction has a bad influence, the protection of the existing tunnel and the measures taken to the side of the new construction are insufficient; measures should be taken to deal with the intermediate strata in order to reduce or eliminate the impact. Generally, methods such as grouting method and freezing method are adopted to strengthen and improve the stratum, and methods of isolating influence may also be adopted, such as underground continuous walls, pipe sheds, steel pipe piles, and arch protection.

Under normal circumstances, the above measures are not completely independently used, and comprehensive application can achieve better control effect. In addition, construction monitoring is an important part of underground engineering. Due to the high risk, complexity, and unpredictability of adjacent construction, monitoring is even more significant. Monitoring data must be fed back in time to guide follow-up construction.

Based on the influence degree of the adjacent construction, the risk level is judged, and then the actual situation of the site is combined to take corresponding measures. For the grade I risk level, the existing structure, the new structure, and the surrounding rock between the existing structure and the new structure should be taken reinforcement measures. For the grade II risk level, the new structure and the surrounding rock between the existing structure and the new structure should take reinforcement measures, and the monitoring of existing structures is needed. For the grade III risk level, the surrounding rock between the existing structure and the new structure should take reinforcement measures, and the monitoring of the existing and the new structures are needed. For the grade IV risk level, the monitoring of the existing and the new structures are needed so as to grasp the impact of construction in real time.

7. Conclusions

In this paper, based on the Hoek–Brown failure criterion, the theoretical analysis is conducted on the underwater tunnel adjacent construction. In consideration of seepage flow, the influence zone of underwater tunnel adjacent construction is divided, then influence zones are matched with risk levels in the code. The following conclusions can be obtained.

( 1 ) The elastic–plastic stress state of the surrounding rock of underwater tunnel is different from that of ordinary tunnel. The seepage action can cause a certain degree of deterioration to the surrounding rock, which cannot be ignored.

( 2 ) According to the proposed zoning method for the influence zone of underwater tunnel adjacent construction in this study, the influence zone can be divided into extensively strong, strong, fairly strong, weak, and no-impact zones. During the construction, measures will be taken on the basis of the different influence zone.

( 3 ) The influence zones are matched with the risk levels, and different countermeasures are taken according to different risk levels, which have guiding significance for the site construction.

Symbols σ 1 :

The maximum principal stress, MPa

σ 3 :

The minimum principal stress, MPa

σ c :

The uniaxial compressive strength of complete rock mass, MPa

m :

The dimensionless parameter

S :

The dimensionless parameter

α :

The dimensionless parameter

D:

The disturbance factor that reflects the type of rock mass

GSI:

The integrity of rock mass

a :

The radius of the tunnel, m

h a :

R p :

The surrounding rock plastic zone radius, m

h p :

P p :

The elastic–plastic contact interface at the normal stress, N

R e :

R L :

h 0 :

The water head far enough, m

P 0 :

The original rock stress of the surrounding rock, N

P a :

The lining support force, N

σ r :

σ θ :

The tangential stress, MPa

γ w :

The unit weight of water

ξ :

The equivalent pore water pressure coefficient.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.