A New Mining Scheme for Hanging-Wall OreBody during the Transition from Open Pit to Underground Mining : A Numerical Study

School of Resources and Civil Engineering, Northeastern University, Shenyang 110819, China School of Manufacturing Science and Engineering, Southwest University of Science and Technology, Mianyang 621010, China State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology, Xuzhou 100083, China Shock and Vibration of Engineering Materials and Structures Key Laboratory of Sichuan Province, Southwest University of Science and Technology, Mianyang 621010, China School of Environment and Resource, Southwest University of Science and Technology, Mianyang 621010, China


Introduction
e hanging-wall ore-body is the ore-body which residues inside open-pit slope after the normal surface mining.Figure 1 shows a mode of the occurrence of hanging-wall ore-body at the Dagushan Iron Mine in China.According to statistical data, the residual ore-body outside the open-pit boundary accounts for 5-16% or even more of the total reserve until the open pit is closed [1].
In order to exploit sufficiently mineral resource and realize the smooth transition for the ore output, it is necessary to extract the hanging-wall ore-body during the transition from open pit to underground mining.However, it is highly challenging to operate simultaneously normal open-pit mining and mining for hanging-wall ore-body.Traditional mining methods for hanging-wall ore-body include surface mining and back-fill stoping method such as that in the Fenghuangshan Iron Mine of China and the Kidd Creek Zinc-copper Mine of Canada, the sublevel caving method such as that in the Yeshan Iron Mine and the Longshou Mine of China, and the sublevel open stoping and caving mining method such as that used by Liu et al. [1].Unfortunately, all these mining methods cannot terminate the serious mutual restriction between the hanging-wall orebody mining and the normal open-pit mining during the transition stage.For instance, when using the open stoping and back-fill stoping method to extract hanging-wall orebody, it is needed to limit the stope scale or adopt additional reinforcement of the open-pit slope, so that the mining efficiency is restricted; a large scale of slope instability could occur at the initial stage of mining when using the sublevel caving method to extract hanging-wall ore-body [2]; therefore, it has to limit the starting time of hanging-wall ore-body mining to ensure the safety of open-pit operation.Generally, these traditional mining methods fail to achieve large-scale mining of hanging-wall ore-body during the transition stage, which causes ore output decline for most mines in this period.Although there are a lot of studies on opencast-underground combined mining [3][4][5], it is rarely mentioned how to realize simultaneously the effective and safe operation of the open-pit mining and inner-slope mining during the transition stage.
In traditional mining methods, the slope usually needs to be stable when open-pit mining and inner-slope mining are operated simultaneously, which limits the stope scale and restricts the mining efficiency of the hanging-wall ore-body.In the present study, a new mining scheme for hanging-wall ore-body during the transition stage by employing a newly improved underground mining method, namely, the induced caving mining method, is proposed.In this scheme, the slope is allowed to collapse to a certain degree and the collapsed slope body is supposed to be absorbed by the mined-out area generated by the planned mining of the hanging-wall orebody. is mining scheme could avoid the influence of the inner-slope mining to the normal open-pit mining and thus guarantee the mining efficiency during the transition stage.Apparently, the slope instability behavior, which is greatly controlled by the rock mass structure of the slope and the exact mining plan (e.g., mining sequence), has a significant impact on the application of the proposed mining scheme.erefore, numerical simulations will be carried out to investigate this problem in the present study.
Numerical modeling is an important modern research approach, and it is expected to reproduce the displacement, stress, and strain distribution pattern of rock mass during mining thorough numerical simulations.Continuum-based numerical simulation methods, such as the finite difference method (FDM) and the finite element method (FEM), were widely applied to study the patterns of the movement and failure of the strata and surface induced by underground mining.For example, Zhao et al. [6] used the FDM software FLAC 3D to obtain the responding laws of the deflection and horizontal thrust of the roof rock beam under shallow mining conditions.Lin et al. [7] used the FLAC 3D to analyze the surface response and interaction of two layers of superthick coal seam mining.Li et al. [8] used the FLAC 3D to determine the minimum thickness of the safe mining roof of submarine gold deposits.Meanwhile, the FEM method was applied to assess the stability of debris and rock slopes [9], to determine the elastic compliance tensor of fractured rock masses [10], and to simulate rock cutting and its fragmentation process [11].It has certain defects to regard the rock mass as a continuous material when using the FDM and FEM to simulate mining problems.Discontinuous deformation and movement and large displacement are greatly involved in the surrounding rock for underground mining and in the slope rock for open-pit mining when slope failure takes place.
On the other hand, discontinuum-based numerical methods, such as the discrete element method (DEM) [12] and the discontinuous deformation analysis (DDA) [13] method, are by nature capable of simulating discontinuous and large displacement of rock masses.For example, the DEM-based software UDEC has been used to simulate the slope instability caused by the hanging-wall ore-body mining with the sublevel caving method [2].e DEM-based software PFC 3D has been used to analyze the failure process of soil-rock mixtures [14].
e DDA has been used to simulate the excavation of the tunnel and its reinforcement in columnar jointed basalt [15].In addition, some researchers used discrete element and finite element coupling methods to simulate large-scale underground excavation [16,17].Another unified continuous-discontinuous numerical method, namely, the numerical manifold method (NMM), has been used to analyze the stability and failure characteristics of the footwall slope [18].
In the present study, the DDA is adopted as the numerical method to do numerical simulations of the discontinuous deformation/movement and large displacement of the open-pit slope caused by inner-slope mining with the proposed mining scheme of hanging-wall ore-body.Key factors affecting the implementation of the mining scheme, such as the rock mass structure and the mining sequence, are deeply investigated through DDA simulations of instances based on practical engineering settings.
e proposed mining scheme is proved to be feasible, and some theoretical and technical guidance for the application of the mining scheme in practical mining operation are provided.

Induced Caving Mining Method.
e so-called induced caving mining method is a newly improved underground mining method as a combination of the block caving mining method (Figure 2(a)) and the sublevel caving mining method (Figure 2(b)) essentially.It contains the ore-body breakage process by underground pressure in block caving, and the drilling, blasting, and drawing process in sublevel caving.Figure 3 shows a typical induced caving mine layout. 2 Advances in Civil Engineering e ore-body is partitioned vertically downwards into three subareas with di erent mining characteristics, called the induced caving subarea, the sublevel caving subarea, and the bottom recovery subarea, respectively.e mining layer placed in the bottom of the induced caving subarea is regarded as the undercut level, and an extensive horizontal panel (undercut) beneath the mined block is formed by fan-shaped hole blasting.It should be noted that in the undercutting process, only about 30-50% of the volume of the blasted ore is extracted after each blasting, so as to provide a swell volume for the roof caving, while the remaining blasted ore temporarily stays at the undercut as a bu er layer to prevent the impact harm during the roof caving.Stress redistribution and gravity combine to trigger the roof progressive fracturing and caving into the undercut, and the caved ore is extracted in the mining process of sublevels below.Fracturing and caving extend progressively upwards as the caved ore is extracted, resulting in signi cant surface depression.In the sublevel caving subarea, usually at least two sublevels are placed to make sure a su cient extraction of the caved ore above.e mining operation in the sublevel caving subarea is exactly

Advances in Civil Engineering
the same as that in the sublevel caving mining method.e bottom recovery subarea provides the last chance for recovering the ore in the stope.erefore, dense production drifts in this mining subarea are employed to improve the ore recovery rate.It can be seen that the induced caving combines the advantages of large capacity, low cost of the block caving, and simple layout, flexible application of the sublevel caving.At present, the induced caving mining method is being applied to several mines in China, such as the Xiaowanggou Iron Mine, the Yanqianshan Iron Mine, the Beiminghe Iron Mine, and the Paishanlou Gold Mine, etc.

e Induced Caving Mining of Hanging-Wall Ore-Body.
Successful application cases of the induced caving mining method improve the understanding of the mechanism of roof caving and surface subsidence.On this basis, a new mining scheme for hanging-wall ore-body via induced caving during the transition from open pit to underground mining, as illustrated in Figure 4, is proposed to avoid the influence of inner-slope mining on the normal open-pit mining and balance the risks and benefits of hanging-wall ore-body mining in the transition stage.Due to the effects of weathering, surface blasting vibration, and mining unloading, etc., the hanging-wall orebody usually becomes favorable to be caved.Based on the preexisting development engineering of the open pit, it is convenient to implement the underground operation of the hanging-wall ore-body mining.As shown in Figure 4, an undercut level is placed within the ore-body to induce the upper ore to cave.While caving extends progressively upwards as the broken ore is extracted gradually, the slope will collapse when the caving reaches to a certain degree.It is expected that with a rational mining strategy, mainly including the design of mining span, height, and sequence, the occurrence time and the scale of the slope instability can be controlled, ensuring that there is enough mined-out area to absorb the collapsed slope body when the slope instability occurs.In such a way, the threat of slope instability to the normal open-pit mining can be eliminated, and the influence of the hanging-wall ore-body mining on the open-pit mining can be avoided.With perfect implementation, the proposed mining scheme can protect the normal open-pit mining from the threat of the slope instability throughout the whole transition stage.Moreover, some monitoring techniques, such as GB-InSAR, LiDAR, etc, can be used to forecast the slope instability to improve the safety of this mining scheme.In terms of mining strategy planning, both the block caving and sublevel caving can provide a lot of practical experiences [19,20].In the following, DDA simulations based on practical engineering settings will be carried out to examine the feasibility of the proposed mining scheme and the key influencing factors.

Engineering Background of Numerical Investigation
e Hainan Iron Mine is an open-pit mine located in Changjiang Li Autonomous County, Hainan, China.It is a large sedimentary metamorphic iron deposit.e ore-body is chiefly composed of hematite and its surrounding rock mainly consists of dolomite, phyllite, and quartz-schist.
e main ore-body situates at an elevation of −601.91 m to 202.44 m, with a length of 1670 m east-west, and a horizontal projection width of 350 m. e open pit with the closed loop at 168 m is being developed at an elevation of 24 m, and as planned, it will turn into underground mining below an elevation of 0 m. e south and north slopes of the open pit have hanging-wall ore-bodies below, which are primarily branches of the main ore-body and also individuals around.According to a series of geological surveys, the rock mass in the slope of the open pit mainly consists of two joint sets with good persistence.One set is approximately horizontal with evenly spaced bedding planes with a dip angle varying slightly with the survey sites (basically maintained within 0 °∼5 °) and the interlayer spacing is approximately 1.0 m, whereas the other set points to the foot of the north slope and inwards to the south slope with the dip angle and interlayer spacing varying significantly.e two main joint sets intersect to cut the rock into a jointed rock mass.Figure 5 shows the distribution of the geological survey sites and the sectional views of the south and north slopes.
e induced caving mining will be considered for the hanging-wall ore-bodies in the south and north slopes in Figure 5.According to the joint set orientations and some literatures of slope instability studies [21,22], due to the mining of the hanging-wall ore-bodies via induced caving, the north slope may take place sliding failure, whereas the south slope may take place topping failure.Determining slope instability under various mining conditions and the internal mechanisms and controlling factors is obviously impossible merely on a visual and empirical basis.In the next two sections, DDA simulation study on induced caving mining of the hanging-wall ore-bodies will be conducted, and two major factors, namely, the rock structure and the mining sequence will be considered.In the six simulation models constructed, the slope rock mass all contains two joint sets, marked as j 1 and j 2 , respectively.

Effect of Rock Mass Structure
e dip angle of j 1 maintains 0 °consistently, representing the approximate horizontal bedding plane in practice.e dip angle of j 2 is changed from 50 °to 150 °, and the value is taken every 20 °, representing the other major joint set with the dip angle varying with the survey sites.When the dip angles of the j 2 are 50 °, 70 °, and 90 °, the j 2 points to the toe of slope, corresponding to the north slope, while the j 2 are inclined at 110 °, 130 °, and 150 °, the j 2 points inwards the slope, corresponding to the south slope.Figure 6 shows the geometry of the model and locations of the measurement points in the simulation.In addition, modeling according to the actual joint spacing will result in 4 Advances in Civil Engineering excessive calculation due to the large number of rock blocks in the model.To reduce the calculation amount, the spacing of both joint sets in the model is enlarged with the same ratio based on that in practice and the number of blocks in each model is maintained around 3.0 × 10 3 .e exact number of rock blocks in the six models is shown in Table 1.Although such a simpli cation will reduce the reliability of the simulation results in certain aspects, e.g., the volumetric swelling of collapsed rock mass, it should have no decisive impact on the results of the slope instability, the main purpose of the simulation.
Because both joint sets have good persistence in practice, the two joint sets are set to penetrate the models.Moreover, the hanging-wall ore-body is divided into three sections in the models, with the top section representing the actual induced caving subarea and the middle and bottom sections representing the actual sublevel caving subarea.
e numerical method used to do the simulations is the DDA method.It is an implicit method with good convergence to compute the static and dynamic behaviors of discrete systems of deformable blocks [23]. is method has been successfully applied to solve many rock mass dynamics problems [24].In the DDA, each block can be of any shape and takes six variables, namely, , where the former three are for block displacements and the later three for block strains.e displacement increment in each time step at any point in the block is (uv) T Τ i D i , where Individual blocks are connected and form a block system by contacts between blocks and constraints on single blocks.Assuming there are n blocks in the de ned block system, the simultaneous equilibrium equations have the following form: where F i is the loading on block i distributed to the six deformation variables.Submatrix K ii depends on the material properties of block i and K ij where i ≠ j is de ned by the contacts between block i and block j. e physical and numerical control parameters used in the DDA simulations in this section are listed in Table 2. Generally, the physical parameters are based on some laboratory tests, such as density test, deformation test, compressive strength test, tensile strength test, and shear test, etc.
It is assumed that strength degradation will take place to a joint [25], so the cohesion, tensile strength, and the friction angle values of the joint are reduced after failure, as shown in Table 2. e numerical control parameters are determined based on extensive trial simulations to achieve high computational e ciency with a small number of open-close iterations in each time step.A good review of the openclose iteration in the DDA to treat contact between blocks with the penalty method can be found in [26].

Excavation Time and Sequence.
Using DDA static calculation, the initial stress equilibrium state of the model is obtained by applying gravity force to the blocks in the model.
ereafter, the excavation of the hanging-wall ore-body is executed, and the calculation turns to dynamic. Figure 7 shows the Y-direction stress evolution curves at measurement point MP1 (Figure 6) in the J1∼J6 models.It is shown that the six models reach initial stress equilibrium at t 20 s after approximate 20,000 time steps of DDA calculation.To  Advances in Civil Engineering simulate actual mining operations, the hanging-wall orebody is excavated by two steps, with the first step to excavate the induced caving subarea and the second step to excavate the sublevel caving subarea.e second excavation is performed at the end of the slope instability caused by the first excavation, and the whole calculation finishes when the slope instability caused by the second excavation completes.Figure 8 takes the J4 model (0 °/110 °) as an example to show the horizontal displacement evolution curves of the four measurement points MP2-MP5 locating on the top surface of the slope.In the J4 model, the first excavation is executed at the 20,000 time step (t � 19.999100 s), whereas the second  Advances in Civil Engineering excavation at the 570,000 time step (t 74.839971 s). e whole calculation nishes at 1,200,000 time step (t 137.759331 s).A series of trial calculation results show that the excavation time steps in the J4 model are also applicable to other ve models.is group of simulations simpli es the mining excavation into two steps to save computer time with a focus to investigate the e ect of the rock structure on the slope instability behavior in the proposed mining scheme of hanging-wall ore-body via the induced caving mining method.e impact of the mining sequence will be speci cally studied in Section 5.

Analysis of Simulation Results. Figures 9(a)-9(e)
show the slope states of the J1-J5 models at the 570,000 and 1,200,000 time step, which represent the slope instability morphologies after the rst and second excavation, respectively.Figure 9(f ) shows the slope states of the J6 model at the 100,000, 200,000, 570,000, and 1,200,000 time step, to demonstrate the whole failure process of the slope in more detail.

Slope Instability Mode.
Figure 9 shows that di erent slope instability modes are observed under various rock mass structures.e instability modes can be categorized into sliding and toppling as predicted in Section 3. e sliding instability appears in the J1∼J3 models, representing the conditions in the north slope, while the toppling instability appears in the J4∼J6 models, representing the conditions in the south slope in practice, respectively.
In the J1∼J3 models, joint set j 2 points to the toe of the slope.After the excavation of the hanging-wall ore-body, the resistance of the rock discontinuities becomes weaker than the sliding force of the rock under the slope free face; thus, the slope is destabilized and sliding failure takes place.In the J4∼J6 models, j 2 points inwards the slope.Once the hangingwall ore-body is excavated, the rock above in the slope will deform under the bending e ect and topple into the minedout area after the failure strength of the rock joints are reached.In the J1∼J5 models, the deformation and movement of the adjacent failed strata seems continuous.Comparatively, in the J6 model, the rock falls o as sparse discrete clusters into the mined-out area.e entire failure process involves the separation, falling, rolling, and bouncing of the rock mass.
It can be found that joint set j 2 plays a dominant role in the slope failure to control the failure mode.A dip angle smaller than 90 °is likely to induce sliding failures, while a value lager than 90 °is likely to induce toppling failures.As for the J6 model, in addition to the gentle inclination of the j 2 , the interlaced, long, and narrow rock blocks with sharp angles account for the special failure phenomena.e interlaced diamond-shaped blocks are bene cial for the slope to achieve stability during the toppling failure process.

Risk Assessment in Open-Pit Mining.
Using the mined-out area generated by the planned mining of the hanging-wall ore-body to absorb the collapsed slope rock mass is the core idea of the proposed mining scheme.e above DDA simulation results show that whether the collapsed rock will endanger normal open-pit mining is closely related to the slope rock mass structure.Figure 10 shows the nal pro les of collapsed slopes for the J1∼J6 models.
As for the J1∼J3 models, Figure 10 shows that in the J1 model, part of the collapsed blocks rush to the bottom of the open pit, threatening the open-pit mining below, whereas, the pro les of the collapsed slopes in the J2 and J3 models are obviously lower than the original slope pro les, which indicates that the mined-out area in these two models can su ciently absorb the collapsed rock; thus, the open-pit mining is safe.From these three models, it is concluded that a gentle inclination of the predominant joint set indicates that the collapsed rock mass is more likely to slide downwards to the bottom of the open pit to threaten open-pit mining.

Advances in Civil Engineering
For the J4∼J6 models, a large number of collapsed blocks rush to the bottom of the open pit in the J5 model.It is worth noting that although a bulgy slope profile is formed in the J5 model after the first excavation, there are no blocks rushing to the bottom of the open pit until the second excavation, as shown Figure 9(e).It shows that as the mining of the hanging-wall ore-body proceeds, the risk may increase.In the J4 model, although no bulk rock blocks rush towards the bottom of the open pit, the profile of the final collapsed slope has a high degree of coincidence with the original slope profile before hanging-wall ore-body mining.
e above conclusion from the J5 model indicates that there exist potential threats to the open-pit mining if further mining of the hanging-wall ore-body in model J4 is carried out.e profile of the collapsed slope in the J6 model is not as full as that of the J5 and J4 models; however, individual collapsed blocks bounced out to the bottom of the open pit, which also poses a threat to the open-pit mining.is is because the long and narrow blocks in model J6 fall off sporadically and get quite large kinematic energy after the hanging-wall orebody is excavated.
Additional measures, such as retaining dams and arresting barriers, can be used to bar rolling rocks.However, such practice is generally passive and unreliable in the case of a large amount of rolling stones.erefore, the present study proposes to eliminate the possible threat of the collapsed rock mass to the open-pit mining.
As indicated by the simulation results of the J1∼J6 models, the mode of slope instability due to induced caving of the hanging-wall ore-body, as well as the possible threat For a given geological condition of hanging-wall ore-body, the feasibility of the application of the proposed mining scheme can be expected to be achieved by adjusting the exact mining plan, e.g., optimizing the mining sequence, as investigated below in Section 5.

Influence of Mining Sequence
For the Hainan Iron Mine, mining the hanging-wall orebody under the south slope during the transition stage is more suitable because the open-pit haulage road is situated on the north slope.Figure 11 shows the production condition of the open pit.According to the previous numerical simulation results, using the induced caving mining method to extract hanging-wall ore-body under the south slope will trigger slope topping failure.Topping failure is highly likely to cause massive slope destruction; therefore, numerous collapsed rock blocks that rush from the slope to the bottom of the open pit may harm normal surface mining (as illustrated by Model J5).Here, the mining of hanging-wall ore-body under the south slope is taken as an example, and four simulation scenarios with di erent mining sequence are designed to explore the e ect of mining sequence on the development of slope instability and the mining safety, so as to determine a highly reliable mining sequence to guide the actual mining design.

Model Establishment.
e layout of the rock structure is provided in Figure 5(a).e ore-body under the slope has two joint sets, one of which is approximately horizontal and the other is approximately vertical, both with a spacing around 1.0 m.By combining the strength, mode of occurrence, depth, and location of hanging-wall ore-body with the actual mining experience, the initial mining level (undercut level) is tentatively set at an elevation of 60 m, with the second mining level at 45 m and the third mining level at 30 m. Figure 12(a) shows the layout of the stope structure on the A-A′ section (Figure 5(a)).Figures 12(b) and 12(c) illustrate the model geometry and the corresponding DDA model, respectively.Di erent from the abovementioned six models in Section 4, the present model not only sets the undercut level but also partitions the 2∼4 mining levels into independent blocks.Excavation of each independent block in the 2∼4 mining levels represents the actual mining of a production drift; thus, di erent mining sequences may be set within a mining layer.Table 3 lists the physical parameters of the ore-body.e physical parameters of the rock and the DDA numerical control parameters are the same as that in Table 2. e joint sets of the ore-body and rock mass in the model are also set as penetrating joints, and the joint spacing of the ore-body and rock is both increased by 5 times to 5.0 m, producing 3,072 blocks eventually in the DDA model.

Mining
Sequence.Similar to the previous simulation examples, solving the initial stress equilibrium of the model is necessary prior to excavating hanging-wall ore-body.Figure 13 plots the evolution curve of the stresses of measurement point MP1 (Figure 12(b)).It is shown that the model reaches the initial stress equilibrium after the 20,000 time step (t 19.86700 s); then, the undercut level is excavated.ereafter, the roof under the combined action of the stress redistribution and the gravity cave into the undercut, and the slope su ers topping instability simultaneously.However, the failure and movement of the slope rock mass are limited due to the supporting role of the caved ore, as shown in Figure 14(a).en, the caved ore is removed, representing the ore-drawing operation in practice.
e removal of the caved ore is accompanied by a further massive topping instability of the slope, as shown in Figure 14(b).At this time, a bulgy loose slope is formed.However, no bulk rock blocks rush in the bottom of the open pit, which indicates that setting the undercut level at 60 m is reasonable.
Based on the results of the J5 model above, the risk of bulk rock rushing downwards to the open pit from the slope may increase with the further mining progress.erefore, to seek for the following optimally safe mining plan, four scenarios with di erent mining sequences for the 2∼4 mining levels in the model, as shown in Figure 15, are considered.Scenario 1∼3 starts from the outside slope inwards, the inside slope outwards, and center to both sides, respectively.Scenario 4 has similar mining sequence with scenario 3; however, a safety ore-pillar of a certain size is reserved near the hanging-wall of the ore-body, and it is excavated in the end.e size of the safety ore-pillar is not designed based on its stability analysis but on the assumption that it does not fail during the mining.

Conclusions
In this study, a mining scheme for hanging-wall ore-body via the induced caving mining method during the transition from open pit to underground mining is proposed.e numerical simulation in this study proves the feasibility of the proposed mining scheme and explored its mechanism and some critical impacting factors including the rock structure and the mining sequence.
e numerical simulation results show that the mode of slope instability triggered by the induced caving mining of the hanging-wall ore-body is closely correlated with the rock mass joint orientations of the slope.When the predominant joint set enables sliding instability of the slope, a more gentle inclination of the predominant joint set means a higher possibility of the rushing out of the collapsed rock mass to the bottom of the open pit as a threat.Comparatively, when the predominant joint set hints topping instability of the slope, the relationship between the inclination of the predominant joint set and the threat of the collapsed slope rock to the open-pit mining becomes further complicated because certain joint set interactions turn to against the failure of the slope.Moreover, the mining sequence of hanging-wall ore-body has a signi cant impact on the development of slope instability and its threat to the open-pit mining.A reasonable mining sequence of the hanging-wall ore-body can help to e ectively eliminate the threat that slope instability poses to the normal open-pit mining.ese conclusions are instructive to the mining design of actual engineering.
Generally, the proposed mining scheme for hanging-wall ore-body via induced caving breaks the traditional concept that the slope must remain stable during the transition stage, which avoids the mutual restriction of hanging-wall orebody mining and normal open-pit mining, and helps to realize smooth transition for the ore output during the transition from open-pit to underground mining.is mining scheme is hopefully to be applied and validated in practical mining engineering in the future.
Data Availability e data used to support the ndings of this study are available from the corresponding author upon request.

Figure 1 :
Figure 1: A mode of the occurrence of hanging-wall ore-body at the Dagushan Iron Mine in China.

4. 1 .
Model Establishment and Simulation Scenarios.Based on the engineering background of Hainan iron mine, six numerical simulation instances are carried out to evaluate the significance of rock mass structure in the development of the slope instability and mining hazard assessment.

Figure 4 :
Figure 4: Illustration of mining scheme for hanging-wall ore-body via the induced caving mining method.(a) Mine layout.(b) Caving after undercutting.(c) Slope failure due to mining.

Figure 6 :
Figure 6: Geometric dimensions of the model and distributions of the measurement points (MPs).

Figure 10 :
Figure 10: Pro les of the collapsed slopes for the J1∼J6 models.

Figure 11 :Figure 12
Figure 11: Production condition of the open pit.(a) Aerial view, (b) south slope with the sectional view shown in Figure 4(b), and (c) broken outcrop.

Figure 12 :Figure 14 :Figure 13 :Figure 15 :
Figure 12: Stope layout and the DDA Model.(a) Layout of the stope structure on the A-A′ section.(b) Model geometry and distribution of measurement points.(c) DDA model.

Table 1 :
Rock joint and block number information in the six models.

Table 2 :
Mechanical and numerical control parameters of the models.