Stability Assessment of the Excavated Rock Slope at the Dagangshan Hydropower Station in China Based on Microseismic Monitoring

A high-resolutionmicroseismic (MS) monitoring systemwas implemented at the right bank slope of the Dagangshan hydropower station in May 2010 to analyse the slope stability subjected to continuous excavation. ,e MS monitoring system could real-time capture a large number of seismic events occurring inside the rock slope. ,e identification and delineation of rock mass damage subject to excavation and consolidation grouting can be conducted based on the analysis of tempospatial distribution of MS events. However, how to qualitatively evaluate the stability of the rock slope by utilizing these MS data remains challenging. A damage model based on MS data was proposed to analyse the rock mass damage, and a 3D finite element method model of the rock slope was also established. ,e deteriorated mechanical parameters of rock mass were determined according to the model elements considering the effect of MS damage. With this method, we can explore the effect of MS activities, which are caused by rock mass damage subjected to excavation and strength degradation to the dynamic instability of the slope. When the MS damage effect was taken into account, the safety factor of the rock slope was reduced by 0.18 compared to the original rock slope model without considering the effect. ,e simulated results show that MS activities, which are subjected to excavation unloading, have only a limited effect on the stability of the right bank slope. ,e proposed method is proven to be a better approach for the dynamical assessment of rock slope stability and will provide valuable references for other similar rock slopes.


Introduction
e stability of rock slopes subject to excavation is a major challenge in geotechnical engineering during both the design and construction stages. Many achievements have been made to evaluate the excavation influence of rock slope [1][2][3]. For instance, Jin proposed that steep terrain had a significant influence on the dynamic behaviour of the slopes under blasting [4]. Guo et al. investigated the relationship between slope type and slope effect value [5]. Wayne discussed the mechanical processes in the rock mass around a kimberlite volcanic excavation that occurred as the rock responded to the formation of the excavation [6]. Lu et al. introduced a modified tensile-compressive damage model to the simulation of the excavation damaged zones of the Xiluodu high rock slope [7]. is model reproduced the excavation damaged zones for high rock slope excavation with the dynamic finite element method. Yang et al. analysed the dynamic behaviour of a high road cutting rock slope under blasting excavation by FLAC 3D [8]. Sato et al. described the results of an in situ excavation disturbance experiment using mechanical excavation with a boom header at the Tono mine and compared the findings with the results obtained using blasting excavation [9]. Kentli and Topal conducted field and laboratory researches to determine the engineering geological properties of rocks exposed along motorways and to assess the excavation behaviour and stability of the cut slopes [10].
Unlike the conventional monitoring methods, which focus only on the information regarding surface movements, the microseismic (MS) monitoring data would be sufficient for us to understand the microscale changes occurring inside the rock mass [11,12]. Application of the MS technique provides in-depth knowledge of the mechanism associated with rock slope failure through the results of excavated rock mass. e MS monitoring technique is a useful tool for detecting unstable areas of rock slopes and assessing the corresponding hazards. Considerable researches have been conducted on the application of MS/acoustic emissions (AEs) to predict rock slope failures. For example, Ma et al. proposed the stability analysis and reinforcement effect of the slope based on MS energy density and the magnitudefrequency relation [13]. Kolesnikov et al. proposed seismo-AE for monitoring slope stability in the north by filtering and polarization analysis [14]. Senfaute et al. examined the MS signals recorded 15 h prior to a rock fall in Mesnil-Val, France, which offered the possibility of using the MS monitoring system to monitor high-risk sections of coastline and understand the mechanism of cliff failure [15]. Amitrano et al. identified a power law acceleration of seismicity rate and energy before a cliff collapse [16]. Eberhardt et al. presented preliminary MS monitoring results recorded through the Randa rockslide laboratory and combined several conceptual and numerical models in terms of the tensile and shear mechanisms of slope failure [17]. Shiotani introduced a prospective method to install the AE sensor in a rock slope to evaluate the slope stability [18]. Spillmann et al. established a seismic monitoring to detect active faults of landslides [19]. Numerous studies have investigated the relationship between rock rupture events and mechanical parameters based on the notion that the AEs and damage variables are consistent. Several researchers introduced damage variables into damage mechanics to derive the relationship between the number of AEs and the damage variable and to study the mechanical properties of rock materials during the progression from microscopic damage to macroscopic rupture [20]. However, a comprehensive approach that could combine the effect of MS damage and numerical modelling to analyse the rock slope stability is rarely discussed.
is study presents a practical application of MS monitoring for stability assessment of the excavated rock slope at the Dagangshan hydropower station, Southwest China. e brief description of the right bank slope of the hydropower station is presented in Section 2. In Section 3, the MS activity characteristics of the right bank slope are described, including the MS monitoring system and tempospatial distribution of MS events. Section 4 focuses on the stability analysis of rock slope considering the effects of MS data. e conclusions are followed in Section 5.

Brief Description of the Right Bank Slope of the Dagangshan Hydropower Station
e right bank slope of the Dagangshan hydropower station is up to 700 m high and is located on the Dadu River near Ya'an city (29°58′N, 103°01′E), approximately 350 km southwest of Chengdu, Sichuan province, China. e project has a double-curvature arch dam with a maximum height of 210 m and a total installed capacity of 2,600 MW, which controls a drainage area of approximately 62,727 km 2 or 81% of the total drainage area of the Dadu River. Meanwhile, it is one of the highest excavated slopes in hydroelectric constructions that are currently used along the main stream of the Dadu River. e excavation height exceeds 300 m. e layouts of the natural slopes and the right bank slope after excavation are shown in Figure 1. e excavated right bank slope is characterized by complex geological structures, large-scale excavation, a high stress level, and strong rock unloading. In particular, the numerous faults, dikes, unloading fracture zones, and joints, such as β 85 , β 4 , c L5 , c L6, f 231 , and two deep unloading belts of XL-316 and XL-915, form the major potential unstable blocks in these areas. A typical transverse section VI-VI of the right bank slope is shown in Figure 2 to provide an in-depth understanding of the complexity of the geological conditions. e unstable blocks and large scale of the excavation may result in the potential instability of the rock slope. In fact, some tensile fractures occurred at a shallow depth. is phenomenon has rarely been encountered in other hydropower projects elsewhere.
Numerous monitoring data recorded by conventional measurement apparatuses, such as multipoint extensometers and convergence gauges, were obtained at different elevations of the right bank slope. Several macrocracks occurred on the surface of the right bank slope during the continuous excavation, as illustrated in Figure 3. e similar increasing trend of the deformation of the rock slope at an elevation of 1,101 m was obtained by using the convergence gauges and multipoint extensometer in July 2010. e same growth shown in Figures 4 and 5 is attributed to the continuous excavation of the rock slope.

MS Monitoring
System. Macrocracking appears on the slope surface due to the accumulation of microscopic fractures in the surrounding rock mass of the excavated slope. e MS monitoring system installed on the right bank slope of the Dagangshan hydropower station was manufactured by the ESG Corporation, Canada. e system is comprised of a Hyperion digital signal processing system, Paladin digital signal acquisition system, and twenty-four uniaxial acceleration transducers installed in the boreholes drilled from the sidewalls of the tunnels inside the right bank slope, which contain sophisticated hardware and powerful software. e real-time MS monitoring system layout installed on the right slope was developed based on network technology and built on the distributed data acquisition units with central access to the online data processing system. e twenty-four acceleration sensors covering the dimensions of 600 m (length), 600 m (width), and 400 m (depth) were connected to four Paladins acquisition systems using the copper twisted-pair cables ( Figure 6). e four substations were connected in series to the Hyperion processing system with optical fibre cables. e twisted-pair cable, fibre optics, and radio telemetry were typically configured during the installation process for the MS monitoring systems. Once the amount of sensors to be 2 Advances in Civil Engineering installed in the right bank slope was determined, a reasonable arrangement of the sensors was implemented to make full use of the MS monitoring system. e ultimate optimal arrangement of the sensor array is shown in Figure 7.
All signals were transmitted to the data analysis system installed on the right bank slope in the container at an elevation of 1,035 m due to geographical limitations. e data were digitized by the processing system with a sampling frequency of 10 kHz, and preliminary event detection was performed when the signals recorded by the four substations exceeded a given threshold, which was identified by the short-time average versus long-time average (STA/LTA) algorithm. e location of a seismic event wasdetermined by using the constant P-and S-wave velocities, which were initially estimated from artificial blasting tests.
Different wave velocities would result in different first arrival times on elastic waves to sensors; therefore, the wave velocities set in the MS monitoring system exert a great influence on positioning of MS events. e wave velocities in rock mass within the scope shall be calibrated prior to the positioning error test of the MS monitoring system. According to the provided rock sample data, the measured wave velocities were inverse calculated and calibrated using data analyse software from the ESG MS monitoring system. e average wave velocities in rock mass are determined as P wave 4400 m/s and S wave 2540 m/s. e detailed specification can be found in [21]. An event can be located by the first P-wave arrivals at 4-24 sensors and S-wave arrivals at 1/4-1/2 based on the simplex algorithm [22]. Moreover, the system can evaluate seismic source parameters, such as the energy index, apparent stress, magnitude, and apparent volume.

Tempospatial Characteristics of Microseismicity.
e MS monitoring system was established on the right bank slope of the Dagangshan hydropower station to analyse the cause of microcracks in the rock slope subject to excavation. MS monitoring is a 3D and real-time monitoring technique, and it has been improved with the rapid progress of computer and data acquisition technology compared with traditional instrumentations. e MS monitoring system installed at the Dagangshan hydropower station was first adopted to investigate the macrocracks subject to the continuous excavation of the rock slope [21]. After more than one year of continuous monitoring and analysis, a total of 1,337 seismic events with M w values ranging from −2.0 to −0.2 were detected along the right bank slope. Figure 8 shows the correlation between the accumulative energy release and occurrence frequency of MS events. e daily frequency of MS events ranges from 1 to 10, with small bursts every 3-5 days. e mean frequency between February 25 and October 10, 2011, was approximately 10 events per day. e majority of the events were recorded during excavation. On-site observation and analysis indicated that MS activities were caused by stress redistribution in the deep slope subject to continuous excavation. e frequency of the MS events increased drastically after blasting excavation. In addition, there was no MS signal during the period from December 1, 2010, to February 20, 2011, because there was no excavation activity during the Chinese Spring Festival. After the break period, the recovery of excavation on the rock slope resulted in an increased number of MS events. Many MS events were recorded in May 2011 with continuous excavation of the right bank between 1,135 m and 1,240 m. e number of recorded MS events decreased significantly after the excavation ceased in September 2011. e histogram of recorded MS events illustrates that the MS quantity is clearly attributable to the excavation unloading. e stress redistribution during the rock slope excavation is the primary source of the rock slope stability. e accumulative trend of seismic events is in good agreement with the rock slope excavation, as illustrated in Figure 8. Advances in Civil Engineering Figure 9 shows the spatial distribution of MS events recorded between May 4, 2010, and September 30, 2012. e figure highlights the distribution and migration characteristics of microfracture clusters in the deep rock mass of the right bank slope. Figure 9(a) shows that the majority of MS events occurred along faults XL-316 and f 231 and at elevations between 1,070 and 1,255 m, thus confirming the relationship between MS events and geological structures. Furthermore, MS event clusters were mainly concentrated in the preexisting geological structures. Figure 9(b) shows that the seismic deformation has a close correlation with the preexisting geological structures, such as faults XL-316 and f 231 . e maximal seismic deformation is mainly located along the aforementioned faults. e seismic deformation is prominent in the areas of preexisting faults and man-made structures. Moreover, continuous excavation may result in energy unbalance and stress redistribution in the deep rock mass of the rock slope, which will cause the movements of weak structures and slippage of the faults. Figure 9(c) shows the energy release density of MS at the rock slope subject to excavation. e maximal density of energy dissipation was located along the faults (i.e., faults XL-316 and f 231 ) in the deep rock mass, which is consistent with the spatial distribution of the seismic deformation contours (Figures 9(a) and 9(b)). us, the occurrence of MS event can be used to identify and delineate the activated geological structures. e energy release mainly concentrates on the potential sliding surface, approximately 80 m away from the excavated slope surface. e excavation-induced deformation can be estimated through the location, time, energy release, and scalar seismic moment magnitude of MS events. e monitoring results show that there is a close correlation between MS frequency and its accumulated energy release, which can be used to quantify the extent of damage to rock mass inside the rock slope. is in turn allows the identification and determination of the potential fracture regions and sliding surface in the excavated slope and thereby forecasts the slope failure.
Clear unloading effects are associated with stress redistribution, which results in the deterioration of the mechanical properties of the rock mass. erefore, the strain increasingly adjusts due to rock mass degradation, which in turn generates macrocracks on the rock slope surface. Numerous investigations and laboratory tests have demonstrated that the mechanical parameters (i.e., elastic modulus, cohesion, and internal friction angle) of the rock mass in the damage zone show a marked decrease compared with the macromechanical properties before excavation [23,24]. Energy release inevitably occurs in the form of seismic signals before its irreversible damage. Seismic signals, which include abundant information about the internal changes of rock mass, can be captured and analysed. erefore, the MS monitoring system is used to determine the spatial location and magnitude of MS events through seismic wave information and to estimate the development of nonlinear strain regions induced by MS activity. e results can explain the macrocracking phenomena on the excavated slope surface. Furthermore, they can provide a new approach to investigate spatial fracture patterns induced by excavation and forecast the stability of the excavated slope. In addition, the MS monitoring results can be used for an in-depth analysis of the failure characteristics of the excavated slope. Excavation-induced MS damage leads to the formation of cracks on the rock slope surface. Figure 10(a) shows that the MS events with a large magnitude and high energy mainly focus on the locations near the large cracks on the slope surface. In contrast, the MS events with a small magnitude and low energy typically concentrate along the small cracks. With further excavation, the scope of the damage zone will gradually expand. e macrocracks on the slope surface are formed by the accumulation of tiny fractures in the deep rock mass of the rock slope, as shown in Figure 10(b). ose small fractures are mainly caused by excavation unloading. e successive appearances of the macrocracks indicate the gradual evolution of damage in the rock slope. A series of tiny fractures gradually evolve and finally lead to macroscopic damage and deformation, as illustrated in Figure 10. e contours of seismic deformation in the selected periods along with the statistical analysis results can provide some information to understand the fracturing behaviour inside the rock slope (Figures 9 and 10). us, the MS monitoring system can offer a 3D visualization of rock mass damage and provide valuable information to understand the microcracking behaviour in the deep rock mass of the rock slope.  elastic material containing an initial random distribution of microfractures [25,26]. e mechanical parameters of the model elements are assumed to follow the Weibull distribution, which includes the statistical variability of the bulk failure strength in the RFPA code. In the elastic damage mechanics, the elastic modulus of the element may decrease gradually as the damage expands until the stress increases to the threshold value. e elastic modulus of the damaged element can be expressed as follows [26]:

Stability Analysis of the
where ω is the damage variable, E is the elastic modulus of the damaged material elements, and E 0 is the elastic modulus of the undamaged elements. E, E 0 , and ω are all scalars [27,28].

RFPA 3D -Centrifuge Method.
Recently, an equivalent method for mechanical analysis in geotechnical engineering was provided by centrifuge modelling. Some achievements have been obtained in mining and civil engineering projects [29,30]. To develop a powerful and effective alternative approach, we have implemented the fundamental principles of centrifuge modelling in the RFPA 3D program, thus forming the RFPA 3D -centrifuge code.
In this code, the stress field in the mesoscopic elements is obtained by adopting the FEM. To simulate centrifuge loading, the mass density of mesoscopic elements is increased gradually. e constitutive law of the mesoscopic elements is described via elastic and elastic damage mechanics. e ultimate tensile strain and envelope in the Mohr-Coulomb criterion are used as failure thresholds [28].
In centrifuge experiments, volumetric force, which is an artificial simulation of gravitational force, is added to   e weight of the original body is defined as c p � ρg, where g represents the gravitational acceleration and ρ the density of rock mass. e specific weight for the model is c m � ρa, in which a represents the overall acceleration vector. In accordance with the condition that the model and its original body should experience the same stresses, that is, σ p � σ m , it conforms to the equation, which is ρgh p � ρah m . en,

Advances in Civil Engineering
In the present study, to establish the model, the aforementioned fundamental principles of centrifugal loading are introduced to represent the failure process and safety factors. e self-weight of the basic mesoscopic elements increases gradually and linearly according to a specified step size. e balance between internal and external forces is identified by performing an iteration calculation at each step. With selfadaptive loading, the failure analysis will not stop until the occurrence of macroscopic failure in the rock slope. When the number of failure elements reaches the maximum, the numerical calculation is able to identify the sliding failure surface. e ratio of the self-weight of an element upon failure to its initial weight is defined as the safety factor K [24]: where S step is the loading step, where the number of the failed elements reaches the limit, Δ g represents the centrifuge coefficient, and c is the material unit weight.

Rock Mass Degradation Constitutive Considering the Effect of MS Damage.
Energy dissipation is a vital characteristic of rock failure. Rock mass damage and strength degradation are typically related to energy dissipation, which can reflect the reduction in rock mass strength. e cohesion of the rock mass element is assumed to be lost completely provided that the damage quantity inside a rock mass element reaches the maximum after specific energy dissipation. Rock mass deformation and failure are related to the synthetic result of energy dissipation and energy release, which is the stress state of the rock mass element and lithology. Energy release is the intrinsic cause of failure of the rock mass, whereas energy dissipation induces rock mass damage and causes the degradation of material properties and strength [31].
In an MS event, the ratio of seismic energy (ΔU) to the total energy release (ΔU) during the microfracturing of the rock mass can be defined as the seismic efficiency η, that is, where η is typically calculated from the ratio of radiation energy recorded by a MS monitoring system to the total energy of rock microfracture, equalling 0.1 and 1, respectively. e seismic efficiency at the right bank slope of the Dagangshan hydropower station can be obtained by comparing the seismic wave energy with the total energy released by explosive charges during the blasting tests. For simplicity, the radiation energy and loss of heat energy are ignored. e energy released from the blasting site is transformed completely into elastic wave energy in the rock mass.
Based on the aforementioned rock mass mechanical parameters revised by the energy damage criterion, 1,587 rock blasting events were selected from the beginning of MS monitoring on May 4, 2010, until the end of the slope excavation on September 2, 2011. Statistics of the seismic source parameters of the blasting events were recorded by the MS monitoring system, as shown in Table 1. e total energy of each blasting event is 216,000 kJ using #2 rock emulsion explosive with a charge of 72 kg for each blast [32]. e ratio of the mean seismic radiation energy (647 J) recorded by the MS monitoring system to the total energy was determined to be approximately 0.003%, which represents the equivalent seismic efficiency on the right bank slope of the Dagangshan hydropower station. e releasable strain energy U e of a rock mass element can be calculated as follows:

Advances in Civil Engineering
A rock mass damage variable D can be defined as the ratio of the assigned energy ΔU of the rock mass element based on the seismic radiation energy to the total releasable elastic strain energy U e [31]: Substituting (5) into (6) yields where E 0 and v 0 represent the initial elastic modulus and Poisson's ratio of the rock mass element, respectively. σ i (i � 1, 2, 3) represents the three principal stresses acting on the rock mass element.
erefore, the corresponding cohesion c r and elastic modulus E r of the rock mass element, considering the seismic energy dissipation, can be expressed as where c 0 and E 0 represent the initial cohesion and elastic modulus of the rock mass element, respectively. Numerical simulation of the rock slope excavation can be undertaken to consider the effects of MS damage by degrading the mechanical properties of the rock mass according to the damage variable D.
To determine the effect of MS damage on rock slope stability, the following assumptions must be made regarding the damage to the rock mass element [24], which are also the limitation of this study: Based on these assumptions, a flowchart is generated to illustrate the "feedback analysis" of the rock slope stability considering the effect of MS damage, as shown in Figure 11.

Numerical Modelling.
e 3D model simulates the topography and prominent structures of the two large unloading fissure zones XL-316 and XL-915 and fault f 231 . e model covers 690 m across the river, 400 m along the river, and a height of 1,060 m from the upstream area to the downstream area, as shown in Figure 12(a). Different colours represent different mass models, as shown in Figure  12(b). e mechanical parameters are listed in Table 2 [33]. e model implemented in RFPA 3D uses a bottom-up technique that subdivides the 3D model into 1,118,264 hexahedron elements and 1,153,843 nodes, as shown in Figure 12(c). e model conducts encryption processing on the shallow slope. e number of elements in the shallow slope comprises up to approximately 93% of the 3D model. Normal constraints are applied in the computation field. Fixed hinges are used for the bottom boundary, and the slope surface is free. e centrifugal load coefficient in the calculation is 0.005. e finite element model entities and mesh subdivision are presented in Figure 12. e C programming language is used to convert the interface routine data by extracting the nodes and units from the ANSYS model and inserting them into the RFPA 3D code. Figure 13 shows the correlation between the MS energy density of the Dagangshan right bank slope under excavation unloading obtained from the MS system and numerical simulation results and the AE distribution calculated by RFPA 3D . e MS energy density can be used to characterize the location and strength of MS events in rock mass, which can be used to delineate potential danger zones in the bank slope. e formation of potential slip surfaces during the slope excavation reveals that the slope will fail by overall buckling and confirms that the weak structural surface is one of the main controlling factors of the rock slope failure mode, as shown in Figure 13(a). e AEs during excavation process based on calculations from RFPA 3D are mostly distributed along the unloading fissure zone XL-316 and the fault f 231 , forming a potential slip surface as shown in Figure 13(b). Weakening of the dominant structure plane from the slope top downward during excavation process is the main cause of potential slope instability. e AE accumulation zone obtained from numerical solution and the distribution of high MS density zone are highly coherent, which showed that the evolution simulation of Dagangshan right bank slope stability under excavation using RFPA 3D did match the site monitoring results and established the model foundation of RFPA 3D slope stability analysis with MS data feedback.

Input and Analysis of MS Parameters.
e MS parameters from the monitoring system are fed into the above numerical model by a RFPA-MMS program. e actual spatial coordinates of the MS events are used to couple the model. According to the MS damage, the rock mechanical parameters are degraded to even the distribution of the seismic energy within the grid unit. Consequently, the seismic energy of each damaged rock element is characterized in space sphere units after the automatic allocation ( Figure 14). Automatic searching and positioning of MS damage was performed, and the mechanical parameters of the rock elements degraded correspondingly for the numerical model. e revised mechanical parameter of the numerical model rock unit (damaged elasticity modulus E r ) is reassigned based on the MS damage effect; thus, the degradation of mechanical parameters of rock slope and the stability analysis with the consideration of MS damage feedback were carried out.

Feedback Analysis Considering MS Data.
e computational progressive failure process in consideration of MS damage is shown as Figure 15. Upon the completion of bank slope excavation, the computational bank slope which considered the MS damage had its first damage in the uppermiddle zone along the unloading fissure zone XL-316, as shown in Figures 15(a) Finally, macrodamage appears on the side slope until collapse occurs and deformation along the valley direction occurs at the toe of the excavation slope. e subsidence towards the bottom of the rock slope takes place, and macroscopic failure begins and continues gradually in Step 7 (Figure 15(g)). e upper slope mainly shows tensile failure, and the base of the slope mainly shows shear failure with slight tensile failure. e blue spheres represent tensile failure, and the red spheres show shear failure (Figure 15(h)). e instability mechanisms of the rock slope are predominantly the initiation, propagation, evolution, and failure of the potential fracture surface along the main faults or weak zones. e calculated safety factor of the rock slope is 1.03. e safety factor of the rock slope is reduced by 0.18 (from 1.21 to 1.03) compared with the case without considering the MS-determined damage [34]. erefore, the effect of rock mass damage induced by excavation is significant in the stability analysis of the rock slope. Nevertheless, the safety factor still satisfies the requirements of hydraulic design specification in China. e MS monitoring and 3D simulation results indicate that the main faults play a significant role in rock slope instability. e initiation, propagation, evolution, and aggregation of microcracking inside the rock slope can be identified and delineated in real time by such an integrated approach. e rock mass parameters are updated in real time in the 3D FEM model considering the effect of MS damage. e numerical calculation of the excavated rock slope is performed through the RFPA 3D centrifugal loading method. e simulated results include the real-time safety factors considering the effect of MS damage and reproduce the 3D progressive failure processes of the rock slope. e results reveal the damage evolution mechanism of rock slope instability and visually reflect the potential sliding surface of the slope determining the rock slope failure mode.

Conclusions
MS monitoring provides an approach for quantifying rock mass degradation. A 3D FEM model was established based on the characteristics of MS data recorded at the right bank slope of the Dagangshan hydropower station. According to the proposed rock mass degradation criterion, the mechanical parameters of the rock mass were modified as inputs to update the properties of elements in the numerical model. e stability of the excavated rock slope was evaluated considering the MS damage effect, which is the contribution to the theory of this study. e following conclusions can be drawn: (1) e spatial distribution of MS delineated the weak geological structures and identified the potential sliding surfaces inside the rock slope. ere may be a risk of deformation on the right bank slope. e main cause of fissures in the rock slope surface was local small block deformation due to the existence of unloading fracture zone XL-316 and fault f 231 . (2) e 3D numerical simulation of the failure processes of the right bank slope revealed that the deep faults or weak structural planes in the rock slope played a significant role in the instability mechanism. e effects of rock mass damage and energy dissipation can be used to describe the constitutive equations and strength criterion of a rock mass. Based on the results of blasting testing events, the equivalent seismic efficiency of MS on the right bank slope of the Dagangshan hydropower station was calculated as 0.003%. e safety factor of the rock slope obtained by the RFPA 3D -centrifuge method was 1.03, which was 0.18 lower than the safety factor obtained without account for MS damage. is result is a direct quantitative verification of the rock damage caused by MS events in the process of rock slope excavation. (3) According to the numerical modelling and MS monitoring results, more attention should be paid to the junction of unloading fracture zone XL-316 and fault f 231 for the convergence and deformation situation. A basis for the safety control assessment of the dam construction and water storage phases is provided.   Step 46

Advances in Civil Engineering
Data Availability e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.