Study on optimization of layout parameters of high-level boreholes in Pingdingshan coal mine

To reduce the gas disaster of high gas coal seam and improve the efficiency of gas extraction by high drilling, the layout parameters of drilling holes in Pingdingshan coal mine are optimized. Based on the analysis and calculation of the "three zones" of the movement towards the overly strata of No.10 coal in Pingdingshan coal mine, the height of caving zone and fissure zone in 24,130 working face are 10.06–14.46 m and 38.75–49.95 m respectively. The elevation angle, azimuth angle and the length of high-level boreholes are studied and analyzed by COMSOL numerical simulation software. The simulation results show that the optimum layout parameters of high-level boreholes are as follows: The elevation angle of borehole should be controlled at 9°–12°, the azimuth angle should be 30°–45°, and the length of borehole should be 150 m. Then the optimum layout parameters of high-level boreholes are determined for engineering application of 24,130 working face. Borehole data onto actual mine show that the optimum layout parameters of high-level boreholes were elevation angle between 8°and 11°, azimuth angle between 30° and 42°, and length of boreholes between 145 and 155 m. The simulation results are basically consistent with the measured data. The maximum gas concentration in working face, upper corner and return air roadway is stably controlled below 1%. The safe mining of 24,130 working face is ensured, which provided a certain reference value of gas control in the goaf of Pingdingshan mine and adjacent mines.

measures to improve gas drainage efficiency by examining the factors influencing gas migration in goaf areas, without further discussing the impact of of the difference of the Angle and location of the actual borehole on gas drainage from fractured zones in the goaf.
Therefore, this study focuses on the actual situation of the 10# coal seam in Pingdingshan coal mine, using numerical simulation software to determine the range of high-level borehole arrangement based on comprehensive parameters including elevation angle, azimuth angle, and borehole length.Field validation was conducted in the roadway of the 24,130 working face in Pingdingshan Coal Mine, and the measured data from the field were found to be consistent with the simulation results.Finally, the optimal layout parameters for high-level boreholes in the 24,130 working face of Pingdingshan coal mine were determined.

Prediction of the height of "three zones" of overlying strata fractures
Using the empirical formula method, the calculation of the height of the caving zone for the 24,130 fullymechanized mining working face is performed 26 .Based on the geological survey data of Pingdingshan coal mine, the immediate roof of the 10# coal seam mainly consists of mudstone and sandy mudstone, with some areas comprising siltstone and fine to medium-grained sandstone.According to the measured results of rock mechanics parameters, the average compressive strength of the mudstone is 27.5 MPa, and the average coal thickness is 5.5 m.Therefore, the calculation is performed based on the assumption of a medium-hard rock layer: where M is the cumulative thickness of the mined coal seam, m (it is taken as 5.5 m); H m is the Height of the caving zone, m.
To calculate the height of the fractured zone within the 24,130 gas drainage working face, the immediate roof of the 10# coal seam is typically composed of K2 limestone.Therefore, it can be treated as a hard rock layer for calculation purposes.

Construction of numerical calculation model
To determine the impact of parameters such as elevation angle, azimuth angle (angle with respect to the tunnel wall), and borehole length on gas drainage effectiveness, this study utilizes the COMSOL numerical simulation software based on the specific conditions of the site.A coupled model is constructed considering the coal seam occurrence conditions and extraction parameters.The variation in gas concentration is used as the indicator for optimal borehole arrangement.By studying the various influencing factors, the optimal layout parameters for high-level boreholes are determined.
Based on the actual parameters of the high-level boreholes in the 24,130 working face of Pingdingshan coal mine, this study selects the following parameters: inclined length of the working face is 186 m, goaf length is 200 m, dimensions of the intake and return airways are 3 × 4 m, height of the working face is 3 m, width is 5 m, height of the caving zone is 15 m, height of the fractured zone is 45 m, length of unmined solid coal behind the working face is 80 m, and coal thickness in the solid coal region is 11.5 m.The overlying strata above the coal seam in the unmined area are considered as rock layers.The geometric model is shown in Fig. 1.The meshing of the model is shown in Fig. 2. The elevation angle of the borehole is denoted as a 1 , the angle between the borehole and the coal seam is denoted as a 2 , and the length of the borehole is L. Due to the significant computational workload of three-dimensional models, this study only simulates individual boreholes and does not consider the interaction between boreholes.The main focus is to optimize the borehole elevation angle, azimuth angle, and length as key technical indicators through simulation. (

Numerical model parameters and boundary conditions
The numerical simulation in this paper is based on the relevant physical parameters of Pingdingshan coal mine.The inlet speed is controlled at 2.5 m/s, the air volume is 1800 m 3 /min, the return air outlet is controlled by the pressure difference between ends 110 Pa, the high level borehole is set as the mass outflow boundary, the extraction flow rate is q ch4 , the solid coal is set according to the original gas content, and the goaf is set according to the proportion of residual gas.The height of the fracture zone is 45 m, and the height of the collapse zone is 15 m.The relevant parameters are shown in Table 1.

Analysis of simulation results of influence of high level borehole parameters on extraction effect
The effect of elevation angle of borehole The elevation angle of a single borehole is set as a 1 = 6°,9°,12°, and 15°, borehole length L = 150 m, angle between borehole and coal wall a 1 = 30°, extraction flow q ch4 = 3 m 3 /min, and extraction time is 240 h. Figure 3 shows the change of gas concentration at 1#(93,50,30), 2#(93,100,30) and 3#(93,150,30) monitoring points in the goaf.Figure 3 shows that the gas concentration at 1# starts to decrease significantly after 12 h of extraction, and gas extraction is completed after 50 h.Among the tested elevation angles, the lowest residual gas concentration is observed at an elevation angle of a 1 = 12°, which is nearly equivalent to that at a 1 = 9°.The highest residual gas concentration is observed at a 1 = 6°, followed by a 1 = 15°.The gas concentration at 2# begins to decline notably after 40 h of extraction, and gas extraction is completed after 140 h.Among the tested elevation angles, the lowest residual gas concentration is observed at a 1 = 9°, while the highest concentration is observed at a 1 = 6°, and a 1 = 15° yields the second-highest concentration.The gas concentration at 3# starts to decrease significantly after 100 h of extraction, and there is still a certain concentration of residual gas after 240 h.Among the tested elevation angles, the lowest residual gas concentration is observed at a 1 = 9°, followed by a 1 = 12°, while the highest concentration is observed at a 1 = 6°.
Figure 4 shows that the gas concentration is relatively high at a 1 = 6°, 15°, highest when a 1 = 6°, and secondhighest when a 1 = 15°.The concentration is relatively lower at a 1 = 9° or 12°, lowest when a 1 = 9°, and secondlowest when a 1 = 12°.This suggests that the optimal range for the elevation angle a 1 is greater than 6° and should be maintained between 9°and 12° for better gas extraction effectiveness.5 shows the change of gas concentration at 3 monitoring points of 1#-3# in the goaf.Figure 5 shows that the gas concentration at 1# began to decrease significantly after 6 h of extraction, and the gas extraction is completed after 80 h, in which the residual gas concentration is the lowest when the azimuth Angle a 2 = 30° and 45°, the highest when a 2 = 20°, and the second when a 2 = 60°.The gas concentration at 2# began to decrease significantly 42 h after extraction, and is completed 180 h after extraction.The lowest concentration of residual gas is found at a 2 = 30°, the highest concentration is found at a 2 = 20°, followed by a 2 = 45 and 60°.After 100 h of extraction, the gas concentration at 3# began to decrease significantly, and after 240 h of extraction, there is still a certain concentration of gas.The highest concentration of residual gas  www.nature.com/scientificreports/ is found when a 2 = 20°, the lowest concentration is found when a 2 = 30°, followed by a 2 = 45 and 60°, indicating that the azimuth between 30° and 45° is the best effect.Figure 6 displays the gas concentration distribution in 240 h extraction under different a 2. The figure shows that the residual gas content is lower in the intake air side and higher in the return air side.Additionally, when a 2 is between 30° and 45°, the residual gas content throughout the goaf is relatively low, indicating a more significant gas extraction effect and better gas extraction efficiency.The effect of borehole length For each individual borehole, the lengths are set as L = 120 m, 140 m, 160 m, and 180 m, with an elevation angle of a 1 = 9°, azimuth angle of a 2 = 30°, gas extraction flow rate of q ch4 = 3m 3 /min, and a total extraction time of 240 h. Figure 7 illustrate the gas concentration variations at three monitoring points.Figure 7 demonstrates that the residual gas content in the shallow part of the goaf has a low correlation with the length of the borehole, as varying the borehole length does not result in significant differences in residual gas content.At 2#, the highest residual gas content in the goaf is observed when L = 180 m, while the residual gas content gradually decreases for L = 140 m, 160 m, and 120 m.In the deeper part of the goaf, there is minimal variation in the trend of residual gas content for L = 120 m, 140 m, 180 m, and 120 m.In conclusion, under unchanged conditions, the borehole should not be too short and should be greater than 120 m to ensure a significant influence range in the horizontal direction.Controlling the length within the range of 140-160 m leads to lower residual gas content in the goaf and relatively better gas extraction effectiveness.
Figure 8 presents the gas concentration distribution in the goaf after 240 h of extraction under different borehole lengths.The figure illustrates that the larger the borehole length, the lower the residual gas content in the goaf.However, when the length reaches 180 m, the decrease in residual gas content in the goaf is not significant.

Test site profile
The 24,130 working face in Pingdingshan coal mine is located at an elevation of + 860 to + 936 m.It has a strike length of 1113 m and a dip length of 220 m, covering an area of 244860 m 2 .The mine has an absolute gas emission rate of 57.26 m 3 /min and a relative gas emission rate of 26.8 m 3 /min.The coal seam in this working face has a dip angle ranging from 0° to 5°, with an average thickness of 3.5 m.The coal seam has a bulk density of 1.46 t/m 3 and a permeability of 0.42 m 2 /(MPa 2 d), indicating that it belongs to a high-gas-emitting coal seam.

Test drilling site and drilling construction design
According to the theory of mining ellipsoid zone in overburden rock, when the fully mechanized caving face is fully influenced by the initial and periodic mining processes, the middle of the caving face is compacted, and a certain width mining fracture zone is formed on both sides, and the mining fracture field shows a parabolic distribution in the inclined direction.The layout of high directional long boreholes should be considered from the horizontal and vertical distance, and should be arranged in the fracture zone of high gas enrichment area with dense distribution of horizontal and vertical fractures, so as to ensure the extraction effect.
In general, the horizontal distance from the compaction area above the goaf to the roadway is less than the initial pressing step.Based on the research on the change rule of overlying rock cracks and gas migration rule in the goaf and the analysis of numerical simulation, three sets of high level boreholes are finally arranged to further improve the extraction effect of high level boreholes.
Determining drilling field spacing is an important means to ensure effective length lap between drilling holes, which can ensure efficient extraction, effective utilization, length and drilling field spacing of drilling holes.The calculation process is shown in Formula (3-5).
(3) According to the caving zone and fracture zone calculated in "Introduction" section and formula (3-5), it is determined that the single-hole pumping flow rate is 4 m 3 /min, four high-level drilling holes are constructed in each high-level drilling field, and three test high-level drilling fields are arranged in the working face of 24,130.The distance between the first high-level drilling field and the cutting hole is.The spacing between the drilling fields of the high level drilling holes is to ensure that the drilling stubble is above.Figure 9 shows the layout diagram of high extraction borehole.Detailed drilling layout parameters of the three high-position drilling fields are shown in Table 2.

Analysis of gas extraction effect
The proposed scheme was implemented for the later-stage gas extraction in the 24,130 working face, and the volume fraction of extracted gas was continuously monitored for a period of 3 months.The changes in gas volume fraction at the working face, return airway, and upper corner are summarized in Fig. 10.
According to Fig. 10, with the continuous advancement of the working face, the net amount of gas extracted from the upper borehole shows a trend of first increasing and then decreasing.By comparing the extraction concentrations of the four boreholes in Fig. 10a, it can be found that the average gas concentrations of boreholes  www.nature.com/scientificreports/2# and 3# are 11.68% and 12.72% respectively, while the average values of boreholes 1# and 4# are 5.89% and 5.26% respectively, and the average gas concentrations of boreholes 2# and 3# are 6.63% higher than those of boreholes 1# and 4#.It can be seen that when the Angle and length of the drilling hole remain unchanged, the dip Angle is 8°-11°, and the pumping effect is better when the Angle is between.It can be seen from Fig. 10b that the average gas extraction values of boreholes 5# and 7# are 10.92% and 13.24% respectively, while the average values of boreholes 6# and 8# are 6.21% and 6.13% respectively, and the average gas concentration values of boreholes 5# and 7# are 5.91% higher than those of boreholes 6# and 8#.The gas extraction effect is better when the Angle of borehole is between 30°and 42° when the inclination and length of borehole remain unchanged.
As can be seen from Fig. 10c, the average gas extraction values of boreholes 9# and 11# are 13.05% and 15.47% respectively, while the average values of boreholes 10# and 12# are 6.02% and 6.79% respectively, and the average gas concentration values of boreholes 9# and 11# are 7.86% higher than those of boreholes 10# and 12#.The gas extraction effect of boreholes 9# and 11# is obviously greater than that of boreholes 10# and 12#.It can be seen that when the Angle and inclination of boreholes remain unchanged, the extraction effect is better when the length of boreholes is between 145 and 155 m.Based on the above analysis, the optimal layout parameters for high-level boreholes are an inclination angle between 8° and 11°, an included angle between 30° and 42°, and a length between 145 and 155 m.
As can be seen from Table 3, numerical simulation results show that the optimal extraction borehole parameters are: dip Angle 9°-12°, included Angle 30°-45°, and borehole length 140-160 m.Field experimental results show that the optimal extraction borehole parameters are: The dip Angle is 8°-11°, the included Angle is 30°-42°, and the drilling length is 145-155 m.It can be seen that the field measured results are basically consistent with the numerical simulation data, which verifies the validity of the numerical model.By conducting real-time observations of the gas volume fractions in the mining face, upper corner angle, and return airway during gas extraction at drill sites 1#, 2#, and 3#, the changes in gas volume fractions in the mining face, upper corner angle, and return airway are shown in Fig. 11.As can be seen from Fig. 11, the gas extraction effect of the high drilling is significant, and the average gas concentration of the working face is 0.54% and the maximum gas concentration is 0.61%.The average gas concentration in the upper corner is 0.69%, and the maximum gas concentration is 0.8%.The average gas concentration in the return air roadway is 0.41%, and the maximum gas concentration is 0.55%, which is far lower than the maximum gas volume fraction 1% stipulated in the Coal mine Safety Regulations, and ensures the efficient production during the working face.

Result and discussion
Based on the analysis of the scope of the hidden area of gas accumulation in the upper part of 24,130 working face (the height of roof caving zone and caving fracture zone), the layout scheme of gas extraction borehole in the upper part of the working face is optimized from the Angle, Angle and length of borehole, which provides the test data basis for the design of gas extraction borehole in the upper part of the working face.The results of numerical simulation and field experiment show that the optimal layout parameters of the high drilling hole are the inclination of 8°-11°, the included Angle of 30°-42°, and the drilling length of 145-155 m.After the optimization of drilling fields 1#, 2# and 3#, the average gas extraction concentration of drilling holes is significantly increased, the gas concentration of working face, upper corner and return air roadway is significantly reduced, and the gas treatment effect of goaf is significantly improved.

Conclusion
(1) The analysis and calculation have determined that the height of the goaf zone in the 24,130 mining face is between 10.06 and 14.46 m, and the fracture zone height is between 38.75 and 49.95 m.This provides a theoretical basis for the establishment of a numerical simulation model for high-level borehole gas extraction (2) Based on the relevant physical parameters of Pingdingshan coal mine, a numerical simulation of in-situ borehole gas extraction was conducted, and the simulation results were analyzed.The results indicate that when the inclination angle of the borehole is controlled between 9° and 12°, the azimuth angle is controlled between 30° and 45°, and the borehole length is 150 m, the gas extraction effect of the fracture zone in the goaf is optimal (3) By investigating the extraction volume at different layouts of high-level boreholes and combining it with field measurements, the optimal placement position of the borehole's terminal point was determined.According to the field data, the optimal borehole layout parameters are an inclination angle between 8° and 11°, an azimuth angle between 30° and 42°, and a borehole length between 145 and 155 m, which is consistent with the simulation results (4) Through the monitoring of gas conditions in the mining face, upper corner angle, and return airway, it has been found that the maximum gas concentration in the 24,130 mining face and associated areas remains stable and controlled below 1%.The high-level borehole gas extraction has achieved the expected effect in the goaf gas, ensuring the safe and smooth mining of the working face.

Data availability
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Figure 2 .
Figure 2. Meshing of the model.

Figure 3 .
Figure 3. Variation curve of gas concentration at each monitoring point under different elevation angle.

Figure 4 .
Figure 4. Gas concentration distribution in 240 h extraction under different a 1 .

Figure 5 .
Figure 5. Variation curve of gas concentration at each monitoring point under different borehole azimuth.

Figure 6 .
Figure 6.Gas concentration distribution in 240 h extraction under different a 2 .

Figure 7 .
Figure 7. Variation curve of gas concentration at each monitoring point under different borehole length.

Figure 8 .
Figure 8. Gas concentration distribution in 240 h extraction under different L.

Figure 9 .
Figure 9. Schematic diagram of the drilling construction at the test site.

Figure 10 .
Figure 10.Variation curve of gas concentration in drilling site.

Table 1 .
Physical parameters of porous media.

Table 3 .
Comparison of results.