Abstract

The basic process of cut blasting is to break rock, throw fragments, and form a cavity. Based on the characteristics of cut blasting and the combined effect of stress waves and detonation gas, the evolution process of wedge cut blasting is divided into two stages, and a theoretical model is proposed to investigate the cavity formation mechanism by theoretical analysis and field tests. In phase one, rock breaking is caused by stress waves. By considering the dynamic strength of the rock, a computational model is built for the rock failure zone derived from the coupled cylindrical charge explosion. In phase two, the driving force of the detonation gas overcomes the total resistance of the surrounding rock mass, accelerates fragments, and then throws fragments to form a cavity. The criterion of cavity formation is established on the basis of the quasi-static loading of the detonation gas. The theoretical model provides an overall interpretation of the cavity formation mechanism, in which stress waves break rock and detonation gas throws fragments. A specific case indicates that the range of the failure zone is approximately 18 times the borehole radius in granite and that the hole-bottom spacing of the wedge cut can be designed as 50 cm; in addition, detonation gas is sufficient to overcome the total resistance, accelerate rock fragments, throw fragments, and form a cavity. Field tests present favourable blasting results, with a high utilization rate of boreholes and uniform fragment sizes. Therefore, the model could provide theoretical support and technical guidance for wedge cut blasting in hard rock.

1. Introduction

With the rapid progress of modernization, the construction of high-speed railways, freeways, and underground rock projects has grown rapidly. Currently, drilling and blasting methods are still the main methods of tunnel and underground rock space excavation. More than 95% of mountain tunnels are formed using drilling and blasting in China [1]. Cut blasting is key to creating another free surface for subsequent blasting and influences the overall blasting procedure [2]. Wedge cuts have the advantages of few boreholes, easy rock casting, and low drilling accuracy requirements; thus, these cuts have been widely used in hard rock excavation. Langefors and Kihlstrom [3] studied cut blasting earlier and proposed some cut patterns. Shapiro [4] compared wedge cuts with other cuts and believed that wedge cuts ensured maximum efficiency in soft rock with relatively shallow blast holes, i.e., shallower than 2.5 m. Chakraborty et al. [5] proved that wedge cuts were more productive than parallel cuts in small drivages and developed an empirical relation for the precise prediction of the powder factor and for the specific drilling with wedge cuts by considering rock strength, Barton’s rock mass quality, and hole depth. Cardu and Seccatore [6] offered statistics of industrial trends in which wedge cuts appeared to be adopted in a wider variety of applications.

Due to the complexity of the mechanism of wedge cut blasting, there are few theoretical studies on the cavity formation of wedge cuts and most of the published literature is about engineering applications. Dai and Du [7] fully considered the failure modes of different sides of cut rock mass and provided a qualitative explanation for the cavity formation mechanism of wedge cuts, believing that the bottom, left, and right sides had been broken under the action of explosion and that shear failure occurred on the upper and lower sides. Based on the quasi-static theory of detonation gas and the limit equilibrium theory of rock mechanics, Shan et al. [8] analysed the generation process of a cavity by considering that upper, lower, left, and right sides were the result of shear damage and that the underside was due to tensile damage. Furthermore, numerical simulation is used extensively in cut blasting due to its rapid implementation and low cost [9]. Using the dynamic analysis software LS-DYNA, Wang and Konietzdy [10] simulated stress and crack evolution in an infinite rock mass with a single free face during blasting. Yang et al. [11] performed a study of duplex wedge cut blasting with 12 holes via numerical simulation, and the results revealed the propagation of stress waves in hollow cavities and the formation process of cavities at the bottom of the blastholes. Xie et al. [12] adopted the RHT model in LS-DYNA to analyse damage mechanisms of cut blasting under high in-situ stresses and proposed a modified design method applicable to deep rock mass according to numerical optimization simulation. Hu et al. established a model of wedge cut blasting by combining Auto-CAD with ANSYS and employed LS-DYNA to simulate cavity extension. In terms of temporal and spatial evolution with multidirection cutaway views, the applicability of wedge cut blasting for roadway excavation was studied [13]. Although considerable efforts have been made in past decades to understand the cavity formation mechanism of cut blasting, few studies have considered the whole process of cut blasting from charge initiation to cavity formation. Compared with open-pit blasting, cut blasting is more difficult due to the lack of free surface. The basic process of cut blasting is to break rock, throw fragments, and form a cavity. The hypotheses about cavity formation in wedge cut blasting in the literature are not very reasonable. In addition, the whole process of wedge cut blasting has not been comprehensively considered. Numerical simulation is commonly limited to the research on stress propagation and rock failure, and the consideration about fragments casting is insufficient. Therefore, by considering the characteristics of cut blasting and the effects of stress waves and detonation gas, it is meaningful to investigate the whole process of wedge cut blasting, which includes rock breaking and fragments casting.

Research on postblast fractures suggests that stress waves from the detonating explosive contribute to rock fracturing, while the detonation gas is responsible for moving or throwing fractured rock [14]. Recent studies have revealed that stress waves are responsible for the initiation of the crushed zone and surrounding radial fractures, while the gas pressure further extends the fractures [15]. In this paper, based on the theory of the coaction of stress waves and detonation gas, the evolution process of wedge cut blasting is divided into two stages. In the first stage, assuming that rock breaking caused by stress waves occurs, the objective is to establish a computational model of the rock failure zone by considering the rock dynamic strength. In the second stage, a criterion is proposed for cavity formation based on the quasi-static loading of detonation gas, assuming that the driving force from the detonation gas needs to overcome the total resistance from the reserved surrounding rock mass, accelerate fragments, and then throw fragments to form a cavity. The cavity formation mechanism is investigated via modelling and the theoretical analysis from the initiation of the cylindrical charge to the formation of the cavity.

2. Failure Zone from Coupled Cylindrical Charge Explosion

After charge detonation, the rock mass around a borehole is usually divided into three areas based on the failure degree of the rock: the crushed zone, the fractured zone, and the elastic vibration zone. In the near area, the violent explosion shock wave causes a small range of rock crushing, i.e., the crushed zone. Then, the shock wave rapidly attenuates into stress wave, which generates cracks in the middle area, that is, fractured zone I. The stress wave later weakens into seismic wave, which arouses elastic vibration in the far area, known as the elastic vibration zone. Although the detonation gas penetrates into cracks and elongates them, forming fractured zone II, the detonation gas plays a much weaker role in rock fragmentation than the stress waves. In cut blasting, we focus on the crushed zone and fractured zone I caused by the shock wave and the stress wave, respectively, and the combination of these two zones is called the failure zone. The combined effects of explosion shock wave and stress wave are collectively called stress waves in this paper. The partitions of rock blasting around a borehole are shown in Figure 1 [16].

2.1. Blasting Load

After detonation of a coupled cylindrical charge, the blasting impact load directly exerts a force on the hole wall. Because of the lack of feasible measuring methods, empirical formulas are usually adopted to estimate the detonation pressure [17]. Based on the equation of the state of ideal gases, Fickett and Davis [18] and Henrych et al. [19] proposed a simplified formula to calculate the Chapman–Jouguet (CJ) detonation pressure from a cylindrical charge:where is the CJ detonation pressure, is the explosive density, is the detonation velocity, and is the specific heat ratio. This equation gives the well-known expression when assuming [20].

The impact load arouses shock wave in the vicinity of the borehole, according to acoustic approximation theory [21]:where is the initial pressure of the shock wave on the borehole wall, is the longitudinal wave velocity of the rock, and is the rock density.

The shock wave causes rock crushing near the charge, which consumes a large amount of energy; thus, the wave rapidly attenuates into stress wave and spreads to the surrounding rock mass, which is subjected to complex stress loading. Since the longitudinal (i.e., axial) dimension of the borehole is significantly larger than the transverse (i.e., radial and circumferential) dimensions and all sections perpendicular to the axis of the borehole are the same, both the radial and circumferential stresses of the rock mass and the constraint conditions hardly change along the axis at the moment of explosion. Thus, this situation can be simplified to a plane strain problem, and the stress at a point in the rock mass can be expressed aswhere , , and are the radial stress, circumferential stress, and axial stress, respectively; is the distance from the point to the axis of the cylindrical charge; is the borehole radius; is the coefficient of attenuation, taken as in the shock wave zone and in the stress wave zone; and is the coefficient of the lateral confining pressure, taken as , where is the dynamic Poisson’s ratio of the rock [21].

The dynamic Poisson’s ratio is an important parameter in the numerical analysis of the rock dynamic response. The dynamic Poisson’s ratio is related to the load, moisture content, loading rate, and stress state. In some rock engineering applications with limited field data, a value between 0.2 and 0.3 is a common estimate for Poisson’s ratio [22]. Generally, the static Poisson’s ratio of intact rocks can be directly determined in the laboratory using static tests, and because the dynamic Poisson’s ratio is slightly lower than the static Poisson’s ratio, the value of the dynamic Poisson’s ratio can be estimated by the empirical equation , where is the static Poisson’s ratio of the rock [21, 23, 24].

2.2. Rock Failure Criterion

During blasting, the rock around the borehole is in a state of three-dimensional stress mixed with tension, compression, and shear. This complex stress state of the surrounding rock can be expressed at any point by the equivalent stress [21]:

When blasting, the rock is subjected to instantaneous impact load. Only when the dynamic strength of the rock is adequately considered can the characteristics of rock blasting be truly reflected. The crushed zone is caused by rock dynamic compression, while fractured zone I is the result of tensile failure. According to the failure characteristics and the von Mises criterion, the rock failure criterion is represented aswhere and are the dynamic compressive strength and the dynamic tensile strength of the rock, respectively.

The strain rate effect of rock-like material strength has been generally recognized and accepted by researchers in related fields. According to a study of the strain rate threshold value [25], rock strength presents significant rate-dependent characteristics when the strain rate exceeds 5 × 10−4 s−1. The strain rate range of engineering blasting is approximately 100–103 s−1 [26], and the blasting load belongs to dynamic load with a high strain rate. The dynamic strength of granite increases with the strain rate under a high strain rate, and the rock dynamic strength has a linear relationship with the 1/3 power of the strain rate [27, 28], as follows:where and are the static compressive strength and static tensile strength of the rock, respectively, and is the strain rate, which is higher in the shock wave zone () and lower in the stress wave zone ().

2.3. Failure Zone of Blasting

The crushed zone generated by shock wave iswhere is the radius of the crushed zone and .

The fractured zone I is outside the crushed zone, and the interface between these two zones is in a critical state, divided into compression failure and tensile failure. The stress on the interface [21] can be described bywhere is the radial stress on the interface between the crushed zone and fractured zone I. Therefore, the fractured zone I resulting from the effect of the stress wave iswhere is the radius of fractured zone I.

The rock failure zone is generated by a coupled cylindrical charge and contains the crushed zone and fractured zone I; thus, it is easy to know that the radius of the failure zone is equal to .

3. Cavity Formation in Cut Blasting

In phase one, rock is subjected to the dynamic loading induced by stress waves, which typically has a very short duration of a few milliseconds. In phase two, detonation gas penetrates into the broken rock mass, has work done on rock fragments, and throws them. The cavity formation mechanism in wedge cut blasting is investigated by means of a three-dimensional simplified model of the cut rock mass [8].

3.1. Model of Wedge Cut Blasting

A model of wedge cut blasting is shown in Figure 2, where numbers 1–6 represent two rows of inclined boreholes arranged symmetrically; four capital letters constitute a surface, e.g., ACFE represents the free surface and CJNF represents the right lateral side; is the hole length; is the angle of inclination; is the hole spacing of the same row; is the space at the top of symmetrically inclined holes, called hole-top spacing for short; and is the space at the bottom of symmetrically inclined holes, named hole-bottom spacing for short.

The dynamic loading of stress waves causes rock mass damage in the first stage of cut blasting, and the rock mass is separated into rock fragments of various sizes. Therefore, we treat broken rock as a granular media material that has lost cohesion but reserved internal friction in the subsequent analysis of cavity formation. When detonation gas works in the second stage, it needs to overcome the resistance from the surrounding rock mass, accelerate fragments, and throw fragments to form the cavity. Then, the Mohr–Coulomb strength criterion can be adopted to calculate the shearing resistance on the four lateral sides. A schematic of the simplified force analysis for the 3D cut entity is shown in Figure 3. The cut entity is excavated in a semi-infinite geologic body, which is subjected to an in situ stress field of (vertical component) and (horizontal component). The shearing resistance on each of the two lateral sides, AIJC and EMNF, is expressed as ; expresses the shearing resistance on each of the two lateral sides, AIME and CJNF. expresses the net force of quasi-static pressure from inclined boreholes, which is perpendicular to the free surface pointing outwards.

3.2. Cavity Formation Mechanism

According to the Mohr–Coulomb strength criterion, the shearing resistance on each of the two lateral sides, AIJC and EMNF, can be estimated bywhere is the internal friction angle; is the vertical stress, taken as , where is the volume weight of overburden, and is the buried depth; and , where is the gravitational acceleration constant. Similarly, the shearing resistance on each of the two lateral sides, AIME and CJNF, is given bywhere is the number of boreholes and is the horizontal stress, taken as [29].

It should be noted that cohesion appears in neither of these two equations for the calculation of shearing resistance because the broken rock is regarded as a granular media material that has lost cohesion but reserved internal friction.

After the explosion, boreholes are full of detonation gas. According to the quasi-static theory, the force exerted on the rock mass [8] follows the equationwhere is the perpendicular force exerted on the lateral side, is the quasi-static pressure, and is the charging length. The theoretical borehole pressure is approximately 45% of the CJ detonation pressure, assuming complete reaction at the front [17]. According to Persson et al. [20], the borehole pressure with a coupled cylindrical charge can be estimated by . The net force from inclined boreholes perpendicular to the free surface pointing outwards is shown below:

In the second stage of cut blasting, rock casting and cavity formation occur when the loading force of the detonation gas overcomes the total resistance from the reserved surrounding rock mass and accelerates fragments, i.e., subject to the following constraint:

4. Application to a Project

4.1. Computational Solution

A large-scale connection tunnel in hard rock, with a buried depth of 230 m, is investigated as an example. The physical and mechanical properties of the rock are listed in Table 1. According to the equations mentioned in Section 2, the crushed zone and fractured zone are obtained in granite caused by No. 2 rock emulsion explosive with and . These two zones extend approximately 2.9 times and 18 times the borehole radius; i.e., their scope is approximately 1 to 2.9 times and 2.9 to 18 times the borehole radius, respectively.

The borehole is 45 mm in diameter; i.e., the borehole radius is 22.5 mm. Therefore, the crushed zone and fractured zone caused by the cylindrical charge explosion extend approximately and , respectively, and the range of the rock failure zone is approximately 405 mm. The hole spacing of the same row is designed on the principle of connecting the failure zones of adjacent boreholes, and the spacing should satisfy the condition , taken as here. The determination principle of hole-bottom spacing is usually based on connecting the crushed zones of symmetrical boreholes, and often, the value is no more than 30 cm. Nevertheless, the hole-bottom spacing is too small to cause the bottom of boreholes to overlap, and rock is excessively crushed in engineering applications. When determining the hole-bottom spacing, it is necessary to ensure that the failure zones connect, i.e., , taken as here. Figure 4 shows the borehole layout of the wedge cut. The charging length is 200 cm, and the blasting parameters are presented in Table 2.

When the relevant parameters are substituted into the equations in Section 3.2, the loading force and resistance are derived. The total force on the four lateral sides perpendicular to the free surface pointing inwards is approximately 109.09 MN, which resists rock casting and cavity formation. The net loading force from 8 inclined boreholes perpendicular to the free surface pointing outwards is approximately 488.94 MN. Thus, equation (14) is shown to be useable.

4.2. Field Test

Field tests were performed at the Huizhou underground water-sealed oil storage caverns project, which is located in the coastal region of southern China, as shown in Figure 5. There are no active and regional faults near the project site. According to geologic survey, field drilling, and laboratory tests, the strata are mainly composed of medium-fine grained granite with good rock integrity and nondeveloped joint fissures. The physical and mechanical parameters of the granite are shown in Table 1. The connection tunnel has a width of 12 m, height of 6 m, and length of 150 m, and its cross section has a straight-wall arch profile with an area of approximately 62 m2; the borehole layout and blasting parameters are identical to those introduced in Section 4.1.

Excavation is carried out according to the design. A jumbo driller with three arms is adopted and assisted by manual operation. The cut holes are charged with No. 2 rock emulsion explosive, and a continuous charging structure is used. The blasting parameters are listed in Table 2. Adequate stemming could greatly improve the efficiency of blasting, and a water-soil composite stemming material is used. Figure 6 shows the relevant materials and instruments of the field tests.

Fine mechanical construction guarantees the quality of the field tests, and the task of obtaining a series of statistics has been completed. The utilization rate of boreholes from all twelve field tests is more than 89%, and the average utilization rate is 91.6%, as summarized in Table 3. The advance per round maintains a high level, which is approximately 2.7-2.8 m. The tunnelling face presents a decent flatness after blasting, and the rock has no noticeable outward bulge. Because a detailed fragment size distribution in underground blasting is difficult to obtain, we have chosen three reference values for the fragments to obtain a general fragment size distribution of the broken rock, as shown in Figure 7. The fragments from 20 to 50 cm account for approximately 60% of fragments. The large block ratio is close to 2%, and the largest fragment size is no more than 80 cm. Figure 8 shows some details about the field tests. These results indicate that wedge cut is good for enhancing the advance per round, accelerating excavation, and achieving good blasting performance.

4.3. Discussion

Although the coalescence of stress waves and detonation gas is widely accepted to cause rock fracture and fragmentation during rock blasting, the features of rock breaking still vary considerably because of significant differences in the physical and mechanical properties and structural characteristics of engineering rock mass. Stress wave loading plays a dominant role in rock fracturing and is key to rock crushing and initial crack formation. Related numerical results showed that a model that considers only stress waves could make reasonable predictions of the rock mass response to blasting load [15, 3035]. The effect of rock fragmentation caused by stress waves is more obvious in hard rock.

The crushed zone and the fractured zone appear around the borehole when a coupled charge explosion occurs in an infinite rock, and the scopes of these two zones are approximately 3–7 times and 8–150 times the borehole radius, respectively [36]. Many models have been proposed to estimate the scope of the zones around a borehole. Hustrulid suggested that the did not exceed 3–5, and Esen et al. indicated that the range of was between 1.3 and 6.6 [17]. According to Ханукаев, the crushed zone was approximately 2-3 times the charging radius and the fractured zone was approximately 10–15 times. Dai suggested that the crushed zone and the fractured zone were approximately 1.6–2.7 times and 13–17 times the charging radius, respectively, and these two zones caused by No. 2 rock emulsion explosive in granite were approximately 3.7 times and 24 times, respectively [21]. Tang et al. suggested that the crushed zone and the fractured zone were approximately 2.5–6.0 times and 16–26 times the charging radius, respectively, and these two zones caused by TNT explosive in granite were approximately 1.7 times and 12.9 times, respectively [37]. In this paper, the crushed zone and the fractured zone caused by No. 2 rock emulsion explosive in granite are approximately 2.9 times and 18 times the borehole radius, respectively; that is, the scope of the two zones is approximately 1–2.9 times and 2.9–18 times the borehole radius, respectively. Although the crushed zone and the fractured zone derived from the proposed model are only considered stress waves, these two zones are consistent with the results from the models in the literature. Although this choice may cause the calculations of the failure zone to be slightly low, the impact is very small. Thus, it has been verified that the proposed model shows good applicability for predicting failure zones. The model could predict the blast-induced failure zone caused by the coupled cylindrical charge well in wedge cut blasting. In addition, this model provides theoretical support to determine cut parameters, such as hole spacing and hole-bottom spacing.

The hole-bottom spacing is taken as 50 cm, which is larger than conventional values. As spacing increases, more energy is used to generate and extend cracks; thus, a wider fractured zone can be acquired rather than a wider crushed zone. The increase helps avoid borehole overlap during drilling and excessive rock crushing in blasting. With a larger spacing, the number of boreholes decreases for the same cross section, which saves time and reduces resource consumption. Therefore, although the value of hole-bottom spacing is larger than conventional values, we could achieve better blasting results and greater efficiency.

Some previous studies indicate that detonation gas does not contribute significantly to fracture propagation and is responsible for moving or throwing rock fragments [14, 36, 38]. In the second stage of blasting, detonation gas penetrates into the broken rock mass, swells, and has work done on fragments. The quasi-static force of detonation gas is as high as 489 MN, while the resistance is approximately 109 MN. Therefore, this force can sufficiently overcome the resistance from the reserved surrounding rock mass, and fragment casting occurs after fragments speed up. Then, a cavity forms that provides another free surface and rock swelling space for subsequent blasting. However, if the detonation gas is strong enough that excess energy remains after overcoming the resistance, the gas may have too much work done on the fragments, leading to a higher velocity, farther casting distance and more difficult mucking. Past investigations on the cavity formation mechanism of cut blasting frequently paid more attention to the propagation law of stress waves, while neglecting fragments casting and cavity formation. In this paper, when analysing the formation of a cavity, the broken rock mass is regarded as granular media that have lost cohesion but reserved internal friction under the impact of stress waves. The formation mechanism of the cavity is reasonably explained from a mechanical perspective of detonation gas.

The average utilization rate of boreholes from twelve field tests is as high as 91.6%, and the volume of the cut cavity is approximately 6.5 m3. Thus, the specific charge is approximately 2.9 kg·m−3, which means that slightly more explosives are consumed. The rock is well fragmented, with uniform fragment sizes and hardly any large blocks, and the cavity is well formed. These results mean that the rock mass has been broken previously because of stress wave loading before the detonation gas throws fragments to form the cavity. Field tests have verified the feasibility of the cut design based on rock breaking and fragment casting caused by stress waves and detonation gas, respectively, and the blasting results indicate the validity of modelling and analysing the cavity formation mechanism of wedge cut blasting.

5. Conclusions

In this paper, by considering the characteristics of cut blasting and the combined effects of stress waves and detonation gas, a theoretical model of the whole process from charge initiation to cavity formation is proposed and then applied to engineering practice. The conclusions can be summarized as follows:(1)Dynamic loading from stress waves plays a dominant role in rock breaking. The scope of the rock failure zone caused by the coupled cylindrical charge in granite is approximately equal to 18 times the borehole radius according to theoretical calculation. The hole-bottom spacing of the wedge cut can be designed as 50 cm, which is significantly larger than the conventional value.(2)Detonation gas is responsible for rock fragment casting. When the driving force from quasi-static pressure overcomes the total resistance from the reserved surrounding rock mass and accelerates fragments, fragment casting occurs, leading to cavity formation.(3)Considering the coalescence of stress waves and detonation gas together during rock blasting, the evolution process of wedge cut blasting has been divided into two stages, which contains the compression damage and the tensile damage caused by stress waves and the fragment casting caused by detonation gas. For wedge cut blasting, the theoretical model makes an overall interpretation of the mechanism.

This study is based only on theoretical analysis and a specific case, and more blasting tests are needed to provide theoretical support and technical guidance for wedge cut blasting in hard rock.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This research was financially supported by the National Natural Science Foundation of China (Nos. 51674062 and 51474049), the National Science Foundation for Young Scientists of China (No. 51604061), and the Basic Scientific Research Operating Expenses of Central University (No. N160104009). We wish to thank Guangdong People Blasting Engineering Co. Ltd. and The First Engineering Co. Ltd. of CTCE Group for supporting the field tests and assisting with data collection. Special thanks are due to American Journal Experts for its professional English editing service.