A Study on the Impact of Longwall-Mining Operation on the Stability of a Slope-Pillar Structure

'e slope-pillar structure is defined as an artificial construction shaped by a surface highwall mining and an underground longwall mining. In order to study the impact of stress variations (induced by the longwall-mining operations) on the stability of a slopepillar structure, No. 30101-1 working face of Nanliang Coal Mine was selected as the study subject. A series of analysis including physical similar simulation test and the theoretical analysis were conducted to study the movement pattern of the overlying strata in the working face of the slope and the loading and the stability of the pillar column revealing the mechanism of dynamic load instability of nonuniformly distributed load “slope pillar.” A “slope-pillar” reverse sliding structural model was proposed along with a newly established “slope-pillar” structural mechanical model and a formula to calculate the boundary width of the protection for the pillar under the extreme balance conditions. According to the study, the width of the protective pillar is designed at 20m, and the roadway deformation during the working face and the stability of the pillar can satisfy the safety requirement concerning the working face, which further validates the theoretic deduction formula.


Introduction
e Shenfu-Dongsheng Coal Field, as an important coal energy base in China, is featured with shallow coal seams and diverse topography. Since the early 1990s, scholars have conducted a series of studies on shallow coal seam mining. e authors of [1][2][3][4][5][6] conducted a large number of on-site monitoring and tests to determine the mine pressure, the fracture law of the overlying strata, and strata control. Due to the staggered surface ravines in the western region, the drastic variance in the occurrence of slopes, and the interaction between the special topography and the shallow coal seam mining, a series of new features such as the appearance of underground pressure and the movement of the overburden at the working face were identified. Regarding the mining pressure and strata control under the special topography, Xufeng et al. adopted the roof structural study to the coal seaming under the slope of bedrock and sandy soil, resulting in the establishment of the depth of gully slope and the impacts of the slope angle on the mining pressure of the working face [7]. Taking the coal seam conditions under the Shenfu-Dongsheng Coal Field into account, Zhang et al. conducted simulation and concluded that a greater load from the working force was imposed on the upper slope as the depth of the gully slope increased [8][9][10]. Hu et al. suggested that several factors including the existence of the key strata, the depth of trench, and slope angle tended to influence the generation of dynamic loading pressure on the working face [11]. Wang et al. established a "support-surrounding rock" model under the gully in combination with the geological conditions of the mine and identified the intensity of the mineral pressure following the pattern of "back ditch section > ditch section > ditch bottom section > normal mining section" [12]. Zhang et al. established the structural mechanics model of the "unevenly distributed load" of the key layer under the valley terrain and obtained the calculation formula of the "slide rotation load" caused by the instability of adjacent key blocks during the movement of the overlying strata [13]. e aforementioned studies concerning the movement of the overlying strata and mining pressure law under the special topography during the shallow coal mining mainly focus on the naturally developed terrain and trenches [14]. However, due to the upper thin overlying soil layer of the shallow coal, coal outcrop often occurs around the well zone and the gully region, which is often treated with excavation and refilling, resulting in manmade slope around the boundary. Current studies regard the man-made slopes and mining under the man-made slope center around the sliding and instability mechanism and theories during the coordinated open-pit mining. Sun et al. studied the sliding mechanism and deformation mechanism of the combined mining slope of underground and open-pit mining [15,16]. Ding et al. proposed four stages of slope deformation and failure evolution under the coordinated mining of open wells and revealed the deformation failure mode and instability mechanism of mining slopes under the coupling effect of "key hard rock layers" and "key weak layers" [17,18]. However, the existing research studies discussing slope slippage, overlying rock movement, and coal pillar stability involved in underground mining of slopes were mainly focused on single-object issues, and few discussions and studies have been performed to explore the coupling interaction among all three factors during the mining process. rough the theoretical analysis and similar simulations as well as the field tests, this paper discusses the interaction of coal pillar, overlying rocks, and slopes in the same space during the mining process of the manmade slope formed after the excavation and backfilling of the corner coal treatment area.
is study is expected to provide some references for mining under similar special mining conditions [19].

The Movement Law of the Overlying
Strata under the Slope 2.1. Project Background. As a typical loess hilly and gully landform, the Nanliang Coal Mine is featured with a complex topography, some vertical and horizontal gullies, deep ridges and steep ditches, gentle stratum occurrence, shallow coal seams, and simple geological structure. e working face #30101 locates at the southwest of the mining field. e depth of the coal seam ranges from 80.43 to 112.01 m, with a strike length and dip length of 1398.32 m and 280 m, respectively. Besides, the average thickness of the coal seam is about 2 m. e transportation of the working face is close to the boundary of the minefield along the south side of the channel, which is labelled as the corner coal of 3-1 coal seam. e overlying rock layer is thin with well-developed sandstone fissures, featured with excessive illegal exploration. In order to prevent the surface collapse, the ground fissures from further development, and the illegal exploration in this area, which tends to trigger a spontaneous combustion of the residual coal and 3-1 coal seam goaf, a comprehensive treatment area has been established with a total of 5.05 × 10 −2 km 2 bottom area through the excavation treatment, as shown in Figure 1.
As shown in Figure 2, the bottom of the comprehensive treatment area of 3-1 coal seam and No. 30101-1 working face formed a slope, which was perpendicular to the working face. e slope angle is 53°.

e Experimental Model Design.
During the study of the interaction between the overlying strata, the slope, and the protective pillar during the advancement of the working face under the shallow coal seam slope terrain, a physical simulation method was adopted to simulate the mining process of the 30101-1 working face. In order to minimize the influence of model size, the experiment adopted a geometric similarity ratio of 1 : 100 and established a model of 3 m in length, 0.2 m in width, and 0.9 m in height, as shown in Table 1. Considering that the main material at the pit bottom of the comprehensive treatment area was dominated with the loose loess and gangue, some strata were removed from the comprehensive treatment area and used as the backfill material in the model.
As demonstrated in Figure 3, the spontaneous loads of the overlying strata and slope instability during mining were monitored and recorded via the CL-YB-152 dynamic load data acquisition system which was numbered from left to right starting from 1. A total of 9 displacement detectors were installed above the 3-1 coal seam to monitor the movement of overlying rock and slopes during mining. A number of percentage meters were installed on the surface and in backfill area to monitor surface movement changes during the mining. e development of cracks and fissures during the mining was closely monitored via a BJQF-1 crack width gauge.

2.3.
e Movement Features of the Slope during Mining Process. During the simulation, before the mining working force approached the slope, the movements observed during the mining were similar to the movements of the working face during the regular shallow coal mining. e mining triggered the initial fracture in the basic roof and a series of cyclic fractures. e fractured rocks were limited by the lateral horizontal force from the surrounding strata, resulting in a relatively stable hinged structure. No interaction was observed between the overlying strata and the slope during the movements under the mining. e loess layer on the surface was subject to the impacts of the mining face. e fractures caused by the advancement of exploration continued to expand to the surface, forming a fracture face on the vertical profile. e fracture face continued to develop cyclically, approaching to the slope. Simultaneously, the loess on the surface under the advancement of the exploration and the additional stress corresponding to the horizontal deformation caused by the exploration overcome the cohesion or the tensile strength of the soil along with greater deformation in the surface than the depth of the loess. As a result, the fracture on the surface reached the ultimate equilibrium of the shear failure than that on the deep layers [20]. Details can be viewed in Figure 4.
As the working face approached close to the slope, the overlying strata triggered by the working face impacted the stability of the slope significantly. In simulation, when the distance between the working face and the slope foot was reduced from 20 cm to 18 cm, the slope body was under the impacts of the mining working face. e slope experienced rotation and sliding, with the vector direction pointing to the goaf of the working face. Considering that the backfilling zone was dominated with the loose materials whose resistance to the pressure was nearly zero, the backfilling zone was pressured and squeezed under the structural rotation of the slope, resulting in the heave phenomenon at the bottom of the backfilling zone. e maximum displacement observed at the bottom of the backfilling zone was with a length of 20 mm. e joint between the backfilling zone and slope surface started to develop fissures, resulting in a massive   through crack with an average width of 1.60 mm, allowing the sliding of the slope. As the simulation indicated, under the circumstance that the working face advanced within the impact range of the slope mining, the slope body tended to develop fissures, slide, and generate bigger cracks under the mining [21]. Along with the cracks before mining, a sliding rotation load structure known as "slope pillar" was formed with the pillar as the main load carrier, as shown in Figure 5. e protective pillar for the slope needs to bear not only the weight of the overlying strata but also the dynamic loading produced during the rotation of the slope toward to the goaf, known as "sliding rotation load" [13]. In the simulation test, when the width of the protective pillar was below 20 cm, the instability of the slope body tended to occur, with the crack surface as the boundary and the inner side of the pillar (close to the working face) as the supporting point to facilitate the process of sliding and rotation. During the movement, the hinged pieces experienced instability and sliding, resulting in the "sliding rotation load" transferred from the back side of the strata to the pillar, which led the load to the inner side of the pillar and dynamic mining pressure [22].

e Stress Monitoring of the Slope Pillar under Mining.
In order to further study the stress variance pattern of the protective coal pillar in different widths, the stress data of the protective pillar collected from the monitoring point 4 to 8 when the working surface was less than 100 m from the slope body were sorted and analyzed (the position of the measuring points is shown in Figure 3). Based on the collected vertical stresses born by the floor collected from the monitoring point 4 to 8 via CL-YB-152 dynamic load data acquisition system, a stress curve is developed and presented in Figure 6. According to Figure 6, as the width of the pillar decreased, the stress monitored at point 4 gradually increased to 3.6 MPa. Under the slope instability, the stress decreased to 3.2 MPa. e stress trend observed at the monitoring point 5 followed a similar pattern. A limited stress increase was observed at the point 5 during the slope instability. e monitoring point 6 was located below the inner side of the slope pillar. As the mining advanced closer to the pillar of the slope, the stress increased accordingly, which continued to rise as high as 5.1 MPa during the slope instability. As the working face advanced to 18 m away from the slope, a hinged rock beam was formed above the working face, resulting in a minimum stress at monitoring point 7.
e monitoring point 8 was located behind the working face with a consistent vertical stress, which was still higher than the original stress of the rock. e analysis of the vertical stresses has suggested that when the width of the protective pillar of the slope was smaller than 20 m, driven by the self-weight, the slope tended to rotate and slide toward the goaf in the same space, resulting in instability. e strata fissures above the goaf zone generated less and less horizontal pushing force, resulting in the occurrence of dynamic load when the slope slid and rotated. e dynamic load was observed at the monitoring point 6 installed in the internal side of the slope pillar, featured with drastic stress increase.   e static load of the protective pillar in the slope area demonstrates a nonuniform trend. Based on previous studies [7,8,13,23], assuming that the structure of the slope is single and the slope angle remains unchanged, following the relationship identified between the protective pillar and the slope projection position, two types of side slope-pillar mechanical models have been established, as shown in Figure 7. In Figure 7, F 1 refers to the supporting force of the external side of the pillar to the overlying strata, F 2 refers to the supporting force of the internal side of the slop pillar to the overlying strata, and F 3 refers to the supporting force for broken block. e width of the slope pillar is equal to or greater than the horizontal projection length of the slope, as shown in Figure 7(a). According to the structural mechanics, taking the supporting point of the outside of pillar as the point O as the research object, it is concluded that Taking the internal side of the pillar as the supporting point, the following equation can be obtained: In formula (3), P a refers to the nonuniform distribution load, and P b stands for the uniform distribution load, which can be calculated according to formula (4).
ρg refers to the average body force of the overlying rock layer of the key layer, KN/m 3 ; v a indicates the volume of the overburden with nonuniform load, v a � 0.5a 2 tanα (α the slope angle of toe); v b is the volume of overlying strata with uniform load; and a suggests the width of the slope pillar, m; e width of the slope pillar is less than the horizontal projection length of the slope, as shown in Figure 7    Shock and Vibration 5 Taking the internal side of the pillar as the supporting point, the following equation can be obtained: When the broken block slips and loses stability, F 1 and F 2 can be calculated by formulas (1), (2), (5), (6).
When the slope pillar is unstable, based on the analysis of the movement features of the slope during mining process and the stress monitoring results of pillar, the protective pillar for the slope needs to bear both the weight of the overlying strata and the instantaneous slipping and rotating load F h produced during the rotation of the slope toward to the goaf.
In formula (9), φ refers to the instantaneous impact coefficient; δ indicates the rotation angle; and G B is the overburden load above the slop pillar.
According to (7) and (8), under the action of the static load of the side slope, the load on the internal side of the pillar is greater than that on the external side of the pillar. When the pillar is not affected by the instantaneous slip and rotation load, the pillar only bears the load of the overlying rock layer and the instability failure of F 1 > R coal pillar under static load (R is the ultimate strength of the coal pillar). When the pillar is affected by the instantaneous slipping and rotating load, the pillar needs to bear not only the static load of the overlying rock but also the instantaneous slipping and rotating load, resulting in F 1 + F h < R, which suggests that the pillar is likely to undergo instantaneous slipping load instability and failure.

The Stability Analysis of the Protective Pillar of Slope
According to the "slope-pillar" unevenly distributed load structure dominated by the pillar as the main load bearer, obtained from the simulation (combined with the research results of reference [16,24,25]), a "slope-pillar" structural mechanical model was established to analyze the stability of the pillar of the slope, as demonstrated in Figure 8. As demonstrated in Figure 8, the fissure face AB and premining fissure BE jeopardized the integrity of the slope, resulting in the development of characteristics of double slip surfaces in the slope and stability reduction. Due to the loess and gangue used in the backfilling zone III, the slope reinforcement was limited. Without considering the hydraulic effect and the reinforcement of backfilling zone to the slope foot, the length of the slip surface formed by the width of the pillar and the length of the broken section determine the slope stability [26,27].
According to the slope structure, the length, inclination, cohesion, and internal friction angle of AB are designated by L 1 , β 1 , c 1 , and φ 1 , respectively. e length, inclination, cohesion, and internal friction angle of BC are designated by L 2 , β 2 , c 2 , and φ 2 , respectively. DC refers to the slope surface with a slope angle α.
Assuming after point B, the inclined angle θ between the shear surface of the soil body of the mining and the surface of the maximum principal stress created a boundary surface BE which divided the slope body into part I and II, corresponding to weight W 1 and W 2 . e detailed calculation of W 1 and W 2 is listed in the following formula: Due to the instability and sliding of the fissure blocks in the overlying strata during the slope instability, the horizontal pushing force toward the slope can be overlooked [25]. Assuming that no additional external force working on the slope, the friction from the bedrock imposed on the slope part I can be calculated as follows: P 1 and P 2 refer to the interaction force between part I and part II, respectively, which equals to each other with opposite direction. P 1 can be calculated as follows: e slope safety coefficient F is defined as the ratio between the sliding resistance force of the slope structural face and sliding force, which means F is the sliding resistance force/sliding force [25], which leads to the calculation formula of slope safety coefficient F 2 for part I.
Because β 2 � 0°, the calculation formula of slope safety coefficient F 2 for part II is developed and presented in the following formula: According to the force analysis conducted previously, the safety coefficient of part I was obtained under the extreme equilibrium. Assuming the slope body is a rigid structure as a whole, F � F 1 � F 2 � 1 can be obtained.
According to formula (15), the width of pillar of the slope under the slope extreme equilibrium can be calculated based on the following formula: Since P 1 � P 2 , According to the corresponding conditions of Nanliang No. 30101-1 working face, the following data were obtained: Applying the corresponding data listed above into formulas (10), (11), (12), and (17), the width of the pillar is calculated as follows: According to the actual dimensional requirements in the field, the width of the retention pillar cannot be less than 20 m.

Engineering Validation
To  Figure 8: e schematic diagram of structural stability analysis of "slope pillar." Shock and Vibration measurement stations 1 to 5 were set in the lower section of the nonslope roadway, and the measurement stations 6 to 10 were set in the lower section of the slope roadway. e measurement station layout is illustrated in Figure 9. e comprehensive measuring station adopts the cross-point method, and the CXK-12 mine intrinsically safe borehole imager method was adopted to monitor the approaching amount of the two sides and the amount of roof subsidence of the 30101-1 auxiliary transportation roadway and the degree of damage to coal pillar when the working face passes the slope and the nonslope during the mining process.

D E A
After analyzing the collected data of the surrounding rock displacement of the roadway, the curves between the working face advancing distance and the subsidence of the roof of the mining roadway and the displacement of the both sides of roadway, respectively, were obtained, as shown in Figure 10. Comparing the roadway roof subsidence and the displacement in both sides of roadway between the slope section and nonslope section, it is revealed that the deformation of the surrounding rock in the slope section is larger than that in the corresponding nonslope section. e maximum difference of the roadway floor subsidence  Figures 11 and 12, the inner wall of the borehole was relatively intact, and the borehole did not show any severe deformation, block loss, hole collapse, and so on. e borehole wall of No. 2 was mainly developed with annular cracks, and the borehole wall of No. 7 was dominated with the radial cracks. e crack depth was greater than the depth of No. 2 borehole. e obtained results of the boreholes at each station were collected to develop the broken line diagram of the crack development position, as shown in Figure 13. Before the mining face advanced to the slope, the maximum depth of cracks in the borehole was 85 cm with an average depth of 76.8 cm. After the mining face advanced into the slope, the maximum depth of cracks in the borehole reached 95 cm with an average depth of cracks of 88 cm, demonstrating increased crack depths with limited variances in fragmentation. e analysis reveals that when the width of the protective pillar of the slope reaches 20 m, the deformation of the roadway surrounding rocks under the influence of the mining within the slope tends to increase at an acceptable safety level, with deeper cracks developed in the pillar. However, the fragmentation in the pillar is limited without jeopardizing the stability of the pillar and no dynamic load instability occurred in the coal pillar. e result can meet the safety requirements on-site.

Conclusion
(1) e simulation has revealed that as the working face advanced and the width of the slope protective pillar reached the boundary width, the "slope-pillar" structure formed by the slope, pillar, and fissure face tended to experience rotation, sliding, and instability. e concentrated stress points were identified at the slope foot. Eventually, the cracks grew, forming a massive through fissure. e backfilling zone was pressured and squeezed by the "slopepillar" structure, resulting in an increased stress and the rib spalling on the pillar. (2) rough similarity simulation on the movement of slope strata and the monitoring on pillar stress during mining process, a mechanical model of nonuniform loading of "slope pillar" was established, which demonstrated that the rotation slip of the overlying strata of the slope under the limit equilibrium condition was a contributing factor to the dynamic load instability of the "slope pillar." e criterion for judging the dynamic load instability of the "slope pillar" was given as F 1 + F h < R.
(3) A "slope-pillar" mechanical model was established along with the boundary width calculation formula for the "slope-pillar" structure under extreme equilibrium. e formula was used to calculate the minimum 20-meter width of the protective pillar for No. 30101-1 working face in Nanliang mining zone.  Shock and Vibration 9 (4) e engineering validation indicated that when the width of the protective pillar reached 20 meters, the average cracks depth inside the pillar was 88 cm, no dynamic load instability occurred, and the deformation of the roadway surrounding rock was all within the controllable range, which also further reflected the validity of the boundary width formula for the pillar under the extreme balance conditions.

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 they have no conflicts of interest.