Evaluation of the Total Quality of Tunnel Contour Using Projection Pursuit Dynamic Cluster Method

School of Civil Engineering and Architecture, Zhejiang University of Science and Technology, Hangzhou 310023, China Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Beijing 100083, China College of Water Resource and Hydropower, Sichuan University, Chengdu 610065, China China Railway No. 2 Engineering Group Co., Ltd., Chengdu 610031, China


Introduction
It is well known that the quality of tunnel contour after smooth blasting always correlates with the tunnelling safety and construction cost, as well as the penetration rate [1][2][3][4][5]. Consequently, appropriate evaluation of the total quality of tunnel contour is of vital importance to the advancing of drilling-blasting tunnelling. In general, the contour quality depends on the complexity of geological conditions and the drilling-blasting (D&B) operations which might cause adverse effects on the integrity and stability of the surrounding rock mass [2,6,7]. In view of applicability, the quality of tunnel contour is now mainly determined by the overbreak in quality rating [8,9].
Many attempts have been made to predict the scale of over/underbreak and the range of blast-induced damage, for optimizing the blast design/operation and improvement of the contour quality in D&B tunnelling [10][11][12][13][14][15]. Because of the uncertainty of site geology and the diversity of influential factors of D&B operation, however, most of the investigations on evaluation of contour quality mainly focused on the overbreaks at specific locations of the contour. As of now, very few efforts were devoted to quantify the contour condition of the entire tunnel after a round of D&B operation. Yangkyun Kim [16] proposed an index named Tunnel Contour Quality Index (TCI) to quantify the general geometric features of the tunnel contour after smooth blasting, considering the average overbreak depth and average ratio of actual contour length to the planned contour length (RCL) in two cross-sections of a round and the average of absolute value of difference in overbreak depth between two neighboring sections in longitudinal direction. Elisa Costamagna et al. [2] conducted an assessment of overexcavation and TCI and then compared it with the Q-system values. ey suggested that the TCI index can be selected as a quantitative and applicable tool to evaluate the tunnel contour quality to some extent. Following their previous work, Yangkyun Kim and Amund Bruland [8] also established another index named Tunnel Contour Quality Index for Construction Cost (TCIC) to evaluate the contour quality and the expenditure. Although the indices of TCI and TCIC offered a way to quantify the tunnel contour quality after smooth blasting, only the overbreaks in some specific linear directions in selected cross/longitudinal sections were taken into account. In fact, the total quality of tunnel contour not only depends on the depth of over/ underbreaks at some specific positions but also is determined by the total scale features in terms of area and volume of over/underexcavations. Koopialipoor et al. [17] analyzed the overbreak using both the artificial neural network (ANN) and hybrid genetic algorithm (GA)-ANN and found that the GA-ANN was more reliable to predict the overbreak. Considering the multiplicity of causing factors of overbreak and the complexity of interactions among these factors, Koopialipoor et al. [18] attempted to develop a new overbreak prediction method by using artificial neural network. Navarro et al. [19] reported an engineering tool to predict high risk of overbreak zones based on the MWD data. Day [20] indicated that the brittle overbreak in the heterogeneous rock mass containing hydrothermal veins and breccia can behave quite differently than the homogeneous rock mass, resulting in the traditional prediction tolls of depth of brittle overbreak which will be ineffective. Jang et al. [21] defined an overbreak resistance factor to develop an empirical method for overbreak estimation. Foderà et al. [22] proposed an operative methodology to estimate overbreak volumes and to distinguish the technical overbreak based on the tunnel laser scanning surveys and high-resolution images of the excavated surfaces, and the influence of influential parameters is assessed by the multiple regression analysis. Delonca and Vallejos [23] presented a generalized failure criterion which includes the scale effect for prediction of the stress-induced overbreak for all excavation diameters. It could be found that great attention has been paid to the estimation of overbreak depth and damage zone prediction. However, the contour quality does not only depend on the depth of overbreak at some specific locations.
In engineering applications, the majority of quality rating methods of tunnel contour after smooth blasting are based on empirical or semiempirical formulas, engineering analogies, and field tests [2,24,25]. Recently, some data mining approaches have been tried to estimate the overbreak and damage zone of surrounding rocks, including adaptive finite element method, fractal theory, image-processing technique, wavelet analysis, model testing, artificial neural network, and multiple regression analysis. But most of these methods are generally qualitative and based on single evaluation index, which might be only applicable under some specific conditions [7,[26][27][28][29]. A comprehensive index system including multiple indices which can be used to quantitatively characterize the total quality of contour is still unavailable. In the reality of tunnel construction using D&B method, different departments might use various indices to judge the quality of tunnel contour according to their concerns. For example, the over/underbreak depth is often selected as the index to estimate the contour quality in the Chinese code for construction on tunnel of railway [30] and the code for construction technical specifications for highway tunnels [31]. Other factors including the area and volume of over/underbreaks are neglected in practical applications. However, effects of these factors on the quality estimation of tunnel contour cannot be ignored in many cases. For rational understanding and assessing the construction safety and time, therefore, development of a comprehensive evaluation method for the quantitative assessment of the total quality of a tunnel contour after smooth blasting is very promising in D&B tunnelling.
is paper targets developing a quantitative evaluation approach with multiple indices to fully understand the total quality of tunnel contour based on the projection pursuit dynamic cluster (PPDC) method. Both of the engineering practical experiences, site measurements, technical standards, experts' knowledge, and relevant literatures are included in the proposed evaluation approach. e influences of over/underexcavations in 3 dimensions including the depth, area, and volume are all counted into the assessment to fully understand the total quality of tunnel contour after smooth blasting.

Selection of Evaluation Indices.
e indices representing the quality of tunnel contour might differ from each other, according to the primary purpose of the evaluation. As is known, it is often impossible to involve all factors into the analysis of complicated problem. Appropriate determination of the evaluation indices virtually plays a very important role in the reliability of the assessment. Identification of the proper indices is the first but also the most important step. In this work, determining the indices follows the following principles: high correlation to the contour quality, being applicable to be obtained, considering both qualitative and quantitative factors, and being easy to use.

Structure of the Multiple Indices
System. In order to fully understand the total quality of the contour, the multiple indices system should be established. Following the concept of the analytic hierarchy process (AHP) method for complicated issues [32], a hierarchical structure model of the proposed multiple indices system for evaluating the tunnel contour quality is utilized, as shown in Figure 1. It can clearly show the influential factors involved in this investigation and correlations among various factors.

Data Collection
e data means the values of potential factors that may affect the total quality of tunnel contour. Identification of these data is performed by using site investigation, expert consultation, and experimental study. According to their utilization frequencies in existing cases, these potential indices representing the contour quality can be grouped into two types: the qualitative indices and quantitative indices. For a qualitative index, grading its effect on the contour quality and accordingly giving a score as its values are employed to quantify the influence of this index. On the other hand, the quantitative index should be converted into dimensionless form, to reduce the influences of variations of absolute values of these factors.

Identification of Potential Indices.
e potential indices can be obtained from literatures, engineering experiences, expert knowledge, and technical codes in railway, highway, coal mine, hydropower industries, and other underground constructions. In this work, 83 cases of rock tunnelling using D&B method in China, 20 overseas cases, and 16 technical codes are examined and investigated. In this work, we define the utilization rate to represent the frequency of the index used in quality rating. Following the definition in literature [33], the utilization rate is defined as the ratio of the utilization times of the index to the total number of existing methods. e adoption rate of the index utilized in traditional quality rating is statistically analyzed, as shown in  It could be found that the cast factor, blast hole utilization factor, advance per round, and overbreak can be considered as the main indicators representing the contour quality of tunnel smooth blasting in China from different aspects. However, the overbreak is often selected as the major index to indicate the contour quality in other countries. Statistically, 32 potential indices can be identified to represent the tunnel contour quality after smooth blasting, which should be taken into account in the proposed comprehensive quality rating.

Screening of Qualitative
Indices. In this work, whether an index should be counted into the evaluation depends on its utilization rate in existing methods; that is, the adoption rate should be no less than 50%. Qualified qualitative indices can be identified via on-site investigation and expert consultation. According to the importance degree in the evaluation of contour quality, the influence level of qualitative index is grouped into five levels. Specifically, level 1 means very important factor which will be given a score of 5 points. Similarly, level 2 represents relatively important factor (4 points), level 3 indicates that its importance degree is fair (3 points), level 4 is less important factor (2 points), and level 5 is unimportant (1 point). In this work, the qualitative indices include the economic benefits, stability of the surrounding rock mass, influence of fly-rock, rock supporting factors, vibration effect, cracking of retaining wall of surrounding rock, and difficulty in operation of construction method. According to the result of grading and score, the very important indices are the economic benefits, stability of the surrounding rock mass, rock supporting factors, and difficulty in operation of construction method, as shown in Figure 5. Only one index shows relative importance to the contour quality, that is, the fly-rock.

Collection of Quantitative Indices.
e quantitative indices can be identified via on-site investigation, expert consultation, and experimental study. Similarly, a quantitative index is qualified if its utilization rate is no less than 50%. Based on site investigation and expert consultation, 25 quantitative indices are collected, as shown in Table 1.
According to the utilization rate, the indices which are very important (5 points) to the contour quality include overbreak, underbreak, overbreak rate, average overbreak, maximum size of step between two guns, explosive dosage, and the drilling footage per unit area, as shown in Figure 6. Accordingly, the indices with relative importance include average block size, disturbance depth of the surrounding rock, explosive charge ratio, irregularity degree, blasting vibration velocity, advance per round, drilling footage per unit area, and retaining the thickness of surrounding rock damage zone (4 points), as shown in Figure 7.
e experimental studies including field testing and numerical modelling also can provide some valuable information on quantitative indices. e indices investigated at the site include measured cross-sectional area, overbreak area, underbreak area, maximum overbreak, maximum underbreak, and average overbreak, by using a tunnel profilometer.
e first five indices also can be identified based on numerical simulations.

General Method.
e multiple indices system is considerably critical to the reliability of the quality evaluation of tunnel contour after smooth blasting. Appropriate indices system can provide convincing result of contour quality evaluation. e determination of the multiple indices system is based on the following considerations: (1) Selecting a serial of indices that can generally indicate the contour quality. e index used in most evaluation methods can be considered as a qualified one, which will represent the characteristics of the contour quality to some extent. Indices which are adopted frequently in most existing evaluations should be involved in this work.
(2) Ranking the selected indices based on their utilization rate. e more frequently the index is employed, the more important it is. Meanwhile, the contribution of each index to the entire contour quality is quantitatively evaluated by the principal component   analysis (PCA) method and the gray correlation analysis (GCA). In this case, the importance degree of an index can be indicated by the weight. e framework of determining the multiple indices system is shown in Figure 8.

Primary
Screening of the Indices. Up to date, there is still no universal index system for the evaluation of the contour quality. e indices utilized in one method might differ from those of another method. Admittedly, each index has its value but the limitations cannot be ignored. For example, the average overbreak and maximum overbreak which are widely utilized in many evaluations can only represent the contour quality at one specific point, while the over/underbreaks in 2 or 3 dimensions cannot be assessed. In order to fully understand the contour quality in 3 dimensions, other indices representing the total area and volume of over/underbreaks should be counted in. erefore, 21 indices with at least 50% of utilization rate are primarily selected as parts of the multiple indices system. As shown in Table 2

Optimization of the Indices.
In general, a large number of indices are often not applicable to perform the evaluation work due to both the computing efficiency and applicability. In this case, the selected 21 indices should be further optimized according to their importance degree. e coupling method of principal component analysis (PCA) method and gray correlation analysis (GCA) is employed to determine the importance degree.

Optimization of Point Indices.
In this section, X 1 , X 2 , X 3 , and X 4 denote the maximum overbreak, maximum underbreak, average overbreak, and average underbreak, respectively. eir weights are determined by coupling  method of the GCA and PCA. e GCA is often utilized to determine the correlation degree among different factors in the system on the basis of the similarity or dissimilarity of the development trend [34]. e correlations of factors can be defined by the similarity of the geometric shapes of the sequence curves. e correlation degree is significant if the geometric shapes of the curves are approximate. To perform the GCA, the reference sequence can be defined To reduce the effect of deviation of factor unit, the sequence should be nondimensionalized by the mean value processing method. e reference sequence denoted as x 0 ′ can be expressed as Equation (1), while the comparative sequence indicated as x 0 ′ can be described as Equation (2).
e correlation coefficient is defined as in which w is the resolution coefficient, Δ min is the minimum difference between the two sequences, and Δ max is the maximum difference between the two sequences. erefore, the correlation degree can be determined by In this work, the relational order of the point indices is shown in Table 3. e correlation degree for each point index is shown in Table 4.
It could be found that the maximum overbreak X 1 has the strongest correlation with the tunnel contour quality, while the average underbreak X 4 shows the lowest correlation with the contour quality.
On the other hand, the contribution of each index to the total contour quality also can be determined using the PCA method. As is known, the PCA is probably the oldest and best known of the techniques for multivariate analysis [35].  [35,36]. Following Eriksson et al. [37], the PCA model is in which X is the input matrix (n × p) which represents original parameters, T is the score matrix (n × k) which represents the relationship between original parameters, P is the loading matrix (p × k) which represents the contribution of original parameters, and E is the residual matrix (n × p) which represents the uncaptured variance. n is the number of time steps, p is the number of original parameters, and k is the number of principal components. On the basis of the maximum variance criterion, each principal component obtained from the PCA model is irrelevant to another one. e first principal component has the highest variance which also reveals that it contains the most information, while the second principal component will capture the next highest variance which has already removed the information of the first principal component [36]. And the rest can be deduced in the same manner. For the kth principal component for the point index, it can be denoted as follows: PC k � a 1k P intakepointindice + a 2k P dischargepointindice + a 3k P intakepointindice + a 4k P dischargepointindice+...+a pk P vibration .
e first three principal components have the highest variance. It means that most of the information in the original parameter is visualized by the first three principal components. e normalized eigenvectors of the point indices are shown in Table 5, and the corresponding eigenvalues are given in Table 6. It can be seen that the first principle component involves the average overbreak X 3 and the maximum overbreak X 1 , the second principle component is the maximum underbreak X 2 , and the third principle component is the average underbreak X 4 .
Combining both the results of PCA and GCA, the optimized point indices are shown in Table 7. In point level, the maximum overbreak X 1 , average overbreak X 3 , and the maximum underbreak X 2 are selected as the optimum point indices in sequence.

Optimization of Linear
Indices. Y 1 -Y 5 represent the disturbance depth of the surrounding rock, maximum size of step between two guns, irregularity degree, average block size, and advance per round, respectively. According to the results of PCA and GCA, the optimized linear indices are shown in Table 8. e maximum size of step between two guns Y 2 is determined as the optimum linear index.

Optimization of Surface
Indices. W 1 -W 6 denote the measured cross-sectional area, area of overbreak, area of underbreak, rate of overbreak, rate of underbreak, and half cast factor, respectively. e results of optimized surface indices are given in Table 9. Area of overbreak W 2 , measured area W 1 , and rate of overbreak W 4 are selected as the optimum surface indices.

Optimization of Volumetric
Indices. Z 1 -Z 6 indicate the volume of overbreak, volume of underbreak, blast hole utilization factor, explosive dosage, economic benefits, and stability of the surrounding rock mass, respectively. e optimized volumetric indices are shown in Table 10, based on the analysis of PCA and GCA. It could be found that the volume of overbreak Z 1 and blast hole utilization factor Z 3 can be determined as the optimum volumetric indices.

Development of the Multiple Indices System.
On the basis of the two rounds of screening, the optimal indices that will be used to establish the multiple indices system can be figured out clearly. According to the hierarchical structure, the proposed multiple indices system is given in Table 11. It includes four levels under the criterion layer and in total 9 indices representing the contour quality for the index layer.

Contour Quality Evaluation Using the PPDC Method
e proposed multiple indices system has 9 indices to represent the contour quality from different aspects. As is known, it is hard to quantitatively analyze the multipledimensional data. us, the projection pursuit (PP) method is used to reduce the dimensions of the data in this work, and the dynamic clustering (DC) is employed to establish the rating of the contour quality. As a result, the total contour

Projection Pursuit (PP) Method.
For multiple factors evaluation problems, mapping of multivariate data onto low-dimensional manifolds for visual inspection has been considered as a common technique in data analysis. e projection pursuit (PP) algorithm is such a tool to solve the multiple factors evaluation problems [38]. It can convert high-dimensional data into low-dimensional space by mapping in a projection direction [39]. In low-dimensional space, the salient features of the original data can be analyzed conveniently by the human gift for pattern recognition. Specifically, the PP algorithm provides a linear mapping method that uses interpoint distance as well as the variance of the point swarm to pursue optimum projections [38]. e major task of PP analysis is particularly a dimension-     10 Advances in Civil Engineering reducing mapping issue, to determine the optimal projection direction which can provide more valuable information from the original high-dimensional data. e procedure of PP mapping is given as follows: (1) Normalization of the original data. In order to reduce the effect of different magnitudes of various factors, the original data should be first normalized. Suppose j) is the value of jth factor after normalization, n is the number of samples, and p is the number of factors. For factors that are better with bigger values, the normalization formula is given in Equation 7. For others, the original data could be normalized by Equation 8.
(2) Establishing the function of linear projection. Suppose a � a(1), a(2), . . . , a(p) is a p-dimensional unit vector and z (i) is the projection characteristic value of x * (i, j). e linear projection value can be obtained from (3) Optimization of the projection index function. e optimal projection direction can be determined by the maximum value of the projection index, in which the genetic algorithm (GA) could be employed to perform the optimization. e projection index Q(a) is a function of the projection axis a → , which can be defined as follows: where S z denotes the spread of the data which can be calculated by Equation (12) and D z is the local density of the projection values on the projection direction a which can be defined as Equation (13).
in which z(i) is the mean value of z (i) along the projection axis a → and R � S z × 10% which means the cutoff radius. (4) Calculation of the projection value. After the determination of the optimal projection direction a → , the projection value of samples can be calculated by Equation (9). As a result, a characteristic projection value which represents the total quality of contour can be obtained.

General Method.
In order to qualify the quality grade of the tunnel contour using the projection value of the multiple indices, the rating of contour quality should be developed in advance. In this study, dynamic clustering (DC) algorithm is used to quantitatively determine the quality grade with 5 levels. e concept of classification is to partition a set of data described by frequency distributions   [40]. It is a general partitioning algorithm of a set of objects in K clusters, which is based on the definition of a criterion of the best fitting between the partition of a set of individuals and the representation of the clusters of the partition [41]. e partitioning operation is described as follows: (1) Initial partitioning: M selected samples α i (i � 1, 2, . . . , m) are randomly divided into f initial groups, that is, F z (z � 1, 2, . . . , f). (2) Calculating the distance: the mean vector z z of each initial group F z could be calculated by where H z and n z are the sum of vectors and number of vectors, respectively. Consequently, the distance between the value of ith sample α i and mean vector z z can be calculated by (3) Establishing the first group: the sample value α i is modified by Equation (16), and thus the first group can be obtained.
(4) Repeating the calculation steps of (2) and (3) with modified sample values to get the ranges of other groups.

Rating of Each
Index. e rating of each evaluation index in the proposed multiple indices system is performed by DC algorithm and the confidence interval estimation. e importance degree of each index is grouped into 5 levels. e original data of each index and related effects on contour quality are collected from the huge number of existing projects. Taking the average overbreak as an example, the calculated result is shown in Table 12. e up/ down limits of each level could be determined by the significance level.
Accordingly, the rating of average overbreak can be obtained as shown in Table 13.
Similarly, the rating of other 8 indices involved in the proposed multiple indices system could be obtained, as shown in Table 14.

Contour Quality Rating Based on Projection Value.
e comprehensive contour quality can be determined by the eigenvalues of the data set of the multiple indices. Based on Table 14, the optimal projection direction of the multiple indices can be obtained by using the MATLAB code. e computer configuration for performing the calculation is shown in Table 15.
Calculated optimal projection direction a → is a * � (0.3154 0.3663 0.1952 0.2508 0.5209 0.2346 0.3747 0.2091 0.3934), and the computed eigenvalues using the PP method for each level are given in Table 16. As shown in Table 16, the total quality of tunnel contour can be classified into 5 levels, according to the eigenvalues.    12 Advances in Civil Engineering

Project Background.
A high-speed railway tunnel from Chengdu to Chongqing excavated by the drill-blast method is selected as the site of field study. e tunnel is located at the mountainous area in west of Chongqing. During tunnelling, the smooth blasting was conducted to excavate the rock formations. e designed length of the tunnel is up to 5054 m, and the height of the test cross section is 11.08 m while for width it is 14.90 m. According to the Chinese technical standard TB10003 [42], the grade of the surrounding rock is classified to level IV, meaning the integration and stability of the surrounding rock are poor. e main lithology exposed during excavation is argillaceous sandstone with age of Jurassic (J 2 s). e properties of the surrounding rocks are shown in Table 17.

Measured Data.
e values of the multiple indices including the maximum overbreak, average overbreak, maximum underbreak, maximum size of step between two guns, area of overbreak, measured area, overbreak rate, volume of overbreak, and blast hole utilization factor are measured/ collected at 3 cross-sections, that is, at DK247 + 033.5 m, DK247 + 047.5 m, and DK247 + 049.5 m, as shown in Table 18.

Calculated Eigenvalues.
Based on the measured data, the optimal projection of the 9 indices representing the tunnel contour quality can be obtained by using MATLAB programming; that is, a * � (0.1841 0.4406 0.0007 0.0034 0.5045 0.4941 0.5099 0.1154 0.0001). Substituting the obtained optimal projection direction into the linear projection function, that is, Equation (9), the calculated projection values of comprehensive quality of the 3 measure contours are Z * � (0.57 0.81 2.17).

Rating of Contour Quality.
Consequently, the rating of contour quality can be obtained by comparing the calculated eigenvalues with the developed rating (Table 16). e result of quality rating of the 3 measured tunnel contours is shown in Table 19. e actual profiles of the 3 contours are shown in Figures 9-11.
From the site pictures of the contours, it could be found that the rating of total contour quality defined by the proposed multiple indices-based evaluation method is much different from that obtained using single index. If the single index is selected to describe the contour quality, for example, the magnitude of the maximum overbreak or underbreak of  Note. e area deviation is the difference between the measured cross-sectional area and the designed area.  contour at DK247 + 049.5 m is relatively higher than the other 2 contours, which will lead to a poor contour quality. However, the total volume of over/underbreaks in contour at DK247 + 049.5 m is much lower than that of others, and the observed maximum over/underbreak only occurs in very few positions. e relative low value of over/underbreaks of the contour at DK247 + 049.5 m also can be examined by the measured area and volume of overbreak. In fact, the contour at DK247 + 049.5 m shows much better quality than that of other contours because of low volume of overbreak, which will reduce the expenditure on backfill and rock support. Obviously, the result from the proposed contour quality evaluation method is more accordant with the fact of the site, which will benefit the cost-effective decision-making by the site engineers.

Conclusions
e tunnel contour quality after smooth blasting is considerably critical to the estimation of construction cost and time, as well as the safety of tunnelling and the stability of the surrounding rock. Comprehensive understanding of the contour quality in 3 dimensions is of vital importance to reasonably evaluate the construction time and cost. Compared with single index, the proposed multiple indices can better represent the total quality of tunnel contour. is paper developed a multiple indices-based evaluation approach, and thus the projection pursuit dynamic cluster (PPDC) method is utilized to qualify the quality grade of the tunnel contour. Some conclusions can be drawn as follows: (1) Based on statistical analysis of engineering cases, site investigation, literatures, expert knowledge, and the technical codes, 7 qualitative indices and 25 quantitative indices representing the contour quality are effectively figured out. According to the utilization rate, 21 indices are preliminarily selected to fully describe the contour quality in 3 dimensions. (2) To identify the most influential factors affecting the total contour quality, the coupling method of the GCA and PCA is utilized to determine the multiple indices system. e selected 9 indices include the maximum overbreak, average overbreak, maximum underbreak, maximum size of step between two guns, area of overbreak, measured area, rate of overbreak, volume of overbreak, and blast hole utilization factor.
(3) To rationally determine the quality grade of contour, the projection pursuit (PP) method is used to process the high-dimensional data of the multiple indices, and the dynamic cluster (DC) analysis is utilized to establish the rating of total contour quality. e total quality of tunnel contour in 3 dimensions is grouped into 5 levels, according to the eigenvalues. (4) If the measured values of the 9 indices are available, the optimal projection direction and eigenvalues can be obtained by using the PP algorithm. In this case, the rating of contour quality can be determined accordingly. Compared with the traditional method based on overbreak, the result of contour quality using the proposed method is more accordant with the fact of the site. It is believed that the developed PPDC-based evaluation method will offer a promoting alternative to reasonably describe the total quality of tunnel contour after smooth blasting.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper. Acknowledgments e financial supports for this research project by the National Natural Science Foundation of China (nos. 41602308