Search Documents

By Fangyu Liu

In the grand tapestry of scientific inquiry, there are individuals whose contributions resonate far beyond the confines of laboratories and lecture halls. It is with profound admiration and respect that we present this special collection of Mining, Metallurgy & Exploration, dedicated to celebrating the remarkable career of Dr. Patrick Taylor. Dr. Taylor, recently retired from the Colorado School of Mines as the Emeritus Ansell Distinguished Professor of Chemical Engineering, possesses over 45 years of expertise in mineral processing and extractive metallurgy. His remarkable journey has seen him traversing the globe, from Canada to Egypt, contributing to pyrometallurgy, hydrometallurgy, and mineral processing. Dr. Taylor’s portfolio encompasses lab work, pilot-plant projects, research, and process development across various metals, with an impressive array of authored papers, presentations, and nine patents to his name. A soughtafter consultant for over 20 companies and an esteemed expert witness, Dr. Taylor also led the Kroll Institute for Extractive Metallurgy at Colorado School of Mines. Under his guidance, more than 100 graduate students and post-docs have undertaken transformative research, while Dr. Taylor’s 45-year tenure as a university lecturer and the creation of 10 industryfocused short courses underscore his dedication to knowledge dissemination. His contributions extend to crafting professional engineering exam questions for a quarter-century and active participation in esteemed organizations like SME, TMS, ASM, and MMSA.

Sep 16, 2023

By Ewan Wingate, Ranjit Prasad, Yanjun Liu

The use of phosphorus by mankind is long established. From use in agriculture, foods, high tech electronics, and more recently in EV battery cathode production, one cannot escape its impact on today’s society. This paper will review and describe the circular journey of phosphorus through its value chain from the mining operation of phosphate ore through beneficiation into downstream chemicals production and finally into the EV battery cathode production space. This will include a review of phosphoric acid production from phosphate concentrate using wet process route and thermal route. The wet process route discussed will be via dissolution of phosphate concentrate by sulphuric acid. The thermal route will discuss thermal decomposition of phosphate concentrate using submerged arc furnace to produce yellow phosphorus followed by its combustion to produce phosphorus pentoxide and subsequent hydration to produce phosphoric acid. It will also discuss briefly other phosphate derivatives for specialty chemicals and recovery of phosphorus units from wastewater treatment plants, sewage sludge and incinerator ash. Finally, the paper will discuss the use of phosphorus in EV batteries in the form LFP cathode materials and its production. It will also briefly discuss the LFP battery recycle.

Sep 11, 2023

By Xianyang Qiu, Xiuzhi Shi, Chen Xie, Xiaofeng Huo, Yanhai Wang

In upwards medium-deep hole cut blasting, many confinement effects exist due to the small number of free surfaces, which is mainly provided by an empty hole. If the cut blasting effect cannot meet expectations, then subsequent blasting construction will be seriously affected, and large security risks to operators in the stope can be induced. The number of free surfaces is directly related to the layout of the cutting borehole. Therefore, it is necessary to optimize the layout of the first round of cutting boreholes. Theoretical analysis, numerical simulations, and field tests were used to study the formation process of cut blasting and optimize the layout of the first round of cutting boreholes. Based on the blasting parameters in the field, LSDYNA software was used to establish cut blasting models with four kinds of borehole layouts in the first round. The RHT constitutive model was used to simulate the rock-breaking effect. The numerical results show that the cut blasting effect is optimal when Nos. 3, 5, 7, and 9 are charge holes and Nos. 1, 2, 4, 6, and 8 are empty holes. In addition, the greater the number of empty holes, the larger the release space and the better the cut blasting effect. The uniformity of the empty hole distribution, namely, the uniformity of the release space, also affects the cut blasting effect. The field results indicate that a good cut blasting effect can be achieved by optimizing the borehole layout.

Sep 7, 2023

By Andrea Brickey, RICHARD AMOAKO, Akshay Chowdu

Production forecasts derived from the medium-term schedule represent the best path forward for underground mining operations to achieve corporate and operational goals. Unforeseen circumstances, such as equipment breakdown, can cause significant deviations from established production forecasts in the short term. Such situations require timely mitigation of the effects of the deviations in order to stay as close as possible to the medium-term plan. This paper presents a deterministic mathematical formulation for short-term scheduling that addresses this challenge. The formulation is able to make up for deficits resulting from operational disruptions by minimizing deviations from medium-term production goals and mediumterm activity start dates. We demonstrate this in a case study by simulating disruptions involving reduction in ore extraction capacity and development advance rate. Aside from the rescheduling capability, this formulation can help mine planners conduct analysis on the impact of various disruptions on the schedule. The formulation can help identify those conditions or disruptions that are more detrimental to the achievement of forecast targets, or otherwise. The formulation also affords mine planners the flexibility to prioritize one production goal over the other, depending on operational conditions at the mine. In the future, we seek to incorporate in the formulation respirable dust, diesel particulate matter, and ventilation constraints in order to create more realistic schedules and enhance miners’ health and safety.

Sep 6, 2023

By Yuting Xue, Lihong Zhou, DAVOOD BAHRAMI

The paper discusses the use of machine learning (ML) to develop a predictive model for determining the location and size of underground mine fires using simulated ventilation data. The study highlights the importance of prompt fire location determination for effective firefighting strategies and reducing injury risks. The ML model was trained with simulated data and achieved an accuracy score of 0.920, which improved to 0.962 by grouping closely connected airways.

Sep 1, 2023

By Yang Song, Li Shaokang, Wang Zhiqiang, Wang Jinjun, Li Guofeng, Su Meixuan

The paper investigates the triboelectric charging properties of calcium carbonate (CaCO3) particles under different conditions using a cyclone tribocharger. It explores the effects of airflow velocity, tribocharge material, and particle moving distance on the charge-to-mass ratio of CaCO3 particles. The study aims to provide guidance for the design of triboelectric separation systems for magnesite.

Sep 1, 2023

By Brendan Demich, Brendan Macdonald, Timothy J. Orr, Jonathan K. Hrica, CASSANDRA L. HOEBBEL, Jason Navoyski, Jennica L. Bellanca

The paper describes the process of collecting detailed accounts of near-miss incidents and turning them into impactful visual stories using graphic simulation. The aim is to create high-quality training materials to reduce haul-truck-related accidents and fatalities in mining.

Sep 1, 2023

By Dion Weatherley, Raúl Castro, RENE GOMEZ, Pablo Cid

The study examines the impact of different drawbell geometries on hang-ups during ore extraction in underground mining. Using the discrete element method, the research compares three drawbell geometries: drawer, circular, and El Teniente. The findings indicate that the El Teniente geometry experienced more hang-up events, while the circular geometry had fewer.

Sep 1, 2023

By Amy McBrayer, Andrea Brickey

The paper discusses the mine planning methods used in the Powder River Basin (PRB), highlighting the differences from open-pit metal mines. It emphasizes the need for agile mine operations due to changing commodity markets and the challenges faced by the coal industry, such as volatile prices and reduced engineering manpower. The study includes a survey of mine engineers to document current practices and identify potential improvements.

Sep 1, 2023

By Yongheng Jiang, Wentong Wang, Haoming Luo, Chuanju Liu, Qianjun Zhang, Sha Guo

Rock burst has always been a major problem in deep underground engineering with high stress, and rock burst strength evaluation has become an important research topic. To effectively predict the rock burst hazard in underground rock mass engineering, a cloud model (CM) rock burst intensity evaluation method based on the CRITIC method and order relation analysis method (G1) was proposed in this paper. First, a rock’s uniaxial compressive strength σc, tangential stress σθ, uniaxial tensile strength σt, ratio of uniaxial compressive strength to tensile strength σc/σt (brittleness coefficient), ratio of tangential stress to uniaxial compressive strength σθ/σc (stress coefficient), elastic deformation energy index Wet, and depth of cover H were selected as evaluation indices of rock burst intensity. Ninety-five groups of rock burst measured data at home and abroad were selected, and the objective weight and subjective weight of each index were calculated by using the CRITIC method and G1 method, respectively. The comprehensive weight was determined according to the combined weighting method of game theory, and the sensitivity of each evaluation index was analyzed. By utilizing a forward cloud generator, the membership degrees of different rock burst grades were calculated, and then the rock burst intensity grades of the samples were evaluated and compared with the evaluation results of the CRITIC-CM method and G1-CM method and the actual grades. Finally, the rock burst classification ability of the model was analyzed. To better verify the accuracy and reliability of this model, the rock burst case of the W39 line in the Chengchao Iron Mine was analyzed by using this model. The research results show that the rock burst evaluation results based on CRITIC-G1-CM are basically consistent with the actual rock burst grade, and the rock burst intensity grade evaluation model has good practicability and reliability.

Sep 1, 2023

By Amar Prakash, V. M. S. R. Murthy, Ajoy Kumar Singh, Avinash Paul, John Loui Porathur

Weak, layered, and fragile rock mass formation, if not supported properly, is subject to roof failure thereby affecting safety and productivity in underground coal mines. Though Central Mining Research Institute – Indian School of Mines Rock Mass Rating (CMRI-ISM RMR) based well-defined support design guidelines are established, still occurrence of roof failures in underground coal mines is a real matter of concern for mining engineers and researchers. Numerical modeling techniques are successfully used by several researchers by simulating rock mass condition for stability assessment of mine openings. The analysis becomes crucial in case of weak and fragile rock formation. The present research envelops the determination of 31 cases of rock load by CMRI-ISM RMR under different geo-mining conditions followed by the development of modified Rock Mass Rating (RMR), i.e., rock mass rating dynamic ( RMRdyn) by incorporating P-wave velocity as a new parameter. The rock load determined using RMRdyn and numerical models was correlated and found in close agreement. In addition, the deviation in rock load determined by all the three approaches, i.e., CMRI-ISM RMR, numerical modeling, and RMRdyn, was compared with the actual field data. The percentage deviation obtained in RMRdyn and numerical modeling is less compared to CMRI-ISM RMR.

Sep 1, 2023

By Andres Salazar, Carlos Delgado Robles, Richard K. Preece

Minera Escondida, located in the Antofagasta Region of northern Chile, is the world’s largest producer of copper as concentrate and SX-EW cathode. Projects to increase mill and flotation capacity began soon after start-up in 1990 to off-set decreasing copper grade. Subsequent projects had the objective to increase copper production through heap leach and SX-EW production and a second concentrator. The mine footprint concurrently expanded peripherally from the central high-grade enriched core of Escondida and then into the separate Escondida Norte deposit. Early in the mine life, metallurgical forecasting that was based on metallurgical composites of presumed ore types and historical plant production was adequate for operational planning, plant design, and project valuation. This approach became unreliable, however, as the increase in mine footprint was accompanied by an increased variability in key ore body characteristics such as hardness, copper-iron sulfide mineralogy, and sulfide liberation textures that were not reflected in the geological criteria used for ore type designation. The project study teams turned to small-scale physical and chemical tests to better characterize the spatial distribution of these key ore characteristics. Programs associated with capital growth projects were carried out through the mid-2000s to sample and construct a three-dimensional spatial model of hardness, copper sulfide mineralogy, and rougher flotation recovery from tests on 15-m drill core composites. Mathematical models were developed from the respective ore characteristic data to calculate throughput, flotation recovery, and concentrate copper grade. The spatial models were integrated with the concurrent mineral resource block model using standard geostatistical methods and the metallurgical process calculations were applied to the block estimates. The geometallurgical model, originally constructed on a project basis, is now managed by a dedicated team reporting to the manager of Geosciences & Resources Definition in the Escondida organization. The model has evolved with increasing complexity of the mine and processing plants and currently includes seven forecasting areas for three concentrators and four areas for two leach plants. Metallurgical forecasts of mill throughput, final flotation recovery, and concentrate copper grade still have the highest relevance to mine planning and optimization. Sample preparation, test procedures, and quality control methods are summarized for the ore characteristic tests important to these forecasts. Analytical results from recent annual sampling programs are presented for the hardness indices, rougher flotation recovery, and partial extraction copper analyses. Geostatistical analysis of the SAG mill hardness index is summarized as an example for the spatial estimation process. Mathematic models are presented for throughput (based on hardness indices), final flotation recovery (based on rougher recovery and plant calibration studies), and concentrate copper grade (based on sulfide mineralogy and empirical pyrite rejection studies).

Aug 31, 2023

By Jiayao Wang, Jian Chen, Leiting Shen, Hongyang Wang, Zeping Cao, Ziru Wang

Large amounts of by-products, including coal gangue, circulating fluidized bed (CFB) slag, and coal fly ash (CFA), are generated during mining, washing and burning of high alumina coal (HAC). The stockpile of these by-products of HAC will cause the pollution of air, water and soil. Meanwhile, the unbalanced regional distribution of bauxite and alumina production has been a critical challenge for the sustainability of world alumina refineries, especially for China. Therefore, economical extraction of alumina from by-products of HAC not only resolves the environmental pollution problems caused by stockpile but also decreases the dependence degree of import bauxite in China. Considering the fact that powder alumina can only be produced by alkali methods, the alkali extraction of alumina from by-products of HAC is summarized in this paper. Based on the formation process of by-products of HAC, the principle, research status, and problem of existing methods were analyzed and discussed, as well as the forecast of research trend. The analyses indicate that simplex recovery of alumina from byproducts of HAC generates large amounts of leaching residue (red mud) with difficulty for disposing, while comprehensive extraction of alumina and silica meets the research trend for clean and high-value utilization.

Aug 25, 2023

By Shilpa S. Haramkar, Bhaskar N. Thorat, Sachin V. Jadhav

The filtration of phosphogypsum is a ubiquitous step in the manufacturing of phosphoric acid. It is highly desired to achieve the minimum concentration of P2O5 in the phosphogypsum cake. This work is aimed at addressing P2O5 losses in phosphogypsum (PG) along with the recovery of P2O5 from slurry by using constant pressure filtration at 1 bar pressure. Laboratory tests were carried out in the phosphoric acid production plant situated in Paradeep, India. Recycled water source near acid plant (i.e., pond water) was used for cake washing to reduce the burden on fresh water stock. The findings indicated that pressure filtration offered the most effective P2O5 recovery (nearly 93.48%). The moisture was reduced to 14–16%, and watersoluble P2O5 in the gypsum cake was found to be negligible. A narrow particle size distribution was observed, which offered effective filterability with average cake resistance of 2.38 × 1010 m/kg, indicating a moderately separable cake.

Aug 23, 2023

By Ping Li, Baokun Zhou, Lei Yang, Changfeng Li, Erhui Zhang

A uniaxial compression experiment was conducted on standard sandstone specimens to explore the sandstone’s mutation characteristics and early warning interval during deformation and failure. Meanwhile, acoustic emission (AE) and digital image correlation (DIC) techniques were introduced to realize multi-field analysis and cross-validation from the stress field, deformation field, and elastic wave field. The mutation characteristics of rock failure were investigated based on the cusp catastrophe theory. The results showed four mutation points in the destabilization process of sandstone: compressional density, yielding, destruction, and residual retention. As an internal activity characteristic of the sandstone failure process, AE directly reflected the dynamic failure process of sandstone. The stress variation was closely related to the evolution of AE characteristics. Sandstone loading failure was a process of the mutations of releasable elastic strain energy, dissipative energy, and releasable elastic strain energy. The mutation of dissipative energy was used as the starting mark of the sandstone warning state, and the second mutation of releasable elastic strain energy was used as the end mark of the sandstone warning state. Besides, the warning range corresponded to the loading time of 178 to 208 s. Finally, the sandstone instability warning interval was re-determined from deformation field evolution characteristics in the sandstone damage process, with a loading time of 180–204 s. This was consistent with the warning interval specified using the cusp catastrophe theory. The combination of AE, DIC, and cusp catastrophe theory revealed the mutation characteristics and early warning interval of sandstone failure. Hence, the results provide significant importance to predicting rock failure.

Aug 23, 2023

By Yuanbo Zhang, Yikang Tu, Zijian Su

Stainless steel swarf is a typical by-product produced from stainless steel processing, which contains a significant level of valuable metals, such as iron, chrome, and nickel. About 1 million tons of stainless steel swarf are generated in China annually, but as mixed with impurities such as Al and Si, stainless steel swarf is not suitable to be directly smelted into stainless steel through refining process. In this study, stainless steel swarf was applied as a strengthening agent during limonitic laterite sintering process. The effects of adding stainless steel swarf on green bed permeability, the yield and quality indexes of the sinters, and the microstructure of the sinters were evaluated systematically. The results showed that the concentrated distribution degree of granules and green bed permeability were improved after adding stainless steel swarf. As the additional amount of stainless steel swarf increased from 0 to 3%, the tumble index and yield of the sinters rose from 50.07 to 54.69% and from 62.04 to 67.35%, respectively. The addition of stainless steel swarf would not change the phase compositions, and had little adverse impact on the metallurgical properties of the limonitic laterite sinters.

Aug 23, 2023

By İbrahim Alp, Emine Yoğurtcuoğlu

This article includes the application of roasting pretreatment by assisting the cyanidation process in the recovery of Au–Ag from Bolkardağ Gossan refractory ore and the characterization studies of the ore, roasting sample, and postcyanidation waste sample. While ore contained many metals along with gold and silver, the environment in which these metals were found was carbonate, sulfated, silicated, and iron oxide-containing deposits that have undergone oxidation. In mineralogical studies, there were Fe oxy/oxyhydroxides, Pb/Zn carbonates, and jarosite-bearing structures such as beudantite, calcite, dolomite, clay, and quartz in the ore. The sample used in the experiments was d80 = 29.6 μm, and the test parameters were determined by examining the TG analysis made to the ore. Preroasting experiments were carried out at a roasting temperature range of 200–800 °C, during a roasting period of 10.0 min. Three different samples were characterized by TG for temperature changes, FT-IR for bond structures and permeability, SEM/EDX for visual inspection of structural changes, and XRD Rietveld’s analysis for mineral % content differences. As a result of the cyanidation experiments carried out after the pretreatment, gold was recovered with 80–90% yields and silver with 40–48% yields, respectively.

Aug 21, 2023

By Rodolfo Lagos, Nicos Pfeiffer, Ben Harding, Clayton V. Deutsch

This paper documents a probabilistic resource modeling workflow for the Round Mountain epithermal gold deposit. The set of categorical and continuous variable simulated realizations provides a model of variability and uncertainty given the current orebody knowledge. These realizations can be used to quantify and understand expected deviations in reconciliation during mining and provide a more accurate measure of expected gold than an estimated model. Gold mineralization at Round Mountain is primarily contained within tertiary tuff units, which dip gently to the west and thicken towards the southern part of the orebody. Domain boundaries are modeled as combinations of lithology and alteration; alteration domains are constrained to be within the tuffs. Large-scale geologic uncertainty is strongly influenced by parameter uncertainty, particularly the input histogram. The spatial bootstrap characterizes prior global uncertainty in the input histogram for categorical and continuous variables. Categorical models are generated using hierarchical truncated pluri-Gaussian (HTPG) simulation. Continuous variable simulation is carried out within the simulated domain boundaries using trend models to account for non-stationarity.

Aug 19, 2023

By Lei Zhang, Zhibo Ma, Pingye Guo

The mechanical behavior of hard roofs often leads to significant deformation and strong mine seismicity in gob-side entry retaining (GER) of the longwall mining, affecting the mining safety and efficiency of coal mines. Field investigations indicate that the failure process of the retained entry by hard roofs can be divided into three stages. In this paper, theoretical analysis and numerical simulations are first used to evaluate the damage process of the retained entry under the hard roof. The results show that the instability of the retained entry in the conventional GER method results from the combined static-dynamic loads (in stage II and stage III) from the hard roof. Then, an innovative method, through fracturing roofs to maintain entry (FRME), based on GER with the hard roof is introduced to solve this problem. In this method, the roof of the mining area is pre-fractured, and when the coal seam is mined, the direct hard roof will collapse along the fractured surface due to its deadweight and mining pressure. The hard roof control effect on the static and dynamic loads, taking into account the innovative method, is verified with simulations and theoretical analyses. The results are compared with those of the conventional GER method. Furthermore, based on the numerical simulation study of dynamic responses, it is found that compared to the GER method, the dynamic response of the retained entry using the innovative method is more moderate and more conducive to the stability of the entry. Field application also validated the effectiveness of the innovative method to control the mechanical behavior of the “massive sudden collapse” of the hard roof.

Aug 14, 2023

By Min Chen, Qinghua Liang, Weishan Yan, Yuan Feng, Shigang Wang

Various coal mine robots (CMRs) and unmanned aerial vehicles (UAVs) are implemented to explore unknown mines for improving the safety and efficiency of mining. It is challenging for CMRs and UAVs to achieve accurate positioning due to the absence of GPS, poor lighting conditions, and similar geometric features in complex mine scenes. LiDAR-based localization and mapping methods are admittedly more accurate than others, while long-time running in large-scale scenarios will also introduce non-negligible cumulative errors. This study presents a semantic-aided LiDAR simultaneous localization and mapping (SLAM) with loop closure, which leverages the uniqueness of mine intersection structure to establish stable semantic loop closure. Specifically, we propose a semantic intersection descriptor of translation and rotation invariance, which encodes 3D point clouds of the same intersection from different positions and viewpoints into a unified image. By using the semantic descriptor, we can construct a constant loop constraint when the same intersection is revisited from different directions to reduce cumulative drift. We provide experimental validation using large datasets collected in two large underground mines, namely, a simulated Edgar mine deployed in ROS Gazebo and a public underground mine dataset provided by ETH. Experimental results show that the proposed method has higher localization accuracy and outperforms the existing LiDAR-based SLAM strategies.

Aug 10, 2023