# Numerical Optimization of Drilling Parameters for Gas Predrainage and Excavating-Drainage Collaboration on Roadway Head

Authors: Mingkun Luo; Lei Yang; Haiou Wen; Dongsheng Zhao; Kun Wang
Journal: Geofluids
 (2022)

---

## Abstract

In order to improve efficiency of gas drainage, the reasonable layout parameters for borehole on roadway head are investigated. Assuming that the coal mass has a dual pore structure with fractures and pores, a multiphysical field coupling mathematical model is proposed for gas drainage. The governing equations of gas adsorption, seepage, and diffusion are considered. The process of gas preextraction and excavating-extraction collaboration on roadway heading face in 2606 haulage roadway of Zhangcun mine is modeled by finite element (FE) method. The effects of gas drainage under different drilling arrangements are analyzed. The results show that the fastest decreasing rate of gas content occurs at the beginning of drainage, and decreasing rate of gas content tends to be stable in the later stage. After removing the predraining on roadway head, the gas content at the center of coal wall on the roadway heading face will rebound. The included angle between the borehole and the roadway middle line determines the spatial area of gas drainage. Too small included angle will reduce the area of gas content reduction, while too large included angle will result in a blind area on both sides of the roadway. The change in borehole spacing affects the decreasing rate of gas content in the rock and coal around the roadway. Large spacing may cause an inadequate decreasing rate in the middle of the roadway, while small spacing may cause an inadequate decreasing rate in both sides of the roadway. In the field application, the pure gas flow and gas concentration of drainage show a downward trend in whole. The gas concentration decreased from 13.2% to 10.8%, and the pure gas flow decreased from 10.6 m3/min to 8.15 m3/min. Research achievements can provide a basis for gas drainage in underground roadway.


---

## Body

## 1. Introduction

Coal mine gas is the main harm of coal production, it is also a kind of clean energy [1]. As the depth of mining increases, the permeability of coal seams decreases gradually [2, 3]. In this case, the danger and intensity of coal and gas outburst are increasing [4]. Underground gas drainage is the main method to reduce gas content and prevent gas disasters [5, 6]. Therefore, studying the seepage law and reasonable drilling layout parameters during gas drainage is of vital meaning for effective gas predrainage and reducing engineering costs.Predecessors have done a lot of research on the seepage law during gas drainage. Wu [7] proposed a gas-coal coupling model based on the interaction among gas diffusion, seepage, and coal deformation. Yin et al. [8, 9] considered the Klinkenberg effect and ab/desorption induced strain effect to establish a gas-solid mathematic model and perform numerical simulations with the aid of finite elements. Zhang et al. [10] obtained the porosity and permeability equations with adsorption caused deformation. According to the gas seepage characteristics and dynamic change in coal permeability, Liang et al. [11] studied the gas seepage processes in fracture, microporous adsorption gas desorption and diffusion, and deformation of coal and rock mass and established a couple seepage model. Lin et al. [12] established a multiphysics coupled model of diffusion, stress, and seepage fields and introduced the dynamic diffusion coefficient to explore the seepage evolution during gas drainage. In addition, some scholars have studied the desorption and adsorption behavior of coal at the microscopic level [13, 14] and the influence of polymer on coal desorption and adsorption [15]. For the borehole spacing, too small spacing will increase the amount of work in the drilling construction and rise the instability risk of drilling such as pierce or collapse holes resulting from the pressure relief damage of coal and rock, while too large spacing will leave blind area of drainage, which not only can hardly eliminate the outbursts but also aggravate the hidden dangers of gas accidents [16]. To this end, scholars have launched a study on the extraction radius. Zhang et al. [17] used a self-developed device to conduct experimental research on different drilling arrangements and explored the extraction radius when the drilling holes interact with each other. Lin et al. [18] studied the effect of the adsorption constant on the extraction radius through numerical simulations. Qin et al. [19] obtained reasonable drainage parameters through the study of drainage radius. These provide a reference for numerical simulation and engineering practice of underground gas drainage.However, there are few studies for gas predrainage and excavating-drainage collaboration on roadway head. This paper is based on the assumption of dual porosity and single permeability and comprehensively considers the gas slippage effect, effective stress, and deformation characteristics of coal seam, and a multiphysics coupled model composed of fluid migration, coal deformation under stresses, and gas diffusion controlling equations is established. This mathematical model was embedded in the FE software-COMSOL Multiphysics, to simulate gas drainage law of the 2606 haulage roadway in Zhangcun Mine. The processes of gas predrainage and excavating-drainage collaboration on roadway head under different schemes were simulated. The research results can provide guidance for optimizing of drainage boreholes layout.


## 2. Coal Seam Gas Drainage Mathematical Model

### 2.1. General Hypothesis

Based on the general law of gas transport within coal seam and related literatures [20–23], some hypotheses are proposed: (1) coal and rock mass has a dual pore structure composed with pores-fractures. (2) Gas diffusion in the pore system of coal matrix follows Fick’s law, and gas seepage in fracture system follows Darcy’s law. (3) Water only moves in the fractures. (4) Ignore the influence of coal seam temperature on gas drainage. (5) Gas is regarded as an ideal gas. (6) The tensile stress is positive.


### 2.2. Permeability Evolution Model

According to hypothesis, the matrix porosity model can be obtained [24]:
(1)φm=1+S0φm0+αmS−S01+S,where φm0 is the original matrix porosity; αm=1−K/Ks is the Biot’s coefficient for matrix; K is the bulk modulus of coal, GPa; Ks=Es/31−2v is the bulk modulus of coal skeleton, GPa; Es is the elastic modulus of coal skeleton, GPa; v is the Poisson’s ratio of coal; S=εv+pmg/Ks−εa is the matrix porosity strain; εv is the volumetric strain; the term “0” stands the original value.The relationship between adsorbed gas strain of coal skeleton and the adsorbed amount can be expressed as follows [25]:
(2)εa=αsgVsg,where αsg is the sorption induced strain factor, kg/m3; Vsg is the adsorbed gas content, m3/kg.Based on the modified Langmuir equation, the adsorbed gas content may be expressed by the following [26]:
(3)Vsg=VLpmgPL+pmg,where VL and PL are the Langmuir constants; pmg is gas pressure in coal matrix, MPa.The fracture porosity model is as follows [27]:
(4)φf=φf0−3φf0φf0+3Kf/Kεa−εa0−εv−εvo,where φf0 is the original fracture porosity; Kf=qKn is the effective stiffness, GPa; q is the fracture width, m; Kn is fracture stiffness, GPa/m.The permeability of fracture can be obtained by cubic law [28]:
(5)k=k0⋅φfφf03=k01−3φf0+3Kf/Kεa−εa0−εv−εvo3,where k0 is the origin permeability, m2.The relative permeabilities are set to illustrate the gas and water transport speed in two-phase flow, which are derived as follows [29]:
(6)krg=krg01−sw−swr1−swr−sgr21−sw−swr1−swr2,krw=krw0sw−swr1−swr4,where krg0 is gas endpoint permeability; sw is water saturation; swr is bound water saturation; sgr is residual gas saturation; krw0 is water endpoint permeability.


### 2.3. Seepage Field Controlling Equation

#### 2.3.1. Gas Transport in Matrix

Gas mass in one unit coal matrix contains both free and adsorbed components [30]:
(7)mm=φmρg+Vsgρcρgs,where ρs is the density of coal skeleton, kg/m3; ρgs is gas density under standard conditions, kg/m3.Satisfying the ideal gas law, the gas density can be obtained as follows:(8)ρg=MgRTp,where Mg is the gas molar mass, g/mol; R is the gas mole constant, J/(mol·K); p is the pressure, MPa; T is coal temperature, K.Due to the effect of drainage breaking the dynamic balance in gas ad/desorption in coal seams, the free gas in matrix pores diffuses into coal fracture caused by the concentration difference. The gas balance in coal matrix can be derived from Fick’s diffusion law. For [28],
(9)∂mm∂t=−MgτRTpmg−pfg,where τ is the gas sorption time, s; pfg is the fracture pressure, MPa.Combining Equations (3) and (7)–(9), the controlling equation for gas migration in coal matrix is gained.
(10)∂∂tφmMgRTpmg+VLpmgpL+pmgρcMgRTpa=−MgτRTpmg−pfg.


#### 2.3.2. Fluid Migration within Fractures

The modified Darcy’s law is adopted to govern the fluid flow of gas and water in fractures [31]:
(11)qg=−kkrgμg1+bpfg∇pfg,qw=−kkrwμw∇pfw,where b is the Klinkenberg factor, MPa; krw is water relative permeability; krg is gas relative permeability; μw is the dynamic viscosity of water, Pa·s; μg is the dynamic viscosity of gas, Pa·s.The gas-water balance equations in the fracture are expressed as follows [32]:
(12)∂∂tsgφfMgRTpfg−∇⋅MgRTkkrgpfg+bμg∇pfg=1−φfMgτRTpmg−pfg,∂swφfρw∂t‐∇⋅ρwkkrwμw∇pfw=0,where sg means gas saturation and ρw means water density, kg/m3.


### 2.4. Stress Field Controlling Equation

The deformation of coal seam concludes solid stress caused strain, gas and water pressure caused strain, and gas sorption induced strain. The controlling equation for coal stress-strain field is as follows [33]:
(13)Gui,jj+G1−2vuj,ji−αmpm,i−αfpf,i−Kεa,i+Fi=0,where G is the shear modulus of coal, GPa; pf is the fluid pressure in fractures, Pa; Fi is the volume force, GPa; αf is the Biot’s coefficient for fracture.


## 2.1. General Hypothesis

Based on the general law of gas transport within coal seam and related literatures [20–23], some hypotheses are proposed: (1) coal and rock mass has a dual pore structure composed with pores-fractures. (2) Gas diffusion in the pore system of coal matrix follows Fick’s law, and gas seepage in fracture system follows Darcy’s law. (3) Water only moves in the fractures. (4) Ignore the influence of coal seam temperature on gas drainage. (5) Gas is regarded as an ideal gas. (6) The tensile stress is positive.


## 2.2. Permeability Evolution Model

According to hypothesis, the matrix porosity model can be obtained [24]:
(1)φm=1+S0φm0+αmS−S01+S,where φm0 is the original matrix porosity; αm=1−K/Ks is the Biot’s coefficient for matrix; K is the bulk modulus of coal, GPa; Ks=Es/31−2v is the bulk modulus of coal skeleton, GPa; Es is the elastic modulus of coal skeleton, GPa; v is the Poisson’s ratio of coal; S=εv+pmg/Ks−εa is the matrix porosity strain; εv is the volumetric strain; the term “0” stands the original value.The relationship between adsorbed gas strain of coal skeleton and the adsorbed amount can be expressed as follows [25]:
(2)εa=αsgVsg,where αsg is the sorption induced strain factor, kg/m3; Vsg is the adsorbed gas content, m3/kg.Based on the modified Langmuir equation, the adsorbed gas content may be expressed by the following [26]:
(3)Vsg=VLpmgPL+pmg,where VL and PL are the Langmuir constants; pmg is gas pressure in coal matrix, MPa.The fracture porosity model is as follows [27]:
(4)φf=φf0−3φf0φf0+3Kf/Kεa−εa0−εv−εvo,where φf0 is the original fracture porosity; Kf=qKn is the effective stiffness, GPa; q is the fracture width, m; Kn is fracture stiffness, GPa/m.The permeability of fracture can be obtained by cubic law [28]:
(5)k=k0⋅φfφf03=k01−3φf0+3Kf/Kεa−εa0−εv−εvo3,where k0 is the origin permeability, m2.The relative permeabilities are set to illustrate the gas and water transport speed in two-phase flow, which are derived as follows [29]:
(6)krg=krg01−sw−swr1−swr−sgr21−sw−swr1−swr2,krw=krw0sw−swr1−swr4,where krg0 is gas endpoint permeability; sw is water saturation; swr is bound water saturation; sgr is residual gas saturation; krw0 is water endpoint permeability.


## 2.3. Seepage Field Controlling Equation

### 2.3.1. Gas Transport in Matrix

Gas mass in one unit coal matrix contains both free and adsorbed components [30]:
(7)mm=φmρg+Vsgρcρgs,where ρs is the density of coal skeleton, kg/m3; ρgs is gas density under standard conditions, kg/m3.Satisfying the ideal gas law, the gas density can be obtained as follows:(8)ρg=MgRTp,where Mg is the gas molar mass, g/mol; R is the gas mole constant, J/(mol·K); p is the pressure, MPa; T is coal temperature, K.Due to the effect of drainage breaking the dynamic balance in gas ad/desorption in coal seams, the free gas in matrix pores diffuses into coal fracture caused by the concentration difference. The gas balance in coal matrix can be derived from Fick’s diffusion law. For [28],
(9)∂mm∂t=−MgτRTpmg−pfg,where τ is the gas sorption time, s; pfg is the fracture pressure, MPa.Combining Equations (3) and (7)–(9), the controlling equation for gas migration in coal matrix is gained.
(10)∂∂tφmMgRTpmg+VLpmgpL+pmgρcMgRTpa=−MgτRTpmg−pfg.


### 2.3.2. Fluid Migration within Fractures

The modified Darcy’s law is adopted to govern the fluid flow of gas and water in fractures [31]:
(11)qg=−kkrgμg1+bpfg∇pfg,qw=−kkrwμw∇pfw,where b is the Klinkenberg factor, MPa; krw is water relative permeability; krg is gas relative permeability; μw is the dynamic viscosity of water, Pa·s; μg is the dynamic viscosity of gas, Pa·s.The gas-water balance equations in the fracture are expressed as follows [32]:
(12)∂∂tsgφfMgRTpfg−∇⋅MgRTkkrgpfg+bμg∇pfg=1−φfMgτRTpmg−pfg,∂swφfρw∂t‐∇⋅ρwkkrwμw∇pfw=0,where sg means gas saturation and ρw means water density, kg/m3.


## 2.3.1. Gas Transport in Matrix

Gas mass in one unit coal matrix contains both free and adsorbed components [30]:
(7)mm=φmρg+Vsgρcρgs,where ρs is the density of coal skeleton, kg/m3; ρgs is gas density under standard conditions, kg/m3.Satisfying the ideal gas law, the gas density can be obtained as follows:(8)ρg=MgRTp,where Mg is the gas molar mass, g/mol; R is the gas mole constant, J/(mol·K); p is the pressure, MPa; T is coal temperature, K.Due to the effect of drainage breaking the dynamic balance in gas ad/desorption in coal seams, the free gas in matrix pores diffuses into coal fracture caused by the concentration difference. The gas balance in coal matrix can be derived from Fick’s diffusion law. For [28],
(9)∂mm∂t=−MgτRTpmg−pfg,where τ is the gas sorption time, s; pfg is the fracture pressure, MPa.Combining Equations (3) and (7)–(9), the controlling equation for gas migration in coal matrix is gained.
(10)∂∂tφmMgRTpmg+VLpmgpL+pmgρcMgRTpa=−MgτRTpmg−pfg.


## 2.3.2. Fluid Migration within Fractures

The modified Darcy’s law is adopted to govern the fluid flow of gas and water in fractures [31]:
(11)qg=−kkrgμg1+bpfg∇pfg,qw=−kkrwμw∇pfw,where b is the Klinkenberg factor, MPa; krw is water relative permeability; krg is gas relative permeability; μw is the dynamic viscosity of water, Pa·s; μg is the dynamic viscosity of gas, Pa·s.The gas-water balance equations in the fracture are expressed as follows [32]:
(12)∂∂tsgφfMgRTpfg−∇⋅MgRTkkrgpfg+bμg∇pfg=1−φfMgτRTpmg−pfg,∂swφfρw∂t‐∇⋅ρwkkrwμw∇pfw=0,where sg means gas saturation and ρw means water density, kg/m3.


## 2.4. Stress Field Controlling Equation

The deformation of coal seam concludes solid stress caused strain, gas and water pressure caused strain, and gas sorption induced strain. The controlling equation for coal stress-strain field is as follows [33]:
(13)Gui,jj+G1−2vuj,ji−αmpm,i−αfpf,i−Kεa,i+Fi=0,where G is the shear modulus of coal, GPa; pf is the fluid pressure in fractures, Pa; Fi is the volume force, GPa; αf is the Biot’s coefficient for fracture.


## 3. Physical Model and Solution Settings

### 3.1. Studied Case

Zhangcun Coal Mine is situated at the eastern edge of Qinshui Basin in Shanxi Province. The regional stratigraphic division belongs to the stratigraphic community in the southern section of Taihang Mountain in the Shanxi stratigraphic division. The stratum generally moves in the NNE direction, dipping gently to W, with an inclination angle of 5-15°. Lane 2606 is located on the north side of lane 26, 300 m away from lane 2605 in the east. The ground elevation of the excavation area is about+972~+992 m. The cross-section of the roadway is 5.8×4.0 m (width×height) rectangle with a cross-sectional area of 23.2 m2. The mined seam belongs to No. 3 coal seam in the Shanxi Formation. The coal seam is stable, with less variation in coal thickness and simpler coal seam structure. The buried depth is about 537 m, and the thickness is 5.33 m-6.19 m, and the average thickness is 5.86 m. The original gas content during the 2606 roadway driving period is 8.5-10.0 m3/t, and the coal seam temperature is 15°C-17°C. During the excavating process, the mining depth continued to increase, and the gas content for 3# coal seam showed an increasing trend.Therefore, gas preextraction and excavating extraction collaboration were adopted to ensure the safe operation of the 2606 roadway during excavation. During the excavating period, the gas predraining is carried out every 130 m. The gas predraining is arranged head-on with 5 drilling holes on the working face. Drilling length is 150 m, and the safety distance is 20 m. The method of excavating extraction is to arrange the stepping drill field in the roadway at a step distance of 50 m (the same side drill field spacing is 100 m). The stepping drilling field is arranged, and each stepping drilling field is arranged with 6 holes, and the drilling depth is 130 m. The layout of the 2606 haulage roadway of gas preextraction and excavating extraction is shown in Figure1.Figure 1
Schematic diagram of the layout of the 2606 haulage roadway of gas preextraction and excavating extraction collaboration.
(a)
Top view(b)
Front view


### 3.2. Physical Model and Definite Solution Conditions

According to 2606 roadway gas extraction scheme, the physical geometric model as shown in Figure2 is established. And through the simulation results obtained after changing spacing and angle of the boreholes, determine the best scheme for gas drainage on 2606 haulage roadway of Zhangcun Coal Mine. The reference line A-B in the physical model is set to see the changes in gas content. The initial gas pressure is 0.8 MPa, the overburden gravity of 14.85 MPa in the z direction is applied of the model, and roller bearing boundary conditions are set around the model. The bottom boundary is fixed. The outer boundaries of the model are impermeable. The diameter of the borehole is 94 mm, and the negative pressure for gas drainage is 20 kPa. The other used parameters are shown in Table 1, which are derived from literature [33] and field data.Figure 2
Numerical simulation of physical geometry model.Table 1
Numerical simulation parameters.
ParameterValueParameterValueInitial CH4 pressure (p0, MPa)0.80Gas sorption time CH4 (τ, d)4.34Young’s modulus of coal seam (E, MPa)3500Coal density (ρc, kg·m-3)1380Young’s modulus of skeleton (Es, MPa)8469Klinkenberg factor (b, MPa)0.62Poisson’s ratio of coal (v)0.30Original temperature (T, K)289.15Langmuir volume constant of CH4 (VL, m3·kg-1)0.0323Matrix initial porosity (φm0)0.04Langmuir pressure constant of CH4 (PL, MPa)2.0833Initial porosity of fracture (φf0)0.018Initial permeability (k0, m2)2.56×10−17Residual gas saturation (sgr)0.15Dynamic viscosity of water (μw, Pa·s)1.01×10−3Initial water saturation (sw0)0.6Dynamic viscosity of CH4 (μg1, Pa·s)1.03×10−5Irreducible water saturation (swr)0.42


## 3.1. Studied Case

Zhangcun Coal Mine is situated at the eastern edge of Qinshui Basin in Shanxi Province. The regional stratigraphic division belongs to the stratigraphic community in the southern section of Taihang Mountain in the Shanxi stratigraphic division. The stratum generally moves in the NNE direction, dipping gently to W, with an inclination angle of 5-15°. Lane 2606 is located on the north side of lane 26, 300 m away from lane 2605 in the east. The ground elevation of the excavation area is about+972~+992 m. The cross-section of the roadway is 5.8×4.0 m (width×height) rectangle with a cross-sectional area of 23.2 m2. The mined seam belongs to No. 3 coal seam in the Shanxi Formation. The coal seam is stable, with less variation in coal thickness and simpler coal seam structure. The buried depth is about 537 m, and the thickness is 5.33 m-6.19 m, and the average thickness is 5.86 m. The original gas content during the 2606 roadway driving period is 8.5-10.0 m3/t, and the coal seam temperature is 15°C-17°C. During the excavating process, the mining depth continued to increase, and the gas content for 3# coal seam showed an increasing trend.Therefore, gas preextraction and excavating extraction collaboration were adopted to ensure the safe operation of the 2606 roadway during excavation. During the excavating period, the gas predraining is carried out every 130 m. The gas predraining is arranged head-on with 5 drilling holes on the working face. Drilling length is 150 m, and the safety distance is 20 m. The method of excavating extraction is to arrange the stepping drill field in the roadway at a step distance of 50 m (the same side drill field spacing is 100 m). The stepping drilling field is arranged, and each stepping drilling field is arranged with 6 holes, and the drilling depth is 130 m. The layout of the 2606 haulage roadway of gas preextraction and excavating extraction is shown in Figure1.Figure 1
Schematic diagram of the layout of the 2606 haulage roadway of gas preextraction and excavating extraction collaboration.
(a)
Top view(b)
Front view


## 3.2. Physical Model and Definite Solution Conditions

According to 2606 roadway gas extraction scheme, the physical geometric model as shown in Figure2 is established. And through the simulation results obtained after changing spacing and angle of the boreholes, determine the best scheme for gas drainage on 2606 haulage roadway of Zhangcun Coal Mine. The reference line A-B in the physical model is set to see the changes in gas content. The initial gas pressure is 0.8 MPa, the overburden gravity of 14.85 MPa in the z direction is applied of the model, and roller bearing boundary conditions are set around the model. The bottom boundary is fixed. The outer boundaries of the model are impermeable. The diameter of the borehole is 94 mm, and the negative pressure for gas drainage is 20 kPa. The other used parameters are shown in Table 1, which are derived from literature [33] and field data.Figure 2
Numerical simulation of physical geometry model.Table 1
Numerical simulation parameters.
ParameterValueParameterValueInitial CH4 pressure (p0, MPa)0.80Gas sorption time CH4 (τ, d)4.34Young’s modulus of coal seam (E, MPa)3500Coal density (ρc, kg·m-3)1380Young’s modulus of skeleton (Es, MPa)8469Klinkenberg factor (b, MPa)0.62Poisson’s ratio of coal (v)0.30Original temperature (T, K)289.15Langmuir volume constant of CH4 (VL, m3·kg-1)0.0323Matrix initial porosity (φm0)0.04Langmuir pressure constant of CH4 (PL, MPa)2.0833Initial porosity of fracture (φf0)0.018Initial permeability (k0, m2)2.56×10−17Residual gas saturation (sgr)0.15Dynamic viscosity of water (μw, Pa·s)1.01×10−3Initial water saturation (sw0)0.6Dynamic viscosity of CH4 (μg1, Pa·s)1.03×10−5Irreducible water saturation (swr)0.42


## 4. Analysis of Numerical Simulation

### 4.1. Different Scheme Design

To determine a reasonable drainage drilling arrangement scheme and obtain the best drainage effect. By changing the borehole spacing and the angle, the following five schemes were formulated up for the 2606 haulage roadway gas drainage in Zhangcun Coal Mine, and using the FE software COMSOL Multiphysics, the gas drainage effect of these schemes is analyzed. The drainage schemes are shown in Table2.Table 2
Drilling layout parameters of different schemes.
Borehole orderPredrainageExcavating drainage1#2#3#4#5#1#2#3#4#5#6#Scheme oneAngle between borehole and middle line of the roadway/°00000012467Drilling spacing/m1.80.5Scheme twoAngle between borehole and middle line of the roadway/°00000000000Drilling spacing/m1.80.5Scheme threeAngle between borehole and middle line of the roadway/°000001357911Drilling spacing/m1.80.5Scheme fourAngle between borehole and middle line of the roadway/°00000012467Drilling spacing/m1.30.1Scheme fiveAngle between borehole and middle line of the roadway/°00000012467Drilling spacing/m2.30.9


### 4.2. Evolution of Gas Content under Different Schemes

In order to ensure the extraction progress of 2606 haulage roadway, the head predrainage after 10 days of drainage will be canceled. In Figure3, gas content on the reference line A-B under different borehole arrangement angles shows decrease trend as extraction time prolongs. The fastest decrease rate of gas content occurs in the initial stage of drainage. As gas pressure in coal seam decreases, the gas flow decreases, so the decrease rate of gas content tends to be stable in the later stage. Comparing Figures 3(a) and 3(b), the decrease rate of gas content in schemes 1 and 2 is basically the same near on the roadway middle line. However, the area of gas content decrease is different. Scheme 1 has a greater impact within 22 m around the roadway middle line, and scheme 2 has a greater impact within 9 m around the roadway middle line. The borehole is drilled at an angle to the roadway middle line in scheme 1, which covers wider area on both sides of the roadway. And the included angle between the borehole and the roadway middle line is 0° in scheme 2. Therefore, the spatial area of influence of the first scheme is larger than that of the second scheme. From Figure 3(c), compare scheme 3 and scheme 1, although the angle between the borehole and the roadway middle line becomes larger, and gas content decrease area is increases, but the overall decrease rate of gas content in scheme 3 is slowed down. The most important thing is that with the extension of the boreholes, the angle between the boreholes in scheme 3 gradually increases, and drainage blind areas appear, causes the gas content to decrease slowly (as can be seen in the 7 m area on the left side of the middle line of the roadway in Figure 3(c)), and the ideal drainage effect cannot be achieved. After removing the predrainage on roadway head, the gas content at the middle line of coal wall on the roadway heading face will rebound. The gas content in the roadway middle line is relatively lower than the gas content of the roadway two flanks after 10 d for gas predrainage and excavating-drainage collaboration on roadway head. When there is a pressure gradient, the gas in the roadway two sides will move to the direction of the roadway middle line. At this time, the gas content in the roadway middle line will rise.Figure 3
Evolution curve of gas content on reference line AB under different arrangements of included angle.
(a)
Scheme two(b)
Scheme one(c)
Scheme threeAs can be seen from Figure4. The change of the borehole spacing has a weak influence on the drainage area. However, this affects the gas content decreasing rate in the rock and coal around the roadway. Comparing Figures 4(a) and 4(b), in scheme 4, the peak gas content within 5 m on the roadway middle line after 10, 30, 60, 90, and 120 days of drainage is 8.3 m3/t, 7.6 m3/t, 7.3 m3/t, 6.9 m3/t, and 6.4 m3/t, respectively. The peak gas content within 5 m on the roadway middle line in scheme 1 is 7.8 m3/t, 7.0 m3/t, 6.8 m3/t, 6.4 m3/t, and 5.9 m3/t, respectively. Compared with scheme 4, the gas content of roadway two sides decreases faster in scheme 1, while the gas content of the roadway both sides decreases slowly, and as a result, the gas content of the scheme 4 rises more than the scheme 1 of the coal wall on the roadway heading face. From Figures 4(b) and 4(c), in the case of 10 days of gas predrainage and excavating drainage collaboration on roadway head, the gas content peak on the roadway midline line in scheme 5 is 6.6 m3/t, and the gas content peak on the roadway midline line the in scheme 1 is 6.1 m3/t. The gas content of the scheme 1 drops faster than that of the scheme 5 on the roadway middle line in the early stage.Figure 4
Evolution curve of gas content on reference line A-B under different borehole spacing.
(a)
Scheme four(b)
Scheme one(c)
Scheme fiveFrom the above analysis, the too small included angle between the borehole and the roadway middle line results the spatial area of gas content decrease, it not only can hardly achieve the desired drainage effect but also may increase engineering costs. When the too large angle between the borehole and the roadway middle line will result a blind area for drainage on the roadway of two sides, spacing of borehole is too small, and the results lower the decrease rate of gas content on the roadway two sides. Spacing of borehole is too large, which caused an inadequate decreasing rate of gas content on the roadway middle line, and the drainage effects of both are poor. Therefore, in general, among the five gas drainage schemes designed for the 2606 haulage roadway of Zhangcun Coal Mine, scheme 1 is the best drainage effect.


## 4.1. Different Scheme Design

To determine a reasonable drainage drilling arrangement scheme and obtain the best drainage effect. By changing the borehole spacing and the angle, the following five schemes were formulated up for the 2606 haulage roadway gas drainage in Zhangcun Coal Mine, and using the FE software COMSOL Multiphysics, the gas drainage effect of these schemes is analyzed. The drainage schemes are shown in Table2.Table 2
Drilling layout parameters of different schemes.
Borehole orderPredrainageExcavating drainage1#2#3#4#5#1#2#3#4#5#6#Scheme oneAngle between borehole and middle line of the roadway/°00000012467Drilling spacing/m1.80.5Scheme twoAngle between borehole and middle line of the roadway/°00000000000Drilling spacing/m1.80.5Scheme threeAngle between borehole and middle line of the roadway/°000001357911Drilling spacing/m1.80.5Scheme fourAngle between borehole and middle line of the roadway/°00000012467Drilling spacing/m1.30.1Scheme fiveAngle between borehole and middle line of the roadway/°00000012467Drilling spacing/m2.30.9


## 4.2. Evolution of Gas Content under Different Schemes

In order to ensure the extraction progress of 2606 haulage roadway, the head predrainage after 10 days of drainage will be canceled. In Figure3, gas content on the reference line A-B under different borehole arrangement angles shows decrease trend as extraction time prolongs. The fastest decrease rate of gas content occurs in the initial stage of drainage. As gas pressure in coal seam decreases, the gas flow decreases, so the decrease rate of gas content tends to be stable in the later stage. Comparing Figures 3(a) and 3(b), the decrease rate of gas content in schemes 1 and 2 is basically the same near on the roadway middle line. However, the area of gas content decrease is different. Scheme 1 has a greater impact within 22 m around the roadway middle line, and scheme 2 has a greater impact within 9 m around the roadway middle line. The borehole is drilled at an angle to the roadway middle line in scheme 1, which covers wider area on both sides of the roadway. And the included angle between the borehole and the roadway middle line is 0° in scheme 2. Therefore, the spatial area of influence of the first scheme is larger than that of the second scheme. From Figure 3(c), compare scheme 3 and scheme 1, although the angle between the borehole and the roadway middle line becomes larger, and gas content decrease area is increases, but the overall decrease rate of gas content in scheme 3 is slowed down. The most important thing is that with the extension of the boreholes, the angle between the boreholes in scheme 3 gradually increases, and drainage blind areas appear, causes the gas content to decrease slowly (as can be seen in the 7 m area on the left side of the middle line of the roadway in Figure 3(c)), and the ideal drainage effect cannot be achieved. After removing the predrainage on roadway head, the gas content at the middle line of coal wall on the roadway heading face will rebound. The gas content in the roadway middle line is relatively lower than the gas content of the roadway two flanks after 10 d for gas predrainage and excavating-drainage collaboration on roadway head. When there is a pressure gradient, the gas in the roadway two sides will move to the direction of the roadway middle line. At this time, the gas content in the roadway middle line will rise.Figure 3
Evolution curve of gas content on reference line AB under different arrangements of included angle.
(a)
Scheme two(b)
Scheme one(c)
Scheme threeAs can be seen from Figure4. The change of the borehole spacing has a weak influence on the drainage area. However, this affects the gas content decreasing rate in the rock and coal around the roadway. Comparing Figures 4(a) and 4(b), in scheme 4, the peak gas content within 5 m on the roadway middle line after 10, 30, 60, 90, and 120 days of drainage is 8.3 m3/t, 7.6 m3/t, 7.3 m3/t, 6.9 m3/t, and 6.4 m3/t, respectively. The peak gas content within 5 m on the roadway middle line in scheme 1 is 7.8 m3/t, 7.0 m3/t, 6.8 m3/t, 6.4 m3/t, and 5.9 m3/t, respectively. Compared with scheme 4, the gas content of roadway two sides decreases faster in scheme 1, while the gas content of the roadway both sides decreases slowly, and as a result, the gas content of the scheme 4 rises more than the scheme 1 of the coal wall on the roadway heading face. From Figures 4(b) and 4(c), in the case of 10 days of gas predrainage and excavating drainage collaboration on roadway head, the gas content peak on the roadway midline line in scheme 5 is 6.6 m3/t, and the gas content peak on the roadway midline line the in scheme 1 is 6.1 m3/t. The gas content of the scheme 1 drops faster than that of the scheme 5 on the roadway middle line in the early stage.Figure 4
Evolution curve of gas content on reference line A-B under different borehole spacing.
(a)
Scheme four(b)
Scheme one(c)
Scheme fiveFrom the above analysis, the too small included angle between the borehole and the roadway middle line results the spatial area of gas content decrease, it not only can hardly achieve the desired drainage effect but also may increase engineering costs. When the too large angle between the borehole and the roadway middle line will result a blind area for drainage on the roadway of two sides, spacing of borehole is too small, and the results lower the decrease rate of gas content on the roadway two sides. Spacing of borehole is too large, which caused an inadequate decreasing rate of gas content on the roadway middle line, and the drainage effects of both are poor. Therefore, in general, among the five gas drainage schemes designed for the 2606 haulage roadway of Zhangcun Coal Mine, scheme 1 is the best drainage effect.


## 5. Field Applications

According to the simulation of gas combined drainage on 2606 haulage roadways in Zhangcun Coal Mine. It is preliminarily determined that the spacing of the preextraction boreholes at the front end of the 2606 haulage roadway is 1.8 m. The distance between the excavating extraction boreholes is 0.5 m, and the angles between borehole and middle line of the roadway are 0°, 1°, 2°, 4°, 6°, and 7°. The boreholes are connected to the pipe network for drainage. The gas drainage concentration and drainage volume monitored during the period are shown in Figure5.Figure 5
2606 roadway gas concentration and pure gas flow change.After nearly 2 months of drainage, the gas concentration and drainage of pure quantity around the 2606 haulage roadway have significantly decreased trend, and the drainage effect is ideal. The gas concentration decreased from 13.2% to 10.8%, a decrease of 18.2%. And the pure gas flow decreased from 10.6 m3/min to 8.15 m3/min, a decrease of 23.1%. Results show that it is feasible to adopt scheme 1 as the gas drainage plan during the 2606 roadway excavation period.


## 6. Conclusions

(1)
Considering that the coal mass is a dual porous structure, a multiphysics coupling mathematical model of coal seam gas drainage was established involving gas sorption, diffusion, and seepage. The COMSOL Multiphysics software is used to simulate gas predrainage and excavating-drainage collaboration on roadway head(2)
The process of gas drainage with different borehole arrangements in 2606 roadway in Zhangcun Coal Mine is simulated. The results show that the fastest decreasing rate of gas content occurs in the initial stage of drainage, and the decreasing rate of gas content tends to be stable in the later stage. After removing the predraining on roadway head, the gas content at the center of coal wall on the roadway heading face will rebound(3)
The included angle between the borehole and roadway middle line determines the spatial area of gas drainage. Too small included angle will reduce the area of gas content reduction, while too large included angle will result in a blind area on both sides of the roadway. The change in borehole spacing affects the decreasing rate of gas content in the coal wall around the roadway. Large spacing may cause an inadequate decreasing rate in the middle coal wall of the roadway, while small spacing may cause an inadequate decreasing rate in both sides of the roadway(4)
Gas drainage is carried out in 2606 haulage roadway of Zhangcun Coal Mine. After nearly 2 months of drainage, gas concentration decreased from 13.2% to 10.8%, with a decrease of 18.2%. And the pure gas flow decreased from 10.6 m3/min to 8.15 m3/min, with a decrease of 23.1%. The pure gas flow and gas concentration of drainage show a downward trend in whole which matches the results of numerical simulations



---

*Source: 3241211-2022-05-20.xml*