Control of the Chamber in Extra-Thick Coal Seam Based on Fracture Characteristics of Coal Rock and Mass

As the development characteristics of cracks of the chamber are closely related to the bolt and cable anchoring parameters, the original rock stress of the surrounding rock, and fracture degree, this paper proposes the control method for chamber in extrathick coal seam by the crack development characteristics. Based on the numerical model analysis, the chamber generates a large amount of shear cracks.-en, these cracks further extend and expand.When the cracks penetrate each other, the chamber bearing capacity decreases sharply, and the corresponding displacement also increases significantly. Based on the above rules, a chamber control method by high-strength anchor bolts and grouting reinforcement is proposed to improve the rock shear strength. By optimizing the anchor bolt parameters and grouting parameters, such as spacing, length, and the grouting pressure, the cracks are significantly reduced and the cracks are prevented from penetrating. Finally, an industrial test was conducted in the chamber at Datong Coal Group, and the chamber deformation was well controlled.


Introduction
With the development of the coal industry, several largescale coal production bases have successively appeared in China, forming many modern high-yield and high-efficiency mines with a capacity of 10 million tons [1]. For example, the Datong Coal Group can produce more than 0.13 billion tons coal in one year, while the output of the mines with a capacity of 10 million tons accounts for about 70% [2]. e main coal seam of those modernized mines exploits generally thick and extra-thick coal seams [3,4]. Due to the large coal output of these mines, the sections of underground electromechanical chambers, gangue washing chambers, and assembly chambers are usually above 70 m 2 [5]. When the coal seam is thick and the strength is high, the chamber is preferentially put in the coal seam. e support theory for chambers in such extra-thick coal seams mainly includes the loosened ring theory, reinforcement arch theory, and prestressed support theory. e loosened ring theory generally emphasizes that the anchor ends of the anchor bolt are arranged outside the loosened ring, which provides a force point for the pretightening force application [6]. is theory provides theoretical guidance for the design of lengths of anchor bolt but does not consider some support parameters such as pretightening force (PTF) and anchor bolt spacing [7]. In addition, the loosened ring of the chamber is usually replaced by a plastic zone. As shown in Figure 1(a), the chamber in coal seams has a large span and height, and the strength of the coal seam is lower than that of the rock. erefore, the plastic zone is generally large. Although the anchoring end of the anchor cable can generally be placed in an intact zone outside the loosen ring, the loosen ring thickness is often greater than the bolt length [8]. erefore, according to the loosen ring theory, it is impossible to explain the support effect for the bolt in the fracture zone of a chamber in the extra-thick coal seam. At the same time, it cannot provide a theoretical basis for the anchor bolt spacing.
Since the tensile strength of coal is lower than its shear strength, tensile damage is formed on the surface of the chamber, as shown in Figure 1(b). Based on this feature, the reinforcement arch theory and prestressed support theory are proposed. e reinforcement arch theory and prestressed support theory optimize the anchor bolt support parameters to make the formation of compressive stress layer as much as possible on the surface of roadways or the chambers [9]. Both of them can provide basis for the support parameters, such as anchor bolt spacing, length, and pretightening force of the chamber, but they ignore the original rock stress and fracture degree of the chamber.
After the formation of chambers, its deformation and failure are mainly the process of initiation, expansion, and penetration of cracks [10][11][12]. e chamber will inevitably produce many cracks, as shown in Figure 2. e crack development is directly connected with the anchor bolt support parameters, and it is also closely related to the original rock stress and the crack development degree of the chamber. erefore, this paper proposes the control method of the chamber of extra-thick coal seam based on the fracture development characteristics. e research of this subject can provide the basis for the selection of support parameters for anchor bolt and cable. At the same time, it can consider the original rock stress and fracture degree of the surrounding rock of the chamber.
After introducing the chamber at Datong Coal Group, this paper uses numerical simulation to establish a numerical model for the development of cracks in the chamber, which is used to reveal the relationship between the fracture development characteristics and the chamber deformation. On this basis, the corresponding control technology of the chamber is proposed, and finally the industrial test is conducted in the chamber at Datong Coal Group.

Background
Tashan Coal Mine is the biggest mine in the Datong Coal Group. It is currently mining carboniferous coal seams, and its annual production of coal per working face is 15 million tons. Due to the large coal output, the mine has high requirements for the section of the chamber. According to the requirements of the electromechanical equipment for the installation space, as shown in Figure 3(a), the chamber size is 9.5 m × 11.85 m. Due to the high strength of the #3-5 coal seam, the chamber is put in the #3-5 coal seam. From the drill core near the chamber, the buried depth of the coal seam is 528.5 m and its thickness is 18.1 m. e floor and roof are all mudstone, and the rock stratum is shown in Figure 3(b).

Constitutive Model of Surrounding Rock of
the Chamber e choice of the constitutive model for rock in the numerical model is crucial. is chapter first analyzes three typical constitutive models of the rock. en based on the brittleness characteristic of coal rock mass, the appropriate constitutive model is chosen.

Constitutive Model of the Rock.
In order to get the constitutive model for the rock, domestic and foreign scholars have given three typical constitutive models of rock based on a large number of compaction tests, as shown in Figure 4 [13]. When the rock reaches the peak stress, three different paths appear: Path ①: the ideal elastoplastic model, that is, the stress remains unchanged after the coal rock mass reaches the peak stress. Path ②: when the rock reaches the peak stress, the stress falls along a certain curve and decreases to a certain stress level. Path ③: after the rock reaches the peak stress, the stress suddenly drops to a certain level. e stress-strain curve of the rock in the traditional numerical model often adopts Path ①, which considers the rock as an ideal elastoplastic model. However, a large number of experiments and engineering practices have shown that the rock shows a strain-softening phenomenon after reaching the peak, namely Path ② and Path ③. When the brittleness of the rock is strong, the stress-strain relationship after peak is close to Path ③. e coal from Tashan Coal Mine has typical brittle characteristics, so this paper adopts the constitutive model of Path ③ [14].  2 Advances in Civil Engineering

Strain-Softening Model Based on the Hoek-Brown
Criterion. In order to consider native joints or cracks in the chamber, the Hoek-Brown strength criterion is used in this paper. Its specific expression is as follows [15]: where σ 3 and σ 1 are the third and first principal stresses at failure, respectively; m b , s, and a are constants; σ c is the uniaxial compressive strength. e expressions of m b , s, and a are as follows [16]: where GSI is the geological strength index. e parameter D is the "disturbance factor." m i is the material constant for the intact rock. e GSI, m i , and D values are the key parameters for calculating m b , s, and a; the RocData software is used to find the m i values for different lithologies. e disturbance factor D can be calculated from the following:

Advances in Civil Engineering
where λ is the ratio of elastic wave velocity before and after excavation. e GSI value is evaluated based on the rock structure and the weathering degree on the surface of rock. e accuracy of GSI evaluation depends mainly on the experience of geologists. In order to facilitate the use of the GSI system and improve its accuracy, Cai et al. [14] gave a method for quantifying the GSI value based on the rock mass V b and the joint surface coefficient J c . e specific equation is given as follows: V b and J c in the above equation can be expressed as follows: where s 1 , s 2 , and s 3 are the joint spacing of the rock in three directions. c 1 , c 2 , and c 3 are the joint angles of the joint faces in three directions of the rock block. J W , J S , and J A are macroscopic fluctuation coefficients, joint microsmooth coefficients, and joint alteration coefficients, respectively. As shown in Figure 5, the joint spacing of the rock in three directions of the rock is obtained by field statistics.
For J W , J S , and J A , the specific values are as follows: e J W values are presented in Table 1. According to the joint macroscopic fluctuation of the rock in Tashan Coal Mine, the value here is 1.00.
e J s values are presented in Table 2. According to the joint smoothness of the rock in Tashan Coal Mine, the value here is 0.75. e J A values are presented in Table 3. According to the alternation of the rock in Tashan Coal Mine, the value here is 0.75.
According to equations (4)∼(6), the GSI value is given. e values of m b , s, and a can be calculated using the GSI, m i , and D values, as shown in Table 4.
To simulate the strain-softening behavior of coal rock after peak stress, Cai et al. proposed a strain-softening model based on the Hoek-Brown strength criterion, which is embodied by the residual GSI r . After the peak stress, the strain-softening process of rock mass is simulated by reducing the initial GSI. e residual GSI r value can be calculated from initial GSI value as follows [14]:

Slip and Open Model for Crack Development.
Many cracks will occur during the destruction process of the chamber. In order to simulate the crack development, the chamber is divided into the Voronoi elements. e slip and open model between the Voronoi elements are shown in Figure 6 [17]. e equations for K n and K s in the model are as follows [18]: where ΔZ min is the smallest width of the zone. G is the shear modulus and K is the bulk modulus of the rock. ey can be expressed as follows [19]: where E is the elastic modulus and μ is Poisson's ratio, respectively. To further calibrate K n and K s , the uniaxial compression is simulated to obtain the relationship of stress and strain, which is then adjusted with the test. After multiple adjustments, more accurate K n and K s are obtained.

Rock Strength Reduction.
e chamber contains native cracks and joints. erefore, the experimentally determined rock mass mechanical parameters are usually higher than those of the rock from the stope. erefore, the strength of the rock in the laboratory is reduced by the RQD (rock quality designation) value in engineering, which can then be used as the mechanical parameter of the rock from the stope. e elastic modulus E is from the laboratory rock, while the elastic modulus E m is from the stope rock. ey have the following relationship [20]: e σ cm is the uniaxial compressive strength from laboratory rock, while σ cm is the uniaxial compressive strength from the stope. Similarly, σ cm /σ c has the following relationship with E m /E [21]:

Advances in Civil Engineering
where q is a key parameter. According to reference [22], it is taken as 0.7. e uniaxial tensile strength σ tm is about 0.1 times the uniaxial compressive strength σ cm [23]. From the experimentally obtained rock parameters, the mechanical parameters in the stope are calculated by equations (10) and (11); see Table 5 for details.

Model Establishment.
e specific modeling process of the model is as follows: according to the geological conditions and size of the chamber, as shown in Figure 7, a 2D numerical model with a size of 150 m long and 54.43 m width is established. e mesh size of the model is 0.3∼1.5 m. e strainsoftening model based on the Hoek-Brown criterion is employed for coal rock mass. e horizontal displacement on around model is constrained, and the model bottom is constrained on the vertical displacement [24]. e model top is applied with a 12.5 MPa vertical load [25].

Verification of the Numerical Simulation.
To verify the rationality and accuracy of the constitutive model and mechanical parameters of the numerical model, a      Advances in Civil Engineering verification study is performed by comparing the crack development characteristics and the surrounding rock deformation characteristics. As shown in Figure 7, the integrity of the rib at 2, 5, and 10 m is obtained using the borehole television. From Figure 8, a large amount of damage is generated in the coal body at 2 m, forming a fracture zone. At 5 m, the coal body has many cracks, while at 10 m, there are only a few cracks in the coal body. In this position, the numerical simulation of the fracture development is in good consistent with the observation results of the borehole television.
erefore, crack development characteristics verify the reliability of the simulation model.
As shown in Figure 9, a multipoint displacement meter is placed in the dome of the chamber to monitor the amount of roof separation and compare it with numerical simulation. From Figure 8, the numerical simulation results of the chamber roof separation are in good consistent with the field measurement results in terms of both value and trend. e chamber deformation verifies the reliability of the numerical simulation.

Crack Development Characteristics of the Chamber and Its Relationship with Deformation
After the excavation of the chamber, the chamber changes from the original triaxial stress state to the biaxial stress state, and even to a uniaxial stress state. Many cracks are produced in the chamber. e fractures initiate, expand, and eventually penetrate, which will inevitably destroy the integrity of the surrounding rock and reduce its resistance to deformation. is section studies the relationship between fracture development characteristics and deformation and reveals the mechanical mechanism of surrounding rock fracture development.

Relationship between Fracture Development
Characteristics and Deformation. As shown in Figure 10, in order to reveal the relationship between the fracture development characteristics and the surrounding rock deformation, cracks in the both ribs and the lateral displacement of the ribs are separately extracted. As shown in Figure 10(a), many cracks are produced in the chamber. In the range of 0∼7.4 m, the cracks penetrate each other. In the range of 7.4∼12.7 m, the cracks initiate and develop, while in the range of 12.7∼20 m, there are basically no cracks in the both ribs.
In this paper, the surrounding rock is divided into three areas, which are intact area, crack initiation and development area, and crack penetration area. e above three areas are mainly divided according to the development of cracks in the numerical simulation. e cracks penetrate each other in the fracture penetration area, while in the crack initiation and development area, although there are cracks, the cracks do not penetrate each other. And, there are no cracks in the intact area, as shown in Figure 11. e following can be seen from Figure 10 (1) At the range of 0∼7.4 m (I area), the lateral movement of the chamber coal rib is obvious. is zone is the crack penetration area of the ribs. e lateral  displacement is largest at the surface of the chamber ribs, the lateral displacement of which is 0.43 m. In Figure 10(a), the closer to the surface of the chamber ribs, the larger the lateral displacement variation, which indicates the more serious damage and weaker bearing capacity. Under the effect of mine pressure, serious lateral displacement occurs in this penetration area. (2) At the range of 7.4∼12.7 m (II area), the variation of the lateral displacement in chamber ribs is significantly reduced. In Figure 10(a), this area is the initiation and development area of the rib cracks. No crack penetration occurs in this area, and therefore, the surrounding rock has certain bearing capability and the lateral displacement is significantly reduced. (3) At the range of 12.7∼20 m (III area), the lateral displacement of the chamber ribs is basically 0.05 m.
In Figure 10(a), this area is the intact area of the ribs. e lateral displacement of the surrounding rock is caused by the surrounding rock movement.
According to the above rules, many cracks are produced in the chamber due to the effect of mining pressure. After the cracks penetrate each other, the chamber bearing capacity decreases sharply, which is manifested by a sharp increase in the displacement of the chamber in the crack penetration area. To reduce the occurrence of cracks and inhibit the penetration of cracks, it is necessary to further reveal the reason of crack development in the chamber.

Mechanism of Crack Development.
According to the above laws, the deformation of surrounding rock is due to the crack development, and the deformation of chamber is  Advances in Civil Engineering large in the crack penetration zone. To reveal the reason of the crack development in the chamber, as shown in Figure 12, the mechanism of crack generation is given.
As shown in Figure 12, due to mining pressure and mining disturbance, many cracks are generated in the chamber. e main mechanism for crack generation is due to the fact that shear stress in the chamber is over strength limit. As a result, shear failure occurs. Compared with shear failure, tensile failure is relatively less. e percentage of cracks from shear failure and tensile failure of the chamber are 98.2% and 19.3%, respectively. erefore, the shear failure zone of the chamber occurs substantially simultaneously with the tensile failure. e shear failure causes the friction angle and cohesive force of the chamber to decrease sharply and the tensile strength also drops. erefore, the shear failure zone also undergoes tensile failure at the same time.

Control of the Chamber in Extra-Thick Coal Seam
From the relationship between crack development characteristics and deformation of the chamber and the mechanism of the crack development, it is gotten that increasing    Advances in Civil Engineering the shear strength limit of the chamber and reducing the crack penetration is an effective method to control the chamber deformation.

Grouting Reinforcement.
After the excavation of the chamber, cracks are initiated in the chamber. Due to the mine pressure, the cracks continue to expand and eventually penetrate. e grouting technology injects the slurry into the cracks of the coal rock mass, making it a new structure with high strength and high stability, thus achieving a good support effect. At present, the grouting method commonly used in coal mines is the pressure injection method. In this method, the chemical slurry is injected into the coal rock mass by a pressure injection pump and the loose coal rock mass is connected into an organic whole. From the analysis of coal rock body failure mechanism and mechanical principle, the main effects of slurry on coal rock mass can be summarized as network skeleton effect, bond reinforcement effect, filling compaction effect, and failure transformation mechanism [26].
(1) Network skeleton effect. During the reinforcement construction process, the slurry penetrates into cracks of the coal rock mass. After reacted and consolidated, a new network skeleton structure is formed in the coal rock body. When the external load exceeds the rock strength limit, it will be destroyed and deformed greatly. en, the reinforcement material network will play a skeleton role with its good toughness and bond strength performance, thereby improving the shear strength limit of the rock [27]. (2) Bond reinforcement effect. e cohesion and tensile strength on the crack of coal rock mass are almost zero. After reinforcement, the reinforcement material can bond the fractured coal rock mass, which is already in a plastic state, maintaining its cohesive force and friction angle [28]. e shear strength limit of rock can be improved, according to the Mohr-Coulomb criterion.
(3) Filling compaction effect. Under the influence of the pumping pressure and the expansion pressure formed by the self-reaction, the slurry can fill the gaps of the interconnected coal rock mass and can also fill and compact the coal rock mass as a whole. is effect can significantly increase the shear strength of coal rock mass. (4) Failure transformation mechanism. When there is a crack in a continuous medium, a strong stress concentration is formed during the load-bearing process with the maximum stress being concentrated at the tip of the crack. After reinforcement, the crack is filled with reinforcing material, and the bonding effect of the material on the crack will greatly weaken or eliminate the stress concentration at the crack tip, thereby avoiding rock shear failure [29].
e above-mentioned effects of the slurry increase the shear strength of the rock and reduce the occurrence of shear failure.

High Prestressed Anchor.
Long anchor cable or lengthened anchor bolt can provide axial and shear forces simultaneously. First, the force of the surrounding rock is changed by the axial force, and the anchor is transformed from the biaxial stress state to the triaxial stress state. Second, through the transverse connection with the anchor, the shear stress and bending distance are born [30]. e anchor bolt can improve the rock cohesive force. e specific formula is as follows: where c and φ are the cohesive force and friction angle before the anchor failure, respectively. c 0 (MPa) is the rock cohesion without the anchor bolt. n is the number of anchor arrangements in the anchor body. c m (MPa) is the additional cohesion provided by a single anchor and is expressed as In the above equation, F bmax is the maximum shear force that the bolt can withstand under pure shear conditions. Its expression is as follows: where σ s is the bolt yield strength, MPa; d is the bolt diameter, m; and s is the free area of the test body, m 2 .
In addition, the anchor can also improve the friction angle of the surrounding rock φ [31]. According to the Mohr-Coulomb strength criterion, the shear strength limit of the rock can be increased to some extent. e anchor cable can provide a high pretightening force (PTF) to reinforce the chamber, which effectively connects Advances in Civil Engineering the broken surrounding rock together to form a unified composite and thus significantly improves the shear strength of the crack [32,33].

Control Technology of Surrounding Rock.
Based on theoretical analysis, after the anchor support and grouting reinforcement, the shear strength limit of the chamber can be improved. However, the method for determining the anchor parameters and grouting parameters is not yet given. e anchor parameters and grouting parameters will be optimized in this section based on the crack development characteristics by means of numerical simulation.

Anchor Support.
e support parameters of the anchor include spacing, length, PTF, etc. According to the theoretical analysis, the anchor can improve the shear strength limit of the rock, inhibit and reduce the surrounding rock shear failure, and decrease the possibility of the cracks interpenetrating. Based on this theory and taking the PTF of the bolt as an example, the method of determining the pretightening force is given. As shown in Figure 13, the crack expansion characteristics of the chamber ribs are given when the PTF of the bolt is 25, 50, 100, and 125 kN. e following can be seen from Figure 13: (1) When the PTF of the anchor bolt is 25 kN, the surrounding rock of the chamber ribs forms a crack penetration area, a crack initiation and development area, and an intact area. e range of the crack penetration area is 0∼6.7 m. In this area, the surrounding rock is seriously damaged, resulting in a drastically decreased bearing capacity and serious deformation.
(2) When the PTF of the anchor bolt is 50 kN, the crack expansion of the chamber ribs is similar to that at 25 kN, but the crack penetration area is significantly reduced. (3) When the PTF of the anchor bolt is 100 kN, the surrounding rock of the chamber ribs forms a crack initiation and development area and an intact area, but there is no crack penetration area. e chamber ribs have a certain load-bearing capacity, and thus the deformation and damage are well controlled. (4) When the PTF of the anchor bolt is 125 kN, the chamber ribs also only have crack initiation and development area and intact area. Moreover, the crack initiation and development area are reduced, while the intact area is increased, indicating that the deformation and damage of the chamber are well controlled.
From the above result, when the PTF is 100 kN, the surrounding rock of the chamber ribs only has a crack initiation and development area and an intact area, and the surrounding rock deformation is well controlled. erefore, the PTF of the anchor bolt is determined to be 100 kN. Other parameters of the anchor are determined by optimization methods similar to that of the anchor bolt PTF. e final support parameters are determined: the spacing of anchor bolts is 900 mm × 900 mm, the PTF is 100 kN, the length of the dome anchor bolt is 2500 mm, and the anchor bolt length at two ribs is 1800 mm. e spacing of anchor cables is 2700 × 1800 mm, the PTF is 150 kN, and the length is 8300 mm.

Grouting Reinforcement.
In the grouting reinforcement, the spacing of the grouting pipe and the grouting pressure are the core parameters. e grouting spacing is determined by the diffusion radius of the slurry, which mainly depends on the grouting pressure, the degree of crack development, the rock properties, etc., and the flow dynamics parameters of the slurry. e diffusion radius is closely related to the grouting pressure under certain surrounding rock conditions and grouting materials. A numerical simulation is built, and the relationship of the diffusion range and grouting pressures is obtained. From the simulation, while the grouting pressure reaches 5.0 MPa, the slurry diffusion radius is basically stable at 1800 mm. erefore, the anchor bolt spacing is 900 mm, and the grouting pressure is 5 MPa.

Fracture Control Effect.
From the above analysis, it is found that reasonable grouting and anchor support parameters can inhibit the penetration of cracks and control the deformation of the chamber. As shown in Figure 14, the crack development characteristics in the chamber under the combined effect of these two supporting methods are further analyzed. en, the surrounding rock control effect is evaluated. e following can be drawn from Figure 14: (1) By combined support, there are only a few cracks in the chamber dome, while cracks are relatively more on the dome sides. However, in general, the fractures significantly decrease and there is no crack penetration zone. erefore, the cracks in the chamber dome are well controlled. (2) By combined support, there are a few cracks in the chamber ribs, but the cracks do not penetrate each other.
Based on the above analysis, after the combined support of grouting and high-strength anchor, the cracks in the chamber only have crack initiation and development area and intact area. According to the relationship between the crack development characteristics and deformation of the chamber, it can be drawn that the chamber deformation can be well controlled.

Engineering Application
7.1. Support Scheme. From the above laws, grouting equipment with a grouting pressure of 5.0 MPa is placed every two anchor bolts. e anchor support is applied after grouting, as shown in Figure 15. Specific supporting parameters are given in Section 6. Figure 16, deformation at different chamber depths are monitored by the multipoint displacement meters. ey are put at depths of 2, 4, 6, 8, and 10 m.  I  II  III   I  II  III   II  III   II  III   I-crack penetration area  II-crack initiation and development area  III-intact area   0  4  8  12 16 20 Figure 13: Optimization for the PTF for anchor bolt.  It can be seen from Figure 16(a) that the maximum value of lateral displacement is 18 mm and the area is the crack initiation and development area. e lateral displacement of the surrounding rock beyond 6 m is basically kept at 6 mm, and the area is the intact area.

Chamber Deformation Measurement. As shown in
It can be seen from Figure 16(b) that the maximum value of roof separation is 12 mm and the area is the crack initiation and development area. e roof separation beyond 4 m is basically kept at 9 mm, and the area is the intact area.
From the above laws, the chamber deformation is well controlled.

Conclusion
After the introduction of geological conditions of the chamber, this paper reveals the relationship between the crack development characteristics and deformation of surrounding rock, as well as the mechanical mechanism of crack generation based on numerical simulation analysis. e following results can be drawn: (1) e drawbacks of the loosened ring theory, reinforcement arch theory, and prestressed support theory are analyzed. e control theory of surrounding rock in a large-section chamber of extra-thick coal seam based on chamber crack development characteristics is proposed. (2) After driving the chamber, much shear damage is generated in the chamber ribs due to the effect of the mine pressure and excavation disturbance, resulting in the formation of many cracks. ese cracks further extend and expand. When the crack penetrates, the bearing capacity of the chamber ribs decreases sharply and the deformation of the surrounding rock increases significantly. (3) From crack development characteristics and the surrounding rock deformation, as well as the mechanism of crack generation, the support technology combining high-strength anchor with the grouting reinforcement is proposed. By optimizing the above supporting and grouting parameters, the generation of cracks in the chamber is controlled, and crack penetration is suppressed. Finally, an industrial test was conducted in the chamber at Datong Coal Group, which showed that chamber deformation and failure were well controlled.

Data Availability
Part of the data used to support the findings of this study are included within the article, and the rest of the data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.