Goaf Gas Control Improvement by Optimizing the Adjacent Roadway Large-Diameter Boreholes

This study introduced gas control technology in goaf using adjacent roadway large-diameter (550 mm) boreholes to control gas accumulation in the upper corner of a fully mechanized working face in high-gas coal seams. The gas control process in the upper corner and gas interception in goaf by large-diameter boreholes was analyzed using the CFD model of the gas flow in goaf. The latter considered the control equation of gas flow, the established permeability model of goaf, and the gas emission law in goaf. Using the 2-105 working face of the Tenghui Coal Mine, Shanxi Province, China, as a case study, the distribution patterns of gas concentration and flow field in the goaf for various extraction flow parameters and different positions of boreholes were numerically simulated. The dependences between various locations, drainage flows, and the gas concentration in the upper corner were determined and fitted by engineering equations. The evolution pattern of the spontaneous combustion zone in the goaf under the drainage conditions was also analyzed. The optimal borehole configuration parameters ensuring the extraction flow rate exceeding 3 m·s−1 and the effective gas control in the upper corner of the working face at a distance of 5 m–15 m behind the working face were identified. The engineering practice proved the feasibility of gas control in the goaf using the adjacent roadway large-diameter borehole. The gas concentrations in the return airflow and the upper corner of the working face were kept below 0.65 and 0.8%, respectively, to ensure production safety and improve the gas utilization efficiency.


Introduction
Despite the global trend of green energy promotion, coal resources remain vital to industrial development, accounting for over half of China's power production. With increased mining intensity and burial depth, the coal seam's gas content also keeps increasing. Various low-gas mines have been converted into high-gas ones, with large gas emissions from the working face, and the gas accumulation in the upper corner often exceeds the limit. e gas emission is hazardous and restricts safe coal production. Meanwhile, mining gas utilization via an extraction system is lucrative from the green energy standpoint [1].
To mitigate the mine gas emission and utilization problems, great efforts have been invested in researching the goaf's gas flow. Majdi et al. studied the calculation method of the caving zone's height in the goaf [2]. Guo et al. proposed the annular fault zone formed above the goaf, constituting the gas accumulation and migration channels [3]. Zhang and Huang derived gas diffusion law in the goaf by studying the gas flow field in the goaf [4]. Zhou et al. applied the twodimensional model to simulate the fracture's gas flow under mining [5]. Karacan and Esterhuizene developed a comprehensive "dynamic" gas reservoir model by coupling the permeability field with the mechanical model to simulate the effects of gas emission and extraction on gas production in the goaf [6]. Wang et al. proposed a method to reduce the gas content and prevent gas outbursts using the large-diameter borehole gas extraction [7]. Yuan et al. introduced a theory and technology of using large-diameter surface boreholes to control the gas in the goaf, which was successfully implemented in the Huainan Mining Area of China [8]. However, theoretical studies [2][3][4][5][6] on the gas distribution and flow evolution in the goaf provided no practical solutions to the gas problems arising in the coal mining process. On the other hand, available recommendations for gas control in the goaf proposed in [7,8] envisaged long construction periods, lacked flexibility and mobility, and failed to solve the problem of gas overrun in the upper corner of the mining face caused by abnormal changes in geological conditions or airflow and abnormal gas emission. Based on the available theoretical achievements, this paper proposes the gas control in the goaf using a large-diameter borehole constructed in the adjacent roadway of the working face. is technology, characterized by high flexibility and productivity, adopts an improved geological drilling machine, producing a largediameter (550 mm) extraction borehole. e drilling parameters can be easily adjusted, according to particular geological conditions, changes in airflow, and abnormal gas emission. In combination with existing gas control measures, the proposed approach mitigates the problem of local gas accumulation in the working face of high-intensity mining quickly, efficiently, and accurately.

Background and Aims of the Proposed Approach.
Under the "U-shaped" ventilation mode, the upper corner of the fully mechanized working face in the high-gas coal seam is close to the coal wall and the goaf side; the airflow velocity is very low, and some areas are in the eddy current state. e gas emitted from the goaf circulates near the upper corner and accumulates in the eddy current area. At the same time, the continuous emission of high-concentration gas in the goaf due to air leakage in the working face provides a high amount of gas in the upper corner, easily reaching an overrun state [9]. At present, the common measures taken to mitigate this problem are as follows. e first one is to change the ventilation mode [10] to prevent the gas from forming the eddy current in the upper corner of the working face. e second one envisages the promoted preextraction of the coal seam [11] and reduces gas emissions from the mining seam. e third one is to control the gas in the goaf via pipe burying or spline extraction in the goaf, high-level drilling or roadway extraction, surface drilling extraction [12], etc. ese measures are characterized by a long construction period and high costs. Especially in the case of sudden local gas accumulation, a long time is needed for adjustment. erefore, it is necessary to use rapid, efficient, flexible, and low-cost techniques to reduce local gas accumulation.

Advantages of the Adjacent Roadway Large-Diameter
Extraction.
e large-diameter gas extraction boreholes of adjacent roadway construction at the air-return side of the working face shown in Figure 1 can change the eddy current formed at the upper corner of the working face into the original flow field of the goaf, introduce new disturbance, extract the goaf gas in the form of large flow and low negative pressure, and intercept the gas from the goaf due to the air leakage of the working face to solve the problem of local gas overrun.
is technology features mechanized operation, simple process, and low engineering complexity, which can effectively solve the problem of connection and balance of mine "mining and extraction" in terms of time. e borehole spacing and the number of boreholes can be adjusted at any time with the working face advance. With strong flexibility, it can exert an effective interception effect on the emitted goaf gas. e key issue of this technology application is to determine the optimal drilling configuration parameters. e computational fluid dynamics (CFD) numerical simulation can be used to determine the gas concentration field and flow field distribution pattern of the goaf at various large-diameter boreholes' parameters. en, the drilling configuration design can be optimized to get the best extraction parameters. e technology of large-diameter drilling in the adjacent roadway to control the goaf is illustrated in Figure 1. A gas drainage process is schematically presented in Figure 1(a), where the extracted gas (i.e., coalbed methane or CBM) is transported via pipelines, with an improved utilization rate of energy resources and environmental protection. Figures 1(b)-1(d) show the specific construction equipment and technology of gas drainage.

eories of Goaf Gas Migration.
Under the mining and the overlying rock pressure, the size of voids in the goaf exhibits certain trends in the vertical and strike directions. In the vertical direction, the caving zone, fracture zone, and bending subsidence zone are formed [13]. In the strike direction, with different accumulation states and pressure on the coal and rock mass, three areas are formed: the natural accumulation area, the caving load area, and the recompaction area [14]. During the coal seam mining process, the goaf loses its supporting force, and a considerable number of secondary fissures are formed in the caving zone and the fracture zone, which become the main gas migration channels. e coal and rock mass in the caving zone is filled with a considerable amount of broken coal and rock mass, resulting in larger porosity. e fracture zone has a smaller porosity due to the support of the caving zone. e goaf's porosity shows nonuniform continuous distribution characteristics [15] and presents an O-shaped distribution, as shown in Figure 2. e height of the fracture and caving zones in the goaf can be obtained through physical analog modeling, numerical simulation, field observations, and theoretical calculations [16]. Determination of these values is conducive to analyzing the physical properties and gas flow characteristics of the goaf. It can provide basic support for the establishment of CFD numerical analysis. According to the law of movement and failure of the overlying coal and rock mass in the goaf, the caving zone and the fracture zone in the goaf can be calculated as follows: where H k is the caving zone height, H is the mining height, W is the main roof subsidence, K pb is the crushing expansion coefficient of the caving rock, and α is the dip angle of the coal seam, while a, b, c are the coefficients related to the roof properties, which can be determined according to the empirical formula in [5]. According to the "O"-ring theory, the goaf's porosity distribution gradually increases from the goaf's edge to the central compaction area. Based on the mathematical model of a continuous surface of the subsidence of the stope overburden, the porosity distribution model of the caving zone in the goaf has been theoretically deduced [17] in the following form: where n is the porosity of any point, c is the volume weight of the caving rock, α is the dip angle of the coal seam, l y is the dip width of goaf, β is the regression coefficient, h d is the thickness of the immediate roof, and l is the length of the broken rock block of the main roof. For the convenience of calculation, the following assumptions are made for the goaf in this study [18]: (1) e shearer and all kinds of equipment of the hydraulic prop in the fully mechanized top coal caving face are ignored; only the air leakage from the working face into the goaf is calculated for the goaf gas, and the return air roadway and extraction drilling on gas emission and concentration distribution are studied.  (2) e working face, the inlet, and return air roadways are assumed to be cuboids filled with an incompressible ideal gas. (3) e goaf is assumed to be porous and isotropic [19]. erefore, the goaf's gas mixture flow satisfies the mass equation and the momentum equations [20]. e mass equation has the following form: e momentum equations can be expressed as follows: e turbulent kinetic energy equation is e transport equation of the turbulent dissipation rate is given below: where μ t is the turbulent viscosity, μ i is the time-averaged velocity, ρ is the density, G k is the generic term of turbulent kinetic energy k caused by the velocity gradient, S ε and S k are user-defined source terms, and C 1ε , C 2ε , and C 3ε are empirical constants.
If the goaf is regarded as the porous medium, it possesses the additional (viscous and inertial) resistance to the gas mixture in the goaf. By adding the loss term in the momentum equation, including viscous resistance and inertial resistance, the following equation is obtained: where k is the goaf permeability (the reciprocal of viscous resistance), μ is the dynamic viscosity, v i is the velocity component in direction i, and C 2 is the inertia resistance coefficient. Due to the small air leakage in the goaf, the inertia resistance is usually neglected. As shown in equation (7), permeability is the key parameter for simulating the goaf gas migration law. According to Darcy's law [21], the permeability model can be represented by the Kozeny-Carman equation [22,23] in the following form: where D p is the average particle size ranging from 0.014 to 0.016 m [24] and n is the porosity of the porous medium.
Equation (8) constitutes the governing equation of the gas flow in the goaf. For the particular boundary and initial conditions, the gas flow pattern can be determined.

Gas Source Analysis.
e gas emission in a fully mechanized working face can be calculated using the different-source prediction method, with the coal seam gas content used as the basic parameter [25]. It follows the principle of establishing an empirical formula based on years of detection data according to the gas content of coal seam, the gas gushing pattern of gas gushing source, the occurrence conditions, and mining technical conditions of coal seam [26]. e gas emission of each part is calculated to achieve the gas emission calculation.
Gas emission from a fully mechanized working face includes gas emissions from the mining seam, from the coal remaining in the goaf, and from the adjacent seam, as shown in Figure 3. e calculation formula of gas emission in the mining seam Q 1 is as follows: and that in the adjacent seam Q 2 is e calculation formula of gas emission from the coal remaining in the goaf Q 3 is (11) where λ is the gas emission coefficient, W 0 is the original gas content, W c is the residual gas content, m is the mining seam thickness, i is the parameter of the i-th adjacent seam, l is the perimeter of the roadway fault face, v 0 is the advancing speed of the working face, t is the exposure time of the coal seam (which is taken as six months), l 1 is the strike length of the working face, l 2 is the strike length of the goaf, and q 0 is the gas emission intensity of the coal wall. e latter parameter can refer to the calculation of coal wall gas emission in "Prediction method of mine gas emission" (AQ1018-2006): where V f is the volatile content in the coal, taken as 18.79%.

Basic Parameters of the Working Face.
e Tenghui Coal Mine is located in Xiahua Township, Hejin City, Shanxi Province, China, with the mine field located on the southern edge of the Hedong Coalfield. e strike longwall mining method is adopted for the working face. e average thickness of the coal seam is 5.2 m, the widths of the inlet and return air roadway of the working face are 4 m, the height is 3 m, the length of the working face is 185 m, the mining height is 3 m, the mining width is 5 m, the roof periodic weighting step is 15 m, and the buried depth of the coal seam is 368-480 m. e U-shaped ventilation mode is adopted. e lithology of the roof is sandy mudstone and mudstone, and the lithology of the floor is mudstone and sandy mudstone, which belongs to a medium-hard roof.
According to the previous basic theoretical analysis, the geological conditions of the Hedong Coalfield, the actual onsite test data of the Tenghui Coal Mine working face, the experimentally measured data, and the empirical parameters' values, the calculation formula parameters were obtained and are listed in Table 1. e physical parameters of the numerical model are summarized in Table 2.

Establishment of a Geometric Model of a Fully
Mechanized Working Face. According to the actual situation in the coal mining face, the size and boundary conditions were set as follows: (1) e size of inlet and return roadways was 10 × 4 × 3 m 3 , the size of the working face was 185 × 5 × 3 m 3 , and space was filled with fluid. (2) e goaf size was 185 × 200 × 57 m 3 , in which the fracture zone was 40 m, the caving zone was 17 m, and space was a porous medium. (3) Except for the inlet and return air roadways of the working face, all the outer boundaries of the extraction drilling were set as the "wall," and the interface between the working face and the goaf was set as the "interior." It was assumed that the stope had no material exchange or conditional reaction with the outside environment. e established geometric model is shown in Figure 4(a), and red arrows represent the model's corresponding coordinate system in Figure 2.

Mesh Generation and Model Realization. (1) Mesh
Generation. e established physical model's calculation domain was meshed using a tetrahedral mesh, as shown in Figure 4(b). e inlet and return air roadways' mesh sizes were both 0.5 m, and those of the caving and fracture zones in the goaf were 2 m each. When the large-diameter borehole was connected to the goaf, spherical densification was carried out around the goaf, and the densified mesh size was 0.1 m, as shown in Figure 4(c).
(2) Solver Setting. e gas concentration and flow field distribution in the goaf were numerically simulated at the borehole distances of 0, 5, 10, 15, and 20 m from the working face under the extraction-free and extraction conditions. It was assumed that the flow rate of the extraction borehole remained unchanged. e simulation conditions were as follows: the inlet velocity of 2.5 m·s − 1 was set at the inlet of the inlet roadway, the "outlet flow" was set at the outlet of the return roadway, the inlet velocity of the large-diameter extraction borehole was − 5 m·s − 1 , the borehole diameter was 550 mm, and the vertical distance from the floor was 2 m.
When the goaf was treated as the porous medium, permeability was a parameter to be determined. In the fracture zone, the permeability variation range is usually between 10 − 6 and 10 − 11 m 2 [27]. In this study, the porosity of the fracture zone was taken as a fixed value of 0.03 [28]. Substitution of this value into equation (8) yields the permeability range between 10 − 6 and 10 − 11 . e porosity of the caving zone of the actual working face was substituted into equation (3) to get the porosity distribution, which had a "shovel" shape.
is implies large porosity values of the shallow part and two roadway sides and small ones of the middle and inner parts, as shown in Figure 4(d). e porosity value was substituted into equation (8) to obtain the permeability parameters.
A user-defined function (UDF) was employed during the CFD numerical simulation to compile the caving zone's permeability.
e "simple" algorithm was selected in the pressure-velocity coupling algorithm, and the renormalization group (RNG) K-epsilon (k-ε) turbulence model [29] was adopted to simulate the gas flow mixture in the stope to improve the simulation accuracy.

Analysis of Gas Migration Law in Goaf under Drainage
Conditions. A coordinate system was established according to Figure 4, with the depth direction, the inclination direction, and the vertical direction of the goaf taken as the X-,

Advances in Civil Engineering
Y-, and Z-axes, respectively. e range of the X-axis [− 2 m, 0] was regarded as the upper corner of the working surface. According to the coal mine safety regulations, 0.8% of the gas concentration in the upper corner is the critical value, while that below 0.8% is considered to be safe, implying the draining effectiveness. e numerical simulation results of gas migration under the conditions of Z � 2 m and no extraction and large-diameter borehole extraction in the goaf area are illustrated in Figure 5. It can be observed from Figure 5(a) that due to air leakage from the working face, the airflow enters the goaf to bring the high-concentration gas in the goaf to the return air side and finally dilutes into the return airflow. Due to the coal wall effect, the air current forms a vortex near the upper corner of the working surface, causing local overrun, with the gas concentration exceeding 6%. e deeper the goaf, the higher the gas concentration, up to 70%, making the goaf a gas reservoir. erefore, extraction measures must be taken to restrain the upper corner gas from exceeding the limit.
In this paper, simulation schemes with different flow rates and locations listed in Table 3 were designed to explore the impact of large-diameter boreholes on gas migration in goaf. e gas distribution law was obtained by comparison. According to the actual working conditions of the on-site extraction pump, the maximum extraction flow rate of the borehole was 100 m 3 ·min − 1 , indicating that the maximum flow rate of the mixed gas in the borehole was 7 m·s − 1 .
Other boundary conditions and the calculation process of the solver were the same as before. Besides, the gas concentration values in the large-diameter borehole and the upper corner were taken as dependent variables. e numerical simulation results are presented in Figure 6.
Best fitting of the numerical simulation results was performed, yielding the following formula: where c denotes the gas concentration in the upper corner, v is the extraction rate, and y is the distance from the working face, while R 2 is the correlation coefficient, which the value is 0.958 that to indicate a close correlation. As seen from the fitting formula, when the drilling position is fixed, the higher the extraction rate, the lower the gas concentration in the upper corner. When the extraction rate is constant, the farther the drilling position from the working surface, the higher the gas concentration in the upper corner.
As illustrated in Figure 6(a), with an increase in the extraction rate, the concentration in the gas pipeline gradually grows. e farther the borehole from the working    surface, the higher the gas concentration in the borehole. e gas concentration was in the range from 20 to 60%. Simultaneously, the gas concentration gradient gradually dropped.
In Figure 6(b), the blue area represents the zone where the gas concentration in the upper corner was below 0.8%. It can be observed that the gas concentration in the upper corner gradually decreased as the extraction rate increased. When the borehole extraction rate was lower than 3 m·s − 1 , the upper corner's gas concentration decreased faster than that when the extraction rate exceeded 3 m·s − 1 . When the borehole was located 10 m behind the working face, the gas concentration in the upper corner was the lowest, and the extraction effect was the best. e farther the borehole from the working surface, the weaker the ability to control the upper corner. When the borehole was 20 m away from the working surface, the extraction rate exceeded 6 m·s − 1 , ensuring that the upper corner gas concentration was less than 0.8%. e analysis of Figures 5 and 6 and fitting formula (14) revealed that the negative pressure of the borehole changed the distribution of the gas flow field in the goaf. Part of the gas that was directly pumped out contributed to intercepting the gas gushing out of the goaf. Simultaneously, it disturbed the eddy currents formed in the upper corners of the working surface and prevented the local gas accumulation under certain conditions of extraction flow and drilling positions. However, the negative pressure also attracted high-concentration gas in the goaf to move closer to the working face, leading to increased air leakage from the   Advances in Civil Engineering working face. is increased the gas concentration gradient on the return air side, leading to a higher risk of spontaneous combustion in the goaf.

Influence of Large-Diameter Borehole Extraction on Spontaneous Combustion Zone in Goaf.
When the gas in the goaf is subjected to high-flow extraction, the air leakage state in the goaf will change, varying the width of the oxidation zone and the air leakage in the goaf. In this paper, the oxygen concentration in the goaf was used as a zoning criterion. e area with the oxygen concentration exceeding 18% in the goaf was identified as a heat dissipation zone, while that with the oxygen concentration in the goaf between 8 and 18% was regarded as an oxidation heating zone. e oxygen concentration in the goaf area in the breath zone was less than 8%. To study the distribution law of the spontaneous combustion zone in the goaf, the distribution variation of the spontaneous combustion zone was compared when the borehole was located at different positions under the condition that the extraction rate of the large-diameter borehole was 6 m·s − 1 . In addition, the distribution variation of the spontaneous combustion zones in the goaf under the condition of different extraction rates of the large-diameter borehole was compared when the borehole was located at 10 m.
At Z � 2, there was no extraction; when the borehole was located at 10 m, the extraction rate was 6 m·s − 1 . When the borehole was located at 20 m, the oxygen concentration distribution cloud diagram was obtained, as shown in Figure 7.
According to the numerical simulation results on oxygen distribution in Tables 4 and 5 and the cloud diagram in Figure 7, larger drilling depths in the goaf corresponded to larger extraction rate, wider spontaneous combustion zone of the goaf, and higher risk of spontaneous combustion. e disturbance effect of the extraction of large-diameter boreholes caused the spontaneous combustion zone of the goaf to widen and get deeper in the goaf. erefore, the risk of spontaneous combustion in the goaf increased. e simulation demonstrated that when the borehole was located at a position of 20 m, the mixed gas extraction rate of the borehole was 6 m·s − 1 , the maximum width of the spontaneous combustion zone in the goaf was 30 m, and the shortest spontaneous combustion period of the coal seam of the Tenghui Coal Mine was 35 d. is implied that the safe advance rate of the working face was no less than 0.86 m per day.

Field Test Results
According to the numerical simulation results, when the flow rate of the extraction borehole was no less than 3 m·s − 1 , the borehole located in the range of 5 m-15 m behind the working face inhibited the gas accumulation in the upper corner of the working face. At the flow rate of the extraction borehole no less than 7 m·s − 1 , the borehole located [5 m, 20 m] from the working face could inhibit the gas accumulation in the upper corner of the working face. However, the flow rate of 7 m·s − 1 was the maximum limit value that could be assigned to a large-diameter drilling pipeline of the extraction system. e 2-105 working face was located southwest of the mine field. As the second working face of the first mining area of the mine, it was arranged in the west wing of the first mining area, east-west direction. On the east wing of the working face was the return air roadway, belt transportation lane, and auxiliary transportation lane of the first mining area. e south and west wings were the boundaries of the coal mine field, and the north wing was the 2-104 working face. Combined with the simulation results and the actual conditions of the Tenghui Mine's 2-105 working face, the extraction borehole was designed. e construction site was located in the 2-105 adjacent roadway, which was the preparation roadway for the next working face. e dimensions of the adjacent roadway of the working face were 3 m × 3 m × 4 m, and the roadway length was 800 m. e engineering design implied a borehole interval (spacing) of 10 m and a borehole distance of 2 m from the floor, considering the production operability and safety. When the extraction borehole entered the goaf and was located 5 m from the working face, the control valve was opened for gas extraction, while the borehole located 15 m from the working face was closed. e control valves of subsequent boreholes were operated in turn, and the flow rate of the extraction borehole was no less than 3 m·s − 1 . e boreholes were numbered in increasing order with the distance from the working face, as shown in Figure 8.
With the advance of the working face, when the 1# borehole connected to the goaf entered the goaf and was 5 m from the working face, the 1# borehole control valve was opened for gas extraction. When the 2# borehole entered the goaf and was 5 m from the working face, the 2# borehole control valve was opened for gas extraction, and the 1# borehole control valve was closed simultaneously. e alternate switching of boreholes was similar to the subsequent drilling.
According to the performed investigation of three boreholes, the gas concentration in the upper corner of the working face and the return air roadway was within the controllable range, the gas concentrations in the upper corner and the return air roadway ranged between 0.5 and 0.8% and between 0.45 and 0.61%, respectively. As the working face advanced, with an increase of the distance between the working face and the large-diameter borehole, the evolution curves of gas concentration in the upper corner of the working face and the return air roadway were "V"-shaped. When the borehole was 5 m away from the working face, the gas concentrations in the upper corner and the return air roadway reached their maximal values of 7 and 0.6%, respectively. When the borehole was about 12 m away from the working face, the borehole had the strongest control over the upper corner, the lowest gas concentration in the upper corner was about 0.55%, and the gas concentration in the return air roadway was about 0.45%. As the working face continued to advance, the gas concentration in the upper corner and the return air roadway gradually increased and changed periodically with the opening of the borehole. e variation trend of the gas concentration in the upper corner was consistent with the trend of the simulated data-fitting equation.
e gas concentration in the large-diameter borehole gradually increased. e deeper the borehole into the goaf, the higher the gas concentration in the borehole. When the borehole was located in the shallow part of the goaf, i.e., about 5 m behind the working face, the negative pressure around the borehole was relatively low, and the gas injection in the goaf was relatively small. e gas concentrations in the pipeline and the upper corner were about 15 and 0.6%, respectively. When the borehole was located 12 m behind the coal mining face, the gas concentration reached the highest value of about 37%. At this time, the gas concentration in the upper corner was maintained at about 0.55%, and then the gas flow in the pipeline tended to be stable. Meanwhile, the gas concentration in the borehole changed periodically with the conversion of the borehole. e comparative analysis of Figures 8(a) and 8(b) revealed that there was a certain relationship between gas concentrations in the extraction pipeline, the upper corner, and the return air roadway, i.e., when the gas concentration in the extraction pipeline increased, the gas concentration in the upper corner and the return air roadway tended to decrease. e extraction flow and borehole spacing determined by the numerical simulation were verified in practice And, when the drilling hole was 10m-12m behind the working face, the extraction effect had been the best effect. During the actual advancement of the working face, the corresponding relationships between the advancing distance and the gas concentration in the upper corner, gas extraction volume, and gas concentration in the return air roadway were observed, and the operation model's specific parameters in this particular mine were derived. e main innovation of this paper is that it substantiated and successfully implemented a large-diameter (550 mm) borehole in the roadway adjacent to the working face to control the gas accumulation. To the best of the authors' knowledge, boreholes with a diameter not exceeding 300 mm have been applied yet for this purpose. According to the relevant definition [30], the underground large-diameter hole generally refers to the near horizontal hole with the final    hole diameter of more than 200 mm, which is used for pipeline laying (including ventilation pipe, gas drainage pipe, and drainage pipe), through tunnel rescue (for provision, drainage, escape, etc.) and gas drainage. us, the adopted borehole diameter of 550 mm corresponds to largediameter definition, while larger diameter values are envisaged to be tested in the follow-up study.

Conclusions
is study analyzed the gas accumulation patterns caused by the development of fissures in the goaf of the fully mechanized working face, gas storage, and the formation of eddy current in the upper corner. It introduced a technical scheme of gas control in the goaf, which used the adjacent roadway large-diameter boreholes, analyzed gas migration law in the goaf, and determined the optimal configuration parameters.
(1) Within the CFD framework, firstly, a borehole gas flow simulation scheme was designed under the conditions of different borehole locations and extraction rates. en, the optimal position of largediameter drilling was determined through simulation. (2) Next, the influence of large-diameter borehole drainage on the spontaneous combustion zone of the goaf was analyzed. erefore, the optimal drilling configuration parameters based on the upper corner gas concentration in the fully mechanized mining face of the Tenghui Mine were determined: the velocity of the borehole was no less than 3 m·s − 1 , the effective range of the borehole placement was 5 m-15 m behind the working face, and the advance rate of the working face was no less than 0.86 m per day.
(3) e engineering parameters were designed according to the numerical simulation results. By investigating the three boreholes' extraction parameters, it was determined that this process could effectively mitigate the local gas overrun in the target mine. e gas concentration in the upper corner was controlled below 0.8%. e gas concentration in the return air roadway was controlled below 0.65%, which was in line with the coal mine safety regulations' relevant provisions. (4) A part of the emitted gas entered the extraction system, which ensured the comprehensive utilization of gas as a clean energy source and provided a solution for similar mines to mitigate gas emission problems.

Data Availability
All data and models used to support the findings of this study are available from the corresponding author upon request.

12
Advances in Civil Engineering