Deformation Analysis of Large-Scale Rock Slopes Considering the E ﬀ ect of Microseismic Events

: To research the macroscopic deformation of rock microseismic damage, a high-precision microseismic monitoring system was established on the left bank slope of the Baihetan hydropower station in Southwestern China. Based on the microseismic monitoring and ﬁeld deformation data, the seismic source radius was applied to characterize the rock fracture scale. Numerical simulations introduced the rock micro-fracture information into the three-dimensional numerical model of the left bank slope and established the damage constitutive model. The unloading deformation process of the left bank abutment rock mass is described by numerical calculations. The feedback analysis method considering the e ﬀ ect of microseismic damage is preliminary exploratory research, which provides a new idea for the stability analysis of similar high rock slopes. quasi-dynamic source The MS data were read into the corresponding rock unit in the model, and the source radius was entered. The numerical calculation can obtain the fracture density and fracture damage variable of the corresponding rock unit and the strain increment in the rock unit under the fracture damage. The stress and strain values of this element in the damage status could be obtained from the FLAC3D calculation and circulation. After all the MS data during the selected period were input into the numerical model, the stress and strain ﬁeld in the damage status of the rock slope could be obtained. distribution of and


Introduction
Slope stability analysis is an old and complex subject. Slope stability analysis methods mainly include the engineering geological analysis, model testing, limit equilibrium, limit analysis and numerical simulation methods. A single evaluation cannot accurately assess the rock state of a slope. Many scholars combined different methods and considered different factors to analyze the surrounding rock stability. Conforti and Ietto [1] conducted an integrated method to research slope stability. Slope stability is affected by many factors. The slope security and stability show varying degrees of attenuation under external disturbances, such as rainfall, earthquakes, weathering and excavation. In particular, Ietto et al. [2] studied the effect of weathering characteristics on the stability of rock slopes. Li et al. [3] combined numerical simulation and on-site monitoring methods to evaluate the stability of a rock mass. Based on the gravity increasing method (GIM), Li et al. [4] established a micromechanical model for studying slope the stability and failure processes. Liu et al. [5] applied a combination of monitoring technology and construction dynamics to study the stability of surrounding rock deformation. Xu et al. [6] introduced microseismic monitoring technology to the stability analysis of the surrounding rock mass. Zhang et al. [7] conducted numerical simulations of the damaged area. The physical and mechanical properties of the internal weak structural plane tend to deteriorate due to these disturbances, which are the fundamental causes of slope instability. Therefore, the subsequent behavior of the slope is analyzed using parameters corresponding to the degree of rock failure; reasonable predictions are rendered via this analysis. However, the mechanical parameters of rock slopes are affected by factors such as various cracks in the internal rock mass structure, joints, weak interlayers and rock mass size effects. Based on laboratory measurements and field measurements alone, it is difficult to properly derive rock mass mechanical parameters under excavation disturbances. data. The stress-strain conditions of the rock in fracture and damage status were better simulated. The functions of the FLAC3D software were expanded, providing a new idea for the quantitative evaluation of rock mass damage during the excavation of rock slopes.

Damage Constitutive Model Considering the Rock Fracture Scale
MS data obtained by the monitoring systems contain abundant source fracture information. Based on extensive research on MS waveform data, scholars have obtained source damage characteristic parameters related to seismic energy accumulation and release, such as the seismic moment, seismic magnitude, seismic energy, source radius, apparent stress, stress drop and seismic deformation. Cracks, joints and faults in a rock mass are discontinuities that greatly decrease the mechanical properties of the rock mass. In the present study, the source radius is used to represent the rock fracture scale, thus establishing a direct relationship between the rock fracture scale and rock failure. By considering the size and number of fractures, it is possible to determine an efficient method to estimate the compliance matrix of fractured rock mass based on micro-mechanics.
The constitutive relationship of the fractured rock mass mainly contains the intact rock strain and the fracture induced strain [22]. To study the fractured rock mass more conveniently, it can be divided into many representative volume elements (RVEs) as study objects [22]. Horii and Nemat-Nasser [23] proposed an overall estimate method for the instantaneous modulus based on the loading conditions when studying a linear, elastic and brittle solid with possible fractional sliding.
In two dimensions, the flexibility matrix of the rock mass with randomly distributed joints can be written in Equation (2) [24]: where E is the elastic modulus of the intact rock mass; ν is the Poisson s ratio; a 0 is the damage variable of the rock mass fractures and can be expressed as follows: where ω is the fracture density of the rock mass and can be expressed as follows: where N is the number of micro-fractures in the unit volume; V is the volume of the unit; r 0 is the MS source radius. In elastic mechanics, the displacement in the length direction can be considered equal to 0 if the length of the subject is much larger than the length in the other two directions, and the size and shape of each part are largely the same. This kind of displacement represents the plane strain problem. In planar problems, the stress, strain and displacement are functions that are independent of the length direction. As the profile of the left bank slope model is selected, it can be approximately simplified to a plane strain problem.
Based on Equation (2), Equation (5) is transformed to obtain the damage constitutive elastic matrix of the fractured rock mass under three-dimensional plane strain as follows:

Project Overview
The Baihetan hydropower station, located on the convergence of Ningnan County in Sichuan Province and Qiaojia County in the Yunan Province of China, is a large project recently constructed on the Jinshajiang River. Figure 1 shows the location of the hydropower station. The project has a double-curvature arch dam with a height of 289 m and a total installed capacity of 16,000 MW. The underground powerhouse's cavern systems are located in both the left and right bank mountains, 350 m downstream from the dam, with horizontal and vertical depths of 110-380 m and 160-420 m, respectively [25,26]. The left bank is the southeast slope of the Daliang Mountain region, with a peak height of 2600 m, and it has a sloped terrain inclined towards the Jinshajiang River. The right bank is the west slope of the Yao Mountain, with a peak height of 3000 m, and its landform is characterized by steep and gentle slopes. Figure 2 shows the high slopes of the spandrel groove after excavation. Hydro China Huadong Engineering Corporation is responsible for slope excavation design work. Through the experimental data of Hydro China Huadong Engineering Corporation, we produced Tables 1-4. We have considered the main geological conditions that are conducive to the numerical simulation analysis. Figure 3 shows a typical section of the left bank slope along the dam-arch axis. The left bank abutment has an approximate north to south direction, inclining to the east with a slope of 60 • . The stratum lithology of the left dam foundation is mainly composed of P 2 β 3 -P 2 β 4 laminar basalt. The first type of basalt with columnar joints is mainly distributed at an elevation of 700-600 m-generally with a length of 3 m and a diameter of 13-25 cm. The left bank slope is a bedding slope with a shape of N30 • -50 • E and SE∠15-25 • . There are many interlayer staggered zones (i.e., the appearance of C3 is N40 • -55 • E, the appearance of C3-1 is SE, 15 • -20 • ) and the faults with high deep angles are in the NE-direction (i.e., F17, F16 and F14). The characteristics of the main interlayer rock mass and faults are listed in Tables 1 and 2, respectively. According to in situ stress measurements, the maximum principal stress is σ 1 = 8-11 MPa, the medium principal stress is σ 2 = 7-9 MPa, and the minimum principal stress is σ 3 = 6-8 MPa [25][26][27].

Project Overview
The Baihetan hydropower station, located on the convergence of Ningnan County in Sichuan Province and Qiaojia County in the Yunan Province of China, is a large project recently constructed on the Jinshajiang River. Figure 1 shows the location of the hydropower station. The project has a double-curvature arch dam with a height of 289 m and a total installed capacity of 16,000 MW. The underground powerhouse's cavern systems are located in both the left and right bank mountains, 350 m downstream from the dam, with horizontal and vertical depths of 110-380 m and 160-420 m, respectively [25,26]. The left bank is the southeast slope of the Daliang Mountain region, with a peak height of 2600 m, and it has a sloped terrain inclined towards the Jinshajiang River. The right bank is the west slope of the Yao Mountain, with a peak height of 3000 m, and its landform is characterized by steep and gentle slopes. Figure 2 shows the high slopes of the spandrel groove after excavation. Hydro China Huadong Engineering Corporation is responsible for slope excavation design work. Through the experimental data of Hydro China Huadong Engineering Corporation, we produced Tables 1-4. We have considered the main geological conditions that are conducive to the numerical simulation analysis. Figure 3 shows a typical section of the left bank slope along the dam-arch axis. The left bank abutment has an approximate north to south direction, inclining to the east with a slope of 60°. The stratum lithology of the left dam foundation is mainly composed of P2β3-P2β4 laminar basalt. The first type of basalt with columnar joints is mainly distributed at an elevation of 700-600 m-generally with a length of 3 m and a diameter of 13-25 cm. The left bank slope is a bedding slope with a shape of N30°-50° E and SE∠15-25°. There are many interlayer staggered zones (i.e., the appearance of C3 is N40°-55° E, the appearance of C3-1 is SE, 15°-20°) and the faults with high deep angles are in the NE-direction (i.e., F17, F16 and F14). The characteristics of the main interlayer rock mass and faults are listed in Tables 1 and 2, respectively. According to in situ stress measurements, the maximum principal stress is σ1 = 8-11 MPa, the medium principal stress is σ2 = 7-9 MPa, and the minimum principal stress is σ3 = 6-8 MPa [25][26][27].        Table 2. Characteristics of the faults.

F14
Downstream of the arch dam Table 3. Intraformational disturbed belts.

No. Occurrence Width (cm) Condition
The thickness is 10-40 cm, and the occurrence varies greatly Breccia and cuttings are the main ones, and the thickness varies greatly The breccia and the cuttings are mainly distributed, and the top and bottom are distributed with a small amount of mud The breccia and the cuttings are mainly distributed, and the top and bottom are distributed with a small amount of mud

Failure Characteristics of the Left Bank Slope after Excavation
The As shown in Figure 4c, the expansion width is approximately 1 cm at the top and approximately 0.3 cm at the toe. Therefore, the slope excavation was suspended immediately on 29 December 2014. Related supporting measures were completed on February 2015 [28]. After approximately one year, the excavation of the slope restarted on 11 February 2016. Some similar sprayed concrete cracked at the peripheries of the entrances to WML1 at an elevation of 605 m on 25 May 2016, as shown in Figure  4d.

Model Setup
The geological profile of the slope along the arch axis at the elevations ranging from 550 m to 880 m was selected to build the numerical model, as shown in Figure 5. The software ANSYS was used to generate grids. The grid node data were imported to FLAC3D via a plug-in. The model

Model Setup
The geological profile of the slope along the arch axis at the elevations ranging from 550 m to 880 m was selected to build the numerical model, as shown in Figure 5. The software ANSYS was used to generate grids. The grid node data were imported to FLAC3D via a plug-in. The model consisted of 246,075 hexahedron elements. The numerical model had a bottom elevation of 480 m and a maximum elevation of 880 m. The model was 570 m long, perpendicular to the river, and 50 m along the river. The model mainly includes nine weak structural planes, as shown in Figure 3 (i.e., the intra formational disturbed belts LS337, LS331, LS3318, and LS3319, interlayer disturbed belts C3-1 and C3, and faults on the posterior border F17, f110, and f108). Table 3 shows the characteristics of the intra formational disturbed belts. The positive direction of the x-axis points to the right bank slope (SE 50 • ); the positive direction of the y-axis points to the river upstream (SW 40 • ); the positive direction of the z-axis is vertically downward. The x-direction normal constraints are imposed on both sides of the calculation area; the y-direction normal constraints are imposed on its front and back; the x-, y-, and z-direction constraints act on some of the key nodes at the bottom, and the free faces are free. The rock mass in the slope is classified into four categories according to the GB50287-2006 Code for Hydropower Engineering Geological Investigation in China: Class II, class III, class IV and class V, with the lithology of basalt. Moreover, the classes II, III, IV and V correspond to the values of 90-75%, 75-50%, 50-25% and below 25%, of the RQD criterion (rock quality designation) respectively [29].
Based on the experimental results listed in Table 4, the mechanical parameters of the rock mass and structural plane in the model were set. After applying the boundary conditions to the model, the C language was applied to compile the damage constitutive model, considering the effect of the rock fracture scale. When FLAC3D software runs, the complied information is directly loaded and attached to all rock elements in the model as a dynamic link library. The excavation of the spandrel arch of the foundation slope began on 5 September 2013. The layered and segmental blasting excavation method was adopted, with each layer approximately 10 m thick. The excavation progress is shown in Figure 6. To reflect the surrounding rock mass deformation response mechanism of the slope under layered excavation, the reinforcement and excavation were considered as truly as possible in the numerical simulation. A numerical model considering anchor cable reinforcement is shown in Figure 7. The anchor cables were adopted on-site with their prestressing force of 3000 kN. The length of the anchor cable was 30 m to 70 m (the length depends on the geological exploration conditions) and the arrangement spacing was 4 m × 4 m. The anchoring angle was perpendicular to the slope foundation surface. Figure 8 illustrates the anchoring diagram in the numerical model. and class V, with the lithology of basalt. Moreover, the classes II, III, IV and V correspond to the values of 90%-75%, 75%-50%, 50%-25% and below 25%, of the RQD criterion (rock quality designation) respectively [29]. The breccia and the cuttings are mainly distributed, and the top and bottom are distributed with a small amount of mud Based on the experimental results listed in Table 4, the mechanical parameters of the rock mass and structural plane in the model were set. After applying the boundary conditions to the model, the C language was applied to compile the damage constitutive model, considering the effect of the rock fracture scale. When FLAC3D software runs, the complied information is directly loaded and attached to all rock elements in the model as a dynamic link library. The excavation of the spandrel arch of the foundation slope began on September 5 2013. The layered and segmental blasting excavation method was adopted, with each layer approximately 10 m thick. The excavation progress is shown in Figure 6. To reflect the surrounding rock mass deformation response mechanism of the slope under layered excavation, the reinforcement and excavation were considered as truly as possible in the numerical simulation. A numerical model considering anchor cable reinforcement is shown in Figure 7. The anchor cables were adopted on-site with their prestressing force of 3000 kN. The length of the anchor cable was 30 m to 70 m (the length depends on the geological exploration conditions) and the arrangement spacing was 4 m × 4 m. The anchoring angle was perpendicular to the slope foundation surface. Figure 8 illustrates the anchoring diagram in the numerical model.

MS Parameter Embedding
The MS monitoring system produced by the Engineering Seismology Group, Canada was adopted in the project. The system was installed completely on 10 November 2014 [27,30]. The multipoint extensometers were successfully installed and commissioned on 14 December 2014. The slope was excavated near the 680 m elevation. Therefore, there is a lack of comparison between MS data above 680 m and field data, which has a certain impact on the simulation results. To build the quantitative relationship between rock mass damage and MS activities of the rock slope, the spandrel arch area of the slope was chosen as the research object to analyze the MS events during the selected period. MS events in two periods were selected. The first period was from 1  172.7, 60, 405) was selected. Based on the preceding coordinates, the coordinate transformation relationship was established to input the MS data in the selected area to the numerical model. During the excavation period, there were a total of 309 MS events. The damage constitutive model considering the rock fracture scale was adopted in all rock elements of the model. It was assumed that the MS monitoring system adopted an annular fault plane. The source fracture zone was equivalent to a disk with a radius of r, which referred to the source fracture scale according to Bruner's model calculation [31]. The formula is as follows: where K c refers to the angular frequency of the P-wave or S-wave; K s = 2.34 for instantaneous stress relief; K p = 2.01 and K s = 1.32 for a quasi-dynamic rounded source [32]. The MS data were read into the corresponding rock unit in the model, and the source radius was entered. The numerical calculation can obtain the fracture density and fracture damage variable of the corresponding rock unit and the strain increment in the rock unit under the fracture damage. The stress and strain values of this element in the damage status could be obtained from the FLAC3D calculation and circulation. After all the MS data during the selected period were input into the numerical model, the stress and strain field in the damage status of the rock slope could be obtained. Figures 9 and 10 show the spatial distribution of MS events and the relationship between MS information and the model excavation process. The MS monitoring system produced by the Engineering Seismology Group, Canada was adopted in the project. The system was installed completely on 10 November 2014 [27,30]. The multipoint extensometers were successfully installed and commissioned on 14 December 2014. The slope was excavated near the 680 m elevation. Therefore, there is a lack of comparison between MS data above 680 m and field data, which has a certain impact on the simulation results. To build the quantitative relationship between rock mass damage and MS activities of the rock slope, the spandrel arch area of the slope was chosen as the research object to analyze the MS events during the selected period. MS events in two periods were selected. The first period was from 1 Based on the preceding coordinates, the coordinate transformation relationship was established to input the MS data in the selected area to the numerical model. During the excavation period, there were a total of 309 MS events. The damage constitutive model considering the rock fracture scale was adopted in all rock elements of the model. It was assumed that the MS monitoring system adopted an annular fault plane. The source fracture zone was equivalent to a disk with a radius of r, which referred to the source fracture scale according to Bruner's model calculation [31]. The formula is as follows: where refers to the angular frequency of the P-wave or S-wave; = 2.34 for instantaneous stress relief; = 2.01 and = 1.32 for a quasi-dynamic rounded source [32]. The MS data were read into the corresponding rock unit in the model, and the source radius was entered. The numerical calculation can obtain the fracture density and fracture damage variable of the corresponding rock unit and the strain increment in the rock unit under the fracture damage. The stress and strain values of this element in the damage status could be obtained from the FLAC3D calculation and circulation. After all the MS data during the selected period were input into the numerical model, the stress and strain field in the damage status of the rock slope could be obtained. Figures 9 and 10 show the spatial distribution of MS events and the relationship between MS information and the model excavation process.    Figure 11 demonstrates the displacement nephogram of the rock slope during the excavation of the spandrel arch slope considering the MS event damage effect. Figure 12 shows the displacement vector of the rock slope at excavation to 550 m. The △Mzjc-3 monitoring point of the multipoint extensometer was chosen as the displacement monitoring point of the numerical simulation. Figure  13 indicates the relationship between the horizontal and total displacement of the displacement monitoring point and the excavation process. (1) The excavation-caused slope deformation was composed of the vertical resilience displacement and horizontal displacement. The slope deformed upwards to the valley. However, the degree of deformation of the slope surface was weakened due to the additional anchoring effect. The maximum deformation of the slope was less than 30 mm. (2) The maximum displacement of each layer after excavation occurred near the excavation face, especially at the foot of the weak structural plane. When the excavation reached the 720 m elevation, the largest displacement appeared in the C3 and C3-1 exposed areas; when it reached the 680 m elevation, the largest displacement appeared in the F17 exposed area; and when it reached the 650 m elevation, the largest displacement appeared in the LS3319 exposed area. (3) The influence of the horizontal unloading effect on the slope increased constantly with the lower excavation plane, thus increasing the proportion of the horizontal displacement in the total. When excavation reached the 720 m elevation in the middle and upper part, the horizontal displacement accounted for approximately 21% of the total displacement; when excavation reached the 550 m elevation in the lower part, the horizontal displacement accounted for approximately 41% of the total displacement.  Figure 11 demonstrates the displacement nephogram of the rock slope during the excavation of the spandrel arch slope considering the MS event damage effect. Figure 12 shows the displacement vector of the rock slope at excavation to 550 m. The Mzjc-3 monitoring point of the multipoint extensometer was chosen as the displacement monitoring point of the numerical simulation. Figure 13 indicates the relationship between the horizontal and total displacement of the displacement monitoring point and the excavation process. (1) The excavation-caused slope deformation was composed of the vertical resilience displacement and horizontal displacement. The slope deformed upwards to the valley. However, the degree of deformation of the slope surface was weakened due to the additional anchoring effect. The maximum deformation of the slope was less than 30 mm. (2) The maximum displacement of each layer after excavation occurred near the excavation face, especially at the foot of the weak structural plane. When the excavation reached the 720 m elevation, the largest displacement appeared in the C3 and C3-1 exposed areas; when it reached the 680 m elevation, the largest displacement appeared in the F17 exposed area; and when it reached the 650 m elevation, the largest displacement appeared in the LS3319 exposed area. (3) The influence of the horizontal unloading effect on the slope increased constantly with the lower excavation plane, thus increasing the proportion of the horizontal displacement in the total. When excavation reached the 720 m elevation in the middle and upper part, the horizontal displacement accounted for approximately 21% of the total displacement; when excavation reached the 550 m elevation in the lower part, the horizontal displacement accounted for approximately 41% of the total displacement. Figure 14 shows the shear strain rate of the rock slope considering the MS damage effect during the excavation process of the spandrel arch on the left bank slope. The red area indicates the largest shear strain rate. After the slope excavation began, the stress concentrated near the weak structural plane due to excavation disturbance, resulting in a large shear strain rate. Shear strain occurred and became more obvious at a lower excavation elevation. The unloading depth of the surrounding rock mass also increased. When excavation reached an elevation of 550 m, the shear strain rate increased dramatically in the wedge area formed by F17 and LS3319, indicating that excavation in the lower part of the slope will affect the deformation on the upper structural plane. The large shear deformation occurred in the interlayer disturbance zone LS3319, mainly due to the presence of large ground stress and weak layers on the upper disc wall. The increase in horizontal unloading after excavation expands the horizontal displacement, resulting in tensile shear deformation of the surrounding rock. Moreover, the shallow surrounding rock mass in this area was columnar-jointed basalt with poor resistance to excavation disturbance. Therefore, to avoid dangerous slip blocks formed by slip surface cutting, the columnar jointed basalt at altitudes of 660 m to 628 m should pay attention to the excavation process and support over time.          ground stress and weak layers on the upper disc wall. The increase in horizontal unloading after excavation expands the horizontal displacement, resulting in tensile shear deformation of the surrounding rock. Moreover, the shallow surrounding rock mass in this area was columnar-jointed basalt with poor resistance to excavation disturbance. Therefore, to avoid dangerous slip blocks formed by slip surface cutting, the columnar jointed basalt at altitudes of 660 m to 628 m should pay attention to the excavation process and support over time.  Figure 15 shows the deformation comparison of rock slopes under different approaches during excavation process. The △Mzjc-3 monitoring point of the multipoint extensometer was selected to record the excavation processes data, both on the feedback analysis model and the model without importing the MS data. Excavation data at elevations from 650 m to 550 m during the MS monitoring period between 1 March 2016 and 31 December 2016, and the corresponding results of the multipoint extensometer, were selected. It was found that the displacement of the numerical simulation considering MS damage was slightly less than the displacement of the field monitoring. This finding might be attributed without considering the effect of the columnar basalt in the numerical simulation. Compared with the model without considering MS damage, the displacement of the slope considering the MS damage effect was closer to the field monitoring data. However, the increment was not obvious, being only approximately several millimeters. Some interpretations might be that the selected research period was short and that the MS events in the selected period had a small source radius and were distributed. The monitoring displacement of the field multipoint extensometer in the early stage grew faster than that of numerical simulation. This result was mainly because the multipoint extensometer was installed on 14 December 2014 and the displacement during the early stage of installation would clearly grow.  Figure 15 shows the deformation comparison of rock slopes under different approaches during excavation process. The Mzjc-3 monitoring point of the multipoint extensometer was selected to record the excavation processes data, both on the feedback analysis model and the model without importing the MS data. Excavation data at elevations from 650 m to 550 m during the MS monitoring period between 1 March 2016 and 31 December 2016, and the corresponding results of the multipoint extensometer, were selected. It was found that the displacement of the numerical simulation considering MS damage was slightly less than the displacement of the field monitoring. This finding might be attributed without considering the effect of the columnar basalt in the numerical simulation. Compared with the model without considering MS damage, the displacement of the slope considering the MS damage effect was closer to the field monitoring data. However, the increment was not obvious, being only approximately several millimeters. Some interpretations might be that the selected research period was short and that the MS events in the selected period had a small source radius and were distributed. The monitoring displacement of the field multipoint extensometer in the early stage grew faster than that of numerical simulation. This result was mainly because the multipoint extensometer was installed on 14 December 2014 and the displacement during the early stage of installation would clearly grow.

Conclusions
Based on the damage-constitutive relationship considering the scale of rock failure, the deformation feedback analysis of the surrounding rock of the left bank slope of Baihetan hydropower station was carried out. The following conclusions can be drawn. Deformation comparison of the rock slope using different approaches during excavation processes.

Conclusions
Based on the damage-constitutive relationship considering the scale of rock failure, the deformation feedback analysis of the surrounding rock of the left bank slope of Baihetan hydropower station was carried out. The following conclusions can be drawn.
Compared with the traditional constitutive relationship, the extended damage constitutive relationship, considering the rock fracture scale can directly reflect the damage-effect of the microcrack within the rock. The damage constitutive model can select the characterization parameters of the rock fracture scale according to the calculation requirements.
The proposed method integrated MS data with the FLAC3D simulation. In particular, this method applied the source radius to characterize the rock mass fracture scale and analyze the surrounding rock deformation characteristics. Compared to the traditional method, this method can better invert the mechanical properties of the rock mass under construction disturbance, and realize the simulation analysis of the slope stability, considering the effect of MS damage. This method provided a new idea for the quantitative evaluation of surrounding rock failure and the analysis of surrounding rock feedback.
The numerical results indicated that the displacement of the model, considering the impact of MS damage, is closer to the field monitoring data than the model without considering the influence of MS damage. The increment was not obvious; i.e., only a few millimeters. This finding was mainly due to various reasons and limitations. The selected research period was short; the MS event source radius was small; and the distribution was scattered. Moreover, the slope deformation in the simulation considering MS damage was slightly less than that of field monitoring. This reduction may be because the numerical simulation did not consider the influence of the columnar basalt with the obvious deformation characteristics in the area. Therefore, the feedback analysis of rock mass mechanical properties is critical for the stability evaluation of high rock slopes under complex conditions.