# Cocompetitive Characteristics and Quantitative Design of Engineering Parameters for Coal Gas Predrainage Boreholes

Authors: Tong-qiang Xia; Jiao-fei He; Zi-long Li; Xiao-yan Sun; Dun-shuai Sun; Jian-hang Lu; Hao-jie Cui
Journal: Geofluids
 (2022)

---

## Abstract

Quantitative design of the optimal gas extraction engineering parameters has always been a key scientific problem to meet the predrainage target of the coal seam, which is paid great attention to and urgently expected to solve for the gas extraction engineering of coal mine. In this paper, the governing equation of the evolution of coal seam gas content is established under the interference and cooperation of gas extraction by boreholes, and the calculation formula of gas cooperative efficiency is defined. Based on gas-air mixed flow, a multifield coupling 3d model of coal gas extraction using boreholes was developed. The content, method, and principle of optimal spacing parameter design of extraction borehole were proposed, and the constraint conditions and objective function of optimal design were constructed. The flow characteristics of gas extraction using boreholes in the coal seam area and the influence of different parameters on gas extraction effect are studied. The cooperative and competitive behaviors of gas extraction by boreholes in coal seam are revealed. The results show that borehole spacing, seam fracture permeability, and gas adsorption constant have significant effects on the synergistic and competitive effects of gas extraction, while the negative pressure of gas extraction and coal matrix permeability have little effects, but the negative pressure of gas extraction has great effects on gas extraction concentration. Finally, Comsol With Matlab is used to solve the multiconstraint nonlinear optimization problem of gas drainage engineering parameters along coal seam in 18401 working face of Tunlan Coal Mine. The effect of actual gas extraction investigation is very remarkable, which provides a quantitative calculation basis for scientific gas extraction in coal mine.


---

## Body

## 1. Introduction

Coal is still one of the most important energy sources in the world for a long time, and safety and efficiency in coal mining are particularly important [1, 2]. Gas accident is the “first killer” threatening the safety production and the biggest hidden danger of coal mine to achieve the goal of “zero gas overrun, zero coal seam outburst, and zero enterprise death” [3, 4]. Order [2015] No. 82 of the State Administration of Work Safety issued 10 regulations on strengthening coal mine gas control, such as “pumping before excavation, pumping before mining, and pumping up to standard.” At present, in the technology of gas accident prevention and gas control, in addition to ensuring good ventilation network and environment, gas drainage using coal seam borehole is one of the most commonly used and effective methods [5, 6]. A coal seam borehole is an artificial channel established between coal seam gas and drainage pipeline system. The scientificity and rationality of borehole design directly affect the difficulty and quantity of drilling construction, as well as the exertion of later drainage efficiency [7]. However, the design of coal seam gas drainage mainly depends on experience, resulting in a tight gas predrainage period and conflict with mining operation, so that the unsafe mining phenomenon generally exists when the drainage effect is not up to standard [8]. How to scientifically predict the effect under different engineering design parameters of gas drainage (such as borehole spacing, diameter, and negative pressure) and then quantitatively design the best engineering parameters to meet the predrainage target of coal gas is a key scientific problem of gas drainage engineering. Domestic and foreign scholars have paid great attention to and urgently expect to solve it in recent years.Numerous studies have shown that gas extraction is a complex multifield coupled process, which involves gas transport and solid deformation of the coal body, with the various processes interlocking and influencing each other [9–13]. In the experimental and theoretical studies on the multifield coupled flow of coal gas, Yang et al. introduced the coupled equations of the evolution of stress, damage, and permeability coefficients induced by the mining disturbance influence; then, a model for the coupled solid-gas action of the coal rock fracture process was established [14]. Lu et al. studied the effective radius of borehole gas drainage by the numerical simulation method [15]. Li et al. proposed to determine the relationship between the borehole spacing and the effective extraction radius by using the method of gas content in the coal body around boreholes at a specific moment [16]. Guo et al. established a gas-solid coupling model based on the characteristics of double porosity medium structure by considering the influence of double porosity medium and gas slippage effect of coal body, carried out numerical calculation and analysis by the finite element method, and gave the reasonable extraction spacing of boreholes [17]. Hu et al. considered the Klinkenberg effect of gas flow in low-permeability coal medium and the dynamic change of the coal pore structure with mining stress and gas pressure, respectively, obtained the dynamic evolution model of coal permeability and porosity, and deduced the solid-gas coupling model of low-permeability coal seam [18]. Aiming at the reasonable layout of gas drainage boreholes in coal seam, Wang et al. analyzed the influence radius of single-borehole gas drainage by establishing the seepage field control equation and coal seam deformation field control equation of borehole gas drainage [19]. Hao et al. established a seepage-stress coupling model considering the rheological characteristics, dynamic changes of permeability, and adsorption characteristics of coal and determined the reasonable drainage radius of coal seam [20]. Liang et al. established a solid-gas coupling model to simulate the effective drainage radius of No. 4 coal seam in Shaqu mine, based on the processes of fracture gas seepage, micropore adsorption, desorption and diffusion, and coal-rock mass deformation [21]. We firstly consider the objective gas leakage characteristics during the gas extraction process and coupled the multiphysical processes under the action of gas drainage such as coal deformation, matrix gas seepage, and fracture air-gas mixed flow, and the 2d simulation of gas drainage was carried out by using a single borehole [22].To sum up, the research of domestic and foreign scholars on the gas drainage fluid structure coupling model enriches the gas multiphysical coupling theory but rarely considers the engineering practical problem of drainage concentration attenuation caused by fracture leakage around the borehole. At the same time, the simulation of the gas drainage effect of 3d boreholes is rarely considered, while the gas drainage effect of single borehole simulation has a great deviation from the engineering practice. In order to solve the above problems, based on the multiphysical coupling model of gas drainage developed in the early stage [22], this paper puts forward the content, principle, and method of parameter design of gas drainage engineering and focuses on the characteristics of gas flow in coal seam area and the influence of different gas storage characteristics and drainage parameters on gas drainage effect. On this basis, the 18401 working face of Tunlan Coal Mine, which belongs to Xishan Coal Electricity, is quantitatively designed, which provides theoretical and technical support for the parameter design of gas drainage engineering.


## 2. The Extraction Characteristics and Quantitative Design Model of the Boreholes

### 2.1. Analysis on Interference Characteristics of Borehole Extraction

Coal seams are a typical pore-fracture system, and gas exists mainly in adsorbed and free states in the pore fracture of the coal body, of which about 90% of the gas is adsorbed in the micropores of the coal matrix and its surface and less than 10% of the gas is free in the fracture and macropores [23, 24]. In the process of gas extraction, free gas from the coal matrix and fracture system as well as outside air from the wall of the roadway seeps into the coal seam borehole due to the drainage effect of negative extraction pressure. At the same time, the gas desorption-diffusion in the matrix causes the coal deformation and then affects the permeability of the coal matrix and fracture system [25]. The distribution of residual gas content around single and double boreholes is shown in Figure 1 [25]. In the figure, r0 indicates the radius of the borehole, L is the effective extraction length of the predrainage borehole, W0 is the original gas content of the coal seam, and Wc is the gas content at the boundary of the coal seam borehole. W is the gas content at distance r from the center of extraction borehole; re1andre2 are the radius of influence of single-hole extraction of Z1andZ2, respectively; Lz is the spacing between Z1andZ2 boreholes; W1andW2 are the residual gas content in the area around the borehole after single-hole extraction of Z1andZ2, respectively; ΔW1 is the gas extraction capacity per ton of coal at radial distance r from single-hole Z1; and ΔW2 is the gas extraction capacity per ton of coal at radial distance Lz‐r from single-hole Z2.Figure 1
Evolution diagram of residual gas content around single and double boreholes [25].
(a)
Single drilling(b)
Adjacent drillingThe relationship between total extraction volumeQ and radial gas content at a specific point in the borehole extraction can be expressed as [26]
(1)Q=∫r0reW0−WrρL2πrdr,where Q is the gas extraction volume of the borehole (m3) and ρ is the apparent density of coal (m3/t).In the process of gas extraction of the boreholes, if the sum of the influence radius of adjacent holes is greater than the distance between holes, the gas extraction effect of the boreholes will be affected by the superposition effect of adjacent holes, and the pressure and content of coal seam gas extraction will decay faster. At the same time, because the pressure and content decay faster, there is a certain interference in gas extraction between boreholes. Then, the gas content atr is
(2)Wr=W0−∆W1+∆W2+∆W′,Lz≤re1+re2,where ΔW1 is the gas extraction capacity per ton of coal at radial distance r of the single-borehole Z1, ΔW2 is the gas extraction capacity per ton of coal at radial distance Lz‐r of the single borehole Z2, and ΔW′ is the relative attenuation of gas extraction capacity per ton of coal at radial distance r of the borehole Z1, which is caused by borehole interference.ΔW1 and ΔW2 in the above equation can be expressed as
(3)∆W1=W0−W1,∆W2=W0−W2.If the sum of the radius of influence of adjacent boreholes is less than the spacing of boreholes, it means that there is no superposition effect of extraction effect between extraction boreholes, and the actual effect of gas extraction is the same as that of single-borehole extraction. That is,(4)Wr=W1,Lz>re1+re2.Substituting equation (3) into equation (2) yields the relative gas extraction capacity attenuation per ton of coal due to borehole disturbance:
(5)∆W′=Wr+W0−W1−W2,Lz≤re1+re2.To quantitatively analyze the cooperative and competitive characteristics between boreholes, for a certain extraction time, the cooperative efficiency of gas extraction at a distancer from a borehole is defined as the ratio between the reduction in gas content extracted by the actual boreholes and the reduction in gas content without considering the interference between boreholes, that is,
(6)ηr=0,Lz>re1+re2,W0−Wr∆W1+∆W2,Lz≤re1+re2.


### 2.2. Calculation Model of Design Parameters

The flow of gas in the coal body under the action of gas extraction is the result of the coupling of multiple processes such as coal deformation, matrix gas seepage, and fracture air-gas mixing.(1)
Deformation control equations for dual media coal under ground stress and gas adsorption/desorption stress effects [10]:(7)Gui,kk+G1−2υuk,ki−α+K·εLPLPL+pm12pm1,i−βpf1+pf2,i+fi=0,where ui is the displacement component (m), G is the shear modulus of coal (MPa), ν is Poisson’s ratio of coal, K is the bulk modulus of coal (MPa), α and β are the Biot coefficients of matrix and fracture systems, respectively, εL is gas Langmuir volumetric strain constant, PL is the gas Langmuir pressure constants (MPa), fi is the component of physical force (N), pm1andpf1 are the matrix and fracture gas pressure of the system (MPa), and pf2 is the pressure of the air in the fracture system (MPa)
(2)
Control equations for gas flow transfer in matrix and fracture systems [22](8)A∂pm1∂t−∇·kmμpm1+φmDm1∇pm1=−81+2a2kmμpm1−pf1−B·pm1,φf∂pf1∂t−∇·kfμpf1+φfDf1∇pf=81+2a2kmμpm1−pf1−C·pf1,where pm1andpf1 are the gas pressure in the matrix and fracture system, respectively (MPa); kmandkf are the permeability of the coal matrix and fracture (m2); φm and φf are the porosity of the coal matrix and fracture; Dm1andDf1 are the diffusion coefficients of the matrix and fracture gases (m2/s); μ is the gas dynamic viscosity coefficient (N·s/m2); t is time (s); a is the side length of the cubic coal matrix block (m); and the physical expressions of A, B, and C are
(9)A=φm+ρspaVLPLPL+pm12−φm−αKεLPLpm1PL+pm121K+b0a0Kf−1,B=φm−αK1K+b0a0Kf−1∂εV∂t,C=φf0Kf1K+b0a0Kf−1∂εV∂t−εLPLPL+pm12∂pm1∂t,where εV=ε11+ε22+ε33 is the volume strain. The subscript, 0, denotes the initial value of the variable
(3)
Crack system and air flow control equations [26]:(10)φf∂pf2∂t−∇·kfμpf2∇pf+φfDf1∇pf2=−C·pf2,where pf2 is the air pressure of the matrix system (MPa)
(4)
Cross-coupled variable control equationsThe permeability of the coal matrix and fractures can be expressed using the cubic law as [22]
(11)kmkm0=φmφm03,kfkf0=φfφf03.The porosity of the coal matrix and fractures can be expressed as [22]
(12)φm=α+φm0−αexp1K1K+b0a0Kf−1ΔεV−Δεs,φf=φf0+φf0Kf1K+b0a0Kf−1ΔεV−Δεs.Equations (7)–(12) comprise a coal gas extraction coupled flow-solid model. Combining equations (1)–(6), the dynamic evolution of coal seam gas flow parameters, gas pressure, and content can be solved, and the borehole coordination and interference characteristics of coal seam prepumping can be revealed, which provides a theoretical basis for the quantitative design of coal seam gas drainage engineering parameters.


### 2.3. Optimization of Design Constraints

Gas extraction is a fundamental measure to manage gas disasters, and its primary purpose is to reduce the pressure and content of coal seam gas; eliminate the risk of coal and gas outburst, gas overrun, gas flare, and other risks in the coal mining process; and realize safe simultaneous extraction of coal and gas. The quantitative design of coal seam gas extraction mainly refers to the quantitative design of engineering parameters such as the borehole diameterd of gas drainage, the effective extraction length of borehole L, the negative pressure of borehole ∆p, and the borehole spacing D. The main idea of quantified gas extraction design is to determine the gas prepumping time t0, according to the mining succession arrangement first. Then, take the key indicators of gas extraction (gas pressure p, content m, and extraction rate η) and gas extraction concentration c (CH4) as constraints. At last, based on the mathematical model of gas extraction flow, the most reasonable engineering parameters of coal seam gas extraction (drilling diameter d, effective extraction length L, negative pressure ∆p, and spacing D) are quantitatively determined. Thus, the constraints for the quantitative design of gas extraction engineering parameters are
(13)cCH4>cCH4LmCH4<mCH4LpCH4<pCH4Lη=QmcQmc+Qmf>ηLt=t0⟶S.t∆p,d,LMaxD,where c is the gas extraction concentration (%), m is the coal seam gas content (m3/t), p is the coal matrix gas pressure (Pa), η is the gas extraction rate (%), t is the gas extraction time (d), t0 is the gas prepumping time (d), the subscript “L” is the expected minimum target value for extraction, S.t is the target constraint, and MaxD is the maximum drilling spacing (m), in which the gas extraction safety meets the standard within the expected period of gas extraction.


### 2.4. Solution of Nonlinear Constrained Optimization Problems

This is a constrained nonlinear programming problem consisting of multifield coupled computational equations, objective function, and constraint equations, which have great complexity. Comsol With Matlab has powerful numerical analysis and nonlinear optimization compilation functions, which provide a powerful computational support platform for solving the optimization problem of equations (1)–(13).The main purpose of the quantitative design of coal seam gas extraction is to achieve the harmonization of safety and economy and to achieve the maximum arrangement of drilling spacing under the premise of satisfying safety; the basic steps are shown in Figure2.
(1)
Determine the prepumping timet0 of the borehole and the minimum borehole spacing D0 (Dmin) based on empirical trials.(2)
The basic parameters of borehole extraction as well as the basic occurrence parameters of coal seam gas are substituted into the multifield coupled calculation model. The values of each parameter are solved when the prepumping timet0 is reached and to determine whether the constraints are satisfied.(3)
The conditions for the constraint of the drilling extraction parameters can be determined based on the gas source to meet the parameters and their critical values. If the constraint conditions are met, increase the drilling spacing and continue to calculate according to the second step until the constraint conditions cannot be met.(4)
If the constraints are not met, a smaller spacing of boreholes is selected to carry out the cycle above.(5)
Compare the boreholes that meet the borehole drainage constraints in different borehole spacing and determine the maximum gas drainage spacing and output.Figure 2
Borehole spacing optimization design process.


## 2.1. Analysis on Interference Characteristics of Borehole Extraction

Coal seams are a typical pore-fracture system, and gas exists mainly in adsorbed and free states in the pore fracture of the coal body, of which about 90% of the gas is adsorbed in the micropores of the coal matrix and its surface and less than 10% of the gas is free in the fracture and macropores [23, 24]. In the process of gas extraction, free gas from the coal matrix and fracture system as well as outside air from the wall of the roadway seeps into the coal seam borehole due to the drainage effect of negative extraction pressure. At the same time, the gas desorption-diffusion in the matrix causes the coal deformation and then affects the permeability of the coal matrix and fracture system [25]. The distribution of residual gas content around single and double boreholes is shown in Figure 1 [25]. In the figure, r0 indicates the radius of the borehole, L is the effective extraction length of the predrainage borehole, W0 is the original gas content of the coal seam, and Wc is the gas content at the boundary of the coal seam borehole. W is the gas content at distance r from the center of extraction borehole; re1andre2 are the radius of influence of single-hole extraction of Z1andZ2, respectively; Lz is the spacing between Z1andZ2 boreholes; W1andW2 are the residual gas content in the area around the borehole after single-hole extraction of Z1andZ2, respectively; ΔW1 is the gas extraction capacity per ton of coal at radial distance r from single-hole Z1; and ΔW2 is the gas extraction capacity per ton of coal at radial distance Lz‐r from single-hole Z2.Figure 1
Evolution diagram of residual gas content around single and double boreholes [25].
(a)
Single drilling(b)
Adjacent drillingThe relationship between total extraction volumeQ and radial gas content at a specific point in the borehole extraction can be expressed as [26]
(1)Q=∫r0reW0−WrρL2πrdr,where Q is the gas extraction volume of the borehole (m3) and ρ is the apparent density of coal (m3/t).In the process of gas extraction of the boreholes, if the sum of the influence radius of adjacent holes is greater than the distance between holes, the gas extraction effect of the boreholes will be affected by the superposition effect of adjacent holes, and the pressure and content of coal seam gas extraction will decay faster. At the same time, because the pressure and content decay faster, there is a certain interference in gas extraction between boreholes. Then, the gas content atr is
(2)Wr=W0−∆W1+∆W2+∆W′,Lz≤re1+re2,where ΔW1 is the gas extraction capacity per ton of coal at radial distance r of the single-borehole Z1, ΔW2 is the gas extraction capacity per ton of coal at radial distance Lz‐r of the single borehole Z2, and ΔW′ is the relative attenuation of gas extraction capacity per ton of coal at radial distance r of the borehole Z1, which is caused by borehole interference.ΔW1 and ΔW2 in the above equation can be expressed as
(3)∆W1=W0−W1,∆W2=W0−W2.If the sum of the radius of influence of adjacent boreholes is less than the spacing of boreholes, it means that there is no superposition effect of extraction effect between extraction boreholes, and the actual effect of gas extraction is the same as that of single-borehole extraction. That is,(4)Wr=W1,Lz>re1+re2.Substituting equation (3) into equation (2) yields the relative gas extraction capacity attenuation per ton of coal due to borehole disturbance:
(5)∆W′=Wr+W0−W1−W2,Lz≤re1+re2.To quantitatively analyze the cooperative and competitive characteristics between boreholes, for a certain extraction time, the cooperative efficiency of gas extraction at a distancer from a borehole is defined as the ratio between the reduction in gas content extracted by the actual boreholes and the reduction in gas content without considering the interference between boreholes, that is,
(6)ηr=0,Lz>re1+re2,W0−Wr∆W1+∆W2,Lz≤re1+re2.


## 2.2. Calculation Model of Design Parameters

The flow of gas in the coal body under the action of gas extraction is the result of the coupling of multiple processes such as coal deformation, matrix gas seepage, and fracture air-gas mixing.(1)
Deformation control equations for dual media coal under ground stress and gas adsorption/desorption stress effects [10]:(7)Gui,kk+G1−2υuk,ki−α+K·εLPLPL+pm12pm1,i−βpf1+pf2,i+fi=0,where ui is the displacement component (m), G is the shear modulus of coal (MPa), ν is Poisson’s ratio of coal, K is the bulk modulus of coal (MPa), α and β are the Biot coefficients of matrix and fracture systems, respectively, εL is gas Langmuir volumetric strain constant, PL is the gas Langmuir pressure constants (MPa), fi is the component of physical force (N), pm1andpf1 are the matrix and fracture gas pressure of the system (MPa), and pf2 is the pressure of the air in the fracture system (MPa)
(2)
Control equations for gas flow transfer in matrix and fracture systems [22](8)A∂pm1∂t−∇·kmμpm1+φmDm1∇pm1=−81+2a2kmμpm1−pf1−B·pm1,φf∂pf1∂t−∇·kfμpf1+φfDf1∇pf=81+2a2kmμpm1−pf1−C·pf1,where pm1andpf1 are the gas pressure in the matrix and fracture system, respectively (MPa); kmandkf are the permeability of the coal matrix and fracture (m2); φm and φf are the porosity of the coal matrix and fracture; Dm1andDf1 are the diffusion coefficients of the matrix and fracture gases (m2/s); μ is the gas dynamic viscosity coefficient (N·s/m2); t is time (s); a is the side length of the cubic coal matrix block (m); and the physical expressions of A, B, and C are
(9)A=φm+ρspaVLPLPL+pm12−φm−αKεLPLpm1PL+pm121K+b0a0Kf−1,B=φm−αK1K+b0a0Kf−1∂εV∂t,C=φf0Kf1K+b0a0Kf−1∂εV∂t−εLPLPL+pm12∂pm1∂t,where εV=ε11+ε22+ε33 is the volume strain. The subscript, 0, denotes the initial value of the variable
(3)
Crack system and air flow control equations [26]:(10)φf∂pf2∂t−∇·kfμpf2∇pf+φfDf1∇pf2=−C·pf2,where pf2 is the air pressure of the matrix system (MPa)
(4)
Cross-coupled variable control equationsThe permeability of the coal matrix and fractures can be expressed using the cubic law as [22]
(11)kmkm0=φmφm03,kfkf0=φfφf03.The porosity of the coal matrix and fractures can be expressed as [22]
(12)φm=α+φm0−αexp1K1K+b0a0Kf−1ΔεV−Δεs,φf=φf0+φf0Kf1K+b0a0Kf−1ΔεV−Δεs.Equations (7)–(12) comprise a coal gas extraction coupled flow-solid model. Combining equations (1)–(6), the dynamic evolution of coal seam gas flow parameters, gas pressure, and content can be solved, and the borehole coordination and interference characteristics of coal seam prepumping can be revealed, which provides a theoretical basis for the quantitative design of coal seam gas drainage engineering parameters.


## 2.3. Optimization of Design Constraints

Gas extraction is a fundamental measure to manage gas disasters, and its primary purpose is to reduce the pressure and content of coal seam gas; eliminate the risk of coal and gas outburst, gas overrun, gas flare, and other risks in the coal mining process; and realize safe simultaneous extraction of coal and gas. The quantitative design of coal seam gas extraction mainly refers to the quantitative design of engineering parameters such as the borehole diameterd of gas drainage, the effective extraction length of borehole L, the negative pressure of borehole ∆p, and the borehole spacing D. The main idea of quantified gas extraction design is to determine the gas prepumping time t0, according to the mining succession arrangement first. Then, take the key indicators of gas extraction (gas pressure p, content m, and extraction rate η) and gas extraction concentration c (CH4) as constraints. At last, based on the mathematical model of gas extraction flow, the most reasonable engineering parameters of coal seam gas extraction (drilling diameter d, effective extraction length L, negative pressure ∆p, and spacing D) are quantitatively determined. Thus, the constraints for the quantitative design of gas extraction engineering parameters are
(13)cCH4>cCH4LmCH4<mCH4LpCH4<pCH4Lη=QmcQmc+Qmf>ηLt=t0⟶S.t∆p,d,LMaxD,where c is the gas extraction concentration (%), m is the coal seam gas content (m3/t), p is the coal matrix gas pressure (Pa), η is the gas extraction rate (%), t is the gas extraction time (d), t0 is the gas prepumping time (d), the subscript “L” is the expected minimum target value for extraction, S.t is the target constraint, and MaxD is the maximum drilling spacing (m), in which the gas extraction safety meets the standard within the expected period of gas extraction.


## 2.4. Solution of Nonlinear Constrained Optimization Problems

This is a constrained nonlinear programming problem consisting of multifield coupled computational equations, objective function, and constraint equations, which have great complexity. Comsol With Matlab has powerful numerical analysis and nonlinear optimization compilation functions, which provide a powerful computational support platform for solving the optimization problem of equations (1)–(13).The main purpose of the quantitative design of coal seam gas extraction is to achieve the harmonization of safety and economy and to achieve the maximum arrangement of drilling spacing under the premise of satisfying safety; the basic steps are shown in Figure2.
(1)
Determine the prepumping timet0 of the borehole and the minimum borehole spacing D0 (Dmin) based on empirical trials.(2)
The basic parameters of borehole extraction as well as the basic occurrence parameters of coal seam gas are substituted into the multifield coupled calculation model. The values of each parameter are solved when the prepumping timet0 is reached and to determine whether the constraints are satisfied.(3)
The conditions for the constraint of the drilling extraction parameters can be determined based on the gas source to meet the parameters and their critical values. If the constraint conditions are met, increase the drilling spacing and continue to calculate according to the second step until the constraint conditions cannot be met.(4)
If the constraints are not met, a smaller spacing of boreholes is selected to carry out the cycle above.(5)
Compare the boreholes that meet the borehole drainage constraints in different borehole spacing and determine the maximum gas drainage spacing and output.Figure 2
Borehole spacing optimization design process.


## 3. Model Verification and Simulation

### 3.1. Model Parameters and Initial Boundary

The strike length of 8# coal seam in Tunlan 18401 working face is 1760 m, the dip length is 235 m, the average thickness of coal seam is 3 m, the dip angle is 4°, and the average mining depth is 750 m. The original gas content and gas pressure of the coal seam are 13.4 m3/t and 1.9 MPa, respectively, and the permeability coefficient of the coal seam is 3.63 m2/(MPa2·d). As shown in Figure 3, the physical model of gas extraction from a three-dimensional hole cluster and its boundary conditions are illustrated. The basic parameters of gas drainage calculation are shown in Table 1.Figure 3
Schematic diagram of numerical model of gas extraction borehole.Table 1
Basic parameters of gas extraction simulation.
ParameterValueParameterValueYoung’s modulus of coal (E, MPa)2813Initial permeability of coal fracture (kf0, m2)4.0×10−18Young’s modulus of coal grains (Es, MPa)8439Normal atmospheric pressure (pa, MPa)0.101325Fracture strength of coal (Kn, MPa)4800Langmuir pressure constant (PL, MPa)0.922Poisson’s ratio of coal (υ)0.339Langmuir volume constant (VL, m3/kg)0.023Density of coal (ρs, kg/m3)1250Langmuir volumetric strain constant (εL)0.01266Initial porosity of the coal matrix (φm)0.05Dynamic viscosity of methane (μ, N·s/m2)1.84×10−5Initial permeability of coal matrix (km0, m2)8.9×10−20Density of methane at standard condition (ρs, kg/m3)0.717Initial gas pressure of coal seam (p0, MPa)1.9Air leakage (N0, m/s)4×10−4


### 3.2. Model Verification and Result Analysis

Figure4 shows the distribution of coal seam gas content after 60 days of drilling. It can be seen from the figure that the influence range of gas drainage of the array boreholes gradually extends to the axial center line of the two boreholes and symmetrically distributed on both sides. The gas content at the center line of the two boreholes is the largest, and that at the borehole is the smallest. Extraction volume around single-borehole ∆W1=∆W2=W0−W1=W0−W2 is significantly less than the borehole collaborative extraction volume W0−Wr, and the superposition effect of two single-borehole extraction ∆W1+∆W2 is greater than the actual extraction content of the borehole extraction W0−Wr, which verifies that there are cooperative and competitive extraction behaviors among boreholes.Figure 4
The gas content distribution of coal seam extracted by the boreholes.
(a)
The change of gas content after 60 days of gas extraction(b)
Residual gas content between point (0, 1201, 1.5) and point (15, 120, 1.5)Figure5 shows the field example matching gas drainage results of the single borehole and boreholes. It can be seen from the figure that under the action of gas drainage, the net flow and concentration of borehole gas drainage decay exponentially. The gas flow concentration under borehole drainage can match the actual situation well. Under the single-borehole gas extraction, the pure flow and concentration of gas extracted from the borehole are relatively high for the same simulation conditions, as the hole cluster effect is ignored. Therefore, in the parameter design of on-site gas drainage engineering, the spacing of gas drainage is often too large, and the extraction effect is far from reaching the expected effect. So, it is necessary to fully consider the gas drainage effect of the boreholes to improve the accuracy of parameter design of gas drainage engineering.Figure 5
Matching gas drainage results of the single borehole and boreholes in the field.
(a)(b)


## 3.1. Model Parameters and Initial Boundary

The strike length of 8# coal seam in Tunlan 18401 working face is 1760 m, the dip length is 235 m, the average thickness of coal seam is 3 m, the dip angle is 4°, and the average mining depth is 750 m. The original gas content and gas pressure of the coal seam are 13.4 m3/t and 1.9 MPa, respectively, and the permeability coefficient of the coal seam is 3.63 m2/(MPa2·d). As shown in Figure 3, the physical model of gas extraction from a three-dimensional hole cluster and its boundary conditions are illustrated. The basic parameters of gas drainage calculation are shown in Table 1.Figure 3
Schematic diagram of numerical model of gas extraction borehole.Table 1
Basic parameters of gas extraction simulation.
ParameterValueParameterValueYoung’s modulus of coal (E, MPa)2813Initial permeability of coal fracture (kf0, m2)4.0×10−18Young’s modulus of coal grains (Es, MPa)8439Normal atmospheric pressure (pa, MPa)0.101325Fracture strength of coal (Kn, MPa)4800Langmuir pressure constant (PL, MPa)0.922Poisson’s ratio of coal (υ)0.339Langmuir volume constant (VL, m3/kg)0.023Density of coal (ρs, kg/m3)1250Langmuir volumetric strain constant (εL)0.01266Initial porosity of the coal matrix (φm)0.05Dynamic viscosity of methane (μ, N·s/m2)1.84×10−5Initial permeability of coal matrix (km0, m2)8.9×10−20Density of methane at standard condition (ρs, kg/m3)0.717Initial gas pressure of coal seam (p0, MPa)1.9Air leakage (N0, m/s)4×10−4


## 3.2. Model Verification and Result Analysis

Figure4 shows the distribution of coal seam gas content after 60 days of drilling. It can be seen from the figure that the influence range of gas drainage of the array boreholes gradually extends to the axial center line of the two boreholes and symmetrically distributed on both sides. The gas content at the center line of the two boreholes is the largest, and that at the borehole is the smallest. Extraction volume around single-borehole ∆W1=∆W2=W0−W1=W0−W2 is significantly less than the borehole collaborative extraction volume W0−Wr, and the superposition effect of two single-borehole extraction ∆W1+∆W2 is greater than the actual extraction content of the borehole extraction W0−Wr, which verifies that there are cooperative and competitive extraction behaviors among boreholes.Figure 4
The gas content distribution of coal seam extracted by the boreholes.
(a)
The change of gas content after 60 days of gas extraction(b)
Residual gas content between point (0, 1201, 1.5) and point (15, 120, 1.5)Figure5 shows the field example matching gas drainage results of the single borehole and boreholes. It can be seen from the figure that under the action of gas drainage, the net flow and concentration of borehole gas drainage decay exponentially. The gas flow concentration under borehole drainage can match the actual situation well. Under the single-borehole gas extraction, the pure flow and concentration of gas extracted from the borehole are relatively high for the same simulation conditions, as the hole cluster effect is ignored. Therefore, in the parameter design of on-site gas drainage engineering, the spacing of gas drainage is often too large, and the extraction effect is far from reaching the expected effect. So, it is necessary to fully consider the gas drainage effect of the boreholes to improve the accuracy of parameter design of gas drainage engineering.Figure 5
Matching gas drainage results of the single borehole and boreholes in the field.
(a)(b)


## 4. Influencing Factors of Gas Drainage Effect Using Boreholes

### 4.1. Influence of Borehole Spacing on Drainage Effect

Figure6 shows the simulation results of gas drainage using boreholes when the borehole spacing D is 2 m, 4 m, and 6 m, respectively. It can be seen from the figure that the smaller the borehole spacing, the stronger the competition and synergy between boreholes, the smaller the gas flow and concentration, the faster the attenuation of coal seam gas pressure, and the relatively high gas extraction rate. For example, if the drilling spacing D is 2 m, the gas flow and extraction rate are 0.020 m3/min and 37.3%, respectively, after 100 days of extraction, while the gas flow and extraction rate are 0.028 m3/min and 21.9%, respectively, when D is 6 m.Figure 6
The influence of the borehole spacing on gas drainage effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


### 4.2. The Influence of Drainage Negative Pressure on Drainage Effect

Figure7 shows the simulation results of gas drainage using boreholes when the negative pressure Δp is 13 kPa, 20 kPa, and 25 kPa, respectively. It can be seen from the figure that the negative pressure of extraction has no obvious impact on the competition and synergy between boreholes but has a significant impact on gas drainage concentration. The negative pressure has little effect on the net amount of gas drainage, the change of coal seam gas pressure, and the extraction rate but has a greater impact on the concentration of gas drainage. The greater the negative pressure, the greater the gas leakage of the borehole, resulting in the faster attenuation of the concentration of gas drainage. Such as Δp=13kPa and Δp=25kPa, it takes 205.8 days and 152.6 days for gas concentration to decay from the initial state to less than 30%, respectively.Figure 7
Influence of different negative pressure on gas drainage effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


### 4.3. Influence of Permeability on Drainage Effect

Figure8 shows the simulation results of gas drainage using boreholes under different initial permeability of coal matrix km0 and initial permeability of coal fracture kf0. It can be seen from the figure that the change in initial permeability of coal matrix km0 has no significant impact on the competition and synergy between boreholes. While the greater the initial fracture permeability kf0, the stronger the gas migration capacity of the fracture system, the stronger the competition effect between boreholes, the higher the purity and concentration of gas drainage, and the greater the gas drainage rate. For example, when the initial permeability of the coal matrix is the same value km0=8.9×10‐20m2, within 100 days of gas drainage time, the corresponding gas extraction rates of different coal fracture initial permeability kf0=4.0×10‐18 and kf0=4.0×10‐17 are 20.2% and 45.8%, respectively. Under the same extraction time and the same coal fracture initial permeability kf0=8.9×10‐18m2, the extraction rates of different coal matrix initial permeability km0=4.0×10‐20 and km0=4.0×10‐19m2 are 20.2% and 20.4%, respectively. This shows that the transmission of gas in the coal seam is mainly determined by the fracture system, and the matrix system mainly determines the diffusion process of coal, which has no obvious effect on the improvement of gas transmission capacity.Figure 8
Influence of different permeability of coal seam on gas drainage effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


## 4.1. Influence of Borehole Spacing on Drainage Effect

Figure6 shows the simulation results of gas drainage using boreholes when the borehole spacing D is 2 m, 4 m, and 6 m, respectively. It can be seen from the figure that the smaller the borehole spacing, the stronger the competition and synergy between boreholes, the smaller the gas flow and concentration, the faster the attenuation of coal seam gas pressure, and the relatively high gas extraction rate. For example, if the drilling spacing D is 2 m, the gas flow and extraction rate are 0.020 m3/min and 37.3%, respectively, after 100 days of extraction, while the gas flow and extraction rate are 0.028 m3/min and 21.9%, respectively, when D is 6 m.Figure 6
The influence of the borehole spacing on gas drainage effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


## 4.2. The Influence of Drainage Negative Pressure on Drainage Effect

Figure7 shows the simulation results of gas drainage using boreholes when the negative pressure Δp is 13 kPa, 20 kPa, and 25 kPa, respectively. It can be seen from the figure that the negative pressure of extraction has no obvious impact on the competition and synergy between boreholes but has a significant impact on gas drainage concentration. The negative pressure has little effect on the net amount of gas drainage, the change of coal seam gas pressure, and the extraction rate but has a greater impact on the concentration of gas drainage. The greater the negative pressure, the greater the gas leakage of the borehole, resulting in the faster attenuation of the concentration of gas drainage. Such as Δp=13kPa and Δp=25kPa, it takes 205.8 days and 152.6 days for gas concentration to decay from the initial state to less than 30%, respectively.Figure 7
Influence of different negative pressure on gas drainage effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


## 4.3. Influence of Permeability on Drainage Effect

Figure8 shows the simulation results of gas drainage using boreholes under different initial permeability of coal matrix km0 and initial permeability of coal fracture kf0. It can be seen from the figure that the change in initial permeability of coal matrix km0 has no significant impact on the competition and synergy between boreholes. While the greater the initial fracture permeability kf0, the stronger the gas migration capacity of the fracture system, the stronger the competition effect between boreholes, the higher the purity and concentration of gas drainage, and the greater the gas drainage rate. For example, when the initial permeability of the coal matrix is the same value km0=8.9×10‐20m2, within 100 days of gas drainage time, the corresponding gas extraction rates of different coal fracture initial permeability kf0=4.0×10‐18 and kf0=4.0×10‐17 are 20.2% and 45.8%, respectively. Under the same extraction time and the same coal fracture initial permeability kf0=8.9×10‐18m2, the extraction rates of different coal matrix initial permeability km0=4.0×10‐20 and km0=4.0×10‐19m2 are 20.2% and 20.4%, respectively. This shows that the transmission of gas in the coal seam is mainly determined by the fracture system, and the matrix system mainly determines the diffusion process of coal, which has no obvious effect on the improvement of gas transmission capacity.Figure 8
Influence of different permeability of coal seam on gas drainage effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


## 5. The Influence of Gas Adsorption Characteristics on Drainage Effect

Figure9 shows the simulation results of gas drainage using borehole pumping with different volume constants VL and pressure constants PL. It can be seen from the figure that the gas adsorption constant of coal has a significant impact on the competition and synergy between boreholes. The smaller the pressure constant PL and the larger the volume constant VL, the higher the net amount and concentration of gas drainage, the faster the attenuation rate of coal seam gas pressure, and the greater the gas drainage rate. For example, when the volume constant VL is 0.023 m3/kg, the corresponding drainage gas flow of different pressure constants PL is 0.922 MPa and 6.922 MPa for 60 days is 0.0322 and 0.0171 m3/min, respectively. The corresponding gas drainage flow of different volume constants VL=0.023 and 0.043 m3/kg is 0.0322 and 0.0443 m3/min, respectively.Figure 9
Influence of gas adsorption constant on gas extraction effect.
(a)
Gas drainage flow rate(b)
Gas drainage concentration(c)
Cooperative efficiency of gas extraction(d)
Gas extraction rate


## 6. Example of Optimized Design of Gas Drainage

The coal seam structure of 18401 in Tulan coal mine is complex, and there is a risk of coal and gas protrusion. In order to prevent the occurrence of gas overrun accidents during backstopping, gas extraction and management measures must be used to ensure the safe production of the working face. According to the working face drainage standard conditions, the residual gas pressurepCH4L should be lower than 0.74 MPa, the residual gas content mCH4L should be lower than 8 m3/t, and the gas extraction rate η should be higher than 40% are selected as the target constraints of the parameters. In the design of the drill hole spacing for the 18401 working face, the range of drill hole spacing D is calculated from 2 m to 6 m, the calculation time t0 is 900 days, the initial spacing D0 is 2 m, and the calculation step is 1 m. At each calculation step, the software automatically judges whether the extraction index meets the standard according to the calculation result until the end, and the maximum drill hole spacing D that meets all indexes is preferably selected from the range of drill hole spacing. The simulation results of different drill hole spacing are shown in Table 2.Table 2
Prediction of gas drainage effect at different spacing of boreholes.
IndexExpected valueBorehole spacing2 m3 m4 m5 m6 mPredicted valuePredicted valuePredicted valuePredicted valuePredicted valueResidual gas volume (m3/t)<84.386.146.997.768.66Gas extraction rate (%)>4067.354.247.842.135.4Maximum coal matrix gas pressure (MPa)<0.740.0970.1190.1360.1520.169Are all constraints met?✓✓✓✓✓According to the simulation prediction results in Table2, the residual gas content and gas extraction rate of coal can no longer meet the extraction target when the drill hole spacing is greater than 6 m, and the maximum drill hole spacing to meet all constraints is 5 m. Therefore, in order to ensure the safe production of mine, the drill hole spacing should be less than 5 m when designing drill holes. In consideration of the actual excavation schedule of the working face, the drilling spacing inward the second setup room of the track roadway and at 765 m inward from the alley entrance of the air-return measure roadway of the belt roadway are designed according to 3 m, and it is designed as 5 m at 38 m inward from the alley entrance of the return air measures of track roadway. The specific layout is shown in Figure 10 and Table 3.Figure 10
Extracting seam drilling construction layout.Table 3
Parameters of the extracting seam drilling construction layout.
LocationHole depthDip anglePositionApertureNumber of boreholesTotal footageOpening heightSpacingTrack roadwayAir-return measure roadway 38 meters inward220 m1~2°90°113 mm11825960 m1.5 m5 mInward from the second setup room150 m1~2°90°113 mm30746050 m1.5 m3 mBelt roadwayAir-return measure roadway 765 meters inward90 m-1~-2°90°113 mm30727630 m1.5 m3 mFigure11 shows the gas extraction data of 18401 working face during 2017.06.07~2017.12.31. From the figure, it can be seen that the cumulative volume of gas extraction in the track roadway is 2,709,400 m3, the average pure volume of extraction is 9.32 m3/min, and the average extraction concentration is 31.53%; the cumulative volume of gas extraction in the belt roadway is 205,100 m3, the average pure volume of extraction is 0.69 m3/min, and the average extraction concentration is 28.24%. After extraction, the coal seam gas content of 18401 working face is obviously reduced, and the extraction effect is relatively satisfactory.Figure 11
Variation of gas flow and concentration in track and belt roadway systems.
(a)
Track roadway(b)
Belt roadway


## 7. Conclusions

In this study, the governing equation of coal seam content evolution under the cooperative-competition of gas extraction by the boreholes was established, the calculation formula of gas cooperative efficiency was defined, and the cooperative and competitive behavior of gas extraction by the boreholes was revealed. The results show that drilling spacing, seam fracture permeability, and gas adsorption constant have significant effects on the synergistic and competitive effects of the boreholes, while the negative pressure of gas extraction and coal matrix permeability have little effects. But the negative pressure of gas extraction has great effects on gas extraction concentration.A multifield coupling 3d model of coal seam extraction of the boreholes was developed based on gas-air mixed flow developed. In the model, we put forward the content, method, and principle of the parameter design of the optimal borehole drainage spacing and constructed the constraint conditions and objective function of the parameter optimization design of gas drainage engineering. It provides a theoretical basis for the quantitative design of coal seam gas drainage engineering parameters.Comsol With Matlab was used to solve the multiconstraint nonlinear optimization problem of gas extraction engineering parameters along the coal seam group in 18401 working face of Tunlan Mine, and the optimal drilling layout parameters were determined. From July 6, 2017, to December 31, 2017, the cumulative capacity of gas extracted by the extraction system of track roadway and belt roadway was 2,709,400 m3 and 205,100 m3, respectively, which has achieved an ideal drainage effect.



---

*Source: 5943321-2022-06-28.xml*