Effect of Shaft Pillar Extraction on Stability of Main Shaft: A Case Study at Xincheng Gold Mine, China

Mining of ore body in the vicinity of a shaft has a significant influence on its stability. The in situ monitoring and numerical simulation are employed to analyze the effect of shaft pillar extraction at Xincheng Gold Mine. The XI# ore body is recently found around and beneath the shaft, and mining in this area may be detrimental to the shaft. Firstly, on the base of geological survey and in situ displacement monitoring, mechanical parameters of rock mass are obtained and the displacement around the shaft is measured. Secondly, the sensitivities of fivemain factors that may affect the shaft displacement are analyzed bymeans of orthogonal experiment according to the numerical simulation with FLAC3D. Finally, a numerical model is established according to the in situ condition; in order to forecast the shaft displacement induced bymining activities of XI# orebody, theMining Priority Index (MPI) is put forward and used to select the optimal mining sequence. Based on the comparison between the numerical results with the monitoring data, it is determined that the ore within 100m from the shaft is not suggested to be extracted until the last period of the shaft life.


Introduction
Shafts are vital access ways for many underground mines, providing passageways for personnel, materials, and ventilation air.The shaft must be designed as a long-term workable and operational structure.As we know, a series of adverse factors threaten the stability of the shaft, such as faults, high geological stress, and mining activities.The conventional method of protecting vertical shafts and their associated service excavations from the mining disturbance is to leave a sufficient large pillar in the vicinity of the shaft.However, with the mining operation going deeper, sometimes fresh ore bodies may be found near the shaft; this is the case that happened at Xincheng Gold Mine, China.In this regard, the stability of shaft will be mainly concerned during the mining operations when the ore body near the shaft is going to be mined out.Furthermore, delineating a shaft pillar in conformity to a certain limit angle is the most common approach to be adopted for this issue.As a result, the amount of mineral resources cannot be mined out around the shaft, and the unworkable mineral will increase with the depth of the shaft.Therefore, this method is not applicable when the shaft depth is greater than 2.5 to 3.0 km in South Africa [1].
Recently, the downturn of the mining commodities price forces mining companies to extend the possible productivities, such as the ore body in the shaft pillar.As for the Xincheng Gold Mine in Shandong Gold Group, China, in recent years, the XI# ore body is found around the main production shaft that is currently used to hoist minerals from the V# ore body, Figure 1.In this regard, it is essential to keep the shaft workable and operational for a long time if the XI# ore body is intended to be mined out.
Many studies have been conducted on the mining operation of shaft pillars and a great number of experiences have been accumulated.Effects of the extraction of tabular deposits around vertical shafts in deep-level mines were numerically studied by Budavari and Croeser [1]; however, due to the limitation of the computing power and the difficulty in quantifying the rock mass parameters and specifying in situ stress conditions, the calculation results usually deviate from the in situ reality to some degree.When extracting the reef around the shaft at the start of mining operations, several  alternatives to protect shaft pillars for deep vertical shaft systems were put forward [2], including the use of backfill, suitable mining layouts, and delaying shaft steelwork and lining installation until residual deformations are tolerable.A combination of strategically placed satellite pillars, yield pillars, and backfill method was chosen to protect the shaft from possible damaging stress changes and to minimize displacement in the shaft-reef area at Beatrix number 3 vertical shaft [3].In Sumitomo Akabira Colliery, changing the displacement angle of the shaft pillar from 60 ∘ to 80 ∘ , the shaft damage was limited when the main fault crossing the shaft was considered as special attention [4].Mining of the Ross pillar at Homestake Mine indicated that unexpected shaft wall movement could occur at a distance of 450 m from the current mining area [5].A waste-cut across the shaft was proposed to extract the entire shaft pillar at the number 4 shaft of the African Rainbow Minerals (ARM) mine, which is located in highly stratified and poor-quality rock masses [6].A shaft destressing slot is excavated to protect the main production shaft at Newmont Canada's Golden Giant Mine from potential damage from the high mininginduced stresses, and it allows the mine to extract a lot of high grade ore [7].Degradation occurs to the X41 shaft at Copper Mine because of two major geological structures, the W41 and W42 faults, which intersect the shaft at two distinct locations [8,9].By virtue of the high geotechnic stress, serious deformation occurs during the sinking of a main shaft at Jinchuan Mine number 3; a tolerable deformation of 95 mm before lining is reasonable in order to maintain the shaft stability, as suggested by Zhao et al. [10].In addition, it is noted that the dynamic interaction between the conveyances and the conveyance guiding system of aging mine shafts should be analyzed in order to examine the workability of shaft [11].At Chengchao Iron Ore in China, the shaft exhibited evident degradation because of mining activities, Huang et al. [12] thought that the east shaft can continue to be used despite the surface subsidence rate of 0.314 mm/d and 0.144 mm/d in horizontal and vertical directions, respectively.Besides, some researches presented that lengthening or shortening the shaft lining by no more than 3 mm/m or 1 mm/m, respectively, generally does not cause sufficient damage to threaten the safety of shaft operations [13], and the shaft vertical displacement of 35 mm can be used as a threshold for lining damage of the X41 shaft at Copper Mine [8].
All of the case studies and above conclusions indicate that the extraction of shaft pillar is feasible if the sufficient measures such as backfilling, suitable mining layouts, and displacement monitoring are taken to mitigate the effect of mining operation on shaft displacement, even though in some cases shaft degradation occurs because of geological structures that are not clearly understood.
The theoretical analysis, in situ monitoring, and numerical methods are generally employed to examine the mininginduced responses to the shaft.Specifically, the theoretical model is of poor generality to represent the in situ condition of underground mining.The in situ monitoring could provide a lot of data, but they are sometimes difficult to be interpreted because dozens of uncertainties exist.On the other hand, the numerical simulation could give the prediction by using sophisticated numerical models; however, the rock mass parameters and in situ stress data are difficult to ensure accuracy and reliability.Therefore, the issue of shaft displacement and damage is complex, which should be solved by combining the performance of different methods.The solution strategy in this paper includes geological survey, in situ monitoring, and numerical simulation, hoping to achieve a better prediction of shaft deformation and to provide a reasonable suggestion to the mining operations.

Description of the Mine Site
The Xincheng Gold Mine is operated by Shandong Gold Group Co., Ltd, located in Shandong Province, China,   The shaft at Xincheng Gold Mine is 499 m deep, and the ground surface elevation is +33 m.According to the latest exploration near the main shaft, a new ore body numbered XI# that consists of 7 ore veins at −380 m level and 2 ore veins at −530 m level are found, Figure 1.Distribution of the XI# ore veins is scattered, the average strike is NE57 ∘ , and the dip of the veins is in the 25-30 ∘ range, which is parallel to the stratigraphic dip, and detailed information of the ore veins is listed in Table 1.In the beginning, a shaft pillar (Figure 3) has been delineated according to the respective rock movement angle of different strata, Table 2.Under these circumstances, more than 1140 kt ore cannot be mined out, which accounts for 67.4% of the XI# ore body.The first mining areas are designed in the XI-1 ore vein at −380 m level and in the XI-8 and XI-9 ore veins at −530 m level, and 280 kt ores have been mined out up to March, 2017.In order to protect the shaft, the following stopping sequence needs to be optimized and the stopping boundary needs to be determined.
The upward horizontal cut and backfill stopping method has been adopted for many years at Xincheng Gold Mine, and a wealth of experience has been accumulated.Compared with the open stopping method and stope caving method, the stopping-and-filling method is better in terms of maintaining the surrounding rock mass stable and preventing the surface subsidence.Therefore, this method is also adopted to excavate the XI# ore body in field production.

Geological Survey
The detail rock mechanical parameters of the Xincheng Gold Mine were obtained during the shaft sinking which are listed in Table 2.The thickness of the quaternary surface soil is thin, while the bedrock is mainly magmatic rock and metamorphic rock with high mechanical strength, which are relatively intact, stable.Besides, the major geological structures in this area include the Jiaojia Fault, the Houjia Fault, and the Hexi Fault, as shown in the 155# vertical cross-section map, Figure 4.The Hexi Fault especially is close to the bottom of the shaft and overlapped with the ore veins XI-8 and XI-9 at −530 m level; thus, it needs to be given special attention.
In fact, the three faults in this area are crushed zones with thickness of dozens of meters, rather than the general faults with the main fracture surfaces.Former studies have shown that the horizontal geotechnic stress is significant at Xincheng Gold Mine [14], the  max and  min are both horizontal at the −280 m level and below, and the direction of  max is perpendicular to the strike of the strata, while the direction of  min is parallel to it contrastively.So, it can be further inferred that the faults in this area tend to be reverse faults, and this deduction provides a foundation for the follow-up analysis of the shaft stability.Thickness of the shaft liner is 0.3 m, and the net diameter is 5 m.Inside the shaft are a cage compartment, a ladder compartment, a skip compartment, and a pipe line compartment.Arrangement of the shaft section is depicted in Figure 5.
According to the Chinese relevant specification for protecting the mine shaft, when the rollers sliding along the steel shaft guides are used, the clearance of sliding guide groove on each side should be kept less than 10-15 mm in order to maintain the regular operation of the hoister in the shaft [15].So the displacement of the main production shaft should not be larger than this threshold in order to guarantee the normal operation of the shaft.Currently, the hoister of the main production shaft runs well, indicating that the  clearance does not exceed the above-mentioned threshold.Most importantly, considering that the shaft will be used longer than 10 years, the 10 mm is selected as the upper limit for the shaft displacement, which is much less than the critical allowable displacement of the concrete liner.

In Situ Displacement Monitoring
In order to trace the response of the main production shaft at Xincheng Gold Mine during mining of the XI# ore body, it is essential to monitor the shaft displacement.In situ displacement monitoring can also provide data for validating the numerical simulation so as to improve the reliability of prediction.However, it is almost impossible to install monitoring equipment directly inside the shaft, because it may disturb the daily lifting operation.On the other hand, the shaft deformation is dependent on its adjacent rock stratum movement during underground excavation, and the method of monitoring the surrounding rock mass around the shaft is applied extensively to evaluate the displacement and stability of the shaft.

Installation of the Monitoring Equipment. The Sliding
Micrometer, a high-precision strain meter (±0.002 mm/m),  is applied to monitor the rock mass at −380 m level near the shaft.This monitoring system consists of the probe, cable, rod, readout unit, data-processing unit, and calibration device.It is noted that the probe uses the ball-and-cone positioning principle in the measuring marks of the measuring tube, and it needs regular calibration before and after each measurement.
4 Sliding Micrometer pipes, M1, M2, M3, and M4, were installed in the monitoring boreholes and then were grouted with cement mortar which should be cured for about 24 hours before measurement operate, in case the cement hardening may have detrimental influence on the monitoring displacement.On-site measurement diagram is shown in Figure 6.Over more than the past three years, it is found that the max cumulative displacements at monitoring boreholes M1, M3, and M4 are 3.5 mm, −1.4 mm, and −0.3 mm, respectively, while M2 is up to 28.9 mm, because the borehole M2 was located in an isolated pillar after one slice of ore was mined out, Figure 7.However, as this mined-out area was backfilled, the displacement monitored at this borehole becomes stable, which indicates that the noticeable displacement monitored at M2 results from the mining activities.However, this mining-induced displacement in surrounding rock mass has little influence on the stability of shaft because the stope is far enough away from the protected shaft.
Correspondingly, the monitoring displacements of M3 and M4 (−1.4 mm and −0.3 mm) reflect the shaft displacement more directly; it indicates that the shaft displacement   is much smaller than the threshold of 10 mm.Although a certain degree of fluctuation exists among the data during this monitoring period, according to the latest data, the max lateral displacement of shaft is 1.4 mm.

Sensitive Analysis of Factors Affecting the Shaft Displacement
The deformation mechanism of shaft can be very complicated; it is necessary to find the dominant factors affecting the shaft displacement based on a concept model.The influential factors include lithology of the rock mass, in situ stress state, mined-out area shape in terms of the ore body shape, location of the shaft, mining method and mining intensity, dynamic disturbances such as blasting vibration, special geological structure such as faults, groundwater, characteristics of the lining, service life, and maintenance.The concept of orthogonality is important in design of experiments because it says something about independence.Experimental analysis of an orthogonal design is usually straightforward because you can  estimate each main effect and interaction independently.In this section, orthogonal experimental design was used to find the dominant one that affects the shaft displacement remarkably.

Orthogonal Experimental Design.
The orthogonal numerical experiment has been done in virtue of the calculation of different scenarios with FLAC 3D software, as shown in Figure 9, and the numerical model is simplified to be a conceptual model in plain strain condition, whose size is 1000 m * 800 m * 20 m, and the 3rd dimension is just 20 m, in which direction the displacement is not concerned in this simplified model; in this respect, a 2D model may be also acceptable for this case.And the shaft is not considered in the grid of FLA3D in order to eliminate the error that may be induced because its diameter of 5 m is too small to be considered.In this respect, the displacement of the shaft is interpolated from the nodes nearby.The maximum and minimum principal stresses occur in the  and  directions while the gravitational stress is applied along the  direction, specifically,   = H,   = H ( = 1.0, 1.5, and 2.0, as shown in Table 4), and  ℎ = 0.5H.The stress boundary conditions are adopted for the two numerical models in this paper, and when the calculation reached an equilibrium status, the initial displacements are set to zero, and then, the new displacements are considered to be caused by mining activities, because this numerical model is only a simplified concept model, primarily aimed at finding the dominant factors affecting the shaft displacement in general situations and providing useful references to the researchers.Referring to the rock mechanics parameters of Xincheng Gold Mine, the orebody and surrounding rock mass are of the similar lithology in this model, as listed in Table 3.
Selecting the lateral pressure coefficient of the in situ stress state, distance between the shaft and mined-out area, location of the shaft (at the footwall or hanging wall of ore body), dip of ore body, and thickness of ore body as controllable factors in the orthogonal experiments and the levels of these factors are designed and listed in Table 4.

Orthogonal Experimental Results. The orthogonal table
L 18 (2 × 3 7 ) is applied and the result is shown in Table 5.Through the analysis of Fisher's  distribution, numerical results indicate that dip of the ore body, lateral pressure coefficient, and distance between the shaft and mined-out area are significant factors affecting the shaft displacement.In contrast, thickness of the ore body and the location of the shaft are indistinctive.On condition that the shaft is located at the hanging wall, the deformation region is more extensive than that at the footwall, as shown in Figure 10.
Based on this orthogonal experiment, the dip of the ore body, lateral pressure coefficient, and distance between the shaft and mined-out area are specially considered in the later numerical simulation.In particular, the mining activity is most closely relevant to the distance between the mining stope and the shaft.

Modelling Strategy.
A bigger numerical model of the Xincheng Gold Mine was established subsequently, whose size is 1000 m * 1200 m * 1000 m, Figure 11;  is the strike direction,  is the dip direction, and  is the gravity direction.The Hexi Fault is adjacent to the bottom of the shaft and overlaps with the XI-8 and XI-9 ore veins at −530 m level, so it is meshed intensively.As a matter of fact, the Hexi Fault is a crushed zone with thickness of dozens of meters, rather than a general fault with the main fracture surfaces.Thus, the "interface" element is not adopted in this model; instead, it is embodied through a series of weaker mechanical parameters; as shown in Table 6, the buried depth from 465.50 to 529.74  is the location of Hexi Fault.The same method is also applied to the other two faults.
The ideal elastic-plastic model with Mohr-Coulomb yield surface in the FLAC 3D software is adopted, that is, where  1 and  3 are the maximum and minimum principal stress, respectively,  is internal cohesive force, and  is angle of internal friction.When  > 0, failure occurs.

Mechanical Parameters of Rock Mass. The modified
Hoek-Brown criterion with the RocLab software [16] is applied to obtain the reasonable rock mass parameters, which is expressed as  where m  is a reduced value of the material constant m  as given by and s and  are constants for the rock mass given by is a factor which depends upon the degree of disturbance, it varies from 0 for undisturbed in situ rock masses to 1 for much disturbed rock masses,   is a material constant according to the lithology, and geological strength index (GSI) was introduced to quantify the rock mass and surface conditions of the discontinuities.In this study, rock masses are treated as undisturbed materials; namely,  = 0 is set for all of the rock strata around the shaft,   is the uniaxial compressive strength (UCS) of the intact rock, and   can be selected according to Table 2.The values of geological strength index GSI cannot be obtained directly; they must be obtained based on in situ investigation of the rock mass around the shaft.However, it is scarcely possible to scan the joints and fissures on the surface of the shaft.Therefore, the in situ survey work was conducted in the drifts near the shaft.
The rock joints and interfaces at the drifts surface of all of the strata are scanned by using the ShapeMetriX3D system (Austrian Startup Company, 2008), and then, a 3D digital image of the surface containing the joints information can be built, as shown in Figure 12.Based on this real 3D rock surface, some visible rock mass features, such as the orientations, traces, areas, lengths, and position of joints, are retrieved.In addition, according to the field investigation, the GSI values of the rock mass around the shaft are estimated based on the geological descriptions in [17].The mechanical parameters of all the rock strata are listed in Table 6.Due to the difference of strata, the premining stress distribution depends on both depth and rock properties [18], as shown in Figure 13.Obviously, by using the stress boundary conditions, when the calculation reaches to an equilibrium status, the strata with higher Young's modulus hold higher stress than the rest of the geological formations.

Simulated Mining Sequence.
Limited by the computing ability of computer, this big numerical model, Figure 11, is difficult to be established so meticulously.As a result, it aims to give advices on the mining sequence from a macroscopic view at present, and the simulative excavation scope of each step is large.
In the numerical simulation, scenario I is approaching to the shaft gradually, which is advantageous for stopping mining according to the shaft displacement, and it consists of 6 steps, Figure 14 (a).Firstly, the XI-1, XI-2, XI-3, and XI-4 veins at −380 m level are mined in 4 steps towards the main shaft until a distance of 100 m; and then, the XI-8 and XI-9 veins at −530 m level are extracted; at last, the rest of the orebodies within 100 m distance from the shaft, including XI-5, XI-6, and XI-7 ore veins and a small part of the XI-1 vein,   are extracted.The output ore at each step is listed in Table 7. Besides, backfilling is simulated to be completed before the next step of extraction, whose mechanical parameters are listed in Table 8.

Numerical Simulation Results
. When the ore is mined out completely, the contour of displacement field is illustrated in Figure 15.Rock mass tends to move towards the mined-out area, and the maximum displacement occurs in the hanging wall, reaching to 56.7 mm.Because the shaft locates at the footwall of the XI# ore body, the excavation leads to a limited influence on it.However, due to high horizontal in situ stress and the mining activities, the Hexi Fault distorts seriously; as a result, the displacement of the shaft bottom is up to 12 mm.
Figure 16 shows the displacement along the lateral (strike, dip) and vertical direction, respectively.It is shown that the lateral displacement accounts for a significant proportion in the total displacement of the shaft.Lateral distortion along the strike direction is no more than 3 mm even when the XI# ore body is excavated completely.In contrast, the lateral distortion along the dip direction is much more serious, which exceeds 10 mm at the shaft bottom.In the vertical Just as the definition implies, the higher the MPI is, the smaller the incremental shaft displacement occurs when the ore production remains the same, or more ores will be mined out if the incremental shaft displacement remains constant.Substantially, the MPI value is a comprehensive evaluation index considering both the economic benefit and the safety of shaft.Obviously, the overall MPI value of the XI# ore body is 1 (100% to 100%).If the MPI is larger than 1, excavating this part of ore body is relatively better than the others.Otherwise, the mining conditions may be worse, especially, if the MPI value is negative, indicating that this mining step is not reasonable; namely, spring-back shaft displacement occurs.
Table 7 shows information of the scenario I that the MPI values declines gradually with mining operation, indicating that the overall mining sequence, excavating approaching to the shaft, is reasonable.However, the negative MPI value which appeared at step 5 indicates that the excavating of XI-8 and XI-9 veins should be coordinated in order to avoid this kind of excavation.
The modified mining sequence scenario II is to extract the XI-8 and XI-9 veins (Step 5) together with step 3 and step 4, Figure 12(b), while other steps stay the same as the scenario I; thus, only 5 steps are left in the scenario II.As for the scenario II, the orebodies XI-1 and XI-2, XI-3, XI-4, XI-8, and XI-9 at two sides of the shaft are mined simultaneously, which may be better for balancing the possible tilt of shaft.The relevant displacement results are shown in Figures 17(b) and 18, and the maximum displacement (at the shaft bottom) does not exceed 10 mm if all the orebody is mined out for the scenario II.
For the scenario II, the deformation patterns along the shaft depth change a little but the maximum shaft displacement decreases evidently.The MPI values are listed in Table 9; negative values do not emerge; correspondingly, displacement of the shaft also increases steadily in a small scale, and there is no spring-back displacement of shaft.
By comparison, scenario II is better in controlling the shaft deformation than scenario I, and the maximum shaft displacement is 9.7 mm when the XI# ore body is excavated completely.Although the maximum displacement does not exceed the safety threshold 10 mm, the MPI value at step 5 is too low, which means that the ore production within the distance of 100 m from the shaft is small but the influence on the shaft displacement is significant; as a result, this part of ore is suggested to be left as the shaft pillar.Compared to the original shaft pillar delineated according to the rock movement angles, more than 1000 kt ore can be mined out.According to the numerical simulation, the max shaft displacement at −380 m level is 1.6 mm at this time, and the monitoring shaft lateral displacement is 1.4 mm.In this regard, we deem that the numerical simulation results are well in agreement with the monitoring data; of course, it needs to be validated further against the in situ monitoring for a long time.Even so, the main value of this crude numerical model mainly comes from the comparison with itself; by this way, a better mining sequence and a proper mining boundary are obtained.
area at −380 m level at −530 m level

Figure 1 :
Figure 1: The overall layout of Xincheng Gold Mine.

Figure 2 :
Figure 2: Location map of Xincheng Gold Mine in Shandong Province, China.

Figure 3 :
Figure 3: The XI# ore body and shaft pillar.

Figure 2 .
Figure 2. The construction of Xincheng Gold Mine started in 1975, and it was completed and put into production in 1980.The shaft at Xincheng Gold Mine is 499 m deep, and the ground surface elevation is +33 m.According to the latest exploration near the main shaft, a new ore body numbered XI# that consists of 7 ore veins at −380 m level and 2 ore veins at −530 m level are found, Figure1.Distribution of the XI# ore veins is scattered, the average strike is NE57 ∘ , and the dip of the veins is in the 25-30 ∘ range, which is parallel to the stratigraphic dip, and detailed information of the ore veins is listed in Table1.In the beginning, a shaft pillar (Figure3) has been delineated according to the respective rock movement angle of different strata, Table2.Under these circumstances, more than 1140 kt ore cannot be mined out, which accounts for 67.4% of the XI# ore body.The first mining areas are designed in the XI-1 ore vein at −380 m level and in the XI-8 and XI-9 ore veins at −530 m level, and 280 kt ores have been mined out up to March, 2017.In order to protect the shaft, the following stopping sequence needs to be optimized and the stopping boundary needs to be determined.The upward horizontal cut and backfill stopping method has been adopted for many years at Xincheng Gold Mine, and a wealth of experience has been accumulated.Compared with the open stopping method and stope caving method, the stopping-and-filling method is better in terms of maintaining

Figure 6 :
Figure 6: On-site measurement using the Sliding Micrometer system.

Figure 7 :
Figure 7: Arrangement of the 4 monitoring boreholes of Sliding Micrometer at −380 m level.
The 4 monitoring boreholes are horizontal and about 10-12 meters deep; boreholes M1 and M2 are near the XI# orebody and their borehole mouths are 110 m and 130 m far away from the shaft, respectively.Boreholes M3 and M4 are near the shaft and their bottoms are 7 m and 8 m from the shaft wall, Figure 7.

4. 2 .
Monitoring Results.The initial monitoring values were collected on October 26, 2013, and the deformation of rock mass was obtained by subtracting the initial data from the subsequent monitored values.Excavating of the XI# ore body at −380 m level started on January 23, 2014.As of the latest collecting, March 27, 2017, more than 280 kt ores were mined out.The cumulative displacements of the 4 monitoring boreholes are shown in Figure 8. Positive value means tensile strain and the negative one stands for compressive strain.

Figure 9 :
Figure 9: Numerical model of FLAC 3D of the concept model.

Figure 10 :
Figure 10: Shaft displacement for different depths, the legend "20 ∘ , 10 m, 50 m, 1.0" means the dip and thickness of orebody are 20 ∘ and 10 m, the distance between the ore and shaft is 50 m, and the lateral pressure coefficient is 1.0.

Figure 11 :
Figure 11: The overall numerical model of Xincheng Gold Mine.

Figure 12 :
Figure 12: Joint sets (indicated in red and green colors) retrieved with ShapeMetriX3D.

Figure 13 :
Figure 13: Contour of initial maximum principal stress.

Figure 16 :
Figure 16: Shaft displacement during mining with scenario I.

Table 1 :
Information of the XI# ore veins.

Table 2 :
Main mechanical parameters of the rock samples retrieved around shaft.

Table 3 :
Mechanical parameters of the numerical orthogonal experimental.

Table 4 :
Design of factors and levels for the orthogonal experimental design.

Table 5 :
Variance analysis of the orthogonal experiment.△" indicate that the mean sum of square of this factor is smaller than that of error term, being incorporated into the error term. "

Table 6 :
Mechanical parameters of rock mass calculated according to Hoek-Brown criterion.Intact rock core cannot be obtained in the stratum from 0 to 33.8 m beneath the surface. *

Table 7 :
Details of scenario I.

Table 8 :
Backfilling mechanical parameters in the numerical simulation.

Table 9 :
Details of scenario II.