Study on the Process and Mechanism of Slope Failure Induced by Mining under Open Pit Slope : A Case Study from Yanqianshan Iron Mine , China

Mining under an open pit slope results in the collapse and slide of the slope. In this paper, a combination of methods including Google Earth and field investigations is applied to investigate the process of eastern slope failure induced by underground mining in the Yanqianshan IronMine over five years. According to the observed ground deformation features, the geomorphic zone of the eastern slope can be divided into four regions (caved rock zone, cracking zone, toppling zone, and sliding zone). Break angles and fracture initiation angles at different times are counted separately. Based on the above work, the process of initiation and development of slope failure is studied. +e analysis results show that the process of slope failure could be chronologically divided into three stages. First, a collapse pit, caused by the falling of the overlying strata above the goaf, appeared on the eastern slope. +en, the rock mass around the collapse pit slid into the pit to form a small landslide. Finally, because of mining disturbances and rock creep, a large landslide occurred on the northeastern phyllite slope. +e control mechanisms of each failure stage are discussed separately. Finally, the RFPA3D code is employed to simulate the slope failure process under the influence of underground mining. +e results are consistent with the field observations, which provided information on deformation failure and the mechanics of the slope that could not be directly observed in the field and deepened the mechanism analysis.


Introduction
In order to maximize mineral resource exploitation, the transition from open pit to underground mining is a feasible development for more open pit mines.During the transitional period, the residual orebody under the final slope should be exploited first.Mining under slopes will result in the movement and deformation of the rock stratum of pit slopes, which may further damage infrastructure and buildings on the surface.It not only threatens the safety of production but also endangers the environment.erefore, increasing attention has been paid to research on the process and mechanism of slope failure caused by underground mining.
Even to this day, surface deformation caused by underground mining is still one of the most challenging problems in mining engineering.Research on coal mining subsidence prediction has been performed for nearly two hundred years.According to a large number of measurements of mining subsidence, empirical formulas, curves, graphs, and influence functions [1,2] have been developed.However, most empirical and analytical methods for predicting subsidence are only effective for simple geometries, which are generally inadequate for complex surface geometries.
Compared with coal mining, the shape and occurrence of a metal orebody are more complex, so the geometric configuration of the mining-out area changes rapidly.Meanwhile, the geological structure of metal mines is more complex, and the surface movement is significantly restricted by geological features such as faults and joints [3].
e above studies mainly focus on mining under relatively flat surfaces.ere are limited cases in which mining under hillsides or slopes can trigger stratum collapses and landslides in these specific areas [4][5][6].
Advances in ground deformation monitoring contribute to the study of ground deformation.In recent decades, remote sensing technologies that have been successfully applied to surface monitoring include LiDAR laser scanning [7,8], photogrammetry [9], and differential interferometric synthetic aperture radar (DInSAR) [10][11][12].
Combining historical aerial photographs with digital photogrammetric modeling and point cloud analysis techniques as well as geomorphological mapping, Clayton et al. [13] studied the failure process and mechanism of the Mitchell Creek Landslide (MCL) in northern British Columbia.
Compared with the expensive and sophisticated monitoring techniques described above, Google Earth offers a more accessible way to observe surface damage by providing free three-dimensional high-resolution satellite images [14][15][16].
In contrast to the expensive high-tech observation methods, field investigation has the defect of low measurement accuracy, but long-term field investigation and records of the surface failure characteristics remain valuable data for studying surface collapse.According to long-term periodic mapping of ground cracks at the Kiirunavaara Mine, Villegas et al. summarized three different stages of surface deformation [17].
A numerical simulation can effectively analyze the stresses, strains, and displacement distributions of rock in the mining process.Some researchers have used numerical simulation methods based on continuous medium mechanics, such as the finite element method (FEM) and finite difference method (FDM), to study the movement and failure modes of strata and surfaces caused by underground mining [18][19][20][21][22][23][24][25][26].
However, different from conventional underground excavation [27], the mechanical behavior of strata in response to underground mining is very complex.On the one hand, roof caving and surface collapse caused by underground mining turn complete rock masses into loose blocks; therefore, the mechanical properties of rock masses before and after collapse are significantly different [28][29][30].Underground mining often leads to large deformations, discontinuous movements, and failure of the ground.Numerical simulation methods that rely on grids, such as FEMs and FDMs, will encounter the problem of grid deformation.
ese issues hinder the application of continuous medium numerical simulation in the study of surface failure caused by caving mining.
e distinct element method (DEM) has the intrinsic ability to simulate the discontinuous and large displacement of rock masses.Some researchers have used UDEC software, used to apply the DEM, to simulate roof failure and collapse during underground coal mining [31][32][33][34].Salmi et al. studied a coal mining-induced landslide in Nattai North, Australia, using the DEM-based software UDEC [6].Svartsjaern et al. utilized the DEM-based software PFC2D to simulate footwall rock mass damage induced by sublevel caving at the Kiirunavaara Mine [35].Using the DEM-based software 3DEC, Xu et al. simulated the influence of underground mining at the Yanqianshan Iron Mine on the eastern slope and predicted the damage range on the slope [36].However, when DEM is used to carry out numerical simulation, if the number of blocks in the model is too large, the computational efficiency will be very low.Meanwhile, there is a risk that the calculation results will not converge.
In addition, some researchers have used a finite element/discrete element approach to simulate the influence of large-scale underground mining on the failure of open slopes [29].
RFPA3D is a numerical software package that can simulate completely progressive failure and macroscopic mechanical behavior of rock materials.Readers are referred to more detailed information about the RFPA2D model by Li et al. [37][38][39].Liang et al. developed the 2D numerical code (RFPA2D) into a 3D code (RFPA3D) to investigate complex fracture growth in three dimensions [40].
In recent years, RFPA3D has been successfully applied to the mechanical analysis of engineering cases.Xu et al. used the creep calculation method of RFPA3D to study the mechanism of the Jiweishan landslide in the Sichuan province, China [41].Zhang et al. explored the mechanism of the southern slope landslide of the Fushun West Open pit coal mine in China [42].
Based on Google Earth satellite photos and records of in-site surface failure characteristics of the eastern slope in the Yanqianshan Iron Mine for more than four years, this paper analyzes the overall process of slope failure caused by underground mining, illustrates the failure mechanism, and finally simulates the failure process using RFPA3D code.e research area is located on the eastern slope of the Yanqianshan open pit, as shown in the red box in Figure 1.In this area, the main iron orebody is located in the middle, strikes about east-west, dips to the northeast at approximately 70-88 °, and has a width from 60 m to 200 m in the south-north direction.Carbonaceous phyllite and diorite are exposed in the north of the research area.ere are fragmented structures and strongly weathered zones in the carbonaceous phyllite.Diorite is on the northwestern part of the phyllite, and a faulted contact is found between the carbonaceous phyllite and the diorite.A wall-like branched 2

Geological Settings
Advances in Civil Engineering orebody, paralleling the strike of the main orebody, is inserted in the phyllite.e southern region of the study area is mainly distributed with granite mixed rocks.ere are three dominant groups of joints in each type of rock mass in the study area.e distributions of those three dominant sets of joints are quite similar in terms of occurrence, so there are three main sets of joints that can represent the structural surface distribution of the entire study area.e average orientations and dip angles of the three dominant sets of joints are 40 °∠75 °, 220 °∠20 °, and 300 °∠80 °, respectively [36].
After the end of open pit mining, the residual orebody under the eastern final slope began to be exploited.A sublevel caving method (Figure 2 e horizontal mining range of each layer is approximately 100 m to 220 m long from east to west and approximately 110 m to 130 m wide from south to north.e orebody is mined from the first to the fifth layer, and the mining direction for each layer is from east to west.
During the period from 2013 to 2017, underground mining was gradually promoted.By December 2017, the fifth layer had been mined.Continuous underground mining induced strata movement and caving, led to surface subsidence, and finally triggered landslides.e corresponding development state of underground mining and surface subsidence in each working period is described in detail in the next section.

Investigation and Analysis of the Geomorphological Evolution of the Eastern Slope Caused by Underground Mining
e field investigation and the Google Earth images are combined to identify and summarize the typical geomorphic deformation features.Detailed engineering geomorphological maps were made to record the surface deformation features on the eastern slope of the Yanqianshan Iron Mine for five years (2013-2017).It can make a fundamental contribution to further mechanism analysis.

Classification of Surface Deformation Features.
Figure 3 shows a series of visible deformation features (including scarps, counterscarps, and open cracks) emerging on the eastern slope.Based on field investigations, a classification standard of typical morphological deformation features was defined (Table 1).

Comprehensive Survey Method of the Geomorphological
Evolution: Google Earth and Field Investigations.Since the end of 2013, the surface of the eastern slope began to deform and generated a large number of cracks.e measurement of such cracks on-site represents a major challenge because of the large study area and varied topography.us, limited by labor costs and safety factors, it is very difficult to determine the morphology of the main cracks using only the artificial field measurements.
Google Earth (GE), launched by Google USA, provides a cost-effective and efficient way to survey surface Advances in Civil Engineering deformation.GE is a virtual global tool that documents the Earth's surface by integrating images from satellites, aerial photography, and geographic information systems (GISs) into a 3D global model.ree advantages of GE enable its use to conduct longterm accurate and cost-free monitoring of geomorphological changes over a wide range of areas [43,44].
(1) GE offers high-resolution satellite photos for free, 46-60 cm per pixel, with excellent spatial resolution (2) GE provides multiple layers for users, among which the topographic map layer can use DEM data collected by NASA's Space Shuttle Radar Topography Mission (SRTM) to observe topography at any angle in 3D space (3) GE has a time slider function, which adds a timecontrol function to the virtual earth view to provide satellite images of the same location at different times.
Because of the large extent of the open pit, the GE monitoring technique has its own advantages, as well as a number of shortcomings, as follows.

Advances in Civil Engineering
First, the poor weather in the area can affect the resolution of satellite images, resulting in faulty images.Second, although ground damage traces can be mapped by GE, the visible deformation features cannot be identified and classified reliably (scarps, counterscarps, or open cracks).erefore, it is still necessary to combine the GE and field investigations in our study to document the progress of surface deformation with underground mining.
e process of surface deformation as determined by the combined approach is shown in Figure 4.
(1) Since surface damage on the eastern slope caused by underground mining first appeared at the end of October 2013, the GE image taken in October 2012 was selected as the reference photo before the surface deformation.Take the surface deformation in October 2016 as an example.Based on the mark points and coordinates, the two GE images (2012-10 and 2016-10) were processed, and the plane images with the same area and the same proportion were output.After comparison, the surface deformation traces are identified and depicted in the treated plane image (2016-10).(2) We regularly conducted a field survey on the ground deformation features of the eastern slope at intervals of approximately 3-4 months.e main work includes locating the features by measuring the coordinates of a point on the trace with a portable GPS device, recording the trend and identifying the destruction type (scarps, counterscarps, and open cracks).After updating the GE satellite photos, a supplementary survey should be organized as soon as possible.According to the survey information before and after the satellite photos were taken and with the help of the virtual three-dimensional space provided by GE, the damage traces depicted on the satellite photos were reviewed to determine their accuracy and the types of the ground deformation features were identified and categorized (scarps, counterscarps, and open cracks).We found that the traces on the satellite photos could be verified by the surface survey, thereby demonstrating that the abovementioned techniques can complement and verify each other.
(3) According to the geomorphic deformation and underground mining in October 2016, the engineering geomorphological map illustrates the surface deformation caused by underground mining at the time.
Six GE images taken from 2013 to 2017 were selected to study the surface deformation process during underground mining.In this way, the geomorphological failures from 2013 to 2017 were plotted separately to study the evolution process of the surface deformation with the advance of mining.

Investigation Results of Geomorphological Deformation.
In this section, based on GE images, field investigations, and the classification of geomorphological deformation, engineering geomorphological maps, including surface failure features and the underground mining area, were plotted.
e geomorphological change process of the slope caused by continuous underground mining is described.

Investigation Results in November 2013.
e orebody under the eastern slope of the Yanqianshan Iron Mine began to be exploited in September 2012.By the beginning of November 2013, the −69 m level had been mined.en, a mined-out area, 136 m long in the east-west direction and 115 m wide in the south-north direction, was generated below the eastern slope (Figure 5(c)).
Despite the formation of the goaf, no subsidence was found on the eastern slope.Only one crack emerged along the western boundary of the horizontal projection of the goaf (Figure 5(b)).ere were some rolling stones on the slope.
Based on Figure 5, the engineering geomorphological map on November 8, 2013, was drawn (Figure 6).On February 26, 2014, the first ground collapse on the eastern slope occurred (Figure 7(d)).e failure appeared on the transport highway at an elevation of 16 m, and the rock mass suddenly fell down to form a well-shaped collapse pit. Figure 7(a), taken on March 12, 2014, clearly shows the topographic destruction on the eastern slope after the initial surface collapse.It caused a small landslide of fragmented rock around the collapse pit to roll into the pit.A series of tension cracks appeared at the back of the landslide head scarp (Figure 7(c)), and others developed along the dominant joint set (average occurrence orientation/dip:317 °∠41 °).It is obvious that the discrete rock blocks separated from the original rock mass by cracks tend to slide towards the collapse pit.

Investigation Results in
According to Figure 7, the engineering geomorphological map on March 12, 2013, was drawn (Figure 8).e survey results showed that, compared to March 2014, the range of the landslide increased significantly with the expansion of the collapse range (Figure 9).e head scarp of the small landslide developed up to +78 m.Surrounding tensile failures occurred at the back of the head scarp.

Investigation Results in
e northern boundary of the small landslide was defined by a side scarp trending in the east-west direction near the outcropped zone of the branched orebody.Some Advances in Civil Engineering tensile cracks, ranging in length from 9 m to 20 m, were distributed behind the side scarp on the northern phyllite slope.e trending directions of these cracks were basically consistent with the trending direction of the dominant joint set (average occurrence orientation/dip: 40 °∠75 °).ese cracks were controlled by the dominant structure planes.
e side scarp defining the southern boundary of the small landslide was consistent with the horizontal projection of the southern boundary of the goaf.In contrast to the phyllite in north, the rock mass quality of granite mixed rocks on the southern slope is better.Only an 18 m long tensile crack was found at +36 m.
According to Figure 9, the engineering geomorphological map on October 12, 2015, was drawn (Figure 10).Compared with the previous period (October 2015), the Google Earth image (Figure 11(a)) of February 2016 shows that the range extension of the small landslide was limited.
e landslide head scarp did not move upward, but the tension cracks in the rear of the scarp continued to develop.Existing cracks gradually penetrated, and new cracks were constantly produced.New tension cracks, also trending along the dominated joints (average occurrence orientation/ dip: 40 °∠75 °), emerged on the northern phyllite slope.is Step 1: by comparing the photos before and a er the ground deformation, the surface deformation traces are identified and mapped on the plane satellite photo  After eight months of mining, Figure 13(a) shows that the geomorphology of the eastern slope had changed significantly.Tensile cracks behind the small landslide head scarp developed and connected with each other to separate rock mass, resulting in the head scarp rising up to +124 m. e location of the northern boundary of the small landslide did not change.However, with phyllite cut by the cracks behind the northern side scarp bending and toppling downward into the collapse pit, the existing tensile cracks were transformed into counterscarps (Figure 13(c)).Simultaneously, new tension fissures were found in the rear of the counterscarps.
e influence of the creep characteristics of phyllite on the eastern slope was fully revealed in the last eight months.Although the horizontal projection of the mining area did not change significantly during this period, a large landslide occurred on the northern phyllite slope.e new large landslide head crown, located at the northeast corner of the open pit, was identified by an arcshaped scarp (Figure 13 According to Figure 13, the engineering geomorphological map on October 17, 2016, was plotted (Figure 14).(Figure 15(h)).It was clear that the sliding body of the large landslide began to slide southwestward.e sliding body still did not slip into the collapse pit at this time.However, many existing tensile cracks had developed into counterscarps, which can prove that the toppling failures of the phyllites intensified under the squeezing action of the landslide body (Figures 15(c) and 15(g)).
According to Figure 15, the engineering geomorphological map on December 27, 2017, was plotted (Figure 16).en, as the range of the collapse pit expanded, a small landslide gradually occurred around the pit.Finally, under the influence of underground mining disturbance and the creep deformation of phyllite, a large landslide, sliding southwestward, occurred on the northern phyllite slope.(4) e location and scope of the previous small landslide were controlled by the mined-out area.e location and scope of the later large landslide was controlled by the contact fault between the diorite and phyllite, the dominant joint sets, and the potential sliding surface.

Surface Failure Zoning.
e surface failure caused by underground mining is generally considered to be divided into three regions: caved rock zone, fractured zone, and continuous subsidence zone (Figure 17) [45].
e angles between the lines connecting the boundary of the deepest mining area to the boundaries of the three regions and horizontal lines are, respectively, identified as the angle of break, angle of fracture initiation, and angle of subsidence.
ese data are powerful tools for studying the mechanics of surface failure caused by underground mining.
Because the time of occurrence and the range and extent of the eastern slope failure are difficult to predict, it is impossible to arrange a long-term observation network on the slope due to safety concerns.On the contrary, the remote monitoring methods of surface subsidence (such as DIn-SAR) are costly and complicated.Considering storage safety and cost, large-scale subsidence monitoring has not been carried out on the eastern slope.As a result, the continuous subsidence zone cannot be identified.However, the information from the caved rock zone and fractured zone can still be obtained through Google Earth images and field investigations, which can also provide sufficient information for the study of slope failure mechanism [17].
To further study the failure mechanism, the fractured zone of the eastern slope was further subdivided into three subregions (cracking zone, dumping zone, and sliding zone) based on the investigation results of the geomorphic ere are no obvious linear failure characteristics, such as fissures, in the region.However, shear slip deformations, such as scarps and cracks, occur along the edges of the region.e rock mass in the region presents an overall sliding tendency.
Based on the engineering geomorphological maps in different periods (Figures 6, 8, 10, 12, 14, and 16), the variation of the surface failure partition on the eastern slope of the Yanqianshan Iron Mine is drawn to illustrate the failure process and explore the failure mechanism (Figure 18).e underground mining area lies on the south side of the east-west axis of the eastern slope.In the early stage of mining, the failure range can also maintain north-south symmetry (Figure 18(b)).However, since the end of 2015, the failure area of the northern phyllite slope started to be larger than that of the southern granite mixed rock.By the end of 2017, the farthest distance from the east-west axis of the mining area to the boundary of the failure zone on the northern phyllite slope was approximately 320 m, which is much larger than approximately 90 m of the southern granite slope (Figure 18(f )).On the contrary, the surface failure mode of the northern phyllite slope was more complicated than that of the south, and all three subzones of the failure appeared and evolved on the northern slope.
In this case, although the angle of subsidence cannot be obtained due to the lack of displacement monitoring data, the failure mechanism can still be studied through the break angles and fracture initiation angles [17].
According to the two sections (profiles I and II) in Figure 18(a), the changes of the break angle (Figure 19) and Advances in Civil Engineering the fracture initiation angle (Figure 20) of five periods (2014-3, 2015-10, 2016-2, 2016-10, and 2017-12) were statistically analyzed (Figure 21).e failure characteristics of granite mixed rock are similar to traditional surface collapse.After the collapse pit was formed, tension cracks appeared around the collapse pit.With the expansion and extension of the mining area, the break and fracture initiation angles increased accordingly.
Unlike the angle change rule in granite mixed rock, after dropping from 100 (March 2014) to 73 (February 2016), the break angle of the northern phyllite remained constant at 73 from February 2016 to December 2017.However, at the same time, the corresponding fracture initiation angle kept decreasing.
e following conclusions can be obtained according to the subdivision of the surface failure area and the break and fracture initiation angles: (1) Underground mining has the greatest impact on the northern phyllite slope, which eventually leads to a large landslide.(2) Before February 2016, the break angle is greater than 90 °. e reason for this is that the stratum rock above the goaf is arched with a cantilever structure.e 3.6.Analysis of the Mechanism of Slope Failure. Figure 22 shows the failure mechanism of the eastern slope as follows: (c) ere is a special mechanics structure on the northern slope, in which the wall-like branched orebody inserted into the phyllite is similar to a retaining wall. is structure prevents the collapse from propagating north further.However, due to the disturbance caused by continuous mining and the creep characteristics of phyllite, the superficial phyllite has a tendency to slide and topple.e expansion of the underground mining range and depth causes constant disturbance to the slope.e phyllite rock blocks cut by tension cracks begin to bend and topple into the collapse pit, and the range of cracking and toppling is constantly expanding.At the same time, the potential sliding surface of the shallow phyllite is continuously developed, which extends to the northeast and produces shear failure away from the collapse pit.(d) Under the comprehensive influence of continuous underground mining disturbance, creep deformation of rock mass, and pushing of phyllite, the branch orebody is damaged.is means that the retaining wall near the surface fails.
e phyllite moves downward along the potential sliding surface, resulting in a large landslide.

A Brief Introduction to RFPA3D.
e RFPA3D model is based on the following assumptions: (1) e rock mass is assumed to be elastic-brittle with a residual strength at the elemental scale.e mechanical behavior of the rock is subject to an elastic damage constitutive law.
e residual strain/ deformation of the rock mass under unloading is not considered.
(2) If the minimum stress of an element exceeds the tensile strength, the element is considered to fail in the tensile mode.When the shear stress of an It is well known that acoustic emission (AE), generated by rapid microcrack growth, is a particular phenomenon correlated with the release of elastic energy produced by brittle fracture within the deformation of brittle materials.It provided abundant information in relation to the failure process in heterogeneous brittle rock.Monitoring AE or microseismic events generated during brittle material deformation is an effective way to detect the origination and the evolution of cracks [46,47].In engineering cases, such as underground mining and landslides, AE or microseismic events are good indicators for assessing rock mass failure.A single AE event representing a microcracking event can indirectly estimate the damage evolution.erefore, the AE counts are accounted for by the number of damaged elements, and changes in the number of AE events can reflect the degree of material destruction.

Numerical Modeling.
Because of the imperfection of the preprocessing function of RFPA3D, a modeling method presented by Hu et al. is adopted to build a more accurate geological numerical model.Figure 23 shows the modeling process.First, the specialized mining software package SURPAC is used to establish a geological solid model, and then the solid model is discretized into a block model  Advances in Civil Engineering composed of hexahedral elements.Subsequently, the information of the block model is dealt with by the conversion code named stor.mprogrammed in MATLAB to generate data files that can be read by RFPA3D.Finally, the data files are imported to the RFPA3D code to build the RFPA3D calculate model.Based on the location and boundary of the mining area and the open pit, it is determined that a reasonable numerical simulation research area is 1.5 to 2 times the size of the largest mining area.As shown in Figure 1, the region is 900 m in the east-west direction, 1100 m in the north-south direction, and approximately 500 m in the vertical direction.
e model contains four types of rock masses, including orebody, phyllite, mixed granite, and diorite.e macromechanical parameters of the rock masses are listed in Table 2.
e model consists of 8,97,471 hexahedral elements and 9,40,713 nodes.e side boundaries of the calculation domain are restrained in the normal direction, the bottom boundary is fixed, and the surface is free.e creep calculation method [41] was used to simulate the failure of the eastern slope.In the simulation, orebodies in three mining zones are mined to simulate the mining activities that took place from 2013 to 2017 in three stages of mining excavation, i.e., 2013∼March 2014, March 2014∼February 2016, and  After the first stage of mining excavation (2013-March 2014), the associated stress field in the rock mass was redistributed under the combined influence of the mining excavation and gravity loading.
e high tension stress concentration, corresponding to the green part in the tensile stress evolution in Figure 25(a), occurred in the overlying rock of the mined-out area.At the same time, there was ground surface subsidence on the slope (Figure 26(a)).
e tensile stress along the surface projection boundary of the goaf grew constantly, and finally, the surface rock mass cracked under tensile force (Figures 24(b   After the third stage of mining excavation (February 2016 ∼ December 2017), cracks appeared between the phyllite and diorite rock mass on the northern slope, and these cracks eventually developed to form a continuous fracture zone.On the contrary, due to creep deformation, the engineering quality of rock masses continually declined.Finally, under the action of gravity, a slip surface of the phyllite slope was formed by shear slip, and the phyllite rock mass, which was separated by the fracture zone and the slip surface, slid southwestward (Figures 24(f           Advances in Civil Engineering e simulated failure pattern is in good agreement with the observed in situ sliding failure pattern, as shown in Figure 26. Figure 27 shows the numerically simulated AE events or associated damage event count sequences during mining.It can be shown that there are three peak periods of AE events, respectively, in steps 32∼36, 50∼54, and 87∼90.Each peak corresponds to the large-scale destruction after each excavation, which proves the driving effect of the excavation under the slope on the failure of the slope and verifies the reliability of the numerical simulation.

Conclusions
In this paper, Google Earth images and field investigations were used to analyze the process of fracture propagation on the eastern slope under the influence of underground mining.Based on the site investigations, the failure mechanism and the damage process of the eastern slope were investigated by numerical models.According to the above comprehensive research method, the following conclusions can be obtained:  5) e location and scope of the previous small landslide were controlled by the mined-out area.e location and scope of the later large landslide was controlled by the contact fault between the diorite and phyllite, the dominant joint sets of the phyllite, and the potential sliding surface.(6) RFPA3D is used to simulate the process of slope failure with underground mining.e stress analysis reveals that the slope failure is mainly caused by tensile failure and shear slip.Displacement analysis and AE event statistics establish three stages of slope failure: collapse, a small landslide caused by the collapse, and a large landslide on the northern phyllite slope.e simulation results are consistent with the field observations.
(a)), with horizontal slice mining from top to bottom, was employed.e orebody beneath the eastern slope is divided from top to bottom into five horizontal layers: the elevation ranges for the first to the fifth layer are −51 m to −69 m, −69 m to −87 m, −87 m to −105 m, −105 m to −123 m, and −123 m to −141, respectively (Figure 2(b)).

Figure 1 :
Figure 1: e location and geological model of the Yanqianshan Iron Mine.

Figure 2 :
Figure 2: Illustration of (a) the typical sublevel caving method (image copyright: Atlas Copco) and (b) the sublevel caving under the slope.

Figure 3 :
Figure 3: Illustration of the classification of surface deformation features.
March 2014.By the end of February 2014, the mining of the −69 m level had been completed, and the −87 m level was being mined.e goaf, 145 m long in the east-west direction and 120 m wide in the south-north direction, is shown in Figure 7(b).
October 2015.By the beginning of October 2015, the −69 m and −87 m levels had been mined out, and the mining of the −105 m level had started.e horizontal projection of the mining area was 180 m long from east to west and 130 m wide from north to south (Figure 9(b)).

3. 3 . 4 .
Investigation Results in February 2016.By the end of February 2016, the −105 m and −123 m levels were in production.e horizontal plan of the mining area was 200 m long in the east-west direction and 130 m wide in the south-north direction (Figure 11(b)).
based on the on-site survey information, with the help of the 3D virtual space provided by GE, the traces mapped on the plane image are reviewed and classified (scarp, counterscarp and open crack) e trace lines of ground deformation features Step 3: according to the situation of ground deformation and underground mining at this time, the engineering geomorphologic map is made 2016-10 N Six GE images are taken from 2013 to 2017 (the shooting dates are marked with red flags)

3. 3 . 5 .
Investigation Results in October 2016.By October 2016, the −105 m and −123 m levels were being mined, and the horizontal projection of the mining area remained unchanged (Figure 13(b)).
(d)) that developed from open cracks and slip failures (Figures 11(c)-11(e)) six months ago.e western boundary of the new large landslide was identified by an approximately 170 m long side scarp trending in the northeast-southwest direction.Coincidentally, the boundary coincided with the trace of the faulted contact between phyllites and diorites.e eastern boundary of the large landslide was made up of a series of discontinuous cracks that were 8-10 m long and trended along the dominated joints (300 °∠80 °) (Figures 13(e) and 13(f )).

Figure 7 :
Figure 7: e survey of geomorphological changes on the eastern slope of the Yanqianshan Iron Mine (March 2014).(a) 3D view from Google Earth (12/3/2014).(b) e location of the mining area, which is observed through the same viewing angle as that of the Google Earth image (March 2014).(c) Tension cracks on the back of the small landslide head scarp.(d) e initial collapse pit (26/2/2014).

8
Advances in Civil Engineering 3.3.6.Investigation Results in December 2017.By the end of 2017, the main orebody had been mined, including each orebody of the −51, −69, −87, −105, and −123 m levels, and the −141 m level was being mined.e horizontal projection of the mined-out area was approximately 220 m long in the east-west direction and approximately 132 m wide in the south-north direction (Figure 15(b)).

Figure 15 (
a) shows the further development of the new large landslide.e head scarp of the landslide was still located at +150 m in the northeast corner of the open pit.Discontinuous fractures along the eastern boundary developed into continuous scarps (Figures 15(e) and 15(f )), and the scarps were connected with the southern boundary of the old small landslide.Despite the lack of underground borehole data, strong evidence inferred the existence of a sliding surface below the northern phyllite slope because of the cracking (Figure 15(d)) and sliding near the head crown, shear failure along the eastern and western boundaries of the large landslide, and the main set of existing joints (220 °∠20 °)

Figure 9 :
Figure 9: e survey of geomorphological changes of the eastern slope of the Yanqianshan Iron Mine (October 2015).(a) 3D view from Google Earth (12/10/2015).(b) e location of the mining area, which is observed through the same viewing angle as that of the Google Earth image (October 2015).(c) e small landslide on the eastern slope (east view).(d) e small landslide on the eastern slope.

Figure 10 :
Figure 10: e engineering geomorphological map and mining progress on October 12, 2015.(a) Distribution of ground destruction.(b) At profile I. (c) At profile II.
stratum collapse process caused by underground mining is the progress of repeated failure and reformation of the stress balance arch in the overburden strata from a mechanical point of view.(3) With the extension of the mining area, the fractured zone developed northward.However, the northward extension of the caved rock zone stopped near the wall-like branched orebody.It is different from the traditional mechanics of surface subsidence in which the caved rock zone and the fractured zone expand together.e branched orebody intruded into the northern phyllite and acted as a retaining wall, which prevented the deep phyllite from sliding towards the collapse pit.Under the combined influence of mining disturbance, rock mass creep, and the pushing of phyllite at the back of the branched orebody, the part of the branched orebody near the surface gradually undermined and eventually failed.It brought about cracking and toppling in the phyllite slope.

Figure 11 :
Figure 11: e survey of geomorphological changes of the eastern slope of the Yanqianshan Iron Mine (February 2016).(a) 3D view from Google Earth (20/2/2016).(b) e location of the mining area, which is observed through the same viewing angle as that of the Google Earth image (February 2016).(c) Scarp away from the small landslide.(d) Tension cracks away from the small landslide.(e) Tension cracks behind the side scarp (northern boundary) of the small landslide.

Figure 12 :
Figure 12: e engineering geomorphological map and mining progress on February 20, 2016.(a) Distribution of ground destruction.(b) At profile I. (c) At profile II.

Figure 13 :
Figure 13: e survey of geomorphological changes of the eastern slope of the Yanqianshan Iron Mine (October 2016).(a) 3D view from Google Earth (17/10/2016).(b) e location of the mining area, which is observed through the same viewing angle as that of the Google Earth image (October 2016).(c) Counterscarps showing rock toppling failures on the phyllite slope.(d) Scarps away from the small landslide.(e, f ) Tension cracks along the eastern boundary of the large landslide.

February
2016∼December 2017.e three stages of mining excavation were arranged successively at Step 2, Step 40, and Step 80.

4. 3 .
Numerical Results.Figures 24 and 25 present the compressive and tensile stress evolution of the rock mass; the corresponding deformation displacements in the z-direction of the rock mass are shown in Figure 26.

Figure 14 :
Figure 14: e engineering geomorphological map and mining progress on October 17, 2016.(a) Distribution of ground destruction.(b) At profile I. (c) At profile II.

Figure 15 :
Figure 15: Survey of geomorphological changes of the eastern slope of the Yanqianshan Iron Mine (December 2017).(a) 3D view from Google Earth (27/12/2017).(b) e location of the mining area, which is observed through the same viewing angle as that of the Google Earth image (December 2017).(c) Tension cracks north of the counterscarps.(d) Tension cracks near the large landslide head crown (August 2017).(e, f ) e eastern boundary of the large landslide (August 2017, October 2017).(g) Counterscarps behind the northern boundary of the small landslide.(h) e large landslide (view east).

16
Advances in Civil Engineering rolled and slipped into the pit with the mined-out area enlarging, and the collapse pit expanded outwards (Figures 24(d), 24(e), 25(d), and 25(e)).

Figure 16 :
Figure 16: e engineering geomorphological map and mining progress on December 27, 2017.(a) Distribution of ground destruction.(b) At profile I. (c) At profile II.

Figure 20 :Figure 21 :Figure 22
Figure 20: e change in fracture initiation angles over five periods.(a) At profile I. (b) At profile II.

Figure 23 :
Figure 23: e modeling process.(a) e solid model in SURPAC, (b) the block model in SURPAC, and (c) the RFPA3D model.

Figure 22 :
Figure 22: An illustration of the failure mechanism of the eastern slope.

( 1 )
e combined survey method of Google Earth and field investigations can effectively implement longterm and low-cost monitoring of surface deformation.(2) ree typical characteristics of surface deformation were found in the failure process of the eastern slope of the Yanqianshan Iron Mine: open cracks, scarps, and counterscarps.Open cracks, developed along the existing dominant structure planes, were mainly distributed on the phyllite slope north of the mining area.e open cracks divided the phyllite rock mass into discrete subvertical layered rocks.With layered rocks bending and toppling, open cracks developed into counterscarps.Scarps are divided into two categories: some around the mining area caused by ground collapse and the others far away from the collapsed open pit caused by massive rock mass sliding.(3) e process of eastern slope failure can be divided into three stages.In the first stage, roof failure of the mined-out area brought about surface collapse.In the second, it resulted in a small landslide in which the rock mass around the collapse pit slid and fell into the open pit.Finally, mining disturbance and creep deformation of the phyllite slope caused a large landslide on the northern phyllite slope.(4) Based on the characteristics of ground deformation, the destruction area of the eastern slope can be divided into four zones: caved rock area, cracking zone, toppling zone, and sliding zone.Combined with the analysis of break and fracture initiation angles, it is shown that the wall-like branched orebody, which is inserted into phyllite, acts as a retaining wall in the transition stage from the early small landslide to the final large landslide.After the collapse pit was formed, the wall-like orebody restricted the southward displacement of phyllite to the north of the mining area.With the disturbance of underground mining and the pushing of the potential sliding body on the northern phyllite slope, the retaining wall-like structure gradually trended towards instability.It resulted in a dynamic and cyclic change pattern of geomorphic destruction features on the phyllite slope: first, open cracks were formed; then, the existing open cracks developed into counterscarps; finally, new open cracks were generated which could become new counterscarps.(

e
Yanqianshan Iron Mine, which was first constructed in 1960, is located in Anshan city, Liaoning Province, China.e open pit mining was finished in September 2012.After more than half a century of mining, a huge open pit had formed, which was 1400 m long from east to west and 600 m to 700 m wide from north to south.e elevation of the closed loop of the pit is approximately 90 m and that of the ultimate bottom is −183 m. e southern, northern, western, and eastern slope angles are approximately 35 to 38 °, 32 to 35 °, 30 °, and 43 °, respectively.

Table 1 :
Description of the deformation feature classification.

Table 2 :
e macromechanical parameters of the rock masses.