Abstract

The macroscopic deformation and failure of engineering rock mass may occur as a result of evolution and breakdown of its internal microfracture. Therefore, the macroscopic state of rock mass can be obtained from fracture scale of microfracture in real time. To assess instability and predict macroscopic deformation and failure of engineering rock mass, a time-frequency analysis technique based on S transform was proposed to investigate microseismic waveform and reveal the correlation between macroscopic deformation failure and microseismic frequency characteristics of engineering rock mass in combination with fracture scale. To minimize the influence of external factors on parameters calculated, a significant amount of microseismic data from three large-scale hydropower projects in southwestern China was collected as the statistical sample. The analysis of correlation between fracture scale and frequency characteristics of microseismic events was carried out based on the statistical sample. Combining with microseismic data and multipoint extensometers in the underground powerhouse of the Houziyan hydropower station, engineering verification was conducted. The result shows that the high-frequency components decrease and microseismic signals display low-frequency characteristic as the fracture scale increases; the microseismic high-frequency components decreased at first and then increased during the deformation process of surrounding rock mass, and the frequency of microseismic events shifts from high band to a lower one before deformation.

1. Introduction

The occurrence and development of rock microfracture release energy in the form of elastic waves, that is, acoustic emissions (AEs) or microseismicity (MS). Researches show that each AE or MS contains abundant seismic source information. The analysis of AE/MS can reveal the evolution of the cracks occurred inside the rock masses and then serve as the judgment basis for the internal damage of the rock masses [13].

The analysis methods of AE/MS include the parameters analysis method and waveform analysis method [4]. Parameters analysis method analyzes the source characteristics by statistical parameters of the signals, such as event rate, number of events, energy rate, and b-value, which are widely used in the experimental study of AE. Zhang et al. [5] studied the AE characteristics of rock failure under uniaxial multilevel loadings by statistical parameters of the signals, such as AE event rate, number of events, and energy rate. Li et al. [6] revealed the relationship among number of AE events, event rate stress, and the time of uniaxial compression experiment. Li et al. [7] found the changing trend of the b-value of AEs and their spatial distribution of the fractal value with stress level in the processes of rock failure under uniaxial loading. The researches above have made significant progresses in the mechanism of rock failure, but these parameters only describe a single characteristic of the signal, providing insufficient information about the source characteristics and making great casualness in the choice of parameters [4, 8]. Waveform analysis method uses signal processing methods to analyze the recorded waveform in time domain and reveals the seismic source information. The AE/MS signal can be defined as an instantaneous elastic wave [5]. Waveform analysis method can thus reflect the failure information of rock mass in real time. In the propagation process, the amplitude and duration of the waves are uncertain due to the attenuation, superposition, and interference. However, the frequencies of the waveforms still keep intrinsic characteristics and uniqueness [9] and are more stable [10]. Different frequency characteristics can reflect the different failure processes within the rock mass [11]. Some scholars have conducted extensive studies on the frequency characteristics of AE signals of rock. The indoor rockburst test of limestone has been conducted by He et al. [12] to study the frequency characteristics of AE signals before the occurrence of rockbursts. Lu et al. [13] analyzed the frequency-spectrum evolvement rule of AE signals in the entire processes from deformation to fracture by testing the compound coal-rock samples obtained from the Sanhejian coal mine. Li et al. [14] studied the AE frequency characteristics of four kinds of rock under seepage and without seepage by the indoor uniaxial compression experiment. The frequency characteristics of signals in earthquake and MS monitoring are also concerned by researchers. Neuberg et al. [15] proposed the use of low-frequency earthquakes for volcanic eruption prediction by studying low-frequency earthquakes in volcanic earthquakes and their induced mechanisms. By analyzing the MS monitoring signals prior to and after rockburst, Lu et al. [16] pointed out that the frequency spectrum of the picked-up signals changed from multipeak to single peak before rockburst and the low-frequency component of signals increased significantly. It is also suggested that the low-frequency phenomenon can be used to predict rockburst. The time-frequency analysis of rockburst and blasting vibration signals in mine MS monitoring have been conducted by Zhao et al. [17] to identify the different waveforms. The researches above mainly focus on single-frequency characteristics. However, few studies have been done on seismic source mechanism combining seismic source parameters and frequency characteristics of MS events obtained from engineering rock mass.

The present study adopts the time-frequency analysis technique based on S transform to process a significant amount of pickup waveforms of MS events obtained from three hydropower projects in southwestern China. The correlation between source fracture scale and frequency characteristics of MS events was analyzed. The evolution regularity of MS signals frequency during the surrounding rock mass deformation was carried out compared with the results of the multipoint extensometers in the region of MS events concentration. Some reference can be provided for the stability analysis and hazards prediction of engineering rock mass.

2. Waveform Analysis Method

The waveform analysis methods such as fast Fourier transform (FFT) [18], short-time Fourier transform (STFT) [19], and wavelet transform (WT) [20] are commonly used in the signals of indoor AE and on-site MS monitoring [4]. FFT is a conversion tool to transform signals from time domain to frequency domain, but it does not have localization ability in time domain and cannot accurately describe the frequency characteristics of nonstationary signals. The STFT can analyze the signals in the time-frequency domain under the selected window function and time width, but its window width is fixed. Therefore, it is still not able to overcome the limitation of FFT in the nonstationary signals analysis, and the choice of window function and analysis time width plays a significant role on the result of signals analysis. The WT is to decompose the signals at different scales by the selected wavelets function, and it can decompose nonstationary signals. However, the WT decomposes one-dimensional signals into time-scale domains rather than time-frequency domains, which cannot provide intuitive time-frequency characteristic. Moreover, the selection of wavelet function and analysis scale has a great influence on the analysis result, resulting in the limitation of WT in time-frequency analysis of MS signals.

S transform (ST) is a time-frequency analysis method proposed by geophysicist Stockwell in 1996 [21]. It can be defined aswhere is the time, which controls the position of the window function on the time axis, is the analysis signals, is the frequency, and is the ST-time spectrum matrix after transformation.

ST can be written in the form of Fourier spectrum H(f):

With regard to discrete signals, Fourier spectrum can be obtained by the following formula:where k = 0, 1, …, N − 1 represents discrete time points, is the discrete signal length, and is the sampling time interval.

For discrete signals, assuming , the S transform of the final discrete signal can be expressed as

Compared with STFT, ST overcomes the defect of window time width invariability and can self-adaptively adjust the analysis time width according to variation of frequency. ST can provide intuitive time-frequency characteristics without selecting window functions and analysis scales [1921]. Owing to these advantages, ST has been widely used in time-frequency analysis of signals. In the present study, ST time-frequency analysis technology is adopted for MS signals waveform analysis, which can provide a new idea for waveform analysis.

The MS signals are transformed by ST as shown in Figure 1, in which the amplitude is normalized amplitude. Figure 1(a) shows the time-domain waveform of MS signals. Figures 1(b) and 1(c) are the time-frequency spectrum and time-frequency amplitude spectrum of the signals, respectively. The MS signal frequency parameters have been marked in Figure 1(b). The dominant frequency is the frequency band which high-amplitude part of the time-frequency distributes, and the center of the dominant frequency band is the center frequency.

3. Correlation between Fracture Scale and Frequency Characteristics of MS Signals

3.1. Processing of MS Monitoring Data

The elastic waves released from rock fractures will suffer from attenuation due to the complexity of the media in the propagation processes. Thus, there are some deviations in the calculation of seismic source parameters of single MS event. However, when the number of statistical samples increases, the differences caused by such attenuation and manual handling will be minimized as much as possible [22]. In order to minimize the influence of attenuation in the propagation of MS signals, the waveforms with the largest amplitude channels are selected for waveform analysis. According to the above reasons, the relationship between fracture scale and frequency characteristics of picked-up signals is statistically analyzed based on the abundant MS data obtained from the left bank slope of the Jinping first-stage hydropower station during excavation [23] and impoundment [24], the underground powerhouse of the Houziyan hydropower station [25, 26], and the underground powerhouse of the Baihetan hydropower station [27]. The numbers of statistical samples for each research object are shown in Table 1.

3.1.1. Fracture Scale

The MS monitoring systems produced by ESG Corporation, Canada, are adopted for each research object in Table 1. Based on the assumption of ring fault plane, the source fracture zone is equivalent to a disk with a radius of r, which represents the rupture scale at the focal point, namely, seismic source radius [28, 29]. MS monitoring system can effectively capture microfractures prior to macroscope deformation and failure of rock mass. The seismic source radius of MS events will show an increasing trend before the rock mass deformation and failure [30]. Therefore, the variation of seismic source radius can reflect some characteristics of macroscope deformation and failure of rock mass. However, the calculation of seismic source radius may be inaccurate due to attenuation and superposition of signals and artificial waveform processing, in particular, when the seismic source radius of MS events in short-time periods and small areas are selected for predicting the deformation and failure of rock mass. More stable frequency parameters of MS signals are used as early warning indicators in the present study. The prediction may be realized through establishing the relationship between signal frequency and deformation and failure of rock mass based on the seismic source radius.

The Madariaga model [29] is used to calculate the seismic source radius parameters. Field observations at the Strathcona mine in Canada show that compared to the Brune model [28], the calculated value using the Madariaga model is more suitable to practical engineering [31]. Its expression iswhere is a constant determined by the seismic source model and is a function of the angle between the perpendicular of the fault plane and the source-sensor ray. When the fault plane is undetermined, for the P wave and for the S wave. is the velocity of the P wave or S wave in the seismic source region, and is the corner frequency. Its formula [32] iswhere and are the velocity power spectrum and displacement power spectrum of the particle, respectively.

3.1.2. Waveform Processing

The ST time-frequency analysis technique can convert the one-dimensional signals into the time-frequency spectrum in the time-frequency domain and describe the temporal variation rule of the signal frequency. The energy of the MS signal is usually distributed in multiple frequency bands, and the ST technique can accurately distinguish the frequency components contained in the MS pickup signals. For instance, a MS event and its time-frequency spectrum based on the ST technique at the left bank slope of the Jinping first-stage hydropower station during impoundment at 0:30 on August 29, 2014, are shown in Figures 1(a) and 1(b). The four highlighted areas in the time-frequency spectrum reflect that the waveform mainly contains four frequency components. The main frequency ranges of the four regions are calibrated in the time-frequency domain, and the center frequency of each region are calculated as shown in Table 2.

The MS signals of each engineering project in Table 1 are processed by the abovementioned method. Seven hundred and ninety-eight, 641, 855, and 495 frequency samples are obtained for the Jinping first-stage project during excavation and impoundment, Houziyan project, and Baihetan project, respectively.

3.2. MS Monitoring Data Analysis

As we know, AE/MS signals with different fracture scales have different frequency ranges [33]. With regard to natural earthquake, large-scale deformation and failure occur in large-scale rock mass under the action of the geological structure. The frequencies of seismic signals are distributed in low-frequency bands. Under laboratory conditions, AE frequencies of rock samples are usually distributed in high-frequency bands. Ohnaka [34] and Haskell [35] found the relationship between the frequency of AE signals and its crack length as follows:where is the time in which the crack closes, opens, or slides; is the length of crack; is the speed in which the crack generates; is the corner frequency of the elastic wave in which the crack generates; and is the angle between the propagation direction of the signal and the crack surface. When is π/2, (7) can be simplified as

It can be found that large-scale crack development will generate AE signals with low frequency. Meanwhile, Ohnaka and Mogi [36] analyzed the frequency characteristics of AE generated by both rocks under uniaxial compression, and they found that the proportion of low-frequency AE signals obviously increased when the rock was near destruction. They believed that this was related to the low-frequency signal generated by the large cracks that occurred in the rock when rock was adjacent to failure. The similar AE regularities were also found in many rock [9, 12, 37] and concrete [38] laboratory tests.

To analyze the frequency characteristics of MS signals in the engineering rock mass, the center frequencies of MS signals in Table 1 are selected as the frequency parameters. A statistical analysis of seismic source radius and center frequency of MS signals in the projects above is conducted. The statistical results are shown in Figure 2. Figure 2(a) shows that the crack scales of MS events occurred in the left bank slope of the Jinping project during excavation range from 2.9 to 7.3 m and mainly concentrated between 3.3 and 4.5 m. Figure 2(b) shows that the crack scales of MS events in the left bank slope of the Jinping project during impoundment range from 1.9 to 5.2 m and concentrate between 2.2 and 4.4 m. Comparing Figures 2(a) and 2(b), it can be deduced that the deformation of the left bank slope of the Jinping project during impoundment may be smaller than that during excavation. Compared with the excavation period, the frequency of MS signals during impoundment is distributed in a wider range. Figures 2(c) and 2(d) show the relationships between source radius and center frequency of MS signals obtained from the Houziyan project and Baihetan project during the statistical periods, respectively. Compared with slope engineering, the crack scales of MS events in underground powerhouses have a wider range. The source radiuses of MS events in the Houziyan project range from 2.3 to 38.7 m and concentrate between 3 and 11 m. The source radiuses of MS events in the Baihetan project range from 4.3 to 39.7 m and concentrate between 5 and 15 m.

Comparing the relationship between source radius and center frequency of each research object in Figure 2, the signal frequency has a wider distribution range when the fracture size of MS events is small in both slope engineering and underground powerhouse engineering. As the size of the rupture increases, the range of frequency distribution of MS signals gradually decreases. The high-frequency MS signals also gradually decrease, and the MS signals are characterized by low-frequency signals. These results are consistent with the phenomena observed in indoor AE tests by Ohnaka and Mogi [36].

The results above indicate that the frequency of AE signals is clearly related to the deformation of rock mass. Namely, as the fracture scale increases, the high-frequency component of AE signals decreases and the signals are characterized by low frequency. The failure of rock is the process of initiation, developing, transfixion, and forming a large-scale fracture of internal microfractures, eventually leading to macrodeformation and failure [12]. Prior to the rock failure, rock fracture is characterized by single-microscale random fractures and AE signals show high-frequency characteristics. When the rock is adjacent to failure, the microfractures continuously develop and interpenetrate to form large-scale fractures, and the signal frequency is characterized by low frequency. Therefore, the deformation behavior of rock mass with the scope of interest can be investigated through analyzing the variation of frequency characteristics of MS events in the selected period, which will provide some references for stability evaluation of engineering rock mass.

4. Case Study

4.1. Project Overview

The Houziyan hydropower station, located in Kangding County, Ganzi Prefecture in Sichuan Province, is the ninth cascade hydropower station on upstream of the Dadu River. The underground plant system is located in the right bank slope. The dimensions of the main powerhouse are 219.50 m in length, 29.20 m in width, and 68.70 m in height. The elevation of the top arch is 1732.225 m. The underground powerhouse is mainly buried horizontally at depths of 280 to 510 m and vertically at depths of 400 to 660 m. It is a typical deep underground powerhouse with complex lithology and large-scale faults [25, 26].

Due to great depth, high geostress, and complex geologic structures, the problem of the multicavity effect is prominent. A large amount of rib spalling and rockbursts occurred subject to excavation unloading in the underground caverns of the Houziyan hydropower station (Figure 3). A real-time MS monitoring system was thus adopted to assess the stability of underground caverns during excavation. The installations of the MS monitoring system and previous achievements about the underground powerhouse of the Houziyan project can be referred to Li et al. [25, 26].

4.2. Correlation between MS Events Accumulation and Surrounding Rock Mass Deformation

Figure 4 shows the absolute displacement curves of the multipoint extensometers M6XZ1-2 and M4CF3-8 installed between omnibus bar caves #1 and #3 in the downstream of main powerhouse. More than 10 mm displacement of multipoint extensometers M6XZ1-2 and M4CF3-8 occurred on the beginning of December 2013 and January 2014, respectively. Combined with the construction dynamic condition, the phenomenon can be interpreted that the deformation was primarily attributed by the excavation of layer V and VI of the main powerhouse. The excavation exposure is shown in Figure 5.

The MS activities characteristics of the main powerhouse between November 26 and December 14, 2013, and between December 23 and January 10, 2014, are shown in Figure 6, respectively, which corresponds to the two stages of surrounding rock mass deformation in Figure 4(b). It can be found that MS events concentrated near the region of multipoint extensometer M4CF3-8 obviously and showed great consistence in temporal distribution. The phenomenon demonstrates that MS activity can reflect the visual deformation measured by multipoint extensometers accurately. Namely, the deformation of surrounding rock mass usually lags behind the concentration of MS events. Therefore, the variation of MS events can be regarded as the basis for judging the deformation of surrounding rock mass.

4.3. Frequency Evolution Characteristics of MS Signals during Deformation Processes of Surrounding Rock Mass

In the deformation period of surrounding rock mass above, the MS events accumulated in surrounding rock mass near multipoint extensometers are selected as a research object. The MS signals are processed based on the ST waveform processing method, and the frequency variation characteristics during the deformation processes of multipoint extensometers are also analyzed. Figure 7 shows temporal distribution regularities of the center frequency of MS events concentrated near multipoint extensometers, and the displacement process curve of multipoint extensometers M4CF3-8, respectively. From the displacement curve of multipoint extensometer M4CF3-8 in Figure 7(b), surrounding rock mass deformation occurred two times in the selected period, namely, between November 28 and December 12, 2013, and between January 2 to 9, 2014, respectively. Figure 7(a) shows that the high-frequency components of MS signals underwent two ascending phases and three dropping phases during the selected period. Combined with the results of surrounding rock mass deformation, the high-frequency components of center frequencies of MS signals in the concentrated zone began to decrease on November 28, 2013. The results in Section 3.2 show that the decrease of high-frequency signals indicates the increase of source rupture scale. Coincidentally, surrounding rock mass began to deform on November 28, 2013. Afterwards, several high-frequency signals appeared, the ratio of high-frequency components began to increase, but there were still a large number of high-frequency components missing. Since December 8, 2013, the high-frequency components of MS signals in the concentrated zone reduced again and then kept constant. The high-frequency components reduced again on December 31, 2014, which indicated a further increase in the source rupture scale. Coincidentally, obvious deformation occurred on January 2, 2014, as shown in Figure 7(b). Since January 6, 2014, the high-frequency components of MS signals increased sporadically, and the center frequency tends to rise. Correspondingly, the displacement curve of multipoint extensometer flattens out after January 9, 2014.

During the two deformation processes of surrounding rock mass as shown in Figure 8(b), the high-frequency components of MS signals all went through the evolutionary process of reducing first and then increasing. This indicates that the frequency variation of MS signals corresponds well with surrounding rock mass deformation. Therefore, the deformation process of surrounding rock mass can be revealed by analyzing the frequency variation characteristics of MS signals in the concentrated zone in the selected period.

The distribution range of the dominant frequency of MS signals in the concentrated zone during the two deformation stages in Figure 6 is equally divided into several bands. The percentage of MS signals components in different ranges is statistically analyzed. The accumulated variation regularities of the dominant frequency in the two deformation stages are obtained as shown in Figures 8 and 9, respectively. It can be found from Figure 8 that the low-frequency components of MS signals obviously increased between November 28 and 30, 2013, and the phenomenon of frequency transferring to low frequency occurred, compared to that between November 26 and 27, 2013. Correspondingly, an obvious deformation of surrounding rock mass appeared during the selected period as shown in Figure 8(b). The components of MS signals were further transferred to lower frequency band between December 1 and 4, 2013. However, the high-frequency components of MS signals increased between December 5 and 14, 2013. Coincidently, the surrounding rock mass deformation continued to develop on December 5, 2013, and ended on December 12, 2013, as shown in Figure 7(b). It can be found from Figure 9 that the low-frequency components of MS signals increased and the high-frequency components decreased, accompanied by the frequency changing from high frequency to low frequency. This phenomenon occurred during the period between December 27, 2013, and January 1, 2014, before an obvious deformation appeared. The MS signals with low frequency continued to occur during the period between January 2 and January 6, 2014. Then, the high-frequency components of MS signals increased between January 7 and 10, 2014. These phenomena can verify the results of Section 3.2 very well. Consequently, prior to obvious deformation of surrounding rock mass, the frequency of MS signals usually shift from high frequency to low frequency. The decrease of MS signal frequency in the selected concentrated zone can be thus regarded as a precursor for deformation and failure of surrounding rock mass.

5. Conclusions

(1)Compared with the conventional methods, the ST method can provide intuitive time-frequency characteristics without selecting the window function and analysis scale. Moreover, ST can self-adaptively adjust the analysis time width according to the signals and can accurately describe the temporal variation characteristics of nonstationary signals such as MS signals.(2)The correlation between source radius and frequency characteristics of MS signals is investigated based on numerous MS data from three large-scale hydropower projects. The correlation between fracture scale of rock mass and frequency of MS signals is also revealed. The results demonstrate that high-frequency components of the MS signal decrease, and the signals are characterized by low frequency when the rupture scale increases.(3)A case study is conducted based on the results of MS monitoring and multipoint extensometers results obtained from the underground cavern of the Houziyan project. The results demonstrate that the frequency variation of MS signals corresponds well with the deformation process of surrounding rock mass. The high-frequency component of MS signals decreases at first and then increases during the deformation processes of surrounding rock mass. Prior to obvious deformation of surrounding rock mass, the phenomenon of MS signals shifting from high frequency to low frequency usually occurs. Therefore, the decrease of MS signal frequency can be regarded as a precursor of macroscopic deformation failure of surrounding rock masses, which can provide some references for stability evaluation and macroscopic deformation failure forecasting of surrounding rock mass.

Data Availability

The authors declare that all the data supporting the conclusions of the present study can be obtained from the corresponding author.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

Nuwen Xu proposed the idea of correlation analysis between fracture scale and frequency of MS signals. Peiwei Xiao, Peng Jiang, and Biao Li processed and analyzed the abundant MS data. Peiwei Xiao, Bo Qian, and Nuwen Xu wrote the paper.

Acknowledgments

The authors are grateful for the financial support from the National Natural Science Foundation of China (no. 51679158) and the National Program on Key Basic Research Project (no. 2015CB057903).