Spatial andTemporalDistribution Law and Influencing Factors of the Mining-Induced Deformation and Failure of Gas Boreholes

Because gas boreholes are easy to damage by integrated coal mining and gas exploration, based on the practice of relieving pressure in deep thin coal seams in the Huainanmining area, a multidimensional coupling numerical simulationmethod was used to reveal the space-time evolution characteristics and influence factors of fracture deformation of gob-side gas boreholes. Results indicate that the danger zone for borehole fractures is primarily between 5 and 12m above the roof of the roadway. +e final-hole position has little effect on the stability of boreholes, and migrating the open-hole position to the entity coal side and roadway roof side can improve the stability of the borehole. +e initial failure of the borehole occurs at a distance of 10m behind the coal face. +e failure of the borehole is largely stable at a distance of 100 to 120m behind the coal face. With the increase in mining height, which leads to an increase in the movement of strata and an increase in pressure relief range, the shear stability of the borehole is reduced, and the extrusion stability of the borehole is improved. A hard roof condition promotes borehole shear stability, while a weak roof condition promotes borehole extrusion stability. +is change can decrease the maintenance difficulty associated with “minor supports” in boreholes to a certain extent by reinforcing the support strength of “primary supports” in roadway retaining walls. +e simulation results are consistent with observed results for the 11 test boreholes, and the accuracy of the numerical simulation is verified.


Introduction
Gas as a resource associated with coal is not only a type of clean energy but is also a dangerous source of energy that threatens safe mining practices in coal mines.How to extract gas safely and efficiently has always been a pertinent issue in coal mine safety in China [1].e key technology of integrated pillarless coal production and methane extraction is a new gas extraction technology that was developed in recent years on the basis of pressure relief mining technology.According to the occurrence conditions of coal seams, the key layer is first exploited to carry out pressure relief, and then the gas extraction roadways arranged in the top/bottom rock stratum of the key layer were replaced by the gob-side retaining roadway [2].e traditional U-shaped ventilation system was changed to a Y-type ventilation system, in which the pressure-relieved gas that was induced by the mining disturbance was continuously extracted by the drainage boreholes that were drilled behind the face in the gob-side retaining roadway.
e technology accounts for the safety and technical problems, including gas control, coal mining, roadway supports, ground temperature, and ground pressure control, and can realize safe and efficient integrated pillarless coal production and methane extraction [2].However, due to the movement of strata and adjustments to stress in the rock strata above and below the panel after mining, the inclined cross-measure boreholes drilled in the retaining roadway suffered severe mining effects such that the extraction boreholes tend to fracture, and the gas extraction rate is poor.Dense extraction boreholes must be used in order to meet the requirements of gas extraction.erefore, how to maximize the drainage effect and reduce the drilling volume and construction cost is an urgent need to further achieve safe and efficient integrated coal production and methane extraction [3].
To improve the effect of borehole drainage, studies on the permeability of rock masses during mining and the seepage law of pressure-relieved gas have been conducted by various researchers, who have obtained abundant research results [4][5][6][7][8][9][10].However, there has been less research focusing on the stability of boreholes, another factor that affects the effect of gas extraction, and existing literature primarily focuses on the mechanisms of surface venthole stability and protection technology.Whittles et al. investigated the effects that changes in the roof geology and roadway support system may have on deformation and closure of the gas drainage boreholes that are drilled from the tailgate across the gob of an active long wall panel.An analytical model was developed to estimate the bending deformation, axial strain, and the rupture of borehole casing due to rock shear and the optimum roadway support system, and the spacing between boreholes was also evaluated [11].Liu et al. [12] conducted a stability analysis of surface ventholes based on the theory of surface subsidence and rock mechanics, which were also used to determine the borehole and casing diameters and the properties of casing fill materials.Sun et al. [13][14][15] established a shear and tensile failure model through a mechanical analysis and three-dimensional numerical simulation method, and the impact of the radius and thickness of the well casing, cement ring, and key stratum on borehole stability was investigated and verified in the field test.Chen et al. [16] analyzed the deformation and failure of a surface venthole casing under compression, tension, and shear conditions using numerical simulations.Périé [17] studied the failure characteristics of the surrounding rock of a gas extraction well using a field test, and it was demonstrated that the intermediate principal stress has an important influence on the stability of the gas well wall.Liu et al. [18] proposed that the structural design of the drilling well should be based on the properties of the formation.It is best to maintain a certain distance between the casing and the shaft wall. is approach not only can provide space for vertical compression and bending deformation of gas extraction wells but can also buffer the damage when horizontal displacement or deformation of the strata occurs, and the screen section of the gas drainage well should be strengthened.
e above studies used a variety of methods to analyze the stability of the surface ventholes from different angles, but due to significant differences in the layouts of inclined cross-measure boreholes that are drilled in the retaining roadway and surface ventholes, it is necessary to study the space-time evolution characteristics and main influence factors of fracture deformation of the gob-side borehole according to its layout characteristics in order to provide references for the design and protection of underground drainage boreholes.

Engineering Background
e Zhuji coal mine is located in the Huainan Panxie mining area and has an annual production capacity of 4 Mt/a.e coal mine has many complicated geological characteristics, such as large changes in the loose layer thickness, it is deeply buried, has a high gas content, a high in situ gas pressure, and a high geotemperature.At present, the coal mine mainly exploits the 11-2 and 13-1 coal seams.e 13-1 coal seam is a coal and gas outburst coal seam with an average gas content of 6.98 m 3 /t, and gas content with such a high value can be seen in many literatures regarding gas drainage problems in the Huainan mining area, China [19,20].To eliminate the outburst danger of the 13-1 coal seam, the upper mining method was adopted to exploit the 11-2 coal seam firstly.e 1111(1) working panel is in the 11-2 coal seam.It has an average thickness of 1.26 m and an average dip angle of 3 °, ranging from 1 °to 5 °. e mining height for the coal is 1.8 m, and the face elevation is -877.6 to -907.0 m with an average mining depth of 910 m. e panel length is 1612 m, and the panel width is 220 m. e working seam is overlain by seam 13-1, which is an outburst coal seam with a thickness of 4 m at a distance of 68 m.As shown in Figure 1, the face was mined using retreat longwall (LW) mining with full seam extraction, and the authors have conducted laboratorial trials regarding strengthening borehole configuration under identical engineering background [21].It was ventilated using a Y-type ventilation system with the maingate and inlet part of the tailgate acting as the intake roadways and the outlet part of the tailgate and the floor return airgate acting as the return roadway.ere were seven connection roadways (CR1∼CR7) from the start line to the stop line between the tailgate and the floor return airgate.A 3 m wide concrete wall was constructed behind the LW face to replace the mined-out roadway wall along the gob side for building the retaining roadway.
During the mining of panel 1111 (1), in order to prevent pressure-relieved gas of overlying coal seam 13-1 from entering into the working face and to eliminate the outburst danger of coal seam 13-1, a group of two inclined crossmeasure boreholes were constructed at intervals of 25 m in the tailgate to extract the overburden pressure-relief gas.However, due to the redistribution of mining-induced stress and the overlying strata movement caused by the coal mining, deformation and failure associated with shear deformation, distortion, and extrusion appeared in the inclined cross-measure boreholes.
e single-hole gas extraction flow was less than 0.3 m 3 /min, and the effect of pressure-relief gas extraction could not be fully realized.
In view of the aforementioned engineering background, it is necessary to study the space-time evolution characteristics and main influence factors determining failure/deformation of the gob-side boreholes, and then some effective protection measures can be adopted to increase the stability of the borehole and extend the effective extraction time of the boreholes.Finally, the safe and efficient coal and gas simultaneous extraction can be realized.ree zones with different degrees of deformation were formed over the gob after mining (namely, the caved zone, the fractured 2

Geological Model and Simulation Program
Advances in Materials Science and Engineering zone, and the bending zone, as shown in Figure 1(a)), which were formed due to deformation of the overlying strata that occurs after the removal of panels.e deformation and failure modes of the gas drainage boreholes can be classified into two main categories: compressional failure resulting from concentrated compressive stresses both in the vertical and horizontal, and shear failure due to the dissimilarity of rock properties and incongruity of horizontal movement (shear offset) and longitudinal distortion between adjacent layers, as shown in Figure 2.

Rationale for Compressional Failure Mode.
For the compressional failure mode, set borehole and strata as an integral research object (Figure 3), and assume that both of them are in elastic state, then in situ stress along the cross-section direction of borehole is presented as follows based on elastic mechanics theory: where α and β are azimuth angle and inclination of the borehole, respectively; σ V , σ H , and σ h are vertical stress, maximum horizontal stress, and minimum horizontal stress of strata, respectively.e extrusion force applied on the borehole can be represented by the sum of even load and uneven load, that is: where P 1 � max σ x , σ y   and P 2 � min σ x , σ y  .Afterwards, set casing as research object, based on the Von Mises yield criterion and adopt inverse formulation in elastic state, and then formula for ultimate extrusion strength of casing can be calculated by where where σ s is yield strength of casing, K p is uneven coefficient of stress, and r and t are outer diameter and wall thickness of casing.
Based on the aforementioned theory, the safety factor of compressional failure of the borehole can be calculated according to (5).When the calculated value is less than 1, the borehole can be considered to be undergoing compressional failure: where f c is the compressional safety factor of the borehole.). 4 Advances in Materials Science and Engineering induced failure, and thus it can be roughly represented by sine function as [11] u(y where y is the distance along the borehole from shear plane (see Figure 4); u(y) is the displacement at location y perpendicular to the y axis; A is amplitude of displacement; and a is wavelength (twofold value of the width of the shearing zone), which is closely correlated with mechanical properties and stress environment of strata.
In view of bending deformation principle in mechanics of materials, the bending-moment equation for the borehole casing under shear-slip deformation state is based on following equation: where E is elasticity modulus of casing.en the shear stress on the cross section of an annular girder is expressed by where Q is shear force acted on the cross section of the girder; S * z is static moment versus neutral axis for a section beyond x on the cross section, as shown by the shaded area in Figure 5, ]; I z is inertia moment of cross section versus neutral axis; r and t are outer diameter and wall thickness of the annular girder, respectively, Equation (8) shows shear stress of the borehole casing under shear-slip deformation state, therefore, the maximum shear stress on the cross section of casing can be expressed by In the same measure as mentioned above for compressional case, the shear failure safety factor of the borehole casing can be calculated according to (10).When the calculated value is less than 1, the borehole can be considered as undergoing shear failure: where f s is the shear safety factor of the borehole, τ lim is the shear strength of the casing, and τ max is the maximum shear stress of the the casing section.Detailed calculations for each parameter can be found in the literature [22].

Multidimensional Coupled Numerical Simulation
Method.To obtain the characteristics of the macroscopic mining-induced stress distribution and stratal movement after coal mining, a 3D numerical model is usually adopted, note that there also have 2D FSDT and HSDT solutions besides 3D solutions.However, accounting for the accuracy of the simulation and the calculation of costs, as shown in Figure 6, the multidimensional coupled numerical simulation method was used to investigate the following.First, a large-scale 3D numerical model (larger grid size without considering the supporting effect) was constructed.e basic distribution law of the macroscopic mining-induced stress field and displacement field and the compaction rate of the gob at different places behind the coal face were obtained by 3D numerical simulation.Next, a 2D numerical model was used to further refine the grid and investigate the supporting effect.e distribution characteristics of the mining-induced  e thickness of the coal seam, the roof strata, and the floor strata of the panel were 2.0 m, 98 m, and 20 m, respectively.e distribution of strata is shown in Table 1.Because of the symmetry of stratal movement in the working face, the 1/2 working face width (110 m) was simulated, and the simulated mining distance was 150 m.During the process of mining, a gob-side entry retaining roadway was implemented.e size of the retaining roadway was 4.0 m wide and 3.0 m high.e dimension of the concrete wall was 3 m wide and 3 m high.e model used a nonuniform grid.Small grids were used near the top and bottom floors of the coal seam.Large grids were used away from the mining area.e minimum grid size was 1 m × 1 m × 2 m, and the maximum grid size was 3 m × 4 m × 2 m. e entire model was divided into 750,000 units and 776,286 nodes.e upper boundary of the model was the stress boundary, at which an equivalent load of 20 MPa was applied.
e remaining boundaries were displacement boundaries, and a normal displacement constraint was imposed.

Two-Dimensional Numerical
Model.Similar to the 3D numerical model, only the size of the Y direction was set to 0.4 m. e mesh density was encrypted based on the 3D model.e minimum mesh size of the encrypted 2D model was 0.25 m × 0.25 m × 0.1 m. e total 2D model was divided into 367,820 units and 464,420 nodes.e model boundary conditions were the same as those used in the 3D model.

Constitutive Model and Rock Mass Parameters.
Considering the simulation accuracy and computational cost, the 3D numerical model adopted the Mohr-Coulomb Model [23], whereas the 2D numerical model adopted the  6 Advances in Materials Science and Engineering Strain-Softening Ubiquitous-Joint Model, which could not only consider the stress-strain relationship of the material after peak but could also meet the needs of the analysis of joints and weak surfaces.Detailed rock mass parameters are shown in Table 2.

Simulation Program.
e coal seam excavation was excavated along the coal seam direction (Y direction) step-bystep in the 3D numerical simulation.Each excavation distance was 10 m, and a total of 150 m was excavated.After excavation, the monitoring lines were arranged along the Y direction at x � 200 and z � 29 to extract the compaction rate of the gob area at different distances behind the coal face, which provided an inversion basis for the 2D model to analyze the deformation and failure rules of drainage boreholes at different distances behind the coal face.
e coal seam excavation was excavated in one step in the 2D numerical simulation.To obtain the spatial distribution law of the deformation and failure of the drainage boreholes, as shown in Figure 8, a total of nine sets of drilling paths were designed.e initial location of the 1# through 3# drilling paths was x � 80 m and z � 32 m, and the final location was x � 110 m, 120 m, 130 m, and z � 97 m, respectively.e final location of the 4# through 9# drilling paths was x � 120 m and z � 97 m, and the initial location was x � 80 m, z � 37 m, 42 m, 47 m and x � 83 m, 86 m, 77 m, and z � 32 m, respectively.At the same time, according to the compaction rate of the gob area at different distances behind the coal face obtained from 3D simulation, the dynamic evolution rule of the drainage borehole deformation and failure behind the coal face was retrieved from the subsidence of the roof node of the caving zone.
To obtain the influence rule of different mining heights and roof properties on drainage borehole stability, three mining heights (2 m, 3 m, and 4 m) and three roof conditions (a common roof with 5 mudstone and 7 m sandstone, a hard roof with 12 m sandstone, and a soft roof with 12 m mudstone) were designed.
To obtain the influence law of different support conditions on the stability of the drainage boreholes, as shown in Table 3, a total of three types of supports were designed.
e cable unit in FLAC3D was used to simulate rock bolts and cables.e size of the bolt and cable was Φ22 × 2000 mm and Φ21.8 × 6000 mm, and the mechanical parameters of the bolt and cable are shown in Table 4. e concrete wall was simulated with three kinds of strength, and the strength parameters are shown in Table 5.
Considering the aforementioned factors, a total of seven sets of simulation schemes were designed, and the detailed simulation schemes are shown in Table 6.

Compaction Rate of Gob at Different Distances behind the Coal Face.
e subsidence ratio of the gob roof at different distances behind the coal face is shown in Figure 9. e transverse coordinates in Figure 9 indicate the distance from the coal face.A positive value represents an area behind the coal face, and a negative value represents the front of the coal face.Due to the fragmental dilatancy of caving rock in the gob, the compaction rates of the caving rock in the gob at different distances behind the coal face are different; therefore, the ratios of the subsidence of the gob roof are different.erefore, the ratio of subsidence of the gob roof can be used to characterize the compaction rate of the caving rock in the gob behind the coal face.It can be concluded from Figure 9 that the compaction rate of the caving rock at the position of the coal face is approximately 24%, and with an increase in the distance behind the coal face, the compacting rate of the caving rock in the gob also increases.e compacting rate of the caving rock is nearly 100% at a position 120 m behind the coal face.It can be concluded that the mining influence will tend to be stable behind the coal face at 120 m, which is consistent with the field measured value of panel 1111(1) in the Zhuji coal mine.erefore, the correctness of the numerical simulation is verified.

Spatial Distribution Law of Drainage Borehole Deformation and Failure.
Assuming that protective casing of the drainage boreholes is V89 * 8 mm N80 oil casing, the critical shear slip of the borehole casing can be calculated to be 0.01 m according to the literature [24].e displacement and stress of each node in the three-dimensional model can be extracted, and the shear slip and the compressional safety factors at different locations of each drilling path can be obtained.
e shear slip and compressional safety factor distribution curves of the 1# through 9# drilling paths are shown in Figure 10.e red dashed line represents the critical shear slip of the drainage boreholes, and the blue dashed line represents the critical compressional safety factor of the drainage boreholes.
Comparing the distribution of shear slip and compressional safety factors of the 1# through 3# drilling paths in Figure 10, the effect of the final location on the stability of the drainage boreholes can be obtained.e dangerous regions of shear and compressional failure of the drainage boreholes are both mainly located at a position near the final location of the boreholes.e dangerous regions of the shear failure of the drainage boreholes are mainly located in strata above the roadway roof between 1-12 m, 18-19 m, 21-27 m, and 30-32 m. e interlayer shear slip in the strata above the roadway roof at 1-5 m is the most dramatic with a maximum shear slip of 0.08 m, and the shear slip in the region above the roadway roof at 32 m is smaller.e dangerous regions of the compressional failure of the drainage boreholes are mainly located in strata above the roadway roof at 5-11 m.With the migration of the final location to the gob side, the dangerous regions and the risk degree of the compressional failure are gradually reduced.erefore, the final location of the drainage boreholes has little influence on borehole stability.e final location of the drainage boreholes can be migrated to the gob side in the drilling design so that the borehole can pass through the high gas penetration zone to achieve the best extraction effect.
Comparing the distribution of shear slip and compressional safety factors of the 4# through 6# drilling paths in Figure 10, the effect of the vertical distance between the  8 Advances in Materials Science and Engineering initial location of the boreholes with the roadway roof on the stability of the drainage boreholes can be obtained.Due to the upward shift of the initial location of the drainage boreholes, severe shearing and sliding in the strata above the roadway roof at 1-5 m is avoided.e shear slip of strata is significantly reduced with a maximum shear slip of 0.03 m, and the risk degree of borehole shear failure is significantly reduced.With an increase in the vertical distance between the initial location of the boreholes with the roadway roof, the dangerous area and the risk degree of the boreholes shear failure are reduced.Similarly, with an upward shift of the initial location of the drainage boreholes, severe compression in the strata above the roadway roof at 5-11 m is avoided.e dangerous area and the risk degree of borehole compressional failure are also reduced.

Advances in Materials Science and Engineering
Comparing the distribution of shear slip and compressional safety factors of the 7# through 9# drilling paths in Figure 10, the effect of the horizontal distance between the initial locations of the boreholes with the filling wall on the stability of the drainage boreholes can be obtained.With a decrease in the distance between the initial locations of boreholes with the filling wall, the dangerous area of the borehole shear failure is basically unchanged, but the risk degree of shear failure increases significantly.When the horizontal distance between the initial locations of boreholes with the filling wall is 10 m, the maximum shear slip of boreholes is only 0.026 m.When the horizontal distance between the initial location of boreholes with the filling wall is 4 m, the maximum shear slip of boreholes reaches 0.17 m, which is 6.5 times that of the former.erefore, the dangerous area of the borehole shear failure is mainly located in the 5 m of strata above the roadway roof, and the lower the layer, the greater the shear slip.As the horizontal distance between the initial locations of boreholes with the filling wall decreases, the danger area and risk degree of the borehole compression failure are slightly reduced, which indicates that the dangerous area of the borehole compressional failure is primarily located in the region near the sidewall on the nonmining side.

Temporal Distribution Law of the Deformation and
Failure of Boreholes.According to the subsidence ratio of the gob roof at different distances behind the coal face in the 3D model, the evolution rule of the dangerous area and risk degree of the borehole with the advance of the coal face can be retrieved by monitoring the subsidence value of the gob roof node in the 2D numerical model.e evolution law of the shear and compressional stability with advance of the coal face of the 2# drilling path is shown in Figure 11.From Figure 11, it can be concluded as follows: (1) e shear slip of the borehole increases approximately linearly with the increase in distance behind the coal face.e shear slip and growth rate of lower strata are significantly higher than those of upper strata, especially in the strata above the roadway roof at 5 m. e time of borehole shear failure initially occurred at a distance of 10 m behind the coal face, and with continuous advance of the coal face, the dangerous area and risk degree continued to increase until it was basically stable at a distance of 120 m behind the coal face.(2) e safety factor of borehole compressional failure gradually decreases with the advance of the coal face as a whole.e mining-induced stress adjustment of lower strata is less than that of upper strata, and the time to reach stability is also shorter than for upper strata.e mining-induced stress adjustment at a depth of 6 m along the borehole is basically stable at a distance of 60 m behind the coal face.e mining-induced stress adjustment at a depth of 60 m along the borehole still has greater adjustment at a distance of 60 m behind the coal face. is shows that the stress adjustment is lagging in the upper strata.e time of borehole compressional failure initially occurred at a distance of 10 m behind the coal face, and compressional failure occurred again in the upper strata at a distance of 70 m behind the coal face.e compressional failure of the borehole was basically stable at a distance of 100 m behind the coal face.

Effect of Mining Height on the Stability of Boreholes.
e variation rules of the borehole stability of the 2# drilling path with mining heights of 2 m, 3 m, and 4 m are shown in Figure 12.It is concluded from the diagram as follows: (1) With an increase in mining height, the dangerous area of the borehole shear failure also increases and is primarily concentrated in the upper strata at 20 m   is indicates that the increase in mining height mainly results in an increase in the displacement of upper strata, whereas lower strata are less affected by the mining height due to the supporting effect of the filling wall.
(2) As the mining height increases, the dangerous area and risk degree of borehole compressional failure reduced. is is due to the increase in mining height, which results in an increase in the pressure relief range that is caused by mining.

Effect of Roof Lithology on Borehole
Stability.e variation rules of the borehole stability of the 2# drilling path with the roof lithologies including a common roof, a hard roof, and a soft roof are shown in Figure 13.We conclude the following from the diagram: (1) e hard roof conditions are most favorable to shear stability of the borehole, and the soft roof conditions are the second most favorable.When the roof lithology is a hard roof, the subsidence of the lateral roof obviously decreases due to the high strength and stiffness of the roof, especially in lower strata, such that the shear slip of a hard roof is lower than that of other types of roofs.When the roof lithology is a soft roof, stratal movement is lower than for the common roof condition due to the dilatancy of the soft roof, and the shear slip of the soft roof is lower than that of the common roof.

Advances in Materials Science and Engineering
(2) e soft roof conditions are most favorable to compressional stability of the borehole, and the hard roof is the second most favorable.When the roof lithology is a hard roof, the concentration of horizontal stress in the lateral lower strata of the gob is higher than that of the common roof, such that it is easy to form a dangerous area of compressional failure in the lower strata.When the roof lithology is a soft roof, the mining disturbance effect is lower than for the common roof, and the dangerous area of the compressional failure is significantly reduced.12 Advances in Materials Science and Engineering 4.6.Effect of Support Scheme on the Stability of Borehole.e variation rules of the borehole stability of the 2# drilling path with support schemes 1, 2, and 3 are shown in Figure 14.We conclude the following from the diagram: (1) It can reduce the dangerous area and risk degree of borehole shear failure to strengthen the support strength of the retaining roadway, especially in lower strata.is is mainly because the subsidence of the roadway roof on the filling wall side is controlled by improving the supporting strength.erefore, the interlayer shear slip is reduced.It can be concluded that the reinforcing support strength of "primary supports" in roadway retaining walls can decrease the difficulty associated with maintenance of "minor supports" of boreholes, to some extent.
(2) Strengthening the support strength also has a certain effect on the stability of the borehole compressional failure, but the effect is not significant with regard to the stability of borehole shear failure.

Field Test
As shown in Figure 1, we drilled a total of 11 drainage boreholes in three groups, sequentially, in CR 4, CR 3, and CR 2. Different protection measures were adopted in the drainage boreholes.To grasp the law of deformation and failure of drainage boreholes, a borehole camera exploration device was used to measure borehole stability during the advance of the face.After five months of monitoring, boreholes 1 through 8 were found to be broken, whereas boreholes 9 through 11, which had a dual thick-walled combined casing, did not fracture during the extraction.Details of the borehole failure position along the casing, the height into the roof, and the distance behind the coal face where the casing failed are shown in Figure 15.
From Figure 15, we conclude that boreholes were broken at a distance of 2-14 m behind the coal face, and 87.5% of that occurred over a distance of 9-14 m behind the coal face.e failure position of boreholes occurred at depths of 5-10 m along the casing, which were all in strata above the roadway roof at 4-9 m. e observed results are consistent with the numerical simulation, indicating that the proposed multidimensional coupling numerical simulation method can be used as an effective method for borehole design and protection.

Conclusions
Gas boreholes that penetrate mineable coal seams may be subject to distress caused by integrated mining coal and gas.
e concentrated compressive stresses and the lateral shear offsets are especially important because they are the most damaging for boreholes.To reveal the space-time evolution characteristics and influence factors of fracture deformation of the gob-side gas boreholes, a multidimensional coupling numerical model that considers the combined effects of topography, various mining heights, and support schemes is presented based on the practice of deep thin coal seam pressure relief in the Huainan mining area.rough this study, the following conclusions are drawn:  Advances in Materials Science and Engineering (2) With an increase in mining height leading to an increase in the movement of strata and an increase in the range of pressure relief, the shear stability of the borehole is reduced and the extrusion stability of the borehole is improved.A hard roof condition promotes borehole shear stability; a weak roof condition promotes borehole extrusion stability.It can decrease the maintenance difficulty associated with "minor supports" in boreholes to some extent by reinforcing the support strength of "primary supports" in roadway-retaining walls.ese results are consistent with the results of the numerical simulation analysis.

Figure 1 :
Figure 1: Panel 1111(1) and testing boreholes.(a) Plan view.(b) A-A section.e figure was reproduced from Xue et al. [21] (under the Creative Commons Attribution License/public domain).

Figure 3 :
Figure 3: Sketch for mechanical model of casing under compression in strata.

3. 3 .
Model Design and Rock Parameters 3.3.1.ree-Dimensional Numerical Model.To eliminate the boundary effect and combine the engineering geology and mining technical conditions, as shown in Figure 7, the 3D model was scaled to 200 m wide, 200 m long, and 120 m high.

Figure 6 :Figure 7 :
Figure 6: Schematic diagram for the multidimensional coupling numerical simulation method.

Figure 8 :
Figure 8: Schematic diagram for the location of the drainage boreholes in the 2D model.Coordinates attached on the right section indicate initial location and final location for each drilling path.

Figure 9 :
Figure 9: Roof subsidence curve at different distances behind the coal face.

Figure 11 :Figure 12 :
Figure 11: Evolution laws for the stability of the borehole with advance of the coal face.(a) Evolution of the dangerous location and shear slip of the borehole with advance of the coal face.(b) Evolution of the dangerous location and safety factor of borehole compressional failure with advance of the coal face.

( 1 )
Results indicate that the danger zone for borehole fracturing is mainly in the scope of 5-12 m above the roof of the roadway.e final-hole position has little effect on the stability of boreholes, and migrating the open-hole position to the entity coal side and roadway roof side can improve the stability of the borehole.e initial failure of the borehole occurs at a distance of 10 m behind the coal face.e failure of the borehole is basically stable at a distance of 100-120 m behind the coal face.

Figure 13 :
Figure 13: Effect of roof lithology on the stability of the borehole.(a) Effect of roof lithology on the shear stability of the borehole.(b) Effect of roof lithology on the compressional stability of the borehole.
(3) e tracing observation of the borehole stability during advance of the 11 test boreholes shows that borehole destruction occurs at a distance of 9-14 m behind the coal face, and the failure position of boreholes occurs at depths of 5-10 m along the casing.

Figure 15 :Figure 14 :
Figure 15: Data for borehole failure position along the borehole, height into the roof, and distance behind the coal face where the borehole failed.

Table 1 :
Distribution of strata and their parameters in the basic model.

Table 3 :
ree different support schemes.

Table 4 :
Material parameters for bolts and cables.

Table 5 :
Parameters for backfilling wall materials.

Table 2 :
Mechanical parameters of strata.