Numerical Analysis and Deformation Mechanism Study on an Excavated High-Steep Slope of a Hydropower Station

State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining & Technology (Beijing), Beijing 100083, China School of Mechanics and Civil Engineering, China University of Mining and Technology (Beijing), Beijing 100083, China China ,ree Gorges Projects Development Co., Ltd., Chengdu, Sichuan 610041, China Changjiang Survey, Planning, Design and Research Co., Ltd., Wuhan, Hubei 430010, China ,ree Gorges High-Tech Information Technology Co., Ltd., Beijing 100033, China


Introduction
e instability of engineering rock mass often causes losses of human lives and infrastructure damage, imposing large financial costs to societies and many potential secondary disasters [1][2][3][4][5][6], especially the instability of slope engineering. As an example, more than 5000 landslides were reported in China in 2015, resulting in hundreds of deaths and billions of economic losses. erefore, it was especially important to analyze and evaluate the stability of the high-steep slopes. Among the stability analysis methods for geotechnical engineering, the development of slope stability analysis method is proved to be highly recognized [7,8].
In recent years, an increasing number of scholars have adopted the numerical simulation methods to analyze slope stability, while the typical numerical simulation method has been the FLAC3D strength reduction method, which was based on finite difference [9][10][11].
rough this method, stress deformation characteristics and the change of plastic area during slope excavation could be effectively analyzed, while slope support measures reasonably would be simulated. is would lead to a comprehensive study on slope stability. Consequently, this method has been widely utilized in both engineering research and design field [12][13][14][15][16][17]. e stability of high rock slopes constitutes the most important engineering geological and geotechnical problem for large-scale water conservancy and hydropower projects [18][19][20][21]. e stability control has become a significant key technology for the success or failure of water conservancy and hydropower projects, both affecting and restricting the development of hydropower resources as well as the construction of hydropower projects [22,23]. However, the appropriate mechanical parameters of slope rock mass constitute the foundation, upon which, the smooth progress of slope stability calculation could be ensured. rough the geological information analysis of the slope in Xiaowan Hydropower Station and according to the deformation and damage in the excavation of the accumulation slope, the shear strength parameters of the accumulation slope were inverted through the rigid limit equilibrium method. Moreover, the deformation parameters of the accumulation slope were inverted based on the monitoring and calculated displacements towards the slope direction regarding the accumulation slope key point [24].
In addition, a high number of scholars have analyzed the stability analysis of slopes of multiple hydropower stations.
rough FLAC numerical simulation analysis of the large excavation slope of Liyuan Hydropower Station, the deformation and failure characteristics of the slope under the conditions of excavation and normal water level were obtained, which provided a theoretical basis for the construction and operation of the hydropower station [25]. Based on the slope around the flood discharge tunnel at Yangqu Hydropower Station, a complicated three-dimensional (3D) numerical simulation and safety assessment were implemented in FLAC3D, while 3D slide arcs of good shape were obtained [26]. Taking the right bank slope in Dagangshan Hydropower Station in Sichuan province of China as the research background, the FLAC3D numerical simulation method was utilized to analyze its stability [27][28][29].
rough the FLAC3D program utilization as a computational analysis tool, the deformation mechanism of the bank slope of the Haizhou Open Pit Coal Mine Hydropower Station lower reservoir was analyzed. e results demonstrated that cyclical fluctuations and speeds of the reservoir water level affected its bank slope stability [30]. In order to study the deformation mechanism and control countermeasures of the high and steep slope at the Huangjinping Hydropower Station on the Dadu River in Sichuan province, the FLAC3D numerical simulation software was used to carry out a comparative study on the distribution of plastic zones prior to and following the slopes reinforcement [31]. In order to comprehensively understand the influence of deep-seated discontinuities on the left slope of Jinping I Hydropower Station and its global stability, both Sarma's method and the strength reduction method based on FLAC3D were adopted to simulate the plastic zone distribution. Also, the results indicated that the deep-seated cracks played a significant role in the slide surface location designation [32]. e time-dependent behavior of the left bank abutment slope at Jinping I Hydropower Station had a major influence on the normal operation and long-term safety of the hydropower station. A geomechanical model containing various faults and weak structural planes was established, while a numerical simulation was conducted under normal water load condition through FLAC3D processing to analyze the slope stability [33][34][35][36][37][38]. To research the macroscopic deformation of the left bank slope of the Baihetan Hydropower Station in Southwestern China, the unloading deformation of the left bank abutment rock mass was studied through FLAC numerical calculations [39].
In this paper, the left bank slope of the discharge tunnel of Wudongde Hydropower Station was taken as the research object. Based on the recommended values of rock mechanics parameters and slope safety monitoring results, the FLAC3D numerical model was used to inversely analyze the mechanical parameters of the slope rock mass. In addition, based on the results, the stability of the current and excavation slopes was studied, while the displacement field, the stress field, and the plastic zone distribution of four typical sections of the slope were analyzed. Combined with the geological characteristics and the numerical analysis results, the internal and external factors that affected the deformation law of the high-steep rock slope were described in detail. Furthermore, the corresponding treatment measures were put forward to ensure the stability of the slope of Wudongde Hydropower Station.

Topography and Geomorphology.
e Wudongde Hydropower Station dam site was located at the lower reaches of the Jinsha River at the junction of Huidong County in Sichuan Province and Luquan County in the Yunnan Province. It constitutes the first cascade hydropower station in the four cascades of the lower reaches of the Jinsha River ( Figure 1). e outlet of the flood discharge tunnel of the hydropower station was located at the left side of the left bank diversion tunnel and of the tail water of the power station. e elevation of the flip bucket was 860 m. A plunge pond was arranged downstream of the flip bucket to dissipate energy, while the plunge pool was formed by the Xiaohuashan platform excavation. e rock mass was discontinued and the impact resistance was poor. e plunge pool was of closed drainage structure type, and the foundation corridor was set to encircle the plunge pool. e bottom of the pond was 795 m high, 245 m long, and 110 m wide.
First, the upper cover of the high slope on the left side of the exit of the flood discharge tunnel was removed. Following, it was excavated to form an artificial slope. e maximum elevation of the overburden was approximately 1070 m, while the highest open line elevation of the slope was approximately 970 m. Also, the excavation elevation of the bottom plunge pool foundation corridor was 786.5 m, and the maximum height of the artificial slope was 183.5 m. e slope above the elevation of 925 m was dolomite excavation. Moreover, the strike of rock stratum was nearly N-S, the dip angle ranged from 28°to 35°, and the dip angle of the slope was approximately 80°. e artificial slope was a gentle slope in angle. e artificial slope below 925 m was the bottom layer of thin dolomite fold, the strike of rock stratum was nearly E-W, the dip angle ranged from 70°to 85°, and the dip angle of the slope was approximately 73°. e artificial slope was steep.

Formation
Lithology. e strata exposed before the excavation were the quaternary Huashangou cumulative body, the folded basement snowfall formation, the sedimentary cover Guanyinya formation, and the Dengying formation: (1) e total thickness of quaternary Huashangou deposit was approximately 50 m-70 m, which could be divided into three layers: the upper part was the collapse deposit gravel layer with a thickness of approximately 20 m; the middle part was the alluvial silty clay layer with a thickness of 10 m-26 m; and the lower part was the collapse and alluvial accumulation layer of erratic boulder and gravel with a thickness of 20 m-25 m. (2) e stratigraphic lithology of the Luoxue formation on the folded basement comprised of limestone and dolomite. With F 6 of Huashangou fault as the boundary, the upper side was thin interbedded with the 10th Luoxue formation (Pt 10 2l ) or very thin limestone. Also, the rock mass was hard, while the lower side was with thin interbedded or very thin dolomite of the Luoxue formation (Pt 2l ), while the rock was medium-hard. Influenced by the Huashangou fault, wrinkles were developed in the rock mass, while the properties of rocks varied. e strike direction was generally 250°∼290°, the tendency was S, and the dip angle was approximately 70°∼85°with local antidip.
(3) e Guanyinya formation (Z 2g ) of sedimentary caprock was mainly composed of gray thin interbedded or very thin dolomite. At the bottom of Guanyinya formation near fault F 6 , approximately 3 m∼5 m of thick carbonaceous shale was found, which was basically composed of medium-hard rock, with poor to fragmented rock mass. e Dengying formation (Z 2d ) was a shallow gray thick and medium-thick dolomite with generally exposed hard to medium hard rock, relatively complete. Sedimentary caprock strata strike was approximately 340°∼10°, while the tendency was E and the dip angle was 28°∼35°.

Geological Structure
e left-side slope of the flood discharge tunnel was penetrated by a large-scale fault F 6 . e fault F 6 was a normal fault, the vertical offset was approximately 420 m, the fault strike was 70°∼90°, the tendency was SE, the inclination angle was approximately 70°, and the fault width was generally 3 m∼8 m, mainly composed of broken rock. e slope fissures were generally undeveloped but locally developed, mainly constituting the surface fissures. e lengths of the fissures were generally 5 m∼10 m. e fissures were generally flat and rough, mostly closed or slight. A total of 213 blocks were found in the slope, constituting 153 blocks less than 1000 m 3 , accounting for 71.9%, while 54 blocks were from 100 m 3 to 1000 m 3 , accounting for 25.3%. Moreover, 6 blocks were from 1000 m 3 to 10000 m 3 , accounting for 2.8%. e maximum bulk volume was approximately 4691 m 3 . Among these blocks, 12 were basically stable blocks, 187 blocks were Advances in Civil Engineering of poor stability, and 14 blocks were of worst stability. e depth was generally approximately 2 m∼5 m, and the deepest point was at approximately 29 m. In summary, the overall quality of the rock mass on the left-side slope of the exit of the flood discharge tunnel was poor. e construction site of the left-bank slope project is presented in Figure 2.
Since the left-side slope of the flood discharge tunnel had sustained deformation during the slope excavation, the deformation law and stability of the high-steep slope during the long-term construction of the Wudongde Hydropower Project were the main concerns. erefore, the stability study of the left-side high slope of the flood discharge tunnel was of high significance for the safe construction and longterm operational safety of the Wudongde Hydropower Project.

Geological and Numerical Models
Based on the engineering geological conditions exposed by the excavation, an overall 3D geological generalization model of the slope was established based on the geological data. Based on the slope safety monitoring and geophysical prospecting results, the mechanical parameters of slope rock mass were inversely analyzed through a FLAC3D threedimensional numerical model. On this basis, the characteristics of displacement field, stress field, and plastic zone of the left-side slope at the outlet of spillway tunnel under excavation and unloading conditions were studied, while the overall stability of the slope was thoroughly analyzed.

Establishment of the Numerical Calculation Model.
According to the actual engineering geological conditions, combined with the distribution of four grades of rock mass and fault, the engineering geological model of the left-bank high slope of the discharge tunnel of Wudongde Hydropower Station was simplified.
A simplified 3D numerical model was constructed based on the actual geological conditions of the left-bank slope of the Wudongde Hydropower Station. e calculation range was 1650 m × 1000 m × 1550 m (X × Y × Z) and the highest elevation was 1760 m, where the X-axis was perpendicular to the inner side of the mountain, the Y-axis was positive with the parallel slope, pointing to the upstream side, while the Zaxis was positive with the vertical direction (right-handed system). Based on the requirements of practical problems and the computational capabilities, the elements of the faults, unloading relaxation zones, and different rock grades in the engineering area should be refined. e physical model of the generalization of the left-side slope is presented in Figure 3. e model adopts a meshing method combining tetrahedral elements and hexahedral elements. At the bottom and top of the model that are less affected by fault and slope excavation, the hexahedral element grid is used. In complex terrain areas, a tetrahedral cell grid is adopted, which is distributed in a gradual manner. In key research areas such as excavation affected areas and faults, the grid is finely divided, and the grid size is also determined according to the size of the volume unit. e model consisted of 753,839 elements and 147,449 nodes. e grid division of the model and the main stratigraphic lithology and structural plane distribution in the study area are shown in Figure 4. For the layered rock mass, a jointed elastic-plastic layer model was adopted, through which the shear failure and tensile failure along the layer and the rock mass could be simulated. In this model, the rock mass strength anisotropy was considered for both directions of the parallel and vertical planes.
Apart from the slope designed as free boundary, the rest of the computation model including bottom and all around is designed as fixed restricting boundary. Under the initial condition and regardless of tectonic stress, only initial stress field produced by weight stress is taken into consideration.
In order to better analyze the evolution law of the rock mass inside the slope, four typical sections were selected for study in the excavation affected area. e typical cross section position is shown in Figure 5.

Model Mechanical Parameters.
Based on the analysis of the internal and external deformation monitoring data of the slope rock mass, combined with the recommended values of the rock mass mechanical parameters provided by the geological survey, the incremental displacement intelligent inversion method was utilized to obtain the mechanical parameters of the rock mass and structural plane of the slope.
e inversion method to obtain the mechanical parameters is as follows: (1) select a number of typical geological profiles on the engineering slope and lay surface and deep displacement monitoring points near the profiles; (2) analyze the on-site displacement monitoring data and select the monitoring data with application value as the inversion data; (3) give the recommended values of the mechanical parameters of the rock mass and determine the parameters to be inverted, namely, the modulus of deformation, cohesion, and coefficient of friction; (4) according to different slope positions, determine different inversion key parameters; (5) set up different comparative working conditions according to the actual project; (6) establish the finite element model of the key section and take the optimization In view of the relaxation effect of rock mass with poor quality caused by excavation of high and steep slope, the quality of rock mass apparently decreased. During the mechanical parameters inversion, the cohesion of relaxed rock mass with large deformation was reduced by 20%-50% of non-unloaded rock mass, while the friction coefficient was reduced by 50%-60% of non-unloaded rock mass. In the model, the deformation modulus, friction coefficient, and cohesion of the rock masses of grades IV 2 , IV 1 , and III 2 for P 10 t21 and P t21 were mainly inverted. e mechanical parameters of the slope rock mass obtained by inversion are presented in Table 1.

Displacement Field.
As it can be observed from Figures 6 and 7, as the slope was continuously excavated, the height of the slope increased. Particularly, with the excavation of the middle and lower rock masses, the excavation and unloading relaxation of the middle and lower slopes weakened the supporting effect on the upper rock masses, while the magnitude of the slope displacement gradually increased. e deformation vector of the slope rock mass appeared to be outward and demonstrated a downward deformation tendency, while the horizontal displacement outside the slope exceeded both the vertical displacement and the upward and downstream displacements. In addition, the shallow surface layer of the slope was loosened due to excavation and unloading, while the rock mass displacement in this area exceeded the deep displacement, which gradually decreased within the slope. e part with large deformation of the slope was located at the hanging wall and foot wall of fault F 6 and the rock mass of grade IV 2 was of poor quality below 830 m. Among these, the fault F 6 steeply leaned to the slope interior, while the deformation of the upper rock mass was significantly higher compared to the lower rock, which exerted a squeezing effect on the fault and the lower rock mass, whereas the lower rock slope provided the reverse resistance. However, when the middle and lower sides of the slope were weakened, the poor rock mass was excavated and the slope was unloaded. e rock mass and mechanical parameters were reduced, while the upper part of the slope was deformed by the upper part pull and the upper part displayed an outward sinking deformation. e deformation of grade IV 2 rock mass in the hanging wall and footwalls of fault F 6 on the slope was generally 200 mm∼350 mm. Certain footwalls were severely affected by fault F 6 and fault fracture zone, while the regional deformation was of 350 mm∼380 mm. e slope, in an entity, was mainly deformed towards the horizontal direction (Xdirection), while it was horizontally directed to the slope exterior. e area with a high amount of deformation was on both sides of fault F 6 , in between 910 m and 800 m in elevation, while the deformation was of 250 mm to 350 mm. Among these areas, the deformation area of the upper wall of fault F 6 was relatively high, mainly because the rock mass on the upper wall of fault F 6 was mainly of grade IV 2 . Also, the distribution area was wide, while following the excavation and unloading of the slope, stress relaxation of the existing Advances in Civil Engineering slope rock mass occurred along with horizontal deformation outward of the slope. e amount of deformation in the direction of the slope parallel to the slope (Y-direction) was −50 mm∼68 mm. e deformation area was mainly distributed on the hanging wall of fault F 6 , which demonstrated that the hanging wall rock of F 6 was deformed towards the downstream side. Among these, the maximum deformation was 60 mm∼68 mm, located at an elevation of 850 m∼865 m, mainly because the rock mass in this area was of poor quality, while it was constrained by the foot wall and the rock mass on the front side of the spillway tunnel. Consequently, after the slope rock mass was unloaded, the parallel slope could only deform downstream. For the footwall rock of F 6 , the overall deformation in the direction of the parallel slope was low, while only a portion of the grade IV 2 rock mass was deformed upstream by the gravity and deformation of the overlying rock. In addition, the height drop the deformation was higher, while the maximum deformation at the elevation of 806 m was approximately 42 mm. For the deformation in the vertical direction (Z-direction) of the slope, the subsidence deformation from the elevation of 910 m to the rock mass at the top of the slope was mainly 40 mm∼100 mm; the settlement gradually decreased from top to bottom, while the deformation in this area was mainly affected by the gravity of the slope rock mass subsequently to excavation. e middle and lower slopes were characterized by upward deformation and mainly located at the lower plate of fault F 6 . e corresponding deformation was approximately 50 mm∼170 mm. Among these, the intersection between the berm at elevation of 813 m and the fault had sustained the highest upward  Advances in Civil Engineering deformation of 240 mm in length. e deformation of the rock mass in this area occurred mainly due to the deformation rebound caused by the excavation of the overlying rock mass. Figure 7 presents the displacement cloud diagrams of 4 typical sections following the slope excavation. rough the displacement cloud diagrams comparison of the rock masses of sections 1 and 2, the following conclusions were drawn: the maximum displacement of section 1 was approximately 280 mm, located at an elevation of 925 m∼865 m. Near the hanging wall and foot wall of the fault F 6 , the maximum displacement of section 2 was approximately 360 mm, located near the hanging wall and foot wall of the fault F 6 with an elevation of 865 m∼850 m, while the two parts with high deformation were all of grade IV 2 rock masses. It could be observed that although the geological structure, stratum lithology, and rock mass quality of the two areas were similar, the maximum deformation amount varied highly, due to different elevations, which indicated that the excavation depth had a significant influence on the unloading deformation of the slope under similar conditions. e higher the excavation depth was, the higher the amount of rebound deformation generated by the slope was. rough the displacement cloud diagrams comparison of sections 3 and 4 rock masses, it could be concluded that the rock masses of the two sections were not penetrated by fault F 6 , while the rock mass was good. Consequently, the overall deformation of the slope was relatively low. e maximum deformation of section 3 was approximately 240 mm, located near the elevation of 830 m. e mass of this part was IV 2 and the depth of the slope was high. erefore, this part of the rock body sustained high rebound deformation subsequent to excavation and unloading. e maximum deformation of section 4 was approximately 150 mm, located near the elevation of 876 m. e mass of this part of the rock mass was grade IV 1 , but because the excavation depth of this part of the slope was shallow and restricted by the downstream mountain mass, the overall deformation of this section was low. In summary, the overall deformation of sections 1 and 2 of rock mass was significantly higher compared to other sections, while the area with high deformation was mainly located near the hanging wall and footwall of the fault F 6 , which indicated the grade IV 2 rock in the hanging wall and footwall areas of the F 6 slope fault. e body was subjected to disturbances under the action of excavation unloading, while the deformation of the rock mass was also significant. In addition, the shallow surface layer of the slope was loosened due to excavation unloading, while the rock mass deformation in this area exceeded the deep deformation. Moreover, the deformation tended to decrease as the slope depth increased.

Stress Field Analysis.
e geological conditions of the left-side slope at the outlet of the spillway tunnel of Wudongde Hydropower Station were complex. Consequently, the stress field of the slope was relatively complex. e natural stress field of the slope was mainly controlled by the gravity of the slope and the geological structure. After the slope excavation, the slope was unloaded and relaxed, while the relaxation range was approximately distributed along the slope surface. e general depth of stress disturbance was 30 m to 60 m. After slope excavation, a certain degree of compressive stress concentration existed at the foot and at the junction of the slope. e stress magnitude was approximately 2 MPa to 5 MPa, while the local junction of the two slopes reached 12 MPa (Figures 8-10). Figure 8 presents the stress field characteristics of the left-side slope of the outlet of the spillway tunnel under natural stress. e initial stress field had the following characteristics: (1) Maximum principal stress: it could be observed from the cloud diagram that the maximum principal stress field of the slope increased with depth from the surface to the interior of the slope in the natural state. e maximum principal stress ranged from 0 to 40 MPa and uniformly varied with depth. At the top and foot of the slope, due to the influence of weathering and unloading, a few tension stress zones appeared, but the influence depth was relatively shallow, which indicated that the slope was in a state of stress stability under natural conditions. At the foot of the slope, due to the influence of fault F 6 , a local stress mutation occurred, but its influence depth was shallow and the scope was small, which did not affect the stability of the slope in its natural state (2) Minimum principal stress: it could be obtained from the cloud diagram that in natural state, the minimum principal stress increased as the elevation decreased, while the overall change trend was more uniform. Its magnitude range was 1 MPa to11 MPa. e minimum principal stress in the shallow part of the slope mostly ranged between 0 and 1 MPa, which indicated that the surface of the slope was in tension stress state, due to weathering and unloading. Similarly, in the fault F 6 area at the foot of the slope, a small range of stress mutation existed, but it did not affect the overall stability of the slope.
After the slope was completely excavated, the slope was unloaded and relaxed, while the relaxation range was roughly distributed along the slope. e general stress disturbance depth was 30 m to 60 m. A certain amount of compressive stress concentration existed at the foot and boundary parts of the slope. e stress level was approximately 2 MPa to 5 MPa, and the junction of the two slopes reached 12 MPa (Figure 9). A large range of tension stress area existed within the slope, which mainly occurred on all levels of slope surface, while the tension stress of the slope was generally 0 MPa to 0.4 MPa. Also, the maximum tension stress value was 1 MPa, while the larger part of tension occurred near the slope opening line and the junction area, where the lithology highly changed. e maximum and minimum principal stresses of the four typical sections are presented in Figure 10.

Strain Zone Analysis.
e scope and depth of the plastic zone of the slope were relatively high. e plastic zone of the rock mass gradually increased with excavation. e plastic Advances in Civil Engineering zone was mainly distributed on the excavated slope surface. e plastic zone was concentrated on the fault F 6 and the tectonic influence zone. e rock mass basically entered the plastic state, mainly as compression-shear failure. Tensionshear failure occurred on the surface of faults and slopes and local rock masses, as presented in Figure 11.

Deformation Mechanism of Slope and Reinforcement Suggestions
6.1. Slope Deformation Evolution. e process of excavation and unloading of engineering slope is always accompanied by stress adjustment, deformation response, and damage of rock mass. On the one hand, the rock mass of the slope was subjected to dynamic damage caused by blasting. On the other hand, excavation changed the boundary conditions of the slope, and the rock mass sustained unloading deformation. During the long-term geological evolution, strain energy accumulated within the rock mass of the slope, while the slope was in a relatively stable state prior to engineering activities. e excavation disturbed the initial state of the slope, and the rock mass was adjusted accordingly. A certain amount of strain energy, accumulated previously, was released, whereas another amount was transferred to the deep part of the rock mass, resulting in stress redistribution and accompanied by the deformation of the slope. e excavation resulted in strain energy release and unloading relaxation. e instantaneous deformation caused by unloading was rebound, which produced static and dynamic disturbance to the slope and rock mass, causing the opening and dislocation of the original structural planes in rock mass, as well as the expansion of primary joints or new cracks,   especially near the geological weak planes, such as faults and lithology, soft rock, and strong unloading relaxation zone. e quality, the mechanical parameters, the bearing capacity, and the stability of retained rock mass might be significantly reduced. e evolution of slope deformation was microdamage and rupture development of rock mass, along with continuous decline of macromechanical strength, accompanied by the formation, development, and penetration of slip surface ( Figure 12). For slowly changing rock slopes with time-dependent deformation characteristics, the sliding surface of such slopes gradually evolved over a long period of time. For slopes with progressive deformation caused by excavation and unloading, if no controllable structural surface existed, the development of deep sliding surfaces would end in the early or middle stage of development, without entering the final stage of large-scale cumulative damage. Conversely, if the deformation was fully developed, the sliding surface evolution would have entered an irreversible state. At present, the left slope of the spillway tunnel exit is still in the stage of slow deformation, while no uniform potential sliding surface was formed. e characteristics of the left-side slope of the Wudongde spillway tunnel exit included (1) soft foundation effect, (2) antitilt cutting of fault F 6 , (3) high and steep slope formed by excavation, (4) high unloading relaxation depth, (5) relaxation deterioration of "locking section" at the slope middle, (6) nonuniform differential deformation, and (7) lower deformed rock mass pulling the upper slope to settle.

Analysis of Deformation Mechanism.
e excavation mechanical response model of the left-side slope at the outlet of the spillway tunnel was mainly controlled by factors, such as the lithology of the slope, rock mass structure, excavation slope shape, excavation unloading scale, and construction disturbance effect. According to the geological conditions revealed by the excavation of the slope, the upper part was ZSSR (e -6)  thin dolomite interspersed with sandy mudstone, shale, and thick dolomite, while the lower part was thin dolomite and limestone. A large-scale fault F 6 existed within the slope, and the trend was NNW160°∼180°. Plastic shale debris was often found near the main section. e mass of the slope rock mass was mainly IV 2 and IV 1 , while the low amount was III 2 , III 1 , and V. e slope unloading relaxation zone was generally 21.4 m deep, and the deepest part was at 29.2 m. e main causes of slope deformation included the following: (1) Internal cause 1: due to the soft foundation effect, the rock mass had poor lithology, as well as a wide range of unloading relaxation, while the quality of rock mass was apparently reduced. e rock mass quality of slope was mainly IV 2 and IV 1 (Figure 13), a portion of which was III 2 . Geophysical detection of slope demonstrated that the acoustic wave in relaxation zone was 2.5 km/s to 3.2 km/s, the acoustic wave in nonrelaxation zone was 3.0 km/s to 4.5 km/s, and the local area was 3.0 km/s to 5.0 km/s. A "soft foundation effect" existed within the lower part of the slope to a certain extent, while similar deformation features often appeared on the slope of the "upper hard and softer" structure (2) Internal cause 2: cutting and weakening effects of fault F 6 and caprock Z 2g on slope. e fault strike of F 6 intersected with the slope strike at a high angle, steeply leaning into the slope, demonstrating an antidipping rock mass structure. F 6 was a normal fault, the caprock was staggered, the vertical offset was approximately 420 m, and the main fault zone was composed of mud-clastic. Field monitoring demonstrated that the highly deformed part of the slope and the large stress part of the anchor cable were located at the hanging wall and foot wall area of fault F 6 . e numerical calculation demonstrated that along with the gradual exposure of F 6 on the slope surface, the downstream side slope of the fault and the upstream side slope (rock mass connected with front slope of spillway tunnel) were separated to form a weak boundary of the downstream side slope of the fault, which weakened the three-dimensional side restraint effect of the front slope of the spillway tunnel on the entire left side of the slope. e deformation of downstream side slope of fault F 6 extruded the fault and its footwall slope, demonstrating the self-compression of the weak part of the fault and the deformation slip towards the outside side of the slope. Also, the phenomena of high lateral deformation and low internal deformation of the unloading deformation of the hanging wall and foot wall itself were apparent. e upstream side slope bore the compressive load transmitted through the Advances in Civil Engineering fault, while the slope connected with the front slope of the spillway tunnel acted as support and impedance to the fault expansion and its downstream side slope. In addition, the Guanyinya formation was a rock mass with poor integrity and low deformation modulus, while the joint weakening with the F 6 fault also affected the slope deformation of local convex terrain near the slope opening line (3) External cause: the scope and depth of excavation unloading relaxation were high, the scope of slope excavation unloading disturbance was high, and the unloading speed was high. e wave velocity decreased along the hole depth, and, in part, continuous low-wave-velocity zones existed in the interior, except the surface relaxation zone of the slope. e strong unloading relaxation depth of the slope was high, while the relaxation zone was generally 21.4 m and the deepest point was at 29.2 m. A high amount of new cracks occurred within the deep rock mass, and the local continuous distribution led to the wave velocity reduction zone, which caused the overall wave velocity to be low and the relaxation zone to be high. With time stress adjustment, crack initiation, evolution, and continuous propagation were caused. e relaxation zone extended towards the depth with excavation, while the quality and mechanical parameters of rock mass were deteriorated.
For example, the long hole acoustic wave test at an elevation of 853m shows that between June 2016 and November 2016, the test relaxation depth is about 2.5 m, and the average acoustic wave velocity of the rock mass in the relaxation zone is 3100 m/s∼3300 m/s. After the slope was excavated from 830 m elevation to 806 m elevation, the relaxation depth became significantly larger. e relaxation depth tested on May 6, 2017 reached 11.8 m, and the average acoustic wave velocity of the rock mass in the relaxation zone was 2500 m/s∼2900 m/s. e relaxation depth is shown in Figure 14. During the blasting operation on the slope, the instantaneous strong dynamic load disturbance caused the reduction of rock mass properties, especially the faults and areas with poor lithology quality, which caused the sudden increase of deformation in these areas in a short time. As shown in Figure 15, when performing blasting operations at an altitude of 830 m, the M14 multipoint displacement gauge is located at an altitude of 850 m, which is close to the blasting operation point and fault, and the rock mass quality level of this area is low, so it caused a sudden change in slope displacement in this area. However, the M13 multipoint displacement gauge is located at an elevation of 880 m, relatively far from the blasting operation point and fault, and the rock mass quality level in this area is relatively high, so although the increase in displacement occurs, the increase in displacement is relatively low. At the same time, after the lower rock mass is excavated, the support effect of the lower rock mass on the upper rock mass is significantly reduced, which will also increase the overall deformation of the upper rock mass.   In summary, the geological conditions of the left-side slope at the outlet of Wudongde spillway tunnel were complex, and the rock mass quality was poor. Under the condition of rapid excavation and unloading, the relaxation deformation depth of rock mass exceeded the general rock slope relaxation deformation depth. e slope deformation and the stress increment of anchor cables were generally high, which had an adverse effect on the slope stability. With the excavation of slope engineering, corresponding measures should be taken to improve the overall stability of the left-side slope.

Reinforcement and Treatment.
With the continuous excavation of the later slope engineering, the stability of the slope would still be reduced under the influence of engineering disturbance under the current conditions. erefore, corresponding measures should be taken to improve the overall stability of the slope. It was suggested that unbonded prestressed anchor cables should be added in the area, where the slope deformation and anchor cable overload could occur within the range from 850 m to the opening line of engineering slope to strengthen anchorage. Moreover, overload anchor cables should be loosened, and the deformation of slope under continuous excavation could be reduced. e concrete should be sprayed upon the slope near the F 6 fault, while the cracked and broken rock mass should be cleaned up, while the steel mesh should be hang and the concrete should be sprayed. Anchor rods of the system were arranged near the fault, perpendicularly to the surface of the slope, which were tied tightly with the reinforcement of the hanging net, while the drainage holes of the system increased. A certain number of anchor holes should be selected for borehole acoustic wave and borehole color TV detection in the construction of new anchor cables in the engineering slope area. In addition, a certain number of anchor cable dynamometers should be arranged. Also, the constant resistance and high deformation monitoring anchor cables should be added to the large deformation area to realize the intelligent and dynamic monitoring of the anchor cable force within the slope.

Conclusions
Based on the above analysis, a few conclusions can be drawn: (1) Combined with the recommended values of rock mass mechanical parameters and slope safety monitoring results, the rock mass mechanical parameters were obtained through inversion analysis. (2) e numerical simulation results demonstrated that after the excavation is completed, the larger slope deformation area is distributed in the hanging and foot wall of fault F 6 and the grade IV 2 rock mass below 830 m elevation, and the deformation of the hanging wall rock mass is significantly greater than that of the foot wall, and the maximum deformation range is 350 mm∼380 mm. e depth of stress disturbance is 30 m∼60 m, and there is a certain degree of stress concentration at the foot of slope; the stress level is about 2 MPa∼5 MPa. e plastic zone of the slope is concentrated in the fault F 6 fracture zone and the area with poor rock mass, and the rock mass of the slope basically entered into the plastic state, mainly demonstrating compressive shear failure. (3) e internal cause of the deformation mechanism is the overall poor quality of the slope rock mass, which is represented by the structure of "upper hard and lower soft." At the same time, fault F 6 and stratum Z 2g have the effect of cutting and weakening the slope. e external factors are the large unloading relaxation range, depth and disturbance range, and the fast unloading speed. (4) As the slope project continues to advance, the stability of the slope will still decrease, so the corresponding treatment measures such as slowing slope ratio, consolidation grouting, and deep anchor cable support are recommended below 806 m. At the same time, it is proposed to use the remote sliding force monitoring system to monitor the sliding force at the deep part of the slope in real time, which will provide an effective reference for the stability assessment and treatment of the slope of the hydropower station.
Data Availability e data are all available and have been explained in this article. Readers can access the data supporting the conclusions of the study.

Conflicts of Interest
e authors declare that no conflicts of interest exit regarding the submission of this paper.