Rational Boreholes Arrangement of Gas Extraction from Unloaded Coal Seam

College of Mining, Liaoning Technical University, Fuxin, Liaoning Province, China State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology, Xuzhou, Jiangsu Province, China State Key Laboratory Cultivation Base for Gas Geology and Gas Control, Henan Polytechnic University, Jiaozuo, Henan Province, China Energy and Mineral Engineering, Pennsylvania State University, University Park, PA, USA


Introduction
For coal seams with low gas permeability and high gas, gas extraction measures must be taken before mining to reduce the coal seam gas content.However, in many cases, although a large number of projects have been taken, the gas extraction effect is not ideal.Protective layer mining can effectively reduce the stress and increase the permeability, and the permeability coefficient is hundreds to thousands of times [1,2].Drilling borehole in the floor roadway is one of the most commonly used methods in the protected coal seam gas extraction.And reasonable arrangement of boreholes is key to the efficient gas extraction of coal seam [3,4].
Since 1960s, scholars at home and abroad have carried out various research studies on the technology of protective layer mining and pressure relief gas extraction.Li et al. [5] simulated stress distribution and deformation law of the protected layer in lower protective layer mining and determined protection scopes of the trend and tendency of the protective layer.Zhang et al. [6] made a similar material physical model and analyzed the movement law of overlying strata and effect of the expansion and deformation of the protected layer.Wang et al. [7] studied the evolution law of fracture of coal and rock strata in the upper protective layer and considered that gas extraction should consider the factors of spatial and temporal effects of the evolution of the mining fracture.Wang et al [8] and Liu [9] tested the permeability change of the protected coal seam after protective layer mining and established the coupling dynamic model of the coal rock mass deformation and gas flow.
Zhang et al. [10][11][12][13] combined with eld practice, respectively, used drilling borehole in the oor roadway, highlevel drilling, and ground drilling and combined the method to extract the released gas from the protected layer and achieved good gas control e ect.Wang et al. [14] and Hu et al. [15] studied the gas emission and the pressure relief law of the protected layer in the mining face of the close and steep inclined and steep slope protection layer.On this basis, the parameters of pressure relief gas extraction were optimized.Ma et al. [16,17] investigated the in uence of stress rate and grain size on gas seepage characteristics in crushed coal particles.
Based on the multiphysics eld theory, the hydraulicmechanical coupling model of gas extraction in the unloading coal seam was established, and the gas extraction process by drilling borehole in the oor roadway in the overburden of Panyi Coal Mine 1551 (1) panel was simulated.e in uence of di erent drilling arrangements on the gas extraction e ect was analyzed.

Hydraulic-Mechanical Coupling Model for
Unloaded Gas Extraction 2.1.Basic Hypothesis.After mining the protective layer, the overlying strata will form the caving zone, fracture zone, and subsidence zone (vertical three zones) in the vertical direction, which will cause the pressure relief of the protected coal seam, as shown in Figure 1 [18].e subsidence and deformation of the protected coal seam make the fracture expand and perforate, and the permeability along the bedding direction is increased, and the gas owability is obviously enhanced.Meanwhile, the changed gas pressure gradient will break the original gas adsorption/desorption equilibrium state and transfer partially adsorbed gas into free gas.Under the action of the gas pressure gradient and extraction negative pressure of borehole, gas ows from high-pressure area to low-pressure area.As gas is pumped out, the pressure and content of coal seam gas are gradually reduced.And the outburst danger is eliminated.
According to the law of gas occurrence and migration in the pressure-unloading coal seam, the following hypothesis is given [18][19][20][21][22][23]: (a) the coal seam is a porous elastic continuous medium.(b) Gas is transported to the pores and fractures of coal seams by adsorption and free state.(c) e gas migration process of pressure-unloading coal seam is as follows: rstly, gas is desorbed from the pore surface of the coal seam.Secondly, gas di uses into ssures.irdly, the gas ows from the crack to the borehole (Figure 2).(d) Gas is ideal, neglecting the action of gravity.(e) e tensile stress is positive, and the compressive stress is negative.

Governing Equation of the Hydraulic Field.
e gas in the coal seam is composed of adsorption and free state.And the gas content per unit volume of coal is de ned as follows [19]: where ϕ is the coal seam porosity, ρ g is the gas density (kg/m 3 ), ρ c is the coal seam density (kg/m 3 ), and ρ gs is the gas density under standard condition (kg/m 3 ).
According to the equation of state of ideal gas, the gas density is de ned as follows: where M g is the molar mass of gas (kg/mol), R is the gas mole constant (J/(mol•K)), p is the gas pressure (MPa), and T is the temperature of the coal seam (K).e volume of adsorbed gas can be expressed by the Langmuir equation [19]: where V L is the volume constant of Langmuir (m 3 /kg), p L is the pressure constant of Langmuir (Pa), and p g is the gas pressure (MPa).
According to the conservation of mass, the sum of the change amount of the adsorbed gas and the free gas in the coal seam is equal to the amount of gas owing into the borehole.Considering the e ect of gas and gas slipping, combined with Darcy's law, it can be obtained as follows: where t is the time (s), T s is the temperature under standard condition (K), k is the permeability (m 2 ), μ g is the dynamic viscosity of gas (Pa•s), and b 1 is the Klingberg coe cient (Pa).Substituting Equations ( 1)-(3) into Equation ( 4), the governing equation of the gas seepage eld in coal seam is obtained as follows [20]: 2 Advances in Civil Engineering

Governing Equation of the Mechanical Field.
Coal is an elastic porous continuous medium.And its total strain is the sum of the solid stress, gas pressure, and the strain caused by gas adsorption/desorption.erefore, we can obtain the following equation [21]: where δ ij is the Kronecker symbol; G is the coal shear modulus (Pa), G E/2(1 + υ); K is the bulk modulus of coal (Pa), K E/3(1 − 2υ); E is Young's modulus of coal (Pa); υ is Poisson's ratio of coal; α is the Biot coe cient; and ε s is the gas adsorption strain.rough experimental research, Levine [24] has found that the coal strain caused by gas adsorption accords with the equation of Langmuir form, which can be expressed as follows: where ε max is the maximum adsorption strain and P 50 is the adsorption pressure constant, which refers to the gas pressure when strain is 50% of the maximum adsorption strain (Pa).e geometric and static equilibrium relations of coal are as follows: where F i is the volume force (MPa) and u i is the displacement in the i direction (m).i, j x, y, z.Combining Equations ( 6) and ( 8), the modi ed Navier equation considering pore pressure and adsorption can be obtained: e Mohr-Coulomb criterion is used as the criterion for judging the strength failure of coal rock: where σ 1 is the maximum principle stress (MPa), σ 3 is the minimum principle stress (MPa), φ is the fractural angle ( °), and f c0 is the uniaxial compressive strength (Pa).

Porosity and Permeability of the Unloading Coal
Seam. e permeability of coal seam is closely related to the stress state, mechanical property, and adsorption property of coal.And it can be expressed as follows [25,26]: (11) where k 0 is the initial permeability of coal seam (m 2 ), C p is the compression coe cient of coal seam (Pa −1 ), and c is the adsorption strain coe cient.e subscript "0" is the initial value.
e relation between permeability and porosity is described by the cubic law, which can be expressed as follows [27][28][29]: where φ 0 is the initial porosity of coal seam (m 2 ).Substituting Equation (10) into Equation ( 11), the dynamic evolution equation of coal seam porosity is obtained as follows: e uid-solid coupling model of coal seam gas extraction is obtained by the combination of Equations ( 5), ( 9), (11), and (13).We can use Comsol Multiphysics to carry out

Advances in Civil Engineering
numerical simulation, to study the law of gas extraction in the unloading coal seam of protective layer mining, and to determine the reasonable drilling arrangement scheme.

Numerical Simulation of Rational Holes Arrangement in Gas Extraction
3.1.Geometric Model and Fixed Solution Condition.e Panyi Coal Mine in the Huainan mining area is an outburst mine, and coal and gas outburst accidents have occurred many times, which seriously restricts the safe and efficient production of the mine [9].C13 coal seam is the main mining seam in coal mine; the average thickness of the coal seam is 5 m, the gas pressure is 4.2 MPa, and the permeability is 2.25 × 10 −17 m 2 , which belongs to the high-gas, lowpermeability, and dangerous coal seam.B11 coal seam is located at the bottom of C13 coal seam; the thickness of the coal seam is 2 m, the average distance between the two layers of coal is 64 m, and the risk of B11 coal seam outburst is relatively small.erefore, B11 coal seam is used as the lower protection layer to shield the mining of C13 coal seam.
e tilt width of 1551 (1) working face in B11 coal seam is 220 m, the strike length is 1460 m, and the average seam inclination angle is 5 °. e tilt width of 1551 (3) working face in C13 coal seam is 200 m and the strike length is 1400 m.Limited by the memory performance of the computer, we only simulate the gas extraction in unloading coal seam when the 1551 (1) working face is mined to 800 m, and the geometric model size is 1000 m × 400 m × 95 m (Figure 3).For the physical model, the bottom boundary is fixed and a slip boundary was applied for other sides.e load from the weight of overlying strata with an average pressure of 10 MPa is applied to the top boundary.All the external boundaries for gas are adiabatic and impermeable boundary.Under the initial condition, the coal seam is in the free stress state, and the mechanical parameters of overlying strata are shown in Table 1 [11].
e diameter of gas extraction borehole is 113 mm, and the suction negative pressure is 25 kPa.And other related parameters are shown in Table 2.

Simulation Result Analysis.
Vertical stress distribution of the protected-layer C13 coal seam is shown in Figure 4 after the protective-layer B11 coal seam is mined to 800 m. e stress concentration is formed on both sides of the gob in B11 coal seam, and stress transfer results in stress concentration in the corresponding area of C13 coal seam.e maximum vertical stress is 15 MPa, and the stress concentration coefficient is 1.5.e C13 coal seam area corresponding to the B11 coal seam gob produces pressure relief, the minimum vertical stress is 7 MPa, and the pressureunloading ratio is 0.32.
According to Equation (10), permeability of the coal seam relief zone rises and permeability of the coal seam stress concentration zone is reduced.As shown in Figure 5, the permeability of the protected coal seam is divided into the enhanced zone, the reduced zone, and the original zone.e maximum permeability is 2.85 × 10 −14 m 2 , which is 1266 times as high as 2.25 × 10 −17 m 2 of the initial permeability in coal seam, while the minimum permeability is 3.65 × 10 −18 m 2 , which is 0.162 times the initial permeability in coal seam.
In order to determine the reasonable spacing of boreholes, two types of gas extraction drilling layouts are studied, as shown in Figure 6.
Original arrangement: the boreholes are evenly distributed in the coal seam.Borehole spacing is 20 × 20 m.
Optimized arrangement: the boreholes are nonuniformly distributed in the coal seam.Borehole spacing of the core enhancement zone is 40 × 40 m, that of the general penetration zone is 20 × 20 m, and that of the permeability reduction zone is 5 × 20 m.
From Figure 6, we can see that the layout of boreholes near the open-off cut and working face is relatively dense.If we increase the strike length of the working face, there will be a small difference in the number of boreholes in the two layout types.
By using the fluid-solid coupling model of coal seam gas extraction, the law of gas extraction under the above two types is simulated so as to determine the reasonable way of borehole distribution.Figure 7 is the distribution map of gas pressure in the unloaded coal seam with gas extraction up to 180 days.Obviously, the initial gas pressure of the coal seam is 4.2 MPa; when it is reduced to 0.74 MPa, the coal seam gas extraction rate is more than 30%.erefore, the reduction of gas pressure to 0.74 MPa can be considered to be the standard.
As shown in Figure 7(a), under the uniform distribution mode (original arrangement), the gas pressure is reduced rapidly in the core enhancement zone, and it can be lower than 0.74 MPa after the extraction of 30 days.In the permeability reduction zone, the gas pressure is reduced slowly, the extraction of 180 days is still greater than 0.74 MPa, and the extraction blind area is 4.1 × 10 4 m 2 .However, as shown in Figure 7(b), under the optimized nonuniform hole layout mode (optimized arrangement), the coal seam gas pressure decreases rapidly.After 180 days of extraction, the effective extraction area basically covers the working face of the protected seam, and the blind area is only 0.36 × 10 4 m 2 , which is 0.0878 times that of the original arrangement.
Figure 8 displays the gas pressure of the cutoff line CD in the middle coal seam after 180 days of extraction.Generally, in the middle area of coal seam 150-250 m, gas pressure of original arrangement is almost equal to that of optimized arrangement.In the areas of 100-150 m and 250-300 m, the gas pressure of original arrangement is larger than that of optimized arrangement.It can be found that an extraction blind area at both sides of the researched coal seam will occur if we use the original arrangement.e range of blind area at each side is about 25 m.
Gas extraction capacity of two types of borehole layouts is compared and analyzed, as shown in Figure 9.In the early stage of extraction, gas extraction capacity is increased rapidly by the promotion of high gas pressure gradient.In 30 days of extraction, gas extraction capacity of the two kinds of cloth hole schemes reached 4.77 × 10 6 m 3 and 5.84 × 10 6 m 3 , respectively.With the increase of extraction time, gas pressure of the coal seam decreases, the gas pressure gradient Advances in Civil Engineering

Conclusions
(1) e uid-solid coupling model of gas extraction in the unloading coal seam and the dynamic evolution equation of permeability and porosity of the coal seam are constructed.Numerical simulation of the gas extraction process in the unloading coal seam of the oor roadway is carried out.(2) e stress on both sides of the gob of the protective layer is higher than that of the protected coal seam zone, and the permeability of coal seam decreases.In the middle of the gob, the stress of the coal seam zone is reduced and the permeability is increased.e protected coal seam has the characteristics of permeability zoning, which can be divided into the permeability-enhanced zone, the permeabilityreduced zone, and the original permeability zone.
(3) In the 1551 (3) working face of C13 coal seam, when the uniform borehole distribution method is adopted, the gas pressure in the permeability reduction zone is lower than 0.74 MPa after 180 days of extraction, and the extraction blind area is 4.1 × 10 4 m 2 .e e ective extraction area after 180 days of extraction is basically covered with the protected layer when the nonuniform distribution method is used for the permeability zoning.
Data Availability e data used in the article are from the simulation results by Comsol Multiphysics.And we wrote the equations as Comsol codes independently to conduct these simulations.erefore, you can contact Chaojun Fan (chaojunfan@139.com) if you want to learn more about our work.

Conflicts of Interest
e authors declare that they have no con icts of interest.

Figure 1 :
Figure 1: Gas extraction in the unloaded coal seam after protective layer mining.

Figure 4 :
Figure 4: Vertical stress distribution of the protected coal seam.

Figure 8 :Figure 9 :
Figure 8: Gas pressure of the cuto line CD after 180 days of extraction.

Figure 7 :
Figure 7: Gas pressure in the coal seam after 180 days of extraction.(a) Original arrangement.(b) Optimized arrangement.

Table 1 :
Mechanical properties of overlying strata.

Table 2 :
Parameters of gas extraction in the coal seam.