Search Documents

By D. F. Malan, J. A. L. Napier

A number of shallow coal or hard rock mines employ pillar mining systems as a strategy for roof failure control. In certain platinum mine layouts, pillars are designed to ‘crush’ in a stable manner as they become loaded in the panel back area. The correct sizing of pillars demands some knowledge of the pillar strength and the overall layout stress distribution. It is particularly important to understand the impact of the layout geometry on the effective regional ‘stiffness’ of the rock mass around each pillar. An important design strategy is to model relatively detailed layout configurations which include a precise representation of the local pillar layout geometry and to analyse multiple mining scenarios and extraction sequences to select optimal pillar sizes and barrier pillar spacing. Although computational solution techniques are now impressive in terms of run time efficiency, a major difficulty is often encountered in assigning suitable material properties to the pillars and in devising an effective material description of the layered rock strata overlying the mine excavations. This paper outlines an efficient numerical strategy that can be used to assess large-scale pillar layout performance while retaining the ability to modify individual pillar constitutive behaviour. The proposed method is applied to selected layouts to compare estimated average pillar stress values against values determined by detailed modelling and against observed behaviour.

May 3, 2023

By D. F. Malan, D. Ile

This study explores the use of backfill in hard rock bord-and-pillar mines to increase the pillar strength and extraction ratio at depth. The use of backfill will also minimize the requirement for tailings storage on surface and the risk of environmental damage. A literature survey indicated that backfill is extensively used in coal mines, but rarely in hard rock bord-and-pillar mines. To simulate the effect of backfill confinement on pillar strength, an extension of the limit equilibrium model is proposed. Numerical modelling of an actual platinum mine layout is used to illustrate the beneficial effect of backfill on pillar stability at greater depths. The magnitude of confinement exerted by the backfill on the pillar sidewalls is unknown, however, and this needs to be quantified using experimental backfill mining sections equipped with suitable instrumentation.

May 3, 2023

By D. F. Malan, P. M. Couto

We propose a layout design for the UG2 Reef where weak geological alteration layers are present. The collapse of the Everest platinum mine in South Africa indicated that these layers substantially weaken the pillars. The popular Hedley and Grant pillar strength formula cannot be used where these alteration layers are present. Underground investigations at Everest mine and numerical modelling of the layout were conducted using the TEXAN code and a limit equilibrium model. Simulations of a collapsed area and an intact area allowed for a preliminary calibration of the model. This was subsequently used to explore modified layouts for these ground conditions. An alternative is to compartmentalize the blocks of ore using barrier pillars. The numerical modelling predicted that the barrier pillars appear to remain stable even in the case of large-scale collapses, provided their width exceeds 25 m. Main access routes into the mine can be protected by a double row of pillars at least 15 m wide to provide a safe travelling way. As a cautionary note, these conclusions are based on the model calibration and this needs to be refined in future. Calibration of the limit equilibrium model remains a challenge owing to the large number of parameters involved.

May 3, 2023

By J. A. Maritz, D. F. Malan

Pillar strength is affected by pillar shape, but this has largely been ignored in past research studies. Bord-and-pillar layouts are typically designed using empirical strength equations developed for square pillars. Owing to the poor quality of pillar cutting, many hard-rock pillars have an irregular shape and it is not clear how this affects pillar strength. Furthermore, the strength of rectangular pillars in comparison with square pillars is also difficult to quantify. The ‘perimeter rule’ is widely adopted for rectangular pillars, but its applicability for pillars with irregular shapes has never been tested. We used numerical modelling in this study to investigate the effect of pillar shape on strength. An analytical limit equilibrium model of a square and a strip pillar also provided useful insights. For slender pillars, the strength of a long rib pillar is essentially similar to that of a square pillar. In contrast, for rib pillars with a large width to height ratio, there is a substantial increase in strength. The study found that the perimeter rule should not be used for irregularly shaped pillars. Displacement discontinuity modelling, using a limit equilibrium approach, is proposed as an alternative to determine the strength of these pillars.

May 3, 2023

By F. J. N. Bruwer, T. R. Stacey

The research described considers whether the variability in coal material strength, as derived through a series of uniaxial compressive strength (UCS) tests, could be used to indicate the variability in coal pillar strength. The aim is to be able to use a distribution of UCS tests as input into the coal pillar strength calculation. This will allow the pillar design to be expressed in terms of a probability of failure rather than as the commonly used safety factor. To achieve this, the bulk strength factor associated with commonly used pillar strength formulae was replaced with a distribution of UCS results divided by an adjustment factor. The factor was determined so as to ensure that the resulting bulk strength does not deviate from the statistically determined bulk strength published in the original formulae. This approach enabled pillar strength distributions to be obtained using industry-accepted strength formulae, subsequently allowing for a probability of failure to be calculated for a specific pillar design. Using a regional coal material strength curve as a baseline, coalfields in which the coal is stronger than the regional mean can be identified and the pillar designs optimized. This is based on the stronger coals achieving lower probabilities of failure at similar safety factors. The research has considered actual UCS data from multiple mines in the Mpumalanga coalfields of South Africa, and has proved that the variability in material strength between coalfields could allow for some optimization using the proposed approach. Based on the data used in the study, a 2.78% increase in extraction could be achieved. However, further research will be required to validate the results of the study in an underground environment.

May 3, 2023

By D. F. Malan, R. P. Els

The limit equilibrium model, used in displacement discontinuity codes, is a popular method to simulate pillar failure. This paper investigates the use of physical modelling to calibrate this model. For the experiments, an artificial pillar material was prepared and cubes were poured using the standard 100 mm × 100 mm civil engineering concrete moulds. The friction angle between the cubes and the platens of the testing machine was varied by using soap and sandpaper. Different modes of failure were observed depending on the friction angle. Of interest is that significant loadshedding was recorded for some specimens which visually remained mostly intact. This highlights the difficulty of classifying pillars as failed or intact in underground stopes where spalling is observed. The laboratory models enabled a more precise calibration of the limit equilibrium model compared to previous attempts. Guidelines to assist with calibration of the model are given in the paper. The limit equilibrium model appears to be a useful approximation of the pillar failure as it could simulate the stress-strain behaviour of the laboratory models.

May 3, 2023

By D. F. Malan, T. E. Oates

A recent experimental pillar extraction project at a UG2 bord-and-pillar mine presented a unique opportunity to compile a new pillar database. Currently, the South African hard rock bordand- pillar mines are designed using the Hedley and Grant formula with a modified K-value. This empirically derived formula was developed for uranium mines in the Elliot Lake district of Canada. The use of this formula for the design of pillars in South Africa is questionable. Very few pillar failures have nevertheless been observed and its current calibrations for the various reef types are possibly too conservative. A new UG2 pillar database of 66 pillars, of which seven are classified as failed, was compiled by the authors. This enabled a revised ‘first-order’ calibration of the K-value for the Hedley and Grant formula. The new estimated value for the UG2 is K = 75 MPa. This gives a pillar strength that is more conservative than the PlatMine formula. This work should nevertheless be considered as only a preliminary calibration as the database was small. Further work is also required to determine whether the exponents in the formula for the width and height parameters are appropriate for UG2 pillars.

May 3, 2023

By D. F. Malan, J. C. Esterhuyse

In this study we explore the use of a novel numerical modelling approach to study the effect of pillar reinforcement on pillar stability. Case studies in the literature indicate that tendons, strapping of the pillars, and shotcrete or thin spray-on liners are commonly used to reinforce pillars. No clear methodology exists to select the type of support or to design the capacity of the support required, however. This has led to ongoing collapses in some mines in spite of heavy support being used to reinforce unstable pillars. A limit equilibrium model with confinement on the edge of the pillar was used to simulate the interaction of the support with the failing pillar. The model correctly predicts that an increase in confinement will lead to a decrease in the extent of pillar failure. As the displacement discontinuity boundary element method allows for the efficient solution of largescale bord-and-pillar layouts, the effect of pillar confinement can now be studied on a mine-wide scale. Accurate calibration of the limit equilibrium model is, however, required before this method can be used for the design of effective pillar support.

May 3, 2023

By J. van der Walt, W. Spiteri

Flyrock remains a significant threat to the health and safety of mine employees and integrity of infrastructure, as well as to the safety of the neighbouring communities and their property. This investigation was motivated by the general lack of fundamental research and mathematically quantifiable data in the literature regarding the relationship between blast design parameters and their impact on flyrock. The focus was to develop a concept that can be used to mathematically quantify the trajectory of flyrock resulting from a blast, which can be used for future research. The ultimate goal for this technique, once it has been fully developed, is to: ➤ Enable mining operations to generate a database with accurate historical flyrock measurements resulting from their blasting operations ➤ Allow research teams to conduct scientific investigations into flyrock and the impact of various blast design parameters ➤ Generate point-cloud data to visualize blasts and flyrock in a virtual reality environment for training and education purposes. This paper summarizes a conceptual technique and preliminary fieldwork that was carried out to determine the technique’s feasibility and motivate further development. The results show conclusively that a modified photogrammetric technique is capable of capturing flyrock data for further processing and analysis. The data acquisition procedure can, at this point, be used to meet the first aim of the project, namely to gather a field database of historical flyrock generation. Further development of the technique is ongoing and it is envisioned that the scientific-based technique will provide a method whereby future flyrock studies will be comparable and that assumptions will be limited.

Apr 5, 2023

By S. Schutte, J. Pelders, F. Magweregwede, S. Ngobese, W. Botha

A people-centric, systematic approach that involves different stakeholders is considered necessary to successfully address changes associated with modernizing the South African mining industry. In this paper we discuss the inclusion of various South African mining industry stakeholders, particularly the employees, in equipment design and development processes. Data was acquired through a literature review and inputs from South African mining industry stakeholders through focus group discussions, interviews, and an online survey. It was suggested that some of the areas that could be improved by the inclusion of employees in the process are design quality, ergonomics, equipment acceptance by the workforce, safety, efficiency, skills, insight into context of use, and early design iterations and identification of problems. Among the issues considered to be barriers to inclusion were lack of buy-in, lack of resources, difficulty in manufacturers accessing the mines, difficulties in involving too many people, intellectual property concerns, and unavailability of mine employees due to their key responsibilities. Identified critical factors in the equipment design and development process include iterative design, stakeholder participation throughout the process, needs analysis through consultation, and change management. A guideline was developed for worker inclusion in equipment design and development that could be used by the South African mining industry. The study identified several economic and social benefits of including stakeholders in the early stages of design and development. It is recommended that the South African mining industry considers using the developed guideline.

Apr 5, 2023

By J. W. van der Merwe, A. M. Maleka, V. J. Matjeke

The 10.5–12.5% chromium, titanium-stabilized ferritic stainless steel (FSS) plates used to fabricate coal line jumbo wagons is susceptible to undesirable chromium depletion during welding with 309L electrodes. In this work, non-stabilized FSS plates were welded using 308L welding electrodes. The objective was to weld non-stabilized FSSs using gas-shielded metal arc welding without causing sensitization. The welding process was optimized by lowering the heat input. Sensitization was assessed by rigorous etching techniques and anodic polarization scans. In addition, the welds were subjected to tensile testing, and their fracture surfaces were examined for intergranular cracking. Non-stabilized FSS plates were successfully welded using gas metal arc welding.

Apr 5, 2023

By M. Handley

The Witwatersrand Basin contains by far the most gold ever found, and has hosted mining from its discovery in 1886 to the present. For many years, South Africa was the world’s largest producerof gold, nearly all of which came from the Witwatersrand. Since 2000, South Africa has fallen back several positions because of declining gold output. There are many complex and varied reasons for this; however, declining gold resources in the Witwatersrand Basin are not one of them. As far as the author knows, there are no qualified estimates in the literature of gold remaining in the Witwatersrand Basin. This paper collects mining data from several sources and then, using elementary gold accounting and error analysis, attempts to estimate the amount of gold remaining in the Witwatersrand Basin. It includes gold underground in existing and defunct mines, in evaluated ore resources, and dumps. Compiled data from 1887 to 2019 shows that 50 200 t of gold were produced by Witwatersrand mines, while it is estimated that a further 48 100 t remainunderground, both within and outside of mine leases, a further 1 600 t on the surface in tailings,and an unknown amount in rock dumps. Nearly all of this gold will remain inaccessible with current mining methods, and major technical developments in mining will be necessary before any of the gold can be categorized as a code-compliant resource or reserve. To win this prize, the mining industry will have to rethink its approach to mining, both in old mining leases and in the unexploited ground. It will also have to find effective means of preventing gold theft and informal mining, which are on the rise.

Apr 5, 2023

By A. Mainza, A. Mabentsela

A finite element model of a 5 m radius DC arc ilmenite furnace in idling mode was used to test the notion that the slag solidifies when it comes in to contact with colder pig iron, thus constituting the initial step in the formation of solid slag at the slag-pig iron interface. It was found that a slag that is 150°C hotter than the pig iron does not solidify at the interface. The 150°C temperature difference between the slag and pig iron is a result of solid slag at the slag-pig iron interface, not the other way around as suggested in the literature. Calculations show that the thickness of the frozen slag at the slag-pig iron interface is 1.7 cm for the furnace used. It is proposed that slag solidification begins with the slow co-current flow of molten slag and pig iron in the outer parts of the furnace. This provides enough time for molten slag to interact with molten pig iron without solidifying. As the reduction products form due to reduction of the slag by carbon in the pig iron, the slag solidifies. This approach negates the need for the slag to solidify by merely coming into contact with an inherently colder pig iron. Making use of a low thermal contact conductance between the slag and pig iron was found to be sufficient to numerically capture the presence of solid slag at the slag-pig iron interface and to preserve the 150C difference between the slag and pig iron phase.

Mar 30, 2023

By T. Nghipulile, G. Dzinomwa, J. Kurasha, M. Amwaama, B. Mapani, K. Maweja, T. E. Moongo

The effect of mineralogy on the grindability was investigated using three copper ores – two sulphides and one oxide. The dominant copper minerals were identified by optical microscopy and mineral chemistry derived from SEM-EDS analysis. The sample designated sulphide 1 was borniterich, sulphide 2 ore was mainly chalcopyrite, and the oxide ore was predominantly malachite and minor azurite. The gangue minerals were identified using semi-qualitative XRD analysis. Sulphide 1 contained more than 80% (w/w) of quartz compared to about 70% in the other two ores. The Bond work indices were 13.8, 21.6, and 17.3 kWh/t for sulphide 1, sulphide 2, and oxide ore respectively. This suggested that the chalcopyrite-rich ore is the hardest, while the malachite-rich ore has intermediate hardness, and the bornite-rich ore is the softest. The brittleness indices of the ores were calculated using the chemical composition of the gangue, and a good correlation between brittleness indices and Bond work indices was observed, which highlights the importance of the gangue composition in determining the fracture behaviour of the ores. There is scope for further investigation into the relationship between ore mineralogy and comminution behaviour using other breakage characterization techniques.

Mar 30, 2023

By R. Nyirenda, B. Besa, N. Dzimunya

The accurate estimation of peak particle velocity (PPV) is crucial during the design of bench blasting operations in open pit mines, since the vibrations caused by blasting can significantly affect the integrity of nearby buildings and other structures. Conventional models used to predict blast-induced vibrations are not capable of capturing nonlinear relationships between the different blasting-related parameters. Soft computing techniques, i.e., techniques that are founded on the principles of artificial intelligence, effectively model these complexities. In this paper, we use the random forest (RF) algorithm to develop a model to predict blast-induced ground vibrations from bench blasting using 48 data records. The model was trained and tested using WEKA datamining software. To build this model, a feature selection process using several combinations of Attribute Evaluators and Search Methods under the WEKA Select Attributes tab was performed. The correlation coefficient of the actual data and RF model-predicted data was 0.95, and the weighted average of the relative absolute error (RAE) was 10.9%. The RF model performance was also compared to the equivalent-path-based (EPB) equation on the testing data-set, and it was seen that the RF model can effectively be used to predict PPV. The study also demonstrates that the EPB equation is a suitable empirical method for predicting PPV.

Mar 30, 2023

By R. H. Matjie, J. R. Bunt, A. C. Collins, J. C. van Dyk, C. A. Strydom

South African coal discards derived from feed coal for thermoprocesses pose human health and environmental problems and incur high disposal costs. These issues need to be resolved. Coal fines and a coal fines/ K2CO3 blend were combusted at 700°C in a laboratory rotary kiln to produce ashes containing metakaolinite (Al2O3.2SiO2) and illite/muscovite/orthoclase amorphous materials. The blend ash and coal fines ash, containing 18% K2O and <1% K2O respectively, were leached with water and with 1 M and 8 M NaOH at 80°C for 4 hours using a solid to liquid ratio of 1:5. The 8 M NaOH leach yielded 17% sodalite (Na8Si6Al6O24(OH)2) formation through SiO32-, Al(OH)4-, and Na+ precipitation in the leached ashes. The 1 M NaOH leach yielded no sodalite. X-ray fluorescence analysis of the original ashes and leached ashes showed low Al and Si dissolution efficiencies due to the sodalite formation. Sequential NaOH leaching resulted in the highest potassium dissolution efficiency of 89% for the ashes of fines/K2CO3 blend and 59% for the fines. It may be possible to increase the ash-derived zeolite concentration by using even higher NaOH concentrations (12–20 M NaOH) during leaching. Sodalite derived from coal ash could possibly be utilized for water purification in industrial applications, or as molecule separators. The potassium-containing liquid could be used in fertilizer manufacture.

Mar 30, 2023

By S. M. H. Ahmadi, K. Barani, F. Moradpouri, R. Ghaedrahmati

As exploration is time-consuming, costly, and risky, determination of the erosion surface of a metalliferous deposit before geophysical surveying and exploration drilling might be very helpful. Geochemical haloes can be used to determine whether the erosion surface is supra-ore or sub-ore and thus reduce the risk of exploration operations. The aim of this investigation is to determine the erosion surface of the North ROK porphyry deposit (NRPD) in northwestern British Columbia in Canada using linear productivity (LP), which is the content of an element defining the halo multiplied by the width of the halo. A total of 2045 soil samples from the B horizon were analysed using ICP-MS for 36 elements, including Cu, Mo, Pb, Zn, Au, As, Ag, Ni, Co, Fe, and Mn. The data-set was snalysed to obtain the statistical parameters and the elements Cu, Mo, Pb, and Zn were chosen to calculate the linear productivity and the total linear productivity. These four elements were modelled using probability plots to identify and separate subpopulations in terms of anomalous haloes and background, including the threshold values of each subpopulation. The results of the probability plot modelling and thresholds values were then used to map the distribution of each element in a GIS to calculate the linear productivity. The total linear productivity indicated that the erosion surface is supra-ore. Finally, a 3D orebody model of the Cu, Mo, Pb, and Zn distributions was constructed using borehole data and used to validate the results.

Feb 6, 2023

By A. M. Venter, C. W. Siyasiya, N. K. K. Arthur, M. G. Willemse, D. Marais

Although additive manufacturing is fast gaining traction in the industrial world as a reputable manufacturing technique to complement traditional mechanical machining, it still has problems such as porosity and residual stresses in components that give rise to cracking, distortion, and delamination, which are important issues to resolve in structural load-bearing applications. This research project focused on the characterization of the evolution of residual stresses in Ti-6Al-4V extra-low interstitial (ELI) additive-manufactured test samples. Four square thin-walled tubular samples were deposited on the same baseplate, using the direct energy deposition laser printing process, to different build heights. The residual stresses were analysed in the as-printed condition by the neutron diffraction technique and correlated to qualitative predictions obtained using the ANSYS software suite. Good qualitative agreement between the stress measurements and predictions were observed. Both approaches revealed the existence of large tensile stresses along the laser track direction at the sections that were built last, i.e., centre of the top layers of the samples. This in addition leads to large tensile stresses at the outer edges (corners) which would have the effect of separating the samples from the baseplate should the stresses exceed the yield strength of the material. Such extreme conditions did not occur in this study, but the stresses did lead to significant distortion of the baseplate. In general, the microstructures and spatial elemental mapping revealed a strong correlation between the macro-segregation of elemental V and the distribution of the β-phase in the printed parts.

Feb 6, 2023

By C. Ndou, T. Mothupi, H. Mandende

The technological advances in efficient, rapid, and non-destructive hyperspectral core logging systems for systematic mineral mapping have led to the discovery and exploitation of new mineral deposits Hyperspectral imaging in the long-wave infrared range has been recently used successfully to identify various phosphatebearing minerals (monazite, xenotime, and britholite), with limited work on apatite associated with mafic-ultramafic layered intrusions. In this study we investigate the effectiveness of a hyperspectral imaging (HSI) system with long-wave infrared (LWIR) bandwidthsto identify apatite in the Upper Zone of the Bushveld Complex. The accuracy of the HSI results was validated by mineralogical and geochemical data. The two apatite-enriched zones detected by HSI suggesting widespread development of apatite throughout the uppermost 600 m of the Upper Zone. The lower apatite-enriched zone is approximately 40 m thick, while the upper apatite-enriched zone is about 23 m thick, consistent with previous thickness determinations by traditional logging and analytical methods. Spectral mixing observed in the response of apatite is ascribed either to the common association of apatite and olivine in these rocks, or to differences between the spatial resolution of the hyperspectral image and the size of apatite grains. The VNIR-SWIR wavelength region did not show prominent spectral features of apatite. Nonetheless, HSI in the LWIR range is effective in mapping apatite and should therefore be considered as an exploration tool. This research advances our knowledge of the reflectance spectroscopy of REE-bearing minerals, which makes it easier to detect, identify, and quantify REE-bearing silicate minerals by HSI.

Feb 6, 2023

By A. McQuillan, N. Bar

Geotechnical models developed during the planning stages of open pit mines are three-dimensional so as to capture the spatial variation in lithological, structural, hydrogeological, and geomechanical conditions. Geological models that describe the lithological and structural (faulting and folding) characteristics of a deposit are always 3D. Likewise, boreholes and piezometers used to develop geomechanical properties and groundwater models are drilled at spatial offsets across the deposit to understand the lateral and vertical characteristics. Yet when geotechnical analysis is completed, often the three-dimensional geological, hydrogeological, and structural models as well as geometrically complex 3D mine designs for optimizing economic mineral recovery and overburden removal are simplified to two-dimensional sections. In this paper we demonstrate that this simplification can lead to the wrong failure mechanism being identified, analysed, and/or a conservative factor of safety being calculated and hence an over-estimation of slope stability. Through case studies we show how three-dimensional analysis methods are more suited to rock slopes, particularly those with anisotropic material strength, when singularities such as geological faults are present, and nonlinear slope geometry. When the same slopes are analysed in two dimensions, the failure mechanism calculated is often fundamentally incorrect. The case studies further reveal that the factor of safety calculated in three dimensions is not always higher than the two-dimensional factor of safety.

Feb 6, 2023