# A Case Study on the Gas Drainage Optimization Based on the Effective Borehole Spacing in Sima Coal Mine

Authors: Ming Ji; Zhong-guang Sun; Wei Sun
Journal: Geofluids
 (2021)

---

## Abstract

Based on the dynamic expressions of permeability and porosity of the coal seam derived in the paper, a multiphysical field coupling numerical model of gas migration under the interaction of stress field and seepage field was established. The gas drainage project #3 Coal Seam operated by Sima Coal Industry Co., Ltd., was selected as the study object. Taking different drainage time periods in various positions of drainage holes into consideration, combined with the advance situation of the 1207 working face in the Sima Coal Mine, a mixed layout gas drainage scheme featured with the effective borehole spacing was obtained through the COMSOL multiphysics simulation. In addition, a series of field industrial tests were performed to validate the research result, revealing that comprehensively considering the extraction time of coal and optimizing the layout of extraction boreholes can effectively improve the engineering economic benefits.


---

## Body

## 1. Introduction

Coal seam gas drainage is one of the important measures to control mine gas [1, 2] However, a majority of researches take the effective drainage radius as the basis of drainage hole layout, overlooking the influences of the superimposed drainage and fracture expansion without taking the influence of different drainage times of drainage holes in different positions of the coal seam into account. Additionally, most researches adopt a single borehole layout, which often causes some common issues such as a waste of resource, uneven gas drainage, and failure of meeting relevant standards. Many existing studies fail to consider the establishment of a multiphysical field coupling model of gas migration under the interaction of the stress field and seepage field and overlook the dynamic nature of the permeability and porosity of the coal seam. All the aforementioned which are lacking in previous studies have made it necessary to conduct an in-depth optimization research on coal seam gas drainage.In this paper, the gas drainage optimization of #3 Coal Seam in Sima Coal Industry Co., Ltd., of Lu’an Group was studied. Based on the changes of coal seam stress, gas pressure, and gas adsorption and desorption which tend to exert an important impact on the coal seam permeability and porosity, the dynamic expressions of coal seam permeability and porosity were deduced which was used to build the numerical model of a multiphysical field coupling for migration under the interaction of stress field and seepage field. Additionally, taking different drainage time periods in various positions of drainage holes into consideration, combined with the advance situation of the 1207 working face in the Sima Coal Mine, a mixed layout gas drainage scheme featured with the effective borehole spacing was obtained by the COMSOL multiphysics simulation. In addition, a series of field industrial tests were performed for the purpose of validating the research.


## 2. Gas Solid Coupling Analysis of Coal Seam Gas Migration

### 2.1. The Migration Law of Coal Seam Gas

Two types of gas flow in coal seams have been identified including the diffusion flow and the laminar flow targeting different types of coal seams [3–5]. Generally speaking, in the coal seam with a microporous structure or low permeability, the diffusion flow plays a critical role while in the coal seam with microporous structure or high permeability, the laminar flow dominates [6–8]. In this study, the gas flow in the process of drainage is treated as laminar flow, following Darcy’s law.A number of factors can affect the migration of coal seam gas, such as geological conditions of storage, occurrence of coal seam, mechanical properties of coal and rock, permeability and porosity of coal seam, gas adsorption capacity of coal seam, and gas pressure. Among those factors, gas pressure, gas adsorption capacity of coal seam, porosity, and permeability of coal seam are considered as the primary ones. Gas pressure is the driving force of gas migration while the gas adsorption capacity of coal seam decides the capacity of the coal seam to store gas. Meanwhile, the porosity and permeability are the indicators of smooth gas migration in the coal seam [9–11].The gas pressure of the #3 Coal Seam in Sima mine was measured at 0.31 MPa which was used as the basic gas pressure of this study, despite that the gas pressure measured stayed below the critical value of 0.74 MPa. The gas adsorption constants a and b of the #3 Coal Seam were identified at about 17.52 m3/t and 0.75 MPa-1, separately, which were considered low, indicating that the gas adsorption capacity of the #3 Coal Seam was weak. During the field drainage, several measures can be adopted to facilitate the drainage and extraction including the selection of the appropriate negative pressure of the drainage hole, the optimization of the layout of the hole, the promotion of the desorption of gas, and the increase of the gas diffusion and migration speed.


### 2.2. The Derivation of Dynamic Expressions of Porosity and Permeability of Coal Seam

During the gas extraction, the gas tends to be desorbed and migrated, exposing the coal skeleton to stress changes, resulting in small deformation, and the porosity and permeability change accordingly. According to the existing research and the field conditions of Sima mine and the 2D model in the numerical simulation calculation, the plane strainεsis introduced, leading to a modified dynamic expression of coal seam porosity [12–14].
(1)φ=1−1−φ01+εs1−∆pks,where φ0 is the initial porosity of coal, εv is the volume strain, ∆P is the variation of gas pressure, ∆p=P−P0, P is the real time gas pressure, P0 is the initial gas pressure, and ks is the skeleton modulus of the coal and rock mass.According to the Kozeny-Carman formula, combined with the previous dynamic expression of porosity, the modified dynamic expression of permeability can be obtained as follows:(2)ke=k01+εs1+εs+∆p/ks1−φ0φ03,wherek0 is the initial permeability of coal and εs is the plane strain.


### 2.3. The Multiphysical Field Coupling Model of Coal Seam Gas Migration

This study focuses on the deformation of the coal seam due to the forces imposed and the gas flow [15–18]. Therefore, the following assumptions have been made:
(1)
The flow field formed by the gas drainage hole under a certain negative pressure is the radial flow field(2)
The adsorption content of the gas is described by the Langmuir equation(3)
Gas flow is described by Darcy’s law(4)
The seepage process of the gas is regarded as an ideal isothermal process(5)
The coal seam is isotropic with only a small linear elastic deformation(6)
The influence of water in the coal seam on gas drainage is overlookedTaking the assumptions listed above into consideration, the seepage movement of coal seam gas as the ideal gas should follow the control equations of gas flow listed as follows [13, 14, 19–25].


#### 2.3.1. Law of Conservation of Mass

Gas in the coal seam often takes two forms including the free state and the adsorption state. Despite the statuses, the source and total amount of gas in the coal seam remain unchanged. In other words, the gas migration should first follow the law of conservation of mass:(3)∂m∂t+∇∙ρgqg=Qp,where m is the gas content (kg/m3), ρg is the gas density (kg/m3), qg is the Darcy seepage velocity of gas (m/s), Qp is the source or sink items (kg/(m3·s)), and t is the time variable (s).


#### 2.3.2. The State Equation of Gas

The gas in seepage movement is regarded as an ideal gas, which satisfies the equation of state of ideal gas:(4)ρg=MgPRT,where Mg is the molecular weight of gas (kg/kmol), P is the gas pressure (Pa), R is the ideal gas constant (kJ/(kmol·K)), and T is the absolute temperature (Ko).


#### 2.3.3. Langmuir Equation

The coal seam is a porous medium with many pores and fissures, in which gas coexists in the adsorption state and the free state. The Langmuir equation describing the gas content defines the relationship between gas adsorption and desorption with gas pressure:(5)m=MgRTφP0+abρs1+bPP2,where φ is the porosity, P0 is the atmospheric pressure (Pa), a is the Langmuir constant (m3/kg), b is the Langmuir constant (Pa−1), and ρs is the coal seam density (kg/m3).


#### 2.3.4. Darcy’s Law

When a pressure gradient is identified in the distribution of gas pressure in the coal seam, the gas migration will occur, which is often a linear seepage flow and follows Darcy’s law:(6)qg=−keμg∇p+ρgg∇z2,where qg is the Darcy flow velocity (m/s), ke is the permeability of the coal seam (m2), μg is the gas dynamic viscosity (Pa·s), and g is the acceleration of gravity (m·s−2).The process of gas extraction is featured with two physical phenomena including the seepage movement of gas as fluid and the microdeformation of coal. The stress and deformation of the coal body cause changes to the pore characteristics and then affect the gas seepage movement. Consequently, the gas seepage tends to change the gas pressure within the coal body, resulting in different effective stresses. The interaction of two phenomena acts like a chain reaction featured with mutual constraint and influence.The dynamic expression of coal seam porosity is substituted into equation (5), then
(7)m=MgRT1P0−1−φ0P01+εs1−∆Pks+abρs1+bPP2.The dynamic expression of coal seam permeability is substituted into equation (6), then
(8)qg=−k0∇p+ρgg∇z2μg1+εs1+εs+∆p/ks1−φ0φ03.By substituting formula (4) and derived formulas (7) and (8) into formula (3), the following results are obtained:
(9)∂Mg/RT1/P0−1−φ0/P01+εs1−∆P/ks+abρs/1+bPP2∂t−∇·MgPk0∇p+ρgg∇z2RTμg1+εs1+εs+∆p/ks1−φ0φ03=Qp.Formula (9) is the gas-solid coupling model finally derived in this paper. The model is embedded into the COMSOL multiphysics software for gas drainage optimization.


## 2.1. The Migration Law of Coal Seam Gas

Two types of gas flow in coal seams have been identified including the diffusion flow and the laminar flow targeting different types of coal seams [3–5]. Generally speaking, in the coal seam with a microporous structure or low permeability, the diffusion flow plays a critical role while in the coal seam with microporous structure or high permeability, the laminar flow dominates [6–8]. In this study, the gas flow in the process of drainage is treated as laminar flow, following Darcy’s law.A number of factors can affect the migration of coal seam gas, such as geological conditions of storage, occurrence of coal seam, mechanical properties of coal and rock, permeability and porosity of coal seam, gas adsorption capacity of coal seam, and gas pressure. Among those factors, gas pressure, gas adsorption capacity of coal seam, porosity, and permeability of coal seam are considered as the primary ones. Gas pressure is the driving force of gas migration while the gas adsorption capacity of coal seam decides the capacity of the coal seam to store gas. Meanwhile, the porosity and permeability are the indicators of smooth gas migration in the coal seam [9–11].The gas pressure of the #3 Coal Seam in Sima mine was measured at 0.31 MPa which was used as the basic gas pressure of this study, despite that the gas pressure measured stayed below the critical value of 0.74 MPa. The gas adsorption constants a and b of the #3 Coal Seam were identified at about 17.52 m3/t and 0.75 MPa-1, separately, which were considered low, indicating that the gas adsorption capacity of the #3 Coal Seam was weak. During the field drainage, several measures can be adopted to facilitate the drainage and extraction including the selection of the appropriate negative pressure of the drainage hole, the optimization of the layout of the hole, the promotion of the desorption of gas, and the increase of the gas diffusion and migration speed.


## 2.2. The Derivation of Dynamic Expressions of Porosity and Permeability of Coal Seam

During the gas extraction, the gas tends to be desorbed and migrated, exposing the coal skeleton to stress changes, resulting in small deformation, and the porosity and permeability change accordingly. According to the existing research and the field conditions of Sima mine and the 2D model in the numerical simulation calculation, the plane strainεsis introduced, leading to a modified dynamic expression of coal seam porosity [12–14].
(1)φ=1−1−φ01+εs1−∆pks,where φ0 is the initial porosity of coal, εv is the volume strain, ∆P is the variation of gas pressure, ∆p=P−P0, P is the real time gas pressure, P0 is the initial gas pressure, and ks is the skeleton modulus of the coal and rock mass.According to the Kozeny-Carman formula, combined with the previous dynamic expression of porosity, the modified dynamic expression of permeability can be obtained as follows:(2)ke=k01+εs1+εs+∆p/ks1−φ0φ03,wherek0 is the initial permeability of coal and εs is the plane strain.


## 2.3. The Multiphysical Field Coupling Model of Coal Seam Gas Migration

This study focuses on the deformation of the coal seam due to the forces imposed and the gas flow [15–18]. Therefore, the following assumptions have been made:
(1)
The flow field formed by the gas drainage hole under a certain negative pressure is the radial flow field(2)
The adsorption content of the gas is described by the Langmuir equation(3)
Gas flow is described by Darcy’s law(4)
The seepage process of the gas is regarded as an ideal isothermal process(5)
The coal seam is isotropic with only a small linear elastic deformation(6)
The influence of water in the coal seam on gas drainage is overlookedTaking the assumptions listed above into consideration, the seepage movement of coal seam gas as the ideal gas should follow the control equations of gas flow listed as follows [13, 14, 19–25].


### 2.3.1. Law of Conservation of Mass

Gas in the coal seam often takes two forms including the free state and the adsorption state. Despite the statuses, the source and total amount of gas in the coal seam remain unchanged. In other words, the gas migration should first follow the law of conservation of mass:(3)∂m∂t+∇∙ρgqg=Qp,where m is the gas content (kg/m3), ρg is the gas density (kg/m3), qg is the Darcy seepage velocity of gas (m/s), Qp is the source or sink items (kg/(m3·s)), and t is the time variable (s).


### 2.3.2. The State Equation of Gas

The gas in seepage movement is regarded as an ideal gas, which satisfies the equation of state of ideal gas:(4)ρg=MgPRT,where Mg is the molecular weight of gas (kg/kmol), P is the gas pressure (Pa), R is the ideal gas constant (kJ/(kmol·K)), and T is the absolute temperature (Ko).


### 2.3.3. Langmuir Equation

The coal seam is a porous medium with many pores and fissures, in which gas coexists in the adsorption state and the free state. The Langmuir equation describing the gas content defines the relationship between gas adsorption and desorption with gas pressure:(5)m=MgRTφP0+abρs1+bPP2,where φ is the porosity, P0 is the atmospheric pressure (Pa), a is the Langmuir constant (m3/kg), b is the Langmuir constant (Pa−1), and ρs is the coal seam density (kg/m3).


### 2.3.4. Darcy’s Law

When a pressure gradient is identified in the distribution of gas pressure in the coal seam, the gas migration will occur, which is often a linear seepage flow and follows Darcy’s law:(6)qg=−keμg∇p+ρgg∇z2,where qg is the Darcy flow velocity (m/s), ke is the permeability of the coal seam (m2), μg is the gas dynamic viscosity (Pa·s), and g is the acceleration of gravity (m·s−2).The process of gas extraction is featured with two physical phenomena including the seepage movement of gas as fluid and the microdeformation of coal. The stress and deformation of the coal body cause changes to the pore characteristics and then affect the gas seepage movement. Consequently, the gas seepage tends to change the gas pressure within the coal body, resulting in different effective stresses. The interaction of two phenomena acts like a chain reaction featured with mutual constraint and influence.The dynamic expression of coal seam porosity is substituted into equation (5), then
(7)m=MgRT1P0−1−φ0P01+εs1−∆Pks+abρs1+bPP2.The dynamic expression of coal seam permeability is substituted into equation (6), then
(8)qg=−k0∇p+ρgg∇z2μg1+εs1+εs+∆p/ks1−φ0φ03.By substituting formula (4) and derived formulas (7) and (8) into formula (3), the following results are obtained:
(9)∂Mg/RT1/P0−1−φ0/P01+εs1−∆P/ks+abρs/1+bPP2∂t−∇·MgPk0∇p+ρgg∇z2RTμg1+εs1+εs+∆p/ks1−φ0φ03=Qp.Formula (9) is the gas-solid coupling model finally derived in this paper. The model is embedded into the COMSOL multiphysics software for gas drainage optimization.


## 2.3.1. Law of Conservation of Mass

Gas in the coal seam often takes two forms including the free state and the adsorption state. Despite the statuses, the source and total amount of gas in the coal seam remain unchanged. In other words, the gas migration should first follow the law of conservation of mass:(3)∂m∂t+∇∙ρgqg=Qp,where m is the gas content (kg/m3), ρg is the gas density (kg/m3), qg is the Darcy seepage velocity of gas (m/s), Qp is the source or sink items (kg/(m3·s)), and t is the time variable (s).


## 2.3.2. The State Equation of Gas

The gas in seepage movement is regarded as an ideal gas, which satisfies the equation of state of ideal gas:(4)ρg=MgPRT,where Mg is the molecular weight of gas (kg/kmol), P is the gas pressure (Pa), R is the ideal gas constant (kJ/(kmol·K)), and T is the absolute temperature (Ko).


## 2.3.3. Langmuir Equation

The coal seam is a porous medium with many pores and fissures, in which gas coexists in the adsorption state and the free state. The Langmuir equation describing the gas content defines the relationship between gas adsorption and desorption with gas pressure:(5)m=MgRTφP0+abρs1+bPP2,where φ is the porosity, P0 is the atmospheric pressure (Pa), a is the Langmuir constant (m3/kg), b is the Langmuir constant (Pa−1), and ρs is the coal seam density (kg/m3).


## 2.3.4. Darcy’s Law

When a pressure gradient is identified in the distribution of gas pressure in the coal seam, the gas migration will occur, which is often a linear seepage flow and follows Darcy’s law:(6)qg=−keμg∇p+ρgg∇z2,where qg is the Darcy flow velocity (m/s), ke is the permeability of the coal seam (m2), μg is the gas dynamic viscosity (Pa·s), and g is the acceleration of gravity (m·s−2).The process of gas extraction is featured with two physical phenomena including the seepage movement of gas as fluid and the microdeformation of coal. The stress and deformation of the coal body cause changes to the pore characteristics and then affect the gas seepage movement. Consequently, the gas seepage tends to change the gas pressure within the coal body, resulting in different effective stresses. The interaction of two phenomena acts like a chain reaction featured with mutual constraint and influence.The dynamic expression of coal seam porosity is substituted into equation (5), then
(7)m=MgRT1P0−1−φ0P01+εs1−∆Pks+abρs1+bPP2.The dynamic expression of coal seam permeability is substituted into equation (6), then
(8)qg=−k0∇p+ρgg∇z2μg1+εs1+εs+∆p/ks1−φ0φ03.By substituting formula (4) and derived formulas (7) and (8) into formula (3), the following results are obtained:
(9)∂Mg/RT1/P0−1−φ0/P01+εs1−∆P/ks+abρs/1+bPP2∂t−∇·MgPk0∇p+ρgg∇z2RTμg1+εs1+εs+∆p/ks1−φ0φ03=Qp.Formula (9) is the gas-solid coupling model finally derived in this paper. The model is embedded into the COMSOL multiphysics software for gas drainage optimization.


## 3. The Introduction of the Calculation Model

### 3.1. The Calculation Model Parameter Setting

Before establishing the numerical model, the determination of the relevant parameters involved is vital. Based on relevant literature and geomechanics tests, after testing and adjusting the coefficient of the COMSOL multiphysics 5.0 built-in module, the main parameters used in the model were obtained, as shown in Table1.Table 1
Model parameter.
ItemsSymbolValueGas density (kg/m3)rhog0.716Gas dynamic viscosity (Pa·s)vl1.08e-5Equivalent compressibility of matrix (1/Pa)χp2.18e-3Fluid compressibility (1/Pa)chif1.3942e-5Coal density (kg/m3)rhoc1470Initial porosity of coalφ00.0318Poisson’s ratio of coalmu0.33Initial permeability of coal (m2)k00.907e-15Young’s modulus of coal (MPa)E3000Initial pressure of coal gas (MPa)P00.31Standard atmospheric pressure (Pa)Pa1.01e5Biot-Willis coefficientαB0.801


### 3.2. The Model Establishment

The original gas pressure of the #3 Coal Seam in the Sima mine was 0.31 MPa. Due to the large error of gas concentration measurement, the original gas pressure reduced by 30% was selected as the standard for the simulation, leading to setting the gas drainage standard pressure at 0.22 MPa.In order to provide direct guidance to the field construction, the effective borehole spacing was introduced in the simulation. The gas pressure of the whole seam can be reduced to the maximum borehole spacing value within the standard pressure value after drainage for a certain period of time. In other words, through the numerical simulation, the most reasonable drilling spacing value can be determined, which not only meets the need of reducing gas pressure but also reduces the difficulty and cost of on-site construction.The method similar to the “bucket theory” was used to calculate the effective drilling spacing, aiming to find out the area with the least superimposed effect of borehole pumping in the simulation calculation model and ensure that the pressure within the area drops right below the standard value.The porous elastic module and Darcy’s law module [13, 14, 26] in COMSOL were used in this simulation, and the required variables and parameters were inputted to the software. According to the actual drilling gas drainage process, an idealized two-dimensional calculation model with length×height of 16m×6.6m was established as shown in Figure 1. The drainage boreholes were arranged in the middle of the model with the bottom as a fixed constraint. Meanwhile, the upper, left, and right boundaries have a boundary load of 8.25 MPa as the actual confining pressure condition. All sides of the model were airtight with an initial gas pressure of 0.31 MPa in the model.Figure 1
The calculation model.As demonstrated in Figure1, a drill hole was arranged in the 2D calculation module. Following the single variable method, the effective extraction radius corresponds to various diameters of extraction drill holes (95 mm, 113 mm, and 133 mm), extraction pressures (20 kPa, 25 kPa, 30 kPa, 35 kPa, 40 kPa, and 45 kPa), and extraction durations (90 d, 180 d, and 360 d). The detailed calculations are presented in Table 2.Table 2
The simulation results of single-bore model.
Single boreEffective radius of drainage (m)Hole diameterNegative pressure90 d180 d360 d95 mm20 kPa0.450.570.7295 mm25 kPa0.460.570.7395 mm30 kPa0.460.580.7495 mm35 kPa0.470.590.7595 mm40 kPa0.470.600.7695 mm45 kPa0.480.610.77113 mm20 kPa0.480.600.76113 mm25 kPa0.480.610.77113 mm30 kPa0.490.620.78113 mm35 kPa0.490.620.79113 mm40 kPa0.500.630.80113 mm45 kPa0.510.640.81133 mm20 kPa0.500.630.80133 mm25 kPa0.510.640.81133 mm30 kPa0.510.650.82133 mm35 kPa0.520.660.83133 mm40 kPa0.530.670.84133 mm45 kPa0.530.670.85According to Table2, as the extraction duration, pressure, and diameter of the extraction drill hole increased, the effective radius of the single drainage hole increased, demonstrating a positive correlation.Currently, the commonly adopted extraction drill holes are sized at113 mm in diameter and 133 mm in diameter, with the extraction pressure between 25 kPa and 40 kPa. The extraction often consumes substantial time, more than 1 year and even 2 to 3 years. Due to the relatively thicker coal seal in the Sima mine, multiple extraction holes in multiple layers were adopted. The extraction holes are often arranged in three ways including the three-flower style, four-flower style, and five-flower style, as demonstrated in Figure2.Figure 2
Three arrangements of drainage holes.A few simulations with various conditions were performed including the extraction drill holes of 113 mm and 133 mm in diameter under the extraction pressure of 25 kPa, 30 kPa, 35 kPa, and 40 kPa. The extraction duration was set at 360 d. The calculation of the effective extraction radius is demonstrated in Figure3 below.Figure 3
The effective distances between bores under different suction pressures.Figure 4
The sections divided in the roadway for extraction.According to Figure3, the drainage hole of 133 m in diameter demonstrated better extraction. The negative pressure of the extraction hole was 35kpa in the coal mining practice. The five-flower layout of the extraction holes was superior than the three-flower layout while the three-flower layout is superior the four-flower layout.


## 3.1. The Calculation Model Parameter Setting

Before establishing the numerical model, the determination of the relevant parameters involved is vital. Based on relevant literature and geomechanics tests, after testing and adjusting the coefficient of the COMSOL multiphysics 5.0 built-in module, the main parameters used in the model were obtained, as shown in Table1.Table 1
Model parameter.
ItemsSymbolValueGas density (kg/m3)rhog0.716Gas dynamic viscosity (Pa·s)vl1.08e-5Equivalent compressibility of matrix (1/Pa)χp2.18e-3Fluid compressibility (1/Pa)chif1.3942e-5Coal density (kg/m3)rhoc1470Initial porosity of coalφ00.0318Poisson’s ratio of coalmu0.33Initial permeability of coal (m2)k00.907e-15Young’s modulus of coal (MPa)E3000Initial pressure of coal gas (MPa)P00.31Standard atmospheric pressure (Pa)Pa1.01e5Biot-Willis coefficientαB0.801


## 3.2. The Model Establishment

The original gas pressure of the #3 Coal Seam in the Sima mine was 0.31 MPa. Due to the large error of gas concentration measurement, the original gas pressure reduced by 30% was selected as the standard for the simulation, leading to setting the gas drainage standard pressure at 0.22 MPa.In order to provide direct guidance to the field construction, the effective borehole spacing was introduced in the simulation. The gas pressure of the whole seam can be reduced to the maximum borehole spacing value within the standard pressure value after drainage for a certain period of time. In other words, through the numerical simulation, the most reasonable drilling spacing value can be determined, which not only meets the need of reducing gas pressure but also reduces the difficulty and cost of on-site construction.The method similar to the “bucket theory” was used to calculate the effective drilling spacing, aiming to find out the area with the least superimposed effect of borehole pumping in the simulation calculation model and ensure that the pressure within the area drops right below the standard value.The porous elastic module and Darcy’s law module [13, 14, 26] in COMSOL were used in this simulation, and the required variables and parameters were inputted to the software. According to the actual drilling gas drainage process, an idealized two-dimensional calculation model with length×height of 16m×6.6m was established as shown in Figure 1. The drainage boreholes were arranged in the middle of the model with the bottom as a fixed constraint. Meanwhile, the upper, left, and right boundaries have a boundary load of 8.25 MPa as the actual confining pressure condition. All sides of the model were airtight with an initial gas pressure of 0.31 MPa in the model.Figure 1
The calculation model.As demonstrated in Figure1, a drill hole was arranged in the 2D calculation module. Following the single variable method, the effective extraction radius corresponds to various diameters of extraction drill holes (95 mm, 113 mm, and 133 mm), extraction pressures (20 kPa, 25 kPa, 30 kPa, 35 kPa, 40 kPa, and 45 kPa), and extraction durations (90 d, 180 d, and 360 d). The detailed calculations are presented in Table 2.Table 2
The simulation results of single-bore model.
Single boreEffective radius of drainage (m)Hole diameterNegative pressure90 d180 d360 d95 mm20 kPa0.450.570.7295 mm25 kPa0.460.570.7395 mm30 kPa0.460.580.7495 mm35 kPa0.470.590.7595 mm40 kPa0.470.600.7695 mm45 kPa0.480.610.77113 mm20 kPa0.480.600.76113 mm25 kPa0.480.610.77113 mm30 kPa0.490.620.78113 mm35 kPa0.490.620.79113 mm40 kPa0.500.630.80113 mm45 kPa0.510.640.81133 mm20 kPa0.500.630.80133 mm25 kPa0.510.640.81133 mm30 kPa0.510.650.82133 mm35 kPa0.520.660.83133 mm40 kPa0.530.670.84133 mm45 kPa0.530.670.85According to Table2, as the extraction duration, pressure, and diameter of the extraction drill hole increased, the effective radius of the single drainage hole increased, demonstrating a positive correlation.Currently, the commonly adopted extraction drill holes are sized at113 mm in diameter and 133 mm in diameter, with the extraction pressure between 25 kPa and 40 kPa. The extraction often consumes substantial time, more than 1 year and even 2 to 3 years. Due to the relatively thicker coal seal in the Sima mine, multiple extraction holes in multiple layers were adopted. The extraction holes are often arranged in three ways including the three-flower style, four-flower style, and five-flower style, as demonstrated in Figure2.Figure 2
Three arrangements of drainage holes.A few simulations with various conditions were performed including the extraction drill holes of 113 mm and 133 mm in diameter under the extraction pressure of 25 kPa, 30 kPa, 35 kPa, and 40 kPa. The extraction duration was set at 360 d. The calculation of the effective extraction radius is demonstrated in Figure3 below.Figure 3
The effective distances between bores under different suction pressures.Figure 4
The sections divided in the roadway for extraction.According to Figure3, the drainage hole of 133 m in diameter demonstrated better extraction. The negative pressure of the extraction hole was 35kpa in the coal mining practice. The five-flower layout of the extraction holes was superior than the three-flower layout while the three-flower layout is superior the four-flower layout.


## 4. The Simulation and Analysis of the 1207 Working Face of Sima Coal

### 4.1. The Simulation Process Analysis

The strike length of the 1207 working face in the #3 Coal Seam of the Sima mine is 220 m with an inclined length of 1092 m. During the production, in order to ensure the safety of production, gas extraction was carried out. However, this kind of pumping scheme is single without considering the specific circumstances.As mentioned previously, the drainage time greatly affects the gas drainage. Different drainage time periods should adopt different drainage hole layout schemes, such as different borehole spacings or patterns. As the no. 1207 working face advanced, different extraction durations were identified in different positions of extraction boreholes in the mining roadway. Under such circumstance, different sections should be divided based on the different extraction durations, allowing individual extraction design corresponding to different field conditions.According to the comprehensive analysis of the driving speed, working face layout time, and mining speed of the mining roadway in the no. 1207 working face, and taking the drainage effectiveness into consideration, the coal seam gas drainage in the no. 1207 working face can be divided into three sections in simulation, as shown in Figure4. The average extraction time of the three sections is 110D, 210D, and 310D, respectively.The average extraction time, negative pressure, and hole diameter of the three sections were identified, and the gas extraction optimization of the three sections was also studied in detail. The calculation model of the optimization of the extraction changes from the single-row hole layout to the five-row pattern hole layout until the optimal result was obtained. Under the extraction of the experimental scheme, the goal of the whole layer of the #3 coal extraction up to 0.22 MPa was achieved.


### 4.2. The Analysis of Simulation Results

#### 4.2.1. The Analysis of the First and Middle Sections

The extraction time of the first section and the middle section was set at 110 days and 210 days, respectively. According to the simulation calculation, the coal seam gas pressure can be reduced to the standard value within a fixed time only when the five-flower layout was adopted. Due to the different extraction durations, the effective drilling spaces of the two sections varied. The gas pressure cloud diagram in the simulation process is shown in Figures5(a) and 5(b).Figure 5
The gas pressure cloud diagram from simulation.
(a)
The first section(b)
The middle sectionAccording to Figure5, with more gas extraction, the pressure of the gas in the coal seam experienced constant changes. The gas pressure closer to the drainage hole dropped most. A further distance from the drainage hole leads to less decrease in the gas pressure until the threshold distance was reached where the gas pressure stayed at the initial value, indicating that the coal outside the impact boundary was not subject to the impacts of the extraction. In order to describe the gas pressure distribution in the five-flower arrangement, the structure line of the value line graph shown in Figure 6 was inserted into the simulated gas pressure cloud graph, which can vividly illustrate the distribution of the gas pressure value in the area wrapped by the five-flower hole. The value line graph is shown in Figure 7.Figure 6
The structure line of the value line graph.Figure 7
The value line graph.According to Figure5, the maximum gas pressure in the wrapping area of the five-flower hole is 0.22 MPa, and the gas pressure in other places is lower than the set value. At this time, the distance between two adjacent drainage holes at the same level is concluded as the effective drilling spacing.According to the simulation results, the effective drilling spacing values of the two sections in the five-flower arrangement are obtained, as shown in Table3.Table 3
Effective borehole spacing values in different sections.
Roadway segmentationExtraction timeHorizontal effective drilling spacing valueLongitudinal row spacingFirst segmentation110 d1.15 m2.20 mMiddle segmentation210 d2.45 m2.25 m


#### 4.2.2. The Analysis of the Results of the End Sections

The extraction time of the end section was set at 310 days. The simulation suggests that the coal seam gas pressure can be reduced to the standard value within a fixed time when the three-flower arrangement was adopted. The gas pressure cloud diagram in the simulation process is shown in Figure8.Figure 8
The gas pressure cloud diagram.Similarly, from Figure8, being closer to the drainage hole leads to higher pressure drops. The gas pressure reaches the peak value of 0.22 MPa at the triangle center of gravity. In order to better describe the distribution of gas pressure in the three-flower hole, the construction line of the value-taking line diagram shown in Figure 9 was inserted into the simulated gas pressure cloud diagram, which can vividly illustrate the distribution of the gas pressure value in the area wrapped by the three-flower hole, and the value-taking line diagram is shown in Figure 10.Figure 9
The structure line of the value line graph.Figure 10
The value line graph.According to Figure10, the maximum gas pressure in the wrapping area of the three-flower hole is 0.22 MPa, and the gas pressure in other places is lower than the set value. At this time, the distance between two adjacent extraction holes at the same level is the effective drilling spacing.The simulation has suggested that the horizontal effective spacing between boreholes is 2.00 m with the vertical row spacing of 2.00 m when the three-flower pattern layout is adopted in the end section.


### 4.3. A Proposal of a Mixed Layout Scheme

According to the previous calculation results, a mixed layout scheme was developed specifically for the gas drainage in the no. 1207 working face of the Sima mine. The details are presented as follows.The first section: a five-flower pattern layout was adopted with the horizontal drilling spacing of 1.15 m. The upper row of drilling is 5.50 m away from the floor longitudinally with the middle row of drilling 3.30 m away from the floor, and the bottom row of drilling1.10 m away from the floor.The middle section: a five-flower pattern layout is adopted with the horizontal drilling spacing as 2.45 m. The upper row of drilling was 5.55 m away from the floor longitudinally with the middle row of drilling 3.30 m away from the floor, and the bottom row of drilling1.05 m away from the floor.The end section: a three-flower pattern arrangement was adopted with the horizontal drilling spacing as 2.00 m. The upper row of drilling was 4.92 m away from the floor longitudinally with the bottom row of drilling 1.68 m away from the floor.


## 4.1. The Simulation Process Analysis

The strike length of the 1207 working face in the #3 Coal Seam of the Sima mine is 220 m with an inclined length of 1092 m. During the production, in order to ensure the safety of production, gas extraction was carried out. However, this kind of pumping scheme is single without considering the specific circumstances.As mentioned previously, the drainage time greatly affects the gas drainage. Different drainage time periods should adopt different drainage hole layout schemes, such as different borehole spacings or patterns. As the no. 1207 working face advanced, different extraction durations were identified in different positions of extraction boreholes in the mining roadway. Under such circumstance, different sections should be divided based on the different extraction durations, allowing individual extraction design corresponding to different field conditions.According to the comprehensive analysis of the driving speed, working face layout time, and mining speed of the mining roadway in the no. 1207 working face, and taking the drainage effectiveness into consideration, the coal seam gas drainage in the no. 1207 working face can be divided into three sections in simulation, as shown in Figure4. The average extraction time of the three sections is 110D, 210D, and 310D, respectively.The average extraction time, negative pressure, and hole diameter of the three sections were identified, and the gas extraction optimization of the three sections was also studied in detail. The calculation model of the optimization of the extraction changes from the single-row hole layout to the five-row pattern hole layout until the optimal result was obtained. Under the extraction of the experimental scheme, the goal of the whole layer of the #3 coal extraction up to 0.22 MPa was achieved.


## 4.2. The Analysis of Simulation Results

### 4.2.1. The Analysis of the First and Middle Sections

The extraction time of the first section and the middle section was set at 110 days and 210 days, respectively. According to the simulation calculation, the coal seam gas pressure can be reduced to the standard value within a fixed time only when the five-flower layout was adopted. Due to the different extraction durations, the effective drilling spaces of the two sections varied. The gas pressure cloud diagram in the simulation process is shown in Figures5(a) and 5(b).Figure 5
The gas pressure cloud diagram from simulation.
(a)
The first section(b)
The middle sectionAccording to Figure5, with more gas extraction, the pressure of the gas in the coal seam experienced constant changes. The gas pressure closer to the drainage hole dropped most. A further distance from the drainage hole leads to less decrease in the gas pressure until the threshold distance was reached where the gas pressure stayed at the initial value, indicating that the coal outside the impact boundary was not subject to the impacts of the extraction. In order to describe the gas pressure distribution in the five-flower arrangement, the structure line of the value line graph shown in Figure 6 was inserted into the simulated gas pressure cloud graph, which can vividly illustrate the distribution of the gas pressure value in the area wrapped by the five-flower hole. The value line graph is shown in Figure 7.Figure 6
The structure line of the value line graph.Figure 7
The value line graph.According to Figure5, the maximum gas pressure in the wrapping area of the five-flower hole is 0.22 MPa, and the gas pressure in other places is lower than the set value. At this time, the distance between two adjacent drainage holes at the same level is concluded as the effective drilling spacing.According to the simulation results, the effective drilling spacing values of the two sections in the five-flower arrangement are obtained, as shown in Table3.Table 3
Effective borehole spacing values in different sections.
Roadway segmentationExtraction timeHorizontal effective drilling spacing valueLongitudinal row spacingFirst segmentation110 d1.15 m2.20 mMiddle segmentation210 d2.45 m2.25 m


### 4.2.2. The Analysis of the Results of the End Sections

The extraction time of the end section was set at 310 days. The simulation suggests that the coal seam gas pressure can be reduced to the standard value within a fixed time when the three-flower arrangement was adopted. The gas pressure cloud diagram in the simulation process is shown in Figure8.Figure 8
The gas pressure cloud diagram.Similarly, from Figure8, being closer to the drainage hole leads to higher pressure drops. The gas pressure reaches the peak value of 0.22 MPa at the triangle center of gravity. In order to better describe the distribution of gas pressure in the three-flower hole, the construction line of the value-taking line diagram shown in Figure 9 was inserted into the simulated gas pressure cloud diagram, which can vividly illustrate the distribution of the gas pressure value in the area wrapped by the three-flower hole, and the value-taking line diagram is shown in Figure 10.Figure 9
The structure line of the value line graph.Figure 10
The value line graph.According to Figure10, the maximum gas pressure in the wrapping area of the three-flower hole is 0.22 MPa, and the gas pressure in other places is lower than the set value. At this time, the distance between two adjacent extraction holes at the same level is the effective drilling spacing.The simulation has suggested that the horizontal effective spacing between boreholes is 2.00 m with the vertical row spacing of 2.00 m when the three-flower pattern layout is adopted in the end section.


## 4.2.1. The Analysis of the First and Middle Sections

The extraction time of the first section and the middle section was set at 110 days and 210 days, respectively. According to the simulation calculation, the coal seam gas pressure can be reduced to the standard value within a fixed time only when the five-flower layout was adopted. Due to the different extraction durations, the effective drilling spaces of the two sections varied. The gas pressure cloud diagram in the simulation process is shown in Figures5(a) and 5(b).Figure 5
The gas pressure cloud diagram from simulation.
(a)
The first section(b)
The middle sectionAccording to Figure5, with more gas extraction, the pressure of the gas in the coal seam experienced constant changes. The gas pressure closer to the drainage hole dropped most. A further distance from the drainage hole leads to less decrease in the gas pressure until the threshold distance was reached where the gas pressure stayed at the initial value, indicating that the coal outside the impact boundary was not subject to the impacts of the extraction. In order to describe the gas pressure distribution in the five-flower arrangement, the structure line of the value line graph shown in Figure 6 was inserted into the simulated gas pressure cloud graph, which can vividly illustrate the distribution of the gas pressure value in the area wrapped by the five-flower hole. The value line graph is shown in Figure 7.Figure 6
The structure line of the value line graph.Figure 7
The value line graph.According to Figure5, the maximum gas pressure in the wrapping area of the five-flower hole is 0.22 MPa, and the gas pressure in other places is lower than the set value. At this time, the distance between two adjacent drainage holes at the same level is concluded as the effective drilling spacing.According to the simulation results, the effective drilling spacing values of the two sections in the five-flower arrangement are obtained, as shown in Table3.Table 3
Effective borehole spacing values in different sections.
Roadway segmentationExtraction timeHorizontal effective drilling spacing valueLongitudinal row spacingFirst segmentation110 d1.15 m2.20 mMiddle segmentation210 d2.45 m2.25 m


## 4.2.2. The Analysis of the Results of the End Sections

The extraction time of the end section was set at 310 days. The simulation suggests that the coal seam gas pressure can be reduced to the standard value within a fixed time when the three-flower arrangement was adopted. The gas pressure cloud diagram in the simulation process is shown in Figure8.Figure 8
The gas pressure cloud diagram.Similarly, from Figure8, being closer to the drainage hole leads to higher pressure drops. The gas pressure reaches the peak value of 0.22 MPa at the triangle center of gravity. In order to better describe the distribution of gas pressure in the three-flower hole, the construction line of the value-taking line diagram shown in Figure 9 was inserted into the simulated gas pressure cloud diagram, which can vividly illustrate the distribution of the gas pressure value in the area wrapped by the three-flower hole, and the value-taking line diagram is shown in Figure 10.Figure 9
The structure line of the value line graph.Figure 10
The value line graph.According to Figure10, the maximum gas pressure in the wrapping area of the three-flower hole is 0.22 MPa, and the gas pressure in other places is lower than the set value. At this time, the distance between two adjacent extraction holes at the same level is the effective drilling spacing.The simulation has suggested that the horizontal effective spacing between boreholes is 2.00 m with the vertical row spacing of 2.00 m when the three-flower pattern layout is adopted in the end section.


## 4.3. A Proposal of a Mixed Layout Scheme

According to the previous calculation results, a mixed layout scheme was developed specifically for the gas drainage in the no. 1207 working face of the Sima mine. The details are presented as follows.The first section: a five-flower pattern layout was adopted with the horizontal drilling spacing of 1.15 m. The upper row of drilling is 5.50 m away from the floor longitudinally with the middle row of drilling 3.30 m away from the floor, and the bottom row of drilling1.10 m away from the floor.The middle section: a five-flower pattern layout is adopted with the horizontal drilling spacing as 2.45 m. The upper row of drilling was 5.55 m away from the floor longitudinally with the middle row of drilling 3.30 m away from the floor, and the bottom row of drilling1.05 m away from the floor.The end section: a three-flower pattern arrangement was adopted with the horizontal drilling spacing as 2.00 m. The upper row of drilling was 4.92 m away from the floor longitudinally with the bottom row of drilling 1.68 m away from the floor.


## 5. The Industrial Test

### 5.1. The Detailed Test Plan

In this industrial practice, the air-return roadway of the no. 1207 working face was selected. Without affecting the original production and extraction, the extraction practice test was carried out at 450 m~485 m in the middle section and 750 m~785 m in the end section of the air-return roadway. According to the previous calculation results, the specific scheme is proposed as follows:(1)
At 450 m~485 m in the middle section of the air-return roadwayThe five-flower hole layout was adopted along with a negative pressure of 35 kPa, a diameter of 133 mm, an average length of 161 m, and a horizontal drilling spacing of 2.45 m. Longitudinally, the upper row of holes was 5.55 m away from the floor, with the middle row of holes 3.30 m away from the floor, and the bottom row of holes 1.05 m away from the floor.(2)
At 750 m~785 m in the end section of the air-return roadwayThe three-flower hole layout was adopted with a negative pressure of 35 kPa, a diameter of 133 mm, an average length of 161 m, and a horizontal drilling spacing of 2.00 m. Longitudinally, the distance between the upper row and the floor was 4.92 m, with the distance between the bottom row and the floor at 1.68 m.


### 5.2. The Analysis of the Field Test

The industrial test scheme for gas drainage revealed that the drainage boreholes cover the coal seam effectively and evenly. The residual gas content in the middle and end sections of the air roadway was reduced by 31.24% and 32.16%, respectively, with the analytic gas content reduced by 46.86% and 49.81%, respectively, and the coal seam gas pressure reduced by 45.03% and 48.33%, respectively.


## 5.1. The Detailed Test Plan

In this industrial practice, the air-return roadway of the no. 1207 working face was selected. Without affecting the original production and extraction, the extraction practice test was carried out at 450 m~485 m in the middle section and 750 m~785 m in the end section of the air-return roadway. According to the previous calculation results, the specific scheme is proposed as follows:(1)
At 450 m~485 m in the middle section of the air-return roadwayThe five-flower hole layout was adopted along with a negative pressure of 35 kPa, a diameter of 133 mm, an average length of 161 m, and a horizontal drilling spacing of 2.45 m. Longitudinally, the upper row of holes was 5.55 m away from the floor, with the middle row of holes 3.30 m away from the floor, and the bottom row of holes 1.05 m away from the floor.(2)
At 750 m~785 m in the end section of the air-return roadwayThe three-flower hole layout was adopted with a negative pressure of 35 kPa, a diameter of 133 mm, an average length of 161 m, and a horizontal drilling spacing of 2.00 m. Longitudinally, the distance between the upper row and the floor was 4.92 m, with the distance between the bottom row and the floor at 1.68 m.


## 5.2. The Analysis of the Field Test

The industrial test scheme for gas drainage revealed that the drainage boreholes cover the coal seam effectively and evenly. The residual gas content in the middle and end sections of the air roadway was reduced by 31.24% and 32.16%, respectively, with the analytic gas content reduced by 46.86% and 49.81%, respectively, and the coal seam gas pressure reduced by 45.03% and 48.33%, respectively.


## 6. Conclusion

As the extraction duration, the negative pressure of the extraction, and the extraction drill hole diameter increase, the effective extraction radius of the single drainage hole increases, demonstrating a positive correlation.Based on the simulation, a proposal including an effective spacing among extraction holes was developed. In the proposal, the extraction holes are arranged based on the effective drainage hole distances, taking various extraction durations corresponding to different locations in the roadway into consideration.According to the gas pressure cloud diagram, a greater pressure drop is identified closer to the extraction hole. The gas pressure stays stable after reaching the threshold distance.The gas pressure decreases with the increase of extraction time and finally grows stable. Some superposition effects have been identified between two drainage holes with a more intensive gas pressure drop in the superposition area.



---

*Source: 5510566-2021-09-08.xml*