Search Documents

By Wenkai Wan, Charlie C. Li

Strain burst can occur in underground hard rock excavation when the in-situ stresses are high enough. The prerequisite for rock burst is that the total energy released from the rock mass is higher than the energy needed for rock fracture. The total released energy is contributed by both the burst rock and the surrounding rock. The part of the energy from the burst rock is constant per unit volume of the rock, but the part from the surrounding rock is associated with the burst volume. A greater amount of energy will be released from the surrounding rock in the case of a larger volume of burst rock. The postulation for the energy conversion in strain burst is verified through a FLAC3D model in the paper. It is assumed that the strain burst occurs in a circular tunnel that is located in the hydrostatically loaded rock mass in the model. Six pairs of notches with different depths are created in the opposite walls of the tunnel to simulate the bursts. The energies released from the notch rocks and from the surrounding rock are calculated during the burst. The simulation results show that the percentage of the energy released from the surrounding rock in the total released energy increases with the notch depth. The numerical modelling confirms the postulation that the surrounding rock contributes more energy than the burst rock itself in large-scale strain bursts

Apr 26, 2022

By M Ali, E Satriawan, M Napitupulu, A Munandar

The Deep Mill Level Zone (DMLZ) panel cave mine is one of three active caving mines at PT Freeport Indonesia’s Grasberg Underground Mining Complex. The DMLZ mine is located at a depth of approximately 1,500-1800 m below the surface, making it one of the world’s deepest panel cave mining operations. This depth results in both high in-situ and induced stress that acts on the hard and massive rock mass, triggering rockburst conditions. As a result of the bursting conditions, ground support and maintaining the support capacity are an essential part of operating the DMLZ mine in an efficient and safe manner.

Apr 26, 2022

By Patrick Linzer, Lindsay Linzer

By its nature, seismology has been a data intensive field since the first digital computers were introduced in the 1950s. The first empirical models were developed in the same period and are still widely used today, for instance, the Gutenberg-Richter frequency-magnitude relationship (Gutenberg and Richter, 1956). Machine learning (ML) methods make it possible to redefine these empirical models using large datasets and examining considerably more parameters, where their significance to the problem being addressed is methodically and repeatably evaluated using robust statistical methods. In this paper we develop a simple logistic model for the short-term prediction of seismic events above a threshold magnitude from a mine’s event catalogue comprising more than a million events with 14 basic parameters from which a pool of over 830 derived parameters was developed, covering a series of 11 precursory time windows. The logistic model was fitted stepwise using SAS with α = 0.025 for both entry and acceptance. 25 variables from the pool met the acceptance criterion. The findings were surprising in that parameters a rock engineer or seismologist might have expected to be the most significant in the logistic model were not; instead, it was an increase in standard deviation of the immediately-prior seismic magnitudes that carried the most weight, followed by the average maximum source displacement and median apparent volume change. These findings make intuitive sense because the increasing scatter of the magnitude shows an increasing number of large events being recorded that could be indicative of instability in the system that could lead to large events. The model is rudimentary and enormous scope exists for more exploration in terms of data pre-processing, model types, and response variable definitions.

Apr 26, 2022

By Zong-Xian Zhang, Charlie C. Li

It has been observed that the rock debris of a strain burst event are often in the mode of tensile fracture. Seismic waves could create tensile fractures in the material nearby reflecting surfaces, such as at the boundary of the fracture zone and the rock surface of an underground space. The fractured rock could be ejected to the space when the seismic waves are powerful enough. The average ejection velocity of the ejected rock is a measure of the magnitude of a burst event. The ejection velocity is directly associated with the released energy from the seismic waves and with the strain energy stored in the rock prior to the burst. In some cases, the seismic wave induced tensile fracture and the energy conversion occur in successive moments during bursting. In this paper, tensile fracture caused by seismic waves is investigated through an 1D model. The relationship of the average ejection velocity of the ejected rock in a rockburst event is analytically established with the energy of the seismic waves and the strain energy stored in the rock.

Apr 26, 2022

By B. Simser, T. Butler, A. Hall

The longterm future of Glencores’ mining operations in the Sudbury basin area of Ontario, Canada likely depends on mining deposits at depths below two kilometres. To gain a better insight into the rockmass response to mining, a geological scouthole program has been implemented at Glencores’ Nickel Rim Deep and Onaping Depth Project to gain insight into rockmass conditions ahead of the devel-opment face. The scouthole program involves drilling long (up to 600 m) subhorizontal boreholes drilled ahead of the development end and collecting data about the rockmass conditions from various geophysi-cal and mechanical tools (in addition to conventional core logging). Once all data has been collected and the borehole is available, seismic sensors are installed up to 450 m down the hole as part of the micro-seismic monitoring program with the leading sensor always located ahead of the development end. The recorded seismic data is used to provide realtime insight into the rockmass response to mining. On May 20, 2017 a flurry of seismicity with a maximum moment magnitude (Mw) of 0.7 occurred causing damage to the ramp heading. Data from the microseismic monitoring system had indicated that conditions at the development end were changing prior to the rockburst. A back analysis of the scouthole data indicated a clear anomaly in the rockmass at the location of the rockburst. The scouthole data also provided sufficient information to suggest that once this anomalous zone was exited, mining conditions would improve. With the microseismic monitoring system providing realtime information from the rockmass, the area was re-habilitated and ramp development continued successfully. Details of the rehabilitation and adjustments to the ramp development are discussed. Some correlations between the scouthole data and recorded micro-seismicity are also presented.

Apr 26, 2022

By Liliana Rybarska-Rusinek, Aleksandr M. Linkov, Ewa Rejwer-Kosińska

Joined modeling of stressed state and induced seismicity is a means to join geomechanics with seismic observations. When applied to quantifying rockburst hazard, it provides, on one hand, estimations in terms of mechanics, such as stress concentration, and, on the other hand, in statistical terms of the risk analysis developed in mining seismology. The modeling became possible due to progress in computers and numerical methods. However, despite tremendous increase of computational abilities, the joined modeling is still quite time expensive. To make it available for wide use and real time simulations, there is need to reduce the time cost of its key stages: calculation of stresses, caused by mining, and simulation of seismic events, induced by these stresses. The present work aims to decrease the time cost. The stages are time costly because of the same reason, which is a very large number N (up to millions) of elements, representing rock structure, openings and flaws. This leads to immense computational complexity (roughly of order N2) of conventional accounting for pairwise interactions between the elements. The complexity may be notably reduced, when employing the fact that interactions decrease with growing distance between elements. Specifically, the fast multipole methods (FMM) avoid the complexity N2, by employing translations of far-fields from whole collections of sources. This reduces the complexity to the order 22, where p is the number of terms involved in a single translation, is the number of translations, and is the dimension of the quantity translated. Still, for conventional FMM, the complexity remains too high to routinely perform joined modeling of stresses and seismicity.

Apr 26, 2022

By Stephanie Singarayer, Damodara Reddy Chinnasane, Mike Yao, Danielle Pelletier

Mining at Vale’s Creighton Mine has been active for nearly 120 years in Sudbury, Ontario, Canada. Being the deepest mine in the Sudbury region, mining induced seismicity is a common occurrence at Creighton Mine. Over the past five years, 26 large magnitude seismic events ranging from 2.5 Mn to 3.7 Mn (Nuttli) have occurred. Out of these large magnitude events, three 3.0+ Mn events occurred between April 2020 and January 2021. The majority of the 3.0+ Mn events are believed to be the fault-slip type that could be linked to the known major geological structures and/or long-term effect of mining over many years. The recent large magnitude events are believed to be triggered by high induced stresses while mining the first few stopes of a new lower horizon in the top-down centre out mining sequence.

Apr 26, 2022

By Christina Dahnér, Kristina Jonsson, Savka Dineva

The iron-ore underground Kiirunavaara mine (Sweden) has been seismically active since 2008. In 2021 the production is on levels 1079 and 1108 m, at the deepest. The seismic activity has changed in character over time with one step change in 2015 when the number of large (ML>2) seismic events increased, and another one when the mine experienced a seismic event with Mw 4.2 in 2020. In 2017, the mine evaluated Hydraulic Fracturing (HF) as a possible mitigation technique for larger seismic events and a one-hole trial was performed. The trial proved promising, and it was decided to go forward with a large-scale trial covering an entire production area in the mine. The large-scale trial consists of 25 fracturing holes, 50 – 100 m long, running along 500 m out of the ~4500 m long orebody and spanning between production levels 1108 m and 1252 m. The fracturing process continued for a year, resulting in ~38,000 seismic events with −3.3<ML<1.4. The breakdown pressure varied between fractures in the same hole as well as between holes in the same area. Average breakdown pressure per hole was between 12 MPa and 50 MPa. The seismic response to HF in 2021 varied greatly from hole to hole and between fractures in the same hole. The variation in response could be related to the geological conditions and the infrastructure, such as large caving ore passes. In this paper some details about the experiments and broad preliminary analysis of the effect of the HF by using seismic source parameters of the induced seismic events are presented: 1) Distribution of seismic activity, 2) the ‘equivalent’ magnitudes for each hole calculated from the cumulative seismic energy and seismic potency released during the fracturing, compared with the local geological conditions; 3) the predominant failure mechanism of the seismic events caused by HF was defined. A difference in the behaviour between larger and smaller events was noticed, as ML>0 events seem to be compression type failures and ML<0 events seem to be related to slip on pre-existing structures. No seismic events related to creation of a new fracture have been identified yet.

Apr 26, 2022

By Ehsan Moosavi, Kamyar Tolouei, Mehran Gholinejad

The literature defines rockbursts in slightly different ways based on their generation mechanisms, field phenomena, and degree of damage, but in a way that is almost consistent. Generally, a rockburst refers to seismic events that mainly happen in underground excavations, for purposes such as underground mining, road and railway tunnels, nuclear power, etc., and that are followed by the collapse and outburst of rocks. Long-term rockburst prediction can be abstracted as a classification problem mathematically, in which we can introduce deep learning to solve it. Further, using deep learning in rockburst prediction can overcome drawbacks such as subjectivity and inconsistency brought from previous traditional approaches. A rockburst data-set, was employed to evaluate the current method for predicting rockburst grade, and the good results of overall success rate were obtained. The results indicated that the proposed algorithms can speed up parameter optimization search, the proposed method is robust model and might hold a high potential to become a useful tool in rockburst prediction research.

Apr 26, 2022

By E Marambio, R Brändle, G Fischer, G von Rickenbach, M Fuentealba, L Burgos, R Luis Fonseca

In underground excavations, under high stress conditions, ground support systems must provide safe and effective response to dynamic events. These systems must be capable to resist dynamic impacts and yielding during the loading process. In this context dynamic testing of the reinforcement and retaining elements that compose the ground support system are required to study and improve the behaviour of these elements under dynamic load events. During the last years, Geobrugg has been working on the improvement of retaining products by testing them in a large-scale impact test facility located at Walenstadt, Switzerland. The test facility is composed of a double level platform with a square-shaped pyramidal trunk geometry, in the upper level housing a loading mass that drop from a height up to 5m. The loading mass is guided by one central steel pipe, and the impact occurs in the sample to be tested located at the lower level in a slab with an area of 3.6m × 3.6m where the ground support system is installed. During the last years, this innovative facility has been used to test several configurations of ground support systems. The results of these tests have enabled to improve the understating of the behaviour of ground support systems under dynamic loads. In this manuscript, the back analysis through numerical methods (central difference and Newmark’s method) of a ground support system tested with a double dynamic impact in 2019 is presented.

Apr 26, 2022

By Jake Beale, Joe Wickline, M. Van Dyke, Ted M. Klemetti

A longwall coal mine in southwestern Virginia mining the Pocahontas #3 Coal Seam experienced a history of mining-related seismic events. The events were never large enough to cause concern until three major events occurred recording 3.4, 4.3, and 3.4 ML in 2005, 2006, and 2007 respectively. Two of the seismic events generated an air blast which damaged stoppings which fueled a mine fire in 2005 and 2007. The mine responded to the events by narrowing the longwall widths, changing ventilation controls, installing a 7-sensor seismic system around the active panel, and creating a mining hazard map to avoid similar geology in the future. All of the actions taken have prevented seismic activity from impacting mining operations, but it did not prevent other large seismic events from occurring (Van Dyke et al., 2018). In 2016, another large seismic event, measured at 3.7 ML, occurred in the gob that was reported by the local media and recorded by the United States Geological Survey (USGS). While the event did not directly impact mining operations, mine management wanted to know where and under what conditions could large seismic events occur in the future. The National Institute for Occupational Safety and Health (NIOSH) became involved to help understand what conditions precipitated these events so that future mine design could more accurately assess seismic hazards and reduce the risk to mine worker safety. Researchers at NIOSH and mine management created a seismic forecasting map that combined three factors: (1) overburden thickness, (2) sandstone thickness, and (3) the expected caving height of the gob. The criteria were developed based on the geologic data the mine collected from core drilling, e-log analysis, and fiberscope observations. The forecast map was able to provide an accuracy of 72%–85% of events 1.5 ML or greater and 100% accuracy of events that were 3.0 ML or greater based on historical data. The seismic forecasting map was applied to the entire future mining area to show areas where seismic potential could be the greatest based on current mine design planning, ultimately improving the ability of the mine operator to assess the seismic hazards surrounding active and future mining areas (Van Dyke et al. 2018). Following this research, the seismic forecasting map was implemented into the mine’s hazard recognition program and is continuously updated and evaluated prior to mining each longwall panel. This paper will cover in more detail how the seismic and lithologic data was processed and how the mine implemented a new seismic mitigation program based on the forecasting results.

Apr 26, 2022

By M. Sears, W. H. Su, J. Wickline, M. A. Van Dyke, G. Esterhuizen, Z. Khademian

It is commonly recognized that seismicity can be induced by rock failures or slip along geologic structures around mining excavations. The unpredictable nature of these events and their proximity to the working area may raise safety concerns in underground mining operations. Modeling can be a useful tool for estimating seismic potential for a given mine layout, but the challenge is to detect a seismic event and estimate its intensity in the model. One well-documented approach is to monitor the variation in the energy content of the model and record potentially damaging seismic events when there is a significant release of energy. Researchers from the National Institute for Occupational Safety and Health (NIOSH) employed this approach using 3DEC modeling software to simulate the anomalous seismic events that occurred in a deep longwall mine in Virginia in 2005. A four-panel model of the mine was constructed utilizing a 2D cross section at the location of the mining face at the time of the event. Seismic potentials associated with the development of each panel are evaluated by recording the seismic energy released during the simulation. The results demonstrate the occurrence of two relatively significant events during the mining of the last panel, which are consistent with the historical observations at the mine.

Apr 26, 2022

By Kai Koosmen, Mehdi Serati, Tengzhou Wang

Spalling and rockburst are unwanted brittle failures at the boundaries of underground excavations that are commonly observed in hard rock and high-stress environments at significant depths, although, spalling has also been observed in comparatively shallow soft rocks. The crack initiation (CI) stress threshold (representing the stress level at which new micro cracks begin to form in an intact rock) is widely used to establish the in-situ spalling strength, and is generally regarded as the lower limit of the long-term rock strength. This study analyses the evolution of CI and other mechanical properties for a high-porosity sandstone with an average UCS of 38 MPa under unconfined cyclic compressive loading. The results indicate that rock properties are variably affected by cyclic loading with; no significant effect on Youngs Modulus; a moderate increase in Poisson’s ratio and fracture angle; a moderate decrease in UCS; and a major reduction in the CI stress threshold. All of these changes were generally amplified by cyclic loading at higher stresses. From the results, it is concluded that rockburst and spalling limit in underground structures may reduce over time if the rock is subject to repetitive cyclic loading, and that the reductions will be greater for higher cyclic stresses. This may be a significant factor in the design and operation of some deep underground structures.

Apr 26, 2022

By P. K. Kaiser, M. Cai

When designing an efficient and effective support system for burst-prone ground, it necessary to select the most appropriate, i.e., the ‘right’ support components, including areal and reinforcement components. For this paper, it is assumed that the areal support is not the weakest link and the support performance is controlled by the rock reinforcement and holding components (bolts). We discuss the loading modes of rockbolts under static and dynamic loadings and analyze published data from laboratory and field tests. The mechanisms of energy absorption of dynamic rockbolts are presented and rockbolts are categorized based on factors such as energy absorption capacity, energy absorption rate, direct and indirect energy absorption capacities, rock mass bulking or displacement control, and shearing resistance capacity. Guidelines for selecting the ‘right’ rockbolts to build an effective rock support system for burst-prone ground are provided.

Apr 26, 2022

By Ping Zhang, Shahram Mozaffari, Andreas Ylmefors

The risk associated with mining induced seismicity and the related rock bursts has become one of the major threats to the safety and sustainability of mining at the Kiirunavaara mine. The causes of mining induced seismicity at the Kiirunavaara mine are very complex and have affected our understanding on the source mechanisms of seismicity and further the mitigation of the seismic risk. Detailed analysis of historical mining induced seismicity at the Kiirunavaara mine was conducted. The focus of this analysis was placed on the classification of large seismic events and the evaluation of the influential factors on the classified seismic events. 145 events with local magnitude larger than 1.5 that occurred between January 1, 2011 and March 31, 2021 were analyzed. For this analysis, seismic events located in the hanging wall were not included, as these have not been observed to cause large damages to the mine itself. It was found that the magnitude of the events had good correlation with the stress change in relation to the production location. The seismic events were classified as mining-front associated, ore-pass associated, mining-front and ore-pass associated and others, according to their hypocenter location. Most of the mining-front associated events, concentrated around the production levels and slightly above and below them, generally correlated with shear cracks, and most of them occurred when the ore flow between two production levels started. The numerical modelling showed that slip of the geological discontinuity, due to unclamping between hanging wall and footwall, could be the major causative mechanisms, which further supports this finding. The ore-pass associated events are located close to the ore passes. Large seismic events located between ore passes have the source mechanism of both fault-slip and non-shear, but events located at the outer side of ore passes were all fault-slip type. By checking with the ore pass scanning data, it was concluded that the most of the events located between ore passes could be pillar bursts. An elastic core concept has been then used to define a critical pillar width. The paper is not intended to be prescriptive but rather to initiate discussion of the possible classification of seismic events at the Kiirunavaara mine and to present new perspectives of the causes of the seismicity for improving the present seismic risk management.

Apr 26, 2022

By P K. Kaiser, D A. Malovichko

Highly stressed excavations in hard brittle rock may fail by self-initiated, triggered or dynamically loaded strainbursts (Moss and Kaiser 2021). The low-frequency part of seismic radiation associated with strainbursts is described by a crush-type seismic source mechanism with a significant negative (implosive) isotropic component and a ‘pancake-shape’ CLVD deviatoric component. The seismic source is co-located with the excavation damage site. For such strainbursts, the ejection velocity is primarily related to the depth of strainbursting, the rock mass bulking characteristics, the duration of the bulking process, and the effectiveness of the support system. For dynamically loaded strainbursts, the displacement and energy demand may be augmented by the impact of the stress wave from a distant seismic event. This article describes the failure process and factors contributing to rock deformation and energy demands for support design. Case examples are used to demonstrate the dependence of rock displacements, rock ejection and energy demand on the factors described above. A seismic source model is proposed to assess the depth of strainbursting and the duration of the related bulking process by using seismic monitoring data. The implications for support design are briefly discussed.

Apr 26, 2022

By Peter Zhang, Mark Van Dyke, Ted Klemetti, Morgan Sears

Room-and-pillar mining is a mining method used in the United States to extract coal from underground, and pillar recovery is widely practiced for maximizing coal recovery and mining economics. Under certain conditions, room-and-pillar mining can create high stresses in coal pillars resulting in a coal outburst. These outbursts represent one of the more severe safety hazards in underground coal mines, potentially causing injuries and fatalities. Over the past ten years, coal outbursts have resulted in three fatalities in U.S. room-and-pillar mines. Though coal outbursts generally occur under deep cover, they are more likely to occur during pillar recovery and under multiple-seam mining conditions. Coal outbursts are generally associated with strong roof and floor and heavy loading of coal pillars, but it is difficult to predict whether, where, and when they will occur. In this respect, risk management provides a proactive tool to prevent coal outbursts in room-and-pillar mining by identifying mining areas with high risk of outburst and mitigating the risk accordingly.

Apr 26, 2022

By B. Ollila, L. Brown

As mining fronts converge in sill pillar extraction, the mine seismic hazard can be significantly altered and large seismic events have been known to occur. This study used common seismological tools to analyze the extraction of highly stressed sills from two deep Canadian mines with the objective of characterizing the rock mass failure process and relation to large seismic events. The approach treated a sill pillar as a mine-scale failure process, meaning the temporal and spatial analysis scales needed to exist on a similar scale to a significant part of the mine. The temporal analysis scale overlapped with the entire mining sequence, included several years of data to establish long term change in the characteristics of rock mass failure. The spatial analysis scale included the rock mass surrounding a large part of the convergent mining sequence. To accommodate the large spatial and temporal scales, the highest magnitude 2% of events in the dataset were analyzed. Similar results were shown for sill pillars in both mines. Changes in event location, time, magnitude, energy index and apparent volume were similar to expected behavior of rock failure. Indications of strain hardening and strain softening were observed in the sills. The hardening period was characterized by a high rate of small events, high energy index and low apparent volume, giving an indication of high stress and low rock yield. The softening phase was characterized by a decrease in the rate of small events concurrent with an increase in the rate of large events, low energy index and high apparent volume, indicative of release of stress and increase in deformation. Large magnitude seismicity was strongly related to the softening phase, with several Nuttli magnitude 3.0+ events occurring shortly after the transition from hardening to softening failure regimes. Sill failure was indicated by consistently low EI and low levels of seismicity and seismic hazard in the sill pillar.

Apr 26, 2022

By Rodrigo Estay Huidobro, Andrés Martinoli, Claudia Pavez-Orrego

The effects of induced seismicity caused by mining activities are one of the main problems that the mining industry faces. Although it has been documented by countries with a great mining tradition (United States, Canada, Australia and Chile) for over a century, the problem has been intensified throughout the years due to the greater depths that mineral extraction reaches, where the stress levels and environmental conditions increase the possibility of a rockburst occurrence. Seismic events inside the mine cause important losses of human lives and infrastructure damage. To deal with this problem it is necessary to implement several kinds of safety protocols. For example, seismic indicators allow to identify, through the analysis of a set of seismic parameters -magnitude, energy, time, decay rate, among others-, periods of time when a seismic event of considerable magnitude is most likely to occur. The objective of this work is to use the frictional energy generated by a microearthquake inside the mine as a seismic indicator, linking it with the occurrence of large magnitude seismic events which are outer the background activity. Thus, we seek to determine the set of seismic parameters which allow the best performance of the frictional energy/heat indicator. In the present study, we considered magnitude, frictional energy, and time as seismic parameters to study the efficiency of the proposed indicator. The optimal set of parameters was determined using two different methodologies: accumulated thermal energy and alert bands. The Piece Skill Score and the Accuracy indexes were used to evaluate each set of parameters. Both approaches show a high rate of true positive and true negative hits. Furthermore, it was concluded that both methods show complementary results in terms of the detection of critical events (Mw > 3.0).

Apr 26, 2022

By E Mutaz, M Bahaaddini, D J. Williams, M Serati

The contribution of the mining industry to the global economy has reached (and gone beyond) $1,800 billion in recent years that necessitates (more than before) the need to ensure workability, viability, health, and safety of mine sites and workers to meet the ever-increasing global demand for mineral resources. Amongst the many challenges facing the excavation industries, rockburst and spalling failures in deep mines in hard and brittle rock formations are yet happening that are still deemed as unsolved problems. During unconfined compressive testing of brittle rocks, microcracks initiate at approximately 0.3–0.5 of the peak strength (which is commonly known as the Crack Initiation stress level) mainly in the direction of the major principal stress, then grow unstably and coalesce to form macro fractures at approximately 0.7–0.8 of the peak load at a so-called Crack Damage stress threshold. From a stress-strain perspective, the crack damage stress threshold () is associated with the reversal point of the volumetric strain, which is also proved to be strongly affected by the intermediate principal stress. This paper extends the previous studies on the impact of true triaxial loading conditions on the CD stress level in high porosity rocks, and introduce a new Crack Mode-Changing Stress (CMCS) defined as the minimum principal stress required to change rock fracturing from splitting to a sliding failure mode. A large number of experiments on sandstone specimens were carried out at different true triaxial loading conditions and the effect of principal stresses on the CI and CD as well as CMCS were studied. The results show that not only the intermediate principal stress, but also the minor principal stress can significantly affect the CD and CMCS thresholds.

Apr 26, 2022