Blasting Vibration Control Using an Improved Artificial Neural Network in the Ashele Copper Mine

Blasting is currently the most important method for rock fragmentation in metal mines. However, blast-induced ground vibration causes many negative effects, including great damage to surrounding rock masses and projects and even casualties in severe cases. +erefore, prediction of the peak particle velocity (PPV) caused by blasting plays an important role in reducing safety threats. In this paper, a genetic algorithm (GA) and an artificial neural network (ANN) algorithm were jointly used to construct a neural network model with a 4-5-1 topology to predict the PPV. For this model, the ANN parameters were optimized using the GA, and the deviating direction, horizontal distance, vertical distance, Euclidean distance, explosive type, burden, hole spacing, and maximum charge per delay were used as input information. Moreover, principal component analysis (PCA) was used to extract the first four principal components from the eight input factors as the four inputs of the ANNmodel. +e model was successfully applied to protect an underground crushing cave from blasting vibration damage by adjusting the blasting parameters. Compared with several widely used empirical equations, the GA-ANN PPV prediction model produced significantly better results, while the Ambraseys–Hedron method was the best of the empirical methods. +erefore, the improved GA-ANN model can be used to predict the PPV on site and provide a reference for the control of blasting vibration in field production.


Introduction
Blasting is currently the most economical and widely used method for artificial caving in underground mines. However, the utilization efficiency of the energy released by explosive charges is very low. According to statistics, only 20-30% of the explosive energy is used for rock fragmentation, and the rest of the energy is lost in other forms, resulting in some negative effects, such as surrounding rock vibration, flyrock, shock waves, tremendous noise, and a large amount of released heat [1][2][3][4]. e negative effects result in a large amount of energy loss and cause severe damage to the surrounding rock masses, which easily leads to rock bursts, large-scale collapses, and other disasters, further threatening the safety of underground personnel and equipment [5][6][7][8]. In particular, as the mining depth increases, the geostatic stress also gradually increases, and the issue of blast-induced rock bursts becomes increasingly prominent. According to statistics, approximately twothirds of rock bursts in mines occur as a result of blastinduced disturbances [9]. e roadway of the −647 m middle section of the Hongtoushan copper mine experienced severe damage over a distance of tens of meters. According to the literature, the peak particle velocity (PPV) is the most commonly used predictor to evaluate blasting vibration [10]. erefore, conducting studies on PPV prediction would help control the blasting vibration and reduce the corresponding hazards.
According to the literature, the PPV is influenced by many factors, including the maximum charge per delay, burden, hole spacing, total charge, number of delays, rock mass quality, and stress [11,12]. Many scholars have studied the relationship between these influencing factors and the PPV. Based on these relationships, empirical methods, numerical simulation methods, linear regression methods, and other methods have been used to predict the PPV, and many meaningful achievements have been obtained [4]. To date, there is still no widely accepted method. In recent decades, with the development of computer technology, neural networks have been increasingly applied in PPV prediction due to the advantages of self-adaptivity and nonlinearity [13,14].
In a study on the blasting parameter-based prediction of blasting vibration using a neural network approach, Singh and Singh [15] found that the neural network approach accurately predicted ground vibration, and the error did not increase as the number of inputs and nonlinearity increased. Khandelwal and Singh [16] used an artificial neural network (ANN) approach to predict the blast-induced ground vibration of a surface coal mine and compared it with the conventional statistical analysis method. Khandelwal and Singh [17] established an ANN model to predict blast-induced ground vibration. e model contained 10 input parameters, 15 hidden neurons, and two output parameters. A comparative analysis demonstrated that the model has a higher prediction accuracy than conventional predictors and multivariate regression analysis. Monjezi et al. [18] established a four-layer neural network model with a 4-1085-1 topology to predict the PPV and determined the distance from the blasting face to be the parameter that had the greatest impact on the PPV by sensitivity analysis. Amnieh et al. [19] used an ANN to establish the relationship between the PPV, block volume, explosive types, and blasting parameters and used this relationship to determine the appropriate burden, hole spacing, and maximum charge corresponding to the PPV values for a safe blasting pattern. Alvarez-Vigil et al. [20] divided the influencing factors into two categories, viz., environmental conditions and waveform parameters, proposed an ANN model for the prediction of ground vibration caused by blasting in an open-pit mine, and verified the effectiveness of this method by comparing it with multiple linear regression. Jayawardana et al. [21] used a neural network to study the relationship between the vibration mitigation effects of dual in-filled trenches under different conditions and the control parameters and established a vibration prediction model based on the relationship. Some scholars combined the neural network method with other methods for the prediction of blasting vibration. Verma and Singh [22] established a PPV predictor using a genetic algorithm (GA) based on 127 data sets. ey compared its predictors with traditional vibration predictors, multivariate regression analysis, and ANNs and discussed the superiority of the GA network. Mohamed [23] separately used an ANN and fuzzy logic analysis to predict the PPV and found that the fuzzy model had a slightly higher prediction accuracy than the ANN model. Armaghani et al. [24] proposed a flyrock and PPV prediction model based on the particle swarm optimization (PSO) algorithm and an ANN and input the top 10 influential blasting parameters to estimate the flyrock and PPV. ey found that the maximum charge per delay was the blasting parameter that influenced the flyrock distance and PPV the most. Hajihassani et al. [25] used a hybrid ANN and PSO algorithm to predict the air blast overpressure and obtained a training correlation coefficient of 0.94, indicating that the PSO-based ANN model has a high prediction accuracy. Hajihassani et al. [26] established a hybrid algorithm model for PPV prediction (an ANN algorithm and an imperialist competitive algorithm (ICA)), the ICA-ANN model, based on the blasting data of a granite quarry in Malaya, and the application results showed that the ICA-ANN model has a high prediction accuracy. Azimi et al. [27] used the monitoring data of blast-induced ground vibration from the Sungun copper mine to establish a GA-ANN prediction model and found that the prediction accuracy was higher when the horizontal and radial distances were replaced with the modified radial distance. Nguyen et al. [28] analyzed and affirmed the prediction accuracy of blast-induced PPV using a hybrid model based on the hierarchical k-means (HKM) clustering algorithm and an ANN (an HKM-ANN model) and a hybrid model based on fuzzy c-means clustering (FCM) and support vector regression (SVR) and concluded that the HKM algorithm could effectively improve the prediction accuracy of ANN models.
Based on the monitoring data of ground vibration caused by underground blasting in the Ashele copper mine, the present study established a PPV prediction model by combining the GA with the ANN. In this model, the ANN parameters were optimized by the GA, and the optimal number of nodes and connection weights of the hidden layers in the neural network were obtained. e model was applied to provide guidance for optimizing the blasting parameters and controlling blasting vibration.

Combined Method Based on a Neural
Network and GA e backpropagation (BP) network is currently one of the most widely used neural network models. e BP network generally consists of an input layer, hidden layers, and an output layer. e learning algorithm of the BP network can be summarized as the process of iteratively performing forward computation, reverse computation, and weight and threshold updating according to a given sample set. e training weight coefficient is reversed. e error of the middle layer is obtained based on the error of the output, and then the weight coefficient of the input is obtained. is process is repeated until the correct output is obtained. e BP network is essentially the input-output mapping and does not require any precise mathematical expression between the input and the output. As long as a known model is used to train the BP network, the network has the capability of mapping between input and output pairs. It actually transforms the input/output problems of a sample set into a nonlinear optimization problem and uses the most common gradient descent optimization method and iterative operation to solve for the weights. e neural network uses a gradient descent-based BP algorithm and its variants, and its dependence on the initial weight selection inevitably leads to the problem of a local optimum, resulting in a decrease in the reliability of the neural network approach. In view of the shortcomings of the BP neural network, this paper adopts an improved GA-BP neural network to overcome the shortcomings of the BP neural network. e GA simulates the evolutionary phenomena of biological inheritance, genetic variation to adapt to the environment, and survival of the fittest in nature [29]. It treats each solution as a chromosome, each element of the solution as a gene, and all chromosomes as a population. e chromosomes are evaluated using a predetermined objective function to calculate the fitness. e algorithm initially randomly generates some chromosomes, performs genetic operations such as genetic inheritance and variation on the chromosomes based on the fitness, selects chromosomes with a high fitness, and conducts the next generation of genetic operations to obtain new populations. Because each iteration selects the elites from the population of the previous generation, the new generation inherits the superior genes of the previous generation and reorganizes in time to produce better chromosomes and eliminate defective new individuals. erefore, through repeated iterations, the solution evolves toward the optimal solution until the termination condition is met. Due to the high level of parallelism, the GA has been widely applied in neural networks.

Production and Monitoring of the Ashele Copper Mine.
e Ashele copper mine, located in Xinjiang Province in China, has a copper reserve of 920,000 tons, with an average grade of 2.4%. e main ore body has a strike length of 700 m, an average thickness of 20 m to 45 m, and an inclination of 55°to 85°. e deposit is mainly formed by volcanic ejection and deposition and is then remolded by hydrotherm fluids. e average rock quality designation (RQD) of the surrounding rock in the hanging wall is 81.2%, with 41.5% of the surrounding rock in the footwall and 70.9% of the ore body. e current mining depth exceeds 1000 m, with 1243 m of the shaft length. e hanging wall of the ore body is mainly composed of dacite porphyry, and the footwall of the ore body is mainly basalt, with chalcopyrite as the main mineral. e results of an in situ stress test showed that the maximum horizontal principal stress at the maximum mining depth was close to 45 MPa, with an inclination of 18°, which was almost perpendicular to the strike of the ore body, and the minimum principal stress was approximately 20 MPa, with an inclination of 64°. e Ashele copper mine mainly employs the large-diameter deep-hole mining method, which allows a high production capacity. Ores were gathered from the bottom using trenches. After the formation of the trenches, the vertical crater retreat (VCR) method was adopted to form cutting grooves at the mining face, and then large-hole lateral caving was used for mining.
e Ashele copper mine produces 6,000 tons of ore daily; the charge for a single blast is as high as 4 tons, and the maximum charge per delay reaches 400 kg. Frequent blast and stress transfers due to excavation disturb the surrounding rock mass of roadways, in some cases causing severe damage and requiring significant effort to maintain and rehabilitate the secondary support. Stope collapses have occurred at two different levels of the Ashele copper mine, which were attributed to the blasting disturbance of adjacent stopes after microseismic activity analysis [30]. e NUBOX-6016 intelligent vibration monitor and the NUBOX-6016 portable data acquisition device used in this study were designed based on optimization tests for blasting, vibration, shock, and noise on site (see Figure 1). e instruments can digitally convert and store the dynamic and static analog signals generated by sensors (including the acceleration, pressure, strain, and temperature) and have a triggering mechanism to ensure the correct recording of signals only for the features concerned. ey allow data acquisition for up to 2048 segments to achieve continuous automatic recording of multisegment vibration signals. e maximum sampling frequency is 200 kHz.
e signal-tonoise ratio (SNR) is ≥62 dB. e minimum sampling frequency is 5 Hz. is equipment is equipped with a TP3V-4.5 sensor, which has the advantages of easy installation, high robustness, reliability, compactness, high measurement precision, and strong anti-interference. e stope of the Ashele copper mine adopted retreats mining from the hanging wall to the footwall, and during the blasting vibration field test, the monitoring points were placed in the footwall of the ore body to avoid measuring the blasting vibration waves that propagated through the goaf.
Every two or three levels are used as a mining combination in the Ashele copper mine. e mining sequence is from bottom to top in the mining combination. At present, the first stope at a depth of 1012 m is to be mined according to the production plan. Unfortunately, a crushing chamber is located near the stope, which is used to crush the ore produced below 600 m depth. Apparently, it has a long service cycle. erefore, it is necessary to control the PPV of the stope blast to protect the crushing chamber.

Factors
Controlling the PPV. Samples of the inputs and outputs should be offered to train and test the blasting vibration prediction model based on the neural network. e related factors concerning the blasting vibration should be selected carefully. en, the inputs of the ANN model, including blasting parameters and field conditions, were determined after investigation. Field conditions mainly include the buried depth, maximum principal stress, deviating direction, horizontal distance, vertical distance, and Euclidean distance. e blasting parameter mainly contains the explosive type, burden, hole spacing, delay number, total charge, and maximum charge per delay.

Weak Effect of Burial Depth and Maximum Stress.
ere are few studies considering the impact of the buried depth and the maximum stress on blasting vibration. Little correlation between the buried depth, the maximum stress, and the PPV was found before. In this paper, the blasting vibration at a deeper level was studied based on the shallow level, so the two impactors were taken to analyze the impact on the PPV. e blasting vibration data were all collected at 562 m depth, 762 m depth, and 912 m depth, which were the main mining levels in the Ashele copper mine. e PPV values of different depths and maximum principal stresses Shock and Vibration 3 are shown in Figure 2.
e PPV tends to increase with increasing burial depth. It can be predicted that the influence of burial depth on PPV will become increasingly prominent with increasing mining depth. According to the results of the in situ stress test of the Ashele copper mine, the maximum principal stress of the in situ stress is near horizontal and approximately perpendicular to the ore body. In addition, there is little change in the direction of the maximum principal stress. Hence, only the value of maximum principal stress is considered for the model. e PPV grows with increasing maximum horizontal principal stress. However, the correlation coefficients between the buried depth, maximum stress, and PPV are small. It can be concluded that there is no significant correlation between them. erefore, the buried depth and the maximum stress are not taken as the inputs of the PPV prediction. At the same time, the correlation coefficients between the delay number, total charge, and PPV are small, illustrating that the two factors have little impact on the PPV. erefore, the delay number and the maximum stress are eliminated.

Deviating Direction.
e scatterplot matrix of the PPV database with correlation is shown in Figure 3. e deviating direction is the angle between the stope centerline and the line linking the blasting center and the monitoring point. e deviating direction is a positive number when the stope centerline rotates clockwise. Otherwise, the deviating direction is a negative number. Apparently, the PPV is the smallest when the deviating direction is 0. With the increase in the absolute value of the deviating direction, the PPV rises gradually. e correlation coefficient between the deviating direction and PPV is 0.307, and a significance test finds significance at the level of 0.05. It can be concluded that the deviating direction has some influence on the blasting vibration.

Horizontal Distance, Vertical Distance, and Euclidean Distance between the Blasting Source and Monitoring Point.
Most empirical formulas incorporate the influence of propagation distance on blasting vibration. e horizontal distance is the horizontal projection of the distance between the blasting center and the monitoring point. Clearly, there is a high correlation between the horizontal distance and PPV. As we thought, the PPV decreases as the horizontal distance increases. Similar to the horizontal distance, the vertical distance, which is the vertical projection of the distance between the blasting center and the monitoring point, has an obvious influence on blasting vibration. However, there are some differences between the influences of horizontal distance and vertical distance on blasting vibration. e vertical propagation distance should be considered when studying the propagation of blasting vibration at different levels. With the increase in the Euclidean distance between the blasting center and monitoring point, the PPV decays gradually. In other words, with increasing horizontal distance, vertical distance, and Euclidean distance, the PPV decreases gradually. e correlation coefficients of the horizontal distance, vertical distance, and Euclidean distance are 0.607, −0.587, and −0.73, respectively. ey are significantly correlated with the PPV at the 99% confidence level. erefore, they are all incorporated into the PPV prediction model.

Factor of Explosive Type.
e blasting speed, power and other parameters vary greatly for different explosive types. In this paper, two types of explosives are mainly used to break the rock. ey are emulsion explosive and ammonium oil explosive, which are represented by 1 and 2, respectively. e blasting vibration generated by blasting with emulsion explosive is significantly greater than that of blasting with ammonium oil explosive. e correlation coefficient between the explosive type and PPV is −0.392,  which is significant at the 0.01 level. According to the theory proposed by Cohen [31], an absolute value of coefficient greater than 0.3 can be regarded as moderate correlation. Consequently, the explosive type is taken as one of the inputs of the model.

3.2.5.
Burden. e burden is the distance from the explosive cartridges to the free surface. When the burden is less than 2.5 m, the PPV increases with decreasing burden. Moreover, the PPV increases with the growth of the burden when the burden is more than 2.5 m. e correlation coefficient between the burden and PPV is 0.324, which is relatively small. However, a significance test identifies significance at the level of 0.05. erefore, the burden can be considered to have an impact on the PPV and can be used to predict PPV.

Hole Spacing.
With increasing hole spacing, the mean values and the minimum values of the PPV present an increasing trend. When the hole spacing is 2.0 m, the PPV values were relatively concentrated. e value range of the PPV is large at most hole spacings. e correlation coefficient between the hole spacing and PPV is 0.350, which is significant at the 0.05 level. As a result, the hole spacing is selected as the input of the model.

Maximum Charge per Delay.
Similar to the distance, the maximum charge per delay often appears in empirical formulas. When the maximum charge per delay equals 220 kg, the PPV is the smallest according to the fitting curve.
e PPV increases gradually with the increase in the maximum charge per delay after exceeding 220 kg. e PPV also increases with the decrease in the maximum charge per delay when it is less than 220 kg. e correlation coefficient between the maximum charge per delay is 0.381, which is significant at the level of 0.01 (see Figure 4). Hence, the maximum charge per delay is selected for model input.
As a result, the deviating direction, horizontal distance, vertical distance, Euclidean distance, explosive type, burden, hole spacing, and maximum charge per delay are selected to predict the PPV.

Standardization of the Data.
e values of all kinds of information vary greatly. e ratio of the largest to the smallest is more than 1000. It is necessary to standardize the input data to eliminate the influence of dimension and variation size. In this paper, the Z-score method, which is the most widely used standardized method at present, is used to standardize the input data of the model. e Z-score method is shown in equation (1). After establishment of the standardized model, all the prediction sample data need to be standardized according to the same model of the training samples.

Principal Component Analysis (PCA) of Input
Information. e input data of the ANN model contain eight kinds of variables according to Section 3.1. e input data contain considerable repetitive and redundant information, and an obvious correlation is found between each factor and PPV. For example, the layout characteristics of PPVs with different horizontal distances and Euclidean distances are similar. e reason for this phenomenon is that the horizontal distance has a high correlation with the Euclidean distance in the Ashele copper mine. e correlation coefficient between them is 0.930. In other words, some information overlap exists.
ere are no apparent direct connections between the blasting parameters, such as the hole spacing and burden. If the explosive consumption is constant, the hole spacing is restricted by the burden. e hole spacing and burden are not completely independent. e existence of this repetitive and redundant information directly affects the prediction accuracy of PPV, which requires some data processing before use for the ANN model. e method of PCA extracts the key information to establish new variables that are unrelated to each other. At the same time, it deletes redundant information while making the best effort to preserve the original information. In this paper, PCA is used to eliminate the influence of the repetitive and redundant information included in the data of eight factors after standardization. e error contribution of PCA is shown in Figure 5. e larger the variance percentage is, the greater the contribution of the principal component to the output. e variance contributions of the first four principal components are 39.20%, 31.75%, 14.43%, and 7.78%. e cumulative variance contribution of the first four principal components accounts for 93.16%. In this paper, the first four principal components are selected to predict PPV. e eigenvalues of the four principal components are 3.136, 2.540, 1.154, and 0.622.

Establishment of the PPV Prediction Model.
Generally, the training error decreases rapidly in the training process, while the testing error sometimes decreases with relatively narrow fluctuations. Referring to the early stopping method, the maximum number of learning sessions was 200,000. en, the best number of learning steps was determined to reduce the overfitting problem according to the error of learning and testing. In addition, a simple structure of the neural network model is helpful to avoid overfitting problems. Consequently, a topology with one hidden layer was used. e hidden layer had a node search range of 5-15. e neural network model had four inputs and one output. e population size was 200, the total number of generations was 30, the code string length was 15, the crossover probability was 0.8, the mutation probability was 0.2, the inversion operator was 0.3, and the initial weight was from −0.5 to 0.5. e individual error and test target error were controlled at 10 −6 and 10 −5 , respectively. e inertia coefficient was set to 0.5, and the generation with a stable fitness Shock and Vibration 5 value was set to 5. In this paper, the mean squared error of the outputs from the measured results was chosen to calculate the fitness value. e above parameter setting was used for the nine generations of genetic evolution. e best fitness value for each generation is shown in Figure 6. e optimal fitness value of the initial population was 1.711. e best fitness value of the first-generation population was 1.700, indicating a significant decrease. en, a slight decrease in the best fitness value was found in the second generation. ere was no significant change in the third and fourth generations, indicating that in the current genetic operation process, no chromosomes were better than those of the second generation. Afterwards, the best fitness value of the fifth generation was 1.207, showing a significant decrease and indicating that during the evolution from the fourth generation to the fifth generation, a better chromosome was generated, so the best fitness value was greatly reduced. Until the ninth generation, none of the individuals were better than those of the fifth generation. At this point, the generation with the stable best fitness value was set to 5, and the genetic evolution was terminated. In this case, the best fitness values for the second to fourth generations were the same, as were the corresponding weights and topological parameters of the neural network. erefore, the best fitness value and the corresponding weights and topological parameters of the neural network of any generation from the fifth generation to the ninth generation can be applied to the PPV neural network. e fitness values of each individual in the ninth generation and their corresponding topological parameters in the neural network are shown in Figure 7.   ere were 50 individuals in the population. e best fitness value of the ninth generation corresponded to individual no. 31, and the topological parameter of its neural network was 5. Finally, the number of nodes in the hidden layer of the neural network was 5. Based on this, a PPV prediction model with a 4-5-1 topology was established. is model had one input layer, one hidden layer, and one output layer. e input layer contained four inputs, the hidden layer contained 5 nodes, and the output layer contained one output. In addition, there were 31 initial weights. e frame diagram of the three-layer neural network is shown in Figure 8.
e outputs of this model during learning and testing are shown in Figure 9. During learning and testing, the output values of the GA-ANN model were very close to the measured values. e errors of all learning and testing cases were within an acceptable range.

Blasting Vibration Control for Crushing Cave Protection.
ere is a crushing cave near the blasting stope at the 1012 m mining depth of the Ashele copper mine, which is shown in Figure 10. Two blasting vibration monitoring points were arranged near the stope to collect the PPV generated by the blasts. e blasting parameters and field conditions are shown in Table 1. e PPV values at points 1 and 2 induced by large-hole blast were 116.1 mm/s and 166.2 mm/s according to the PPV prediction model and blasting design, respectively. e PPV at monitoring point 2 caused by a large-hole blast is close to the minimum safety value recommended in the blasting safety regulation. e blasting parameters were adjusted to control the blasting vibration to avoid damage to the crushing cave caused by blasts near the stope. e burden and hole spacing were changed to 3.5 m and 2.5 m from 3.3 m and 3.4 m, respectively. Moreover, the maximum charge per delay was reduced from 400 kg to 300 kg, while the delay number increased from 10 to 13. e corresponding delay times were 0, 50 ms, 110 ms, 150 ms, 200 ms, 250 ms, 310 ms, 380 ms, 460 ms, 550 ms, 650 ms, 760 ms, and 880 ms. After the adjustment of the blasting parameters, the predicted PPVs at monitoring points 1 and 2 were 55.60 mm/s and 87.80 mm/s, respectively. Finally, the blast of the near stope did not cause apparent damage to the crushing cave. e measured PPVs at points 1 and 2 caused by the large-hole blast and medium-deep hole blast were 56.28 mm/s and 96.99 mm/s, respectively.

Discussion
According to previous studies, the PPV is influenced by many factors, including environmental parameters such as the stress and rock mass quality and blasting parameters such as the maximum charge per delay, total charge, Euclidean distance, and number of delays. Most previous studies have considered the maximum charge per delay and the Euclidean distance to be the two most important parameters affecting the PPV. For example, empirical PPV predictors, such as those of the United States Bureau of Mines [32], Langefors and Kihlstrom [33], the Bureau of Indian Standards criteria [34], Ambraseys and Hendron [35], and the Central Mining Research Institute (CMRI) predictors [36], consider only the influences of the maximum charge per delay and the Euclidean distance on the PPV. ese empirical formulas are shown in Table 2.
Using the learning and testing samples, the above five empirical models were used to establish the PPV prediction equations to predict the PPV values at different sites during the two blasts. e results predicted by the empirical equations were relatively poor (see Table 3).
ere is no significant difference between the PPV values predicted by the five empirical equations. is is because the same factors are taken into account by the five empirical equations with different expressions. e Ambraseys-Hendron equation was the best among these empirical equations for PPV prediction, and its prediction errors were approximately 13.36% and 87.81%. As a result, the PPV values predicted by the GA-ANN algorithm were the closest to the measured values, and the prediction errors were less than 10%, indicating that the prediction accuracy of the GA-ANN algorithm was the highest (see Figure 11). It is believed that the PPV prediction accuracy will be further improved with the accumulation of data. In general, the PPV should be predicted according to the actual factors rather than simply using empirical equations.   at two different monitoring points are very similar. In addition, the maximum charge per delay and Euclidean distance are found to significantly influence the PPV through correlation analysis. In fact, the influence of deviating direction, hole spacing, burden, and explosive type on the PPV cannot be ignored according to the data in this paper.
Coupled with the advantages of neural networks in solving nonlinear problems, these reasons make the predicted PPV using the model in this paper more accurate than the empirical formulas. is can also be seen from the PPV prediction results of the two monitoring points. Compared with the traditional empirical formula, the predicted PPV values of the neural network at the two monitoring points with different deviating directions have marked differences.

Conclusion
is paper proposed a GA-ANN PPV prediction method. e input information of the PPV prediction method contains field conditions and blasting parameters, such as the deviating direction, horizontal distance, vertical distance, Euclidean distance, explosive type, burden, hole spacing, and maximum charge per delay. PCA was used to eliminate redundant and repetitive information of the field conditions and blasting parameters. e mean squared error of the outputs from the measured results was chosen as the fitting function to search for the best node number of the hidden layer. In this way, the GA-ANN PPV prediction model with a 4-5-1 topology was established for the Ashele copper mine. In the learning and testing data, the predicted results were in good agreement with the measured values. is model was used to predict and control the blasting vibration of a deep stope located at a 1012 m buried depth in the Ashele copper mine.
rough the adjustment of blasting parameters, such as the maximum charge per delay, delay number, burden, and hole spacing, the blasting vibration was successfully reduced, and the crushing cave was protected from damage. By comparing the prediction results of the GA-ANN PPV prediction model with those of traditional empirical PPV prediction models, the superiority of the GA-ANN PPV prediction model was verified. Among the traditional PPV prediction methods, the Ambraseys-Hendron method was the best. e proposed GA-ANN PPV model can be used to guide the design of blasting parameters for field production.
Data Availability e datasets used in the present study are available from the corresponding author upon reasonable request.