Abstract
The accuracy of the drainage radius plays a vital role in the gas drainage effect, and the establishment of the drainage model and numerical calculation of the model has an essential value for the accurate determination of the drainage radius. Based on the elastoplastic constitutive model of mining coal and rock mass, effective stress equation, gas content equation, gas flow continuity, and control equation established by predecessors, the gas-solid coupling model of coal-rock deformation and pressure-relief gas flow in protective layer mining was established in this paper. The numerical simulation prototype was based on the mining engineering practice of a short-distance and underprotective steep seam mining in a mine in western Henan of China. COMSOL Multiphysics coupling software was used to numerically calculate the coupling model, and the influence radius of different mining pressure-relief area and different drainage time was obtained. The results show that under the same geological conditions and mining conditions, the impact range of drilling in the pressure-relief area of the protected layer is larger than that of the nonrelief area. As the working face of the protection layer advances, the pressure-relief area of the protected layer gradually increases, the influence range of the drainage borehole increases, the drainage borehole has undergone a process of initial stress-stress concentration-stress reduction-stress recovery successively, and the strike drainage radius is larger than the inclined.
1. Introduction
Mineable reserves of steeply inclined coal seams account for approximately 17% of the total coal reserves. Out of these reserves, 62% are distributed in north China and the rest in south China [1]. The steeply inclined coal seam is recognized as a difficult coal seam in the mining industry. Due to its special occurrence, the steeply inclined coal seam is severely damaged by extrusion, and most of the steeply inclined coal seams are high gas or outburst coal seams. A large number of studies have shown that the mining of protective layer can release the elastic potential of coal seam, increase the permeability of coal seam, facilitate the gas flow and desorption of the protected layer, and reduce the internal energy of coal seam gas. The mining of protective layer combined with pressure-relief gas drainage is considered to be the safest and most effective regional antioutburst measure [2–5]. Protective seam mining technology involves two coal seams: the protective seam and the protected seam. Generally, the protective seam contains no potential outbursts, and the protected seam contains potential outbursts. Many studies in Chinese coalfields have focused on protective seam mining, leading to a large amount of data and achieving good effects on gas extraction [6–9].
Coal and rock mass is a multifractured porous medium composed of solid skeleton and fractured pores, whose pores and fissures are the channels for gas migration and storage [10]. In a pore-fractured coal seam, the pore system absorbs gas and diffuses, while the fracture system seeps free gas. Under normal conditions, the diffusion and seepage are a parallel mass transfer process with the synchronous mass exchange. In the original coal body, the free gas and the adsorbed gas are in a dynamic equilibrium state, and there is constant heat exchange between them and no material transfer on the macroscopic level. The equilibrium state is mainly controlled by the environmental temperature and gas pressure of the coal body. When the coal seam is affected by mining activities, the coal-rock mass is deformed, the pore pressure changes, the original equilibrium state of gas adsorption/desorption is broken, and the coal seam gas adsorption/desorption develops in a single direction, thus forming concentration gradient and pressure gradient in the coal seam, and the resulting gas flow is called pressure-relief gas flow [11–13]. After the protective layer is mined, the overlying rock moves and deforms, the original stress balance of coal and rock is broken, and the original stress balance of the coal-rock mass is broken, which in turn causes the volumetric strain and porosity of the coal-rock mass to be changed, and the coal seam gas pressure and gas flow are changed finally [14, 15]. The change of gas pressure changes the stress state of coal-rock mass and also causes the desorption of adsorbed gas, which further changes the mechanical properties of coal [16].
This study uses the elastoplastic constitutive model of mining coal and rock mass, effective stress equation, gas content equation, gas flow continuity, and control equation established by predecessors, and the gas-solid coupling model of coal-rock deformation pressure-relief gas flow in protective layer mining was established. Furthermore, COMSOL Multiphysics coupling software was used to solve the coupling model of pressure-relief gas drainage in underprotective steep seam mining.
2. Gas-Solid Coupling Model
2.1. Basic Assumption
The research object of the gas-solid coupling model of coal-rock deformation and pressure-relief gas flow is the coupling effect between adsorption/desorption, diffusion, seepage movement of pressure-relief gas, and deformation movement of coal and rock mass. The study of gas-solid coupling law is a complicated problem, involving many disciplines such as geotechnical mechanics and seepage mechanics [17, 18]. In order to facilitate numerical solution and analysis, the following basic assumptions are put forward: (1)The gas-solid coupling system comprises single-phase gas and solid-phase coal-rock. The stress-strain relationship conforms to the mining elastic-plastic constitutive model, and the mining coal-rock mass is regarded as a fractured pore medium(2)The coal body cannot be hardened and is subject to correlation and regular yield loading and unloading criteria. M-c criterion is used in this paper. According to this criterion, when the pressure is not large (generally <10 MPa), the strength criterion can be expressed as follows:where is shear stress (MPa); is cohesion (MPa); is normal stress (MPa); and is internal friction angle. (3)The coal seam contains free gas and adsorbed gas, and the relationship between gas content and pressure obeys the Langmuir equation [19]:
The free gas flow state obeys the ideal gas state equation, namely: where is the mass of adsorbed gas per unit volume of coal seam, that is, the mass concentration of adsorbed gas diffusion in the coal seam (g/mL); is the mass concentration of free gas (g/mL); is the limit gas adsorption capacity of coal (mL/g); is the adsorption constant (MPa-1); is the coal seam gas pressure (MPa); is the coal bulk density (g/mL); is the gas molar mass (CH4 is 16 g/mol); is the gas under standard conditions molar volume (22.4 L/mol); is the coal porosity (%); and p0 is the standard atmospheric pressure (0.1 MPa). (4)The adsorption of gas is nonequilibrium adsorption(5)The diffusion of gas in the coal seam conforms to Fick’s law of diffusion, and the diffusion coefficient remains unchanged [20]:where is the diffusion flux of gas through unit area [kg/(s.m2)]; is the diffusion coefficient of coal seam pore system (m2/s); is the mass of adsorbed CBM contained in unit volume of coal seam, namely the diffusion mass concentration of coalbed methane in the adsorbed state in the coal seam (kg/m3); is the Hamiltonian operator; and is the unit vector in the coordinate directions, respectively. (6)The seepage of gas in fractures conforms to Darcy’s law, and the permeability of coal and rock mass varies with the porosity and mining stress of coal and rock mass:where is the seepage velocity vector of the free gas (volume flux vector m/s); is the permeability of the coal seam fracture system (m2); is the gas dynamic viscosity (CH4 is ); p is the free gas in the fracture system state CBM pressure (Pa); and is the unit vector in the coordinate direction, respectively. (7)At the same level, the original gas content and gas pressure of the coal seam are the same(8)Gas is a positive pressure fluid, and its mass density is only related to its own pressure
According to the continuum theory, the establishment of the gas flow coupling model of mining pressure-relief needs to be based on the gas flow control equation, the gas flow continuity equation, the gas state equation of gas flow, the gas content equation of coal seam, and the coal-rock mass deformation equation [21].
2.2. Gas Flow Continuity Equation
2.2.1. Continuity Equation of Diffusion Motion of Pore System
In the pore system, take a microelement, and in time , according to the principle of mass conservation, the mass of the inflowing microelement minus the sum of the outflow mass plus the generation amount of the mass source should be equal to the unit time. The mass change of the microelement, for the pore system, q is a negative exchange negative mass source, kg/(m3.s), and then the continuous equation of the diffusion motion of the pore system is [22, 23]
2.2.2. Continuity Equation of Seepage Movement in Fracture System
Taking a microelement in the fracture system, after time , the gas density and porosity will change, which will cause the control of the fluid quality in the body to change. According to the principle of mass conservation, the mass change of the microelement includes the mass of the seepage flowing into the microelement minus the sum of the outflow mass and the input amount. The continuity equation of the seepage motion of the fracture system is
2.2.3. Coupling Effect between Gas Diffusion and Seepage
The above research is carried out by dividing the flow of coal seam gas into diffusion and seepage as two open systems. The adsorption/desorption of coal seam gas and gas diffusion are regarded as a system. Gas desorption converts the adsorbed gas on the inner wall of micropores into free gas, which diffuses into the seepage system along with the pores, and the two systems complete the mass transfer through the mass exchange source. Considering gas diffusion and seepage as a semiclosed system, this system has only output, and gas diffusion and seepage will inevitably decrease concentration. Gas diffusion and seepage are controlled by concentration, influence each other, and control each other. Therefore, there is a coupling effect between gas diffusion and seepage. Since the gas adsorption and desorption only change the material form, it does not cause the change of the gas quality in the microelement. Therefore, in order to facilitate the research, analysis, and solution, the gas flow mass conservation equation is obtained by adding equations (8) and (9):
2.3. Adsorption Expansion Strain
Most of the gas in the coal seam is in the adsorption state, and the size of the pore surface area in the coal mass and the strength of the attractive force between the coal molecules and between the coal molecules and the gas molecules directly affect the surface tension of the coal body and the gas adsorption capacity [24]. When the gas molecules are adsorbed by the coal, the surface tension of the coal decreases, and the volume expands and deforms. The volume expansion rate of the coal is proportional to the amount of adsorption, and the coal generates expansion stress due to the adsorption expansion under the constraint conditions [25]. Assuming that coal is an elastic body, the mechanical properties and adsorption properties are the same in all directions, the adsorption properties are not affected by external forces, the adsorption deformation of the coal skeleton contact point is in a unidirectional compression state, and because the pore pressure is also the same in all directions, the coal expands in all directions According to the relationship between the adsorption expansion stress-strain and the adsorption thermodynamic parameters, the expansion stress-strain can be expressed by the following formula [26–28]: where is the adsorption expansion stress (MPa); is the Poisson’s ratio; is the adsorption expansion strain; is the adsorption expansion strain increment; is the universal gas constant, ; is the absolute temperature (K); and is bulk modulus. The other symbols represent the same meaning as the above.
2.4. Relationship between Effective Stress and Pore Pressure and Adsorption Expansion Stress
Due to adsorption expansion stress, the effective stress calculation changes. The common structural model, adsorption deformation mode, coal mass deformation mode, and coal seam stress state of gas and coal all have an impact on the connection between the two [29, 30]. The relationship between the effective stress of gas-bearing coal body, pore pressure, and adsorption expansion stress is as follows: where is the effective stress of gas-bearing coal body; is the total stress; is the pore pressure; and is the adsorption expansion stress.
The existence of pore pressure and adsorption expansion stress can reduce the effective stress of coal, but there are essential differences between the two [31]. The adsorption expansion stress is caused by the expansion deformation caused by the decrease of the surface energy of the coal skeleton after adsorbing the gas. The force transmission must rely on the coal skeleton, and the existence of cracks causes the reverse stress difference in the crack space due to the lack of skeleton filling [32, 33]. Even when the effective stress remains unchanged, the adsorption expansion stress still acts on the fracture surface, which leads to the closure of the fracture and affects the permeability characteristics of the fractured coal body. On the one hand, the adsorption expansion stress reduces the effective stress and promotes permeability. On the other hand, the adsorption expansion stress will lead to the closure of cracks, decrease permeability, and inhibit gas flow.
2.5. Equations of Porosity and Permeability of Gas-Bearing Coal
The predecessors have carried out many studies on the evolution models of porosity and permeability under adsorption expansion deformation and pore pressure. The representative models of porosity and permeability are as follows: (1) the Seidle-Huitt (S-H) model only considers the shrinkage effect of the coal matrix [34]; (2) the Palmer-Mansoori (P-M) model assumes that the deformation of the coal body is a uniaxial strain condition, which is consistent with the (SH) model, and coal permeability and porosity adopt the cubic law [35]; (3) in the Shi-Durucan (S-D) model, (SD) model and (P-M) model assume the same [36]; (4) Cui-Bustin (C-B) model deduces the change of porosity with stress and gas pressure according to the relationship between pore volume and stress [37]; (5) Robertson-Christiansen (R-C) model assumes that the coal body is under isobaric conditions, and the cubic model; on this basis, it is deduced that the permeability varies with gas pressure [38]. Among them, the most common one is the (P-M) model. The model parameters are relatively simple and easy to obtain. Under the same assumptions, the (C-B) model can also be further simplified to the (P-M) model. During the mining process of the protective layer, the movement and deformation of the overlying rock and the drainage of the pressure-relief gas from the protected layer will disturb the protected coal seam, and the original stress and pore pressure of the protected layer will change. Then the effective stress of the coal and rock mass will change. Combined with PM model and C-B model, according to the principle of effective stress of coal body, the influence of coal matrix adsorption expansion/desorption and shrinkage deformation and mining stress change, the evolution model of porosity and permeability can be obtained as follows: where is the Lame constant, ; is the bulk modulus of pores and fissures, usually 0.6~0.8 times the bulk modulus of coal; is the initial porosity; and is the initial permeability.
2.6. Deformation Equation of Mining Coal and Rock Mass
It is assumed that the stress state of any microelement coal-rock mass is shown in Figure 1.
For any microelement as shown in Figure 1, the stress balance equation in the three directions of , , and can be expressed as [39]
Equation (17) can be expressed as a tensor:
By substituting the effective stress equation (14) into equation (18), the equilibrium differential equation for the elastic-plastic deformation of the coal-rock skeleton is obtained as
The geometric equation of coal and rock mass reflects the relationship between strain and displacement, and the tensor is in the form of where is the strain tensor and is the displacement tensor.
In the study of gas-solid coupling in porous media, the deformation of coal and rock mass is controlled by effective stress, and the effective stress-strain constitutive relationship of coal and rock mass satisfies the constitutive equation of generalized Hooke’s law [40]: where is the coal-rock skeleton strain and is the elastic constitutive symmetry matrix of the rock. In three-dimensional coordinates, is a square matrix of order 6. For an isotropic orthotropic body, the elastic modulus and Poisson’s ratio can be used to represent
3. Numerical Simulation
3.1. Numerical Simulation Experimental Prototype
In this paper, the numerical simulation is solved with the help of COMSOL Multiphysics software, and the engineering and geological background is mainly based on the mining of No. 1-8 coal protected by the protected layer No. 2-1 coal seam in a mine in west Henan province, China. The average thickness of No. 1-8 coal seam is 1.1 m, and the average thickness of the vertical level is 0.71 m; the average thickness of the lead plumb of the protected layer of No. 2-1 coal seams is 6.2 m, and the average thickness of vertical level is 4.0 m; the average dip angle of coal seams is 50°, the lead plumb layer spacing of No. 2-1 coal and No. 1-8 coal is 14.0 m, and the dip distance between layers is about 9.0 m. The relative layer spacing is about 12 times the mining height of the protected layer. Its coal stratum is dominated by sandstone, mudstone, quartz sandstone, and coal seam, and the distribution of coal stratum and lithology is shown in Figure 2.
3.2. Physical Modeling and Meshing
The numerical simulation of protective layer working face along the inclined layout and along the strike mining establishes the numerical physical model of the inclined and strike, respectively. The length of the working face of the protection layer in strike direction is 270 m, with 20 m coal pillars at each end of the model, the total length of the model is 310 m, the height is 81.2 m, and 117,540 meshes are divided; the length of the working face in inclined direction is 75 m, the width of the model is 88.47 m, the height is 91.43 m, and 37,770 meshes are divided; the physical model and mesh division are shown in Figures 3 and 4.
3.3. Boundary Conditions
The average mining depth of the coal seam is 330 m, the vertical stress is 8.25 MPa, and the horizontal stress is 1.5 times of the vertical stress, which is 12.38 MPa. Only gas source exists inside the coal seam, and there is no gas flow field outside. The coal seam mining area is set as fixed as the first type of boundary conditions and fixed as atmospheric pressure 0.1 MPa, and fixed negative pressure 13 KPa is set in the borehole space.
3.4. Initial Equilibrium of Ground Stress
The basic parameters of the coal-rock body are assigned, and the rock mechanics parameters of the coal-rock body are shown in Table 1. These parameters can present the distribution and equilibrium state of stress in coal strata under the state of original in-situ stress.
3.5. Gas Initial Balance
Combined with the actual mine, the initial pore pressure, permeability, and porosity parameters of the coal seam are initially assigned, and the basic parameters of the initial state gas are shown in Table 2. The stress and gas in Table 2 are the parameters in the initial state.
3.6. Coupling Calculation
The protective layer working face is retrieved according to the actual control time step in the field, and the coupled calculation of gas seepage, stress, and coal-rock body deformation is performed.
4. Results and Analysis
4.1. The Gas Pressure Change Rule of the Protected Layer along the Strike Direction
Due to the large size of the strike model, if the boreholes are arranged exactly according to the actual borehole spacing in the field, the simulation results will be too dense and difficult to distinguish. In order to more intuitively observe the transport of unloading gas in the protected layer under the action of extraction, only 5 groups of drainage boreholes (boreholes are numbered from left to right as 1#~5#) are arranged in the protected layer in the strike profile; the borehole spacing is 60 m, and the borehole diameter is 100 mm. Figure 5 shows the cloud diagram of gas pressure distribution in the protected layer when 30 m, 60 m, 90 m, 120 m, 150 m, and 270 m are, respectively, mined at the working face of the protected layer.
It is evident from Figure 5 that compared with the drainage boreholes that are not subject to mining pressure-relief by the lower protective layer, the drainage influence range of the boreholes within the pressure-relief area of the protected layer is larger, and the gas pressure around the boreholes is smaller. With the advancement of the protected layer’s working face, the protected layer’s unloading range gradually increases, and the influence range of the drainage borehole also increases. In addition, under the influence of mining stress, the drainage borehole undergoes the original stress-stress concentration-stress reduction-stress recovery and finally stabilizes. Figure 6 shows the curve drawn from the gas pressure of the protected seam, and the bottom of the sharp decline of the curve is the position of the extraction borehole on strike.
(a) 30 m
(b) 60 m
(c) 90 m
(d) 120 m
(e) 150 m
(f) 270 m
The law of pressure-relief gas drainage is as follows. (1)When the lower protection layer of No. 1-8 coal is retrieved for 30 m, it can be seen from Figure 6(a) that only borehole #1 is in the pressure-relief area, borehole #2 is in the stress concentration area in front of the working face of the protection layer, and boreholes #3~#5 are in the original stress area; and the horizontal distance of borehole #2 from the working face of the protection layer is 45 m at this time. The effective drainage radius of borehole #1 in the pressure-relief area is 3.8 m, the effective drainage radius of borehole #2 in the stress concentration area is only 1.4 m, and the effective drainage radius of boreholes #3~5 in the original stress area is only 1.8 m(2)From Figure 6(b), it can be seen that when the protective layer is retrieved for 60 m, only borehole #1 is still in the pressure-relief area, borehole #2 is in the area with higher stress concentration, and boreholes #3~5 are still in the original stress area. After 90 days of extraction, the effective drainage radius of borehole #1 in the pressure-relief area is 4.4 m, the effective drainage radius of borehole #2 in the pressure-relief area is reduced to 1.3 m, and the effective drainage radius of boreholes #3-5 in the original stress area is 1.9 m(3)When the protective layer is retrieved for 90 m, borehole #1 is still in the pressure-relief area, borehole #2 also enters the pressure-relief area, borehole #3 enters the stress concentration area, and boreholes #4 and #5 are still in the original stress area, as shown in Figure 6(c). After 105 days of drainage, the effective drainage radii of boreholes 1~5# were 4.8 m, 4.9 m, 1.6 m, 2.0 m, and 2.0 m in order(4)As the protective layer working face advances, the relief area also gradually increases. When the recovery distance increases to a certain extent, the collapsed rocks in the goaf behind the protective layer working face are gradually compacted. The coal body of the protected layer originally in the relief area is recompacted. The stress is gradually restored, at which time the permeability of the coal body of the protected layer decreases compared with the previous relief, and the drainage effect also decreases(5)When the working face is retrieved for 150 m, it can be seen from Figures 5(e) and 6(e) that borehole #2 is exactly in the stress recovery area, and boreholes #1 and #3 are near the boundary of the pressure-relief area, borehole #4 is in the stress concentration area, and borehole #5 is in the original stress area. At this time, the effective drainage radius of boreholes 1~5# is 4.9 m, 3.5 m, 5.7 m, 2.7 m, and 3.0 m in order, and the maximum gas pressure in the stress recovery area has increased compared with the previous one, recovering to 1.27 MPa, which is slightly lower than the original gas pressure. The permeability of the corresponding areas in descending order is pressure-relief area > stress recovery area > original stress area > stress concentration area(6)When the protective layer No. 1-8 coal is retrieved 270 m along the strike, it can be seen from Figures 5(f) and 6(f) that the mining relief range covers the area where boreholes 1#~5# are located. The pressure-relief is the most evident at the two ends of the strike open-cut eye and the stopping line, and the central part is recompacted as the stress recovery area. There is no stress concentration area directly above the mining area without leaving coal pillars in the mining area. In 200 days of drainage, the maximum effective drainage radius of the decompression area at both ends reaches 6.2 m, and the effective drainage radius of the central compaction area is only 3.9 m
4.2. The Gas Pressure Change Rule of Protected Layer along the Inclined Direction
The inclined model size is relatively small, in the inclined profile with a spacing of 15 m the same arrangement of 5 groups of boreholes, and borehole diameter of the same 100 mm. Considering the simulation of the working face along the inclined arrangement, along the direction of recovery, inclined on the protective layer working face once finished, inclined is protected layer gas pressure distribution as shown in Figure 7.
The drainage time of 60 days, 90 days, 150 days, and 200 days was selected to protect the coal seam gas pressure along the inclined direction and plotted as a curve (Figure 8), and the bottom of the curve dropping sharply is the location of the drainage borehole on the inclined.
From Figures 7 to 8, it can be seen that the degree of pressure-relief at the inclined upper part of the tilt is significantly higher than that at the lower part, and the drainage influence range gradually increases with the increase of drainage time. The variation of the pressure-relief drainage radius with drainage time is shown in Table 3, and a logarithmic relationship was found when the two were fitted (as shown in Figure 9).
When the protective layer workings are mined over, the fissures around the protective layer drainage area can be maintained for a long time, which is the full pressure-relief area and can be used as the main action area for gas drainage. After the overlying rock layer is moved and stabilized, a larger stress recovery area will appear in the middle of the protected layer. The lower part of the inclined is less depressurized than the upper part, and the spacing of drainage boreholes in this area should be shortened appropriately. In addition, from the comparison of inclined and strike, the inclined drainage radius is smaller than the strike. When arranging drainage boreholes, the spacing between boreholes in the inclined direction should be smaller than that in the strike direction. Intensive drilling should be carried out in the protection and expansion boundary areas, especially in the inclined lower boundary areas.
To sum up, in the process of protected seam mining, the coal body of the protected seam all undergoes the process of “stress concentration→pressure-relief→stress recovery,” and the space and time for sufficient pressure-relief should be grasped to drain the pressure-relief gas from the protected seam efficiently.
5. Conclusions
(1)On the one hand, the adsorption expansion stress reduces the effective stress and promotes the permeability. On the other hand, the adsorption expansion stress will lead to the closure of cracks, decrease the permeability, and inhibit the flow of gas(2)Compared with the drainage boreholes that are not affected by the mining pressure-relief effect of the protective layer, the drainage boreholes within the pressure-relief area of the protective layer have a larger influence, and the gas pressure around the boreholes is smaller. With the advance of the working face of the protective layer, the pressure-relief area of the protected layer gradually increases, and the influence range of the drainage borehole also increases. In addition, under the influence of mining stress, the drainage borehole has undergone a process of initial stress-stress concentration-stress reduction-stress recovery(3)When coal mining is completed, the pressure-relief is most obvious at both ends of the strike opening hole and stop-mining line, and the central recompaction is the stress recovery area. Under the condition that no coal pillar is left in the goaf, there is no stress concentration area directly above the goaf. After 200 days of drainage, the maximum effective drainage radius of the pressure-relief area at both ends is 6.2 m, and the effective drainage radius of the central compaction area is only 3.9 m. The degree of pressure-relief in the inclined upper part is higher than in the lower part. With the increase of drainage time, the influence range of drainage gradually increases, and the relationship between drainage radius and drainage time is logarithmic. The maximum effective drainage radius of the inclined pressure-relief area is 4.9 m at 200 days of drainage(4)The fissures around the goaf of the protective layer can be maintained for a long time. In order to fully relieve the pressure, it can be used as the main area for gas drainage, and the middle of the protected layer is the stress recovery area. The pressure-relief degree of the inclined lower part is smaller than that of the upper part. From the comparison between the inclined and strike, the radius of the inclined extraction is smaller than that of the strike. When the boreholes are arranged, the spacing of the boreholes in the inclined direction should be smaller than that in the strike direction. Intensive boreholes should be carried out in the protected and extended boundary areas, especially in the inclined lower boundary areas
Data Availability
The test data used to support the findings of this study are included within the article. Readers can obtain data supporting the research results from the test data table in the paper.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (Nos. 52074107 and 51704100); the Doctoral Fund of Henan Polytechnic University (No. B2016-004); the Postdoctoral Research Foundation of Henan Province (Third-Class Financial Support & First-Class Financial Support); the Henan Provincial Key Laboratory of Gas Geology and Gas Control-Open Project of the State Key Laboratory Cultivation Base jointly built by the province and the ministry (No. WS2018B13); the Fundamental and Frontier Technology Research Project of Henan Province (No. 162300410038); and the Key Scientific Research Project of Higher Education Institutions of Henan Provincial Department of Education (No. 15A440007).