Next Article in Journal
A Dynamic Hybrid Cryptosystem Using Chaos and Diffie–Hellman Protocol: An Image Encryption Application
Previous Article in Journal
Performance Analysis of QoS-Oriented OFDMA Protocol Based on IEEE 802.11ax for Cognitive Radio Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction and Optimization of Blasting-Induced Ground Vibration in Open-Pit Mines Using Intelligent Algorithms

1
School of Resources and Safety Engineering, Central South University, Changsha 410083, China
2
Hongda Blasting Engineering Group Co., Ltd., Guangzhou 510623, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(12), 7166; https://doi.org/10.3390/app13127166
Submission received: 29 May 2023 / Revised: 7 June 2023 / Accepted: 10 June 2023 / Published: 15 June 2023

Abstract

:
Prediction and parameter optimization are effective methods for mine personnel to control blast-induced ground vibration. However, the challenge of effective prediction and optimization lies in the multi-factor and multi-effect nature of open-pit blasting. This study proposes a hybrid intelligent model to predict ground vibrations using a least-squares support vector machine (LSSVM) optimized by a particle swarm algorithm (PSO). Meanwhile, multi-objective particle swarm optimization (MOPSO) was used to optimize the blast design parameters by considering the vibration of particular areas and the bulk rate of blast fragmentation. To compare the prediction performance of PSO-LSSVM, a genetic-algorithm-optimized BP neural network (GA-BP), unoptimized LSSVM, and BP were used, by applying the same database. In addition, the root-mean-squared error (RMSE), the mean absolute error (MAE), and the correlation coefficient (r) were regarded as the evaluation indicators. Furthermore, the optimization results of the blasting parameters were obtained by quoting the established vibration prediction model and bulk rate proxy model in MOPSO and verified by field tests. The results indicated that the PSO-LSSVM model provided the highest efficiency in predicting vibrations with an RMSE of 1.954, MAE of 1.717, and r of 0.965. Furthermore, the blasting vibration can be controlled by using the two-objective optimization model to obtain the best blasting parameters. Consequently, this study can provide more specific recommendations for vibration hazard control.

1. Introduction

Blasting is the preferred technique for open-pit mining activities, primarily aimed at achieving rock fragmentation. However, the release of blast energy will have adverse effects, e.g., seismic waves, flying rocks, and dust [1]. Hence, pre-blast prediction and optimization are essential tools to control blasting effects. In addition, research on blasting effects is shifting from traditional empirical formulas to intelligent algorithmic predictions.
To predict and optimize blasting effectiveness, researchers have conducted many studies. The traditional prediction methods are based on empirical models. Zhu et al. [2] modified the Kuz–Ram empirical model to improve the accuracy of blasting block prediction. Matidza et al. [3] compared the performance of the Sadovsky formulation with other empirical models in predicting blast vibration. Recently, intelligent algorithms have gained widespread use in predicting blast-induced outcomes, including the predicting of adverse effects, e.g., backbreak [4], dust emissions [5], and vibration [6], the predicting of direct blasting results, e.g., block degree [7] and throwing distance [8]. Moreover, blasting impacts the operations of various production segments in open-pit mines [9]. For example, bulk ore produced by blasting implies additional loading, transportation, and crushing costs [10]. In addition, blast cracking and bulkiness affect crushing and grinding [11]. Therefore, various indirect effects can be predicted based on blasting, such as the prediction of shovel loading time, truck payload [12], production energy consumption [13], etc. Above all, the objects of prediction and optimization for blasting effectiveness are multifaceted [14]. As the research of intelligent algorithms in mining and natural resources is intensifying, the prediction and optimization of multiple objects by intelligent algorithms are receiving more attention [15,16,17].
The different contributions of explosive energy are essentially responsible for the diversity of blasting results, i.e., a smaller fraction of energy (about 20–30%) is used for rock fragmentation and removing, whereas most of the energy is wasted in various adverse effects [18]. Among the various adverse effects, blasting-induced ground vibrations are considered the most harmful [19], accounting for approximately 40% of the explosive energy [20]. For example, blasting vibrations affect crushers and explosives magazines in open-pit mines [21]. Therefore, when there are important buildings and unstable terrain around the blasting area, blasting-induced ground vibration should be strictly controlled [22]. Peak particle velocity (PPV) is generally considered a crucial indicator of blast-induced ground vibration intensity [23]. Over time, numerous empirical models have been developed and proven effective in predicting PPV [3,24]. However, a literature review found that empirical models are less accurate [25,26]. Therefore, adopting intelligent algorithms and incorporating more factors affecting PPV into prediction models improves prediction accuracy [27,28,29].
Scholars have paid significant attention in recent years to enhancing intelligent algorithms for accurate PPV prediction. For example, Faradonbeh and Monjezi [20] used the cuckoo optimization algorithm (COA) to optimize gene expression programming (GEP), and Zhang et al. [30] optimized XGBoost using the particle swarm algorithm (PSO). Therefore, the hybrid models predict PPV more effectively, e.g., by optimizing the parameter configuration of prediction models. In addition, preprocessing the input data using algorithms (e.g., clustering techniques and feature selection) can enable predictive models to exhibit a dominant performance in a comprehensive comparison [19,31,32]. Comparing the hybrid model with a single model, some scholars found a significant improvement in PPV prediction accuracy with the hybrid approach [33]. However, these situations are common in that the prediction results based on blasting parameters cannot directly give reference values, especially since the adjustment of blasting design parameters still needs to rely on the engineers’ experience.
To improve the prediction accuracy, metaheuristic algorithms are utilized in hybrid models for predicting blast vibration [20,30,33]. Additionally, metaheuristic algorithms are also widely applied in mine blasting optimization, e.g., the grasshopper optimization algorithm (GOA) is employed to identify the blasting solution that minimizes dust generation [34]. Meanwhile, multi-objective particle swarm optimization (MOPSO), developed from PSO, is capable of handling multiple metrics in mining applications [35].
Therefore, establishing a highly reliable prediction model and giving reference values of blasting design parameters is particularly important for engineers in controlling blasting vibration hazards and guaranteeing blasting effectiveness. This study proposes a hybrid model of PSO and least-squares support vector machine (PSO-LSSVM) to predict the PPV of a sand and gravel mine. Furthermore, the LSSVM, GA-BP, and BP were considered for comparison purposes. In addition, the optimization strategy of the blasting design parameters is proposed based on multi-objective particle swarm optimization, and a blast vibration control scheme for sensitive areas is given.

2. Project Overview

This study was conducted at a gravel mine on a sea island in Zhoushan, Zhejiang Province. The mine has an original high–steep slope (Area e), as shown in Figure 1. Being mined from the top of the slope means the steps are negative relative to the elevation of the detonation point, which proves that the propagation of blasting seismic waves has a negative elevation effect under such geological conditions [36]. The mine mines the ore to produce sand and gravel aggregate, with a production scale of 4.74 million m3/year. Notably, due to the topography, the ore-crushing plant (Area c) is adjacent to the slope and is a typical blasting-vibration-sensitive area.
In this mine, the ore is mainly crystalline bosom tuff, rhyolite porphyry, and sunken tuff, with a dense structure and hard rock type. Besides, the fault and joints are not developed in the mine area, and the rock body is generally intact. Therefore, the area with a single lithology (rhyolite porphyry) was used as the test area (Area b). The main physical and mechanical parameters of the rock body are shown in Table 1. Rock blasting was carried out by deep-hole step blasting with a continuous charge structure and a step height of 15 m. The blasting parameters for the mine rock are shown in Figure 2. The arrow (d) in Figure 2 shows the direction of blast throwing, i.e., back from the existing slope.

3. Data Collection

In order to establish the prediction model of PPV and optimize the blasting parameters, a statistical scheme of the PPV data, including the blast design parameters and negative elevation (H), was developed. The values of the maximum explosive charge capacity (W), monitoring distance (R), negative elevation (H), burden (B), spacing (S), subdrilling (H0), powder factor (PF), and PPV were collected and recorded. In this statistical scheme, W, R, and B have been used by other scholars as controllable parameters to predict PPV [37], while S and PF were also shown by Nguyen and Bui [38] to be related to blasting vibration. In addition, negative elevation steps can have a decreasing effect on vibration velocity [36]. Therefore, H and H0 (meaning the maximum depth at which explosives are placed) were counted and used as input variables. The monitoring equipment for PPV caused by blasting vibration was the TC-4850 blasting [39] vibration tester, and the monitoring points were arranged as shown in Figure 3. R and H were determined by considering laser ranging and calibration on the construction drawings. A total of 30 blasting events were investigated. Due to topographic constraints, PPV in high–steep slope and blasting steps were monitored separately for different periods of blasting. Five measurement points were arranged in front of each blasting event. In addition, less-effective recovery data were obtained about the high–steep slope because of the large number of falling rocks. Therefore, a total of 50 sets of PPV data were recorded. The characteristics of the recorded input and output data are shown in Table 2.

4. Methods

4.1. Particle Swarm Optimization Algorithm

PSO belongs to the intelligent swarm algorithm (SI), which is derived from the observation and simulation of insect populations in nature. The basic principle is that the particle dynamically adjusts the movement speed and position according to its optimal position and the globally known optimal position, and the selection of the optimal position is determined by the fitness function [40]. In dealing with the optimization of indicators with multiple parameters, PSO has a powerful function. The dynamic adjustment of the particles is achieved by two core formulas, i.e., the velocity update formula and the position update formula, as follows:
V i k + 1 = ε V i k + c 1 r a n d ( 1 ) ( L P i k P i k ) + c 2 r a n d ( 2 ) ( G P i k P i k )
P i k + 1 = P i k + V i k + 1
Among them, ε is the inertia weight coefficient, i.e., the weight of the updated particle velocity vector; c1 is the individual learning acceleration coefficient; c2 is the global learning acceleration coefficient; rand(1) and rand(2) are random functions, ranging from 0 to 1; LPik and GPik are the local optimal position and the global optimal position of the kth iteration of the ith particle, respectively.

4.2. Least-Squares Support Vector Machine

SVM is an important branch of machine learning algorithms, minimizing structural risk as a fundamental principle, i.e., achieving a balance between empirical error and confidence intervals to improve the generalization ability of the algorithm. LSSVM was improved by SVM [41], which converts the quadratic inequality constraint of the optimization function into the equation constraints of the linear equation system, and the training time of the algorithm is shortened. At the same time, LSSVM inherits the advantages of SVM in dealing with problems such as small samples, nonlinearity, and local minima. The main principles of LSSVM are as follows:
Given a training sample dataset D = { ( x 1 , y 1 ) , ( x 2 , y 2 ) , , ( x N , y N ) } , i = 1 to N, where x i R d , i.e., a total of N d dimensions for input vectors xi, yi, is the output indicator. A nonlinear function maps the low-dimensional samples to the high-dimensional feature space. The nonlinear function estimation modeling is given as
y i = ω φ ( x i ) + b
where ω is the weight vector and b is the bias term.
In order to evaluate and classify the sample set, the following equality constraint was set as an optimization function in LSSVM.
min   J ( ω , e ) = min ( 1 2 ω 2 + 1 2 γ i = 1 N e i 2 ) s . t . y i = ω φ ( x i ) + b + e i
Among them, γ is the regularization parameter, which is used to determine the trade-off between model complexity and accuracy; ei is the error vector, which represents the difference between the actual value of the output indicator and the predicted value.
Solve the above optimization function by constructing the Lagrangian function and finding partial derivatives:
L ( ω , b , e , ξ ) = 1 2 ω 2 + 1 2 γ i = 1 N e i 2 i = 1 N ξ i [ ω φ ( x i ) + b + e i y i ]
where ξ i is the Lagrange multiplier.
Take the partial derivative for the parameters ω , b , e i , and ξ i , and set the partial derivative to zero [42]; we have:
{ L w = 0 ω = i = 1 N ξ i ω φ ( x i ) L b = 0 i = 1 N ξ i = 0 L e = 0 ξ i = γ e i L ξ = 0 ω φ ( x i ) + b + e i y i = 0
By eliminating ω and ei, the four linear equations of Equation (6) are simplified as:
[ 0 E T E Ω + 1 γ E ] [ b ξ ] = [ 0 y ]
Among them, y = [ y i , , y N ] T , ξ = [ ξ i , , ξ N ] T , E = [ 1 , , 1 ] T , and Ω is an N × N symmetry matrix of kernel function, as follows:
Ω = K ( x i , x j ) = φ ( x i ) φ ( x j )
Equations (7) and (8) lead to the decision function of the LSSVM model:
y ( x ) = i = 1 N ξ i K ( x , x i ) + b
where the radial basis function (RBF) is widely used as the kernel function. Furthermore, the RBF kernel has also been widely used and shown to achieve the best prediction performance in LSSVM prediction model studies of open-pit blasting vibration [43,44].
K ( x i , x j ) = exp ( x x i 2 2 σ 2 )
where σ is the key parameter of the kernel function.

4.3. Optimization of LSSVM by Particle Swarm Algorithm

It can be seen from Equations (7) and (8) that γ and σ are the critical parameters in determining the weight of the error vector in the LSSVM model [45]. Thus, improving the ability to learn and generalize is vital when applying the LSSVM model to predict nonlinear data with a small sample. In order to improve the prediction accuracy, this paper introduced the PSO algorithm into the LSSVM. Figure 4 shows the optimization flow of the hyperparameters ( γ and σ ) in the LSSVM model using PSO. The root-mean-squared error (RMSE) was chosen as the fitness function for the optimization of the two hyperparameters γ and σ :
fit ( γ , σ ) = 1 n i = 1 n ( y i y ^ i ) 2
Among them, n is the number of samples; y i is the actual value of the sample; y ^ i   is the predicted value of the sample.
In the model’s training, the particle velocity was updated by Equation (1), and the particle position was updated by Equation (2). The individual and global fitness were updated by calculating and finding the minimum RMSE value in each step. The minimum RMSE was recorded in the iterative process to find the corresponding optimal hyperparameters. The time cost of the algorithm was appropriately improved, while the prediction accuracy of the algorithm was enhanced.

4.4. Multi-Objective PSO Optimization

The MOPSO is recognized as an efficient tool for optimizing the proxy model of study objects, which was proposed by Coello and Lechuga [46] based on the regular PSO. Furthermore, the particles are updated according to Equations (1) and (2). The conventional PSO can only set one objective function (fitness function) to obtain a unique optimal solution. In contrast, the MOPSO can handle multiple functions. Therefore, MPSO can obtain a set of optimal solutions, i.e., Pareto solutions [35,47]. MOPSO can quickly converge to the Pareto front for its practical searchability. Then, the optimal values (function value and variables) are obtained by extracting the Pareto solution [48].
PPV and the bulk rate (Br) were taken as the two optimized objects to control the ground vibration in a sensitive area and ensure the blasting fragmentation effect. Therefore, the developed PSO-LSSVM (PPV) is called in MOPSO as one of the proxy models. Furthermore, an empirical model, which combined the Kuznetsov model with the Rosin-Rammler model (Kuz–Ram model), was considered as the proxy model for Br. The Kuz–Ram model of the Br is described as Equation (12) [49]. Moreover, the process of two-objective optimization by MOPSO is shown in Figure 5. After initializing the population, PSO-LSSVM (PPV) and Kuz–Ram (Br) are invoked through function calls. Then, the non-dominated solution is determined by comparing the fitness values and updated to the repository set. Meanwhile, the local optimum and the global optimum are selected in this update. The particles move according to Equations (1) and (2) until the maximum number of iterations.
{ x ¯ = K P F 0.8 ( 11.4 + H 0 ) 1 6 ( 115 100 ) 19 30 β = ( 2.2 14 B 115 ) ( 1 0.3 B ) ( 0.5 + 0.5 S B ) 11 15 x e = x ¯ ( 0.693 ) 1 6 B r = e ( x 0 x e ) β
where x ¯ is the average block size, β is the uniformity factor, xe is the characteristic block size, and Br is the ratio of ores with a block size greater than x0.

5. Modeling for Prediction and Optimization

5.1. PSO-LSSVM Model Building

In this study, intelligent algorithms were built based on the MATLAB software. Moreover, the RMSE was used as a performance index for modeling. The performance of the prediction and RMSE were negatively correlated. The smaller the RMSE value, the better the prediction effect is. In order to obtain the best prediction performance of the PSO-LSSVM model, the maximum number of iterations (Imax) and the number of particles in the swarm (npop) of PSO were determined according to the experimental method, and the other parameters were taken empirically. In addition, the inertia weight coefficient ( ε ) was iterated according to formula ε ( ε 0.4 ) T / k . As shown in Figure 6, after 100 iterations, there was no significant change in the fitness curve. Among them, the fitness curve with npop = 20 had the smallest RMSE at the final steady state and reached stability at approximately 70 iterations. Therefore, Imax was taken as 100, and npop was taken as 20. The initial parameters of the model are shown in Table 3.
Therefore, the optimal values of γ and σ 2 were 17.315 and 0.917, respectively, when PPV was used as the output indicator. Based on the optimized hyperparameters and training set, the regression fit of the PSO-LSSVM model was obtained, as shown in Figure 7. It can be seen that the PSO-LSSVM model fit the PPV better.

5.2. Parameter Configuration of the Models Used for Comparison

In order to validate the predictive performance of the proposed joint model (PSO-LSSVM), other hybrid models (GA-BP) and single models (LSSVM and BP) were used for comparison purposes. BP neural networks were proposed by Rumelhar et al. [50]. Their structure consists of an input layer, a hidden layer, and an output layer. In this study, a three-layer neural network was used. BP has the disadvantage of high workload and reliance on empirical values when assigning the connection weights and thresholds between neurons [51]. GA belongs to the metaheuristic algorithm [33]. Using the algorithm function of GA (e.g., selection, crossover, and mutation) to optimize BP is considered adequate [52].
According to the method in [53], the LSSVM was trained three times for the training datasets with three different sets of γ and σ2, i.e., (50, 1.9), (10, 0.6), and (4, 0.04). The one with the smallest RMSE was taken to participate in the comparison. The research experience in [52] was introduced to show that the BP and GA-BP methods were applied in the prediction of PPV. The accuracy of the training model was set to 0.0001; the learning rate was 0.1; the maximum number of iterations was 1000. The weight threshold was set in the range of [−2, 2] for the BP neural network without optimization. The multiple predictions (10 times) of BP were performed, and the group with the smallest RMSE was taken to participate in the discussion. All the RMSEs of LSSVM and BP in the modeling process are shown in Figure 8.

5.3. Two-Objective Optimization Model

PSO-LSSVM (PPV) for the sensitive area (c) and Kuz–Ram (Br) in the chosen blasting area (a) were used as a proxy model in the MOPSO. After the site survey, the R of the sensitive area (c) and blasting area was 200 m, and H was 90 m. According to the actual ore, the bulk size x0 was 80 cm. Therefore, the two-objective optimization strategy required considering five variables, i.e., W, PF, S, B, and H0. Therefore, the number of decision variables (nVar) was five. The parameter set is shown in Table 4.
In the MOPSO, all the data were treated the same way as the PPV prediction model, i.e., the problem of different magnitudes was avoided by means of normalization and anti-normalization. The data were manipulated by Equation (13). For example, the random variable (X) assigned to the PSO-LSSVM model was within 0 and 1, while the predicted PPV was anti-normalized to the original magnitudes. Besides, the randomly generated variables were anti-normalized to the values with the range of [xmin, xmax] of measured variables and quoted by Equation (13) to obtain the bulk rate. Therefore, the optimization was conducted on 30 blasting experiences.
X = x x min x max x min

6. Results and Discussion

6.1. Comparison and Evaluation of Prediction Models

To verify the performance of the developed prediction model (PSO-LSSVM), the same training datasets and testing datasets were utilized. Furthermore, the root-mean-squared error (RMSE), mean absolute error (MAE), and correlation coefficient (r) were selected as the assessment indicators. The RMSE and MAE were negatively correlated with the performance of the models, while on the contrary for r. The efficiencies of PSO-LSSVM, GA-BP, the unoptimized LSSVM model, and BP were compared by calculating and evaluating the performance indicators. As shown in Table 5, it can be seen that the hybrid models (PSO-LSSVM and GA-BP) outperformed the single prediction models (LSSVM and BP). In the single model, LSSVM seemed to exhibit a slightly better predictive performance than BP, with smaller error metrics (RMSE and MAE) and larger fit coefficients (r). Importantly, Table 5 shows that PSO-LSSVM had a more reliable performance compared with the remaining models.
Additionally, a PPV profile between the predicted and the measured for training and testing is shown in Figure 9 and Figure 10. Based on the graphical representation of the regression fit in Figure 9 and Figure 10, the predicted values of PSO-LSSVM and GA-BP had a tendency to be close to the fit line. For instance, the predicted values of the single model had more discrete points, with a lower r. Consequently, the PSO-LSSVM model was closest to the measured values with an r of 0.983 and 0.965 for the training and test sets.

6.2. Optimization of Blast Design Parameters

In this section, a strategy of two-objective optimization is presented, i.e., using MOPSO to format the Pareto front by assigning and solving PSO-LSSVM (PPV) and Kur–Ren (Br). Therefore, a stable Pareto front was obtained after 100 iterations. As shown in Figure 11, the Pareto solution meant a set of optimized function values, which were limited to each other. Spontaneously, the corresponding independent variable could be extracted from the repository set. Furthermore, Points A and C were the two endpoints of the Pareto-optimal solution set, which corresponded to the two cases of the minimum PPV and Br, respectively. Point B was the middle point, which had a more balanced consideration between the two objectives.
Thus, the three classical solutions (A, B, and C) are shown in Table 6. It is worth noting that the transformation of PPV in Sensitive Region c was weaker than the Br of the blasting fragmentation when the parameters changed. In order to validate the optimization results, a blasting scheme was designed using the blasting parameters extracted from Point B. Then, limited blasting field tests were carried out, and a total of three set dates were recorded, finally. Therefore, considering the B solution in the blasting scheme, the measured PPV values for Sensitive Area c were 0.752 cm/s, 1.00 cm/s, and 1.01 cm/s, and the Brs of the blasting fragmentation bulk were 11.02%, 9.06%, and 8.42%. Consequently, the actual measured values were similar to the optimized values, indicating that the optimization results were convincing.

7. Conclusions

Vibration hazard control is the same as rock fragmentation and is the goal pursued in open-pit blasting. Meanwhile, reliable and accurate blasting vibration prediction models and referenceable blasting design parameters are important for the control of blasting-induced ground vibration. To accurately predict and control blasting vibrations, this paper proposed the PSO-LSSVM prediction model and the MOPSO optimization model. The following conclusions can be drawn:
  • The PSO-LSSVM model outperformed GA-BP, unoptimized LSSVM, and BP in terms of the prediction efficiency and accuracy of PPV, with better prediction error (RMSE, MAE) and goodness of fit (r). More specifically, the hyperparameter-optimized PSO-LSSVM reduced the prediction error (RMSE, MAE) by more than 55% and improved the goodness of fit (r) by more than 10%, compared with the least-efficient one (unoptimized BP).
  • With the multi-objective optimization function of PSO, a relatively optimal blasting design was found based on existing blasting tests. When considering the PPV of the sensitive area (R = 100 m, H = 90 m) and the Br, the blasting design parameters for Zhoushan mine were taken as W = 173.68 kg, PF = 0.38 kg/m3, S = 5.96 m, B = 3.19 m, and H0 = 1.91 m, corresponding to a PPV of 0.893 cm/s and Br of 9.8%. Therefore, it can provide engineers with a more specific and reliable blasting vibration control reference.
  • In order to predict and control the PPV more accurately, more experimental data and influencing factors should be considered. Furthermore, reliable agent models are a requirement for more accurate multi-objective optimization results.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/app13127166/s1.

Author Contributions

All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by P.L. The first draft of the manuscript was written by P.Z. and J.G. and all authors commented on previous versions of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Project supported by the Central South University School-Enterprise Joint Project “Analysis of Blasting Production Data of Dahuang Mountain Tuff Mine and full operation chain optimization study” (2022XQLH080); supported by the Central South University-Hongda Blasting Engineering Group Postgraduate Joint Training Base (2020pyjd91).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available in the Supplementary Materials.

Acknowledgments

The team of authors express their gratitude to the editors and reviewers for valuable recommendations that have been taken into account to improve significantly the quality of this article. We also thank Xiaobing Yan and Yi Liu for providing the field guidance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khandelwal, M.; Singh, T.N. Prediction of blast-induced ground vibration using artifificial neural network. Int. J. Rock Mech. Min. Sci. 2009, 46, 1214–1222. [Google Scholar] [CrossRef]
  2. Zhu, Z.; Yin, Y.X.; Chen, M.; Wei, D.; Lu, W.B.; Liu, J.C. Study on the prediction of blasting block size in Changjiushenshan chert mining. Blasting 2021, 38, 17–23+31. [Google Scholar]
  3. Matidza, M.I.; Jianhua, Z.; Gang, H.; Mwangi, A.D. Assessment of blast-induced ground vibration at jinduicheng molybdenum open pit mine. Nat. Resour. Res. 2020, 29, 831–841. [Google Scholar] [CrossRef]
  4. Khandelwal, M.; Monjezi, M. Prediction of backbreak in open-pit blasting operations using the machine learning method. Rock Mech. Rock Eng. 2013, 46, 389–396. [Google Scholar] [CrossRef]
  5. Bakhtavar, E.; Hosseini, S.; Hewage, K.; Sadiq, R. Green blasting policy: Simultaneous forecast of vertical and horizontal distribution of dust emissions using artificial causality-weighted neural network. J. Clean. Prod. 2021, 283, 124562. [Google Scholar] [CrossRef]
  6. Mostafa, T.M. Artificial neural network for prediction and control of blasting vibrations in Assiut (Egypt) limestone quarry. Int. J. Rock Mech. Min. Sci. 2009, 46, 426–431. [Google Scholar]
  7. Wang, R.C.; Wu, S. LM-based nerve network model for predicting burst degree. J. Hydropower Gener. 2019, 38, 100–109. [Google Scholar]
  8. Yu, Z.; Shi, X.; Zhou, J.; Chen, X.; Miao, X.; Teng, B.; Ipangelwa, T. Prediction of blast-induced rock movement during bench blasting: Use of gray wolf optimizer and support vector regression. Nat. Resour. Res. 2019, 29, 843–865. [Google Scholar] [CrossRef]
  9. Kinyua, E.M.; Jianhua, Z.; Kasomo, R.M.; Mauti, D.; Mwangangi, J. A review of the influence of blast fragmentation on downstream processing of metal ores. Miner. Eng. 2022, 186, 107743. [Google Scholar] [CrossRef]
  10. Bakhtavar, E.; Sadiq, R.; Hewage, K. Optimization of blasting-associated costs in surface mines using risk-based probabilistic integer programming and firefly algorithm. Nat. Resour. Res. 2021, 30, 4789–4806. [Google Scholar] [CrossRef]
  11. Nielsen, K.; Lownds, C.M. Enhancement of taconite crushing and grinding through primary blasting. Int. J. Rock Mech. Min. Sci. 1997, 34, 226.e1–226.e14. [Google Scholar] [CrossRef]
  12. Ohadi, B.; Sun, X.; Esmaieli, K.; Consens, M.P. Predicting blast-induced outcomes using random forest models of multi-year blasting data from an open pit mine. Bull. Eng. Geol. Environ. 2020, 79, 329–343. [Google Scholar] [CrossRef]
  13. Silwa, J.D.; Amaya, J.G.; Asso, F. Development of a predictive model of fragmentation using drilling and blasting data in open pit mining. J. S. Afr. Inst. Min. Metall. 2017, 117, 1087–1094. [Google Scholar]
  14. Bakhtavar, E.; Abdollahisharif, J.; Ahmadi, M. Reduction of the undesirable bench-blasting consequences with emphasis on ground vibration using a developed multi-objective stochastic programming. Int. J. Min. Reclam. Environ. 2017, 31, 333–345. [Google Scholar] [CrossRef]
  15. Jang, H.; Topal, E. A review of soft computing technology applications in several mining problems. Appl. Soft Comput. 2014, 22, 638–651. [Google Scholar] [CrossRef]
  16. Wang, L.; Yao, Y.; Wang, K.; Adenutsi, C.D.; Zhao, G.; Lai, F. A novel surrogate-assisted multi-objective optimization method for well control parameters based on tri-Training. Nat. Resour. Res. 2021, 30, 4825–4841. [Google Scholar] [CrossRef]
  17. Sacchelli, S.; Bernetti, I. Integrated management of forest ecosystem services: An optimization model based on multi-objective analysis and metaheuristic approach. Nat. Resour. Res. 2019, 28 (Suppl. 1), 5–14. [Google Scholar] [CrossRef]
  18. Armaghani, D.J.; Hajihassani, M.; Mohamad, E.T.; Marto, A.; Noorani, S. Blasting-induced flyrock and ground vibration prediction through an expert artifificial neural net work based on particle swarm optimization. Arab. J. Geosci. 2014, 7, 5383–5396. [Google Scholar] [CrossRef]
  19. Bui, X.; Choi, Y.; Atrushkevich, V.; Nguyen, H.; Tran, Q.; Long, N.Q.; Hoang, H. Prediction of blast-induced ground vibration intensity in open-pit mines using unmanned aerial vehicle and a novel intelligence system. Nat. Resour. Res. 2020, 29, 771–790. [Google Scholar] [CrossRef]
  20. Faradonbeh, R.; Monjezi, M. Prediction and minimization of blast-induced ground vibration using two robust meta-heuristic algorithms. Eng. Comput. 2017, 33, 835–851. [Google Scholar] [CrossRef]
  21. Bakhtavar, E.; Yousefi, S. Analysis of ground vibration risk on mine infrastructures: Integrating fuzzy slack-based measure model and failure effects analysis. Int. J. Environ. Sci. Technol. 2019, 16, 6065–6076. [Google Scholar] [CrossRef]
  22. Nguyen, H.; Drebenstedt, C.; Bui, X.N.; Bui, D.T. Prediction of blast-induced ground vibration in an open-pit mine by a novel hybrid model based on clustering and artifificial neural network. Nat. Resour. Res. 2019, 29, 691–709. [Google Scholar] [CrossRef]
  23. Hasanipanah, M.; Golzar, S.B.; Larki, I.A.; Maryaki, M.Y.; Ghahremanians, T. Estimation of blast-induced ground vibration through a soft computing framework. Eng. Comput. 2017, 33, 951–959. [Google Scholar] [CrossRef]
  24. Murmu, S.; Maheshwari, P.; Verma, H.K. Empirical and probabilistic analysis of blast-induced ground vibrations. Int. J. Rock Mech. Min. Sci. 2018, 103, 267–274. [Google Scholar] [CrossRef]
  25. Ghasemi, E.; Ataei, M.; Hashemolhosseini, H. Development of a fuzzy model for predicting ground vibration caused by rock blasting in surface mining. J. Vib. Control. 2013, 19, 755–770. [Google Scholar] [CrossRef]
  26. Ainalis, D.; Kaufmann, O.; Tshibangu, J.P.; Verlinden, O.; Kouroussi, G. Assessing blast source pressure modelling approaches for the numerical simulation of ground vibrations. In Proceedings of the 23rd International Congress on Sound and Vibration, Athens, Greece, 10–14 July 2016. [Google Scholar]
  27. Shirani Faradonbeh, R.; Jahed Armaghani, D.; Abd Majid, M.Z.; Md Tahir, M.; Ramesh Murlidhar, B.; Monjezi, M.; Wong, H.M. Prediction of ground vibration due to quarry blasting based on gene expression programming: A new model for peak particle velocity prediction. Int. J. Environ. Sci. 2016, 13, 1453–1464. [Google Scholar] [CrossRef] [Green Version]
  28. Zhou, J.; Shi, X.; Li, X. Utilizing gradient boosted machine for the prediction of damage to residential structures owing to blasting vibrations of open pit mining. J. Vib. Control. 2016, 22, 3986–3997. [Google Scholar] [CrossRef]
  29. Armaghani, D.J.; Momeni, E.; Abad, S.V.A.N.K.; Khan, D.M. Feasibility of ANFIS model for prediction of ground vibrations resulting from quarry blasting. Environ. Earth Sci. 2015, 74, 2845–2860. [Google Scholar] [CrossRef] [Green Version]
  30. Zhang, X.; Nguyen, H.; Bu, X.N.; Tran, Q.H.; Nguyen, D.A.; Bui, D.T.; Moayedi, H. Novel soft computing model for predicting blast-induced ground vibration in open-pit mines based on particle swarm optimization and XGBoost. Nat. Resour. Res. 2019, 29, 711–721. [Google Scholar] [CrossRef]
  31. Xue, X. Neuro-fuzzy based approach for prediction of blast-induced ground vibration. Appl. Acoust. 2019, 152, 73–78. [Google Scholar] [CrossRef]
  32. Zhou, J.; Asteris, P.G.; Armaghani, D.J.; Pham, B.T. Prediction of ground vibration induced by blasting operations through the use of the Bayesian Network and random forest models. Soil Dyn. Earthq. Eng. 2020, 139, 106390. [Google Scholar] [CrossRef]
  33. Yang, H.; Hasanipanah, M.; Tahir, M.M.; Bui, D.T. Intelligent Prediction of blasting-induced ground vibration using ANFIS optimized by GA and PSO. Nat. Resour. Res. 2020, 29, 739–750. [Google Scholar] [CrossRef]
  34. Hosseini, S.; Monjezi, M.; Bakhtavar, E. Minimization of blast-induced dust emission using gene-expression programming and grasshopper optimization algorithm: A smart mining solution based on blasting plan optimization. Clean Technol. Environ. Policy 2022, 24, 2313–2328. [Google Scholar] [CrossRef]
  35. Sadrossadat, E.; Basarir, H.; Karrech, A.; Durham, R.; Fourie, A.; Bin, H. The optimization of cemented hydraulic backfifill mixture design parameters for different strength conditions using artifificial intelligence algorithms. In Proceedings of the 28th International Symposium on Mine Planning and Equipment Selection—MPES 2019, Perth, WA, Australia, 2-4 December 2019; Topal, E., Ed.; Springer: Berlin/Heidelberg, Germany, 2019; pp. 219–227. [Google Scholar]
  36. Zhang, X.; Shi, J. Prediction of vibration velocity of bench blasting reflecting negative elevation effect. Geofluids 2021, 2021, 6662809. [Google Scholar] [CrossRef]
  37. Yu, Z.; Shi, X.Z.; Zhou, J.; Chen, X.; Qiu, X.Y. Effective assessment of blast-induced ground vibration using an optimized random forest model based on a harris hawks optimization algorithm. Appl. Sci. 2020, 10, 1403. [Google Scholar] [CrossRef] [Green Version]
  38. Nguyen, H.; Bui, X.N. Predicting blast-induced air overpressure: A robust artificial intelligence system based on artificial neural networks and random forest. Nat. Resour. Res. 2018, 28, 893–907. [Google Scholar] [CrossRef]
  39. Guo, J.; Zhang, C.; Xie, S.; Liu, Y. Research on the prediction model of blasting vibration velocity in the dahuangshan mine. Appl. Sci. 2022, 12, 5849. [Google Scholar] [CrossRef]
  40. Kennedy, J.; Eberhart, R. Particle swarm optimization. In Proceedings of the ICNN’95—International Conference on Neural Networks, Perth, Australia, 27 November–1 December 1995. [Google Scholar]
  41. Suykens, J.A.K.; Vandewale, J. Least squares support vector machine classifiers. Neural Process. Lett. 1999, 9, 293–300. [Google Scholar] [CrossRef]
  42. Wu, Y.Y.; Li, N.N.; Xue, W.; Xu, L. Improved PSO optimized LSSVM for liquid rocket engine fault detection. Comput. Simul. 2020, 37, 49–51. [Google Scholar]
  43. Komadja, G.C.; Rana, A.; Glodji, L.A.; Anye, V.; Jadaun, G.; Onwualu, P.A.; Sawmliana, C. Assessing ground vibration caused by rock blasting in surface mines using machine-learning approaches: A comparison of CART, SVR and MARS. Sustainability 2022, 14, 11060. [Google Scholar] [CrossRef]
  44. Nguyen, H.; Choi, Y.; Bui, X.-N.; Nguyen-Thoi, T. Predicting Blast-induced ground vibration in open-pit mines using vibration sensors and support vector regression-based optimization algorithms. Sensors 2020, 20, 132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Yue, C.; Wu, Y.X.; Wei, Z.; Wang, G.; Wang, Y.; Li, X.; Zhou, Y.S. Prediction of blasting vibration effects in open pit mine based on PSO-LSSVM model. Eng. Blasting 2020, 26, 1–8. [Google Scholar]
  46. Coello, C.; Lechuga, M.S. MOPSO: A proposal for multiple objective particle swarm optimization. In Proceedings of the 2002 Congress on Evolutionary Computation. CEC’02 (Cat. No.02TH8600), Honolulu, HI, USA, 12–17 May 2002. [Google Scholar]
  47. Wang, D.; Tan, D.; Liu, L. Particle swarm optimization algorithm: An overview. Soft Comput. 2018, 22, 387–408. [Google Scholar] [CrossRef]
  48. Qing, S.; Tao, X.; Tatsuo, Y.; Nan, S. Multi-objective optimization of gas metal arc welding parameters and sequences for low-carbon steel(Q345D) T-joints. J. Iron Steel Res. Int. 2017, 24, 544–555. [Google Scholar]
  49. Cai, J.D.; Zheng, B.X.; Wang, X.G.; Li, P.F. Research on block size prediction and blasting control technology for multiple-size stone mining. J. Rock Mech. Eng. 2012, 31, 1462–1468. [Google Scholar]
  50. Rumelhart, D.E.; Hinton, G.E.; Williams, R.J. Learning representations by back-propagating errors. Nature 1986, 323, 533–536. [Google Scholar] [CrossRef]
  51. Yu, J.C.; Jin, A.Y.; Pan, J.W.; Wang, J.T.; Zhang, C.H. Seismic vulnerability analysis of arch dams based on GA-BP neural network. J. Tsinghua Univ. (Nat. Sci. Ed.) 2022, 62, 1321–1329. [Google Scholar] [CrossRef]
  52. Guo, Q.P.; Yang, S.J.; Zhu, Z.H.; Xiang, Z.B.; Zhang, Z.H.; Hu, G.Q. Blast vibration velocity prediction using GA-BP neural network. Blasting 2020, 37, 148–152. [Google Scholar]
  53. Zhao, D.Q.; Li, C.X.; Lan, S.N. A time lag control algorithm for structural seismic response based on least squares support vector machine. Vib. Shock. 2013, 32, 165–172. [Google Scholar] [CrossRef]
Figure 1. A view of Zhoushan tuff mine.
Figure 1. A view of Zhoushan tuff mine.
Applsci 13 07166 g001
Figure 2. Rock blasting parameters.
Figure 2. Rock blasting parameters.
Applsci 13 07166 g002
Figure 3. Monitoring point arrangement (a is the area for the verification of the optimized blasting parameters scheme).
Figure 3. Monitoring point arrangement (a is the area for the verification of the optimized blasting parameters scheme).
Applsci 13 07166 g003
Figure 4. PSO optimization LSSVM model process.
Figure 4. PSO optimization LSSVM model process.
Applsci 13 07166 g004
Figure 5. Two-objective optimization process of MOPSO.
Figure 5. Two-objective optimization process of MOPSO.
Applsci 13 07166 g005
Figure 6. Fitness curve of PPV.
Figure 6. Fitness curve of PPV.
Applsci 13 07166 g006
Figure 7. PSO-LSSVM regression fitting to PPV.
Figure 7. PSO-LSSVM regression fitting to PPV.
Applsci 13 07166 g007
Figure 8. RMSEs of LSSVM and BP in the modeling.
Figure 8. RMSEs of LSSVM and BP in the modeling.
Applsci 13 07166 g008
Figure 9. Regression curve of predicted and measured PPV for training by different methods.
Figure 9. Regression curve of predicted and measured PPV for training by different methods.
Applsci 13 07166 g009
Figure 10. Regression curve of predicted and measured PPV for testing by different methods.
Figure 10. Regression curve of predicted and measured PPV for testing by different methods.
Applsci 13 07166 g010
Figure 11. Pareto solution for two-objective optimization.
Figure 11. Pareto solution for two-objective optimization.
Applsci 13 07166 g011
Table 1. Mechanical parameters of rocks in the test area.
Table 1. Mechanical parameters of rocks in the test area.
Rock
Type
Density
(Kg/m3)
Cohesion
(MPa)
Angle of Internal Friction (°)Compressive Strength
(MPa)
Tensile Strength
(MPa)
Modulus of Elasticity (GPa)
rhyolite porphyry26800.94020640.50.26
Table 2. Input and output parameter characteristics.
Table 2. Input and output parameter characteristics.
IndexesW
(kg)
R
(m)
H
(m)
B
(m)
S
(m)
H0
(m)
PF
(kg/m3)
PPV
(cm/s)
minimum1561.820350.50.250.14
mean174.48124.4431.573.136.131.370.328.21
maximum1951256543.36.52.590.4433.58
standard
deviation
12.58034.530137.3800.1470.3400.6010.0489.250
Table 3. The initial parameters of PSO-LSSVM.
Table 3. The initial parameters of PSO-LSSVM.
ModelnVar *c1c2 ε 0 * npoppopmax *popmin *vmax *vmin *Imax
PSO-LSSVM21.51.70.94010000.011−1100
* nVar: number of decision variables; ε 0 : initial inertia weight coefficient; popmax: optimization parameter maximum value; popmin: optimization parameter minimum value; vmax: maximum update speed; vmin: maximum update speed.
Table 4. Parameters setting of PSO.
Table 4. Parameters setting of PSO.
ModelnVarc1c2 ε 0 npopnrep *nGrid *Gamma *Beta *mu *
MOPSO51.51.70.92001007220.1
* nrep: repository size; nGrid: the grid dimension of the decision variable; gamma: elimination factor of non-dominated solutions; beta: selection factor for non-dominated solutions; mu: mutation rate.
Table 5. Performance indicators of the prediction models.
Table 5. Performance indicators of the prediction models.
ModelsTraining DatasetsTest Datasets
RMSEMAErRMSEMAEr
PSO-LSSVM1.9121.2620.9831.9541.7170.965
LSSVM3.4092.2770.9393.7403.2610.887
GA-BP2.7231.9890.9612.9072.0200.906
BP3.4382.6420.9354.1272.6610.878
Table 6. Three sets of classical solutions for the Pareto front.
Table 6. Three sets of classical solutions for the Pareto front.
SolutionW (kg)PF (kg/m3)S (m)B (m)H0 (m)PPV (cm/s)Br (%)
A169.950.296.333.211.260.6830.196
B173.680.385.963.191.910.8930.098
C180.320.415.023.042.351.1090.041
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, J.; Zhao, P.; Li, P. Prediction and Optimization of Blasting-Induced Ground Vibration in Open-Pit Mines Using Intelligent Algorithms. Appl. Sci. 2023, 13, 7166. https://doi.org/10.3390/app13127166

AMA Style

Guo J, Zhao P, Li P. Prediction and Optimization of Blasting-Induced Ground Vibration in Open-Pit Mines Using Intelligent Algorithms. Applied Sciences. 2023; 13(12):7166. https://doi.org/10.3390/app13127166

Chicago/Turabian Style

Guo, Jiang, Peidong Zhao, and Pingfeng Li. 2023. "Prediction and Optimization of Blasting-Induced Ground Vibration in Open-Pit Mines Using Intelligent Algorithms" Applied Sciences 13, no. 12: 7166. https://doi.org/10.3390/app13127166

APA Style

Guo, J., Zhao, P., & Li, P. (2023). Prediction and Optimization of Blasting-Induced Ground Vibration in Open-Pit Mines Using Intelligent Algorithms. Applied Sciences, 13(12), 7166. https://doi.org/10.3390/app13127166

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop