^{1}

^{*}

^{1}

^{1}

^{1}

^{2}

^{1}

^{1}

^{1}

Regardless of beneficial associated with internal waste dump (IWD) method, practices of this method within boundaries of pit-slope have some serious problems on stability issues due to this area is zone of potential failure. This zone is known as dynamic reactive zone which is easy to deform by external force, and inherent dangers of failure posing a threat to slope. Therefore, it is paramount to study the induced shear stress behavior in this zone particularly when IWD method is adopted within this zone. In this paper, a numerical study for investigating IWD-induced shear stress behavior has been carried out using Finite Element Method (FEM) with Strength Reduction approach. Different scenarios as per pit-slope depths, IWD heights and buffer zone lengths have been accounted and simulated using PHASE 2 to understand changes in induced shear stress imposed on the pit-slope. It is found that shear stress imposed on pit slope seems change dramatically with increasing IWD height for case of buffer zone length is less than 100-m-long.

Internal waste dump (IWD) within boundaries of pit-slope is one of solutions to cut waste material hauling cost, particularly when the waste contains low grade valuable mineral which is planned to proceed further in future. However, practices of this method is challenging considering to area within pit-slope boundaries are a potential failure zone (Hustrulid et al., 2013), as illustrated in

deform by external force. Thus, when an IWD material is dumped in this zone, the slope surface will easy to distress, and shear stress imposed on the pit-slope surface will swiftly increase. If this shear stress is not controlled, tensile stress-induced cracks can open with a subsequent reduction in the cohesive and friction forces restraining pit-slope in place. This situation can reduce the effective normal force on potential failure planes. When the depth of pit-slope increases, shear stress on the pit-slope surface will extensive due to increase in load gravity, the extent of this distressed zone will more radical accordingly. This condition brings the shear stress imposed on pit-slope is larger than the strength of the pit-slope, and it can cause the slope failure (Cho & Song, 2014; Scott et al., 2007). Furthermore, it will get worse when poor blasting practice-induced tension cracks occur in this slope (Khandelwal & Singh, 2009; Persson et al., 1994; Calder, 1977).

Theoretically, increase in pit-slope depth will be followed by increasing shear stress due to increasing gravity load imposed on the slope. And the pit-slope will be more distressed when an IWD is dumped at top of the pit-slope due to increasing induced shear stress. This shear stress on the pit-slope will decrease when the IWD is moved away from boundaries of the pit-slope. This study investigates shear stress behavior on pit-slope when those three parameters, i.e. pit-slope depth, IWD height and buffer zone length, are changed. In order to investigate and explore shear stress behavior on pit-slope in regards with intercorrelation between pit-slope depth, IWD height and buffer zone length, a study of numerical analysis for IWD-induced shear stress imposed on pit-slope has been carried out using Finite Element Method (FEM) with Strength Reduction approach to understand changes in stress reduction factor (SRF) with variation pit-slope depth and IWD height under different buffer zone length.

In numerical simulation, stability of pit-slope is typically expressed with factor of safety i.e. ratio between resisting and driving forces of the rock mass (Wesseloo & Read, 2009). It can also be defined with a ratio between the tangent stress related to the shear strength of rock mass and shear stress induced weight and other factors influencing the rock mass behavior. The equation is generally expressed as follows:

Factorofsafety = τ r / τ s . (1)

where, τ_{s} is shear stress, while τ_{r} is maximum shear resistance estimated by Mohr-Coulomb Failure Criterion. This maximum shear resistance can be calculated by equation:

τ r = c + σ tan φ . (2)

where, c is cohesion, σ is normal stress, φ is friction angle.

Considering to a decision-making process of numerical method that leads to a realistic explanation particularly in case of slope design (Rocscience, 2005; Itasca, 2004; Rockfield, 2001), the continuum model of PHASE2 is adopted in this study. The PHASE2 is an extremely versa-tile two-dimensional elastoplastic finite element stress analysis program. One of the major features of this software is the using of shear strength reduction (SSR) method that can be set fully automated and can be used with either Mohr-Coulomb or Hoek-Brown strength parameters (Rocscience, 2005). In slope stability analysis, the SSR technique is used to determine a stress reduction factor (SRF) or factor of safety value that brings a slope to the verge of failure (Hammah et al., 2005). In this study, SRF 1.25 is considered as standard of safety (Wasseloo & Read, 2009).

Two key parameters that are considered in slope stability analysis are cohesion and angle of friction. A slope is considered unstable in the SSR technique when cohesion and friction angle in the shear strengths of slope successively diminishes until the finite element model does not converge to a solution (Hammah et al., 2005). These conditions are expressed by equations:

c m o b = c / FactorofSafety . (3)

φ m o b = φ / FactorofSafety . (4)

where, c_{mob} is mobilized cohesion and φ_{mob} is mobilized friction angle.

A case from a surface porphyry copper and molybdenum (Mo) deposit mine where IWD is taking place has been considered. This studied mine currently is at level of 240-m-deep. The studied mine has 3 IWD sites where located 50 to 600-m-long from boundaries of pit-slope. Current situation of IWD of the studied mine is illustrated in

designed well, the pit-slope will be largely distressed by IWD-induced shear stress. Furthermore, if this situation is not properly controlled, slope failure may occur as demonstrated in

In this study, the numerical model was built, and set-up based on actual condition of the studied mine such as: 1) Lithology consists of three main layers i.e. oxidation zone, transition zone and fresh rock; 2) Rock properties of main layers and IWD material; and 3) Current and long-term mine planning.

The numerical modeling of this study is performed by adopting continuum mechanics code PHASE2 which developed by Rocscience. The model geometry is created generally based on actual field geometry of the studied mine. The investigation area is 4-km-long and 2-km-wide (

In the numerical modeling, the selection of parameters is often considered from the most difficult element in the generation of a model due to the high uncertainty in the parameter data-base, in which the geotechnical parameters were conventionally derived from laboratory tests. The tests were conducted for three sample types i.e. oxidation rock sample, transition rock sample and fresh rock sample. The rock properties from the laboratory test are summarized in

According to the literatures (Li et al., 2011; Jimenez-Rodriguez & Sitar, 2007; Low, 2007; Jimenez-Rodriguez et al., 2006), cohesion (c) and friction angle (φ) is key parameters in reliability analysis of slope stability. Based on the properties of rock samples given in

In order to investigate the shear stress behavior on safe distance limit of buffer zone for a safe pit-slope and IWD stabilities, the models were simulated for different geometry of pit-slope depth from 240 to 440-m-deep, IWD height from 0 to 150-m-high, and buffer zone length from 0 to 200-m-long. The concept of the numerical simulation can be seen in

Material | Unit weight, MN/m^{3} | Friction angle, degree | Cohesion, MPa | Tensile strength, MPa | Young’s modulus, GPa | Poisson’s ratio |
---|---|---|---|---|---|---|

Oxidation material | 0.023 | 26.5 | 0.5 | 0.6 | 1.5 | 0.254 |

Transition material | 0.026 | 32.4 | 0.6 | 0.8 | 2 | 0.254 |

Fresh rock | 0.027 | 36.6 | 0.8 | 1 | 3 | 0.332 |

Fresh rock Waste material | 0.025 | 35 | 0.24 | 0.05 | 0.05 | 0.3 |

There were in total 48 simulations that have been run to investigate the shear stress behavior on pit-slope for differences scenarios of pit-slope depth, IWD height and buffer zone length. The simulation results are summarized in

Case 1: Effect of IWD on shear strain behavior on pit-slope.

In this case, the study was done based on simulation result for model of 240-m-deep pit-slope which illustrates actual pit-slope depth of the studied mine. The simulation was done for two conditions i.e. without and with IWD at top of the pit-slope. The simulation results can be seen in

Case 2: Effect of increasing pit-slope depth on shear stress behavior.

Pit-slope depth, d_{s} (m) | IWD height, h_{d} (m) | Buffer zone length, l_{b} (m) | ||||
---|---|---|---|---|---|---|

0 | 50 | 100 | 150 | 200 | ||

240 | 0 | 1.68 | - | - | - | - |

50 | 1.33 | 1.49 | 1.54 | 1.56 | 1.62 | |

100 | 1.14 | 1.31 | 1.41 | 1.48 | 1.54 | |

150 | 1.06 | 1.19 | 1.34 | 1.35 | 1.36 | |

340 | 0 | 1.46 | - | - | - | - |

50 | 1.22 | 1.26 | 1.33 | 1.36 | 1.42 | |

100 | 1.12 | 1.17 | 1.27 | 1.28 | 1.3 | |

150 | 0.94 | 1.07 | 1.19 | 1.15 | 1.2 | |

440 | 0 | 1.36 | - | - | - | - |

50 | 1.15 | 1.21 | 1.26 | 1.29 | 1.32 | |

100 | 1.07 | 1.13 | 1.16 | 1.19 | 1.22 | |

150 | 0.87 | 0.9 | 1.12 | 1.15 | 1.19 |

in gravity load degrades mobilized cohesion and mobilized friction angle step by step until the slope reaches to a quasi-static state. This situation can lead hazardous condition when an IWD is dumped just next to boundaries of the 440-m-deep pit-slope. This dumping work can trigger the pit-slope more distressed due to increase in IWD-induced shear stress. This cohesion and friction degradation process will drive the slope model to deform until the bond material break. The bond breaks due to shear failure will occur at the toe of the pit-slope. After the gravity load accumulates at toe of the pit-slope, tensile failure will occur at the crack tip. If this tensile failure propagates toward surface, it can bring the slope collapse (Nishimura et al., 2010; Kaneko et al., 1997).

Case 3: Effect of increasing IWD height on shear stress behavior on pit-slope and IWD.

Case 4: Effect of buffer zone on shear strain behavior on pit-slope with IWD.

The increase in IWD-induced shear stress on pit-slope can be restrained by introducing buffer zone between IWD and pit-slope.

In this section, the discussion is constructed by numerical evidence i.e. shear stress behavior for different geometry of pit-slope, IWD and buffer zone. The governing configuration is derived for three different pit-slope depths integrated with three IWD heights and five different buffer zone lengths.

three IWD heights (i.e. 50, 100, 150-m-high) for 0-m-long buffer zone. The figure shows that the SRF decreases with increasing IWD height. The SRF also decreases with increasing pit-slope depth. A similar characteristic with this correlation as shown in

S R F = k e − b h d . (5)

where, h_{d} is IWD height, and b and k are a constant. Considering to

The constants b and k that illustrates correlation between SRF and IWD height for any depths of pit-slope and length of buffer zone are given in

b = − 2 × 10 − 5 l ′ b + 0.003 . (6)

where l ′ b is buffer zone length in a range of 0 to 100-m-long.

Constant k which is physically representing SRF looks like controlled by pit-slope depth and buffer zone length. In order to confirm it, constant k is

Buffer zone length, l_{b} (m) | Slope height, d_{s} (m) | Constant | Corr. Coef. | |
---|---|---|---|---|

b | k | R | ||

0 | 240 | 0.002 | 1.4697 | 0.998 |

340 | 0.003 | 1.4108 | 0.998 | |

440 | 0.003 | 1.3522 | 0.996 | |

50 | 240 | 0.002 | 1.6582 | 1 |

340 | 0.002 | 1.3708 | 1 | |

440 | 0.003 | 1.4407 | 0.995 | |

100 | 240 | 0.001 | 1.6407 | 0.999 |

340 | 0.001 | 1.4105 | 1 | |

440 | 0.001 | 1.3259 | 0.997 | |

150 | 240 | 0.001 | 1.688 | 0.999 |

340 | 0.002 | 1.4905 | 0.999 | |

440 | 0.001 | 1.3557 | 0.997 | |

200 | 240 | 0.002 | 1.7899 | 0.997 |

340 | 0.002 | 1.5426 | 1 | |

440 | 0.001 | 1.3778 | 0.996 |

plotted as function of pit-slope depth for any length of buffer zone (i.e. 0, 50, 100, 150, 200-m-long), as shown in

more than 100-m-long, constant k decays in similar rate of 0.001. The correlation between k for buffer zone length less than 100-m-long can be expressed in equation as follows:

k = c e − y d s . (7)

where d_{s} is pit-slope depth, and y and c are a constant. Constants y and c are given in

y = 6 × 10 − 6 ( l ′ b ) + 0.0004 . (8)

Different with constants b and y which is less influenced by buffer zone when the length is more than 100-m-long, the constant c in Equation (7) is certainly influenced by buffer zone without any conditions.

c = 3.8 × 10 − 3 l b + 1.6595 . (9)

where l_{b} is buffer zone length from 0 to 200-m-long.

Buffer zone length, l (m) | Constant | Corr. Coef. | |
---|---|---|---|

y | c | R | |

0 | 0.0004 | 1.6247 | 0.999 |

50 | 0.0007 | 1.886 | 0.712 |

100 | 0.001 | 2.0873 | 0.972 |

150 | 0.001 | 2.1851 | 0.997 |

200 | 0.001 | 2.4357 | 0.997 |

By substituting Equations (8) and (9) to Equation (7) to get constant k and substituting the constant k and Equation (6) to Equation (5), the sensitivity chart of IWD-induced shear stress behavior on a pit-slope under stability level SRF 1.25 was constructed. The sensitivity chart is given in

the data which plotted in the graph is for buffer zone length 10 to 100-m-long considering to a significant change in SRF under this buffer zone range of length.

It is observed in

Moreover, there are two convergence characteristics in interrelationship between pit-slope depth and IWD height for increasing buffer zone length. A slight convergence is found for case of buffer zone length less than 50-m-long. While a considerable convergence is found for case of buffer zone length more than 50-m-long. A slight convergence indicates a considerable sensitivity of increasing IWD height on safe deep limit of pit-slope. In this case, buffer zone is considered not effective to control IWD-induced shear strain imposed on pit-slope. It may due to the area within 50-m-long from the boundaries of pit-slope is zone of potential failure which is sensitive responsive and responsiveness to external force. While, a considerable convergence indicates less shear stress behavior on pit-slope due to substantial effectivity of buffer zone to control IWD-induced shear strain.

Internal waste dump (IWD) within boundaries of pit-slope is one of solutions to cut the waste material hauling cost owing to location of IWD close to the pit. However, this method also can introduce a serious problem on stability of the pit-slope considering area within the boundaries of pit-slope is a zone that is easy to deform. Therefore, it is paramount to investigate shear stress behavior on pit-slope with IWD case to prevent the slope failures. In this study, IWD-induced shear strain that imposed on pit-slope is simulated by using Finite Element Method (FEM) with strength reduction approach.

Simulation result shows that shear stress on pit-slope changes with changing geometry of pit-slope depth, present of IWD on top of pit-slope as well as IWD height, and buffer zone length. The change occurs particularly at toe region of pit-slope. If this induced-shear stress is not controlled, tensile failure will occur at crack tip accordingly. If this tensile failure propagates toward surface, the slope is very possible to collapse. Thus, control shear stress imposed on pit-slope is key on implementation of IWD.

Simulation results show that the induced-shear stress increases with increasing pit-slope depth and increasing IWD depth but decreases with increasing buffer zone length. This study exhibits the shear strength dependency of shear stress by sensitivity analysis. The sensitivity analysis is done based on performance of geometrical combinations analysis between pit-slope depth, IWD height and buffer zone length by means of numerical simulation.

According to the sensitivity analysis, which is concluded from 48 simulations, safe deep limit of pit-slope is considerably influenced by IWD height for buffer zone length less than 100-m-long. The shear stress on this case of pit-slope seems changed dramatically with increasing IWD height. It may due to the area within less than 100-m-long buffer zone is potential failure zone which largely distressed. Considering to these analyses results, a geo-technical engineering can have a predictive capability for design of IWD, and that predictive capability can only be achieved. Moreover, considering this study was conducted at low rainfall infiltration condition, further analysis of saturated stability analyses should be conducted to predict shear stress behavior on more critical conditions.

The authors express their gratitude to Mining Research and Design Center, Mongolian University of Science and Technology for providing the field data. The paper is the effect of the scientific project between Mongolian University of Science and Technology and Kyushu University.

The authors declare no conflicts of interest regarding the publication of this paper.

Wahyudi, S., Shimada, H., Sasaoka, T., Hamanaka, A., Amarsaikhan, T., Mao, P., Dintwe, T. K. M., & Moses, D. (2020). Study of Internal Waste Dump-Induced Shear Stress Behavior on Pit-Slope. Journal of Geoscience and Environment Protection, 8, 71-86. https://doi.org/10.4236/gep.2020.84006