Application of Explosive Equivalency Approach in Blast-Induced Seismic Effect Prediction Using EXPLO5 Thermochemical Code

: Blasting is a key process that plays a significant role in various industries, including mining and construction. To measure the effectiveness and potential impact of a blast generated by different explosives, industry professionals use a widely accepted parameter known as TNT (trinitrotoluene) equivalent. This manuscript provides an overview of the approach based on the application of the explosive equivalency principle in the prediction of the seismic effects caused by the detonation of different explosives. The explosive equivalents of studied explosives are derived from the results of thermochemical calculations using the EXPLO5 code and compared to field tests. The results have demonstrated that the equivalency approach can potentially be a useful tool in the assessment of blast-induced seismic effects.


Introduction
Explosives have been utilized for centuries to achieve construction projects, clear land, and break rock.However, the use of explosives can have significant negative effects on the environment, i.e., blast-induced seismic effects, fly-rock, and air blasts.Among these effects, seismic effects are considered the ones with the highest importance [1].Thus, it is crucial to control these negative effects to comply with existing regulations and laws.
The impact of blast-induced seismic effects on buildings involves three essential steps: estimation of ground motion, analysis, and the establishment of permitted PPV (peak particle velocity) limits [2].Today, due to technological progress, only approach to achieving results differs.Many scholars have developed equations to predict ground vibration, which are summarized by Kumar et al. [3].Some of them are adapted to specific geological conditions and areas.For example, Ozer proposed different equations for different zones and geological conditions within the "Istanbul, Kadikoy-Kartal Railway Mass Transport System" project [4].ISEE gave the equation for the 95% confidence line for typical data from quarry blasting [5].A similar equation for the 95% confidence line has been given by Ak et al. for measured PPV during blasting works in magnesite mines [6].Nicholson proposed a prediction equation developed from measurements during Bengal Quarry blasting [7].Badal, in his thesis "Blast vibration studies in surface mines", proposed an equation for the "Jindal Power Open Cast Coal Mine" [8].Mesec conducted his research on numerous blasting sites in deposits with different GSI (Geological Strength Index) [9].Moreover, damage characteristics in a rock mass have been extensively studied [10].Agrawal modified the scaled distance equation using the superimposition factor of vibration waves [11].
Most of the published articles about blast-induced seismic effects are based on a significant amount of measured values gathered during a particular project or quarry development.These data are then analyzed using statistical analysis or ANN (Artificial Neural Network) [12,13].Considerable quantities of open-pit blast-induced seismic effects measurements can be collected to give a regression curve.Civil works usually do not have this quantity of measured data; hence, the execution of the initial test blast is fundamental for developing an equation for the safe execution of blasting works [14].Guidelines exist for measuring instruments and installation methods during measurement [15].
To increase the precision of blast-induced seismic effect measurements, we established a research program and have performed several studies according to it.A greater number of monitoring points were used in a single measurement line [16], positioning of measurement instruments at relevant predetermined distances [17], and detecting and excluding atypical values to achieve more accurate input data and, thus, consequently, more precise results [18].
For contractors and owners, ensuring the safety of surrounding buildings during blasting works is of utmost concern [19].In this regard, the power of an explosive is a key factor in determining its effectiveness and safety.To allow for a standard comparison of explosive power, the concept of TNT equivalent has been introduced, which enables professionals to make informed decisions in blasting operations.
The TNT equivalence, which is a common engineering tool, is one of the potential tools for improving the effectiveness of blasting and the precision of blast-induced seismic measurement results that has been widely investigated [20].It is used to relate the effects of the yield of a certain explosive to that of TNT (or another reference explosive).It is established mainly empirically by different types of experiments, which often do not yield the same TNT equivalence values.The most used tests are the Trauzl test, plate dent, ballistic mortar, sand crush, and air blast test.Predicting blast effects on structures with TNT equivalence is certainly usable.Formby and Wharton conducted a ballistic mortar test on several commercial explosives to obtain their TNT equivalence [21].Analysis of criteria to determine TNT equivalence has been performed using several different calculation methods for different explosive materials and charge shapes [22].Attempts by several authors to correlate the TNT equivalent and different detonation parameters showed that there is no certain method for determining the TNT equivalent [23,24].Most regularly, an equivalent TNT mass links the mass of a certain explosive to the equivalent mass of TNT using the ratio of their detonation heat [25].Locking used several different equations with different detonation parameters to theoretically estimate the TNT equivalent for several standard high explosives [23].Numerical studies of TNT equivalent were performed for various high explosives for "incident and normally reflected peak overpressure and impulse" [26].Air blast TNT equivalence has been achieved for a variety of commercial explosives using ballistic mortar tests with 12 piezo-electric gauges, with two different calculating methods [27].
The TNT equivalence can also be estimated from the thermochemical equilibrium code results [28].Thirty years ago, Copper [20] claimed that the thermochemical computer codes could not correctly predict blast waves due to an absence of "chemically reactive equations of state".However, thanks to the progress in the field of thermochemical calculations and the development of more accurate equations of state, today's computer codes are capable of predicting the detonation performance and TNT equivalence quite accurately [24].Using thermochemical calculations, Jeremic and Bajic derived the TNT equivalent for a larger number of high explosives [29].
This study, as a part of the above-mentioned research program, aims to improve further the precision of blast-induced seismic effect calculation using the TNT equivalence approach, whereby TNT equivalents are calculated using the detonation parameters calculated by the thermochemical code EXPLO5.By doing so, we strive to contribute to the development of safer and more effective blasting techniques in mining projects, construction, and other relevant fields.

Materials and Methods
For this study, the original experimental data of Škrlec et al. [30] were taken in electronic form.To determine the blast-induced seismic effect of different explosives, the authors varied the mass of the charge weight per delay while the depth of the borehole and the size of the charge itself (cartridge) were kept constant.The test blast was carried out in the diabase rock, which is characterized by high compressive strength (Table 1).The explosives tested included 5 boreholes filled by EME (Emulsion explosive) charges and 5 boreholes for each EMM/EPS (Emulsion matrix/Expanded Polystyrene) mixture.The boreholes were drilled with 32 mm diameter and 240 mm depth.The charges consisted of PVC confinement 25 mm in diameter.Each explosive charge was set off separately using electric detonators of energy as reference detonator No. 4.
The ground oscillation velocities were measured for each borehole.The measurement setup is shown in Figure 1.Instantel BlastMate III and Instantel MiniMate Plus units were used to measure the ground oscillation velocity.A total of six units were used, with three units per line spaced 2 m apart.
development of safer and more effective blasting techniques in mining projects, construc-tion, and other relevant fields.

Materials and Methods
For this study, the original experimental data of Škrlec et al. [30] were taken in electronic form.To determine the blast-induced seismic effect of different explosives, the authors varied the mass of the charge weight per delay while the depth of the borehole and the size of the charge itself (cartridge) were kept constant.The test blast was carried out in the diabase rock, which is characterized by high compressive strength (Table 1).The explosives tested included 5 boreholes filled by EME (Emulsion explosive) charges and 5 boreholes for each EMM/EPS (Emulsion matrix/Expanded Polystyrene) mixture.The boreholes were drilled with 32 mm diameter and 240 mm depth.The charges consisted of PVC confinement 25 mm in diameter.Each explosive charge was set off separately using electric detonators of energy as reference detonator No. 4.
The ground oscillation velocities were measured for each borehole.The measurement setup is shown in Figure 1.Instantel BlastMate III and Instantel MiniMate Plus units were used to measure the ground oscillation velocity.A total of six units were used, with three units per line spaced 2 m apart.The maximum PPV from each measurement is a function of the scaled distance (SD).The general form of dependence of PPV and SD is given by the equation [31]: where SD is the scaled distance (m/kg 1/2 ), H is the constant in blast design, and β is the attenuation constant.The constants H and β are defined by test blasts.The scaled distance is given by Equation (2).
where R is the distance (m) and W is the charge weight per delay (kg).The maximum PPV from each measurement is a function of the scaled distance (SD).The general form of dependence of PPV and SD is given by the equation [31]: where SD is the scaled distance (m/kg 1/2 ), H is the constant in blast design, and β is the attenuation constant.The constants H and β are defined by test blasts.The scaled distance is given by Equation (2).
where R is the distance (m) and W is the charge weight per delay (kg).
It should be noted that, in the case of using different types of explosives for blasting works within the same project, the relationship of the effect of detonation properties of diverse explosives on blast-induced seismic effect, i.e., PPV, can be achieved via the equivalent mass approach.This method involves relating the effects of the yield of a certain explosive to that of a reference explosive (usually TNT): where W e (kg) present the TNT equivalent weight, i.e., the weight of TNT that produces the same effect at the same distance as an actual explosive, W (kg) is the weight of the actual explosive, and TNT e is the TNT equivalent for the actual explosive expressed in percentage.

Thermochemical Calculations
The detonation properties of analyzed explosives are computed using the EXPLO5 thermochemical code [32].The code predicts detonation properties (e.g., detonation temperature, detonation velocity, pressure, the energy of detonation, heat, etc.) by applying the "chemical equilibrium steady-state Chapman-Jouguet (C-J) detonation model".The equilibrium composition of detonation products the code calculates by mathematically describing the state of equilibrium in a multicomponent and multiphase system, applying the free energy minimization technique and mass balance principle.This method was originally developed by White et al. [33] and later adapted for computer application by Mader [34].The method is based on the fact that in chemical equilibrium, the chemical potential of reaction products is equal to the chemical potential of reactants, i.e., the fact that in the equilibrium state, the free energy of products has the minimum value.
EXPLO5 has the potential to use a number of equations of the state of gaseous detonation products.For this analysis, the modified Becker-Kistiakowsky-Wilson (MBKW) equation of state is used, given that earlier research demonstrated that it best explains the behavior of low-density explosives [35].
Different approaches have been used to estimate the strength of an explosive relative to TNT based on the results of thermochemical calculations.Basically, all the approaches are based on the interpretation of the relative strength in terms of certain detonation properties [23,24].Most frequently, the heat of detonation (Q) is used to estimate the TNT equivalent [23]: In [24], the authors analyzed the accuracy of several approaches for the evaluation of relative strength in terms of the Trauzl test and found that the best concurrence between experimental results and calculation is obtained using the product between the heat of detonation (Q) and the square root of the volume of detonation products (V 0 ): In this study, we also used the detonation energy, E 0 (i.e., maximum available energy for doing mechanical work), the detonation pressure (p CJ ), and velocity of detonation (D) to estimate the TNT equivalent: × 100 ( 6) × 100 ( 7) The detonation energy in Equation ( 6) is derived from the expansion isentrope of detonation products as described in [36].

Results and Discussion
The parameters of detonation for explosives used in this study are calculated using the EXPLO5 thermochemical code, applying the Chapman-Jouguet detonation model.The calculated detonation parameters, along with the experimentally determined detonation velocities, are given in Table 2.It can be seen from Table 2 that experimental and calculated detonation velocities for EME differ by less than 1%, while in the case of EMM/EPS mixtures, the difference increases as the amount of EPS increases, reaching ~25% for the mixture containing 80% of EPS.As shown in our previous studies [35,37,38], EXPLO5 can predict the detonation velocity of a large series of different explosives with an error of less than 5%, so a large discrepancy observed between experimental and calculated detonation velocities for mixtures containing a larger amount of EPS points to the conclusion that these explosives behave non-ideally.Among other things, the non-ideal behavior manifests in a decrease in detonation velocity with an explosive charge diameter decrease, and ideal detonation codes always predict larger detonation velocities than the experimental ones.
Recently, Dobrilović et al. [24] have demonstrated that ideal detonation codes can be used to predict the TNT equivalent and relative strength of explosives from the heat of detonation and the volume of detonation products for both ideal and nonideal explosives.In light of this finding, and aware of the limitation of ideal detonation codes when it comes to accurate prediction of D and p of commercial explosives, in this research, we used the ideal detonation code, but we calculated the TNT equivalent using different detonation parameters.
Although, in most studies, TNT is used as a reference explosive when comparing the effect of the yield of a certain explosive, in this study we have chosen EME as a reference explosive.The main reason for choosing EME is that both EME and EMM/EPS mixtures belong to the same class of explosives-commercial explosives.This is important because, according to Marshall [39], the equivalency principle can provide suitable results for an equal class of explosives, i.e., explosives that perform in the same manner.Compared to TNT, which behaves ideally, explosives used in this study exhibit non-ideal behavior and have significantly lower values of detonation parameters.
With regards to the above, Equations ( 4)-( 8) change in such a way that the detonation parameters for EME are used instead of those for TNT, and the TNT equivalent (TNTe) is replaced by EME equivalent (EMEe).The EME equivalents for EMM/EPS mixtures are calculated based on the detonation parameters given in Table 2 and using modified Equations ( 4)- (8).The results are given in Table 3.
It is evident from Table 3 that Equations ( 4)-( 6) give fairly similar EMEe values (difference up to 7.3%), while Equation (8) gives significantly lower values of EMEe.On the other hand, the EMEe values calculated by Equation ( 7) greatly deviate from the values obtained by other equations.Which of the equations given in Table 3 gives better results, i.e., which detonation parameter enables the most accurate estimation of the peak particle velocity by applying the equivalence principle, will be discussed below.
In the experiment [30] carried out using EME explosives, the authors measured PPV as a function of charge weight and distance from the blast (Figure 1).For analysis, only measurements from Line 1 were used since the PPV values for Lines 1 and 2 are similar, with small discrepancies due to microgeological conditions and different directions.The maximum PPV from each blast and monitoring position (MO) for EME are given in Table 4.The dependence of experimentally measured PPV on the explosive mass and distance from the blast is expressed through the scaled distance (Equation (1)).Taking the logarithm of both sides in Equation (1) yields the following equation: When plotted in log-log coordinates, Equation ( 9) gives a linear PPV-SD relationship of the form: where y = log(PPV), a = log(H), x = log(SD), and b = β.One can easily determine the constants a and b in Equation (10) and the constants H and β in Equation ( 9) through the linear regression analysis.The slope of line (b) is computed from the expression [40]: and intercept a is computed from the expression: The coefficient of correlation r is obtained by evaluating the expression [40]: and the sample standard deviation is evaluated from: where n is the number of points, x i and y i are the ith variable in the set i = 1 to n, and x and y are the means of each of the variables in the given set [40].
Based on the experimental data given in Table 4, the fitting constants a and b (Equation ( 10)) for EME are derived using the Instantel Blastware application.The application works by importing data directly from records and by manually entering the values of the actual charge weight per delay (W) and the distance of the monitoring points from the explosion (R) for each MO into the application.The values of PPV are recorded by the instruments at the monitoring points, based on which the application calculates the fitting constants (best-fit line) and additionally calculates and presents a 95% confidence line with equation (upper limit).The output results for monitoring points in Line 1 for EME are shown in Figure 2.
and the sample standard deviation is evaluated from: where n is the number of points, xi and yi are the ith variable in the set i = 1 to n, and x and y are the means of each of the variables in the given set [40].
Based on the experimental data given in Table 4, the fitting constants a and b (Equation ( 10)) for EME are derived using the Instantel Blastware application.The application works by importing data directly from records and by manually entering the values of the actual charge weight per delay (W) and the distance of the monitoring points from the explosion (R) for each MO into the application.The values of PPV are recorded by the instruments at the monitoring points, based on which the application calculates the fitting constants (best-fit line) and additionally calculates and presents a 95% confidence line with equation (upper limit).The output results for monitoring points in Line 1 for EME are shown in Figure 2.
95% confidence line, upper limit: 95% confidence line, lower limit: Once the fitting constants for the reference explosive EME are determined, the equivalence principle can be applied to predict the blast effect, i.e., PPV value at a specified distance and charge weight, for any explosive, provided its EMEe is known.For Since the Instantel Blastware application does not give the fitting constants a and b for the best-fit line, the same is calculated by CurveExpert-curve fitting and data analysis code.The constants H and β in Equation ( 9), calculated from constants a (a = log(H)) and b (b = β), and final equations for the best-fit line and 95% confidence lines are: Best-fit line: PPV EME = 395.64•SD−1.161 (σ = 7.306 mm /s, r = 0. 8357) 95% confidence line, upper limit: 95% confidence line, lower limit: Once the fitting constants for the reference explosive EME are determined, the equivalence principle can be applied to predict the blast effect, i.e., PPV value at a specified distance and charge weight, for any explosive, provided its EME e is known.For illustration, let us assume we want to predict the PPV, which will be generated by the detonation of 15.5 g of an explosive X at a distance of 6 m.The EME e of explosive X equals 38.5%.The equivalent weight of EME (W e ), which gives the same effect as 15.5 g of explosive X, can be obtained by Equation (3) (W e = W•EME e /100 = 5.97 g).The scaled distance for EME equals SD(W e , R) = R/ √ W e = 6/ √ 5.97/1000 = 77.7 m/kg 0.5 .By substituting the value of SD (W e , R) in Equation ( 15), we obtain PPV = 2.53 mm/s.Therefore, 15.5 g of explosive X will produce PPV = 2.53 mm/s at a distance of 6 m, the same as 5.97 g of EME at the same distance.Similarly, we can predict a mass of mixture X that will give the same value of PPV, at 6 m distance, as 72 g of reference explosive EME.In this case W e = 72 g, SD(W e, R) = 22.3 m/kg 0.5 , and the calculated PPV for that SD (Equation ( 15)) will be 10.7 mm/s.The mass of mixture X, which will produce PPV = 10.7 mm/s at R = 6 m, equals W = W e /(EME e /100) = 187 g.
To validate the above-described equivalency approach to predict blast-induced seismic effects, we used the experimentally measured PPV vs. charge weight and distance data from Škrlec et al. [30] for four mixtures based on an emulsion matrix (EMM) and expanded polystyrene (EPS).It is well known that the performance (or power) of explosives determined by different tests does not follow the same ranking order.For example, for blasting work, the detonation heat and the amount of detonation products are the most significant parameters.However, if an explosive performs strong disintegration work, then the most important parameters are the velocity of detonation and pressure [34].Maienschein [28] claims that the detonation energy calculated by thermochemical codes allows assessment of the TNT equivalency with respect to peak pressure, while the heat of detonation allows estimation of TNT equivalency with respect to quasi-static pressure at a long time.However, we could not find any previous research on which detonation parameter is best to use to derive TNT or any other explosive equivalence concerning PPV values.
Considering the above-mentioned, the first step in the validation of the approach was to find which detonation parameter, i.e., which equation for the calculation of EME e (Equations ( 4)-( 8)), correlates best with the measured PPV.The analysis showed that EME e derived from Q, E 0 , and QV 0 0.5 (Equations ( 4)-( 6)), give similar results and all of them overpredict PPV, EME e derived from p CJ (Equation ( 7)) greatly underpredicts PPV, while EME e derived from detonation velocity (Equation ( 8)) best reproduce experimental PPV-SD data for all studied explosives.For illustration, the results of the analysis for EMM/EPS-40/60 mixture are given in Figure 3. Considering that, further validation is achieved using EME equivalents calculated by Equation (8).detonation of 15.5 g of an explosive X at a distance of 6 m.The EMEe of explosive X equals 38.5%.The equivalent weight of EME (We), which gives the same effect as 15.5 g of explosive X, can be obtained by Equation ( 3 .By substituting the value of SD (We, R) in Equation ( 15), we obtain PPV = 2.53 mm/s.Therefore, 15.5 g of explosive X will produce PPV = 2.53 mm/s at a distance of 6 m, the same as 5.97 g of EME at the same distance.Similarly, we can predict a mass of mixture X that will give the same value of PPV, at 6 m distance, as 72 g of reference explosive EME.In this case We = 72 g, SD(We,R) = 22.3 m/kg 0.5 , and the calculated PPV for that SD (Equation ( 15)) will be 10.7 mm/s.The mass of mixture X, which will produce PPV = 10.7 mm/s at R = 6 m, equals  =  /( 100) ⁄ = 187 g.To validate the above-described equivalency approach to predict blast-induced seismic effects, we used the experimentally measured PPV vs. charge weight and distance data from Škrlec et al. [30] for four mixtures based on an emulsion matrix (EMM) and expanded polystyrene (EPS).It is well known that the performance (or power) of explosives determined by different tests does not follow the same ranking order.For example, for blasting work, the detonation heat and the amount of detonation products are the most significant parameters.However, if an explosive performs strong disintegration work, then the most important parameters are the velocity of detonation and pressure [34].Maienschein [28] claims that the detonation energy calculated by thermochemical codes allows assessment of the TNT equivalency with respect to peak pressure, while the heat of detonation allows estimation of TNT equivalency with respect to quasi-static pressure at a long time.However, we could not find any previous research on which detonation parameter is best to use to derive TNT or any other explosive equivalence concerning PPV values.
Considering the above-mentioned, the first step in the validation of the approach was to find which detonation parameter, i.e., which equation for the calculation of EMEe (Equations ( 4)-( 8)), correlates best with the measured PPV.The analysis showed that EMEe derived from Q, E0, and QV0 0.5 (Equations ( 4)-( 6)), give similar results and all of them overpredict PPV, EMEe derived from pCJ (Equation ( 7)) greatly underpredicts PPV, while EMEe derived from detonation velocity (Equation ( 8)) best reproduce experimental PPV-SD data for all studied explosives.For illustration, the results of the analysis for EMM/EPS-40/60 mixture are given in Figure 3. Considering that, further validation is achieved using EME equivalents calculated by Equation (8).As a part of the validation, we compared the experimental PPV-SD(W) data given in Table 5 against the data calculated by Equation ( 15) using EME equivalent weight (PPV-SD(We)).The best-fit curve for EME explosives is also given for the comparison (Figure 4).As a part of the validation, we compared the experimental PPV-SD(W) data given in Table 5 against the data calculated by Equation ( 15) using EME equivalent weight (PPV-SD(W e )).The best-fit curve for EME explosives is also given for the comparison (Figure 4).
It can be observed from Figure 4 that a satisfactory agreement between the experimental and predicted PPV-SD data exists for all explosives.It should be noted that the experimental PPV-SD data for EME (Figure 2) exhibit a significant scattering of the results, especially in the area of lower SD, i.e., closer to the blast point.In this region, PPV exponentially increases with the decrease of SD, so a small variation in measured distance or explosive weight may cause a significant change in PPV, which results in larger experimental errors.For example, for SD ≈ 10 m/kg 0.5 , the PPV ranges between 48.5 mm/s and 27.3 mm/s (the difference is 21.2 mm/s, almost double).For EMM/EPS mixtures, for SD = 10-15 m/kg 1/2 , the minimum and the maximum PPV differ between 5 and 13.5 mm/s.With that in mind, it can be concluded that the predicted PPV values for EME/PPS mixtures are within experimental error (Figure 4).It can be observed from Figure 4 that a satisfactory agreement between the experimental and predicted PPV-SD data exists for all explosives.It should be noted that the experimental PPV-SD data for EME (Figure 2) exhibit a significant scattering of the results, especially in the area of lower SD, i.e., closer to the blast point.In this region, PPV exponentially increases with the decrease of SD, so a small variation in measured distance or explosive weight may cause a significant change in PPV, which results in larger experimental errors.For example, for SD ≈ 10 m/kg 0.5 , the PPV ranges between 48.5 mm/s and 27.3 mm/s (the difference is 21.2 mm/s, almost double).For EMM/EPS mixtures, for SD = 10-15 m/kg 1/2 , the minimum and the maximum PPV differ between 5 and 13.5 mm/s.With that in mind, it can be concluded that the predicted PPV values for EME/PPS mixtures are within experimental error (Figure 4).
To further validate the equivalency approach for blast-induced seismic effect prediction, we analyzed how accurately the weight of explosives that will generate a specified PPV value at a specified distance from the blast point can be predicted.For this, we took the PPV limit value from the HR DIN 4150:2011 standard [41] (Table 6) as the maximum permitted value for "dwellings and buildings of similar design and/or occupancy" (Line 2), which equals 20 mm/s.Table 6.Guideline values for vibration velocity to be used when evaluating the effects of short-term vibration on structures as per HR DIN 4150:2011 [41].To further validate the equivalency approach for blast-induced seismic effect prediction, we analyzed how accurately the weight of explosives that will generate a specified PPV value at a specified distance from the blast point can be predicted.For this, we took the PPV limit value from the HR DIN 4150:2011 standard [41] (Table 6) as the maximum permitted value for "dwellings and buildings of similar design and/or occupancy" (Line 2), which equals 20 mm/s.Table 6.Guideline values for vibration velocity to be used when evaluating the effects of short-term vibration on structures as per HR DIN 4150:2011 [41].The value of SD at which a specified PPV of 20 mm/s is reached for individual explosives is determined from experimental PPV-SD data as illustrated in Figure 5.The best-fit equations for EMM/EPS explosives are obtained by the regression analysis of experimental PPV-SD data given in Figure 4.The results of the regression analysis are summarized in Table 7.The value of SD at which a specified PPV of 20 mm/s is reached for individual explosives is determined from experimental PPV-SD data as illustrated in Figure 5.The best-fit equations for EMM/EPS explosives are obtained by the regression analysis of experimental PPV-SD data given in Figure 4.The results of the regression analysis are summarized in Table 7. From the best-fit line equations for individual explosives (Table 7), and with a PPV limit value of 20 mm/s taken from the standard, the SD is calculated according to the equation:

Line
The results of the calculation are summarized in Table 7. From the best-fit line equations for individual explosives (Table 7), and with a PPV limit value of 20 mm/s taken from the standard, the SD is calculated according to the equation: The results of the calculation are summarized in Table 7.
As expected, the SD shifts to lower values (i.e., closer to the explosion point) with a decrease in detonation performance.From the known values of SD corresponding to PPV = 20 mm/s, one can calculate the maximum permitted weight of individual explosives at a specified distance from the blast point using Equation (2).The results of the calculation are given in Table 8.
However, if experimental PPV-SD data are not available, one can apply the equivalency principle to estimate the maximum permitted weights at specified distances.This can be achieved in the following way: first, using the best-fit equation for reference explosive, the SD value at which the specified PPV value is reached is calculated using Equation (18).Then, from the known SD value (in this case for PPV = 20 mm/s, SD = 13.07 m/kg 1/2 for EME as reference explosive), the maximum permitted weight at a specified distance is calculated by Equation (2).Then, the maximum permitted weights of any other explosive can be calculated by Equation (3), applying the equivalency principle and knowing their EMEe.As visible from Table 8 and Figure 6, the weights predicted using the equivalency principle (WEP) are lower than those determined from the best-fit equation of experimental PPV-SD data.The difference between predicted and experimental PPV values goes from 17.6% for EMM/EPS-40/60 to 69.7% for EMM/EPS-20/80.In our opinion, there are two key reasons for such large differences.The first is generally a large scattering of measured PPV values, which could primarily be attributed to the small amount of explosive used in the experiment (15-42 g for EMM/EPS mixtures).

Conclusions
The possibility of the application of the equivalence approach in the prediction of blast-induced seismic effects is studied in this work.After a thorough analysis of the results, the following conclusions were drawn: (a) The results have shown that the equivalency approach is a potential tool in the assessment of the blast-induced seismic effect.It is especially useful in the situation when the type of explosive is changed within a single project; (b) It was demonstrated that the best results are obtained when the equivalent weight of the studied explosive is calculated as a ratio of detonation velocities of studied and The second reason lies in the fact that the PPV range in the experiments was significantly below the PPV limit value of 20 mm/s.For example, the maximum mean PPV value for EMM/EPS-20/80 mixture at SD ≈ 16 m/kg 1/2 equals only 6.9 mm/s, for EMM/EPS-30/70 mixture at SD ≈ 13 m/kg 1/2 equals 11.8 mm/s, etc.That means that SD which corresponds to PPV = 20 mm/s (Equation ( 18)), calculated by the best-fit equation, is

Figure 2 .
Figure 2. The best-fit line and corresponding 95% confidence line for EME.Since the Instantel Blastware application does not give the fitting constants a and b for the best-fit line, the same is calculated by CurveExpert-curve fitting and data analysis code.The constants H and β in Equation (9), calculated from constants a (a = log(H)) and b (b = β), and final equations for the best-fit line and 95% confidence lines are: Best-fit line:  = 395.64•  .(σ= 7.306 mm/s, r = 0.8357)(15)

Figure 2 .
Figure 2. The best-fit line and corresponding 95% confidence line for EME.

Figure 3 .
Figure 3.Comparison of experimental and predicted values of PPV for EMM/EPS-40/60 mixture.

Figure 3 .
Figure 3.Comparison of experimental and predicted values of PPV for EMM/EPS-40/60 mixture.
Note: SD (W, R)-scaled distance for specified actual weight (W) and distance R, SD (W e , R)-scaled distance for equivalent weight (W e ) and distance R.

Figure 4 .
Figure 4. Comparison of predicted and experimental PPV values for EMM/EPS mixtures.(Note: Red dashed curves represent the best-fit curve for EME calculated by Equation (15) using actual weight, solid blue curves represent the best-fit curve for EMM/EPS mixtures calculated by Equation (15) using EME equivalent weight, and solid green curves represent the best-fit of experimental data).

1 BuildingsFigure 4 .
Figure 4. Comparison of predicted and experimental PPV values for EMM/EPS mixtures.(Note: Red dashed curves represent the best-fit curve for EME calculated by Equation (15) using actual weight, solid blue curves represent the best-fit curve for EMM/EPS mixtures calculated by Equation (15) using EME equivalent weight, and solid green curves represent the best-fit of experimental data).

Figure 5 .
Figure 5. Graphical presentation of determination of SD at which a specified PPV value is reached.

Figure 6 .
Figure 6.Comparison of predicted and experimental W-R curves for fixed PPV value of 20 mm/s.

Figure 6 .
Figure 6.Comparison of predicted and experimental W-R curves for fixed PPV value of 20 mm/s.

Table 2 .
Detonation parameters of used explosives.D expt experimental detonation velocity measured in steel tubes 21.5 mm inner diameter, 150 mm long, 2.7 mm wall thickness, D, p CJ , Q, E 0 , and T are calculated detonation velocity, pressure, heat, energy, and temperature, respectively. *

Table 3 .
Calculated EME equivalents for tested explosives.

Type of Structure Vibration at the Foundation at a Frequency of Vibration at Horizontal Plane of Highest Floor at All Frequencies 1 Hz to 10 Hz 10 Hz to 50 Hz 50 Hz to 100 Hz *)
At frequencies above 100 Hz, the values given in this column may be used as minimum values.
At frequencies above 100 Hz, the values given in this column may be used as minimum values.

Table 7 .
Best-fit line equations and scaled distance corresponding to PPV = 20 mm/s for studied Graphical presentation of determination of SD at which a specified PPV value is reached.

Table 7 .
Best-fit line equations and scaled distance corresponding to PPV = 20 mm/s for studied explosives.

Table 8 .
Experimental and predicted maximum permitted weights at different distances that generate PPV = 20 mm/s.: W fit , W EP -maximum permitted weights determined from experimental data using best-fit curve and estimated applying the equivalency principle, respectively. Note