## 2 Materials and Methods  
For this study, the original experimental data of [PERSON] 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]:  
\[PPV=H(SD)^{-\beta} \tag{1}\]  
where SD is the scaled distance (m/kg\({}^{1/2}\)), H is the constant in blast design, and \(\beta\) is the attenuation constant. The constants H and \(\beta\) are defined by test blasts. The scaled distance is given by Equation (2). \[SD=\frac{R}{\sqrt{W}} \tag{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):  
\[W_{e}=W.\frac{TNT_{e}}{100} \tag{3}\]  
where W\({}_{\text{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\({}_{\text{e}}\) is the TNT equivalent for the actual explosive expressed in percentage. \begin{table}
\begin{tabular}{c c c} \hline \hline  & **Diabase Rock** & \\ \hline Compressive strength & 200–400 & MPa \\ Density & 2.85–3.15 & g/cm\({}^{3}\) \\ Water absorption & 0.2–1 & \% mass \\ Spatial mass & 2.8–3.1 & g/cm\({}^{3}\) \\ Porosity & 0.1–1 & \% volume \\ \hline \hline \end{tabular}
\end{table} 
Table 1: Characteristics of diabase rock [30]. Figure 1: Ground oscillation velocity measurement setup for Line 1 [30]. 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 [PERSON] et al. [33] and later adapted for computer application by [PERSON] [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.