Next Article in Journal
Enhanced Flotation Recovery of Fine Molybdenite Particles Using a Coal Tar-Based Collector
Next Article in Special Issue
Discrimination of Different AE and EMR Signals during Excavation of Coal Roadway Based on Wavelet Transform
Previous Article in Journal
Recent Advancements in Metallurgical Processing of Marine Minerals
Previous Article in Special Issue
Investigation on the Failure Mechanism of Weak Floors in Deep and High-Stress Roadway and the Corresponding Control Technology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Investigation on Influence of Two Combined Faults and Its Structure Features on Rock Burst Mechanism

1
School of Mines, China University of Mining & Technology, Xuzhou 221116, China
2
Jiangsu Engineering Laboratory of Mine Earthquake Monitoring and Prevention, China University of Mining & Technology, Xuzhou 221116, China
3
Xuzhou Wushuo Information Co., Ltd., Xuzhou 221116, China
4
School of Computer Science & Technology, China University of Mining & Technology, Xuzhou 221116, China
*
Authors to whom correspondence should be addressed.
Minerals 2021, 11(12), 1438; https://doi.org/10.3390/min11121438
Submission received: 8 November 2021 / Revised: 15 December 2021 / Accepted: 16 December 2021 / Published: 19 December 2021

Abstract

:
With the increase in coal mining depth, engineering geological conditions and the stress environment become more complex. Many rock bursts triggered by two combined faults have been observed in China, but the mechanism is not understood clearly. The focus of this research aims at investigating the influence of two combined faults on rock burst mechanisms. The six types of two combined faults were first introduced, and two cases were utilized to show the effects of two combined faults types on coal mining. The mechanical response of the numerical model with or without combined faults was compared, and a conceptual model was set up to explain the rock burst mechanism triggered by two combined faults. The influence of fault throw, dip, fault pillar width, and mining height on rock burst potential was analyzed. The main control factors of rock burst in six models that combined two faults were identified by an orthogonal experiment. Results show that six combinations of two faults can be identified, including stair-stepping fault, imbricate fault, graben fault, horst fault, back thrust fault, and ramp fault. The particular roof structure near the two combined faults mining preventing longwall face lateral abutment pressure from transferring to deep rock mass leads to stress concentration near the fault areas. Otherwise, a special roof structure causing the lower system stiffness of mining gives rise to the easier gathering of elastic energy in the coal pillars, which makes it easier to trigger a rock burst. There is a nonlinear relationship between fault parameters and static or dynamic load for graben faults mining. The longwall face has the highest rock burst risk when the fault throw is between 6 and 8 m, the fault dip is larger than 65°, the mining height is greater than 6 m, and the coal pillar width is less than 50 m. The stair-stepping, imbricate, horst, and ramp fault compared to the other fault types will produce higher dynamic load stress during longwall retreat. Fault pillar width is the most significant factor for different two combined faults, leading to the rise of static load stress and dynamic proneness.

1. Introduction

As a typical dynamic disaster phenomenon, rock bursts have emerged as the most critical issue for underground engineering, which usually leads to casualties, equipment damage, and roadway destruction with the release of lots of elastic energy [1,2,3,4,5,6,7,8,9]. The statistics of a large number of rock bursts show that the area affected by geological fault structure is more prone to dynamic disasters [10,11,12,13,14]. For example, by the late 1970s, there had been at least a thousand fault-induced rock bursts at the Witwatersrand gold mine in South Africa [15]. In Longfeng Coal Mine, faults induced rock burst 36 times (50 times in total), 62% of which occurred when the roadway was close to the faults [8]. 3 November 2011 witnessed that 10 people were killed and 75 were trapped underground in the Qianqiu Coal Mine due to the fault activation [16]. Otherwise, multiple rock bursts triggered by faults were recorded in many other countries, including the USA, Poland, Australia, South Africa, Germany, Russia, and Canada [17,18,19,20,21,22,23].
Considerable efforts have been expended to obtain a comprehensive understanding of rock burst mechanisms triggered by faults or fault slipping. Keneti et al. [11] investigated the influence and mechanism of rock burst corresponding to geological discontinuities (e.g., faults, geological contacts, etc.). The nature of a deformation zone and fault rock related to a rock burst in South Africa was investigated by Stewart, et al. [24]. Uskov et al. [23] investigated the relationship between the tectonic features of faults and rock burst risk in the Norilsk Region, Russia. Li, et al. [25] introduced the concept of the fault coal pillar and suggested that fault-induced rock burst be categorized as fault slip, coal pillar instability, and their coupling type. Cai, et al. [26] verified two types of fault rock burst through similar simulation experiments: mining-induced quasi-static stress-dominated and seismic-based dynamic stress-dominated. Wang, et al. [27] investigated the stress distribution characteristics and instability mechanism of fault surfaces under mining conditions. However, the present works only are focused on the relationship between fault slip and rock burst. There is less attention paid to combined geological structures, such as combined faults. Recently, severe several rocks bursts triggered by two combined faults have been reported in China, such as “2.22” rock burst in the Xinjulong Coal Mine (four people killed), “8.2” rock burst in the Tangshan Coal Mine (seven people killed and five people injured), “6.9” rock burst in the Longjiapu Coal Mine (nine people killed and 12 people injured), “10.20” rock burst in the Longyun Coal Mine (21 people killed and four people injured), and “3.27” rock burst in the Qianqiu Coal Mine (six people killed and 13 people injured) [28]. However, the influence of two combined faults and its structure features on the rock burst mechanism is not understood. A systematical investigation on the rock burst mechanism caused by two combined faults is indispensable to conduct the targeted prediction and prevention.
The physical experiment is a powerful tool for assessing the rock burst potential. However, due to the limitations of model size, complex mechanical properties of coal and rock mass, and excavation sequence of a longwall face, a large-scale parametric investigation cannot be carried out. In comparison, the numerical technique has proved to be a powerful tool for studying the mechanics of underground rock [29]. A large-scale discrete model was established by Wei, et al. [30] to assess the influence of faults on peak abutment stress, the distribution characteristics of strain energy of coal seam, and the energy released by fault slip. An explicit Finite Element Method (FEM) was developed by Manouchehrian and Cai [6] to perform a parametric study with different fault positions, length, and dip relating to the deep tunnel. An FLAC3D model was built up by Jiang et al. [31] to analyze the influence of fault parameters on fault behavior and the possibility of rock burst induced by fault slip. Therefore, it is accepted that numerical simulation has unique advantages in large-scale excavation simulation and rock burst risk assessment [32,33].
In this paper, a numerical investigation is conducted to understand the rock burst mechanism and parameter-influenced mechanical response triggered by two combined faults. The types of two combined faults were first introduced based on the types of two embedded faults and the movement direction of the upper and lower walls of the faults. Two FLAC3D models, i.e., the base model and two combined fault models, were set up to investigate the mechanical response difference between them and the rock burst mechanism triggered by combined faults. Subsequently, a series of parametric tests (including fault throw, dip, mining height, and fault pillar width) were performed to research the influence of different parameters of two combined faults on rock burst potential. In addition, the orthogonal experiment method with four factors and three levels (34) was adopted to identify the main control factors of rock burst triggered by combined faults.

2. Classification of Two Combined Faults under Mining Background

2.1. Classification of Two Combined Faults

A slip fault is a geological phenomenon characterized by a fracture zone existing between two rock blocks, in which two blocks can slip along the fracture [34]. The angle between the fault plane and horizontal surface (known as the dip) and the direction of slip along the fault are used by earth scientists to classify faults. It is widely accepted that faults could be classified as dip-slip and strike-slip types depending on whether they move along the dip direction or move horizontally. Especially for dip-slip faults, two types are commonly adopted: (1) normal fault—The hanging wall of the fault moves downward relative to the footwall of the fault; (2) thrust fault—The hanging wall of the fault moves upward relative to the footwall of the fault.
However, in reality, a fault does not exist independently in nature, which means that a fault combined with other faults or folds could form new composite structures [35,36]. Particularly for two adjacent faults, several two combined faults types were systematically proposed by HU [35]: The scheme of the different two combined faults is depicted in Figure 1, and the detailed identification is shown in Table 1.

2.2. Influence of Two Combined Faults Types on Coal Mining

Over the decades, with the coal resources in China rapidly getting exhausted, the mining depth has increased substantially, e.g., the Huafeng Coal Mine, a typical rock burst coal mine in eastern China, has reached −1300 m below the surface [37]. It was observed that the increasingly complicated engineering geological conditions had posed a great threat to the safety production of Chinese coal mines. Considerable efforts have been expended to acquire a comprehensive understanding of mechanism of the dynamic failure induced by different factors, including roof rupture [38,39,40], fault slip [8,26,41], fold structure [42,43], coal pillar [7,44,45], etc. However, more complex geotechnical conditions are characterized by composite geological structures that have been threatening the safety in underground mining.
Recently, especially for the aforesaid two combined faults, several catastrophic accidents have occurred in China. For example, the “2.22” rock burst disaster (occurred on 22 February 2020) of the Xinjulong Coal Mine in Shandong Province, China, trapped a total of 14 people, four of whom died [46]. Meanwhile, the investigation report [46] indicated that the wedge-shaped graben fault structure (GFS) is one of the direct causes of rock burst.
Figure 2a shows the panel layout and seismic events distribution (captured and calculated by seismic monitoring system) from the time the longwall started to retreat until rock burst occurred on longwall (LW) 2305S in the Xinjulong Coal Mine. To facilitate the understanding of roof strata of LW2305S, the A-A cross-section in Figure 2a is shown in Figure 2b, and the stratigraphic column of LW2305S is shown in Figure 2c. Two faults, i.e., FD6 and FD8 run through LW2304S and extend into LW2305S. As shown in Figure 2b, the two faults form a typical GFS. A tremendous amount of seismic events occurred near the LW2305S between the fault FD6 and FD8. Particularly, there were several high-magnitude events that exceeded 105 J, distributing at the area between two faults. It means that the area between the two faults has a significantly higher dynamic disaster risk than the other areas. On the one hand, the existed GFS will lead to specific roof structures (seen in Figure 2b) and further abnormal roof activity. Moreover, the GFS will reduce mine system stiffness and result in unstable failure of the mine system [6].
Similarly, there are stair-stepping faults detected in Xuzhuang Coal Mine, which is a typical burst-prone coal mine in Jiangsu Province, China [8]. As shown in Figure 3a, two faults, the F113 fault and F115 fault, existed on both sides of the retreat direction of LW7335. Moreover, the dip cross-section of coal mass formed a stair-stepping structure (seen in Figure 3b). Since 2014, the Xuzhuang Coal Mine was equipped with a seismic monitoring system called “SOS” developed by the Central Institute of Mining in Poland to continuously monitor seismicity related to mining activities. Lots of seismic events were detected during the LW7335 retreat between 2018/09/08 and 2019/01/31. As depicted in Figure 3d, many seismic events were distributed near the solid coal area and F115 fault rather than the goaf area. An extraordinary phenomenon also indicated that the special roof with stair-stepping faults could lead to an increasingly dynamic disaster potential.
It follows from the above statement that the two combined faults have been a common geological structure and threaten the safety of coal mines. When there is longwall face mining near the faults, the interaction between faults and coal extraction is a key factor leading to the fault reactivation [26]. Otherwise, in order to minimize the disturbance of operation to the faults, longwall panels are usually separated from the faults by fault pillars. Especially for two combined faults (as shown in Figure 2 and Figure 3), the advance of the longwall face is usually consistent with the strike of the fault to reduce the disturbance to the faults. However, in this case, the fault parameters (i.e., throw, dip, and the width of fault pillars) are usually distinct everywhere due to the complex geological surroundings. Especially rock burst could be triggered by the instability of the fault pillar and disturbance from fault reactivation proposed by Cai et al. [26]. Therefore, the influence of two combined faults and their structure features on the mechanism of rock burst should be further systemically investigated in order to guide the monitoring and prevention of two combined faults-triggered rock bursts.

3. Influence of Two Combined Faults on Rock Bursts

Obviously, understanding the triggering mechanism of two combined faults on rock bursts is an essential task for the improvement of mine safety. In recent years, considerable hypotheses have been expended to acquire a comprehensive understanding of rock burst mechanisms, such as strength [47], rigidity [48,49], burst liability [50,51], stability [52,53], and energy theory [54,55], in which stress, strain, and energy are mainly considered. To further understand the links among the three parameters and explain the rock burst mechanism, a theory that integrated the stress, energy, and displacement (i.e., stiffness) of the roof–coal–floor system was proposed by Zhang, et al. [56]. Therefore, it is critical to systematically research the triggering mechanism of two combined faults on dynamic hazards according to stress, energy, and displacement.
In this section, numerical modeling was performed to investigate the mechanism of rock burst triggered by two combined faults by simplifying to some extent. However, it is not straightforward to recognize stable and unstable rock failure (i.e., rock burst) in numerical modeling. Usually, the judgment of unstable failure in numerical modeling is performed via comparing the model response with a known model response, which is regarded as a stable failure [6]. In longwall excavation modeling, the failure mode in a model without geological structures can be considered stable (base model). On this basis, the stress, displacement, and energy responses of geological models are compared with those of non-geological models to investigate the triggering mechanism of rock burst.

3.1. Model Setup and Material Properties

Taking the graben fault as an example, FLAC3D numerical modeling was performed to investigate the mechanical response caused by two combined faults. In the process, two numerical models (a base model and a graben fault model (GFM)) were established. The only difference between the base model and GFM is the presence of faults, as shown in Figure 4. Figure 4a,c show the numerical model without faults (base model) and with faults (GFM model), respectively. To show the longwall layout more clearly, the grids hyalinized over the longwall face in Figure 4a,c are shown in Figure 4b,d, respectively. The strata were simplified into floor, coal seam, immediate roof, main roof, and overburden according to the classical mining overburden movement theory [57] due to the model not being targeted to a specific engineering context.
Figure 4 displays the numerical models with the size of 500 m × 400 m × 130 m. This paper aims to obtain the mechanical response characteristics of a specific fault structure rather than to do a parametric study on the rock strata and faults. Therefore in this process, the mechanical parameters of strata and faults are obtained according to reliable literature. The origin (i.e., elastic) mechanical parameters of strata are referred from the literature [58], as seen in Table 2. In addition, to obtain the stress distributions in the model, the static displacement boundary is applied around the models. The horizontal displacements around the model and the vertical displacements at the model’s bottom are fixed in the normal direction. At the top of the models, a vertical load (P = γ H = 15.0 MPa, approximately equal to the hydrostatic pressure when the coal seam buried 600 m) was applied to simulate the in situ stress, while stress coefficients along the x- and y-directions were set to 0.8 and 1.5, respectively.
In addition, to simulate the excavation of the longwall face, the coal blocks in the mined area were set to the null model. Moreover, the boundary between the longwall panel and the model boundary should be 50 m to eliminate the boundary effect during the excavation of the longwall face. On this basis, the strain-softening model was adopted to simulate the weakening effect of rock mass after mining. The softening coefficient method suggested in the FLAC3D manual [59] was adopted due to the absence of reliable measurements for the strata after excavation. Moreover, the fault was simulated using the interface elements (seen in Figure 4d). Since the mechanical parameters of the faults cannot really be measured, they were determined by trials in a reasonable range. Thereby, the common values of interfacial cohesion and interfacial friction angle were adopted herein as referred to Karakaş [60]; the normal stiffness (kn) and shear stiffness (ks) suggested in the FLAC3D manual were adopted (seen in Table 2). Among these, kn and ks were primarily estimated to be greater than ten times the equivalent stiffness of the stiffest neighboring zone, which can be expressed as:
k n = k s 10 max [ K + 4 G / 3 Δ Z min ]
where K and G are the bulk and shear moduli, respectively, and Δ Z min is the smallest width of a zone in the normal direction of the interfacing surface.

3.2. Mechanical Response of Two Combined Faults to Mining

To investigate the mechanical response characteristic and triggering mechanism of two combined faults to mining, the difference including mining stress, displacement (i.e., stiffness), and elastic energy between graben fault model (GFM) and the base model was systematically analyzed in this section.
Figure 5 denotes the vertical stress distribution difference between the graben fault model and base model (the cross-section of 10 m ahead of the longwall face is selected since the abutment pressure peak area is about 10 m ahead of the longwall face, which means a higher rock burst potential over there). The vertical stress peak of the base model and GFM is 20.93 MPa and 27.70 MPa, respectively, which also indicates that the vertical stress peak of GFM is approximately 32% (6.77 MPa) higher than that of the base model. Moreover, abnormal stress transitions occur near faults in GFM, as seen in Figure 5b. It can be concluded that the stress transfer on both sides of the coal pillar of the fault is obstructed due to the existence of faults, which results in stress concentration near the fault area. When the fault pillar encounters high stress and a high-stress difference state at the same time, it is easy to produce high-magnitude events, which may lead to fault activation and then trigger rock bursts [26].
Figure 6 shows the vertical displacement distribution difference between GFM and base model (the cross-section of 30 m behind the longwall face is selected, since the first weighting distance of the longwall can vary 20 to 60 m [61], which is accompanied by severe roof collapse and slip of goaf). The vertical displacement peak of the base model and GFM is 212 mm and 302 mm, respectively, which also indicates that the vertical displacement peak of GFM is approximately 42% (90 mm) higher than that of the base model. In addition, the mining-affected span and height of GFM are much larger than that of the base model (compare Figure 6a with Figure 6b). Especially as seen in Figure 6b, the roof displacement field presents a trend of migration from faults to the middle of the longwall face, which means that the existence of faults reduces the stiffness of the roof–coal–floor system.
In numerical modeling, elastic energy density W E can be calculated as follows [58]:
W E = σ 1 2 + σ 2 2 + σ 3 2 2 μ σ 1 σ 2 2 μ σ 2 σ 3 2 μ σ 1 σ 3 2 E
where μ and E are the Poisson ratio and elastic modulus respectively, and σ 1 , σ 2 , and σ 3 are the maximum, intermediate, and minimum principal stresses of the element body, respectively.
Figure 7 depicts the elastic energy distribution difference between the GFM and base model. The elastic energy peak of the base model and GFM is 5 × 105 J and 7.5 × 105 J, respectively, which also indicates that the elastic energy peak of GFM is approximately 42% (90 mm) higher than that of the base model. In addition, not only the energy concentration area of GFM is larger than that of the base model, but also the energy peak area of GFM is closer to the coal mass surface (Compared Figure 7a with Figure 7b). It is not hard to find that rock bursts are easily triggered in elastic energy concentration areas once the roadway and coal mass encounter instability, which may lead to catastrophic casualties and property loss.

3.3. Mechanism of Rock Burst Triggered by Two Combined Faults

It follows from the above argument that the stress, displacement, and energy field of a longwall face mining near two combined faults are obviously different from those without faults. To explain the difference and mechanism of rock burst triggered by two combined faults, a theoretical model was established here.
According to the mechanical behavior of the roof, floor, and coal, the roof and floor can be assumed as linear elastomers, and coal can be assumed as a nonlinear material with strain-softening behavior [26]. When longwall mining is without faults, the stress of the coal can be expressed as σ s , whose stress–strain relationship is shown on the right in Figure 8, and the corresponding roof and floor deformation behavior is shown on the left in Figure 8.
As shown in Figure 8, when the coal mass is under quasi-static loading, the strain is generated after the force is applied, and the strain generated by the cooperative deformation of roof and floor satisfies the relation:
k 1 Δ ε 1 = k 2 Δ ε 2
where k 1 is the loading stiffness of the roof and floor, and k 2 is the unloading stiffness of the coal. The strain of the “roof–coal–floor” system can be expressed as follows:
Δ ε = Δ ε 1 + Δ ε 2 = k 1 + k 2 k 1 Δ ε 2 Δ ε 2 Δ ε = 1 1 + k 2 / k 1
It can be seen that when k 1 + k 2 = 0 , Δ ε 2 / Δ ε . It means that the “roof–coal–floor” structure is in an extremely unstable state, and it is easy to induce rock burst when disturbed by the dynamic load. Furthermore, with the increase in deformation and energy release of the “roof–coal–floor” system, the whole system gradually tends to be stable (corresponding to point S), and the whole process corresponds to the complete process of rock burst. The whole system in this process follows the first law of thermodynamics. In other words, the elastic energy released from the roof and floor (E1) is equal to the energy consumed by coal deformation (E2) plus the elastic energy released outward (E3), from which it can be found that the elastic energy released outward (E3) directly determines the intensity of rock burst. The relationship between them can be expressed as:
E 3 = E 1 E 2 k 1 + k 2 2 k 1 ( σ b + σ a ) Δ ε 2 = k 1 + k 2 2 k 1 k 2 ( σ b 2 σ a 2 )
where σ b and σ a are the stress of coal before and after rock burst (corresponding to points S1 and S in Figure 8, respectively). It can be anticipated that when k 1 + k 2 < 0 , E 3 > 0 , which is consistent with the conclusion of strength and stiffness theory [62]. When the “roof–coal–floor” is disturbed by dynamic load (roof breakage, fault slip, etc.), it can be considered that external energy E4 is applied to it. At this time, the stress state of coal rises from point S1 to point S2, which is closer to the peak strength of coal. More energy is released when instability occurs, and the strength is greater when a rock burst occurs. In addition, when the stress of coal is located at point S 2 , although it does not exceed the peak strength at this time, due to the strain-softening behavior of the coal mass or the dynamic load disturbance, the deformation increases with time and enters the deformation stage after the peak strength (reaching point S1 or S2), which will induce rock burst [63].
It can be seen from Section 3.2 that the stress of coal pillars on both sides of the longwall significantly increases due to the presence of faults, which results in a significant increase in the stress level of the “roof–coal–floor” system (as shown by the green curve in Figure 8). In addition, the existence of combined faults reduces the overall stiffness of “roof–coal–floor” (seen in Section 3.2), which means that the stiffness of roof and floor and the unloading stiffness of coal increase significantly (roof and floor stiffness rise from point k 1 to point k 1 , and coal unloading stiffness rises from point k 2 to point k 2 ). According to Equation (5), when the loading stiffness of the roof and floor and the unloading stiffness of coal increase simultaneously, the elastic energy released during the instability of “roof–coal–floor” will greatly increase (it is fully proved by the energy field distribution in Section 3.2), which could lead to the rock burst and more release of elastic energy.
Therefore, the occurrence of rock burst triggered by two combined faults is the result of the joint action of the “roof–coal–floor” system. The special roof structure near the two combined faults mining prevents longwall face lateral abutment pressure from transferring to deep rock mass, which leads to stress concentration near the fault areas. Otherwise, a special roof structure causing the lower system stiffness of mining gives rise to the easier gathering of elastic energy in coal pillars, which is higher than those without faults and means higher rock burst risk.

4. Effects of Two Combined Faults on Rock Burst Potential

It follows from the above statement that the existence of two combined faults leads to a significant increase in the rock burst potential of the longwall face. Thereby, it is necessary to carry out a systematic investigation, including the influence of different parameters of two different combined faults on rock burst potential.

4.1. Definition of Rock Burst Potential Assessment Indexes

Many case studies have proven that the occurrence mechanism of rock burst should be accounted for the static load, dynamic load, and superposition effect of them [64,65,66,67,68,69]. Moreover, when longwall face mining near the faults, fault pillar instability and fault activation are two main causes of rock burst [25,26]. Therefore, two rock burst potential assessment indexes were proposed in this section considering aforesaid static and dynamic load to investigate the response characteristics of different parameters to mining near two combined faults.
There are three types of principal stress transformation of coal and rock mass under load (as seen in Figure 9): (1) the Mohr circle remains the same size and shifts its position, i.e., | σ 1 σ 3 | remains constant; (2) the Mohr circle changes in size while its position remains constant, i.e., σ 1 and σ 3 get bigger or smaller simultaneously and | σ 1 σ 3 | makes changes; (3) the Mohr circle changes in size and position, i.e., | σ 1 σ 3 | makes changes. In this process, according to the Mohr–Coulomb strength criterion of single discontinuity [70], rock mass failure occurs when the Mohr circle intersects or tangles the strength curve, as shown in Figure 9. It is further deduced that the smaller the distance between the Mohr circle and strength curve is, the worse the stability of the rock mass is and the more rock burst risk.
On this basis, PANG, et al. [71] considering the transformation of three principal stresses proposed the criterion of instability of complete rock mass, which can be expressed as:
k 3 d = | σ 1 σ 3 σ 1 σ 3 | 2 c cot φ + σ 1 + σ 3 2 c cot φ + σ 1 + σ 3
where c and φ represent the cohesion and internal friction angle of coal and rock mass, respectively. It can be inferred that the larger the k 3 d is, the closer to or beyond the strength limit of coal and rock mass, the stronger the disturbance state of surrounding rock, leading to the higher the rock burst risk.
In addition, it can be deduced that in the process of normal stress transformation near the fault, when the principal stress is in the nearly same direction as the fault dip, the fault is easier to be activated (which means that the stiffness of the “roof–coal–floor” system is lower). On this basis, a disturbance instability index (DII) was proposed here as an improvement on k 3 d to assess the rock burst potential caused by static load transformation, which can be expressed as:
D I I = k 3 d f sin
where f sin is the sine value of the direction of maximum principal stress at a certain point in the surrounding rock and the normal vector of the fault plane.
Furthermore, in order to evaluate risk caused by fault activation, maximum fault slip rate (Fv) is proposed to assess the rock burst potential caused by dynamic load (the faster fault slip velocity can lead to greater dynamic load and more rock burst potential [32,72]. Especially in FLAC3D, the Fv of the fault plane can be recorded directly using embedded FISH language.

4.2. Influence of Different Parameters of Two Combined Faults on Rock Burst Potential

When mining longwall faces near the faults, several parameters, including fault throw, fault dip, mining height, and fault pillar width are the essential factors influencing the stability of longwall faces [26]. Thereby, in this section, a series of numerical modeling was performed to systematically research the influence of different parameters of two combined faults on rock burst potential. In this process, the control variable method was adopted. For example, when studying the effect of fault drop changes on the longwall face, five numerical models were set up, in which model parameters and the excavation mode of the longwall face were consistent with the model (as shown in Figure 4c,d), except that the fault throw was successively changed to 4 m, 6 m, 8 m, 10 m, and 12 m.
Figure 10a shows the distribution of DII and Fv with different fall throws during longwall retreat. It can be seen that DII fluctuates as the fault throw increases, and DII reaches the peak when the fault throw is between 6 and 8 m. In addition, Fv firstly increases with the increase in fault throw and then gradually stabilizes in a relatively stable range, which means that the increase in fault throw results in a limited disturbance of dynamic load. Comprehensively considering the response of DII and Fv to the increase in fault throw, the longwall has a relatively high rock burst potential when the fault throw varies from 6 to 8 m.
Figure 10b depicts the distribution of DII and Fv with different fall dips during longwall retreat. It is noted that DII fluctuates as the fault dip increases when fault dip is less than 65°, and DII increases rapidly when fault dip is larger than 65°. Fv reaches the peak when the fault dip is equal to 55°. Moreover, when the fault dip reaches 65°, the rapid rise of DII means that the increase in fault dip will lead to a sharp increase in the stress concentration of the fault pillar.
Figure 10c displays the distribution of DII and Fv with different mining heights during longwall retreat. It is found that both DII and Fv show an obvious upward trend with the increase in mining height. It means that the increase in mining height will lead to the increase in the static and dynamic load of the longwall face, further resulting in higher rock burst risk, which is unfavorable to risk management. In addition, DII and Fv show a more obvious upward trend when the mining height is larger than 6 m. Therefore, the mining height should be reduced to below 6 m as far as possible in order to effectively reduce the rock risk potential.
Figure 10d shows the distribution of DII and Fv with different fault pillar widths during longwall retreat. Both DII and Fv increase first and then decrease with the increase in pillar width, and DII and Fv reach the peak when the pillar width is equal to 30 m and 40 m, respectively. Additionally, both DII and Fv are at extremely high levels when the fault pillar width is between 20 and 50 m. Consequently, the fault pillar width should be more than 50 m to effectively reduce the burst risk.
According to the analysis of the influence of different parameters of two combined faults on rock burst potential, it can be concluded that the change of fault throw, fault dip, mining height, and fault pillar width will make the static and dynamic load change significantly. Furthermore, none of them are linearly related to the level of static or dynamic load. It is observed that the longwall face has the highest rock burst risk when the fault throw is between 6 and 8 m, the fault dip is larger than 65°, the mining height is greater than 6 m, and the coal pillar width is less than 50 m.

4.3. Identification of Main Control Factors of Rock Burst

It has been widely acknowledged that the monitoring, warning, and prevention of rock burst must be targeted [73]. As for mining near the combined faults, the main control factors should be identified for further targeted prevention and control measures. The influenced parameters that triggered rock burst risk such as fault throw, dip, mining height, and fault pillar width have been investigated in Section 4.2. However, mining height can be adjusted flexibly when it leads to a high burst risk. Thereby, in this section, the fault throw, dip, and fault pillar width were tested to identify the main control factors of rock burst in different two combined faults models. Moreover, the DII and IV were used to evaluate the static and dynamic burst risk, respectively.
As for the fault throw, dip, and fault pillar width, if each factor contains four parameters and the experiment is conducted using the control variable method, there will be a total of 64 schemes, which will result in a huge amount of computation and data processing and is unfavorable for data analysis. To solve the problem, the orthogonal experiment method [74] is widely adopted to study this approach of multiple factors and levels (i.e., parameters) problems. Due to three factors being considered here, an orthogonal table with four factors and three levels (34) is adopted, as shown in Table 3. In addition, the difference between the different parameters of each factor should be as large as possible for the more obvious results difference and facilitating the analysis (seen in Table 3). Table 4 shows the orthogonal test schemes of four factors and three levels (34). On this basis, several numerical models (as shown in Figure 11) corresponding to the six aforesaid two combined faults types were set up to identify the main control factors of rock burst, in which the model setups and material properties are the same as those mentioned in Section 3.1.
Figure 12a shows the DII distribution of different orthogonal experiment results. For six different combined faults, the distribution of DII is relatively high and close to each other except for some special cases, which indicates that the excavation of the longwall face in the middle of the faults leads to lateral abutment pressure, and the existence of faults on both sides leads to stress concentration near the faults. Furthermore, the similar DII distribution of different combined faults indicates that the correlation between stress concentration and faults type is weak. Additionally, for a special two combined faults type, the difference of DII among different experimental schemes was also small except for some cases, which means that the longwall face will always be affected by high static stress when mining near the combined faults and targeted prevention and control measures should be implemented.
Figure 12b shows the IV distribution of different orthogonal experiment results. It can be seen that the IV distribution varies among different combined faults, in which all experimental schemes have relatively low IV values in back thrust fault and graben fault models and the rest of the models have relatively high IV values. Furthermore, it is noticed that the maximum IV values all appear in scheme 9 of the other four combined faults models, in which the fault throw and dip are the maximum and the pillar width is the minimum in all experimental schemes (seen in Table 4). The larger the fault throw and dip, the smaller the pillar width when mining near the combined faults, which will lead to higher dynamic risk.
In orthogonal experimental analysis, the sample range (R) is usually used to determine the order of influence factors, which can be expressed as:
R = max { K 1 , K 2 , , K i } min { K 1 , K 2 , , K 3 }
where i represents the number of influencing factors, and K i represents the sum of the corresponding test results when the parameters of any column of orthogonal test table are the same (as shown in Table 4).
Figure 12c denotes the sample range of DII relating to different orthogonal experimental schemes. It can be seen that the fault pillar width is the primary control factor in the stair-stepping fault and back thrust fault model. Fault throw is the main control factor in the graben fault model, and there are few differences between the three factors in the imbricate fault model, horst fault model, and ramp fault model. It can be concluded that the main control factors are variable when mining near the different combined faults. It warns us that the exploration of fault occurrence must be strengthened, especially for the actual mining process, as faults are variable.
Figure 12d presents the sample range of IV corresponding to different orthogonal experimental schemes. It can be seen that fault pillar width is the primary control factor in all fault models. Especially for the stair-stepping fault, imbricate fault, and ramp fault model, the sample range of IV of the fault pillar width is much greater than that of the fault throw and dip, which means inappropriate fault pillars should be avoided primarily to avoid triggering significant activation of the fault.
It follows from the above statement that the presence of two combined faults will result in high static load stress, and specific combined fault types will produce higher dynamic load stress (e.g., stair-stepping, imbricate, horst, and ramp fault model) during longwall retreat. As for different combined faults, the main control factors leading to the rise of static load stress of longwall are different, in which fault pillar width is the most significant factor. In addition, the fault pillar width is the primary control factor in all fault models, leading to dynamic proneness.

5. Conclusions

In underground burst-prone coal mines, the combined faults have become the critical factors that induced rock burst potential. The types of two combined faults were first classified in this work. On this basis, the rock burst mechanism and control factors were systematically investigated. The following conclusions can be drawn:
(1)
According to types of two embedded faults, six types of two combined faults can be summarized, including stair-stepping fault, imbricate fault, graben fault, horst fault, back thrust fault, and ramp fault.
(2)
The particular roof structure near the two combined faults mining preventing longwall face lateral abutment pressure from transferring to deep rock mass leads to stress concentration near the fault areas. Otherwise, a special roof structure causing the lower system stiffness of mining gives rise to the easier gathering of elastic energy in coal pillars, which makes it easier to trigger a rock burst.
(3)
The analysis of the influence of different parameters of graben faults on rock burst potential shows that none of them are linearly related to the level of static or dynamic load. It is noticed that the longwall face has the highest rock burst risk when the fault throw is between 6 and 8 m, the fault dip is larger than 65°, the mining height is greater than 6 m, and the coal pillar width is less than 50 m.
(4)
This approach of combined faults (i.e., stair-stepping, imbricate, horst, and ramp fault model) will produce higher dynamic load stress during longwall retreat. The main control factor for different two combined faults leading to the rise of static load stress and dynamic proneness is different in which fault pillar width is the most significant factor.

Author Contributions

Conceptualization, A.C. and Y.L.; methodology, Y.L.; software, Y.L.; validation, S.J., Q.H. and Y.P.; formal analysis, Y.L. and S.J.; investigation, A.C. and Y.L.; resources, A.C. and Y.L.; data curation, X.Y.; writing—original draft preparation, A.C., Y.L., X.B. and X.Y.; writing—review and editing, A.C. and Y.L.; visualization, A.C. and Y.L.; supervision, Y.L.; project administration, A.C.; funding acquisition, A.C. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was carried out by the funded projects: National Natural Science Foundation of China (51734009), Shandong Provincial Department of Science and Technology (2019SDZY02) and Graduate Research and Practice Innovation Program of Jiangsu Province (KYCX21_2350).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declared that they have no conflict of interest to this work.

References

  1. Akdag, S.; Karakus, M.; Nguyen, G.D.; Taheri, A.; Bruning, T. Evaluation of the propensity of strain burst in brittle granite based on post-peak energy analysis. Undergr. Space 2021, 6, 1–11. [Google Scholar] [CrossRef]
  2. Cao, A.; Dou, L.; Cai, W.; Gong, S.; Liu, S.; Jing, G. Case study of seismic hazard assessment in underground coal mining using passive tomography. Int. J. Rock Mech. Min. Sci. 2015, 78, 1–9. [Google Scholar] [CrossRef]
  3. Faradonbeh, R.S.; Taheri, A.; Karakus, M. The propensity of the over-stressed rock masses to different failure mechanisms based on a hybrid probabilistic approach. Tunn. Undergr. space Technol. 2022, 119, 104214. [Google Scholar] [CrossRef]
  4. Jiang, L.; Wu, Q.; Wu, Q.; Wang, P.; Xue, Y.; Kong, P.; Gong, B. Fracture failure analysis of hard and thick key layer and its dynamic response characteristics. Eng. Fail. Anal. 2019, 98, 118–130. [Google Scholar] [CrossRef]
  5. Khan, M.; He, X.; Farid, A.; Song, D.; Li, Z.; Tian, X.; Ni, M. A novel geophysical method for fractures mapping and risk zones identification in a coalmine, Northeast, China. Energy Rep. 2021, 7, 3785–3804. [Google Scholar] [CrossRef]
  6. Manouchehrian, A.; Cai, M. Numerical modeling of rockburst near fault zones in deep tunnels. Tunn. Undergr. Space Technol. 2018, 80, 164–180. [Google Scholar] [CrossRef]
  7. Vardar, O.; Zhang, C.; Canbulat, I.; Hebblewhite, B. Numerical modelling of strength and energy release characteristics of pillar-scale coal mass. J. Rock Mech. Geotech. Eng. 2019, 11, 935–943. [Google Scholar] [CrossRef]
  8. Wang, C.; Cao, A.; Zhu, G.; Jing, G.; Li, J.; Chen, T. Mechanism of rock burst induced by fault slip in an island coal panel and hazard assessment using seismic tomography: A case study from Xuzhuang colliery, Xuzhou, China. Geosci. J. 2017, 21, 469–481. [Google Scholar] [CrossRef]
  9. Wang, H.; Shi, R.; Lu, C.; Jiang, Y.; Deng, D.; Zhang, D. Investigation of sudden faults instability induced by coal mining. Saf. Sci. 2019, 115, 256–264. [Google Scholar] [CrossRef]
  10. Dou, L.; He, X. Theory and Technology of Rock Burst Prevention; China University of Mining and Technology Press: Xuzhou, China, 2001. [Google Scholar]
  11. Keneti, A.; Sainsbury, B.-A.; Dargaville, R. Consideration of Strain-Bursting Phenomena Associated with Large-Scale Discontinuities: A Numerical Study. Pure Appl. Geophys. 2021, 178, 3581–3600. [Google Scholar] [CrossRef]
  12. Sainoki, A.; Schwartzkopff, A.K.; Jiang, L.; Mitri, H.S. Numerical Modeling of Complex Stress State in a Fault Damage Zone and Its Implication on Near-Fault Seismic Activity. J. Geophys. Res.: Solid Earth 2021, 126, e2021JB021784. [Google Scholar] [CrossRef]
  13. Wang, P.; Jiang, L.; Jiang, J.; Zheng, P.; Li, W. Strata Behaviors and Rock Burst–Inducing Mechanism under the Coupling Effect of a Hard, Thick Stratum and a Normal Fault. Int. J. Geomech. 2018, 18, 04017135. [Google Scholar] [CrossRef]
  14. Wang, P.; Jiang, L.-s.; Zheng, P.-q.; Qin, G.-p.; Zhang, C. Inducing mode analysis of rock burst in fault-affected zone with a hard–thick stratum occurrence. Environ. Earth Sci. 2019, 78, 467. [Google Scholar] [CrossRef]
  15. Pang, Y. Study on Rock Burst Initiation and Failure Propagation; Tsinghua University: Beijing, China, 1999. [Google Scholar]
  16. Feng, J.; Wang, E.; Ding, H.; Huang, Q.; Chen, X. Deterministic seismic hazard assessment of coal fractures in underground coal mine: A case study. Soil Dyn. Earthq. Eng. 2020, 129, 105921. [Google Scholar] [CrossRef]
  17. Keneti, A.; Sainsbury, B.-A. Review of published rockburst events and their contributing factors. Eng. Geol. 2018, 246, 361–373. [Google Scholar] [CrossRef]
  18. Ortlepp, W.D.; Armstrong, R.; Ryder, J.A.; O’Connor, D. Fundamental Study of Micro-Fracturing on the Slip Surface of Mine-Induced Dynamic Brittle Shear Zones. In Proceedings of the 6th International Symposium on Rockburst and Seismicity in Mines, Perth, Australia, 9–11 March 2005. [Google Scholar]
  19. Hofmann, G.F.; Scheepers, L.J. Simulating fault slip areas of mining induced seismic tremors using static boundary element numerical modelling. Min. Technol. 2011, 120, 53–64. [Google Scholar] [CrossRef]
  20. Sainoki, A.; Mitri, H.S. Back analysis of fault-slip in burst prone environment. J. Appl. Geophys. 2016, 134, 159–171. [Google Scholar] [CrossRef]
  21. Alber, M.; Fritschen, R. Rock mechanical analysis of a Ml= 4.0 seismic event induced by mining in the Saar District, Germany. Geophys. J. Int. 2011, 186, 359–372. [Google Scholar] [CrossRef] [Green Version]
  22. Tajdus, A.; Cala, M.; Tajdus, K. Seismicity and Rock Burst Hazard Assessment in Fault Zones: A Case Study. Arch. Min. Sci. 2018, 63, 747–765. [Google Scholar] [CrossRef]
  23. Uskov, V.A.; Eremenko, A.A.; Darbinyan, T.P.; Marysyuk, V.P. Geodynamic Hazard Assessment for Tectonic Structures in Underground Mining of North Ore Bodies in the Oktyabrsky Deposit. J. Min. Sci. 2019, 55, 77–87. [Google Scholar] [CrossRef]
  24. Stewart, R.A.; Reimold, W.U.; Charlesworth, E.G.; Ortlepp, W.D. The nature of a deformation zone and fault rock related to a recent rockburst at Western Deep Levels Gold Mine, Witwatersrand Basin, South Africa. Tectonophysics 2001, 337, 173–190. [Google Scholar] [CrossRef]
  25. Li, Z.; Dou, L.; Cai, W.; Wang, G.; He, J.; Gong, S.; Ding, Y. Investigation and analysis of the rock burst mechanism induced within fault–pillars. Int. J. Rock Mech. Min. Sci. 2014, 70, 192–200. [Google Scholar] [CrossRef]
  26. Cai, W.; Dou, L.; Si, G.; Hu, Y. Fault-Induced Coal Burst Mechanism under Mining-Induced Static and Dynamic Stresses. Engineering 2021, 7, 687–700. [Google Scholar] [CrossRef]
  27. Wang, H.; Shi, R.; Song, J.; Tian, Z.; Deng, D.; Jiang, Y. Mechanical model for the calculation of stress distribution on fault surface during the underground coal seam mining. Int. J. Rock Mech. Min. Sci. 2021, 144, 104765. [Google Scholar] [CrossRef]
  28. Administration, N.M.S. Government Information Publicity. Available online: https://www.chinamine-safety.gov.cn/ (accessed on 15 December 2021).
  29. Starfield, A.M.; Cundall, P.A. Towards a methodology for rock mechanics modelling. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1988, 25, 99–106. [Google Scholar] [CrossRef]
  30. Wei, C.; Zhang, C.; Canbulat, I.; Huang, W. Numerical investigation into impacts of major fault on coal burst in longwall mining—A case study. Int. J. Rock Mech. Min. Sci. 2021, 147, 104907. [Google Scholar] [CrossRef]
  31. Jiang, L.; Kong, P.; Zhang, P.; Shu, J.; Wang, Q.; Chen, L.; Wu, Q. Dynamic Analysis of the Rock Burst Potential of a Longwall Panel Intersecting with a Fault. Rock Mech. Rock Eng. 2019, 53, 1737–1754. [Google Scholar] [CrossRef]
  32. Bizzarri, A. Rupture speed and slip velocity: What can we learn from simulated earthquakes? Earth Planet. Sci. Lett. 2012, 317–318, 196–203. [Google Scholar] [CrossRef]
  33. Sainoki, A.; Mitri, H.S. Dynamic modelling of fault-slip with Barton׳s shear strength model. Int. J. Rock Mech. Min. Sci. 2014, 67, 155–163. [Google Scholar] [CrossRef]
  34. U.S. Geological Survey. What Is a Fault and What Are the Different Types? Available online: https://www.usgs.gov/faqs/what-a-fault-and-what-are-different-types?qt-news_science_products=0#qt-news_science_products (accessed on 10 November 2021).
  35. Hu, M. Structural Geology; Petroleum Industry Press: Beijing, China, 2015. [Google Scholar]
  36. Stawikowski, W. Structural Geology, 2nd ed.; Cambridge University Press: Cambridge, UK, 2017; Volume 23. [Google Scholar]
  37. Chen, S.; Yin, D.; Liu, H.; Chen, B.; Jiang, N. Effects of coal’s initial macro-cracks on rockburst tendency of rock–coal composite samples. R. Soc. Open Sci. 2019, 6, 181795. [Google Scholar] [CrossRef] [Green Version]
  38. Ma, K.; Yuan, F.; Wang, H.; Zhang, Z.; Sun, X.; Peng, Y.; Wang, H. Fracture mechanism of roof key strata in Dongjiahe coal mine using microseismic moment tensor. Geomat. Nat. Hazards Risk 2021, 12, 1467–1487. [Google Scholar] [CrossRef]
  39. Gao, F.; Stead, D.; Kang, H. Simulation of roof shear failure in coal mine roadways using an innovative UDEC Trigon approach. Comput. Geotech. 2014, 61, 33–41. [Google Scholar] [CrossRef]
  40. Mondal, D.; Roy, P.N.S. Fractal and seismic b-value study during dynamic roof displacements (roof fall and surface blasting) for enhancing safety in the longwall coal mines. Eng. Geol. 2019, 253, 184–204. [Google Scholar] [CrossRef]
  41. Wang, H.; Xue, S.; Shi, R.; Jiang, Y.; Gong, W.; Mao, L. Investigation of Fault Displacement Evolution During Extraction in Longwall Panel in an Underground Coal Mine. Rock Mech. Rock Eng. 2020, 53, 1809–1826. [Google Scholar] [CrossRef]
  42. Cao, A.; Jing, G.; Dou, L.; Wu, Y.; Zhang, C. Statistical analysis of distribution patterns of coal seams in fold zones in Northwest China. Int. J. Min. Sci. Technol. 2018, 28, 819–828. [Google Scholar] [CrossRef]
  43. Chen, G.; Wang, Z.; Dou, L. The Study on Prevention of Rockburst when Mining in Folding Structure Areas. In Proceedings of the Asia-Pacific Power and Energy Engineering Conference (APPEEC 2011), Wuhan, China, 25–28 March 2011. [Google Scholar]
  44. Xue, Y.; Cao, Z.; Shen, W. Destabilization and energy characteristics of coal pillar in roadway driving along gob based on rockburst risk assessment. R. Soc. Open Sci. 2019, 6, 190094. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Cao, A.-Y.; Dou, L.-M.; Wang, C.-B.; Yao, X.-X.; Dong, J.-Y.; Gu, Y. Microseismic Precursory Characteristics of Rock Burst Hazard in Mining Areas Near a Large Residual Coal Pillar: A Case Study from Xuzhuang Coal Mine, Xuzhou, China. Rock Mech. Rock Eng. 2016, 49, 4407–4422. [Google Scholar] [CrossRef]
  46. Platform RBC Investigation Report of “2.22” Ground Pressure Impact Accident of Shandong New Julong Energy Co., LTD. Available online: https://mp.weixin.qq.com/s/cZ65mX6ZYuhs5ndyN6kCrA (accessed on 21 April 2020).
  47. Romashov, A.N.; Tsygankov, S.S. Generalized model of rock bursts. J. Min. Sci. 1993, 28, 420–423. [Google Scholar] [CrossRef]
  48. Chen, Z.H.; Tang, C.A.; Huang, R.Q. A double rock sample model for rockbursts. Int. J. Rock Mech. Min. Sci. 1997, 34, 991–1000. [Google Scholar] [CrossRef]
  49. Ma, T.-H.; Tang, C.-A.; Tang, S.-B.; Kuang, L.; Yu, Q.; Kong, D.-Q.; Zhu, X. Rockburst mechanism and prediction based on microseismic monitoring. Int. J. Rock Mech. Min. Sci. 2018, 110, 177–188. [Google Scholar] [CrossRef]
  50. Cai, W.; Dou, L.; Si, G.; Cao, A.; He, J.; Liu, S. A principal component analysis/fuzzy comprehensive evaluation model for coal burst liability assessment. Int. J. Rock Mech. Min. Sci. 2016, 81, 62–69. [Google Scholar] [CrossRef]
  51. Kidybiński, A. Bursting liability indices of coal. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1981, 18, 295–304. [Google Scholar] [CrossRef]
  52. Cook, N.G.W. A note on rockbursts considered as a problem of stability. J. South. Afr. Inst. Min. Metall. 1965, 65, 437–446. [Google Scholar] [CrossRef]
  53. Linkov, A.M. Rockbursts and the instability of rock masses. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1996, 33, 727–732. [Google Scholar] [CrossRef]
  54. Hudson, J.A.; Crouch, S.L.; Fairhurst, C. Soft, stiff and servo-controlled testing machines: A review with reference to rock failure. Eng. Geol. 1972, 6, 155–189. [Google Scholar] [CrossRef]
  55. Xu, S.; Yao, Z.; Zhou, J.; Liu, Y.; Fang, C.; Du, M. Transmittance enhancement of micro-grating structure sapphire with high-refractive index Y2O3 layer. J. Mod. Opt. 2017, 64, 777–780. [Google Scholar] [CrossRef]
  56. Zhang, C.; Jin, G.; Liu, C.; Li, S.; Xue, J.; Cheng, R.; Wang, X.; Zeng, X. Prediction of rockbursts in a typical island working face of a coal mine through microseismic monitoring technology. Tunn. Undergr. Space Technol. 2021, 113, 103972. [Google Scholar] [CrossRef]
  57. Yang, D.; Zhang, Y.; Chen, Z. Analysis on Catastrophe Theory during First Weighting Sliding Instability and Support Crushing of Main Roof with Large Mining Height in Shallow Coal Seam. Appl. Sci. 2020, 10, 5408. [Google Scholar] [CrossRef]
  58. Zhu, G.-a.; Dou, L.-m.; Li, Z.-l.; Cai, W.; Kong, Y.; Li, J. Mining-induced stress changes and rock burst control in a variable-thickness coal seam. Arab. J. Geosci. 2016, 9, 365. [Google Scholar] [CrossRef]
  59. Itasca, I. FLAC3D: Fast Lagrangian Analysis of Continua in 3 Dimension, Version Minneapolis; Itasca Consulting Group: Minneapolis, MN, USA, 2012. [Google Scholar]
  60. Karakaş, A. Practical Rock Engineering. Environ. Eng. Geosci. 2008, 14, 55–57. [Google Scholar] [CrossRef]
  61. Le, T.D.; Bui, X.-N. Effect of Key Parameters on Top Coal First Caving and Roof First Weighting in Longwall Top Coal Caving: A Case Study. Int. J. Geomech. 2020, 20, 04020037. [Google Scholar] [CrossRef]
  62. Cai, W.; Dou, L.; Si, G.; Cao, A.; Gong, S.; Wang, G.; Yuan, S. A new seismic-based strain energy methodology for coal burst forecasting in underground coal mines. Int. J. Rock Mech. Min. Sci. 2019, 123, 104086. [Google Scholar] [CrossRef]
  63. Zubelewicz, A.; Mróz, Z. Numerical simulation of rock burst processes treated as problems of dynamic instability. Rock Mech. Rock Eng. 1983, 16, 253–274. [Google Scholar] [CrossRef]
  64. Cai, W.; Bai, X.; Si, G.; Cao, W.; Gong, S.; Dou, L. A Monitoring Investigation into Rock Burst Mechanism Based on the Coupled Theory of Static and Dynamic Stresses. Rock Mech. Rock Eng. 2020, 53, 5451–5471. [Google Scholar] [CrossRef]
  65. Cao, A.; Jing, G.; Ding, Y.-l.; Liu, S. Mining-induced static and dynamic loading rate effect on rock damage and acoustic emission characteristic under uniaxial compression. Saf. Sci. 2019, 116, 86–96. [Google Scholar] [CrossRef]
  66. Dou, L.-m.; Mu, Z.-l.; Li, Z.-l.; Cao, A.-y.; Gong, S.-y. Research progress of monitoring, forecasting, and prevention of rockburst in underground coal mining in China. Int. J. Coal Sci. Technol. 2014, 1, 278–288. [Google Scholar] [CrossRef] [Green Version]
  67. He, J.; Dou, L.; Gong, S.; Li, J.; Ma, Z. Rock burst assessment and prediction by dynamic and static stress analysis based on micro-seismic monitoring. Int. J. Rock Mech. Min. Sci. 2017, 93, 46–53. [Google Scholar] [CrossRef]
  68. Li, Z.-l.; Dou, L.-m.; Wang, G.-f.; Cai, W.; He, J.; Ding, Y.-l. Risk evaluation of rock burst through theory of static and dynamic stresses superposition. J. Cent. South Univ. 2015, 22, 676–683. [Google Scholar] [CrossRef]
  69. Wang, C.; Cao, A.; Zhang, C.; Canbulat, I. A New Method to Assess Coal Burst Risks Using Dynamic and Static Loading Analysis. Rock Mech. Rock Eng. 2019, 53, 1113–1128. [Google Scholar] [CrossRef]
  70. Jaeger, J.; Cook, N.; Zimmerman, R. Fundamental of Rock Mechanics; Wiley-Blackwell: Hoboken, NJ, USA, 2007. [Google Scholar]
  71. PANG, Y.; WANG, G.; LI, B. Stress path effect and instability process analysis of overlying strata in deep stopes. Chin. J. Rock Mech. Eng. 2020, 39, 682–694. [Google Scholar]
  72. Bhat, H.; Dunham, E. Attenuation of Radiated Ground Motion and Stresses from Three-Dimensional Supershear Ruptures. J. Geophys. Res. (Solid Earth) 2008, 113, B08319. [Google Scholar] [CrossRef] [Green Version]
  73. Dou, L.; He, X. Technique of classification forecasting rock burst in coal mines. J. China Univ. Min. Technol. 2007, 36, 717–722. [Google Scholar]
  74. Box, G.; Hunter, S.; Hunter, W. Statistics for Experimenters: Design, Innovation, and Discovery, 2nd ed.; Wiley-Interscience: New York, NY, USA, 2005. [Google Scholar]
Figure 1. The scheme of the different two combined faults. (a) Plane view represents stair-stepping fault. (b) Plane view represents imbricate fault. (c) Plane view represents graben fault. (d) Plane view represents horst fault. (e) Plane view represents back thrust fault. (f) Plane view represents ramp fault. The red arrows represent the moving direction of the fault block. The black fill represents coal seams and others represent other sedimentary formations.
Figure 1. The scheme of the different two combined faults. (a) Plane view represents stair-stepping fault. (b) Plane view represents imbricate fault. (c) Plane view represents graben fault. (d) Plane view represents horst fault. (e) Plane view represents back thrust fault. (f) Plane view represents ramp fault. The red arrows represent the moving direction of the fault block. The black fill represents coal seams and others represent other sedimentary formations.
Minerals 11 01438 g001
Figure 2. (a) Seismic events distribution (from the time the longwall started to retreat until rock burst occurred) of LW2305S. Purple crosses represent events with energy from 0 to 103 J, blue filled circles represent events with energy from 103 to 104 J, pink unfilled rhombuses represent events with energy from 104 to 105 J, pink filled rhombuses represent events with energy from 105 to 106 J, and red filled circles represent events with energy over 106 J. (b) A-A cross-section of LW2305S in (a). The red arrows represent the moving direction of the fault block. (c) Stratigraphic column of LW2305S. (It is noticed that the stratigraphic column of a longwall face is usually plotted by comprehensively considering the strata distribution around the longwall face rather than for specific sites. The same applies to similar descriptions elsewhere in this work).
Figure 2. (a) Seismic events distribution (from the time the longwall started to retreat until rock burst occurred) of LW2305S. Purple crosses represent events with energy from 0 to 103 J, blue filled circles represent events with energy from 103 to 104 J, pink unfilled rhombuses represent events with energy from 104 to 105 J, pink filled rhombuses represent events with energy from 105 to 106 J, and red filled circles represent events with energy over 106 J. (b) A-A cross-section of LW2305S in (a). The red arrows represent the moving direction of the fault block. (c) Stratigraphic column of LW2305S. (It is noticed that the stratigraphic column of a longwall face is usually plotted by comprehensively considering the strata distribution around the longwall face rather than for specific sites. The same applies to similar descriptions elsewhere in this work).
Minerals 11 01438 g002aMinerals 11 01438 g002b
Figure 3. (a) Seismic events distribution of LW7335 between 2018/09/08 and 2019/01/31. (b) E-E cross-section of the mining coal seam in (a). (c) Stratigraphic column of LW7335. (d) The seismic events frequency distribution on the dip position of the LW7335. Green filled circles represent events with energy from 0 to 103 J, yellow filled circles represent events with energy from 103 to 104 J, and red filled circles represent events with energy over 104 J.
Figure 3. (a) Seismic events distribution of LW7335 between 2018/09/08 and 2019/01/31. (b) E-E cross-section of the mining coal seam in (a). (c) Stratigraphic column of LW7335. (d) The seismic events frequency distribution on the dip position of the LW7335. Green filled circles represent events with energy from 0 to 103 J, yellow filled circles represent events with energy from 103 to 104 J, and red filled circles represent events with energy over 104 J.
Minerals 11 01438 g003aMinerals 11 01438 g003b
Figure 4. Numerical model of base model and graben fault model. (a) Grid model of base model. (b) Schematic of longwall face layout in base model. (c) Grid model of graben fault model. (d) Schematic of longwall face layout in graben fault model. The rock layers above the longwall face in (b,d) were made transparent in order to display the layout of the longwall.
Figure 4. Numerical model of base model and graben fault model. (a) Grid model of base model. (b) Schematic of longwall face layout in base model. (c) Grid model of graben fault model. (d) Schematic of longwall face layout in graben fault model. The rock layers above the longwall face in (b,d) were made transparent in order to display the layout of the longwall.
Minerals 11 01438 g004
Figure 5. Vertical stress distribution curve (cross-section of 10 m in advance of the longwall face). (a) Base model; (b) Graben fault model.
Figure 5. Vertical stress distribution curve (cross-section of 10 m in advance of the longwall face). (a) Base model; (b) Graben fault model.
Minerals 11 01438 g005
Figure 6. Vertical displacement distribution contour (cross-section of 30 m behind the longwall face). (a) Base model; (b) Graben fault model.
Figure 6. Vertical displacement distribution contour (cross-section of 30 m behind the longwall face). (a) Base model; (b) Graben fault model.
Minerals 11 01438 g006
Figure 7. Elastic energy distribution contour (cross-section of 30 m behind the longwall face). (a) Base model; (b) Graben fault model.
Figure 7. Elastic energy distribution contour (cross-section of 30 m behind the longwall face). (a) Base model; (b) Graben fault model.
Minerals 11 01438 g007
Figure 8. Schematic diagram of rock burst mining near the combined faults.
Figure 8. Schematic diagram of rock burst mining near the combined faults.
Minerals 11 01438 g008
Figure 9. Normal stress transformation form based on Coulomb failure criterion where σ 1 and σ 3 are the initial maximum and minimum principal stress of rock mass respectively. σ 1 and σ 3 are respectively the maximum and minimum principal stresses of rock mass after transformation.
Figure 9. Normal stress transformation form based on Coulomb failure criterion where σ 1 and σ 3 are the initial maximum and minimum principal stress of rock mass respectively. σ 1 and σ 3 are respectively the maximum and minimum principal stresses of rock mass after transformation.
Minerals 11 01438 g009
Figure 10. Change of rock burst potential under the influence of different factors based on proposed rock burst potential assessment indexes. (a) Influence of different fault throw on rock burst potential. (b) Influence of different fault dip on rock burst potential. (c) Influence of different mining height on rock burst potential. (d) Influence of different fault pillar width on rock burst potential.
Figure 10. Change of rock burst potential under the influence of different factors based on proposed rock burst potential assessment indexes. (a) Influence of different fault throw on rock burst potential. (b) Influence of different fault dip on rock burst potential. (c) Influence of different mining height on rock burst potential. (d) Influence of different fault pillar width on rock burst potential.
Minerals 11 01438 g010aMinerals 11 01438 g010b
Figure 11. The numerical models of the different two combined faults. (a) Stair-stepping fault model. (b) Imbricate fault model. (c) Graben fault model. (d) Horst fault model. (e) Back thrust fault model. (f) Ramp fault model. The Minerals 11 01438 i001 blocks represent overburdens; the Minerals 11 01438 i002 blocks represent main roofs; the Minerals 11 01438 i003 blocks represent immediate roofs; the Minerals 11 01438 i004 blocks represent coal seams; and the Minerals 11 01438 i005 blocks represent floors.
Figure 11. The numerical models of the different two combined faults. (a) Stair-stepping fault model. (b) Imbricate fault model. (c) Graben fault model. (d) Horst fault model. (e) Back thrust fault model. (f) Ramp fault model. The Minerals 11 01438 i001 blocks represent overburdens; the Minerals 11 01438 i002 blocks represent main roofs; the Minerals 11 01438 i003 blocks represent immediate roofs; the Minerals 11 01438 i004 blocks represent coal seams; and the Minerals 11 01438 i005 blocks represent floors.
Minerals 11 01438 g011aMinerals 11 01438 g011b
Figure 12. Orthogonal experimental results. (a,b) represent DII and IV distribution of different orthogonal experimental schemes, respectively. (c,d) represent the sample range of DII and IV correspond to different orthogonal experimental schemes respectively.
Figure 12. Orthogonal experimental results. (a,b) represent DII and IV distribution of different orthogonal experimental schemes, respectively. (c,d) represent the sample range of DII and IV correspond to different orthogonal experimental schemes respectively.
Minerals 11 01438 g012aMinerals 11 01438 g012b
Table 1. The detailed identification of two combined faults (as referred to HU [35]).
Table 1. The detailed identification of two combined faults (as referred to HU [35]).
Types of Two Combined FaultsTypes of Two Embedded FaultsIdentification
Stair-stepping faultTwo normal faultsThe stair-stepping fault is formed by two parallel normal faults, the hanging walls of which descend in the same direction, resembling a staircase.
Imbricate faultTwo thrust faultsThe imbricate fault is formed by two parallel thrust faults, the hanging walls of which descend in the same direction, resembling a roof tile on the house.
Graben faultTwo normal faultsThe graben fault is composed of two symmetrical and inclined normal faults with fault walls rising on both sides and a common descending wall in the middle.
Horst faultTwo normal faultsThe horst fault is composed of two symmetrically inclined normal faults with fault walls descending on both sides and common uplifting walls in the middle.
Back thrust faultTwo thrust faultsIt is similar to the graben fault, but the back thrust fault consists of two thrust faults.
Ramp faultTwo thrust faultsIt is similar to the horst fault, but the ramp fault consists of two thrust faults.
Table 2. Model mechanical parameters and thickness of the fault and strata. /: The interface elements do not have thickness assignments.
Table 2. Model mechanical parameters and thickness of the fault and strata. /: The interface elements do not have thickness assignments.
ItemsThickness (m)Bulk Module
(GPa)
Shear Module
(GPa)
Cohesion
(MPa)
Internal Friction Angle
(°)
Density
(kg·m−3)
Overburden371.310.420.30332500
Main roof202.301.521.7382700
Immediate roof51.310.420.30332500
Coal seam30.150.080.11251400
Floor651.971.140.40352600
Interface elements/interfacial cohesion = 2.0 MPa, interfacial friction angle = 30°, kn = ks ≈ 2.16 × 1010 Pa/m.
Please note that only the normal stress and tangential stress of the fault are considered here without the effect of pore water pressure, which means that the mechanical parameters here do not consider the effect of pore water.
Table 3. Orthogonal test factor levels and model parameters.
Table 3. Orthogonal test factor levels and model parameters.
LevelABC
Throw (m)Dip (°)Pillar Width (m)
144510
286040
3127570
Table 4. Orthogonal test schemes of four factors and three levels (34).
Table 4. Orthogonal test schemes of four factors and three levels (34).
Level (Scheme)ABCD
Throw (m)Dip (°)Pillar Width (m)/
144510/
246040/
347570/
484570/
586010/
687540/
7124540/
8126070/
9127510/
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cao, A.; Liu, Y.; Jiang, S.; Hao, Q.; Peng, Y.; Bai, X.; Yang, X. Numerical Investigation on Influence of Two Combined Faults and Its Structure Features on Rock Burst Mechanism. Minerals 2021, 11, 1438. https://0-doi-org.brum.beds.ac.uk/10.3390/min11121438

AMA Style

Cao A, Liu Y, Jiang S, Hao Q, Peng Y, Bai X, Yang X. Numerical Investigation on Influence of Two Combined Faults and Its Structure Features on Rock Burst Mechanism. Minerals. 2021; 11(12):1438. https://0-doi-org.brum.beds.ac.uk/10.3390/min11121438

Chicago/Turabian Style

Cao, Anye, Yaoqi Liu, Siqi Jiang, Qi Hao, Yujie Peng, Xianxi Bai, and Xu Yang. 2021. "Numerical Investigation on Influence of Two Combined Faults and Its Structure Features on Rock Burst Mechanism" Minerals 11, no. 12: 1438. https://0-doi-org.brum.beds.ac.uk/10.3390/min11121438

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