Failure Mechanism and Stability Control Technology of Slope during Open-Pit Combing Underground Extraction: A Case Study from Shanxi Province of China

: With the development of society, the demand for mineral resources is gradually increasing, and the current situation of decreasing total resources dictates the inevitable interaction between open-pit combing underground extraction (OPUG) in time and space. In this research, we took the Anjialing coal mine in Shanxi Province of China as a case study, and tested the physical and mechanical properties of coal rocks in the laboratory. The similarity criterion was used to build a similar experimental model for the deformation evolution of the slope of the open-pit mine section; the digital scattering method was used to test the inﬂuence of the underground mining process parameters on the deformation evolution of the open-pit slope. The results showed that there was an obvious distribution of “three zones” above the mining goaf, namely, a collapse zone, fracture zone, and slow subsidence zone. When the mining face was continuously advanced towards the bottom of the open pit, the supporting stress of the mining face transferred to the side of the open-pit slope. Additionally, large displacement and stress concentration were observed on the slope near the stoping line, which caused the slope body to move along the uppermost part of the slope ﬁrst, and thereafter along the lower part. Various techniques for slope stability control are discussed, including the optimization of spatial and temporal relationships between open-pit and underground mining, the optimization of mining plans, and the use of monitoring and early warning systems. The results can provide a guide for slope stability control of similar open-pit mines in the process of mining coal resources.


Introduction
With the continuous advancement of surface mining technology, the proportion of surface coal production has been increasing in the past decades. Restricted by factors such as mining technology and stripping ratio, the mining scope of surface mines is limited to the upper section resources of the entire ore body, and resources not suitable for surface mining at depth are completed by underground mining means [1][2][3]. For horizontal or near-horizontal coal mines with miscellaneous coal seams, it is economically unreasonable to open-pit mine the deeper resources, and a combination of open-pit mining method and underground extraction method is generally used [4][5][6][7]. This form of mining is also known Sustainability 2022, 14, 8939 2 of 21 as the method of open-pit combing underground extraction (OPUG) [8][9][10][11][12]. In the process of OPUG, the effects of underground and opencast mining are both mutual superimposition, forming a large dynamic composite system, so the deformation mechanisms of slopes and rock bodies are more complicated, and the slope already formed within the opencast mine will be affected by the disturbance of the underground mining. Settlement deformation and weak layer activation will continuously occur, posing a serious threat to the stability of the slope and safe production [13][14][15][16].
To take necessary and timely strengthening and prevention measures and enable the mining enterprise to obtain optimal economic benefits, it is essential to obtain knowledge of slope mechanisms and deformation characteristics under the mutual influences of open-pit combined underground mining. For conventional slope stability in actual open-pit combined underground mining engineering, the impact of open-pit mining on slope deformation is often analyzed separately, without taking into account the factor of underground mining. Presently, numerical simulation methods are commonly applied for analysis and evaluation, and remarkably critical results have been obtained [17]. Stead et al. analyzed rock slope failure mechanisms in a complex stress environment by using modified numerical modeling techniques [18]. Scholtès and Donzé simulated the expansion and evolution process of fissure-bearing rock masses using the three-dimensional discrete unit method and conducted stability analysis for a specific fissure-bearing rock slope engineering example [19]. Shen and Karakus proposed the Hoek-Brown criterion-based finite element nonlinear strength reduction method for three-dimensional rocky slope stability analysis using the FLAC3D program [20]. Amini et al. divided the rock masses in the reverse slope into two parts: one part was considered as completely disconnected from the parent rock, and the other part was considered as intact rock masses connected with the parent rock [21]. Johari and Lari divided the damage modes of rocky slopes into four types, explored the influence of different heights and damage modes on the instability probability of slope systems, and established a model to directly assess the instability probability of rocky slopes. The influence of different heights and damage modes on the probability of slope system instability was discussed, and a model for directly assessing the probability of rocky slope instability was established [22]. Zhou et al. used RFPA-2D to study the stability of multiple sections of open-pit intersections containing faults and concluded that the presence of faults changes the stress distribution, damage patterns, and water flow paths [23]. Luo et al. used the stochastic finite element method (FEM) to analyze the probability of rock mass stability of mine slopes and underground quarries, and the numerical simulation results showed that the mining of underground quarries had some influence on the stability of mine slopes [24]. Dagdelen and Traore determined the optimum mining depths for six open-pit and shaft mines based on mixed integer linear programming (MILP) [25]. Gischig et al. investigated the significance of kinematic release mechanisms on planar translational slope failure using the three-dimensional distinct element codes [26]. They studied the three-dimensional block shape and volume of the unstable rock masses simulated with different discontinuity set geometries. Li et al. [27] showed through three-dimensional numerical simulation that plastic damage occurs gradually with increasing depth, and a large amount of tensile damage occurs at the top of the slope when the mining depth reaches the edge determined by the seismic exploration method. These research results can solve the related problems of conventional slope stability analysis, among others. Meanwhile, there are still few studies on the joint mining process. Jiang et al. [28] studied the maximum mining capacity of different shafts in the Shirengou iron ore mine under open-pit combined mining conditions based on ore loss rate and grade. Cao et al. [29] studied the impact of underground mining on the mining environment of opencast coal mines based on a two-dimensional heat transfer equation with a phase change component. Wang et al. [30] concluded that mining of shallow buried coal seams causes shear damage to polygonal rock masses, leading to the instability of articulated structures, and that increasing slope and mining height leads to the more rapid development of fractures within gullies, which has a significant effect on the stability of gully slopes.
Due to the complex engineering environment of open-pit joint mining, the deformation value and deformation range of slope rock have far exceeded the original single mining theory understanding. Some scholars have also used model tests, field monitoring, or integrated means to study the deformation and damage law of slope rock under the combined surface and underground mining conditions to analyze the coordinated interaction between the two mining effects. Li et al. [31] have concluded through field investigations, centrifugal modeling, and numerical simulations that the subsurface mining area causes overturning damage to the upper rock masses, resulting in sliding shear of the underlying rock. Singh investigated the vibration response of roof slabs to adjacent open-pit blasting by investigating seven underground coal mines in India. Kovrizhnykh [32] used an elastic creep model to solve the problem of deformation and damage of rock masses in shaft and surface mining under hydrostatic pressure, determining the time to complete damage under long-term stress. The above studies have shown that the interaction between open-pit and underground extraction poses a threat to the safe and orderly production of open-pit mines and that the slope instability and ecological problems caused by such coordinated open-pit and underground extraction constitute a major potential hazard to large-scale open-pit mining. Although some research work has been conducted on coordinated open-pit and underground extraction, systematic results have not yet been seen, and the prevention and control measures for slope hazards of coordinated open-pit and underground extraction under the action of complex factors are still in the initial stage. It has become a research trend to carry out systematic theoretical research on the safety evaluation, monitoring, prediction, forecasting, and treatment of mining slopes, and to provide safety measures and optimization methods for combined surface and underground mining engineering.
Therefore, the engineering geological conditions and mining conditions of the Anjialing II open-pit mine are presented here. We adopted a combination of theoretical research, model tests, and mechanical parameter testing experiments to investigate the spatial and temporal relationship between the mining face and slope of the open-pit mine in the process of OPUG. The deformation and damage characteristics of the slope of the open-pit mine under the influence of pit mining disturbances, the deformation and damage process and instability mechanisms of the mining slope, and the slope stability control methods were investigated. The research results can provide an important basis for improving the resource recovery rate, reducing resource waste, and optimizing the industrial layout, with broad application prospects and significant economic and social benefits. The Anjialing opencast coal mine located in Shanxi Province of China was selected as the study area ( Figure 1). Founded in 1982, the coalfield is a large coal production and energy export base of 100 million tons. The mining area is bordered by the Taihang and Wutai Mountains to the east and Inner Mongolia to the west and has a cold, dry continental climate with particularly severe sand and wind. The area is in a lower temperature state all year round, with an annual average temperature of 5.4-13.8 • C and a large daily temperature difference of 18-25 • C. The average annual sunshine period is 2693.3 h, with a minimum of 2444.5 h and a maximum of 2883.1 h. The maximum humidity is 80% and the minimum is 0. Precipitation is concentrated between July and September and rainfall is unevenly distributed, with annual average precipitation of 426.7 mm and annual average evaporation of 2006.7 mm.

Location of the Anjialing North and the Well-Working Coal Mining II Mine
At present, the north gangue of Anjialing is mainly affected by the sinking of the 9 # coal seams mining face of well-working coal mining II mine area, which is manifested by the sinking of the slope of the north gangue of the 9 # coal No. 29 process, the overlying coal seams form "three zones", i.e., a collapse zone, fracture zone, and subsidence zone. The direct topping of the 9 # seam coal in the shaft working mine resulted in the formation of funnel-type subsidence in the flat pan above 1405 in the north end of the open pit. The fissure zone and bending zone caused by the mining of the 9 # coal resulted in fissures in the flat pan above 1360 elevation in the north end of the open pit. The impact on the geotechnical body of the north end of the open pit is mainly manifested in the reduction in the physical and mechanical parameters of the geotechnical body of the north gangue caused by the disturbance of the shaft working mine, resulting in the reduction in the physical and mechanical parameters of the geotechnical body of the north gangue. The impact on the rock and soil body of the north end is mainly due to the lowering of the physical and mechanical parameters of the rock and soil body of the north gangue as a result of the underground mining disturbance, which will lead to the development of joints and fissures in the rock and soil body of the north end. In addition, the 1330, 1360, and 1390 elevation pans are the main transportation pans in the open pit. If large fissures and funnel-type subsidence occur, the safety of the open pit will be seriously affected. The spatial location of the north gangue of the Anjialing open-pit mine and the 9 # coal of the well-working coal mining II mine are shown in Figure 2. At present, the north gangue of Anjialing is mainly affected by the sinking of t coal seams mining face of well-working coal mining II mine area, which is manifeste the sinking of the slope of the north gangue of the 9 # coal No. 29209 mining face, No. 2 mining face, and No. 29211 mining face by back mining. During the mining proces overlying coal seams form "three zones", i.e., a collapse zone, fracture zone, and su

Testing of Geotechnical and Physical Properties in the Study Area
The understanding of the rock parameters and geological conditions is beneficial for decreasing the potential for geo-hazards during engineering [25,28,[33][34][35][36]. The physical and mechanical parameters of the rock and soil were tested by collecting soil and rock from the slope of the open pit, making standard soil and rock specimens, using a straight shear meter, a rock shear resistance meter, and a triaxial shear tester. The test specimens, equipment, and test results are shown in Figure 3 and Table 1.

Testing of Geotechnical and Physical Properties in the Study Area
The understanding of the rock parameters and geological conditions is beneficial for decreasing the potential for geo-hazards during engineering [25,28,[33][34][35][36]. The physical and mechanical parameters of the rock and soil were tested by collecting soil and rock from the slope of the open pit, making standard soil and rock specimens, using a straight shear meter, a rock shear resistance meter, and a triaxial shear tester. The test specimens, equipment, and test results are shown in Figure 3 and Table 1.

Testing of Geotechnical and Physical Properties in the Study Area
The understanding of the rock parameters and geological conditions is beneficial for decreasing the potential for geo-hazards during engineering [25,28,[33][34][35][36]. The physical and mechanical parameters of the rock and soil were tested by collecting soil and rock from the slope of the open pit, making standard soil and rock specimens, using a straight shear meter, a rock shear resistance meter, and a triaxial shear tester. The test specimens, equipment, and test results are shown in Figure 3 and Table 1.   Twenty soil samples, numbered ZJ1-ZJ20, were collected ( Table 1). The physical and mechanical properties of the soil samples were tested in the laboratory and the samples were classified and named accordingly. The weathered bedrock specimens collected, numbered KJ-1-KJ10, are shown in Table 2. The physical and mechanical properties of the weathered bedrock specimens were tested in the laboratory and the specimens were named accordingly.

Model Tests
Model tests were conducted using the similarity principle to study the evolution of open-pit slope deformation under different well mining conditions. Using the similarity criterion, a similar experimental model was constructed for the evolution of slope deformation in coordinated open-pit and well mining. Using the digital scattering method, the influence of well mining process parameters on the evolution of slope deformation in the open pit was tested to provide parametric support for the stability of slopes in the process of OPUG.

Test Methods and Equipment
A similar material simulation test has the advantages of intuitive results and repeatable analysis of influencing factors, which plays an important role in guiding mining. A similar material simulation experiment of slope stability in the Anjialing open-pit mine under the effect of coordinated open-pit and shaft mining includes a similar material model, a BW-5 micro pressure box, and an XJTUDP 3D optical photogrammetry system and a digital scatter monitoring system. The digital scatter correlation method is a non-contact photometric test mechanics method that determines the displacement and strain of an object's surface by obtaining information about its deformation. The main equipment used include an image recording computer, a high-resolution CCD camera, an image processor, etc. Figure 4a shows the BW-5 pressure box and Figure 4b shows the YJZ-32A digital strain gauge. Figure 4c shows the XJTUDP 3D optical measurement system, and Figure 4d shows the digital scattering light measurement system, etc.

Test Methods and Equipment
A similar material simulation test has the advantages of intuitive results and repeat able analysis of influencing factors, which plays an important role in guiding mining. A similar material simulation experiment of slope stability in the Anjialing open-pit mine under the effect of coordinated open-pit and shaft mining includes a similar materia model, a BW-5 micro pressure box, and an XJTUDP 3D optical photogrammetry system and a digital scatter monitoring system. The digital scatter correlation method is a noncontact photometric test mechanics method that determines the displacement and strain of an object's surface by obtaining information about its deformation. The main equipment used include an image recording computer, a high-resolution CCD camera, an im age processor, etc. Figure 4a shows the BW-5 pressure box and Figure 4b shows the YJZ-32A digita strain gauge. Figure 4c shows the XJTUDP 3D optical measurement system, and Figure  4d shows the digital scattering light measurement system, etc.

Theory and Methods
Similar materials experiments are conducted using similar materials, following similar theories and guidelines, to build a model similar to the site and then carry out a simulated mining study to observe the movement of the overburden and the distribution of support pressure during mining. The model is then used to observe the movement of overburden and support pressure distribution during mining, and to calculate or extrapolate the pattern of roof movement and support pressure distribution during actual mining under these conditions.

Similar Material Parameters
Using the similarity principle, similar model similarity parameters are obtained based on geometric similarity, kinematic similarity, and dynamic similarity calculations, as shown in Table 3.

Proportioning of Model Materials
A distilled water solution equipped with cement, lime, kaolin, paraffin, and retarder was used as the cementitious material and calcium carbonate. Gypsum and silica sand were used as aggregates and similar models were constructed according to the actual slope range of the studied open-pit mine. According to the mechanical parameters of the different rock layers of the slope, similar materials were made by adjusting the ratios of lime, gypsum, and silica sand, as shown in Table 4.

Point Settings and Model Boundaries
The size of the model table is 5000 mm × 300 mm × 2000 mm. Considering the influence of well mining on the evolution of slope deformation, the displacement monitoring points are laid out according to 200 mm × 100 mm. In the process of model building, the BW-5 pressure box was pre-buried to test the stress value, and the YJZ-32A intelligent digital strain gauge tested the deformation value to achieve real-time monitoring of stress and displacement during the advancement process. Seven measurement lines parallel to the coal seam were laid out on the overlying rock layer at the top of the coal seam, with equally spaced measurement points on each line, making a total of 118 optical displacement measurement points with a horizontal spacing of 200 mm and a vertical spacing of 100 mm. Twelve stress test points were also laid out in the rock layer 8 m from the top of the coal seam with a horizontal spacing of 100 mm.
The model boundary conditions are displacement constraints on the left, bottom, and top edges of the model, free surface on the front of the model, and horizontal displacement constraints on the back of the model. The initial stress state is simulated by the self-gravity of the model seam, and the well mining pressure and slope deformation are simulated by the advancement of the coal seam. Figure 5 shows the arrangement of displacement measurement points and stress tests for similar models.

Point Settings and Model Boundaries
The size of the model table is 5000 mm × 300 mm × 2000 mm. Considering the influence of well mining on the evolution of slope deformation, the displacement monitoring points are laid out according to 200 mm × 100 mm. In the process of model building, the BW-5 pressure box was pre-buried to test the stress value, and the YJZ-32A intelligent digital strain gauge tested the deformation value to achieve real-time monitoring of stress and displacement during the advancement process. Seven measurement lines parallel to the coal seam were laid out on the overlying rock layer at the top of the coal seam, with equally spaced measurement points on each line, making a total of 118 optical displacement measurement points with a horizontal spacing of 200 mm and a vertical spacing of 100 mm. Twelve stress test points were also laid out in the rock layer 8 m from the top of the coal seam with a horizontal spacing of 100 mm.
The model boundary conditions are displacement constraints on the left, bottom, and top edges of the model, free surface on the front of the model, and horizontal displacement constraints on the back of the model. The initial stress state is simulated by the self-gravity of the model seam, and the well mining pressure and slope deformation are simulated by the advancement of the coal seam. Figure 5 shows the arrangement of displacement measurement points and stress tests for similar models.

Experimental Program
According to the actual situation mining conditions, similar material tests were conducted to study the effect of conditions such as continuing mining after the working face of the 4 # and 9 # coal seams in the open pit was advanced to the corresponding stoping line on the stability of the slope. At the same time, the overlying rock stress (initial incoming pressure, periodic incoming pressure step, etc.) and deformation law were obtained to determine the reasonable stoping line location, collapse angle, and the first appearance of ground deformation results. Finally, according to the results, the deformation and damage process of the slope of the north gangue is analyzed, and the interaction law between the reasonable stoping line of the shaft mine and the open pit is studied. Then, the factors of cutting the slope of the north gangue to reduce the weight and the pressure foot of the inner row are considered, and the influence of measures such as the pressure foot of the inner row on the stability of the north gangue is analyzed. The interaction between the reasonable stoping line and the open pit is investigated, and the impact of measures such as internal pressure footing on the stability of the north end is analyzed, taking into account the factors of slope cutting and weight reduction of the north end.  When the mining face advanced to the initial stoping line (i.e., 486.5 m from the foot of the slope), the 9 # coal collapse zone was connected with the 4 # coal. When the mining face advanced to reasonable stoping line 1 and reasonable stoping line 2, the overburden rock damage and deformation became more and more adequate. By the time the mining face advanced to the final stoping line, the overburden rock within the mining area as a whole slipped towards the side slope, and the steps of the open pit side slope at +1363 m level and above were affected.

Analysis of the Formation Mechanism and Rock Movement Pattern of the "Three Zones" in Single-Seam Mining
The results of the collapse and deformation of the overlying rock during the advance of the 29,211 mining face are shown in Figure 7. The 29,211 mining face has similarities with 29,210 in terms of the formation of "three belts" and the rock movement pattern. When 29,211 was advanced to 719 m from the foot of the slope, the basic roof first emerged with a height of 12 m. When the face was advanced to 644 m from the foot of the slope, the overburden damage developed to 75 m from the top of the coal seam, and the overburden was separated from the seam. When the mining face was advanced 594 m from the foot of the slope, the overburden rock fall developed to a height of 92.5 m from the top slab. When the mining face was advanced 544 m from the foot of the slope, the overburden rock damage reached a height of 105 m from the top slab. After the overburden destruction, a trapezoidal distribution is formed in the vertical direction, and the "three zones" of bending and sinking zone, fallout zone, and fracture zone are distributed. At 1 m from the reasonable stoping line, the bending and deformation of the rock layer developed on the surface.

Analysis of the Formation Mechanism and Rock Movement Pattern of the "Three Zones" in Single-Seam Mining after the Footing of the Slope
The overburden rockfall and deformation on the 29,209 mining face are shown in Figure 8. As the 29,209 mining face advances towards the foot of the slope, the pattern of overburden rockfall is similar to that of the 29,210 and 29,211 mining faces. From the physical model test results, when the mining face is advanced to 707 m from the foot of the slope, the first rise in the direct roof occurs, with a height of 5.5 m and a fallen angle of 52 • . Later, the direct roof rises with the mining, with a cycle pressure step of about 22 m. When the mining face is advanced to 656 m from the foot of the slope, the first rise in the basic roof occurs, with a fall height of 13 m. When the mining face is advanced to 227.5 m from the foot of the slope, the overlying rock layer of the coal seam falls. When the mining face was advanced to the foot of the slope, the coal seam overlying rock layer fell to a height of 80 m, and the overlying rock layer was separated from the seam. When the mining face was advanced to 555 m from the foot of the slope, the coal seam overlying rock layer fell to a height of 95 m. When the mining face was advanced to 454 m from the foot of the slope, the coal seam overlying rock layer fell to a height of 100 m, and when the mining face was advanced to 404 m from the foot of the slope, the ground surface was sunken and deformed. When advancing to the initial stoping line, the slope tends to slide, and a pressure square reinforcement plan should be adopted to ensure that mining reaches the stoping line without landslides.
The results show that during a similar simulation of the 29,209 mining face, the initial and periodic pressure on the overlying coal seam is obvious, as the distribution of the "three zones", namely, the bubble zone, the fracture zone, and the slow sinking zone in the mining area. In addition, the mining of the mining face will lead to surface deformation and even crack penetration to the surface. The mining line will have an impact on the stability of the slope.

Analysis of the Evolution of Slope Deformation in the Process of OPUG
The area of monitoring and calculation of the open-pit slope in the similar simulation of open-pit and shaft coordinated mining in profile 1 is shown in Figure 9.
The results of digital scatter photometry of the impact of mining face mining of the 4 # coal seam on the slope of the open pit are shown in Figure 10. When the 24,207 mining face was mined, the slope of the open pit was deformed to some extent, the maximum vertical displacement was about 0.5 m, and the corresponding strain value was small. The main reason for this analysis is that the mining face was far away from the slope location, so the deformation and movement of the overlying rock layer caused by the mining of the mining face had less impact on the rock body of the slope. Therefore, the mining of mining face 24,207 will not pose a hazard to the stability of the slope. When the 24,208 mining face advanced to the designed stoping line, the slope of the open pit was locally deformed and cracks appeared, with maximum vertical displacement and strain of approximately 1.0 m and 0.2, respectively, with vertical displacement and vertical strain mainly concentrated in the 1280 flat pan. The local deformation and cracks in the slope of the open pit do not affect the overall stability of the slope but are mainly due to the deformation and subsidence of the overlying rock layer after the mining of the mining face to form an arch structure, which has a large thrust in the direction of the slope face, resulting in some deformation and horizontal displacement of the slope. Figure 11 depicts the calculation area at mining face 29,210 in the 9 # coal seam.
In the model test, when the mining face advances to near the designed stoping line, the fractures in the rock body of the slope have penetrated, and at the same time, the fractures have extended to the flat plate of the lower slope. If mining continues, the fractures in the rock body are bound to penetrate (see Figure 12). Figure 13 shows the results of digital scattering photometry of the impact of mining the 29,210 mining face on the slope of the open pit. At this time, the maximum vertical displacement value is close to 1.8 m, the strain value increases to about 0.4, and the maximum horizontal displacement is larger in the direction of the proverbial surface. When the mining line continues to advance 60 m to 120 m after the designed stoping line, the maximum vertical displacement and strain value increase significantly. In addition, the horizontal displacement value is large in the direction of the proverbial surface. The crack penetrated to the foot of the slope, the upper slope rock is close to sliding out, and the whole slope may slide.

Analysis of the Evolution of Slope Deformation in the Process of OPUG
The area of monitoring and calculation of the open-pit slope in the similar simulation of open-pit and shaft coordinated mining in profile 1 is shown in Figure 9.  Figure 10. When the 24207 mining face was mined, the slope of the open pit was deformed to some extent, the maximum vertical displacement was about 0.5 m, and the corresponding strain value was small. The main reason for this analysis is that the mining face was far away from the slope location, so the deformation and movement of the overlying rock layer caused by the mining of the mining face had less impact on the rock body of the slope. Therefore, the mining of mining face 24207 will not pose a hazard to the stability of the slope. When the 24208 mining face advanced to the designed stoping line, the slope of the open pit was locally deformed and affect the overall stability of the slope but are mainly due to the deformation and subsidence of the overlying rock layer after the mining of the mining face to form an arch structure, which has a large thrust in the direction of the slope face, resulting in some deformation and horizontal displacement of the slope. Figure 11 depicts the calculation area at mining face 29210 in the 9 # coal seam.     In the model test, when the mining face advances to near the designed stoping line, the fractures in the rock body of the slope have penetrated, and at the same time, the fractures have extended to the flat plate of the lower slope. If mining continues, the fractures in the rock body are bound to penetrate (see Figure 12). Figure 13 shows the results of digital scattering photometry of the impact of mining the 29210 mining face on the slope of the open pit. At this time, the maximum vertical displacement value is close to 1.8 m, the strain value increases to about 0.4, and the maximum horizontal displacement is larger in the direction of the proverbial surface. When the mining line continues to advance 60 m to 120 m after the designed stoping line, the maximum vertical displacement and strain value increase significantly. In addition, the horizontal displacement value is large in the direction of the proverbial surface. The crack penetrated to the foot of the slope, the upper slope rock is close to sliding out, and the whole slope may slide.     Figure 15 shows the results of digital scatter photometry of the impact of mining on the slope of the open pit at 29211, where the overburden rock layer has fallen and sunk to form an arch-like structure, which has a horizontal thrust effect on the slope, resulting in the horizontal displacement of the slope deformation to the adjoining surface, thus producing a tensile effect on the slope and causing sliding. This will cause the slope to pull apart and deform in the settlement, which will lead to the sliding of the slope. When the    Figure 15 shows the results of digital scatter photometry of the impact o the slope of the open pit at 29211, where the overburden rock layer has fallen form an arch-like structure, which has a horizontal thrust effect on the slope, the horizontal displacement of the slope deformation to the adjoining surfac ducing a tensile effect on the slope and causing sliding. This will cause the s apart and deform in the settlement, which will lead to the sliding of the slope  Figure 15 shows the results of digital scatter photometry of the impact of mining on the slope of the open pit at 29,211, where the overburden rock layer has fallen and sunk to form an arch-like structure, which has a horizontal thrust effect on the slope, resulting in the horizontal displacement of the slope deformation to the adjoining surface, thus producing a tensile effect on the slope and causing sliding. This will cause the slope to pull apart and deform in the settlement, which will lead to the sliding of the slope. When the mining face advances to the vicinity of the reasonable stoping line, the maximum vertical displacement and strain are about 0.8 m and 0.3, respectively, and the displacement in the horizontal direction is larger towards the critical plane, with local rupture in the rock body of the slope, which has little impact on the stability of the slope at this time. However, when advancing to the vicinity of the designed stoping line, the internal cracks of the rock body of the slope gradually penetrate towards the bottom of the slope, and the maximum vertical displacement and strain are 1.5 m and 0.4, respectively, with the maximum vertical displacement and strain towards the critical plane. The horizontal displacement increased significantly in the direction of the critical plane. If the mining face continued to advance across the designed stoping line, the rock body of the slope would deform and rupture under the joint action of tension and pressure, and the slope would tend to slide out as a whole, which would seriously affect the stability of the rock body of the slope. body of the slope gradually penetrate towards the bottom of the slope, and the maximum vertical displacement and strain are 1.5 m and 0.4, respectively, with the maximum vertical displacement and strain towards the critical plane. The horizontal displacement increased significantly in the direction of the critical plane. If the mining face continued to advance across the designed stoping line, the rock body of the slope would deform and rupture under the joint action of tension and pressure, and the slope would tend to slide out as a whole, which would seriously affect the stability of the rock body of the slope.

Analysis of Coupling Effects in the Process of OPUG
The study of the spatial and temporal relationship between open-pit mining and shaft mining is a key issue in open-pit coordinated mining. The purpose of open-pit comining is to effectively utilize the complementary roles of open-pit mining and well-working mining while minimizing the coupling effects of both parties. In the process of OPUG, their mutual influence is mainly manifested in three aspects. (1) Shaft mining leads to the weakening of the sloping rock in its influence domain while activating the internal weak structural surface of the slope system, again disturbing the slope system which has been adjusted by the stress state due to the unloading effect of open-pit excavation, easily leading to the occurrence of slope instability disasters. (2) The unloading effect of opencast mining will also have an impact on shaft mining, especially the repeated loading of dynamic loads such as blasting in opencast mines, resulting in damage to the rock body and even loosening of the surrounding rock of the roadway. (3) Due to the overlap of the two influence domains, their mining effects are superimposed, coupled, and amplified, thus leading to disasters.

Analysis of Coupling Effects in the Process of OPUG
The study of the spatial and temporal relationship between open-pit mining and shaft mining is a key issue in open-pit coordinated mining. The purpose of open-pit co-mining is to effectively utilize the complementary roles of open-pit mining and well-working mining while minimizing the coupling effects of both parties. In the process of OPUG, their mutual influence is mainly manifested in three aspects. (1) Shaft mining leads to the weakening of the sloping rock in its influence domain while activating the internal weak structural surface of the slope system, again disturbing the slope system which has been adjusted by the stress state due to the unloading effect of open-pit excavation, easily leading to the occurrence of slope instability disasters. (2) The unloading effect of opencast mining will also have an impact on shaft mining, especially the repeated loading of dynamic loads such as blasting in opencast mines, resulting in damage to the rock body and even loosening of the surrounding rock of the roadway. (3) Due to the overlap of the two influence domains, their mining effects are superimposed, coupled, and amplified, thus leading to disasters.

Optimizing the Spatial and Temporal Relationship in the Process of OPUG
The spatial and temporal relationship between the mining operations of the Anjialing shaft mining face No. 2 mine and the opencast coal mine are two different coal mining methods in the same spatial and temporal context, and the spatial and temporal relationship between shaft mining face overlaid with the mining effects of opencast mining is a complex systemic issue. In terms of the stability and safety of the slopes of the opencast mine, the most favorable spatio-temporal relationship is for shaft mining to lag behind opencast mining, and then proceed with shaft mining when the inner row of the opencast mine is stable. However, due to succession problems, shaft mining currently overtakes opencast mining, and when shaft mining pushes to the designed stoping line, the inner row of the opencast mine cannot follow in time. This is why it is necessary to dynamically adjust the location of the stoping line and reduce the temporal and spatial parameters, such as the shaft mining intensity.

Optimizing the Mining Plan in the Process of OPUG
In the process of advancing the open pit to the east, the mining plan was adjusted to change from parallel advance to U-shaped mining, while accelerating the rate of internal drainage follow-up to ensure the safety of the open-pit slopes. During the mining process, the slope moves back towards the shaft mining void, and the backfill is used to stop the slope slippage trend. The backfilling process both prevents large tension cracks in the slope platform and increases the slip resistance of the slope, avoiding the possibility of landslides under the influence of well mining disturbances. For underground mining, pressurized footing backfill increases the overburden load on the quarry, which affects the stability of the mining tunnel and quarry envelope.

Monitoring and Early Warning Techniques
This project aims to maximize the recovery of coal resources under the premise of safe production in the surface mine. At the time when the 29,210 and 29,211 back mining faces are mined to the designed stoping line, the surface mine 1360 and above are in the shaft mining subsidence zone. The amount of settlement is less than that of 1405 and above the pan, and with effective slope maintenance, transport vehicles can pass normally. The Anjialing open-pit mine has a more advanced monitoring and surveillance system. During the joint mining process of the open-pit mine and shaft mining and after the shaft mining to the designed stoping line, when the internal row of the Anjialing Road open-pit mine has not been followed up, the slope monitoring work is the key to ensuring safe production. Therefore, the monitoring system is given full play, and the combination of full-coverage monitoring and early warning technology, such as geological radar monitoring and global navigation satellite system (GNSS) monitoring of the slope, is used to achieve prediction and forecasting of landslides on the slopes of the open-pit mines to ensure the safety and stability of the slopes under the influence of coordinated open-pit and shaft mining.

Conclusions
(1) The mining of the 24,207 mining face does not affect the stability of the slope. During the advancement of the 24,208 mining face, the distribution of "three belts", namely, the collapse zone, fracture zone, and slow sinking zone, is obvious in the vertical direction after the overlying rock layer emerges, but it does not pose a threat to the stability of the slope. The 29,209 mining face advanced to the initial stoping line. When the 29,209 mining face reached the initial stoping line, the rock deformation developed on the surface, obvious cracks appeared in the flat pan of the slope, and the slope tended to move in the direction of the airside. At this time, the slope reinforcement measures of internal pressure foot were adopted, so that the mining face continued to advance to the designed stoping line without causing sliding damage to the slope.
(2) As the 29,210 mining face continues to advance, the "three zones" of the overlying rock layer become more and more obvious. The angle of rock collapse is about 54 • , the height of collapse is about 22 m, and the cycle pressure step is about 19 m. As the mining face continues to advance 60-120 m, the fracture will soon penetrate to the bottom of the slope, and the upper step may slide out. The 29,210 mining face is suitable to stop mining when it advances to the reasonable stoping line 2 (341 m from the foot of the slope).
(3) As the mining face of 29,211 advanced, the rock collapse angle was about 56 • , the height of rock collapse was about 25 m, and the cycle collapse step was about 17 m. When the mining face advanced to 1360 flat, a fracture zone was formed under the slope. The fracture penetrated the whole slope rock body and the slope of the open-pit mine showed "situational landslide" damage. Therefore, the appropriate stoping position for the 29,211 mining face is when the mining face is advanced to the reasonable stoping line 2 (288 m from the foot of the slope).
(4) The digital scatter calculation results show that as the mining face of the mine continues to advance, the over-supporting stress of the mining face shifts to the side of the open-pit slope, and large displacement and stress concentration occur on the side of the slope near the stoping line. The deformation of the open-pit slope body first appears at the lowermost stratification of the slope. As the mining face advances, the overlying rock layer collapse height increases, and the deformation of the open-pit slope body continues to develop towards the upper part of the slope body. When the mining face gradually advances to the designed stoping line, the slope body is "pushed" along the uppermost part of the open-pit slope, and a fracture along the stratification layer appears below it. If the mining face continues to advance, the slope body may be "pushed" along the next level, and large displacements and stress concentrations may occur on the side of the slope near the stoping line.
(5) Full-coverage monitoring and early warning technology that combines the optimization of the spatial and temporal relationship, the mining plan in the process of OPUG, and the use of a combination of slope geological radar monitoring and GNSS monitoring has been adopted to achieve safe and stable control of slopes in the process of OPUG.