Abstract

Based on the thermodynamic and elastodynamic theories, the controlling equation of temperature-stress coupling action on rocks containing random damage units is established by combining the Mogi-Coulomb damage criterion. And the numerical calculation model of combustion cavity expansion under temperature-stress coupling condition is established by using ABAQUS secondary development program. The fracture field evolution law during the expansion of the gasification cavity was studied under two conditions: perpendicular to (condition 1) and along the intermediate principal stress (condition 2). It is found that under the condition 1, the gasification cavity gradually forms a smaller fracture circle, while the surrounding rock at the floor of the gasification cavity generates a wide range of equivalent damage areas under the condition 2, which is unfavorable for practical engineering. The condition 1 deployment scheme is more practical in terms of the degree of rupture and the subsequent gasification process. The gasification cavity is obviously affected by the horizontal stress. When the horizontal stress is small, the stability of the surrounding rock is seriously damaged. It is necessary to fully consider the influence of in situ stress in the layout of the gasification cavity; at the same time, measures are added to the process to reduce the degree of rupture of the surrounding rock. The evolution law of the temperature field and rupture field of the surrounding pressure in the combustion cavity during the whole process of UCG was numerically simulated. The floor rupture zone develops gradually with the advancement of the working face of the combustion cavity, and deeper rupture zones appear in several areas, which need special attention in the engineering.

1. Introduction

Underground coal gasification (UCG) is a new and environmentally friendly mining method [1], which can produce combustible gas mixture by controlled combustion of coal in situ underground under the action of a series of chemical reactions [2]. It is a mining method that integrates well construction, mining, and gasification, with the characteristics of green mining and clean utilization [3]. In the process of UCG, a key technology is how to effectively evaluate and control the stability of surrounding rocks in the combustion zone [4].

The rock temperature around the combustion cavity formed by UCG is as high as 700~1000° [5]. Under the action of high temperature [6], the surrounding rock of combustion cavity zone generates a large number of microcracks at the boundary of mineral particles due to thermal expansion [7], which seriously affects its stability [8]. Related studies on the evolution of rock fracture damage at high temperatures were also carried out by scholars such as Meng et al. [9], Fan et al. [10], and Ding et al. [11]. Therefore, how to effectively control the stability of the gasification cavity is a key technology in the process of UCG [12].

In essence, the destruction of surrounding rock caused by coal mining is a spatiotemporal evolutionary process [13]. Under the high temperature and in situ stress of coal bed gasification, the surrounding rock around the combustion cavity will form a fracture zone [14]. Once the rupture zone is connected to the underground aquifer, a water permeation accident will be triggered. At the same time, the gas in the gasification cavity may leak or overflow the ground, polluting the environment [15], which will cause the gasification cavity to fail to produce normally or even cause production shutdown accidents [16]. In addition, more serious is the connection between the rupture zone and the underground aquifer that may cause a large amount of groundwater leakage, which will cause great damage to the groundwater resources [17].

As the combustion cavity formed by UCG is often in a three-directional stress state [18], with the expansion of the combustion cavity area, the structural equilibrium of the surrounding rock is broken [19], resulting in the redistribution of the surrounding rock stress field [20]. The combustion cavity is often in a true triaxial stress state, with its stability closely related to the three-dimensional stress environment. Especially, the layout of the gasification cavity will determine the long-term operational stability in the subsequent gasification process. Under the high temperature of coal seam combustion and three-dimensional ground stress, the overlying rocks in the combustion cavity will form fracture zone [21]. When the fracture zone of overlying rocks is connected with the water-bearing layer in the upper part of the coal seam [22], it will cause the accident of water penetration in the roof [23]. At the same time, the gas in the gasification cavity may leak or even spill out of the ground, which may prevent the normal operation of the gasification cavity and even cause a production stoppage [24].

In addition, the research team found that the rock rupture characteristics [25] under the real triaxial stress environment are very different from those of the conventional triaxial [26]. Therefore, the rupture evolution characteristics of combustion cavities under the influence of different stresses are analyzed by combining the three-dimensional ground stress state characteristics [27]. It is of great significance to reveal the rupture mechanism of combustion cavity under the combined effect of high temperature and three-dimensional stress field. Meanwhile, it is important to propose a reasonable gasification cavity layout and construction scheme [28].

In view of this, a controlling equation of rock temperature-stress coupling [29] action is established in this manuscript based on thermodynamic and elastodynamic theories [30]. The numerical calculation model of the combustion cavity expansion under the temperature-stress coupling condition is established by using ABAQUS secondary development program. The law of fracture field evolution during the expansion of the combustion cavity is studied, the influence of different stress fields on the stability of the gasification cavity is analyzed, and a reasonable construction scheme is proposed.

2. UCG Numerical Modeling

A geometric three-dimensional model (Figure 1) was created in ABAQUS to represent the teardrop-shaped characteristics [31] of the combustion cavity under single combustion point conditions [32]. Prabu and Jayanti [33] and Yang et al. [14] found that the gasification cavity formed after combustion was roughly teardrop-shaped; so, this feature will be used for the gasification cavity in this study. The teardrop-shaped combustion cavity is simplified to consist of half a sphere (radius ) and half an ellipsoid (short axis , long axis ). According to the height relationship of the seam, the limit height of the combustion cavity is the full combustion thickness of the coal seam (15 m).

The numerical calculation model uses a rectangular body of , which is simplified to three rock layers: the roof strata, the coal strata, and the floor strata (Figure 1(a)). The boundary conditions of the model are set as follows: the bottom surface is constrained by vertical displacement; the height of the overlying rock layer of the model is 550 m, which is loaded on the upper boundary of the model in the form of a uniform load with a magnitude of 15 MPa (); the maximum and minimum horizontal stresses are 12 MPa and 7 MPa, respectively. Two combustion cavity working conditions are designed: condition 1, perpendicular to the horizontal combustion direction; condition 2, along the horizontal combustion direction. The initial temperature of each rock layer before combustion is 30°C. The temperature of the combustion cavity surface (representing the gasification surface) is raised to 1000°C after the start of gasification.

Natural rocks contain a large number of microdefects such as randomly distributed microfractures inside [34]. A numerical model of rock rupture evolution containing random damage distribution is developed using the secondary development subroutine USDFLD [35]. The total percentage of damage units is , with basic mechanical parameters about 1/2 of the intact units [36]. In addition, considering the influence of temperature, according to the theory of elasticity, the solid equilibrium differential equation that can be expressed by displacement is [37]

where is the displacement component, is the body force component, is the linear thermal expansion coefficient of the rock, is the temperature field, which is determined by the heat conduction control equation, and are the rock damage shear modulus and bulk modulus, and is Poisson’s ratio.

According to the principle of heat balance, it is assumed that the specific heat and thermal conductivity of the rock are constants that do not change with temperature. Then the differential equation of rock heat conduction control can be expressed as [38]

In the formula, is the density of the rock, and are the specific heat capacity and thermal conductivity of the rock, and is the internal heat source of the rock.

Figure 2 shows the numerical calculation process. Since the Mogi-Coulomb strength criterion [39] is more suitable for the true triaxial ground stress environment [40], the strength criterion is used to determine whether the element is damaged [41]. The stiffness degradation of the rock failure unit is treated, and its stiffness is 1/10 of the original initial stiffness [42]. In order to facilitate the calculation, the rock layers are homogenized. The physical and mechanical parameters of each rock layer in the model are shown in Table 1.

Figure 3 shows the fracture characteristics of the surrounding rock of the gasification chamber under the two working conditions. The fracture zone characteristics are similar in both conditions, and the fracture circle is formed gradually near the surface of the gasification cavity. In addition, due to the arched structure of the gasification cavity roof, the fracture zone of the surrounding rock shifts to both sides, resulting in a relatively small fracture zone. Comparing the two working conditions, it can be found that the rupture zone is slightly larger in condition 2 ( along the horizontal combustion direction) than in condition 1.

Since the computational model uses the Mogi-Coulomb intensity criterion, its calculation equation [43] is

Among them, the magnitude of the equivalent effect force has an important pointing significance in judging the rupture trend of the unit; so, the equivalent effect force is used to quantify the rupture trend of the gasification cavity, as shown in Figure 4. The larger the equivalent force is, the more obvious the rupture trend is, and the higher equivalent force (over 10 MPa) is defined as the equivalent damage field in this paper.

Although there is little difference in the range of the rupture zone under the two working conditions, it can be found by comparing the two working conditions: in condition 1, the surrounding rock near the surface of the gasification cavity gradually forms a fracture circle, whose shape characteristics are roughly similar to those of the fracture zone, indicating that the damage characteristics of the surrounding rock can be better quantified by using equivalent forces. In addition, the “” damage field was formed in the two arched shoulders and two bottom corners of the gasification cavity, which poses a potential risk for the continuation of the gasification process.

In condition 2, the surrounding rocks of the gasification cavity produced a wide range of equivalent damage zones, especially the surrounding rocks deep in the floor, which gradually converged with the far-field damage zone formed at the bottom boundary, forming a huge damage penetration zone. This is unfavorable for actual engineering. Once there is an aquifer in the bottom of the floor, it will cause the damage field to be connected to the aquifer, resulting in a substantial water-conducting rupture zone, which will cause a water penetration accident in the floor. The engineering can be avoided by reducing the height of the combustion cavity, the length of the combustion cavity direction, increasing the protection of coal pillars, and other ways to avoid the occurrence of water gushing accident at the bottom plate [44].

3. Horizontal Cross-Sectional Fracture Characteristics of the Combustion Cavity

For a single combustion point, the combustion cavity is approximately teardrop-shaped in the horizontal plane. Since the coal seam has no boundary in the horizontal direction, the combustion cavity expands along the horizontal direction until it reaches the designed combustion cavity width and finally forms a teardrop-shaped combustion cavity. In order to show more clearly the influence of stress field orientation on the rupture characteristics of the combustion cavity, a horizontal section combustion cavity model is established (Figure 5). Comparing the rupture evolution process of the combustion cavity under two working conditions, it can be found in the following:

In condition 1, the rupture zone is asymmetrically distributed near the front and rear ends of the gasification chamber owing to the fact that is perpendicular to the horizontal combustion direction and the asymmetric circular shape of the combustion cavity boundary. As the gasification process proceeds, the rupture zone gradually expands. Since the gasification cavity will continue to retreat in the axial direction of the gasification channel after reaching the combustion boundary, the operating condition can be regarded as an advanced precracking of the surrounding rock of the gasification cavity in the next stage. As long as the fracture range of the coal seams on both sides is effectively controlled, the stability requirements of the surrounding rock of the gasification cavity can be met.

In condition 2, the rupture zone is symmetrically distributed due to along the horizontal combustion direction, mainly concentrated in both sides of the gasification cavity (coal seam). Although the extent of the rupture zone is slightly reduced in this condition compared to condition one, the extent of rupture at the front and rear ends of the gasification cavity is extremely small. There is basically no impact on the next stage of gasification process, which is a waste from the perspective of energy utilization. Of course, it is more advantageous for some projects to use a separate gasification chamber scheme (with lower recovery rate).

4. Fracture Characteristics of the Gasification Cavity in Axial Vertical Section

In order to analyze the influence of different stress states on the fracture characteristics in the axial and vertical section of the gasification cavity, the influence of the horizontal lateral stress on the gasification cavity was studied, and the vertical calculation model shown in Figure 6 was established. The boundary conditions are basically the same as those in Figure 1, where the vertical stress is constant at 15 MPa and the different horizontal stress magnitudes (8,10,12 MPa) are mainly changed. The basic dimensions of the gasification cavity and the physical and mechanical parameters of each rock layer are consistent with Table 1.

Figure 7 shows the effect of different horizontal stresses on the rupture zone of the gasification cavity. When the horizontal stress was small, a wide range of rupture zones appeared near the roof and two bottom corners of the combustion cavity. The stability of the surrounding rock all suffered serious damage, which will have an extremely negative impact on the continued advance of gasification afterwards. With the increase of horizontal stress, this phenomenon is alleviated, especially when the horizontal stress reaches 12 MPa, and the rupture zone of surrounding rock near the combustion cavity is greatly reduced. It shows that the influence of ground stress needs to be fully considered when the gasification chamber is arranged. When the horizontal stress is significantly lower than the vertical stress, a large-scale fracture zone will appear in the surrounding rock near the gasification cavity; so, special attention is needed in the process.

Based on the above analysis results, the stress distribution characteristics of working condition 2 are used to simulate more realistically the development pattern of gasification cavity in the whole process of UCG. In the numerical simulation design, the gasification cavity development is divided into 2 stages according to the thermoelastic instantonal equation and the temperature field control equation (Figure 2). In stage 1, the gasification cavity is small in size and expands equidistantly along the width and height directions in the coal seam until the top of the gasification cavity develops to the roof of the coal seam (1 m/d). The high-temperature gasification surface is gradually transferred to the downstream cavity wall of the gasification cavity in this stage. In stage 2, the gasification working surface is advanced horizontally and uniformly at a constant rate (1 m/d), with a constant height of the gasification chamber.

5. Characteristics of Surrounding Rock Temperature and Rupture Field in the Whole Process of UCG

Based on the above analysis results, the stress distribution characteristics of condition 2 are used to more truly simulate the development law of gasification chamber in the whole process of UCG. In the numerical simulation design, the gasification cavity development is divided into 2 stages according to the thermoelastic instantonal equation and the temperature field control equation (Figure 8). In stage 1, the gasification cavity is small in size and expands equidistantly along the width and height directions in the coal seam until the top of the gasification cavity develops to the roof of the coal seam (1 m/d). The high-temperature gasification surface is gradually transferred to the downstream cavity wall of the gasification cavity in this stage. In stage 2, the gasification working surface is advanced horizontally and uniformly at a constant rate (1 m/d), with a constant height of the gasification chamber.

Figure 9 shows the characteristics of the temperature field distribution around the combustion cavity during the combustion of the coal seam. In stage 1, the space of the combustion cavity is small, and the temperature on the surface of the coal wall and the roof and floor is basically the same, roughly around 950°C. With the burning of the coal wall, the combustion cavity enters stage 2 and advances to the left, with the space gradually increasing. The temperature of the surface of the combustion coal wall basically remains around 950°C during this process, and the temperature of the downstream of the gasification cavity wall gradually decreases, and when the length is 68 m, the temperature of the floor surface downstream has been reduced to around 350°C.

Figure 10 demonstrates the change pattern of temperature of the surrounding rock at different depths of the floor of the coal seam with the UCG process. The temperature on the surface of the floor upstream from the combustion cavity is basically about 950°C, with more drastic temperature changes. With the increase of the vertical distance from the coal seam, the temperature of the rock seam decreases rapidly. For example, the surrounding rock is about 1.5 m away from the floor surface, the temperature field changes gradually, and its maximum peak temperature is roughly half of the floor surface. The temperature field of surrounding rock 3 m from the floor surface has little influence, and the peak temperature is only about 210°C.

Figure 11 shows the evolution of the failure characteristics of the surrounding rock of the combustion cavity during the whole UCG process. In stage 1, the roof of the combustion cavity as well as the two bottom corners generated a wide range of rupture zones, especially the damage of the bottom plate was more obvious. Later, as the gasification process continued, it entered stage 2. A deeper rupture zone appeared near the initial arch on the right side of the gasification cavity, which was prone to roof collapse. The rupture zone of the floor is gradually developed with the advancement of the burning cavity working face, and deeper rupture zones appear in several areas, which need special attention in the engineering to avoid further increase of the rupture zone by leaving coal pillars and other means if necessary.

6. Conclusion

(1)Based on the thermodynamic and elastodynamic theories, the controlling equation of temperature-stress coupling action on rocks containing random damage units is established by combining the Mogi-Coulomb damage criterion. And the numerical calculation model of combustion cavity expansion under temperature-stress coupling condition is established by using ABAQUS secondary development program(2)The fracture field evolution law during the expansion of the gasification cavity was studied under two conditions: perpendicular to (condition 1) and along the intermediate principal stress (condition 2). It is found that under the condition 1, the surrounding rock near the surface of the gasification cavity gradually forms a smaller fracture circle and an “” type damage field appears; under the condition 2, the surrounding rock at the floor of the gasification cavity generates a wide range of equivalent damage areas, which is unfavorable for practical engineering(3)The rupture characteristics of the horizontal section of the gasification chamber under the two working conditions were compared. It is found that the fracture areas of condition 1 are asymmetrically distributed near the front and end of the gasification cavity. In condition 2, the fracture areas are symmetrically distributed and mainly concentrated in the horizontal sides of the gasification cavity. The condition 1 deployment scheme is more practical in terms of the degree of rupture and the subsequent gasification process(4)The gasification cavity is obviously affected by the horizontal stress. When the horizontal stress is small, the stability of the surrounding rock is seriously damaged; with the increase of horizontal stress, the range of the surrounding rock fracture area is greatly reduced. It is necessary to fully consider the influence of in situ stress in the layout of the gasification cavity; at the same time, measures are added to the process to reduce the degree of rupture of the surrounding rock(5)The evolution law of the temperature field and rupture field of the surrounding pressure in the combustion cavity during the whole process of UCG was numerically simulated. The upstream temperature of the gasification cavity shows an increasing trend, and the downstream temperature gradually decreases. A deeper rupture cavity appears near the initial vault on the right side of the gasification cavity, which is prone to roof collapse. The floor rupture zone develops gradually with the advancement of the working face of the combustion cavity, and deeper rupture zones appear in several areas, which need special attention in the engineering

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

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

Acknowledgments

We gratefully acknowledge the support provided by Natural Science Foundation of Jiangsu Province, China (Grant No. BK20200628), China Postdoctoral Science Foundation Funded Project (Grant No. 2020M671649), and the Fundamental Research Funds for the Central Universities (Grant No. 2020QN42).