After the seven models were solved and computed, damage contour maps of the rock mass were obtained, as shown in
Figure 9. Since fracturing is the main form of action that explosives have on the nearby rock mass, and the rock mass damage variable can represent the fracturing state of the rock mass, the damage extent of the rock mass in different models was analyzed to select the optimal drilling scheme. Based on the analysis results and comprehensive research objectives, this study evaluates the blasting effects using two indicators: the level of residual holes at the bottom and the uniformity of the damage volume.
4.3.1. Comparison of Residual Hole Depth for Different Control Ratio R Values
Although post-processing software LS-PREPOST v4.3 can visualize the damage changes and damage propagation during the blasting process, it cannot directly measure the length of residual holes after blasting. To address this, the SECTION function in the MODEL card of the software can be used to slice the blasted rock mass and obtain cross-sectional views at different positions, as shown in
Figure 10a. The cross-sectional views are parallel to the drilling workface, with an interval of 0.5 m between adjacent sections. The depth
H values of the cross-sectional views from the hole top are successively 2.0 m, 2.5 m, 3.0 m, 3.5 m, and 4.0 m.
The cross-sectional view of the rock mass after the first row of easer holes is shown in
Figure 10b, resulting in five damage contour maps at different depths. Different colors in the damage contour maps represent varying degrees of damage, and by comparison with the legend, the damage values of the mesh elements can be visually obtained. The damage value
D indicates the degree of rock fragmentation: when
D = 0, the rock is undamaged; when 0 <
D < 1, the rock is damaged but not completely destroyed; and when
D = 1, the rock is completely destroyed.
In practical blasting, the damage value when the rock is fractured is less than 1. In this study,
D = 0.75 is used as the damage threshold for rock fracture. When
D < 0.75, the rock is considered not fractured; when
D ≥ 0.75, the rock is considered fully fractured [
33].
Using the damage contour map of the rock mass after the detonation of the row 1 easer holes, the black ellipses in
Figure 11a mark the cross-sections at depths of 3.0 m and 4.0 m, intersecting with the adjacent hole axis plane, as shown in
Figure 11b. The damage values of each grid on these intersection lines are statistically analyzed. In
Figure 11c, it can be observed that at a depth of 3.0 m, the damage values of all grids in the cross-section are greater than 0.75, indicating that this section has experienced complete fracturing. At a depth of 4.0 m, less than half of the grids in the cross-section have damage values greater than 0.75, meaning this section has not experienced complete fracturing.
For easier statistical analysis, the maximum penetration depth of a single hole is defined to represent the rock-breaking efficiency of easer holes. The maximum penetration depth of a single hole is the maximum value where the cross-section in the rock mass shows complete fracturing.
The damage contour maps for row 1 of easer holes in different models are analyzed and shown in
Figure 12. This figure presents the damage contour arrays at cross-sections with depths of 2.5 m, 3.0 m, 3.5 m, and 4.0 m. The horizontal axis represents the control ratio
R corresponding to the minimum burden of the easer hole, and the vertical axis represents the cross-section depth
H. An analysis of the damage contour maps for row 1 of easer holes resulted in the following findings:
For R = 1.0 to 1.2, the maximum penetration depth is between 3.0 m and 3.5 m.
For R = 1.3 to 1.6, the maximum penetration depth is between 3.5 m and 4.0 m.
To obtain a more precise maximum penetration depth, the rock mass was further sliced at 10 cm intervals, as shown in
Figure 13. The maximum penetration depths with different
R values are shown in
Table 7.
For quantification and considering potential drilling errors on-site, the left endpoint of the interval (the smaller value) is taken as the maximum penetration depth value for each hole. The difference between the hole length and the maximum penetration depth is approximately considered as the residual hole length after blasting. So, the maximum penetration depths and residual hole lengths can be calculated as shown in
Figure 14.
The final depths of the residual holes for each row of easer holes in the different models are summarized in
Table 8.
From
Table 8, it can be seen that for the average depth of residual holes, when
R ranges from 1.1 to 1.6, it is consistently less than when
R = 1.0. This indicates that the design scheme where the minimum burden at the bottom of the hole increases proportionally has a better effect in reducing the depth of residual holes. Specifically, the average depths at
R = 1.0 and
R = 1.1 are larger, at 46.7 cm and 45.0 cm, respectively, whereas at
R = 1.3 to 1.5, the depths are smaller, measuring 36.7 cm, 36.7 cm, and 35.0 cm, respectively. Connecting the endpoints of residual holes of all schemes in sequence yields the residual hole section, which characterizes the smoothness of the section after blasting. Similarly, the variance of residual hole depths calculated in this paper can be used to characterize the smoothness of the section. From
Table 8, it can be observed that when
R = 1.0 and
R = 1.6, the variance of residual hole depths is larger, at 0.036 and 0.037, respectively, indicating poorer smoothness of the section after blasting, whereas when
R = 1.3 to 1.5, the variance of residual hole depths is smaller, at 0.006 each, indicating better smoothness of the section after blasting. Therefore, in terms of the smoothness of the residual hole section after blasting, the blasting effect is superior when
R = 1.3 to 1.5 compared to schemes with an equally spaced minimum burden.
4.3.2. Relationship between Rock Mass Damage Volume V and Control Ratio R
The damage volume of the rock mass refers to the volume of rock that meets the fragmentation criteria under the action of detonation waves and detonation products after an explosive detonation. This volume indicates the extent and scope of the damage to the rock mass. As shown in
Figure 15, the damage threshold value of the rock mass is set to
D = 0.75 in LS-PREPOST, displaying only the parts of the rock that meet the fragmentation criteria. Using the Measure function, the volume of the rock mass with
D greater than 0.75 in each hole’s burden area is calculated.
(a) Relationship between V and R for each row of easer holes
The statistical results of the damage volumes of the rock mass burdened by each row of easer holes under different conditions are shown in
Figure 15. The horizontal axis represents the control ratio
R for different schemes, while the vertical axis represents the damage volume
V of the rock mass burdened by each hole, with S1 to S6 representing the first to sixth rows of easer holes in each scheme.
From
Figure 16, it can be observed that as
R increases, the damage volume
V of the rock mass burdened by rows 1, 2, and 3 of easer holes (S1 to S3) gradually decreases. When
R = 1.6, the damage volume is the smallest, reduced by about 17% compared to when
R = 1.0. This is because, with the increase in
R, the minimum burden of the hole bottom for rows 1, 2, and 3 of easer holes gradually decreases, resulting in a decreasing proportion of energy from the explosive being used for rock fragmentation.
For rows 4, 5, and 6 of easer holes (S4 to S6), the damage volume V of the burdened rock mass increases with R. However, for S5 and S6, the damage volumes start to decrease when R increases to 1.4 and 1.5, respectively. After the decline trend appears, the damage volume V of S6 decreases more rapidly. This is because, with the increase in R, the minimum burden of the hole bottom for the last three rows of easer holes also increases, resulting in an increasing proportion of energy from the explosive being used for rock fragmentation. However, when R increases to 1.5 or 1.6, the minimum burden of the hole bottom for the last two rows of easer holes becomes too large, exceeding the destructive capability of the holes. At this point, the energy from the explosive is insufficient to achieve continuous fragmentation from the hole to the new free face.
(b) Relationship between total damage volume Vs and R
Figure 16 shows the relationship between the total damage volume vs of the rock mass and the control ratio
R for different models. The horizontal axis represents the control ratio
R for each scheme, while the vertical axis represents the total damage volume vs for each condition, with specific values of total rock damage volume indicated in parentheses for each model. It can be seen from
Figure 17 that as
R increases, the total damage volume vs gradually increases, reaching its maximum value of 4.638 m
3 at
R = 1.3. Beyond this point, the total damage volume vs begins to decrease, reaching its minimum value of 4.276 m
3 at R = 1.6. Therefore, from the perspective of total damage volume, the optimal blasting effect is achieved when
R = 1.3. However, since the total damage volumes for
R = 1.2 and
R = 1.4 are similar to that for
R = 1.3, with differences within 2%, these two conditions can also be considered as alternative solutions.
In summary, a comprehensive analysis of the blasting effects was conducted by comparing different models from three perspectives: post-blast residual hole depth, section flatness after blasting, and damage volume. With the exception of the R = 1.6 condition, the design scheme in which the minimum burden of the easer hole bottom increases proportionally from row to row achieved better blasting effects than the equal burden scheme. This validates the feasibility of the design principles proposed in this paper. Based on the three blasting effect evaluation indicators, the blasting effect is optimal when R = 1.3. Therefore, in the engineering background of this project, the optimal control ratio for the minimum burden of the easer hole bottom should be 1.3.