A Comprehensive Evaluation Method of Bench Blast Performance in Open-Pit Mine

: In blasting operation, some undesirable impacts, such as ﬂy-rock, fragmentation, and back break, are induced. If the blasting design is not optimized, these mentioned impacts would reduce the blasting e ﬃ ciency. To improve and optimize the blast design, blasting e ﬀ ect evaluation is essential. Due to the complexity of interactions among blasting parameters, empirical methods may not be appropriate for blast design optimization. A two-level mathematical model based on fuzzy mathematics, is proposed in this work. In total, 11 typical parameters were chosen and classified into three groups. The blasting e ﬀ ect is evaluated from three aspects, and then the comprehensive evaluation is given. A blasting e ﬀ ect evaluation system was developed based on the mentioned method on the platform of VC ++ . Some other techniques, such as image processing, were integrated into the system, which allowed for obtaining all of the parameters rapidly and conveniently. The system was applied in practical bench blast engineering. The results obtained from the system can provide e ﬀ ective information for the optimization of the next blast design.


Introduction
The bench blast has been widely applied in many fields, such as mining engineering and hydraulic tunnel excavation. The effects of the bench blast would directly affect subsequent operations, such as loading, hauling and crushing. Otherwise, the bench blast would also cause some undesirable impacts, such as fly-rock, vibration and back-break, if the blast design is not optimized [1]. To improve and optimize the bench blast design, a comprehensive, valid and convenient for application evaluation method of blast performance is therefore expected by mining engineers.
As efforts to find more scientific and the reliable methods to evaluate bench blast performance, various research studies were initiated. These works can be divided into three types from different viewpoints: (1) statistical analysis method; (2) harmful effects evaluation method; (3) numerical modeling method.
The statistical analysis method mainly employed the Artificial Neural Network (ANN) [2][3][4][5][6][7]. The ANN is an information processing system that is used for modeling complex relationships between input and output data. An ANN consists of an input layer, hidden layers and an output layer. Once the network has been trained by the previous blasting parameters and results, the results can be predicted by the network according to the new data inputs. However, in hidden layers, all the input parameters are calculated as whole, so this method can only show approximate relationships between parameters and results.
The numerical simulation method has been widely applied to predict the results of the bench blast, such as the Finite Element Method (FEM) [8][9][10], Discrete Element Method (DEM) [11][12][13] and Discontinuous Deformation Analysis (DDA) [14][15][16]. The feedback from simulated results can provide a direction for design. However, the numerical model needs many parameters, such as geological and construction conditions, and some parameters are unavailable or difficult to be determined. The simulated results always have a significant difference with the results measured on site. Therefore, the numerical modeling method is only considered as a reference for bench blast design.
The harmful effects evaluation method mainly focuses on predicting the harmful effects induced by bench blasts [17][18][19][20]. For example, the famous Sadovsky equation [21] was employed to predict the peak particle velocity (PPV) for a certain bench blast design, and some empirical formulas were obtained to evaluate the distance of the fragments flying away from the blast source [22]. In this method, the single factor of bench blast performances was analyzed, instead of the comprehensive evaluation of the bench blast. Although they are still widely used due to their simplicity and relatively good reliability of results, these equations have one-sidedness because they only consider the harmful effects of blasting to evaluate the blasting effect.
To overcome the limitations of the abovementioned methods, a comprehensive evaluation method of bench blast in an open pit mine was proposed in this work. Eleven parameters were determined to evaluate the performance of bench blasts. Then, an evaluation model based on the linear programming and the fuzzy mathematics was proposed. Furthermore, a program code was developed on the platform of VC++. An application case of a bench blast was finally analyzed.

Evaluation Parameters
Bench blasting performance has some undesirable impacts, such as vibration and back break; therefore, effect evaluation for bench blasting is quite complicated. The rationality of the established method for evaluating blasting effect is related to the fact that it could reflect the blasting effect correctly, and then provide information for optimizing the blast design.
In bench blast design, security issues are the key factors to be considered. In bench blasting performance, fly-rock is one of the most important causes of accidents, and blasting vibration is particularly harmful to buildings; therefore, the blasting safety evaluation is basic for the whole blasting effect evaluation. Blasting quality is the purpose of blast design. Blasting quality can directly affect construction efficiency, which would affect the subsequent operations, such as loading, hauling and crushing. Under the premise of ensuring blasting safety, optimizing the blast design is important for maximizing the blasting efficiency. Additionally, in bench blasting performance, energy and economic waste usually exist; therefore, the effective measures to reduce blasting costs and improve economic efficiency should be taken, which is the main factor that a blast design needs to be considered. Therefore, in this study, the blasting effect is evaluated from three aspects: blasting safety, quality and cost-11 typical parameters are chosen and are divided into these three classes. Eleven typical parameters are used to describe the effect of bench blasting. Fly-rock: the shorter the distance of flying stone, the less harm caused by blasting. Vibration: the smaller the peak value of blasting vibration, the smaller the harm caused by vibration. Boulder yield: the block size meets the engineering requirements, and the block rate is low. Muck pile height: the blasting pile is concentrated and has certain looseness, which should not be too high or too low, so as to meet the requirements of efficient shoveling of equipment. Forward distance: using the slag body reserved from the previous blasting, the forward striking distance of blasting pile body is controlled to make the blasted rock mass collide with the residual slag body, reduce the block rate and improve the energy utilization rate. Back break: under excavation, this will cause processing delay or idling, reduce productivity and waste money. Depth of lateral posterior sulcus: in order to avoid the slope collapse, the back gully should be built at an appropriate place outside the slope to remove the water. Bench smoothness: high step flatness is conducive to improve the efficiency of the next drilling and shoveling transportation. Explosive unit consumption: reasonable explosive unit consumption is conducive to improving energy utilization and saving cost. Blasted volume per meter: the higher the blasting quantity, the better the blasting effect. Auxiliary materials: lowering the consumption of auxiliary materials but maintaining the same blasting effect, meaning the cost is lower. Finally, the comprehensive evaluation is given. The parameter classifications are listed in Table 1.

Fuzzy Mathematical Model
Due to the complexity of interactions among blasting parameters, empirical methods may not be appropriate for blast design optimization, but the fuzzy model can cope with the complexity of complicated systems in a flexible and reliable way. Therefore, fuzzy mathematics has been widely applied in mining engineering. In this study, a two-level mathematical model was built based on fuzzy mathematics. In the first-level mathematical model, the subclasses (blasting safety, cost and quality) are evaluation objectives, and the parameters belonging to the subclass are constraint conditions. The first-level mathematical model is described by Equation (1): where A T A R M f E are the evaluated results of safety, quality and cost, respectively, and b j , a j are the weight and membership degrees. For the second-level mathematical model, the comprehensive blasting effect is the evaluation objective, and the subclasses are constraint conditions. The model is described by Equation (2): where B Ni is the comprehensive evaluation, and α, β, γ are the weight of safety, quality and cost, respectively.

Membership and Weight Determination
All evaluation parameters are described by linguistic variables. Then linguistic variables are translated to memberships, respectively. The linguistic variables and memberships are shown in Table 2. Weight determination is the most important factor that affects the accuracy of fuzzy comprehensive evaluation. Due to the different focus of each engineering assessment, the fixed weight method may not be suitable. In this study, the weight adjustment method is used. The initial weight is given by expert valuation for some special engineering, so if there are some unreasonable weight values, the weight can be adjusted. If the adjusted value ε is given, the new weight value W k and other new weight values W i can be recalculated by Equation (3): After the adjusted weight value is known, the other weight W i values can also be recalculated by Equation (5):

Program Development
In this study, a blasting effect evaluation system (BEES) was developed based on the aforementioned method. The detailed algorithm was designed, and the corresponding code was developed on the platform of VC++ (The version is 6.0, Microsoft). The program structure is shown in Figure 1.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 12 Weight determination is the most important factor that affects the accuracy of fuzzy comprehensive evaluation. Due to the different focus of each engineering assessment, the fixed weight method may not be suitable. In this study, the weight adjustment method is used. The initial weight is given by expert valuation for some special engineering, so if there are some unreasonable weight values, the weight can be adjusted. If the adjusted value ε is given, the new weight value ( 1 ,2 , , 1 , 1 , , ) After the adjusted weight value is known, the other weight ' i W values can also be recalculated by Equation (5):

Program Development
In this study, a blasting effect evaluation system (BEES) was developed based on the aforementioned method. The detailed algorithm was designed, and the corresponding code was developed on the platform of VC++ (The version is 6.0, Microsoft). The program structure is shown in Figure 1. Boulder yield is calculated from rock fragmentation size distribution. The rock fragmentation size distribution is an important parameter for blasting effect evaluation, but the manual measurement is tedious and time-consuming. Today, software is available for use with any camera or photo type, or cameras with inbuilt software that automatically gives rock fragmentation size distribution both in table view and graphically. So, in this system, the image analysis technique is used for rock fragmentation size statistics. The boulder yield can be calculated easily by setting the standard value. The algorithm of fragment size statistics is shown in Figure 2. The muck pile height can also be obtained by the method of image processing-the algorithm is shown in Figure 3a, and its explanation on how muck pile height is determined is shown in Figure 3b. The vibration can be measured by vibration measuring instrument. The distance meter can measure the distance of flying rock. The forward distance, under excavation, depth of side and back ditch and the flatness of terrace Boulder yield is calculated from rock fragmentation size distribution. The rock fragmentation size distribution is an important parameter for blasting effect evaluation, but the manual measurement is tedious and time-consuming. Today, software is available for use with any camera or photo type, or cameras with inbuilt software that automatically gives rock fragmentation size distribution both in table view and graphically. So, in this system, the image analysis technique is used for rock fragmentation size statistics. The boulder yield can be calculated easily by setting the standard value. The algorithm of fragment size statistics is shown in Figure 2. The muck pile height can also be obtained by the method of image processing-the algorithm is shown in Figure 3a, and its explanation on how muck pile height is determined is shown in Figure 3b. The vibration can be measured by vibration measuring instrument. The distance meter can measure the distance of flying rock. The forward distance, under excavation, depth of side and back ditch and the flatness of terrace can be obtained by Appl. Sci. 2020, 10, 5398 5 of 12 field observation and image processing. The blasting cost could be calculated by inputting the bench blasting engineering information.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 12 can be obtained by field observation and image processing. The blasting cost could be calculated by inputting the bench blasting engineering information. In this system, the default data range and weight values are given by the expert evaluation method. That is to organize the experts in the field of Open-pit Bench Blasting Engineering, use professional knowledge and experience, and through intuitive induction, give a strong professional can be obtained by field observation and image processing. The blasting cost could be calculated by inputting the bench blasting engineering information. In this system, the default data range and weight values are given by the expert evaluation method. That is to organize the experts in the field of Open-pit Bench Blasting Engineering, use professional knowledge and experience, and through intuitive induction, give a strong professional In this system, the default data range and weight values are given by the expert evaluation method. That is to organize the experts in the field of Open-pit Bench Blasting Engineering, use professional knowledge and experience, and through intuitive induction, give a strong professional judgment and evaluation (see Tables 3 and 4). However, the data ranges and weight values may be different for different construction environments and blast designs, it is not suitable for evaluation with the fixed data ranges and weight values. So, the data ranges and weight values should be changed with engineering conditions as it will be easier for deciders to obtain the advantages and disadvantages of blast design, and to better optimize the next blast design. In the system, the normal user and administrator models are provided. The administrator can adjust the data ranges and weight values according to the actual situation. The adjusted data are saved to the text in binary form and can be recalled next time.   After all the information input is completed, all the parameters are calculated and classified according to subclasses. The BEES provides dialogs for users to observe the evaluation of each parameter (see Figure 4). According to all the parameter evaluation results, the blasting effect evaluation can be obtained-the algorithm is shown in Figure 5-and the blasting effect evaluation results are displayed by the form of a graph (They can be seen in the next two figures, Figures 6 and 7).

Application
To verify the evaluation and optimization method proposed in this paper, a practical bench blast engineering in a copper deposit was selected as an application case. The photo shot on site is shown in Figure 6. This copper mine is mined by open-pit deep-hole step blasting. A No. 2 rock explosive is selected, and the structure of continuous charge is adopted. The total charge amount in the blasting

Application
To verify the evaluation and optimization method proposed in this paper, a practical bench blast engineering in a copper deposit was selected as an application case. The photo shot on site is shown in Figure 6. This copper mine is mined by open-pit deep-hole step blasting. A No. 2 rock explosive is selected, and the structure of continuous charge is adopted. The total charge amount in the blasting   Second, the BEES shows the comprehensive evaluation result-the result is shown by histogram-as shown in Figure 9, the red column represents the comprehensive evaluation result.
The blasting effect evaluation of this bench blast is between fair and good. The blasting cost evaluation is similar to the blasting effect evaluation. The blasting quality evaluation is between good and very good, but the blasting safety evaluation is poor. So, this blast design is not optimized. The next blast design should focus on burden and row spacing, etc., to control the blasting fly-rock and vibration.

Conclusions
Blasting effect evaluation is very necessary for blast design. In this paper, a method of blasting effect evaluation was proposed and an evaluation system was developed. The traditional evaluation methods are tedious and time-consuming. However, the developed evaluation system integrates several techniques, such as image analysis and parameter classification. All the selected parameters can be obtained easily to achieve rapid evaluation. The evaluation results are divided as subclass evaluation results and comprehensive evaluation results. The system has been used in practical bench blast engineering. The results are easier for deciders to obtain the advantages and disadvantages of blast design and provide the effective information for the optimization of the next blast design. The comprehensive evaluation system of blasting effect developed is not only convenient, accurate and objective, but also beneficial to enrich and perfect the fine management system of digital blasting management systems. The use of this software greatly reduces the burden of engineers and technicians and realizes the rapid statistics and evaluation of blasting effect.

Application
To verify the evaluation and optimization method proposed in this paper, a practical bench blast engineering in a copper deposit was selected as an application case. The photo shot on site is shown in Figure 8. This copper mine is mined by open-pit deep-hole step blasting. A No. 2 rock explosive is selected, and the structure of continuous charge is adopted. The total charge amount in the blasting area is 4392 kg. Detonators are used to detonate explosives. The latest Series IV advanced vibration and overpressure monitors are used for the blasting vibration measurement.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 12 area is 4392 kg. Detonators are used to detonate explosives. The latest Series IV advanced vibration and overpressure monitors are used for the blasting vibration measurement. In general, there is visual error between photo measurement and site measurement; therefore, the correction coefficient is necessary. The camera is used to measure the image of the explosive pile particles. The visible area of the blasting pile is selected to shoot in an open direction. The area of each shot is about 40-30 m. The camera should be located in the front of the shooting area, and the elevation angle of the line of sight is within 35 degrees. Insert a pole, which should be located in the middle of the field of view, with an effective visible length of 1-2 m. The image scale coefficient is determined by the size of the post on the image and the actual size of the post. Then, the error of the image measurement method is corrected according to the measured particle size of the explosive pile sample, and the correction coefficient between the image measurement method and the geometric method is obtained. In this study, the coefficient is 0.8 and the rock fragment distribution is shown in Figure 9.
With the different construction environments, the data ranges and weight values should be changed according to the construction condition. The BEES provides the default data ranges, and also provides the permission to modify the data. In this study, the data ranges are listed in Table 5, and the weight values were not modified (see Table 4).  The comprehensive blasting effect evaluation has two steps. First, in this study, a representative blasting operation was selected and 11 representative parameter values were obtained (see Table 6). According to the values of the 11 parameters, the BEES evaluates blast effect from blasting safety, quality and cost, respectively, and shows the evaluation results by histogram-as shown in Figure  8, the red column represents the evaluation result of the current index.  In general, there is visual error between photo measurement and site measurement; therefore, the correction coefficient is necessary. The camera is used to measure the image of the explosive pile particles. The visible area of the blasting pile is selected to shoot in an open direction. The area of each shot is about 40-30 m. The camera should be located in the front of the shooting area, and the elevation angle of the line of sight is within 35 degrees. Insert a pole, which should be located in the middle of the field of view, with an effective visible length of 1-2 m. The image scale coefficient is determined by the size of the post on the image and the actual size of the post. Then, the error of the image measurement method is corrected according to the measured particle size of the explosive pile sample, and the correction coefficient between the image measurement method and the geometric method is obtained. In this study, the coefficient is 0.8 and the rock fragment distribution is shown in Figure 9.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 12 blast design and provide the effective information for the optimization of the next blast design. The comprehensive evaluation system of blasting effect developed is not only convenient, accurate and objective, but also beneficial to enrich and perfect the fine management system of digital blasting management systems. The use of this software greatly reduces the burden of engineers and technicians and realizes the rapid statistics and evaluation of blasting effect.  With the different construction environments, the data ranges and weight values should be changed according to the construction condition. The BEES provides the default data ranges, and also provides the permission to modify the data. In this study, the data ranges are listed in Table 5, and the weight values were not modified (see Table 4).  The comprehensive blasting effect evaluation has two steps. First, in this study, a representative blasting operation was selected and 11 representative parameter values were obtained (see Table 6). According to the values of the 11 parameters, the BEES evaluates blast effect from blasting safety, quality and cost, respectively, and shows the evaluation results by histogram-as shown in Figure 6, the red column represents the evaluation result of the current index. Second, the BEES shows the comprehensive evaluation result-the result is shown by histogram-as shown in Figure 7, the red column represents the comprehensive evaluation result.
The blasting effect evaluation of this bench blast is between fair and good. The blasting cost evaluation is similar to the blasting effect evaluation. The blasting quality evaluation is between good and very good, but the blasting safety evaluation is poor. So, this blast design is not optimized. The next blast design should focus on burden and row spacing, etc., to control the blasting fly-rock and vibration.

Conclusions
Blasting effect evaluation is very necessary for blast design. In this paper, a method of blasting effect evaluation was proposed and an evaluation system was developed. The traditional evaluation methods are tedious and time-consuming. However, the developed evaluation system integrates several techniques, such as image analysis and parameter classification. All the selected parameters can be obtained easily to achieve rapid evaluation. The evaluation results are divided as subclass evaluation results and comprehensive evaluation results. The system has been used in practical bench blast engineering. The results are easier for deciders to obtain the advantages and disadvantages of blast design and provide the effective information for the optimization of the next blast design. The comprehensive evaluation system of blasting effect developed is not only convenient, accurate and objective, but also beneficial to enrich and perfect the fine management system of digital blasting management systems. The use of this software greatly reduces the burden of engineers and technicians and realizes the rapid statistics and evaluation of blasting effect.