Next Article in Journal / Special Issue
Determination of Seismic Safety Zones during the Surface Mining Operation Development in the Case of the “Buvač” Open Pit
Previous Article in Journal
Prismatic to Asbestiform Offretite from Northern Italy: Occurrence, Morphology and Crystal-Chemistry of a New Potentially Hazardous Zeolite
Previous Article in Special Issue
Research on Power Plant Ash Impact on the Quality of Soil in Kostolac and Gacko Coal Basins

Article

# Models of Determining the Parameters of Rock Mass Oscillation Equation with Experimental and Mass Blastings

Faculty of Mining and Geology, Djušina 7, 11060 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Minerals 2018, 8(2), 70; https://doi.org/10.3390/min8020070
Received: 21 December 2017 / Revised: 1 February 2018 / Accepted: 7 February 2018 / Published: 16 February 2018

## Abstract

The explosion caused by detonation of explosive materials is followed by release of a large amount of energy. Whereby, a greater part of energy is used for rock destruction, and part of energy, in the form of seismic wave, is lost in the rock mass causing rock mass oscillation. Investigations of the character and behavior of the pattern of seismic wave indicate that the intensity and nature of the seismic wave are influenced by rock mass properties, and by blasting conditions. For evaluation and control of the seismic effect of blasting operations, the most commonly used equation is that of M.A. Sadovskii. Sadovskii’s equation defines the alteration in the velocity of rock mass oscillation depending on the distance, the quantity of explosives, blasting conditions and geological characteristics of the rock mass, and it is determined based on trial blasting for a specific work environment. Thus, this paper offers analysis of the method for determination of parameters of the rock mass oscillation equation, which are conditioned by rock mass properties and blasting conditions. Practical part of this paper includes the experimental research carried out at Majdanpek open pit, located in the northern part of eastern Serbia and the investigations carried out during mass blasting at Nepričava open pit, located in central Serbia. In this paper, parameters n and K from Sadovskii’s equation were determined in three ways—models in the given work environment. It was noted that, in practice, all three models can be successfully used to calculate the oscillation velocity of the rock masses.

## 1. Introduction

When the seismic wave comes across a certain point in the rock mass, it then disturbs the balance of rock mass particles at that point. The particles oscillate around their balance position for a certain period until they return to a complete standstill. Oscillation of the rock mass particles is something that is manifested and felt as shock, i.e., rock mass vibration. Shock or rock mass oscillation can be a consequence of natural processes or human activities (in the performance of mining). Measurements and monitoring of seismic processes can be done in different ways, depending on the cause of the shock .
Intensity of the shock may be determined by measuring one of the three basic parameters that characterize oscillation of the aroused rock mass, and those parameters are: displacement of the rock mass particles, oscillation velocity of the rock mass particles and acceleration of the aroused environment. Sizes of these three parameters point to the intensity of the force which caused them, as well as to the degree of danger caused by the shock. The parameter that is most often applied to estimate seismic effect of the blasting is the oscillation velocity of the rock mass. It is considered that this is the parameter which best defines danger from shock and damages it may cause. Due to that, adequate standards for the protection from shock are based on data that refer to oscillation velocity [2,3].
As the relation between the rock mass oscillation velocity and basic parameters affecting its magnitude, being: the amount of explosive, the distance from the blast site, characteristics of the rock material and the type of blasting, the equation of M.A. Sadovskii, where the oscillation velocity v is given in the form of the function, is most frequently used:
$v = K · R − n ,$
here R represents reduced distance, meaning the distance from the blast site to the monitoring point r, reduced to the used explosive amount Q. K and n parameters conditioned by soil characteristics and blasting conditions, thereby v is the decreasing convex function of the variable R.
By applying the equation of rock mass oscillation while blasting, the determination of the rock mass oscillation velocity is made possible for each blast operation in advance, thus blasting is, as regards seismic effect, under control, giving opportunity to plan the magnitude of shock waves for each future blast operation . In this way adverse blasting effects are reduced. Adverse effects of blasting imply, in addition to the seismic ones, those of air blast waves, fly rock, etc. Thus, production efficiency is increased and, at the same time, construction and mining facilities, as well as the environment near the blast site are protected.

## 2. Recording Rock Mass Oscillation Velocity

Intensity of the shock is measured by instruments—seismographs, in the field. The main parts of a seismograph are sensors—geophones by way of which rock mass oscillation velocity is being recorded. Measuring of the rock mass oscillation velocity needed for this paper has been done by the instrument type Vibraloc, produced by the Swedish firm ABEM. Detector for detecting oscillation contains three seismometers placed in a common housing oriented in directions X, Y and Z from the rectangular coordination system. Seismograph is constructed to be able to record oscillation velocity, at the same time it is possible to calculate acceleration and displacement of the rock mass. The instrument also records frequency value. Recording of air shocks may be done on the fourth channel using a microphone that is connected to the external connector. Basic characteristics of the seismograph Vibraloc are:
• Type and designation of instrument: Vibraloc
• Manufactured: in Sundbyberg, Sweden
• Measuring possibilities: speed, acceleration displacement
• Number of channels: three channels—speed, fourth channel—air channel
• Number of components: three—transversal, vertical, side
• Frequency range: 2–250 Hz
• Location possibility: flat floors, boards, foundations, ground
• Trigger level, V, L, T (trigger levels): 0.1–200 mm/s
• Channel trigger level A (microphone): 2–150 Pa
• Sampling: 1000, 2000 or 4000 Hz
• Recording length: 1–100 s or automatic length
• Data transfer: Vibraloc PC software
• Data analysis: UVSZ and UVSZA software (Vibraloc, Sundbyberg, Sweden)
At the monitoring point, waves that arise as the consequence of blasting, appear almost simultaneously. They also cause the movement of rock mass particles in different directions. Material rock mass particles oscillate along a very complex path. To define intensity and direction of the action of shocks at a certain point of time, it is necessary to record them in three directions, and that by way of:
• horizontal component vt, normal to the direction of wave spreading (transversal component)
• horizontal component vl, in direction of wave spreading (longitudinal component) and
• vertical component vv.
Based on three oscillation velocity components resulting rock mass oscillation velocity may be determined $v r e z = v t 2 + v l 2 + v v 2 ,$ in the equation of Sadovskii: oscillation velocity.
In Figure 1 is shown the recording of oscillation velocity of the ground intended for blasting No. I—measuring point No. 3 at PK Nepričava.

## 3. Equation of Rock Mass Oscillation

To establish the correlation between the oscillation velocity and three basic parameters affecting its size: the explosive quantity, properties of rock material and the distance, there have been developed several mathematical models in the world. One of the most frequently used models, i.e., equations, is the equation of Sadovskii defining the law on velocity alteration of soil oscillation depending on the distance, the explosive quantity, and the way of blasting [5,6]. The equation defined in this way offers the possibility to determine the seismic effect of blasting towards a facility or a settlement, whereby the connection, between the velocity of soil oscillation and consequences that can affect facilities, is used. The equation of M.A. Sadovskii  is given in the form:
$v = K · R − n = K · ( r Q 3 ) − n$
where there are:
v—the velocity of rock mass oscillation (cm/s)
K—the coefficient conditioned by soil characteristics as well as blasting conditions where the explosive amount is given by way of the volume. K is being determined by terrain surveying,
n—the exponent, conditioned by soil properties and mining conditions and determined by field measurements as well,
r—the distance from the blast site to the monitoring point (m),
Q—the amount of explosive (kg),
R—reduced distance, expressed as $R = r Q 3$.

#### 3.1. Derivation of Equation of Rock Mass Oscillation

Sadovskii’s equation has been derived from the following condition: if the radius of the charge r0 and the distance from the blast site to monitoring point r increase in the same or approximately the same ratio, the soil oscillation velocity v remains the same. Sadovskii equation, given in Equation (1) may be expressed in the form :
$v = K v · ( r 0 r ) n$
which denotes:
v—the velocity of rock mass oscillation (cm/s)
Kv—the coefficient conditioned by rock mass characteristics as well as blasting conditions where the explosive amount is given through radius of the explosive charge. K is being determined by terrain surveying,
r—the distance from the blast site to the monitoring point (m),
r0—the radius of the explosive charge,
n—the exponent, conditioned by rock mass properties and mining conditions and determined by field measurements as well.
The radius of the explosive charge r0 and the amount of explosive Q are related by the equation:
By replacing the value r0 from the Equation (4), in the Equation (3), we obtain:
$v = K v ( 3 Q 4 π 3 r ) n = K v · ( 3 4 π 3 ) n · ( Q 3 r ) n = K v · K 1 · ( r Q 3 ) − n = K · ( r Q 3 ) − n = K · R − n ,$
whereas:
whereby:
R—the reduced distance is the distance from blasting point to observation point reduced to a quantity of explosive and given in the following form $R = r r 0$.
Thus, we have obtained the oscillation law of rock mass, i.e., the equation of Sadovskii.

#### 3.2. The Property of the Rock Mass Oscillation Equation

The rock mass oscillation equation, given in the Equation (1) may be written in the form:
$v ( R ) = K · R − n .$
Then from the Equation (7) for the derivation v′(R), we obtain:
$v ′ ( R ) = − n · K · R − n − 1 ,$
which may be written in the form:
$v ′ ( R ) = − n R · K · R − n .$
Having in mind Equation (7), from the Equation (9), we obtain:
$v ′ ( R ) = − n · v ( R ) R .$
If Lagrange’s equation is applied to the function v(R) from (7), in the final growth increment, we obtain:
For small absolute values h, from (11), we obtain the relation:
$v ( R + h ) ≈ v ( R ) + h · v ′ ( R ) .$
Having (10) in mind, relation (12) becomes:
$v ( R + h ) ≈ v ( R ) − n · h · v ( R ) R .$
A significant property of the equation of rock mass oscillation velocity depending on the reduced distance which is obtained from the relation (13) is as follows: if the reduced distance R from any level is increased (decreased) by 1%, the oscillation velocity of the soil v will decrease (increase) by approximately n% .
Indeed, for $h = ± R 100 = ± 1 % R$, from (12) we obtain:

#### 3.3. Models to Determine the Parameters of Equation of Rock Mass Oscillation Velocity

In the Equation (1) two parameters appear, K and n, which need to be determined for a specific work environment and specific blasting conditions. Given the property of the rock mass oscillation equation, it is possible to determine the parameters K and n in many ways, i.e., models, using the values obtained by experimental measurements.

#### 3.3.1. Model 1—Determining the Parameters by Applying the Least Square Method

The Least Square Method is mainly used to obtain the parameters K and n which represents a common model [9,10].

#### 3.3.2. Model 2—Determining the Parameters through the Application of Quotient of Products of the Same Number of Reduced Distances and Oscillation Velocities

According to this model, determination of the parameter n is based on the quotient of the product of equal number of experimental data of the rock mass oscillation velocity v, using the oscillation velocity equation. In this way, first, we determine the parameter n, and then with parameter n found, we determine the parameter K from the rock mass oscillation equation, which will be shown cases .
Let v1, v2, …, vN be the rock mass oscillation velocity values obtained for corresponding reduced distances R1, R2, …, RN. Proceeding from the equation for the rock mass oscillation velocity, given by the Equation (1) we can assume, that for pairs (Rm, vm) the following relations are valid:
from where we obtain the relation:
$v 1 · v 2 · … · v s v s + 1 · v s + 2 · … · v 2 s = K s · ( R 1 · R 2 · … · R s ) − n K s · ( R s + 1 · R s + 2 · … · R 2 s ) − n ,$
i.e., relation:
$v 1 · v 2 · … · v s v s + 1 · v s + 2 · … · v 2 s = ( R 1 · R 2 · … · R s R s + 1 · R s + 2 · … · R 2 s ) − n , 2 s ≤ N .$
From relation (17) by the logarithm operation we obtain:
$− n = log v 1 · v 2 · … · v s v s + 1 · v s + 2 · … · v 2 s log R 1 · R 2 · … · R s R s + 1 · R s + 2 · … · R 2 s .$
In this way we obtained the value of the parameter n, which we will include into the Equation (1). Now we have the Equation (1) in which there is the unknown parameter K.
If we insert the values of pairs (Rm, vm), m = 1, 2, …, N, into the Equation (1), we will obtain one after the other:
$K 1 = v 1 · R 1 n , K 2 = v 2 · R 2 n , … K N = v N · R N n ,$
wherein for the parameter K, we take the mean value, i.e.,
$K = K 1 + K 2 + … + K N N .$
With the values of n and K obtained in this manner, the rock mass oscillation equation is reduced to: $v = K · R − n$.

#### 3.3.3. Model 3—Determining the Parameters Using the Linear Connection between the Logarithms of Reduced Distances and the Logarithms of the Oscillation Velocity

To obtain a linear connection between reduced distance R and rock mass oscillation velocity v, the Equation (1) was converted to logarithm. In the observed experiment, in the coordinate system logR-0-logv, certain points (logRi, logvi) will be below logR axis. In this case, we will increase the values of the logv ordinates by a positive number, so that all points are above the abscissa axis logR .
Such a procedure offers possibility to determine of parameters K and n from the equation for rock mass oscillation in a completely new way.
If the Equation (1) is converted to logarithm, we will obtain equation:
$log v = log K − n · log R .$
By introducing the change:
Equation (21) becomes a line equation:
$y = a x + b .$
We will increase ordinate values $y = log v$ by one positive number c, so that all points $M i ( x i , y i + c )$ are above axis $x = log R$, i.e., so that $y i + c = loh v i + c > 0$. In this manner, the Equation (23) becomes equation:
$y + c = a x + b + c .$
In practice, for c the following value can be taken:
$c = 3 .$
In this case Equation (24) becomes equation:
$y + 3 = a x + b + 3 .$
Let us write Equation (26) in the form:
$y + 3 = A x + B ,$
where:
$B = b + 3 .$
If we connect points $M i ( x i , y i + 3 )$ by line segments $M i M i + 1$(i = 1, 2, ..., N − 1) we will get one broken line. In so doing, the area bounded by x-axis, line segments $x 1 M 1$ and $x N M N$ and broken line, consisting of (N − 1) trapezoid [12,13], has the value:
$S T = 1 2 [ ∑ i = 1 N − 1 ( y i + 3 + y i + 1 + 3 ) · ( x i + 1 − x i ) ] .$
Let’s mark the area bounded by line (27), lines and x-axis with S. Its value is:
$S = ∫ x 1 x N ( A x + B ) d x = A 2 ( x N 2 − x 1 2 ) + B ( x N − x 1 ) .$
Let us assume that:
$S = S T ,$
which means that:
$A 2 ( x N 2 − x 1 2 ) + B ( x N − x 1 ) = S T .$
In this manner we have obtained one link between coefficients A and B of the equation of line (27).
From Equation (32), it’s:
$B = S T x N − x 1 − A ( x N + x 1 ) 2 .$
By replacing the value for B from the Equation (33) into Equation (27), we obtain:
$y + 3 = A ( x − x N + x 1 2 ) + S T x N − x 1 .$
From Equation (34), we have:
$A = y + 3 − S T x N − x 1 x − x N + x 1 2 .$
Based on obtained experimental data of oscillation velocities v1, v2, ..., vN and appropriate reduced distances R1, R2, ..., RN and having in mind relations (22) and (27), we will calculate the value of ST according to Equation (29). For each pair (xi, yi + 3) from Equation (35), we will determine the parameters Ai, i = 1, 2, ..., N, and in so doing, for parameter A, we will take the arithmetic mean value of parameters Ai:
$A = 1 N · ∑ i = 1 N A i .$
If the value of parameter A obtained according to Equation (36), is replaced in the Equation (33), we will obtain the value of parameter B. Based on the obtained values of parameters A and B and considering the relations (22), the Equation (27) becomes:
$y + 3 = A x + S T x N − x 1 − A ( x N + x 1 ) 2 ,$
i.e., equation:
$log v + 3 = A log R + S T x N − x 1 − A ( x N + x 1 ) 2 ,$
where from we obtain:
$log v = A log R + S T x N − x 1 − A ( x N + x 1 ) 2 − 3 .$
If we put that:
$S T x N − x 1 − A ( x N + x 1 ) 2 − 3 = log K ,$
Equation (39) maybe written in the form:
i.e., in the form:
$log v = log K · R A ,$
wherefrom we obtain:
$v = K · R A .$

#### 3.4. Defining Statistical Criteria

For the above-mentioned models 1, 2 and 3, based on experimental data, we have obtained equations, which make possible to determine the oscillation velocities of the rock mass v depending on the reduced distance R. To assess the degree of connection between v and R, we have used the curved line dependency index ρ [14,15].
The evaluation of the relationship degree of two variables  to values of the curved line dependency index ρ is given in the following survey:
• 0.0 < ρ < 0.2—none or highly poor correlation,
• 0.2 < ρ < 0.4—poor correlation,
• 0.4 < ρ < 0.7—significant correlation,
• 0.7 < ρ < 1.0—strong or highly strong correlation.
In order to make sure that the conclusions we reach based on the sample concerning the characteristics of the basic data set are credible, the sample needs to be representative. The sample will be representative, if it is, by its structure, like the basic set .

## 4. Methodology and Results of Research

Practical part of this paper includes:
• Experimental research carried out in Majdanpek open pit, situated in northern part of east Serbia and,
• Investigations carried out during mass blasting in Nepričava open pit, located in central Serbia.

#### 4.1.1. General Characteristics of the Work Environment in Majdanpek Open Pit

In Majdanpek open pit we performed experimental investigations on copper ore . The basic physical and mechanical properties of the rock mass were determined on samples in the laboratory. By examining the physical and mechanical properties of the work environment, the following values were obtained:
• compressive strength, σp = 109 MPa
• bending strength, σs = 19 MPa
• tensile strength, σI = 13 MPa
• cohesion, C = 7.5 MPa
• strength coefficient, f = 11
• volume density, γ = 26 kN/m3
• angle of internal friction, φ = 50°
• porosity, p = 3.6%

#### 4.1.2. Method of Blasting in Majdanpek Open Pit

The measurements of shocks resulting from blasting in Majdanpek open pit were performed during blasting using half-second electric detonators. Delay time of initiation between boreholes was 0.5 s, which led to ten explosions and appropriate rock mass oscillation velocities.
As explosive powdered ammonium nitrate for general purposes Amoneks-I manufactured by Trayal from Kruševac, Serbia was used. The holes were arranged in a single line while one cartridge of explosive was placed in each hole. Diameter of explosive cartridge was 28 mm, cartridge length 0.15 m, and cartridge weight 0.1 kg. Electric capacitor was used for initiation of the explosive .
During experimental investigations performed in Majdanpek open pit, the following blasting parameters were specified:
• depth of borehole: 0.5 m,
• weight of explosive charge per borehole: 0.10 kg,
• number of boreholes: 10,
• distance between boreholes: 1.0 m,
• distance between measuring point and the first borehole: 5.0 m,
• time delay between initiation of boreholes: 0.5 s.

#### 4.1.3. Calculation of Parameters of Rock Masses Oscillation Equation in Majdanpek Open Pit

The values of distance from the blasting point to the point of observation r, quantity of explosive Q, calculated values of reduced distances R, registered values of rock mass oscillation velocities by components vt, vv, vl and resulting rock mass oscillation velocities vrez for a total of ten explosions, are given in Table 1.
Based on data given in Table 1, the equation of rock mass oscillation is calculated by the Equation (1)—as per models 1, 2 and 3. The calculation of the curve was carried out for values of reduced distances from R = 10.7722 to R = 30.1621. Thus, curve parameters were calculated enabling us to determine the equation of rock mass oscillation in the form of:
● Model 1
$v 1 = 146.7047 R − 1.8326$
Graphic representation of the rock masses oscillation equation in Majdanpek open pit is given in Figure 2.
● Model 2
$v 2 = 158.5234 R − 1.8568$
● Model 3
Table 2 shows the data for calculating the parameters of rock masses oscillation equation in Majdanpek open pit according to Model 3.
According to Equation (29), we obtain:
$S T = 2.787876$
According to Equation (36), we find:
$A = − 1.71420 ,$
while according to Equation (33), we get:
$B = 7.66477 .$
If the obtained values of parameters A and B from (47) and (48) are put into Equation (27), we get:
$y + 3 = − 1.171420 x + 7.66477 .$
Since y = logv, x = logR, then the Equation (49) becomes to equation:
$log v + 3 = − 1.71420 log R + 7.66477 .$
Figure 3 graphically represents the values of logarithms of the rock mass oscillation velocities and the logarithms of the reduced distances, whereby the logarithms of the rock mass oscillation velocity are increased by number 3.
In diagram shown on Figure 3, the appropriate line is represented by Equation (50). From the Equation (50), we obtain:
$log v = − 1.71420 log R + 4.66477 .$
Since:
$4.66477 = log 106.1412 ,$
then Equation (51) can be written in the form:
$log v = log 106.1412 R − 1 , 7142 ,$
where from we obtain:
$v = 106.1412 R − 1.7142 .$
Based on the obtained Equations of soil oscillation (44), (45) and (54), it is possible to calculate values of rock mass oscillation velocities for corresponding reduced distances for models 1, 2 and 3.
In Table 3, is presented the survey of reduced distances R, recorded oscillation velocities vr, calculated oscillation velocities vi1, vi2, vi3 as well as the difference between recorded and calculated soil oscillation velocities for models 1, 2 and 3.
Based on the data in Table 3, a statistical analysis was carried out and the following values were obtained:
• for model 1:
The curved line dependency index ρ1 between the reduced distance R and rock mass oscillation velocity v is:
ρ1 = 0.9730 (there is a strong correlation between R and v, given in the Equation (44)),
• for model 2:
ρ2 = 0.9694 (there is a strong correlation between R and v, given in the Equation (45)),
• for model 3:
ρ3 = 0.9789 (there is a strong correlation between R and v, given in the Equation (54)).

#### 4.2.1. General Characteristics of Nepričava Open Pit

The Nepričava open pit is a part of the municipality of Lajkovac. It is situated on the left bank of the Kolubara River in Central Serbia, and belongs to Kolubara District.
The limestone deposit Nepričava consists of rock masses of sedimentary origin, not entirely homogeneous in composition. Limestones are the most important lithological members from both economic aspect and aspect of mining exploitation technology.
Physical and mechanical properties of limestone in the deposit are spatially very uneven. Examination of physical and mechanical properties of the working environment, submitted the following results:
• Compressive strength 21.1–93.0 MPa
• Volume weight 17.1–26.3 kN/m3
• Density 0.97–0.99 kN/m3
• Porosity 0.0004%–0.026%
• Water suction 0.08%–0.23%

#### 4.2.2. Blasting Method

Measurements of seismic shocks at the Nepričava open pit were performed during blasting, conducted for deposit exploitation . Three blasting operations were performed.
Balkanit 60/1000 and Videksil 65/2000 (Videks Company, Sofia, Bulgaria) were used as explosive materials. The activation of explosives in boreholes was performed using detonating fuse C-12, while initiating of detonating fuse was performed with blasting cap No. 8. decelerators between boreholes had the interval delay of 20 ms.
Basic data related to the number of boreholes Nb, the overall explosive amount Quk, the maximal explosive amount by deceleration interval Qi, overall borehole depth Luk and average stemming length Lpc, are presented in Table 4.

#### 4.2.3. Calculation of Rock Mass Oscillation Equation Parameters

Values of distances from blast sites to monitoring points r, the amount of explosive Q, calculated values of reduced distances R, recorded values of rock mass oscillation velocities by components vt, vv, vl and resulting oscillation velocities vrez for blasting from I to III of totally five measuring points MM are given in Table 5.
Based on data given in Table 5, the equation of rock mass oscillation is calculated by the Equation (1)—by the models 1, 2 and 3. The calculation of the curve was carried out for values of reduced distances from R = 15.2808 to R = 45.0858. Thus, curve parameters were calculated enabling us to determine the equation of rock mass oscillation in the form of:
● Model 1
$v 1 = 34.0616 R − 1.2663$
Graphic representation of the rock masses oscillation equation in Nepričava open pit is given in Figure 4.
● Model 2
$v 2 = 50.4220 R − 1.3730$
● Model 3
Table 6 shows the data for calculating the parameters of rock masses oscillation equation in Nepričava open pit according to Model 3.
According to Equation (29), we obtain:
$S T = 2.546514 .$
According to Equation (36), we find:
$A = − 0.8891 ,$
while according to Equation (33), we get:
$B = 5.2588 .$
If the obtained values of parameters A and B from (58) and (59) are put into Equation (27), we get:
$y + 3 = − 0.8891 x + 5.2588 .$
Since y = logv, x = logR, then the Equation (60) comes down to equation:
$log v + 3 = − 0.8891 log R + 5.2588 .$
Figure 5 graphically represents the values of logarithms of the rock mass oscillation velocities and the logarithms of the reduced distances, whereby the logarithms of the rock mass oscillation velocity are increased by number 3.
In diagram shown on Figure 5, the appropriate line is represented by Equation (61). From the Equation (61), we obtain:
$log v = − 0.8891 log R + 2.2588 .$
Since:
$2.2588 = log 9.5729 ,$
then Equation (63) can be written in the form:
$log v = l log 9.5729 R − 0.8891 ,$
where from we obtain:
$v = 9.5729 R − 0.8891 .$
Based on the obtained equations of rock mass oscillation (55), (56) and (65), it is possible to calculate values of rock mass oscillation velocities for corresponding reduced distances for models 1, 2 and 3.
In Table 7, there is presented the survey of reduced distances R, recorded oscillation velocities vr, calculated oscillation velocities vi1, vi2, vi3, as well as the difference between recorded and calculated rock mass oscillation velocities for models 1, 2 and 3.
Based on the data in Table 7, a statistical analysis was carried out and the following values were obtained:
• for model 1:
The curved line dependency index ρ1 between the reduced distance R and rock mass oscillation velocity v is:
ρ1 = 0.7786 (there is a strong correlation between R and v, given in the Equation (55)).
• for model 2:
ρ2 = 0.7831 (there is a strong correlation between R and v, given in the Equation (56)).
• for model 3:
ρ3 = 0.6937 (there is a significant correlation between R and v, given in the Equation (65)).

## 5. Conclusions

To establish the relationship between the oscillation velocity of the rock mass and basic parameters affecting its magnitude, being: the amount of explosive, the distance from the blast site, characteristics of the rock mass and the type of blasting, it is the equation of M. A. Sadovskii that is used most commonly.
Significant property of the equation of rock mass oscillation velocity depending on the reduced distance, obtained in this paper, is as follows: if the reduced distance R from any level is increased (decreased) by 1%, the oscillation velocity of the rock mass v will decrease (increase) by approximately n%.
In this paper, parameters n and K from Sadovskii’s equation were determined in three ways—models in the given work environment. In addition to the usual least squares method, two new models were applied to determine the parameters. In model two, the parameters of the rock mass oscillation equation were determined through the application of quotient of products of the same number of reduced distances and oscillation velocities. With model three, the parameters of the rock mass oscillation equation were determined using the linear connection between the logarithms of reduced distances and the logarithms of the oscillation velocity. In so doing, appropriate functions were obtained which showed the rock mass oscillation velocity depending on reduced distance.
Comparing values of the recorded oscillation velocities of the rock mass with the corresponding calculated ones, they are approximately then calculated corresponding indices of curvilinear correlation show that between the reduced distance and the rock mass oscillation velocity there is a significant and very strong curvilinear relationship expressed by the functions obtained, so we conclude that all three models can be used to calculate the rock mass oscillation velocity.

## Acknowledgments

This paper resulted from the Project No. 33025 and 33029 financed by the Ministry of Education of the Republic of Serbia realized by the Faculty of Mining and Geology in Belgrade.

## Author Contributions

Suzana Lutovac and Branko Gluščević performed the modelling, optimization and analysis of results. The paper was written by Suzana Lutovac with the support of Rade Tokalić and revised by Čedomir Beljić. Field measurement was performed by Jelena Majstorović.

## Conflicts of Interest

The authors declare no conflict of interest.

## References

1. Borla, O.; Lacidogna, G.; di Battista, E.; Niccolini, G.; Carpinteri, A. Electromagnetic Emission as Failure Precursor Phenomenon for Seismic Activity Monitoring, Fracture, Fatigue, Failure, and Damage Evolution; Springer: Cham, Switzerland; Heidelberg, Germany; New York, NY, USA; Dordrecht, The Netherlands; London, UK, 2015; Volume 5, pp. 221–229. [Google Scholar]
2. Kahriman, A. Analysis of Parameters of Ground Vibration Produced from Bench Blasting at a Limestone Quarry, Soil Dynamics and Earthquake Engineering 24; Elsevier: Amsterdam, The Netherlands, 2004; pp. 887–892. [Google Scholar]
3. Sambuelli, L. Theoretical Derivation of a Peak Particle Velocity–Distance Law for the Prediction of Vibrations from Blasting, Rock Mechanics and Rock Engineering; Springer: New York, NY, USA; London, UK, 2009; Volume 42, pp. 547–556. [Google Scholar]
4. Savić, M. Blasting at of the Open Pits, Monograph; Cooper Institute Bor: Bor, Serbia, 2000; pp. 317–319. [Google Scholar]
5. Trajković, S.; Slimak, Š.; Lutovac, S. Technique of Blasting and Shock Waves; Faculty of Mining and Geology: Belgrade, Serbia, 2005; p. 199. [Google Scholar]
6. Savić, L.; Ristović, I.; Trajković, S.; Savić, L. 2010: Experimental Research and Defining of the Optimum Bit Blade Sharpening Angle for Drilling Rocks and Ores; Technics Technologies Education Management (TTEM): Sarajevo, Bosnia and Herzegovina, 2010; Volume 5, pp. 205–217. [Google Scholar]
7. Medvedev, S.V. Seismics of Mountainous Explosions; Nedra: Moscow, Russia, 1964; pp. 42–43. [Google Scholar]
8. Lutovac, S. Model for Determining the Parameters of the Law of Oscillation of the Rock Mass in Blasting. Ph.D. Thesis, Faculty of Mining and Geology, Belgrade, Serbia, 2010; pp. 35–38. [Google Scholar]
9. Simeunović, D. Mathematics; Faculty of Mining and Geology, Mining Department: Belgrade, Serbia, 1985; pp. 101–103. [Google Scholar]
10. Lutovac, S.; Trajković, S.; Katona, O.; Savić, L.; Leković, B. Parameter Determination of Soil Oscillation in Limestone; Technics Technologies Education Management; DRUNPP: Sarajevo, Bosnia and Herzegovina, 2012; Volume 7, pp. 376–383. [Google Scholar]
11. Rakić, A. Seismic of Blasting, the Society of Engineers and Technicians; NIS Naftagas: Novi Sad, Serbia, 2005; pp. 54–60. [Google Scholar]
12. Kudryavtsev, V.A.; Demidovich, B.P. A Brief Course of Higher Mathematics; MIR Publishers: Moscow, Russia, 1984; pp. 293–295. [Google Scholar]
13. Lutovac, S.; Medenica, D.; Gluščević, B.; Tokalić, R.; Beljić, Č. Some models for determination of parameters of soil oscillation law during blasting operations. Energies 2016, 9, 617. [Google Scholar] [CrossRef]
14. Petrović, Ž. Business Statistics, Megatrend; University of Applied Sciences: Belgrade, Serbia; Faculty for Management: Zaječar, Serbia, 2004; p. 182. [Google Scholar]
15. Vukadinović, S. Elements of Probability and Mathematical Statistics; Technique for Copying Bureau of Construction: Belgrade, Serbia, 1970; pp. 329–331. [Google Scholar]
16. Žižić, M.; Lovrić, M.; Pavličić, D. Methods of Statistical Analysis; Faculty of Economics in Belgrade: Belgrade, Serbia, 1996; pp. 141–150. [Google Scholar]
17. Trajković, S. Area of Application of Laws of Rock Mass Oscillation Caused by Blasting, with Reference to Mining Works. Ph.D. Thesis, Faculty of Mining and Geology, Belgrade, Serbia, 1993. [Google Scholar]
18. Trajković, S.; Tokalić, R.; Lutovac, S. Study of Seismic Impaction the Surrounding Building Structures as the Result of Blasting at Nepričava Open Pit—Lajkovac; Faculty of Mining and Geology: Belgrade, Serbia, 2004; pp. 1–47. [Google Scholar]
Figure 1. Recording of oscillation velocity of the ground intended for blasting No. I—measuring point No. 3 at PK Nepričava.
Figure 1. Recording of oscillation velocity of the ground intended for blasting No. I—measuring point No. 3 at PK Nepričava. Figure 2. Graphic representation of rock masses oscillation velocity in Majdanpek open pit.
Figure 2. Graphic representation of rock masses oscillation velocity in Majdanpek open pit. Figure 3. Graphic representation of values logR and logv + 3 for calculating the parameters in Majdanpek open pit
Figure 3. Graphic representation of values logR and logv + 3 for calculating the parameters in Majdanpek open pit Figure 4. Graphic representation of rock masses oscillation velocity in Nepričava open pit.
Figure 4. Graphic representation of rock masses oscillation velocity in Nepričava open pit. Figure 5. Graphic representation of values logR and logv + 3 for calculating the parameters in Nepričava open pit.
Figure 5. Graphic representation of values logR and logv + 3 for calculating the parameters in Nepričava open pit. Table 1. Survey of blasting parameters and measurement results in Majdanpek open pit.
Table 1. Survey of blasting parameters and measurement results in Majdanpek open pit.
No.Blastingr (m)Q (kg)Rvt (cm/s)vv (cm/s)vl (cm/s)vrez (cm/s)
1I5.00.110.77220.7000.9001.1501.6194
2I6.00.112.92660.3000.9251.1501.5060
3I7.00.115.08100.2550.5400.8501.0388
4I8.00.117.23550.2600.3500.6500.7827
5I9.00.119.38990.2450.3250.5500.6842
6I10.00.121.54430.1550.2400.4050.4956
7I11.00.123.69880.2050.2750.4400.5579
8I12.00.125.85320.1100.2000.2850.3651
9I13.00.128.00770.0850.1900.2300.3102
10I14.00.130.16210.0900.1300.2000.2550
Table 2. Representation of data for calculating the parameters in Majdanpek open pit according to Model 3.
Table 2. Representation of data for calculating the parameters in Majdanpek open pit according to Model 3.
Rixiviyiyi + 3hi = xi+1 − xi(yi + 3 + yi+1 + 3)hi
10.77222.3769691.61940.4820563.4820560.1823181.256447
12.92662.5592871.50600.4094573.4094570.1541490.993879
15.08102.7134361.03880.0380663.0380660.1335350.773576
17.23552.8469710.7827−0.2450062.7549940.1177810.633130
19.38992.9467520.6842−0.3795052.6204950.1053590.518209
21.54433.0701110.4956−0.7019862.2980140.0953130.449347
23.69883.1654240.5579−0.5835762.4164240.0870100.388613
25.85323.2524340.3651−0.0075841.9924160.0800450.305922
28.00773.3324790.3102−0.1705381.8294620.0741070.256630
30.16213.4065860.2550−0.3664921.633508--
Table 3. Survey of recorded and calculated rock mass oscillation velocities for models 1, 2 and 3.
Table 3. Survey of recorded and calculated rock mass oscillation velocities for models 1, 2 and 3.
No.Rvr (cm/s)vi1 (cm/s)vi2 (cm/s)vi3 (cm/s)vr − vi1vr − vi2vrvi3
110.77221.61941.88211.92001.8043−0.2627−0.3006−0.1849
212.92661.50601,34751.36861.32000.15850.13740.1860
315.08101.03881.01591.02791.01350.02290.01090.0253
417.23550.78270.79540.80220.8061−0.0127−0.0195−0.0234
519.38990.68420.64100.64460.65880.04320.03960.0254
621.54430.49560.52840.53010.5499−0.0328−0.0345−0.0543
723.69880.55790.44370.44410.46700.11420.11380.0909
825.85320.36510.37830.37780.4023−0.0132−0.0127−0.0372
928.00770.31020.32670.32570.3507−0.0165−0.0155−0.0405
1030.16210.25500.28520.28380.3089−0.0302−0.0288−0.0539
Table 4. Survey of blasting parameters in Nepričava open pit.
Table 4. Survey of blasting parameters in Nepričava open pit.
BlastingNbQuk (kg)Qi (kg)Luk (m)L (m)
I8185.02680.03.0–3.2
II7135.020.070.03.0–3.2
III7105.019.061.03.0–3.2
Table 5. Survey of blasting parameters and measurement results in Nepričava open pit.
Table 5. Survey of blasting parameters and measurement results in Nepričava open pit.
No.Blas.MMr (m)Q (kg)Rvt (cm/s)vv (cm/s)vl (cm/s)vrez (cm/s)
1IMM–1256.90185.045.08580.1170.1360.1100.2104
2IMM–2220.81185.038.75210.1990.2960.2270.4228
3IMM–3224.58185.039.41370.1740.3130.1600,3922
4IMM–4161.78185.028.39230.2670.4170.2370.5490
5IMM–5179.10185.031.43200.2170.2110.1240.3271
6IIMM–1184.41135.035.94790.1950.1710.1100.2817
7IIMM–2151.52135.029.53650.3410.2760.2530.5064
8IIMM–3152.79135.029.78400.3720.4660.2210.6359
9IIMM–487.57135.017.07041.1700.8210.6151.5560
10IIMM–599.91135.019.47590.3700.4800.3040.6780
11IIIMM–1137.43105.029.13080.2850.2540.1800.4221
12IIIMM–2137.33105.029.10960.3000.3760.2460.5403
13IIIMM–372.09105.015.28081.0300.5650.3391.2227
14IIIMM–477.99105.016.53140.3370.2560.3150.5276
Table 6. Representation of data for calculating the parameters in Nepričava open pit according to Model 3.
Table 6. Representation of data for calculating the parameters in Nepričava open pit according to Model 3.
Rixiviyiyi + 3hi = xi+1 − xi(yi + 3 + yi+1 + 3)hi
15.28082.7265971.22270.2010623.2010620.0786650.437507
16.53142.8052620.5276−0.6394172.3605830.0320840.186174
17.07042.8373461.55600.4421183.4421180.1318320.798046
19.47592.9691780.6780−0.3886082.6113920.3769401.889123
28.39233.3461180.5490−0.5996572.4003430.0249500.119379
29.10963.3710680.5403−0.6156312.3843690.0007280.003292
29.13083.3717960.4221−0.8625132.1374870.0138310.061646
29.53653.3856270.5064−0.6804282.3195720.0083460.040609
29.78403.3939710.6359−0.4527142.5472860.0538550.238567
31.43203.4478260.3271−1.1174891.8825110.1342450.485376
35.94793.5820710.2817−1.2669151.7330870.0751140.290859
38.75213.6571850.4228−0.8608562.1391440.0169280.071151
39.41373.6741130.3922−0.9359832.0640170.1344540.471298
45.08583.8085670.2104−1.5587451.441255--
Table 7. Survey of recorded and calculated rock mass oscillation velocities for models 1, 2 and 3.
Table 7. Survey of recorded and calculated rock mass oscillation velocities for models 1, 2 and 3.
No.Rvr (cm/s)vi1 (cm/s)vi2 (cm/s)vi3 (cm/s)vr − vi1vr − vi2vr − vi3
115.28081.22271.07841.19340.84770.14430.02930.3750
216.53140.52760.97621.07120.7904−0.4486−0.5436−0.2628
317.07041.55600.93731.02510.76820.61870.53090.7878
419.47590.67800.79320.85530.6832−0.1152−0.1773−0.0052
528.39230.54900.49210.50980.48870.05690.03920.0603
629.10960.54030.47680.49260.47790.06350.04770.0624
729.13080.42210.47640.49210.4776−0.0543−0.0700−0.0555
829.53650.50640.46810.48290.47180.03830.02350.0346
929.78400.63590.46320.47740.46830.17270.15850.1676
1031.43200.32710.43270.44330.4464−0.1056−0.1162−0.1193
1135.94790.28170.36500.36870.3962−0.0833−0.0870−0.1145
1238.75210.42280.33190.33260.37060.09090.09020.0522
1339.41370.39220.32490.32490.36510.06730.06730.0271
1445.08580.21040.27400.27020.3239−0.0636−0.0598−0.1135
{"ft_lang_label":"__label__en","ft_lang_prob":0.8686849,"math_prob":0.9847587,"size":38764,"snap":"2022-40-2023-06","text_gpt3_token_len":10419,"char_repetition_ratio":0.19011869,"word_repetition_ratio":0.2186438,"special_character_ratio":0.28905687,"punctuation_ratio":0.17558467,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.99745345,"pos_list":[0,1,2,3,4,5,6,7,8,9,10,11,12],"im_url_duplicate_count":[null,null,null,2,null,2,null,2,null,2,null,2,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2022-09-29T08:30:02Z\",\"WARC-Record-ID\":\"<urn:uuid:d4d03bb3-fd82-4a5d-8c38-1ecfadb21379>\",\"Content-Length\":\"323807\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:69e03a4a-b31a-46a0-984d-ed643f4d44f3>\",\"WARC-Concurrent-To\":\"<urn:uuid:a0d239fa-7969-4090-9131-3d783a82c2f3>\",\"WARC-IP-Address\":\"104.18.24.151\",\"WARC-Target-URI\":\"https://www.mdpi.com/2075-163X/8/2/70/htm\",\"WARC-Payload-Digest\":\"sha1:NZ65X4UDWIBHBGKIQ673HZCXUWH5PGCU\",\"WARC-Block-Digest\":\"sha1:LYA522LX3LYUE24HMJLLDA2FVGQPKBS3\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2022/CC-MAIN-2022-40/CC-MAIN-2022-40_segments_1664030335326.48_warc_CC-MAIN-20220929065206-20220929095206-00184.warc.gz\"}"}