Optimization of Genetic Algorithm through Use of Back Propagation Neural Network in Forecasting Smooth Wall Blasting Parameters

: With the continuous development in drilling and blasting technology, smooth wall blasting (SWB) has been widely applied in tunnel construction to ensure the smoothness of tunnel proﬁle, diminish overbreak and underbreak, and preserve the tunnel’s interior design shape. However, the complexity of the actual engineering environment and the deﬁciency of current optimization theories have posed certain challenges to the optimization of SWB parameters under arbitrary geological conditions, on the premise that certain control targets are satisﬁed. Against the above issue, a genetic algorithm (GA) and back propagation (BP) neural network-based computational model for SWB design parameter optimization is proposed. This computational model can comprehensively reﬂect the relation among geological conditions, design parameters, and results by training and testing the 285 collected sets of test data samples at different conditions. Moreover, it automatically searches optimal blasting design parameters through the control of SWB targets to acquire the optimal design parameters based on speciﬁc geological conditions of surrounding rocks and under the speciﬁed control targets. When the optimization algorithm is compared with other current optimization algorithms, it is shown that this algorithm has certain computational superiority over the existing models. When the optimized results are applied in practical engineering, it is shown that in overall consideration of the geological conditions, control targets, and other inﬂuencing factors, the proposed GA_BP-based model for SWB parameter optimization has high feasibility and reliability, and that its usage can be generalized to analogous tunneling works.


Introduction
Drilling and blasting method is the main technique for tunnel construction [1]. With the progress in construction techniques, the SWB technique has been widely accepted in drilling and blasting operations. SWB can effectively control the blasting effect, and hence diminish overbreak and underbreak, keep the tunnel's rocky wall smooth, and further maintain the stability of surrounding rocks, lessen the supporting workload, reduce the support materials and engineering costs required, and shorten the construction period [2]. The core issue of SWB lies in the control of overbreak and underbreak [3]. Mahtab [4] believes that the combination of traditional blasting methods with simulation technology can assist in the further evaluation of overbreak and underbreak through the tunnel while the tunnel is drilled and blasted forward; the research outcomes can predict the total

Methods
The traditional BP neural network is prone to parameter underfitting due to improper selection of initial parameters while training with small and medium samples. To address this problem, the genetic algorithm (GA) is combined with the BP neural network. The preferable weights and thresholds of the initial network are obtained by GA, and thus the fitting precision of the BP neural network is improved [24][25][26][27].

Genetic Algorithm (GA)
The implementation of GA includes the following 5 steps: [28] (a) Population initialization.
Individuals are encoded by the real coding method. Each individual is a real string composed of 4 components: weight of connection between the input layer and the hidden layer, threshold of the hidden layer, weight of connection between the hidden layer and the output layer, and threshold of the output layer. The individuals comprise all weights and thresholds of the neural network. Provided that the network structure is known, a network with a definite mapping structure, number of nodes, weights, and thresholds can be constructed.
According to the initial values of the BP neural network obtained by the individuals, the system output is predicted after training the BP neural network using the training data, and the absolute value of the error and the variance E between the predicted output and the expected output as individual fitness F is taken, as calculated by the Equation (1): where n is the number of the network's output nodes; y i is the expected output of the ith node of the BP neural network; o i is the predicted output of the ith node; k is the coefficient for normalization; in this paper, k = 1.
The selection operation in GA is based on the selection strategy of fitness proportion. The selection probability, p i , of each individual i is: where F i is the fitness value of individual i. As it is preferred that fitness be as small as possible, the fitness value shall be inverted prior to individual selection. k is the coefficient with the same value as Formula (1); N is the number of individuals in the population.
As individuals are encoded by real coding, the crossover operation is performed by a real number crossover method. The crossover operation on the kth chromosome a k with the lth chromosome a l at position j is performed by the following method: where: b is a random number within (0, 1).
The jth gene of the ith individual, a ij , is selected to undergo mutation, as operated by the following method: a ij = a ij + a ij − a max × f (g) r > 0.5 a ij = a ij + a min − a ij × f (g) r ≤ 0.5 (5) Mathematics 2022, 10, 1271 where a max is the upper bound to gene a ij ; a min is the lower bound to gene a ij ; f (g) = r 2 (1 − g/G max ) 2 ; r 2 is a random number; g is the count up to the current iteration; G max is the maximum evolution count; r is a random number used to judge the mutation operation within (0, 1), which is automatically generated when selecting a ij .

BP Neural Network
The BP neural network is a multilayer feedforward neural network that can be regarded as a nonlinear function, whose independent and dependent variables are the network's input value and predicted value, respectively. When the number of input nodes is n and the number of output nodes is m, the BP neural network expresses the function mapping relation from the n independent variables to the m dependent variables. The BP neural network shall be trained prior to prediction so that it is endowed with associative memory and predictive ability. The training process of BP neural network includes the following 7 steps: [29] (a) Network initialization. According to the system's input and output sequences (X, Y), determine the number of nodes, n, at the network's input layer, the number of nodes, l, at the hidden layer, the number of nodes, m, at the output layer, initialize the weights w ij and w jk of connections between neurons at the input layer, hidden layer, and the output layer, respectively, and the threshold, a, of the hidden layer and the threshold, b, of the output layer, and give the learning rate and the neuron excitation function.
(b) Hidden layer output calculation.
According to the input variable X, weight w ij of connection between the input layer and the hidden layer, and threshold a of the hidden layer, calculate the hidden layer output H.
where l is the number of nodes at the hidden layer; f is the excitation function of the hidden layer, which can be expressed in many ways. Here, the excitation function is selected as: (c) Output layer output calculation.
According to the hidden layer output H and the weight of connection w ij and threshold b, calculate the predicted output O of the BP neural network.
According to the network's predicted output O and expected output Y, calculate the network's predicted error e.
According to the network's predicted error e, update the weights w ij and w jk of network connections.
According to the network's predicted error e, update the thresholds a and b of network nodes.
(g) Judge whether the algorithm iterations come to an end; if not, return to Step (b).

GA_BP Neural Network
The parameter fitting calculation of the GA_BP neural network falls into two components: the BP neural network and GA optimization. The calculation flowchart is shown in Figure 1.
According to the network's predicted error e, update the weights wij and wjk of network connections.
where η is the learning rate.
According to the network's predicted error e, update the thresholds a and b of network nodes.

GA_BP Neural Network
The parameter fitting calculation of the GA_BP neural network falls into two components: the BP neural network and GA optimization. The calculation flowchart is shown in Figure 1.  Where the BP neural network is composed of two parts: One is the part of the BP neural network structure determination, namely to determine the results of the BP neural network, and thereby the length of GA individuals, according to the number of input and output parameters of the fitting function; the other is the part of the BP neural network parameter fitting, which is responsible for fitting the neural network's input and output parameters, determining the nonlinear relations between parameters, and predicting the function output after network training. The GA optimization component uses the weights and thresholds of the genetic optimization algorithm of the BP neural network. Each individual in the population comprises all weights and thresholds of a network. Individuals' fitness is calculated via the fitness function, in order for the genetic algorithm to find the individual with the optimal fitness value through selection, crossover, and mutation operations, and Mathematics 2022, 10, 1271 6 of 21 thus to optimize the initial weights and thresholds of the BP neural network. The concrete procedure of GA_BP neural network goes as follows: (a) Determine the input and output parameters of GA_BP neural network; (b) initialize the weights and thresholds between the initial parameters of the BP neural network; (c) optimize the above weights and thresholds by GA and select the optimal ones; (d) fit the input and output parameters of GA_BP neural network; (e) error test as to whether the requirement is met, if yes end the calculation, or else return to Step (c).

Parameter Optimization
The calculation for the SWB parameter optimization is to determine the nonlinear relations between input parameters and output parameters by fitting all these parameters via the GA_BP neural network, and to implement the prediction of the output results under the condition of input parameters; next, through the control of the output results, it automates the optimization model to search for the optimal solutions among the input parameters. The calculation process falls into two parts: One is the fitting of GA_BP neural network parameters, and the other is the calculation for parameter optimization. The concrete calculation model is shown in Figure 2.
parameter fitting, which is responsible for fitting the neural network's input and output parameters, determining the nonlinear relations between parameters, and predicting the function output after network training. The GA optimization component uses the weights and thresholds of the genetic optimization algorithm of the BP neural network. Each individual in the population comprises all weights and thresholds of a network. Individuals' fitness is calculated via the fitness function, in order for the genetic algorithm to find the individual with the optimal fitness value through selection, crossover, and mutation operations, and thus to optimize the initial weights and thresholds of the BP neural network. The concrete procedure of GA_BP neural network goes as follows: (a) Determine the input and output parameters of GA_BP neural network; (b) initialize the weights and thresholds between the initial parameters of the BP neural network; (c) optimize the above weights and thresholds by GA and select the optimal ones; (d) fit the input and output parameters of GA_BP neural network; (e) error test as to whether the requirement is met, if yes end the calculation, or else return to Step (c).

Parameter Optimization
The calculation for the SWB parameter optimization is to determine the nonlinear relations between input parameters and output parameters by fitting all these parameters via the GA_BP neural network, and to implement the prediction of the output results under the condition of input parameters; next, through the control of the output results, it automates the optimization model to search for the optimal solutions among the input parameters. The calculation process falls into two parts: One is the fitting of GA_BP neural network parameters, and the other is the calculation for parameter optimization. The concrete calculation model is shown in Figure 2. The concrete steps of calculation for the SWB parameter optimization go as follows:  Initialize the GA_BP neural network and determine the optimal neural network parameters;  Fit the input parameters and output parameters of SWB among the training samples of the GA_BP neural network; The concrete steps of calculation for the SWB parameter optimization go as follows: • Initialize the GA_BP neural network and determine the optimal neural network parameters; • Fit the input parameters and output parameters of SWB among the training samples of the GA_BP neural network; • After fitting the parameters of the GA_BP neural network, conduct an error analysis into the GA_BP neural network using the test samples. If the requirement is met, go to the next step, else return to the first step and modify the basic parameters of the neural network; • When the fitting results meet the requirement of test error, the nonlinear relations between the input and output parameters are reflected, and thus the prediction of results is implemented; • According to the purpose of practical engineering and the parameter design requirement (input parameter control), determine the calculation parameters when the event counter P = 1; • Try figuring out the SWB calculation results under the feasible condition when the event counter P = 1, using the nonlinear relations reflected in the GA_BP neural network.
• After going through the case when the event counter P = 1, determine the calculation parameters when the event counter P = 2, and return to the previous step for calculation until covering all event counters for P = P + 1; • Control the calculation results and thus automatically search for the optimal values of them under this trial condition, so as to derive the optimal design parameters by inversion and implement the optimization of SWB design parameters.

Parameters
In practical tunneling, SWB involves numerous data. Taking the parameters into consideration of the calculation model for SWB parameter optimization is bound to the problems of data redundancy and computational complexity. Therefore, it is necessary to screen the important parameters in SWB and perform the corresponding simplification of the concrete fitting parameters. In general, the SWB effect is subject mainly to geological conditions and blasting parameters. Accordingly, the selected parameters should include factors in three aspects: geological conditions, blasting design parameters, and blasting result parameters. Considering the practical engineering application, and for the convenience of the uniform measurement of test data, the SWB parameters are simplified with overall consideration of the blasting design standard and scholars' research outcomes [8,30,31]. Figure 3 is a schematic diagram of the SWB design parameters. Figure 4 is a schematic diagram for the calculation of average linear overbreak and underbreak. results is implemented;  According to the purpose of practical engineering and the parameter design requirement (input parameter control), determine the calculation parameters when the event counter P = 1;  Try figuring out the SWB calculation results under the feasible condition when the event counter P = 1, using the nonlinear relations reflected in the GA_BP neural network.  After going through the case when the event counter P = 1, determine the calculation parameters when the event counter P = 2, and return to the previous step for calculation until covering all event counters for P = P + 1;  Control the calculation results and thus automatically search for the optimal values of them under this trial condition, so as to derive the optimal design parameters by inversion and implement the optimization of SWB design parameters.

Parameters
In practical tunneling, SWB involves numerous data. Taking the parameters into consideration of the calculation model for SWB parameter optimization is bound to the problems of data redundancy and computational complexity. Therefore, it is necessary to screen the important parameters in SWB and perform the corresponding simplification of the concrete fitting parameters. In general, the SWB effect is subject mainly to geological conditions and blasting parameters. Accordingly, the selected parameters should include factors in three aspects: geological conditions, blasting design parameters, and blasting result parameters. Considering the practical engineering application, and for the convenience of the uniform measurement of test data, the SWB parameters are simplified with overall consideration of the blasting design standard and scholars' research outcomes [8,30,31]. Figure 3 is a schematic diagram of the SWB design parameters. Figure 4 is a schematic diagram for the calculation of average linear overbreak and underbreak.   between the input and output parameters are reflected, and thus the prediction of results is implemented;  According to the purpose of practical engineering and the parameter design requirement (input parameter control), determine the calculation parameters when the event counter P = 1;  Try figuring out the SWB calculation results under the feasible condition when the event counter P = 1, using the nonlinear relations reflected in the GA_BP neural network.  After going through the case when the event counter P = 1, determine the calculation parameters when the event counter P = 2, and return to the previous step for calculation until covering all event counters for P = P + 1;  Control the calculation results and thus automatically search for the optimal values of them under this trial condition, so as to derive the optimal design parameters by inversion and implement the optimization of SWB design parameters.

Parameters
In practical tunneling, SWB involves numerous data. Taking the parameters into consideration of the calculation model for SWB parameter optimization is bound to the problems of data redundancy and computational complexity. Therefore, it is necessary to screen the important parameters in SWB and perform the corresponding simplification of the concrete fitting parameters. In general, the SWB effect is subject mainly to geological conditions and blasting parameters. Accordingly, the selected parameters should include factors in three aspects: geological conditions, blasting design parameters, and blasting result parameters. Considering the practical engineering application, and for the convenience of the uniform measurement of test data, the SWB parameters are simplified with overall consideration of the blasting design standard and scholars' research outcomes [8,30,31]. Figure 3 is a schematic diagram of the SWB design parameters. Figure 4 is a schematic diagram for the calculation of average linear overbreak and underbreak.   Average overbreak and underbreak are calculated by the formulae: for average linear overbreak, and: S q = S q1 + S q2 + . . . + S qn I q1 + I q2 + . . . + I qn for average linear underbreak, where I c and I q correspond, respectively, to the arclength corresponding to each overbreak area and underbreak area.
For optimization calculation, the input parameters are the ones that affect the blasting results, whereas the output parameters are the blasting results. Moreover, in practical engineering, due to the influences of geological conditions, boring equipment, and tunneling purpose, the parameters already have fixed values or designed values and cannot or need not be optimized further. Therefore, the core purpose of SWB parameter optimization is to search for the optimal solutions of spacing between lines of least resistance W, auxiliary hole pitch E b , and ambient hole pitch E a , under the condition of controlling the average linear overbreak and average linear underbreak.

Determination of GA_BP Neural Network Topology and Basic Parameters
Through the above description, the number of input parameters of the GA_BP neural network can be determined to be 11, whereas the number of output parameters is 2. The neural network adopts a three-layer topological structure. The excitation function is selected as a sigmoid function.
Generally, the number of nodes in the hidden layer is calculated by the following empirical formulas [29]: where n is the number of input layer nodes, l is the number of hidden layer nodes, m is the number of output layer nodes, and a is a constant between 0 and 10.
In this paper, the number of nodes in the input layer is n = 11 and the number of nodes in the output layer is m = 2. Therefore, according to the above formula, it is considered that the value range of l is between 4 and 20. Figures 5 and 6 show the convergence speed during sample training and the error rate after sample fitting under different hidden layer nodes during trial calculation. Considering the results of Figures 1 and 2, it is most reasonable to set the number of nodes of the hidden layer as l = 13.    Relevant references [32,33], are selected to determine the network's initial basic parameters: The number of nodes at the hidden layer of GA_BP neural network is 13, with 11 × 13 + 2 × 13 = 169 weights and 13 + 2 = 15 thresholds; individual encoding length in GA is 169 + 15 = 184, population size is 20, evolution count is 100, crossover probability is 0.94, and mutation probability is 0.2. The final GA_BP neural network topology is shown in Figure 7. Relevant references [32,33], are selected to determine the network's initial basic parameters: The number of nodes at the hidden layer of GA_BP neural network is 13, with 11 × 13 + 2 × 13 = 169 weights and 13 + 2 = 15 thresholds; individual encoding length in GA is 169 + 15 = 184, population size is 20, evolution count is 100, crossover probability is 0.94, and mutation probability is 0.2. The final GA_BP neural network topology is shown in Figure 7.

Control Targets of SWB Parameter Optimization
The purpose of parameter optimization is to search for the optima involves numerous parameters, and, without a unique evaluation index results, the evaluation is a multitask and multipurpose problem. At the cur multitask and multipurpose optimization using neural network is implem by two methods: One is to figure out the multiple goals into mutually inde solutions by the method of Pareto solutions, and the other is to turn the

Control Targets of SWB Parameter Optimization
The purpose of parameter optimization is to search for the optimal results. SWB involves numerous parameters, and, without a unique evaluation index of the blasting results, the evaluation is a multitask and multipurpose problem. At the current stage, the multitask and multipurpose optimization using neural network is implemented mainly by two methods: One is to figure out the multiple goals into mutually independent target solutions by the method of Pareto solutions, and the other is to turn the multiple goals into single goals by some calculation model. Combining the practical engineering, the two target parameters, the average linear overbreak and average linear underbreak, are difficult to become mutually independent target solutions. Therefore, the latter approach can be more effective. Moreover, the optimal solutions for the results of SWB typically need to meet two requirements: One is that the tunnel profile shall be as smooth as possible; the other is that the contour line of the practically blasted tunnel shall be as designed as possible. Therefore, reflected in the calculation model for optimization, the control targets of the output parameters are:

•
The minimum of "average linear overbreak S c + average linear underbreak S q "; • The minimum of average linear overbreak S c .

Project Overview
This paper relies on the East Tianshan tunneling project for the Barkol-Hami (firstclass) highway G575, which is located in the area of Hami, Xinjiang. The two-way separated superlong tunnel covers a full length of 11.767 km of the sinistral tunnel and a full length of 11.776 km of the dextral tunnel; the maximum burial depth is 1225 m and the average burial depth is 706 m. Near the tunnel, one can find the typical western mountainous area with obvious underwater seepage. Affected by the burial depth and seepage, the overall strength of the tunnel's surrounding rocks is low, and some construction sections are even vulnerable to the gushing of a large amount of underground water after being tunneled. The schematic diagram of tunnel construction is shown in Figure 8, and the schematic diagram of the project overview is shown in Figure 9.

Sample Collection
While 20 sets of SWB test data are collected from this project, 265 sets of the SWB test data are also collected from other tunneling projects. Figure 10 provides the sources of other SWB data. The input and output parameters of SWB are fitted with the 265 sets of data from other projects as the training samples of the GA_BP neural network and the 20

Sample Collection
While 20 sets of SWB test data are collected from this project, 265 sets of the SWB test data are also collected from other tunneling projects. Figure 10 provides the sources of other SWB data. The input and output parameters of SWB are fitted with the 265 sets of data from other projects as the training samples of the GA_BP neural network and the 20

Sample Collection
While 20 sets of SWB test data are collected from this project, 265 sets of the SWB test data are also collected from other tunneling projects. Figure 10 provides the sources of other SWB data. The input and output parameters of SWB are fitted with the 265 sets of data from other projects as the training samples of the GA_BP neural network and the 20 sets of field test data as the test samples.

Sample Collection
While 20 sets of SWB test data are collected from this project, 265 sets of the SWB test data are also collected from other tunneling projects. Figure 10 provides the sources of other SWB data. The input and output parameters of SWB are fitted with the 265 sets of data from other projects as the training samples of the GA_BP neural network and the 20 sets of field test data as the test samples.  According to the BQ method in the "Standard for Engineering Classification of Rock Masses" (GB 50218-94) of China and the BQ-RMR relation in (17), the parameter surrounding rock rank pertains to geological conditions among the sample data. Limited by space, this paper only displays the surrounding rock mass ranks under preliminary investigation in the Xinjiang-based project, as shown in Table 1. RMR =(BQ − 80.786)/6.0943 (17)  Still limited by space, only the 20 sets of blasting test sample data in the Xinjiang-based project for testing the fitting precision of the GA_BP neural network are presented in Table 2, where the average overbreak is expressed in positive values and the average underbreak is expressed in negative values.

Sample Training Test
GA_BP neural network is used to train and test the above-collected samples. Figure 11 presents the fitness curves of the GA-optimized initial parameters of the BP neural network. Figure 12 presents the overall prediction error curves at distinct iteration counts when GA_BP neural network is fitting the input and output parameters of SWB. From Figures 11 and 12, both the fitting precision and error can meet the practical requirements when the GA_BP neural network is fitting the input and output parameters of SWB, without the problem of underfitting or overfitting, and the nonlinear relations between input and output parameters of SWB can be achieved satisfactorily, hence the prediction of SWB results can be implemented. input and output parameters of SWB can be achieved satisfactorily, hence the prediction of SWB results can be implemented. Figure 11. GA optimized fitness curves. Figure 11. GA optimized fitness curves.

The Accuracy of the Prediction in Both Training and Testing Data
Tables 3 and 4 list the measured value, predicted value, and error rates for 20 sets of training data and 20 sets of testing data. It can be considered that, whether it is training data or test data, the error between the predicted value and the measured value is within 10%. Considering drilling work, this 10% error can be neglected. In general, the accuracy of the prediction meets the actual engineering requirements.

The Accuracy of the Prediction in Both Training and Testing Data
Tables 3 and 4 list the measured value, predicted value, and error rates for 20 sets of training data and 20 sets of testing data. It can be considered that, whether it is training data or test data, the error between the predicted value and the measured value is within 10%. Considering drilling work, this 10% error can be neglected. In general, the accuracy of the prediction meets the actual engineering requirements.

Parameter Input for Optimization Calculation
After the prediction of blasting results is implemented using the GA_BP neural network to complete SWB parameter fitting, the SWB design parameters can be optimized by controlling the predicted blasting results. Relying on the practical engineering project, this paper has designed five sets of different input parameters, among which the fixed parameters are shown in Table 5, and the ranges of values of the parameters to be optimized are shown in Table 6. As the values of W, E b , and E a are of little significance below the mm level in practical engineering, the parameters are optimized only at the cm level.  Table 7 displays the optimization calculation results at different design numbers and the predicted overbreak and underbreak values at the corresponding optimal values. The optimization calculation results in Table 6 are applied in the engineering practice, and the overbreak and underbreak values in real blasting are compared with the predicted overbreak and underbreak values with the results shown in Table 8. The actual effects of engineering blasting are shown in Figures 13 and 14.    The optimization calculation results in Table 6 are applied in the engineering practice, and the overbreak and underbreak values in real blasting are compared with the predicted overbreak and underbreak values with the results shown in Table 8. The actual effects of engineering blasting are shown in Figures 13 and 14.   From the comparison between the predicted results of optimization calculation and the practical engineering results in Table 8, it can be found that the relative errors between the predicted overbreak/underbreak values and the actual values in engineering are small and within 10%, demonstrating that the method of optimization calculation is scientific and effective and that the predicted results of optimization are realistic and From the comparison between the predicted results of optimization calculation and the practical engineering results in Table 8, it can be found that the relative errors between the predicted overbreak/underbreak values and the actual values in engineering are small and within 10%, demonstrating that the method of optimization calculation is scientific and effective and that the predicted results of optimization are realistic and scientific. From Figures 13 and 14, it can be found that applying the optimized results of parameters in practical engineering can effectively enhance the accuracy of SWB during blasting tunneling and further ensure the security and stability of the surrounding rocks.

Optimization Calculation Results and Engineering Application
The accomplishment in engineering application demonstrates that the GA_BP neural network-based calculation model for SWB parameter optimization can be applied to not only effectively implement the prediction of SWB effect but also upgrade the SWB effect and the overall safety in tunnel construction.

Discussion of the Superiority of GA_BP Neural Network Algorithm
To verify the superiority of the GA_BP neural network algorithm in data processing, a comparison is made between the predicted results under the test samples after the fitting calculation of its training samples and the actual values of blasting, the predicted values of the BP neural network, and the predicted values by the GA_ISVR algorithm. The predicted values of the three are shown in Table 9, and the algorithm average training errors are shown in Figures 15 and 16. a comparison is made between the predicted results under the test s ting calculation of its training samples and the actual values of blas values of the BP neural network, and the predicted values by the G The predicted values of the three are shown in Table 9, and the algor ing errors are shown in Figures 15 and 16.    From Figures 15 and 16 and Table 9, one can find that the GA_BP neural network has a significant upgrade over the BP neural network in the accuracy of overbreak or underbreak prediction after the GA is introduced, and that the error of SWB parameter fitting by the GA_BP neural network is significantly smaller than by the GA_ISVR algorithm after the sample size is enlarged. Evidently, the GA_BP neural network can be applied for excellent fitting and prediction of SWB parameters. In practical engineering, the GA_BP neural network outperforms the BP neural network and GA_ISVR algorithm in predicting SWB results.
In addition, when discussing the superiority of the GA-BP algorithm, there is still a key problem, that is, the inaccuracy and uncertainty of the data. Therefore, it is also necessary to verify and analyze the accuracy and certainty of the data, so as to avoid calculation errors or analysis errors. In future work, based on References [34][35][36], the accuracy and certainty of the data will be investigated through the "prediction system based on fuzzy logic".

Discussion of Optimization Results of SWB Parameters
First, we define an indicator K called relative change rate: where a and b represent the upper and lower bounds of the change interval, respectively. K represents a percentage. According to Table 6, Scheme 1, Scheme 2 and Scheme 4 are selected. The sensitivity calculation of the optimization parameter index is shown in Table 10, and the analysis results are shown in Table 11.  After normalizing the optimization parameters "auxiliary hole pitch E b , spacing between lines of least resistance W and ambient hole pitch E a " according to the "change rate of 10%", it can be found that in each scheme, the relative change of over and underbreak caused by the change of "auxiliary hole pitch E b " is the largest, its maximum value reaches 63.35%, and the average change rate of the three schemes is 31.32%. In contrast, the maximum relative change rate of overbreak and underbreak caused by the change of "spacing between lines of least resistance W" is 43.96%, and the average change rate of the three schemes is 23.00%. The maximum relative change rate of overbreak and underbreak caused by the change of "ambient hole pitch E a " is 43.61%, and the average change rate of the three schemes is 23.67%. It can be observed that for the optimization parameter index, "auxiliary hole pitch E b " is the most sensitive factor affecting the value of overbreak and underbreak.
In addition, according to Table 11, the overall relative change rate of Scheme 1 is 25.58%, that of Scheme 2 is 19.47%, and that of Scheme 4 is 33.42% when all optimization parameters change by 10%. Referring to Table 6, the σ c of rock in Scheme 2 is the largest and its relative change rate is the smallest, while the σ c of rock in Scheme 4 is the smallest and its relative change rate is the largest. Therefore, it can be considered that this change law may be related to rock properties, which need further study and analysis.

Conclusions
The combination of genetic algorithm and the BP neural network can improve the generalization degree and calculation accuracy of the prediction model, so as to solve the problem of insufficient accuracy caused by insufficient data quantity. Based on this, in this work, a new algorithm based on the GA_BP neural network is proposed, which is applied to the prediction and optimization of the tunnel SWB parameters. By training the input data (geological conditions and SWB parameters) and output data (overbreak/underbreak), the algorithm model builds the nonlinear relationship and realizes the prediction of the SWB effect. Moreover, based on the control of the prediction, the optimization of the design parameters of the SWB is realized.
In addition, through the analysis of the optimization of tunnel SWB parameters, it is believed that the "E b (auxiliary hole pitch)" has the greatest impact on overbreak and underbreak. Therefore, this index should be given priority when determining parameters.
The application results demonstrate that the algorithm model is effective and feasible to predict and optimize the parameters of SWB, and the results can meet the requirements of practical engineering.