Research on Coal Acoustic Emission Characteristics and Damage Evolution During Cyclic Loading

The acoustic emission, energy, and damage evolution of coal samples for three kinds of uniaxial cyclic loading and unloading are deeply analyzed in this study. The evolution of total absorption energy, elastic strain energy, and dissipated energy of coal samples is related to the stress path, and the increasing amplitudes cycle loading has an obvious damage effect on coal samples. During the loading stage, the acoustic emission phenomenon is most active when loading is increasing and the Felicity and post-Kaiser phenomena appear. The acoustic emission phenomenon during constant loading does not obviously change, but rather becomes active with the increase of the equivalent load. The damage to the coal sample shows nonlinear change increasing loading and unloading and shows linear change for other stress paths. Compared to waveforms with stepwise increasing amplitudes cyclic loading, the failure process of the coal sample is more closely related to the size of the external load, which indicates that reasonable hydraulic design is beneficial to the stability of the confining pressure in the chamber of an underground pumped storage hydropower station.


INTRODUCTION
In the current energy transition from fossil energy to environmental energy, the proportion of energy from new sources, such as solar, wind, and hydropower, continues to hit new highs. However, these new types of energy generation are greatly affected by the weather, and so cannot provide energy consistently. Energy storage technology is one solution to the instability of these new energy sources (Bilgili et al., 2015). At present, energy storage technologies are mainly divided into chemical battery storage, compressed air energy storage (CAES), and underground pumped storage hydropower stations (UPSHs) (Rahman et al., 2020;Fan et al., 2020a). These technologies are applicable to different times and spaces. CAES requires strict geological conditions to meet the requirements of impermeability for high pressure air, which is why it is mainly used in rock salt mine caves (Fan et al., 2019). In contrast, a UPSH is more flexible in terms of technology and has lower environmental requirements (Gaudard and Madani 2019). China in particular has long relied on coal as its main energy source, and it has many abandoned mines, some of which are located near the country's main wind and solar power sources. Using these abandoned mines to build UPSHs can not only solve the problem of abandoned mines destroying the environment, but also contribute to the sustainable development of local economies (Menendez et al., 2019a;Menendez et al., 2019b). A UPSH consists of two reservoirs. The upper reservoir is located on the surface of the ground, which is an industrial square (an area protecting the ground factories of coal mines) for coal mining, and the lower reservoir is located underground. The lower reservoirs are usually built 1) to excavate and secure additional caverns, 2) to make use of existing roadways, or 3) to excavate new roadways (Madlener and Specht 2013). During the operation of a coal mine pumped storage power station, the discharge and impounded processes of water are repeated, which causes periodic loading and unloading of the surrounding rock inside the power station, leading to a change in the mechanical properties of the mine's surrounding rock. Determining how to maintain the stability of the rock surrounding the underground reservoir during the operation of an energy storage power station is an urgent problem that needs to be solved.
In order to solve the problem of surrounding rock stability, scholars have performed a great amount of research. Menéndez et al. (2019c) studied the influence of the pumping and storage process on the stability of the rock surrounding a UPSH in northern Spain. The cyclic loading and unloading pressures generated by the confining pressure of the UPSH chamber had an irreversible impact on the stability of the chamber's surrounding rock. Fu et al. (2020) studied the characteristic stress and mechanical properties of marble cracks for three different cyclic loading stress paths, and found that the crack initiation stress of the marble first increased, and then decreased, with the increasing number of cycles. Based on these test results, Fu put forward the CVFE constitutive model and applied it to Realistic Failure Process Analysis software (RFPA2D), and then studied the stability of the surrounding rock in the chamber. From their research, it was apparent that the load types of the surrounding rock of the CAES or UPSH were mainly cyclic. These cyclic loads could be classified into four classes according to the shape of the load-time curve (Heap et al., 2009;Jia et al., 2018;Liu et al., 2018;Peng et al., 2020): constant amplitude cyclic loading, incremental cyclic loading, variable amplitude cyclic loading, and random cyclic loading.
Using different circulation paths, scholars have carried out laboratory experiments on the mechanical properties of different lithological rocks; usually granite, marble, sandstone, rock salt, and coal (Wang et al., 2021a;. Tang et al. (2020) used acoustic emission monitoring technology to study the crack development law of the post-peak of granite under a cyclic load and found that the strain concentration zone expanded centrally along both sides of the crack, the crack opening degree increased with the increasing number of cyclic loads, and the acoustic emission B value increased first and then decreased. Liu et al. (2021) studied the acoustic emission characteristics of weakly cemented rock mass under different triaxial stress paths are studied. Liu et al. (2012) studied the influence of cyclic frequency on the damage degree of sandstone, and found that with an increase in frequency, the initial stiffness of the sandstone was higher and the number of cycles it withstood before failing was also higher. Rocks with different lithologies have different brittlenesses. Nejati and Ghazvinian (2014) studied the influence of rock brittleness on fatigue characteristics and found that in the crack damage stage, microcrack density increased with increasing rock brittleness while fatigue life decreased. Salt rock is a good gas storage medium with low porosity, high ductility, stable chemical properties, and strong selfhealing abilities. With this in mind, Roberts et al. (2015) studied the law of the cyclic loading creep and the conventional static loading of salt rock, and found that there was no significant difference in the creep characteristics of salt rock for different stress paths. Wang et al. (2021b) studied the creep characteristics of salt rock for five different long-term cycles and found that the damage process of salt rock could be divided into the attenuation, stability, and acceleration stages, and that the larger the stress amplitude was, the larger the cumulative irreversible deformation was.
However, the surrounding rock of an abandoned coal mine energy storage power station is usually dominated by coal, which is a soft rock and has different mechanical properties than hard rock . Zhong et al. (2019) studied the crack development of coal samples containing saturated water for different cyclic frequencies, loading rates, and amplitudes, and they found that pore water destroyed the original cracks and lead to tensile damage. In terms of water permeability, Yang et al. (2018) studied the change of coal permeability for a cyclic loading stress path and found that with the increasing number of cycles, coal permeability declined irreversibly, which was caused by irreversible deformation. However, there are still few research studies on the damage and fracture development of coal during a cyclic load. In this research, acoustic emission monitoring is used to study the change of the fracture development damage of coal under this type of load. The surrounding rock of the UPSHs is in different stress environment, the research results of this paper will be beneficial to UPSHs to choose a reasonable working condition, effectively reduce the destruction of the chamber surrounding rock, reduce the operation cost of the pumped-storage power station.

Specimen Preparation
The coal block gathered in the field is black and bright with parallel beddings, good bonding, and strong brittleness, as shown in Figure 1A. According to the requirements of the International Rock Mechanics Test Regulations, coal samples are cut perpendicular to the coal bedding to produce 20 samples containing horizontal bedding, with a size of 100 mm × 100 mm×100 mm, as shown in Figures 1B,C. Previous research results show that the rock strength and the elastic modulus decrease exponentially with the increase of the scale. When the cross-sectional area of the sample is larger than 50 mm × 50 mm, the uniaxial compressive strength of the cylindrical rock samples with a slenderness ratio between 1 and 3 changes only slightly, regardless of whether the cross-section is round or square. Ultrasonic testing equipment is used to test the uniformity of the rock standard samples. After the test, the p-wave velocity of coal samples is mostly near 1370 m/s, and 22 usable coal samples are screened out. The samples with large deviations are removed according to the test results.

Test Instrument
In this uniaxial compression test, a CSS-44300 servo electronic testing machine purchased by the State Key Laboratory of Coal Resources and Safe Mining of the China University of Mining and Technology is used, as shown in Figure 2A. This testing machine can synchronously collect the loading time, axial load, and axial deformation of a test block.
In this research, a PCI-2 acoustic emission monitoring system from the American Acoustical Company is used to record the acoustic emission activities of rock during the loading process. The system is composed of an acoustic emission host, preamplifier, and acoustic emission probe, as shown in Figure 2C. The acoustic emission preamplifier can provide gains of 20 dB, 40 dB, and 60 dB in three different gears. The type of acoustic emission probe used is Beijing Soft Lead RS-54A, with a diameter of 8 mm and a monitoring frequency range of 100-900 kHz, as shown in Figure 2D. In order to ensure the  Frontiers in Earth Science | www.frontiersin.org January 2022 | Volume 10 | Article 818452 3 accuracy of sampling in this research, acoustic emission probes are arranged according to the linear positioning method, and the lead breaking method is used to determine that the lead breaking amplitude monitored by the two probes exceeds 95 dB, which meets the requirements of the acoustic emission monitoring test. In each test, the ambient sound is collected and analyzed, the preamplifier gain is set to 40 dB, and the sampling frequency is set to 3 MHz.

Test Scheme
There are three different cyclic loading modes designed for the uniaxial compression test of coal samples, as shown in Figure 3, which are increasing cyclic loading amplitude, constant amplitude cyclic loading, and graded cyclic loading. Firstly, in order to determine the stress threshold of uniaxial cyclic loading, five rock samples are selected for the static uniaxial compression tests to obtain the uniaxial compressive strength and the Young's modulus of the coal samples. Then, the following three cyclic loading tests are conducted: 1) A uniaxial increment cyclic loading amplitude test is conducted. The load control loading method is adopted for loading. The loading rate is controlled at 0.05 MPa/s (0.5 kN/ s), and the cyclic stress gradient is set to 3 MPa (30 kN). At the end of six cycles, the specimen is compressed and destroyed. 2) A uniaxial constant amplitude cyclic loading test is conducted using the load control loading method for loading and unloading, and the loading and unloading rate is controlled at 0.05 MPa/s (0.5 kN/s). The first load is 90% of the uniaxial strength of the sample, then the sample is unloaded to 0.5 MPa, and then the sample is cyclically loaded and unloaded six times in this interval. Finally, the sample is compressed and destroyed. 3) Cyclic loading tests are carried out in stages. The loading rate is controlled at 0.05 MPa/s (0.5 kN/s). In the first stage, the loading is conducted at 0.05 MPa, and then the cyclic loading is carried out six times within the range of 0.05 MPa-30% of the uniaxial strength. Then, in the second stage, the cyclic loading is carried out six times within the range of 0.05 MPa-60% of the uniaxial strength. In the final stage, the sample is compressed and destroyed.

ANALYSIS OF THE ACOUSTIC EMISSION CHARACTERISTICS OF THE COAL SAMPLES FOR CYCLIC LOADING AND UNLOADING Test Results and Analysis of the Mechanical Properties of the Coal Samples
The mechanical properties of the coal samples are shown in Table 1, and the stress-strain curve with uniaxial compression is shown in Figure 4. It can be seen from the figure that the Frontiers in Earth Science | www.frontiersin.org January 2022 | Volume 10 | Article 818452 maximum uniaxial compressive strength of the coal sample is the same as the strain at failure, and the physical properties of the four blocks are consistent, so these properties can be used for subsequent tests. In the process of loading, the coal sample first appears to have a splitting crack. Then the sample is destroyed in the splitting and shear mixed fracture mode, and a large number of small splashes appear when the destruction occurs.

Coupling Analysis of the Mechanical Properties and Acoustic Emission of Coal Samples During Cyclic Loading and Unloading
The deformation and the failure of rocks are always accompanied by the release of strain energy in the form of elastic waves, which is called the acoustic emission phenomenon (Liu et al., 2020a;Liu et al., 2020b;Jiang et al., 2016). The AE of a rock mass with the action of loading stress can be divided into two categories. One category is the acoustic emission caused by severe elastic vibration due to the change of the elastic strain of the rock, and the other category is the acoustic emission caused by the sliding friction between microcracks (Chmel and Shcherbakov 2013). Four coal samples are taken for each cyclic loading path. For better research, only representative test results are selected for analysis. In order to comprehensively observe and analyze the relationship between the stress on the coal samples and the acoustic emission, these data are placed on a graph. As shown in Figure 5, there are three curves and one group of scatter points in the figure, which represent the curves of the cumulative counts and the times of the stress, strain, and acoustic emission events of the sample. A scatter plot of the acoustic emission ring count rate and the time is created.
It can be seen from Figure 5 that the strain under uniaxial compression is divided into four stages: the compaction stage, the elastic deformation stage, the plastic deformation stage, and the post-peak stage. Additionally, during the loading process of sample 1-1, the acoustic emission process is divided into three stages. The first stage is the acoustic emission contingency stage in which only sporadic acoustic emission signals appear and the acoustic emission count rate and cumulative count values are very low because the sample is in the compaction stage, and there are basically no new cracks in the interior. The second stage is the AE development stage, in which the AE count rate increases gradually but still remains at a low level. The cumulative count of AE events increases gradually and presents a step shape, indicating that new cracks begin to appear inside the test and form cracks with the expansion and convergence of new cracks. The final stage is the peak stage of acoustic emission, and the peak stage occurs in a sample near destruction. For this stage, the acoustic emission event count rate increases by an order of magnitude, and the number of accumulated acoustic emission events dramatically increases. This is due to the fast expansion of cracks and broken blocks.
For uniaxial cyclic loading and unloading, the relationship between the acoustic emission ringing count rate, the cumulative count of the acoustic emission events, the loading stress, the  sample strain, and the monitoring time of the rock specimens is shown in Figures 7-9. In order to facilitate analysis, the AE ring count rate in the figure is enlarged.
In the process of increasing cycle loading and unloading, an obvious Kaiser effect appears; that is, when the stress value reaches the peak value of the previous loading, the bearing rock presents an obvious acoustic emission phenomenon (the purple dotted line in Figure 6. When the stress value does not reach the peak value of the last loading stress, acoustic emission also occurs, which is called the Felicity phenomenon. In each cycle, with the increasing load, the phenomenon of weak acoustic emission first appears, and then the acoustic emission frequency increases rapidly. When the stress value exceeds the maximum value of the last cycle, the acoustic emission time increases rapidly and reaches the maximum value at the stress peak of this cycle. At the unloading stage of each cycle, a large number of acoustic emission phenomena appear at the initial unloading stage and then disappear with the decrease of the stress value to a certain extent. However, as the stress value at the unloading point of each cycle continues to increase, the acoustic emission phenomenon becomes active again in the unloading process, and the load of this phenomenon is close to the peak value of the last cycle, which is the post-Kaiser phenomenon. The appearance of the Kaiser effect and the Felicity effect indicates that in a certain unloading process, the rock not only recovers its elastic deformation, but also changes its plastic deformation Kang et al., 2021).
With the process of cyclic loading and unloading, the maximum emission count rate of each cycle increases from 400 times to about 1,250 times before failure. The maximum emission rate in each cycle is greater than that in the previous cycle; specifically, at the moment before failure, the maximum emission count rate reaches 1,500. This shows that with the increase of the number of cycles, the peak value of the cyclic load increases and the damage degree of the rock intensifies. As can be seen from the cumulative count curve, the cumulative AE count increases significantly by three times over a span of time, and it then reaches the maximum value near the failure of the sample. The significant increases corresponded to 9.5, 16.47, and 17.28 MPa, respectively. This shows that these are the stages at which a large number of fissures are produced or at which the fissures are rapidly expanded and transfixed.
In the loading process of constant amplitude cyclic loading and unloading, as shown in Figure 7, the AE phenomenon appears periodically, and the cumulative amount of sound emission gradually increases, increasing especially rapidly in the first cyclic loading stage and the final failure stage. In each cyclic loading process, the AE counting rate value increases gradually and reaches the extreme value at the loading stress maximum. The maximum AE counting rate in the first cycle is 94 times, and the maximum AE counting rate in the following cycles is between 50 and 70 times. At each cycle unloading stage, the AE activity weakens rapidly, and when the stress reaches the lowest point, the AE phenomenon basically disappears. In the constant amplitude loading and unloading cycle, the maximum stress value is close to the uniaxial compressive strength of rock, so an obvious acoustic emission phenomenon appears in each cycle, the maximum emission count rate fluctuates within a certain range, and the acoustic emission activity does not decline significantly. It can be seen that for the loading close to the maximum uniaxial compressive strength of the sample, new cracks are constantly generated and developed inside the coal sample. The coal sample is in a state of continuous damage, and each time the coal sample is in the active stage of acoustic emission, this indicates that the coal sample is about to experience yield failure.
For the condition of piecewise cyclic loading and unloading, as shown in Figure 8, the AE phenomenon appears periodically, and the cumulative number of the AE increases gradually, and increases especially rapidly for the increasing amplitude and the final failure stage. In the process of cyclic loading and unloading in the first stage, the maximum AE count rate in each cycle reaches the maximum of 208 times in the first cycle, 66 times and 74 times in the following two cycles, and then stabilizes within the range of 21-34 times, presenting a trend of fluctuation and overall decline. With the process of cyclic loading and unloading, the AE activity of the coal samples decreases gradually. In the process of cyclic loading and unloading in the second stage, the maximum AE count rate reaches a maximum value of 234 times in the first cycle and then fluctuates. Although the trend of the acoustic emission count rate of the coal samples is similar to that in the first stage, the acoustic emission activity in each cycle is FIGURE 6 | Relationship between axial stress, axial strain, and acoustic emission for cyclic loading with a uniaxial increment. higher than that in the first stage. This indicates that the upper limit of the load has a significant impact on the acoustic emission of the coal samples. The greater the load is, the more active the acoustic emission is. When the load in the second stage exceeds that in the first stage, the AE events increase rapidly, indicating that the Kaiser effect of the acoustic emission activity of the coal samples is an inherent attribute of the coal samples. The AE activity increases sharply when the coal sample is near failure, and then decreases sharply and disappears after failure.
To summarize, whether in uniaxial loading or in each cycle of loading and unloading, the acoustic emission phenomenon of the coal samples exists in three stages: the calm stage, the increasing stage, and the active stage. In different types of cyclic loading and unloading, the cumulative value of the AE keeps increasing until the sample is destroyed, and the AE events are mainly caused by the propagation of cracks, indicating that new cracks are constantly generated and developed inside the coal sample, and the damage degree of the coal sample keeps deepening. In different unloading stages of the cycle, acoustic emission phenomena also occur, indicating that the coal samples also have damage changes in the unloading process, but the degree is relatively small. Compared with cyclic loading and unloading at a low stress level, the high stress cycle close to uniaxial compressive strength continuously promotes the formation and development of cracks in the coal samples, so the acoustic emission is always in an active stage.

ANALYSIS OF ENERGY AND DAMAGE VARIATION OF COAL SAMPLES FOR DIFFERENT CYCLIC LOADING AND UNLOADING STRESS PATHS Evolution Law of Coal Sample Energy for Different Cyclic Loading and Unloading Paths
With regards to energy, during the process of rock deformation and failure, the rock always exchanges energy with the outside world. Part of the energy absorbed by rock is used for the internal damage and plastic deformation of rock, namely the dissipated energy, while the other part of the energy is stored in the form of elastic energy and released in the process of external load reduction. Thus, the energy transferred to the rock sample can be regarded as the sum of the elastic energy and the dissipated energy: where U is the energy stored by the test loaded sample, U e is the elastic energy stored in the sample during test loading, and U d is the energy dissipated by the sample during the test loading process. These terms can be obtained with the following formula: where u is the energy storage density, which is represented by the area surrounded by the loading curve and the strain axis, u e is the elastic energy density, which is represented by the area surrounded by the unloading curve and the strain axis, and u d is the area surrounded by the loading curve, the unloading curve, and the strain axis, as shown in Figure 9.    In summary, the elastic energy density and the dissipated energy density of the loaded rock mass can be obtained by calculating the area indicated in Figure 9.
According to the above calculation method for the dissipated energy density, the stress-strain curves of the coal samples for different cyclic loading and unloading modes are integrated to obtain the total energy, elastic energy, and dissipated energy corresponding to each stress cycle, as shown in Table 2 and Figures 10, 11.
For the three cyclic loading and unloading stress paths, the elastic energy accounts for most of the total energy absorbed by the rock, indicating that the energy transferred to the coal sample is mainly stored in the form of elastic energy before rock failure. The energy evolution of the coal sample is related to the external load value and the number of cycles, which is most affected by the external load value. For the same load level, the elastic energy of the coal sample remains the same, but the dissipation energy undergoes variation. For the cyclic increase condition, the total energy, elastic energy, and dissipated energy show nonlinear increasing trends, and the proportion of dissipated energy becomes larger, indicating that the upper limit of the load increases, the cracks in the sample continue to expand, and damage development accelerates. During constant amplitude cyclic loading and unloading, the elastic energy shows an overall increasing trend, indicating that with the increasing number of cycles, the inner rock gradually hardens. The dissipated energy decreases gradually, then suddenly decreases after the first cycle, and then decreases after the second cycle, indicating that the primary pores, cracks, and secondary cracks of the sample close during cyclic loading. The energy of each section is similar to that of the constant amplitude. When the loading process enters a new stage, the energy increases rapidly, which confirms that new damage occurs in the coal sample for the new load level.
As can be seen from Figure 12, regardless of the loading mode, the dissipated energy ratio reaches the maximum value between 0.14 and 0.22 in the first cycle, then decreases slowly in the subsequent cycle, and tends to be stable between 0.07 and 0.09. This shows that in the first cycle, the energy is due to the original crack closure and the friction effect of consumption. In the subsequent cycle, the elastic energy stored in the coal sample increases continuously, and the elastic property will be released immediately after the coal sample is damaged; this is why the coal sample damaged piece splash phenomenon occurs. During the cyclic loading and unloading at constant amplitude, the dissipated energy fluctuates, which is because the upper limit of the stress is close to the uniaxial failure strength of the coal samples, and cracks are connected in the coal body during cyclic loading. In the piecewise loading process, the proportion of elastic energy increases with the increase of the upper limit of stress, which is similar to the change of the proportion of dissipated energy in the increasing cycle of loading and unloading, indicating that when the external load increases, the energy is always stored as elastic energy first. Since the upper limit of the stress is stable, the decreasing trend of the second stage is similar to that of constant amplitude cyclic loading and unloading and the dissipated energy proportion increases, indicating that there is a phenomenon of unstable development of cracks in the loading process of the coal samples.

Coal Sample Damage Evolution Law for Different Cyclic Loading and Unloading Paths
In solid materials, damage refers to the existence and development of internal cracks. In the material, these defects are discontinuous, but at larger sizes, they are considered continuous. There are many ways to define damage variables. One method involves variables based on continuum damage theory (CDT), such as the residual axial  Frontiers in Earth Science | www.frontiersin.org January 2022 | Volume 10 | Article 818452 strain, residual volume strain, and elastic modulus. Another method involves variables based on energy methods, such as dissipated energy, acoustic emission, and ultrasonic speed. Regardless of the parameters used, the damage evolution must be consistent with microcrack initiation, stability, and unstable propagation (Xiao et al., 2010). A large number of researchers have applied the dissipative structure theory to the field of rock mechanics. They believe that the damage and failure process of rock is an irreversible process of energy dissipation, and the change of rock mechanical properties can be accurately reflected by defining specific damage variables. In this research, the ratio of the accumulated dissipated energy of a certain cycle to the accumulated total strain energy is defined as the damage variable of this cycle, and the total strain energy is the sum of the subelastic strain energy and the accumulated dissipated energy, as shown in Formula (8). Then, the coal sample damage variable for each cycle is calculated as follows: where i is the ith cycle, j is the total number of cycles, D i is the damage change of the ith cycle, U d i is the accumulated dissipated energy of the ith cycle, U e j is the elastic energy of the jth cycle, and U d is the total dissipated energy after the rock is destroyed.
It can be seen from the above formula that when the sample is not loaded or in a fully elastic strain state, the dissipated energy U d 0 of the sample is 0, and the damage of the sample is 0. When the coal sample is damaged during the loading process, the dissipated energy U d i > 0 and the damage is irreversible. Then, the damage variable D is defined as: According to Formula Eq. 9, when the rock is completely destroyed, i = j, D = 1, and there is an important loading and unloading process. The range of D is between 0 and 1, which conforms to the defined rules of the damage variables. According to this rule, the damage variable values of the coal samples for various cyclic loading conditions are shown in Table 3.
It can be seen from Table 3 that for different cyclic loading and unloading paths, the damage degree is positively correlated with the loading and unloading times and the load size; that is, the damage degree of the samples increases with the increase of the cycle times and the load, which conforms to the basic characteristics of rock mechanics. Before loading, the damage degree for all samples is 0, and after loading, the damage variable is about 0.4. In addition, according to previous research experience, when the damage degree of the rock sample exceeds 0.5, the sample is easily destroyed in the subsequent loading process, indicating that it is difficult to prepare coal and rock samples with damage degrees exceeding 0.5. Figure 13 shows that the loading curve and the unloading curve of the coal sample coincide with the constant amplitude cycle of 0.3σ ucs . It can be seen from the figure that the sample damage changes slowly in the process of constant amplitude loading with a small upper limit of the load. With an increase in the load, the damage of the sample increases nonlinearly; that is, the rate of damage increases continuously. However, for a constant amplitude load, the damage to the sample increases linearly within six cycles, and the rate of increase is stable between 0.04 and 0.05.

CONCLUSION
In this research, the acoustic emission and energy evolution characteristics of coal samples are studied for different cyclic loading and unloading conditions. The damage variables of coal samples are defined with the energy method, and the damage evolution characteristics of coal are studied. The main conclusions are as follows: 1) For different stress paths, the AE law of the sample is highly correlated with the size of the external load. The closer the load is to the compressive strength of the coal sample, the more active the AE is, and cracks are constantly generated and expanded. For the low load cycle, the activity of the AE phenomenon decreases with the increase of the number of cycles, but for the cycle close to the compressive strength of the coal sample, the AE phenomenon is always active. The Felicity phenomenon appears in the acoustic emission of the coal samples, and the Kaiser effect appears in the unloading FIGURE 13 | The relationship between the damage and cycle loading and unloading times.
Frontiers in Earth Science | www.frontiersin.org January 2022 | Volume 10 | Article 818452 process, indicating that part of the plastic deformation is restored in the increasing cyclic unloading stage. 2) For the different types of cyclic loading and unloading described in Section 3, the disturbance of the increasing cyclic mode for the coal sample is the most obvious. With this stress path, the total absorbed energy, elastic strain energy, and dissipated energy of the coal increase nonlinearly, and the more elastic strain energy is stored, the more fragments are produced by coal sample failure. The elastic energy changes only slightly with the subsection loading and constant amplitude loading.
3) The damage degree of the coal sample increases with the increasing number of cycles. The smaller the load is, the smaller the damage degree is with the increasing number of cycles. The increase of the cyclic loading and unloading has the greatest influence on the damage to the coal, and the damage degree increases exponentially, which is not conducive to the stability of the coal samples. This shows that the design of a UPSH with an appropriate water level can benefit the stability of the surrounding rock of the chamber.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.