Rock Slope Stability Analysis by Using Integrated Approach

Slope stability assessment is an essential aspect of mining and civil engineering. In this study, Songwe open-pit mine in Malawi was investigated to establish possible pit slope instability. In performing the analysis, an integrated approach entailing rock mass characterisation, kinematic and numerical methods were applied. Based on rock mass classification system, Songwe Hill carbonatite rock mass is characterised as a good rock but still it possesses numerous random discontinuities that present a complex challenge in geo-technical engineering. Dip 6.0 software was used in carrying out kinematic analysis based on the attributes of discontinuities. The results show that there is a 16% likelihood of planar failure in the divided slope sections of the planned pit. Thus, slope angle optimisation to 41˚ has been proposed as a counter-measure to minimise the potential risk of planar failure. At the optimised angle, the risk of planar failure could be reduced by 44%. On the other hand, wedge failure was found to be improbable since no joint intersections were found in the critical zone of potential failure. For numerical analysis, finite element code was applied using FLAC 3D 5.0 application. The results demonstrate that overall slope angle of 41˚ would offer a favourable balance between safety and mining economics as mining operations progress to deeper horizons thereby avoiding a costly push back solution due to instability.


Introduction
Slope stability assessment is an essential aspect of mining and civil engineering.Depending on the analysis technique, the assessment can be aimed at: identifying endangered slope sections, investigating potential failure mechanisms, determining the slope sensitivity to different triggering mechanisms, designing of optimal slopes with regard to safety, reliability and economics and testing the efficacy of different support and stabilisation options [1].Despite advancements in the computer-based slope stability analysis, slope instability still poses a major challenge in large mining operations.Llano-Serna et al. [2] quantify that roughly, two open pit failures occur worldwide annually.It is articulated that large and deep instability typically involves rock and is strongly influenced by the orientation and spacing of discontinuities existent in the rock mass as compared to near surface slope instability which is particularly influenced by soil properties and soil-water interaction [3] [4] [5] [6] [7].Read & Stacey [8] also highlight that as surface mining operations are expanding and getting deeper, the influence of in situ stresses could be of reckon to slope instability.This compounding effect in slope stability has to be accorded attention as the demand for mineral resources increase to cater for the world booming population.
In order to make accurate mine design, a thorough site investigation that includes geological and discontinuity mapping should precede stability studies to supply necessary and accurate input data for the stability analysis.This is because occurrence of slope failures in unstable mine slopes may affect the mining operations due to fatality, machinery damage and render the recovery of ore uneconomical if the mine was being actively mined [7] [8].For instance, in 2012 a fatal slope failure left casualties at Jwaneng large diamond mine in Botswana.
The incident lead the company to resort to a push back and dip slope mining solution after recognising the structural complexity of the country rock [9] [10].
The failure at the mine was apparently triggered by a planar release along a structural bedding that dips 20˚ -30˚ northwest which is the only section of the mine with a layered (anisotropic) rock mass [11] [12].Thus, from the initiation of the mining projects, it is crucial to undertake a robust stability evaluation and design to better understand the rock mass, and methods of analysis that can model the rock mass and take into account variability in all of the geotechnical parameters [13].
In the case of the study area, the discovered rare earth elements (REE) deposit has a surface disseminated mineralisation style that will inevitably require surface mining method, specifically open-pit mining.The geological and structural mapping of the carbonatite complex at the Hill has shown that the rock mass possesses numerous random discontinuities which present a complex challenge in geotechnical engineering.These discontinuities may adversely affect the rock slope stability, potentially culminating in failures.In addressing such a geotechnical engineering problem, the study employs a combination of methods to evaluate stability conditions of the mine.Kinematic analysis, as a conventional method, graphically examine modes of slope failures that can possibly occur in a rock mass with respect to an existing or proposed rock slope [5]  The method analyses the displacement of geological materials along structural planes with no reference to the forces that induce the movement.Thus, kinematic analysis, which is capable of handling large discontinuity data sets, is undertaken first to assess possible modes of failure.On the other hand, numerical modelling complements in understanding the stability conditions by simulating the effect of stress distribution on the slope.
A good number of cases have been investigated applying kinematic analysis [3] [17] [18] [19] [20].For instance, in a study by [18] found planar sliding failures were predominant in some of benches with high dip angle (45˚ ≤ α ≤ 55˚) which he described as not ideal for mine stability.Furthermore, the rock mass was prone to sliding along the foliations particularly at low angles (31˚ -33˚).
The condition was argued to be due to greater size of the component force in the direction of sliding of the foliation planes with an angle of 31˚ than the ones with higher dip angles as well as their surfaces with nearly planar-slightly rough surfaces and mica flakes positioned in such a manner that is parallel to the foliation surfaces.However, the common trend in studies that apply kinematic analysis, use the method independently without evaluating the stress conditions on the slopes.Therefore, this study takes an integrated approach in the assessment of pit slope instability on a planned Songwe open-pit mine for optimal design.

Location and Geology of Study Area
Songwe hill is situated in Phalombe district southern Malawi (Figure 1).The Hill has a North-South diameter of 800 m and measures 450 m East to West.
The planned open-pit mine at the site is anticipated to cover almost the whole stretch of the hill.The overall dimensions of the pit will be approximately 650 m North to South and 400 m East to West.Regarding pit height, it is expected that the deepest level will be 300 m when measured from the highest reduced level (RL) of the pit on the southern side of the pit (Figure 1(c)).The current, prefeasibility studies show that the mine will produce 8.536 million tonnes of ore and 26.930 million tonnes of waste at a 3.15 strip ratio [21].The mine life is projected to take 18 years of mining activity.
In terms of regional geology, the area lies within the Chilwa Alkaline Province (CAP) (Figure 1).Fenites form an aureole around the carbonatite intrusion.Large blocks of fenite show evidence of being in situ.They are interpreted as fractured blocks from the margins, or the roof of the carbonatite.It is postulated that the carbonatite intrusion never reached the surface since the fenite is continuous with only rare carbonatite veinlets [24].In terms of the texture, the fenites display a coarsegrained equigranular igneous texture, strongly suggesting an igneous protolith.
The fenitisation at Songwe is predominantly potassic with a composition of orthoclase and minor aegirine.
Secondary geological units at the site consist of veins of manganese and iron (Mn-Fe) and phonolite dykes.Broom-Fendley, et al. [24] presented new petrographic and geochemical observations following thin sections where he uncovered veins of flourite, apatite, calcite/ferroan-calcite, quartz and barite with accessories that include xenotime, zircon, rutile, hematite, synchesite and parasite.Regarding structural geology, the site is within the active tectonic environment of the Malawi Rift system (MRS) which is part of the main East African Rift System (EARS).Thus, faulting and development of joints may not be an uncommon phenomenon.However, the structural disruption at the carbonatite complex is not vividly reflected on a large scale except for a major fault on the eastern side of the hill that runs NW to SE (Figure 1(b)).The subtle evidence of structural deformation is revealed in sharp lithological breaks across the hill.
One of the structural patterns in the area is an apparent offset of fenite breccia body.The fenite-brecciated section is characterised by a depression between the apparent ridges on the north facing section of the hill as it can be observed in Figure 1(a).
Witley et al. [25] demonstrated that the lithological breaks correspond well with the faults interpreted from the ground magnetic survey although it is admitted that the fault traces were considered as an approximation (Figure 2).The low confidence in the magnetic data is premised on the fact that the resolution of the magnetic image was low at the scale of the geological map.
The surficial predominating structural features at the study site are joints.Surface geological mapping conducted revealed that the joints have different orientations.It was also discovered that most of the joints are healed by rough cementing calcite infill (Figure 3).This type of infill in the joints could generate enough impact to minimise chances of sliding compared to soft clay infill.This is because the binding effect of calcite could enhance cohesion between rock blocks.
Figure 2. Geological map with faults (long dashed lines) interpreted from a ground magnetic survey [25].World Journal of Engineering and Technology

Methods and Materials
The approach to the study is summarised in Figure 4.The input data for the study was obtained from field investigation and laboratory experiments.Field investigations involved geological and structural mapping and geotechnical logging.Geotechnical log sheets provided discontinuity information; rock quality designation (RQD), joint spacing and joint conditions carefully recorded from diamond drill core.The data was logged at a regular interval of a 3 meter run, as such an average value was computed to characterise the rock on each drill hole.
Persistence of joints was determined through scan line mapping of rock cuts created when making drill pads.The length of majority of the joints extended up to 3 meters.
To determine rock material properties, laboratory tests were conducted on representative samples collected from drill holes.The samples were prepared according to International Society of Rock Mechanics (ISRM) standards [26].
In order to comprehend the characteristic behaviour and rock conditions of the carbonatite complex, rock mass classification was performed based on a quantitative Rock Mass Rating (RMR) system and a correlation to qualitative system was done.Among the available classification systems, RMR system was selected based on its easiness and versatility in engineering practices, involving tunnels, mines, slopes, and foundations.In addition, the method has stood the test of time with a high global acceptance owing to its varied applications amounting to over 351 case histories [27].As applied in slope analysis the method makes use of the five basic parameters to characterise a rock mass.The parameters include: UCS of rock material, RQD, discontinuity spacing ( )   = ∑ [27]   (1 where p R is the rating of each of the five parameters. To perform the kinematic analysis, Dip 6.0 software [14], was used in plotting the discontinuities based on the dip and dip direction of the joints.Dip 6.0 is an application designed to analyse features related to engineering analysis of rock structures.The program permits the analysis and visualisation of structural data employing similar techniques applied in manual stereographic projections.The parametric setting of the analysis entailed friction angle set at 35˚ based on shear strength test results.This value is represented by a plane circle termed friction cone in the stereonet determined by counting the great circles of the stereonet from outside inwards.The dip of the overall pit slope was set at the maximum possible angle of 45˚and the direction was varied with respect to the slope face as 360˚ or 0˚, 90˚, 180˚, and 270˚.The slope angle (α) of the slope face is represented by curved plane in the stereonet.Lastly, dip of the discontinuities (β) and dip directions are represented by vector points since the number of joints is high.The analysis was carried out in a pole vector mode and discontinuities from drill hole logging were assumed continuous and through-going since there is no excavated surface to ascertain continuity.
For numerical analysis, finite difference method (FDM) using FLAC 3D 5.0 software was applied in simulating the potential slope instability of the pit slope.
The FDM approach was preferred because the constitutive models are handled with no adjustment to the solution algorithm as such the modelling is justifiably more accurate for plastic collapse loads and plastic flow [28].
where: τ is shear stress, c is cohesion, φ is internal friction angle, σ is normal stress, and F is slope safety factor.
The model was governed by the following boundary conditions; roller boundaries were assigned in the xand y-direction thereby fixing the boundary planes in the xand y-direction respectively, pinned boundary (i.e., constrained in the x-, yand z-directions) was applied at the bottom of the model and an uncon- ( ) ( ) ( ) ( ) Two scenarios were considered in the modelling.In the first scenario, shear strain behaviour on the pit-slope was investigated at the current planned depth of 250 m.The planned depth is within the bounds of proven ore reserve hence the geological confidence is high.The second scenario is for the ultimate depth of 300 m.At this depth, the geological confidence is relatively low since less than 10% of the drilled holes reached 300 m.In both cases, the analysis was conducted at different overall slope angle (OSA).The stability of the pit walls was also assessed by varying the rock mass parameters based on the GSI values that fall in the rock mass class.The adjusted parameters (UCS, cohesion, friction angle and tensile strength) were derived by using Rocklab 1.0 software [29].The application converts the Hoek-Brown classifying parameters into a Mohr-Coulomb fit applicable to FLAC 3D .

Rock Mass Classification
The major aim of rock mass classification is to provide a common basis for effective communication concerning geomechanical problems that could be associated with the rock mass condition.As such, the representative RMR value established in the study and its interpretation adhere to a set of the system's guidelines.The summary of the results is presented in Table 2. From the results, the rock mass at the study site has an RMR value of 62 which can be described as a good rock with a rock class of "II" [30].
The evaluated quantitative RMR value was correlated to qualitative rock mass World Journal of Engineering and Technology characterisation for validation.Geological Strength Index (GSI) system, introduced by Hoek in 1994, is a widely accepted qualitative scheme for describing the rock mass reliably fast and simple [31].The method is a result of a combination of observations of the rock mass conditions (Terzaghi's descriptions) with the relationships developed from the extensive experience gained using the RMR-system.Thus, the relationship between rock mass structure and rock discontinuity conditions is used to estimate an average GSI value represented in the form of diagonal contours.The GSI value was empirically correlated to RMR based on the relationship formula by Singh & Goel [30] as: The GSI value of the Songwe carbonatite is computed to be 57.Thus, the qualitative interpretation of the rock mass based on computed GSI value implies that the rock mass is a good rock with a structure of multi-faceted angular blocks generated by three or more joint sets with a surface quality of being rough, slightly weathered and iron stained surface.The description is quite in agreement with the quantitative RMR interpretation of the rock mass characterisation.The GSI value is significant in this study for assessing rock strength condition at the modelling stage.

Planar Failure
For planar failure to occur, the joint plane has to dip downward in almost the same direction as the slope face at a less steep angle.Where the joint plane is exposed in the slope face, the plane is regarded to "daylight", a condition conducive for rock mass sliding.Hoek E. & Brady [5], Kliche [32] and Goodman, 1989 [4] proposed that the following conditions must be fulfilled for planar sliding to occur on a discrete plane.
1) The plane on which sliding occurs must strike parallel or nearly parallel (within approximately +/− 20 degrees) to the slope face.
2) The failure plane must "daylight" in the slope face.This implies that the dip of the slope face must be steeper than the apparent dip of the failure plane.
3) The dip of the failure plane must be greater than the angle of friction of this plane.
The performed analysis shows that the risk of potential planar failure is relatively high on the north facing section of the Songwe mine.The kinematic analysis indicates that 83 joints in the north facing section would be critical in influencing planar slope failure of the planned pit (Figure 5).This represents a 17% probability of failure for a sectional north dipping joints and a crude 6% probability with respect to all the plotted joints.In the east and west facing section of the planned pit, the risk of failure is relatively low with 55 and 51 joint planes plotting in the critical zone of potential failure respectively (Figure 6 & Figure 7).Statistically, the probability of failure for the slopes is at 17% in both cases when considering the discontinuities in the respective sections.However, the crude probability of failure is at 4% in the eastern section and 3.7% in the western section.
The section with the lowest risk of failure is the south facing section, which has 36 joint planes plotting in the critical zone of potential failure (Figure 8).
This represents a 13% sectional probability of failure and a crude probability of 2.6%.Thus, the sectional overall probability of planar failure of the entire mine is 16% which is high against the acceptable stability criteria empirically determined by [8] [17] [33].
To minimise the potential risk of planar failure identified in the analysis, slope angle optimisation would be crucial.In mining, exploitation of ore in open pit mine involves relating slope angle design to a pair of contradiction between resource recovery and safety.The heightening and steepening of excavation bench slopes may affect the security of the overall slope thereby inducing instability.Thus, Slope angle optimisation helps in striking a balance between safety and mining economic profits.Based on the current analysis, an optimised slope angle could be adopted within an acceptable angle between 40˚ and 45˚.

Wedge Failure
Wedge failure can only occur on two intersecting planes of the joints.Generally, the trend of the intersection line that satisfies wedge failure must be orientated within 90˚ of the dip and the plunge of the intersection line must "daylight" in the slope face implying that the dip of the slope face has to be greater than the plunge of the intersection line [4] [5] [6].In this kinematic analysis, all sections of the planned open pit registered no critical intersections in the potential failure zones both in the primary and secondary zone (Figures 14-17).Within the primary critical zone, any intersection plotted could directly satisfy the frictional and kinematic conditions for wedge failure.On the contrary, the intersection in the secondary critical zone would be such that the sliding would occur on a joint plane with a dip greater than the friction angle and the other joint acting as a release plane for wedge failure.Since no intersections plotted in the critical zones of potential failure, wedge failure could be an unlikely phenomenon unless artificial propagation of joint planes under active mining operations.

Pit Height and Slope Stability
In slope stability, FOS is used to determine the stability conditions.FOS is a ratio between supply (shear strength) and demand (shear stress).The minimum requirement for stability is that FOS should be at unity.However, in mines, the standard used criterion for safety assurance is for FOS to be ≥1.2 [8] [33].Figure 18 presents the stability conditions evaluated in case 1 where the height is at 250 m.As aforementioned, the OSA was varied from 45˚ to 40˚.At 45˚ the FOS is     wall stability at 45˚ and 43˚ (Figure 19).On the other hand, it can be noted that the FOS values are still reliable at OSA of 41˚ and 40˚.The FOS at 40˚ which is 1.26 is more reliable than the FOS at 41˚ (1.22).Nevertheless, in mining, efforts are made to balance between safety and mining economic profits.That entails making slope angle steeper and high to avoid excavating much waste without effecting instability.In the event that the mine extends to 300 m, then it would be recommendable to optimise the OSA to 41˚.The proposed optimisation option at the initiation of the mining operation could assist in evading costs due to push backs or onset of instabilities.For instance, a study by Bye and Bell (2001) [34], demonstrated that the optimization of slope design at Sandsloot Platinum mine extended the life of mine allowing for more two benches after the 300m planned [36] where decreasing displacement with increasing depth of excavation was noticed and supported by physical models.
This phenomenon can be attributed to stress distribution on the slope face.
The high maximum principal stress concentration observed at the toe region due to gravity loading and angle of pit slope is reduced towards the crest as a result World Journal of Engineering and Technology of thrusting from the bottom where the minimum principal stress is at low concentration.Similarly, maximum shear stress concentration reduces towards the crest of the slope.Tensile stress, which is not significant in as far as failure is concerned, is developed on the upper portion of the slope face.These stress interactions transform the confined state of stress into unconfined state due to the excavation of the pit.Furthermore, the length of the slope face provides for stress release in the redistribution of stresses.This could account for the small displacement values in the pit wall at increased excavation depth.

Significance of GSI in Slope Stability Analysis
Since rock condition has a direct impact on rock behavior, the values of GSI describing the rock mass state were applied in numerical modelling to establish the strength conditions of the rock slope.The results of the analysis are graphed in Figure 21.The GSI values appear to have a direct proportionality to material stability on the pit slope.At a lower GSI value of 55, the safety factor is 3.12, 3.08, 3.0 and 2.9 at 40˚, 41˚, 43˚ and 45˚ respectively.When the GSI value is increased to the upper limit of the rock class, the safety factor increases to 5.75, 5.67, 5.56 and 5.47 at 40˚, 41˚, 43˚ and 45˚ in the order of OSA.The values of rock slope stability for the currently established GSI value are shown in the red outline (Figure 21).The trend corroborates a non-failure condition with all values above unity as observed in the FOS using unadjusted parameters (see Figure 18 and Figure 19).However, the stability condition wane with increase in OSA.Sonmez and Ulusay (1999) [37] in five well-documented slope instability examples that they studied, they included the GSI values in the back analysis of the failed slopes.The outcome indicated that the use of GSI value confirmed the limit equilibrium analysis outcome for the failed slopes in all the cases.Thus, it can be concluded that the use of GSI in stability analyses gives an insight of the slope strength conditions of rock mass to failure.

Conclusions
In this paper, potential slope failure modes and stability conditions of the pit slope were investigated using kinematic and numerical approaches.
The kinematic analysis of the planned mine at Songwe carbonatite complex revealed a potential planar slope instability at a high possible OSA of 45˚.The

Figure 3 .
Figure 3. Predominating calcite infill in joints (Pictures taken by first author).

sJ
, condition of discontinuity surface (joint length ( ) l J , joint aperture ( ) a J , joint roughness ( ) r J , joint infill ( ) i J and joint weathering ( ) w J ) and groundwater conditions (GW).RQD is derived as a percentage of the sum of intact core pieces longer than 10cm divided by the total length of core run.The RMR method can be presented statistically as: D. Moses et al.DOI: 10.4236/wjet.2020.83031411 World Journal of Engineering and Technology
The analysis was performed in elasto-plastic state in finding the solution with Mohr-Coulomb constitutive model and failure criterion.To determine the safety factor of the stability, Shear Strength Reduction (SSR) technique is used in FDM code.The SSR involves a systematic iterative search of the factor of safety (FOS) value that stretches the slope to limits of failure.The Mohr-Coulomb materials were factored or reduced to determine the safety factor as presented in Equation (1): D. Moses et al.DOI: 10.4236/wjet.2020.83031412 World Journal of Engineering and Technology

Figure 5 .
Figure 5. Planar slope stability condition in north facing pit slope at an angle of 45˚.

Figure 6 .
Figure 6.Planar slope stability condition in east facing pit slope at an angle of 45˚.

Figure 7 .
Figure 7. Planar slope stability condition in west facing pit slope at an angle of 45˚.

Figure 8 .
Figure 8. Planar slope stability condition in south facing pit slope at an angle of 45˚.

Figure 9 .
Figure 9. Planar slope stability condition in north facing pit slope at an angle of 41˚.

Figure 10 .
Figure 10.Planar slope stability condition in east facing pit slope at an angle of 41˚.

Figure 11 .
Figure 11.Planar slope stability condition in south facing pit slope at an angle of 41˚.

Figure 12 .Figure 13 .
Figure 12.Planar slope stability condition in west facing pit slope at an angle of 41˚.

1 .
22 which is reasonably satisfactory for mine stability.When the OSA is successively reduced to 43˚, 41˚ and 40˚, in turn, FOS increases from 1.29, 1.32 and 1.34 respectively.Thus, OSA in inversely proportional to FOS.

Figure 14 .
Figure 14.Wedge slope stability condition in north facing pit slope at an angle of 45˚.

Figure 15 .
Figure 15.Wedge slope stability condition in east facing pit slope at an angle of 45˚.

Figure 16 .
Figure 16.Wedge slope stability condition in south facing pit slope at an angle of 45˚.

Figure 17 .Figure 18 .
Figure 17.Wedge slope stability condition in west facing pit slope at an angle of 45˚.

Figure 19 .
Figure 19.Pit wall stability condition at 300 m height.

Figure 20 .
Figure 20.Horizontal displacement of rock mass on the pit wall.
boundary at the top of the model.The material properties used in the study are presented in Table1.Rock density (γ), UCS, Young's modulus (E),

Table 1 .
Mechanical properties of the rock units.
*Computed as the sum of the rock mass characterisation parameters.World Journal of Engineering and Technology