Abstract
Brittleness is one of the most significant properties of geomaterials. However, very few studies have been conducted on factors influencing the rock brittleness indices. In this paper, conventional triaxial compression tests were carried out to investigate the effects of confining pressure and bedding angle on the brittleness of slate. From the perspective of energy, brittleness is an index that could reflect the release rate of energy that accumulated in the slate under the effect of external energy after reaching peak strength. Therefore, a new brittleness index of slate based on postpeak energy release is proposed herein. The applicability of this index is illustrated by comparing with other five existing brittleness indices. The following results can be obtained. (1) The confining pressure exerts a great influence on the brittleness of slate. With the increase of confining pressure, the brittleness of slate decreases significantly. The dispersion of brittleness values of slate declines with increasing confining pressure. (2) There is a parabolic relationship between slate brittleness and bedding angle. As bedding angle increases, the brittleness is intensified and reaches its maximum at a bedding angle of about 45° and then decreases gradually. (3) In contrast to the previous indices, the brittleness index proposed in this paper can describe the whole process of the postpeak stage through an index of the energy release, which makes this measure more suitable for rock that has the characteristics of step-drop or bench-drop at the postpeak stage.
1. Introduction
Brittleness is one of the most significant mechanical properties of geological materials [1]. On the one hand, brittleness indices are widely used in rock engineering. For instance, rock brittleness is used to characterize the rock drilling performance during the tunnelling excavation process using a TBM [2, 3]; shale brittleness is utilized to evaluate the mechanical properties and fracturing effects of shale gas reservoirs [4, 5]. In addition, some brittle rock, such as granite, has outstanding applications in geological engineering due to its stability [6]. The thick granite layer in Beishan, China, was chosen as the underground nuclear waste disposal base [7, 8], and the underground laboratory of Atomic Energy Canada Limited was also constructed in a granite layer [9, 10]. On the other hand, the failure process of brittle rock generally releases a lot of energy and leads to dynamic disasters in rock engineering. For instance, the failure of brittle rock in deep underground excavation projects regularly triggers rockbursts [11, 12], and the mining process in deep coal mines with hard surrounding rock might induce coal bumps [13, 14].
The research on the rock brittleness has significance in rock engineering. Generally speaking, rock brittleness is a concept contrary to rock plasticity, referring to rock fractures at or slightly beyond the yield stress. To characterize rock brittleness, a number of brittleness indices have been proposed. These main brittleness indices can be divided into the following three types, respectively, based on the strength parameters [15–19], deformation parameters [15, 20–22], and energy parameters [23–26]. Among these indices, brittleness indices based on strength and deformation parameters have a shortcoming, i.e., these indices fail to define the rock scale from brittleness to ductility [26]. In contrast, brittleness indices based on energy parameters are gaining popularity for the conciseness. In general, the aforementioned criteria have observably promoted research into rock brittleness; however, research into the factors influencing rock brittleness indices is insufficient. In respect to confinement, some scholars [27, 28] characterized the strength of intact brittle rock considering confinement using an s-shaped failure criterion highlighting the transition of failure modes, while the research [29] revealed the brittle behavior of soft calcarenites under low confinement regimes [30]. The drained triaxial tests were conducted on Beaucaire marl under a wide range of confining pressure and the brittle-to-ductile transition was described [30]. A study [31] investigated the effect of bedding structure on the brittle failure of deep shale. In addition, an evaluation method for the strength and failure of brittle rock containing initial cracks was proposed [32]. Besides, [33] the effect of porosity on rock brittleness was investigated.
The statement above indicates that there are very few studies considering the factors influencing the rock brittleness indices. There are bedding planes in rock mass related to geotechnical engineering operations and there are different confining pressures acting on the surrounding rocks in underground engineering works. Laboratory tests, including triaxial compression tests [34], acoustic emission tests [35], indentation tests [36], and numerical simulations [37, 38], show that bedding angle and confining pressure directly affect the mechanical behavior of rock, such as its strength and failure mode. The study of the influence of bedding angle and confining pressure on brittleness has certain significance for supporting design of rock slopes and underground tunnels, quantitative evaluation of surrounding rock stability, and installation of pipelines and installation of pipelines and conduits in shallow formations [38, 39]. Bedding angle and confining pressure have a significant influence on the brittleness of rock; however, there is a lack of research into the brittleness of rocks due to bedding angle and confining pressure. Therefore, it is valuable to study the influence of bedding angle and confining pressure on brittleness.
In this paper, conventional triaxial compression tests were conducted to investigate the effects of confining pressure and bedding angle on the brittleness of slate. In the second section, the preparation of the specimens and the test processes are introduced. In the third section, the experimental results under different confining pressures and bedding angles are demonstrated. In the fourth section, a new brittleness index of slate based on postpeak energy release is presented. The effects of confining pressure and bedding angle on the brittleness of slate are analyzed by using the new brittleness index in this section. In the fifth section, the applicability of the index proposed by this paper is illustrated by comparison with five existing brittleness indices. The new brittleness index based on energy release is of great significance to the quantitative evaluation of surrounding rock stability.
2. Experimental Setup
2.1. Site Sampling
In order to demonstrate the influences of bedding and confining pressure on the brittleness of rock, the selected rock specimens should have a developed bedding structure and relatively high brittleness. So, slate was selected as the research sample in this case. According to the sampling situation in the field (Figure 1), the sample is a sandy slate. The selected sample is composed of detritus and limestone. In the sample, there is a porous form of cementation, which is mainly composed of siliceous matter, calcium carbonate, clay, iron oxide, calcium sulphates, etc. Testing showed that the slate has a loose structure and poor weathering and water resistance. Through scanning electron microscopy and by extracting salient microcrack parameters of specimens, it is found that the microcrack area ratio of these specimens in a direction parallel to the bedding plane is between 0.024 and 0.036, and the microcrack area ratio normal to the bedding plane is between 0.034 and 0.051 [40], which indicates that the microcrack density dispersion of these specimen is small and has good consistency.

(a)

(b)

(c)
2.2. Test Scheme
To study the influences of bedding and confining pressure on the brittleness of slate, it is necessary to obtain the complete stress-strain curves of slate under different bedding angles and confining pressures. The confining pressures were set to 0 MPa, 10 MPa, 20 MPa, 40 MPa, and 60 MPa: the bedding angles were set to 0°, 15°, 30°, 45°, 60°, 75°, and 90°. The experimental scheme is illustrated in Figure 2.

2.3. Preparation of Specimens with Different Bedding Angles
The specific preparation steps of slate specimens with different bedding angles are as follows:(1)Firstly, the bulk rock specimens were placed on a workbench, and then the dip direction and dip angle (φ) of the sample were measured with a geological compass. After that, the slate specimens were correspondingly rotated in the vertical plane where the dip direction was located (the rotation angle is −φ). Then the rotated rock specimens were fixed by using clamping device; then the bottom surface of the rock sample was flattened by a mechanical grindstone. Then the dip direction and dip angle (φ) of the sample were measured again with the geological compass. The specimens were rotated and flattened again if the bedding angle is greater than 1°. Thus, slate specimens with horizontal bedding were obtained.(2)The slate specimens obtained in Step (1) were placed on the working platform of drilling machine. As shown in Figure 3(a), the bedding angle of the sample was adjusted so that the true dip angle is equal to the designed angle β, as shown in Figure 2.(3)The rotated rock specimens were then fixed using a clamping device to ensure that there was no sliding in the sample-drilling process. The specimen preparation and triaxial compression tests should indeed refer to the ISRM standard (φ 50 mm × 100 mm); however, due to the significant influence of the bedding planes during slate processing, the specimens are often badly damaged. According to the results of Yang et al. [41], the triaxial compression strength of rock specimens measuring φ 37 mm × 74 mm is slightly different from that at φ 37 mm × 74 mm; therefore, a core-drill with an internal diameter of 37 mm was used to drill the specimens as shown in Figure 3(b). A number of cylindrical cores with a diameter of 37 mm and bedding angle β were obtained as shown in Figures 3(c) and 3(d).(4)The true dip angle of the bedding of sample was adjusted and Steps (2) and (3) were repeated. Finally, all the specimens with bedding angles of 0°, 15°, 30°, 45°, 60°, 75°, and 90° were acquired: the number of cores with each bedding angle was no less than 12.(5)These cores were trimmed and sliced to cylinders with a height of 74 mm, as shown in Figure 3(e). According to the ISRM suggested method for preparing specimens for conventional triaxial compression testing, the two ends of each specimen were ground to ensure that both planes are horizontal within accuracy of ±0.05 mm and perpendicular to the longitudinal axis within ±0.25°.

(a)

(b)

(c)

(d)

(e)
2.4. Experimental Apparatus
The triaxial quasistatic compression test was carried out using an MTS815.03 testing machine (MTS Corporation, USA). It is mainly used for mechanical tests of materials such as rock and concrete. It is equipped with a servocontrolled, fully automatic, triaxial pressurization and measurement system (Figure 4). An axial force up to 4600 kN and a confining pressure of up to 140 MPa can be applied.

The experimental process is as follows:(1)The prepared rock specimens are numbered one by one with the format sample of “β–n,” where β represents the bedding angle of specimen, n denotes the number of specimens.(2)After measuring the size of the specimens, it was placed in the loading chamber. The axial displacement sensors and the circumferential displacement sensors were installed.(3)The control mode is displacement-controlled: specimens were loaded at a rate of 0.01 mm/s. The predetermined initial confining pressure (0, 10, 20, 40, and 60 MPa) was applied progressively as a static hydraulic pressure. During the loading process, the data sampling frequency was 2 Hz. The specimens were loaded until they were destroyed. During the experiment, the data were collected automatically and the complete stress-strain curves were drawn.(4)After the specimen was destroyed, the test was stopped, the specimen was removed, and the form of destruction of the specimen was recorded and described.
3. Test Results
Some 74 slate specimens were successfully tested in conventional triaxial compression experiments. A three-dimensional scatter plot of the relationship between compressive strength, confining pressure, and bedding angle is shown in Figure 5: the compressive strength increases linearly with increasing confining pressure. The typical stress-strain curves of these specimens are shown in Figures 6–10.


(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)

(a)

(b)
It can be seen from Figures 6–10, that the stress and strain curves of slab rock are affected by confining pressure and bedding angle. The compressive strength is enhanced with the increase in confining pressure. The slate stress-strain curves show different ways of falling at the postpeak stage. Some curves exhibit direct drops, as shown in specimen 0-3 (bedding angle, 0°): some curves drop by one level, as shown in specimen 30-1 (bedding angle, 0°). There are also specimens showing multiple steps, as shown in specimen 75-5 (bedding angle, 75°). It can be seen that the shape of the postpeak curve varies under different bedding angles. In summary, the slate stress-strain curve is affected by confining pressure and bedding angle both at the prepeak stage and the postpeak stages. The brittleness is also affected by confining pressure and bedding angle. Therefore, the study on the influences of confining pressure and bedding angle on the brittleness of slate is of significance during the excavation of this type of chamber in practical engineering works.
4. Effects of Confining Pressure and Bedding Angle on the Brittleness of Slate
4.1. A Brittleness Index Based on Postpeak Energy Release for Slate
4.1.1. Energy Accumulation at the Prepeak Stage
According to the law of conservation of energy, the amount of work done by an external force on an object is equal to the amount of energy change in the object:where W is the amount of the work done by the external force on the object, and ∆U represents the amount of the energy change in the object.
Energy accumulation mainly occurs in the prepeak stage of the stress-strain curve of slate subjected to conventional triaxial compression:where is the energy accumulated at the axial prepeak stage in Figure 11, refers to the energy accumulated at the circumferential prepeak stage in Figure 11, is the increment of the axial strain energy at the prepeak stage in Figure 11, , represents the increment of the circumferential strain energy at the prepeak stage in Figure 11, V is the volume of the specimen, is the axial stress-strain curve for slate, is the confining pressure, and is the circumferential strain on the specimen at the peak stress.

Most of the work done by the external load before the peak strength is converted to strain energy of the specimens, and a small percentage of the work is converted to the energy dissipated in crack propagation and so on. Since slate is a hard rock, the uniaxial compressive strength of slate is usually over 100 MPa. According to the analysis of the stress-strain curves of the slate specimens, the stage of nonlinear elastic deformation and the plastic stage of the slate are both small. Therefore, the energy dissipation in the nonlinear elastic deformation and the plastic stages is not considered in the analysis of the brittleness index of slate. The total mechanical work done by the mechanical press iswhere is the energy accumulated at the prepeak stage.
4.1.2. Unsteady Dissipation of Energy at the Postpeak Stage
The slate specimens were destroyed at postpeak stage, and the energy was released at this stage. The mechanical work done by the mechanical press was converted into dissipative energy at the postpeak stage. Besides this, the strain energy accumulated before reaching peak strength was also converted into dissipative energy until the residual stage. The mechanical work done by the press iswhere is the work done by press at the axial postpeak stage from the peak strength to point B in the stress-strain curve, denotes the strain at point P, is the work done by press at the circumferential postpeak stage, and refers to the work done by the press in the postpeak stage.
The residual strain energy after the peak iswhere is the residual strain energy after the peak, E is the deformation modulus of the specimen, is the compressive strength, and represents the circumferential confining pressure.
From Equations (3), (5), and (6), the increment in dissipated energy after peak strength iswhere is the mechanical work done by the press at the postpeak stage in Figure 11.
4.1.3. A New Brittleness Index Based on the Unsteady Release of Energy
The brittleness is an index that reflects the release rate of energy accumulated in the slate under the effect of external energy input after the peak strength from the perspective of energy, which iswhere Br is the brittleness of slate, represents the increment of dissipation energy released at the postpeak stage, and denotes the mechanical work done by the press after reaching peak strength.
The brittleness of slate can be obtained from Equations (5), (7), and (8):
4.2. Brittleness Calculation
The brittleness indices of slate could be calculated by using Equation (9). The calculated results are displayed in Table 1.
4.3. Correlation Analysis of Brittleness and Confining Pressure
A three-dimensional scatter diagram showing brittleness, confining pressure, and bedding angle is shown in Figure 12. A fitting curve between brittleness and confining pressure is illustrated in Figure 13.


It can be seen that there is a negative correlation between the brittleness and the confining pressure from Figure 13. With the increase of confining pressure, the brittleness of slate decreases. The conclusion is consistent with the general understanding among current scholars. With the increase of confining pressure, there is a transition in the rock from brittleness to plasticity. As the confining pressure increases, the lower the brittleness of the rock and the less prone it is to brittle failure; therefore, this property is often used in engineering practice. If the exposed rock is subjected to confining pressure, e.g., by shotcreting, the possibility of brittle failure by rock bursting can be reduced.
Furthermore, the dispersion of brittleness of slate decreases with the confining pressure increasing.
Confining pressure exerts a significant influence on the brittleness of slate. This is because the energy accumulated at the axial prepeak stage (Wmf·A) shows significant linear growth with increasing confining pressure, while the residual strain energy at the postpeak stage (WR) and the work done by the press in the axial postpeak stage (Wmb·A) show exponential growth with the increase of confining pressure (Figure 14). Substituting Wmf·A, Wmb·A, and WR into Equation (9), it can be found that the brittleness shows exponential growth with increasing confining pressure.

(a)

(b)

(c)
4.4. Correlation Analysis: Brittleness and Bedding Angle
A fitting curve between brittleness and bedding angle is shown in Figure 15. It can be seen that there is a parabolic relationship between slate brittleness and bedding angle. The brittleness is minimised when the bedding angle is 0°. With increasing bedding angle, the brittleness increases and reaches a maximum when the bedding angle is about 45° and then decreases gradually thereafter. When the bedding angle is 45°, the slate specimen is almost destroyed along the bedding plane according to the Jaeger failure criterion. Due to the weak strength of the bedding plane, the specimen is easily fractured along the bedding plane, showing a relatively strong brittleness. In addition, it can be seen from the comparison that the average value of slate brittleness is 4 to 12 under the influence of confining pressure, and 3 to 6 under the influence of bedding angle. This indicated that the confining pressure has a much greater influence on the brittleness of slate than the bedding angle.

Bedding angle has a remarkable influence on the brittleness of slate. It can be seen, from Figure 16(a), that when the bedding angle increases, the energy accumulated at the prepeak stage (Wmf·A) shows a parabolic relationship whereby it first decreases and then increases; therefore, when the bedding angle is 0°and 90°, the energy accumulated at the axial prepeak stage is larger, while the smallest value occurs at 45°. From Figures 16(b) and 16(c), it can be seen that the residual strain energy at the postpeak stage (WR) and the work done by the press at the axial postpeak stage (Wmb·A) show a parabolic relationship, whereby it first increases, and then decreases, with increasing bedding angle; however, the residual strain energy at the postpeak stage is greater when the bedding angle is 90° than when the bedding angle is 0° and vice versa for the work done by press at the axial postpeak stage. Substituting Wmf·A, Wmb·A, and WR into Equation (9), the brittleness has a parabolic relationship whereby it first increases and then declines with increasing bedding angle.

(a)

(b)

(c)
5. Discussion
With the development of rock mechanics, many kinds of brittleness indices have been proposed based on stress-strain curves for different purposes and research objectives; therefore, in this section, the applicability of the index proposed here is illustrated by comparison with five existing brittleness indices.
The typical brittleness indices are as follows:(1) is a brittleness index based on the area under the stress-strain curvewhere A1 is the area under the line that the slope of it is E50 (E50 is the deformation modulus at 50% compressive strength), and it passes through the point of peak compressive strength and A2 is the area under the stress-strain curves.(2) is a brittleness index based on the axial peak strain and the axial residual strain:where is the axial peak strain and refers to the axial residual strain.(3) is a brittleness index based on compressive strength and residual strengthwhere is the compressive strength and is the residual strength.(4) is a brittleness index based on the relative amount, and absolute rate, of stress drop at the postpeak stagewhere σA·P represents compressive strength, σA·R is residual strength, |KAC| is the slope of the line AC where A is the onset of yielding, and C is the residual point of the yield process.(5) is a brittleness index based on peak strain, residual strain, and the residual strengthwhere is the axial peak strain, represents the axial residual strain, and is the strain at the prepeak stage such that the stress is equal to the residual strength.
It can be seen from Figure 17 that the elastic modulus, peak strength, and residual strength of specimens 30-13 and 60-3 are similar before their peak compressive strength was mobilized. In the earlier stage of the postpeak stage, the stress-strain curves of the specimens are similar: in the postpeak stage, the stress on 30-13 drops at a constant rate, while there are multiple stress-drops seen on specimen 60-3. The mechanical work done by the press, on the two specimens, before their peak strength is similar but there is more energy released from sample 60-3. As a result, the brittleness of specimen 30-13 should be greater than that of 60-3. Therefore, a good brittleness indicator can not only reflect the conclusion that brittleness of 30-13 is greater than that of 60-3 in a qualitative sense, but also could distinguish the brittleness of the two specimens quantitatively.

The brittleness values of specimens 30-13 and 60-3 were calculated by using the proposed brittleness index and the five existing brittleness indices mentioned above (Table 2).
It can be seen from Table 2 that the brittleness values of samples 30-13 and 60-3 calculated by the brittleness index proposed by this paper are 2.948 and 2.372, indicating that there is a significant difference between the brittleness of the two specimens. Therefore, the brittleness index proposed by this paper can distinguish the brittleness of the two specimens; however, the brittleness values calculated using the five brittleness indicators mentioned above were very similar. Some results show that the brittleness of sample 30-13 calculated by the five brittleness indicators is less than that of specimen 60-3. This is contrary to the facts; therefore, the new brittleness index was deemed more suitable for this particular slate.
To investigate the applicability of the five existing brittleness indicators mentioned above for slate, the brittleness of the slate specimens was calculated using the five brittleness indicators. The relationship between the magnitudes of brittleness and confining pressure was obtained and the results are demonstrated in Figure 18.

(a)

(b)

(c)

(d)

(e)
It can be seen from Figure 18 that the brittleness values of the specimens calculated by different indices all decrease with increasing confining pressure. The conclusion is consistent with the facts as observed experimentally; however, the brittleness values calculated by indices (1), (2), (4), and (5) all decrease within a small range with increasing confining pressure, which cannot reflect the effect thereof. For index (3), although effect of the confining pressure on the brittleness calculated by index (3) is obvious, the brittleness values of different specimens under different confining pressures are all within 0 to 1: this indicates insufficient ability of index (3) to distinguish the brittleness of different specimens. The relationship between the brittleness calculated by the brittleness index proposed by this paper and confining pressure is shown in Figure 18: the effect of confining pressure is obvious, and with the increase of confining pressure, the variation in brittleness decreases and the difference between different specimens is obvious.
To sum up, no matter whether it is a reflection of the stress behavior in the postpeak stage of a single specimen, or the rules of the effect of confining pressure of all specimens, the brittleness index proposed by this paper is more suitable for use on slate samples. This is because most previous indices are based on data at peak compressive strength and residual strength. From the conventional triaxial compression tests of slate (Section 3), it can be seen that the stress-stain curves of slate indicate different types of postpeak behaviors; however, considering only the two points of compressive strength and residual strength, it is difficult to describe the brittleness of slate, for two points cannot reflect the true postpeak behavior with sufficient accuracy. In contrast to the previous indices, the brittleness index based on postpeak energy release proposed in this paper can describe the whole postpeak stage by considering the energy released in the postpeak stage. Therefore, for rock exhibiting characteristic step-drop, or bench-drop, in the postpeak stage, the brittleness index based on postpeak energy release proposed in this paper is more suitable.
6. Summary and Conclusions
Brittleness is an important indicator of failure in geological materials. The brittleness is an index that could reflect the release rate of energy that accumulated in the slate under the effect of external energy after the peak strength from the perspective of energy. Therefore, a brittleness index based on postpeak energy release for slate was constructed in this paper, which is defined as the increment of dissipated energy released at the postpeak stage divided by the mechanical work done by the press tester after reaching the peak strength. Then conventional triaxial compression tests of 74 specimens were carried out to reveal the effects of confining pressure and bedding angle on the brittleness of slate. The following results were obtained:(1)The confining pressure exerts a significant influence on the brittleness of slate. With increasing confining pressure, the brittleness of slate shows a rapid decrease. The dispersion of brittleness values of slate decreases with increasing confining pressure.(2)There is a parabolic relationship between slate brittleness and bedding angle. The brittleness is the smallest when the bedding angle is 0°. With the increase of bedding angle, the brittleness increases and reaches its maximum when the bedding angle is about 45°, and it then decreases gradually thereafter.(3)Most previous indices are based on data at compressive strength and residual strength points; however, considering only these two points makes it difficult to describe the brittleness of slate, for use of only two points cannot reflect postpeak stage behavior with sufficient accuracy, especially for the rock that undergoes different types of postpeak stress decrease, such as slate.(4)The brittleness index based on postpeak energy release proposed in this paper can describe the whole of the postpeak stage through an index based on the energy released in the postpeak stage. Therefore, for a rock that exhibits characteristic step-drop, or bench-drop, in its postpeak stage, the brittleness index based on postpeak energy release proposed in this paper is more suitable.
Abbreviations
D: | Diameter of specimens |
H: | Height of specimens |
: | Axial strain of the specimens at the peak strength |
: | Compressive strength |
: | Axial strain at the eventual residual stress |
: | Circumferential strain of the specimens at the peak stress |
: | Residual strength |
: | Circumferential strain at the eventual residual stress |
E: | Deformation modulus of the specimens |
ES: | Secant modulus |
W: | Amount of the work done by the external force on the object |
∆U: | Amount of the energy change in the object |
: | The energy accumulated at the axial prepeak stage |
: | Increment of the axial strain energy at the prepeak stage |
V: | Volume of the specimens |
: | Axial stress-strain curve |
: | Confining stress-strain curve |
: | The energy accumulated at the circumferential prepeak stage |
: | Increment of the circumferential strain energy at the prepeak stage |
: | Confining pressure |
: | The energy accumulated at the prepeak stage |
: | The work done by press tester at the axial postpeak stage |
: | Residual strain energy |
: | Increment of dissipation energy released at the postpeak stage |
: | The work done by the mechanical press at the circumferential postpeak stage |
: | The work done by the mechanical press at the postpeak stage |
: | The brittleness of slate |
A1: | The area under the line such that its slope is E50 |
A2: | The area under the stress-strain curves |
E50: | The deformation modulus at 50% compressive strength |
: | Bedding angle |
: | Strain at the prepeak stage that the stress is equal to the residual strength |
Sig.: | Significance |
F: | F-test value. |
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by National Natural Science Foundation of China (grant no. 51804309), the National Key R & D Programme of China (2016YFC0801400 and 2016YFC0600700), and the State Key Laboratory of Coal Resources and Safe Mining (China University of Mining and Technology) under grant nos. SKLCRSM16KFB07, SKLCRSM16DCB01, and SKLCRSM17DC11, funded by the Key Laboratory of Coal Mine Safety and High-Efficiency Mining as coestablished by the Province and the Ministry (JYBSYS2018201), supported by the Yue Qi Distinguished Scholar Project, China University of Mining and Technology (Beijing), funded by Young Elite Scientists Sponsorship Program by CAST (2017QNRC001), and was also supported by Open Fund of State Key Laboratory of Water Resource Protection and Utilization in Coal Mining (grant no. SHJT-17-42.10). All of these supports are gratefully acknowledged.