DynamicResponseCharacteristics andFailureMechanismofCoal Slopes with Weak Intercalated Layers under Blasting Loads

Rock slopes with weak intercalated layers could experience disturbance from various deep mining activities; however, their dynamic stability has not been thoroughly investigated. In this paper, the dynamic response characteristics and failure mechanism of the coal slopes with weak intercalated layers under blasting loads were studied by means of numerical analysis, shaking table tests, and field tests. /e effects of dynamic loads with different frequencies on the dynamic response of the slope were analyzed, and the natural frequency of the slope was also determined. /e results show that the dynamic amplification effect of the slope is smaller than that of the homogeneous slope, and weak layers weaken the wave propagation in the rock mass. Both experimental and field investigation results show that the slope’s natural frequency was approximately 35Hz. /e slope deformation decreased with the distance of the blasting source. Cracks appear along the weak interlayer firstly under the action of horizontal vibration; then, longitudinal cracks occur at the slope crest. With the increase of dynamic loads, cracks continue expanding, deepening, and penetrating in the main controlled weak interlayer; then, the sliding body presents tensile shear failure along the sliding surface. /is study could provide insights into the understanding of the coal slope instability and failure mechanism; this could benefit the blasting operation of the coal slope in fields.


Introduction
According to the International Energy Agency (IEA, 2019), open-pit mining has made a significant contribution to the industry, producing around 17 billion tonnes of mineral fuels and industrial minerals annually [1]. With the sustainable demand for mineral resources and increasing consumption of shallow resources, most open-pit mines reach their designed ultimate depth. Under such a circumstance, deeper open-pit mining is performed on the higher and steeper slope [2][3][4][5][6], and good slope stability of open-pit mining is required to achieve its safety and financial targets [7,8]. Also, underground mining is combined with open-pit mining to produce more mineral resources.
Rock slopes could experience the disturbance from various deep mining activities in the process of combined open-pit mining and underground mining. In such deep mining process, the high-stress hard rock experiences various extreme events such as blasting operations, mechanical shock earthquakes, excavation unloading, and other dynamic disturbances around the deep underground rock [9][10][11][12][13][14][15][16], causing the significant change of the stress field and displacement field in the slopes and thus worsening the instability of slopes and landslide [17][18][19][20]. For example, on May 13, 2011, the north side of Shengli open-pit mine in Inner Mongolia, China, experienced a catastrophic landslide, which had a length of about 1.3 km, an area of 0.886 km 2 , and a volume of approximately 85 million m 3 [21]. Note that there also have been four small-and medium-scale landslides reported in this area during the past five years. Considerable economic and colliery losses have been caused. rough field investigation and disaster exploration, we found that the landslide could be triggered by the initiation and development of joints and structural planes in the slope and the subsequent blasting dynamic disturbance. Special attention should be thus paid to understanding the dynamic stability of coal slopes subjected to blasting dynamic disturbance. e rock mass of the deep concave slope in the open-pit and underground combined mining is usually composed of multilayer soft and hard interlayer structural planes, of which the mechanical properties influence the slope deformation characteristics [19,21,22]. A number of studies have shown that structural planes are the main controlling factor in controlling coal slopes' local instability. For instance, the dynamic stability of coal slopes could be significantly affected by the mechanical behavior, structural form, and spatial distribution of structural planes [23,24]. e shear failure of structural planes of coal slopes has been reported to cause severe slope instability problems [18]. Weak structural surfaces could be formed by multitudes of joints, and rock mass integrity could be impaired if they are widely distributed throughout rock slopes [25][26][27]. Given the uncertain geometric and mechanical parameters between them, such as uncertain directions, dip angles, and locations, the distribution of joints could be complex and difficult to understand [26,27]. erefore, under blasting dynamic disturbance, the dynamic mechanical properties of coal slopes with weak intercalated layers could be rather complicated.
e dynamic response characteristics and stability of rock slopes, particularly those disturbed by excavation blasting, have been studied by several researchers. e rock stress measurements were performed by using techniques such as the compact conical-ended borehole overcoring (CCBO) technique, developed by Obara et al. [28]. eir measurement and 3D finite element analysis (FEA) results showed that the horizontal stress component is important in estimating the stability of rock slopes, even the mine is located near the top of a mountain. As the joints in rock could weaken structural surfaces, practical and realistic simulations of rock mass stability should include the influence of joints [29]. For soft-rock slopes, their instability and failure are mainly caused by the expansibility, dispersivity, extremely low shear strength levels, softening behavior, and preferential flows of soft-rock materials [30]. Given the stochastic behavior of the geomechanical properties (e.g., cohesion, friction coefficient, and unit weight) of the slope rock masses as well as the hydrogeological medium properties (e.g., hydraulic conductivity), probabilistic analysis has been applied to rock slope stability [31][32][33]. e risk models could be used to optimize design variables such as the slope angle through Monte Carlo simulation to fulfill the objective function that is formulated regarding the excavation cost [33].
To further understand the real dynamic response of rock slopes, shaking table tests have been used in the lab to reveal the variations of the dynamic response of rock slopes and their failure modes under seismic excitation [34][35][36]. Various types of lithology combinations and structures, slopes with discontinuous joints [37], and the slope with weak intercalation [37,38] have been included in the lab test. In large-scale shaking table tests, the model size, the ground motion input parameters, and waveforms can be controlled based on the case studies to enable the slope dynamic responses to be measured to greater accuracy [39][40][41][42][43][44]. e testing results generally show that the slope angle, elevation, and surface microtopography have large impacts on the slope stability. e damage evolution process of the slope during earthquakes includes stages from the formation of tensile cracks in the strata and their propagation to the crest of the slope and sliding failure [45,46]. However, the above review shows that few studies focus on the understanding of the dynamic response characteristics of slopes under blasting loads; in particular, little research has been found on the analysis of the dynamic response of coal slopes with weak intercalated layers due to blasting loads.
In this paper, numerical simulations, shaking table tests, and field tests were performed to explore the dynamic stability of a coal rock slope with weak bedding layers under blasting loads. A case study on an actual engineering project was carried out. A mechanical slope model was first established based on LS-DYNA software platform to simulate the dynamic deformation characteristics of the slope with weak layers under blasting loads. e shaking table test was performed to investigate the dynamic failure evolution process of the slope under blasting dynamic loads in terms of the evolution rules of horizontal displacement and shear strain nephogram. Moreover, field tests and field investigations were used to verify the dynamic failure mechanism and deformation characteristics of the slope obtained from the numerical and shaking table test results. Our research aims to provide insights into the understanding of the dynamic response of coal slopes with weak intercalated layers due to blasting loads.

Case Study
As a typical coal slope with weak intercalated layers, Datang Shengli coalfield was selected for the case study. It is located in the middle and eastern part of the Inner Mongolia plateau, northern China (Figure 1), which belongs to the hilly plateau terrain with a height of 970.0-1,326.6 m. e geomorphology of the slope is formed by structural denudation, mainly consisting of four geomorphic units: structural denudation topography, denudation accumulation topography, erosion accumulation topography, and karst platform [47].
Shengli east No. 2 open-pit mine is located in the transition zone of denudation accumulation terrain and low hilly terrain in the east of Shengli coalfield, with an elevation of 990.0-1,122.0 m and a relative elevation difference of 143.0 m, shown as the "Study area" in Figure 1(a). e coalbearing strata in this area are lower coal-bearing member (K1bx12) of Xilin formation of Bayanhua group of lower cretaceous and upper coal-bearing member (K1bsh1) of Shengli formation. e coal seam of Shengli group is well developed, which is the primary coal-bearing section of this section. e geological condition of the deformation area is rather complicated, considering the fact that there are several weak mudstone intercalations with low strength from top to bottom (i.e., four weak intercalations with gentle dip and one weak intercalation with steep dip). Figure 2 shows the geomorphology of the slope. Table 1 shows the parameters of the corresponding rock samples. e mechanical properties of these rocks were obtained according to a series of laboratory tests, mainly including triaxial test, uniaxial compression test, and direct reduction test. Due to the rheology of the weak mudstone layer and the impact of blasting load, the shear strength of the coal rock slope is reduced. e weak mudstone layer is inclined, forming a natural potential slip surface, which is detrimental to the stability of the slope. In case of instability, the slope would slip along the weak mudstone layer. Hence, it is necessary to study the stability of the Datang south slope to better analyze the disaster mechanism of the slope with weak layers under the influence of dynamic loads. Table Tests. As the monitoring points on the slope are 200-400 m away from the blasting location in our study, the slope is mainly subjected to the earthquake waves. When the seismic wave generated by blasting acts on the slope, the vibration wave propagates in the form of longitudinal and transverse waves, but the peak value of the horizontal vibration displacement is much larger than that of the vertical vibration displacement; that is, the slope vibration is mainly horizontal under the earthquake waves. According to the above-mentioned blasting characteristics, the shaking table test applies a simple vibration in the horizontal direction to simulate the vibration wave. e peak amplitude is 10 mm, and the vibration frequency is between 10 and 40 Hz. e testing results are used to generally understand the dynamic response characteristics and failure mechanism of coal slopes with weak intercalated layers under blasting loads.

Testing Models and Measurement
System. e shaking table test was performed in the Institute of Mechanics, Liaoning Technical University, China. e bidirectional electric servo shaking table is 1.5 m × 1.0 m. Waves can be used as the input motions, and the effective frequency range is 0.5-50 Hz. A maximum horizontal acceleration of 1.5 g and a maximum vertical acceleration of 1.2 g can be obtained. To simulate the semi-infinite foundation of the model, the effect of the boundary should be minimized when the experimental tank is designed. A rigid, sealed experimental tank comprising organic glass and carbon steel plates was developed for the tests, and its inner dimensions are 1.5 m × 0.3 m × 0.453 m ( Figure 3). e slope is the prototype of the field slope, as shown in Figure 2. e materials and test parameters of the model were deduced based on the Buckingham π theorem of similarity [48]. e acceleration a, mass density ρ, and dimension L were selected as the fundamental parameters, with similarity ratios of C a � 1, C ρ � 1, and C L � 25, respectively. From the similarity calculation, the parameters of the rock mass simulation materials in the tests are listed in Table 1. e simulation materials were mainly formed from silica sand, calcium carbonate, and gypsum at different ratios, according to a series of triaxial tests. Cement, lime, kaolin, paraffin, and distilled water solution with retarder were selected as cementing materials. After stirring evenly, filling modeling was carried out according to the size of the model. e apparent density of similar materials is 1.5 g/cm 3 . Mica powder is used for layering between layers. e thickness of each loading is generally 5-10 mm.
In order to test the dynamic response of the slope when subjected to blasting vibrations with different frequencies, six identical models were prepared. To track the displacement of the slope subjected to seismic loads, the optical measurement technique analysis was employed. e displacement data were collected by using XTDIC measurement system, including a camera with two million pixels, and its acquisition frequency is 50 Hz, as shown in Figure 3. e camera was used to collect the displacement data in xand z-directions.
e highest acquisition frequency of Advances in Civil Engineering XTDIC measurement system is 340 Hz, and its measurement error within the range of 5 m is about 0.05-0.1 mm. e twodimensional displacement of every marked point can be obtained in the tests, and the layout scheme of the marked points is shown in Figure 3. e XTDIC measurement system is based on binocular stereo vision technology, using the high-speed camera to collect the real-time images in all stages of the deformation object. e accurate recognition of marker points or digital speckle was used to realize stereo matching, including coding and noncoding marker points, to calculate the deformation of the object. is method has characteristics of high precision, high speed, and easy operation and noncontact. Figure 3 illustrates the setup for measuring displacement. e cameras were stationary after calibration; a few states were collected as a benchmark when the model was static. When the shaking table began to shake, the cameras began to collect displacement data. Finally, the spatial coordinates of each state of the mark points were calculated to analyze displacements of the marker points; i.e., the displacement data were obtained based on the analysis of the change of coordinate values. Note that the data analysis of XTDIC measurement system supports coordinate transformation.
Blasting vibration is mainly the horizontal displacement produced by the longitudinal wave. us, the experiment mainly simulates the deformation and failure of the slope  e weak interlayer and the vicinity of the fault were slightly staggered with vibration, and debris fell with vibration. Figures 5 and 6 show that when the vibration frequency is 35 Hz, the cracks appear and penetrate rapidly at the fault outcrop and weak interlayer junction. At the same time, some small cracks were generated on the slope surface and the top of the slope, and the cracks gradually began to develop. With the increase of vibration time, the particle displacement at about 0.5 s was amplified due to resonance, and the slope was suddenly seriously damaged and collapsed. Figure 5(d) shows that the shear strain first appears in the weak structural plane inside the slope. With the continuous increase of vibration time, vertical cracks are formed at the top of the slope of each step, and the internal dislocation of the slope along the direction of shear belt is more and more obvious, resulting in the shear failure in other places in the direction of the shear belt that has not been damaged, and the shear belt looks more continuous. Finally, a shear outlet appears at the junction of the weak layer at the bottom of the slope and the fault, and at the same time, it is connected with the fractures in the vertical direction.
To further study the dynamic response characteristics of the slope, take a point at the top of the slope as an example, its displacement-time history curves under the action of vibration waves of different frequencies are shown in Figure 7. e results show that with the increase of the input frequency, the displacement of the measuring point also increased. However, when the frequency was 35 Hz, there was an anomaly and the displacement of the measuring point suddenly increased up to 18.9 mm, indicating that the shape of slope or the parameters of rock and soil material change greatly, and the slope is destroyed dynamically.
e above analysis shows that the dynamic damage evolution process of the slope can be summarized as follows: first of all, shear bands are formed in the weak layers due to horizontal vibration. Secondly, vertical tension cracks are generated at the top of the slope, and the two gradually develop with the action of vibration until they are broken through, forming a fracture plane composed of tension shear. It is not a traditional shear slip plane, nor a tension crack generated by acceleration amplification at the top of the slope, but a combination of the two. erefore, our shaking table test shows that when the slope model is subjected to vibration loads with different frequencies, the integrity is still good, although there are slight degrees of damage. When the vibration frequency is 35 Hz, the slope suddenly collapses; this could be caused by resonance. erefore, it can be determined that the natural frequency of the slope model is about 35 Hz.

Numerical Model.
Finite element method (FEM) dynamic analyses were then used to simulate the wave propagations through rock masses; thus, the internal, continuous stress development in rock masses can be characterized and understood. Due to the capabilities of transient dynamics simulation, the LS-DYNA program was adopted as a solver to conduct the FEM dynamic analyses.
is program, which mainly uses Lagrange, ALE, and Euler algorithms, is fully functional with geometric nonlinearity (e.g., large displacement, large rotation, and large strain), material nonlinearity, and contact nonlinearity. It not only provides an explicit solution but also generates implicit solutions concurrently. e rock and explosive elements were simulated by the ALE algorithm of the 4-node solid element PLANE62, the explosive element, and the exploded element share nodes.
is setting could have the fastest calculation speed, and it can continue the calculation even if the contact element has been broken. ere are two kinds of rock and soil parameters in this study: weak-layer visco-elastoplastic materials and other rock mass plastic materials. e constitutive model of weak-layer viscoelastic materials is the linear viscoelastic model proposed by Herrmann and Peterson [49].
For other rock mass plastic materials, * MAT_ PLASTIC_KINEMATIC material was selected as the simulated material, and the yield condition is where P and C are input constants and _ ε ij and ε   Advances in Civil Engineering 7 follows: when the explosion happens, the detonation wave front satisfies the following equation: where p D , ρ D , μ D , and C D are the pressure, density, particle velocity, and sound velocity on the detonation wave front, respectively; k, D, and ρ e are the polytropic exponent, explosive detonation rate, and explosive density, respectively. After the initiation of high-energy explosives, the pressure p of the inner element of the explosive was calculated by the JWL equation, which describes the change of chemical energy when the explosive is exploded. In order to simulate the dynamic failure mechanism of the slope with weak intercalated layers under blasting dynamic loads, one numerical model with the homogeneous slope (i.e., Model 1) and the other one with the slope having weak interlayers (i.e., Model 2, as shown in Figure 8 , in section shaking table tests,  blasting waves with a series of vibration frequencies including 15 Hz, 20 Hz, 25 Hz, 30 Hz, 35 Hz, 40 Hz, and 45 Hz were used to test the dynamic response characteristics of the slope in our experiment. e testing results showed that, when the frequency was 35 Hz, the slope collapsed due to the resonance (also the maximum horizontal displacement of the monitoring points of 18.29 mm occurred), indicating that the natural frequency of the slope model is about 35 Hz. Hence, the blasting wave with 35 Hz frequency was used in our numerical model to simulate the dynamic failure mechanism of the slope with weak intercalated layers under blasting dynamic loads.

Numerical Results.
e equivalent stress nephograms of homogeneous slope (Model 1) and slope with weak intercalated layers (Model 2) under blasting loads are, respectively, shown in Figures 9(a) and 9(b). Figure 9(a) shows that the detonation wave took the explosion source as the center and spread outwards in a spherical manner. e failure area of the explosion source gradually increased, forming the fracturing and crushing area, and the rock mass was completely destroyed. With the distance from the explosion source, the shock wave decayed into a stress wave and failed to crush the rock mass directly. However, the tangential tensile stress of the rock mass exceeded its tensile strength, and tensile failure occurred in the surrounding rock. e red region shows that the maximum tensile stress was first reached and the rupture occurred. With the increase of the sphere radius, the reflection wave was generated when it reached the boundary transits.
Compared to Model 1 results, the weak interlayer in Model 2 had an obvious barrier effect on the stress wave propagation when the explosion wave reached the weak interlayer. When the stress wave propagated from the hard rock to the weak interlayer, the wave resistance of the weak interlayer was less than that of the hard rock. When the stress wave reached the interface of the two media, part of the wave energy entered the soft medium, leading to the attenuation of the stress wave. When the stress wave propagated from the weak interlayer to the hard rock, a large amount of stress wave energy was reflected at the interface of the two media, and only a small number of waves penetrated the interface and entered the hard rock. is phenomenon shows that the weak interlayer not only effectively weakened the stress of the slope under blasting loads but also changed the mechanical properties of the weak interlayer, making the slope easier to produce creep deformation. e displacement-time history and velocity-time history of the three typical monitoring points (i.e., 2640, 6248, and 11707) of the two models in the x-and y-directions are shown in Figures 10 and 11. e changes of peak ground acceleration (PGA) and peak ground displacement (PGD) at the points are shown in Figure 12. As shown in Figure 10, the horizontal and vertical displacements of the monitoring points of the slope with weak interlayers under blasting loads   Figures 10-12 also show that the blasting wave amplitude was extremely large in the area close to the explosion source. As the explosion source distance increased, the maximum value of the vibration amplitude of the particle decreased and tended to be stable with the increase of the distance. From the results of Figures 10(a)-10(d) and Figures 11(a)-11(d), we obtained the PGA and PGD curves of two models, which are shown in Figure 12. Figure 11 shows that the horizontal vibration velocity amplitude of the particle at the same point was larger than the vertical vibration velocity,    and the horizontal PGA and PGD were approximately twice as much as the vertical vibration amplitude. Both the horizontal and vertical acceleration-time history curve decayed rapidly. After about 3,000 microseconds, the vibration amplitudes of the particle decayed to an extremely small level; thus, they are not shown in Figure 11. After that, the particle would make a slight vibration near the equilibrium position until the energy dissipation was complete. e particle vibration amplitude of Model 2 is significantly less than that of Model 1, which is due to the damping effect of weak interlayers on the propagation of the stress waves. Figure 13 shows the blasting vibration monitors that include the NUBOX-6016 intelligent vibration monitor and the supporting speed sensor produced by Sichuan Top Measurement and Control Technology Co., Ltd. ey are mainly used for recording signals such as blasting, vibration, and noise. Deep hole blasting was adopted in our tests. According to the vibration measurement data in the trial and error tests, the frequencies of blasting vibration were kept at 7-18 Hz; hence, when judging the impact of blasting vibration on the slope and newly poured mass concrete, the frequency should be <10 Hz or 10-50 Hz [49]. e displacement sensor is required to be rigid and fixed with the binding surface. As the open-pit mine is a soft-rock slope, the sensor area is small, and the mass is light, and it is difficult to form a rigid fixation with the contact surface of loose soft rock. In order to fix the sensor, the bedrock was first bonded with the cement column by using gypsum. e sensor was then bonded to the top surface of cement column by using gypsum.

Blasting Vibration Monitoring Test.
During the blasting operation of the slope, the displacement sensor was fixed on the side slope for data monitoring, as shown in Figure 13, to analyze the dynamic deformation characteristics of the slope under the action of blasting. e relative maximum displacements of the monitoring point at the top of the slope during the two blasting operations in May and June are shown in Figure 14.
e results show that, under the influence of blasting vibration, the cumulative maximum creep displacement of the slope monitoring point is approximately 57 mm and 48 mm. e displacement variations follow the typical development of three stages of creep: initial creep, stable creep, and accelerated creep. Additionally, to analyze the impact of blasting wave frequency on the dynamic response characteristics of the slope, the displacement development of a monitoring point subjected to different blasting loads with frequencies from 25 to 45 Hz is presented in Figure 15. e field test results show that the displacement of the monitoring points under the action of blasting waves with different frequencies is different at the early age of around 3 seconds. e maximum horizontal displacement of the monitoring points is 18.29 mm when the frequency of the main vibration is 35 Hz, indicating that the natural frequency of the slope is about 35 Hz.
us, the slope is most prone to instability failure under the action of this vibration frequency, which is consistent with the results of the numerical simulation and the shaking table test.

Deformation
Mechanism of the Landslide. Besides the above blasting vibration monitoring test performed in the field, the field investigation was also carried out to observe the deformation of the landslide and understand its failure mechanism.
Our numerical and shaking table test results show that the factors triggering landslides could include bad geological conditions and improper blasting loading. is could be verified by our field investigation. For example, there are many weak mudstone intercalated layers with low strength in the south side of the first mining area of the open-pit mine. e weak intercalated layers have great rheology to cause their deformation to increase with time. Also, the weak layer is inclined, forming a natural potential slip surface, which is detrimental to the stability of the slope. If this instability occurs, the slope could slip along the weak mudstone layer. Our comprehensive study also shows that the influence of blasting vibration on the slope stability of the open-pit mine is also an important external factor of slope failure and landslide. As an example, Figure 16 presents a slide shear outlet of the landslide that occurred in the platform at the elevation of 912-948 m on August 22, 2011. Different from traditional landslides, the landslide in this study presents several unique deformation and failure characteristics during the deformation evolution process, which are highlighted in the following three aspects. e slope deformation duration is long, and its deformation rate is large. Although a series of landslides have occurred since 2011, the slope has been still in a state of constant deformation. e monitoring points between fracture zones have obvious migration at the beginning, and the deformation time has been as long as half a year; particularly, the cumulative displacement can reach tens of meters. e maximum deformation rate is always maintained at more than  tens of millimeters per day, and the rate during the severe deformation period can reach several meters per day. For example, the maximum deformation rate at CL26 and CL112 monitoring points can reach 1.08 m per day. e landslide has a long time, high rate of deformation, and evolution process, which is not the traditional landslide, as shown in Figure 17. e failure mode of the landslide is the progressive and ulcerated sliding failure. e landslide began to damage near the F68 fault (Figures 18 and 19), and two arc-shaped fracture zones were successively formed, as shown in Figure 19. e failure area of the slope extends outwards step by step, presenting a gradual and ulcerated failure mode. According to the displacement monitoring curves (Figure 17), there are two groups of measuring points, including near fracture belt points (CL109, CL115, CL117, and CL135) and remote points (CL136, CL142, CL156, and CL104). e obvious difference between the displacement and velocity of the two groups of monitoring points can be found. erefore, although these two groups of points are in the deformation zone, they cannot slide along a certain slip surface as a whole, but they just present multistage ulceration sliding deformation. e failure mode of the landslide is very different from the traditional slide failure mode of traction or push-back. e deformation of the landslide is characterized by near-horizontal creep slip of the loose slide. During the deformation and failure process of the landslide, in the hanging wall area of the F68 fault below the elevation of 924 m, loose slides constantly appear on the slope and eventually form a large area of loose slides accumulation area (see Figures 18 and 20). e gradient of the loose sliding body is similar to the horizontal level, but the sliding body is still able to move northward continuously, causing the exposed coal area at the bottom of the pit to be covered by the sliding body, which is in a wavy shape, as shown in the elliptic region in Figures 20(a) and 20(b). In addition, in November 2011 and October 2012, the slope foot pressing measures were implemented in the flat plate area at an elevation of 936 m and 924 m, respectively, which not only failed to achieve east-west closure but also led to a large number of foot material migration to the depth of the slope (Figure 20). e movement of nearhorizontal loose slide and presser foot material is very different from the deformation and failure of the traditional landslide.

Discussion
e dynamic response characteristics of a coal rock slope in the open-pit mine were studied by numerical simulations, model tests, and field monitoring. e numerical simulation is used to study the wave propagation characteristics in the coal rock slope under the action of the blasting wave, and the influence of the weak intercalation on the wave propagation characteristics in the slope is also analyzed. It is difficult to analyze the wave characteristics of the slope in shaking table tests and field monitoring tests. However, the numerical simulation adopts the finite element method, which belongs to the study of dynamic response characteristics in the elastic domain of small deformation, and does not involve the study of the process of large deformation failure, deformation response characteristics, or failure modes. Hence, the shaking table test was used for further analysis. e purpose of the shaking table test is to study the deformation characteristics, failure process, and failure mode of the slope under the blasting wave and to explore the influence of blasting waves with different vibration frequencies on the dynamic deformation response of the slope. Moreover, the conclusion of the shaking table test is further verified by field monitoring.
Additionally, our comprehensive study shows that the natural frequency of the coal slope in this case study is 35   14 Advances in Civil Engineering Advances in Civil Engineering 15 estimated by monitoring the maximum slope deformation when the slope is under blasting vibration with a series of frequencies. is procedure could thus be potentially used in the field to identify the natural frequency of coal slopes. Compared to the slope with homogeneous materials, the slope with weak intercalated layers has significant dynamic behaviors because the weak layers weaken the blasting wave propagation in the rock mass. e different mechanical properties of weak interlayer and the bulk slope zone induce a huge attenuation and interaction of blasting waves at the interface zone, as shown in Figure 9. However, the shaking table test and field investigations show that the weak intercalated layer could be the location where cracks initiate.
is could be because the weak properties of the interlayer fail to resist large blasting wave energy. Further dynamic loading would trigger the propagation of cracks, causing the tensile shear failure along the sliding surface. erefore, it is of great importance to identify the location, dimension, and structural properties of weak intercalated layers before the mining operation in the field. Studies such as [47] have been performed in this research topic.
Blasting loads are investigated in our study to benefit the engineering practice of underground mining or open-pit mining. While the dynamic failure modes and mechanisms of coal slopes are also paramount to its aseismic design, future studies could focus on the seismic response and deformation mechanism of coal slopes with weak intercalated layers under a strong earthquake.

Conclusions
Numerical methods, shaking table tests, and field tests were performed to study the dynamic response characteristics of a coal rock slope with weak interlayers under blasting loads. e following conclusions can be made: (1) e weak interlayer has a damping effect on the propagation of blasting waves in the slope. e analyses of equivalent stress show that weak interlayer effectively improves the stress state of the slope under blasting loads, and the weak interlayer makes the slope easier to produce creep deformation.
(2) e shaking table test and field monitoring test results identify the natural frequency of the slope to be approximately 35 Hz. When the vibration frequency of the blasting wave is 35 Hz, the sliding body will slide along the weak bedding layer. e slope's natural frequency can be potentially determined by monitoring the maximum slope deformation when the slope is under blasting vibration with a series of frequencies.
(3) e specific deformation and failure process of the coal slope is shown as follows: first of all, cracks were formed in weak interlayers by blasting vibration, and the cracks expanded to form shear zones, which then led to vertical tension cracks on the slope crest. With the increase of vibration load, the cracks gradually penetrated and formed a sliding surface composed of tension and shear. e slip surface is formed by the combination of the shear slip surface of the weak bedding interlayer and the tension crack generated by the acceleration amplification at the top of the slope. (4) Slope dynamic deformation is negatively correlated with the distance of the blasting source. e farther the distance from the blasting source is, the less the displacement, PGA, and PGD of the slope gradually decrease (5) Future studies could focus on the seismic response and deformation mechanism of coal slopes with weak intercalated layers or other complex geological structures under a strong earthquake.

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 that there are no conflicts of interest regarding the publication of this paper.