Abstract

This work aims to characterize the in situ stress field along the Lijiang to Shangri-La railway and identify possible engineering geological problems when constructing tunnels along this railway on the margin of the Tibetan Plateau. The in situ stress measured at 76 points in 12 boreholes by the hydraulic fracturing method was analysed. A rose diagram of the maximum principal stress direction was plotted based on the measured in situ stress data. The results show that the measured in situ stress is mainly horizontal stress, corresponding to a strike-slip fault-related tectonic stress field with a moderate to high in situ stress value. The main stress values have a good linear relationship with the burial depth, and the maximum horizontal principal stress (σH) increases by 1.1–8.8 MPa per 100 m, with an average gradient value of 3.6 MPa per 100 m. The maximum and minimum horizontal principal stresses and the stress differences increase with depth, and the lateral pressure coefficient (σH/) is generally 1–1.5. The ratio of the maximum and minimum effective stresses is less than the threshold at which faulting occurs, resulting in faults that are relatively stable at present. The direction of the maximum horizontal principal stress is oriented at a small angle to the axial direction of the deeply buried tunnel along the railway line; therefore, the tunnel sidewalls could readily deform during the construction process. Rock bursts are expected to occur in strong rock masses with high risk grades, whereas slight to moderate deformation of the rock surrounding the tunnel is expected to occur in weak rock masses. This study has significance for understanding the regional fault activity and issues related to the construction of deeply buried tunnels along the Lijiang to Shangri-La railway.

1. Introduction

In situ stress is a fundamental parameter in a wide range of endeavours in rock mechanics and engineering construction, and it is related to many engineering problems, such as rock bursts, rock slack, and rock deformation [14]. Since A. Heim (1912) proposed the concept of in situ stress in an Alpine tunnel project [5], the in situ stress theory has been widely applied in underground tunnel projects. In general, the in situ stress can be measured in the field or simulated by finite element models or other model codes. A series of in situ stress measurement techniques have been developed to interpret stress in different geological conditions at a given point [1, 2]. In situ stress measurement methods include the hydraulic fracturing method [2, 69], stress relief method [710], flat jacking method [9], borehole breakout method [9], drilling-induced tensile fracture method [11], acoustic emission method [12], strain recovery method [9], differential strain curve analysis method [9, 13], and geophysical method [9]. Among them, the hydraulic fracturing method and the stress relief method with overcoring are usually employed for in situ stress measurement in projects with boreholes [610]. Additionally, numerical simulations of the in situ stress field have become an efficient approach to study the stress distribution on a regional scale [1416]. These studies can be used to understand the different measurement approaches and analyse the in situ stress features from a 2- or 3-dimensional simulation perspective. As China’s western development strategy is further implemented, the number of deep and large road and railway tunnel projects will increase in the future. Unlike shallow engineering projects, deeply buried tunnel projects feature high sidewalls, large spans, high in situ stresses, and complex geological conditions [1216]. The deformation or damage of the rocks surrounding deeply buried tunnels is distinctly influenced by the in situ stress field. Here, in situ stress not only determines the stability of surrounding rock masses but also acts as a load that can cause deformation and damage in tunnel projects. Therefore, determining the regional in situ stress characteristics is a prerequisite for analysing the stability of a surrounding rock mass, realizing the excavation design of underground projects, and conducting scientific decision-making.

As the region with the most intense tectonic activity on the margin of the Tibetan Plateau, the region between Lijiang and Shangri-La is characterized by large-scale active faults with a dense distribution. These long faults exhibit frequent, large-magnitude seismic activity and large displacements, which leads to a very complex in situ stress field. The railway from Lijiang to Shangri-La features numerous tunnels, and the stability of the surrounding rocks is closely related to the stress state of the rock masses (Figure 1). Many scholars have explored the region’s tectonic stress field and other crustal dynamics through the analysis of focal mechanism solutions, the fault fracture modes, the characteristics of crustal deformation, inversion of the drainage pattern, stress measurements, and physical model tests in the Lijiang and Shangri-La regions [1924]. However, the in situ stress cannot be effectively verified from the previous literature because of the limited number of measurement points. On this basis, a systematic in situ stress test using the hydraulic fracturing method was applied to boreholes that were drilled to survey the site of a proposed deeply buried tunnel project along the Lijiang to Shangri-La railway. The in situ stress was characterized from the geomechanics perspective, and the effects of fault activity and construction of a deeply buried tunnel on rock bursts and large deformations were also evaluated. This work not only provides insight into the in situ stress characteristics on the margin of the Tibetan Plateau but is also useful for tunnel axial arrangement and lining design for the Lijiang to Shangri-La railway.

2. Geological Setting

The region between Lijiang and Shangri-La in Yunnan, China, is located in the rhombus-shaped Sichuan-Yunnan fault block on the southeastern margin of the Tibetan Plateau (Figure 1(a)). This area is part of the middle section of the Hengduan Mountains, representing the frontal part of the suture zone between the Indian Plate and the Eurasian Plate. The study area is adjacent to the active faults along the block boundary to the southwest and is within the Zhongdian-Dali seismic belt in the southern part of the North-South seismic belt. The study area is experiencing uplift and tectonic movement to the southeast, driven by the indention of the Indian Plate and extrusion of the Tibetan Plateau (Figure 1(a)). Therefore, the structural pattern of the study area has been strongly transformed and deformed by neotectonic movement [2527], and the study area exhibits a NNW-SSE stress direction according to the World Stress Map (Figure 1(b)) [18, 28]. There are many deep and large active faults within the region, including the Zhongdian-Hailuo fault, the Xiaozhongdian-Daju fault, the Xiaojin River-Lijiang fault, the Zhongdian-Longpan-Qiaohou fault, and the Daju-Lijiang fault. These faults have accommodated intense active tectonic movement around the Tibetan Plateau in China. The railway from Lijiang to Shangri-La is oriented in the NNW direction and crosses the Xiaojin River-Lijiang fault, the Daju-Lijiang fault, and the Zhongdian-Longpan-Qiaohou fault, as shown in Figure 1(c). In the Quaternary period, the regional stress field was oriented in the NNE direction during the early Pleistocene period, shifted to the E-W direction in the late early Pleistocene to the early late Pleistocene period, and finally shifted to the NNW direction in the middle period of the late Pleistocene. The present crustal stress direction continues to be oriented in the NNW direction according to the World Stress Map (Figure 1(b)) [18, 26]. In particular, the stress fields with NNE and NNW directions were associated with the development phase of a fault basin in northwestern Yunnan, while the stress field with the E-W direction coincided with the highly differentiated uplift and subsidence phase [29].

3. In Situ Stress Measurement Method

The in situ stress measurement was conducted based on the boreholes drilled as part of the deeply buried tunnel survey. Considering the feasibility of measuring the in situ stress in deep holes and the working conditions in situ, the hydraulic fracturing technique was chosen as the main method to measure the in situ stress along the Lijiang to Shangri-La railway. However, there are still some limitations because the surface measurements should be done in vertical or subvertical boreholes. Here, 12 deep boreholes with 76 in situ measurement points in the proposed deeply buried tunnel were selected for the in situ stress test based on the lithological characteristics, geomorphological features, and fault distribution along the Lijiang to Shangri-La railway. The hydraulic fracturing method uses high-pressure pumps to generate high water pressure that either forms new fractures or reopens preexisting fractures.

A schematic diagram of the hydraulic fracturing process and a field photograph are shown in Figure 2. The hydraulic fracturing equipment was mainly composed of two 120 cm long sealing packers, an impression packer, a high-pressure pump, a chart recorder, and a flow meter (Figure 2). While the flow of high-pressure fluid between the oil pump and the sealing packers occurred via a pressure-resistant hose, the flow between the water pump and the isolated test section of the borehole occurred via drilling rods, which enabled the packers to be hoisted easily. This generates tensile stresses in the borehole wall. Pressurisation continues until the borehole wall ruptures through tensile failure and hydrofracturing is initiated. The hydraulic fracturing pressure and calculation of the in situ stress are described in detail by Haimson and Cornet [6]. The hydraulic fracturing method allows a direct measurement of the minimum stress in the plane perpendicular to the borehole axis, which is normally the minimum horizontal principal stress (σh), and the accuracy is good (<∼ ±3%). The maximum horizontal stress is calculated from equations that include failure criteria and parameters evaluated from the field pressure data. The accuracy is lower for the maximum horizontal principal stress (σH), with an error of ∼ ±3–8%. The orientation of the maximum horizontal principal stress is accurate to less than ±5°, corresponding to the A-quality class as defined as within ±15° by Heidbach et al. [18]. In the Lijiang to Shangri-La railway, 76 data points were obtained from the studied boreholes, as shown in Table 1.

4. Results of In Situ Stress and Faults

4.1. Type of In Situ Stress Field

Among the 76 in situ measured points in the 12 drilling boreholes, 64 data points have a maximum horizontal principal stress (σH) larger than the vertical principal stress (σv), accounting for 84.2% of the total measured points. The other 12 data points of maximum horizontal principal stress are less than the vertical principal stress, accounting for 15.8% of the total measured points. Therefore, the horizontal in situ stress is the dominant stress in the rocks along the Lijiang to Shangri-La railway, and this stress predominantly results from tectonic stress. This in situ stress field represents a typical tectonic stress field. Among the 76 in situ measured points, there are 10 measured points with σH > σh > sv or σv > σH > σh, accounting for 13.2% of the total measured points. However, there are still 56 points with σH > σv > σh, accounting for 73.6% of the total measured points and showing that the in situ stress field is conducive to the formation of strike-slip faults. The western boundary of the Zhongdian-Longpan-Qiaohou fault (No. 23 in Figure 1(c)) presents dextral strike-slip offset features [30], and the horizontal displacement is larger than the vertical displacement, as determined with the Global Positioning System (GPS) technique. This GPS observation result is consistent with the results of the in situ stress measurement. From the geodynamics perspective, the horizontal stress of the in situ stress field dominates the tectonic stress field and is mainly produced by the eastward movement of the Tibetan block due to the collision of the Indian Plate and the Eurasian Plate (Figure 1(a)) [17, 31]. The rhombus-shaped block of the western Zhongdian-Longpan-Qiaohou fault rotates clockwise under the horizontal stress oriented in the SE and SSE directions. Consequently, this rhombus-shaped block is characterized by greater displacement in the south and west than in the north and east.

4.2. Characteristics of the Horizontal Principal Stress Orientation

The rose chart of the horizontal principal stress orientation and focal mechanism solution is illustrated in Figure 3. According to the statistical analysis, N10°–60°W is the dominant orientation of the maximum horizontal principal stress (σH), which results in a predominant orientation of N15°W (Figure 3(a)). Moreover, the principal stress orientation and predominant orientation of the modern tectonic stress field are N35°–5°W and N15°W, respectively, according to the analysis of focal mechanism solutions in the northwestern Yunnan region (Figure 3(b)) [32]. However, the World Stress Map of the region (Figure 1(b)) suggests that the dominant orientation is N-S with an error of ±15°–25°. Compared with that of the World Stress Map, the horizontal principal stress orientation from on-site measurements is accurate to less than ±5°, indicating that on-site measurements more accurately reflect the in situ stress field along the Lijiang to Shangri-La railway. From the structural trace perspective, fault scars suggest that the regional principal stress field has been oriented to the NNW from the middle period of the late Pleistocene to the present [33]. The result from structural traces is generally consistent with the horizontal stress orientation from the on-site measurement. However, some on-site points reflect local stress characteristics, which are controlled by the local landforms, lithology, crustal structure, etc. In particular, borehole No. 1 exhibited a local stress field related to its position near a deep gorge, and the principal stress was rotated towards the river valley, with the maximum horizontal stress oriented in the direction of N42°–53°E. Meanwhile, the principal stress of borehole No. 9, which is located near the Chongjiang River fault (striking N32°W), is affected by the fault and is deflected in the direction of N50°–60°W.

4.3. In Situ Stress Characteristics Relationship with Depth
4.3.1. Variation Pattern of the Principal Stress with Depth

Statistical analysis of the data in Table 1 shows that the maximum horizontal principal stress is 3.94–25.09 MPa, and the stress gradient is 1.1–8.8 MPa/100 m, with an average of 3.6 MPa/100 m. The minimum horizontal principal stress is 2.43–18.2 MPa, and the stress gradient is 1.3–5.1 MPa/100 m, with an average of 2.6 MPa/100 m. A scatter diagram of the in situ stress relationship with burial depth is plotted in Figure 4. The relationship between the horizontal principal stress and the burial depth was determined by the least square linear regression method. The maximum horizontal principal stress is σH = 1.853 + 0.0266H, with a correlation coefficient R2 = 0.702, and the minimum horizontal principal stress is σh = 1.043 + 0.0182H, with an R2 = 0.627. As shown in Figure 4, the maximum and minimum horizontal principal stresses basically increase with burial depth and exhibit good linear relationships. However, the horizontal principal stresses still exhibit discrete phenomena when the burial depth is less than 600 m, with the largest difference of 5–6 MPa at a burial depth of approximately 380 m.

4.3.2. Variation Pattern of the Lateral Pressure Coefficient with Depth

The lateral pressure coefficient (LPC, σH/) is the ratio of the maximum horizontal principal stress and the vertical stress, which is illustrated in Figure 5. The LPC generally ranges from 0.75 to 1.91, with an average value of 1.21. The LPC is higher than 1 at 64 on-site measured points, accounting for 84% of the total measured points. The LPC ranged between 1.0 and 1.5 at 54 of the 64 points, accounting for 71% of the 64 points. However, the LPC at the other 10 points ranged from 1.5 to 2.0, accounting for 13%. This result indicates that regional stress is dominated by tectonic stress. Here, the in situ stress in our study area can reach a moderate to high level, which is based on the criterion of the in situ stress field from Xue et al. [34]. The LPC shows great discreteness when the burial depth is shallower than 400 m. However, the LPC gradually approaches 1 as the burial depth increases (Figure 5), indicating that the horizontal stress or tectonic stress becomes equal to the gravitational stress.

4.3.3. Variation Pattern of the Horizontal Stress Difference with Depth

The regional horizontal shear stress can be indirectly reflected by the difference between the maximum and minimum horizontal principal stresses [18, 33]. The scatter diagram of horizontal stress differences (ΔσHh) with burial depth is illustrated in Figure 6. The relationship between ΔσHh and burial depth (H) is ΔσHh = 0.0084H + 0.8103, R2 = 0.6301, based on linear least square regression. The scatter diagram shows that ΔσH–h increases with burial depth, which results in the gradual increase in ΔσHh from 1 MPa at approximately 200 m depth to 7.4 MPa at approximately 850 m (Figure 6).

4.4. Fault Characteristics along the Lijiang to Shangri-La Railway

The Lijiang to Shangri-La railway crosses the Xiaojin River-Lijiang fault, the Daju-Lijiang fault, and the Zhongdian-Longpan-Qiaohou fault from south to north, as illustrated in Figure 1(c) and Table 2. The Xiaojin River-Lijiang fault zone does not show obvious signs of Quaternary deformation and is mainly associated with erosional linear depressions and drainage systems eroding the preexisting thrust belt. The Xiaojin River-Lijiang fault was an active thrust fault with a left-lateral component in the Neogene–early Pleistocene (NQ1) period, with a maximum historical earthquake magnitude of approximately M 6.1 (Table 2). The Lijiang-Daju fault controls the fault depressions in the subsided Lijiang-Daju basin with a maximum depth of >1200 m and was active in the late Quaternary (Q3Q4) period [23, 24]. The Lijiang-Daju fault is a normal fault with a left-lateral component and a maximum earthquake magnitude of M 7.0 in 1996. The Lijiang-Daju fault is characterized by fault triangles, fault scarps, and escarpments composed of dislocated late Quaternary strata, and broken strata with traction deformation are observed. The Longpan-Qiaohou fault shows obvious Quaternary activity and great deformation in the fault zone, resulting in an inconsistent distribution of strata on either side and large differences in occurrence. The Longpan-Qiaohou fault is a thrust fault with a right-lateral component and a historical maximum earthquake magnitude of M 6¾ in 1751 (Table 2). Along the Longpan-Qiaohou fault, Quaternary fault basins, such as the Jiuhe, Jianchuan, and Shaxi basins, have developed, and the maximum deposition thickness has reached >600 m in the Jianchuan Basin. The Zhongdian fault is located on the northern side of the railway and was active in the Quaternary, specifically during the late Pleistocene to Holocene, and it is responsible for the formation of the Zhongdian and Xiaozhongdian depression basins [23, 24]. The Zhongdian fault is a right-lateral strike-slip fault with a maximum earthquake magnitude of 6¼ in 1933.

5. Analysis and Discussion

5.1. Activity of the Fault

Fault activities mainly include stick-slip movement and creep movement. Creep movement gradually releases crustal energy, which causes many consecutive earthquake events and crustal deformation. However, stick-slip movement is associated with abrupt changes and results in displacement, crustal rupture, and surface deformation and is always associated with an earthquake event with a number of aftershocks after a long quiescent period [1922]. According to the frictional sliding criterion of faults, the relationship between the crustal stress state and fault activity follows Coulomb’s law of friction. Assuming that the cohesion (intrinsic strength) of a fault is 0 and considering the concepts of effective stress, average stress, and maximum shear stress, the formula of seismic activity can be obtained as follows [35]:

In the formula, σ1 and σ3 are the maximum and minimum principal stresses on the periphery of the fracture. For strike-slip fault activity, σ1 and σ3 are equal to the maximum horizontal stress σH and the minimum principal stress σh, respectively, but for reverse fault activity, σ1 and σ3 are the maximum horizontal stress σH and the vertical principal stress σv, respectively. U0 is the pore pressure, and μ is the friction coefficient of the fault. The friction coefficient (μ) was equal to 0.5–0.7 if the pore water pressure was ignored [3638]; however, μ should be 0.6–1.0 to evaluate the shallow fault activity [37, 38].

If the ratio of the maximum and minimum effective principal stresses (σ1/σ3) is less than or equal to µm, from equation (1), the faults are stable. However, if this ratio (σ1/σ3) is larger than µm, faulting will occur. If the pore water pressure is ignored and μ = 0.5 is the critical friction coefficient of a fault instability, a value of μm = 2.6 can be inferred. After applying the in situ stress data from Table 1 in equation (1), the K value is in the range of 1.2–1.9, which is less than 2.6, indicating that the ratio of the maximum and minimum effective stresses (σ1/σ3) in the study area is less than the critical value at which faulting will occur and that the faults are relatively stable at present. This result may be affected by the Ms 7.0 Lijiang earthquake and its aftershocks in 1996, which resulted in the release of crustal energy to a certain extent and, consequently, a reduction in the stress level to some extent.

5.2. Stability of a Rock Mass in Tunnel Cavern

The principal stress along the Lijiang to Shangri-La railway is dominated by intense compression in the horizontal direction. The Lijiang to Shangri-La railway line is oriented in the north-south direction at a small angle to the maximum horizontal principal stress direction. Additionally, the shear failure of a rock mass mainly depends on the principal stress difference. The larger the burial depth of the underground excavation is, the larger the difference between the maximum and minimum principal stresses and the lower the stability of the surrounding rocks. According to the energy dissipation and energy release principles [39], as shown in equation (2), the maximum release rate of the destructive energy of the surrounding rocks appears in the direction of the minimum principal stress, which is unfavourable for the stability of the sidewalls of a tunnel. When the surrounding rocks are damaged by in situ stress, the energy release rate in the principal stress direction is proportional to the stored elastic strain energy and principal stress difference (Δσ). The energy release rate increases with increasing burial depth and principal stress difference (Δσ). The energy release equation is as follows [39]:

In the formula, K3 is the material constant, Ue is the elastic strain energy, and Gc is the critical value of the maximum energy release rate.

A rock mass composed of hard rock can store a very high elastic strain energy under high in situ stress conditions because it has high compressive strength and strong elasticity and brittleness parameters. The principal stress forms a stress concentration area in the surrounding walls when hard rock is excavated during the construction of a tunnel. When the principal stress of the surrounding rocks exceeds the ultimate compressive strength of the rock mass, the high strain energy of the rock mass can be released instantaneously, causing a rock burst, which typically includes splitting, shearing, and ejection phenomena. The uniaxial compressive strength of the limestone and basalt in the tunnel along the Lijiang to Shangri-La railway ranges from 70 to 100 MPa. The maximum horizontal principal stress (σH) along the railway is 25.09 MPa. The ratio of the tangential stress of the surrounding rocks and the compressive strength is approximately 0.6–0.8 MPa along railway line. Based on the Russenes criterion in Table 3 [40], there is a high probability of strong rock burst phenomena in the hard rock surrounding the tunnel project along the Lijiang to Shangri-La railway.

The weak surrounding rock has low compressive strength and high plasticity. The principal stress will readjust and redistribute in the surrounding walls when excavating this weak rock during the construction of the tunnel. When the stress exceeds the resistance of the surrounding rocks to bending, the surrounding rocks above and below the tunnel could become loose, producing vertical displacement and flexural phenomena. Soft rocks, including phyllite, shale, and slate, are widely distributed along the path of the tunnel project of the Lijiang to the Shangri-La railway. This soft rock has a uniaxial compressive strength of 8–25 MPa. The ratio of the maximum horizontal stress and compressive strength along the railway line is 0.16–0.50. Based on the strength-stress criterion in Table 4 [41], there is a high probability of slight to moderate deformation in the soft rock (e.g., phyllite, slate, and shale) surrounding the tunnel along the Lijiang to Shangri-La railway.

6. Conclusion

Based on the results of this study, the following conclusions can be drawn:(1)The in situ stress field along the Lijiang to Shangri-La railway is dominated by horizontal stress, which is described by σH > σv > σh and belongs to the strike-slip fault-based tectonic stress field. The maximum horizontal principal stress orientation is N10°–60°W, and the predominant direction is N15°W.(2)The in situ stress along the Lijiang to Shangri-La railway is classified as moderate to high. The maximum and minimum horizontal principal stresses and the horizontal stress differences (ΔσH-h) linearly increase with increasing burial depth. Although the lateral pressure coefficient (LPC, σH/) displays a certain discreteness in the range of 1.0–1.5, it exhibits a decreasing trend with increasing burial depth and gradually approaches 1 asymptotically.(3)The ratio of the maximum and minimum effective stresses along the Lijiang to Shangri-La railway is less than the critical value, indicating that the faults are still relatively stable at present. The risk probability of earthquakes is relatively low.(4)A small angle is observed between the maximum horizontal principal stress and the axial direction of the tunnel, which is unfavourable for the stability of the sidewalls of a tunnel. A high probability of a strong rock burst exists in the intact hard rock mass. On the other hand, slight to moderate deformation is expected to occur in soft rock masses (e.g., phyllite, slate, shale) in the tunnel project along the Lijiang to Shangri-La railway.

Data Availability

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

Conflicts of Interest

The authors declare that they have no known conflicts of interest or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

The authors thank colleagues of China Railway Eryuan Engineering Group Co. Ltd (CREEC) and Southwest Jiaotong University for kind support in the field. This work was funded by the Key Research and Development Program of China Railway Eryuan Engineering Group Co. Ltd (CREEC) (No. KYY2017006-(17-19)) and Fundamental Research Funds for the Central Universities (No. 2682020CX11).