{'fraction_non_alphanumeric': 0.02597151095236046, 'fraction_numerical': 0.011081464905108232, 'mean_word_length': 4.749010390895596, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Experimental Investigation on AE Precursor Information of Rockburst

This chapter investigated the AE precursor information of rockburst. Firstly, the characteristics of rock cracks initiating, propagating and linking was analyzed with the spatial-temporal distribution characteristics of AE, AE evolution pattern and the precursor feature of limestone failure were found. Secondly, the characteristic point of the relatively quiet period was defined as a characteristic precursor point of rock failure. The corresponding time was suggested. When AE energy release rate reached at the maximum, as the critical point for predicting rock failure. Thirdly, the ideal precursor parameters to predict rock failure were suggested after studying the variation regulation of AE energy parameters for rock failure. The spatial-temporal-energy evolution model of AE events on the fracture surface was established during rock failure, as well as the characteristics of fracture surface and the energy evolution.

with dissipation, Carpinteri and Lacidogna (2006) analyzed evolving relation between AE energy and dissipation energy under the concrete beam loading using the three point bending test.
In a word, many scholars conducted lots of AE experiments; most of their work was focused on the relationship between AE parameters with rock fracture and failure. However, in this chapter, the precursor of rock failure was utilized to predict rockburst using the distribution feature and evolution pattern of AE/MS in the deep mining.

Spatio-temporal Evolution Pattern of Rock Failure
Under the external loading, the virgin rock cracks were initiated, propagated, and linked. At the same time, AE events were released in the form of elastic stress waves. Therefore, with the spatial-temporal distribution characteristics of AE source 3D-location, the characteristics of rock cracks initiating, propagating and linking was analyzed, AE evolution regular and the precursor feature of limestone failure were revealed finally.

AE Experiment
In this experiment, 10 limestone specimens were divided into two groups. During the testing, the axial equivalent displacement loading method was used to load the specimens, and AE monitoring system was synchronously operated. The loading speed was controlled at 0.5 mm/min. In the AE monitoring testing, data collection was performed six channels, and the sampling frequency was set at 1 MHz. The threshold voltage value was set at 100 mv and the threshold amplitude value was set at 45 dB. AE signals would be filtered for removing noise signals, and a high signal/noise ratio could be achieved. During the test, six nano30 sensors were used to monitor AE signals and the pre-amplifier gain was set at 40 dB. The sensors were fixed on rock surface by adhesive band and petroleum jelly. In order to reduce the influence of the end effect, there should be 20 mm distances between the end of samples and the sensors. The layout of AE sensors were illustrated in Fig. 6.1.

Principle of AE Event Location
In order to locate the micro-cracks inside rock sample by AE monitoring, the sensors were arranged in a geometric pattern to form the sensor array on the surface of rock sample. The location of AE source was then determined based on AE signal parameters and location methods.
For the sudden type of AE signals, time difference location method and regional location method were mainly used to achieve the goal. Combined with the spatial location of sensor, P wave velocity and other related parameters, the former method was usually use as the measurement of parameters, such as the time difference among AE channel signals reaching the sensors to determine the location of AE source by the simplex algorithm or the Geiger algorithm. This location method was widely used in small rock samples and large components because of its higher accuracy, though the calculation was more complex than the latter one. Due to the simpler location operation and higher processing speed, the regional location method was generally applicable to situations where the location accuracy was low and the time difference could not be located.
The simple location algorithm and Geiger location algorithm (Miller et al. 1999) were two commonly used algorithms of AE location. The principle of these two algorithms was to invert the location of AE source by picking up the P-wave arrival time difference through different position sensors. Geiger algorithm was one of Gauss-Newton minimum fitting function applications, which was suitable for AE/ MS and earthquake events in small areas. In this chapter, Geiger algorithm was used to locate the AE source based on the time difference between the arrivals of P wave. In order to determine the real AE source coordinates, this algorithm was done by iterating over a hypothetical AE source. In each iteration process, a correction vector Dh was calculated by least squares method, and added it to the previous iteration result to obtain a new AE source location coordinate, then put the new coordinate parameters into the time-distance equation (6.1) to judge whether the new coordinate satisfies the equation. If the requirement was met, the new coordinate was the location of AE source; otherwise, it was need to continue iterating until the time distance equation was satisfied. In the AE location test, sensors were arranged in accordance with Fig. 9.3, and the location of AE events could be determined by measuring the relative time difference between each sensor picking up the P wave. The time-distance discriminant equation for the iteration results is as follows: where x; y; z are the initial assumed AE source coordinates; t is the time at which the initial AE event occurs; x i ; y i ; z i are the spatial coordinates of ith sensor, respectively; v p is the velocity of P wave; t i is the time that the ith sensor picks up the P wave.
The time at which the P wave propagates to the ith sensor is t 0;i , The first-order Taylor expansion of P wave arrival time can be calculated using the initial assumed AE source coordinates: where, where t c;i is the time that the P wave spread from the assumed AE point coordinates to the ith sensor. There were six sensors in this experiment for three-dimensional location. Therefore, six equations were listed in the form of a matrix as follow: @t 1 @x @t 1 @y @t 1 @z 1 @t 2 @x @t 2 @y @t 2 @z The modified vector Dh in Eq. (6.3) was solved by Gauss elimination method, which was Dh ¼ Dx; Dy; Dz; Dt ½ . After the modified vector Dh was calculated for each possible AE source coordinate matrix, the iteration was continued with ðh þ DhÞ as the new coordinates until the Eq. (6.1) was satisfied. At this time, the satisfied coordinates could be used as the real location.

Spatio-temporal Distribution Characteristics of Rock Failure
As the primary micro-cracks closed, expanded and linked, each sample generates lots of AE events in the entire loading process. The collecting system could be manipulated the AE date and obtained the visualized location images. Take only X2 sample date for example, as shown in Figs. 6.2 and 6.3. As shown in Fig. 6.2, when the stress was at low level (\0:2r c ), AE activity was active and gathered round middle part and top part of sample. This phenomenon could be explained as follows: in one hand, the original developmental cracks of samples were gradually closed under compression loading; on the other hand, a small quantity of AE events were occurred because of friction between samples and the pressure plate of machine. In the elastic deformation and stable growth crack stage (0:4À0:6r c ), AE counts increased steadily and accompanied with the development of microfracture. AE events dispersed gradually from the middle part to around part. When stress increased at the 0:8r c , AE events increased quickly and gathered round the main fracture, which presented obvious zonal distribution. The sample surface began showing the macrocrack, and AE event increased. In the post-failure stage, AE activities did not stop, but stress dipped quickly.
As the spatial distribution of AE events change with the loading time, we have analyzed the feature of spatio-temporal evolution in the different loading stage. Though the distribution was very heterogeneous, we still found two main phenomenons after summarizing all sample results, as shown in Figs. 6.2 and 6.3: (1) Intensive phenomenon: It was easy to be found that AE events occurred mainly in the middle part of X2 sample. However, when the stress increased to the 40, 60 and 80% of peak stress, AE events were always around the main plane of fracture, and only few events was in others parts. Therefore, it was an effective way to predict roughly the final position and attitude of main fracture plane using the spatial distribution of AE event.
(2) Zonation: AE events internal sample presented the obvious zonate shape, and the location of shape was basically consistent with the final macro-fracture. The main fracturepart from the zonation was shown.
6.3 Characteristics of Relatively Quiet Period for Rock Failure

Deformation and Failure Characteristics of Loaded Rock
Stress-strain curves of limestone samples were shown in Fig. 6.4, the deformation trend of a specimens were similar during the entire loading process. There are four basic stages of variation trend, including the original crack compaction, the stable crack growth, the unstable crack propagation and the rock failure stage (Wang et al. 2015a, b). In the chapter, only one sample was choosed to analyze and discuss the mechanical properties and AE characteristics of the loaded rock. As shown in Fig. 6.6, the sample No. 3 was appeared the original crack compaction stage after the beginning of uniaxial compression. With the stress increasing, the sample came into elastic deformation and the curve remained linear variation because of stable crack growth, it spended for a long time. Cracks initiated, extended and developed into the unstable propagation. In the plastic deformation stage, curve started to change non-linearly and rock volume also increased slightly. It was thought that many micro-fractures propagated and jointed one another which caused the volume increase of mesostructures. When the microcracks jointed and formed into macroscopic cracks, rock relative stability was destroyed, and stress reached the peak 102.89 MPa (Wang et al. 2015a, b). Stress didn't dipped sharply after the peak, and rock still had certain supporting ability. Therefore, there was secondary stress drop in the post-failure stage. Due to heterogeneity of rock, there was different destruction style, such as oblique shear failure and tensile failure. In the same way, the different mechanical properties of samples was because the size and the heterogeneity of rockmass. The analyzation of rock mechanical properties was important to understand the AE phenomenon (Wang et al. 2015a, b).

Relationship Between the Constitutive Characteristics and the Accumulative AE Counts
We could clearly find the relationships of sample No. 3 between stress-strain curve and the accumulative AE counts in Fig. 6.5. AE activity was constantly changed during the whole uniaxial loading (Wang et al. 2015a, b). However, there were still several changing characteristics in the accumulative AE counts. Before stress reached 20% of the peak stress, the original cracks were gradually closed and the friction occurred between samples and the pressure plate which generated a quantity of AE events occurred, the slope of this stage was larger than the slope of the elastic stage in the accumulative count curve. In the elastic deformation stage, cracks were in stable growth, and AE activity appeared obviously changing trend at the 30% of the peak stress, which we could regard as the Kaiser Effect point. From 30 to 60% of the peak stress, AE accumulative counts were linearly increased, which meanted AE activity stably developing. From Fig. 6.5, it was easy to discover that AE activity was unstable by comparison, which laterally reflected the heterogeneity of rock properties. When stress reached 60% of the peak stress, rock came into the plastic deformation stage and AE activities were enhanced. When stress was approximately at 75% of the peak stress, there were almost no AE activity and this phenomenon lasted 5 s. In other words, the sample deformation was in a relatively  Fig. 6.5. Then AE activity generated actively. When stress passed through the yield-point to peak stress, AE activities were rapidly enhanced, and the parameters of AE events were also increased, such as the amplitude, energy, frequency, duration, et al, which implied the magnitude of AE event stronger, and the size of fracture bigger. Therefore, in author's opinion, we could use the relatively quiet period as the precursor of rock failure to provide a predicting signal for rockmass instability. Similarly, we could use this opinion to analyze the micro-seismic events of seismicity to evaluate the rockmass stability in mining engineering (Wang et al. 2015a, b). After the peak stress, rock damage was in the post-failure, but AE activities did not stop, and intermittent emission constantly appeared. With stress declining, AE count increased periodically.

Relationship Between AE Rate and Constitutive Characteristics
Due to the samples inhomogeneous, mechanical and AE characteristics are nonlinear under uniaxial compression. Meanwhile, AE rate represents the intensity of AE activity and the developing speed of rock microcrack. Therefore, AE parameters of time domain can clearly reflect the real-time damage evolution, such as AE rate, AE energy rate and AE amplitude. AE duration and AE frequency reflect the properties of single AE event. In order to legibly understand the characteristics of AE and find the feature of AE relatively quiet period, we need to analyze AE parameters respectively. AE rate is a sensitive parameter about real time, which can reflect the frequency of AE activities and evaluate the damage status. The relationships between the constitutive characteristics and AE rate about the sample NO. 3 are shown in Fig. 6.6.
From the above, it was easy to discover that AE activity was always active during the whole uniaxial compression. As shown in Fig. 6.6, AE accumulative count appeared a small increase at 15% of the peak stress. We also found AE rate sharply increased to 800 at 25 s which was caused by inhomogeneous development of microcrack in the sample. In the initial deformation stage, AE rate increased firstly which implied the initial microcavities closed gradually. After 10 s, the rate became dipping trend and appeared a sharp increase around the time 25 s. Then AE activity maintained steady change in the elastic rock deformation, and the rate also appeared an obvious increase at the Kaiser Effect point. In the plastic deformation damage, the rate increased gradually and stress appeared the linear increase. In this stage, cracks propagation was unstable, and micro-cracks merged progressively. Before 70 s, the rate appeared a sharp increase, but stress curve didn't appear any dip trend. However, AE rate decreased sharply afterwards and there was a relatively quiet period appeared and lasted 5 s before the stress reached the peak stress. The loaded sample quickly passed through the quiet period before the arrival of the peak stress. AE rate increased rapidly when stress reached the peak. In the post-failure damage, AE rate became intermittent and several periods appeared with the process of stress decline. From the change trend of AE activity in the failure stage, we could easily find how strong the intensity of rock damage (Wang et al. 2015a, b).

Evolution Characteristics of AE Amplitude
in the Time Domain Figure 6.7 was a scatter diagram in which each scatter point represents a single AE event. AE amplitude was the maximum amplitude value of each AE wave which could be gained by Fourier Transform, which reflected the magnitude and intensity of AE activity. In view of the effectivity of AE extraction and the environment noise, we set the threshold at 45 dB. Therefore, the minimum amplitude of AE signals was 45 dB. As shown in Fig. 6.7, most amplitude of AE signals was below 60 dB, and few of AE signals were beyond 80 dB. As the amplitude reflected AE event energy directly, majority of AE events were had low-energy. On the contrary, only 0.7% of the total events were high energy. In other words, local slight fractures generated in the samples during the entire loading process. In the initial loading stage, most of low-energy AE events were generated, and some of amplitudes began distributing beyond 80 dB, which represented the completion of carities closing and the beginning of microcracks propagation. In the elastic deformation, the AE amplitudes were steady distributed in the time domain. But in the plastic deformation stage, the amplitudes presented the ascent distribution. Then in the unstable crack propagation stage, lots of high-amplitude AE signals appeared rapidly. When stress reached 75% of the peak stress, a relatively quiet period appeared as shown in Fig. 6.7. After the peak stress, rock strain began softening, and many high-amplitude and high-energy AE events were generated in the same time, which reflected the high energy was releasing with rock failure instability (Wang et al. 2015a, b).

Analyzation of Precursory Characteristics of Rock Failure
During the testing, a relatively quiet period of AE generation was found before rock failure during the test. To some extent, AE activities implied stress state of micro-fracture and reflected the damage evolution in the rock. But there must be some factors which caused the AE intermittent existed in the whole compress process. On the one hand, with uniaxial compressive stress increasing, strain energy of rock system transformed as deformation energy and dissipated energy. In the different deformation stages, energy percent of deformation and dissipation was relatively changing, which could be due to rock inhomogeneity, deformation energy could be also divided into volume energy and fracture energy, as the dissipated energy could be divided into the thermal energy, AE energy and radiant energy. The dissipated energy could be regarded as energy release with system energy import. Therefore, there was a balanced time domain in the process of energy import and energy release, which caused the propagation of micro-cracks, and energy release reach dynamic and temporary balance before the peak stress. This is one of reasons we can see the relatively quiet period in the loaded time domain.
On the other hand, high-damage intensity caused the sharp attenuation of AE propagation before the peak stress. There were many microstructural cracks existed in sample, which caused inhomogeneity of rock properties. Therefore, the propagation of AE waves was badly effected and the quiet period appeared. With stress increased, the state of balance was broken, and macro-fracture was formed causing the stress dip. The relatively quiet period might be mainly caused by inhomogeneity of samples (Wang et al. 2015a, b).
Therefore, we can use the AE quiet period as precursory information for prediction of rock failure and the practical engineering applications.

Variation Regulation of AE Energy
Parameters for Rock Failure

Experimental Setup and Procedure
In the experiment, 10 granite specimens were tested, approximately 50 mm in diameter and 100 mm in height. The loading machine and AE testing system were operated at the same time to obtain the mechanical behaviors and AE dates of specimens under uniaxial compression. The loading speed was 0.005 mm Â s −1 . In addition, AE sampling frequency was set at 1 MHz. Six AE sensors were pasted on the surface of sample to detect AE waves. The pre-amplifiers were set at 40 dB. Moreover, eliminating the influence of background noise, the threshold value of PCI-II acquisition system was set at 45 dB.

Energy Parameters Variations in Each Channel
As six sensors were installed on the sample surface, it was meant that there were six channels could detect the same AE waves. Due to the inhomogeneity of rock properties, the attenuation of AE signal might be different in local positions.
Therefore, one AE signal might be detected by the different position sensors, and some of sensors might not receive the signal, which caused AE energy parameter variations were different for each channel. Variation trend of AE energy was shown in Fig. 6.8, energy of 6 channels almost had the same variation trend. In the initial loading stage, AE energy increased firstly and then decreased subsequently. Energy of channel 1 and channel 2 was relatively insensitive, because the attenuation of AE energy was seriously when wave reached at the sensors. During the medium stage, AE energy kept plain evolution. After 70 s, AE energy increased sharply and reached at the peak, which lasted in the almost the entire failure process. With the destruction of rock supporting stability, AE energy also sharply dipped. The reasons for this phenomenon could be explained as follow. The original inner cavities were compacted with lots of energy releasing. While in the elastic stage, rock deformation and cracks propagation were generally stable. Thus, AE energy release was also stable. However, as rock in the plastic stage, micro-fractures were in unstable expansion. Then, microcracks jointed macrocrack which led to rock failure while high energy fiercely released (Wang et al. 2015a, b).

Relationship of Energy and Stress-Strain Curve on Time Domain
As shown in Fig. 6.9, there was a fine relation between AE energy release and constitutive characteristics of loaded rock. It was obvious to find that most of AE signal energy was 300 MV/s, and high energy event appeared in high stress stage, which directly reflected the stability of damage evolution. AE energy had different stages of relevant changing characteristics, describing as follow: the original crack compaction stage, the elastic deformation stage, the plastic deformation stage and the failure stage. During the first stage and the second stage, rock got the compression, but rock kept strong elastic property, so stress can be linearly developing.
In this stage, a lot of AE waves had been generated with the cavities compaction  Wang et al. (2015a, b) and cracks initiation. In the third stage, rock sample was in plastic deformation and cracks developing were relatively stable. Before the yield stress, there was relatively stable AE energy release which reflected the balance of energy transformation. After 65 s, stress curve bagan changing in the unstable slope, directly reflecting the strain unstability. Meanwhile, many high energy signals were generated and implied the intense cracks expansion (Wang et al. 2015a, b). Therefore, variation characteristics could be applied to replace the field stress monitoring as the evaluation factors of rockmass stability. As AE signal energy increased abruptly, it could be thought as the precursor for rock rupture.

Analysis of Damage Features
As shown in Fig. 6.9, high AE energy began appearing after 65 s, which meant high-intensity damage in the sample. Therefore, it was very important to analyze the energy evaluating characteristics of this period. In this chapter, we selected the summation energy, average energy and maximum energy to analyze as shown in Figs. 6.10, 6.11, and 6.12. Three kinds of energy are defined as follow: Summation of energy was the cumulative value of AE energy detected in all channels at the same time; the average of energy was the value as the total energy divided by the number of channels at a certain time; the maximum of energy was the maximum energy among each channel detected at a certain time. As the above, there should be pointed out that a channel could detect multiple energy signals in the same time, but each channel could not detect one event energy at the same time, which caused three values different at the same time point (Wang et al. 2015a, b).
As shown in Figs. 6.10, 6.11, and 6.12, three initial indexes were presented the same variation trend from 69 to 99 s, but the values jumped in the time domain. So we had fitted the initial curve plotted in Figs. 6.10, 6.11, and 6.12. From the fitted Fig. 6.9 Relationship between AE and stress with time. From Wang et al. (2015a, b) curves, we could find that the energy values reached the maximum before 90 s, indicating the maximum release of AE energy was before rock failure. Therefore, we estimated rock stable status by observing the variation of three indexes simultaneously. Compared with single energy index, we could more exactly judge rock energy release. Corresponding the fitted equations f(t), summation of energy was obtained based on the experimental data as presented in Eq. (6.4). From derivation, the time of extreme values in the parabolic curves are t 1 , t 2 and t 3 , which could be used as the critical times (Wang et al. 2015a, b). The fitted curve of energy summation could be described as follow: f ðtÞ ¼ À147:2t 2 þ 24,330t À 971,000 ð6:4Þ The time of the maximum value in energy summation curve could be obtained by the derivative of Eq. (6.5): where f 0 t ð Þ ¼ 0, so t 1 ¼ 82:64s. Based on the calculations above, the time of average energy, and the maximum energy could be obtained as t 2 = 82.89 s and t 3 = 81.57 s. Calculated the average of t 1 , t 2 and t 3 was 82.37 s. It was obvious to find that the times of three indexs almost reached at the same time. Meanwhile, it was thought that the time of most intense AE energy release occurred around 82.37 s and regarded the point of time as a critical feature to predict rock failure (Wang et al. 2015a, b).

Energy Release Rate of Rock Damage
As shown in Fig. 6.9, as lots of scattered AE signals released, AE energy was always in dynamic evolution during the entire loaded process. Thus, it was hard to exactly describe the real energy release and the actual damage level. Then the concept of energy release rate was proposed, which referred to AE energy released per unit time, and eliminated the instability of AE release. Thereby, it was significative to further study the mechanism of energy release rate and analyze its reasonability. Being the numerical signal, AE energy is a quantitative description of continuous AE signals. From electronics theory, AE energy release rate is proportional to mean square voltage of AE signals, which means that the voltage can replace AE energy to describe the damage variation. Mean square voltage to calculating AE energy can be formulated as follows: Furthermore, the energy release rate can be calculated as: and the mean square voltage is calculated as where V(t) is voltage of AE signals; DT is mean time interval. Energy release rate is expressed by the mean square voltage, as shown in Eq. (6.8).
Thus, the damage stability of rock can be judged by analyzing the energy release rate. Meanwhile, there is a positive correlation between mean square voltage and energy release rate. Therefore, the evolution process can be directly described by the variation of voltage to provide information or predicting of rock failure (Wang et al. 2015a, b).

Evolution Patterns of Spatial-Temporal-Energy on Rock Fracture Surface
Macroscopic failure of rock materials is the result of initiation, propagation and coalescence of rock cracks with a rapid release of transient elastic wave (AE). Much work of rock mechanics behavior and damage evolution were studied, and several valuable experimental results were obtained (Vidya Sagar et al. 2013). However, it is not comprehensive to study the evolution of rock failure process only from the features and roughness of fracture surface, rock failure process is a change status by energy release. Therefore, some work about energy point of rock failure were studied and several valuable experimental results were obtained (Peng et al. 2015). Wasantha et al. (2014) investigated the pattern of energy release during deformation of dry sandstone and saturated sandstone with different bedding directions under uniaxial compression. Carpinteri et al. (2006) studied the energy and dissipation energy of AE in the loading process of three point bending test of concrete beams in detail. Meanwhile, as we know, fracture surface is the main carrier of deformation and failure process; there are few researches which consider the relationship between fracture surface and energy of AE. Therefore, in this chapter, failure process of rock was studied through analyzing the changing pattern of time, space and energy of AE event on the fracture surface. In this chapter, the spatial-temporal-energy evolution model of AE events on the fracture surface was established during rock failure, and the characteristics of fracture surface and the energy evolution were obtained during failure process, it provided us with scientific evidence for predicting rockburst.

Experimental Setup and Procedure
In this chapter, limestone samples were taken from the roof of coal seam at 1212 m depth in Xiezhuang Mine in Shandong province. These samples were processed as standard cylindrical rock samples of 50 mm in diameter and 100 mm in length. The experiment was adopted the axial displacement loading method, and the loaded rate was controlled at 0.005 mm/s. The loading machine and AE motoring system was operated synchronously.
The preamplifiers of AE were set at 40 dB before the samples were loaded. However, reducing the influence of environment noise, the threshold value of PCI was set at 45 dB, and sampling rate was set at 1 MHz. AE testing system with 6 sensors for realizing real-time 3-D positioning of limestone samples during the entire deformation and fracture process. To reduce the influence of end-effects, 6 AE sensors were fitted on the limestone specimen 20 mm away from top and bottom of sample. Petroleum jelly and the adhesive band were used to fix the sensors on specimen surface. After the experiment, the relevant parameters of each event were analyzed during the entire process by Geiger algorithm, such as the AE rate, AE energy, AE position, and AE frequencies.

Relationship Between Spatial Distribution of AE Events and Rock Fracture
To study the evolution of spatial-temporal-energy evolution in the course of limestone failure, a series of AE tests were carried out for limestone failure under uniaxial compression. Stress strain data and AE data were obtained. Geiger algorithm was used to analyze AE events during loading process, and obtain the parameters of each event (time, position, energy). The corresponding parameters of AE events during the loading process of X3 rock sample were illustrated in Table 6.1.
A large number of AE events were obtained immediately, due to the internal deformation and cracks propagation. AE events reached several thousand or even tens of thousands during a loading process, if all the events were analyzed, then the workload was too great. Many small energy events could be ignored in the rock failure process. As a dominant reason for the ultimate failure of rock, which could only be large energy events. So, a large number of AE events were screened for energy, and finally the appropriate threshold of energy was determined. AE events over threshold were defined as a strong AE event. Comparison of AE spatial distribution and actual rupture propagation before and after the screening of X3 rock sample was demonstrated in Fig. 6.13.
Comparison of AE events spatial distribution and actual fracture propagation before and after the screening of X3 rock sample was shown in Fig. 6.13. As shown in Fig. 6.13b, AE events could be better reflected the evolution characteristics of internal rupture. There were many AE events in Fig. 6.13a, which leaded to the chaotic distribution of AE events in rock sample, and it was difficult to determine the specific location of main fracture surface. As shown in Fig. 6.13b, compared with the actual fracture propagation, AE events were mainly distributed near the main fracture surface of rock sample, which could clearly reflect the initiation, propagation and coalescence process of main fracture surface. Therefore, it is reasonable to study the spatial-temporal-energy evolution characteristics of AE in the loading process using energy screening method.

Spatial-Temporal Evolution of Strong AE Events on Rock Fracture
It was found that AE spatial distribution contained full information about internal deformation characters and cracks evolution of rock failure. Therefore, it provided a basis for studying the propagation process and spatial morphology of rock during loading process. AE localization of sample X3 in different loading period under uniaxial compression was shown in Figs. 6.14 and 6.15, where the color represented the magnitude of events energy value.
Analyzing the spatial-temporal-energy evolution characteristics of strong AE events, the spatial-temporal-energy evolution model of rock fracture surface was established. In this chapter, using the energy screening methods, it was showed that   6.5 Evolution Patterns of Spatial-Temporal-Energy on Rock Fracture Surface most of strong AE events were distributed on rock fracture surface. As a matter of fact, the weak AE events had little effect on rock failure, and the occurrence of strong AE events corresponded to damage deformation. Therefore, it could be better analyzed the spatial position, propagation direction and spatial shape of cracks propagation in rock samples, and it could be more accurately determine the location of fracture surface, so as to provide predictions for the final rupture of rock condition.

Spatial-Temporal-Energy Evolution Model of Strong AE Events in Limestone Fracture
As shown in Fig. 6. 16, it was the relationship curves between cumulative energy of strong AE events and time in the fracture surface of rock samples during the loading process. Take X3 rock sample for example, when stress was less than 20% of the peak stress, in the fracture compaction stage, more low energy events were monitored in the lower part of rock sample, and the middle energy events were relatively large in the lower part, which may be the internal joint fracture. The lower part was more developed, under the action of external load joints began to close, resulting in AE events, and activities were more intense. With the stress increased, rock sample entered the cracks initiation stage, and a small amount of low energy AE events were observed at the top and bottom of rock sample. The fracture surface did not form at this stage. When stress reached at 40-60% of the peak stress, rock specimen entered cracks in the stable development stage, and AE events were significantly increased in the lower part of rock sample, event energy was larger, forming the strip shape with the pre-generated event. When stress reached at 80% of the peak stress, rock sample entered the unstable development stage of cracks, AE events increased rapidly and energy increased significantly. AE events on the strip began to be concentrated, and then internal cracks began to pass through, forming macro-cracks, and fracture surface was formed. In the later of this stage, number of AE events and the accumulated energy hardly increased, and there was a brief silence zone, namely, AE "quiet period", which could be used as a predicting point of rock destabilization failure (Wang 2015b). When stress exceeded the peak stress and entered the failure stage, AE event was obviously reduced compared to the previous stage. The macro fracture connection and coalescence of rock sample form the main fracture surface, and the rock sample was failure. From the above analysis, it could be seen found that the formation of fracture surface mainly concentrated on the plastic stage and the failure stage. Therefore, the pattern of energy evolution pattern of rock fracture surface should begin with the plastic stage, which could be better revealed the evolution pattern of propagation and coalescence of internal cracks in rock samples.
The cumulative energy of plastic section of each rock sample, the abscissa was as the time, and the ordinate was selected as the cumulative energy. The Gauss Amp function was used to fit, and the fitting function was y ¼ y 0 þ Ae À xÀxc ð Þ 2 2x 2 , the dependent variable y was used to represent the cumulative AE energy, x c was used to represent the critical moment, y 0 was used to represent cumulative AE energy in the plastic stage, A was used to represent the amplitude of cumulative AE energy, x was the relevant parameters, and it could be satisfied the equation 2x ¼ FWHM=sqrt ln 4 ð Þ ð Þ, which could calculate the value of x. So the original formula changes to E ¼ E 0 þ Ae À tÀtc ð Þ 2 2x 2 . The obtained AE data were analyzed in the experiment, and the cumulative energy-time curve of strong event in fracture surface was plotted. Then the above functions were used to fit, the fitting results were shown in Fig. 6.16, and the related parameters were shown in Table 6.2.  As shown in Fig. 6.17, besides the rock samples X3 and X5 (the correlation coefficient of X3 is 0.93 and X5 is 0.91), the fitting correlation coefficients of other 3 rock samples were all above 0.97. The results showed that Gauss Amp function fitted well, and the parameters E 0 , A, t c and x could be obtained by fitting the results.
During the loading process, the input energy was mainly converted into two parts: one part was dissipated in the form of damage energy and plastic deformation energy of rock specimen. The other part was accumulated inside the rock sample in the form of elastic strain energy and could be released during unloading. The ratio of elastic energy and dissipated energy to input energy had a great influence on rock deformation and failure modes. As a part of dissipation energy, AE energy might be related to the external input energy. The area between stress-strain curve and transverse axis was the strain energy density. The related parameters of samples energy were shown in Table 6.3.
As shown in Table 6.3, vary of strain energy density and cumulative AE energy was showed at the plastic stage. However, when strain energy density was large,  cumulative AE energy corresponding was correspondingly large. Its proportion in a certain range ranged from 0.5 to 0.7%.
In conclusion, it could be concluded that cumulative AE energy increased with time during the failure process. Especially, after the plastic stage, cumulative AE energy of strong events in the fracture surface conformed to Gauss Amp function distribution. Therefore, the spatial-temporal-energy evolution model of limestone on the fracture surface was proposed: where k is a material constant, and its value is between 0.5 and 0.7%; t c is the critical moment; t is the loading time, and it ranges from the initial stage of plastic stage to the peak stress; e c is strain corresponding to peak stress; x is the coefficient in Gauss Amp function; A is the increment of cumulative AE energy from the beginning of plastic phase to the peak stress; E is cumulative AE energy of strong events.
Energy of AE events directly characterized the intensity of internal mesoscopic damage, and the weak event had little effect on rock failure, and the ultimate failure of rock was dominated by strong event. On the basis of predecessors, energy was taken into account in the spatial-temporal evolution analysis of AE events, so as to analyze the spatial-temporal-energy evolution of AE events on the fracture surface. The proposed spatial-temporal-energy evolution model combined the spatial-temporal events and energy of fracture surface in the process of limestone rupture and failure, and it described the energy evolution pattern in the formation of fracture surface under the uniaxial compression. The model was used to describe the energy change process of fracture surface systematically, and it provided a theoretical basis for predicting rockburst. However, this model is only applicable to the uniaxial compression test of limestone, and further study is needed for other rock and other stress conditions.

Spatial Fractal Dimension Evolution Model of AE Events on Fracture
In 1973, mathematician Mandelbrot first proposed the idea of using fractional and fractal geometry to describe irregular objects and chaotic phenomena in nature (Mandelbrot 1983). Analyzing AE signal produced in the propagation of micro-cracks, it was found that AE parameters followed not only fractal characteristics in time, but also fractal characteristic in spatial distribution. The balls covering method or the circle covering method after projection were widely used for three-dimensional distribution of AE events. For the cylindrical rock specimens, column covering method was more appropriate. In this method, the spatial distribution characteristics of AE events were apt to be described exactly. Due to the shape of rock specimens, it was very convenient for the method to calculate the fractal dimension of spatial AE event distribution according to the coordinates of AE events. Therefore, the method could be used to determine the critical value of rock failure, predict the rockbursts, and evaluate the stability of rockmass. When AE events were covered by a cylinder with radius r and height h, the relationship between AE event number and the radius of cylinder was: where C1 represents height-radius ratio. Equation (6.10) could be further expressed as: where r is a given radius, M(r) is a value which can be obtained through Eq. (6.11).
In a double logarithmic coordinate system, point series (log r, log M(r)) corresponding to given radii can be portrayed. The spatial AE event distribution has fractal characteristic in the range of a given scale if those points approximately distribute along a straight line. The slope of line is the exact fractal dimension of spatial AE event distribution. Based on the above-mentioned column covering method, using the different radius r of the cylinder to cover the strong AE events, the relationship between number of strong events M(r) and radius was obtained. In a double logarithmic coordinate system, point series (log r, log M(r)) corresponding to given radii could be portrayed. Then the linear point of data points was fitted, and the obtained slope was the fractal dimension of spatial distribution of large events in rock sample under the stress level. The curve between the spatial fractal dimension and the stress level of strong AE events of limestone rock samples was shown in Fig. 6.18. When stress was at a low level (below 0.2r c ), the spatial fractal dimension, close to 3, was relatively large indicated that AE events in this stage were more evenly distributed in the interior space of rock samples. With the increase of stress, rock sample entered the cracks initiation stage and cracks stabilization expansion stage (0.2-0.7r c ). Because of less AE events, D value of this stage decreased slowly. With the increase of stress, rock specimen entered the unsteady expansion stage (0.7-0.9r c ), the internal damage increased and the fractal dimension decreased continuously. When stress value reached about 80% of the peak stress, D value decreased rapidly, indicated that the internal events began to accumulate near the macroscopic main fracture surface. When stress reached the peak stress, the macro-fracture surface was formed gradually, and D value dropped to the minimum value, which was close to 2.2.
By analyzing the relationship curve between the spatial fractal dimension and the stress level of AE event, the spatial fractal dimension decreased linearly with the increase of stress level. Therefore, the least squares method was used to fit the curve, and the fitted curve was shown in Fig. 6 where D is the spatial fractal dimension, and its value is between 2 and 3; r c is the peak stress; r is the stress value during loading process.

Conclusion
Rock failure was accompanied by changing of AE parameters as well. In view of this issue, this section was built on the laboratory test, explored and quantized the feature of failure precursory information, and obtained the following conclusions: (1) Though the distribution was very heterogeneous, through analyzed the feature of spatio-temporal evolution in the different loading stage, we still found two main phenomenons after summarizing all sample results: intensive phenomenon and zonation. (2) The characteristic point of the relatively quiet period was defined as a characteristic precursor point of rock failure. The ratio of stress at the relatively quiet period was from 88 to 98% of the peak stress in the unstable propagation stage.
(3) In the plastic stage, AE energy rapidly increased and remained at high values indicating intense energy release before rock failure. In the failure stage, AE energy sharply decreased and rock instability occurred. The corresponding time, when AE energy release rate reached at the maximum, can be suggested as the critical point for predicting rock failure. Meanwhile, AE energy release rate could be considered as an important rock failure index, which might provide physical basis for predicting rockburst.
(4) Strong AE events could fully reflect the evolution of internal damage and deformation, and were concentrated in the vicinity of main fracture surface. Through the spatial distribution analysis of strong AE events, the position of rupture surface could be determined more accurately. (5) The range of fractal dimension of strong AE event was from 2 to 3, the value of the initial stage was close to 3, and the value of final failure stage was close to 2. The spatial fractal dimension showed a sharp decrease near the yield point (about 0.8r c ), indicating that rock was about to lose stability and could be used as a precursory feature of rockburst. Moreover, there was a linear relationship between fractal dimension and stress level. (6) The fracture surface was mainly produced in the plastic stage and the failure stage under uniaxial compression. Cumulative AE energy offracture surface conformed to Gauss Amp function distribution. Therefore, the spatial-temporal-energy evolution model in the fracture surface for strong AE event was established.