Coal and Gas Comining in the High-Gas Detong Mine, China: Analytical and Numerical Parameter Optimization of Gas Extraction Boreholes

. The upper corner gas overrun caused by goaf gas gushing in the high-gas mine can be mitigated by roof positioning of long boreholes in the goaf fracture zone, creating an artiﬁcial gas migration channel and achieving safe comining of coal and gas. This study numerically simulated the fracture evolution in the coal rock overlying goaf in workface #2-104 of the Detong Mine, China. Using the O-shaped circle theory and the FLAC 3D commercial software package, the range parameters of goaf fracture zone gas-concentrated areas were obtained by deﬁning the pressure relief coeﬃcient, which provided a basis for the arrangement of holes in the goaf located at 15m from the workface goaf collapse zone to the workface ﬂoor and 15–60m from the structure zone to the workface. According to the obtained goaf parameters, the distribution and migration law in the goaf under the conditions of initial state and extraction with diﬀerent roof borehole parameters were simulated by the FLUENT software. The optimal borehole location was 30 m from the ﬂoor and 25m from the return airway. The ﬁeld experiment with three boreholes arranged in the optimal extraction position provided the extraction concentration in the borehole above 30%, the gas extraction volume of 12 m 3 · min − 1 , and the concentration at the upper corner below 1%, which ensured safe coal and gas comining.


Introduction
While coal mining intensity and depth keep growing annually with increased coal demand, mine productivity is restricted by gas outburst risks.Since mining gas is also a clean energy resource, the realization of coal and gas comining is conducive to improving the utilization rate of resources and reducing greenhouse gas emissions [1].Influenced by the mining and overburden, the fracture zone generated by the bending and separation of the goaf roof becomes a high-concentration gas area, which is the key area to control goaf gas and prevent the overrun of the upper corner gas of the mining face [2].A recent theoretical study comprehensively integrated the O-shaped circle theory, pressure relief mining and pumping, and coal and gas comining [3].Based on the subsidence of a coal mine in Britain, Fang et al. [4] reported that the compaction and stress recovery of rock mass in goaf were mainly related to mining height, buried depth, and other parameters.ey performed the numerical regression analysis of the goaf compaction mechanical properties via the FLAC 3D software.Cheng et al. [5] concluded that the mining fracture of the overburden had the characteristics of ∩-shaped high cap.Qian and Xu [6] revealed a circular fault zone developed above the goaf, constituting a gas accumulation zone and migration channel.Yavuz [7] analyzed the elliptical paraboloid zone's pressure relief gas extraction mechanism and proposed the corresponding coal and methane comining technology.Karacan et al. [8] investigated the effect of the coal mining fracture field on mining stress distribution and coal seam gas pressure.Guo et al [9] coupled the permeability field with the results of the mechanical model and proposed a comprehensive "dynamic" gas reservoir model to simulate the effect of gas gushing and extraction on goaf gas production.Li et al. [10] analyzed the evolution and development law of the goaf gas macroflow channel based on the scale and time effect of the rock mass.Although the above studies provide a theoretical basis for gas accumulation in the goaf and the realization of coal and gas comining, no practical control measures have been proposed.e method proposed in Li et al. [10] can solve the problem of gas accumulation and lacks operability in practical applications.According to the evolution of goaf overburden fracture under mining, the extraction height of roof HLB and its optimal extraction parameters were determined through numerical simulation.

Fracture Development of Overburden and Gas Accumulation in the Upper Corner in Goaf
After mining the coal seam, the coal rock mass accumulates in goaf.According to the evolution of coal rock mass, accumulation state, and gas flow characteristics, the mining space can be subdivided into collapse zone, fracture zone, and bending subsidence zone [11].e coal seam mining reduces the goaf support, triggering numerous secondary fractures in the collapse and fracture zones, which turn into the main migration channel for the gas.With the workface advance, macrocracks developed in the overburden in the middle of goaf are gradually compacted and closed by the ground pressure.e porosity in goaf exhibits nonuniform and continuous distribution patterns, resulting in an O-shaped structure zone [12].e gas in goaf floats up and accumulates in the O-shaped structure zone due to the buoyancy, forming gas concentration areas.e schematic diagram of the three-band distribution has been shown in Figure 1.Due to goaf air leakage, gas flows into the upper corner, resulting in gas overrun [13].Under U-shaped ventilation conditions, the air flow in the goaf flows from the transportation lane of the working face to the cut-off hole of the working face, a small part of which flows to the goaf, and most of the working face flows.
e air leakage in the goaf presents a parabolic shape, which brings out the gas in the deep part of the goaf, where the upper corners of the working surface meet, so that the upper corner of the working surface has a higher gas concentration.At the same time, because the upper corner of the coal mining face is close to the coal wall and the goaf, when the wind flow passes the end of the working face, the wind speed near the coal wall is reduced due to the sudden vertical turn of the roadway, and vortex phenomenon appears in the upper corner of the working face.ere is a phenomenon of air circulation in the vicinity, so that the gas in the goaf and working face is not easily taken away by the wind, and the gas in the upper corner is prone to accumulate.erefore, besides controlling the air leakage in goaf, one has to provide the goaf gas control and upper corner gas overrun prevention.To this end, a continuous artificial negative pressure channel should be set in the gas concentration area in the fracture zone, overlying the goaf to drain the gas in goaf by the pressure gradient to prevent gushing out of high-concentration gas.

Determination of the High-Efficiency
Pumping Area of Goaf Gas

Evolution of Goaf Overburden Simulated by FLAC 3D .
To investigate the evolution of overburden in goaf after coal seam mining and determine the optimal position of roof HLB, this study used the FLAC 3D software package for simulating the plastic stress variation and distribution of overburden after coal seam mining, which provided a basis for studying goaf gas movement law and determining the optimal extraction parameters.Detong Coal Mine is located in Linfen City, Shanxi Province, and the mine field is located in the north of Hedong Coal Field, the coal mining method of the working face has adopted the strike longwall mining method, the average thickness of the coal seam is 5 m, the mining height is 4 m, and the coal seam is buried at a depth of about 390 m.
e "U" type ventilation method is used and the roof lithology is sandy mudstone and mudstone, and the floor lithology is mudstone and sandy mudstone.Based on workface 20104 of the Detong Mine, distance a model with a length of 430 m, a width of 295 m, and a height of 382.3 m was established in this paper.According to the field measurement results, the model's vertical and horizontal stresses (σ z and σ x ) were 20 and 18 MPa, respectively.
e main mechanical parameters of numerical simulation were obtained based on the geological data and geological histogram.
e Mohr-Coulomb criterion was adopted for stresses in the simulation [14].e particular strata of the coal rock mass are described in Table 1, and the established model is shown in Figure 2(a).
e length of a unit square in Figure 2 in the direction of the coal seam advance was 5 m.When the excavation length reached 120 m, the plastic failure area continued to expand along the advancement direction but remained stable and unchanged in the vertical direction.e range of the plastic failure area was 60 m from the floor, and its upper limit could be considered the upper bound of the fracture zone [15], which was consistent with the calculation results.us,   Advances in Civil Engineering it could be concluded that the upper limit of the fracture zone was 60 m. e goaf collapse zone could be determined according to the value and direction of the maximum and minimum horizontal principal stresses σ H and σ h after mining.At σ H > 0 and σ h > 0, it was a bidirectional tensile stress zone mainly distributed in the goaf collapse zone rock mass [15].
In this way, the distributions of the maximum and minimum horizontal principal stresses were obtained through FLAC 3D numerical simulation, as shown in Figure 3.
According to the simulation results, coal seam mining led to stress redistribution in the surrounding rock mass.In the process of stress balancing, a series of mechanical effects such as overburden caving, collapse, separation, and displacement of the overlying coal seam strongly manifested themselves.With the workface advance, the goaf behind the workface was gradually compacted under the overburden weight.As observed, the goaf caving zone was mainly distributed within 15 m from the floor，so the distance from the collapse zone to the floor in further calculations was taken as 15 m.

Distribution of Gas Infiltration Zone in Goaf
Undermining.In the process of gas adsorption-desorption dynamic equilibrium, if the gas is in a saturated adsorption state, the adsorbed gas will be desorbed as the coal body is unloaded; if the gas is in an unsaturated adsorption state, the coal body needs to reach a certain value when the pressure is relieved.e pore gas pressure is reduced to the critical gas desorption pressure and the gas is desorbed.is pressure relief value can be called the effective pressure relief coefficient for gas desorption and drainage.In other words, the lower the saturation of gas adsorption, the higher the degree of pressure relief required for large-scale gas desorption.To observe the coal rock mass evolution after coal seam excavation, the pressure relief coefficient r was defined as where σ Z is the vertical pressure after coal seam excavation and σ 0 is the initial pressure.A previous study [16] revealed that the effective extraction range of gas required a great pressure relief degree of coal rock.Based on the research results, the effective pressure relief coefficient can be preset at 0.8 (i.e., σ z � 0.2σ 0 ), which corresponds to the efficient extraction range.e numerical simulation performed via FLAC 3D software predicted the distribution of workface vertical pressure.e simulation results were sliced according to the distribution range of the collapse and fracture zones, yielding the stress distributions in goaf at different heights, as shown in Figure 4.
According to the pressure relief distribution of the vertical slice, the O-shaped compaction stress rising area was observed in the overlying strata of the workface, and a "vertical stress relief ring" was observed near the O-shaped    Advances in Civil Engineering compaction stress rising area.With an increase in the slice height, the stress in the pressure relief area of the overburden increased gradually and returned to the initial rock stress value; the stress variation rate in the surrounding rock gradually decreased, which was in line with the O-shaped circle theory [3].As observed, a large pressure relief area was observed near the back of the workface.With an increase in goaf height, the range of the pressure relief area gradually reduced, and this area should be avoided in the extraction drilling arrangement to prevent the borehole cut-off, which would deteriorate the gas extraction efficiency.e simulation results revealed that the area with an initial vertical pressure of 20 MPa and pressure relief coefficient of 0.8 was within the range of 0-45 m from the floor height.To avoid the drilling being cut off, the drilling should avoid the area with the pressure relief coefficient >0.8. e optimal height range of the pressure relief channel was 15-45 m.According to the distribution range of pressure relief coefficient and trigonometric function, the effective pressure relief angle of gas efficient extraction range was 76 °in the strike and 85 °in the dip directions.e area with a pressure relief coefficient of 0.2 was a low-permeability gas area, which could be used as a compaction area.According to the simulation results, the distance between the boundary of the compaction area and the workface was 20 m.

Gas Migration Law in Goaf Pressure
Relief Area e migration of gas mixture in goaf was comprehensively affected by air leakage, emission amount of workface gas, and the pressure of coal rock mass in goaf.A further numerical simulation was performed for the case study of workface #2-104 of the Detong Mine, China.

Goaf Gas Migration Control Equation and Assumptions.
For the convenience of calculation, the go extraction area (goaf ) was simplified as follows: (1) e influence of shearer, hydraulic prop, and electromechanical equipment on the fluid in coal mining workface was ignored (2) e mining workface, intake airway, return airway, and goaf areas were assumed to be cuboids, and the gas was treated as an incompressible ideal gas (3) e goaf was assumed to be porous and isotropic [17] (4) e gas flow was regarded as isothermal motion, and the influence of water vapor and other gases in goaf was ignored According to the above assumptions, the goaf can be geometrically modeled and meshed, as shown in Figure 5.
According to the above assumptions, the gas mixture flow in goaf satisfied the mass and momentum equations [18].
e fluid flow in the goaf followed Darcy's law, which could be expressed as where s i is the source term of the i th momentum equation; μ is the viscous resistance coefficient, Pa•s; α is the permeability coefficient of porous media, m 2 ; ] j is the velocity component.e flow model that could deal with the turbulent state of gas flow was k − ε turbulent model, derived by a rigorous statistical technique, and the equations of turbulent kinetic energy k and dissipation rate ε were as follows [19]: where μ t is the turbulent viscosity, μ i is the average velocity, ρ is the density, G k is the turbulent kinetic energy generation term caused by velocity gradient, S ε /S k are the custom source terms, and C 1ε /C 2ε /C 3ε are the empirical coefficients.

Gas Sources in Goaf Analysis and Calculation.
e gas sources of goaf mainly included the gas gushing out of workface (q 1 ), of adjacent layer (q 2 ), and out of goaf (q 3 ).Under the action of buoyancy and air leakage, part of the gas gushed out with airflow, and the other part floated up and was stored into the cracks induced by mining.
ese processes are expressed as follows [20]: where q 1 is the amount of gas gushing from the mining layer, q 2 is the amount of gas gushing from the adjacent layer, q 3 is the amount of gas gushing from the residual coal in the goaf, λ is the gas gushing coefficient, W 0 and W c are the initial and Advances in Civil Engineering residual gas contents, respectively, m is the mining layer thickness, i is the parameter of the i th adjacent coal seam, V 0 is the mining velocity of workface, t is the coal seam exposure time, l is the perimeter of roadway section, l 1 is the workface length, l 2 is the strike length of goaf, and q 0 is the gas gushing intensity of the coal wall.

CFD Numerical Simulation of Mixture Gas Flow in Goaf.
e parameter values of each calculation formula shown in Table 2 were obtained based on the field measured data and the empirical formula.erefore, parameters q 1 , q 1 , and q 3 were taken as 2.97 × 10 −5 , 1.85 × 10 −9 , and 2.6 × 10 −8 kg•(m 3 •s) −1 , respectively.

Geometric Model.
According to the actual dimensions of the site and used scale, the geometric model was set as follows: (1) e coal mining face's length, width, and height dimensions were 185 × 5 × 4 m, and the space type was fluid (2) e scale of the intake and return airway was 10 × 4 × 4 m, and the space type was fluid (3) e scale of goaf was 185 × 200 × 60 m, the fracture and collapse zones were 45 and15 m from the floor, respectively, and the space type was porous media

Boundary Conditions.
e return airway was set as "Out-flow" in the simulations, while the intake airway was assigned a "Velocity-in" value of 2.5 m•s −1 .
Since the goaf was regarded as porous media, its permeability was required for calculations.erefore, the following permeability model based on sigmoid function was used to simulate the permeability parameter evolution [21].
where α 1 � (10/x 0 ), x 0 is the strike compaction boundary (preset at 70 m); α 2 � (10/y 0 ), y 0 is the inclination compaction boundary (preset at 20 m); l y is the workface incline length; k min is the porosity of crushed rock mass after compaction, ranging from 0.1∼0.17[19]; k max is the maximum porosity at the workface, which could be obtained by measurement and is commonly ranged from 0.31∼0.4[22].Substituting equations ( 5) and ( 6) into the Kozeny-Carman function [23], one can obtain the goaf permeability distribution [24]: where z is permeability at maximum porosity, k is porosity, and D p is the average particle size, ranging from 0.014 to 0.016 m.
e permeability was compiled with User Defined Functions (UDF) and Renormalization Group (RNG) k − ε model [25] was adopted to simulate the mixture gas flow in stope considering the turbulence.e distribution of gas in goaf was obtained.
Figure 6 shows the gas concentration distribution in the return airway side section and the goaf plan at different heights.As observed, the gas concentration distributions of the intake and return airways were quite different.e gas concentration of the intake airway was lower, and the lowconcentration area was larger than that of the return airway.In the vertical direction, significant stratification of gas was observed near the workface due to air leakage.e closer it 6 Advances in Civil Engineering was to the mining area, the lower the gas concentration was.e upper and deep parts of the goaf were high-concentration gas accumulation areas.Meanwhile, the vortex was developed near the return airway in the workface under the action of negative pressure and coal wall, and the airflow velocity decreased, resulting in gas accumulation and overrun, jeopardizing production safety.
e specific distribution was as follows.
Along the strike direction of goaf, the gas concentration in the intake airway was nearly zero due to the influence of fresh airflow.However, the amount of fresh air entering the goaf decreased with the deepening of goaf, and the gas concentration grew from 1 to 60%. e gas concentration near the workface on the return air side was no less than 6%, with a large gradient and high concentration.In the vertical direction, the gas concentration in the return airway was high, with up to 30% being concentrated in the range of 0-30 m in the vertical height and 0-20 m behind the workface.
With no countermeasures, the gas in goaf would get into the return airway under air leakage in the workface.At the same time, due to the action of negative pressure and coal wall, the vortex develops near the upper corner, resulting in the gas overrun in the upper corner, jeopardizing production safety in goaf under HLBs.
For the optimal extraction position obtained above, i.e., 15-45 m in the vertical direction, with the respective support of coal wall and overburden pressure, the optimal position in the horizontal direction can be determined as follows [26]:

Advances in Civil Engineering
where α is the collapse angle of fault zone on the return airway side, set at 85 °; β is the dip angle of the coal seam, set at 0 °; L is the advancing distance of the workface, set at 200 m; c is the modification coefficient, set at 1. e compaction area width was preset at 20 m, while a 21.3 m horizontal distance from the return airway was adopted.e borehole type of the roof's HLB diameter extraction was defined as "Velocity-in," the velocity was −5 m•s −1 , and the diameter was 203 mm.
Figure 7 shows the gas concentration distribution of goaf for the HLB (borehole) located 15 m from the floor and 21.3 m from the return airway.As observed, the HLB played a key role in gas drainage in the goaf; the gas in the goaf was discharged by the negative pressure of the HLB and the ventilation airflow.e HLB presence changed the gas flow field at the return air side and the fracture zone.e gas converged to the HLB.e gas concentration at the return air side was lower than that without extraction.
Cloud diagrams of the plane distribution of gas concentration in the goaf, as shown in Figure 8, indicate that the gas concentration in the upper corner decreased to 4% (which was significantly lower than that without extraction), while the low-concentration area at the return air side increased.A negative pressure zone was developed between the gas near the extraction port and the fracture zone, which forced the high-concentration (up to 30%) gas extraction from the fracture zone.
To further study the HLB location effect on the gas distribution in goaf, the gas concentration distribution and extraction volume with the HLBs in different positions were simulated.e preset HLB distances from the floor were 20, 30, and 40 m, respectively, and those from the return airway were 25, 30, and 35, respectively.e results are listed in Table 3.
As observed in Table 3, the optimal extraction parameter corresponded to the borehole located 30 m from the roof and 25 m from the return airway.At this HLB location, the high-concentration gas in the fracture zone could be extracted without affecting the original ventilation airflow.e gas concentration at the upper corner was significantly reduced, and the extraction concentration reached the maximum value.However, simulation results implied that a single borehole could not control the gas concentration at the upper corner within the safe range (>1%).
According to the data in Table 4, the relationship between the HLB parameters and the gas concentration in the drainage pipeline and the gas concentration in the upper corner is obtained through data fitting.
e relationship between the HLB and the gas concentration of the drainage pipeline [27]: e relationship between the HLB and upper corner gas concentration: where W c is the gas concentration in the drainage borehole, W u is the gas concentration in the upper corner, x is the distance from the return air tunnel, and h is the height from the floor.Hence, the number of boreholes should be increased to the number determined by the following equation: where Q j and Q f are the absolute gas emissions of the workface and air row, respectively, K is the disbalance coefficient ranging between 1.2 and 1.7 6 , V is the gas flow rate in the sealed pipe, and C is the gas volume fraction in the extraction pipe (%).e actual values of the respective parameters were substituted into equation (11), yielding the required number of boreholes N � 3.
According to observation Figures 9 and 10, gas in the fracture zone was accumulated in the borehole under the negative pressure of the borehole; a large amount of highconcentration gas was extracted, and as the high-concentration gas in the goaf moved to the deep part of the goaf as a whole, the gas emission amount was reduced.
e gas concentration slice of Z � 2 m in the goaf showed that the gas brought out by the airflow entering the goaf was reduced; the upper corner gas concentration was less than 1%, the production safety was ensured, and the safe comining of coal and gas could be realized.

Analysis of Drilling and Extraction Effect in Roof Strike
5.1.HLB Application Site Design.According to the simulation results, three directional HLBs were drilled in the roof of workface # 2-104 of the Detong Mine; the drilling yard was located at the paralleled return airway (see Table 4 and Figure 11).After the completion of drilling construction, continuous observation was performed on the variation of parameters in the process of the workface advance to 200 m, the variation laws of gas concentration in the upper corner and HLB, and the extraction volume with the advance of workface during mining.

Field Application Analysis.
As observed in Figure 12, the variation of gas concentration in the HLB borehole could be subdivided into three stages.In the first stage, the workface advanced by 20 m in the initial mining; fissures above the goaf were rarely developed, and there was no gas concentration area in the goaf.In this case, the gas concentration and net gas extraction amount were almost zero.In the second stage, the workface advanced by 30-75 m; the goaf roof began to collapse under the overburden pressure, the "three-zone" of goaf began to emerge, and the fractures in the fracture zone formed a gas migration channel.At this stage, the gas concentration and volume in the boreholes increased.In particular, gas concentrations in the three    Advances in Civil Engineering concentration areas in the deep part of the goaf.en, the gas floated up to the fracture zone, and the gas concentration in the borehole had a stable and high level.Additionally, a periodic pressure variation of goaf was observed; the gas concentration of the three boreholes was 30-38%, and the gas extraction volume was maintained at 12 m 3 •min −3 .
In the first and second stages, gas concentration in the upper corner was consistent with the gas extraction volume of the borehole.In a rising trend, the maximum gas concentration in the upper corner reached 0.8%.In the third stage, when the gas concentration and gas extraction volume in the borehole reached a high level, the gas concentration at the upper corner was relatively low, and vice versa.e gas concentration in the upper corner was maintained within the safe range of 0.5-0.6%, and roof HLB parameters obtained by the numerical simulation were quite reliable.
Comparing the actual results on site with the numerical simulation results, the changes in borehole concentration and the upper corner gas concentration were basically consistent with the simulation results.e actual field data is    Advances in Civil Engineering higher than the simulation result.In the following research, the numerical model should be further optimized to avoid errors.

Conclusions
Coal mining induces deformation and damage of goaf overburden, resulting in many secondary fractures and gas concentration areas.e gushing of high-concentration gas seriously restricts mining production safety.erefore, the roof-drilled highly located boreholes (HLBs) are constructed to create an artificial channel to drain the high-concentration gas stored in the fracture zone and improve the production safety with coal and gas comining.eoretical analysis and FLAC 3D numerical simulation were applied to the case study of workface #2-104 of the Detong Mine (China), yielding the goaf overburden stress evolution pattern.rough analysis, it was found that the range of goaf collapse zone was located 15 m from the floor, and the range of fracture zone was 15-60 m from the floor.By defining the pressure relief coefficient, the optimal extraction range of the fracture zone was assessed as 15-45 m; the effective pressure relief angles were 76 °in the strike and 85 °in the dip.FLUENT software was used to provide the numerical simulation of gas migration law in goaf under initial and different extraction parameters.e results showed that the optimal position of roof-drilled HLB was 30 m from the floor and 25 m from the return airway.It was verified that the problem of the upper corner gas overrun could be mitigated by drilling three directional HLBs.Engineering practice revealed that when three roof-drilled HLBs were located 30 m from the workface floor and 25 m from the return airway, the gas concentration in the boreholes exceeded 30%, while the gas concentration in the upper corner did not exceed 0.8%.e accuracy of the simulation results was verified, which provided a reference for solving gas overrun problems in similar coal mines and realizing effective coal and gas comining.

Figure 1 :
Figure 1: Distributions of three fracture development zones of overburden in the goaf.(A) Collapse zone.(B) Fracture zone.(C) Bending subsidence zone.

Figure 2 :Figure 3 :
Figure 2: Numerical calculation model (a) and stress state distribution map (b) of the goaf.

Figure 5 :
Figure 5: Numerical simulation model.(a) Geometric model of working face.(b) Mesh generation of working face.

Figure 8 :
Figure 8: Cloud diagrams of the plane distribution of gas concentration in the goaf: (a) 2 m from the floor; (b) 15 m from the floor.

Figure 10 :Figure 11 :
Figure 10: Plane distribution cloud diagram of gas concentration in goaf during drainage of three boreholes: (a) 2 m from the floor; (b) 15 from the floor.

1 )Figure 12 :
Figure 12: Dynamic change of gas in drainage borehole and upper corner with the movement of working face.

Table 1 :
Physical and mechanical properties of rock strata.

Table 2 :
Calculation formula parameter value.

Table 3 :
Gas concentrations at nine different HLB locations.

Table 4 :
Construction parameters of roof HLBs.