Next Article in Journal
Hybrid Interfacial Transition Zone between Wet—On—Wet Casted Concrete—Microstructure and Mechanical Properties
Previous Article in Journal
Controlling the Oxygen Defects Concentration in a Pure BiFeO3 Bulk Ceramic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical and Field Investigations of Acoustic Emission Laws of Coal Fracture under Hydro-Mechanical Coupling Loading

1
Key Laboratory of Deep Coal Resource Mining, Ministry of Education, School of Mines, China University of Mining and Technology, Xuzhou 221116, China
2
Xiaoyun Coal Mine, Jining Energy Group, Jining 272000, China
*
Author to whom correspondence should be addressed.
Submission received: 28 July 2022 / Revised: 25 August 2022 / Accepted: 26 August 2022 / Published: 20 September 2022
(This article belongs to the Topic Fluid Mechanics)

Abstract

:

Abstract

Taking coal under hydro-mechanical coupling as the research object, the discrete element software PFC3D (particle flow code) was used to analyze the relationships among the force, acoustic emission (AE), and energy during coal fracture. Based on the moment tensor (MT) inversion, we revealed the AE event distribution and source type during crack initiation and propagation until the final failure of coal. Meanwhile, we examined the relationships among the stress, number and type of cracks, magnitude, KE, and b value of AE under different water and confining pressures. The results show that the numerical simulation can effectively determine the microscopic damage mechanism of coal under different conditions. Moreover, the rupture type of the numerical simulation is consistent with the field investigations, which verifies the rationality of the simulation. These research results can provide reference for safety production evaluation of water inrush mines.

Highlights

  • The relationships among stress, acoustic emission (AE), and energy during coal fracture under hydro-mechanical coupling loading were analyzed.
  • Moment tensor reveals event distribution, source type change, and b value of AE during crack initiation and propagation in coal.
  • The relationships among stress, number and type of cracks, AE, number of contacts and KE were revealed under different water and confining pressures.

1. Introduction

With the continuous increase in the mining depth of coal mines, geological conditions have become particularly complex [1]. Discontinuous joints such as microcracks in a rock stratum provide necessary locations for the storage and migration of groundwater [2,3]. In the seepage process, the generated water stress affects the stress field of a coal and rock mass [4]. Under the superposition of multiple stress fields, the frequency and intensity of catastrophic accidents such as rock bursts and water inrush remarkably increase [5,6]. When a water inrush accident occurs, the overall strength and plastic deformation of a coal pillar significantly change, which has a notable impact on the support capacity and resistance of the roadway confining pressure system, thereby affecting the width of the coal pillar [7].
Studies on the mechanical properties of coal and rock under hydro-mechanical coupling have mainly focused on the interaction between the seepage and stress fields, and the research methods have primarily included laboratory experiments and numerical simulations. For example, Hui et al. [8] conducted a series of laboratory acoustic emission (AE) tests to study the influence of fluid-filled rock joints on the propagation and attenuation of stress waves. Gardner et al. [9] investigated the interaction relationship between a rock mass and a fluid based on the three-dimensional discrete element and lattice Boltzmann methods. For the purpose of improving the calculation efficiency and accuracy, an extensive discussion on aspects that still need to be optimized was conducted. Chen et al. [10] studied the characteristics of the permeability evolution during sandstone failure under triaxial compression by experiments and simulations. Cai et al. [11] analyzed the evolution of the crack growth, coordination number, volumetric strain, and permeability under hydro-mechanical coupling in sandstone and established a synergistic mechanism among the microstructure, accumulated damage, and permeability evolution. Liu et al. [12] analyzed various parametric changes of a rock mass with hydro-mechanical coupling under unloading conditions. Wang et al. [13] investigated the mechanical characteristics and permeability evolution of red sandstone under hydro-mechanical coupling. Yao et al. [14] analyzed the mechanical properties and failure precursor characteristics of a coal and rock system caused by water content using AE and infrared radiation. Wang et al. [15] studied the deformation characteristics of granite under hydro-mechanical coupling and analyzed the changes in the permeability before and after expansion under different confining pressures. Zhang et al. [16] studied the influence of water and confining pressures on the permeability of fractured sandstone. Many scholars have studied the changes of rock mechanical parameters under the hydro-mechanical coupling loading in the laboratory, but there are few studies on the influence of seepage under different water and confining pressure [17,18]. In addition, few scholars have focused on the source rupture type in the process of rock seepage through AE.
In recent years, AE has been commonly used in the analysis of coal and rock failure. As an accompanying phenomenon in this process, it contains substantial information about the internal damage process of rocks, which helps to understand the mechanism of rock failure [19,20,21]. For instance, Shimizu et al. [22] discussed the effects of different fluid viscosity and particle size distributions on fracture propagation. Jia et al. [23] analyzed the relationship between AE and spatial fractal dimensions under different water soaking times. Ye et al. [24] obtained the relationships among AE events and the joint angle, water injection pressure, stress, and deformation amount in the process of a rock slip instability. Makhnenko et al. [25] observed that AE activity coincides with the onset of inelastic response in a fluid-saturated rock and explored the relationship between the clustering of AE events and inhomogeneous deformation. Li et al. [26] applied AE to dynamic disaster monitoring of coal and rock systems. Li et al. [27] concluded that there is a good positive correlation between the change in coal AE events and the water pressure curve under the action of triaxial hydraulic coupling. Chitrala et al. [28] reported the changes in the AE during hydraulic pressure damage of sandstone under different applied stresses. Chen et al. [29] used AE to analyze the evolution characteristics of granite damage during uniaxial compression and established a relationship between the permeability and the confining pressure. Makoto et al. [30] studied the relationship between the AE and fluid pressure during the failure of the Eagle Ford Shale, and analyzed the main types of AE. Lu et al. [31] successfully reproduced the stress redistribution, AE event evolution, permeability change, and formation of water inrush channels during a mining process. Jiang et al. [32] analyzed the AE event, energy, peak frequency, and crack type under the action of hydraulic fracturing in a layered rock. Most scholars widely use AE to locate the crack position of rock, and then study the force of the specimen [33]. However, there are few studies on the focal mechanism of AE, and analysis of the rupture form of rock/coal mass under different stress conditions is still lacking.
In summary, most studies mainly focused on the basic mechanical parameters of coal and rock masses under the action of hydro-mechanical coupling. However, they did not conduct detailed investigations on the focal mechanism of this action, and rarely used an energy method to evaluate the impact tendency of coal under hydro-mechanism coupling [34,35]. To better understand the characteristics of AE and its impact tendency under hydro-mechanism coupling, this study used the PFC3D software to analyze the relationships among the force, AE distribution (type), energy, and b value in the failure process of a coal sample under different water pressures. In addition, it is consistent with the results of our field investigation. These results are highly significant for guiding the safe production of coal mines.

2. Damage Mechanism of Coal under Hydro-Mechanism Coupling

To ensure no shear failure occurs during the loading of a rock, the limit value of the rock force can be obtained by the Mohr–Coulomb criterion [36] as follows:
τ = C + γ s σ
where τ is the shear stress of the rock, MPa, C is the internal cohesion, Pa, and γ s is the internal friction coefficient, which can be also expressed as
γ s = tan φ
It is assumed that the minimum principal stress on the rock is σ 3 . When the rock is in the critical limit state, the maximum principal stress, σ 1 , can be calculated. Owing to the existence of water pressure in the rock pores, a part of the confining pressure and vertical stress is offset; therefore, the water pressure eliminates the confining pressure (reducing the servo stress). The triaxial stress state of the rock under water pressure is shown in Figure 1.
When the Mohr–Coulomb circle is tangent to the envelope curve, the length of the right-angle side of the grey triangle area is λ 1 σ 1 λ 2 σ 3 / 2 , and the hypotenuse length is C / tan φ + λ 1 σ 1 + λ 2 σ 3 / 2 . Based on trigonometric functions,
λ 1 σ 1 λ 2 σ 3 / 2 = C / tan φ + λ 1 σ 1 + λ 2 σ 3 / 2 sin φ
Combined with Equation (2), the maximum principal stress under the critical stress state of the sample can be expressed as
σ 1 = 2 C + λ 2 γ s + γ s 2 + 1 σ 3 λ 1 γ s 2 + 1 γ s
In the experiments, when the internal cohesion and the friction angle have certain values, because λ 1 and λ 2 are less than 1, the stress–strain curve of the sample under water pressure is above that of the sample without water pressure; however, the peak stress in the former is smaller. This is attributed to the effect of the pore water pressure reducing the effective average stress of the rock, making it easier for the rock to reach the ultimate strength when the deviatoric stress remains constant.

3. Modeling Methodology

3.1. Engineering Background

A water inrush accident occurred in the Xiaoyun coal mine in Shandong province on 11 September 2018. The coal body of the 1318 working face was severely affected by the water pressure of proximately 3 MPa, and the soaking time was 100 days. A sample was taken from the soaking water coal body of the above working face, which has a simple structure and is relatively stable. To make the sample moisture non-volatile, it was stored in plastic wrap before experiments and processed into a cylinder with a diameter of 50 mm and a length of 100 mm. An MTS Landmark 370.50 rock testing machine was used for loading.

3.2. Model Setup

In the discrete element model, the core of the fluid flow algorithm is the topological structure formed by pipe domains. As shown in Figure 2, many domains that can store the pore water pressure are formed between the particles, and these domains are represented by red squares. A fracture between adjacent particles is regarded as a seepage pipe to simulate the flow of water. A pipeline is represented by a yellow line. When there is a pressure difference between adjacent domains, the fluid flows through the pipe to them.
To study the mechanical characteristics and instability mechanism of coal under the action of seepage, a model was set up, which is shown in Figure 2. For the model, the standard size (Φ50 mm × 100 mm) was adopted, and a total of 3117 particles were generated. The radius was 2–2.5 mm. The linear contact bond model was used between particles, and the servo mechanism was used to apply different confining pressures around the model. When the confining pressure reached the target stress, the particles in contact with the wall were fixed. A high-pressure water area was applied within the range of 0–10 mm in the z-axis direction at the bottom of the model. Under the action of the pressure difference, a water body seeped from the bottom to the top of the model. At this time, the water pressure is close to linear distribution along the longitudinal direction, which can be regarded as the equilibrium state of hydro-mechanical coupling. After the balance was reached, the upper and lower walls were controlled to apply an axial pressure at a speed of 0.03 m/s, and the loading was stopped when the residual stress reached 70% peak stress.
The model parameter calibration was based on the uniaxial compressive strength, strain, and failure form of the specimen. Based on Figure 3, the peak stresses from the experiments and numerical simulations are 6.09 MPa and 6.05 MPa, respectively, and the corresponding strains are 1.93% and 1.95%. The difference between them is relatively small, and the failure modes are also similar. From the variation trends of the AE events, the microscopic properties inside the coal body are changed owing to the softening effect of water. When the strain is approximately 5%, the number of AE events shows remarkable fluctuations, which is followed by a quiet period. When the strain reaches 1.6%, the AE events abruptly increase, subsequently cracks gradually penetrate to form macroscopic cracks, and, finally, the number of AE reaches the maximum. Based on the above analysis, the rationality of the microparameter values can be verified, which are listed in Table 1.

3.3. MT Calculation in the PFC

In seismology, MT is a mathematical representation of fault motion and an important tool for source characterization. We use M to represent the source MT, which is a 3 × 3 symmetric matrix with nine couples and six independent components, as shown in Figure 4.
Figure 4 shows the force couples in the Cartesian coordinate system, which is related to source strength and fault direction. We can further decompose the MT into elementary sources, e.g., compression, shear, and tensile. Then, we can intuitively obtain the rupture type of A.E. Hudson’s source type plot [37,38], which is a way to vividly represent decomposition of an MT into isotropic, compensated linear vector dipole (CLVD) and double-couple (DC) components. The MT decomposition is shown in the following Equations (5) and (6):
M = M D C + M i s o + M v c l v d = M 1 ϕ S + M 2 ϕ D + M i s o + M v c l v d
M 11 M 12 M 13 M 21 M 22 M 23 M 31 M 32 M 33 = 1 2 M 22 M 11 M 12 0 M 21 1 2 M 22 M 11 0 0 0 0 + 0 0 M 13 0 0 M 23 M 31 M 32 0 + 1 3 M 11 + M 22 + M 33 1 0 0 0 1 0 0 0 1 + 1 3 1 2 M 11 + M 22 M 33 1 0 0 0 1 0 0 0 2
The PFC 5.0 software (Itasca Consulting Group, Inc., Minneapolis, MN, USA) can directly obtain the force magnitude, direction, and displacement of particles, and then obtain the MT according to the definition of AE. The calculation process of the AE in PFC is shown in Figure 5.
An AE event is composed of particles around microcracks in PFC 5.0 software. When the contact between two particles is broken, a crack will occur. Two particles in contact with the crack and particles in contact with these two particles are marked as an AE event [39,40]. From Figure 5, (a) is the AE start stage, (b) is the AE end stage, and the range delineated by the red circle is the AE event. Through MT decomposition, tensile, shear, and compression failure mode can be obtained. The AE rupture mode corresponding to PFC 5.0 numerical simulation software is shown in Figure 6.
Similar to earthquakes, the intensity of AE can also be expressed by magnitude. In PFC software, the AE intensity is calculated according to the strain energy (Ek) change before and after the event [41]. For each step of calculation, the Ek of all particles in the AE event can be written as Equation (7):
E k = i = 1 n 1 2 F n l 2 k n + F s l 2 k s
where F n l , F n s , k n , and k s are normal stress, shear stress, normal stiffness, and shear stiffness, respectively; n is the number of contacts within the AE range. For an AE event, its magnitude can be measured by the change of Ek at the beginning and end of AE [42]:
M e = 2 3 log Δ E k 6.0
Based on the knowledge of seismology, the relationship between the magnitude and the frequency is
log N = a b M e
where Me is the magnitude, N is the frequency, and a and b are constants. The b value can be used to evaluate the damage degree of a rock mass and the size of a crack. A large b value implies a large proportion of small-sized cracks in the rock mass and a slow initiation and development of new cracks.

4. Mechanical Characteristics and Failure Mechanism of Water-Soaked Coal

4.1. AE and Stress Changes during Failure

In the simulations, the AE and the energy were monitored to investigate the characteristics of coal failure. To study the coal failure mechanism under hydro-mechanical coupling, the simulations described in Section 4.1 were taken as an example, and the confining and water pressures were set as 1 MPa and 2.5 MPa, respectively. The microscopic characteristics of the AE distribution (type) and energy change during the failure process were analyzed in detail.

4.1.1. Relationship between AE and Stress

To study the coal failure characteristics, the fish language and the history module were used to monitor the vertical stress, strain, AE, energy, etc. Simultaneously, the monitoring circle module was used to record the parameter, and the results are shown in Figure 7 and Figure 8.
The simulated stress–strain curve, cracks, and locations of the AE sources during uniaxial compression of water-soaked coal are shown in Figure 7. The coal failure process can be divided into four stages. The first stage is the compaction stage, in which the stress gradually increases without cracks. The second stage is the crack initiation stage, and the crack initiation stress is 5.33 MPa. The crack distribution is mainly observed in the lower part of the model, and the increase in the number of shear cracks is significantly greater than that in the number of tensile cracks. The third stage shows a notable increase in the cracks. In this stage, the increase rate of the shear cracks is much larger than that of the tensile cracks, and microcracks gradually develop upward. However, the damage degree of the lower part of the entire sample is much greater than that of the upper part, and the peak stress reaches 7.8 MPa. The fourth stage is the post-peak stage, in which the number of cracks sharply increases, the cracks gradually penetrate, and, finally, the entire instability occurs.
By depicting the distribution of the micro contact force in the particle system, the macro mechanical mechanism can be revealed. Figure 8 shows the contact force between particles in the model, which can better reflect the force in different directions of the model. From Figure 8, as the axial stress increases, the internal contact force shows a gradually increasing trend. The shape of the contact force gradually evolves from “spherical” to “elliptical”, the stress levels at different positions gradually change, and at the final stage of failure, the local stress abruptly increases. The vertical stresses at monitoring points F1–F4 are consistent with the stress–strain curve. Except for at point F1, closeness to the water injection area implies a high vertical stress. Based on the change in the yy-stress, when the strain is 1.8%, the horizontal stress remarkably changes, and the model shows an instability at this stage.

4.1.2. AE Event Distribution and Field Verification

In the loading process, a total of 509 AE events were recorded. Figure 9 shows the T–K diagram of the MT obtained using the standard inversion method.
From Figure 9, at the four stages of loading, total of 15, 77, 65, and 417 AE events are recorded, respectively. The distribution of the AE events is relatively scattered, and most of them are distributed in the shear area. At stages (a) and (b), 15 and 77 events are generated, respectively, with magnitudes of between −6.68 and −5.3. The corresponding stresses at the end of these two stages are 5.32 MPa and 7.4 MPa, which are approximately 68% and 94% of peak stress, respectively. The AE events during this period are mainly dominated by the water pressure. At stages (c) and (d), the axial stress gradually reaches the peak strength of the coal. During this period, the effect of the axial stress is far greater than that of the water pressure, and the number of AE events sharply increases (mainly shear failure), which leads to instability of the sample.

4.1.3. Proportions of Source Type

To study the types of AE events in coal under the action of hydro-mechanical coupling, the number and proportion of tensile, shear, and implosion source types at different failure stages were quantitatively counted. The results are shown in Figure 10.
From Figure 10, as the axial stress increases, the numbers of implosions and shears sharply increase; however, their proportions first increase and subsequently decrease. The entire damage is dominated by shear cracks. In stages (a) and (b), the proportions of shear and tension source types become larger, from 69% and 25% to 70% and 29%, respectively. Tensile failure increases in a slightly larger proportion than shear cracks. Observing the AE distribution shows that the damage is mainly concentrated in the lower part of the model. The particles in the water pressure area are subjected to tension. Simultaneously, tensile failure is prone to occur under axial stress. With the increase in the axial stress, the influence of the water pressure gradually weakens. Therefore, in stages (c) and (d), the proportion of tensile source type slightly decreases, whereas that of implosion source type increases to 8%. In summary, the water pressure mainly affects the early stage of loading and is primarily influenced by the axial stress in the later failure process.

4.1.4. Energy Change

Coal failure process is accompanied by generation, accumulation, dissipation, and release of energy, and the evolution of energy is closely consistent with the expansion of internal cracks in coal. Using the energy module of the PFC3D, the energy variation was recorded during coal failure process.
From Figure 11, the kinetic energy tends to increase in the initial stage of loading. The main cause is the force unbalance between some particles under the action of the axial pressure, which destroys the linear parallel bonds between them. There is no increase in the slip energy within a certain period. When the vertical stress of the coal body reaches approximately 50% peak stress, the kinetic and slip energies gradually increase. When the strain reaches 2%, the model slips overall. At this time, the slip and kinetic energies gradually reach maximum. The total strain energy is consistent with the overall change trend of the stress–strain curve.

4.2. Effects of Water and Confining Pressures on Coal

The confining and water pressures were set as 0.5, 1.5, and 2.5 MPa, respectively. The specific settings are summarized in Table 2.

4.2.1. Failure Characteristics

A total of 12 models were established, and each model number and its corresponding parameters are listed in Table 2. To analyze the internal damage of coal under the effects of different water and confining pressures, the AE events distribution of the different models were examined, and the results are shown in Figure 12.
From Figure 12, it is seen that the confining and water pressures significantly affect the failure pattern of coal. When the model is damaged, the AE events are mainly concentrated in the water pressure area. Based on the comparative analysis of P-1–P-4, P-5–P-8, and P-9–P-12, as the water pressure increases with the same confining pressure, the internal damage and high-energy events increase, and most of the AE events are mainly concentrated in the lower part of the model. When the water pressure is kept constant, the model accumulates more elastic energy before it reaches failure with increasing confining pressure. When the model is damaged, more high-energy events are generated.
From Figure 13, the force of the model has a large anisotropy under the effects of the confining, water, and axial pressures. As the water pressure increases, the number of contacts in the lower part is much smaller than that in the upper part, mainly owing to the greater damage in the lower part. As the confining pressure increases, the number of contacts increases significantly. The maximum number of contacts is increased from 28 to 29, and the influence of the water pressure is gradually weakened, indicating that the strengthening of a support can help protect the stability of the roadway.

4.2.2. Relationship between AE and Stress–Strain Curve

To analyze the effects of different water and confining pressures on the overall strength of the model, the stress–strain curves of the P-1–P-12 models were examined, and are shown in Figure 14.
Based on Figure 14, when the water pressure remains constant, a large confining pressure implies a high peak strength of the sample and a large deformation before failure. When the confining pressure is constant, a large water pressure implies a low peak strength and a small axial deformation before the peak; however, the slope increases, which is consistent with the analysis results of Equation (3).
Figure 15 shows the changes in the numbers of shear and tensile cracks during the failure of different models. The number of shear cracks is approximately twice that of tensile cracks. As the water pressure increases before uniaxial compression, more cracks are generated during the seepage process, and the number of shear cracks is greater than that of tensile cracks. In the uniaxial compression process, when the confining pressure is small, the water pressure has a significant effect on the number of AE events (NP-2 ≈ 1.23NP-3 ≈ 2.01NP-4). When the confining pressure is high, the effect of the water pressure significantly weakens (NP-10 ≈ 1.05NP-11 ≈ 1.4NP-12), and the degree of weakening significantly strengthens.

4.2.3. Frequency–Magnitude Curve and b Value

Smaller earthquakes occur more frequently than large earthquakes [43]. This trend can be expressed by the magnitude–frequency relationship and evaluated by b value.
Figure 16 shows that there is a good linear relationship between the water and confining pressures and the b value. If the confining and water pressures are high, the value of b is small, which proves that the large-scale cracks inside the model gradually increase and rapidly and unsteadily expand, and the AE distribution becomes more diffusing. Therefore, as the water and confining pressures increase, the internal stress concentration of the model rises. When cracks occur inside the model, microcracks are more likely to grow and penetrate large-scale cracks. Simultaneously, a high confining pressure implies a large number of AE events inside the rock mass. Thus, the water and confining pressures have remarkable impacts on the AE of the coal body.

4.3. Burst Tendency of Coal under Hydro-Mechanism Coupling

There is a certain relationship between the dissipated and strain energies during the uniaxial compression of coal, which is shown in Figure 17.
A schematic of the energy burst tendency is shown in Figure 17, where F1 represents the difference between the elastic energy accumulated before the peak and the dissipated energy consumed during the plastic deformation and crack development of the coal. F2 represents the energy consumed when the coal is damaged after the peak, i.e., F1/F2 is an energy burst tendency indicator of the coal. In the energy module part of the PFC3D 5.0 software (Itasca Consulting Group, Inc., Minneapolis, MN, USA), the total strain energy includes the strain energy ( E c ) and the bond strain energy ( E p b ), and the total dissipated energy includes the frictional energy ( E μ ) and the kinetic energy ( E k ). Among them, the slip energy remains equal to zero until the bond breaks. It can be seen that the particle flow can monitor the change in the energy in the deformation and failure process of the model in real time. Therefore, in the full stress–strain curve, the elastic and dissipation energies at any time can be expressed in PFC3D as
U i e = E p b + E c U i d = E f + E k
After deducting the energy consumed by the plastic deformation of the rock and crack development, the burst tendency of the coal body can be determined using the energy method as follows:
F 1 = U p e a k e + U p e a k d F 2 = U p e a k e + U r e s i d u a l e K E = F 1 / F 2
Figure 18 shows the calculation results of the burst tendency of coal under different confining and water pressures. When the water and confining pressures are 0 MPa and 1.5 MPa, respectively, the coal body has the lowest burst tendency, and KE is 2.4. When the water and confining pressures are 1.5 MPa and 0 MPa, respectively, the coal body has the highest burst tendency, and KE is 2.67. From the KE distribution, when the water pressure remains constant, as the confining pressure increases, KE becomes smaller. When the confining pressure remains constant, as the water pressure increases, KE becomes larger. Therefore, the burst tendency under different conditions can be accurately determined.

5. Onsite Investigations

The 1314 working face is located at the level of −790 m; the mining seam is the #3 coal with the thickness of 1.5–3.7 m, and the buried depth is 648–700 m. The roadway is arranged along the roof of the #3 coal seam. The coal seam structure of the working face is simple, with an inclination angle of 13°–26°.
On 11 September 2018, a water inrush accident occurred in this working face, and the coal was immersed in a water pressure of approximately 3 MPa for approximately 100 days, which had a certain impact on the stability of the roadway under the influence of water. By counting MS events during soaking, the Hudson’s source type is obtained according to focal mechanism inversion, as shown in Figure 19. There are a wide range of failure types in the coal, mainly shear and tensile cracks, and its failure mainly occurs in the type of shear cracks. The field observation results are consistent with the simulation conclusions.
In terms of forces, a shear crack can be expressed by two perpendicular force dipoles with zero angular momentum. Therefore, the shear failure source is usually called DC. In order to intuitively show the internal damage form of the coal, several MS events around the working face were selected for detailed analysis. Three types of beachballs, full, deviatoric, and DC, were drawn, as shown in Figure 20. It is found that the shear component accounts for a large proportion in the above MS events, and it is concluded that the force of the shear component is the main factor leading to the failure of the coal. In summary, the failure of coal under water pressure is mainly dominated by shear source, mixed with a small amount of other components, which is strictly consistent with the simulation conclusions. From Figure 21, several fault plane solutions were solved, including two normal fault, three reverse fault, and one normal oblique slip, most of which face the interior of the working face.

6. Conclusions

(1)
Affected by hydro-mechanical coupling, the damage degree of coal in the flooded part is relatively large, and its damage mainly occurs as shear cracks; this is verified by field observation. The closeness of coal to the water injection area implies a high vertical stress.
(2)
Monitoring of the energy changes shows that the kinetic and slip energies increase slightly in the initial stage of loading. When the peak stress is greater than 50%, a sharp energy increase occurs, and the increase rate becomes increasingly high.
(3)
There is a good linear relationship between the water pressure, confining pressure, and b value. With increasing water and confining pressures, the damage degree and AE energy inside the model increase, the value of b reduces, and most of the AE events are mainly concentrated in the water injection area. When the confining pressure is low, the water pressure has a significant effect on the number of AE events; otherwise, the effect of water pressure gradually weakens.
(4)
When the water pressure is constant, as the confining pressure increases, KE decreases. When the confining pressure is kept constant, as the water pressure increases, KE increases.

Author Contributions

Conceptualization, Y.-M.W.; Data curation, J.-H.W.; Investigation, Z.-W.Z. and H.-F.C.; Methodology, C.-P.L.; Software, J.-F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (51574225).

Informed Consent Statement

Written informed consent has been obtained from the patient(s) to publish this paper.

Data Availability Statement

The data in this manuscript are available from the authors.

Acknowledgments

We gratefully wish to acknowledge the collaborative funding support from the National Natural Science Foundation of China (51574225). In addition, the data in this manuscript are available from the authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

  1. Ma, D.; Miao, X.; Bai, H.; Huang, J.; Pu, H.; Wu, Y.; Zhang, G.; Li, J. Effect of mining on shear sidewall groundwater inrush hazard caused by seepage instability of the penetrated karst collapse pillar. Nat. Hazards. 2016, 82, 73–93. [Google Scholar] [CrossRef]
  2. Yu, Y.; Cheng, W.; Zhou, G. Numerical simulation study on permeability properties of coal seam during high pressure water injection process. Electron. J. Geotech. Eng. 2014, 19, 10315–10326. [Google Scholar]
  3. Kou, M.M.; Liu, X.R.; Wang, Z.Q.; Tang, S.D. Laboratory investigations on failure, energy and permeability evolution of fissured rock-like materials under seepage pressures. Eng. Fract. Mech. 2021, 247, 107694. [Google Scholar] [CrossRef]
  4. Grgic, D.; Amitrano, D. Creep of a porous rock and associated acoustic emission under different hydrous conditions. J. Geophys. Res. Solid Earth 2009, 114, B10201. [Google Scholar] [CrossRef]
  5. Li, H.; Li, H.; Gao, B.; Jiang, D.; Feng, J. Study of Acoustic Emission and Mechanical Characteristics of Coal Samples under Different Loading Rates. Shock. Vib. 2015, 2015, 1–11. [Google Scholar] [CrossRef]
  6. Wen, Z.; Wang, X.; Chen, L.; Lin, G.; Zhang, H. Size Effect on Acoustic Emission Characteristics of Coal-Rock Damage Evolution. Adv. Mater. Sci. Eng. 2017, 2017, 1–8. [Google Scholar] [CrossRef]
  7. Chen, X.; Yu, J.; Tang, C.A.; Li, H.; Wang, S. Experimental and Numerical Investigation of Permeability Evolution with Damage of Sandstone Under Triaxial Compression. Rock Mech. Rock Eng. 2017, 50, 1529–1549. [Google Scholar] [CrossRef]
  8. Hui, Y.; Hfda, B.; Jz, C. Effects of filling fluid type and composition and joint orientation on acoustic wave propagation across individual fluid-filled rock joints. Int. J. Rock Mech. Min. Sci. 2020, 128, 104248. [Google Scholar]
  9. Gardner, M.; Sitar, N. Modeling of Dynamic Rock–Fluid Interaction Using Coupled 3-D Discrete Element and Lattice Boltzmann Methods. Rock Mech. Rock Eng. 2019, 52, 5161–5180. [Google Scholar] [CrossRef]
  10. Chen, Z.; Wang, M. Pore-scale modeling of hydromechanical coupled mechanics in hydrofracturing process. J. Geophys. Res. Solid Earth 2017, 122, 3410–3429. [Google Scholar] [CrossRef]
  11. Cai, Y.; Chen, X.; Yu, J.; Zhou, J. Numerical Study on the Evolution of Mesoscopic Properties and Permeability in Sandstone Under Hydromechanical Coupling Conditions Involving Industrial Internet of Things. IEEE Access 2018, 6, 11804–11815. [Google Scholar] [CrossRef]
  12. Liu, S.; Zhu, Q.; Shao, J. Deformation and mechanical properties of rock: Effect of hydromechanical coupling under unloading conditions. Bull. Eng. Geol. Environ. 2020, 79, 1–18. [Google Scholar] [CrossRef]
  13. Wang, Q.; Hu, X.; Zheng, W.; Li, L.; Zhou, C.; Ying, C.; Xu, C. Mechanical Properties and Permeability Evolution of Red Sandstone Subjected to Hydro-mechanical Coupling: Experiment and Discrete Element Modelling. Rock Mech. Rock Eng. 2021, 54, 2403–2405. [Google Scholar] [CrossRef]
  14. Yao, Q.; Wang, W.; Zhu, L.; Xia, Z.; Tang, C.; Wang, X. Effects of moisture conditions on mechanical properties and AE and IR characteristics in coal–rock combinations. Arab. J. Geosci. 2020, 13, 1–15. [Google Scholar] [CrossRef]
  15. Wang, Z.H.; Ren, W.G.; Tan, Y.L.; Konietzky, H. Experimental and Numerical Study on Hydromechanical Coupled Deformation Behavior of Beishan Granite considering Permeability Evolution. Geofluids 2020, 2020, 1–14. [Google Scholar] [CrossRef]
  16. Zhang, S.; Zhang, D.; Wang, Z.; Chen, M. Influence of Stress and Water Pressure on the Permeability of Fissured Sandstone Under Hydromechanical Coupling. Mine Water Environ. 2018, 37, 774–785. [Google Scholar] [CrossRef]
  17. Martyushev, D.A.; Galkin, S.V.; Shelepov, V.V. The Influence of the Rock Stress State on Matrix and Fracture Permeability under Conditions of Various Lithofacial Zones of the Tournaisian–Fammenian Oil Fields in the Upper Kama Region. Mosc. Univ. Geol. Bull. 2019, 74, 573–581. [Google Scholar] [CrossRef]
  18. Riabokon, E.; Turbakov, M.; Popov, N.; Kozhevnikov, E.; Poplygin, V.; Guzev, M. Study of the Influence of Nonlinear Dynamic Loads on Elastic Modulus of Carbonate Reservoir Rocks. Energies 2021, 14, 8559. [Google Scholar] [CrossRef]
  19. Zhang, Q.; Wang, L. Research Status and Prospect of Rock Acoustic Emission Technology. Appl. Mech. Mater. 2012, 170, 461–464. [Google Scholar] [CrossRef]
  20. Zhang, Q.; Wang, L. Research Status and Prospect of Rock Mass Acoustic Emission Instrument. Adv. Mater. Res. 2013, 671, 280–283. [Google Scholar] [CrossRef]
  21. Zhang, Y.H.; Li, Y.; Zhang, Y. The Research of Shaoguo Tunnel Failure Characteristics of Triaxial and Acoustic Emission Characteristics. Appl. Mech. Mater. 2014, 501, 1757–1760. [Google Scholar] [CrossRef]
  22. Shimizu, H.; Murata, S.; Ishida, T. The distinct element analysis for hydraulic fracturing in hard rock considering fluid viscosity and particle size distribution. Int. J. Rock Mech. Min. Sci. 2011, 48, 712–727. [Google Scholar] [CrossRef]
  23. Jia, Z.; Ren, L.; Liu, Q.; Peng, Y.; Xu, D.; Zha, E. Influence of water-soaking time on the acoustic emission characteristics and spatial fractal dimensions of coal under uniaxial compression. Therm. Sci. 2017, 21, 327–334. [Google Scholar] [CrossRef]
  24. Ye, Z.; Ghassemi, A. Investigation of Microseismicity and Permeability Evolution in Shale Fractures during Stimulation. SPE Prod. Oper. 2020, 35, 0797–0808. [Google Scholar] [CrossRef]
  25. Makhnenko, R.Y.; Ge, C.; Labuz, J.F. Localization of deformation in fluid-saturated sandstone. Int. J. Rock Mech. Min. Sci. 2020, 134, 104455. [Google Scholar] [CrossRef]
  26. Li, J.G.; Liu, H. Application Conditions on Acoustic Emission (AE) Technique Monitoring Coal and Rock Dynamic Disasters in Mines. Adv. Mater. Res. 2012, 413, 235–240. [Google Scholar] [CrossRef]
  27. Li, N.; Sun, W.; Huang, B.; Chen, D.; Zhang, S.; Yan, M. Acoustic Emission Source Location Monitoring of Laboratory-Scale Hydraulic Fracturing of Coal Under True Triaxial Stress. Nat. Resour. Res. 2021, 30, 2297–2315. [Google Scholar] [CrossRef]
  28. Chitrala, Y.; Moreno, C.; Sondergeld, C.; Rai, C. An experimental investigation into hydraulic fracture propagation under different applied stresses in tight sands using acoustic emissions. J. Pet. Sci. Eng. 2013, 108, 151–161. [Google Scholar] [CrossRef]
  29. Chen, L.; Liu, J.F.; Wang, C.P. Characterization of damage evolution in granite under compressive stress condition and its effect on permeability. Int. J. Rock Mech. Min. Sci. 2014, 71, 340–349. [Google Scholar] [CrossRef]
  30. Makoto, N.; Chen, Y.; Yamamoto, K.; Morishige, Y.; Imakita, K.; Tsutumi, N.; Kawakata, H.; Ishida, T.; Tanaka, H.; Arima, Y.; et al. Tensile-dominant fractures observed in hydraulic fracturing laboratory experiment using eagle ford shale. Geophys. J. Int. 2020, 222, 769–780. [Google Scholar]
  31. Lu, Y.; Wang, L. Numerical simulation of mining-induced fracture evolution and water flow in coal seam floor above a confined aquifer. Comput. Geotech. 2015, 67, 157–171. [Google Scholar] [CrossRef]
  32. Jiang, Z.; Li, Q.; Hu, Q.; Liang, Y.; Xu, Y.; Liu, L.; Wu, X.; Li, X.; Wang, X.; Hu, L.; et al. Acoustic emission characteristics in hydraulic fracturing of stratified rocks: A laboratory study—ScienceDirect. Powder Technol. 2020, 371, 267–276. [Google Scholar] [CrossRef]
  33. Stanchits, S.; Surdi, A.; Gathogo, P.; Edelman, E.; Suarez-Rivera, R. Onset of Hydraulic Fracture Initiation Monitored by Acoustic Emission and Volumetric Deformation Measurements. Rock Mech. Rock Eng. 2014, 47, 1521–1532. [Google Scholar] [CrossRef]
  34. Li, B.; Liu, J.; Bian, K.; Ai, F.; Hu, X.; Chen, M.; Liu, Z. Experimental study on the mechanical properties weakening mechanism of siltstone with different water content. Arab. J. Geosci. 2019, 12, 1–14. [Google Scholar] [CrossRef]
  35. Chen, Y. Permeability Evolution in Granite Under Compressive Stress Condition. Geotech. Geol. Eng. 2017, 36, 641–647. [Google Scholar] [CrossRef]
  36. He, M.; Li, Q.; Li, X.; Xu, L.; Kühn, M. Hydromechanical behaviors of andesite under different stress states during fluid injection. J. Rock Mech. Geotech. Eng. 2021, 13, 727–744. [Google Scholar] [CrossRef]
  37. Aki, K.; Richards, P.G. Quantitative Seismology, 2nd ed.; University Science Books: Sausalito, CA, USA, 2002; ISBN 0-935702-96-2. [Google Scholar]
  38. Hudson, J.A.; Pearce, R.G.; Rogers, R.M. Source type plot for inversion of the moment tensor. J. Geophys. Res. 1989, 94, 765–774. [Google Scholar] [CrossRef]
  39. Zhou., Z.; Zhou, J.; Zhao, Y.; Chen, L.; Li, C. Microscopic Failure Mechanism Analysis of Rock Under Dynamic Brazilian Test Based on Acoustic Emission and Moment Tensor Simulation. Front. Phys. 2021, 8, 592483. [Google Scholar] [CrossRef]
  40. Zhao, Y.; Zhao, G.; Zhou, J.; Ma, J.; Cai, X. Failure mechanism analysis of rock in particle discrete element method simulation based on moment tensors. Comput. Geotech. 2021, 136, 104215. [Google Scholar] [CrossRef]
  41. Itasca Consulting Group Inc. PFC2D (Particle Flow Code in 2 Dimensions); ICG: Minneapolis, MN, USA, 1999. [Google Scholar]
  42. Hanks, T.C.; Kanamori, H. A moment magnitude scale. J. Geophys. Res. 1979, 84, 2348–2350. [Google Scholar] [CrossRef]
  43. Shearer, P.M. Introduction to Seismology; Birkhauser: Basel, Switzerland, 1979. [Google Scholar]
Figure 1. Schematic of critical stress state under influence of water pressure based on Mohr–Coulomb criterion.
Figure 1. Schematic of critical stress state under influence of water pressure based on Mohr–Coulomb criterion.
Materials 15 06510 g001
Figure 2. Schematic of incorporation of fluid coupling.
Figure 2. Schematic of incorporation of fluid coupling.
Materials 15 06510 g002
Figure 3. Variation curves of vertical stress and AE impact number of coal sample during uniaxial compression. The red and black lines represent the AE number and the stress-strain curve, respectively. The upper left corner figure is the failure figure of laboratory and simulated respectively.
Figure 3. Variation curves of vertical stress and AE impact number of coal sample during uniaxial compression. The red and black lines represent the AE number and the stress-strain curve, respectively. The upper left corner figure is the failure figure of laboratory and simulated respectively.
Materials 15 06510 g003
Figure 4. The nine force couples representing the components of MT.
Figure 4. The nine force couples representing the components of MT.
Materials 15 06510 g004
Figure 5. Generation process of an AE event: (a) is the start of an event and (b) is the finish state of this event.
Figure 5. Generation process of an AE event: (a) is the start of an event and (b) is the finish state of this event.
Materials 15 06510 g005
Figure 6. Schematic diagram of AE event failure mode in particle model.
Figure 6. Schematic diagram of AE event failure mode in particle model.
Materials 15 06510 g006
Figure 7. Stress–strain curve and number of cracks; upper snapshots show source distributions at different stages. I, II, III, and IV are the stages of crack generation, which are pore compaction stage, fracture expansion stage, peak stress stage and post peak stage, respectively.
Figure 7. Stress–strain curve and number of cracks; upper snapshots show source distributions at different stages. I, II, III, and IV are the stages of crack generation, which are pore compaction stage, fracture expansion stage, peak stress stage and post peak stage, respectively.
Materials 15 06510 g007
Figure 8. Changes in model internal stress and contact force. Panel (A) shows the internal contact force at different stages of the model. The warmer the color and the thicker the network, the greater the contact force. Panel (B) is the stress change at the monitoring point under the different stage corresponding to (A). (a–d) are the stages of crack generation, which are pore compaction stage, fracture expansion stage, peak stress stage and post peak stage, respectively.
Figure 8. Changes in model internal stress and contact force. Panel (A) shows the internal contact force at different stages of the model. The warmer the color and the thicker the network, the greater the contact force. Panel (B) is the stress change at the monitoring point under the different stage corresponding to (A). (a–d) are the stages of crack generation, which are pore compaction stage, fracture expansion stage, peak stress stage and post peak stage, respectively.
Materials 15 06510 g008
Figure 9. Hudson T–K locations of source models from (ad).
Figure 9. Hudson T–K locations of source models from (ad).
Materials 15 06510 g009
Figure 10. Nature of source versus time in uniaxial compression tests: (A) number and ratio of sources and (B) proportion of sources. (a–d) are the stages of crack generation, which are pore compaction stage, fracture expansion stage, peak stress stage and post peak stage, respectively.
Figure 10. Nature of source versus time in uniaxial compression tests: (A) number and ratio of sources and (B) proportion of sources. (a–d) are the stages of crack generation, which are pore compaction stage, fracture expansion stage, peak stress stage and post peak stage, respectively.
Materials 15 06510 g010
Figure 11. Various energy changes in coal failure process.
Figure 11. Various energy changes in coal failure process.
Materials 15 06510 g011
Figure 12. Source distributions of different models at different stages.
Figure 12. Source distributions of different models at different stages.
Materials 15 06510 g012aMaterials 15 06510 g012b
Figure 13. Number of contacts in different orientations.
Figure 13. Number of contacts in different orientations.
Materials 15 06510 g013
Figure 14. Stress–strain curves of different models.
Figure 14. Stress–strain curves of different models.
Materials 15 06510 g014
Figure 15. Development of cracks in coal body failure process in different models.
Figure 15. Development of cracks in coal body failure process in different models.
Materials 15 06510 g015aMaterials 15 06510 g015b
Figure 16. Frequency–amplitude curve and b value.
Figure 16. Frequency–amplitude curve and b value.
Materials 15 06510 g016aMaterials 15 06510 g016b
Figure 17. Relationship between dissipated and strain energies on stress–strain curve during uniaxial compression of coal. O, A and D correspond to the starting, peak stress and post peak value of the stress-strain curve respectively. B is the strain value of the model after unloading. C and D are the strain values corresponding to the peak and post peak stresses, respectively.
Figure 17. Relationship between dissipated and strain energies on stress–strain curve during uniaxial compression of coal. O, A and D correspond to the starting, peak stress and post peak value of the stress-strain curve respectively. B is the strain value of the model after unloading. C and D are the strain values corresponding to the peak and post peak stresses, respectively.
Materials 15 06510 g017
Figure 18. Change in KE under different water and confining pressures.
Figure 18. Change in KE under different water and confining pressures.
Materials 15 06510 g018
Figure 19. Hudson’s source type of coal failure under water pressure.
Figure 19. Hudson’s source type of coal failure under water pressure.
Materials 15 06510 g019
Figure 20. Beachballs in the 1314 working face.
Figure 20. Beachballs in the 1314 working face.
Materials 15 06510 g020
Figure 21. Fault plane solutions of the 1314 working face.
Figure 21. Fault plane solutions of the 1314 working face.
Materials 15 06510 g021
Table 1. Microparameters of numerical model.
Table 1. Microparameters of numerical model.
Mechanical Parameters of Coal
Particle radius/mm2.0–2.5Elastic modulus/GPa0.5
Density/kg·m31250Friction angle/°45
Cohesive strength/MPa4.5Porosity0.18
Parameters for Flow Model
Ap_zero/mm1.3 × 10−3Bulk_W/MPa2.2 × 103
Flow_perm/mm/s7.0 × 10−2Flow_dt/s1.0 × 10−3
P_give/MPa2.5Gap_mul0
Table 2. Setting of confining and water pressures.
Table 2. Setting of confining and water pressures.
ModelP-1P-2P-3P-4P-5P-6P-7P-8P-9P-10P-11P-12
Confining pressure/MPa0.50.50.50.511111.51.51.51.5
Water pressure /MPa00.51.52.500.51.52.500.51.52.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, J.-F.; Lu, C.-P.; Zhan, Z.-W.; Cui, H.-F.; Wang, Y.-M.; Wang, J.-H. Numerical and Field Investigations of Acoustic Emission Laws of Coal Fracture under Hydro-Mechanical Coupling Loading. Materials 2022, 15, 6510. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15196510

AMA Style

Song J-F, Lu C-P, Zhan Z-W, Cui H-F, Wang Y-M, Wang J-H. Numerical and Field Investigations of Acoustic Emission Laws of Coal Fracture under Hydro-Mechanical Coupling Loading. Materials. 2022; 15(19):6510. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15196510

Chicago/Turabian Style

Song, Jie-Fang, Cai-Ping Lu, Zhao-Wei Zhan, Hai-Feng Cui, Yan-Min Wang, and Jian-Hua Wang. 2022. "Numerical and Field Investigations of Acoustic Emission Laws of Coal Fracture under Hydro-Mechanical Coupling Loading" Materials 15, no. 19: 6510. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15196510

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop