Method and apparatus of controlling drilling for rock burst prevention in coal mine roadway

Abstract

A method for controlling drilling for rock burst prevention drilling in a coal mine roadway is provided. The method comprises: acquiring rock mechanical parameters of coal mass in surrounding rock of a roadway to be subjected to burst-preventing drilling construction, and calculating a surrounding rock critical softening depth, a critical ground stress and a critical mining peak stress for rock burst initiation in the roadway; calculating a critical mining-induced stress index of the roadway to realize quantification of burst risk; then determining critical conditions for drillhole burst and a quantitative relationship between the critical conditions for drillhole burst and for roadway rock burst initiation; quantitatively determining construction parameters of burst-preventing drillholes according to the surrounding rock critical softening depth, a critical plastic softening zone radius for drillhole burst, and the critical mining-induced stress index; and controlling a drilling machine to operate according to the determined construction parameters.

Claims

1. A method for controlling drilling for rock burst prevention in a coal mine roadway, comprising the following steps: S1: acquiring rock mechanical parameters of coal mass in surrounding rock of a roadway to be subjected to burst-preventing drilling construction; S2: calculating a surrounding rock critical softening depth L.sub.pcr, a critical ground stress P.sub.cr and a critical mining peak stress P.sub.mcr of a surrounding rock stress concentration zone for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction; S3: acquiring a mining peak stress P.sub.m of the coal mass in the surrounding rock of the roadway to be subjected to burst-preventing drilling construction, optimizing the critical mining peak stress P.sub.mcr of the surrounding rock stress concentration zone for roadway rock burst initiation, and calculating a critical mining-induced stress index K.sub.cr of the roadway to be subjected to burst-preventing drilling construction; $\begin{array}{cc}{K}_{\mathrm{cr}}=\frac{P_m}{P_{\mathrm{mcr}}^{\star }}& \left(1\right)\end{array}$ wherein P.sub.mcr* is the optimized critical mining peak stress of the surrounding rock stress concentration zone for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction; S4: determining critical conditions for drillhole burst occurrence; calculating a critical fracture zone radius r.sub.dcr, a critical plastic softening zone radius r.sub.pcr and a critical stress P.sub.hcr for drillhole burst occurrence; S5: determining a relationship between the critical conditions for drillhole burst occurrence and critical conditions for roadway rock burst initiation; S6: quantitatively determining a drillhole diameter, a drillhole depth L.sub.drill and drillhole spacing D.sub.drill of burst-preventing drillholes according to the surrounding rock critical softening depth L.sub.pcr for roadway rock burst initiation, the critical plastic softening zone radius r.sub.pcr for drillhole burst occurrence and the critical mining-induced stress index K.sub.cr; and S7: controlling a drilling machine to operate according to the determined drillhole diameter, drillhole depth L.sub.drill and drillhole spacing D.sub.drill of the burst-preventing drillholes.

2. The method according to claim 1, wherein the rock mechanical parameters in S1 comprise uniaxial compressive strength σ.sub.c, elastic modulus E, a burst modulus index K=λ.sub.1/E, residual modulus reduction λ.sub.2, and a residual strength coefficient ξ, wherein λ.sub.1 is post-peak softening modulus.

3. The method according to claim 2, wherein the critical fracture zone radius r.sub.dcr, the critical plastic softening zone radius r.sub.pcr and the critical stress P.sub.hcr for drillhole burst occurrence calculated in S4 are as shown in the following formula: $\begin{array}{cc}{r}_{\mathrm{dcr}}=r_0\sqrt{\frac{\alpha }{\beta }}& \left(2\right)\end{array}$ $\begin{array}{cc}{r}_{\mathrm{pcr}}=r_0\sqrt{\frac{\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{E}{{\lambda }_1}+1}& \left(3\right)\end{array}$ $\begin{array}{cc}{P}_{\mathrm{hcr}}=\frac{m+1}{2}{\sigma }_c\left(\frac{p_{\mathrm{dcr}}}{{\sigma }_c}+\frac{1+{\lambda }_1/E}{m-1}\right){\left(\left(1-\xi \right)\frac{1}{K}+1\right)}^{\frac{m-1}{2}}-\frac{{\sigma }_c{\lambda }_1/E}{2}{\left(\left(1-\xi \right)\frac{1}{K}+1\right)}^{\frac{m+1}{2}}-\frac{1+{\lambda }_1/E}{m-1}{\sigma }_c& \left(4\right)\end{array}$ wherein $p_{\mathrm{dcr}}=\left(\frac{\alpha }{1-q}\right)\left[1-{\left(\frac{r_{\mathrm{dcr}}}{r_0}\right)}^{q-1}\right]+\left(\frac{\beta }{1+q}\right)\left[1-{\left(\frac{r_{\mathrm{dcr}}}{r_0}\right)}^{q+1}\right]$  is an acting stress of a surrounding rock fracture zone on a plastic softening zone when a drillhole burst occurs; $\alpha ={\sigma }_c\left[\frac{{\lambda }_2}{E}+\frac{{\lambda }_2}{{\lambda }_1}\left(1-\xi \right)+\xi \right],\beta ={\sigma }_c\left[\frac{{\lambda }_2}{E}+\frac{{\lambda }_2}{{\lambda }_1}\left(1-\xi \right)\right],m=\frac{1+\sin \phi }{1-\sin \phi },$  φ is an internal friction angle of a coal rock medium in the plastic softening zone of the roadway surrounding rock; $q=\frac{1+\sin {\varphi }^{\prime }}{1-\sin {\varphi }^{\prime }},$  φ′ is an internal friction angle of the coal rock medium in the fracture zone of the roadway surrounding rock; and r.sub.0 is a drillhole radius or a drill bit cutting radius.

4. The method according to claim 3, wherein the relationship between the critical conditions for drillhole burst occurrence and the critical conditions for roadway rock burst initiation determined in S5 meets the following relational expression:
P.sub.mcr*>P.sub.cr>P.sub.hcr  (5).

5. The method according to claim 4, wherein the drillhole diameter in S6 is determined according to the arrangement mode of the burst-preventing drillholes in the rock burst roadway and the self-conditions of the mine; the drillhole depth L.sub.drill is determined based on the surrounding rock critical softening depth L.sub.pcr for roadway rock burst initiation, and is as shown in the following formula:
L.sub.drill=η.sub.dη.sub.LL.sub.pcr  (6) wherein η.sub.d is a correction coefficient for coal seam thickness; when the coal seam thickness is greater than 0 m and less than 4 m, the value range of η.sub.d is 0.8≤η.sub.d≤0.9, when the coal seam thickness is 4-8 m, the value range of η.sub.d is 0.9<η.sub.d≤10; when the coal seam thickness is greater than 8 m, the value range of η.sub.d is 1.0<η.sub.d≤1.2; η.sub.L; is a burst-preventing safety coefficient for the drillhole depth; two determination methods are provided for η.sub.L: one is to determine η.sub.L according to the critical mining-induced stress index K.sub.cr of burst risk evaluation, namely η.sub.L=0.85+0.5K.sub.cr; and the other method is to determine η.sub.L according to a burst risk level obtained by burst risk evaluation based on a comprehensive index method; the drillhole spacing D.sub.drill is determined based on the critical plastic softening zone radius r.sub.pcr for drillhole burst occurrence, and is as shown in the following equation:
D.sub.drill=2η.sub.pcrr.sub.pcr  (7) formula (3) is combined with formula (7) to further obtain $\begin{array}{cc}{D}_{\mathrm{drill}}={\eta }_{\mathrm{pcr}}d\sqrt{\frac{\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{1}{K}+1}& \left(8\right)\end{array}$ wherein η.sub.pcr is a burst-preventing safety coefficient for the burst-preventing drillhole spacing, d is the diameter of a drill bit in burst-preventing drilling construction, d=2r.sub.0; two determination methods are provided for η.sub.pcr: one is to determine η.sub.pcr according to a critical stress index method of burst risk evaluation, namely η.sub.pcr=2.325−1.75K.sub.cr; and the other method is to determine η.sub.pcr according to the burst risk level obtained by the burst risk evaluation method based on the comprehensive index method.

6. An apparatus for controlling drilling for rock burst prevention in a coal mine roadway, the apparatus comprising a memory and a processor that is configured to perform the method according to claim 1.

Description

BRIEF DESCRIPTION OF DRAWINGS

(1) FIG. 1 is a flowchart of a method for controlling drilling for rock burst prevention in a coal mine roadway provided by an embodiment of the present invention;

(2) FIG. 2 is a schematic diagram of a mechanical model of roadway rock burst initiation provided by an embodiment of the present invention;

(3) FIG. 3 is a schematic diagram of a mechanical model of drillhole burst occurrence provided by an embodiment of the present invention;

(4) FIG. 4 is a schematic diagram of the relationship between drillhole burst occurrence and roadway rock burst initiation provided by an embodiment of the present invention;

(5) FIG. 5 is a flowchart of design of drilling burst-preventing key parameters provided by an embodiment of the present invention;

(6) FIG. 6 is a schematic diagram of the influence of a drillhole diameter on the pressure relief effect in the thickness direction of coal mass provided by an embodiment of the present invention;

(7) FIG. 7 is a schematic diagram of a drilling depth of drillholes for preventing and controlling rock burst initiation in a roadway provided by an embodiment of the present invention; and

(8) FIG. 8 is a comparison result diagram of the amount of drilling cuttings per meter of a typical drillhole before and after burst-preventing drilling construction in a roadway provided by an embodiment of the present invention.

(9) In the figures: 1, mining-induced stress concentration zone; 2, disturbance stress wave; and 3, drillhole.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(10) The specific implementation of the present invention will be further described in detail below in combination with the accompanying drawings and embodiments. The embodiments below are adopted to illustrate the present invention, but not to limit the scope of the present invention.

(11) This embodiment takes the main 5# coal seam of a mine in Hebei as an example, and by the adoption of a method for controlling drilling for rock burst prevention in a coal mine roadway provided by the invention, drilling parameters of burst-preventing drillholes in a rock burst roadway of the coal seam can be determined, and a drilling machine is controlled to operate according to the drilling parameters.

(12) A burst-preventing drilling parameter determination method for a rock burst roadway in a coal mine, as shown in FIG. 1, comprises the following steps:

(13) S1, acquiring rock mechanical parameters of coal mass in surrounding rock of the roadway to be subjected to burst-preventing drilling construction, the rock mechanical parameters comprising uniaxial compressive strength σ.sub.c, elastic modulus E, a burst modulus index K=λ.sub.1/E, residual modulus reduction λ.sub.2, and a residual strength coefficient ξ, wherein λ.sub.1 is post-peak softening modulus;

(14) in this embodiment, the coal seam has an average thickness of 7.03 m, an inclination angle of 13 degrees, and an average buried depth of 984 m. The average uniaxial compressive strength of the coal mass is 10 MPa. The coal seam has a weak burst tendency, a roof has a weak burst tendency, and a floor has no burst tendency. The physical parameters of coal rock, support strength and geometric characteristic parameters of the roadway are shown in Table 1 for details;

(15) S2: calculating a surrounding rock critical softening depth L.sub.pcr, a critical ground stress P.sub.cr and a critical mining peak stress P.sub.mcr of a surrounding rock stress concentration zone for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction;

(16) S2.1: acquiring a support stress p.sub.s of the roadway to be subjected to burst-preventing drilling construction;

(17) S2.2: calculating a critical fracture zone radius ρ.sub.fcr and a critical softening zone radius ρ.sub.cr for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction, as shown in the following formulas:

(18) $\begin{array}{cc}{\rho }_{\mathrm{fcr}}={\rho }_0\sqrt{\frac{p_s\left(q-1\right)+\alpha }{\beta }}& \left(1\right)\end{array}$$\begin{array}{cc}{\rho }_{\mathrm{cr}}={\rho }_0\sqrt{\frac{p_s\left(q-1\right)+\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{E}{{\lambda }_1}+1}& \left(2\right)\end{array}$

(19) wherein ρ.sub.0 is a roadway radius after the roadway to be subjected to burst-preventing drilling construction is equivalent to a homogeneous, continuous and isotropic circular roadway;

(20) $q=\frac{1+\sin {\varphi }^{\prime }}{1-\sin {\varphi }^{\prime }},$
φ′ is an internal friction angle of a coal rock medium in a fracture zone of the roadway surrounding rock,

(21) $\alpha ={\sigma }_c\left[\frac{{\lambda }_2}{E}+\frac{{\lambda }_2}{{\lambda }_1}\left(1-\xi \right)+\xi \right],\beta ={\sigma }_c\left[\frac{{\lambda }_2}{E}+\frac{{\lambda }_2}{{\lambda }_1}\left(1-\xi \right)\right];$

(22) S2.3: calculating the surrounding rock critical softening depth L.sub.pcr and the critical ground stress P.sub.cr for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction, as shown in the following formulas:

(23) 0$\begin{array}{cc}{L}_{\mathrm{pcr}}={\rho }_0\sqrt{\frac{p_s\left(q-1\right)+\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{1}{K}+1}-\frac{B}{2}& \left(3\right)\end{array}$$\begin{array}{cc}{P}_{\mathrm{cr}}={\sigma }_c\left\{\frac{m+1}{2}{\left(\frac{p_{\mathrm{fcr}}}{{\sigma }_c}+\frac{1+{\lambda }_1/E}{m-1}\right)\left[1+\left(1-\xi \right)\frac{E}{{\lambda }_1}\right]}^{\frac{m-1}{2}}-{\frac{{\lambda }_1/E}{2}\left[1+\left(1-\xi \right)\frac{E}{{\lambda }_1}\right]}^{\frac{m+1}{2}}-\frac{1+{\lambda }_1/E}{m-1}\right\}& \left(4\right)\end{array}$

(24) wherein B is the width of the roadway to be subjected to burst-preventing drilling construction,

(25) $m=\frac{1+\sin \phi }{1-\sin \phi },$
φ is an internal friction angle of the coal rock medium in a plastic softening zone of the roadway surrounding rock, and p.sub.fcr is an acting stress of the surrounding rock fracture zone on the plastic softening zone when a rock burst is initiated in the roadway to be subjected to burst-preventing drilling construction, as shown in the following formula:

(26) $\begin{array}{cc}{p}_{\mathrm{fcr}}={p_s\left(\frac{{\rho }_{\mathrm{fcr}}}{{\rho }_0}\right)}^{q-1}+\left(\frac{\alpha }{1-q}\right)\left[1-{\left(\frac{{\rho }_{\mathrm{fcr}}}{{\rho }_0}\right)}^{q-1}\right]+\left(\frac{\beta }{1+q}\right)\left[1-{\left(\frac{{\rho }_{\mathrm{fcr}}}{{\rho }_0}\right)}^{q+1}\right]& \left(5\right)\end{array}$

(27) S2.4: calculating the critical mining peak stress P.sub.mcr of the surrounding rock stress concentration zone for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction, as shown in the following formula:

(28) $\begin{array}{cc}{P}_{\mathrm{mcr}}=2P_{\mathrm{cr}}-\frac{2P_{\mathrm{cr}}-{\sigma }_c}{m+1}& \left(6\right)\end{array}$

(29) S3: acquiring a mining peak stress P.sub.m of the coal mass in the surrounding rock of the roadway to be subjected to burst-preventing drilling construction, optimizing the critical mining peak stress P.sub.mcr of the surrounding rock stress concentration zone for roadway rock burst initiation, and calculating a critical mining-induced stress index K.sub.cr of the roadway to be subjected to burst-preventing drilling construction, and thus achieving the quantification of burst risk, wherein the critical mining-induced stress index K.sub.cr quantitatively represents the possibility degree of rock burst occurrence in the roadway to be subjected to burst-preventing drilling construction;

(30) firstly, acquiring the mining peak stress P.sub.m of the coal mass in the surrounding rock of the roadway to be subjected to burst-preventing drilling construction, and optimizing the critical mining peak stress of the surrounding rock stress concentration zone for rock burst initiation in the roadway to be subjected to burst-preventing drilling construction according to the section shape of the roadway to be subjected to burst-preventing drilling construction to be P.sub.mcr*=n.sub.1×P.sub.mcr, wherein n.sub.1 is a correction coefficient for the section of the roadway to be subjected to burst-preventing drilling construction; when the section of the roadway to be subjected to burst-preventing drilling construction is rectangular, trapezoidal, straight wall arched or circular in shape, n.sub.1 is 0.89, 0.92, 0.95 and 0.98, respectively;

(31) then calculating the critical mining-induced stress index K.sub.cr of burst risk in the roadway to be subjected to burst-preventing drilling construction, as shown in the following formula:

(32) $\begin{array}{cc}{K}_{\mathrm{cr}}=\frac{P_m}{P_{\mathrm{mcr}}^{\star }}& \left(7\right)\end{array}$

(33) S4: determining critical conditions for drillhole burst occurrence;

(34) S4.1: obtaining a drilled surrounding rock system equation from radial stress continuing conditions of each sub-zone of the drilled surrounding rock by combining the Mohr-Coulomb yield criterion and boundary conditions σ.sub.r(r.sub.0)=0 of a radial stress of the surrounding rock at the drillhole wall according to a coal rock equilibrium differential equation, a geometric equation, a constitutive equation and a coal rock damage evolution equation under uniaxial compression from comparison between a mechanical model of the roadway rock bursts shown in FIG. 2 and a mechanical model of drillhole burst occurrence shown in FIG. 3;

(35) the coal rock equilibrium differential equation:

(36) $\begin{array}{cc}\frac{d{\sigma }_r}{\mathrm{dr}}-\frac{{\sigma }_{\theta }-{\sigma }_r}{r}=0& \left(8\right)\end{array}$

(37) the geometric equation:

(38) $\begin{array}{cc}\begin{array}{c}{\epsilon }_r=\frac{\mathrm{du}}{\mathrm{dr}}\\ {\epsilon }_{\theta }=\frac{u}{r}\end{array}\}& \left(9\right)\end{array}$

(39) wherein r is the radius of the drilled surrounding rock, and while taking different values, r represents different positions of the drilled surrounding rock; ε.sub.r is radial strain of an elastic zone of the drilled surrounding rock, ε.sub.θ is circumferential strain of the elastic zone of the drilled surrounding rock, and u is radial displacement of the drilled) surrounding rock; σ.sub.r(r.sub.0) is the radial stress of the surrounding rock at the drillhole wall, r.sub.θ is a drillhole radius or a drill bit cutting radius, and σ.sub.θ σ.sub.r are a tangential stress of the elastic zone of the drilled surrounding rock and the radial stress of the surrounding rock respectively;

(40) the constitutive equation:

(41) the constitutive relation in the elastic zone of the drilled surrounding rock meeting:

(42) $\begin{array}{cc}a.\begin{array}{c}{\sigma }_r=\stackrel{\_}{E}\left({\epsilon }_r+{\stackrel{\_}{v\epsilon }}_{\theta }\right)\\ {\sigma }_{\theta }=\stackrel{\_}{E}\left({\epsilon }_{\theta }+{\stackrel{\_}{v\epsilon }}_r\right)\end{array}\}& \left(10\right)\end{array}$

(43) wherein

(44) $\stackrel{\_}{E}=\frac{E\left(1-v\right)}{\left(1+v\right)\left(1-2v\right)},\stackrel{\_}{v}=\frac{v}{\left(1-v\right)},$
ν is a Poisson's ratio;

(45) (2) the constitutive relation in the plastic softening zone of the drilled surrounding rock meeting:

(46) $\begin{array}{cc}\frac{{\sigma }_{\theta }}{1-D}=m\frac{{\sigma }_r}{1-D}+{\sigma }_c& \left(11\right)\end{array}$

(47) (3) the constitutive relation in the fracture zone of the drilled surrounding rock meeting:

(48) 0$\begin{array}{cc}\frac{{\sigma }_{\theta }}{1-D}=q\frac{{\sigma }_r}{1-D}+{\sigma }_c& \left(12\right)\end{array}$

(49) the coal rock damage evolution equation:

(50) $\begin{array}{cc}\begin{array}{cc}D=1-\left(1-\frac{r_d^2}{r^2}\right)\gamma -\frac{\xi r_d^2}{r^2}& \left(r<r_d\right)\\ D=\frac{{\lambda }_1}{E}\left(\frac{r_p^2}{r^2}-1\right)& \left(r_d<r<r_p\right)\\ D=0& \left(r>r_p\right)\end{array}\}& \left(13\right)\end{array}$

(51) wherein D is a damage variable of the coal rock medium in the drilled surrounding rock, γ=λ.sub.2/E+(1−ξ)λ.sub.2/λ.sub.1+ξ, r.sub.d is the radius of the fracture zone of the drilled surrounding rock, and r.sub.p is the radius of the plastic softening zone of the drilled surrounding rock;

(52) obtaining the drilled surrounding rock system equation from the radial stress continuing conditions of each sub-zone of the drilled surrounding rock by combining the Mohr-Coulomb yield criterion and boundary conditions σ.sub.r(r.sub.0)=0 of the radial stress of the surrounding rock at the drillhole wall, as shown in the following formula:

(53) $\begin{array}{cc}{P}_h=\frac{m+1}{2}{\sigma }_c\left(\frac{p_d}{{\sigma }_c}+\frac{1+{\lambda }_1/E}{m-1}\right){\left(\left(1-\xi \right)\frac{1}{K}+1\right)}^{\frac{m-1}{2}}-\frac{{\sigma }_c{\lambda }_1/E}{2}{\left(\left(1-\xi \right)\frac{1}{K}+1\right)}^{\frac{m+1}{2}}-\frac{1+{\lambda }_1/E}{m-1}{\sigma }_c& \left(14\right)\end{array}$

(54) wherein, P.sub.h is an stress of the drilled surrounding rock, namely a roadway mining-induced stress, and

(55) $p_d=\left(\frac{\alpha }{1-q}\right)\left[1-{\left(\frac{r_d}{r_0}\right)}^{q-1}\right]+\left(\frac{\beta }{1+q}\right)\left[1-{\left(\frac{r_d}{r_0}\right)}^{q+1}\right]$
is the acting stress of the drilled surrounding rock fracture zone on the plastic softening zone;

(56) S4.2: obtaining the critical fracture zone radius r.sub.dcr, the critical plastic softening zone radius r.sub.pcr and the critical stress P.sub.hcr for drillhole burst occurrence according to a disturbance response instability criterion

(57) $\frac{\mathrm{dr}}{{\mathrm{dP}}_h}=\infty$
for burst initiation, as shown in the following formulas:

(58) $\begin{array}{cc}{r}_{\mathrm{dcr}}=r_0\sqrt{\frac{\alpha }{\beta }}& \left(16\right)\end{array}$

(59) $\begin{array}{cc}{r}_{\mathrm{pcr}}=r_0\sqrt{\frac{\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{E}{{\lambda }_1}+1}& \left(17\right)\end{array}$$\begin{array}{cc}{P}_{\mathrm{hcr}}=\frac{m+1}{2}{\sigma }_c\left(\frac{p_{\mathrm{dcr}}}{{\sigma }_c}+\frac{1+{\lambda }_1/E}{m-1}\right){\left(\left(1-\xi \right)\frac{1}{K}+1\right)}^{\frac{m-1}{2}}-\frac{{\sigma }_c{\lambda }_1/E}{2}{\left(\left(1-\xi \right)\frac{1}{K}+1\right)}^{\frac{m+1}{2}}-\frac{1+{\lambda }_1/E}{m-1}{\sigma }_c& \left(18\right)\end{array}$

(60) wherein

(61) $p_{\mathrm{dcr}}=\left(\frac{\alpha }{1-q}\right)\left[1-{\left(\frac{r_{\mathrm{dcr}}}{r_0}\right)}^{q-1}\right]+\left(\frac{\beta }{1+q}\right)\left[1-{\left(\frac{r_{\mathrm{dcr}}}{r_0}\right)}^{q+1}\right]$
is the acting stress of the surrounding rock fracture zone on the plastic softening zone when the drillhole burst occurs;

(62) S5: determining a relationship between the critical conditions for drillhole burst occurrence and critical conditions for roadway rock burst initiation;

(63) drillhole burst occurrence and roadway rock burst initiation have the same occurrence mechanism, that is, under the high stress condition, coal rock in a softening zone and coal rock in an elastic zone of the roadway or drilled surrounding rock form an unstable balance system, the boundary of the surrounding rock plastic zone generates great nonlinear expansion under external disturbance, and a series of macroscopic responses are triggered. However, for surrounding rock burst-preventing drillholes of a specific roadway, the axial direction of the roadway is perpendicular to the axial direction of the drillhole, as shown in FIG. 4. From the analysis of FIG. 4, it can be seen that the reason for roadway rock burst initiation is that the mining-induced stress of the roadway reaches the critical mining peak stress P.sub.mcr* for roadway rock burst initiation, and the reason for drillhole burst occurrence is that the mining-induced stress of the roadway reaches the critical stress P.sub.hcr for drillhole burst occurrence. Therefore, the relationship between drillhole burst occurrence and roadway rock burst initiation is specifically embodied in the following aspects: 1) both drillhole burst occurrence and roadway rock burst initiation have the same disturbance response instability mechanism, that is, the drillhole can be regarded as a circular roadway without support stress; 2) for the specific roadway and the surrounding rock drillholes thereof, the surrounding rock has the same physical and mechanical parameters; 3) in a spatial position, the axis of the roadway is perpendicular to the axes of the drillholes; and 4) driving stress sources for drillhole burst occurrence and roadway rock burst initiation are the same, that is, both stresses are roadway mining concentrated stresses.

(64) In conclusion, according to the critical stress P.sub.hcr for drillhole burst occurrence and the optimized critical mining peak stress P.sub.mcr* for roadway rock burst initiation, the relationship between the critical conditions for drillhole burst and the critical conditions for roadway rock burst initiation is determined to meet the following relational expression:
P*.sub.mcr>P.sub.cr>P.sub.hcr  (19)

(65) From the relationship between the critical conditions for the drillhole burst and the critical conditions for roadway rock burst initiation shown in equation (19), it can be seen that under the driving of the certain mining-induced stress, the critical stress for drillhole burst occurrence are less than the critical stress for roadway rock burst initiation, that is, drillhole burst occurrence is easier than roadway rock burst initiation, which reveals the phenomenon that due to drilling, the drillhole burst occurs but the roadway burst is not initiated in engineering. Therefore, once the critical conditions for drillhole burst occurrence are destroyed, the critical conditions for roadway rock burst initiation can be destroyed, and thus the rock burst is prevented and controlled. Therefore, a quantitative theoretical basis is provided for determining drilling construction parameters for the purpose of preventing and controlling roadway rock burst initiation.

(66) S6: quantitatively determining construction parameters of burst-preventing drillholes according to the surrounding rock critical softening depth L.sub.pcr for roadway rock burst initiation, the critical plastic softening zone radius r.sub.pcr for drillhole burst occurrence, and a roadway burst risk characterization parameter, namely, the critical mining-induced stress index K.sub.cr.

(67) S7: controlling a drilling machine to operate according to the determined drillhole diameter, drillhole depth L.sub.drill and drillhole spacing D.sub.drill of the burst-preventing drillholes.

(68) In the method provided by the present invention, the construction design principle of burst-preventing drilling is as follows:

(69) (1) the critical plastic softening zone radius r.sub.pcr for drillhole burst occurrence, a surrounding rock critical softening depth L.sub.pcr for roadway rock burst initiation and the critical mining-induced stress index K.sub.cr of roadway burst risk are taken as data bases for quantitatively determining the burst-preventing drillhole construction parameters;

(70) (2) the relationship between the critical conditions for drillhole burst occurrence and the critical conditions for roadway rock burst initiation is taken as a theoretical basis for quantitatively determining the burst-preventing drillhole construction parameters;

(71) (3) in the aspect of drillhole depth, the surrounding rock critical softening depth L.sub.pcr for roadway rock burst initiation is taken as a calculation basis for determining the drillhole depth; and it is guaranteed that the drilling depth reaches and goes beyond a mining-induced stress concentration zone when the roadway rock burst is initiated, and the key is to calculate the surrounding rock critical softening depth L.sub.pcr for roadway rock burst initiation;

(72) (4) in the aspect of drillhole spacing, the critical plastic softening zone radius r.sub.pcr for drillhole burst occurrence is taken as the calculation basis for determining the drillhole spacing; and it is guaranteed that the drillhole spacing is enough to destroy the critical plastic softening radius conditions for drillhole burst, and the key is to calculate the critical plastic softening zone radius r.sub.pcr for the drillhole burst;

(73) (5) based on the actual situation of the mine, an inner space for deformation and instability in the drillhole can be formed in the coal mass through the determined drillhole diameter, a deformation absorption space is continuously provided for deformation of the surrounding rock under load, and the burst-preventing effect is strengthened;

(74) (6) the arrangement mode of the burst-preventing drillholes is determined according to the coal seam thickness and the Poisson's effect.

(75) Based on the above-mentioned burst-preventing design principle, as shown in FIG. 5, a specific method for quantitatively determining the drillhole diameter, drillhole depth L.sub.drill and drillhole spacing D.sub.drill of burst-preventing drillholes provided by the present invention is as follows:

(76) I, determining the drillhole diameter according to the arrangement mode of the burst-preventing drillholes in the rock burst roadway and the self-conditions of a mine.

(77) While the construction depth and spacing of the burst-preventing drillholes are quantitatively determined, in the thickness direction of the coal seam, considering that the burst-preventing drillholes in the rock burst roadway are generally arranged in a row or triangle shapes, the mine should adopt large-diameter drillholes for burst preventing of coal mass to the greatest extent on the basis of self-conditions. The influence of the drillhole diameter on the burst-preventing effect in the thickness direction of the coal mass is shown in FIG. 6. In the figure, l.sub.1 and l.sub.2 are the vertical distances from the burst-preventing boundary of the drillhole to the roof and the floor respectively. When the coal seam is thick and the burst-preventing effect in the thickness direction of the coal seam is limited, the triangular arrangement mode should be considered.

(78) Theoretical calculation shows that as the diameter of the drillhole is increased, the damage range of the drilled surrounding rock is increased, and the critical softening zone radius for drillhole burst occurrence is increased. At present, the maximum drilling diameter of a mining roadway drilling machine is about 0.4 m, and the common diameter is 0.05 m to 0.2 m. Therefore, increasing the drillhole diameter is conducive to increasing the single-hole burst-preventing effect, correspondingly increasing the drillhole spacing and improving the efficiency of drillhole construction. The drillhole diameter mainly depends on the power of the mine drilling machine. The influence factor is single and easy to determine. Therefore, the determination of the drillhole diameter is the premise for further determining the drillhole spacing parameter.

(79) II, using the mining-induced stress concentration zone on the roadway side as a limit equilibrium zone of the roadway surrounding rock, and also as a roadway rock burst initiation zone. This zone is a driving stress source for drillhole burst occurrence and roadway rock burst initiation. The acting main object of the drilling burst-preventing is this roadway rock burst initiation zone, as shown in FIG. 7. Therefore, the depth L.sub.drill of burst-preventing drillholes should not only penetrate through the mining-induced stress concentration zone of the current roadway, but also penetrate through the surrounding rock critical softening depth L.sub.pcr for rock burst initiation.
L.sub.drill=η.sub.dη.sub.LL.sub.pcr  (6)

(80) wherein η.sub.d is a correction coefficient for coal seam thickness; when the coal seam thickness is greater than 0 m and less than 4 m, the value range of η.sub.d is 0.8≤η.sub.d≤0.9; when the coal seam thickness is 4-8 m, the value range of η.sub.d is 0.9<η.sub.d≤1.0; when the coal seam thickness is greater than 8 m, the value range of η.sub.d is 1.0<η.sub.d≤1.2; the specific value of η.sub.d in each value range is determined according to the actual construction working condition; is a burst-preventing safety coefficient for the drillhole depth, and the value of the safety coefficient is associated with the burst risk of the zone to be subjected to drilling construction, so that the determination of the drillhole depth is related to the stress of the roadway. Two determination methods are provided for η.sub.L: one is to determine η.sub.L according to the critical mining-induced stress index of burst risk evaluation, namely η.sub.L=0.85+0.5K.sub.cr, and this method has the advantages that burst risk characterization adopts a continuously quantified numerical value interval; and the other method is to determine η.sub.L according to a burst risk level obtained by burst risk evaluation based on a comprehensive index method commonly used at present, generally, η.sub.L is 1.3 in a strong burst risk zone, η.sub.L is 1.2 in a moderate burst risk zone, and is 1.1 in a weak burst risk zone.

(81) III, determining the drillhole spacing D.sub.drill based on the critical plastic softening zone radius r.sub.pcr for drillhole burst occurrence, as shown in the following equation:
D.sub.drill=2η.sub.pcrr.sub.pcr  (21)

(82) combining formula (21) with formula (17) to further obtain

(83) $\begin{array}{cc}{D}_{\mathrm{drill}}={\eta }_{\mathrm{pcr}}d\sqrt{\frac{\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{1}{K}+1}& \left(22\right)\end{array}$

(84) wherein η.sub.pcr is the burst-preventing safety coefficient for the burst-preventing drillhole spacing, and d is a drill bit diameter in burst-preventing drilling construction, d=2r.sub.0.

(85) The value of the burst-preventing safety coefficient η.sub.pcr for the burst-preventing drillhole spacing is associated with the burst risk of the zone to be subjected to burst-preventing drilling construction, so that the determination of the drillhole spacing is related to a current environmental load. Two determination methods are provided for η.sub.pcr: one is to determine η.sub.L according to the critical stress index method of burst risk evaluation, namely η.sub.pcr==2.325−1.75K.sub.cr and this method has the advantages that burst risk characterization adopts a continuously quantified numerical value interval; and the other method is to determine η.sub.L according to a burst risk level obtained by the burst risk evaluation method based on the comprehensive index method commonly used at present, generally, η.sub.L is 0.75 in a strong burst risk zone, η.sub.L is 1.10 in a moderate burst risk zone, and is 1.45 in a weak burst risk zone.

(86) In the calculation determination formula (20) of the drillhole spacing,

(87) $\sqrt{\frac{\alpha }{\beta }}\sqrt{\left(1-\xi \right)\frac{1}{K}+1}$
expresses a coal mass property factor, η.sub.pcr expresses a stress concentration degree factor, namely, the burst risk, and d expresses a geometric dimension factor of the drillhole diameter.

(88) In this embodiment, the burst-preventing drilling means is utilized to actively prevent and control the rock burst for stoping of 394 working face in the 5# district of the mine, and burst-preventing drilling is implemented 200 m ahead of the two stoping roadways in the working face. The drillhole diameter is 150 mm. Specially, for the zone with strong rock burst risk, the burst-preventing drillhole depth is 15 m, the hole spacing is 1.2 m, the drillholes are arranged perpendicular to the axial direction of the roadway, and the drillhole is 0.5-1.5 m away from the floor. When the working face is enabled to enter the strong burst risk zone of 340 m-487 m by pushing mining, and the burst risk of the surrounding rock is detected through a drilling cuttings method, dynamic phenomena of in-hole bursts, ultra-high drilling cuttings amount and suction and sticking occur in a drilling cuttings amount detection hole for many times. Such phenomena indicate that under the current construction parameters, burst-preventing drillholes fail to destroy the critical conditions for drillhole burst occurrence so as to achieve the purpose of preventing roadway rock burst initiation. As shown in FIG. 8, the maximum drilling cuttings amount of a single hole per meter is 70.0 kg/m, far exceeding a rock burst warning value of 4.3 kg/m.

(89) In order to enhance the burst-preventing effect of the drillholes, in this embodiment, based on the burst-preventing drillhole parameter determination method for the rock burst roadways in the coal mines, the burst-preventing drillhole depth is obtained through optimization calculation, namely, 30.67 m for the strong burst risk zone, 28.31 m for the moderate burst risk zone, and 25.95 m for the weak burst risk zone; the burst-preventing drillhole spacing is obtained through optimization calculation, namely, 1.08 m for the strong burst risk zone, 1.58 m for the moderate burst risk zone and 2.08 m for the weak burst risk zone. See table 1 for details.

(90) In this embodiment, according to the optimization design results of the burst-preventing drillhole parameters, when the drillhole spacing is adjusted to be 1.08 m and the drillhole depth is adjusted to be 30.67 m in the strong burst risk zone of 340 m-487 m of the working face, the pulverized coal amount of the drilling cuttings amount detection hole is reduced to 3.2 kg/m, the situations that the drilling cuttings amount is ultrahigh, suction power appears and the like do not occur, and the burst-preventing effect is improved to a large extent.

(91) Table 1 Main parameters of roadway and burst-preventing drillhole, critical value of bursts and determination results of burst-preventing drillhole parameters

(92) TABLE-US-00001 TABLE 1 Main parameters of roadway and burst-preventing drillhole, critical value of bursts and determination results of burst-preventing drillhole parameters Serial Parameter Name of main control number type parameters Symbol Unit Value 1 main burst modulus index K — 1.30 2 parameters uniaxial compressive σ.sub.c MPa 11.12 of roadway strength of coal rock 3 and elastic modulus of coal rock E Gpa 3.58 4 drillhole internal friction angle of φ degrees 30 coal rock 5 residual modulus reduction λ.sub.2 MPa 8.00 6 residual strength coefficient ξ — 0.20 7 support stress p.sub.s MPa 0.34 8 equivalent radius of roadway ρ.sub.0 m 2.37 9 drillhole diameter d m 0.15 10 critical critical mining peak stress for P.sub.mcr MPa 50.98 conditions roadway rock burst initiation 11 for roadway optimized critical mining P*.sub.mcr MPa 45.37 and peak stress for roadway rock drillhole burst initiation bursts critical stress for drillhole burst P.sub.hcr MPa 28.75 12 surrounding rock critical L.sub.pcr m 23.59 softening depth for roadway rock burst initiation 13 critical plastic zone radius for r.sub.pcr m 0.72 drillhole burst occurrence 14 optimization coal seam thickness correction η.sub.d — 1.0 results of coefficient burst- burst-preventing safety η.sub.L — 1.10, 1.20, preventing coefficient for drillhole depth 1.30 drillhole (classified as weak, moderate parameters and strong burst risk levels according to the comprehensive index method) 15 burst-preventing safety η.sub.pcr — 0.75, 1.10, coefficient for drillhole spacing 1.45 (classified as weak, moderate and strong burst risk levels according to the comprehensive index method) 16 drillhole depth strong L.sub.drill m 30.67 moderate 28.31 weak 25.95 17 drillhole spacing strong D.sub.drill m 1.08 moderate 1.58 weak 2.08 Note: “strong” represents that the roadway to be subjected to burst-preventing drilling construction has the strong rock burst risk, “moderate” represents that the roadway to be subjected to burst-preventing drilling construction has the moderate rock burst risk, and “weak” represents that the roadway to be subjected to burst-preventing drilling construction has the weak rock burst risk.

(93) Note: “strong” represents that the roadway to be subjected to burst-preventing drilling construction has the strong rock burst risk, “moderate” represents that the roadway to be subjected to burst-preventing drilling construction has the moderate rock burst risk, and “weak” represents that the roadway to be subjected to burst-preventing drilling construction has the weak rock burst risk.

(94) In addition, the present application further provides apparatus for controlling drilling for rock burst prevention in a coal mine roadway, the apparatus comprising a memory and a processor, the memory storing a program, and the program being executed by the processor to perform the method for controlling drilling for rock burst prevention in a coal mine roadway. The memory may comprise a volatile memory in a computer-readable medium, a random access memory (RAM) and/or a non-volatile memory, etc., such as a read-only memory (ROM) or a flash memory (flash RAM), and the memory comprises at least one memory chip. The memory is an example of a computer-readable medium.

(95) Finally, it should be noted that the above embodiments are only utilized to illustrate the technical solutions of the present invention and not to limit the same. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that the technical solutions described in the foregoing embodiments can be modified or some or all of the technical features thereof can be equivalently replaced. These modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope defined by the claims of the present invention.