Abstract
The quality assessment of surrounding rocks has tremendous influences on engineering construction, so its investigation is significant. The uniaxial compressive strength of rocks , the integrity coefficients of rock mass , the conditions of the structural plane , the seepage amount of groundwater , and the intersection angle between the axis line and structural plane are selected as the quality assessment index of surrounding rocks at first. Then, the entropy weight-variable fuzzy sets model is introduced. Second, the relative membership degree matrix about the variable fuzzy sets is established, and the weight coefficients are determined based on the entropy weight method. Finally, the fracture behavior level of surrounding rocks is determined using the mean ranking feature value. The conclusions demonstrate that the results obtained based on the variable fuzzy sets method are consistent with the current specification; its accuracy arrives at 80%. Compared with the traditional way, its assessment process has higher reliability and efficiency. So, it could provide a new alternate way to assess the risk level in civil engineering accurately.
1. Introduction
Many factors affect the stability of surrounding rocks in reality [1, 2], for example, the construction level, timbering opportunity, and the geometric deformation of cross-section [3, 4]. In essence, these factors are relevant to the quality of surrounding rocks [5, 6], so the quality of surrounding rocks during the construction process of the tunnel is adopted to assess the stability of the tunnel. Moreover, it is also used to select the engineering structure parameters and evaluate the economic benefit [7]. So, the quality assessment of surrounding rocks in the tunnel has great significance.
Many researchers in many countries have provided many methods to assess the quality of surrounding rocks [8] in recent years. The fractal theory is supplied by Du et al. [9] and Point et al. [10] to assess the surrounding rocks’ quality based on the fractal characteristics of the rock structure face. The fractal dimension about over and under excavation of rock mass is calculated by Boadu, etc., [11] to seek the relationship between structural features and the category of surrounding rocks. The probability and statistical method for rock mass classification are suggested by Boadu and Long [12] to determine the rock strength and deformation parameters based on the classification index of surrounding rocks. Gao et al. [13] perform an artificial intelligence expert system to classify the tunnel’s surrounding rock level. Then, a fuzzy information analysis model is suggested by Zhou et al. [14] to establish the classification of surrounding rocks in the road tunnel. The Fuzzy comprehensive evaluation method is provided by Bieniawski and Balkema [15] to classify the level of the rock mass. Zhang et al. [16] suggested the hierarchy analytical method to determine the weight distribution of different influential indices about the surrounding rocks. The grey optimization theory about the stability classification of surrounding underground rocks was established by Bai et al. [17] based on the correlated grey analysis. Ruan and Wei [18] and He and Chen [19] showed the BP neural network model to deal with the nonlinearity and uncertainty of the classification of the surrounding rocks. Li et al. [20] solved the actual category of surrounding rocks in the hydraulic tunnel based on the artificial neural network theory. The matter-element model of surrounding rock evaluation in underground engineering is applied by Ye et al. [21] according to the procedure of extension method. The support vector machine theory is used to assess the quality of rock mass by Zhao et al. [22], and its higher accuracy is improved; besides, the distance discriminant analysis method is applied to determine the classification of surrounding rock quality in the tunnel by Gu et al. [23] and Chen and Zhou [24]. Wong et al. apply the efficacy coefficient method to destabilization risk assessment of surrounding rocks in the tunnel [25, 26].
Although the above method has promoted the development of the assessment theory of surrounding rock quality, some shortcomings still exist [27–29]. For example, the calculative process is complex, and the assessment process in many methods is often quantitative or qualitative. To overcome the shortcomings of the above methods, the entropy weight-variable fuzzy sets are introduced to assess the quality level of surrounding rocks in the tunnel. It has many virtues, such as the preciseness of the algorithm and operability in practice, and it can well solve the grading standards, which are interval forms. The method is a tremendous improvement on the traditional fuzzy sets model.
The paper is organized as follows: in Section 2, the study overview is introduced at first; in Section 3, the theory and methodology based on the entropy-weight variable fuzzy sets model are introduced; in Section 4, the assessment model about the surrounding rock quality is established, and the assessment results of the proposed model are compared and discussed; in Section 5, conclusions are drawn.
2. Study Area
The Mandarin Duck intersection tunnel is located between the Mandarin Duck intersection village and Dong Niu Po village in Suozhou city, Shanxi province, China. Its locations are plotted in Figure 1. It belongs to the Shanyin to Pinglu districts of National Highway in Shanxi province. It is a two-way and four-lane separated tunnel. The space between the left and right lines is 40 m. The import mileage at the left line is LK178 + 102; the exit mileage is LK182 + 830, and its total length is 4728 m. The longitudinal slope is 1.931%. The import mileage at the right line is RK178 + 044; its total length is 4787 m, and the longitudinal gradient is 1.98%. The maximum depth of the bottom plate in the tunnel at the left line is 316 m, and one at the right line is 317 m; it belongs to the extralong tunnel. The inclined shaft is placed in the middle of the tunnel. The tunnel passes through the tectonic and denudation strip in the mountain; its altitude height is 1400–2000 m, the cutting depth is 100–300 m, the slope angle is significant, and the mountain peak is steep. The lithology is mainly carbonate rock; the Archean metamorphic rock is beneath the carbonate rock. The different thicknesses of loess are covered at its upper. The maximum depth is 317 m. Its geomorphology is the bedrock ridge and gully slope.

The cross-section LK178 + 713∼733, LK182 + 085∼115, K0 + 320∼350, and RK182 + 017∼047, are selected as the classic surrounding rock sections. According to the geological survey reports, the correlated geological parameters are determined. The geological conditions of the classical cross-section are plotted in Figure 2.

(a)

(b)

(c)

(d)
3. Theory and Methodology
3.1. The Basic Principle
It is assumed that is a fuzzy variable in the domain U, for any , a determined number can be found in the closed interval. The relation between and is called as the absolute membership relationship, its mapping is [30]:
It is called as the absolute membership degree function of .
In the domain , is a variable of , , two opposite fuzzy parameters in the variable are, respectively, , , for the variable , two number , can be determined at any position in two continuous closed intervals, they are called as the relative membership degree relation of to and , the mapping is as follows:
It is called as the relative membership degree of and .
The schematic diagram of dynamic variable about any number in the closed interval can be plotted in Figure 3.

Two fuzzy variables and are, respectively, within two closed intervals, their relative membership degree can be defined as and , and , , .where, is defined as the opposite fuzzy sets. It is plotted in Figure 4.

Similarly, for any variable , there exists, respectively, attractive and repelled sets and , thenwhere when , ; when , ; when , . is called as the relative difference function of to . The mapping is as follows:
It is called as the relative difference function of to . It is plotted in Figure 5.

3.2. The Determination of the Relative Membership Degree
To assess the quality of surrounding rocks, the sample sets can be expressed as follows:where denotes the eigenvalue of the index of sample , . means level of index; the attractive domain is expressed as follows:
When the sets is enlarged based on the upper and lower bound of its adjacent intervals, it is depicted as follows:
According to the physical meaning of assessment index, the level standard of index can be determined by the matrix :
Based on and , the parameter is shown as follows:where when , ; when , then ; when , then .
It is assumed that is the attractive domain of fuzzy variable sets , namely, , is included in the upper and lower domain intervals of , it is plotted in Figure 6 as follows.

According to relevant physical analysis, their relative membership degree can be expressed as follows:
3.3. The Determination of Index Weights
(1)Assuming that there are cases of surrounding rock and assessment indices, so the original matrix can be expressed as follows:(2)Normalization treatment is conducted for the main indices Xij. The positive indicator is as follows: The negative indicator: where is assessment scheme; is assessment index; and is the corresponding magnitude of the assessment index in the scheme.(3)The determination of proportion about the evaluation index in the scheme, it is expressed as follows:(4)The entropy of the evaluation index is shown as follows:(5)The weight of the evaluation index is depicted as follows:3.4. The Determination of the Assessment Level
Based on (11), (12) and (18), according to the relevant references [31], a synthetic membership degree can be obtained: is normalized, then the magnitude of normalized synthetic membership degree can be expressed as follows:where
Finally, the assessment level can be determined based on the value
3.5. The Assessment Steps
(1)Based on the relevant data and assessment criterion, the eigenvalue matrix and classification matrix can be determined(2)Determining the attractive domain , range matrix and point value matrix (3)According to (11) and (12), the relative membership degree can be obtained(4)Determining the weights of different indices of the surrounding rock quality based on the entropy-weight method(5)Computing the level eigenvalues according to (19)–(22), then the assessment level of surrounding rock quality can be determined on the basis of the magnitudes of , if , then the result is level ( is positive integer)4. The Establishment of Assessment Model about the Surrounding Rocks’ Quality
4.1. The Establishment of Assessment Index Systems
The uniaxial compressive strength of rocks , the integrity coefficients of rock mass , the conditions of the structural plane , the seepage amount of groundwater , and the intersection angle between the axis line and structural plane are, respectively, selected as assessment index. These indices are all qualitative indices. The five risk assessment indices are divided into five levels: very good (I), good (II), common (III), bad (IV), and very bad (V) as shown in Table 1. Based on reference [32], the monitoring value of the assessment index of ten cross-sections can be shown in Table 2.
4.2. The Construction of the Assessment Frame
The assessment of surrounding rock quality not only influences the construction progress but also affects the life safety of construction staff. Hence, it is significant to assess the quality level of surrounding rocks.
A new quality assessment method of surrounding rock quality is suggested based on the variable fuzzy sets theory as presented in Figure 7. At first, to evaluate the quality level of surrounding rock quality, a complete assessment index system is established. Second, the weight of each assessment index are determined according to an entropy weight theory. Third, the relative membership degree are determined using the variable fuzzy sets theory. Then, the magnitudes of synthetic certainty degree are determined; finally, the quality level of surrounding rock is determined.

4.3. The Determination of Quality Level about the Surrounding Rock
4.3.1. Determining the Attractive Domain, Range Matrix, and Point Value Matrix
Based on Table 1, according to (7), the attractive domain can be expressed as follows:
Based on (8), the matrix can be expressed as follows:
Based on (10) and (11), the point value matrix can be depicted as follows:
4.3.2. The Determination of the Relative Membership Degree Matrix
In the first procedure, according to the monitoring value in Table 2, and in combination with the (11) and (12), it should be determined whether the assessment index in Table 2 is located on the left or right of the point , the data of cross-section 1 are selected as an example, when , then , and the point value can be expressed as follows:
When , , , , , and , is located in the out of intervals, so ; when , , ,, , is located in the out of intervals, so ; when , , , , , is located in the right of , it belongs to , based on (12), the relative membership degree can be obtained .
Similarly, the relative membership degree matrix of cross-section 1 can be obtained as follows:
4.3.3. The Determination of Weight Coefficients of Different Indices
According to Table 2 and (16), the specific gravity matrix can be shown in Table 3.
According to Table 3 and (17), the entropy matrix can be shown in Table 4.
Based on (18), the weight coefficients can be calculated in Table 5.
4.3.4. Determining the Comprehensive Relative Membership Degree Vector
According to (19), and in combination with the matrix of , the comprehensive relative membership degree matrix can be obtained in Table 6.
According to (20) and (21), the normalized comprehensive relative membership degree vector can be obtained in Table 7.
4.3.5. Determining the Level of the Surrounding Rock Quality
According to (22), and in combination with Table 7, the ranking feature value of cross-section 1 can be shown in Table 8.
Similarly, the feature value of cross-section 2–10 can be shown in Table 9, respectively.
To test the assessment results of surrounding rock quality, the results in the paper are compared with the other methods, they are shown in Table 10.
The variable fuzzy set assessment model is applied to assess the surrounding rock quality; the complete results are shown in Tables 9 and 10. It can be found in Table 10 that the quality levels of surrounding rock from cross-sections 1 to 10 are different. The surrounding rock quality level at cross-section 1, 2, and 4 are IV; one at the cross-section 5 is II; and one at rest is III. It means that the quality level of surrounding rocks at the cross-sections 1, 2, and 4 is bad. The necessary consolidation measurement should be performed about the surrounding rocks at these cross-sections; one of surrounding rock quality at the cross-section 5 is good, so no measurement needs to be done. Ones at the rest cross-sections are common, so the qualified rate of surrounding rock quality in all cross-sections arrives at 70%.
Based on the comparative results of the assessment model in Table 10, it can be found the results assessed by using the variable fuzzy sets method are consistent with the current specification for ten different cross-sections, except the cross-sections 2 and 9. Its accurate rate arrives at 80% in the text method, which is higher than the results (70%) from the rough set theory [32]. Moreover, in comparison with the rough set methods, the proposed model can accurately convey the risk degree of surrounding rocks. So, the conclusions are drawn that it is feasible to estimate the quality level of surrounding rock by using the entropy weight-variable fuzzy sets model. The method provides accurate results and brings more details about the surrounding rock quality levels. For example, the uniaxial compressive strength of cross-section 6 is 37, which should belong to level IV, according to Table 1. In addition, the degree of membership of the other indices obtained by the variable fuzzy sets model belongs to level III, so the quality level probability of cross-section six at the level III is more extensive than that of grades I, IV, II, and V. Hence, the surrounding rock quality of cross-section six only belongs to level III and almost impossibly belongs to groups I, IV, II, and V. Furthermore, the level of cross-section 6 is more likely to be level III than that of cross-section 7, 8, and 10 because the mean ranking feature value (3.2718) of cross-section 6 for level III is higher than that of cross-section 7 (2.9787), 8 (2.8092), and 10 (2.5486). In total, the results obtained by using the entropy weight-variable fuzzy sets model demonstrate the quality level not only accurately but also further determine the ranking of surrounding rock quality for different cross-sections at the same level.
5. Discussion and Comparative Analysis
5.1. Comparison with Existing Studies
(1)The rough set theory is suggested to evaluate the quality level of surrounding rocks, and good results are obtained. However, relative to the model of rough sets, the proposed model can accurately convey the risk degree of surrounding rocks by adopting the eigenvalue of level H [33]. So, the proposed method is much stricter in the superior grade, and the integrity is improved to assess the quality level of surrounding rocks.(2)In comparison with the other traditional models, the fuzziness and randomness of evaluating index are considered, and interval-oriented evaluation criteria are adopted. So, the proposed method improves the reliability of the assessment process and effectively detects surrounding rock quality status.5.2. The Advantages and Limitations of the Proposed Model
By comparing the appropriate methods, the advantages of the suggested method can be summarized as follows:(1)The proposed method can accurately convey the risk degree of surrounding rocks, so it has higher accuracy.(2)Compared the traditional method, its assessment process has higher reliability and efficiency.
However, the suggested model still exist some limitations. For example, calculation is complicated, and multiple variable parameters required to calculate the degree of difference; these limited its application, but the theory has still great space for the improvement in the future.
6. Conclusions
Considering the uniaxial compressive strength of rocks , the integrity coefficients of rock mass , the conditions of the structural plane , the seepage amount of groundwater , as well as the intersection angle between the axis line and structural plane, a new assessment method is introduced in this paper to assess the surrounding rock quality level based on the fuzzy variable sets theory. The relative membership degree matrix of the assessment sample is determined at first. Then, the weighting coefficients are calculated using the entropy weighting method. Finally, the quality level of surrounding rock is determined using the mean ranking feature value.
The proposed method is applied to assess the surrounding rock quality level. Finally, its result is compared with that of the current specifications and the rough set theory; the results obtained based on the variable fuzzy sets method are consistent with the current specification; its accuracy is 80%. The qualified rate of surrounding rock quality in all cross-sections is 70%. In other words, except the cross-sections 1, 2, and 4, for any different cross-sections, no measurement need be adopted to consolidate the surrounding rock quality. Moreover, the results obtained by using the entropy weight-variable fuzzy sets model demonstrate the quality level not only accurately but also further determine the quality ranking of surrounding rock quality for different cross-sections at the same level. The proposed method can accurately convey the risk degree of surrounding rocks; and in comparison with the traditional method, its assessment process has higher reliability and efficiency. But its calculation is complicated, and multiple variable parameters required to calculate the degree of difference, so the proposed method can still be improved in the future.
In short, the proposed model could provide an alternate way to assess the risk level in civil engineering accurately.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the cross project in Nanyang Institute of Technology (520009); The Start-Up Foundation (510126). Sichuan Science and Technology Program (No: 2020YJ0424), Regional Innovation Cooperation Programs of Sichuan Province (2021YFQ0050), National Natural Foundation in China (no. 41672357), Guizhou Provincial Science and Technology Projects (No. [2020]2004), the Scientific and Technological Projects of Henan Province (No: 222102320296), and the Opening Project of Sichuan Province University Key Laboratory of Bridge Non-Destruction Detecting and Engineering Computing (2022QYJ02).