The Effect of Blast-Hole Arrangement, Delay Time, and Decoupling Charge on Rock Damage and Vibration Attenuation in Multihole Blasting

Rock damage and vibration attenuation are the basis of blasting design, which will aﬀect rock breaking results and the safety of structures (buildings). In this paper, the eﬀect of blast-hole arrangement, delay time, and decoupling charge on the rock damage and the vibration attenuation in multihole blasting were numerically investigated. Through dynamic analysis software ANSYS/LS-DYNA, the JOHNSON_HOLMQUIST_CONCRETE (JHC) rock model and ﬂuid-solid coupling method were used to establish single-hole and multihole rock blasting models. Based on the analysis of single-hole rock damage and vibration, ﬁrstly, the eﬀect of the above three factors on the rock damage characteristics of multihole blasting was analyzed. Then, the extracted peak particle velocity ( PPV ) data of multihole blasting were ﬁtted to the United States Bureau of Mines (USBM) equation to obtain the blasting vibration attenuation parameters (the site constant K and α ). It is found that the blast-hole arrangement, delay time, and decoupling charge have a much smaller eﬀect on α than K. At the same time, a delay-time-dependent model of PPV correction coeﬃcient was proposed, and its rationality was veriﬁed by ﬁeld data. The results obtained in this paper have reference signiﬁcance for optimizing the blasting design and improving the blasting eﬀect.


Introduction
As an important construction method, engineering blasting is widely used in production and construction activities such as mining engineering, hydraulic engineering, tunnel engineering, and nuclear power engineering [1].During the blasting process, a small part of the explosion energy is used to break the rock, causing the fragmentation and cracking of rock mass in the near zone, that is, the rock mass damage.
e range of the near zone of the blasting source is affected by the lithology, explosive performance, and blasting methods, and it is estimated to be approximately 2-33/100 times the charge radius [2].e near-zone damage characteristics caused by the explosion have an important influence on the subsequent grinding, crushing, or slagging efficiency and energy utilization.e remaining part of the explosion energy will be dissipated in the form of shock waves, thermal energy, and blasting seismic waves, accompanied by adverse effects such as vibrations, shock waves, fly rocks, toxic gas, and noise [3].Among such effects, blast-induced vibration is considered the most adverse effect.e blast-induced vibration can reduce the stability and bearing capacity of the rock mass, resulting in damage or even destruction of buildings (structures).In the study of blast-induced vibration effects, usually, the peak particle velocity (PPV) is used as the vibration variable to evaluate the intensity of blast-induced vibration and control blasting safety [4][5][6].e attenuation law of blasting vibration is one of the main bases for blasting design.Researchers from various countries have established some different empirical equations ( 7)- (11), as listed in Table 1.
Due to the anisotropy, nonlinearity, variability of mechanical properties of rock materials, and the coexistence of solid, liquid, and gas phases, the response of rock mass under explosive load is very complex [12].Many scholars and engineers have conducted long-term research on rock blasting damage and blast-induced vibration.Based on the tensile expression method of the classic blasting damage model, Hu et al. [13,14] constructed a new tensile-compressive damage model, which can quantitatively calculate the damage.e model was linked to LS-DYNA to analyze the excavation damage effects of smooth blasting and presplit blasting on high slopes.Numerical simulation results show that the blasting damage of smooth blasting has significant cumulative characteristics.Based on field tests and numerical simulations, Li et al. [15] studied the explosive load characteristics, the blast-induced damage characteristics for rock mass, and the PPV attenuation law and established the relationship between PPV at 30 m away from the charge hole and the damage depth of rock mass.A rock damage control method is proposed to restrict the rock damage depth under blasting excavation according to the safety threshold of vibration velocity for rock mass.Khandelwal and Singh [16] conducted multiple sets of field blasting tests and monitored the ground vibration values and then studied the artificial neural network method (ANN) to predict PPV.By comparing with the prediction results of PPV empirical formulas in various countries, it proved that ANN has a better prediction level.
e research on single-hole blasting law is representative, but the actual engineering blasting is mainly multihole blasting.erefore, the study of the multihole blasting effect is of great significance to engineering practice.In multihole blasting, the proper hole arrangement will make the blasting energy evenly distributed, improve the crushing result of the rock, reduce the large lump rate, and increase the looseness.Based on field survey data, Zhou et al. [17] used numerical simulation and physical model test methods to simulate and optimize the four typical cut hole arrangements (9-hole, single-spiral-hole, double-spiral-hole, 2-empty-hole-diamond) of shaft formation by one deep-hole blasting, and the relevant blasting parameters obtained were successfully applied to field tests.rough numerical simulation of rock slope presplit blasting under in-situ stress, Wang et al. [18] compared the blasting results of rock mass with the square and plum-blossom-shape arrangement of blast-holes.e results show that the plum-blossom-shape arrangement can improve the blasting quality under prestress conditions.Delayed blasting has an important effect on rock damage and blast-induced vibration.With the gradual popularization of electronic detonators, the precise control of delayed blasting has attracted more attention from researchers [19][20][21].Based on field data, Chen et al. [22] studied the influence of millisecond delay time, detonating sequence, and properties of the seed signal on the vibration reduction effect of millisecond blasting.Based on ANSYS/LS-DANY finite element software, Shao et al. [23] established a three-dimensional tunnel blasting model.e effects of different delay times on the dynamic response of tunnel lining are quantitatively studied, and combined with the PPV safety criterion, a reasonable millisecond delay time is obtained.Qiu et al. [24] theoretically analyzed the formation mechanism of blasting craters in the near-field and the mechanism of vibration reduction in the far-field under the condition of millisecond blasting.A comparison of field tests was carried out, and it was found that only when the delay time is less than the formation of a new free surface, it is possible to form a normal blasting crater.Compared with simultaneous blasting, delayed blasting can effectively reduce the PPV in the near-field.Using LS-DYNA, Hashemi and Katsabanis [25] simulated delayed blasting in the rock mass and explored the influence of initiation time on damage generation and morphological distribution.It is found that the delay time provides favorable conditions for the growth of cracks in the rock mass adjacent to the blast hole.Johansson and Ouchterlony [26] conducted small-scale tests on short-delay blasting.e results show that when the delay time is in the time range of the stress wave interactions, there is no distinct difference in the rock fragmentation compared with no shock wave interactions.Tatsuya et al. [27] simulated the relationship between PPV at certain monitoring points and delay time through computer programs and selected the best delay time.A design concept based on combined delayed blasting is proposed, which can effectively reduce the PPV of important monitoring points.
e decoupling charge is a charge structure often used in controlled blasting.Blasting with decoupling charge can greatly reduce underexcavation, overexcavation, and side pumice, improve construction efficiency, and have an important effect on rock damage and blast-induced vibration.rough numerical simulation, Shao et al. [28] studied the law of rock crack propagation in smooth blasting.e calculation formula of charge decoupling coefficient in smooth blasting is proposed, which can reduce the damage and improve the stability of surrounding rock in smooth blasting.Qi [29] used software ANSYS to simulate rock blasting with decoupling charges and analyzed the corresponding stress response values and nephogram under different radial decoupling coefficients.It is found that the rock mass blasting stress response is the largest in the case of a coupled charge and gradually decreases with the increase of the decoupling coefficient, indicating that the air medium can buffer the effect of detonation waves.Wang et al. [30] Table 1: Different conventional blast vibration equations [7].

√
) α Ambraseys-Hendron [10] PPV � K(W 1/3 /R) α Bureau of indian standards (BIS) [11] PPV � K(W/R 2/3 ) α R is the distance-usually radial-from the measurement point to the blast face (m), W is the maximum charge per delay (kg), and the K and α are nondimensional coefficients related to the topographic and geological conditions between the explosion source and the measuring point, respectively.
2 Shock and Vibration also reached the same conclusions.
rough theoretical analysis and blasting model tests of decoupling charge blasting, Lou et al. [31] studied the formation and propagation of shock waves, as well as the impact pressure characteristics of blast-hole walls.Results show that when the decoupling coefficient is 1.5-3.5, the impact pressure is evenly distributed, the impact time is longer, and the blastinduced vibration is weaker.With the increase of the decoupling coefficient, the initial impact pressure shows a trend of first increasing and then decreasing.
Existing studies have used the methods of theoretical analysis, field test, numerical simulation, and model test to study rock blasting damage and blast-induced vibration.Among them, numerical simulation is a common method for studying blasting, which can easily control blasting conditions such as the number, size, arrangement, and initiation time of blast-holes, and has the advantages that other methods cannot replace.Among the numerical methods, the large-scale general explicit dynamic finite element analysis software LS-DYNA is widely used in the dynamic response of rock engineering under impact loads.LS-DYNA has two methods for blasting analysis: the Lagrange method and the multimaterial fluid-solid coupling method.In the Lagrange method, the grid and the analysis structure are integrated, and the finite element nodes are material points, and clear material interfaces can be obtained.However, serious distortions will occur during the blasting, leading to calculation errors.In the multimaterial fluid-solid coupling method, explosives and other fluid materials (air, water, etc.) use Euler or ALE (Arbitrary Lagrange-Euler) algorithm, and solid materials use the Lagrange algorithm, and the two interact through fluid-solid coupling.In numerical modeling, grid division is an important step.e quantity, uniformity, and division quality of model grids have a great influence on calculation accuracy, scale, and stability.Under the condition that the computer hardware performance can support the calculation, it is necessary to ensure the regular shape of the units, and the parts of concern and with large deformation should be appropriately refined.
e outline of the work given in this paper is as follows.First, we establish single-hole and multihole rock blasting models by software ANSYS/LS-DYNA in Section 2. en, we analyze the single-hole blasting damage and vibration effect in Section 3, followed by the analysis of multihole rock damage in Section 4 and the analysis of multihole vibration attenuation characteristics in Section 5. Finally, we conclude our paper in Section 6.

Models and Materials
e ANSYS/LS-DYNA numerical simulation software was used to simulate the cylindrical charge blasting.To avoid calculation errors caused by the distortion of the model grid, the multimaterial fluid-solid coupling method was adopted.
at is, the explosive and air materials were assigned to the ALE algorithm, and the rock material was assigned to the Lagrange algorithm, and then the two algorithms were coupled.

Model Design.
According to the different analysis focus, the limitation of the model grids number and the highquality grid's requirements, when studying the rock damage in the near zone and the blasting vibration effect in middle and far zone, the analysis models of the near zone and middle and far zone are established, respectively.e size of the near-zone damage analysis model is 2.5 m × 2 m × 6 m (XY × Z direction).To better observe the damaged area, dense meshing is carried out, and the element size is set small.e size of the middle and far zone analysis model is 30 m × 4 m × 6 m (X × Y × Z direction), and the grid is sparser than that in the near area.At the same time, 5 vibration measurement points are set at different horizontal distances (R) from the blast-hole: R � 5 m, 7.6 m, 11.6 m, 17.7 m, and 27 m.In all numerical models, the charge radius is 0.035 m, the blast-hole radius is 0.045 m (adjusted according to the decoupling coefficients in the analysis of decoupling charge blasting).e depth of the blast-hole is 5 m, of which the stemming length is 2 m, and the charge height is 3 m.e detonation points are all set as the geometric center of the explosive.Nonreflecting boundary conditions are applied to the bottom and surroundings of the models to avoid the influence of reflected waves.e upper surface of the model is free.Section a is perpendicular to the Z-axis, and section b is perpendicular to the Y-axis, both passing through the center point of the explosive.e single-hole analysis model considering near zone damage is shown in Figure 1, and the single-hole analysis model considering middle and far zone vibration is shown in Figure 2.
When building the multihole blasting model, based on the single-hole model, only the number of blast-holes, the hole arrangements, the delay times, and the decoupling coefficients were changed.Numerical simulations were then carried out, and the results of single-hole blasting rock damage and vibration were taken as comparison standards to analyze the influence of hole arrangement, delay time, and decoupling charges on the blasting effect.

Explosive.
e explosive material is defined by the MAT_HIGH_EXPLOSIVE_BURN keyword and the Jones-Wilkins-Lee (JWL) state equation are used to define the detonation state, and the keyword is JWL_EOS. is state equation expresses the relationship between the detonation pressure P e and the relative volume V e of the detonation product and the internal energy per unit volume E e : where A e , Be, W, R1, R2 are constants related to the properties of explosives, V e0 is the initial relative volume, and E e0 is the initial internal energy per unit volume.e related parameters of the explosive are listed in Table 2, where VoD represents the detonation velocity, and PCJ represents the Chapman-Jouget pressure.
Shock and Vibration 3  Shock and Vibration

Air.
e air material adopts the MAT_NULL model, which can better simulate the fluid properties.e corresponding state equation is defined by the JWL_LI-NEAR_POLYNOMIAL keyword, as shown below: where P a represents the detonation pressure, V a represents the relative volume, V a0 represents the initial relative volume, E a represents the internal energy per unit volume, E a0 represents the initial internal energy per unit volume, C 0 , C 1 , C 2 , C 3 , C 4 , C 5 and C 6 are coefficients, and μ represents the relative density.e related parameters of the air are listed in Table 3.

Rock.
e rock material adopts the MAT_JOHN-SON_HOLMQUIST_CONCRETE (JHC) constitutive model, which can describe large strains, high strain rates, and high pressures.JHC constitutive model has many parameters.Among them, the density ρ and uniaxial compressive strength f c are based on laboratory test results.Shear modulus G is a combination of elastic modulus E and Poisson's ratio ] and is calculated according to the following formula: e rock Brazilian tensile strength (T) is calculated from f c according to the following formula [33]: e limit surface parameters (A, B, and N) have a great influence on the blast-induced vibration, which is determined by the basic rock parameters according to the following method [34,35].
e JHC normalized equivalent stress σ * is expressed as a function of pressure, strain rate, and damage: where A is the normalized cohesive strength, B is the normalized pressure hardening coefficient, N is the pressure hardening exponent, C is the strain rate coefficient, D is the damage parameter, P * is the normalized hydrostatic pressure, and _ ε * is the dimensionless strain rate.σ * , P * and _ ε * are defined as follows: where σ is the realistic equivalent stress, P is the realistic hydrostatic pressure, _ ε is the realistic strain rate, and _ ε 0 is the reference strain rate.
Without considering the effects of damage and strain rate, (5) is abbreviated as follows: Based on plastic theory, both JHC and Mohr-Coulomb models pass through the following two points on the compression meridian plane: uniaxial compression (σ * � A) and pure shear (σ � c). e following relationship can be obtained: Based on the Mohr-Coulomb model, the maximum principal stress σ 1 has the following relationship with the confining pressure σ 3 : According to the Mohr-Coulomb model, the maximum principal stress σ 1 has the following relationship with the confining pressure σ 3 , the cohesion c, and the internal friction angle φ: σ and P are expressed as follows: Substituting the laboratory data of triaxial compression into (7), the other two limit surface parameters of the JHC model (B � 2.298 and N � 1.0344) were obtained.

Definition of Rock Damage.
e damage variable D of the JHC model (0 ≤ D ≤ 1, D � 0 means intact rock, D � 1 means completely broken rock) is accumulated from the equivalent plastic strain Δε P and the plastic volume strain Δμ P , expressed as follows: where SFmax is the normalized maximum strength.
In addition, the relevant parameters of the JHC model also include the plastic strain before fracture (EFmin), the crushing pressure (p c ), the crushing volume strain (u c ), the compaction pressure p l , the compaction volume strain u l , and three material constants K 1 , K 2 and K 3 .ese parameters are less sensitive to the calculation results, and the initial parameter values of the JHC model can be used.e JHC model is used to simulate the slightly weathered granite and the related parameters of the JHC model are listed in Table 4.

Simulation of the Single-Hole Blasting
e model considering near zone damage (Figure 1) and the model considering middle and far zone vibration (Figure 2) of the single hole were calculated, and the rock damage and vibration attenuation results were obtained.

Rock Damage.
e damage nephogram in the near zone of the single-hole blasting is shown in Figure 3, and different colors represent different damage levels, that is, different D values.It can be seen from the damage nephogram of section a that the damaged area is centered on the charge and extends outward along the radial direction of the blast hole.
e damage level changes from heavier to lighter and finally dissipates in a scattering shape.From the damage nephogram of section b, it can be seen that the damaged area extends outward from the blast hole.e damage level of the initiation point is smaller than the upper and lower parts of the explosive, and the boundary of the damage zone is smoother than section a.
According to the research of Hu et al. [13] on the safety threshold of blasting damage, the critical damage value of retained rock mass is D � 0.19, so the area of D � 0.19-1 is taken as the damaged area.Calculated from the damage contour of the single-hole blasting model, the damaged area of section a is about 35 times of the blast-hole cross section, and the damaged area of section b is about 8 times of the blast-hole longitudinal section.

Blast-Induced Vibration.
e PPVs of the five measurement points of single-hole blasting are shown in Figure 4(a), the horizontal radial direction (X-direction) is the largest, the vertical direction (Z-direction) is second, and the horizontal tangential direction (Y-direction) is the smallest.e horizontal radial velocity-time curve at R � 5 m is shown in Figure 4(b).It can be seen that the vibration duration is about 50 ms, the part with large vibration amplitudes is basically within 20 ms.
Based on the USBM equation (Table 1), the single-hole numerical calculation data is subjected to regression analysis, and the horizontal radial attenuation parameters of the single-hole blasting are obtained.e fitting result is shown in Figure 5, and the attenuation equation is as follows:

Multihole Blasting Rock Damage
Based on the single-hole analysis model considering near zone damage (Figure 1), only the number of blast-holes, the blast-hole arrangements, the delay times, and the decoupling coefficient are changed to establish the multihole blasting models and conduct numerical simulations.

Hole Arrangement.
e hole arrangements often used in multihole blasting can be summarized into two types: rectangle and triangle.e rectangle and triangle arrangements represented by 4 blast-holes are shown in Figure 6.
Figure 7 shows the near-zone damage nephograms of blasting models with different hole arrangements under the same explosive quantity and the same initiation time.It can be seen that, compared with rectangular holes, the size of rock mass fragments under the triangular hole arrangement is more uniform and smaller.
e rectangular hole arrangement has a greater possibility of producing larger-sized fragments, which is not conducive to rock mass breaking and engineering production.
Calculating the damaged area in the range of D � 0.19-1 under the rectangular hole arrangement is 1.543 m 2 , accounting for 34.3% of the total cross-sectional area.e damaged area of the triangular hole arrangement is 1.581 m 2 , accounting for 35.1% of the total cross-sectional area.e damaged area of the triangular hole arrangement is 1.024 times that of the rectangle.at is, the breaking capacity of the triangle on the rock mass is slightly larger than that of the rectangle.Considering that the vibration duration is about 50 ms in single-hole blasting, the rock damage conditions when the delay time Δt is 0 ms, 1 ms, 5 ms, 10 ms, 15 ms, 20 ms, and 50 ms are simulated, respectively.Among them, Δt � 0 ms means that No. 1 and No. 2 explosives detonate at the same time.Figure 9 shows the damage nephogram when Δt is 0 ms, 5 ms, and 50 ms.

Shock and Vibration 7
Define the damage area ratio λ as follows: where   At this time, λ 2 � 1.406, which is an increase of about 40% compared with Δt � 0 ms.It can be considered that λ can increase by about 40% at most compared with simultaneous detonation.According to Stagg [36], the possible reason is as follows: when detonating at the same time, what is produced between explosives No. 1 and No. 2 is the superposition of stress waves, which breaks the rock.During the delay blasting, the No. 1 explosive that detonated first produces detonation gas, which did not dissipate within the short-delay time.e stress wave generated after the detonation of the No. 2 explosive interacts with the detonation gas produced by the No. 1 explosive, resulting in greater crushing power.at is, delayed blasting promotes the interaction between the stress wave induced by the latter blast-hole and the detonation gas generated by the previous blast-hole, and the interaction between the stress wave and the explosive gas can promote rock breaking result more than the superposition of the stress wave.

Decoupling Charge.
e charge structure is one of the factors that have a decisive influence on the blasting effect.In controlled blasting, decoupling charges are usually  Shock and Vibration used.e decoupling coefficient is an important control parameter for the decoupling charge blasting, which is divided into the radial decoupling coefficient and the axial decoupling coefficient.e radial decoupling coefficient refers to the ratio of the blast-hole diameter to the charge diameter, and the axial decoupling coefficient refers to the ratio of the blast-hole length to the charge length in the direction of the blast-hole length.
e construction operability of radial decoupling charges is better than that of axial decoupling charges.erefore, radial decoupling charges are often used in engineering practice.In the radial decoupling charge blasting, the air layer between the explosive and the blast-hole wall can reduce the shock wave pressure generated by the explosion, thereby reducing excessive damage to the rock mass and achieving better crushing results.e radial decoupling coefficient K r is as follows: where d b and d e are the diameter of the blast-hole and the diameter of the charge, respectively.e model is simplified as two blast-holes, and the decoupling coefficient K r is set to 1, 1.5, 2, 2.5, 3, and 5, respectively.K r � 1 is the decoupling charge under the limit condition, which is equivalent to the coupling charge.
e damage nephogram of blasting models with K r � 1, 2, and 5 is shown in Figure 11.According to the rock damage results, with the increase of K r value, the damaged area around and the connecting part of the blast holes gradually decreases.When K r � 1 (coupling charge), the explosive detonation wave acts on the rock directly, and the energy transfer rate is the highest.When K r > 1 and increases, it means that the air layer between the explosive and the rock is getting thicker, and the energy transfer rate, blast-hole pressure, and rock damage area are all reduced.Smooth blasting and presplit blasting use this principle to achieve the expected rock-breaking result.

Multihole Blast-Induced Vibration Attenuation
Based on the single-hole analysis model considering middle and far zone vibration (Figure 2), only the number of blastholes, the blast-hole arrangements, the delay times, and the decoupling coefficient are changed to establish the multihole blasting models considering middle and far zone vibration and conduct numerical simulations.

Hole Arrangement.
e PPVs of the measurement points in the middle and far zone under different hole arrangements are shown in Figure 12(a).ere is almost no difference in vertical vibration PPVs between rectangular and triangular holes.e main differences lie in horizontal radial and horizontal tangential vibrations: (1) In rectangular hole blasting, the horizontal radial PPV is larger than that of the triangular hole, and the difference is about 30 cm/s.( 2) In triangular holes, the horizontal tangential PPV is larger than that of the rectangular hole, indicating that the triangular hole arrangement is more likely to stimulate the generation of shear waves.erefore, in the blasting practice, the blasthole arrangement method should be selected according to the situation of protected items.For structures or instruments that are sensitive to radial vibrations, a triangular blast-hole should be selected, and for tangential vibrations, a rectangular blast-hole can be considered.( 3) At a far distance(≥30 m), the effect of the blast-hole arrangement on vibration tends to be weakened.
According to the USBM equation, the regression curves and optimum equations of horizontal radial blasting vibrations under different hole arrangements are shown in Figure 12(b).It can be seen that the site constants K and α of the rectangular-hole blasting are 60.062 and 1.7886, which are larger than those of the triangular-hole blasting (K � 45.699, α � 1.5883).

Delay Time.
Calculate the vibration PPV when the delay time Δt is 0 ms, 1 ms, 5 ms, 10 ms, 15 ms, 20 ms, and 50 ms.Figure 13 takes the X-direction blasting vibration waveforms at R � 5 m of Δt � 0 ms, 15 ms, and 50 ms as examples.e results show that with the increase of Δt, two peaks gradually appear.When Δt � 50 ms, the two wavebands are independent because 50 ms is close to the duration of a singlehole waveform.
Figure 14 shows the variation of horizontal radial PPVs with different delay times.It can be seen that the delayed blasts produce distinctly smaller PPVs than does the simultaneous blasting: (1) When Δt � 0 ms, that is, simultaneous blasting, the PPV value is the largest, and the vibration superimposition degree is the highest.
e reason is that millisecond blasting can reduce the charge quantity per delay, which allows the energy of the explosion to be released in stages.In the case of simultaneous blasting, the charge quantity is the total amount of 2 holes, while in the delayed blasting, the charge quantity per delay is the amount of 1 hole.e blast-induced vibration increases with the increase of the charge quantity.( 2) When Δt ≠ 0, with the increase of /s, and 3-4 cm/s, the difference is small in the far distance and large in the near.It indicates that the delay time has a small influence on the PPVs of long-distance measurement points but has a large influence on short-distance measurement points.Setting a large delay time is better for vibration control, but it will increase the overall vibration duration, which is unreasonable for engineering practice [23].erefore, the delay time needs to be selected reasonably in engineering practice, and it should not be too large or too small.e horizontal radial PPVs of two holes blasting with different delay intervals and the regression curves and optimum equations based on the USBM equation are plotted in Figure 15.When Δt � 1 ms, 5 ms, 10 ms, 15 ms, 20 ms, 50 ms, the coefficients of determination are 0.99897, 0.99899, 0.9994, 0.99987, 0.99991, 0.99915.When Δt � 0, it means simultaneous blasting, with the largest charge quantity W and the smallest K value.When Δt ≠ 0, the charge quantity W is half of that when Δt � 0. As Δt increases, the K value tends to decrease.e constant α is in a stable range without large fluctuations.

Correction Coefficient.
To quantify the relationship between the change of blast-induced vibration PPVs and the delay times, define ß as the correction coefficient of the radial vibration PPV when the delay time is Δt: where PPV single is the radial vibration PPV of the single-hole blasting.e empirical equations of blast-induced vibration with charge and propagation distance as independent variables have been well-known in the field of blasting engineering and have been widely used.Take the empirical USBM equation as an example, when Δt is constant: where K and α are the site coefficients of single-hole blasting, K Δt and α Δt are the site coefficients of millisecond blasting when the delay time is Δt. e single-hole blasting charge quantity is the same as the maximum charge quantity of delayed blasting, both are W. Substitute into equation (20): It can be seen that ß and R have an inverse proportional function relationship.Let 22) can be written as follows: If the duration of the single-hole blasting waveform is T, let Δt � nT; in this case, T � 50 ms, then n � Δt/50.When the distance R is constant, it can be seen from Figure 16 that as n increases, ß is approximately inversely proportional to n.Assuming that when Δt > 0, ß is in inverse proportion to R and n: where b 1 , b 2 , b 3 are constants.Take the logarithm of both ends of the equation ( 20) to get the following: MATLAB was used to fit equation ( 21) based on the least-square method.b 1 � 1.0666, b 2 � 0.04046, and b 3 � 0.0855 were obtained.
Obtaining equation ( 20) is as follows: e regression coefficient of determination is 0.8866, and the sum of squared error is 0.04431.

Equation Verification.
According to the field tests conducted by Chen et al. [22] when studying the effect of millisecond blasting to reduce vibration, the signal collected from one detonation was taken as the single-hole vibration waveform, and the signal collected from two detonations with a certain delay time was taken as the double-hole vibration waveform.e nonel detonator was used for detonation, and millisecond times of 50 ms and 100 ms were selected, which are the delay times of the two sets of tests.Two measurement points were arranged with a distance of 18 m and 30 m.
e single-hole charge quantity and the double-hole maximum charge quantity are both 10 kg.In the above numerical simulation of blasting, the single-hole charge quantity and the double-hole maximum charge quantity are both 11.6 kg.
e value of the two charge quantities is close.e PPV results of the field tests are listed in Table 5. e waveform duration T of the measured time history curve in single-hole blasting is 200 ms.
According to the data in Table 5, the rationality of equation ( 22) is verified.e distance R, the delay time Δt, and the waveform duration T of each measurement point were substituted into the equation ( 22) to obtain the calculated value of ß. e measured single-hole PPV values and the calculated values of ß were substituted into equation (16), and the predicted values of the double-hole PPV were obtained.
e relevant calculation results are listed in Table 6.e relative errors between the predicted Shock and Vibration values and the measured values are not very small, which is inevitable.Each site has its unique lithology characteristics, and the types and weights of explosives, charging methods, delay accuracy are also different.e correction coefficient ß is affected by these factors, but it can still meet the needs of the project.Δt=1 ms: PPV = 117.42(W 1/2 /R) 1.704  Δt=5 ms: PPV = 101.12(W 1/2 /R) 1.654  Δt=10 ms: PPV = 93.75(W 1/2 /R) 1.643  Δt=15 ms: PPV = 83.36(W 1/2 /R) 1.620  Δt=20 ms: PPV = 80.19(W 1/2 /R) 1.608  Δt=50 ms: PPV = 79.66(W 1/2 /R) 1.641   (b)   17(a): (1) With the increase of K r value, the PPVs decrease greatly, indicating that the decoupling charge can reduce the vibration effect of the explosion.(2) As the propagation distance increases, the influence of K r on PPV gradually decreases, indicating that the vibration reduction effect of the decoupling charge will be affected by the distance.at is, the vibration reduction effect is large near and small at a distance.e regression curves and optimum equations based on the USBM equation of horizontal radial blasting vibrations under different decoupling coefficients are shown in Figure 17(b).e site constants K gradually decreases with the increase of K r , and the rate of decrease gradually slows down.e value of α is still stable.

Comparison of Vibration Attenuation Law with Different
Blast-Holes. Figure 18 is a comparison of horizontal radial PPVs at the 5 m measurement point with one hole, two holes (two holes 1 refers to delay blasting, two holes 2 refers to decoupling blasting), and four holes.e general rule is the blasting cases with more blast-holes produced larger PPVs than did the cases with fewer blast-holes.When the decoupling coefficient of two-hole blasting is large, the PPV will drop to less than that of single-hole blasting.
Based on the site constants K single and α single of singlehole blasting, the changes of K/K single and α/α single with the number of blast-holes under delayed blasting, decoupling charge blasting, and different hole arrangements are shown in Figure 19.e variation range of K is within 0.450-1.507.When W 1/2 /R are the same (the single-hole blasting and the delayed blasting when Δt ≠ 0), the site constant K increases as the number of blast-holes increases.When simultaneous blasting (the single-hole blasting, the delayed blasting when Δt � 0, decoupling blasting, rectangular and triangular hole blasting), as the number of blast-holes increases, W 1/2 /R increases, and the site constant K generally shows a decreasing trend.e reason is that as the number of blastholes increases, the increase in charge quantity is greater than the increase in PPV. e site constants K for the cases of 20 ms and 50 ms are almost the same as those of the single-hole case since the vibration waveforms from the 20 ms and 50 ms blasting can be separated within the interval time.e variation range of α is within 0.932-1.118,and the number of blast-holes has a much smaller influence on α than K.   Shock and Vibration

Conclusions
Based on the ANSYS/LS-DYNA software, single-hole and multihole blasting simulation models have been established.For the damage in the near zone and the vibration effect in the middle and far zone of multihole blasting, the influence of the hole arrangement method, the delay time, and the decoupling charge are studied.Based on the USBM equation, the PPV data are fitted to obtain the blasting vibration attenuation parameters.A calculation equation of PPV correction coefficient with delay time as the independent variable is proposed, and the rationality of the equation is verified by field data.e blasting vibration parameters and PPVs of different numbers of blast-holes are compared.e following conclusions are obtained.
(1) In single-hole blasting, in the cross section, the rock damaged area is about 35 times that of the blast-hole, and in the longitudinal section, it is about 8 times.
e blasting vibration parameters of the single-hole blasting are K � 77.91, α � 1.6.
(2) In multihole blasting, compared with rectangular holes, the size of rock mass fragments under the triangular hole arrangement is more uniform and smaller, which is beneficial to engineering production.e horizontal radial PPVs of rectangular holes are larger than that of triangular holes, but the horizontal tangential PPVs are smaller than that of triangular holes.erefore, in blasting practice, the blast-hole arrangement should be selected according to the vibration-sensitive direction of the protected items.e site constants K and α of the rectangularhole blasting are larger than those of the triangularhole blasting.(3) In the two-hole millisecond blasting, the damaged area of the first blast-hole is not affected by the delay time, and that of the second blast-hole increases with the increase of delay time, but the speed of increase  Shock and Vibration gradually slows down, increasing by about 40% at most.Compared with simultaneous blasting, delayed blasting reduces the charge quantity, thereby effectively reducing PPV, especially in nearby zones.e PPV correction coefficient is inversely proportional to the delay time.When Δt ≠ 0, the K value tends to decrease as Δt increases, and the constant α is in a stable range.(4) With the increase of the decoupling coefficient K r , the rock damage area gradually decreases, and the horizontal radial and vertical PPVs decrease greatly, indicating that the decoupling charge can effectively reduce the vibration of the explosion.e site constants K gradually decreases with the increase of K r , and the value of α is still stable.
(5) e general rule is the blasting cases with more blastholes produced larger PPVs than did the cases with fewer blast-holes.When W 1/2 /R are the same, the site constant K increases as the number of blast-holes increases.When simultaneous blasting, the site constant K generally shows a decreasing trend as the number of blast-holes increases.e number of blast-holes has a much smaller influence on α than K.
It should be mentioned that the simulation in this paper is based on blasting with idealized conditions, and there is a certain gap between the complex blasting conditions in engineering practice.us, readers should be cautious of any extrapolation from the results above.Further investigations are needed to study the rock damage and blasting vibration under more complex blasting conditions through theoretical analysis, numerical methods, and field tests in the future.

λFigureFigure 11 :Figure 12 :
Figure Analysis models with different delay times.

Table 3 :
Related parameters of the air material.Density (kg/m 3 ) C 0 C 1 C 2 C 3 C 4 C 5 C 6 E a0 (GPa) V a0 � T/f c is the normalized maximum Braziliantensile strength (T is the realistic maximum Brazilian tensile strength), D 1 and D 2 are material constants.e damage strength DS is defined as follows: * + T *  * As shown in Figure 8, explosive No. 1 detonates first, and explosive No. 2 detonates after an interval of Δt.
6 Shock and Vibration 4.2.Delay Time.Millisecond blasting is a controlled blasting technique often used in precise blastings such as rock foundation excavation, mineral mining, and demolition of structures.In millisecond blasting, multiple blast-holes in a group of charges are detonated in sequence according to a set time interval.It is an effective method to improve rock fragmentation and reduce blast-induced vibration.e analysis is simplified to a group of two blast-holes.

Table 4 :
Related parameters of the rock material JHC.

Table 6 :
Comparison of equation calculation results and field test data results.