1. Introduction
The energy dissipation and fracture mechanism of rocks under coupled hydro-mechanical unloading is not only one of the research hotspots in geotechnical engineering, but is also one of the key scientific issues to be solved in deep engineering [
1,
2,
3]. In addition, related results showed that rock burst is easy to occur in the water-rich roadway of coal mines, which is closely related to the energy dissipation and fracture mechanism of rocks under coupled hydro-mechanical unloading [
4,
5,
6].
Generally, the joints of rocks develop actively, which show significant anisotropy [
7,
8]. Meanwhile, sandstone, as one of the geological rocks, is common in numerous large-scale engineering [
9,
10,
11]. Further, there are the evolution characteristics of in situ stress on the surrounding rocks of layered water-rich roadways before, in and after excavation in
Figure 1 [
12,
13]. Among them, the surrounding rocks of layered water-rich roadways before excavation are subjected to the long-term comprehensive actions of high geo-stress and strong osmotic pressure, which result in the high-density compaction of layered rock strata, and are accompanied by accumulated energy (see
Figure 1a). The surrounding rocks of layered water-rich roadways during excavation are in the unloading state of horizontal principal stress, and the accumulated energy is released (see
Figure 1b). However, the surrounding rocks of layered water-rich roadways after excavation are subjected to the long-term superimposed actions of the self-weight of the overlying strata (vertical stress) and the mining stress. Subsequently, it easily leads to energy accumulation, and the occurrence of a series of mechanical behaviors such as large deformation or collapsed failure, and even the occurrence of rock burst (see
Figure 1c). Therefore, the energy dissipation law and fracture mechanism of layered sandstones under coupled hydro-mechanical unloading needed to be further studied.
Currently, related research on the energy dissipation law and fracture mechanism of rocks has been conducted by numerous researchers, and relatively abundant results have been obtained. Among them, the energy dissipation law of coal or rocks under a single stress field have both been studied. Subsequently, the energy dissipation law and fracture mechanism of rocks under the single stress field was better revealed [
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26]. In addition, Professor Gong’s team has conducted numerous laboratory mechanical experiments on rocks, where the linear energy evolution law and the fracture mechanism of rocks were further found and revealed [
27,
28,
29,
30,
31]. Meanwhile, numerous researchers have constructed a series of damage evolution models and failure criterions of coal or rocks. For example, the damage constitutive models of brittle rocks according to the energy dissipation law and fracture mechanism were constructed [
32,
33,
34]. Moreover, a series of failure criterions of rocks under different loading modes were also constructed [
35,
36,
37,
38].
Furthermore, the creep mechanical behavior and time-dependent initiation, propagation and interaction mechanism of rock-like mortar with three-dimensional cracks under triaxial H-M coupling were deeply studied [
39,
40]. Subsequently, it was concluded that subcritical crack propagation is the main failure cause of rocks, and the faster the propagation rate of dominant cracks, the more severe the failure rate of rocks. Moreover, the mechanical behaviors of single-fissure mudstone under triaxial compression shear with H-M coupling have been analyzed. Subsequently, it was concluded that the single-fissure mudstone was most prone to failure when the fissure angle was 45° [
41]. Additionally, the mechanical behavior of jointed rocks under H-M coupling was characterized, and its damage constitutive model was constructed [
42]. Moreover, a micromechanical hydro-mechanical-damage coupling model was established based on consideration of the multi-scale structural differences of layered rocks [
43]. In addition, the evolution law of the permeability coefficient and fracture mechanism of rocks under various confining pressures was characterized and revealed, respectively [
44]. The corresponding hydro-mechanical-damage coupling evolution model according to the damage evolution characteristics of the plastic dilatation of anisotropic rocks under H-M coupling have been also established [
45,
46].
This brief literature survey shows that the focus has mainly been on the energy dissipation law and fracture mechanism of rocks under a single stress field. However, there are relatively rare reports of coupled hydro-mechanical unloading. Therefore, the corresponding coupled hydro-mechanical unloading experimentation of layered sandstones has been conducted. Firstly, the anisotropy and competition relationship of energy are obtained. Secondly, the corresponding energy evolution model is established. Meanwhile, the energy evolution and failure mechanism are revealed. Subsequently, the occurrence mechanism of rock burst in a layered water-rich roadway is revealed according to the proposed theory of an “Energy Flow” chain. Meanwhile, the specific discriminant formula with the form of “Energy Flow” of rock burst in a layered water-rich roadway is qualitatively proposed. Finally, the nature of the effective prevention and accurate control of rock burst is obtained from three aspects, which corresponded to that controlling the starting source, cutting off the “Energy Flow” chain, and protecting the target source.
2. Preparation and Test Plans
A series of processing works were conducted on the large, layered sandstones taken from the sampling site, including drilling, cutting and polishing. Subsequently, the standard specimens were obtained.
The TOP INDUSTRIE Rock 600-50 multi-field coupling test system was used for the mechanical tests under H-M coupling. Among them, the stress paths adopted could be seen in
Figure 2. The corresponding test plans were as follows:
(a) Stage Ⅰ: Confining pressure loading stage. At this stage, the confining pressure σ3 was loaded to the corresponding pre-set value a (a = 10, 20, 30 MPa) by the hydrostatic pressure, respectively. Meanwhile, the loading rate adopted at this stage was 1 MPa/min.
(b) Stage ①: Osmotic pressure loading stage. At this stage, the axial pressure and confining pressure were kept constant. Then, the osmotic pressure was loaded to the corresponding pre-set value of 4 MPa with the stress-controlled mode. Meanwhile, the loading rate adopted at this stage was 1 MPa/min.
(c) Stage Ⅱ: Axial stress loading stage. At this stage, the osmotic pressure and confining pressure were kept constant. Then, the axial pressure was loaded until failure in terms of the strain-controlled mode. Meanwhile, the loading rate adopted at this stage was 0.06 mm/min.
Figure 2.
Triaxial loading (unloading) stress paths under H-M coupling: (a) CTLTP, (b) TLUTP.
Figure 2.
Triaxial loading (unloading) stress paths under H-M coupling: (a) CTLTP, (b) TLUTP.
- 2.
Triaxial loading and unloading test under H-M coupling (TLUTP)
(a) Stage Ⅰ: Initial high in situ state recovery stage. At this stage, the confining pressure and axial pressure were loaded to the (corresponding to the buried depth of the sampling site) by the deviatoric stress. Meanwhile, the loading rate adopted at this stage was 1 MPa/min.
(b) Stage ①: Osmotic pressure loading stage. At this stage, the was kept constant. The osmotic pressure was loaded to the corresponding pre-set value of 4 MPa in terms of the stress-controlled mode. Meanwhile, the loading rate adopted at this stage was 1 MPa/min.
(c) Stage Ⅱ: Constant axial pressure–unloading confining pressure stage. At this stage, the osmotic pressure and axial pressure were kept constant. Then, the confining pressure was unloaded to the corresponding pre-set value c (c = 5, 10, 20 MPa) in terms of the stress-controlled mode. Meanwhile, the unloading rate adopted at this stage was 1 MPa/min.
(d) Stage Ⅲ: Axial stress loading stage. At this stage, the osmotic pressure and confining pressure were kept constant. Then, the axial pressure was loaded until failure in terms of the strain-controlled mode. Meanwhile, the loading rate adopted at this stage was 0.06 mm/min.
3. Results
3.1. Stress–Strain Behavior and Fracture Modes
As shown in
Figure 3 and
Figure 4, the global radial deformation of layered sandstones under TLUTP was significantly weaker than that under CTLTP. In addition, the global axial deformation with CTLTP under an identical layered angle and various confining pressures had little differences. However, the global axial deformation with CTLTP under identical confining pressure and various layered angles was quite different. Therefore, the influence of confining pressure was not sensitive to the global axial deformation under CTLTP. However, the bedding effect was sensitive to the global axial deformation under CTLTP.
Nevertheless, the effects of confining pressure and bedding were both sensitive to the global axial deformation under TLUTP. It was also indicated that there were certain errors in estimating the deformation resistance of deep surrounding rocks only according to the results of CTLTP. Therefore, the stress path of TLUTP should be vigorously adopted to conduct a series of related research. Thus, it was expected that the deformation resistance of deep surrounding rocks could be estimated more accurately.
In addition, the peak strength with TLUTP was higher than that with CTLTP under identical working conditions. Therefore, to more accurately estimate the bearing capacity and deformation resistance of deep surrounding rocks, it was necessary to fully rely on the stress path of TLUTP to conduct the laboratory mechanical research.
In addition, as shown in
Figure 3,
Figure 4 and
Figure 5, layered sandstones under CTLTP and TLUTP showed four types of macroscopic fracture modes. Among them, the macroscopic fracture mode with layered angles of 0°, 30° and 90° basically belonged to type Ⅲ. However, the macroscopic fracture mode with a 45° layered angle belonged to type Ⅱ. Meanwhile, the macroscopic fracture mode with a layered angle of 60° basically belonged to type Ⅰ. Among them, the comprehensive effects of stress path CTLTP, bedding angle 0° and confining pressure 20 MPa led to the macroscopic fracture mode belonging to type Ⅳ. It was suggested that the influences of the stress path and confining pressure had little effects on the macroscopic fracture mode, but the bedding effect had significant influences on the macroscopic fracture mode.
In addition, the fracture modes under CTLTP and TLUTP were mainly shear fracture, which suggested that there must be a weak structure with the lowest bearing capacity in layered sandstones [
47,
48]. Under the actions of the effective triaxial loads, the weak structure could gradually develop and finally formed the shear fracture dominated by single- or crossing-fracture planes.
3.2. Strength Anisotropy and Deformation Evolution
There were significant influences of confining pressure on the peak strength under CTLTP and TLUTP (see
Figure 6). Meanwhile, with the increase of confining pressure, the peak strength under CTLTP and TLUTP was raised. However, when
β ∈ (0°, 30°) and
β ∈ (45°, 90°), the peak strength under an identical confining pressure changed little, which suggested that the influenced sensitivity of the layered angles compared to that of the confining pressure on the peak strength under TLUTP was weaker. However, when
β ∈ (30°, 45°), the influence of the layered angles on the peak strength was more sensitive than that of confining pressure.
The peak axial strain with CTLTP was generally higher than that with TLUTP under identical working conditions (see
Figure 7). Corresponding to the deep excavation engineering, there were errors in the estimating results for the deformation resistance of deep surrounding rocks in terms of the stress path CTLTP. Therefore, it is necessary to conduct a series of related studies adopting the stress path TLTUP to more accurately estimate the deformation resistance of deep surrounding rocks. In addition, the influence of confining pressure on the peak axial strain was more significant than that of layered angles.
As shown in
Figure 8, compared with TLTUP, CTLTP had the more significant influence on the peak radial strain, and the anisotropy of peak radial strain under CTLTP was more significant. In addition, the influenced sensitivity of confining pressure on the peak radial strain under CTLTP was significantly stronger than that of layered angles. However, the anisotropy of peak radial strain with TLUTP under low confining pressure was significant.
The influence of confining pressure on the peak volumetric strain under TLTUP was also significant (see
Figure 9). Due to the comprehensive actions of the layered angles, stress path, confining pressure and osmotic pressure, the evolution characteristics of peak volumetric strain under CTLTP were significantly different. In addition, the influences of layered angles were more sensitive on the peak volumetric strain under CTLTP than that of confining pressure. Nevertheless, the influence of confining pressure on the peak volumetric strain under TLUTP was more sensitive than that of layered angles.
3.3. Energy Evolution
Under the comprehensive actions of loading or unloading, the energy evolution under H-M coupling still follows the first thermodynamics law [
49], namely
where
,
and
, respectively, represent the total input energy density, elastic energy density and dissipated energy density under H-M coupling.
In addition, according to Terzaghi’s effective stress principle [
50], the effective axial stress
and the effective confining pressure
can be expressed as:
where
,
and
represent the axial stress, confining pressure and osmotic pressure, respectively.
Meanwhile, combined with
Figure 3 and
Figure 4 and Equation (2), the total input energy density under H-M coupling can be quantitatively calculated:
In addition, similarly with the calculation method of the elastic energy density under the single stress field, the elastic energy density under H-M coupling can be obtained, and the specific expression is as follows:
Combining with Equations (1)–(4), the dissipated energy density under H-M coupling can be obtained, and the specific expression is as follows:
where
E represents the elastic modulus, and
represents the Poisson’s ratio.
and
are the final data points of axial strain and radial strain collected by hydraulic coupling tests, respectively.
Therefore, the energy evolution curves under CTLTP and TLUTP are drawn in
Figure 10 and
Figure 11.
In addition, the deformation of rocks under identical stress paths has a certain “memory effect”. Naturally, the energy evolution of rocks also has certain similarities. Therefore, the energy evolution curve with CTLTP under a layered angle of 0°, confining pressure of 10 MPa and osmotic pressure of 4 MPa was adopted for representative analysis (see
Figure 10a). Similarly, the energy evolution curve with TLUTP under the layered angle of 0°, confining pressure of 5 MPa and osmotic pressure of 4 MPa was adopted for representative analysis (see
Figure 11a).
As shown in
Figure 11a and
Figure 12a, the energy evolution under CTLTP could be divided into three stages, while that under TLUTP could be divided into four stages.
Stage ① corresponded to the initial pore compaction stage (see
Figure 10a). According to Terzaghi’s effective stress principle, the effective stresses at this stage were both small. Therefore,
,
and
were all seen in the initial incubation stage with small values which, respectively, slowly and synchronously raised (0→4.52
, 0→2.77
and 0→1.75
). Stage Ⅰ corresponded to the initial high in situ stress state recovery stage (see
Figure 11a). At this stage,
and
raised quickly (0→42.25
and 0→38.61
), while
raised in a shapely manner (0→3.64
). According to Terzaghi’s effective stress principle, the existence of high effective confining pressure and high effective axial pressure made the compaction effect of the three-dimensional volume significantly enhanced at this stage. It could greatly improve the compaction degree of its internal pores and space, and even made it become rigid. The occurrence of this process would inevitably lead to significant elastic energy accumulation, and then cause
to increase rapidly. However, at this stage,
had almost no change, which was close to the rigid body.
In
Figure 10a, stage ② corresponded to the pre-peak elastic and plastic deformation stage. At this stage, numerous newly formed micro-cracks would gradually derive, develop, nucleate and even coalesce under the actions of external loads, and gradually formed the dominant seepage channels. According to Griffith strength theory [
51], numerous newly developed micro-crack tips were easy to form the high-stress concentration zone, which was accompanied by significant energy accumulation. Therefore,
raised rapidly (2.77
→42.69
). Meanwhile, the gradual formation of the dominant seepage channels would cause the seepage water to absorb the energy, and indirectly replaced the position of dissipated energy. Then, it resulted in
having almost no change at this stage, which steadily and slowly raised (1.75
→10.49
). Therefore,
raised rapidly (4.52
→53.18
). In
Figure 11a,
at stage Ⅱ rapidly decreased (38.61
→20.40
), while
raised rapidly at this stage (3.64
→60.46
). Therefore,
raised steadily (42.25
→80.85
). At this stage, the radial unloading effect would appear to a certain extent, so that the openings of internal original micropores were opened again. Meanwhile, the radial openings of the dominant seepage channels were also opened, which led to
to be rapidly released and the rapid accumulation of
at this stage.
In addition, in
Figure 10a, stage ③ corresponded to the post-peak failure stage. At this stage, under the continuous actions of effective stresses, the layered sandstones accumulated significant damage which would further form the macroscopic shear fracture plane. Meanwhile, the dominant seepage channels would also be completely formed, thus leading to the release of
. The numerical results showed that
decreased rapidly (42.69
→16.61
),
raised rapidly (10.49
→109.99
) and
raised rapidly (53.18
→126.60
). However, in
Figure 11a, stage Ⅲ-① under TLUTP was varied from stage ② under CTLTP. This stage included the secondary compaction sub-stage, linear elastic deformation sub-stage and plastic deformation sub-stage. Due to the combined actions of stages Ⅰ and Ⅱ, the stress–strain curves of the secondary compaction sub-stage and the plastic deformation sub-stage under CTLUP were not significant, and the axial deformation under CTLUP was smaller than that of the whole deformation. However, the energy evolution of stage Ⅲ-① was similar to that of stage ② in
Figure 10a. At this stage, numerous newly formed micro-cracks would gradually derive, develop, nucleate and even coalesce under the actions of external loads, and gradually formed the dominant seepage channels. According to Griffith strength theory, numerous newly developed micro-crack tips can easily form a high-stress concentration zone, which is accompanied by significant energy accumulation. Therefore,
at this stage raised rapidly (20.40
→363.48
). Meanwhile, the gradual formation of the dominant seepage channels would cause the seepage water to absorb part of the energy, and indirectly replace
. The result is that
had almost no change at this stage, which steadily and slowly raised (60.46
→70.66
). Therefore,
at this stage raised rapidly (80.85
→434.14
).
In addition, in
Figure 11a, stage Ⅲ-① was the post-peak failure stage under TLUTP. The energy evolution at this stage under TLUTP was similar to that of stage ③ under CTLTP. At this stage, under the continuous actions of effective stresses, the layered sandstone that accumulated a lot of damage would further form the macroscopic shear fracture plane. Meanwhile, the dominant seepage channels would also be completely formed, thus leading to the release of
. The numerical results showed that
decreased rapidly (363.48
→44.43
),
raised rapidly (70.66
→623.52
) and
raised rapidly (434.14
→667.95
).
3.4. Energy Competition Relationship
There was a significant energy competition relationship. Therefore, the evolution curve of the ratio of elastic energy density to dissipated energy density
under CTLTP and TLUTP with axial strain
was drawn (see
Figure 12 and
Figure 13).
Similarly, the
evolution curve with CTLTP under a layered angle of 0°, confining pressure of 10 MPa and osmotic pressure of 4 MPa was adopted for representative analysis (see
Figure 12a). Similarly, the
evolution curve with TLUTP under a layered angle of 0°, confining pressure of 5 MPa and osmotic pressure of 4 MPa was adopted for representative analysis (see
Figure 13a).
As shown in
Figure 12a and
Figure 13a, the evolution curves of
under CTLTP could be divided into three stages, which showed the inverse “
N” type evolution characteristic first decreasing, then increasing and finally decreasing again. However, the evolution curves of
under TLUTP could be divided into four stages, which showed the “
M” type evolution characteristic first increasing, then decreasing, then increasing again, and finally decreasing again.
Among them, in
Figure 12a, stage ① corresponded to the initial pore compaction stage. At this stage,
dropped rapidly from 1.01 to 0.159, which indicated that compared with
,
at this stage was dominant. It was also suggested that the accumulation rate of
at this stage was weaker than the derivation rate of
, which indirectly indicated that the initial porosity at this stage was not completely compacted. However, in
Figure 13a, stage Ⅰ corresponded to the initial high in situ stress state recovery stage. Under the comprehensive actions of high effective stresses, the initial micropores in layered sandstones at this stage were three-dimensional compacted, even close to rigid, and the compaction degree was greatly and completely improved. Therefore,
at stage I in
Figure 13a raised sharply from 0.1634 to 21.51, which indicated that compared with
,
at this stage was dominant.
In addition, in
Figure 12a, stage ② corresponded to the pre-peak elastic and plastic deformation stage. In this stage,
raised sharply from 0.159 to 4.26, which indicated that
in this stage was gradually displacing
. Finally,
was gradually dominant in this stage. It also indirectly indicated that the derivative rate of
required for the derivation, development and nucleation of pre-peak nascent cracks and the gradual formation of dominant seepage channels was significantly weaker than the accumulation rate of
. Then, it could also verify the rationality of Griffith’s strength theory. However, in
Figure 13a, stage Ⅱ corresponded to the constant axial pressure–unloading confining pressure stage. In this stage,
dropped rapidly from 21.51 to 0.55, and the radial openings of micro-cracks opened rapidly. Then, it could provide the adequate spaces for the release of
accumulated in stage I, which led to the sharply rising growth rate of
at this stage. The dropped characteristic of
was significant, and
was dominant at this stage.
In addition, in
Figure 12a, stage ③ corresponded to the post-peak failure stage. At this stage,
dropped rapidly from 4.26 to 0.146 again. Due to the formation of a macroscopic fracture plane, the
that was accumulated was rapidly dropped, which led to
sharply raising. Compared with
,
was dominant at this stage. In
Figure 13a, stage Ⅲ-① corresponded to the pre-peak elastic and plastic deformation stage. Similarly, compared with the
of stage ② under CTLTP, the
of TLUTP raised rapidly from 0.55 to 5.23 again, and
was dominant. However, compared with the
of stage ② under CTLTP, the increased amplitude of
under TLUTP was gentler, which was the action result of the stage Ⅱ.
Therefore, it was noted that the peak of
at stage Ⅲ-① was labeled as
, and the peak of
at the end of stage Ⅰ was “spurious” peak, which was not included in the subsequent analysis. In addition, in
Figure 14a, stage Ⅲ-② corresponded to the post-peak failure stage. Similar to the
of CTLTP, the
with TLUTP dropped rapidly from 5.23 to 0.21 again in this stage. Compared with
,
was dominant in this stage.
In addition, the
of TLUTP was basically lower than that of CTLTP under identical working conditions (see
Figure 14a–c). Among them, the proportion of
with CTLTP exceeded that with TLUTP, which could reach the highest value of 95.81%. It also indicated that the fracture degree under TLUTP was much weaker than that under CTLTP.
In addition, as shown in
Figure 14d, the criterion for the fracture degree was determined according to the distribution range of
. Among them, when
∈ [0, 3], it was the weak fracture degree; when
∈ (3, 9], it was the medium-strong fracture degree; when
∈ (9, +∞), it was the strong fracture degree. Therefore, the fracture degree under TLUTP was the weak and medium-strong fracture, and the proportion of the weak fracture degree was higher. However, the fracture degree under CTLTP was the medium and strong fracture, and the proportion of the strong fracture degree was higher.
3.5. Energy Anisotropy
In addition, the energy evolution showed significant anisotropy under CTLTP and TLUTP. Therefore, the relationships between
,
and
under CTLTP and TLUTP needed to be further clarified (see
Figure 15).
As shown in
Figure 15a,c,
under CTLTP and TLTUP had a significant linear evolution relationship with the peak stress. However,
was more suitable to estimate the elastic energy storage limit than
.
The peak elastic energy density under CTLTP and TLUTP showed the “W” shaped evolution characteristic with increasing layered angles (see
Figure 16a,c). When the layered angle was 60°,
was basically the lowest; when the layered angle was 30°,
was the highest. It was also noted there was the significant energy anisotropy. However, the influence of the stress path on the evolution characteristics of
was significantly weaker than that of layered angles.
Similarly,
under CTLTP and TLUTP showed the “
N” type evolution characteristic with increasing layered angles (see
Figure 17a,c). When the layered angle was 90°,
under TLUTP was the highest. When the layered angle was 60°,
under TLUTP was the lowest.
Similarly,
under CTLTP showed the “
W”-shaped evolution characteristic with increasing layered angles (see
Figure 18a,c). When the layered angle was 60°,
under CTLTP was the lowest. When the layered angle was 0° or 30°,
under CTLTP was the highest. However,
under TLUTP showed the “
V” shape evolution characteristic with increasing layered angles. When the layered angle was 90° or 0°,
under TLUTP was the highest. When the layered angle was 45° or 60°,
under TLUTP was the lowest.
As the confining pressure increased,
,
and
under an identical bedding angle all showed a significant increasing trend, but the influence of a stress path was not significant (see
Figure 16b,d,
Figure 17b,d and
Figure 18b,d).
In conclusion, the effects of layered angles, H-M coupling and stress paths could influence the energy anisotropy. However, compared with the influences of H-M coupling and stress paths, the influence of layered angles on energy anisotropy was stronger.
3.6. Energy Evolution Model
As shown in
Figure 19 and
Figure 20, the nonlinear evolution relationships between
,
and
, effective for confining pressure and layered angles under CTLTP and TLUTP, were obtained. The specific expression of the energy evolution models was as follows:
.
Based on the above evolution models, the elastic energy storage capacity under CTLTP and TLUTP could be effectively and accurately estimated, which thus provided more accurate guidance and references for field engineering.
3.7. Fracture Mechanism
Combined with
Section 3.1, under the combined actions of H-M coupling, layered angles and stress paths, a series of mechanical behaviors occurred on the weakest structure, and finally the main shear fracture with single or crossing cracks was formed. However, according to Griffith’s strength theory, the energy was easy to accumulate in the tips of the weakest structure, which could be regarded as the starting source.
Therefore, as shown in
Figure 21, according to the concept of “Energy Flow” and combined with the domino effect [
48,
52], the energy accumulated in the high-stress concentration zone at the tips of the weakest structure would continuously impact the low-state damage zone, which could be regarded as the target source here. Furthermore, the gradual derivation and expansion of the weakest structure was constantly driven, so that the local damage zone gradually expanded to the main fracture zone [
47,
48].
4. Sublimation and Application of “Energy Flow” Fracture Mechanism
The proposed “Energy Flow” fracture mechanism of layered sandstone was proposed with a laboratory scale. The “Energy Flow” was a directional medium with a definite starting source and target source, thus forming a complete” Energy Flow “chain. Similarly, as shown in
Figure 22, the fracture position of the confined water key layer could be regarded as a weak structure to a certain extent, which was one of the starting sources (energy accumulation,
U1). Meanwhile, the “Energy Flow” 1 (“Longitudinal Flow”) was generated.
The existence of the goaf and the frontage of the layered roadway would lead to the drainage direction of head distribution in the confined water key layer turn to the layered roadway, layered coal pillar and goaf. Then, they both transferred to form the layered water-rich roadway, layered water-containing coal pillar and water-containing goaf, respectively. Under the combined actions of the osmotic pressure Pg, the dynamic load σd and the mining stress σc, the significant high-stress concentration zone (energy accumulation, U2) was easily formed in the surrounding rocks of the water-rich layered roadway (target source), which was the second starting source. Meanwhile, “Energy Flow” 2 (“Transverse Flow”) was generated. Under the comprehensive impacts of two “Energy Flow” chains of “Transverse Flow” and “Longitudinal Flow”, the phenomenon of rock burst occurred in a water-rich layered roadway. Among them, the occurrence of a sufficient condition of rock burst in a water-rich layered roadway was as follows:
The combined actions of osmotic pressure
Pg, dynamic load
σd and mining stress
σc that led to the superimposed stress
σs in the surrounding rocks of a layered roadway (target source) was prone to form a significant high-stress concentration zone, which must be higher than the critical stress
σ0. The specific expression of superimposed stress
σs was as follows:
The specific expression of the sufficient condition of rock burst in a water-rich layered roadway was as follows:
The specific expression of
U2 accumulated by the corresponding “Transverse Flow” was as follows:
where
Pg was the osmotic pressure,
σd was the dynamic load,
σc was the mining stress and
E was the elastic modulus.
Therefore, the specific energy discrimination expression of rock burst in a water-rich layered roadway was as follows:
Meanwhile, according to the concept of an “Energy Flow” chain [
48], the discriminant converted into the form of “Energy Flow” was as follows:
Namely, at the identical and very short time interval, when the energy gradient tensor of “Longitudinal Flow” and the energy gradient tensor of “Transverse Flow” jointly impacted the layered water-rich roadway directionally, and the sum of them exceeded the critical energy gradient tensor when the rock burst occurred, the phenomenon of rock burst would occur in the layered water-rich roadway, where was the directional vector of “Energy Flow”.
Therefore, the essential attributes of the prevention and control of rock burst could be summarized from three aspects.
Firstly, controlling the starting source; for example, the protective layer mining technology, cantilever beam top-cutting releasing-energy technology and deep cavitation absorbing-energy technology were used to weaken the starting source. Then, it could greatly reduce the potential of rock burst.
Second, cutting off the propagation path (cutting off the “Energy Flow” chain), such as the hydraulic fracturing (cutting off the “Longitudinal Flow” chain) and large-diameter drilling (cutting off the “Transverse Flow” chain), could also greatly reduce the potential of rock burst.
Finally, protecting the target source (namely, adding one or more safety covers to the target source), such as the grouting modification technology on the surrounding rocks, roadway layout technology under goaf, protective layer mining technology, etc., could better reduce the possibility of rock burst.
In conclusion, “Energy Flow”, as a directional medium, could better reveal the occurrence mechanism of rock burst in a water-rich roadway. No matter the laboratory scale or engineering scale, it could provide a new perspective and concept for the effective prevention and control for rock burst in deep engineering.
5. Conclusions
(1) When β ∈ (0°, 30°) and β ∈ (45°, 90°), the influenced sensitivity of layered angles on the peak strength under TLUTP was weaker than that of confining pressure. When β ∈ (30°, 45°), the influence of layered angles on the peak strength was more sensitive than that of confining pressure.
(2) The nonlinear evolution relationships between the peak total input energy, peak elastic energy and peak dissipated energy, effective for confining pressure and layered angles under CTLTP and TLUTP, were obtained. Their specific expression was as follows: .
(3) At an identical and very short time interval, when the energy gradient tensor of “Longitudinal Flow”
and the energy gradient tensor of “Transverse Flow”
jointly impacted the layered water-rich roadway directionally, and the sum of them exceeded the critical energy gradient tensor
when the rock burst occurred, the phenomenon of rock burst would occur in the layered water-rich roadway. Namely, the specific discriminant formula with the form of “Energy Flow” of rock burst in a layered water-rich roadway was as follows:
(4) The nature of the effective prevention and accurate control of rock burst was obtained from three aspects, which corresponded to controlling the starting source, cutting off the “Energy Flow” chain (“Transverse Flow” and “Longitudinal Flow”) and protecting the target source.