Search Documents

By P. Lenegan, R. Masethe, R. I. L. Ferreira

Underground mining in high rock-stress environments can induce unanticipated dynamic rockmass deformation, posing a serious risk to the workforce. Data collected by seismic networks are imperative for the design-as-you-mine process. Often an inappropriate assessment of the seismic hazard is done when unnecessary data are included in analytical techniques, particularly short-term time frames. Depending on the methodology or application of data analysis and interpretation, some of the information may not be relevant and even detract from the objective. By using the spatial event clusters associated with active workplaces to define the shape and size of the polygons for data selection for subsequent analysis, a problematic subjective component can be eliminated. The mining of tabular reefs in South Africa, be they gold- or platinum group metals, is very extensive, and underground workplaces are often widespread. Production stoppages are very costly, especially for marginal mines, so there is just cause not to evacuate an entire operation’s workforce following the occurrence of a large mining-induced seismic event. A rationale based on the damaging peak particle velocities (PPV) of the past seismicity can exclude the unaffected underground personnel and allow a quick orderly withdrawal where required. Micro-deformation may precede rockbursting, which can occur on-shift close to the workface. A mine-wide array of geophones suitable for the accurate location and analysis of seismicity may not be able to detect these higher-frequency emissions. Inexpensive but fit-for-purpose accelerometers, providing real-time, continuous measurement of acoustic rockmass deformation near the workplace, integrated through underground and surface communications networks, are shown to record ‘rock-talk’, clearly distinguishable from the instrument background noise, which can be in the order of hours leading up to an imminent damaging seismic event.

Jan 12, 2025

By T. Matlou, S. Myendeki, V. Midzi, T. Zulu, H. Mabotja, T. Mulabisana, B. S. Zulu, L. Hlatshwayo, T. Makhateng

The Council for Geoscience (CGS) performed a macroseismic investigation following the occurrence of an earthquake in the Boksburg area on 11 June 2023. The event had a magnitude of ML 4.4 and a focal mechanism that showed the event as having strike-slip motion along a northwest-southeast oriented fault. The event resulted in aftershocks, some of which occurred away from the causative fault in the east in an area without a mapped fault. However, the cluster of past events in that area suggests that there is an active fault. Due to the size of the main earthquake, a macroseismic survey was conducted to determine its impact on the community in the region. Analysis of observed data obtained from questionnaires resulted in the creation of 153 intensity data points. Their distribution showed high intensity values close to the earthquake source area. However, the investigation highlighted a number of anomalies, such as the higher than expected intensity values observed far from the event epicentre. It could be that these anomalies are caused by site conditions, for example, amplification of ground motion due to thick soil layers.

Jan 12, 2025

By S. Durapraj, R. T. Masethe, U. Maqina, A. C. Adoko

This study investigates seismicity and rockburst control within the Merensky Reef and UG2 orebody of an intermediate-depth platinum mine located in South Africa's Bushveld Complex. Utilising a mine-wide seismic monitoring system, over 1,900 seismic events were recorded between April 2023 and May 2024, revealing complex seismic patterns driven by mining-induced stress changes. The largest event, with a magnitude of ML 2.0, underscores the seismic risks associated with mining at increasing depths. Key findings include the spatial and temporal variability of seismic events, with high seismic activity observed near pillars and geological structures. A significant portion of events correlated with pillar and foundation failures. Moment tensor analysis identified dominant failure modes, such as stress-driven fracture closures and dynamic pillar failures. To mitigate risks, the study emphasises advanced support system redesign and rockburst risk assessment using the rockburst damage potential (RDP) method. The research highlights the importance of monitoring, geotechnical compliance, and adaptive designs in enhancing safety and reducing seismic hazards in underground mining operations. The study applied the RDP framework, classifying damage across R1 to R5 scales, where most instances fell into low-risk categories (R1 and R3), indicating minor to moderate damage. However, data limitations and the tabular stoping method's complexity led to inconclusive results, highlighting the need for more robust data and customised risk assessment tools. Despite these challenges, the RDP approach provided valuable insights into rockburst mechanisms and the effectiveness of existing support systems.

Jan 12, 2025

By J. Petersen, S. W. Robertson

Heap leaching is used extensively for the processing of low-grade copper, gold, and uranium ores. The hydrodynamics is often the most important factor that limits the rate of metal recovery. Over the past decades the author has performed physical and hydraulic testing of a large number of ore samples considered for heap leaching. The data were reviewed to investigate the effect of physical properties such as particle size distribution (PSD) and bulk density on the hydrology. The unsaturated flow through the bed is traditionally modelled with soil science models such as Brooks-Corey (BC) and Van Genuchten-Mualem (VGM), which relate the hydraulic conductivity to the degree of saturation, also known as the hydraulic conductivity function (HCF). However, most authors found a disconnect between the saturated hydraulic conductivity (KS) and the rest of the HCF, which is explained as resulting from a transition from capillary-to gravity-controlled flow. However, the author found no such transition, and found the flow to be governed by capillary flow throughout. Traditional BC or VGM fits were found to be appropriate, however it was necessary to modify the model to include a point of inflection, above which the hydraulic conductivity remains constant while the moisture content keeps on increasing. It was proposed that this is the result of unconnected or poorly connected pores, which fill with solution but do not create additional flow channels. The point of discontinuity was also found to correspond with the air entry point, measured independently from the air conductivity curves. On average, the air entry point corresponds to the experiential rule that the void saturation should remain below 65% for aerated heaps.

Dec 11, 2024

By J. Petersen, A. E. Lewis, J. Chivavava

Rare earth elements (REEs) have special properties that prompted their extensive use in high-tech applications. The demand for REE has, therefore, increased over the past years, resulting in the supply risk of the materials. Extraction from non-conventional low-grade ores, like easy-to-mine ion adsorption clays (IACs), by hydrometallurgical methods is being explored to supplement the supply of REEs. The extraction of REEs from low tenor IAC leach solutions requires enrichment using cost-effective methods such as precipitation. Therefore, the aim of this study was to understand the recovery of REEs from IAC leach solutions using various reagents for precipitation. In the first part of the study, recoveries of REEs from IAC leach solutions through direct precipitation using H2C2O4 and NH4HCO3, as well as antisolvent crystallization using C2H5OH, were predicted via thermodynamic modelling using OLI Stream Analyzer. Verification experiments were then conducted through the use of an agitated reactor. Simulation results showed that high yields of REEs were possible using each of the reagents, but large quantities of C2H5OH were required. Experimental results confirmed the high yields predicted from simulations. Precipitation with H2C2O4 was found to selectively recover REEs while rejecting Al. However, utilization of NH4HCO3 and C2H5OH resulted in the co-precipitation of Al. Furthermore, the yield of Y from antisolvent crystallization was significantly lower than the theoretical value. The product recovered from antisolvent crystallization consisted of well faceted, large crystalline particles, which filtered faster than the fine, amorphous REE carbonates recovered using NH4HCO3. REE oxalates were crystalline, larger than the REE carbonates and, overall, showed the highest filtration rate.

Dec 11, 2024

By B. Genc, V. J. Mahadeo, S. Ramparsad4, U. Govender

The study was intent on assessing the socio-economic impact of a zinc mine in South Africa through the application of people-centric digital tools using a case study research methodology. The global mining industry is undergoing significant technological shifts to modernize, improve operational sustainability, and align with the environmental, social, and governance (ESG) agenda. While the case study mine is a more modernized mine, continuous improvement around the ‘social’ pillar becomes critical, and to this end application of people-centered digital tools presented as the vehicle to understand community needs and impacts. Despite considerable annual corporate social investment by mines annually, the overall impact or social return on investment is marginally understood. Towards contributing to addressing this gap, digital ESG ‘social-related’ data gathering tools for mines decision-making were developed, including the training-needs assessment tool, community-social needs assessment, small, micro, and medium enterprise (SMME) communications and engagement tool, digital leadership competencies gap assessment tool, and change management blueprint, which further allowed for customization to the mining context. Hence, the novelty of this paper relates to applying digital tools to provide real-time information of the needs of a community, both socially and economically. The findings suggested that these digital tools enabled gathering of information for the establishment of regional and operational baselines towards enabling of social investment returns in relation to socio-economic transformation development interventions. These tools have been applied to a zinc mine context, with lessons and experiences obtained through this study, invaluable to the application of these digital tool in other mines in South Africa. The study was conducted in September 2023 across five host communities, surrounding the zinc mine, in the Northern Cape region, with 120 participants involved in the application of three relevant community-related digital tools, which included the training needs, community-social needs, and SMME communications and engagement tools. It was evident that the digital ESG tools could be used for data gathering for the mine’s informed decision-making and allowed for securing of insights on community needs in near real-time, and further proved practical for community use.

Dec 11, 2024

By F. Ntuli, T. V. Nyangadzayi

Chromium green is an essential metal with various applications in the paints, construction, electronics, and aerospace industries. The main source of chromium green is the carbo-thermic reduction of chromite ores, however, chromite ores are being heavily depleted on a global scale, hence the need to find alternative or secondary sources for the chromium. In order to meet the global demand for chromium there is a need to implement principles of the sustainable circular economy to keep chromium within the material loop. Ferrochrome slag and stainless steel slags are secondary resources that are potential sources of chromium. In this work, recovery of chromium from ferrochrome slag was investigated. The ferrochrome slag was alkaline leached using sodium hydroxide- sodium hypochlorite solution to reduce the use of potent reagents that have a potential negative environmental impact. The leachate produced by oxidative alkaline leaching of ferrochrome slag contained hexavalent chromium, which was used to recover the extracted chromium through a two-stage process involving reduction and precipitation. Reduction of hexavalent Cr to trivalent Cr was investigated using anhydrous FeSO4 at concentrations of 5 g/L, 10 g/L, and 15 g/L. Precipitation was conducted by dosing the reduced leach solution with 1 M sodium hydroxide at dosages of 5 ml/L, 10 ml/L, and 15 ml/L. The experiments were designed using the Taguchi (L9) orthogonal array and analysed according to the ANOVA, general linear model. The optimized conditions for the reduction step were ferrous sulfate dosage of 5 g/L with a pH of 12, and for the precipitation step, 1 M sodium hydroxide at a dosage of 15 ml/L with a setting time of 90 minutes and pH of 7. Physical changes were observed during the two-step chromium recovery, with solution changing from a yellowish colour to dark green colour indicating completion of the reaction. The chromium recovery was found to be 93.3%.

Dec 11, 2024

By J. Petersen, T. Moyo, D. Maharaj

This study investigates the effect of ammonium thiosulfate on the leaching of gold from the contacts of unpopulated cut electronic printed circuit boards (PCBs). Studies in the literature show ammonium thiosulfate to be a viable option in comparison to the commonly used but hazardous cyanide leaching route. However, the simultaneous dissolution of Cu interferes with the Au leaching reaction, often limiting the extraction to moderate levels. To optimize the system, an ammoniumthiosulfate system containing 0.5M ammonia, 0.5M ammonium thiosulfate, and Cu(II) sulfate at varying concentrations was studied as an alternative hydrometallurgical option for the recovery of gold from PCB. Tests were done with unpopulated standardized model PCBs (cut to 15 x 30mm pieces). Gold extraction of 47.3% was achieved after a 24-hour leaching experiment with a 0.02M copper background concentration. Further experiments focused using different initial background Cu concentrations (0M – 0.1M) to determine its effect on both gold and copper extractions. An initial Cu concentration of 0.045M resulted in the highest gold extraction, with 93.7% extraction after 6 hours of leaching. Further analysis of the data shows that the kinetic rate of gold leaching is correlated in a first-order dependence with the actual copper concentration in the solution, which increases with the concomitant dissolution of copper and gold from the PCBs. A lower threshold of 0.017M and an upper threshold of 0.15M copper in solution were determined, in between which gold extraction shows this first-order rate increase. The upper threshold is potentially linked to the onset of copper precipitation, which blinds the metal surfaces, slowing further leaching. The findings of this study suggest that in an industrial recovery operation of gold from waste PCB using the thiosulfate method, copper concentrations need to be carefully monitored, and the process should preferably operate with uncut or coarsely cut untreated boards to minimize exposure to copper circuits.

Dec 11, 2024

By K. C. Sole, M. Makonese

Bushveld Belco is commissioning a plant for production of electrolyte used in vanadium redox flow batteries. Owing to the source material, the electrolyte is contaminated with Fe, at concentrations up to 140 mg/L, which negatively affects battery performance. This study examined the use of ion exchange to reduce the Fe concentration to the target value of < 100 mg/L from a feed solution containing 100 g/L V in 4 M H2SO4. Four resins with possible application were identified: Puromet MTS9570, Puromet MTX7010, Puromet MTS9500, and Puromet MTC1600H. Puromet MTS9570 outperformed the other three resins in all batch experiments with respect to loading capacity and selectivity over V, and was evaluated for Fe removal in fixed-bed column loading tests. Based on these data, preliminary sizing and design of a full-scale column was undertaken. Owing to the difficulty of eluting this resin, a philosophy of discarding the resin on loading to full capacity was considered. These data revealed that, despite the high selectivity of this resin for Fe over V, it would be too expensive and impractical to use this resin for full-scale operation under such conditions. It is recommended to remove iron earlier in the flowsheet prior to generation of the 4 M H2SO4 electrolyte.

Dec 11, 2024

By F. Dinter, C. McGregor, Q. G. Reynolds, L. Hockaday

The necessity of reducing greenhouse gas emissions to limit the impact of anthropogenic climate change means that all industrial processes should investigate their options towards decarbonization. Iron, steelmaking, aluminium, and cement production produce the most industrial emissions. Although most research focuses on decarbonizing these industries, minerals processing industries must be decarbonized to achieve net zero emissions goals by 2050. Incorporating renewable energy sources and avoiding the combustion of fossil fuels are two ways to reduce greenhouse gas emissions. This paper reports on the results from a dynamic process model developed to investigate the feasibility of concentrating solar thermal pretreatment of manganese ores to pretreat carbonate-rich manganese ores for increased ferromanganese smelter productivity and reduced greenhouse gas emissions. The results indicate that pretreated ores reduce the energy requirement for smelting significantly for some ores and lowers total carbon dioxide emissions by 13 to 19% compared to the traditional smelting route. A high-level economic evaluation shows that solar thermal treatment is financially viable for low cryptomelane, high braunite ore that is common in the Kalahari Manganese Field.

Dec 11, 2024

By L. J. Su, J. Y. Zhang, P. W. Gu, J. F. Xu

Surface tension of mold flux is a key factor affecting optimum casting conditions. A thermodynamic model was developed to determine the surface tension of mold flux containing fluoride based on the ion and molecule coexistence theory of slag structure and Butler’s equation. The relationship between composition and surface tension was investigated by this model. Results indicated that the calculated values showed good agreement with literature data, and the average error was 8.19%. The surface tension value of the CaO–Al2O3-based mold flux was larger than that of the CaO– SiO2-based mold flux. The surface tension in a multi-component system decreased with increasing contents of SiO2, CaF2, and Na2O, which are considered surface-active components. This trend was more significant in the CaO–Al2O3-based mold flux.

Dec 11, 2024

By O. S. L. Bruinsma, M. J. Matsela, D. J. van der Westhuizen

The free ligand (FL) model is an intuitive model to calculate the equilibrium concentrations in solvent extraction (SX) based on the distribution coefficient (Kd) at infinite dilution and the fraction of free ligand available for extraction. In this work, the FL model is validated to account for co-extraction of Ni and for designing the extraction and scrubbing sections of Co/Ni separation in a base metal refinery. Cyanex 272 is employed as a selective Co extractant in the process to yield high purity Ni. To enhance the loading capacity of the extractant, 50 mol% of NaOH is pre-loaded into it. Experimental and FL modelling results show that four equilibrium stages are required to purify the 99% Ni feed to a 99.97% Ni raffinate. The 53.59% co-extract with 26.76% free ligand is purified with a reflux scrub of 30 g Co/L in one equilibrium stage to 99.01%, Co in the scrubbed extract. The FL model as well as the Co/Ni/Na-analysis of the aqueous phase confirm that in the scrubbing section, Ni in the extract is mainly replaced by Co. Finally, in the one stage stripper the Co concentration in the scrubbed extract is increased from 5.6 g Co/L in the extract to 30.5 g Co/L that is partially used as reflux scrub liquor.

Dec 11, 2024

By A. Doroudgar, H. Noferesti

Unstable rock failure or rock bursting is one of the most dangerous problems in deep hard rock mines. This type of failure, that occurs uncontrollably, is associated with a significant release of strain energy from the surrounding rock. In this study, the numerical simulation of unstable rock failure was carried out using the three-dimensional explicit Lagranian finite-volume finite difference method to identify the onset and location of unstable failure events in the model. On this basis, a numerical simulation of mine pillars with width-to-height ratios of 0.5, 2 and 3 was performed. A criterion was developed, based on the calculation of the elastic strain energy, to determine the prone areas for the initiation of rockburst in the models. The results of this study showed that, in softer loading systems, due to large amounts of stored elastic strain energy, an unstable equilibrium is established, which leads to the rapid transfer of strain energy to the rock in the form of violent failure. It was also found that the unstable failure starts at the corners and free surface of the model. The new criterion developed in this study may help predict the occurrence of rockbursts in slender pillars and the design of pillars to control human and financial losses caused by this phenomenon.

Nov 27, 2024

By A. Çelik

Wheel loaders are widely used in many sectors due to their mobility, operational flexibility, and low costs. The literature is mainly focused on improving the hydraulic, mechanical, and electronic subsystems of these machines in terms of efficiency or developing strategies for controlling these systems. In this study, the effects of the loader operator and different working conditions on the loading performance were investigated. The study was carried out at the Can Lignite Enterprise (CLI), which is a major coal producer in Turkey. In the study, which lasted approximately four months, a total of 808 measurements were performed to investigate seven different scenarios. The effect of a single parameter on loading performance was investigated in the first six scenarios. The results showed that the biggest impact on loading performance was caused by the operator's experience with 63.4%. On the other hand, the worst and most ideal working conditions determined from the first six scenarios were tested in the seventh scenario. It was tested as a script. The test results showed that the loading performance was 285.5% more effective under the most ideal operating conditions compared to the worst operating conditions. In addition, in the study, the loading performances of five loaders with four different capacities were examined in terms of energy efficiency and it was determined that the Kawasaki 115 work machine, which has a higher model and lower engine operating hours, was more advantageous. In the study, various suggestions were developed in parallel with the research results. The most prominent and easiest to implement of these suggestions is related to the selection of the loader operator.

Nov 27, 2024

By M. Mokgohloa, T. R. Stacey

Boring of tunnels and shafts in hard rock under high stress conditions is becoming more common as mines and tunnels are developed at greater depths. Under the high stress conditions, fracturing of the often-brittle rock occurs in the walls, backs, and ahead of advancing faces of excavations. Fracturing can have a significant impact on boring activities: sidewall spalling, which affects machine gripper capacity, and fractures can also develop in the rock ahead of the face, leading to blocky rock conditions, which may have a significant effect on machine excavation. Fracturing may develop dynamically, leading to rock bursts manifesting in the form of strain bursts, which can be hazardous and destructive. In brittle rock, fracturing is commonly extensional in nature. The focus of this paper is on fracturing in the face of the excavation. Examples of such fracturing behaviour are described briefly in the paper. Numerical analyses were carried out to predict the initiation of extension fractures and their orientations ahead of excavation surfaces, and the resulting formation of rock slabs/plates, and the stability of these plates. The results should be beneficial for the evaluation of the conditions prior to, and during boring operations.

Nov 27, 2024

By M. I. Mabala, B. P. Watson, M. J. Kanda

A failed and stable Upper Group 2 chromitite pillar database was collected in 2007 as part of the PlatMine 1.2 pillar strength project. The database consisted of 177 pillars, of which 33 had failed. Most of the data were collected from various intermediate-depth mines near the town of Thabazimbi. Spans between rows of chain pillars were about 30 m. Pillar loads were estimated using pseudo-3D, elastic, boundary element modelling software. The pillar-load data were evaluated using a maximum likelihood regression analysis to determine strength parameters, following a similar approach to historically acceptable practice. The original formula was named the PlatMine Formula. This paper describes re-evaluations of the database, using a linear formula and several variations of the power formula. The formula providing the lowest standard deviation to the data in the database was considered the best representation of the pillars in the database. Through this rigorous process, a new formula for peak pillar strength was developed for Upper Group 2 chromitite pillars, which is unaffected by volume. The small standard deviation provided by this formula suggests a reliable relationship between strength and w/h ratio, in the w/h ratio range between 1.5 and 4.7. Accuracy, precision, and recall analysis confirmed the efficacy of this equation to the data in the database. Importantly, the new formula provided in the paper does not predict a significantly different strength to the current PlatMine formula, within the confines of the database.

Nov 27, 2024

By D. F. Malan, W. J. Theron

Practical mining aspects should be considered when conducting pillar designs for bord-and-pillar layouts. The current methodology for pillar design will result in increasing pillar sizes with depth. This affects the extraction ratio and will result in onerous ventilation requirements when cutting large pillars. A holistic approach, including all mining engineering requirements, is required to ensure that the rock engineering designs are optimized to ensure efficient mining operations and sustainable production. Bord widths should not only be a function of the rock mass ratings, but should also be selected to fit the specifications of the mechanized equipment. The use of a ‘squat pillar’ formula for hard rock is discussed in the paper and the formula based on the exponents of the Hedley and Grant pillar formula, is explored. The effect of abutments and geological losses on average pillar stress is also explored. These factors must be considered when designing layouts at increasing depths.

Nov 27, 2024

By J. Uys, P. J. Mokganya, J. Olwagen, R. C. W. Webber-Youngman

This article discusses the critical role of leadership in the technology adoption process within organizations, particularly in industries like mining, where innovation is essential for maintaining competitiveness and operational efficiency. It highlights the limitations of traditional technology adoption models, such as the Technology Acceptance Model (TAM), Diffusion of Innovations (DOI), and Unified Theory of Acceptance and Use of Technology (UTAUT), which focus on individual factors and overlook the strategic and organizational roles of leadership. The study proposes the Ambidextrous Leadership-Technology Adoption (ALTA) model, which integrates leadership styles such as transformational, transactional, adaptive, and ambidextrous leadership into the technology adoption process. The model emphasizes the balance between exploration (fostering innovation) and exploitation (ensuring operational efficiency). The findings suggest that transformational leadership fosters innovation by encouraging risk-taking and creativity, while transactional leadership ensures efficient technology integration. The article also identifies key barriers to technology adoption, such as organizational resistance to change, which can be mitigated through adaptive leadership that addresses employee concerns. The research emphasizes the importance of a supportive organizational culture and cross-functional collaboration, which are essential for successful technology adoption. The study concludes by offering practical recommendations for leaders to foster a culture of innovation, manage resistance, and balance short-term efficiency with long-term technological advancement.

Nov 27, 2024

By C. C. E. Pretorius, V. M. Mathoho, A. van Zyl, R. J. Mostert

The in-service degradation of the mechanical properties of steel components through the damage mechanism of high temperature hydrogen attack (HTHA), is a topic of concern in the refining and green hydrogen industry. This damage mechanism occurs in susceptible steels operating in environments at high temperatures and hydrogen pressures. The current investigation deals with the indirect monitoring of mechanical degradation via tracking of the swelling strain in affected structures. An autoclave with an AISI 316 shell was utilized to simulate accelerated HTHA damage at 550 °C and 46 bar for exposure times ranging from zero to 700 hours. The progress of the HTHA damage was tracked using encapsulated high-temperature strain gauges. The correlation between the swelling strain and mechanical property degradation was studied to develop a methodology for the continuous monitoring of embrittlement. The tensile sample orientation of the carbon steel plate was included as a variable, i.e., samples parallel, transverse, or perpendicular to the plate rolling direction were included, since it has been shown that the sample orientation influences the HTHA damage features. For the through-thickness orientation, and upon exceeding a threshold value of exposure time, all tensile properties were severely degraded, with values in the order of a 90 per cent reduction being observed. It was found that the degradation of carbon steel mechanical properties can be correlated to the swelling strain measured during exposure. The critical point for mechanical property degradation in the plate through-thickness orientation, whereafter a severe decrease in the ductility of the material occurs, was found to be in the order of 1% of the total swelling strain measured during exposure, equivalent to 65 microstrain. This threshold was found to be significantly lower than that expected from the literature, where limits of 400 to 1000 microstrain were postulated.

Nov 27, 2024

By D. F. Malan, P. J. Le Roux

Almost no work has been published on the strength of pillars in mines exploiting the LG6/LG6A chromitite bands in the Western Bushveld Complex. The strength of these pillars is unknown, and the hard rock industry still uses the Hedley and Grant formula. Numerical modelling, using inelastic constitutive models, may be of some value in estimating the pillar strength, but this approach is difficult and prone to errors as many assumptions are made. This paper explores the alternative approach of the back-analysis of actual LG6/LG6A mining layouts, using displacement discontinuity codes to simulate the average pillar stresses on a mine-wide scale. From this, and underground observations, a ‘minimum’ pillar strength (K-value in the pillar strength formula) can be estimated. The codes TEXAN and Map3D are both used in the paper to simulate the actual pillar shapes and sizes in an LG6/LG6A mine. The results are encouraging as the two codes produced similar APS values for a complex pillar geometry. An estimated K-value of 77 MPa was determined for the LG6/LG6A pillars in this particular area. This is more or less similar to earlier work done for the UG2 chromitite seam, but higher than the K-value of 39 MPa for a linear formula obtained by other workers for LG6/LG6A pillars. Caution should therefore be exercised before using this higher value for all areas. The effect of element sizes and explicitly simulating the effect of the surface for shallow pillars is illustrated in the manuscript using both codes. The need for actual pillar observations and an iterative design process, which cycles between modelling and observations, is also emphasized.

Nov 27, 2024