Research Article  Open Access
Theoretical Calculation and Experimental Analysis on Initial Shock Pressure of Borehole Wall under Axial Decoupled Charge
Abstract
This paper attempts to calculate the exact initial shock pressure of borehole wall induced by the blasting with axially decoupled charge. For this purpose, Starfield superposition was introduced considering the attenuation and superposition of blasting pressure, and the theoretical solution of initial borehole wall pressure was obtained for the upper and middle airdecked charging structures. Then, the explosive pressure field around the borehole was measured by cement mortar models and a dynamic pressure test system, and the pressures at multiple measuring points were simulated with numerical models established by ANSYS/LSDYNA. The results show that the deviations between simulated and theoretical pressures are smaller than 10%, indicating the reliability of the theoretical formula derived by Starfield superposition. For the upper airdecked charging structure, the initial shock pressure of the charging section followed a convex distribution, with the peak value near the charge centre. With the increase in the distance from the charging section, the borehole wall shock pressure in the air gap underwent a sharp decline initially before reaching a relatively constant level. The minimum pressure was observed at the hole collar. For the middle airdecked charging structure, the pressures at both ends of the charging section obeyed a convex distribution, with the peak value near the charge centre. Finally, the author optimized airdecked charging structure of periphery boreholes within Grade III surrounding rocks of Banjie tunnel, China, and proved the enhancement effect of the theoretical findings on smooth blasting. The research findings provide valuable references to the theoretical and experimental calculation of air column length and other key parameters of airdeck blasting and shed new light on the charging structure determination of smooth blasting and blasting vibration control for the excavation of largesection, deep mining roadways.
1. Research Status
During the excavation of deep roadway or tunnel under high crustal pressure, the stability of the surrounding rock is dependent on the effect of smooth blasting, which itself hinges on the selection of axial decoupling coefficient. Concerning the close ties between the coefficient, the air gap, and the initial shock pressure of the borehole wall, below is a brief review of the existing studies on mine blasting, especially airdeck blasting, that involves these three elements.
In previous research, the initial shock pressure of borehole wall with airdecked charge is considered as 8~11 times the quasistatic pressure of the detonation gas [1], provided that the borehole wall is rigid and subject to orthogonal impact from the detonation wave and that the shock waves which attenuate and superpose along the borehole axis are negligible; i.e., the initial shock pressure of borehole wall is uniformly distributed in the axial direction. This conclusion may apply to short air gap or axially coupled charge but does not suit long air gap or axially decoupled charge. Hence, it is necessary to explore the initial shock pressure of borehole wall under axially decoupled charge.
To disclose the mechanism of airdeck blasting, Kabwe [2] suggested that the surrounding rock can be enhanced by repeated oscillation of shock wave at air space sections and verified the suggestion with the distribution of rock fragments after blasting with top airdecked charge (hereinafter referred to as “top airdeck blasting”) at Chimiwungo Pit of Lumwana, an openpit copper mine. Similarly, Jhanwar and Jethwa [3], through an airdeck blasting at an openpit coal mine in India, discovered that the blasting method can reduce the yield of fines and boulders, and that the length of air space directly bears on the distribution of rock fragments. Fourney et al. [4] concluded that, in top airdeck blasting, the shock wave is reflected upon reaching the plug, and the shock pressure on the surrounding rock lasts 2~5 times that under coupled charge. Moreover, it is critical to determine a proper length of air gap, because the surrounding rock will not be further broken if the pressure of the reflected shock wave fails to surpass the tensile strength of the rock.
Yang et al. [5] numerically simulated the distribution of initial shock pressure at different air gap ratios, pointing out that the optimal directional fracture effect appears at the air gap ratio of 33.3%~50% in slitcharge blasting. Saharan et al. [6] proposed to improve the energy utilization rate of explosives through rational setting of air gap and plug lengths and achieved the expected improvement with different combinations of air gap and plug. Furthermore, Park and Jeon [7] put forward an airdeck blasting method with a thin papertube and quantified its shock absorption effect and tunnelling efficiency through numerical simulation and experiments. Roy et al. [8] applied the mature airdeck blasting technique of openpit mining and tunnelling to an underground coal mine in India and succeeded in enhancing the blasting effect and controlling the deflagration and harmful gases.
Over the years, much research has been done on the pressure distribution of borehole wall. Despite the lack of direct measuring methods, the borehole wall pressure has been estimated with various empirical formulas or detonation theories [9]. Taking borehole wall pressure as the basis of blasting design, Cunningham [10] derived borehole wall pressure via polynomial decay and presented an alternative to the pressure. Otuonye et al. [11] calculated borehole wall pressure through lab tests on stemming movement. Homment et al. [12] determined borehole wall pressure based on the expanding volume of blasting cavity. Swoboba [13] developed a novel numerical model on the propagation of explosion pressure, which describes the crack propagation in the borehole wall and deduces borehole wall pressure based on blasting cavity or crack volume. Likewise, Liu and Katsabanis [14] created a new computing method after exploring the initial shock pressure of borehole wall. Yilmaz and Unlu [15] provided the estimation formula of borehole wall pressure in light of the explosion pressure expressions under decoupled charge in [16, 17].
With the development of computer technology and intelligent testing, a series of new approaches have been introduced to reveal the dynamic effect of rock blasting, e.g., numerical simulation and model test. With the aid of LSDYNA and FLAC3D, Jiang et al. [18] investigated the blasting damage of vertical crater retreat (VCR) mining method, laying the theoretical basis for rational control of stope boundaries and optimization of blasting parameters. Through numerical simulation of bottom airdeck blasting, Wu et al. [19] studied the dynamic pressure features and failure mechanism of nearfield boreholes, compared the dynamic pressure features of different detonation methods, and found that indirect initiation can easily destroy the bottom rock of borehole by enhancing static function intensity.
Focusing on directional pressurerelief blasting, Xiao et al. [20] carried out cement mortar model tests on the dynamic strain of borehole wall and other parts. The test results are in good agreement with those of numerical simulation, indicating that the pressurerelief effect of the borehole can protect the wall materials and the bottom air gap columns. Ma et al. [21] tested 4 types of millisecond blasting models with threecircle boreholes, aiming to find the blast damage on the surrounding rock in vertical shaft excavation and put forward effective damage control measures. Zhao et al. [22] performed a blasting model test with similar materials and dosages by Froude scaling and provided theoretical and empirical references to the selection of materials and charging conditions for blasting test models.
The above studies open a new direction for quantifying the distribution of initial shock pressure in borehole wall. Considering airdeck blasting, Yang et al. [23] probed into the borehole pressure distribution under different charging conditions and disclosed the heterogeneous pressure distribution across borehole wall: the initial shock pressure peaks at the centre of the charging section, far higher than that in the air gap. This means the initial shock pressure of borehole wall distributes unevenly along the axial direction. Nevertheless, there is no report on the mathematical expressions of the uneven distribution pattern.
Under coupled/decoupled charge, Ling [24] tested the initial shock pressure of cement mortar and organic glass and obtained the initial shock pressure at each point of borehole wall. Using manganesecopper piezoresistor, Ni et al. [25] measured the peak initial shock pressure produced by ammonium nitrate explosive, emulsion explosive, and watergel explosive in granite and concrete and conducted a regression analysis of the relationship between borehole wall pressure and the incident angle. Their research shows the correlation between the initial shock pressure and the distance to the explosive source along the borehole axis. Nonetheless, the finding is only suitable for singlestage charging. Further research is needed to ascertain the distribution of borehole wall pressure of double or multistage charging.
2. Introduction
This paper aims to select a proper axial decoupling coefficient for airdeck blasting. To this end, it is necessary to investigate the air gap and the initial shock pressure of the borehole wall under axial decoupled charge. Here, the Grade III surrounding rock of Banjie tunnel, a 4,806mlong deepburied tunnel, is taken as the object. The tunnel is one of the three main tunnels in YongrenGuangdong section of ChengduKunming railway. In the tunnel, the surrounding rock masses are mostly Grade III~V sandstone and sandy mudstone. The tunnel segments through Grades III, IV, and V rocks are, respectively, 4,115m, 530m, and 161m in length. Thus, the object tunnel segment accounts for 85.6% of the total length of Banjie tunnel. The physicalmechanical parameters of Grade III surrounding rock are listed in Table 1.

The previous smooth blasting of Grade III surrounding rock was carried out in two steps with the following parameters: the total sectional area=116.02 m2, the excavation height= 5.81m (upper stage) and 5.52m (lower stage), and the cyclic advancement=3~3.2m. In the upper stage, the construction parameters are as follows: borehole distance= 50cm (periphery holes) and 80cm (auxiliary holes); and borehole depth=3.8m (periphery holes) and 4m (auxiliary holes). For the fivestage double horizontal wedge cutting, the vertical hole is 4m in depth. The density and detonation velocity of 2# emulsion explosive are 1,300kg/m3 and 3,200m/s, respectively. Due to the high detonator cost and complex process, the periphery holes were charged continuously at the bottom (hereinafter referred to as the continuous bottom charging).
The above blasting structure has some problems that led to an unsuccessful tunnelling. After the smooth blasting on the upper stage, the maximum overbreak was as high as 0.4m near the working face (hole bottom) of the hance, and the maximum underbreak near lining (hole collar) stood at 0.25m. At the vault, the postblast outline exhibited as a large flat plate with severe rockfall. In general, the working face had an uneven contour and barely any visible hole profiles (Figure 1).
The unsuccessful blasting in Banjie tunnelling is mainly caused by inappropriate charging. When continuous charging structure is used at the hole bottom, the over concentrated charging in this section is prone to bring about heavy shock pressure to the borehole wall. The uncharged hole collar, however, tends to result in an underbreak due to insufficient shock pressures.
The unsuccessful blasting is mainly attributable to improper charging. With continuous bottom charging, the charged bottom applied a heavy shock pressure onto borehole wall, while the uncharged hole collar generated so few shock pressure as to cause underbreak. In light of these, Starfield superposition was introduced considering the attenuation and superposition of blasting pressure, and the theoretical solution of the initial borehole wall pressure was obtained for top and middle airdeck blasting.
Meanwhile, the author constructed a stable dynamic pressure test system that avoids the interference and signal distortion of the previous test method and validated the theoretical analysis results through model tests. To determine the optimal axial decoupling coefficient, the distribution features of the initial borehole wall pressure were analysed under different axial decoupling coefficients and contrasted with each other by numerical simulation.
Based on the findings, the peripheral hole charging structure was optimized for the blasting of Grade III rocks in Banjie tunnel, aiming to enhance the smooth blasting effect. Then, the damping effect of airdeck blasting was compared with that of the blasting with continuous bottom charging.
3. Theoretical Analysis
3.1. Starfield Superposition Method
Starfield superposition treats the column charge as the superposition of a finite number of spherical charges with equal radius, that is, the equal charging principle. Let d_{c} and r_{e} be the diameter of column charge and the equivalent radius of spherical charge, respectively. Then, the equal charging principle can be expressed as
For simplicity, it is assumed that the peak pressure pattern remains the same independent of the detonation order of spherical charges in each unit, such that the peak pressure at a point in the air gap equals the pressure of the shock wave at the point produced by the blasting of the equivalent unit spherical charge lying the closest to the air gap. This value is close to the actual peak pressure and thus satisfies the computing demand [26].
By Starfield superposition, the time effect of the explosion of unit spherical charge cannot be ignored while computing the blasting effect of the column charge. This is because both the detonation velocity of explosives and the longitudinal wave velocity of rock are both in the 10^{3} m/s order of magnitude. Hence, the following exponential function was adopted to depict the time attenuation of the shock wave induced by unit spherical charge [27]:where A (A=C_{p}(12μ)/[a(1μ)k]) is the pressure attenuation coefficient of equivalent unit spherical charge A=C_{p}(12μ)/[a(1μ)k] [28, 29]; C_{p} is the longitudinal wave velocity of rock; a is a constant related to rock properties; and t is the decay time of shock wave.
For a point in the air gap, each equivalent unit spherical charge has a pressure effect on the time of positive pressure. Without the loss of generality, the peak pressure is assumed to occur when the shock wave induced by a unit spherical charge arrived at a specified point. For example, the pressure of a point in the air gap reaches the peak when the shock wave induced by spherical charge unit k propagates to that point, indicating that the pressure induced by the preceding unit i at this point has attenuated. The decay time can be expressed aswhere l_{i} is the distance from unit spherical charge i to the specified point in the air gap; l_{k} is the distance from unit spherical charge k to the specified point in the air gap; D is the detonation velocity of the explosive.
In accordance with the shock wave theory, the peak pressure of shock wave attenuates with the distance, under a single spherical charge, following the pattern below [30]:where is a scaled distance; B is a constant; α is the attenuation coefficient of peak shock wave pressure along the borehole axis. The scaled distance is calculated as , with R being the distance between the calculation point and the explosive source, and Q being the TNT equivalent of actual charge quantity (Q=Q_{s}, where Q_{s} is the actual charge quantity; is the specific detonation heat of the explosive; is the specific detonation heat of TNT).
According to (1), the actual charge quantity of a single spherical charge can also be expressed as
Thus, (4) can be rewritten aswhere ρ, , B, and α are constants for specific explosive and surrounding rock.
For better accuracy, the attenuation coefficient of shock pressure along the borehole axis is denoted as δ. Assuming that K=0.98α(ρB, the peak pressure of a single equivalent unit spherical charge at the specified point can be obtained as
When ignoring the interaction of explosion for adjacent unit spherical charges, the peak pressure at the specified point of the column charge can be approximately considered as the superposition of the peak stress at the point of all spherical charges.
3.2. Initial Shock Pressure of Borehole Wall with Upper AirDecked Charge
In the upper airdecked charging structure, the length of the air gap is assumed as l_{a}, the charge length at hole bottom as l_{e}, and the charge diameter is d_{c}. As shown in Figure 2, a random point A in the air gap is x (x≤l_{a}) away from the top of the charging section. At the blasting of the unit spherical charge i, the shock pressure of the spherical charge at point A can be expressed as
In this case, the shock pressure produced by the entire section of column charges at point A equals the pressure superposition of n unit spherical charges at that point:
3.3. Initial Shock Pressure of Borehole Wall with Middle AirDecked Charge
In the middle airdecked charging structure, the length of the air gap is assumed as l_{a}, the charge length at hole bottom as l_{e}, and the length of upper charging section and the lower charging section are l_{1} and l_{2}, respectively. As shown in Figure 3, a random point A’ in the air gap is x (x≤l_{a}) away from the top of the lower charging section. At the blasting of the unit spherical charge i, the shock pressure of the spherical charge at point A’ can be expressed as
At the blasting of the unit spherical charge j, the shock pressure of the spherical charge at point A’ can be expressed as
In this case, the shock pressure produced by the upper and lower charging sections at point A’ equals the pressure superposition of all unit spherical charges at that point:where n_{1 }and n_{2} are the number of equivalent unit spherical charges in the upper and lower charging sections, respectively.
The variation curve of the initial shock pressure along the borehole axis of the upper and middle charging structures were derived by (9) and (12) (Figure 4).
4. Experiment Verification
4.1. Pressure Test System
Based on the principle of resistance strain, the pressure test system applies a dynamic pressure on the specimen, and the resulting deformation of the specimen is recorded by the strain gauges at the measuring points. Then, the deformationinduced resistance change is converted into the change of voltage or current, making it possible to deduce the value of deformation. Here, the pressure at each monitoring point is approximated by Hooke’s law, with the aim of reflecting the distribution features of the initial shock pressure of borehole wall along the axial direction.
The test instruments include a BlastUltra multichannel shock tester (Chengdu Tytest Co., Ltd.), a KD6009A strain amplifier (Yangzhou Kedong Electronics Co., Ltd.), etc. The strain gauges are attached to prefabricated strain bricks and embedded in a concrete model to receive the explosive signals. The entire test system is illustrated in Figure 5.
4.2. Preparation of Cement Mortar Models and Strain Bricks
The 40cm×20cm×50cm (L×W×H) cement mortar model was casted with 42.5# ordinary portland cement and screened fine sands (size: <1mm) at the mix proportion of 1:1:0.5 (cement: sand: water). The model was cured for 28 days at room temperature. The holes (depth: 40cm; blasting burden: 5.5cm) were reversed by a solid fiberglass pipe (OD: 12mm) at 10cm away from the front boundary and the back boundary.
For the upper airdecked charging structure, one 2# Nonel rock detonator was installed at the hole bottom, a 28cm air gap was reserved at the upper part, and the axial decoupling coefficient K_{l} was set to 5. For the middle airdecked charging structure, one 2# Nonel rock detonator was installed at the hole bottom and the hole collar, respectively, a 21cm air gap was reserved at the middle part, and the axial decoupling coefficient K_{l} was set to 2.5. In both charging structures, the hole collars were blocked with 5cm long cement plugs (Figure 6).
(a) Upairdeck charge
(b) Middleairdeck charge
The strain bricks are 3cm×3cm×40cm (L×W×H) cuboids. These were prepared with the same mix proportion, seeking to prevent the reflection of instantaneous explosion signals and maintain a uniform wave impedance between the model and the bricks [31, 32]. To capture explosion pressure, two strain sticks were arranged in parallel with the holes by a spacing of 2cm. One of them has 8 measuring points and the other has 11 (Figure 6). On each measuring point, two BX1204AA resistance strain gauges (resistance: 120Ω; sensitivity: 2.08±1%; Yangzhou Kedong Electronics Co., Ltd.) were pasted vertically via halfbridge connection.
Before pasting the strain gauges, the specimen surface was polished with sandpaper at 45° to the axis of the strain brick to remove the sands and gravels. The strain gauges were them pasted with strong glue. The redundant glue must be squeezed out to ensure the good contact between the gauge and the stick.
In the meantime, three standard 5cm×5cm×10cm (L×W×H) specimens were produced and cured for 28 days in the same environment as the models. The physicalmechanical parameters of the models were determined after the curing (Table 2).

4.3. Experimental Result
The experimental parameters are as follows: sampling rate=4 each/min, data collection time=10ms, negative delay=1ms, trigger level=5%, gain=100, bridge voltage=2V, and the lowpass frequency=1kHz. The shock pressures of upper and middle airdecked charging structures are recorded in Table 3.
 
Note. Instead of actual pressure values, the peak pressures obtained by Hooke’s law are a rough representation of the distribution features of initial borehole wall shock pressure in the axial direction. The relative distance from hole bottom is the quotient of the distance from hole bottom and the hole diameter. 
4.4. Discussion
Figures 7 and 8 display the test results of the upper and middle airdecked charging structures, respectively. It is clear that the initial shock pressure along the borehole axis obeys the same distribution pattern as that in Figure 4.(1)For the upper airdecked charging structure, the initial shock pressure of the charging section followed a convex distribution, with the peak value near the charge centre. With the increase in the distance from the charging section, the borehole wall shock pressure in the air gap underwent a sharp decline initially before reaching a relatively constant level. The minimum pressure was observed at the hole collar.(2)For the middle airdecked charging structure, the pressures at both ends of the charging section obeyed a convex distribution, with the peak value near the charge centre. By contrast, the pressure in the air gap was exhibited as a concave distribution and minimized at the middle of the air column.(3)In both the upper and middle airdecked charging structures, the initial borehole wall shock pressure increased with the decrease in the distance from the charging section. When that distance was on the rise, the initial pressure experienced a gradual decrease. The decline rate was fast at the beginning and slow on the later stage. In the end, the borehole wall shock pressure reached a relative stable state.(4)In the middle airdecked charging structure, the pressure values at major monitoring points in the air gap were close to the minimum value. This feature was particularly prominent when the air space was fairly long. Besides, the initial shock pressure in the upper and lower charge sections obeyed the same distribution of that in the upper airdecked charging structure. The initial shock pressure distribution of the two charge sections could be combined into the superposed impact pressure.
To sum up, the axial distribution curve of the initial borehole wall shock pressure in the two charging structures agrees well with that of the theoretical analysis.
5. Numerical Simulation
5.1. Verification of Theoretical Formula for Initial Shock Pressure
5.1.1. Numerical Model and Parameters
Five numerical models were established to verify the accuracy of the theoretical formula derived by Starfield superposition, which calculates the blastinduced initial shock pressure in the air gap. The numerical simulation considers both upper and middle airdecked charging structures, as well as five different charging lengths. In these models, the column charge was 1~4 times longer than equivalent spherical charge. For the middle airdecked charging structure, the lengths of the upper and lower charging sections are both 1~5 times of the diameter of equivalent spherical charge.
As shown in Figure 9, the column charge is 1.2cm in diameter, which puts the diameter of the equivalent spherical charge as 1.46cm according to (1), the length of the reserved air gap is 5cm, and the borehole burden is 5.5 cm. Moreover, the borehole size is Φ1.2 cm×[(1~5)×1.46]cm (the length of column charge)+5 cm (the length of air gap) for the upper airdecked charging structure, and Φ1.2 cm×[2×(1~5) ×1.46]cm (the length of column charge)+5 cm (the length of air gap) for the middle airdecked charging structure. The initiation points were arranged along the charging section at an interval of 1.46cm and designed to initiate simultaneously.
For better accuracy, the rock, explosive, and air were described separately by MATPLASTICKINEMATIC, MATHIGHEXPLOSIVEBURN, and MATNULL models, and the state equations of the explosive and the air were illustrated by EOSJWL and EOSLINEARPOLYNOMAIAL, respectively [33–35]. The parameters of each material model and its state equation are provided in Tables 4–6. Concerning the boundary conditions, the left side and the top of the model were set as free faces, and the other planes as nonreflecting boundaries. Specifically, the explosive and air were processed with ArbitraryLagrangianEulerian (ALE) algorithm, while the rock was treated by Lagrangian algorithm [36, 37].

 
Note: A, B, R_{1}, R_{2}, and ω are material model parameters; R_{1} and R_{2} are nondimensional parameters; ω is Grüneisen constant (change rate of pressure relative to internal energy under constant volume); E_{0} is initial internal energy per unit volume of explosive; is relative volume. 
 
Note: C_{0} ~ C_{6} are multinomial coefficients of the state equation; E_{0} is initial internal energy per unit volume of explosive; is relative volume. 
5.1.2. Result Analysis
The intersection points between the borehole wall and the axial crosssection of air column, which is far away from free faces, were selected to monitor the initial shock pressure of borehole wall. Through simulation, two sets of initial shock pressure values were obtained from the monitoring points for upper and middle airdecked charging structures. Meanwhile, the shock pressures at the specified points were acquired through theoretical analysis.
Figure 10 compares the theoretical and simulated values of initial shock pressures at the specified points under the two charging conditions. It can be seen that the simulated pressures deviated from the results of (9) by 8.82% and those of (12) by 8.96%. Since neither of the deviations surpassed 10%, the simulated results agree well with the conclusion of theoretical analysis, an evidence to the reliability of the theoretical formula derived from the Starfield superposition under both upper and middle airdecked charging conditions.
Furthermore, the theoretical and simulated initial shock pressure at specified points of the air gap increased with the charging length, according to the variation in initial shock pressure with column lengths at specified points. The trend echoes with the classical blasting theory.
5.2. Verification on Axial Distribution Pattern of Initial Borehole Shock Pressure
5.2.1. Numerical Model and Parameters
The author established the numerical models by ANSYS/LSDYNA, using the same dimension, borehole size, and borehole pattern with the test model (Figure 11). For the upper airdecked charging structure, five models were created with different axial decoupling coefficients: K_{l}=2.5, K_{l} =3.5, K_{l} =4, K_{l}=5.0, and K_{l} =6.0. Among them, K_{l} =5.0 corresponds to the charging structure of the test model. In this charging structure, the overall charging length equals the length of the lower charging section.
For the middle airdecked charging structure, five models were created with different axial decoupling coefficients: K_{l}=2.5, K_{l} =3.5, K_{l} =4, K_{l}=5.0, and K_{l} =6.0. Among them, K_{l} =2.5 corresponds to the charging structure of the test model. The charging parameters of these models are shown in Table 7. In this charging structure, the overall charging length equals the sum of the length of the upper and lower charging section.

The material models of the rock, explosive, and air and the relevant state equations (Tables 4–6) are the same with those in the verification of theoretical formula for the initial borehole wall shock pressure. The stemming material is described by the MATSOILANDFOAM model, and its parameters are shown in Table 8.
 
Note: A_{0}, A_{1}, and A_{2} are constants of yield function; PC is truncation pressure of tensile failure; EPS1~ EPS10 are characteristic bulk strains; P_{1} ~ P_{10} are pressures corresponding to characteristic bulk strains. 
5.2.2. Results and Discussion
To capture the initial borehole wall shock pressure, row elements were selected from borehole wall as pressure monitoring points along the axial direction from hole bottom to hole collar. Figures 12 and 13 depict the distribution curves of the initial shock pressure along the borehole axis under different decoupling coefficients, respectively.
As shown in Figures 1213, the simulated pressure distributions of the two charging structures are similar to those obtained by theoretical analysis. In particular, the curves at K_{l }= 5.0 of the upper airdecked charging structure and K_{l }=2.5 of the middle airdecked charging are consistent with the theoretical results and the experimental findings. Suffice it to say that the simulation demonstrates the reliability of theoretical analysis and physical experiments.
6. Case Study and Comparative Model Test of Blasting Vibration
6.1. Case Study
The borehole length in a single blasting row is relatively long in Grade III surrounding rock of Banjie tunnel. To ensure the effect of smooth blasting, the peripheral hole was designed with airdecked charging structure. In the middle airdecked charging structure, the initial borehole wall shock pressure obeyed a wavy distribution (peak value> valley value>0). The peak and valley values, respectively, correspond to the charging section and the air gap. To achieve a consistent impedance of the explosive and rock, it is necessary to appropriate the length of the air column [38].
In Banjie tunnel, the compressive strength and tensile strength of Grade III surrounding rock are 75MPa and 5.6MPa, respectively (Table 1). According to the simulated distribution of initial borehole wall shock pressure in middle airdecked charging structure (Figure 13), the initial shock pressures in the air gap were all higher than the tensile strength of the rock when the axial decoupling coefficient varied in the range of K_{l }=2.5~6. This means blasting cracks formed among the adjacent peripheral holes. However, the shock pressures at major monitoring points in the air gap were above the compressive strength of the rock, when K_{l}=2.5, =3.5, K_{l} =4.0, and K_{l} =5.0. The result shows that the rock was crushed but not sufficient to achieve a good smooth blasting effect.
In contrast, when K_{l} =6.0, the pressures at most monitoring points in the air gap were close to the valley value, except the charging section and a small portion of the air gap, and were below the compressive strength of the rock. In this case, the surrounding rock remained integrated after the blast and kept a halfhole profile, revealing enhanced smooth blasting effect. Therefore, the axial decoupling coefficients of the peripheral hole were designed as K_{l} =6.0 in Grade III surrounding rock of Banjie tunnel. In light of the field conditions, fourstage charges and three airgaps were adopted for tunnelling. The optimized charging structure of peripheral hole in Grade III rock is given in Figure 14.
In Figure 14, the hole profile is clearly visible after blasting, and over 90% of halfholes are reserved. According to the blasting effect in Figure 15, there is barely any overbreak or underbreak after excavation, and the surrounding rock remained stable for a long time.
6.2. Comparative Model Test of Blasting Vibration
In tunnel blasting, the vibration, a key impact factor of rock stability, is the combined effect of the blasting of cutting holes, auxiliary holes, and peripheral holes. However, it is difficult to accurately extract the variation induced by the blasting of peripheral holes in field monitoring. To compare the blasting vibration of continuous bottom charging with airdecked charging structure, two 100cm×20cm×60cm (L×W×H) cement mortar models, denoted as model I and model II, respectively, were prepared for blasting experiment.
The mix proportion, curing environment, and curing time were exactly the same as the cement mortar models of Section 4.2. One borehole (depth: 40cm; blasting burden: 5.5cm) was reversed at 10cm away from the front boundary and the back boundary. With the same explosive payload, Models I and II, respectively, adopt continuous bottom charging and airdecked charging. The blasting was monitored by sensors placed at 10cm, 20cm, 40cm, and 80cm away from borehole. Figure 16 illustrates the cement mortar models and the layout of vibration measuring points.
The monitoring results are shown in Table 9, and the typical waveforms of vibration vector synthesis velocities are given in Figure 17. As can be seen in Table 9, the vibration velocity of model II was much lower than that of model 1 at the same distance from borehole, revealing the good damping effect of airdeck blasting. Compared with continuous bottom charging, the airdecked blasting structure had a 15.05% lower average amplitude under the same explosive payload. Thus, the airdecked blasting structure can reduce the integer blasting vibration velocity during the tunnelling in Grade III surrounding rock of Banjie tunnel and, in turn, assure the stability of the surrounding rock.

7. Conclusions
(1)Based on Starfield superposition, the theoretical solutions of the initial pressure on borehole wall were obtained for upper and middle decked charged structures. Through field monitoring and numerical simulation, it is learned that the initial borehole wall pressure was unevenly distributed at long air columns.(2)For the upper airdecked charging structure, the initial shock pressure of the charging section followed a convex distribution, with the peak value near the charge centre. With the increase in the distance from the charging section, the borehole wall shock pressure in the air gap underwent a sharp decline initially before reaching a relatively constant level. The minimum pressure was observed at the hole collar. For the middle airdecked charging structure, the pressures at both ends of the charging section obeyed a convex distribution, with the peak value near the charge centre.(3)The deviations between simulated and theoretical pressures are smaller than 10%, indicating the reliability of the theoretical formula derived by Starfield superposition. It is concluded that this formula works well in computing the blastinginduced initial shock pressure in the air gap under both upper and middle airdecked charging conditions.(4)Compared with the continuous bottom charging, the airdecked charging is suitable to the actual conditions thanks to the distribution features of initial borehole wall shock pressure. The application of the structure in Grade III surrounding rock of Banjie tunnel led to an expected smooth blasting effect.(5)The research findings provide valuable references to the theoretical and experimental calculation of air column length and other key parameters of airdeck blasting and shed new light on the charging structure determination of smooth blasting and blasting vibration control for the excavation of largesection, deep mining roadways.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research work was supported by the National Natural Science Foundation of China under Grants 51679093 and 51374112.
References
 J.H. Chen, J.S. Zhang, and X.P. Li, “Study of presplitting blasting parameters and its application based on rock blastinginduced damage theory,” Rock and Soil Mechanics, vol. 37, no. 5, pp. 1441–1450, 2016. View at: Publisher Site  Google Scholar
 E. Kabwe, “Improving collar zone fragmentation by top airdeck blasting technique,” Geotechnical and Geological Engineering, vol. 35, no. 1, pp. 157–167, 2017. View at: Publisher Site  Google Scholar
 J. C. Jhanwar and J. L. Jethwa, “The use of air decks in production blasting in an open pit coal mine,” Geotechnical and Geological Engineering, vol. 18, no. 4, pp. 269–287, 2000. View at: Publisher Site  Google Scholar
 W. L. Fourney, D. B. Barker, and D. C. Holloway, “Model studies of explosive well stimulation techniques,” International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, vol. 18, no. 2, pp. 113–127, 1981. View at: Google Scholar
 G. L. Yang, R. S. Yang, C. Huo, and Y. L. Che, “Numerical simulation of airdeck slotted charge blasting,” Advanced Materials Research, vol. 143144, pp. 787–791, 2011. View at: Publisher Site  Google Scholar
 M. R. Saharan, M. Sazid, and T. N. Singh, “Explosive Energy Utilization Enhancement with AirDecking and Stemming Plug, "SPARSH",” in Proceedings of the ISRM European Rock Mechanics Symposium, EUROCK 2017, vol. 191, pp. 1211–1217, Czech, June 2017. View at: Google Scholar
 D. H. Park and S. W. Jeon, “Reduction of blastinduced vibration in the direction of tunneling using an airdeck at the bottom of a blasthole,” International Journal of Rock Mechanics and Mining Sciences, vol. 47, no. 5, pp. 752–761, 2010. View at: Publisher Site  Google Scholar
 S. K. Roy, R. R. Singh, R. Kumar, and U. K. Dey, “Studies into possible use of air decking in solid blasting in Indian underground coal mines,” Transactions of the Institutions of Mining and Metallurgy, Section A: Mining Technology, vol. 117, no. 2, pp. 83–92, 2008. View at: Google Scholar
 S. Esen, I. Onederra, and H. A. Bilgin, “Modelling the size of the crushed zone around a blasthole,” International Journal of Rock Mechanics and Mining Sciences, vol. 40, no. 4, pp. 485–495, 2003. View at: Publisher Site  Google Scholar
 C. Cunningham, “Concepts of blast hole pressure applied to blast design,” Fragblast, vol. 10, no. 12, pp. 33–45, 2006. View at: Publisher Site  Google Scholar
 O. Francis, “Measurements and prsesdictions of borehole pressure variations in model blasting system,” in Proceedings of the 1st International Symposium on Rock Fragmentation by Blasting, pp. 7–20, Lulea, Sweden, 1983. View at: Google Scholar
 P. J. Hommert, J. S. Kuszmaul, and R. L. Parrish, “Computational and experimental studies of the role of stemming in cratering,” in Proceedings of the 2nd International Symposium on Rock Fragmentation by Blasting, pp. 550–562, Colorado, Keystone, 1987. View at: Google Scholar
 G. Swoboba, “The numerical modelling of blast loading,” in Proceedings of the International Symposium on Application of Computer Methods in Rock Mechnics and Engineering, Xian, China, 1993. View at: Google Scholar
 L. Liu and P. D. Katsabanis, “A theoretical approach to the stress waves around a borehole and their effect on rock crushing,” in Proceedings of the 4th International Symposium on Rock Fragmentation by Blasting, pp. 10–13, Vienna, Austria, 1993. View at: Google Scholar
 O. Yilmaz and T. Unlu, “Three dimensional numerical rock damage analysis under blasting load,” Tunnelling and Underground Space Technology Incorporating Trenchless Technology Research, vol. 38, no. 9, pp. 266–278, 2013. View at: Publisher Site  Google Scholar
 S. Nie and M. Olsson, “Study of fracture mechanism by measuring pressure history in blastholes and crack lengths in rock,” in Proceedings of the 27th Annual Conf. Explosives and Blasting Technique, Orlando, Florida, USA, 2000. View at: Google Scholar
 A. Persson, R. Holmberg, and J. Lee, “Rock blasting and explosives engineering,” International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, vol. 32, no. 6, p. 278A, 1995. View at: Google Scholar
 N. Jiang, C. B. Zhou, X. D. Luo, and S. Lu, “Damage characteristics of surrounding rock subjected to VCR mining blasting shock,” Shock and Vibration, vol. 2015, Article ID 373021, 8 pages, 2015. View at: Publisher Site  Google Scholar
 L. Wu, Y. Zhou, and Z. Guo, “The rock failure mechanism of airdecked blasting in hole bottom,” Applied Mechanics and Materials, vol. 201202, pp. 375–378, 2012. View at: Publisher Site  Google Scholar
 D. J. Xiao, B. Li, C. J. Pu, and H. Q. Zhou, “Model test and numerical simulation for directional pressure relief blasting,” Advanced Materials Research, vol. 779, pp. 848–856, 2013. View at: Publisher Site  Google Scholar
 Q. Ma, P. Yuan, J. Zhang, R. Ma, and B. Han, “Blastinduced damage on millisecond blasting model test with multicircle vertical blastholes,” Shock and Vibration, vol. 2015, Article ID 504043, 6 pages, 2015. View at: Publisher Site  Google Scholar
 H. Zhao, X. P. Li, and Q. Dong, “The research of similar material and dosage in blasting model test,” Applied Mechanics and Materials, vol. 724, pp. 152–155, 2015. View at: Publisher Site  Google Scholar
 G.L. Yang, R.S. Yang, and L.L. Jiang, “Pressure distribution along borehole with axial airdeck charge blasting,” Explosion and Shock Waves, vol. 32, no. 6, pp. 653–657, 2012. View at: Google Scholar
 W. M. Ling, “Experimental research on explosion pressure on the wall of a borehole in rock,” Mining and Metallurgy, vol. 13, no. 4, pp. 13–16, 2004. View at: Google Scholar
 Z. F. Ni and Y. M. Li, “Theoretical calculation and regression on initial shock wave pressure of borehole in rock blasting,” Chinese Journal of Rock Mechanics and Engineering, vol. A1, pp. 511–514, 1996. View at: Google Scholar
 X. H. Jin, Y. Sha, and R. Wang, “Analysis on breakage from blasting crater of column charge,” Water Resources and Hydropower Engineering, vol. 41, no. 11, pp. 48–51, 2010. View at: Google Scholar
 W. L. Wang, Drilling and Blasting, China Coal Industry Publishing House, Beijing, China, 1984.
 G. Harries, “The modeling of long cylindrical charge of explosive,” in Proceedings of the 1st International Symposium on Rock Fragmentation by Blasting, pp. 419–431, Lulea, Sweden, 1983. View at: Google Scholar
 W. B. Lu, The Propagation of Stress Wave and Its Effects in Rock Blasting, Wuhan University of Hydraulic and Electric Engineering, Wuhan, China, 1994.
 J. Henrych and G. R. Abrahamson, “The dynamics of explosion and its use,” Journal of Applied Mechanics, vol. 47, no. 1, p. 218, 1980. View at: Publisher Site  Google Scholar
 Q. Zong, Tunneling Blasting Parameters Model Experiment Study on Vertical Well Freezing Soil, University of Science and Technology China, Hefei, China, 2004.
 J. J. Ma, Z. Z. Xiong, W. D. Duan, D. W. Zhong, and L. J. Cai, “Theoretical testing study factors affecting parallel hole cut blasting,” Journal of University of science and Technology (Natural science edition), vol. 24, no. 2, pp. 170–174, 2001. View at: Google Scholar
 Y. F. Chen and G. X. Yi, “Dynamic response analysis of the reinforced concrete column under the effect of explosive impact load,” Advanced Materials Research, vol. 681, pp. 99–104, 2013. View at: Publisher Site  Google Scholar
 I. A. Balagansky and A. A. Stepanov, “Numerical simulation of Composition B high explosive charge desensitization in gap test assembly after loading by precursor wave,” Shock Waves, vol. 26, no. 2, pp. 109–115, 2016. View at: Publisher Site  Google Scholar
 P. Yan, W. Zhou, W. Lu, M. Chen, and C. Zhou, “Simulation of bench blasting considering fragmentation size distribution,” International Journal of Impact Engineering, vol. 90, pp. 132–145, 2016. View at: Publisher Site  Google Scholar
 Y. Han and H. Liu, “Finite element simulation of mediumrange blast loading using LSDYNA,” Shock and Vibration, vol. 2015, Article ID 631493, 9 pages, 2015. View at: Publisher Site  Google Scholar
 D. Su, Y. Kang, D. Li, X. Wang, and F. Yan, “Analysis and numerical simulation on the reduction effect of stress waves caused by water jet slotting near blasting source,” Shock and Vibration, vol. 2016, Article ID 5640947, 18 pages, 2016. View at: Publisher Site  Google Scholar
 W.B. Gu, Z.X. Wang, J.H. Chen, J.Q. Liu, and M. Lu, “Experimental and theoretical study on influence of different charging structures on blasting vibration energy,” Shock and Vibration, vol. 2015, Article ID 248739, 11 pages, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Xiaoming Lou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.