Support Design of Underground Cut and Fill Mine by using Hybrid Numerical Empirical Model

Rock mass classifications are traditionally used to group areas with similar geomechanicalcharacteristics.This classifications is expectedto provide guidelines for stability performance, and to select appropriate support. The rock mass classifications have been successfully applied to many tunnel construction designs. Two widely-used rock mass classifications in underground mines are Rock Mass Rating (RMR) and Norwegian Geotechnical Institute (NGI, Q-system). The main objective of this paper is to predict the required support systems for cut and fill underground mine which focus on decline main ways. In this study, support system requirements were calculated based on rock mass classification systems, and considered of induced stresses. This research is important to propose support systems for underground mine openings on different level activities.


Introduction
Obtaining a stablethe underground mining isone of the main concernsforanunderground mine engineer. A stable underground mining is determined by the suitability and economical support system with a convenient excavation method.
Rock mass classifications are traditionally used to group areas with similar geomechanicalcharacteristics.This classifications is expectedto provide guidelines for stability performance, and to select appropriate support. The rock mass classifications have been successfully applied to many tunnel construction designs. Two widely-used rock mass classifications in underground mines are Rock Mass Rating (RMR) and Norwegian Geotechnical Institute (NGI, Q-system).
The main objective of this paper is to predict the required support systems for cut and fill underground mine which focus on decline main ways. In this study, support system requirements were calculated based on rock mass classification systems, and considered of induced stresses. This research is important to propose support systems for underground mine openings on different level activities.
Kirsch, was the first publisher for the induced stresses distributions of a circular tunnel [1].Hoek and Brown, developedanadditional tool that can be predicting the induced stresses distributions [2]. The impact of induced stresses on stope stability was applied to Mathews empirical stability graph method [3,4]. Rock massclassification, Qsystem is proposed for predicting underground support system [5,6]. The effects of induced stresses are determined by comparing the maximum tangential stresses (σ θ ) to the intact uniaxial compressive strength (σ c ) in a similar manner to the Mathews' method. Wattimena, 2003 applied the Mathews'methodon block caving mines in Australia and proposed a Hybrid NumericalEmpirical Model (HNEM) to predict induced stresses and to calculate the Q system for estimating support requirement on undercut and production levels [7].

Geology
The field site for this study is Cibaliungunderground gold minewhichlocated in the western part of Java Island, Indonesia (Figure1). The resource is estimated to be approximately 1.5 million wmt gold with grade is estimatedto be 9.8 ppm.The type of ore is vein with low sulphidation.Geological mapping and geotechnical description were conducted in the field.

Abstract
A detailed study to propose new design support systems for underground mine was carried out in this paper. The depth was increased to develop the production in accordance with the vein.The cut and fill mining method was selected for this mine in line with the geometry of vein, topographic specifications, and annual production rates.Excavationsize, distance between stope and decline, depth, and the impact of fault governedthe induced stresses. Different depths in multiple excavations were analyzed in this study especially on decline area. The rock properties of investigated areawere defined from the field and laboratory test for each rock type. Support systems for the decline were designed based on the suggestions of the Rock Mass Rating (RMR) system and the Q system for rock mass classification. Induced stress on different levels were calculated to define Stress Reduction Factor (SRF) and applied to predict the support system by Q system. New method namely Hybrid Numerical Empirical Model was developed to support requirement analysis in cut and fill mining method. The support system requirement influenced by depth. The shotcrete will be thicker and rock bolt space more closely as increasing the depth. The host rock of Cibaliung deposit is in Honje Formation consists of andesite, basaltic andesite volcanic and volcanic breccia intercalated with tuffaceous sediment [8]. Andesite rock types, consists of andesite, andesite breccia, polymictic and monomictic breccia. Those host rocks are altered by chlorite-adularia and smectite-illite. The type of ore is a vein with low sulphidation epithermal deposit dominated by quartz vein. This vein follows the Cibaliung fault within NNW-striking/ENEdipping, with two shoot target Cikoneng and Cibitung ore shoot. The Cikoneng ore geometry is 250 m length, 2-10 m width and 200 m depth, whereas Cibitung shoot geometry is 150 m length, 2-15 m width and 300 m depth with dip of ore is 80 o . Some minor faults occur throughout the hanging wall excavation with strike NE-SW dip 75 o . Angel of this fault with axis of tunnel is 21 o .The geology of Cibaliung can be seen in Figure 2 [9].
Three main rock types, quartz vein as ore body, andesite-breccia altered by chlorite-adularia on footwall abundance, and andesitebreccia altered by smectite-illite dominant on hanging wall rock are observed in the study area.

Engineering Properties
The rock mass rating (RMR) system [10,11] and the Qsystem developed by the Norwegian Geotechnical Institute [5] were used to classify rock masses. The physical and mechanical properties of the rocks were investigated based on rock mass classifications from the field and laboratory studies on intact rock samples.
These tests coveran evaluation of uniaxial compressive strength (σ c ), tensile strength (σ t ), Young's modulus (E), Poisson's ratio (v), density (ρ), internal friction angle (ϕ), and cohesion (c). Uniaxial compressive strength, triaxial compressive strength, Brazilian and density tests were conducted in accordance with the ISRM standard. The values of minimum and maximum UCS varies, with the lower mean of UCS is in hanging wall rock. The results of rock properties tests are summarized in Table 1.
In this study, RMR and Q-system classifications were applied rock mass characterization. Some parameters were determined for RMR, includes uniaxial compressive strength (UCS), Rock Quality Designation (RQD), discontinuity spacing, discontinuity conditions, and groundwater conditions. The results of this study are shown in Table 2 to Table 4 for RMR.RQD, number of joint sets (J n ), joint surface roughness (J r ), joint weathering and alteration (J a ), joint water reduction factor (J w ), and stress reduction factor (SRF) are parameters for Q to quantify rock mass. RMR rating can be determined from the classification proposed by Bieniawski, [10]. Meanwhile, the Q value of a rock mass can be calculated by the equation given below, as determined in the classification proposed by Barton et al. [5].

Design of Support Systems
Design of support systems was carried out by considering the in situ stress, stress around the openings, and evaluation of failure conditions. In situ stress was calculated from eq. (2) for different levels within each geological formation. .
. . Hoek and Brown, proposed equations to calculate tangential stress around the roof and the wall of the openings on eq. (3) and eq. (4), respectively. The constants A and B are taken as 3.2 and 2.3 for a horseshoe-shaped gallery [2].
Since the influence of multiple excavations could not be predicted by the abovementionedequations,numerical modeling analysis is a proposed as a solution. ThereforePhase2® software is useful to predict induced stress on two dimensions.
The tangential stresses are determined near the opening on the roof and wall side. Figure 3 shows the initial condition for numerical modeling.The results of tangential stresses by numerical modeling are given in Table 5andFigure4 shows the graphic of induced stress on different levels.     Tangential stresses at the boundaries increase as increasing depth of excavations. In general, the tangentialstresses at the boundaries of both hanging wall and footwall decline excavations in the sidewall decreased with the distance from the decline to the stope. Significant tangential stresses decreases was observed when the distance changed from 5 m to 10 m. The tangential stresses are relatively constant when the distance between declines and stope up to 30 m on the hanging wall decline and 20 m on the footwall decline.
The tangential stresses in the roof increase with increasing of depth. The trend of tangential stress for both hanging wall and footwall declines in the roofwas found to be relatively constant as increasing distance between declines and stope. These results are the basis data to predict damage and calculate the support requirement. Tangential stress on the roof is higher than that on the wall and considered to support analysis.

Prediction of Damage
The relation between intact rock material strength and tangential stress can be used to predict rock burst in underground. A high value for tangential stresses is a good indicator for rock burst. Rock burst potentially damage classified by Grimstad and Barton [6]. The classifications are summarized inTable 6. Calculated values between intact rock and tangential stress to predict rock burst potential are given in Table 7. The intact uniaxial strength of hanging wall decline was 24 MPa and footwall decline was 51 MPa.    Based on the relation between intact rock and tangential stress, the damage potential in the roof is higher than that in the wall sides. Given the condition in the hanging wall,a rock damage will be more likely to occur in the hanging wall when comparedwith the footwall excavations. The results suggest that there will be minor rock spalling when the activities below 300 m depth in the wall sides on hanging wall. However, low to moderate rock spalling will take place in the roof when the mine activities occurs under 150-200 m depth. At the depth lower than 300 m, high rock spalling is predicted to occur and necessitating a heavy support. Although the condition on the footwall was observed to be relative stable, low rock spalling could occur at the depth below 300 m.

RMR Support Design
Bieniawski, proposed guidelines for excavation and support of 10 m span rock tunnels in accordance with the RMR system. The current support system in Cibaliung underground gold mine determined by using RMR. The support systems of individual openings were designed according to the surrounding rock mass characteristic. The equations proposed by Barton et al. [12] and Unal [13], which includes rock mass classifications, excavation width, span excavation, and material rock density was used to determine the splitset number and shotcrete thickness to support requirements [13].The results from eq. 5 to eq.8are showninTable 8 and Table 9.     Based on the equations above, the hanging wall with poor rock mass needs 16 rock bolts and only 12 rock bolts for fair rock mass. In additional, some rock bolt with grouted are also needed to support the systems. Eight grouted rock bolts and six grouted rock bolts are required for poor and fair rock masses, respectively. Shotcrete will be used as the support system. The capacity of shotcrete proposed around 25 MPa for all area withthe minimum thickness of shotcrete is 45 mm for hanging wall poor rock type and 35 mm for fair rock type.
However, the capacity of rock boltandshotcrete for hanging wall and footwall are identical, the numbers of rock bolt on footwall are less than the hanging wall due to the difference rock masses and rock properties. The numbers of rock bolt for footwall are 13 and 10 for fair and good rock mass, respectively. Six grouted rock bolts are needed for fair rock mass, and five grouted rock bolts for good rock mass. The thickness of shotcrete for fair rock mass in footwall is similar with shotcrete on hanging wall around 35 mm, and for good rock mass the thickness of shotcrete around 25 mm for the footwall.

Q System Support Design
In 2003, Wattimena, etalhasproposed a Hybrid Numerical-Empirical model to predict the support requirements in block caving mining method [7]. This model is the combination of numerical modeling and empirical analysis to design the support requirements. From numerical modeling, induced stress on complex design of mine method can be predicted, and the result is applied to the empirical Q system to design support system. To calculate the Q values as eq. (1), induced stress is important as parameter for SRF. To describe SRF, Barton et al. [5] and Grimstad and Barton [6] proposed the relation between tangential induced stress and intact rock uniaxial strength [5,6]. Table 10 shows the classification of SRF, and the ratios between tangential stress and intact rock uniaxial strengthare given in Table 11. By using Q system to design the underground support, the influence of induced stress on differentlevels can be predicted. Based on SRF values, we can obtain minimum and maximum range of SRF as input to calculate the Q values. Based on the values of ratio between tangential stress and intact rock strength, the research area could be classified into four classes of SRF. The classification of the SRF is given in Table 12. To determine Q values, other factorsas seen in Table 13 should be considered as well. The results of Q values are shown in Table  13.When the SRF classification of Barton et al. [5] and Grimstad and Barton [6] was applied, the range of Q values are 0.01 -2.04 and 0.02 -3.06 (extremely poor to poor). Meanwhile, the Q values on footwall are higher than that on the hanging wall. The Q values are between 2.20-4.40 and 3.58 -7.15 which mean poor to fair. Based on Table 13, it can be seen that the lowest Q values based on calculation using eq. (1) is located on the roof. Therefore the roof condition will be considered to analyze the support requirements.
The Qvalue is then used to predict the supportsystem requirement. Two additional parameters should be considered in using the Q system for predicting the support systemfor an underground excavation are Excavation Support Ratio (ESR) and equivalent dimension.      be determined by plotting the D e values against Q that was proposed by Barton et al. [5]. The sample of the support requirement estimationis given in Figure 5.
For the hanging wall decline when the activities up to 88 m, the Q values arebetween 2.04-3.06. As a conservative approach, Q values of 2.04 is used to support system analysis.The support requires includeshotcrete 4-5 cmthickness, and bolt with length 1.7 m and spacing 1.8 m. The support required will be increased with the increasing depth. This can be seen from the bolt spacing. For 150 m depth, the Q values are between 1.02-1.53. Within this range of Q value, the supportsystemrequiresshotcrete with 4-5 cmthickness, and bolt with length 1.7m and spacing 1.7 m. The support system needs    For the footwall decline, the Q values are higher than that on the hanging wall decline. Similar with the hanging wall decline, the highest tangential stresses occur in the roof.The stress was observed to increase with increasingdepth. For example, when the activities took place at depth≤200 m, the Q value is estimated to be in the range of 4.40-7.15 and classified under faircategory. The proposed support system requiresbolt length 1.7 m with space 1.6 m. and shotcreteis not necessary. When the activities increase up to 300 m, the support system includes shotcrete 4-5 cm, bolt length 1.7 m and bolt space 1.8 m. Detail predicted supports are given in Table 14.

Conclusion
Recommended support systems for different levels (88 m, 150 m, 200 m and 300 m depth) were designed using a new method; Hybrid Numerical Empirical Model.
The following results were obtained from damage potential analysis. Damage potential will increase with increasingthedepth. Low to moderate rock spalling will occur for hanging wall on150-200 m depth in the roof, and there will be high rock spalling when the mine activities below 300 m depth. In contrast, footwall will start encounter low rock spalling when the activities below 300 m.
It was recommended by the RMR system classification to use rock bolts and shotcrete. In hanging wall decline, minimum rock bolts for poor rock mass are 16 pieces and8 pieces of rock bolts grouted. The shotcrete is required for this condition with minimum thickness is recommended 45 mm and the capacity minimum 25 MPa.However, footwall classified as good rock mass, the support system is needed. For the good condition of rock mass 10 pieces of rock bolt and 5 pieces of rock bolt grouted are required with minimum 25 mm thickness of shotcrete.
Based on Qsystem, which is analysed by using Hybrid Numerical Empirical Model the requirement of rock bolts length aresame for all rock mass classifications that are 1.7 m. Recommendedsupport system by Q classification influenced by depth. The shotcrete will be thicker and rock bolt space more closely as increasing the depth. When the activities at 88 m depth, the support system requires shotcrete 4 -5 cm thickness and bolt spacing 1.8 m.Meanwhile, the shotcrete is 12 -15 cm thickness and bolt spacing is 1.0 m are required for depth up to 300 m. From this study the advantage to calculate support requirements using SRF is influence of induced on different levels could be predicted.