Analysis and Application on Controlling Thick Hard Roof Caving with Deep-Hole Position Presplitting Blasting

For the thick hard roof (THR) in Datong mining area, mining operations often led to large-scale hanging-roof and frequent and strong strata behavior, threatening mining safety seriously. Based on the instability mechanism, the fracture model for THR was established, including rock blocks articulation and combined cantilever beam, and the limit initial and periodic intervals of THR were determined to be 36.0m and 8.0m, respectively.&e study proposed the deep-hole presplitting blasting (DPB) for weakening THR for mitigating strong strata behaviors. Blasting-induced fracture characteristics were calculated, determining the charging coefficient and holes spacing. LS-DYNA was employed for establishing a DPB model to analyze crack evolution under the synergistic action of blasting stress wave and detonation gas and the attenuation characteristics for rock peak particle velocity, verifying the rationality of blasting parameters. Field measurement analysis indicated that the immediate roof induced a timely collapse to fill the goaf and the THR was effectively cut off near the presplitting line. Meanwhile, the working resistance was utilized with safety allowance. &e field application showed the DPB on controlled THR caving achieved the significant effect.


Introduction
In the Datong Mine area, coal is mined from Jurassic and Carboniferous coal seams at present, which occur below multilayer 10-25 m complex sandstone roofs.Rock mass has high strength and strong integrity, with the linear density of fracture only being 1.1 strip/m.In the process of mining, the thick hard roof (THR) easily forms the large-sized hanging, inducing energy accumulated and abutment pressure extended far [1].When the limit span of THR is reached, the sudden fracturing would lead to complex migration and strong strata behaviors, including the coal wall spalling, roadside deformation and instability, supports breaking-off, and pillars being crushed [2,3].Furthermore, the selection on supports type with matchless performance and low efficiency for controlling THR would worsen the working face condition, threatening the safety and high efficiency.Consequently, presplitting controlling measures and supports selection should be synergistically employed for controlling THR.
To reduce the strong pressure behavior, positive controlling measures should be employed on presplitting and weakening THR in advance for decreasing caving intervals [4,5].e rock blasting and hydraulic fracturing have been applied for weakened THR and dealing with stress concentration [6][7][8].Meanwhile, the deep-hole presplitting blasting (DPB), with convenient operation and huge economic benefits, was widely used in mine-induced stress transferred and rock mass fractured, achieving a remarkable technical effect [9,10].DPB application could actively weaken the integrity of THR between blasting holes at the specific position, form the presetting angle and the interpenetrated fracturing zone, and effectively reduce the length of roof suspension.Near the presplitting line, THR breaks into rock blocks with optimized caving intervals.Meanwhile, blasting the fracturing zone effectively cuts off the transferring way of stress, and achieves stress relief to some extent, which mitigates the strata behavior.Domestic and foreign scholars have carried out many related studies on the structure characteristics of the roof initial caving, the presplitting mechanism and techniques for controlling THR migration.Liu et al. analyzed the support resistance variation under the different combination of layer key-stratums caving and combined selecting reasonable support type by using hydraulic fracturing technology to control THR, achieving good effect on the pressure behavior [11].Yu et al. put forward explosive or liquid carbon dioxide blasting to weaken roofs and large pillar for reducing stress concentration degree through surface drilling [8,12].Wang et al. focused on the initial weighting interval of thick roof in shallow depth seams, employed LS-DYNA3D simulation of DPB to reveal the blasting-controlled roof caving mechanism, and optimized blasting parameters, achieving the expected effect on-site [13].Yang et al. proposed the DPB technology confined blasting in high pressure water-filled medium instead of the conventional air medium, the result showing a significant technical effect on roof presplitting and pressure relief, but the application of the technique was restricted because of the technical complexity and various limitations on the procedure [14].Ning et al. analyzed the mining-induced movement and breaking characteristics of the double-layer THR via microseismic monitoring in longwall top-coal caving working face.After adopting the systematic long-hole blasting, the longwall panel could be smoothly extracted with low-frequency instability [15].
From the viewpoint of the safety controlling roof, the DPB could effectively reduce the size of breaking roof, especially the suspended length of the THR, alleviating the strong strata behaviors.However, the above research generally aimed at the presplitting of THR initial caving, and lacked the study on the synergistic controlling roof between the working resistance of the support and the reasonable periodic caving intervals as well as the optimization of the blasting technical parameters.In the study, according to the THR occurrence characteristics of the 8939 working face in Xinzhouyao Mine of Datong Mining area, theoretical analysis was employed for analyzing the action relation between THR caving intervals and working resistance of supports, obtaining the reasonable blasted horizon, the initial presplitting caving interval, and periodic caving interval with presplitting angle.e characteristics and influencing factors of blasting fracturing were analyzed, and meanwhile, the fracture evolution and the rock peak particle velocity (PPV) attenuation were simulated with LS-DYNA, which discussed and determined the rationality of the optimized blasting technical parameters.Subsequently, the presplitted roof collapse characteristics and support resistance were observed and analyzed through field measurements.Consequently, the rationality and efficiency of the DPB were verified on reducing the caving intervals and controlling the strong strata behavior.

Fracture Characteristics and Caving Intervals
Controlling Mechanism for THR which is characterized by the stable occurrence.e longwall-mechanized top caving is applied for mining along strike, with an advancing length of 1176.0 m and the face length of 104.5 m.Mining height is 3.0 m, with a caving ratio being 1 : 1.4.e layout of working face is shown in Figure 1.Immediate roof is 1.8 m mudstone and 2.7 m siltstone, the main roof above which is composed of 13.7 m medium sandstone and 14.2 m gritstone.e uniaxial compressive and tensile strength of the medium sandstone are 63.5 MPa and 6.05 MPa, respectively, which belongs to typical THR.Lithological characteristics are shown in Figure 2.With no joints and high stability of the rock, the main roof does not break readily as the face advances, inducing the largesized hanging roof and increasing the possibility of strong strata behavior.erefore, the presplitting on THR is expected to form vertical weaken planes, reduce the sizes of caved blocks, and thereby improve mining safety degree [16].

Characteristics of the THR Initial
Breaking.With advancing the working face, when the limit span is reached, the THR is broken up into "V"-shaped structure.For the THR thickness is large, the fractured roof produces horizontal compression and contacting friction after rotary [17,18], which is balanced with the overburden loading under the synergistic action of supports, as shown in Figure 3(a).Where H, l, and θ are the thickness, length, and rotating angle of THR blocks, respectively, a is the height of the articulated surface, Δh is the roof movement subsidence, q is the uniform distributed overburden loading, φ is the internal friction angle of rock mass, and B is the support width.e mechanical model is shown in Figure 3(b).
According to Figure 3(b), the geometric relation could be obtained as follows: For rock blocks remaining stable before contacting gangue, the moment balance is Taking tan φ as the friction coefficient of the THR articulated surface [19], the friction force R 0 between rock blocks is obtained as follows: Considering the impacting action of the THR initial breaking, the top coal is reserved instead of caved for filling 2 Advances in Civil Engineering goaf to act as a cushion layer of the THR at the initial stage of production in the working face.erefore, 13.7 m medium sandstone is determined as the initial broken-articulated strata, and the relation between the support working resistance and the initial intervals is shown in the following equation: With the existing support resistance of 6000 KN, the critical length of the articulated block is determined with 18.8 m.On the basis of safe production on-site, the limit initial caving interval of the roof is set to be L e 36.0 m. erefore, the location of presplitting is determined near the open-o cut, 18.0 m and 36.0 m from the open-o cut, respectively.

Synergetic Controlling on the THR Periodic
Caving and Supports

Characteristics of the THR Periodic Caving Interval.
After the initial collapse of the THR, with advancing the working face, the hanging dimension is large enough, leading to a cantilever structure with the free end in the goaf (Figure 4(a)).erefore, the mechanical model of the cantilever beam with uniformly distributed loading is proposed and employed for analyzing characteristics of the THR periodic breaking [11], as shown in Figure 4(b).
e following stress components of the cantilever beam are obtained [20]: where σ x , σ y , and τ xy are the horizontal, vertical, and shear stress components, respectively, and l z is the half length of the periodic caving interval.
Combined with the stress distribution characteristics of the cantilever beam, the horizontal tensile stress at the xed end (l z , −h/2) reaches the maximum.Consequently, the fracture of rock mass conforms to erefore, the limit periodic caving interval l z could be obtained:

Determination on the ickness of Cantilever Stratication and Loading.
For the large area roof suspended before THR periodic breaking, the caving interval of the upper thin and soft roof is consistent with the lower THR.Considering the fracturing angle of the rock, the hanging roof presents inverted trapezoid.erefore, the "immediate roof" combined cantilever beam structure is put forward [11,12], and the caving height full of the goaf is taken as the structure thickness.e overburden loading is mainly carried by the articulated layer, under which the THR with the cantilever structure acts on the support.e broken and instability structure is shown in Figure 5. e total thickness of the cantilever beam is where H z is the cantilever beam thickness, k i is the coe cient of bulk increase, 1.33, h m is the mining thickness, and h i is the strata layer thickness.e presplitting layer is determined with numbers 1-4, that is, 1.8 m sandy mudstone, 2.7 m nestone, 13.7 m medium sandstone, and 4.2 m nestone are the layers of the  Advances in Civil Engineering cantilever beam.
e vertical height of presplitting holes, consistent with the total thickness of the structure, is not less than 22.4 m, indicating the 13.7 m medium sandstone is completely in the scope of presplitting.
Supposing the volume force is c i and elastic modulus is E i of layered strata, where i 1, 2, m, . . .n. Combined with the key strata theory, the rock-bearing capacity (q n ) m could be obtained [21]: where (q n ) m is the loading exerted by n layer on m layer rock mass.When loading accords with Equation ( 9), q m (q n ) m could be determined.Subsequently, combined with Equation (7), periodic caving intervals of the layered cantilever beam are shown in Table 1.

Determination on the Presplitting Periodic Intervals of THR.
Combined with the parameters of periodic intervals and the cantilever beam shown in Table 1, considering 14.2 m gritstone forming the articulated structure, the support working resistance is determined with the loading of 22.4 m combined cantilever beam and the additional loading applied by articulated strata.According to the mechanical characteristics of the combined cantilever beam, the critical stability condition is obtained: where P t is the top coal unit weight, c is the horizontal distance between the equivalent action point and coal wall, l i is the layer periodic caving dimension, P i c i h i l i , R x is the additional load of articulated blocks acting on the cantilever beam, l x is the acting length of the additional load point on the layered strata, α i is the fracturing angle of roof, simpli ed as 80 °, and θ is the presplitting angle.As shown in Equation (10), the main controlling factors of presplitting are the periodic presplitting intervals and angle.Based on the cantilever beam parameters, the critical working resistance could be obtained: Because of the support working resistance exceeding 25.5 MN, in order to alleviate the strata behavior and ensure  Advances in Civil Engineering support's safe operation, it is necessary to presplit the THR and control migration.After 13.7 m medium sandstone presplitted, the limit length of the cantilever beam is the reasonable THR caving interval.For the 14.2 m gritstone forming the articulated structure, the additional loading decreases with the cantilever length shortening [8,22].e mechanical model after presplitting is shown in Figure 6.After reducing the periodic intervals, the support critical stability condition is shown in Equation ( 12). e variation curve between the presplitting angles and periodic caving intervals is shown in Figure 7.
As shown in Figure 7, the curve between the presplitting angle and the periodic caving interval is approximately linear.
e limit periodic interval is 9.0 m, with the corresponding presplitting horizontal angle is 0 °.
However, the smaller presplitting angle may lead to roof broken and cutting-o near the presplitting position, increasing the management di culty.However, the drilling workload and charging di culty increase with the rise of the horizontal rotation angle, and the location is relatively complicated.Meanwhile, the horizontal rotation angle varying from 5 °to 10 °has little in uence on the presplitting e ecting [1,4], and the corresponding periodic presplitting intervals are 7.5 m-8.4 m. erefore, considering the presplitting angle, the periodic interval, and drilling workload, the optimal presplitting interval and blasting-hole horizontal rotation angle are set to be 8.0 m and 7 °, respectively.

Analysis on Fracturing Mechanism by Blasting
In order to relieve the strata behaviors level, the key is to reduce the caving intervals for cutting o the loading transferring way and realizing the stress relief.DPB, because of the simple construction technology and strong adaptability, is widely used in presplitting the THR as an e ective technology [10].

Fracturing Mechanism of the DPB.
e DPB stimulates overpressure near blasting holes exceeding the dynamic compressive strength of the rock mass and forms a compression crushed zone (Crushed zone I).During the period, the pressure rapidly attenuates to the compression stress wave on the boundary of the crushed zone.After that, the radial fracture is induced by the reversely releasing of compressive stress and the main and wing cracks are interconnected, forming the initial concentric fracture network (Fractured zone, stage II) [23,24].
After the initial crack forming, the detonation product di uses uniformly in the fracture zone, which apply a quasistatic loading to the crack tip and produce the secondary propagation (Fractured zone, stage III).It penetrates through the fracturing zone induced by the detonation gas of the adjacent borehole, which achieves e ective presplitting.
e characteristic partition is shown in Figure 8.

Mechanism of Fracturing Induced by Stress Waves.
e peak intensity of a transmitted shock wave generated by explosion of a cylindrical charge [25], P m , is given by the following equation: where n 0 is the stress intensi cation factor, ρ 0 is the explosive density, D 0 is the velocity of explosive, k r r c /r b , r c , r b are radius of blast hole and charge, respectively, c is the thermal insulation factor, 3, and l c is the axial decoupling charge coe cient.Combined with three-direction stress intensity coe cient C, the radius of crushing zone R c and the radius of the initial crack zone R P are obtained, respectively [26]: where σ c is the static uniaxial compressive strength, ξ is the loading strain rate of the rock, σ td is the dynamic uniaxial tensile strength, b is the side pressure coe cient, α and β are shock and stress wave attenuation coe cients, 2 + b and 2 − b, respectively, and μ d is the dynamical Poisson ratio.
e fracturing characteristic curve of k r (1.0, 1.11, 1.25, 1.5, and 1.875) on main cracks number and (R c + R p ) is shown in Figure 9. e variety on the fracturing length and number of main cracks took k r 1.25 as an in ection point on the whole.With k r increasing (k r ≤ 1.25), the length of main cracks slightly rose, and meanwhile, the number rstly increased and then kept stable.When k r > 1. 25, the e ect of the fracturing sharply dropped with both decreasing linearly.erefore, k r was nally identi ed as 1.25.
Combined with Equation ( 14) and Table 2, R c and R p were set to be 0.13 m and 1.95 m, respectively, obtaining R c + R p 2.08 m.

Mechanism of Crack Propagation Driven by Detonation
Gas.On the basis of the initial cracks, the detonation gas di uses uniformly into the crack tip, resulting in the tensile yield and simulating the cracks secondary propagation.When the quasi-static pressure of detonation gas drops to the critical fracturing value of the brittle rock, the crack reaches the maximum length [14].e ultimate pressure (σ l ) is K IC is the static fracture toughness of rock mass and l k(max) is the maximum length of secondary fracturing induced by detonation gas.
K IC could be measured experimentally as follows [27]: where P max is the peak loading of specimens, a is the crack length, B is the thickness, R is the radius of the disc, and F is the dimensionless stress intensity factor with a low loading rate, F 1.0.Test results are shown in Table 3. e detonation gas fracturing is mainly along the blasthole radial, and the maximum opening of the crack u is in the boundary line of the initial crack zone.Based on the tip, the angle of cracks is 0 °approximately [28], the detonation gas in the isentropic di usion, accords with

Advances in Civil Engineering
where n is the number of main cracks and G is the shear modulus.
Combined with Equation ( 17), the relationship between l k(max) and P x could be obtained.
Based on Figure 9, the development number of main crack is set to be n � 6-8, and meanwhile, the number of main cracks in the mathematical model is 4-8 [29]; therefore, the number of main cracks for detonating gas diffusion is determined to be 4-8.Figure 10 shows the relationship curve between l k(max) and P x .
As shown in Figure 10, l k(max) is monotonically increased with P x and the number of main cracks has a significant effect on the fracturing length induced by detonation gas; that is, under the same fracturing length, the larger the number of main cracks is, the higher the critical gas pressure needed for crack propagation is.Combined with Equation ( 17), the steady pressure of detonation gas in initial cracks is calculated at 58.5 MPa.
e maximum number of main cracks at present (n � 8) corresponds to l k(max) � 0.52 m. erefore, the length of secondary fracturing is l k(max) � 0.52 m.
en, the blast-hole spacing could be determined with L erefore, the optimal blasting holes spacing is set at 5.0 m.

Numerical Model and Constitutive Equation.
e LS-DYNA is employed for establishing the 3-D model, and ALE is applied for analyzing the diffusion of detonation gas and the fracture evolution.Considering the computing time and simulation precision, the model size is set to be 20.0 m × 12.0 m × 15.0 m with nonreflecting boundary condition.Horizontal spacing of blast holes is made at intervals of 5.0 m, with each having a cartridge diameter of 40 mm (k r � 1.25) and depth of 15.0 m through the whole model.Test points arrangement is shown in Figure 11.
According to the field construction and related contemporary research results, the numerical model and field experiment adopt the Class-2 coal mine permissible emulsion explosive [30].e MAT_PLASTIC_KINEMATIC and EOS_LINEAR_ POLYNOMAL keyword offered by LS-DYNA are employed for characterizing the blasting influential process on the rock.e parameters of MAT_HIGH_EXPLOSIVE_BURN and JWL Equation (19) are used for describing the relation between diffusion volume and pressure of crack tips, with explosive parameters shown in Table 4 [31,32].
A, B, R 1 , R 2 , and ω are the performance parameters of explosion, and E 0 and V are the internal energy and volume of detonation gas, respectively.
Table 5 shows the rock mechanical parameters, and Table 6 presents the air status parameters.

Numerical Results and Analysis.
In order to analyze the law of blast fracturing between blasting holes, numerical simulation with the DPB model (L � 5.0 m) was employed for analyzing cracks evolution and rock particle vibration velocity.
(1) Dynamics of Cracks at a Cross Section. Figure 12 shows the dynamic evolution characteristics of the main and wing cracks at a cross section during blasting.
Figure 12 displays the fracture evolution in a profile vertical to the blast-hole axes.Figure 12(a) shows the morphological development of the crushed zone (region A) formed at 149.8 μs, with R c � 0.12 m. e diffusion  When the fracturing length of detonating gas was 0.4 m, the fracture penetration further promoted the development of the wing cracks (R k � 1.94 m) and the rock failed primarily with the direction vertical to the line passing through the holes center (red trajectory in region C).A mount of wing cracks occurred around the symmetry axis between two blast holes and strengthened the development degree of fracturing network and fragmentation lumpiness of the rock.Simulation results confirmed that the structural parameters of charge determined could ensure the favorable presplitting effect on the rock.
(2) Attenuation Law of Blasting Vibration Velocity on the Rock between Blast Holes.After blasting, the different vibration speed corresponds to the different rock damage degree.e higher the PPV of rock mass is, the greater the damage degree of the corresponding rock mass is [33,34].
e PPV critical value of the rock mass is shown in Table 7 [35,36].Based on the integrity of 13.7 m medium sandstone, the demarcation point corresponding to the rock PPV is obtained.
Figure 13 shows the attenuation law of PPV at different positions on the blasting holes profile.In the direction of X, the particle vibration velocity presented the attenuation trend on the whole.e PPV of point A near the blasting hole was up to 38.96 m/s, resulting in the rock mass being approximately broken, and meanwhile, the PPV of point E was 3.35 m/s exceeding the critical value of fracturing.e variation of PPV showed that DPB had a significant effect on the rock presplitted along the line through the holes center and the fracture intensive development in the middle of blast holes.In the direction of Y, the PPV was obviously lower than that of X direction, indicating that the initial particle vibration was dominated by radial compression.After 400.0 μs, the particles C, D, and E showed a slight increase in the vibration velocity, which revealed that the detonation gas    lagging behind blasting stress wave simulated the particle vibration once more, promoted the secondary fracturing of cracks, and formed the interconnecting fractured network between the blasting holes.e PPV evolution indicated that the DPB in 13.7 m medium sandstone achieved the crushing and fracturing of the rock between the blasting holes, which showed the rationality and high e ciency of the blasting parameters.

Technological Parameters of Presplitting Blasting.
Combining theoretical analysis and simulation results, it was determined that the initial and periodic caving intervals of the 8939 working face (key: 13.7 m medium sandstone) were 36.0 m and 8.0 m, respectively.e technological parameters of DPB were determined and implemented.e arrangement of the presplitting blasting holes is shown in Figure 14.
In the process of mining, the initial blasting fracturing of the THR was divided into 3 groups, in which the rst was set at open-o cut and then the presplitting interval was 18.0 m (Twice).After implementing the initial caving presplitting of 13.7 m medium sandstone, the periodic interval of blasting   presplitting was set to be 8.0 m, and the distance between the drilling position and the oor was 1.2 m.In order to avoid the process interferences, serial initiation of two boreholes was generally carried out during the maintenance crew, and the distance between initiation location of the periodic DPB and coal wall was 30 m at least.A reverse charge was used in the blasting hole, in which the detonator and the rst volume explosive were bundled at the bottom of the hole with the cartridge speci cation being 40 mm × 500 mm (Figure 15(a)).e cartridge must be tight and the connection was reliable.e sealing section length was not less than 1/3 of the holes depth, in which the rigid (concrete) was used in conjunction with plasticity (stemming) structure (Figures 15(b) and 15(c)).e stemming in middle cushioned the blasting energy, and the lower concrete resisted the residual energy (Figure 15(d)).
e explosive and the sealing material are shown in Figure 15, and the charge structure is shown in Figure 16.e technical parameters of the roadside and roof presplitting holes are shown in Table 8.After blasting, drilling holes were tested to ensure the blasting e ect.

Advances in Civil Engineering
the presplitting e ecting of the THR was analyzed based on the eld observation, including the working condition of supports and caving morphology.Shapes of the presplitted roof collapse and characteristics of roadside instability in the technical roadway are shown in Figure 17.
e caving morphology of THR in the technical roadway was monitored while mining the working faces.When the working face was advanced near the blasting presplitting line, the presplitted roof was collapsed in time without suspended roof and lagging collapse, and the maximum size of the caved roof fragmentation was 2.5 m × 2.0 m by visual measurement, as shown in Figure 17(a).Figure 17(b) reveals the destroyed degree on the roadside after DPB. e average width of the fractured zone was about 5.7 m with the maximum 8.0 m, and the height was the whole roadside.Meanwhile, the maximum broken depth of the roadside near the blasting hole reached 2.3 m, and the bolts were completely ushed out of the roadsides.e above showed the good technical e ect of DPB.

Working
Resistance of the Supports.ZFS6000/22/35 type caving coal hydraulic support was applied for the 8939 working face mining.e characteristics of the support working resistance were monitored, as elaborated in Figure 18.
Figure 18(a) shows that, during mining, THR breaking led to "Minor-Major periodic weighting," in which a major weighting was usually accompanied by 2-4 minors.13.7 m medium sandstone rstly collapsed with the working face advanced at 36.6 m and the work resistance was up to 5490 KN, which met the supporting strength.en, with working face advanced, the sandstone formed "minor-periodic weighting" and periodic intervals were stable at 7.2 m-10.4 m, coinciding   Advances in Civil Engineering 13 with the theoretical analysis.When the working face was advanced at nearly 100.0 m, 14.2 m fine sandstone firstly collapsed, resulting in local strong strata behavior, and the maximum working resistance reached 6160 KN. en, the periodic breaking interval of 13.7 m medium sandstone was 7.4 m-8.6 m and that of 14.2 m fine sandstone was 30.2 m-38.8 m. e maximum working resistance reached 5860 KN, matching with support type.
Figure 18(b) shows that the working resistance of the support largely located in the interval of 4000 KN-6000 KN and made up 79.69% of the totality, making the support capacity be used effectively.
e time weighted working resistance was 4765 KN, accounting for 79.4% of the rated value (6000 KN).Accordingly, the support had sufficient safety margins in the production process, ensuring the safe production.

Conclusions
On the basis of the controlled THR fracturing in the present study, the following conclusions can be drawn: (1) e initial caving interval of the THR above the coal seam has been large.A large suspended roof gradually has been formed during the initial mining stage, and its rotation and sinking with fracture dynamic load have produced high side abutment pressure on the working face, which has been the main source of strong strata behaviors.(2) e fracture mechanical model on the rock block articulated of the THR initial caving and the periodic caving of the cantilever beam is established.e relationship between the characteristics of the THR firstly fracturing and the support effecting is analyzed, obtaining the limit initial caving interval.Combined with the combined cantilever beam instability and characteristics of articulated rock stratum, the periodic caving interval and presplitting angle of the presplitted roof are obtained, with the optimal L aze � 8.0 m and θ � 0-7 °.
(3) DPB, as an effective presplitting technology, is employed for the THR fracturing.e length and characteristic zoning of DPB stress wave and detonation gas are calculated, respectively, and the influence factors of holes spacing are analyzed.Based on the results, the reasonable value of blasting technical parameters and holes spacing are optimized and determined.(4) LS-DYNA3D was used for revealing the mechanism of crack propagation through evolution and interconnection of main and wing cracks.Meanwhile, the attenuation curve of PPV was obtained for analyzing rock damage scope.e result demonstrated that the rock between holes was successfully prefractured.(5) A plan of blasting technological parameters on DPB was designed and applied in 8939 working face of Xinzhouyao Mine.Field measurement indicated that the blasting presplitted roof has achieved a timely and complete collapse.e DPB reduced the initial and periodic caving intervals effectively and the working resistance of the support was obviously reduced, which showed good control on the strata behavior.e remarkable technical effect has been obtained.
It should be noted that the limit caving intervals were obtained, only considering the static loading action of the articulated-strata and the broken THR.Because the additional dynamic loading induced by the fracturing of the lower THR and overlying strata, imposed on the supports, when the working face is advanced at about 100 m (14.2 m gritstone broken), the working resistance of the supports increases sharply to 6160 KN, exceeding the rated working resistance, which is larger than the theoretical calculation value.
erefore, the future research emphasis is to supplement the calculation of the additional dynamic load applying on the supports, and simulate the characteristics of the dynamic fracturing and the relation between support and surrounding rocks under different presplitting caving intervals and angles for analyzing the THR controlling effect.Meanwhile, auxiliary holes should be added to observe the fracture evolution for analyzing presplitting effect before and after blasting in site.

Figure 3 :
Figure 3: Mechanical analysis on the THR initial breaking.

Figure 4 :
Figure 4: Load-bearing characteristics of the roof structure with cantilever supported condition.

Figure 5 :
Figure 5: Structure characteristics of the combined cantilever beam-hinge structure of overlying strata.

Figure 6 :
Figure 6: Mechanical model of the combined cantilever beam after presplitting.

Figure 7 :
Figure 7: Relationship between the presplitting angle and periodic caving intervals.

Figure 9 :
Figure 9: Effects of k r on the development of main cracks.

Figure 11 :
Figure 11: Model size and layout of observation points.

Figure 10 :
Figure 10: Effect of detonation gas acting on the secondary fracturing length.

Figure 13 :Figure 14 :
Figure 13: Attenuation curve of blasting velocity on the rock elements under di erent directions.(a)Horizontal direction (m/s).(b)Vertical direction (m/s).

Figure 17 :Figure 18 :
Figure 17: Shapes of the presplitted rock collapse and instability.(a) Roof collapse and (b) impact crushing of roadsides induced by blasting.

Table 1 :
Characteristics of strata loading and caving intervals.

Table 2 :
Mechanical parameters of rock mass.

Table 3 :
Fracture parameters of specimens.

Table 8 :
Technical parameters of presplitting blast holes.