Search Documents

By R. C. A. Minnitt, A. Akpabio

This study presents a rigorous comparison between ordinary kriging and commonly used machine learning algorithms, those being, linear regression, support vector regression, decision trees, random forests (RF), and k-nearest neighbours for spatial interpolation of platinum grade estimates in a complex ore body within the Bushveld Igneous Complex. Using only X and Y coordinates as predictors, both ordinary kriging and machine learning models were evaluated at point and block supports under traditional and spatial block cross validation frameworks. While naive validation results suggested superior performance for k-nearest neighbour and random forest (R² = 0.92 and 0.86, respectively), these were revealed to be overly optimistic under spatial dependence. Spatial block cross validation results demonstrated substantial declines in model performance, with R² often falling below zero, particularly for decision trees and k-nearest neighbour, indicating strong overfitting and limited generalisability. Ordinary kriging exhibited more stable, albeit modest, performance under spatial validation, reflecting its strength in geostatistical interpolation when contextual geological variables are unavailable. The study underscores the critical importance of spatially aware validation in resource estimation and highlights that machine learning models constrained to spatial coordinates behave as interpolators rather than true learners of geological variability. Recommendations are provided for future work incorporating geological information to enhance predictive robustness.

Mar 30, 2026

By K. L. Morton, R. .. J. Connelly, P. C. Bouw

Uncertainties in the forecasting of the volume of groundwater likely to enter underground workings present difficulties for mine management in planning and costing the water-related activities of mining. This paper describes a technique that was developed to assess future inflows into a gold mine in the Orange Free State. The approach employed was based on the interactive operation of two separate computer models: a regional-flow model and an inflow model. As a first stage, a regional finite-element model of groundwater flow was set up and calibrated. Predictive runs were made to establish the influence of mining on the water levels from base-line conditions. This step was necessary since observations of the regional water levels were sparse. The predicted water levels were transferred to the inflow model, and the future flowrates were calculated. Both models were calibrated by trial and error until a satisfactory match was obtained between historical records and the values predicted by modeling. In this way, an improved monitoring system to record regional water levels and pressure heads within a mine was established. By the use of information, the initial model predictions can be refined as mining progresses. The models used are based on a non-linear relationship between the area mined and the inflows. A phased behaviour in the rates of inflow increase was noted. In conclusion, the modeling approach can provide mine management with guidelines for dewatering requirements. The confidence levels depend upon the records that have been collected during the mine's history.

Jan 1, 1988

By Q. Fu, Q. Wang, L. G. Wen, M. X. Jia, J. J. Zhao, J. W. Zhou

The BGRIMM process mineralogy analyser (BPMA) is a new automated mineralogy software platform developed in China. This system consists of the BPMA software platform, controlling a combination of commercially available scanning electron microscopy (SEM) and energy dispersive Xray spectrometry (EDS). The BPMA software controls SEM and EDS through their application programming interfaces, to form an auto-SEM and EDS system that automatically acquires BSE images and X-ray spectra, respectively, and classifies minerals by spectrum-matching. After automated measurement and mineral identification, various quantitative process mineralogical data (e.g., mineral composition, mineral abundance, mineral liberation and intergrowth relationships, mineral associations, particle and grain sizes, particle density distributions, calculated assay, elemental deportments, mineral grade recovery and elemental grade recovery, etc.) would be calculated quickly, accurately and automatically. The BPMA could be used for ore characterisation as a new advanced automated mineralogy technology. Keywords: BSE image; X-ray spectrum; automated mineral identification; quantitative mineralogical investigation; automated mineralogy

Jan 1, 2020

By L. Uotinen, M. Rinne, R. Mishra

Underground mining gives rise to geotechnical hazards. A formal geotechnical risk assessment can help to forecast and mitigate these hazards. Frequentist probability methods can be used when the hazard does not have many variables and a lot of data is available. However, often there is not enough data for probability distributions, such as in the case of new projects. The risk assessment is often subjective and qualitative, based on expert judgement. The purpose of this research is to present the use of Bayesian networks (BNs) as an alternative to existing risk assessment methods in underground mines by combining expert knowledge with data as it becomes available. Roof fall frequency forecasting using parameter learning is demonstrated with 1141 sets of roof fall data across 12 coal mines in the USA. The prediction is nearly identical for individual mines, but when multiple mines are evaluated it is difficult to find a single best fit distribution for annual roof fall frequency. The BN approach with TNormal distribution was twice as likely to fit the observed data compared to the Poisson distribution assumed in the past. A hybrid approach using BN combining multiple probability distribution curves from historical data to predict annual roof fall is proposed. The BN models can account for variability for multiple parameters without increasing the complexity of the calculation. BNs can work with varying amounts of data, which makes them a good tool for real-time risk assessment in mines.

Jun 1, 2021

By H. Ramazi, A. Amini

"This article is devoted to the introduction of a new technique of electrical resistivity data processing called residual resistivity (RR). We define RR as measured resistivity minus background resistivity. To determine the background resistivity, the data acquired from electrical resistivity measurements along a given survey line is evaluated, and then an equation is fitted to the data corresponding to a chosen measurement station as a function of current electrode spacing (or array length). The RR technique was applied to several synthetic models to compare the conventional resistivity inversion of each model with its RR-based inversion. A case study was carried out in a karstic area in Zarrinabad, Lorestan Province, western Iran, to detect the location and geometry of probable cavities by conventional resistivity inversion and RR-based inversion. The results showed that the anomalous zones are better highlighted in the RR-based inversion images in comparison with the conventional inversion images. In some cases, anomalies detected by the RR-based images were hidden in the conventional method. IntroductionThe geoelectrical resistivity method has been widely used since the early 20th century, and plays an important role in several fields such as groundwater and subsurface mineral exploration, geotechnical and environmental investigations, and archeological studies (Bayrak and Senel, 2012; Candansayar and Basokur, 2001; Fehdi et al., 2011; Hee et al., 2010; Osella et al., 2003; Wilkinson et al., 2010). The goal of geoelectrical resistivity surveys is to determine the distribution of subsurface resistivity by measuring the current-potential difference on the ground surface (Aizebeokhai et al., 2010). In the electrical resistivity method, an electrical current is injected into the ground by two electrodes, called current electrodes (AB), and the potential difference is measured between another pair of electrodes, the potential electrodes (MN). Several methods of interpretation are available for application to electrical resistivity data, the simplest being the graphical interpretation of the apparent electrical resistivity pseudo-sections along each survey line. A second method, electrical resistivity inversion, has been developed to create relatively accurate two- and threedimensional computational resistivity models of subsurface sections (Loke and Barker, 1996; Niwas and Mehrotra, 1997; Nordiana et al., 2014; Oldenburg and Li, 1994; Papadopoulos et al., 2009; Yilmaz and Narman, 2014). In the inversion method, models of subsurface objectives are mathematically inverted to reach optimal solutions subject to prescribed objective functions, constrains, and convergence criteria. A common problem when using the electrical resistivity inversion method is that the inversion model of an electrical resistivity profile may be non-unique because different inhomogeneities in the investigation media can result in the same electrical response (Smith and Vozzof, 1984)."

Jan 1, 2016

By G. Song, X. Liu

The global coal industry has changed significantly in recent years. In this paper we analyse the coal production and consumption, imports and exports, and coal prices in China and other major coal-producing countries during the past four decades. The analysis shows that global coal production and consumption have declined over the last three years, in striking contrast to the continuous growth for most of the past 35 years. Australia, Indonesia, and Russia have become the world’s major coal providers, while China, India, and the USA, and Germany are the major coal consumers. The world’s, and China‘s, coal industry can be classified into four stages of development: the stable growth phase, stable phase, rapid growth phase, and high-quality phase. Coal prices rose in the world’s four major coal markets in 2016 after four consecutive years of falling prices. The elimination of overcapacity and production-optimizing operations in China could assist in achieving a demand and supply balance.

Aug 8, 2024

By M. Z. Emad, Y. Majeed, G. Rehman

Massive deposits of rock salt are mined in the Salt Range area of Pakistan. The main export products of rock salt include salt lamps, blocks, and tiles, while local consumption includes the chemical industry and domestic use. The profitable production of salt for value addition depends on the quality, size, and shape of the block. Salt processing for value addition is easier when very few cracks are present in the blocks. The current salt mining practice involves drilling and blasting, which introduces undesired fractures in the material. This paper focuses on the extraction of rock salt blocks by applying conventional cuboidshaped block mining methods used in quarries. The techniques include wedges and feathers, expansion chemicals, and blasting with low-yield explosives. Experimental work involved both laboratory-scale and in-situ field testing. Laboratory experiments assisted with the determination of splitting force, loaddeformation curves, and other rock mechanics parameters. Regression analysis proposes a relationship for the estimation of radial strain from the load in the case of the wedges and feathers method. The splitting force obtained from the laboratory tests was used to confirm the accuracy of the already published empirical relationship. The results of the laboratory as well as in-situ tests showed that the wedges and feathers technique is the most suitable method for mining rock salt blocks of the desired size, shape, while minimizing cracking.

Mar 2, 2022

By A. McQuillan, N. Bar

Geotechnical engineers are routinely tasked with advising suitable stand-off distances below highrisk sections of slopes in open-pit mines that are identified to have potential to deform or collapse. Accurate prediction of failed material runout can mean the difference between continuous safe mining and unwanted high-potential incidents that result in loss of production, equipment damage, injury, or loss of life. This paper updates previous empirical relationships presented by the authors for estimating the volume and runout distance of excavated slope failures, in an open-pit mine operation, using slope geometry as the primary predictor. Cases are sourced from varying slope geometries (fall heights up to 385 m, slope angles up to 80°) and a range of commodities (iron ore, coal, nickel, gold, copper, boron, and limestone), excavated in sedimentary, banded-sedimentary, epithermal, and copper-porphyry deposits, across all six inhabited continents. Analysis of these cases identified positive correlations between slope height and runout distance, and slope height and failed material volume. In general, negative correlations were identified between Fahrböschung angle and slope height, and Fahrböschung angle and failed material volume; however, significant scatter is observed in these datasets. A definitive relationship could not be derived comparing Fahrböschung angle with failed material runout. Slope angle was also found to be a poor indicator of runout. Of the parameters analysed, slope height (i.e., fall height) was found to be the simplest and best predictor of runout distance. This paper presents new charts for predicting failed material runout distance for rock slopes. Relationships are defined for structurally and rock-mass driven slope failure mechanisms at average, 75%, and 95% prediction intervals

Jul 10, 2024

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 J. van Zyl, K. B. Le Bron, L. J. Gardner

Pillar design is a crucial aspect of underground mining engineering, as it directly impacts the safety, stability, and overall effectiveness of mining operations. In this paper we present a method to calculate pillar stability based on the strain criterion and pillar load inversion concept, which complements the empirical pillar strength formulae. Those formulae do not account for certain parameters that can affect the stability of the pillars, e.g. weak layers, changes in stress orientations due to mining, orebody dip, mining layout complexity, and the influence of regional stabilizing pillars. The approach presented here is not intended to replace the empirical method; rather, we suggest using it as the initial step in the pillar design process. This should be followed by iterative numerical modelling to study pillar behaviour, define pillar stability, and optimize the pillar design by considering site-specific and representative rock mass properties and criteria.

Jun 10, 2024

By E. R. Leeman

INTRODUCTION The first references in the literature to the so-called C.S.I.R. 'doorstopper' strain cell were made in 1964 1,2,3,4. These were the climax of nearly two years of experimental and developmental work with the instruments and the operating techniques involved. Since that time numerous field measurements have been made with 'doorstoppers' and they have been developed to a stage where they are in commercial production and are being supplied to users throughout the world. No detailed description of the latest 'doorstopper' equipment has, however, appeared in print and this paper is written to bring the literature up to date in this respect. The triaxial strain cell was developed three years ago 5, 6, 7 and was first described at the First Congress of the International Society for Rock Mechanics in Lisbon in 1966 8. Since that date laboratory tests and-particularly important-field tests 7' 9 have demonstrated the validity and practicality of the instruments and techniques. These instruments have also been developed to a stage where they have been released for commercial production. In this form they are very different from those described two years ago and a full description of them at this stage seems to be justified. THE INFORMATION SUPPLIED BY 'DOORSTOPPER' AND TRIAXIAL STRAIN CELLS Both the 'doorstopper' and triaxial strain cells were designed specifically to determine the absolute stress in rock using an overcoring stress relieving technique. They were not designed to measure changes in stress. They may, however, be used for the latter purpose provided the glue which is used to stick the strain gauges to the rock is known to possess sufficiently stable strain-time characteristics as not to affect the accuracy of the results during the time the measurements are made. The 'doorstopper' was designed primarily to measure the major principal stress in situations where its direction and those of the two other principal stresses were either known or could be assumed. For example, at great depth in undisturbed virgin ground, it is reasonable to assume that the maximum principal stress is vertical and that the other two lie in a horizontal plane. The borehole in which the 'doorstopper' is installed would thus be drilled in the horizontal plane parallel to one of the minor principal stresses, the direction of which would either have to be assumed or determined as described later. The stresses measured would be the vertical major principal stress and the principal stress acting in a horizontal plane normal to the borehole axis. It is possible, as will be shown later, to obtained the complete state of stress using 'doorstoppers'. In this event, it is necessary to make measurements in three boreholes drilled in any three known directions to each other and relative to a system of three orthogonal axes. The triaxial strain cell was designed to determine the complete state of stress in a single borehole drilled in any direction in any stress field. At first glance it would appear to supersede the 'doorstopper'. This is not necessarily true since there are

By M. Royle, A. Bloem, R. Uken

Rock mass damage induced by mining increases the fracture frequency and joint aperture width in the pit walls and rock mass behind open pit mine slopes. The depth and intensity of this mechanical damage relative to the pre-mining state is not well understood and yet is a major factor in the development of slope stability models and operational guidelines. The induced damage can be partitioned into domains or zones. The blast damage zone is essentially a ‘free flow’ rock mass volume, which transitions quite rapidly into the excavation damage zone, which extends beyond that. In both zones, transient high in situ water pressures are often linked to unstable slope conditions. Whereas rock damage can enhance passive depressurisation in the rock mass, high infiltration rates on the other hand, as one could expect during rainstorm events, will decrease slope stability within the dilated pit rock mass, especially in the blast damage zone-excavation damage zone, by increasing pore pressures for short periods of time. Ice jacking and freshet are additional factors to consider in (sub) arctic environments and can be exacerbated by unfavorable rock fabric orientation. On a local scale, brittle deformation damage zones, in the form of faults, can contribute to slope instability, or enhanced stability because of slope depressurisation, depending on their location and orientation relative to the pit shape and the water sources. In this study, the shapes of the damage zones and their anticipated impact on infiltration and/or passive depressurisation are considered for different pit morphologies, rock mass conditions, in situ stress, and hydrogeological settings. Some unusual hydrological situations, like mining towards rivers and lakes and the potential effects on pore pressures, are also considered. This investigation has resulted in the development of a method to adjust initial hydraulic conductivity values used to estimate transient pore pressure conditions found in rock masses exposed in different pit morphologies and stress regimes as well as discontinuities’ orientation relative to open faces within the open pit mine. The changes in K values and pore pressures can then be used to determine hydrogeological influence more accurately on future pit slope stability.

Sep 1, 2025

By C. Ntsele

Stirred milling technologies have firmly established themselves in concentrator circuits as energy-efficient alternatives to tumbling ball mills. High and continued increases in electricity costs, coupled with the ever-increasing power consumption by primary grinding circuits due to declining head grades of ore-bodies, provides a compelling case for these technologies. Additionally, the ability to successfully deploy stirred milling in coarse grinding applications extends the range of use in concentrators, thus increasing the opportunity matrix for energy optimization in a flowsheet. This paper therefore details methodologies and strategies that can be applied to evaluate the optimum stirred milling solution for different concentrator applications. These methods and strategies will focus specifically on the two major stirred mill machine types; namely, gravity-induced (VERTIMILL? grinding mill and Nippion-Eirich tower mill) and fluidized (SMD, Isamill?, Knelson-Deswik mill, etc. ). Finally, case studies from circuits across industry and laboratory tests will be used to illustrate the energy benefit from the use of stirred mills.

Jan 1, 2012

By T. Song, X. Wang, X. O. Xia, X. J. Luo, B. Chen

Interparticle breakage is the main crushing mode of an inertia cone crusher according to its working principle and feeding condition. In this paper the breakage behaviour of an inertia cone crusher is investigated, based on energy controlled compression combined with the vibration theory, rather than geometry controlled compression (form conditioned) which plays a major role in conventional eccentric cone crushers. The effect of feeding and operating conditions on breakage behaviour was described quantitatively from the perspective of energy consumption and product size distribution. The influence of feeding conditions, discharge size settings and vibration exciter speeds on the crushing process and product quality has been studied according to the results of laboratory crushing experiments. The results of the experiments show that the particle size of broken products becomes finer with the decrease of a discharge gap within the gap range selected by the experiment. The product size increases significantly when the gap approaches the upper limit of 6 mm. The variation of eccentric speed has a more prominent effect on reducing the product size compared with the change of discharge gap. The change in product size is not apparent when the eccentric speed rises to a certain level. The throughput of an inertia cone crusher is sensitive to the change of discharge gap, which increases linearly with the discharge gap and decreases with the eccentric speed. The change of feeding conditions has little influence on the crusher throughput. Keywords: Inertia cone crusher, interparticle breakage, energy consumption, product size distribution

Jan 1, 2020

By D. N. Esfahani, K. Khazen, M. Mohammadi

Electrowinning is one of the nickel recovery processes from solutions obtained from the leaching of nickel sources. Every year significant amounts of nickel products are used in various industrial processes and daily life, such as catalysts or batteries, which are a valuable source of nickel for recovery. Nickel purity is also very important in the electrowinning process, because the amount of impurities in nickel such as sulphur, if greater than 100 ppm, make it unsuitable for use in catalysts and plating. Because sulfur is intrinsic to the sulphate electrowinning system, we needed to investigate parameters that can affect the quality and purity of nickel. In this regard, experiments were performed in static and circulating electrowinning processes. A monitoring and circulation system was designed and built to measure voltage, current and temperature as well as to circulate the solution in the system at a specified rate. It was then observed that nickel ion concentration, pH, using the circulation system and additives, were factors affecting the purity and quality of nickel, and in the best case the amount of sulphur in nickel was below 100 ppm. Keywords: Nickel, electrowinning, sulphur impurity

Jan 1, 2020

By M. Z. Emad, A. Khan, N. Abbas, K. G. Li

Probability-based empirical methods were employed as an alternative approach to predicting uncertainties associated with rock mass properties. The focus was on developing probabilistic spreadsheets to forecast rock mass classification indexes. Histograms were constructed to describe the best distribution in predicting rock mass properties. The developed models also offer utility in predicting the impact of discontinuities within the rock mass on rock strength and rock mass classification systems. Statistical analyses identified volumetric joint count, joint spacing, joint frequency, and rock strength as the most influential parameters. Moreover, the statistical analysis revealed varying degrees of correlation among different rock mass properties. While some properties demonstrated significant correlations suitable for modelling, others did not align well with any correlation model. The results highlight the need for a comprehensive approach to rock mass characterization, considering multiple factors beyond volumetric joint count. Geological complexities, including tectonic activity and weathering processes, may obscure direct correlations. These results emphasize the importance of empirical modelling and detailed site investigations for accurate assessment of rock mass quality and stability in the Himalaya.

Jun 10, 2024

By A. García-Gutiérrez, A. Aragón-Jerónimo, R. Salazar-Mendoza

A theoretical analysis of a two-region solid-liquid flow system in horizontal pipes is presented, which considers a stationary bed of drill cuttings as a porous medium (-region) below of a two-phase dispersed flow (-region).The method of volume averaging was applied in this study to derive a rigorous mathematical model in which each variable is precisely defined. This model includes the volume-averaged transport equations for both the two-phase dispersed flow and the porous medium regions and terms from a macroscopic forces balance. The solution of the two-region model allowed evaluation of the behaviour of pressure gradient as a function of velocity, total volume fraction of cuttings and the relationship of the height of the stationary bed and pipe diameter. The solution is based on a backward finite difference explicit scheme. The simulated physical system is a pipe of 4.135 m in horizontal length and 0.0508 m in diameter. A one-dimensional, mesh-centred grid, consisting of 100 nodes, is used. The numerical results were compared with experimental data on slurry flows and a good agreement was found. Keywords: two-phase, cuttings transport, solid-liquid flow, drilling modelling, volume averaging.

Jan 1, 2007

By T. C. Maswanganyi, G. L. Smith, R. C. A. Minnitt

This paper identifies the strategic gaps and opportunities for the South African mining industry’s contribution towards advancing sustainable development. Mineral resources companies are confronted by a range of complex social performance challenges, which largely manifest in sporadic community and labour unrest, undermining the sustainability and profitability of these companies. Integrating supply chain management, strategic long-term planning and sustainability planning is identified as a strategic gap that can contribute towards the sustainability of mining companies, if effectively implemented. Contributing towards sustainable development by the South African mining industry is not merely a voluntary requirement on the industry but a legal requirement in terms of the Minerals and Petroleum Resources Development Act No. 28 of 2002. This Act is legislation aimed at advancing socioeconomic development in communities that host mining operations and in areas that provide labour to the South African mining industry. In contributing towards the advancement of sustainable development in mining host communities, not only will mining companies fulfil a legal requirement in terms of the Minerals and Petroleum Resources Development Act, but they will also prevent negative social and economic impacts in those communities. This paper identifies innovative and cost-effective ways of advancing sustainable development by integrating strategic supply chain management, strategic long-term planning, and sustainability planning.

Dec 15, 2025

By S. R. Grobler, S. Bilessuris, M. E. Greenway

Ensuring that all the head ropes on a multi-rope friction winder share the load equally has been a design and maintenance issue since they were first developed. Accurate machining of the rope grooves is necessary to meet this objective. Early analysis of the change of rope loads for ropes in grooves that are mismatched in depth treated the grooves as rigid. Groove depth tolerances deduced from this analysis are stringent and difficult to achieve in practice. Test work indicated that the inherent flexibility of groove lining materials alleviates the load-sharing problem. However, no analysis of rope load changes due to mismatched flexible grooves has been published. In this paper we develop new equations for the rope load variation for flexible grooves. By allowing flexibility to reduce to zero, the corresponding rigid rope groove equations are obtained. New criteria for tolerable groove depth variations are developed from this analysis, which depend on the flexibility of the groove lining material.

Sep 1, 2021

COLLOQUIUM AND GENERAL MEETING Developments in gold recovery processes A General Meeting and Metallurgical Colloquium on the above topic were held on Wednesday, 14th March, 1973, at Kelvin House, Johannesburg. Dr J. P. Hugo (President) was in the Chair. The colloquium was attended by 170 delegates, and was opened by the President at 9.15 a.m.