Next Article in Journal
Robust Current Control of a Small-Scale Wind–Photovoltaic Hybrid System Based on the Multiport DC Converter
Previous Article in Journal
Combining HIIT Plus Cognitive Task Increased Mental Fatigue but Not Physical Workload in Tennis Players
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterizing Rockbursts and Analysis on Hilbert-Huang Transform Spectrum of Microseismic Events, Shuangjiangkou Hydropower Station, Based on Microseismic Monitoring

1
College of Mining, Inner Mongolia University of Technology, Hohhot 010051, China
2
Inner Mongolia Engineering Research Center of Geological Technology and Geotechnical Engineering, Inner Mongolia University of Technology, Hohhot 010051, China
3
Key Laboratory of Geological Hazards and Geotechnical Engineering Defense in Sandy and Drought Regions at Universities of Inner Mongolia Autonomous Region, Inner Mongolia University of Technology, Hohhot 010051, China
4
State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(12), 7049; https://doi.org/10.3390/app13127049
Submission received: 13 May 2023 / Revised: 8 June 2023 / Accepted: 9 June 2023 / Published: 12 June 2023
(This article belongs to the Section Earth Sciences)

Abstract

:
The Shuangjiangkou hydropower station in China has complex geological conditions with high in situ stress. During the tunnel excavation, rockbursts occurred frequently, which seriously affected construction progress. Microseismic (MS) monitoring technology was used to explore rock MS activities to predict rockbursts. The MS monitoring system can capture a large number of MS signals. Based on Hilbert–Huang transform (HHT) instantaneous frequency analysis technology, using MATLAB software (R2022a) to write a program to convert the MS waveform, the frequency and energy characteristics of MS signals at a certain time can be obtained. By analyzing the frequency and energy characteristics of every event, the microseism active areas can be determined, and then rockbursts can be predicted scientifically. This paper selected two different construction sites, which were the main powerhouse and the access tunnel in the main powerhouse, as the research background. Introducing HHT instantaneous time–frequency analysis technology conducted MS event dynamic analysis and predicted rockbursts. The HHT spectrum scientifically and comprehensively displayed MS signal frequency characteristics at a certain time and reflected the change laws of signal instantaneous energy and local abrupt change information. The results indicated that some parameter anomalies in the event spectrum can predict rockbursts. For complex tunnel construction conditions, the HHT time–frequency analysis technology can realize a new idea of using a single-channel signal to predict rockbursts, which was very meaningful.

1. Introduction

Rockburst is a common engineering disaster. In order to more accurately predict rockburst, so as to ensure smooth construction, researchers have used different methods and different angles to study the rockburst process and mechanism. The research methods for rockbursts include theoretical analysis, indoor tests, on-site tests, and numerical simulation. Representative rockburst theories include strength theory, energy theory, stiffness theory, and propensity theory. Researchers use theoretical analysis to study the rockburst mechanism and predict its occurrence. Indoor tests include uniaxial, biaxial, triaxial, and dynamic tests. On-site tests include on-site exploration, stress monitoring, microseismic monitoring, and so on. Researchers use different software to establish different models to analyze rockburst processes. Due to the limitations of experimental conditions and costs, there is relatively little analysis and summary of on-site rockbursts. However, on-site rockburst tests are the most practical and can better reflect the process of rockburst incubation and occurrence. The rockburst mechanism is complex, and rockburst prediction is difficult. Therefore, the prediction and prevention of rockbursts have always been a research hotspot. Naji et al. [1] used empirical approaches to evaluate the proneness of rockburst occurrence. Numerical simulation was used to predict the actual failure zone well. This method can effectively analyze the mechanism of rockburst, but it cannot grasp the internal activity status of surrounding rocks in real time and predict the occurrence of rockburst in a timely manner. Lee et al. [2] summarized the geotechnical issues encountered, including the various modes of rock mass instability and geological hazards, and introduced the current analytical techniques used for prediction in China. This article provided a comprehensive explanation of rock engineering technology. Rehman et al. [3] discussed the major challenges that were encountered during construction herein along with their countermeasures. The research results of this article have an important guiding significance for the excavation of deep-buried tunnels with complex geological environments. Zhang et al. [4] selected four typical rockburst events and described their temporal and spatial characteristics in detail. The research results were of great help in understanding the rockburst mechanism. They had relatively little research on the activity state inside surrounding rocks, which was very important for predicting rockbursts.
During the excavation process, the rock mass will inevitably produce different-scale microcracks, which contain rich information [5,6]. Microseismic monitoring technology is a three-dimensional spatial monitoring method that can form a spatial network structure by arranging a relatively small number of accelerometers, enabling real-time collection of microcrack information inside the rock mass. Then, by analyzing the microcrack information, the microseismic activity laws inside the rock mass can be better understood, revealing the precursor information of rock mass engineering instability and delineating the main damage zones and potential danger zones of the tunnel surrounding rock, which can provide a reference basis for later construction plans. The main characteristics of microseismic monitoring technology include: (1) It can directly determine important parameters such as the time, space, and energy release intensity of microcracks inside the rock masses; (2) Accelerometers can be deployed in areas far away from rock mass that are prone to damage, ensuring the safety of monitoring personnel and equipment, and ensuring that the monitoring system can operate for a long time without being damaged; (3) It can cover a large monitoring area with a limited number of sensors; (4) It can monitor in real time and transmit data remotely.
Microseismic monitoring equipment can capture rock fracture signals in real time. Through in-depth analysis of signal characteristics, we can obtain rich source information, which can be used as the basis for evaluating the stability of rock mass [7]. Compared with other waveform parameters, the frequency-domain characteristics are often unique [8] and relatively stable [9]. The spectrum distribution is an important parameter that can fully reflect the rock failure information carried by the acoustic emission waveform [10]. The MS signal frequency contains rich source information, and the change characteristics of different sequence frequency components can reflect different activity states inside surrounding rocks [11,12]. It has unique advantages and uses the frequency characteristics to conduct rockburst early warning and stability analysis for deep-buried tunnels. Some studies have shown that the high frequency of an MS event corresponds to a small-scale crack [13]. Therefore, the study on the change characteristics of the MS event frequency spectrum has great significance for rockburst early warning and obtaining the MS activity laws.
The spectrum of acoustic emission/MS can better reflect the characteristics of rock damage [14]. An increase in the low frequency of acoustic emission can indicate an imminent rockburst [15,16]. Some studies have shown that lower frequency or a trend of spectrum moving to the low-frequency band is rockburst precursor signals [17,18,19,20,21]. With the rapid development of signal processing technology, a lot of time–frequency analysis theories have been proposed. In recent years, the HHT has been widely used to process nonlinear signals. Huang et al. [22] first applied the HHT method to ocean analysis. Some scholars [23,24,25] applied the HHT to the spectrum research of seismic signals. Khaldi et al. [26] applied the HHT to speech signal denoising with different noise levels, and it was proved that the method was effective. Based on the HHT theory, Mao et al. [27] improved the empirical mode decomposition, improved the decomposition accuracy, and proved the effectiveness of the improved HHT method for low SNR nonstationary signal analysis. Duan [28] applied the HHT to image restoration for the first time and achieved good results. Yu et al. [29] introduced the HHT into gear fault diagnosis for the first time and proposed the local Hilbert energy spectrum concept. Good test results have been obtained. Lian [30] applied the HHT method to seismic phase analysis of seismic signals and achieved good results. Shi et al. [31] proposed a signal analysis method for acoustic nondestructive detection of full-length bonded bolts based on HHT. Yuan et al. [32] used the HHT method to analyze the blasting vibration signal in a lead-zinc mine in South China, and accurately identified the detonator initiation time of seven time delay levels. Veltcheva and Soares [33] applied the HHT method to identify different components of the spectra of real sea waves. Hamdi et al. [34] used the HHT for the extraction of a new relevant damage descriptor to be adopted for acoustic emission (AE) pattern recognition in order to help understand the damage process. Stroeer et al. [35] showed that complementing the HHT with techniques such as zero-phase filtering, kernel density estimation, and Fourier analysis allows it to be used effectively to detect and characterize signals with low signal-to-noise ratios. Litak et al. [36] analyzed the cutting force amplitude and frequency by means of wavelets and HHT. Numerous studies have shown that the HHT method is a relatively scientific signal processing method. Therefore, we can also introduce this method into rockburst prediction.
This paper focuses on rockburst prediction of the tunnels of the Shuangjiangkou hydropower station. The on-site construction environment is very complex, and fracture characteristics inside surrounding rocks are also complex. Therefore, rockbursts occurred frequently and were difficult to predict. The MS signals released contained rich information. The HHT spectrum scientifically and comprehensively displayed MS signal frequency characteristics and reflected the change laws of signal instantaneous energy and local abrupt change information. Compared with the conventional time–frequency analysis method, it can accurately capture the transient characteristics of signals. In this paper, HHT instantaneous frequency analysis technology was introduced. First, HHT instantaneous frequency analysis technology was applied to analyze and identify various on-site signals. Second, through rockburst analysis and prediction for two different sites, it was found that HHT instantaneous frequency analysis technology can achieve a good rockburst prediction effect, which had great significance for on-site construction and safety protection. The HHT time–frequency analysis technology by using a single-channel signal to predict rockbursts can adapt to complex tunnel construction environments.

2. HHT Theory and Signal Recognition

Based on HHT instantaneous frequency analysis technology, using MATLAB software to write a program to convert the MS waveform, the frequency characteristics and amplitude of MS signals at a certain time can be obtained. By analyzing the distribution and change characteristics of frequency and energy over time, the MS activity status inside surrounding rocks was analyzed to predict the possibility of rockburst occurrence.

2.1. HHT Theory

In 1998, Huang et al. [37] proposed a more effective nonlinear signal instantaneous time–frequency conversion technology, namely the Hilbert–Huang transform, which can identify the frequency characteristics of a signal at different time points, and provided a new idea for evaluating the stability of surrounding rock and other equipment failures. The algorithm uses empirical mode decomposition (EMD) signal filtering rules to decompose the signal into several characteristic components in the time domain, and then transform and combine their components. The transformation method can overcome the problem of insufficient signal frequency resolution.

2.1.1. EMD Principle and Algorithm [38]

First, find out all the maximum and minimum points on the time series X(t), and then use the cubic spline function to interpolate the maximum and minimum points of the signal, respectively, to obtain the upper envelope Xmax(t) and lower envelope Xmin(t) of a signal. The mean line m1(t) obtained by connecting the upper and lower envelope lines in turn is as follows:
m 1 ( t ) = [ X max ( t ) + X min ( t ) ] / 2
The first component h1(t) is as follows:
h 1 ( t ) = X ( t ) m 1 ( t )
In the subsequent sifting processes, h1(t) can be treated only as a proto-IMF. Then,
h 11 ( t ) = h 1 ( t ) m 11 ( t )
After repeated siftings in this manner, up to k times, h1k(t) becomes an IMF; that is,
h 1 k ( t ) = h 1 ( k 1 ) ( t ) m 1 ( k 1 ) ( t )
The test is defined as
S D = t = 0 T h 1 ( k 1 ) ( t ) h 1 ( k ) ( t ) h 1 ( k 1 ) 2 ( t ) 2
Huang et al. [24] have proved that when the SD value is between 0.2 and 0.3, it can not only ensure the linearity and stability of IMF, but can also make it have a clear physical significance, and the extraction effect of waveform frequency characteristics is relatively good.
When h1k(t) meets the requirements of SD, the h1k(t) is taken as the first-order IMF.
c 1 ( t ) = h 1 k ( t )
The residual r1(t) is as follows:
r 1 ( t ) = X ( t ) c 1 ( t )
Regard r1(t) as the decomposed signal, and obtain all residuals ri(t) after multiple operations. The result is
r i 1 ( t ) c i ( t ) = r i ( t ) ,   i = 2 , 3 , , n
By summing up (7) and (8), we finally obtain
X ( t ) = i = 1 n c i ( t ) + r n ( t )

2.1.2. The Hilbert Spectral Analysis

The Hilbert spectral analysis method [37] can more intuitively reflect the change characteristics of the instantaneous frequency of nonlinear signals over time. Firstly, every IMF component of the signal is subject to the Hilbert transform, and the corresponding analytic signal A[ci(t)] is obtained as follows:
A [ c i ( t ) ] = c i ( t ) + j [ 1 π c i ( t ) t τ d τ ] = α i ( t ) e j θ i ( t )
where the instantaneous amplitude and instantaneous phase are as follows:
α i ( t ) = c i 2 ( t ) + H 2 [ c i ( t ) ]
θ i ( t ) = arctan H [ c i ( t ) ] c i ( t )
where the instantaneous frequency is as follows:
ω i ( t ) = d θ i ( t ) d t
Then, the IMF component ci(t) can be expressed as:
c i ( t ) = Re [ α i ( t ) exp ( j 2 π f i ( t ) d t ) ]
Show α i ( t ) on the time–frequency plane, and the Hilbert spectrum of ci(t) is obtained.
H i ( t , f ) = α i ( t ) , 0 , f = f i ( t ) f f i ( t )
We finally obtain the analysis signal of the Hilbert spectrum.
x ( t ) = Re [ i = 1 n α i ( t ) exp ( j 2 π f i ( t ) d t ) ]

2.2. Signal Recognition

A drilling and blasting method was adopted at tunnels of the Shuangjiangkou hydropower station, which made the construction noise more complex. The MS system installed can capture a lot of signals. Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5 show some signals in the time domain and their HHT spectrum. The amplitude unit of the amplitude–time curve is set to V (Volt). The vertical coordinate of the HHT spectrogram represents the frequency and the color represents the amplitude. HHT spectrum maps have high temporal resolution and display a lot of effective information.
Figure 1 shows a signal monitored during the drilling hole. Four cyclic waveforms occurred within 800 ms. With the change of frequency over time, there were also four energy concentration areas, which indicated that there were four cyclic waveforms. The frequency distribution (FD) range was about 300–600 Hz. The noise frequency was not obvious in the figure, because the noise amplitude was small, and it was automatically shielded. Due to on-site construction, many drilling signals can be monitored.
An excavation blasting event is shown in Figure 2. The HHT spectrum showed that the rock failure process was extremely complex during excavation blasting. The FD range was about 100–1300 Hz, and the energy distribution of 400 Hz to 900 Hz was relatively concentrated. There was no obvious energy concentration area over time. The frequency points were discrete. Only a part of the blasting excavation signal was intercepted in the figure. There was a clear difference between the blasting signal HHT spectrum and other signals.
Figure 3 shows a signal when a truck whistled. The FD range of the event narrowed over time. The event energy was not obviously distributed below 400 Hz. There was also a small amount of information distribution in relatively high-frequency bands (800–900 Hz). The signal duration is relatively long. The frequency points were discrete.
Figure 4 shows a typical MS signal on 3 January 2018. The information of the MS signal was mainly distributed in 0.06–0.09 s, with a short duration of 0.03 s. The FD range was about 300–600 Hz, the FD had a downward trend over time, and there was an obvious energy concentration area in the figure. The frequency of energy concentration areas was relatively low.
Figure 5 shows the waveform of a rockburst occurred on 13 January 2018. The information of the rockburst signal was mainly distributed in 0.045–0.095 s, with a relatively long duration of 0.05 s. The FD range was about 100–600 Hz, which was relatively wide. In the low-frequency band, the signal information increased significantly; that is, the energy distribution increased significantly. The signal had obvious coda, and its FD range narrowed over time. For rockburst events, the amplitude was larger, the released energy was higher, and the signal component was more complex.

3. Case Analysis

The Shuangjiangkou hydropower station is an important hydropower project located in Sichuan Province, China, as shown in Figure 6. The maximum dam height is 314 m. The hydropower station has complex geological conditions. The diversion and power generation system is arranged on the left bank. The main rock type of tunnels is granite, and the rock uniaxial compressive strength is 60–70 MPa. The maximum principal stress is 37.82 MPa. The total length of the access tunnel in the main powerhouse is 1473.31 m, and cross-sectional dimensions are 11.24 m × 8.7 m (width × height). The powerhouse caverns are mainly buried horizontally at depths of 400 m to 640 m and vertically at depths of 320 m to 500 m. The excavation dimensions of the main powerhouse are 132.56 m in length, 29.30 m in width, and 63.00 m in height. The layout of the main powerhouse and the access tunnel in the main powerhouse is shown in Figure 7. Due to the comprehensive influence of high in situ stress and structural plane during the tunnel excavation, a series of serious damage problems occurred frequently. Among them, the occurrence of rockbursts is the most prominent problem, as shown in Figure 8.
Xue et al. [39] introduced the MS monitoring system of the Shuangjiangkou hydropower station in detail. The main powerhouse (early excavation) and the access tunnel in the main powerhouse are monitored by six accelerometers. At least four accelerometers must be triggered simultaneously to determine the source position. The Paladin data acquisition system has an acquisition frequency of 20 kHz.

3.1. Case Analysis of the Main Powerhouse

On 7 January 2019, two slight rockburst events were monitored at the main powerhouse. The on-site construction progress had slowed down. The HHT spectrum was used to analyze the MS activity laws of the whole rockburst process inside surrounding rocks. The time–frequency information of the rockburst precursor was extracted to provide a scientific basis for efficient construction.
Figure 9a shows an MS signal captured at 00:06:18 on 7 January 2019. The apparent stress value of the event was 16.186 KPa, which was relatively high. The energy released was 170 J. The Es/Ep was 2.72, which indicated that the rock damage type was mainly tensile damage. The FD range of the event was about 200–600 Hz and had a significant downward shift phenomenon over time. The frequency value was relatively low during 0.06–0.08 s. There was some energy concentration in the low-frequency band (200–300 Hz), which indicated there were relatively large-scale cracks inside the rock. The signal lasted approximately 0.08 s.
Figure 9b shows an MS signal captured at 00:08:31 on 7 January 2019. The apparent stress value of the event was 13.394 KPa, which was relatively high. The energy released was 8 J. The Es/Ep was 1.53, which indicated that the rock damage type was mainly tensile damage. The FD range of the event was about 300–700 Hz during 0.01–0.05 s, and was about 300–450 Hz during 0.05–0.065 s. The FD of the event had a downward shift phenomenon over time. The signal lasted approximately 0.055 s.
Figure 9c shows an MS signal captured at 00:12:54 on 7 January 2019. The apparent stress value of the event was 25.308 KPa, which was relatively high. The energy released was 263 J. The Es/Ep was 4.17, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 200–600 Hz during 0.025–0.07 s, and was about 100–300 Hz during 0.07–0.085 s. The FD of the event had a significant downward shift phenomenon over time. The FD range was obviously low, and there was some energy concentration in the low-frequency band. It indicated that the scale of the crack was relatively large. The signal lasted approximately 0.06 s.
Figure 9d shows an MS signal captured at 00:15:04 on 7 January 2019. The apparent stress value of the event was 19.756 KPa, which was relatively high. The energy released was 103 J. The Es/Ep was 3.98, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 300–500 Hz during 0.025–0.07 s, and was about 150–300 Hz during 0.07–0.1 s. The FD of the event had a significant downward shift phenomenon over time. The frequency points had good continuity, and the FD range narrowed. The energy distributions of the event were mainly concentrated within 200–300 Hz. The signal lasted approximately 0.075 s. The frequency spectrum characteristics of the event were typical precursory information of rockburst, indicating that rockburst was likely to occur.
Figure 9e shows a rockburst signal captured at 01:07:33 on 7 January 2019. The apparent stress value of the event was 66.482 KPa, which was high. The energy released was 3060 J. The Es/Ep was 5.67, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 150–400 Hz during 0.025–0.09 s. The frequency energy distribution of the signal had no downward trend over time. The overall FD range of the event was relatively low. The frequency points had good continuity, and the FD range was relatively narrow. The signal lasted approximately 0.065 s. There was a relatively loud “explosion” sound. The time interval between the occurrence of the rockburst event and the previous event was 52 min. It indicated that real-time monitoring had sufficient time to ensure safe construction on site.
Figure 9f shows an MS signal captured at 03:02:48 on 7 January 2019. The apparent stress value of the event was 17.421 KPa, which was relatively high. The energy released was 77 J. The Es/Ep was 5.16, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 200–600 Hz during 0.02–0.065 s, and was about 100–300 Hz during 0.065–0.085 s. The FD range had a downward shift phenomenon over time. The energy concentration was in the low-frequency band. It indicated that after a rockburst occurred, large-scale cracks continued to occur inside surrounding rocks. The signal lasted approximately 0.065 s.
Figure 9g shows an MS signal captured at 03:05:10 on 7 January 2019. The apparent stress value of the event was 11.448 KPa, which was relatively high. The energy released was 52 J. The Es/Ep was 3.45, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 300–600 Hz during 0.025–0.07 s, and was about 100–300 Hz during 0.07–0.085 s. The FD range had a significant downward shift phenomenon over time. There was some energy concentration in the low-frequency band. The signal lasted approximately 0.06 s.
Figure 9h shows an MS signal captured at 03:14:51 on 7 January 2019. The apparent stress value of the event was 5.914 KPa. The energy released was 2 J. The Es/Ep was 11.23, which indicated that the rock damage type was mainly shear damage. The FD range of the event was about 400–700 Hz during 0.025–0.06 s, and was about 100–300 Hz during 0.06–0.07 s. The FD range had a downward shift phenomenon over time. The signal lasted approximately 0.045 s.
Figure 9i shows an MS signal captured at 03:17:13 on 7 January 2019. The apparent stress value of the event was 7.224 KPa. The energy released was 10 J. The Es/Ep was 5.78, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 400–700 Hz during 0.015–0.06 s, and was about 100–300 Hz during 0.06–0.07 s. The FD range had a downward shift phenomenon over time. Some energy concentration was in the low-frequency band. The signal lasted approximately 0.055 s.
Figure 9j shows an MS signal captured at 03:22:21 on 7 January 2019. The apparent stress value of the event was 14.008 KPa, which was relatively high. The energy released was 136 J. The Es/Ep was 1.44, which indicated that the rock damage type was mainly tensile damage. The FD range of the event was about 200–500 Hz during 0.025–0.07 s, and was about 100–400 Hz during 0.07–0.085 s. The FD range had a downward shift phenomenon over time. The overall FD of the event was relatively low. The signal lasted approximately 0.06 s.
Figure 9k shows a slight rockburst signal captured at 03:25:17 on 7 January 2019. The apparent stress value of the event was 44.024 KPa, which was high. The energy released was 1172 J. The Es/Ep was 3.54, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 200–600 Hz during 0.02–0.06 s, and was about 100–300 Hz during 0.06–0.085 s. The FD range had a significant downward shift phenomenon over time. The FD range narrowed and the coda was developed. The signal lasted approximately 0.065 s. There was a relatively loud “explosion” sound.
Before a rockburst, for some events, the FD range had a downward shift phenomenon over time; some energy concentration was in the low-frequency band. The frequency points had good continuity. We need to analyze the HHT spectrum characteristics of every event promptly and in detail. If some events contain rockburst precursor information, it indicates an increased risk of rockburst.

3.2. Case Analysis of the Access Tunnel in the Main Powerhouse

On 13 January 2018, two rockburst events were monitored at the tunnel. There was a relatively loud “explosion” sound on site. The HHT instantaneous frequency analysis technology was used to analyze the rockburst process.
Figure 10a shows an MS signal captured at 07:05:26 on 10 January 2018. The apparent stress value of the event was 77.642 KPa, which was relatively high. The energy released was 154 J. The Es/Ep was 11.19, which indicated that the rock damage type was mainly shear damage. The FD range of the event was about 200–500 Hz during 0.02–0.11 s. The FD range had no downward shift phenomenon over time. The FD range was relatively stable. The energy distribution was relatively dispersed. The frequency points had good continuity. The signal lasted approximately 0.09 s.
Figure 10b shows an MS signal captured at 07:05:36 on 10 January 2018. The apparent stress value of the event was 30.660 KPa, which was relatively high. The energy released was 13 J. The Es/Ep was 6.84, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 300–450 Hz during 0.03–0.12 s. The FD range had no downward shift phenomenon over time. The frequency points had good continuity. The FD range of the event was relatively stable and relatively narrow. The coda was developed, and the energy distribution of the event was relatively dispersed. The signal lasted approximately 0.09 s.
Figure 10c shows an MS signal captured at 08:17:40 on 10 January 2018. The apparent stress value of the event was 59.019 KPa, which was relatively high. The energy released was 102 J. The Es/Ep was 8.06, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 200–600 Hz during 0.025–0.08 s. The frequency points had good continuity. The signal lasted approximately 0.055 s.
Figure 10d shows an MS signal captured at 22:20:58 on 11 January 2018. The apparent stress value of the event was 75.663 KPa, which was relatively high. The energy released was 261 J. The Es/Ep was 5.81, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 250–500 Hz during 0.03–0.09 s. The FD range had no downward shift phenomenon over time. The FD range of the event was relatively stable, and the coda was developed. The signal lasted approximately 0.06 s.
Figure 10e shows an MS signal captured at 02:14:37 on 12 January 2018. The apparent stress value of the event was 62.786 KPa, which was relatively high. The energy released was 79 J. The Es/Ep was 18.25, which indicated that the rock damage type was mainly shear damage. The FD range of the event was about 250–600 Hz during 0.03–0.1 s. The FD range had a downward shift phenomenon over time. The frequency points had good continuity. The FD range of the event narrowed. The coda was developed, and the energy distribution was relatively dispersed. The signal lasted approximately 0.07 s.
Figure 10f shows an MS signal captured at 02:13:49 on 13 January 2018. The apparent stress value of the event was 34.294 KPa, which was relatively high. The energy released was 167 J. The Es/Ep was 4.38, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 100–600 Hz during 0.03–0.13 s. The FD range had a downward shift phenomenon over time. The frequency points had good continuity. The FD range of the event gradually narrowed. The coda was developed, and the energy distribution was relatively dispersed. The signal lasted approximately 0.1 s.
Figure 10g shows an MS signal captured at 03:27:39 on 13 January 2018. The apparent stress value of the event was 24.959 KPa, which was relatively high. The energy released was 43 J. The Es/Ep was 7.74, which indicated that the rock damage type was mainly mixed damage. The FD range of the event was about 200–500 Hz during 0.025–0.13 s. The FD range had a downward shift phenomenon over time. The frequency points had good continuity. The FD range gradually narrowed. The coda was very developed, and the energy distribution of the event was relatively dispersed. The signal lasted approximately 0.105 s.
Figure 10h shows an MS signal captured at 03:38:41 on 13 January 2018. The apparent stress value of the event was 15.936 KPa. The energy released was 64 J. The Es/Ep was 2.06, which indicated that the rock damage type was mainly tensile damage. The FD range of the event was about 200–600 Hz during 0.025–0.07 s, and was about 300–400 Hz during 0.07–0.14 s. The FD range had a downward shift phenomenon over time. The frequency points had good continuity. The FD range gradually narrowed. The coda was very developed. The signal lasted approximately 0.115 s.
Figure 10i shows a slight rockburst signal captured at 04:07:44 on 13 January 2018. The apparent stress value of the event was 58.516 KPa. The energy released was 1154 J. The Es/Ep was 2.91, which indicated that the rock damage type was mainly tensile damage. The FD range of the event was about 100–700 Hz during 0.025–0.065 s. There was obvious energy distribution in the low-frequency band. The signal was relatively complex.
Figure 10j shows a rockburst signal captured at 04:19:06 on 13 January 2018. The apparent stress value of the event was 130.864 KPa. The energy released was 6597 J. The Es/Ep was 1.96, which indicated that the rock damage type was mainly tensile damage. The FD range of the event was about 100–700 Hz during 0.01–0.06 s. There was obvious energy distribution in the low-frequency band. The coda was very developed. The signal was relatively complex.
Before the occurrence of a rockburst, there were several indicators that needed to be noted. The event FD range had a significant downward shift phenomenon over time. The frequency points of some events had good continuity. The energy distribution of some events was in the low-frequency range. The FD range of some events gradually narrowed. Before the rockbursts, only some of these indicators may have occurred. For the comparison of multiple events, there was obvious energy distribution for some events in the low-frequency band over time. The HHT spectrum can well reflect the MS activity characteristics and timely feedback the anomalies inside the surrounding rocks on site.

4. Discussion

The rockburst mechanism is very complex, and it is also very difficult to predict accurately. Therefore, researchers need to explore scientifically by multiple methods and perspectives. Many researchers have found that the Hilbert spectrum can reflect the frequency characteristics and energy characteristics of a waveform at a certain time, and can more comprehensively display the characteristics of an MS waveform to feedback the MS instantaneous activity characteristics inside surrounding rocks.
Rock failure is a very complex process that contains a lot of information, so we need to constantly explore and identify. For MS signals or rockburst signals, due to their different energy released, damage types, and others, there are also some differences in their waveforms and frequency characteristics. Due to different factors inducing rockbursts, the precursor information of rockbursts is also different. According to the field monitoring, many MS events were captured before a rockburst occurred. These MS events contained different information. Some were typical MS events, and some events contained some rockburst precursor information. Therefore, we need to use certain technical means to analyze, identify, and accurately extract the precursor information of rockbursts to accurately predict the occurrence of rockbursts and even predict rockbursts in time and intensity.
The HHT spectrum scientifically and comprehensively displayed MS signal frequency characteristics at a certain time and reflected the change laws of signal instantaneous energy and local abrupt change information. For complex tunnel construction conditions, the HHT time–frequency analysis technology can realize a new idea of using a single-channel signal to predict rockbursts, which was very meaningful. The HHT instantaneous frequency analysis technology was introduced to analyze rockbursts that occurred at the main powerhouse and the access tunnel in the main powerhouse. It was found that before the rockbursts, the frequency of some events had a significant downward trend over time. From the HHT spectrum, the frequency points of some events had good continuity, the energy distribution of some events was in the low-frequency range, the FD range of some events gradually narrowed, and there was some other abrupt information. Rock failure information or rockburst precursor information was very complex. Therefore, rockburst precursor information may be partially or completely abnormal in frequency indicators. For some events, although the frequency has no obvious downward trend over time, the overall FD of events was relatively low. For the longitudinal comparison of multiple events, the energy distribution of some events increased in the relatively low-frequency band over time.

5. Conclusions

The MS system installed at tunnels of the Shuangjiangkou hydropower station can capture a lot of signals. Firstly, MS signal waveforms were recognized by using the HHT. Furthermore, introducing HHT instantaneous time–frequency analysis technology conducted MS event dynamic analysis and predicted rockbursts. The main conclusions are as follows:
(1)
For rockburst events, there was relatively high energy distribution in the low-frequency band, and the signal characteristics were relatively complex. If the energy released by an MS event was relatively high, there was obvious energy distribution for the event in the low-frequency band.
(2)
Before the occurrence of a rockburst, there were several indicators that needed to be noted. The event FD range had a significant downward shift phenomenon over time. The frequency points of some events had good continuity. The energy distribution of some events was in the low-frequency range. The FD range of some events gradually narrowed. Before the rockbursts, only some of these indicators may have occurred, or all of these indicators have occurred. For the comparison of multiple events, the energy distribution of some events increased in a relatively low-frequency band over time.
The HHT spectrum can well reflect the MS activity characteristics and timely feedback the anomalies inside the surrounding rocks on site. In this paper, the HHT instantaneous frequency analysis technology was introduced to reveal the instantaneous frequency change laws of the single-channel signal before rockbursts. A new idea of using a single-channel signal to predict rockbursts was realized.

Author Contributions

Conceptualization, R.X.; methodology, R.X.; software, R.X.; formal analysis, Y.K.; investigation, Y.K.; resources, N.W.; data curation, N.W.; writing—original draft preparation, R.X.; writing—review and editing, N.W.; supervision, N.W.; project administration, R.X.; funding acquisition, R.X. and N.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 42167023, 41941018) and Research Program of Science and Technology at Universities of Inner Mongolia Autonomous Region of China (No. NJZY22363) for which the authors are very grateful.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the anonymous reviewers, who provided valuable suggestions that improved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Naji, A.M.; Emad, M.Z.; Rehman, H.; Yoo, H. Geological and geomechanical heterogeneity in deep hydropower tunnels: A rock burst failure case study. Tunn. Undergr. Space Technol. 2019, 84, 507–521. [Google Scholar] [CrossRef]
  2. Lee, C.; Sijing, W.; Zhifu, Y. Geotechnical aspects of rock tunnelling in China. Tunn. Undergr. Space Technol. 1996, 11, 445–454. [Google Scholar] [CrossRef]
  3. Rehman, H.; Naji, A.M.; Nam, K.; Ahmad, S.; Muhammad, K.; Yoo, H.K. Impact of Construction Method and Ground Composition on Headrace Tunnel Stability in the Neelum–Jhelum Hydroelectric Project: A Case Study Review from Pakistan. Appl. Sci. 2021, 11, 1655. [Google Scholar] [CrossRef]
  4. Zhang, C.; Feng, X.-T.; Zhou, H.; Qiu, S.; Wu, W. Case histories of four extremely intense rockbursts in deep tunnels. Rock Mech. Rock Eng. 2012, 45, 275–288. [Google Scholar] [CrossRef]
  5. Wang, S.F.; Li, X.B.; Yao, J.R.; Gong, F.Q.; Li, X.; Du, K.; Tao, M.; Huang, L.Q.; Du, S.L. Experimental investigation of rock breakage by a conical pick and its application to non-explosive mechanized mining in deep hard rock. Int. J. Rock Mech. Min. Sci. 2019, 122, 104063. [Google Scholar] [CrossRef]
  6. Huang, L.Q.; Xin, W.U.; LI, X.B.; Wang, S.F. Influence of sensor array on MS/AE source location accuracy in rock mass. Trans. Nonferrous Met. Soc. China. 2023, 33, 254–274. [Google Scholar] [CrossRef]
  7. Zhang, Y.P. Rock mass acoustic emission and microseismic monitoring and localizing technology and its application. Eng. Blasting 2002, 8, 58–61. [Google Scholar]
  8. He, M.C.; Zhao, F.; Zhang, Y.; Du, S.; Guan, L. Feature evolution of dominant frequency components in acoustic emissions of instantaneous strain-type granitic rockburst simulation tests. Rock Soil Mech. 2015, 36, 1–8. [Google Scholar]
  9. Ma, C.X.; Zhang, D.L. Acoustic emission Studies on rock cutting. Chin. J. Rock Mech. Eng. 1996, 15, 77–83. [Google Scholar]
  10. Pan, C.L.; Zhu, F.C.; Cao, P.; Feng, T. Characteristics of acoustic emission of bursting-intended rocks under uniaxial compression. J. Cent. South Univ. Technol. 2001, 32, 336–338. [Google Scholar]
  11. Lu, C.P.; Liu, G.J.; Liu, Y.; Zhang, N.; Xue, J.H.; Zhang, L. Microseismic multi-parameter characteristics of rockburst hazard induced by hard roof fall and high stress concentration. Int. J. Rock Mech. Min. Sci. 2015, 76, 18–32. [Google Scholar] [CrossRef]
  12. He, M.C.; Miao, J.L.; Feng, J.L. Rock burst process of limestone and its acoustic emission characteristics under true-triaxial unloading conditions. Int. J. Rock Mech. Min. Sci. 2010, 47, 286–298. [Google Scholar] [CrossRef]
  13. Lu, C.P.; Dou, L.M.; Liu, B.; Xie, Y.S.; Liu, H.S. Microseismic low-frequency precursor effect of bursting failure of coal and rock. J. Appl. Geophys. 2012, 79, 55–63. [Google Scholar] [CrossRef]
  14. Zhu, Z.F.; Chen, G.Q.; Xiao, H.Y.; Liu, H.; Zhao, C. Study on crack propagation of rock bridge based on multi parameters analysis of acoustic emission. Chin. J. Rock Mech. Eng. 2018, 37, 909–918. [Google Scholar]
  15. Armstrong, B.H. Acoustic emission prior to rockbursts and earthquakes. Bull. Seismol. Soc. Am. 1969, 59, 1259–1279. [Google Scholar]
  16. Calder, P.N.; Madsen, D. High frequency precursor analysis prior to a rockburst. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1989, 26, 3–4. [Google Scholar]
  17. Lu, C.P.; Dou, L.M.; Liu, H.; Liu, H.S.; Liu, B.; Du, B.B. Case study on microseismic effect of coal and gas outburst process. Int. J. Rock Mech. Min. Sci. 2012, 53, 101–110. [Google Scholar] [CrossRef]
  18. Lu, C.P.; Dou, L.M.; Zhang, N.; Xue, J.H.; Wang, X.N.; Liu, H.; Zhang, J.W. Microseismic frequency-spectrum evolutionary rule of rockburst triggered by roof fall. Int. J. Rock Mech. Min. Sci. 2013, 64, 6–16. [Google Scholar] [CrossRef]
  19. Lu, C.P.; Dou, L.M.; Zhang, N.; Xue, J.H.; Liu, G.J. Microseismic and acoustic emission effect on gas outburst hazard triggered by shock wave: A case study. Nat. Hazards 2014, 73, 1715–1731. [Google Scholar] [CrossRef]
  20. Meng, X.X.; Liu, W.T.; Zhang, H. A study of rock burst hazard prevention method. Geotech. Geol. Eng. 2018, 36, 2237–2246. [Google Scholar] [CrossRef]
  21. Yuan, R.F.; Li, H.M.; Li, H.Z. Distribution of microseismic signal and discrimination of portentous information of pillar type rockburst. Chin. J. Rock Mech. Eng. 2012, 31, 80–85. [Google Scholar]
  22. Huang, N.E.; Shen, Z.; Long, S.R. A new view of nonlinear water waves: The Hilbert spectrum. Annu. Rev. Fluid Mech. 1999, 31, 417–457. [Google Scholar] [CrossRef] [Green Version]
  23. Huang, N.E. Computer implemented empirical mode decomposition method, apparatus, and article of manufacture for two-dimensional signals. U.S. Patent Pending 6,311,130 B1, 30 October 2001. [Google Scholar]
  24. Huang, N.E.; Wu, Z. A review on Hilbert-Huang transform: Method and its applications to geophysical studies. Rev. Geophys. 2008, 46, 2–4. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, R.R. Characterizing and quantifying earthquake-induced site nonlinearity. Soil Dyn. Earthq. Eng. 2006, 26, 799–812. [Google Scholar] [CrossRef]
  26. Khaldi, K.; Alouane, M.T.; Boudraa, A. Voiced speech enhancement based on adaptive filtering of selected intrinsic mode functions. Adv. Adapt. Data Anal. 2010, 2, 65–80. [Google Scholar] [CrossRef] [Green Version]
  27. Mao, W.; Jin, R.H.; Geng, J.P.; Li, J.Q. A time frequency analysis method for nonstationary signal based on improved Hilbert-Huang Transform and its application. J. Shanghai Jiaotong Univ. 2006, 5, 724–729. [Google Scholar]
  28. Duan, X.T. Research and Application of Hilbert-Huang Transform. Master’s Thesis, The PLA Information Engineering University, Zhengzhou, China, 2005. [Google Scholar]
  29. Yu, D.J.; Cheng, J.S.; Yang, Y. Application of Hilbert-Huang Transform method to gear fault diagnosis. Chin. J. Mech. Eng. 2005, 6, 102–107. [Google Scholar] [CrossRef]
  30. Lian, H.N. Study on application of Fourier, Wavelet and Hilbert-Huang Transform in earthquake Engineering. Master’s Thesis, Institute of Engineering Mechanics, Chinese Seismological Bureau, Harbin, China, 2004. [Google Scholar]
  31. Shi, Z.M.; Liu, L.; Peng, M.; Liu, C.C.; Tao, F.J.; Liu, C.S. Non-destructive testing of full-length bonded rock bolts based on HHT signal analysis. J. Appl. Geophys. 2018, 151, 47–65. [Google Scholar] [CrossRef]
  32. Yuan, H.P.; Liu, X.L.; Liu, Y.; Bian, H.B.; Chen, W.; Wang, Y.X. Analysis of Acoustic Wave Frequency Spectrum Characters of Rock Mass under Blasting Damage Based on the HHT Method. Adv. Civ. Eng. 2018, 2018, 9207476. [Google Scholar] [CrossRef]
  33. Veltcheva, A.D.; Soares, C.G. Identification of the components of wave spectra by the Hilbert Huang transform method. Appl. Ocean Res. 2004, 26, 1–12. [Google Scholar] [CrossRef]
  34. Hamdi, S.E.; Duff, A.L.; Simon, L.; Plantier, G.; Sourice, A.; Feuilloy, M. Acoustic emission pattern recognition approach based on Hilbert-Huang transform for structural health monitoring in polymer-composite materials. Appl. Acoust. 2013, 74, 746–757. [Google Scholar] [CrossRef]
  35. Stroeer, A.; Cannizzo, J.K.; Camp, J.B.; Gagarin, N. Methods for detection and characterization of signals in noisy data with the Hilbert-Huang transform. Phys. Rev. D 2009, 79, 124022. [Google Scholar] [CrossRef] [Green Version]
  36. Litak, G.; Kecik, K.; Rusinek, R. Cutting force response in milling of Inconel: Analysis by wavelet and Hilbert-Huang transforms. Lat. Am. J. Solids Struct. 2013, 10, 133–140. [Google Scholar] [CrossRef]
  37. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.A.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. London. Ser. A Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  38. Huang, N.E. Hilbert-Huang Transform and Its Applications. World Sci. 2014, 16. [Google Scholar]
  39. Xue, R.X.; Liang, Z.Z.; Xu, N.W.; Dong, L.L. Rockburst prediction and stability analysis of the access tunnel in the main powerhouse of a hydropower station based on microseismic monitoring. Int. J. Rock Mech. Min. Sci. 2020, 126, 104174. [Google Scholar] [CrossRef]
Figure 1. Waveform analysis of a drilling hole signal in the time domain and the HHT spectrum.
Figure 1. Waveform analysis of a drilling hole signal in the time domain and the HHT spectrum.
Applsci 13 07049 g001
Figure 2. Waveform analysis of an excavation blasting signal in the time domain and the HHT spectrum.
Figure 2. Waveform analysis of an excavation blasting signal in the time domain and the HHT spectrum.
Applsci 13 07049 g002aApplsci 13 07049 g002b
Figure 3. Waveform analysis of a truck horn signal in the time domain and the HHT spectrum.
Figure 3. Waveform analysis of a truck horn signal in the time domain and the HHT spectrum.
Applsci 13 07049 g003
Figure 4. Waveform analysis of an MS event signal in the time domain and the HHT spectrum.
Figure 4. Waveform analysis of an MS event signal in the time domain and the HHT spectrum.
Applsci 13 07049 g004
Figure 5. Waveform analysis of a rockburst signal in the time domain and the HHT spectrum.
Figure 5. Waveform analysis of a rockburst signal in the time domain and the HHT spectrum.
Applsci 13 07049 g005
Figure 6. Geographical location.
Figure 6. Geographical location.
Applsci 13 07049 g006
Figure 7. Layout of tunnels.
Figure 7. Layout of tunnels.
Applsci 13 07049 g007
Figure 8. Photos of rock (granite) damage: (a) main powerhouse; (b) access tunnel in the main powerhouse.
Figure 8. Photos of rock (granite) damage: (a) main powerhouse; (b) access tunnel in the main powerhouse.
Applsci 13 07049 g008
Figure 9. HHT spectrum of MS events: (a) 00:06, 7 January 2019; (b) 00:08, 7 January 2019; (c) 00:12, 7 January 2019; (d) 00:15, 7 January 2019; (e) 01:07, 7 January 2019; (f) 03:02, 7 January 2019; (g) 03:05, 7 January 2019; (h) 03:14, 7 January 2019; (i) 03:17, 7 January 2019; (j) 03:22, 7 January 2019; (k) 03:25, 7 January 2019.
Figure 9. HHT spectrum of MS events: (a) 00:06, 7 January 2019; (b) 00:08, 7 January 2019; (c) 00:12, 7 January 2019; (d) 00:15, 7 January 2019; (e) 01:07, 7 January 2019; (f) 03:02, 7 January 2019; (g) 03:05, 7 January 2019; (h) 03:14, 7 January 2019; (i) 03:17, 7 January 2019; (j) 03:22, 7 January 2019; (k) 03:25, 7 January 2019.
Applsci 13 07049 g009aApplsci 13 07049 g009bApplsci 13 07049 g009cApplsci 13 07049 g009d
Figure 10. HHT spectrum of MS events: (a) 07:05, 10 January 2018; (b) 07:05, 10 January 2018; (c) 08:17, 10 January 2018; (d) 22:20, 11 January 2018; (e) 02:14, 12 January 2018; (f) 02:13, 13 January 2018; (g) 03:27, 13 January 2018; (h) 03:38, 13 January 2018; (i) 04:07, 13 January 2018; (j) 04:19, 13 January 2018.
Figure 10. HHT spectrum of MS events: (a) 07:05, 10 January 2018; (b) 07:05, 10 January 2018; (c) 08:17, 10 January 2018; (d) 22:20, 11 January 2018; (e) 02:14, 12 January 2018; (f) 02:13, 13 January 2018; (g) 03:27, 13 January 2018; (h) 03:38, 13 January 2018; (i) 04:07, 13 January 2018; (j) 04:19, 13 January 2018.
Applsci 13 07049 g010aApplsci 13 07049 g010bApplsci 13 07049 g010cApplsci 13 07049 g010d
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xue, R.; Kong, Y.; Wu, N. Characterizing Rockbursts and Analysis on Hilbert-Huang Transform Spectrum of Microseismic Events, Shuangjiangkou Hydropower Station, Based on Microseismic Monitoring. Appl. Sci. 2023, 13, 7049. https://doi.org/10.3390/app13127049

AMA Style

Xue R, Kong Y, Wu N. Characterizing Rockbursts and Analysis on Hilbert-Huang Transform Spectrum of Microseismic Events, Shuangjiangkou Hydropower Station, Based on Microseismic Monitoring. Applied Sciences. 2023; 13(12):7049. https://doi.org/10.3390/app13127049

Chicago/Turabian Style

Xue, Ruixiong, Yinghui Kong, and Na Wu. 2023. "Characterizing Rockbursts and Analysis on Hilbert-Huang Transform Spectrum of Microseismic Events, Shuangjiangkou Hydropower Station, Based on Microseismic Monitoring" Applied Sciences 13, no. 12: 7049. https://doi.org/10.3390/app13127049

APA Style

Xue, R., Kong, Y., & Wu, N. (2023). Characterizing Rockbursts and Analysis on Hilbert-Huang Transform Spectrum of Microseismic Events, Shuangjiangkou Hydropower Station, Based on Microseismic Monitoring. Applied Sciences, 13(12), 7049. https://doi.org/10.3390/app13127049

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop