Abstract

Gas drainage through long seam boreholes is an effective method to prevent gas disasters in coal mines. In this paper, a multifield coupling model of gas migration in gas-bearing coal was first established. Then, a quantitative characterization method of gas drainage effect was put forward. Finally, the extraction effect of long boreholes was obtained under different layouts and drilling lengths. The research results show that, under the arrangement of long boreholes along the seam, the gas pressure around the borehole decreases significantly with the extraction time. There is no extraction blank in the middle of the working face. However, it is easy to cause uneven gas drainage in the combined arrangement of the long boreholes along the seam and the penetrating boreholes. Furthermore, it is found that the drainage volume of the long boreholes along the seam is similar to that of the joint layout under the same drainage time. As the length of the borehole increases, the influencing range of gas drainage increases. When the borehole lengths are 150 m and 240 m, the drainage volumes are about 1.31 and 2.50 times that of the 90 m boreholes, respectively. The research achievements could provide a specific reference for the layout of long boreholes along the bedding and the determination of reasonable parameters for gas drainage on site.

1. Introduction

Gas drainage is an effective way to eliminate the danger of coal and gas outbursts in China at present [14]. In recent years, with the development of drilling rig equipment and drilling technology, gas drainage technology with long boreholes has made rapid progress [57]. Gas migration significantly affects the gas drainage effect of long boreholes. Therefore, a correct understanding of the law of coal gas migration is of great significance for preventing and controlling coal mine gas disasters and promoting long drilling technology [812].

At present, many scholars have studied the multifield coupling model of coal seam gas migration. The coal reservoir is a typical double pore system composed of interlaced fractures and coal matrices [1315]. Coal seam gas migration mainly includes two processes: diffusion and seepage. First, the gas is desorbed from the coal matrix and diffused into the fissures under the action of a concentration gradient. This process conforms to Fick’s law of diffusion. Then, the gas in the fissure flows to the borehole under the action of the pressure gradient, and this process conforms to Darcy’s law [16]. Su et al. [17] described a dual pore model that paralleled coal seam gas migration. The diffusion and seepage of gas in gas drainage are simulated, and the influence of diffusion coefficient on gas pressure evolution and permeability in fractures is analyzed. Dong et al. [18] analyzed the mechanism of negative pressure in gas drainage through a gas-solid coupling model and put forward a method to reduce the negative pressure of drainage to improve gas drainage efficiency. Mora and Wattenbarger [19] summarized the correct shape factor formulas considering constant drainage rate from a matrix block and constant pressure in the adjacent fractures through numerical simulation. However, the above research mainly analyzed the gas drainage effect in conventional borehole gas drainage. There are few studies on gas drainage with long boreholes along the seam.

On the other hand, to evaluate the effect of long borehole gas drainage, it is necessary to study the reasonable length of the borehole for long borehole drainage [2022]. If the length of the borehole is too large, due to the attenuation of the negative pressure in the borehole, the gas content at the bottom of the long borehole will be high, making it difficult to reach the drainage index. If the length of the borehole is too short, it will result in low efficiency of gas drainage measures [2325]. At present, the methods to determine the effective gas drainage radius include on-site measurement and theoretical calculation methods [26, 27]. The on-site measurement mainly determines the residual gas content of the gas after the gas drainage. In terms of theoretical calculations, the effective gas drainage radius is determined mainly based on the gas flow theory in the borehole and the characteristics of coal seam gas seepage. [28, 29] The practice of predecessors in gas drainage has achieved remarkable results. In the case of a short borehole length, the negative pressure of gas drainage in the borehole does not decrease significantly, and there is still a large negative pressure for gas drainage at the bottom of the borehole. After a period of time of extraction, there is little difference in residual gas content before and after the borehole. However, for long boreholes along the bedding, when the negative pressure attenuation is apparent, the residual gas content at the bottom of the borehole is still relatively large. Therefore, it is necessary to study the reasonable length of the borehole for long borehole drainage.

This paper first elaborates and establishes a multifield coupling model of gas migration in coal which considers coal matrix gas diffusion, fissure gas seepage, permeability evolution law, and coal deformation law. Then, the change value of the residual gas content per unit time is used to calculate the gas drainage volume and then quantify the gas drainage effect of the long borehole. Finally, the COMSOL software was adopted, and the fully coupled numerical simulation method was used to analyze the drainage effects of different layouts and lengths of long boreholes along the bedding. The research results are expected to provide a theoretical basis for gas drainage from long boreholes along the bedding.

2. Multifield Coupling Model for Gas-Bearing Coal Seams

2.1. Model Assumptions

The multifield coupling model for gas-bearing coal seams is established based on the following assumptions:(a)Coal seam is a homogeneous, isotropic, double-porosity elastic medium(b)The influence of water on gas migration in coal seams is ignored(c)The coal seam is an isothermal system, and gas is an ideal gas(d)The strain of the coal skeleton is infinitesimal [14, 16, 17]

2.2. Effective Stress Principle

In order to better describe the mechanical response of pore and fissure dual media under gas pressure, the effective stress law of dual porous media is introduced when establishing the permeability model [16]:where is the effective stress, MPa; is the total stress, MPa; is the Kronecker delta (1 for i = j and 0 for i ≠ j); and are pressure in the fractures and matrix blocks, respectively, Pa; and βf and βm are the effective stress coefficients for pore and fracture, respectively.

The effective stress coefficients βf and βm of the cracks in (1) can be calculated by the following formulas, respectively:where K is the bulk modulus of the coal body, MPa; Km is the bulk modulus of the coal matrix, MPa; and Ks is the bulk modulus of the coal body skeleton, MPa.

The following formulas can calculate the above three bulk moduli:where E is the elastic modulus of the coal body, MPa; Em is the elastic modulus of the coal matrix, MPa; is Poisson’s ratio of the coal; and is the coal matrix porosity, %.

2.3. Dynamic Evolution Equation of Coal Porosity and Permeability

Palmer and Mansoori [30] proposed a widely used permeability model (PM model) suitable for uniaxial strain conditions. However, the PM model was established based on the coal seam containing only cracks. The coal seam is usually modeled as a dual-porosity medium containing both cracks and pores. Therefore, this study adopts the dual-porosity poroelasticity theory. The revised crack porosity model is shown in the following equation:where is the fracture porosity, %; is the initial fracture porosity, %; is Langmuir volumetric strain constant; is Langmuir pressure constant; M is the limiting axial modulus, MPa; and is the initial pressure in the matrix blocks, Pa.

According to the Kozeny–Carman equation [31, 32], the coal seam permeability k iswhere is the coal seam permeability, m2.

2.4. Diffusion Equation of Gas in Coal

The mass exchange between the matrix and the fissures in the coal body can be expressed by the following equation:where Qm is the mass exchange between coal pores and cracks, kg/(m3·s); is the matrix shape factor, m−2; MC is the molar mass of methane, kg/mol; D is gas diffusion coefficient, m2/s; R is the universal gas constant, J/(mol·K); and T is the coal seam temperature, K.

The following formula can calculate the diffusion coefficient:where ξ is the attenuation coefficient of the dynamic diffusion coefficient, s−1, and is the initial gas diffusion coefficient, m2/s.

According to the Langmuir equation, the gas content per unit volume of coal matrix can be obtained aswhere mp is the gas content per unit mass of the coal matrix, kg/m3; ρc is the apparent density of the coal, kg/m3; is the matrix porosity, %; a is the Langmuir volume constant, m3/kg; b is the reciprocal of Langmuir pressure constant, Pa−1; and Vm is the molar volume of methane under standard conditions, m3/mol.

From the conservation of mass, it can be seen that the relationship between the amount of change in matrix gas content and Qm is

Substituting (6) ∼ (8) into (9), the governing equation for gas diffusion is

2.5. The Seepage Equation of Gas in Coal

According to the law of conservation of mass, the change in gas mass in a fracture is equal to the gas diffused into the fracture by the pore minus the gas flowing into the borehole from the fracture; namely,where is the gas flow velocity in the crack, m/s.

Moreover, the flow of gas in the cracks of the coal body conforms to Darcy’s law; then,where μ is the dynamic viscosity of the gas, Pa·s.

By combining the above equations, the governing equation of gas percolation can be obtained as

2.6. Deformation Equation of Gas-Bearing Coal

The deformation equation of gas-containing coal is composed of the stress balance equation, geometric deformation equation, and stress-strain relationship [33], which can be expressed aswhere G is the coal’s shear modulus, MPa, and ui is the displacement component in the i-direction.

The above control equations constitute a multifield coupling model for gas-containing coal seams that comprehensively consider the coal skeleton compression effect (effective stress effect) and matrix shrinkage effect. The coupling relationship between the physical fields is shown in Figure 1.

3. Geometric Model and Definite Solution Conditions

3.1. Geometric Model

In order to analyze the effect of gas drainage with long boreholes along the seam, two gas drainage schemes under two conditions have been established :(1)Gas drainage with long boreholes along the seamThe geometric model is shown in Figure 2. The 240 m long boreholes along the seam were constructed in the direction of the coal seam’s inclination in the machine roadway or wind roadway in the working face. The borehole spacing distance is 2 m, and ten long boreholes along the seam were constructed.(2)The combination of long boreholes along the seam and penetrating boreholesThe geometric model is shown in Figure 3. The 90 m long boreholes along the seam were constructed in the machine lane and the winding lane in the working face along the direction of the inclination of the coal seam. The borehole spacing distance is 2 m, and 20 long boreholes along the seam were constructed. Moreover, the penetrating boreholes were constructed in the middle of the working face of the bottom-drawing roadway, and the position of the final hole of the penetrating boreholes is shown in Figure 3. The boreholes are spaced 2 m apart in both the strike and inclination directions, and a total of 290 penetrating boreholes have been constructed.The size of the coal seam under these two drainage schemes is 400 × 240 m.

The numerical simulation software used in this study is COMSOL Multiphysics. The simulation uses the PDE module and the solid mechanics module in COMSOL to calculate the model. The parameters used in the numerical simulation are shown in Table 1 [17].

3.2. Model Validation

To verify the correctness of the multi-field coupling model of gas-bearing coal seams established in this paper, the use of mathematical models of Liu et al. [16] coal seam gas pressure data was compared to verify. A model consisting of three horizontal coal layers was established. There is a coal seam with a width of 40 meters and a height of 4 meters between the two rock layers. The drill hole is located in the center of the coal seam. The other parameters used in the model are consistent with those in the Liu model. After COMSOL numerical calculation, the experimental results are shown in Figure 4. It can be seen from Figure 4 that the gas-solid coupling model established in this paper is in good agreement with Liu’s model.

4. Quantitative Characterization of Gas Drainage Effect in Boreholes

It can be seen from Section 1 that, at time t and t+ Δt, the gas content in the coal seam is

In the formula, mt and mt+ Δmt represent the gas content in the coal seam when the drainage time is t and t+ Δt, respectively.

According to the above formula, the gas content extracted from the coal seam after the gas drainage of Δt time is as follows:

Based on the above formula, the gas content extracted within a certain period can be obtained by integration with COMSOL software.

5. Results and Discussion

5.1. Influence of Drilling Layout on the Effect of Gas Drainage

Figure 5 is a cloud map of coal seam gas content under different drainage times in long boreholes along the seam. It can be seen from Figure 5 that, under the action of negative pressure in the borehole, the gas pressure around the borehole is significantly lower than that of the distant coal body. Furthermore, because the borehole spacing distance is only 2 m, the gas pressure between the drill holes is significantly reduced, and there is no apparent blank zone. In addition, the length of the borehole reaches 240 meters and runs through the entire working face. Therefore, the contact area between the borehole and the coal body is large, and the drainage effect is good. On the other hand, as the drainage time increases, the gas content around the borehole gradually decreases, and the influence range of borehole drainage gradually increases. For example, after 360 days of drainage, the gas content of the coal body covered by the borehole was significantly reduced to 3 m3/t, which was significantly lower than the gas content of 5 days of drainage.

At present, due to limited drilling rig equipment in some domestic mining areas, the construction length of long boreholes along the seam is limited, so there must be a blank zone of gas drainage in the middle of the working face. To ensure the safety of the working face during the mining period, penetrating boreholes are usually used to extract the coal seam gas in the middle of the working face. The boreholes arrangement in this section is based on the above facts.

Figure 6 is a cloud map of coal seam gas content under different drainage times in combined drainage with long boreholes along the seam and penetrating boreholes. It can be seen from Figure 6 that the distribution of residual gas content in coal seams is different from that in Figure 5. The gas content in long boreholes along the seam drainage area is significantly lower than that of penetrating boreholes drainage area. For example, at five days of drainage, the gas content in the area covered by the long boreholes along the seam has dropped to 4 m3/t, while the gas content in the penetrating boreholes covering the coal seam is about 8 m3/t. According to the analysis in Section 1, the gas drainage volume per unit time of a borehole is mainly determined by the permeability and the contact area of the borehole. Therefore, the gas drainage effect of long boreholes along the seam is significantly higher than that of penetrating boreholes. In addition, with the increase of the extraction time, the gas content of the coal body around the borehole is consistent with the previous analysis. One point that needs to be particularly pointed out is that, with the increase of the drainage time, the difference between the gas drainage effect of the long boreholes along the seam and the penetrating boreholes gradually decreases.

For example, at 360 days of drainage, the gas content in the long boreholes along the seam area is about 2.5 m3/t, and the gas content in the penetrating boreholes area is also reduced to about 3 m3/t.

It can be seen that the drainage effect of the long boreholes along the seam arrangement is significantly higher than that of the combined long boreholes along the seam and penetrating boreholes. However, it is difficult to intuitively compare the drainage effects of the two drilling arrangements based on the distribution cloud map of the residual gas content after coal seam drainage. To solve the problems mentioned above, this paper proposes a quantitative characterization method for the effect of borehole gas drainage in Section 3. Based on this quantitative characterization method, the gas drainage volume under different borehole layouts can be obtained. This paper takes the combined drainage volume of long boreholes along the seam and penetrating boreholes as the benchmark. It can be concluded that there is an increase in drainage under the condition of long boreholes along the seam relative to the combined drainage, as shown in Figure 7. It can be seen from Figure 7 that the gas drainage volume in the bedding long borehole layout is similar to the gas drainage volume under the combined layout under the same drainage time. For example, in the early stage of drainage, the drainage volume under the long borehole layout is 1.025 times that of the combined drainage, and it gradually decreases in the later stage, down to 0.94 times at the lowest level. In addition, considering the construction cost of tunnel extraction and drilling through layers, it can be considered that the long boreholes along the seam can reduce the pressure and content of coal seam gas faster in the process of gas drainage in coal mines. The danger of gas outbursts in the working face is eliminated, and construction costs are reduced.

Figures 8 and 9 are the streamline diagram of the long boreholes along the seam and the streamline diagram of the combined boreholes through bedding during the 360-day drainage. It can be seen from Figure 8 that the gas migration streamlines of the long boreholes along the seam are perpendicular to the borehole. The gas flows neatly and orderly to the boreholes, with almost no influence on each other, conforming to radial flow. It can be seen from Figure 9 that the gas streamlines in the long boreholes along the seam area still show a radial flow, while the gas flow in the penetrating boreholes area is messy. In addition, the turbulent gas streamlines interfere with each other at the intersection of the long boreholes along the seam and penetrating boreholes. The gas migration between boreholes is subject to drainage forces in multiple directions, not conducive to gas drainage.

The above analysis compares and analyzes the drainage effect of the arrangement of long boreholes along the seam and the combined arrangement of penetrating boreholes and the long boreholes along the seam from the three aspects of gas drainage volume, gas flow line, and drilling construction cost. It can be seen that long boreholes along the seam have advantages in all aspects. Therefore, it is recommended that the arrangement of long boreholes along the seam is used in gas drainage in the working face.

5.2. Influence of Drilling Length on the Effect of Gas Drainage

To study the impact of the length of the long borehole along the bedding on gas drainage, numerical simulations of gas drainage with borehole lengths of 90 m, 120 m, 150 m, 180 m, 210 m, and 240 m were carried out. The geometric model is shown in Figure 10.

When the borehole length is 90 m, the gas content distribution around the borehole at different extraction times is shown in Figure 11. It can be seen from Figure 11 that, under the action of the borehole negative pressure, as the drainage time increases, the gas content around the borehole gradually decreases, and the range of influence of borehole drainage gradually increases.

To quantitatively analyze the distribution of gas content around the borehole under different extraction times, a monitoring line is selected along the direction and inclination of the working face, and the gas content change law on it is analyzed. The monitoring line in the strike direction is 45 m away from the roadway wall, and the monitoring line in the inclined direction is first 5 m away from the borehole. Figure 12 shows the distribution of gas content on the monitoring line under different extraction times. From Figure 12(a), when the borehole length is 90 m, in the strike direction, as the drainage time increases, the gas content of the coal seam gradually decreases. At a distance of 5 m from the borehole, the gas content at a drainage time of 180 days is reduced by 1.0 m³/t and 1.3 m³/t compared with the gas content at a drainage time of 30 days and five days. In the direction of the trend, the gas content of the coal seam gradually decreases with the increase of the extraction time. However, because the length of the borehole is only 90 m, part of the gas in the middle and back of the coal seam has not been effectively drained.

When the borehole length is 150 m, the gas content distribution around the borehole under different extraction time is shown in Figure 13. It can be seen from Figure 13 that, under the action of borehole negative pressure, as the drainage time increases, the gas content around the borehole gradually decreases, and the range of influence of borehole drainage gradually increases. Compared with the 90 m borehole, the 150 m borehole has a broader range of drainage, and the drainage effect is more significant.

In the strike direction, as the drainage time increases, the coal seam gas content gradually decreases. At a distance of 5 m from the borehole, the gas content at 180-day drainage is reduced by 0.94 m³/t and 1.26 m³/t compared with the gas content at the drainage time of 30 days and five days, and the change is small compared to the 90 m borehole, as shown in Figure 14.

When the borehole length is 240 m, the gas content distribution around the borehole under different extraction time is shown in Figure 15. It can be seen from Figure 15 that, under the action of borehole negative pressure, as the drainage time increases, the gas content around the borehole also shows a gradual decrease, and the influence range of borehole drainage further increases.

In the strike direction, as the drainage time increases, the coal seam gas content gradually decreases. In the trend direction, the gas content in the coal seam gradually decreases with the increase of the extraction time, and the gas content at the bottom of the borehole decreases significantly. With the increase of drilling depth, the upward trend value changes less. As the drilling depth increases, the gas at the bottom of the drilling hole can be drained, and the gas content decreases significantly, as shown in Figure 16.

To quantitatively analyze the impact of borehole length on the gas drainage effect, the drainage volume when the borehole length is 90 m is used as the benchmark. The increase in the drainage of the 150 m and 240 m borehole lengths relative to the 90 m borehole can be obtained, as shown in Figure 17. It can be seen from Figure 17 that the longer the borehole length, the better the drainage effect and the greater the impact range of the borehole. When the borehole length is 150 m, the drainage volume is about 1.31 times that of the 90 m borehole; when the borehole length is increased to 240 m, the drainage volume is about 2.50 times that of the 90 m borehole. It can be seen from Figure 17 that the flow line after drilling is perpendicular to the drilling, the gas migration at the tip of the slot conforms to the law of spherical flow, and the gas migration in the middle and rear parts conforms to the radial flow. With the increase of drilling length, the streamline area and the influence range of drilling increase significantly.

6. Conclusions

This paper first elaborates and establishes a multifield coupling model of gas migration in coal. Then, the change value of the residual gas content per unit time is used to quantify the gas drainage effect of the long borehole. Finally, the COMSOL software was adopted to analyze the drainage effects of different layouts and drilling lengths of long boreholes along the bedding. The main conclusions are as follows:(1)A multifield coupling model for gas-containing coal seams was established. The model considers the gas diffusion in the coal matrix, the gas seepage through fractures, the evolution of permeability, and the deformation of coal. At the same time, the dynamic equation of coal diffusion is introduced to make the simulation results more in line with the actual situation.(2)By integrating the space of the coupled model and then making the difference in time, the change value of the residual gas content per unit time can be obtained to calculate the gas drainage volume. Furthermore, the gas drainage effect of long boreholes can be quantified.(3)Under the combined layout of long boreholes along the seam and penetrating boreholes, the gas content of the long boreholes along the seam area is significantly lower than that of the interbedded borehole drainage area.(4)The drainage volume is similar to that under the combined arrangement under the same drainage time. In the early stage of drainage, the drainage volume under the long borehole layout is 1.025 times that of combined drainage, and it gradually decreases in the later stage, down to 0.94 times at the lowest level.(5)As the length of the borehole increases, the scope of influence of borehole drainage increases. When the borehole length is 150 m and 240 m, the drainage volume is about 1.31 and 2.50 times that of the 90 m borehole.

The above results provide a particular reference for the layout of long boreholes along the bedding and the determination of reasonable parameters for gas drainage on site.

Data Availability

All the data, models, and code generated or used during the study appear in the manuscript.

Conflicts of Interest

The author declares no commercial or associative interest that represents conflicts of interest in connection with the work submitted.

Acknowledgments

The research presented in this paper was jointly supported by the Natural Science Foundation of Chongqing (Grant No. cstc2020jcyj-msxmX1013).