Study on the Influence of the Blasting of Preceding Tunnel on the Pre-Penetration Surrounding Rock in Succeeding Tunnel with Small Clear Distance

To study the inuence of the preceding tunnel blasting on the surrounding rock in the pre-penetration area of the succeeding tunnel, combined with the actual monitoring of the Kaifeng Mountain small clear distance tunnel, the relationship between the vibration velocity and the blasting center distance was derived according to the energy attenuation law, predicted the peak velocity of particle vibration on the pre-penetration surrounding rock, came up with the conception of surrounding rock seismic wave vibration velocity ratio, obtained its change situation, and then established a numerical model to further analyze the changes in vibration velocity, plastic strain, and stress. e results show the prediction formula of peak vibration velocity obtained from theoretical analysis of vibration velocity, which can well calculate the peak vibration velocity of some unmeasurable points in tunnel engineering and provide theoretical support for actual tunnel blasting engineering. Under the same clear distance conditions, the front blasting side of surrounding rock in the pre-penetration area of the succeeding tunnel is aected by blasting vibration more than the back blasting side, and the vibration velocity decreases gradually from the arch waist⟶arch shoulder⟶arch foot⟶vault (arch bottom), and the most unfavorable section is that the blasting surface is parallel to the surrounding rock section in the pre-penetration area of the preceding tunnel.e excavation of the lower step has greater destructive force than that of the upper step; with the clear distance that gradually increases, the impact of the preceding tunnel blasting on the surrounding rock in the pre-penetration area of the succeeding tunnel becomes smaller, and D (D is the tunnel span) is the minimum safe clear distance.


Introduction
At present, a small clear distance tunnel is widely used because of the advantages of being convenient for overall road planning, but this type of tunnel usually has a small thickness of the middle rock wall between the two holes, so the construction is di cult. If the vibration velocity generated is too high during blasting excavation of the preceding tunnel, which will cause a very serious impact on the middle rock wall and the adjacent unexcavated succeeding tunnel surrounding rock, change the mechanical properties of surrounding rock and even change the geological structure of surrounding rock, and may easily lead to instability or destruction of the surrounding rock, therefore, it is particularly important to study the vibration velocity prediction in the blasting process of small clear distance tunnel and the in uence of succeeding tunnel in the excavation process. In recent years, many scholars have done a lot of research on the energy attenuation and blasting vibration theory in the process of tunnel blasting, as well as the thickness and the range of plastic area of middle rock wall, the stress change in the surrounding rock, and the construction method in the small clear distance tunnel, and have made some achievements. e blasting of the tunnel is bound to cause the generation, transformation, transmission, and work of energy, and this process is extremely short and may only take tens of microseconds. In addition, due to the uneven internal structure of the rock, the process of energy transmission becomes very complex. Tian et al. [1] found that the blasting vibration energy mainly comes from the explosion in the initial blasting. In the research, many experts and scholars have summarized and deduced particle vibration velocity attenuation equation [2] and blasting vibration energy attenuation equation [3] through engineering measurement and theoretical analysis. On these foundations, Yang et al. [4] identified and separated the waveform of transient unloading induced vibration from the measured vibration signals by means of amplitude spectrum analysis and lowpass filtering and summarized the corresponding attenuation law; Fei et al. and Yao et al. [5,6] studied the energy distribution in different frequency bands, analyzed the attenuation law of energy in seismic wave propagation, and summarized relevant formulas, to describe the variation law of seismic wave energy more quantitatively; Shan et al. [7] used MATLAB software to decompose the vibration signals obtained by blasting of small clear distance tunnel with fivelayer wavelet packet transform, solved the reconstructed signal and frequency band energy of each node, and studied the frequency band energy distribution characteristics of the blasting vibration signal. In addition, by fitting regression analysis of total energy, the corresponding total energy attenuation equation is put forward; Dai et al. [8] calculated the energy contained in the incident wave, reflected wave, and transmitted wave, respectively, based on the SHPB test, and summarized and studied the energy dissipation equation of the corresponding experiment; Liu et al. [9] explored the dynamic characteristics and failure modes of concrete under confining pressure through SHPB dynamic compression experiment and linearly fitted the dynamic growth coefficient (DIF) of concrete with the dynamic confining pressure growth coefficient (DIFc); Ji et al. [10] performed regression analysis on blasting test data based on Sadovsky's equation and studied the blasting vibration velocity propagation law of the surrounding rock in the middle rock in front of and behind the tunnel; Fan et al. [11] aimed at circular tunnel excavation under the condition of in situ stress transient unloading, studied the balance mechanism among surrounding rock kinetic energy, strain energy, and radial stress work, adopted the dimensional analysis method to establish the corresponding vibration velocity attenuation equation induced by the transient unloading of rock strain energy at excavation site, and analyzed the attenuation law that can induce vibration. Many of the above scholars have conducted in-depth studies on the energy form, transmission, attenuation, distribution characteristics, vibration changes, and other mechanisms and laws in the process of tunnel blasting construction and achieved fruitful research results, which provide a solid theoretical basis for exploring the influence of tunnel blasting on surrounding rock.
Many experts and scholars at home and abroad have predicted blasting vibration by various advanced technologies and algorithms and established various prediction models. Advanced even combines neural network with it. Xu et al. [12] studied the blasting flow produced by DLSOS (sublevel open stope method) in the process of ore body mining, analyzed the relationship between ore drawing point spacing and production blasting ring load, and established a back propagation neural network (BPNN) prediction model; Xu et al. [13] established PPV prediction model by combining genetic algorithm with artificial neural network and using mine underground blasting parameters and surface vibration monitoring data; Jayawardana et al. [14] analyzed the relationship between damping effect and control parameters of double-filled grooves by combining neural network and put forward the related vibration prediction model; Guo et al. [15] used genetic algorithm to build a GA-BP neural network prediction model, input the maximum single-stage explosive charge, blasting center distance, and monitoring parameters of the actual blasting project, and realized the prediction of blasting vibration speed; Wang [16] adopted the improved BP neural network algorithm and artificial forest algorithm to analyze the collected actual engineering blasting monitoring parameters and predict the blasting vibration speed and main frequency; in addition, Marilena et al. [17] analyzed 12 experimental cases and summarized the empirical correlation equations of vibration prediction, which provides an excellent basis for evaluating the attenuation of rock mass caused by rock blasting vibration; He et al. [18] carried out dimensional analysis of blasting vibration wave based on field experimental data, decomposed the blasting vibration signal by wavelet packet, researched the peak velocity equation of blasting vibration under the influence of elevation and the attenuation mechanism of energy along slope, and put forward a prediction model of blasting vibration related to its actual engineering; through theoretical analysis and numerical simulation, Jayasinghe et al. [19,20] obtained the vibration velocity waveform of the surface particles caused by the cutting holes blasting and the ground vibration attenuation equation caused by the rock blasting, respectively; Xie et al. [21] determined the pulse function sequence of each blasting hole initiation based on the Anderson linear superposition model, used the convolution theory to calculate the superposition waveform of group hole blasting, and proposed a calculation model of vibration velocity in the vicinity of the tunnel, which can more accurately predict the intensity of blasting vibration; the above experts have done a lot of research on the prediction models of blasting vibration and made some achievements, but the theories and methods based on them are different, and different prediction models are suitable for specific practical projects, but they lack some universal applicability, so the prediction methods of blasting vibration still need further research.
ere are also many scholars who have made in-depth studies on the influence of the construction process of small clear distance tunnel on the lining and surrounding rock. Wu and Jin [22] analyzed the peak vibration velocity under different conditions and determined the maximum effect position by blasting of the preceding tunnel, the reasonable clear distance of the tunnel, and the corresponding control blasting measures; Liu [23] analyzed the surrounding rock stress, displacement, plastic area distribution, ground settlement, and other parameters after tunnel excavation under the conditions of different middle rock thicknesses and obtained the reasonable clear distance, reinforcement method, and effect of middle rock wall, through numerical model calculation; Yang et al. and Guo et al. [24,25] studied the relationship between tunnel excavation and clear distance and found that the stress of rock wall was greatly affected by the clear distance, damage or crack propagation would occur on the front blasting side and back blasting side, the damage degree gradually decreased with the increase in clear distance, and the corresponding minimum safe clear distance is determined; Li et al. [26] simulated the stress characteristics and deformation of small clear distance tunnel structure under different excavation methods and discussed the corresponding reinforcement measures of the middle rock wall; Wang et al. [27] studied the influence law of the excavation process of under passing tunnel on the settlement deformation, vibration acceleration, and vibration velocity of the existing adjacent tunnel lining structure under the conditions of different construction methods, different sandwich thickness, and different train axle load; Guan et al. [28] established three-dimensional numerical models of square, circular, and horseshoe tunnels using the ALE algorithm of LS-DYNA software and compared and analyzed the influence of different types of tunnel distance changes, and it was found that the peak vibration velocity and tensile stress of three types of tunnels increased with the decrease in distance, among which the circular tunnel was the most affected, followed by the square tunnel and horseshoe tunnel; Yang et al. [29] found the weakest position of the preceding tunnel when the succeeding tunnel is excavating and obtained the relationship between the blasting vibration velocity of each particle and time; Jia et al. [30] studied the influence of the blasting excavation of a new tunnel on the vibration of existing adjacent tunnel lining by establishing a numerical model and analyzed the influence of different clear distances, excavation footage, and buried depth; Liu [31] studied the effect of the supporting structure of the preceding tunnel on the succeeding tunnel and the influence of the excavation of the succeeding tunnel on the surface settlement, vertical displacement, lining structure, and surrounding rock stress of the preceding tunnel and put forward the corresponding monitoring and reinforcement measures; the above domestic and foreign scholars have done a lot of research on many aspects, such as rock wall thickness, surrounding rock stress, plastic area distribution, and construction methods during the construction of small clear distance tunnel, and these studies have played an important role in realizing the smooth construction of small clear distance tunnel. e blasting of the preceding tunnel is bound to cause the stress redistribution and strain change in the middle rock wall and the pre-penetration surrounding rock of the succeeding tunnel. However, most of the current researches focus on the analysis of the stress and strain of the lining structure and the middle rock wall of the succeeding tunnel. ere is a lack of research on the change law of the vibration velocity, stress, and strain of the surrounding rock of the preceding tunnel, and it should be further explored.
is article is based on the actual project of the Kaifeng Mountain small clear distance tunnel. Firstly, the relationship between the particle vibration velocity and the burst center distance is derived through the seismic wave attenuation equation and Sadovsky's empirical equation, and then, the peak velocity vibration, frequency, and other data from field monitoring and simulation result are analyzed and fitted to obtain the corresponding coefficients k, α, and β. Finally, the relationship between the vibration velocity and the blasting center distance of tunnel project is used to calculate the seismic wave velocity proportional coefficient V r /V 0 of the surrounding rock in the pre-penetration area, the variation trend of the vibration velocity is judged, and the most unfavorable section in the process of blasting excavation is obtained.
en, the numerical model is established by the finite element software, and further, the changes in rock wall and pre-penetration surrounding rock are analyzed in the succeeding tunnel caused by blasting construction of the preceding tunnel through the vibration velocity, plastic strain, and stress changes. Finally, combined with field monitoring and theoretical derivation, the relevant laws are summarized to provide a reference for the construction design and actual construction of the small clear distance tunnel.

Project Overview and Vibration
Monitoring Scheme e Kaifeng Mountain tunnel is located in Baoji Village, Hongshan Town, Suizhou City, Hubei Province, which is a small clear distance double-hole separated tunnel in the Mazhu Expressway. e biggest difference between the two-hole separated tunnel with small distance and the single-hole tunnel is that the two holes will influence each other in the construction process, resulting in the deterioration of surrounding rock conditions. e blasting construction of the preceding tunnel will increase the vibration influence on the surrounding rock of the succeeding tunnel, which will disturb the surrounding prepenetration rock, worsen the stress conditions of the tunnel, and even lead to the destruction and instability of the surrounding rock. e total length of the tunnel is about 1000 m, the span of the excavation section is 13.12 m, and the thickness of the middle rock wall is 9∼12 m. e tunnel is spread from southeast to northwest direction. It is a doubletrack tunnel with a maximum buried depth of 70 m, located in hilly area, and the surrounding rock is grades IV and V, the geological conditions are complex, and the self-stability is poor. e detailed clearance section design is shown in Figure 1 (all units not shown in the figure are in cm).
e Kaifeng Mountain tunnel is constructed by upper and lower step methods and excavated by drilling and blasting method, and the smooth blasting technology is adopted, and the working faces of the two tunnels are staggered. e explosives are No. 2 rock ammonium nitrate explosive and emulsion explosive, and the blasting method is segmented delayed detonation network and the detonating electric detonator initiation. e total charge of the tunnel face blasting is 78 kg, divided into 8 delays, and the Shock and Vibration maximum charge of a single stage is 15.6 kg. e diameter of cartridge is 0.032 m, the length of cartridge is 0.2 m, and the blasting hole is densely blocked with gun clay, which is composed of three components: clay, medium coarse sand, and water. e stemming length should not be less than 0.4 m, and the charge coe cient should be controlled between 75% and 80%, and intermediate holes (1 # , 3 # , 5 # , 7 # ), cutting holes (9 # , 11 # ), and peripheral holes (13 # , 15 # ) are blasted successively (Figure 2). e actual blasting e ect is excellent, the residual blast hole traces are evenly distributed, the integrity of rock mass is not damaged, there is no obvious blasting crack on the excavation surface, and the blasting indexes meet the design requirements.

Monitoring Program.
When the preceding tunnel is constructed, the original mechanical system of surrounding rock around the tunnel is broken, and the original stress state has changed. With the increase in excavation depth, this change will continue. To monitor the state and mechanical dynamic behavior of the pre-penetrating surrounding rock of the succeeding tunnel, the surrounding rock of the preceding tunnel and the pre-penetrating surrounding rock of the succeeding tunnel can be monitored by means of monitoring and measurement. e results can guide the construction, foresee the dangerous situation of accidents, take timely measures to prevent problems before they happen, and provide an analogy basis for the design and construction of similar projects in the future. e main instruments used in this blasting monitoring system include horizontal and vertical velocity sensors, UBOX-20016 vibration recorder, and computer. e horizontal and vertical sensors are used to collect dynamic analog signals such as vibration velocity and vibration acceleration generated by tunnel blasting. UBOX-20016 portable data acquisition equipment is used for digital conversion and storage of the required signals, and the e ective data obtained will be processed and analyzed by the computer. e speci c working principle is supplemented and shown in Figure 3.
In the process of construction, it is necessary to select representative sections for the blasting vibration test. In the preceding tunnel, several vibration monitoring sensors are arranged at the arch foot and the arch waist on each cross section ( Figure 4) and at di erent distances behind the tunnel face on the vertical section ( Figure 5). e in uence of blasting vibration wave on surrounding rock at di erent distances can be measured. e site layout is shown in Figure 6.
According to geological conditions and construction scheme of Kaifeng Mountain tunnel, combined with the corresponding measurement point layout principles and monitoring requirements, the corresponding measuring points are reasonably arranged in the direction of vertical section. Four monitoring points (1 # , 2 # , 3 # , and 4 # ) are arranged on the wall of the front blasting side of the preceding tunnel, which are 10 m, 30 m, 65 m, and 85 m away from the blasting working face, respectively, and the other two monitoring points (5 # and 6 # ) are, respectively, arranged at 120 m from the explosion source on the back blasting side of the succeeding tunnel and 125 m from the explosion source on the front blasting side. In this study, the in uence of the preceding tunnel blasting on the surrounding rock in the pre-penetration area of the succeeding tunnel is studied, since it is impossible to arrange the monitoring points inside the surrounding rock in the pre-penetration area of succeeding tunnel. erefore, the derivation points are, respectively, arranged on the face and back blasting sides of surrounding rock in the pre-penetration area of the succeeding tunnel: taking the blasting working face as the starting point, the lines of the derivation points extend 90 m and 50 m backward and forward, respectively, and a deduction point is arranged each 10 m. e layout is shown in Figure 5.    Hz. ere is no overlap between the two frequency ranges, and the main vibration frequency is much higher than the natural vibration frequency of the tunnel. So, there will be no resonance amplification in this blasting, and the frequency will not be considered, and only the vibration velocity will be analyzed. During the blasting excavation of this tunnel project, the maximum vibration velocity is 3.564 × 10 −2 m/s, which is far less than the vibration velocity value of 12∼15 cm/s specified in the national standard blasting safety regulations (GB6722-2014) ( Table 2), and meets the blasting safety requirements.

Numerical Simulation Calculation of Blasting Excavation of Small Clear Distance Tunnel
is article takes the Kaifeng Mountain tunnel project as the background, establishes a three-dimensional numerical model using finite element software, and defines the left tunnel as a succeeding tunnel without excavation and the right tunnel as a preceding tunnel under excavation. e influence of blasting excavation of upper and lower steps of the preceding tunnel in the pre-penetration area of the succeeding tunnel is simulated under the conditions of IV type surrounding rock and different tunnel clear distances (0.15D, 0.5D, D, and 2D, while D is the tunnel span 13.12 m), and the changes in vibration velocity, plastic strain, and stress distribution of the surrounding rock in the prepenetration area are analyzed.

Basic Dimensions and Parameters of the Model.
According to existing studies [32]: the influence range of blasting excavation of tunnel is (3-5) D. erefore, to avoid the interference caused by boundary conditions, the model takes 3D as the boundary. is article takes the tunnel clear distance of 0.15D/0.5D/D/2D to model. e span of the tunnel is 13.12 m, the height is 11 m, and the rock layers covered above the tunnel are weathered soil, strongly weathered sandstone, and carbonaceous siliceous shale with thicknesses of 9.6 m, 12.8 m, and 77.6 m, respectively. According to the influence range of rock excavation, the horizontal dimension is 120 m, the vertical height is 100 m, and the depth dimension is 60 m by calculating. e tunnel numerical calculation model is divided into 98573 nodes and 108208 units, as shown in Figure 7. e parameters of tunnel surrounding rock and supporting materials are shown in Table 3.

Determination of Model Boundary Conditions.
e shock wave generated by blasting will be reflected on the surrounding rock boundary, which will cause a big gap between the simulated calculation data and the on-site construction. erefore, it is necessary to set the boundary conditions on the surrounding rock to simulate the continuous state of the rock in the actual project, so that the software calculation results can be more consistent with the actual situation. In the eigenvalue analysis of the model, considering the absorption of scattered wave energy by surrounding rock and its viscoelastic properties [33], the model is set as the viscoelastic boundary condition, to restore the real environmental state of rock as much as possible. Firstly, the boundary conditions of the support are established, and the reaction coefficient of elastic foundation is calculated as follows [34]. e coefficient of vertical foundation reaction is as follows: e coefficient of horizontal foundation reaction is as follows: In the above equations: Secondly, the reflection of wave should be considered in the dynamic analysis, and the damping constant should be manually input on each boundary of the model to generate a damper, and then, the absorption of the scattered stress wave is simulated. e viscous boundary proposed by Lysmer and Kulemeyer [35,36]is used. e equations for calculating the damping ratio in each direction of the model required by the viscous boundary condition are as follows.
P wave is as follows: S wave is as follows: Monitoring of measuring point at the arch foot Figure 6: Field layout of instruments and sensors. 6 Shock and Vibration In the above equations: e above equations can be used to calculate the foundation damping coefficient of each rock layer in the model, as shown in Table 4.

Determination of Blasting Load Parameters.
Before the numerical calculation of tunnel blasting, it is necessary to determine the relevant parameters of blasting loading, including peak load, blasting loading waveform, and loading and unloading time. e calculation of each parameter plays an important role in the accuracy of the model and the degree of conformity with the actual project.
After blasting hole explosive, the formation of blasting stress wave load is evenly distributed on the tunnel lining, and the action direction is perpendicular to the tunnel wall [37]; according to the propagation law of wave, the form of blasting load can be simplified to the curve shown in Figure 8 [38]. P max in the figure is the peak pressure on the blasting hole wall, which is related to rock properties, explosive types, charge forms, and other factors. e explosion stress wave at any point in the rock behaves as load similar to triangle, and the peak pressure decreases sharply after reaching the peak load pressure and keeps spreading forward.
Due to the functional characteristics of the software, it is difficult to establish blasting hole units one by one. erefore, after the blasting holes are gathered together, the pressure is assumed to act vertically on the imaginary surface of the blasting, and then the normal uniform load P max applied to the surrounding rock at the tunnel wall in the numerical simulation is calculated according to the empirical method. at is, in the equation mentioned in the National Highway Institute of the United States, the blasting load per 1 kg is as follows [39][40][41]: In the above equations: � blasting pressure acting on the hole wall (Pa) P de t � explosion pressure per 1 kg explosive (Pa)

Shock and Vibration
sge � specific gravity of explosives (kg/m 3 ) V � blasting speed of explosive (m/s) e loading and unloading time of blasting load can be determined by the following equation [30].
Loading rise time is as follows: Total action time is as follows: In the above equations: K � bulk modulus of compression (Pa); K � E/3(1 − 2μ) Q � blasting hole charge (kg) r � contrast distance (the ratio of the distance r from a point to the blasting center to the radius of the blasting hole r b ) t R � rise time period (s) t s � total action time of triangle waveform (s) μ � Poisson's ratio of rock According to the charging structure, number of blasting sections, charging amount, and data results obtained from blasting monitoring of Kaifeng Mountain tunnel, the total action time of the blasting load applied by the numerical model is determined. It can be seen that the pressure boost time of triangular wave load is 0.008∼0.012 s and the pressure relief time is 0.04∼0.12 s [42]. To fit the actual construction situation as much as possible, after calculation, the numerical simulation adopts the pressure boost time of 0.008 s, the pressure relief time of 0.092 s, and the total time of 0.1 s. e explosive velocity is 3000 m/s, the specific gravity is 1.1 × 10 3 kg/m 3 , adopt uncoupled charge, the cartridge diameter is 0.032 m, and the drilling diameter is 0.04 m. According to (5), the explosion pressure for each 1 kg can be calculated as P B � 1.127 × 10 7 Pa. Each excavation footage of upper step is 2 m, and the maximum charge is 13.2 kg. e blasting pressure can be calculated as 1.488 × 10 8 Pa through (6); each excavation footage of lower step is 2 m, the maximum charge is 26.6 kg, and the calculated blasting pressure is 2.93 × 10 7 Pa.

Load Results
Applied by the Model. According to the parameters of rock and soil in the project, the corresponding indexes in the modeling process are set, the measuring points 1 # ∼6 # at the same position of the model ( Figure 5) are arranged, the vibration velocity time-history curve (Figure 9(a)) of each measuring point is obtained, the peak vibration velocity of each measuring point is extracted, and the numerical simulation vibration data are sorted. e results are shown in Table 5. Compared with the peak value of vibration velocity obtained from the actual engineering monitoring (Figure 9(b)), so as to verify the accuracy of the model, the maximum difference between them is 3.94 × 10 −3 m/s at No. 2 measuring point, and the relative error is 14.7%. e minimum value is 7.8 × 10 −4 m/s at No. 3 measuring point, and the relative error is 4.21%. e average value is 8.23 × 10 −4 m/s, and the average relative error is 10.4%. It can be seen that the error between the peak value of numerical simulation vibration velocity and the peak value of actual monitoring vibration velocity is slight, which indicates that the parameters used in the finite element model are properly selected, and the model can better simulate the reality of the construction site.

Prediction of Pre-Penetrating Surrounding Rock Vibration Velocity in Succeeding Tunnel
When the explosive explodes, the energy is released in the form of shock wave with high temperature (3000°C) and high pressure (10∼100 Gpa) within about 3 ∼ 7 times the radius of the charge from the explosion source firstly [43], as the blasting center distance increases to 120∼150 times the charge radius, and the shock wave transforms into a compressive stress wave. When it exceeds 150 times the charge radius, the compressive stress wave will decay into a periodic vibration seismic wave [44], and the attenuation rate and action range of these three waveforms are different. If the shock wave and stress wave have a large intensity, which leads to plastic damage or even breakage of surrounding rock, it will not cause damage to surrounding rock until it is transformed into seismic wave, but only cause elastic vibration [45]. e surrounding rock in the pre-penetration area of the succeeding tunnel is far away from the blasting center. erefore, the surrounding rock of the tunnel is mainly affected by seismic wave during the blasting process. e vibration velocity of the surrounding rock is proportional to the stress produced by blasting, so the vibration velocity can directly reflect the energy input from the seismic wave to the structure, and it provides a way to measure the stress change in the surrounding rock. e propagation of seismic waves in rock media is often affected by multiple factors, including the damping effect generated by the internal structure of the rock and the scattering reaction caused by its inhomogeneity, and all the above can lead to the consumption of blasting seismic wave and energy; at the same time, the vibration velocity will continue to decay.

4.1.
eoretical Analysis of Vibration Velocity of Tunnel Surrounding Rock. During the tunnel construction, the explosive will disturb the surrounding rocks and soil media at the moment of explosion, and there will be waves radiating in all directions from the blasting point. e blasting seismic waves can be approximately regarded as plane waves at a distance from the disturbance center. e simplest explanation for the problem of seismic wave absorption is assuming that the relative decrease in wave energy is   proportional to the distance traveled by the wave. Under the influence of the internal structure of surrounding rock and other factors, the energy in the process of seismic wave propagation decreases continuously, and the energy consumed is proportional to the propagation distance [3,45,46], namely, In the equation: E � the energy when the seismic wave propagation distance is r (m 2 /s −2 ) E 0 � the initial energy (energy at the source when the propagation distance of seismic wave tends to zero) r � seismic wave propagation distance (m) β � seismic wave absorption coefficient e energy generated by tunnel blasting is proportional to the square of the vibration velocity of the surrounding rock.
When r tends to 0, the seismic wave energy E � cv 0 2 , where v 0 is the initial vibration velocity when the seismic wave is generated (which can be approximately regarded as the vibration velocity of the rock mass at the explosion source), and it is substituted into formula (13), namely, e formula of peak vibration velocity at a certain blasting center distance can be obtained by simultaneously calculating square root on both sides of formula (14): In equation (15), k, α, and β can be derived from the monitoring data of the tunnel blasting project, and the formula can predict the vibration velocity of surrounding rock particles where the actual monitoring cannot be carried out and can be applied to evaluate the blasting effect under known blasting conditions. Fitting calculation of k, α, and β coefficients is as follows.
(1) e coefficients k and α are fitted and calculated by Sadovsky's empirical equation: In the equation: k � coefficient related to propagation medium and blasting conditions Q � the maximum charge of a single stage in delayed blasting (kg) R � distance from blasting center (m) V � blasting particle vibration velocity (m/s) α � vibration attenuation coefficient, related to geological and topographic conditions and the distance from the blasting center (2) e absorption coefficient β is related to the amplitude [49,50]: In the equation: A � amplitude of blasting seismic wave (m) Formula (15) is generally applicable to all kinds of tunnel blasting excavation projects, but the coefficients k, α, and β are related to the actual geological and topographic conditions, explosive charge, blasting center distance, amplitude, and other factors of a certain tunnel project. By fitting and calculating the coefficients k, α, and β according to the actual project monitoring data, the prediction formula of characteristic vibration velocity applicable to this tunnel can be obtained. When applied to this tunnel project, the peak vibration velocity v 0 of the pre-penetration surrounding rock of the succeeding tunnel can be calculated, and then, the influence of the blasting construction of the preceding tunnel on the succeeding tunnel can be judged.

Determination of Peak Vibration Velocity of Tunnel Pre-Penetrating Surrounding Rock.
According to the actual blasting monitoring situation and numerical simulation results of Kaifeng Mountain tunnel (Tables 4 and 5), the corresponding unknowns k, α, and β in the vibration velocity formula (15) are calculated using the above derivation process, and then, the formula is applied and verified.
Combined with Sadovsky's empirical equation (16) in Section 4.1, the data of total charge, blasting center distance, and particle vibration velocity obtained from the actual blasting and simulation result are analyzed. en, the reduced distance Q 1/3 /R in (16) (the maximum charge Q of a single section of this project is 15.6 kg) is taken as the x-axis, the peak vibration velocity v of the monitoring points is taken as the y-axis to establish a coordinate system, and the relationship curve is drawn between the reduced distance and the vibration velocity, as shown in Figure 10(in the Figure: V � peak vibration velocity; Q 1/3 /R � reduced distance).
It can be seen from Figure 10 that the relationship between the reduced distance and the vibration velocity obtained from the actual engineering monitoring is y � 38.339x 0.392 , and the correlation coefficient R 2 � 0.944; the relationship between numerical simulation reduced distance and vibration velocity is y � 39.414x 0.296 , and the correlation coefficient R 2 � 0.942; the difference between the curve and formula fitted by them is very slight, which indicates that the fitting effect is well. erefore, the fitting result obtained from the actual engineering monitoring with relatively high correlation coefficient will be more suitable for this tunnel engineering, with the coefficient k of 38.339 and the vibration attenuation coefficient α of 0.392.
In (15), the seismic wave absorption coefficient β still needs to be solved. (17) is combined to analyze the amplitude and the detonation center distance, and the relevant curve is drawn, as shown in Figure 11 (in the figure: A � amplitude; R � blasting center distance). e x-axis is the distance between the monitoring point and the explosion source (the blasting center distance r). e y-axis is −Ln(A). e fitting formula obtained from actual engineering monitoring is y � 0.0071x + 8.451, and the correlation coefficient R 2 � 0.911; the fitting formula obtained by numerical simulation is y � 0.00691x + 8.454, and the correlation coefficient R 2 � 0.900; and the slope of the fitting curve is the seismic wave absorption coefficient β, and the seismic wave absorption coefficient β is 0.0071 when the fitting value of engineering monitoring is taken.
Combined with the coefficients k, α, and β obtained by fitting the above monitoring data and simulating result, the particle vibration velocity v r in the seismic wave propagation process of the tunnel engineering can be obtained: e theoretical vibration velocity of 1 # ∼6 # measuring points can be calculated by theoretical derivation formula (18) and compared with the actual vibration velocity obtained by monitoring and simulating ( Figure 12) (in the figure: v � peak vibration velocity; r � blasting center distance). It can be seen from the figure that the deviation among the theoretical value, the actual monitoring value, and the numerical simulation value is very slight, and all show a decreasing trend with the increase in blasting center distance, showing the same law basically, so the above theoretical derivation formula can play a certain role in predicting some special unmeasurable points on the surrounding rock. e initial velocity v 0 is the same during blasting, which is difficult to monitor. To study the vibration velocity of each particle in the surrounding rock of the succeeding tunnel, the vibration velocity v r of each particle in the process of seismic wave propagation is divided by the initial velocity v 0 to obtain the ratio χ, which is defined as the seismic wave vibration velocity ratio, namely, According to (19), the seismic wave vibration velocity ratio χ of derivation point can be calculated to reflect the vibration velocity change in surrounding rock in the prepenetration area of the succeeding tunnel, and then, the influence of blasting vibration of the preceding tunnel on the vibration at any position of the succeeding tunnel can be obtained. Combined with the calculated seismic wave velocity ratio χ of each derivation point, a coordinate system with the blasting working face as y-axis and the blasting excavation direction of the tunnel as x-axis ( Figure 5) is established. A curve graph is drawn according to the seismic wave velocity ratio χ of each theoretical calculation point changes with its location, as shown in Figure 13.
From the analysis in Figure 13, it can be seen that the closer the location of the derivation point is to the origin coordinate system, the seismic wave velocity ratio χ is greater; otherwise, it becomes smaller and smaller. When the abscess value is 0, it will reach the maximum value; that is, the surrounding rock section in the pre-penetration area corresponding to the blasting face is the most unfavorable section. It shows that the closer to the blasting working face, the vibration velocity of the surrounding rock and the disturbance are greater. In addition, by comparing the vibration velocity ratio χ of the derivation points on the front blasting side and the back blasting side, the χ value of the front blasting side is larger than the back blasting side. As the distance between the derivation points and the original point becomes further and further, the difference between the two derivation points becomes smaller and smaller. It shows that the impact of vibration velocity on the front blasting side of pre-penetration surrounding rock of the succeeding tunnel is greater than that on the back blasting side, and the influence on both sides gradually becomes consistent with the increase in distance.
Due to the limitation of construction conditions, the monitoring points and monitoring times that can be arranged are limited. At the same time, the stress redistribution and the failure pattern of surrounding rock in the pre-penetration area cannot be obtained completely through on-site monitoring. So, a numerical simulation calculation study is further conducted, which can deeply explore the variation of the plastic zone of the middle rock wall during the preceding tunnel blasting and the stress distribution of the surrounding rock in the pre-penetration area of the succeeding tunnel.

Numerical Simulation Calculation and Analysis of Engineering Construction
In the simulation calculation, the same upper and lower step methods are used as the actual project, and the distance between the upper and lower step working face is 10 meters.
Only the most unfavorable section in the conclusions obtained from the analysis of the theory and monitoring results for analysis are taken, and the model has a total of 8 monitoring points arranged along its pre-excavation contour line. e measuring points at the top and bottom of the arch are 1 # and 8 # , the arch shoulders are 2 # and 3 # , the arch waist is 4 # and 5 # , and the arch feet are 6 # and 7 # . e specific layout is shown in Figure 14.

Vibration Velocity Analysis.
e simulated vibration velocity data are drawn into the maximum vibration velocity envelope diagram (Figure 15), and the vibration velocity comparison diagram of surrounding rock in the pre-penetration area of the succeeding tunnel when the upper and lower steps of the preceding tunnel are excavated is drawn ( Figure 16).
rough Figure 15, it is found that the maximum vibration velocity always appears at the 5 # measuring point, which is the arch waist, and the vibration velocity of other measuring points gradually decreases from arch waist⟶ arch shoulder⟶ arch foot⟶ arch top (bottom). As the clear distance gradually increases, the vibration velocity at the same position gradually decreases, and the area formed by the envelope diagram becomes smaller and smaller. It shows that with the increase in the tunnel clear distance, the impact of blasting vibration on the surrounding rock of the pre-penetration area becomes smaller and smaller.
According to Figure 16, the maximum vibration velocity of the measuring points on the front blasting side (1 # , 3 # , 5 # , 7 # , and 8 # ) and the back blasting side (2 # , 4 # , and 6 # ) is compared and analyzed and found that under the same clear distance condition the vibration velocity of the front blasting  side is greater than that of the back blasting side, and the variation range is larger, so the graph line is saw-toothed, and this is consistent with the law of on-site monitoring data.
With the increase in the clear distance, the graph line gradually tends to a straight line, indicating that the surrounding rock can consume part of the energy generated by   blasting, thus reducing the vibration velocity. Comparing the vibration velocity of the same measuring point under different clear distance conditions, it is found that due to the vibration damping effect produced by the surrounding rock, the larger the clear distance, the smaller the maximum vibration velocity. According to the requirements of "GB 6722-2014 Blasting Safety Regulations," the maximum allowable particle vibration velocity of tunnel blasting is 0.2 m/ s, and when the tunnel clear distance is 0.15D, the measuring points at the side facing the blasting during the excavation of the upper step in the succeeding tunnel and the eight measuring points during the excavation of the lower step all exceed the specification requirements. It shows that if the clear distance is too small, the blasting vibration of the preceding tunnel will easily affect the stability of the surrounding rock in the pre-penetration area of the succeeding tunnel, so when the succeeding tunnel is excavated, the corresponding reinforcement measures should be taken in advance. All the measuring points with the tunnel clear distance greater than or equal to 0.5D meet the specification requirements, so the tunnel clear distance of 0.5D can be used as a dividing line of safety clear distance. Under the same clear distance condition, the vibration velocity of the tunnel at the same position when the lower step is excavated is greater than that when the upper step is excavated.

Plastic Strain Analysis.
When the surrounding rock is destroyed, the plastic deformation and displacement on the failure surface will suddenly change. erefore, the maximum point of plastic strain in each section of the surrounding rock is connected into a line, and the potential failure surface of the surrounding rock can be obtained, to judge the stability of the surrounding rock. Figures 17-20 show the distribution of plastic strain in the rock wall and pre-penetration surrounding rock of the succeeding tunnel when the preceding tunnel is excavated under different clear distance conditions, in which the black line segment is the potential plastic failure surface. rough analysis, no matter what the clear distance of the tunnel is, the plastic strain area generated by the excavation of the upper step in the preceding tunnel tends to develop toward the arch top, and the plastic strain area generated by the excavation of the lower step tends to develop toward the arch bottom. When the clear distance is 0.15D, the middle rock wall and the surrounding rock of pre-penetration area of the succeeding tunnel all enter plastic state, and the potential failure surface of plastic strain is particularly obvious when the lower step is excavated, which is U-shaped and has a large distribution range. e results show that when the clear distance is 0.15D, the pre-penetration surrounding rock and the middle rock wall will be disturbed greatly, so the prereinforcement measures should be taken before the blasting of the succeeding tunnel. With the increase in clear distance, the plastic strain area and potential failure surface gradually become smaller and move away from the middle rock wall; that is, the influence of blasting excavation of the preceding tunnel on the middle rock wall and the surrounding rock of the succeeding tunnel is smaller and smaller. When the clear distance is greater than or equal to D, the surrounding rock in the pre-penetration area of the succeeding tunnel always maintains a state of elastic. So, according to the change law of plastic strain, D can be judged as the minimum safe clear distance.

Stress Analysis.
To further study the stress distribution of the surrounding rock in the pre-penetration of the succeeding tunnel when the preceding tunnel is blasting, the stress distribution during the blasting of upper and lower steps of the preceding tunnel under different clear distance conditions is analyzed, as shown in Figures 21-24.
From Figures 21-24, it is found that the stress of surrounding rock is concentrated mostly around the surrounding rock of the preceding tunnel, and the distribution range expands with the increase in the clear distance of the tunnel, gradually moves away from the middle rock wall and the pre-penetration area of the succeeding tunnel, and tends to a safe and stable state. When the upper step is excavated, the maximum stress of the surrounding rock in the pre-penetration area of the succeeding tunnel is mainly distributed around the            rate of rock blasting. erefore, the static tensile strength of rock mass is used to approximately replace the uniaxial dynamic tensile strength of rock mass in simulation), and the surrounding rock will not cause the type of tensile failure in this condition. When the lower step is excavated, the maximum stress of the pre-penetration surrounding rock is mainly distributed around the arch bottom of the lower step, and the stress distribution range is relatively large. e stress value of the whole middle rock wall and most of the pre-penetration surrounding rock of the tunnel with a clear distance of 0.15D is greater than the tensile strength of the rock, and the surrounding rock is prone to tensile failure. While in the tunnel with a clear distance of 0.5D, the half range of the front blasting side of the middle rock wall is the tensile failure area. When the clear distance of the tunnel is greater than or equal to D, the stress values of all middle rock wall and the surrounding rock of the succeeding tunnel are less than the tensile strength of the rock. e surrounding rock of the tunnel is in a stable state and will not be destroyed, so D is the safe clear distance [57].

Conclusions
(1) According to theoretical derivation and field monitoring analysis, it is concluded that the peak vibration velocity is negatively correlated with the blasting center distance of the preceding tunnel; the larger the distance, the smaller the vibration velocity, the smaller the seismic wave velocity ratio, and the smaller the disturbance to the surrounding rock. e surrounding rock on the front blasting side is more affected by blasting vibration than that at the back blasting side of the succeeding tunnel, and the section of the pre-penetration area of the succeeding tunnel parallel to the blasting working face of the preceding tunnel is the most unfavorable. (2) rough the analysis of the vibration velocity of the numerical simulation, it is found that the vibration velocity of the front blasting side at the same position is greater than that of the back blasting side, and the vibration velocity gradually decreases from arched waist⟶ arched shoulder⟶arched foot⟶arch top (bottom), which is consistent with the law obtained from field monitoring. With the gradual increase in the clear distance, the variation range of vibration velocity between each measuring point decreases, and the impact of blasting vibration on the surrounding rock in the pre-penetration area becomes smaller and smaller.    (3) Under the same clear distance conditions, compared with the stress value and the range of plastic area in the rock wall and the pre-penetration surrounding rock of the succeeding tunnel, it is found that the lower step excavation is larger than that of the upper step excavation of the preceding tunnel, and the disturbance and failure of the pre-penetration surrounding rock and the middle rock wall are stronger. With the increase in the clear distance, the disturbance effect of the blasting of the preceding tunnel on the succeeding tunnel gradually decreases until both the middle rock wall and the pre-penetration surrounding rock reach a stable state. (4) By analyzing the change law of the maximum vibration velocity, it is found that all measurement points meet the specification requirements when the clear distance of the tunnel is greater than or equal to 0.5D. It is found that when the clear distance is greater than or equal to D, the excavation of the preceding tunnel will not affect the middle rock wall and the pre-penetration surrounding rock of the succeeding tunnel by analyzing the stress distribution and plastic strain. erefore, based on the research results of three aspects, it is concluded that the minimum safe clear distance is D. (5) From the prediction formula of vibration velocity, it is found that the vibration velocity of surrounding rock in tunnel blasting construction is mainly related to geological and topographic conditions, explosive charge, blasting center distance, amplitude, and other factors. rough the actual engineering monitoring and the fitting calculation and mutual verification of numerical simulation data, it shows that this theory can be well applied to the corresponding tunnel, can predict the maximum vibration velocity of surrounding rock in the pre-penetration area, and can well guide the engineering construction, which can provide a reference for similar projects.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this study.