Javascript is required
1.
M. Sanei, L. Faramarzi, A. Fahimifar, S. Goli, A. Mehinrad, and A. Rahmati, “Shear strength of discontinuities in sedimentary rock masses based on direct shear tests,” Int. J. Rock Mech. Min. Sci., vol. 75, pp. 119–131, 2015. [Google Scholar] [Crossref]
2.
K. Karaman, A. Kaya, and A. Kesimal, “Use of the point load index in estimation of the strength rating for the RMR system,” J. Afr. Earth Sci., vol. 106, pp. 40–49, 2015. [Google Scholar] [Crossref]
3.
M. C. He, T. H. Xue, and Y. F. Peng, “Study on determination method of mechanical parameters of engineering rock mass,” Chin. J. Rock Mech. Eng., vol. 2, pp. 225–229, 2001. [Google Scholar] [Crossref]
4.
D. M. Ivars, M. E. Pierce, C. Darcel, J. Reyes-Montes, D. O. Potyondy, R. P. Young, and P. A. Cundall, “The synthetic rock mass approach for jointed rock mass modelling,” Int. J. Rock Mech. Min. Sci., vol. 48, no. 2, pp. 219–244, 2011. [Google Scholar] [Crossref]
5.
B. S. Chun, Y. J. Lee, D. D. Seo, and B. S. D. D. Lim, “Correlation deformation modulus by PMT with RMR and rock mass condition,” Tunn. Undergr. Space Technol., vol. 21, no. 3, 2005. [Google Scholar]
6.
L. Zhang, “Evaluation of rock mass deformability using empirical methods–A review,” Undergr. Space, vol. 2, no. 1, pp. 1–15, 2017. [Google Scholar] [Crossref]
7.
L. X. Huang, “Achievements and trends of rock dynamics research,” Rock Soil Mech., vol. 32, no. 10, pp. 2889–2900, 2011. [Google Scholar] [Crossref]
8.
K. Dragomiretskiy and D. Zosso, “Variational mode decomposition,” IEEE Trans. Signal Process., vol. 62, no. 3, pp. 531–544, 2014. [Google Scholar] [Crossref]
9.
P. Y. Li, W. X. Gao, X. J. Zhang, M. L. He, and L. Wang, “Prediction of peak velocity of blasting vibration based on SSA-BP,” BLASTING, pp. 1–12, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20231017.1747.012.html [Google Scholar]
10.
J. Tian, Y. P. Zhang, P. Yan, W. C. Sun, and X. Yang, “Noise reduction smooth model of open-pit deep hole blasting vibration signal based on CEEMD,” BLASTING, pp. 1–14, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20230814.1056.012.html [Google Scholar]
11.
K. Sarkar, V. Vishal, and T. N. Singh, “An empirical correlation of index geomechanical parameters with the compressional wave velocity,” Geotech. Geol. Eng., vol. 30, no. 2, pp. 469–479, 2012. [Google Scholar]
12.
Z. Huang, Q. Gu, Y. Wu, S. J. Li, K. Zhao, and R. Zhang, “Effects of confining pressure on acoustic emission and failure characteristics of sandstone,” Int. J. Min. Sci. Technol., vol. 31, no. 5, pp. 963–974, 2021. [Google Scholar]
13.
D. Xue, Z. Zhang, C. Chen, J. Zhou, L. Lu, X. T. Sun, and Y. T. Liu, “Spatial correlation-based characterization of acoustic emission signal-cloud in a granite sample by a cube clustering approach,” Int. J. Min. Sci. Technol., vol. 31, no. 4, pp. 535–551, 2021. [Google Scholar] [Crossref]
14.
A. I. Lawal, S. Kwon, O. S. Hammed, and M. A. Idris, “Blast-induced ground vibration prediction in granite quarries: An application of gene expression programming, ANFIS, and sine cosine algorithm optimized ANN,” Int. J. Min. Sci. Technol., vol. 31, no. 2, pp. 265–277, 2021. [Google Scholar]
15.
H. M. Zhou, S. C. Zhao, H. Z. Wang, H. Yu, W. H. Li, and X. T. Zhang, “Noise reduction analysis of loose blasting based on LMD-MFE-SVD,” BLASTING, pp. 1–9, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20230424.1538.002.html [Google Scholar]
16.
S. Wang, W. Q. Zhao, S. C. Zhao, G. H. Hao, F. Z. Long, and H. Su, “Noise reduction analysis of tunnel blasting vibration data based on OVMD-MPE algorithm,” BLASTING, pp. 1–9, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20230407.1333.006.html [Google Scholar]
17.
M. Sun, L. Wu, and J. Yang, “Application of time-frequency analysis of blasting vibration in underground caverns based on CEEMDAN-INHT,” BLASTING, pp. 1–9, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20221115.1923.028.html [Google Scholar]
18.
B. Sun, L. Song, S. Zeng, Y. Yi, and Z. Luo, “The open-pit mine blasting vibration signal decomposition method based on FSWT-FDM,” Metal Mine, pp. 1–12, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/34.1055.TD.20230419.1536.002.html [Google Scholar]
19.
L. Liu, S. Zhong, W. Yi, Y. Yang, and Y. Chai, “Study on spectral characteristics of blasting vibration signals at different elevations of open-pit slope based on PEMD,” Metal Mine, vol. 1, pp. 148–153, 2022. [Google Scholar] [Crossref]
20.
L. Qiu and C. Li, “STA/LTA method for picking up the first arrival of natural seismic waves and its improvement analysis,” Prog. Geophys., vol. 38, no. 4, pp. 1497–1506, 2023. [Google Scholar] [Crossref]
21.
T. Diehl, N. Deichmann, E. Kissling, and S. Husen, “Automatic S-wave picker for local earthquake tomography,” Bull. Seismol. Soc. Am., vol. 99, no. 3, pp. 1906–1920, 2009. [Google Scholar] [Crossref]
22.
C. Baillard, W. C. Crawford, V. Ballu, C. Hibert, and A. Mangeney, “An automatic kurtosis-based P-and S-phase picker designed for local seismic networks,” Bull. Seismol. Soc. Am., vol. 1, pp. 394–409, 104AD. [Google Scholar] [Crossref]
23.
J. M. O’Toole, A. Temko, and N. Stevenson, “Assessing instantaneous energy in the EEG: A non-negative, frequency-weighted energy operator,” in 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Chicago, IL, USA, 2014. [Google Scholar] [Crossref]
Search
Open Access
Research article

Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals

yiran yan1,
aobo liu1,
junpeng gai2,
zhenyang xu1*
1
School of Mining Engineering, University of Science and Technology Liaoning, Anshan 114051, Liaoning, China
2
Ansteel Group Corporation Limited, Anshan 114002, Liaoning, China
Acadlore Transactions on Geosciences
|
Volume 2, Issue 4, 2023
|
Pages 215-229
Received: 11-14-2023,
Revised: 12-14-2023,
Accepted: 12-19-2023,
Available online: 12-24-2023
View Full Article|Download PDF

Abstract:

The precision of determining rock mass mechanical parameters is notably impacted by mining blast activities. An advanced method for inverse analysis of these parameters, predicated upon measured blasting vibrations, has been developed. This approach employs a meticulous recognition of initial P-wave and S-wave arrivals within the vibrational energy spectrum. Utilizing principles from elastic wave theory, a novel framework has been established, correlating P-wave and S-wave velocities with dynamic characteristics of rock masses. The efficacy of this method has been substantiated through practical implementation, particularly in the Guanbaoshan Open-pit Iron Mine, Liaoning Province. Here, the derived density ratios were observed to range between 0.98 and 1.01, aligning closely with figures provided by authoritative research institutes. Additionally, the dynamic-to-static Poisson's ratio exhibited variations from 0.85 to 1.03, while the modulus of elasticity ratio dynamically to statically spanned from 2.0 to 2.6. These results, falling within anticipated theoretical ranges, underscore the robust applicability and accuracy of this method. The research contributes significantly to the domain of mining operations, particularly in optimizing blasting processes and enhancing the precision of mechanical parameter acquisition. It presents a pioneering approach, essential for addressing similar challenges in the mining sector.

Keywords: Rock mass mechanical parameters, Blasting vibrations, Initial arrival recognition, Rock mass wave velocity, Elastic wave theory

1. Introduction

The accurate assessment of rock mass mechanical parameters is quite critical for the effective planning and implementation of blasting operations. Inverse analysis methodologies have emerged as a pivotal tool in the computation of dynamic parameters of rock masses in blasting zones, contributing significantly to the refinement of blasting techniques, cost efficiency, and precision in execution.

Presently, the predominant methods for inverse analysis of rock mass mechanical parameters encompass experimental approaches, application of empirical formulas, and ultrasonic testing techniques. Experimental approaches are bifurcated into laboratory tests [1], [2], where rock samples from project sites are examined to ascertain mechanical parameters such as tensile and compressive strength through static experiments. However, this method may not adequately represent the rock structure in blasting areas. Conversely, in-situ tests aim to surmount the limitations of laboratory tests by objectively mirroring the mechanical characteristics of rocks in situ [3], albeit at a higher cost and complexity.

In the realm of theoretical and empirical research, scholars globally have made strides in estimating mechanical parameters of rock masses in blasting zones. Ivars et al. [4] introduced the synthetic rock mass (SRM) technology, facilitating the simulation of jointed rock masses and enabling the creation of models that encapsulate the spatial distribution of rock mass joints for numerical calculation of mechanical parameters. Chun et al. [5] leveraged the rock mass rating (RMR) system to analyze field data for deriving the mechanical parameters of the site-specific rock mass. Furthermore, Zhang [6] formulated equations that establish a correlation between rock quality designation (RQD) and rock mass mechanical parameters, drawing on extensive field data.

Ultrasonic testing methods have been identified as providing rapid inverse calculation of mechanical parameters in blasting areas, with advanced processing of blasting vibration waveforms enhancing signal fidelity. It was found by Huang [7] that acoustic waves, during their propagation, carry substantial information, the alteration of which can reflect the physical and mechanical properties of rock masses. The variational mode decomposition (VMD) algorithm, introduced by Dragomiretskiy and Zosso [8], achieves adaptive signal frequency domain separation, effectively mitigating mode mixing and splitting phenomena, and has gained widespread application in signal decomposition. In pursuit of enhanced prediction accuracy for peak blasting vibration velocity, Li et al. [9] devised a sparrow search algorithm (SSA) optimized back propagation (BP) neural network model, addressing the limitations of traditional BP neural network models in avoiding local optima. Similarly, Tian et al. [10] proposed a signal denoising and smoothing model based on complementary ensemble empirical mode decomposition (CEEMD), establishing an objective function predicated on filtering algorithm similarity and smoothness for optimal solution calculation.

Sarkar et al. [11] conducted measurements of various rock properties, such as density, slake durability index, uniaxial compressive strength (UCS), and P-wave velocity. These measurements underpinned the establishment of a correlation between P-wave velocity and rock mechanical properties, thus facilitating the determination of rock mechanical parameters via ultrasonic testing. Huang et al. [12] explored acoustic emission (AE) effects and failure characteristics of sandstone under varying confining pressures through uniaxial and triaxial compression acoustic emission tests. Their findings indicated a decrease in fractal dimension with increased confining pressure, suggesting a higher orderliness in rock failure. Xue et al. [13] employed a cube-cluster model to effectively characterize the spatial correlation of scattered AE events in granite samples during the rock failure process under triaxial compression conditions.

Further advancements were made by Lawal et al. [14], who proposed predictive models utilizing gene expression programming (GEP), adaptive neuro-fuzzy inference system (ANFIS), and sine cosine algorithm optimized artificial neural network (SCA-ANN) for ground vibrations from blasting in granite quarries. Additionally, Zhou et al. [15] developed a hybrid denoising method combining local mean decomposition (LMD), multiscale fuzzy entropy (MFE), and singular value decomposition (SVD), significantly enhancing denoising efficiency.

Wang et al. [16] employed a novel signal smoothing denoising model, integrating optimal variational mode decomposition (OVMD) with multiscale permutation entropy (MPE). Through the analysis of both simulated and actual engineering signals, it was demonstrated that this model is proficient in accurately identifying signal frequency information. The primary two components were shown to effectively represent the significant elements of the superimposed signal, proving particularly adept for high-precision data sequence analysis, feature extraction, and denoising, with a notable 22.05\% improvement in denoising error over the ensemble empirical mode decomposition (EEMD)-MPE model. Sun et al. [17] introduced the complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) algorithm, which incorporates adaptive noise signals into empirical mode decomposition (EMD). This innovation was found to effectively counter the inherent mode mixing issue inherent in EMD, enhancing the time-frequency analysis precision of seismic wave signals resulting from blast-induced ground vibrations. Moreover, Sun et al. [18] proposed a frequency slice wavelet transform-Fourier decomposition method (FSWT-FDM) for blasting vibration signals. This method, which combines FDM with FSWT, reduces ineffective signal components during decomposition, offering new perspectives for analyzing vibration signals in open-pit mining blasting. Liu et al. [19], in their analysis of open-pit slope blasting vibration tests, utilized the principal empirical mode decomposition (PEMD) method to eliminate noise in vibration signals, effectively revealing the frequency spectrum characteristics of slope blasting vibration signals at various elevations.

In summation, the predominant focus of current research is on ultrasonic testing and empirical formulae for the inversion of rock mass mechanical parameters and the processing of blasting vibration signals. However, there remains a paucity of studies that integrate blasting vibration signals with the mechanical properties of mine rocks for parameter inversion. Most existing methods rely on sample collection for mechanical parameter testing, which may not comprehensively represent regional rock mass characteristics. This study, employing initial arrival recognition methods for P-waves and S-waves within vibration signal energy and integrating elastic wave theory, deduces expressions linking P-wave velocity and S-wave velocity to rock mass dynamic parameters. Consequently, it introduces an inverse analysis method for these parameters based on blasting vibration signals. Applied in the Liaoning Province's Guanbaoshan Open-pit Iron Mine, this method has been proven to yield rational and reliable rock mass dynamic parameters, significantly elevating blasting efficiency and reducing associated costs.

2. Initial arrival Recognition Methods

2.1 Recognition of P-Wave Initial Arrival

The identification of the initial arrival of P-waves in blasting vibration signals is a critical aspect of rock mass dynamic parameter analysis. Commonly implemented techniques, such as the short-term average/long-term average (STA/LTA) method, spectral analysis, and pattern recognition, primarily focus on specific signal aspects. Nonetheless, their efficacy is limited when applied to complex, non-stationary signals like those generated by blasting vibrations. A notable characteristic of these signals is that the energy parameter variations around the P-wave initial arrival are more distinct than amplitude fluctuations, thus providing a more reliable identification marker.

In response to these challenges, a novel approach for P-wave initial arrival recognition, based on the energy of vibration signals, has been developed. This method incorporates the processing of trend items in blasting vibration signals and introduces the envelope derivative operator (EDO). EDO effectively amalgamates frequency domain information with the time-varying envelope of the signal, demonstrating superior capacity in tracking instantaneous changes in the signal. This characteristic renders it particularly effective for accurately demodulating signals, even in environments with substantial noise, thereby making it apt for processing non-stationary blasting vibration signals. The computational steps for the EDO involve the following:

$S[x(t)]=\left|X(t)^2\right|=|x(t)+j H[x(t)]|^2$
(1)

In this equation, $H[x(t)]$ denotes the Hilbert transform of $x(t)$, focusing solely on the instantaneous amplitude without considering the instantaneous frequency. The EDO, akin to the Teager-Kaiser (TK) energy operator method, employs a derivative function as a weighted filter, leveraging the properties of the Fourier transform:

$\begin{gathered}\dot{x}(t) \xrightarrow{\text { Fourier }} j w X(w) \\ X(\omega) \stackrel{\text { Fourier }}{\longrightarrow} x(t)\end{gathered}$
(2)

The EDO denoted as $\Gamma x(t)$ for a continuous signal $x(t)$ is defined as:

$\Gamma[x(t)]=|\dot{x}(t)+\mathrm{j} H[x(t)]|^2=\dot{x}^2(t)+H[\dot{x}(t)]^2$
(3)

where, $\dot{x}(t)=\frac{dx(t)}{dt}$, and $H[x(t)]$ represents the Hilbert transform of $x(t)$. For discrete signals $x(m)$, the TK energy operator typically employs a forward difference approach, while the EDO utilizes the central difference method. The discrete formulation of the EDO is presented as follows:

$\begin{gathered}\Gamma_d[x(m)]=\frac{1}{4}\left[x^2(m+1)+x^2(m-1)+h^2(m+1)+h^2(m-1)\right] \\ +\frac{1}{2}[x(m+1)+x(m-1)+h(m+1)+h(m-1)]\end{gathered}$
(4)

where, $h(m)=H(x(m))h$ represents the $m$-th sampling point within the data window. Subgraph (a) of Figure 1 illustrates the demodulation results of both the TK energy operator and EDO applied to a given signal. Conversely, Subgraph (b) of Figure 1 displays a signal comprising two components, which can be mathematically represented as:

$x(m)=1.5 \cos (\pi m / 32)+3.5 \cos (\pi m / 8)$
(5)
(a)
(b)
Figure 1. Examples of two algorithms for signal demodulation (a) Signal demodulation; (b) Test signal

As depicted in Figure 1, both the TK operator and EDO exhibit proficient tracking capabilities for energy changes within the signal. However, at certain extremal points, the TK operator yields negative values due to its modulation term, leading to singular values in subsequent energy separation processes. In contrast, the EDO demonstrates advantageous non-negative characteristics. The total energy of the blasting vibration wave can be quantified as:

$E=\frac{1}{2} \Delta m \sum\left(x_i^2+y_i^2+z_i^2\right)$
(6)

where, $x_i$, $y_i$, and $z_i$ denote the amplitudes of the recorded blasting vibrations. The time history curves for the signal components $x_i$, $y_i$, and $z_i$, and the total energy, depicted in Figure 2, reveal a pronounced increase in energy coinciding with the initial arrival of the P-wave.

(a)
(b)
Figure 2. Blasting vibration time history curve and energy time history curve

Adaptations of certain principles from the STA/LTA algorithm [20] were made to accentuate changes in energy, thus aiding the identification of the P-wave's initial arrival. The STA/LTA method involves computing the energy of the vibration within designated long and short calculation windows to pinpoint the P-wave's onset. The formula for the STA/LTA calculation is as follows:

$R=\frac{\sum C F_S / S}{\sum C F_L / L}$
(7)

The chosen characteristic function, sensitive to amplitude and frequency changes, is defined as:

$C F=X(i)^2+[X(i)-X(i-1)]^2$
(8)

where, $X(i)$ denotes the value at moment $i$. To enhance accuracy, a weighting coefficient, based on Eq. (8), is introduced:

$w=S T A^{\prime}(i)$
(9)

This coefficient aptly mirrors the signal’s variations in the frequency domain. From Eqs. (7) to (9), it becomes evident that the characteristic function, essential for identifying the moment of P-wave initial arrival, is computed as:

$F=\frac{\sum C F_S / S}{\sum C F_P / L} * w$
(10)

The analysis reveals that upon the P-wave's signal arrival, there is an abrupt shift in energy, characterized by a sharp increase. Consequently, the moment corresponding to the maximum value of Eq. (10) is determined to be the time of the P-wave's initial arrival.

2.2 Recognition of S-Wave Initial Arrival

Following the precise identification of the P-wave's initial arrival, the recognition of S-waves, particularly within an engineering context, often encounters challenges due to the interference of P-wave tails and ambient noise, factors influenced by propagation distance and site conditions. This research advances the understanding of S-wave recognition by comparing the physical properties of P-waves and S-waves to ascertain the initial arrival time of S-waves in measured blasting vibration waves. A critical factor in analyzing P-waves and S-waves is the calculation window width. The relationship between this window width and the dominant frequency of P-waves is generally expressed as:

$l=\frac{1}{f_{\mathrm{P}} \Delta t}$
(11)

where, $l$ denotes the window width length, $f_p$ represents the dominant frequency of P-waves, and $\Delta t$ denotes the instrument's sampling time. The calculation of the angle between the polarization directions of S-waves and P-waves is determined by the following formula:

$M=\left(\begin{array}{ccc}\operatorname{cov}(x, x) & \operatorname{cov}(x, y) & \operatorname{cov}(x, z) \\ \operatorname{cov}(y, x) & \operatorname{cov}(y, y) & \operatorname{cov}(y, z) \\ \operatorname{cov}(z, x) & \operatorname{cov}(z, y) & \operatorname{cov}(z, z)\end{array}\right)$
(12)

where, $M$ symbolizes the covariance matrix of P-waves over the length $l$, with the angle $\alpha$ between its maximum eigenvalue and the polarization direction of P-waves as the deflection angle. The normalization of $\alpha$ is achieved through:

$F_1=\frac{\alpha}{\pi / 2}$
(13)

Furthermore, the formula for calculating the polarization degree of P-waves and S-waves is given as:

$F_2=\frac{\left(\lambda_1-\lambda_2\right)^2+\left(\lambda_2-\lambda_3\right)^2+\left(\lambda_1-\lambda_3\right)^2}{\left(\lambda_1+\lambda_2+\lambda_3\right)^2}$
(14)

where, $\lambda_1, \lambda_2$ and $\lambda_3$ represent the eigenvalues of $M$ corresponding to each respective window. The transformation of the original coordinate system of the vibration signal, specifically radial $(X)$, tangential $(Y)$, and vertical $(Z)$ directions, into the $LQT$ coordinate system is executed using the following formula:

$\left[\begin{array}{l}L \\ Q \\ T\end{array}\right]=\left[\begin{array}{lll}u_{11} & u_{12} & u_{13} \\ u_{21} & u_{22} & u_{23} \\ u_{31} & u_{32} & u_{33}\end{array}\right]\left[\begin{array}{l}X \\ Y \\ Z\end{array}\right]$
(15)

where, $u_{i j}\,\,(j=1,2,3)$ denotes the cosine of the angle between the $i$-th principal direction and the radial $(X)$, tangential $(Y)$, and vertical $(Z)$ directions. The ratio of the transverse energy of P-waves and S-waves to the total energy is calculated as per the subsequent formula:

$F_3=\frac{\sum_{i=1}^l\left(Q_i^2+T_i^2\right)}{\sum_{i=1}^l\left(L_i^2+Q_i^2+T_i^2\right)}$
(16)

The method for identifying the initial arrival of S-waves is expressed through:

$F_{\mathrm{w}}=\sqrt{Q^2+T^2} \cdot F_1^2 \cdot F_2^2 \cdot F_3^2$
(17)

where, $\sqrt{{{Q}^{2}}+{{T}^{2}}}$ represents the transverse energy. Analysis of prior research, alongside Eqs. (11) to (17), leads to the conclusion that the moment yielding the maximum value of Eq. (17) is designated as the initial arrival time of the S-wave. The normalization process for Eq. (17) is formulated as:

$N=\frac{F_{\mathrm{w}}}{\max \left(F_{\mathrm{w}}\right)}$
(18)

3. Relationship Between Rock Mass Wave Velocity and Mechanical Parameters

3.1 Recognition of P-Wave Initial Arrival

The propagation of blasting vibration waves through rock mass media imparts key insights into the dynamic parameters of the rock, such as vibration velocity, frequency, duration, and waveform. These parameters are intricately influenced by the rock's mechanical properties. Additionally, the propagation process itself carries vital physical information pertinent to these mechanical characteristics. An analysis of the measured blasting vibration signals further elucidates the rock mass's elastic modulus, Poisson's ratio, and density. According to elastic wave theory, the correlation between P-wave and S-wave velocities and the dynamic parameters of the rock mass is established as follows:

$\begin{gathered}V_P=\sqrt{\frac{E_d\left(1-\mu_d\right)}{\rho\left(1+\mu_d\right)\left(1-2 \mu_d\right)}} \\ V_S=\sqrt{\frac{E_d}{2 \rho\left(1+\mu_d\right)}}\end{gathered}$
(19)

This equation encompasses longitudinal wave velocity $V_P$ (km/s), transverse wave velocity $V_S$ (km/s), the dynamic elastic modulus $E_d$ (GPa) of the rock mass, dynamic Poisson's ratio ${\mu }_{\text{d}}$, and the density $\rho$ (g/cm³) of the rock mass.

Eq. (19) demonstrates that the dynamic parameters of the rock mass can be deduced from the propagation velocities of blasting vibration waves, as articulated in the subsequent formulas:

$\mu_{\mathrm{d}}=\frac{\left(\frac{V_{\mathrm{p}}}{V_{\mathrm{s}}}\right)^2-2}{2\left[\left(\frac{V_{\mathrm{p}}}{V_{\mathrm{s}}}\right)^2-1\right]}$
(20)
$E_{\mathrm{d}}=\frac{\rho V_{\mathrm{S}}^2\left(3 V_{\mathrm{P}}^2-4 V_{\mathrm{S}}^2\right)}{V_{\mathrm{P}}^2-V_{\mathrm{S}}^2}$
(21)

Rock mass density is recognized as a fundamental physical property. Extensive research by scholars globally [21], [22], [23] has delved into the correlation between density and P-wave velocity. It has been observed that different types of rock masses exhibit diverse functional relationships between density and P-wave velocity, attributable to their distinct properties. Commonly, Gardner's empirical formula, introduced in 1974, is employed in these assessments:

$\rho=1.741 V_{\mathrm{P}}^{0.25}$
(22)

The calculation of P-wave and S-wave velocities within the rock mass at various measuring points and the blast source is formulated as:

$\Delta t=\frac{l}{V_{\mathrm{S}}}-\frac{l}{V_{\mathrm{P}}}$
(23)

where, l signifies the distance (m) between the blast source and the measuring point, and $\Delta t$ represents the time difference $(ms)$ between the initial arrivals of P-waves and S-waves.

From Eqs. (19) to (23), it is inferred that in practical applications, the analysis of the waveform of the measured blasting vibration signal at the measuring points, combined with the recognition of the initial arrivals of P-waves and S-waves, and the computation of P-wave and S-wave velocities between the points, facilitates the determination of the rock mass's elastic modulus, Poisson's ratio, and density in the area between the points. The inversion process of these parameters is delineated in Figure 3.

Figure 3. Inversion flowchart
3.2 Simulation Validation

To authenticate the precision of the method employed for identifying the initial arrival times of P-waves and S-waves based on vibration signal energy, the LS-DYNA dynamic analysis program was utilized to simulate the propagation of blasting vibrations in a rock mass within a predetermined area. The rock mass was represented using 3D Solid 164 elements, featuring dimensions of 20m×10m×15m, a borehole depth of 4m, stemming of 1m, a borehole diameter of 1m, and a charge diameter of 0.8m. Vibration monitoring points were strategically positioned above the model, as depicted in Figure 4.

Figure 4. Numerical model

Figure 5 demonstrates the efficacy of P-wave and S-wave recognition in the blasting vibration waveform at a distance of 15m from the blast center.

(a)
(b)
Figure 5. Numerical model

The rock's physical and mechanical parameters, essential for the analysis, are summarized in Table 1. Utilizing Eq. (23), the P-wave velocity was computed to be 5113m/s, and the S-wave velocity was determined to be 3064m/s.

Table 1. Key parameters of our model

Elastic Modulus $E_d$(GPa)

Poisson's Ratio $(\boldsymbol{\mu})$

Density $\left(\mathrm{g} / \mathrm{cm}^3\right)$

60

0.22

2.62

The proposed method was applied to identify the initial arrivals of P-waves and S-waves. These results were then juxtaposed with the theoretical initial arrival times. Table 2 details the accurately calculated initial arrival times of P-waves and S-waves. Upon comparison with the theoretical initial arrival times, it was found that the differences in the identified initial arrival times of P-waves and S-waves were nominal, with errors remaining below 3%. This outcome validates the efficacy of the recognition method and its compatibility with the requirements of field engineering.

Table 2. Key parameters of our model

P-Wave Initial Arrival Time (ms)

S-Wave Initial Arrival Time (ms)

Theoretical time

0.22

2.62

Experimental time

2.93

4.89

Error

$1.7 \%$

$2.5 \%$

4. Case Study

4.1 Engineering Background and Monitoring Scheme

A field experiment was conducted at the 285m level in the Guanbaoshan project, Liaoning Province. Area A, the focus of the blasting operation, was a rock region comprising 81 blastholes. The holes were spaced 6.0 m apart, with a row spacing of 4.2 m. A delay time of 17 ms was implemented, with hole depths varying between 14.0 m and 15.6 m. The maximum charge per hole was set at 202 kg, resulting in a total designed charge of 15044 kg. The objective of the study was to deduce the challenges associated with blasting in the subsequent area B, through the analysis of vibration conditions in area A. Eight measurement points were strategically placed on-site, distanced 8-10 m apart, along the central line of blasting area B, as illustrated in Figure 6.

(a)
(b)
Figure 6. Test scheme

For the experiment, the Topu integrated iSV-420 and the TC-4850 blasting vibrometer, developed by the Chinese Academy of Sciences, were employed. These instruments are capable of measuring three-dimensional velocity, acceleration, and other vibration parameters, making them apt for the on-site measurement of engineering blasting vibrations. Figure 7 shows the layout of the measurement points, while Figure 8 displays the arrangement of the on-site vibrometers.

Figure 7. Measuring point layout
Figure 8. Layout diagram of on-site vibrometer

The signal measured at point 1\# during the experiment is presented in Figure 9. The blasting vibration waveform exhibited a baseline deviation, necessitating correction before further analysis. The adjusted blasting vibration signal is depicted in Figure 10.

Figure 9. Blasting vibration signal diagram
Figure 10. Corrected blasting vibration signal
4.2 Inversion of P-Waves and S-Waves

The methods outlined earlier for identifying the initial arrivals of P-waves and S-waves were applied to the vibration waveform data collected. Table 3 displays the identification results for the initial arrival times of P-waves and S-waves at each measuring point. Additionally, this paper presents the characteristic function curves for the recognition of P-waves and S-waves in the blasting vibration waveform at a representative measuring point, point 4\#, as illustrated in Figure 11.

Table 3. First arrival recognition of P-waves and S-waves based on blasting vibration waveform

Measuring Point No.

P-Wave Initial Arrival Time (ms)

S-Wave Initial Arrival Time (ms)

$1 \#$

194.68

199.26

$2 \#$

195.13

200.83

$3 \#$

195.23

202.53

$4 \#$

195.25

203.98

$5 \#$

196.13

202.38

$6 \#$

196.25

203.96

$7 \#$

196.28

209.78

$8 \#$

196.24

211.05

(a)
(b)
Figure 11. Recognition effect of blasting vibration waveform: P-waves and S-waves

Upon substituting the initial arrival times of P-waves and S-waves, as identified in Table 3, into Eq. (23), the average velocities of P-waves and S-waves between various measuring points in blasting area B were calculated. The results of these calculations are summarized in Table 4.

Table 4. Inversion results of wave velocities in blast area B

$\boldsymbol{v}_p(\mathbf{k m} / \mathbf{s})$

$\boldsymbol{v}_s(\mathbf{k m} / \mathbf{s})$

From blast area to measuring point 1#

5.152

3.177

From measuring point 1# to 2#

5.227

3.135

From measuring point 2# to 3#

5.017

3.116

From measuring point 3# to 4#

5.272

3.152

From measuring point 4# to 5#

5.442

3.254

From measuring point 5# to 6#

5.328

3.273

From measuring point 6# to 7#

4.993

2.973

From measuring point 7# to 8#

4.856

2.927

To corroborate the reliability and accuracy of the inversion method, ultrasonic testing was performed in blast zone B, with an ultrasonic gauge (Figure 12) measuring the velocities of P-waves and S-waves. The results of these ultrasonic measurements are summarized in Table 5.

Figure 12. Ultrasonic gauge
Table 5. Wave velocity test results

$\boldsymbol{v}_p(\mathbf{k m} / \mathbf{s})$

$\boldsymbol{v}_s(\mathbf{k m} / \mathbf{s})$

From blast area to measuring point 1#

5.337

3.289

From measuring point 1# to 2#

5.432

3.301

From measuring point 2# to 3#

5.288

3.323

From measuring point 3# to 4#

5.520

3.272

From measuring point 4# to 5#

5.639

3.375

From measuring point 5# to 6#

5.617

3.443

From measuring point 6# to 7#

5.300

3.128

From measuring point 7# to 8#

5.075

3.081

The P-wave velocity in blast area B varied from 4.856 km/s to 5.442 km/s, and the S-wave velocity ranged from 2.927 km/s to 3.273 km/s, as indicated in Table 4. Ultrasonic measurements (Table 5) showed P-wave velocities varying from 5.075 km/s to 5.639 km/s, and S-wave velocities from 3.081 km/s to 3.443 km/s. The average errors between the P-wave and S-wave velocities obtained from the measured blasting vibration signal waveforms and those measured by ultrasonics were 4.67% and 4.74%, respectively, as shown in Table 6. This level of accuracy confirms that the method is suitable for practical engineering applications.

Table 6. Key parameters of our model

Relative Error $ \boldsymbol{v}_p(km/s)$

Relative Error $v_s$(km/s)

From blast area to measuring point $1 \#$

$3.60 \%$

$3.52 \%$

From measuring point $1 \#$ to $2 \#$

$3.93 \%$

$5.03 \%$

From measuring point $2 \#$ to $3 \#$

$5.41 \%$

$6.65 \%$

From measuring point $3 \#$ to $4 \#$

$4.70 \%$

$3.80 \%$

From measuring point $4 \#$ to $5 \#$

$3.62 \%$

$3.72 \%$

From measuring point $5 \#$ to $6 \#$

$5.43 \%$

$5.21 \%$

From measuring point $6 \#$ to $7 \#$

$6.15 \%$

$4.80 \%$

From measuring point $7 \#$ to $8 \#$

$4.53 \%$

$5.23 \%$

Average error

$4.67 \%$

$4.74 \%$

4.3 Relative error between ultrasonic test results and inversion results

The velocities of P-waves and S-waves derived from Table 4 were incorporated into Eqs. (20) to (22) for the calculation of the rock mass's mechanical parameters in the experimental area: dynamic elastic modulus, dynamic Poisson's ratio, and density. Table 7 encapsulates these results. Furthermore, Table 8 presents the static parameters of the rock mass recommended by the Anshan Iron and Steel Group Mining Design Institute for blast area B.

Table 7. Inversion results of rock parameters based on blasting vibration signals

Measuring Point Location

Dynamic Elastic Modulus $E_d$ (GPa)

Dynamic Poisson's Ratio $\mu_d$

Density $\left(\mathrm{g} / \mathrm{cm}^3\right)$

From blast area to measuring point $1 \#$

63.142

0.193

2.621

From measuring point $1 \#$ to $2 \#$

63.045

0.219

2.631

From measuring point $2 \#$ to $3 \#$

59.975

0.215

2.604

From measuring point $3 \#$ to $4 \#$

64.012

0.222

2.636

From measuring point 4 to $5 \#$

68.759

0.221

2.657

From measuring point $5 \#$ to $6 \#$

67.794

0.225

2.644

From measuring point $6 \#$ to $7 \#$

56.341

0.197

2.601

From measuring point $7 \#$ to $8 \#$

53.759

0.186

2.583

Table 8. Recommended values of static parameters of rock mass

Rock Mass

Dynamic Elastic Modulus $E_d$ (GPa)

Dynamic Poisson's Ratio $ \mu_d$

Density $\left(\mathrm{g} / \mathrm{cm}^3\right)$

Blast area B

26

0.22

2.62

Figure 13. Comparison diagram of rock mass parameters

Figure 13 depicts a comparison of the rock mass parameters. The dynamic Poisson's ratios, obtained from the inversion of blasting vibration waveforms, varied from 0.186 to 0.226, compared to the institute's provided ratio of 0.22, resulting in a comparative ratio of 0.85 to 1.03. The density derived from the waveform inversion, relative to the density value provided by the Anshan Iron and Steel Group Mining Design Institute, exhibited a comparative ratio of 0.98 to 1.01. The dynamic elastic modulus of the rock mass, ascertained from waveform inversion, ranged from 53 to 69 GPa, in contrast to the institute's recommended static elastic modulus of 26 GPa, yielding a comparative ratio of 2.0 to 2.6. These inversion results fall within the anticipated dynamic-static modulus ratio range of 2.2 to 2.9, substantiating the viability of deducing rock mass mechanical parameters from blasting vibration signal waveforms.

5. Conclusion

This study has addressed the gap in studies focusing on the inversion of mechanical parameters of rock masses within blasting areas using field blasting vibration data. Previously, methods predominantly revolved around sample collection for mechanical parameter testing without an in-depth analysis of the rock mass properties in the measured area. This paper has investigated the disparities in frequency, deflection angle, polarization, and energy between P-waves and S-waves in blasting vibration signals. Expressions for recognizing the initial arrivals of P-waves and S-waves have been derived, and based on elastic wave theory, the dynamic interplay between P-wave velocity, S-wave velocity, and the dynamic parameters of the rock mass has been analyzed. This approach offers an improved representation of the variations in frequency, energy, and deflection angle around the initial arrival time of S-waves. The key findings and prospects of this study are as follows:

Firstly, the average relative error between the P-wave and S-wave velocities, as recognized from the blasting vibration signal waveforms and those measured by ultrasonics in the field, is 4.67\% and 4.74\% respectively. These values fall within the acceptable range for engineering applications, underscoring the method's practical feasibility.

Secondly, by incorporating elastic wave theory, a relationship has been established between P-wave and S-wave velocities from blasting vibrations and the dynamic parameters of the rock mass. The proposed method, premised on the recognition of initial arrivals of P-waves and S-waves, inverts the dynamic parameters of rock from blasting vibration signals. Engineering practice demonstrates that the inverted density values approximate those provided by the research institute, with a comparative ratio of 0.98 to 1.01; the dynamic-static Poisson's ratio comparative ratio is 0.85 to 1.03; and the dynamic-static modulus of elasticity comparative ratio is 2.0 to 2.6. These inversion results align with prior research, validating the reliability and applicability of the dynamic parameters of the rock mass obtained through this method.

Finally, the influence of surface waves on the accuracy of P-wave and S-wave recognition was not examined in this paper, potentially affecting the results. The signal pickup process may inadvertently capture noise, and classifying the initial motion of P-waves and S-waves could be subject to inaccuracies. Additionally, the assumption that the blasting area represents an isotropic continuous medium, without considering the presence of joints, fractures, and other conditions within the rock mass, may introduce errors in the study of inverting the dynamic parameters using P-wave and S-wave velocities.

Funding
This study was funded by the Xingliao Talent Program (XLYC2203173).
Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References
1.
M. Sanei, L. Faramarzi, A. Fahimifar, S. Goli, A. Mehinrad, and A. Rahmati, “Shear strength of discontinuities in sedimentary rock masses based on direct shear tests,” Int. J. Rock Mech. Min. Sci., vol. 75, pp. 119–131, 2015. [Google Scholar] [Crossref]
2.
K. Karaman, A. Kaya, and A. Kesimal, “Use of the point load index in estimation of the strength rating for the RMR system,” J. Afr. Earth Sci., vol. 106, pp. 40–49, 2015. [Google Scholar] [Crossref]
3.
M. C. He, T. H. Xue, and Y. F. Peng, “Study on determination method of mechanical parameters of engineering rock mass,” Chin. J. Rock Mech. Eng., vol. 2, pp. 225–229, 2001. [Google Scholar] [Crossref]
4.
D. M. Ivars, M. E. Pierce, C. Darcel, J. Reyes-Montes, D. O. Potyondy, R. P. Young, and P. A. Cundall, “The synthetic rock mass approach for jointed rock mass modelling,” Int. J. Rock Mech. Min. Sci., vol. 48, no. 2, pp. 219–244, 2011. [Google Scholar] [Crossref]
5.
B. S. Chun, Y. J. Lee, D. D. Seo, and B. S. D. D. Lim, “Correlation deformation modulus by PMT with RMR and rock mass condition,” Tunn. Undergr. Space Technol., vol. 21, no. 3, 2005. [Google Scholar]
6.
L. Zhang, “Evaluation of rock mass deformability using empirical methods–A review,” Undergr. Space, vol. 2, no. 1, pp. 1–15, 2017. [Google Scholar] [Crossref]
7.
L. X. Huang, “Achievements and trends of rock dynamics research,” Rock Soil Mech., vol. 32, no. 10, pp. 2889–2900, 2011. [Google Scholar] [Crossref]
8.
K. Dragomiretskiy and D. Zosso, “Variational mode decomposition,” IEEE Trans. Signal Process., vol. 62, no. 3, pp. 531–544, 2014. [Google Scholar] [Crossref]
9.
P. Y. Li, W. X. Gao, X. J. Zhang, M. L. He, and L. Wang, “Prediction of peak velocity of blasting vibration based on SSA-BP,” BLASTING, pp. 1–12, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20231017.1747.012.html [Google Scholar]
10.
J. Tian, Y. P. Zhang, P. Yan, W. C. Sun, and X. Yang, “Noise reduction smooth model of open-pit deep hole blasting vibration signal based on CEEMD,” BLASTING, pp. 1–14, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20230814.1056.012.html [Google Scholar]
11.
K. Sarkar, V. Vishal, and T. N. Singh, “An empirical correlation of index geomechanical parameters with the compressional wave velocity,” Geotech. Geol. Eng., vol. 30, no. 2, pp. 469–479, 2012. [Google Scholar]
12.
Z. Huang, Q. Gu, Y. Wu, S. J. Li, K. Zhao, and R. Zhang, “Effects of confining pressure on acoustic emission and failure characteristics of sandstone,” Int. J. Min. Sci. Technol., vol. 31, no. 5, pp. 963–974, 2021. [Google Scholar]
13.
D. Xue, Z. Zhang, C. Chen, J. Zhou, L. Lu, X. T. Sun, and Y. T. Liu, “Spatial correlation-based characterization of acoustic emission signal-cloud in a granite sample by a cube clustering approach,” Int. J. Min. Sci. Technol., vol. 31, no. 4, pp. 535–551, 2021. [Google Scholar] [Crossref]
14.
A. I. Lawal, S. Kwon, O. S. Hammed, and M. A. Idris, “Blast-induced ground vibration prediction in granite quarries: An application of gene expression programming, ANFIS, and sine cosine algorithm optimized ANN,” Int. J. Min. Sci. Technol., vol. 31, no. 2, pp. 265–277, 2021. [Google Scholar]
15.
H. M. Zhou, S. C. Zhao, H. Z. Wang, H. Yu, W. H. Li, and X. T. Zhang, “Noise reduction analysis of loose blasting based on LMD-MFE-SVD,” BLASTING, pp. 1–9, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20230424.1538.002.html [Google Scholar]
16.
S. Wang, W. Q. Zhao, S. C. Zhao, G. H. Hao, F. Z. Long, and H. Su, “Noise reduction analysis of tunnel blasting vibration data based on OVMD-MPE algorithm,” BLASTING, pp. 1–9, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20230407.1333.006.html [Google Scholar]
17.
M. Sun, L. Wu, and J. Yang, “Application of time-frequency analysis of blasting vibration in underground caverns based on CEEMDAN-INHT,” BLASTING, pp. 1–9, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/42.1164.TJ.20221115.1923.028.html [Google Scholar]
18.
B. Sun, L. Song, S. Zeng, Y. Yi, and Z. Luo, “The open-pit mine blasting vibration signal decomposition method based on FSWT-FDM,” Metal Mine, pp. 1–12, 2023, [Online]. Available: http://kns.cnki.net/kcms/detail/34.1055.TD.20230419.1536.002.html [Google Scholar]
19.
L. Liu, S. Zhong, W. Yi, Y. Yang, and Y. Chai, “Study on spectral characteristics of blasting vibration signals at different elevations of open-pit slope based on PEMD,” Metal Mine, vol. 1, pp. 148–153, 2022. [Google Scholar] [Crossref]
20.
L. Qiu and C. Li, “STA/LTA method for picking up the first arrival of natural seismic waves and its improvement analysis,” Prog. Geophys., vol. 38, no. 4, pp. 1497–1506, 2023. [Google Scholar] [Crossref]
21.
T. Diehl, N. Deichmann, E. Kissling, and S. Husen, “Automatic S-wave picker for local earthquake tomography,” Bull. Seismol. Soc. Am., vol. 99, no. 3, pp. 1906–1920, 2009. [Google Scholar] [Crossref]
22.
C. Baillard, W. C. Crawford, V. Ballu, C. Hibert, and A. Mangeney, “An automatic kurtosis-based P-and S-phase picker designed for local seismic networks,” Bull. Seismol. Soc. Am., vol. 1, pp. 394–409, 104AD. [Google Scholar] [Crossref]
23.
J. M. O’Toole, A. Temko, and N. Stevenson, “Assessing instantaneous energy in the EEG: A non-negative, frequency-weighted energy operator,” in 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society, Chicago, IL, USA, 2014. [Google Scholar] [Crossref]

Cite this:
APA Style
IEEE Style
BibTex Style
MLA Style
Chicago Style
GB-T-7714-2015
Yan, Y. R., Liu, A. B., Gai, J. P., & Xu, Z. Y. (2023). Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals. Acadlore Trans. Geosci., 2(4), 215-229. https://doi.org/10.56578/atg020404
Y. R. Yan, A. B. Liu, J. P. Gai, and Z. Y. Xu, "Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals," Acadlore Trans. Geosci., vol. 2, no. 4, pp. 215-229, 2023. https://doi.org/10.56578/atg020404
@research-article{Yan2023InverseAO,
title={Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals},
author={Yiran Yan and Aobo Liu and Junpeng Gai and Zhenyang Xu},
journal={Acadlore Transactions on Geosciences},
year={2023},
page={215-229},
doi={https://doi.org/10.56578/atg020404}
}
Yiran Yan, et al. "Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals." Acadlore Transactions on Geosciences, v 2, pp 215-229. doi: https://doi.org/10.56578/atg020404
Yiran Yan, Aobo Liu, Junpeng Gai and Zhenyang Xu. "Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals." Acadlore Transactions on Geosciences, 2, (2023): 215-229. doi: https://doi.org/10.56578/atg020404
Yan Y., Liu A., Gai J., et al. Inverse Analysis of Rock Mass Dynamic Parameters from Blasting Vibration Signals[J]. Acadlore Transactions on Geosciences, 2023, 2(4): 215-229. https://doi.org/10.56578/atg020404
cc
©2023 by the author(s). Published by Acadlore Publishing Services Limited, Hong Kong. This article is available for free download and can be reused and cited, provided that the original published version is credited, under the CC BY 4.0 license.