Underground Copper Mine Research Articles

Geomechanical Analysis of the Main Roof Deformation in Room-and-Pillar Ore Mining Systems in Relation to Real Induced Seismicity

Rockbursts represent one of the most serious and severe natural hazards emerging in underground copper mines within the Legnica–Glogow Copper District (LGCD) in Poland. The contributing factor determining the scale of this event is mining-induced seismicity of the rock strata. Extensive expertise of the copper mining practitioners clearly indicates that high-energy tremors are the consequence of tectonic disturbances or can be attributed to stress/strain behaviour within the burst-prone roof strata. Apparently, seismic activity is a triggering factor; hence, attempts are made by mine operators to mitigate and control that risk. Underlying the effective rockburst control strategy is a reliable seismicity forecast, taking into account the causes of the registered phenomena. The paper summarises the geomechanics analyses aimed to verify the actual seismic and rockburst hazard levels in one of the panels within the copper mine Rudna (LGCD). Two traverses were designated at the face range and comparative analyses were conducted to establish correlations between the locations of epicentres of registered tremors and anomaly zones obtained via analytical modelling of changes in stress/strain behaviours within the rock strata. The main objective of this study was to evaluate the likelihood of activating carbonate/anhydrite layers within the main roof over the excavation being mined, with an aim to verify the potential causes and conditions which might have triggered the registered high-energy events. Special attention is given to two seismic events giving rise to rockbursts in mine workings. Results seem to confirm the adequacy and effectiveness of solutions provided by mechanics of deformable bodies in the context of forecasting the scale and risk of dynamic phenomena and selecting the appropriate mitigation and control measures in copper mines employing the room-and-pillar mining system.

Integrated stochastic optimisation of stope design and long-term production scheduling at an operating underground copper mine

ABSTRACT Conventional underground long-term mine planning is based on a deterministic stepwise framework, which is unable to effectively manage the synergies between the mine planning components or to manage the orebody risk in production schedules and forecasts. The present study enhances prior integrated optimisation approaches to a variation of the sublevel longhole open stoping mining method using backfilling, through a new two-stage stochastic integer programming formulation. A comparison to a sequential stochastic approach shows different stope layouts and extraction sequences for a copper mine with secondary elements being gold and uranium. The integrated approach shows lower horizontal development costs and 6% higher net present value.

Optimizing Retention Bunkers in Copper Mines with Numerical Methods and Gradient Descent

This study examines the optimization of ore receiving bins in underground copper mines, targeting the reduction of rapid wear and tear on bin components. The investigation identifies the primary wear contributors as the force exerted by the accumulated ore and the velocity at which ore particles move. By altering design and operational parameters, the objective is to decrease wear at key points such as transfer areas, thereby improving the efficiency and service life of retention bunkers. A Discrete Element Method (DEM) model of the bin was created and validated against actual mining conditions to study the impact of material flow on wear. The optimization approach used a constrained gradient descent algorithm to minimize factors like particle velocity and pressure force, while maintaining the efficiency of the bin. The findings provide valuable insights for the future design enhancements, potentially improving the operational performance of retention bunkers in the mining industry.

Coprinopsis laanii fruiting on mammal carcasses in an underground mine

ABSTRACT Fungi are important decomposers of organic material, including animal waste. Ammonia and postputrefaction fungi grow in soil enriched in ammonium and nitrogen from carcasses. In 2014, we observed mushrooms fruiting on the flesh of a dead muskrat (Ondatra zibethicus) in an abandoned underground copper mine in southeastern New Brunswick, Canada. We placed an adult beaver (Castor canadensis) carcass near the muskrat to facilitate fungal colonization and fruiting. The beaver carcass was colonized by a variety of molds, especially Acaulium caviariforme. We observed mushrooms of an unidentified copriniid on the flesh 6 years and 9 months after carcass placement. Using morphological and molecular (nuclear internal transcribed spacer [nrITS]) data, we identified the mushrooms as Coprinopsis laanii, a rarely encountered species generally considered lignicolous. We discuss the role of C. laanii, and other postputrefaction fungi, in cave environments.

Integrated Stochastic Underground Mine Planning with Long-Term Stockpiling: Method and Impacts of Using High-Order Sequential Simulations

The integrated optimization of stope design and underground mine production scheduling is an approach that has been shown to effectively leverage the synergies among these two underground mine planning components to generate truly optimal stope layouts and extraction sequences. The existing stochastic integrated methods, however, do not include several elements of a mining complex, such as stockpiles, due to the computational complexity and non-linearity that they might add to the optimization of the long-term mine plan. Additionally, sequential simulation methods that rely on two-point statistics and Gaussian distribution assumptions are commonly used to generate the input realizations of the mineral deposit. These methods, however, are not able to properly characterize complex spatial geometries or the high-grade connectivity of non-Gaussian and non-linear natural phenomena. The present work proposes an extension of previous developments on the integrated stope design and underground mine scheduling optimization through an expanded stochastic integer programming formulation that incorporates long-term stockpiling decisions. An application of the proposed method at an operating underground copper mine compares the cases in which the geological simulated orebody models are based on high-order and Gaussian sequential simulation methods. The extraction sequence and related final stope design are shown to be physically different. It is seen that the optimization process takes advantage of the better representation of high-grade connectivity when high-order sequential simulations are used, by targeting the areas with grades that follow the mill’s blending requirements and by making less use of the stockpiles. Overall, a 4% higher copper metal production and a resultant 6% higher net present value are observed when high-order sequential simulations are used.

The Effect of the Blasthole Diameter on the Detonation Velocity of Bulk Emulsion Explosive in the Conditions of Selected Mining Panel of the Rudna Mine

The blasting technique is currently the basic excavation method in Polish underground copper mines. Applied explosives are usually described by parameters determined on the basis of specific standa ...

The Issue of Determining the Size of Main Excavations Protective Pillars in Deep Underground Copper Mines

The Issue of Determining the Size of Main Excavations Protective Pillars in Deep Underground Copper Mines

Analysis of the Epicenter Location Accuracy for the Local Seismic Network Operated in the Mining Area Towards the Automation of Location Procedures

Legnica-Głogów Copper District (LGCD) is one of the most active seismic regions in Europe. Several thousand seismic events induced by underground copper mines are recorded there each year, with the strongest reaching magnitudes above 4. Seismicity in LGCD is monitored by the LUMINEOS surface seismic network and the mine's underground networks. While the horizontal location uncertainty of the LUMINEOS network is about 300–600 m, the declared epicenter uncertainty of dense mine networks is less than 50 m. It gives us a unique opportunity to test the location accuracy of seismic events recorded by the LUMINEOS network determined with various algorithms and automatic procedure. In our work, we compared the location accuracies of traveltime-based algorithms LocSAT and NonLinLoc as well as waveform-based algorithm BackTrackBB. The iterative, Geiger-type LocSAT algorithm is used in routine daily processing for the LUMINEOS. Its downside is the need to use the S wave onset times, which in the case of the LUMINEOS network are very uncertain. As an alternative, we tested the probabilistic NonLinLoc algorithm, and the waveform based BackTrackBB algorithm. The BackTrackBB algorithm is known to work well with local, high noise seismic networks. We aimed to find out if it could be used routinely with the LUMINEOS data. In addition, we conducted a comparative analysis of the location uncertainty of these algorithms to determine the effectiveness of this parameter in assessing accuracy.

Application of Machine Learning in Forecasting the Impact of Mining Deformation: A Case Study of Underground Copper Mines in Poland

Open access to SAR data from the Sentinel 1 missions allows analyses of long-term ground surface changes. The current data-acquisition frequency of 12 days facilitates the continuous monitoring of phenomena such as volcanic and tectonic activity or mining-related deformations. SAR data are increasingly also used as input data in forecasting phenomena on the basis of machine learning. This article presents the possibility of using selected machine learning algorithms in forecasting the influence of underground mining activity on the ground surface. The study was performed for a mining protective area with a surface of over 500 km2 and located in western Poland. The ground surface displacements were calculated for the period from November 2014 to July 2021, with the use of the Small Baseline Subset (SBAS) method. The forecasts were performed for a total of 22 identified subsidence troughs. Each of the troughs was provided with two profiles, with a total of more than 10,000 identified points. The selected algorithms served to prepare 180-day displacement forecasts. The best results (significantly better than the baseline) were obtained with the ARIMA and Holt models. Linear models also provided better results than the baseline and their performance was very good at up to 2 months forecasting. Tree-based models including their sophisticated ensemble versions: bagging (Random Forest, Extra Trees) and boosting (XGBoost, LightGBM, CatBoost, Gradient Boosting, Hist Gradient Boosting) cannot be used for this type of predictions since Decision Trees are not able to extrapolate and thus are not a valid stand-alone tool for forecasting in this type of problems. A combination of satellite remote sensing data and machine learning facilitated both the simultaneous quasi-permanent monitoring of ground surface displacements and their forecasting in a relatively long time period.

Influence of a Type of Rock Mass on the Stability of Headings in Polish Underground Copper Mines Based on Boundary Element Method

This paper presents the results of a numerical analysis of the impact of rock mass geomechanical parameters on the stability of preparatory headings located within the Legnica-Glogow Copper District. The paper shows the results of numerical calculations prepared for headings driven in two rock mass types with different strength and deformation parameters, which allow illustrating their influence on the safety of mining performed in underground copper ore mines. Numerical modeling was performed using the Examine2D 7.0 software, in the plane strain state. Numerical simulations were performed for an isotropic and for a homogenous medium. The rock medium was described with an elastic model. The parameters of the rock mass for numerical modeling were calculated using the Hoek–Brown classification. The Coulomb–Mohr strength criterion was adopted as a measure for assessing the rock mass effort. Numerical simulations confirmed the dependance between the stability of the analyzed excavations and rock mass geomechanical parameters.

Novel approach for the destress blasting in hard rock underground copper mines

The present study investigates the possibility of developing a novel method for reducing seismicity and rockbursts in deep underground mines based on modifying drilling and blasting patterns. The main goal was to develop and implement firing patterns for multi-face production blasting, which allow increasing the capability of inducing stress relief in the rock mass, manifested in the seismic event. This method may improve stability control in underground workings, and mitigate risks associated with the dynamic effects of rock mass pressure compared with currently used methods. Thus, the seismic energy may be released immediately after blasting in a controlled way. For this purpose, underground tests using modified blasting patterns and precise electronic detonators were carried out. Vibration data recorded from the multi-face blasting in the considered trial panels were assessed in the scope of amplitude distribution. Results of trials have proven that the method is promising and should be further developed to improve the effectiveness of rockburst prevention in deep hard rock mines.

Back-Calculation Method for Estimation of Geomechanical Parameters in Numerical Modeling Based on In-Situ Measurements and Statistical Methods

An important element of numerical modeling for specific mining issues is the selection of model parameters. The incorrect determination of geomechanical parameters can result in significant calculation errors carried throughout the entire problem. This paper presents a method for determining effective geomechanical parameters for technological and residual pillars through the use of numerical modeling, specifically, back-calculation. This is based on the results of numerical simulations, measurement data (e.g., excavation convergence measurements), and statistical methods (a non-linear regression model with “dummy” variables). The result is that appropriate parameters of pillars are set out iteratively so that the displacements of selected points in the numerical model correspond (with some approximation) to the results of mine measurements. The procedure of determining pillar parameters is presented using a case study of one mining field in an underground copper mine, where the deposit is mined using the room and pillar system. Numerical calculations were performed using a Phase2 v. 8.0 program (Rocscience, Toronto, Canada), while statistical calculations used a Statistica computer program. The results of excavation convergence measurements performed in the analyzed mine have been applied. This paper shows that for the presented method, the resulting matching of theoretical values of convergence determined numerically for specified pillar parameters to in-situ results of convergence measurements, is very good (R2 = 0.9896). This work exemplifies a set of the parameters of pillars for an elastic model of rock mass, but this method can also be applied to other models.

Photocatalysis as a method for reduction of ambient NOx in deep underground mines

Air pollution is a critical global issue. The introduction of further restrictions on the emission of harmful gases, such as carbon monoxide and nitrogen oxides (NOx), is posing a growing challenge to various industries. This problem also affects the underground mining sector, where in most cases the entire fleets are powered by diesel engines. It is therefore important, in order to improve the health and safety of underground miners, to develop appropriate methods to reduce exposures of underground miners to harmful gas and particulates. This study investigated the effectiveness of photocatalysis as a method to reduce the concentration of nitrogen oxides in underground copper mines in Poland. The study was based on laboratory measurements using an AVL GAS PEMS gas analyzer. A number of experiments were carried out and the results obtained were subjected to statistical analysis. By fitting appropriate curves – linear and sums of two exponents – an assessment was made of the process carried out comparing the obtained coefficients. This study showed that methodology based on photocatalysis has potential as a strategy for controlling exposures of underground miners to nitrogen dioxide.

Practical Use of Measuring the Deflection of Roof Layers in the Assessment of the Stability of Mining Excavations in the Polish Copper Ore Mine “Polkowice-Sieroszowice”

The paper presents the results of long-term continuous measurements of the deflection of the roof layers in the underground copper mine Polkowice-Sieroszowice, Poland belonging to KGHM Polska Miedź S.A. The measurements were performed with the use of the inclinometric method consisting of continuous registration of changes in the angle of inclination of the roof strata. The measurements were carried out using an inclinometer sensor fixed to the end of a rockbolt in the roof. Measurements presented in the article were made in various regions of the underground mine. The monitoring covered: The exploitation front, machinery chamber and the region of the experimental longwall mining of copper ore. The obtained results proved the usefulness of the developed method in the process of the evaluation of the stability of mining excavations. The sensors were highly sensitive and performed the measurements in a simple way; highly accurate and reliable results were obtained.

Environmental impact improvements due to introducing automation into underground copper mines

A life-cycle assessment (LCA) model was developed to comparatively analyze the use of manual and automated mining equipment in underground copper mine sites. Processes and key variables that were determined to contribute to the environmental impact of operations were identified for six mine sites in a range of geographical locations around the world. Our model successfully calculated carbon dioxide (CO2 eq.) emissions to within 4.9% of the reported annual emissions from the site’s respective companies. The implementation of automation was found to decrease global warming potential by a range of 11.4%–18.0% or 3.9–17.9 kg CO2 eq./t ore. The model was also used to estimate the average reductions across several impact potentials including, acidification (11.9%–17.8%), eutrophication (7.6%–13.7%), and human toxicity (16.0%–20.0%). World-wide the mining industry is moving toward introducing significantly more automation to enhance productivity and safety. This novel work demonstrates an important third dimension that can support this move, reduced environmental impact.

Comparison of H-B and M-C Failure Criterion at Low Confinement and their Impact on the Prediction of Failure Zone Around an Underground Metal Mine

The accurate assessment of rock mass strength and the failure zone near the excavation is important for the rational design of underground structures. To predict the strength and behaviour of the rock mass during excavations, generally, the two most widely used failure criterion namely, linear Mohr-Coulomb and nonlinear Hoek-Brown criteria are used. However, the failure envelopes from the two criteria deviate much, especially, at low confinements. The main motive of this investigation is to understand the stress state around an opening. An underground copper mine of Hindustan Copper Limited (HCL) i.e., Khetri Copper mine in India was considered for this research where the host rock is quartzite schist. The mine is operated at a depth of around 475 m and the ore is extracted using the sublevel stopping method. From the finite element model, it was found that the stress state of the rock mass around the excavation mostly lies in the range between the direct tension and low value of confinements. This is the range where the H-B and M-C failure envelopes differ. The FEM modelling showed a different value of failure zone around the opening, 17.8 m from H-B compared to 8.0 m from the M-C criterion. This research suggested that for the quartzite rock investigated, the support requirements are more if the H-B failure criterion is considered instead of the M-C failure criterion.

Influence of Driving Direction on the Stability of a Group of Headings Located in a Field of High Horizontal Stresses in the Polish Underground Copper Mines

This paper investigates the problem of stability in a group of headings driven in high horizontal stress fields in the copper ore mines of the Legnica-Glogow Copper Belt (LGCB). The headings are protected with the roof bolting system. This problem is of high importance due to special safety regulations which apply in mining workings serving as airways and haulageways. The analysis was performed for a group of four headings driven in the geological and mining conditions of the Polkowice-Sieroszowice mine. The stability of the headings was evaluated with the use of Finite Element Method (FEM). The parameters of the rocks used in the numerical modeling have been determined on the basis of the Hoek–Brown classification, with the use of the RocLab 1.0 software. The parameters of the stress field have been identified on the basis of in situ measurements, which were performed in the Polkowice-Sieroszowice mine in 2012. The measurements were carried out with the use of the overcoring method, which is a stress relief method. A CSIRO HI probe was used as the measuring device. The tests were carried out on three measuring points, on which six successful tests were performed. The measurements confirmed the presence of high horizontal stresses in the rock mass. Numerical modeling was performed using the Phase2 v.8.0 software, in a triaxial stress state and in a plane strain state. The rock mass was described with an elastic-plastic model with softening. Numerical analyses were based on the Mohr–Coulomb failure criterion. It was assumed that the optimal measure of the stability of the group of headings is the range of the formed zone of yielded rock mass in the excavation roof. Numerical simulations have shown that the direction of driving the headings in the field of increased horizontal stresses may be of key importance for the stability of the headings in LGOM mines. The greatest extent of the yielded rock mass zone in the excavation roof occurred when the group of headings was driven in the direction perpendicular to the direction of the maximum horizontal stress component σH. The obtained results served to provide an example of the application of a roof bolting system to protect headings driven in unfavorable conditions in a high horizontal stress field.

Autonomous Loading System for Load-Haul-Dump (LHD) Machines Used in Underground Mining

This paper describes an autonomous loading system for load-haul-dump (LHD) machines used in underground mining. The loading of fragmented rocks from draw points is a complex task due to many factors including: bucket-rock interaction forces that are difficult to model, humidity that increases cohesion forces, and the possible presence of boulders. The proposed system is designed to integrate all the relevant tasks required for ore loading: rock pile identification, LHD positioning in front of the ore pile, charging and excavating into the ore pile, pull back and payload weighing. The system follows the shared autonomy paradigm: given that the loading process may not be completed autonomously in some cases, it takes into account that the machine/agent can detect this situation and ask a human operator for assistance. The most novel component of the proposed autonomous loading system is the excavation algorithm, and the disclosure of the results obtained from its application in a real underground production environment. The excavation method is based on the way that human operators excavate: while excavating, the bucket is tilted intermittently in order to penetrate the material, and the boom of the LHD is lifted on demand to prevent or correct wheel skidding. Wheel skidding is detected with a patented method that uses LIDAR-based odometry and internal measurements of the LHD. While a complete loading system was designed, the validation had to be divided in two stages. One stage included the rock pile identification and positioning, and the other included the charging, excavation, pull back, and weighting processes. The stage concerning the excavation algorithm was validated using full-scale experiments with a real-size LHD in an underground copper mine in the north of Chile, while the stage concerning the rock pile identification was later validated using real data. The tests showed that the excavation algorithm is able to load the material with an average of 90% bucket fill factor using between three and four attempts (professional human operators required between two and three loading attempts in this mine).

Localization of LHD Machines in Underground Conditions Using IMU Sensors and DTW Algorithm

This article presents the concept of using the DTW algorithm to partially solve the problem of locating LHD (load, haul, dump) in an underground mine. The concept assumes the recognition of characteristics—patterns that are hidden in vibrations recorded by vehicles—in segments of the route in the underground excavation, which under appropriate conditions enables the obtainment of information similar to that obtained through the use of RFID gates. With the use of this solution in practice, there are several problems that are addressed in this article. One of the main issues is the different arrangement of the signal fragments resulting from driving along with characteristic parts of the route (bumps, paving damage, lumps of excavated material, etc.) at different driving speeds. This problem was solved by using a combination of the road quality detection algorithm and the DTW algorithm, which estimates the similarity of time series with different lengths. The concept was developed and pre-tested using a test rig and a constructed wheeled robot, and then validated in the conditions of the KGHM underground copper mine in Poland, where the readings from the typical haulage process of an LHD vehicle were analyzed.

Application of Hydraulic Backfill for Rockburst Prevention in the Mining Field with Remnant in the Polish Underground Copper Mines

In the polish underground copper mines owned by KGHM Polska Miedz S.A, various types of room and pillar mining systems are used, mainly with roof deflection, but also with dry and hydraulic backfill. One of the basic problems associated with the exploitation of copper deposits is rockburst hazard. Aa high level of rockburst hazard is caused by mining the ore at great depth in difficult geological and mining conditions, among others, in the vicinity of remnants. The main goal of this study is to investigate how hydraulic backfill improves the geomechanical situation in the mining filed and reduce rockburst risk in the vicinity of remnants. Numerical modeling was conducted for the case study of a mining field where undisturbed ore remnant, 40 m in width, was left behind. To compare the results, simulations were performed for a room and pillar mining system with roof deflection and for a room and pillar mining system with hydraulic backfill. Results of numerical analysis demonstrate that hydraulic backfill can limit rock mass deformation and disintegration in the mining field where remnants have been left. It may also reduce stress concentration inside or in the vicinity of a remnant, increase its stability, as well as prevent and reduce seismic and rockburst hazards. Hydraulic backfill as a local support stabilizes the geomechanical situation in the mining field.