Abstract

Rock burst is a violent release of energy stored in the rock mass, which has been the focus of metal mine engineers. It gradually becomes a serious threat to both the mine’s workers and equipment. Rock burst risk assessment is a key issue in the prevention of rock burst. Based on the microseismic method for rock burst risk assessment that is used in tunnels, this paper presents a method of rock burst risk assessment for metal mines combining a simplified probability function and principal component analysis (PCA). Specifically, the size and details of key monitoring area and risk assessment units are determined according to the distribution of the fracture and driving stress during deep mining process. Based on the correlation analysis, six microseismic parameters are selected to be used for rock burst risk assessment. Simple linear functions are temporarily used to construct probability functions of each parameter before enough rock burst cases are collected. The method has been successfully used in the deep mining section of the Ashele Copper Mine (approximately 1000 m deep). Additionally, among various microseismic indicators, the energy release rate and cumulative energy are the greatest contributors to rock burst of different intensities. Therefore, the proposed method in this paper provides an important reference for the prevention and control of rock burst in metal mines.

1. Introduction

High stress usually occurs in deep underground engineering. When the actual stress exceeds the rock mass strength, some brittle failure may occur, accompanied with the release of a large quantity of energy [13]. This phenomenon is named rock burst. Rock burst have been observed in mines [4, 5], tunnels [6], hydraulic engineering projects [7, 8], and other forms of underground engineering projects and has become a worldwide problem. The number of rock burst cases is reported in the tunnels of Norway [9], in hydraulic engineering projects in China [10], and in the deep mines of China [11], South Africa [12], and Canada [13]. A rock burst may cause equipment damage, construction scheme delay, and even the loss of lives because of its violence [1416]. Professor Fairhurst indicates that understanding and controlling rock burst is one of the greatest challenges of deep mining engineering safety [17]. Suorineni et al. [18] use a cancer metaphor for contemporary geological mechanics in deep mining to describe rock burst. Obviously, rock burst has become a common problem in deep mining.

Some rock burst risk assessment and warning methods have been reported, which include microseismic monitoring [19], infrared monitoring [20], seismic velocity tomography [21], and computed tomography technology [22], and a study was conducted on rock burst tendency [23]. As a result, many great theories and types of risk assessment indicators have been proposed by researchers based on energy, fracture, damage, and nonlinear theory [2427]. In general, rock burst risk assessment methods can be divided into two types: one is performed before mining, and the other is performed during the mining process [28]. The former type pays attention to mining geology, ground stress, and the mechanical properties of the rock mass and provides a theoretical basis for the optimization of mining parameters, a supporting plan, and the prevention and control of rock burst [29, 30]. The latter assesses rock burst risk dynamically according to detailed conditions and information gained from monitoring the mining process and then takes timely rock burst control measures. To some degree, the latter rock burst assessment method is more accurate than the former.

Over the past few decades, microseismic monitoring technology has gradually become the primary safety monitoring and risk assessment tool used for underground engineering projects in the world [3136]. Many scholars presented risk assessment methods and rock burst indicators after studying the law of microseismicity in the process of construction [37]. For example, at the Mt Charlotte Mine, Poplawski [38] proposed departure indexing for rock burst hazard assessment to overcome the disadvantages of traditional method with only using one or two seismic parameters. Alcott et al. [39] discussed a method of using energy, apparent stress, and seismic moment criteria to identify the relevant events for microseismic hazard assessment. Tang and Xia [40] defined the rate of change of the ratio of the seismic stiffness in the nucleation region to that in surrounding region as the indicator of the rock burst prediction. Tang et al. [41] integrated microseismic monitoring and numerical simulation to improve the accuracy of rock burst prediction at Jinping II hydropower station. Zhou et al. [42] adopted heuristic algorithms like genetic algorithm, particle swarm optimization, and support vector machine to establish hybrid models to predict rock burst. Feng et al. [28] suggested a method for predicting various rock burst intensities and successfully targeted measures to reduce rock burst risk at the Jinping II hydropower project. Cao et al. [43] provided a reference for assessing rock burst risk and determining potential rock burst risk areas in the Xuzhuang Coal Mine. Feng et al. [6] established a warning formula to provide some advance warning of rock burst based on the in situ microseismicity and rock burst case studies in tunnels. Xu et al. [35] found the concentration of microseismic events can be used to predict the strainburst location in deep tunnels after analyzing the spatiotemporal evolution of microseismic activity. Dou et al. [44] introduced the R-value scoring method to determine the weights of bursting strain energy, time-space-magnitude independent information, and time-space-magnitude compound information for rock burst warning in Hujiahe Coal Mine. Ma et al. [45] concluded that the geological structure plays a crucial role in the rock burst prediction of deep tunnels at Jinping II Hydropower Station. Additionally, some mathematical modeling approaches and some learning algorithms are used to predict rock burst [42, 46]. In view of the above work, risk assessment and forecasting studies for rock burst have gained some meaningful results, especially the rock burst risk assessment studies for tunnel and hydraulic engineering.

Simple project layout and long liners are the main features of traffic engineering and hydraulic engineering. Considering the production capacity requirement, multiple stopes at different levels need to be mined at the same time, which leads to the coexistence of numerous goafs in a metal mine. In such conditions, it is difficult for engineers to collect useful information during the rock burst development process. To some degree, it also increases the difficulty of providing early rock burst warning. To date, the rock burst risk assessment method used in metal mines is far from perfect. Considering the differences between the metal mines and the tunnels, such as nonlinearity, simultaneous mining in multimiddle section and multiple stopes and complex layout of shaft and roadway, a specific rock burst risk assessment method is attempted in a metal mine. The site chosen for this study is a deep mining area in the Ashele Copper Mine, Xinjiang Province, China, where the rock bursts are considerable threats to the workers and equipment. The spatiotemporal evolution law of microseismicity is summarized based on the microseismic monitoring system. Then, a model is established to provide the spatial-temporal risk assessment of rock burst. It may provide a reference for determining potential rock burst areas.

2. Rock Burst Risk Assessment Method

Microfractures generate in some areas where the stress is concentrated during the mining process [47]. Before the occurrence of a rock burst, microseismicity may increase accompanied with the development of microfracture. A microseismic monitoring system collects seismic waves released by the microfracture using microseismic sensors which are arranged around the target area. Then, the time, location, and parameters of the seismic source can be obtained based on the waves [48]. In general, there are some precursory characteristics of most rock burst in a metal mine. The microseismic method is used to identify the precursory characteristics to provide rock burst risk assessment and management. According to the literature, few rock bursts have no precursory characteristics [6]. Subsequently, all the rock bursts involved in this paper have microseismic precursory characteristics.

2.1. Determination of Risk Assessment Area and Unit

Multistope mining at different levels and a complex distribution of long drifts are the main features of a metal mine, making it difficult for rock burst risk assessment and management to cover all areas. The deep stope mining, especially with large scale blasts, would inevitably cause stress concentration in some areas [49]. Combined with strong blasting disturbance, the microseismicity and rock burst risk would be more serious. Considering this, the potential rock burst areas should focus on the area where the distribution of microseismic events is concentrated. However, there is no need for the miners and equipment to enter all the drifts. For example, people only worked in the drift without entering the goaf when a sublevel and long-hole open stoping mining method was adopted. Therefore, the surrounding rock stability of the drift, used as a walkway for miners and equipment, is the main focus of the rock burst assessment and management. After collecting a large number of rock burst cases and matched microseismicity, we use statistical analysis to fix the dimensions of risk assessment unit.

2.2. Selection of the Microseismic Indicator

According to the literature review, a number of the microseismic parameters, such as cumulative microseismic events, event rate, cumulative energy, energy rate, cumulative apparent volume, apparent volume rate, stress drop, energy index, b value, Schmidt number, and some other parameters, can reflect the microfracture activity and the pregnant process of the rock burst. These parameters can be used as risk assessment indicators individually and can also be used in combination to avoid the randomness of a single parameter [38]. Generally, the selection of the microseismic parameters for rock burst risk assessment mainly relies on empirical analogy. In this paper, the rock burst intensity and the microseismic parameters of the corresponding risk assessment unit are investigated based on the previous studies. And then the correlation analysis of the rock burst intensity and the microseismic parameters are carried out. There will be a significant correlation between the rock burst intensity and a microseismic parameter if their correlation coefficient is big enough. In that case, the microseismic parameter would be used as an indicator of the rock burst risk assessment. On the contrary, it would increase the complexity of the risk assessment model and reduce the risk assessment accuracy. Besides, the indicators used for risk assessment of rock burst should be obtained easily.

According to field destroy situation, rock burst can be classified into intense rock burst, slight rock burst, and none. The main factors considered in the classification are the distance of the flying rock, as well as the depth and volume of the pit, which can be easily obtained. Using the classification criteria, we can collect a database of rock burst cases, which contains intensity and the microseismic parameters of corresponding risk assessment units.

2.3. Probability Distribution Function of Rock Burst of Different Intensities for Each Indicator

It is well known that microseismic indicators reflect a large amount of information on the rock burst development process from different aspects. For each microseismic indicator, there is a relationship between microseismicity and the probability of different rock burst intensities. Generally, the probability distribution function of rock burst of different intensities is established by fitting analysis based on a large number of cases. Here, a linear function is temporarily used to simplify the description of the probability distribution of rock burst due to a lack of cases. We use the average indicator values of rock burst of different intensities as the eigenvalues to conduct the rock burst early risk assessment. The eigenvalue is the central value of the monitored microseismicity for one rock burst intensity. In other words, the probability of one rock burst intensity is 100% when the monitored microseismicity equals the eigenvalue, and the probability is 0 when the monitored microseismicity equals the eigenvalue of an adjacent rock burst intensity. With the increase of deviation between the monitored microseismicity and an eigenvalue, the probability of a corresponding rock burst intensity decreases. If the monitored microseismicity is smaller than the lowest eigenvalue, the rock burst of lowest intensity is likely to occur. Similarly, the probability of an intense rock burst reaches 100% when the monitored microseismicity is bigger than the highest eigenvalue.

2.4. Rock Burst Risk Assessment Formula

As mentioned above, the occurrence of rock burst should be attributed to all the microseismic indicators. Therefore, it should be noted that the rock burst risk assessment should take into account the microseismic indicators according to the correlation analysis. The formula which involves multiple indicators is shown as follows:where i represents rock burst intensity, such as intense, slight, and none, corresponding to 3, 2, and 1, j is microseismic indicator, and Pji is the probability of ith intensity of microseismic indicator j [6].

It is feasible to calculate the weight coefficients of the microseismic indicators in the formula through PCA. Coefficients represent the contribution of the microseismic indicator to the occurrence of rock burst. Perhaps they are not independent. Using PCA, the information contained in the microseismic indicators can be simplified into n new variables, where n is a value is smaller than the number of microseismic indicators. In this way, the ratio of the original information contained in each new variable can be calculated. Then, the weight coefficients of different microseismic indicators can be obtained according to the coefficient in the expressions of the new variables. After normalization processing, the weight coefficients can be taken as the coefficients of the microseismic indicators in the risk assessment formula.

3. Background

The production capability of the Ashele Copper Mine reaches 6000 t per day. The orebody is about 20 to 45 m thick. Its dip ranges from 70° to 80°. The mining depth is close to 1000 m. According to the stress measurement results, the maximum principal stress exceeds 40 MPa at the 0 m level, with a direction that is nearly horizontal. Chalcopyrite is the main mineral, while pyrite and basalt are the surrounding rock. It should be stressed that chalcopyrite and pyrite have a strong rock burst tendency. A large-diameter longhole with the delayed filling open stoping method is dominant in the Ashele Copper Mine. The maximum charge per delay reaches 380 kg, and the maximum charge per blast reaches 5 t. At present, rock burst is a serious hazard in the deep mining area.

A SinoSeism microseismic monitoring system with 28 channels is established in the deep mining area. The microseismic monitoring system is composed of four data acquisition units, nineteen unidirectional sensors, and three triaxial sensors, which are installed from 50 m level to 350 m level. The frequency response range of the sensors is 10∼2000 Hz. When the sampling frequency is 250 Hz, the signal-to-noise ratio is not less than 120 dB. When the sampling frequency is 4000 Hz, the signal-to-noise ratio is not less than 115 dB. Microseismic sensors were installed at the bottom of the hole to avoid the influence of loose zone, which is 6 m deep. The data acquisition unit is connected to a data exchange center by a communication cable, while the data exchange center is connected to the server on the surface by a fiber optic connection. The locating error of the monitoring system is less than 10 m according to the blasting test. Many noise signals are mixed with the rock fracture signals. A neural network model was used to eliminate the blasting vibration signals, mechanical vibration signals, and electrical noise signals [50].

4. Evolution of Microseismic Activity

Microseismic events mainly gathered from 0 m level to 150 m level in the deep mining area of Ashele copper mine, as shown in Figure 1. The size and color represent energy and time of the microseismic events, respectively. Cloud chart distributions of the microseismic event density in the mining area of the Ashele Copper Mine are shown in Figure 2. The events were gathered around the 50 m level in April, 2017, where the roof of the stopes is at the 0 m level. Although some new stopes were added at the 300 m level, there were few changes in the concentrated area and in the clustering degree in September, 2017. The microseismic events appeared in the area ranging from 50 m level to 150 m level and rapidly accompanied the fast growth of the mining stopes growing in December, 2017. The microseismic events continued to be concentrated at the 200 m level in March, 2018. It can be concluded that the mining activities of stopes at the 0 m level and 150 m level have a large impact on the distribution of the microseismic events, and mining activities of stopes at the 300 m level have little impact on it.

The energy index (EI) of a microseismic event is the ratio of the measured radiated microseismic energy generated by the event to the average microseismic energy of all events in the region. It can reflect the stress value indirectly. Cloud chart distributions of lg(EI) based on microseismic events are shown in Figure 3. The area of concentration for lg(EI) is apparently different from that of the event density. The areas of concentration of lg(EI) were distributed from the 50 m level to the 150 m level in April, 2017, where the roofs of the stopes at 0 m level were. As the stopes were mined, the concentration of lg(EI) in the 100 m level to 150 m level range reduced to some extent. Meanwhile, the value of lg(EI) grew rapidly in the 50 m level to 100 m level range, which illustrated that the roof of the stope at the 0 m level was subject to high damage driven stress. The area of concentration continued to decrease, especially near the stope at the 150 m level in December, 2017. Only specific concentration region emerged on the roof of the stope at the 0 m level. The concentration of lg(EI) demonstrated an invisible increase at the 200 m level.

According to the concentrated features of microseismic events and the cloud chart distributions of lg(EI), the high risk area mainly ranges from the 50 m level to the 250 m level. The deep stope deformation monitoring data show that the mining disturbed region is approximately 35 m wide. Microseismic events were concentrated within 30 m from the stope. Based on the projected layout of the stopes, the sublevel haulage drifts on the footwall and hanging wall from the 0 m level to the 350 m level are classified into a group that requires risk assessment work to be covered. The stress reduced area, stress concentration area, and initial stress area are observed around the stope through drift deformation monitoring. In addition, there are many microseismic events in the stress concentration area, where there is less than 30 m from the goaf. Combining the size of the drift, the spatial risk assessment unit, which takes the drift stalk line as the center, is 30 m high, 30 m long, and 30 m wide.

5. Risk Assessment of Rock Burst in Ashele Copper Mine

The rock burst intensity can be divided into three types in the Ashele Copper Mine: intense rock burst (see Figure 4), slight rock burst (see Figure 5), and none rock burst. To consider the temporal evolution process of rock burst, microseismic events rate, microseismic energy rate, and microseismic apparent volume rate use 3 days as the time window, and others use 7 days. The sliding length of the time window is 24 hours. According to the above method, the rock burst with different intensities and the microseismic parameters of the corresponding risk assessment unit are counted. After the correlation analysis between the microseismic parameter and rock burst intensity, six microseismic parameters, such as cumulative event number, cumulative energy, cumulative apparent volume, event rate, energy release rate, and apparent volume rate, are selected to study rock burst risk assessment. The values of the samples can be seen in Figure 6.

The averages of the microseismic indicators for different rock burst intensities are obtained as the eigenvalues. These eigenvalues are the key values for rock burst risk assessment. When the value of a microseismic parameter approaches an eigenvalue of rock burst, the corresponding intensity rock burst is likely to occur. If the value of a microseismic parameter is smaller than the eigenvalue of the lowest level of rock burst, there will be none rock burst. Similarly, the risk of the highest level of rock burst will reach 100% as the value exceeds the highest eigenvalue. In practice, we used the eigenvalue of the adjacent level of rock burst as the boundary for the corresponding rock burst intensity. When the value of Pji is close to 1, the ith intensity rock burst will occur most easily, according to the microseismic parameter j.

In Table 1, Aji represents the eigenvalue of the ith level of rock burst for microseismic indicator j. The probability functions for different rock burst intensities when the value of indicator equals M as follows:

The probabilities of different rock burst intensities for cumulative microseismic events are shown in Figure 7(a). The probabilities for other microseismic parameters can be obtained in the same way.

Before using PCA, we carried out the Kaiser-Meyer-Olkin (KMO) test to judge the feasibility of PCA for the matrix. As a result, the value of KMO is 0.756, which illustrates that PCA can be used to solve the matrix. The total variance of the interpretation is shown in Table 2. Generally, the new principle components of PCA need to contain more than 85% of the original information. In this paper, there are about 91.28% of the original information exiting in the first two principle components. So the first two principle components are selected to calculate the ingredient matrix, which are listed in Table 3. The coefficients of the microseismic parameters list in the new principle components are equal to the ratio of the corresponding ingredient matrix element to the square root of the eigenvalue. Then the two principal component can be expressed as

The proportional contribution of the accumulated microseismic events is calculated as follows:

The microseismic indicator weight model is obtained as follows:

The weight coefficients are determined after normalized processing.

Therefore, the rock burst probability formula for different intensities is as follows:where P1, P2, and P3 represent none rock burst, slight rock burst, and intense rock burst, respectively.

A risk assessment alarm triggered in the north of the footwall haulage drift on 26th March, 2018. The monitored values of the microseismic parameters are shown in Table 4.

The rock burst probabilities for different intensities are calculated as follows:

As a result, the probabilities of the intense rock burst, slight rock burst, and none rock burst in the region are 66.6%, 33.4%, and 0, respectively.

The results show that the probability of the intense rock burst in the region is 66.6% and 33.4% for slight rock burst and 0% for none rock burst. Finally, an intense rock burst occurred in the region on the next day, with a large amount of noise and rock ejection. The result is shown in Figure 8.

The method can predict about 70% of rock burst after two years of application. In practice, there was one rock burst accompanied by no apparent changes in the microseismic indicators, resulting in no warning before its occurrence. In addition, it was difficult to find obvious changes in the microseismic indicators for some rock burst during the establishment of the rock burst database. As a result, we did not pay attention to this type rock burst in this study.

6. Discussion

In the rock burst risk assessment formula, the coefficients of microseismic event rate, cumulative microseismic events, energy release rate, cumulative energy, apparent volume rate, and cumulative apparent volume are 0.180, 0181, 0.191, 0.188, 0.127, and 0.132, respectively, which show the different contributions to the rock burst probability. Obviously, the energy release rate and cumulative energy offer the greatest impact to rock burst from among the microseismic indicators involved in the formula. According to Wang and Park [51], the energy accumulation during the mining process is one of two main impact factors for rock burst. Additionally, Wiles [52] indicated the local energy release rate can be used to predict rock burst because of their strong correlation. However, according to Zhou et al. [36], few approaches based on the single energy index have been found to be particularly successful. The coefficients in the rock burst risk assessment formula may explain the above relationship. The microseismic event rate and cumulative microseismic events are almost equally important, along with the energy release rate and cumulative energy. This is why multiple parameters can be used in combination to avoid the shortcomings of a single parameter.

The application of the proposed method in the Ashele Copper Mine has proven that it is feasible to perform rock burst risk assessment through a combination of a simplified probability distribution function and PCA. With deep mine exploitation, more rock burst cases will be collected, and the eigenvalues of rock burst of different intensities will increase or decrease slightly. By then, a fitting function can be used to replace the liner function to describe the probability of different intensity rock burst. In that case, the successful warning rate would apparently increase. Furthermore, the information contained in the rock burst database should be detailed as much as possible for a future rainy day.

7. Conclusion

The study of rock burst risk assessment and management is important for miners to prevent rock burst. This paper presents a methodology combining a simplified probability function and PCA to allow for early risk assessment of rock burst of different intensities. The microseismic event rate, cumulative microseismic events, microseismic energy rate, cumulative microseismic energy, apparent volume rate, and cumulative apparent volume are taken as microseismic indicators. A liner function was used to describe a simplified probability distribution for rock burst of different intensities for each indicator. Then, the coefficients of the microseismic indicators were calculated by PCA to establish the rock burst risk assessment formula. According to the Ashele Copper Mine rock burst database, the energy release rate, and cumulative energy have the highest contribution to the rock burst probability, followed by the microseismic event rate and cumulative microseismic events, with the apparent volume rate and the cumulative apparent volume last. The application of rock burst risk assessment in the Ashele Copper Mine illustrates that the proposed method is an effective tool for predicting rock burst in metal mines.

Data Availability

The datasets used in the present study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by project 2016YFC0801605 of the National Key Research and Development Plan of China, projects 51704056 and 51974059 of the National Natural Science Foundation of China, and projects N182608003 and N2001001 of the Fundamental Research Funds for the Central University of China.