# The Drainage Horizon Determination of High Directional Long Borehole and Gas Control Effect Analysis

Authors: Yuqi Shang; Guiyi Wu; Qinzhi Liu; Dezhong Kong; Qiang Li
Journal: Advances in Civil Engineering
 (2021)

---

## Abstract

In order to effectively solve the problem of gas concentration overrun in the upper corner of goaf and tailentry during the mining of panel 9303 in Anshun Coal Mine, based on the advantages of controllable trajectory and wide coverage area of directional drilling technology, high directional long boreholes are arranged in tailentry 9303 to extract pressure relief gas. Firstly, the principle of high directional long borehole drainage technology is introduced, and the fracture evolution of overlying strata is obtained through using numerical simulation, theoretical calculation, and field practice, and the fracture evolution range is determined to be 6–12.69 m, and rationality of fracture height obtained by theoretical analysis and numerical simulation is verified by the method of field borehole peep observation. Through the analysis, it is concluded that the best location of the final hole is within the range of 6–12.69 m of the roof of coal seam 9#. The field practice has proved that the final hole position of the high directional long borehole is arranged at 12 m from the roof of coal seam 9#, and the average gas extraction concentration can reach 40%–50% after the borehole enters the stable extraction stage, the purity of gas extraction is up to 8.5 m3/min, and the gas concentration in the upper corner of panel 9303 is stable below 0.5% during mining, which achieves good gas drainage and control effect and provides a new way for gas control under similar geological conditions.


---

## Body

## 1. Introduction

The coal resource in Guizhou Province is very abundant, known as “Jiangnan Coal Sea.” The coal industry is one of the pillar industries in Guizhou Province, but coal seam gas content is generally high in Guizhou Province; coal and gas outburst disasters are serious. In particular, the number of coal layers is large, and the coal seam spacing is close, when the coal seam is mined, the pressure relief gas in the adjacent coal seam is easy to flow into the goaf of the working face through the delamination fracture and broken fracture. With the improvement and popularization of directional drilling technology and equipment in China, the use of high directional long boreholes technology to control coal mine gas problems is widely used in major coal mines.Coal seam mining causes the redistribution of stress field of overlying strata [1, 2], resulting in deformation and failure of overlying strata, thus forming mining-induced fractures. Therefore, study on fracture evolution law of overlying strata on working face is very important to improve gas control ability of high-level directional long borehole technology. To determine the drainage horizon of efficient extraction of high directional long borehole, many studies have analyzed the fracture evolution law of overlying strata in goaf during coal mining. Wang and Mang [3] proposed the gas control mode of goaf by roof directional long borehole instead of high drainage roadway, and it is concluded that compared with the high drainage roadway, the roof directional long drilling can not only achieve the same treatment effect, but also greatly reduce the capital investment and shorten the construction period. Finally, the best drainage position is determined by FLAC3D numerical simulation. Xu et al. [4] on the basis of predecessors, according to the actual working face roof overlying strata characteristics and mining conditions, the establishment of coal mining numerical simulation model, not only according to the results of numerical simulation and theoretical analysis to determine the high directional long drilling layout horizon, but also optimize the drilling structure, effectively reduces the gas concentration of the working face. Using the different-source prediction method, Cai et al. [5] analyzed gas emission source of working face, and the height range of caving zone and fracture zone in goaf is calculated theoretically, after the reasonable high drilling parameters are determined according to the specific situation of the mine. It was found that using roof strike high borehole gas drainage technology can effectively ensure the safe mining of working face. With the deepening of research on the overlying strata of working face, Cheng et al. [6] according to the simulation results of crack field, stress field, and strain field distribution, the range of fracture strengthening zone and compaction zone of goaf roof along the advancing direction of working face are divided. In their study, the borehole part in the fracture strengthening zone is the effective extraction area of boreholes; the length of the effective extraction section of boreholes and the development degree of fractures in the effective extraction section of boreholes jointly determine the extraction efficiency of long boreholes with high strike, revealing the mechanism of gas extraction by long boreholes with high strike in goaf roof. Li and Jin [7] simulate the mining process of panel 3304 in Zaozhuang Coal Mine by UDEC numerical simulation software, and combined with the theoretical analysis of rock fracture, according to the failure characteristics of coal seam roof and the change rule of displacement field in the simulation results, the maximum height position of roof rock caving zone and water flowing fault zone caused by mining is determined. The results show that the results obtained by field test are close to the numerical simulation results.The above research, whether numerical simulation or model test, has carried on the beneficial exploration and the research to the goaf overburden law according to the different actual working condition and the geological condition [8] and has obtained the massive research results. Therefore, on the basis of previous studies, this paper adopts the method of combining numerical simulation, field test, and theoretical analysis to calculate and analyze the distribution range of caving zone and fracture zone, so as to determine the location of holes. Finally, the effect of gas drainage is investigated through field practice to verify the distribution range of fracture zone, and the directional drilling technology is used to accurately control the drilling trajectory, so that the drilling trajectory is always in the fracture zone, so as to intercept and extract the pressure relief gas in the adjacent layer and control the gas in the upper corner of the goaf to achieve the purpose of safe production.


## 2. Engineering Background

The panel 9303 of the Anshun Coal Mine has a buried depth of 398 m and belongs to a high gas mine. The dip angle of the coal seam is near horizontal, which are 8#, 9#, and 10# coal seams (as shown in Figure1), respectively. The coal seam 8# and coal seam 10# are thin and therefore not mined. The mining coal seam 9# is 12.44–12.94 m below the 10# coal seam and 14–16.55 m above the coal seam 8#. The coal seam spacing is relatively close, so it is easy to cause disturbance to the adjacent coal seam due to the excavation of the working face, making the pressure relief gas flows to the working face through the cracks. The panel 9303 adopts the U-type ventilation mode. The goaf is easy to accumulate gas, and there is a risk of exceeding the gas volume fraction in the upper corner.Figure 1
Coal and rock strata histogram.The panel 9303 is located in coal seam 9# of 3rd panel, the panel length is 726 m, and the setup entry length is 173 m. The thickness of coal seam 9# is 1.52 m, and the recoverable reserves of panel 9303 are 290.3 million tons. Coal seam 9# roof is no false roof and direct top of silty mudstone, thickness 5.96 m. The basic roof is limestone with a thickness of 4.79 m. The direct bottom is fine sandstone with a thickness of 3.78 m. The old base is siltstone with a thickness of 8.12 m. The distance between coal seam 9# and coal seam 8# is about 15 m. The coal thickness of coal seam 8# is 0.2–0.58 m, and the average coal thickness is 0.43 m. The maximum gas content of coal seam 8# is 20.83 m3/t.


## 3. Drainage Position Determination

### 3.1. Principle of Gas Drainage Technology with High Directional Long Borehole

Directional drilling technology and equipment reflect the most advanced level of underground drilling in coal mines at this stage [9, 10]. Directional drilling is a kind of drilling method which uses natural bending law of drilling hole or artificial deflecting tool to extend drilling hole to predetermined target according to design requirement [11, 12].The principle of gas drainage by high directional drilling is based on the “three zones” distribution theory and “O” circle theory of overlying strata. When overlying strata move and destroy on the working face, “falling zone, fracture zone, and bending subsidence zone” are formed in the vertical direction, namely, “vertical three zones.” “Coal wall support influence zone, separation zone, and recompaction zone” are formed in the horizontal direction, namely, “transverse three zones” [13, 14]. The separation zone is the mining fracture “O” circle. The “O” circle and the middle and lower parts of the fault zone are the main channels of gas migration and the main places for gas accumulation [15].The horizon determination of high directional long borehole is crucial, which is directly related to the gas drainage effect, based on the distribution of the fracture zone of the overlying strata and “O” circle theory. On the one hand, the high directional long borehole cannot be arranged in the falling zone; otherwise, the rock crushing and caving will directly destroy the borehole, resulting in difficult drainage and even unable to extract. On the other hand, it cannot be arranged in the bending subsidence zone, because of the fracture development is poor in this zone, there is no effective gas flow channel between cracks, and the effect of gas extraction is poor. The best high directional long borehole arrangement horizon is the range of fracture zone in “three zones,” because of the development of fractures in this zone, fractures are connected to each other to form an effective channel for gas flow, and will not cause serious damage to boreholes, which is conducive to gas extraction [16]. However, due to the different geological and mining conditions of working faces in different mining areas, the height of the “three zones” formed by the collapse of the overlying strata is also different. Therefore, in order to reasonably determine the high directional long boreholes arrangement horizon, it is necessary to determine the caving characteristics of overlying strata, the evolution of fractures, and the height of “three zones” of overlying strata.


### 3.2. Theoretical Calculation of “Three Zones” Height of Overlying Strata

The height of overlying strata falling zone and fracture zone is determined by the mining height of working face, coal seam dip angle, roof management method, roof lithology, and other factors. Considering the actual situation of working face, since the roof of coal seam 9# is silty mudstone, it belongs to soft rock. So, the following empirical formula is selected to calculate the height of falling zone and fracture zone [17, 18]:(1)H1=100M6.2M+10.0±2.5,(2)H2=100M3.1M+6.0±6.5.In the formula,H1 is the falling zone height. H2 is the fracture zone height. M is the thickness of coal seam mining; take the average thickness of 1.5 m into formulas (1) and (2), the height of falling zone H1 is 5.27–10.27 m, and the height of fracture zone H2 is 7.58–20.58 m.


### 3.3. Numerical Simulation of Roof Crack Evolution in Goaf during Mining

In order to further accurately determine the fracture development characteristics of pamel 9303 under different advancing distances, based on the coal and rock strata histogram, physical and mechanical properties of coal strata, and panel 9303 mining conditions, the UDEC numerical simulation method is used to establish the evolution model of overlying strata fracture in the process of coal seam mining, so as to study the evolution characteristics of overlying strata fracture under different advancing distances and finally determine the most efficient hole arrangement horizon.


#### 3.3.1. Establishment of Numerical Simulation Model

By cutting the blocks, UDEC separates the deformed materials into discrete block sets to represent the discontinuous medium. The block in the model uses the Mohr–Coulomb criterion to describe the deformation and failure process of rock materials [19], as shown in the following formula:(3)R=c⋅cosφ+12σ1+σ3⋅sinφ.In the formula,φ is the angle of internal friction of rock, c is rock cohesion, σ1 is the maximum principal stress, σ3 is the minimum principal stress, and R is the radius of stress circle, R = (σ1 − σ3)/2. And according to the geological prospecting report of Anshun Coal Mine, the mechanical parameters of the roof-floor of coal seam 9# are obtained, as shown in Table 1.Table 1
The mechanical parameters of major rock mass for roof and floor.
Name of rockDensity (kg·m−3)Bulk modulus (GPa)Shear modulus (GPa)Cohesion (MPa)Angle of internal friction (°)Tensile strength (MPa)Sandstone26005.27112.3394.35Fine sandstone26454.62188.36342.35Siltstone27005.675.3352.86Mudstone25504.164.691.41373.2Coal seam 8#13503.952.22301.04Siltstone27005.63.85.3352.86Limestone25406.851.471.26101.2Silty mudstone213268.138.87424.56Limestone25406.851.471.26101.2Silty mudstone213268.138.87424.56Coal seam 9#13004.335.12381.04Fine sandstone224524.62126.36342.35Siltstone21005.675.3352.86Coal seam 10#13503.952.22301.04Medium sandstone2132268.1311.87424.56Based on this, establish a simplified two-dimensional model by UDEC. The model is 200 m in length and 120 m in height. Displacement constraints are applied to the vertical direction and lower boundary of the model, and equivalent stress loads are applied to the upper boundary. According to the calculation, the equivalent load is 7.11 MPa. The advance length of the model is 120 m for excavation, and the development law of rock strata caving cracks is simulated. 40 m is left at both ends to offset the boundary effect (as shown in Figure2).Figure 2
Numerical calculation model.


#### 3.3.2. Analysis of Roof Fracture Development

In the process of working face advancing, due to the influence of excavation work, the stress field is redistributed, resulting in the continuous evolution of mining cracks with excavation work [20–23]. Due to the supporting effect of coal wall at both ends of working face and the different lithology of rock strata, the collapse and subsidence degree of each rock strata are different, so the separation of rock strata will occur first; thus, the separation cracks will appear first. In general, the gas flow and concentration in borehole extraction will be significantly higher than before after the emergence of separation cracks. With the further excavation of the working face, after the large-scale collapse of the overlying strata, the fracture cracks are formed at both ends of the open cut hole and the working face. With the further increase of the excavation distance, the cracks in the upper central area of the goaf are gradually compacted, and the separation cracks begin to evolve upward. These fractures are interconnected and gradually form fracture channels, resulting in the continuous inflow of free gas from the adjacent layer and surrounding rock into the goaf through this channel, resulting in the gas concentration overrun in the working face.In order to analyze and study the overlying strata fracture evolution law, the fracture evolution diagram of working face advancing to 30 m, 50 m, 80 m, and 110 m is selected for analysis.(1)
When the working face advances 30 m, it can be seen from Figure3 that the immediate roof reaches the bearing limit and begins to collapse due to the excavation work. When the working face is pressed for the first time, the stress of the rock strata is released, and under the action of the dead weight, bending and sinking occur, and the main roof cracks begin to develop. A small amount of delamination fracture appears in the overlying strata and the disturbed coal seam 8#. Broken fissures are formed at both ends of the working face and the open-off cut, and the delamination fracture in the overlying strata shows a parallel distribution. However, the evolution of the fissures has just begun, and the penetrating fissures have not yet occurred.(2)
When the working face advances 50 m, it can be seen from Figure4 that the collapse range of the immediate roof and the main roof increases obviously, and the immediate roof has collapsed, the main roof is fractured, and the working face is under periodic weighting. The upper coal seam 8# is deformed, and the delamination fracture continues to expand upward. The key layer fracture is increasing. The rock strata near the open-off cut side and the working face side are supported by the coal wall to form a shear effect on the upper part, forming a shear failure and then forming a broken fissure.(3)
When the working face advances 80 m, it can be seen from Figure5 that the roof fracture further expands to the overlying strata, and the overlying rock fracture area gradually expands. The delamination fracture and broken fissures continue to penetrate, and the key layer fracture develops well. The compaction area gradually forms in the middle of the overlying strata.(4)
When the working face advances 110 m, it can be seen from Figure6 that the delamination fracture begins to develop laterally, and the fracture in the middle of the goaf gradually disappears. The delamination fracture at the top of the fracture field is gradually replaced by the broken fracture at the bottom. With the continuous advancement of the working face, after 110 m mining in the working face, the fracture distribution form and scope of the overlying strata are gradually stable in the vertical direction. In the horizontal direction, the development of overlying rock fractures moves forward synchronously with the mining of working face, but does not continue to extend upward in the advancing direction. At this time, the fracture of overlying strata develops completely.Figure 3
Evolution diagram of overlying strata fracture when working face advances to 30 m.Figure 4
Evolution diagram of overlying strata fracture when working face advances to 50 m.Figure 5
Evolution diagram of overlying strata fracture when working face advances to 80 m.Figure 6
Distribution of fracture strengthening zone and compaction zone in overlying strata after stable fracture development.In summary, combined with the distribution of the fracture strengthening area and the compaction area of the overlying strata, the fracture strengthening area of the roof of the goaf is the main accumulation area of the pressure relief gas in the upper coal seam 8#. The gas content is high and the concentration is large, which is the best choice for the hole arrangement layer of the drainage borehole. Through the simulation results when the working face is excavated to 110 m, it can be obtained that the best fracture development of the overlying strata is with the height of 12.69 m.The above empirical formula calculation and numerical simulation results show that when the working face advances to 110 m, the working face reaches full mining, the height of the “three zone” of the overlying strata is basically stable, and the maximum fracture zone is 12.69 m. Therefore, the high directional long borehole should be arranged in the position of about 12 m of the coal seam roof.


#### 3.3.3. Mining Fracture Height by Field Test

The borehole peeping method has the advantages of intuitive and convenient, which is widely used in coal mines to detect the geological structure of surrounding rock and the development of fractures [24, 25]. This borehole peeping adopts the intrinsically safe borehole peeper for mining. The minimum diameter of the peeping borehole is 25 mm, and the length of the observable borehole is 30 m. The working principle is to use the probe to convert light into electronic signals, and the signal is transmitted to the image receiver through the cable, so as to observe, record, and store the rock mass structure image in the borehole [26, 27].This study mainly wants to obtain the fracture evolution height of roof strata affected by mining in panel 9303 during mining, so as to verify the rationality of numerical simulation and empirical formula calculation [28]. When arranging the peep boreholes test points, considering that a large amount of dust will be generated during the drilling construction, the test points are selected in the return airway. Combined with the actual position and advancing speed of the working face, two boreholes are constructed at 50 m, 100 m, and 150 m away from the open-off cut as the test points of peep boreholes. The layout of the observation borehole is shown in Figure 7, two boreholes are in one group, and the plane surface of each group is shown in Figure 8. The borehole design of each group is the same. According to the actual situation of Anshun Mine, the vertical height of the final borehole position is determined to be 18 m.Figure 7
Peeping borehole layout schematic diagram.Figure 8
Peeping hole layout plane diagram.In order to analyze the development of fracture at different heights, so as to determine the fracture zone of the overlying strata on the roof, the borehole images of the same borehole with the advance of the working face in different periods and depths are intercepted, as shown in Figures9 and 10.Figure 9
Drilling peep at different depths of 1-2# observation hole when working face advances to 50 m. (a) Final hole height 2 m. (b) Final hole height 6 m. (c) Final hole height 10 m. (d) Final hole height 15 m.
(a)(b)(c)(d)Figure 10
Drilling peep at different depths of 2-2# observation hole when working face advances to 100 m. (a) Final hole height 2 m. (b) Final hole height 6 m. (c) Final hole height 10 m. (d) Final hole height 15 m.
(a)(b)(c)(d)From Figures9 and 10, it can be seen that, from the opening up into the fracture zone, the fracture is more developed, the hole wall is broken until 10 m, and then to the upper rock layer is more complete, but at the end hole height of 15 m, there is a more broken section, which may be broken again after mining stress compression shear failure. The upper hole wall is basically intact, only a few places have slight damage. Therefore, the borehole observation results show that the height of water flowing fractured zone is not more than 15 m.In summary, according to the same drilling hole in different hole depth peep, the analysis found that the fractures are mainly produced in the range of 6 m to 15 m above the roof, and the range changes with the advance of the working face. Through the analysis of images with different advancing distances and different borehole depths, it is found that when the working face advances to 100 m, according to the peep images of 3-1# and 3-2# two boreholes, the fracture development situation shown by the images of 2-2# boreholes is similar. It can be concluded that the fracture development height gradually tends to be stable. Fracture zone range is within 6–15 m.


## 3.1. Principle of Gas Drainage Technology with High Directional Long Borehole

Directional drilling technology and equipment reflect the most advanced level of underground drilling in coal mines at this stage [9, 10]. Directional drilling is a kind of drilling method which uses natural bending law of drilling hole or artificial deflecting tool to extend drilling hole to predetermined target according to design requirement [11, 12].The principle of gas drainage by high directional drilling is based on the “three zones” distribution theory and “O” circle theory of overlying strata. When overlying strata move and destroy on the working face, “falling zone, fracture zone, and bending subsidence zone” are formed in the vertical direction, namely, “vertical three zones.” “Coal wall support influence zone, separation zone, and recompaction zone” are formed in the horizontal direction, namely, “transverse three zones” [13, 14]. The separation zone is the mining fracture “O” circle. The “O” circle and the middle and lower parts of the fault zone are the main channels of gas migration and the main places for gas accumulation [15].The horizon determination of high directional long borehole is crucial, which is directly related to the gas drainage effect, based on the distribution of the fracture zone of the overlying strata and “O” circle theory. On the one hand, the high directional long borehole cannot be arranged in the falling zone; otherwise, the rock crushing and caving will directly destroy the borehole, resulting in difficult drainage and even unable to extract. On the other hand, it cannot be arranged in the bending subsidence zone, because of the fracture development is poor in this zone, there is no effective gas flow channel between cracks, and the effect of gas extraction is poor. The best high directional long borehole arrangement horizon is the range of fracture zone in “three zones,” because of the development of fractures in this zone, fractures are connected to each other to form an effective channel for gas flow, and will not cause serious damage to boreholes, which is conducive to gas extraction [16]. However, due to the different geological and mining conditions of working faces in different mining areas, the height of the “three zones” formed by the collapse of the overlying strata is also different. Therefore, in order to reasonably determine the high directional long boreholes arrangement horizon, it is necessary to determine the caving characteristics of overlying strata, the evolution of fractures, and the height of “three zones” of overlying strata.


## 3.2. Theoretical Calculation of “Three Zones” Height of Overlying Strata

The height of overlying strata falling zone and fracture zone is determined by the mining height of working face, coal seam dip angle, roof management method, roof lithology, and other factors. Considering the actual situation of working face, since the roof of coal seam 9# is silty mudstone, it belongs to soft rock. So, the following empirical formula is selected to calculate the height of falling zone and fracture zone [17, 18]:(1)H1=100M6.2M+10.0±2.5,(2)H2=100M3.1M+6.0±6.5.In the formula,H1 is the falling zone height. H2 is the fracture zone height. M is the thickness of coal seam mining; take the average thickness of 1.5 m into formulas (1) and (2), the height of falling zone H1 is 5.27–10.27 m, and the height of fracture zone H2 is 7.58–20.58 m.


## 3.3. Numerical Simulation of Roof Crack Evolution in Goaf during Mining

In order to further accurately determine the fracture development characteristics of pamel 9303 under different advancing distances, based on the coal and rock strata histogram, physical and mechanical properties of coal strata, and panel 9303 mining conditions, the UDEC numerical simulation method is used to establish the evolution model of overlying strata fracture in the process of coal seam mining, so as to study the evolution characteristics of overlying strata fracture under different advancing distances and finally determine the most efficient hole arrangement horizon.


### 3.3.1. Establishment of Numerical Simulation Model

By cutting the blocks, UDEC separates the deformed materials into discrete block sets to represent the discontinuous medium. The block in the model uses the Mohr–Coulomb criterion to describe the deformation and failure process of rock materials [19], as shown in the following formula:(3)R=c⋅cosφ+12σ1+σ3⋅sinφ.In the formula,φ is the angle of internal friction of rock, c is rock cohesion, σ1 is the maximum principal stress, σ3 is the minimum principal stress, and R is the radius of stress circle, R = (σ1 − σ3)/2. And according to the geological prospecting report of Anshun Coal Mine, the mechanical parameters of the roof-floor of coal seam 9# are obtained, as shown in Table 1.Table 1
The mechanical parameters of major rock mass for roof and floor.
Name of rockDensity (kg·m−3)Bulk modulus (GPa)Shear modulus (GPa)Cohesion (MPa)Angle of internal friction (°)Tensile strength (MPa)Sandstone26005.27112.3394.35Fine sandstone26454.62188.36342.35Siltstone27005.675.3352.86Mudstone25504.164.691.41373.2Coal seam 8#13503.952.22301.04Siltstone27005.63.85.3352.86Limestone25406.851.471.26101.2Silty mudstone213268.138.87424.56Limestone25406.851.471.26101.2Silty mudstone213268.138.87424.56Coal seam 9#13004.335.12381.04Fine sandstone224524.62126.36342.35Siltstone21005.675.3352.86Coal seam 10#13503.952.22301.04Medium sandstone2132268.1311.87424.56Based on this, establish a simplified two-dimensional model by UDEC. The model is 200 m in length and 120 m in height. Displacement constraints are applied to the vertical direction and lower boundary of the model, and equivalent stress loads are applied to the upper boundary. According to the calculation, the equivalent load is 7.11 MPa. The advance length of the model is 120 m for excavation, and the development law of rock strata caving cracks is simulated. 40 m is left at both ends to offset the boundary effect (as shown in Figure2).Figure 2
Numerical calculation model.


### 3.3.2. Analysis of Roof Fracture Development

In the process of working face advancing, due to the influence of excavation work, the stress field is redistributed, resulting in the continuous evolution of mining cracks with excavation work [20–23]. Due to the supporting effect of coal wall at both ends of working face and the different lithology of rock strata, the collapse and subsidence degree of each rock strata are different, so the separation of rock strata will occur first; thus, the separation cracks will appear first. In general, the gas flow and concentration in borehole extraction will be significantly higher than before after the emergence of separation cracks. With the further excavation of the working face, after the large-scale collapse of the overlying strata, the fracture cracks are formed at both ends of the open cut hole and the working face. With the further increase of the excavation distance, the cracks in the upper central area of the goaf are gradually compacted, and the separation cracks begin to evolve upward. These fractures are interconnected and gradually form fracture channels, resulting in the continuous inflow of free gas from the adjacent layer and surrounding rock into the goaf through this channel, resulting in the gas concentration overrun in the working face.In order to analyze and study the overlying strata fracture evolution law, the fracture evolution diagram of working face advancing to 30 m, 50 m, 80 m, and 110 m is selected for analysis.(1)
When the working face advances 30 m, it can be seen from Figure3 that the immediate roof reaches the bearing limit and begins to collapse due to the excavation work. When the working face is pressed for the first time, the stress of the rock strata is released, and under the action of the dead weight, bending and sinking occur, and the main roof cracks begin to develop. A small amount of delamination fracture appears in the overlying strata and the disturbed coal seam 8#. Broken fissures are formed at both ends of the working face and the open-off cut, and the delamination fracture in the overlying strata shows a parallel distribution. However, the evolution of the fissures has just begun, and the penetrating fissures have not yet occurred.(2)
When the working face advances 50 m, it can be seen from Figure4 that the collapse range of the immediate roof and the main roof increases obviously, and the immediate roof has collapsed, the main roof is fractured, and the working face is under periodic weighting. The upper coal seam 8# is deformed, and the delamination fracture continues to expand upward. The key layer fracture is increasing. The rock strata near the open-off cut side and the working face side are supported by the coal wall to form a shear effect on the upper part, forming a shear failure and then forming a broken fissure.(3)
When the working face advances 80 m, it can be seen from Figure5 that the roof fracture further expands to the overlying strata, and the overlying rock fracture area gradually expands. The delamination fracture and broken fissures continue to penetrate, and the key layer fracture develops well. The compaction area gradually forms in the middle of the overlying strata.(4)
When the working face advances 110 m, it can be seen from Figure6 that the delamination fracture begins to develop laterally, and the fracture in the middle of the goaf gradually disappears. The delamination fracture at the top of the fracture field is gradually replaced by the broken fracture at the bottom. With the continuous advancement of the working face, after 110 m mining in the working face, the fracture distribution form and scope of the overlying strata are gradually stable in the vertical direction. In the horizontal direction, the development of overlying rock fractures moves forward synchronously with the mining of working face, but does not continue to extend upward in the advancing direction. At this time, the fracture of overlying strata develops completely.Figure 3
Evolution diagram of overlying strata fracture when working face advances to 30 m.Figure 4
Evolution diagram of overlying strata fracture when working face advances to 50 m.Figure 5
Evolution diagram of overlying strata fracture when working face advances to 80 m.Figure 6
Distribution of fracture strengthening zone and compaction zone in overlying strata after stable fracture development.In summary, combined with the distribution of the fracture strengthening area and the compaction area of the overlying strata, the fracture strengthening area of the roof of the goaf is the main accumulation area of the pressure relief gas in the upper coal seam 8#. The gas content is high and the concentration is large, which is the best choice for the hole arrangement layer of the drainage borehole. Through the simulation results when the working face is excavated to 110 m, it can be obtained that the best fracture development of the overlying strata is with the height of 12.69 m.The above empirical formula calculation and numerical simulation results show that when the working face advances to 110 m, the working face reaches full mining, the height of the “three zone” of the overlying strata is basically stable, and the maximum fracture zone is 12.69 m. Therefore, the high directional long borehole should be arranged in the position of about 12 m of the coal seam roof.


### 3.3.3. Mining Fracture Height by Field Test

The borehole peeping method has the advantages of intuitive and convenient, which is widely used in coal mines to detect the geological structure of surrounding rock and the development of fractures [24, 25]. This borehole peeping adopts the intrinsically safe borehole peeper for mining. The minimum diameter of the peeping borehole is 25 mm, and the length of the observable borehole is 30 m. The working principle is to use the probe to convert light into electronic signals, and the signal is transmitted to the image receiver through the cable, so as to observe, record, and store the rock mass structure image in the borehole [26, 27].This study mainly wants to obtain the fracture evolution height of roof strata affected by mining in panel 9303 during mining, so as to verify the rationality of numerical simulation and empirical formula calculation [28]. When arranging the peep boreholes test points, considering that a large amount of dust will be generated during the drilling construction, the test points are selected in the return airway. Combined with the actual position and advancing speed of the working face, two boreholes are constructed at 50 m, 100 m, and 150 m away from the open-off cut as the test points of peep boreholes. The layout of the observation borehole is shown in Figure 7, two boreholes are in one group, and the plane surface of each group is shown in Figure 8. The borehole design of each group is the same. According to the actual situation of Anshun Mine, the vertical height of the final borehole position is determined to be 18 m.Figure 7
Peeping borehole layout schematic diagram.Figure 8
Peeping hole layout plane diagram.In order to analyze the development of fracture at different heights, so as to determine the fracture zone of the overlying strata on the roof, the borehole images of the same borehole with the advance of the working face in different periods and depths are intercepted, as shown in Figures9 and 10.Figure 9
Drilling peep at different depths of 1-2# observation hole when working face advances to 50 m. (a) Final hole height 2 m. (b) Final hole height 6 m. (c) Final hole height 10 m. (d) Final hole height 15 m.
(a)(b)(c)(d)Figure 10
Drilling peep at different depths of 2-2# observation hole when working face advances to 100 m. (a) Final hole height 2 m. (b) Final hole height 6 m. (c) Final hole height 10 m. (d) Final hole height 15 m.
(a)(b)(c)(d)From Figures9 and 10, it can be seen that, from the opening up into the fracture zone, the fracture is more developed, the hole wall is broken until 10 m, and then to the upper rock layer is more complete, but at the end hole height of 15 m, there is a more broken section, which may be broken again after mining stress compression shear failure. The upper hole wall is basically intact, only a few places have slight damage. Therefore, the borehole observation results show that the height of water flowing fractured zone is not more than 15 m.In summary, according to the same drilling hole in different hole depth peep, the analysis found that the fractures are mainly produced in the range of 6 m to 15 m above the roof, and the range changes with the advance of the working face. Through the analysis of images with different advancing distances and different borehole depths, it is found that when the working face advances to 100 m, according to the peep images of 3-1# and 3-2# two boreholes, the fracture development situation shown by the images of 2-2# boreholes is similar. It can be concluded that the fracture development height gradually tends to be stable. Fracture zone range is within 6–15 m.


## 3.3.1. Establishment of Numerical Simulation Model

By cutting the blocks, UDEC separates the deformed materials into discrete block sets to represent the discontinuous medium. The block in the model uses the Mohr–Coulomb criterion to describe the deformation and failure process of rock materials [19], as shown in the following formula:(3)R=c⋅cosφ+12σ1+σ3⋅sinφ.In the formula,φ is the angle of internal friction of rock, c is rock cohesion, σ1 is the maximum principal stress, σ3 is the minimum principal stress, and R is the radius of stress circle, R = (σ1 − σ3)/2. And according to the geological prospecting report of Anshun Coal Mine, the mechanical parameters of the roof-floor of coal seam 9# are obtained, as shown in Table 1.Table 1
The mechanical parameters of major rock mass for roof and floor.
Name of rockDensity (kg·m−3)Bulk modulus (GPa)Shear modulus (GPa)Cohesion (MPa)Angle of internal friction (°)Tensile strength (MPa)Sandstone26005.27112.3394.35Fine sandstone26454.62188.36342.35Siltstone27005.675.3352.86Mudstone25504.164.691.41373.2Coal seam 8#13503.952.22301.04Siltstone27005.63.85.3352.86Limestone25406.851.471.26101.2Silty mudstone213268.138.87424.56Limestone25406.851.471.26101.2Silty mudstone213268.138.87424.56Coal seam 9#13004.335.12381.04Fine sandstone224524.62126.36342.35Siltstone21005.675.3352.86Coal seam 10#13503.952.22301.04Medium sandstone2132268.1311.87424.56Based on this, establish a simplified two-dimensional model by UDEC. The model is 200 m in length and 120 m in height. Displacement constraints are applied to the vertical direction and lower boundary of the model, and equivalent stress loads are applied to the upper boundary. According to the calculation, the equivalent load is 7.11 MPa. The advance length of the model is 120 m for excavation, and the development law of rock strata caving cracks is simulated. 40 m is left at both ends to offset the boundary effect (as shown in Figure2).Figure 2
Numerical calculation model.


## 3.3.2. Analysis of Roof Fracture Development

In the process of working face advancing, due to the influence of excavation work, the stress field is redistributed, resulting in the continuous evolution of mining cracks with excavation work [20–23]. Due to the supporting effect of coal wall at both ends of working face and the different lithology of rock strata, the collapse and subsidence degree of each rock strata are different, so the separation of rock strata will occur first; thus, the separation cracks will appear first. In general, the gas flow and concentration in borehole extraction will be significantly higher than before after the emergence of separation cracks. With the further excavation of the working face, after the large-scale collapse of the overlying strata, the fracture cracks are formed at both ends of the open cut hole and the working face. With the further increase of the excavation distance, the cracks in the upper central area of the goaf are gradually compacted, and the separation cracks begin to evolve upward. These fractures are interconnected and gradually form fracture channels, resulting in the continuous inflow of free gas from the adjacent layer and surrounding rock into the goaf through this channel, resulting in the gas concentration overrun in the working face.In order to analyze and study the overlying strata fracture evolution law, the fracture evolution diagram of working face advancing to 30 m, 50 m, 80 m, and 110 m is selected for analysis.(1)
When the working face advances 30 m, it can be seen from Figure3 that the immediate roof reaches the bearing limit and begins to collapse due to the excavation work. When the working face is pressed for the first time, the stress of the rock strata is released, and under the action of the dead weight, bending and sinking occur, and the main roof cracks begin to develop. A small amount of delamination fracture appears in the overlying strata and the disturbed coal seam 8#. Broken fissures are formed at both ends of the working face and the open-off cut, and the delamination fracture in the overlying strata shows a parallel distribution. However, the evolution of the fissures has just begun, and the penetrating fissures have not yet occurred.(2)
When the working face advances 50 m, it can be seen from Figure4 that the collapse range of the immediate roof and the main roof increases obviously, and the immediate roof has collapsed, the main roof is fractured, and the working face is under periodic weighting. The upper coal seam 8# is deformed, and the delamination fracture continues to expand upward. The key layer fracture is increasing. The rock strata near the open-off cut side and the working face side are supported by the coal wall to form a shear effect on the upper part, forming a shear failure and then forming a broken fissure.(3)
When the working face advances 80 m, it can be seen from Figure5 that the roof fracture further expands to the overlying strata, and the overlying rock fracture area gradually expands. The delamination fracture and broken fissures continue to penetrate, and the key layer fracture develops well. The compaction area gradually forms in the middle of the overlying strata.(4)
When the working face advances 110 m, it can be seen from Figure6 that the delamination fracture begins to develop laterally, and the fracture in the middle of the goaf gradually disappears. The delamination fracture at the top of the fracture field is gradually replaced by the broken fracture at the bottom. With the continuous advancement of the working face, after 110 m mining in the working face, the fracture distribution form and scope of the overlying strata are gradually stable in the vertical direction. In the horizontal direction, the development of overlying rock fractures moves forward synchronously with the mining of working face, but does not continue to extend upward in the advancing direction. At this time, the fracture of overlying strata develops completely.Figure 3
Evolution diagram of overlying strata fracture when working face advances to 30 m.Figure 4
Evolution diagram of overlying strata fracture when working face advances to 50 m.Figure 5
Evolution diagram of overlying strata fracture when working face advances to 80 m.Figure 6
Distribution of fracture strengthening zone and compaction zone in overlying strata after stable fracture development.In summary, combined with the distribution of the fracture strengthening area and the compaction area of the overlying strata, the fracture strengthening area of the roof of the goaf is the main accumulation area of the pressure relief gas in the upper coal seam 8#. The gas content is high and the concentration is large, which is the best choice for the hole arrangement layer of the drainage borehole. Through the simulation results when the working face is excavated to 110 m, it can be obtained that the best fracture development of the overlying strata is with the height of 12.69 m.The above empirical formula calculation and numerical simulation results show that when the working face advances to 110 m, the working face reaches full mining, the height of the “three zone” of the overlying strata is basically stable, and the maximum fracture zone is 12.69 m. Therefore, the high directional long borehole should be arranged in the position of about 12 m of the coal seam roof.


## 3.3.3. Mining Fracture Height by Field Test

The borehole peeping method has the advantages of intuitive and convenient, which is widely used in coal mines to detect the geological structure of surrounding rock and the development of fractures [24, 25]. This borehole peeping adopts the intrinsically safe borehole peeper for mining. The minimum diameter of the peeping borehole is 25 mm, and the length of the observable borehole is 30 m. The working principle is to use the probe to convert light into electronic signals, and the signal is transmitted to the image receiver through the cable, so as to observe, record, and store the rock mass structure image in the borehole [26, 27].This study mainly wants to obtain the fracture evolution height of roof strata affected by mining in panel 9303 during mining, so as to verify the rationality of numerical simulation and empirical formula calculation [28]. When arranging the peep boreholes test points, considering that a large amount of dust will be generated during the drilling construction, the test points are selected in the return airway. Combined with the actual position and advancing speed of the working face, two boreholes are constructed at 50 m, 100 m, and 150 m away from the open-off cut as the test points of peep boreholes. The layout of the observation borehole is shown in Figure 7, two boreholes are in one group, and the plane surface of each group is shown in Figure 8. The borehole design of each group is the same. According to the actual situation of Anshun Mine, the vertical height of the final borehole position is determined to be 18 m.Figure 7
Peeping borehole layout schematic diagram.Figure 8
Peeping hole layout plane diagram.In order to analyze the development of fracture at different heights, so as to determine the fracture zone of the overlying strata on the roof, the borehole images of the same borehole with the advance of the working face in different periods and depths are intercepted, as shown in Figures9 and 10.Figure 9
Drilling peep at different depths of 1-2# observation hole when working face advances to 50 m. (a) Final hole height 2 m. (b) Final hole height 6 m. (c) Final hole height 10 m. (d) Final hole height 15 m.
(a)(b)(c)(d)Figure 10
Drilling peep at different depths of 2-2# observation hole when working face advances to 100 m. (a) Final hole height 2 m. (b) Final hole height 6 m. (c) Final hole height 10 m. (d) Final hole height 15 m.
(a)(b)(c)(d)From Figures9 and 10, it can be seen that, from the opening up into the fracture zone, the fracture is more developed, the hole wall is broken until 10 m, and then to the upper rock layer is more complete, but at the end hole height of 15 m, there is a more broken section, which may be broken again after mining stress compression shear failure. The upper hole wall is basically intact, only a few places have slight damage. Therefore, the borehole observation results show that the height of water flowing fractured zone is not more than 15 m.In summary, according to the same drilling hole in different hole depth peep, the analysis found that the fractures are mainly produced in the range of 6 m to 15 m above the roof, and the range changes with the advance of the working face. Through the analysis of images with different advancing distances and different borehole depths, it is found that when the working face advances to 100 m, according to the peep images of 3-1# and 3-2# two boreholes, the fracture development situation shown by the images of 2-2# boreholes is similar. It can be concluded that the fracture development height gradually tends to be stable. Fracture zone range is within 6–15 m.


## 4. Parameter Design of Long Height Directional Drilling

Comprehensively considering the results of empirical formula calculation and numerical simulation, and according to the field measurement verification, it can be seen that the hole arrangement layer of gas directional drilling in goaf of panel 9303 is selected at 12 m above coal seam 9#. Due to the first use of directional drilling rig in this mine, in order to simplify the construction technology, the directional high drilling design adopts single-hole parallel arrangement and does not open branch holes.In order to verify the above theoretical research results, the directional drilling rig drilling site is constructed in the lower side of tailentry 9303 (the drilling site position is shown in Figure11). The drilling site specification is wide × deep × high = 10 × 5 × 3 m. The opening position is located at 0.5–1 m of the coal seam roof and the opening spacing is 0.5–0.8 m. The design inclination of the opening section is 10–14°. After drilling the sandstone on the roof of coal seam 9#, the drilling direction is adjusted until the drilling reaches the predetermined position: 12 m above coal seam 9# (as shown in Figure 12), and the final hole spacing is 7 m. The parameters for five specific boreholes are as follows:(1)
Design parameters of 1-1 borehole: opening angle is 14°, opening azimuth is 331.9°, opening height is 0.3 m above coal seam roof, final hole height is 12 m above coal seam, and final design hole depth is 325 m. The final hole is 8 m away from the roadway side.(2)
Design parameters of 1-2 borehole: opening angle is 14°, opening azimuth is 328.9°, opening height is 0.3 m above coal seam roof, final hole height is 12 m above coal seam, and final design hole depth is 335 m. The final hole is 15 m away from the roadway side.(3)
Design parameters of 1–3 borehole: opening angle is 14°, opening azimuth is 325.9°, opening height is 0.3 m above coal seam roof, final hole height is 12 m above coal seam, and final design hole depth is 335 m. The final hole is 22 m away from the roadway side.(4)
Design parameters of 1–4 borehole: opening angle is 14°, opening azimuth is 322.9°, opening height is 0.3 m above coal seam roof, final hole height is 12 m above coal seam, and final design hole depth is 335 m. The final hole is 29 m away from the roadway side.(5)
Design parameters of 1–5 borehole: opening angle is 14°, opening azimuth is 319.9°, opening height is 0.3 m above coal seam roof, final hole height is 12 m above coal seam, and final design hole depth is 335 m. The final hole is 39 m away from the roadway side.Figure 11
Indication diagram of directional long drilling field layout.Figure 12
Design plane surface of directional drilling in return airway.


## 5. Extraction Effect Analysis

During the mining of panel 9303, the drainage effect of five high directional boreholes constructed by directional drilling field is significant. Gas drainage is shown in Figure10. The extraction flow and gas concentration of five boreholes at different advancing distances were counted (as shown in Figure 13). When the mining face is 0∼30 m, the gas concentration of borehole extraction increases significantly, and the pure gas flow rate shows a rising trend to 5.5 m3/min. At this stage, panel 9303 overlying strata has begun to form cracks due to the mining influence, and with the continuous development of cracks, drainage concentration and drainage gas flow are also rising. Subsequently, when the working face advances from 30 m to 110 m, the concentration of gas drainage and the pure flow of gas drainage gradually increase. Since the overlying strata activities of the stope are continuously intensified, this causes disturbance to the upper coal seam 8# and makes the gas in the coal seam 8# continue to flow downward through the intertwined separated fractures and broken fractures. When the working face advances to 30 m, the first weighting of the working face, the expansion deformation of the coal seam 8#, and the gas pressure relief occur. When the working face advances to 50 m, periodic weighting occurs. Before the periodic weighting, the borehole drainage flow and drainage concentration have a significant downward trend. After the periodic weighting, the borehole drainage concentration and drainage flow have a significant upward trend and then tend to be stable until the next periodic weighting. Subsequently, when the working face advances to 110 m, the overall drainage effect of the borehole tends to be stable, entering the stage of efficient and lasting stable drainage. The drainage concentration reaches about 43%, and the maximum pure gas flow rate reaches 8.5 m3/min.Figure 13
Gas extraction by directional long borehole.During the mining period, the gas concentration in the return airway of the working face was 0.25%–0.4%, and the gas concentration in the upper corner was also lower than 0.5%, as shown in Figure14. The high directional long borehole can effectively prevent the pressure relief gas of coal from 8# entering the return air flow through the fracture into the goaf, which leads to the problem of gas concentration overrun in the working face. At the same time, compared with the traditional high pumping roadway, the preparation time of the working face is shortened, the production cost is reduced, and the safe and efficient mining of the working face is ensured.Figure 14
Gas condition of return airway and upper corner in working face.


## 6. Conclusions

(1)
Through theoretical analysis and numerical simulation under the real geological parameters condition, it is obtained that during the mining period of panel 9303 in Anshun Coal Mine, when the working face advances to 110 m position, the working face reaches fully mining, and the height range of “three zones” of overlying strata is basically stable. At this time, the development range of fracture zone is 6–12.69 m, and this area is the best high drainage borehole final hole layer. The rationality of the development height of fracture zone obtained by theoretical analysis and numerical simulation is verified by the method of borehole peeping observation.(2)
By analyzing the field drainage effect of high directional long drilling, after the high directional long drilling in panel 9303 enters the stable and efficient extraction stage, the gas extraction rate is between 40% and 45%, and the pure gas extraction amount is up to 8.5 m3/min. The results show that the high directional long drilling has the advantages of gas concentration, large extraction amount, and stable flow, and the drilling drainage interception and extraction of adjacent layer gas is effective. It can provide reference for the arrangement of gas drainage boreholes in adjacent working faces and similar geological conditions.



---

*Source: 3370170-2021-08-18.xml*