Research Article  Open Access
Xueying Bao, Haiwen Li, "Study on the Evaluation Method of Subgrade Slope Green Protection Effect in DryHot Valley of SichuanTibet Railway", Mathematical Problems in Engineering, vol. 2020, Article ID 7159582, 16 pages, 2020. https://doi.org/10.1155/2020/7159582
Study on the Evaluation Method of Subgrade Slope Green Protection Effect in DryHot Valley of SichuanTibet Railway
Abstract
In view of the harsh environment in the dryhot valley of SichuanTibet railway and the important role of the green protection of the subgrade slope in railway construction and ecological restoration, this paper evaluates the effect of the green protection of the subgrade slope in dryhot valleys of the SichuanTibet railway. First, an evaluation index system for the green protection effect of railway subgrade slopes is established in terms of the soil matrix quality, vegetation community quality, and protection performance of the slopes. Second, game theory is adopted to combine the improvedgroupanalytichierarchyprocess method and the vectoranglecosine method and thus determine the weight of each evaluation index. Moreover, the membershipcloudgravitycenter method is used to evaluate the green protection level of the railway subgrade slope, and an evaluation cloud map is drawn with the help of MATLAB software to further analyze the evaluation. Finally, a section subgrade slope of the DHV area of the SichuanTibet railway is selected as the evaluation object, and the green protection effect is evaluated to verify the applicability and effectiveness of the model. The study provides a theoretical basis for the protection of subgrade slopes and ecological restoration in the DHV area.
1. Introduction
With the rapid development of social economy, the large number of exposed slopes caused by the increasing railway construction not only destroyed the vegetation along the railway but also caused a series of ecological problems and geological disasters [1], such as the sudden loss of slope plant species, slope structure damage, landslides, and debris flow [2]. The subgrade slope, as the foundation to maintain the stability of the railway subgrade, is exposed to the air for a long time, so it is strongly influenced by the natural factors, resulting in the strength attenuation; meanwhile, the rock and soil on the slope surface are ruptured, which can easily aggravate the collapse of the slope and cause the subgrade instability [3]. Geological disasters, such as subgrade collapse, landslides, and spalling of the slope surface, may seriously damage the railway and surrounding environment if the subgrade slope is not properly protected. The green protection (plant protection) of a slope not only prevents these disasters but also plays a role in consolidating the slope, growing vegetation along the subgrade slope, and protecting the ecological environment. The green protection of subgrade slopes has therefore become an indispensable part of railway construction. The implementation of ecological protection measures that take into account the engineering effect and landscape function on the subgrade slope, restore the vegetation on the slope surface, prevent the subgrade from destabilizing, and maintain the health of the slope ecosystem has become a hot issue shared by the majority of soil and water conservation and ecological engineering researchers [4].
In recent years, many experts have analyzed the causes of slope instability [5] and destruction [6] through extensive research and found that the main causes of slope damage are mostly rainfall [7], floods [8], and wind action [9] that eroded the soil of the slope, cracked the subgrade slope, reduced the resistance of the subgrade slope to external forces, and eventually caused the slope to collapse and lose stability. Against the various causes of subgrade slope instability and destruction, many slope protection technologies and methods have been proposed after a large number of engineering tests, which are basically divided into engineering protection technology and plant protection technology. The forms of railway subgrade slope protection [10] include plant protection of hollow brickforeign soil, water interception framework protection, holewindow protection wall, assembled framework protection [11], anchor frame beam protection, prestressing force anchor cable antislide pile supporting technology [12], threedimensional earthwork mesh cushion protection [13], ecological bags, revegetation bags, original ecological tridimensional slope vegetation technology [14], plant fiber blanket road slope protection technology [15], etc. The engineering protection technology is gradually replaced by vegetation slope protection technology because of its high cost, and the protection effect is not as good as vegetation slope protection technology. Vegetation slope protection technology not only saves cost but also green environment and achieves the purpose of slope protection, which is the future development trend.
However, few studies have used mathematical models to evaluate the effect of subgrade slope protection. Moreover, there is no research on the protection of subgrade slopes and related aspects in special areas and special environments, such as the dryhot valley (DHV) areas. Against the above background, the present paper establishes an evaluation index system for the green protection effect of railway subgrade slopes in a DHV area. Game theory (GT) is then adopted to combine the improvedgroupanalytichierarchyprocess (IGAHP) method and the vectoranglecosine method and thus determine the evaluation index weights. Moreover, a membershipcloudgravitycenter (MCGC) mathematical model is proposed to evaluate the green protection effect of a subgrade slope protection project in a DHV area of the SichuanTibet railway. Since the environment of the DHV area of SichuanTibet railway is very fragile, the environmental protection [16] and ecological restoration of the area are the key issues in the construction of the railway. The study not only improves the lack of evaluation of the subgrade slope protection effect but also provides a theoretical basis for the subgrade slope protection and subsequent monitoring and management of the SichuanTibet railway under construction. Moreover, it provides a quantitative evaluation method for the evaluation of subgrade slope protection and also has a certain reference value for slope protection, environmental protection, and ecological restoration in DHV areas.
2. Formation and Distribution of DHV Areas
2.1. Reasons for the Formation of DHV Areas
The DHV areas [17] are formed by the superposition of the original DHV areas and the secondary DHV areas. The original DHV areas already existed, and a special climate was formed due to the combination of special landforms and local microclimates. The area has a tropical and subtropical monsoon climate with a distinct dry and wet season. Moreover, the area has a unique gully system [18], where the foehn effect [19] and local circulation effect of valley wind often occur, causing warming and water loss. Forest vegetation in DHV areas is very difficult to recover. The main reason is that the lack of water makes large areas become barren, especially the loss of surface soil on the slopes of river valleys, which exposes large areas of bare soil and makes the afforestation survival rate low. Secondary DHV areas are formed by the continuous negative disturbance of deforestation. Some studies believe that most of the original vegetation along the banks of the Jinsha River between 3,000 and 5,000 years ago was green broadleaved oak forests. The above factors have combined to create the extremely fragile ecological environment of the DHV areas. The geomorphological characteristics of DHV areas are shown in Figure 1.
2.2. Distribution and Characteristics of DHV Areas
The DHV areas in China are mainly distributed [17] on the slopes of the valleys of the Minjiang River, Dadu River, Yalong River, Jinsha River, Nujiang River, Lancang River, and other rivers in the southwest of China. The areas are mainly located in the arid and semiarid rivervalley areas with latitude of 23° 00′∼28° 10′ north, longitude of 98° 50′∼103° 50′ east, and altitude below 1500 m [20].
The DHV areas [21] are located in a hightemperature valley area, which is rich in solar and thermal energy resources. It has a hot climate with little rain, severe soil erosion, and fragile ecology. Cold, drought, high wind, fire, and other natural disasters are particularly prominent. The dry season in the DHV areas lasts for more than 8 months each year. The average annual rainfall is about 800 mm (rainfall is 550∼650 mm, and the dry season rainfall is 50∼160 mm). However, the evaporation is more than 3 times the rainfall, the relative humidity is small, and the drought index is about 1.5 to 2.5. The highest surface temperature of bare land reaches 75°C, and the water content of topsoil within 40 cm is close to zero, in April and May. Most of the sites in the DHV areas have thick soil, but the vegetation cover is less than 5%. The vegetation in the DHV areas is seminatural savanna, which has simple structure, low species richness, obvious dominance of herb layer, suboptimal shrub layer, and few trees [17].
3. Evaluation Index System
3.1. Analysis of Evaluation Factors
Slope protection engineering involves multidisciplinary and comprehensive system engineering. Therefore, the types and indexes of evaluations of engineering effects must have different properties and characteristics, and problems should be explained at different levels. Meanwhile, such indexes should be interconnected to reflect multiple aspects of the purpose of protection engineering effects comprehensively. From the perspective of system engineering, subgrade slope protection engineering can be regarded as a slopematrixvegetation system [22, 23]. This system and the environment form a unified entity having a certain structure and performance through interaction, intertwining, and infiltration. The construction of a scientifically feasible evaluation index system of the engineering effect must therefore be based on the structural integrity and stability and thus the durability of the engineering system. There are two basic units, namely, the matrix and vegetation, of the structure of the slope protection engineering system. The matrix refers to the basic material (i.e., the soil in the present case) provided by the subgrade slope to meet the normal growth of vegetation, which not only provides a reasonable physical structure for vegetation growth but also serves as a platform and carrier [24] for the transformation and utilization of vegetation moisture and nutrients. Meanwhile, the vegetation refers to the plant community on the slope, which is the main biological factor for systematic development and change. The performance of a slope protection engineering system mainly refers to the protection performance of the subgrade slope, which reflects the ability of the slope to resist external environmental interference. According to the structure and performance of the slopematrixvegetation system, the soil matrix quality, vegetation community quality, and protection performance are taken as firstlevel indexes, and secondlevel evaluation indexes are further analyzed and screened.
The soil in the DHV area is dominated by torrid red soil [25]. This type of soil has a low moisture content, high nutrient loss, and a barren soil layer. Meanwhile, the soil guarantees vegetation survival on the slope in the green protection of the slope. The quality of the soil matrix is therefore a factor that cannot be ignored. The soil matrix quality can be measured by the volumetric weight of soil [26], total soil porosity, soil infiltration rate, soil organic matter and mineral nutrient content, soil expansion and contraction, duration of the soil permanent wilting point (PWP), and the ability of the soil to preserve moisture and fertility [27]. Specifically, the volume weight of soil reflects the comprehensive state of soil organic matter, texture, and compactness; therefore, the volume weight of soil is an important physical characteristic that reflects the quality of the soil matrix. Total soil porosity [25] includes capillary porosity and noncapillary porosity. Capillary porosity refers to the main pore size for soil water retention, and it provides water for plant growth. Noncapillary porosity refers to the main breathable pore size of the soil, and it provides oxygen and nutrients to plant roots in the soil. Moreover, the porosity of the soil affects the infiltration rate [28] of the soil. The infiltration rate affects the circulation, maintenance, and storage of water in the soil. The water infiltrated into the soil is the main source of water for the plant root system. Torrid red soil is compact and the infiltration rate decreases with increasing depth of the soil layer. The infiltration rate is therefore an important factor in evaluating the quality of the soil matrix. Carbonaceous compounds in soil are the main components of soil organic matter [29] and are one of the main sources of nutrients needed for the growth of slope protection vegetation. The other source is mainly mineral nutrients [30] containing trace elements such as total N, total P, available P, and available K [31] in the soil. The content of organic matter and mineral nutrients in the soil is therefore an indispensable factor in measuring the quality of the soil matrix. There are distinct dry and wet seasons in the DHV area, resulting in the large expansion and contraction of soil [32] due to temperature changes. The soil expansion and contraction creates cracks in the slope. In the long run, such cracks cause the slope to collapse, and expansion and contraction damages plant root systems and causes plants to die. Therefore, the root systems of slope plants must have a certain tensile strength. Moreover, the alternation of a dry and wet environment affects the mineralization rate of organic carbon in soil and thus the nutrient content of the soil. The soil PWP [33] in the DHV area lasts for a long time, usually 78 months, and this is the main cause for the sparse vegetation. Slope protection must therefore improve the quality of the soil matrix and thus shorten the duration of the soil PWP to increase the survival rate of plants. Owing to severe soil coarsening in the DHV area, the soil's ability to retain water and fertilizer is reduced, and soil fertility is reduced. This is a major cause of soil infertility, and the soil coarsening resistance and water and fertilizer retention should thus be improved in slope protection.
The quality of a vegetation community is mainly measured by the vegetation coverage, the vegetation survival rate, incidence of vegetation diseases and insect pests, vegetation species richness index, rationality of vegetation spatial collocation, ability of vegetation to resist drought and heat, and tolerance of vegetation to barrenness. Owing to the hot climate in the DHV area, the soil is relatively infertile and lacks nutrients, resulting in extremely scarce vegetation [34] in this area. However, the construction of the railway greatly affected the environment of the area, and plants are therefore planted in the subgrade slope protection project not only to stabilize the subgrade but also to improve the local ecological environment. The numbers and types of slope protection vegetation must meet the requirements of slope protection and environmental improvement in the area, and vegetation coverage [35] can be used as a measure of the number of plants and the area covered by vegetation. Owing to the harsh environment, there are many plant diseases and insect pests in the DHV area, which is a major reason for the low rate of plant survival. The survival rate of plants must therefore be improved by improving the soil and selecting suitable plant species. For slope protection vegetation, shrubs and herbaceous plants that are rich in variety and suitable for DHV growth should be selected. This vegetation has developed roots and strong root meristem, which can further improve the protection ability of slope by improving the soil shear strength. Moreover, vegetation with a strong meristematic ability can improve the coverage of slope vegetation in the short term and enhance the ability of slope stabilization, and vegetation suitable for the area can improve the soil [36] and enhance soil fertility [37]. When planting vegetation, it is also necessary to consider the rationality of the spatial arrangement and distribution. A reasonable structural arrangement can further improve the survival rate of plants and play a positive role in reinforcing the slope. Moreover, the climate of the DHV area is arid and hot; therefore, the selected vegetation should flourish in such a climate.
The protective performance of a slope mainly refers to the slope’s resistance to leaching erosion, antiscouring ability, ability to protect against wind and sand, ability to prevent water and soil loss, and ability to prevent natural disasters. Rainfall is concentrated in the DHV area, but evaporation is generally 3–6 times the rainfall [25]. The amount of water infiltrating the soil is too large in the season of concentrated rainfall, and the infiltration water dissolves and leaches nutrients in the soil, such that the slope resistance to leaching erosion is important. Moreover, a large amount of rainwater will scour the slope surface and base, resulting in water and soil loss. Therefore, the slope must also have a strong ability to prevent water and soil loss and scouring. The land is gradually desertified in the DHV area owing to the coarsening of the soil. In addition, there is often windy weather, which causes severe weather such as sandstorms. Therefore, the ability of a slope to resist wind and sand [38] is also important. Moreover, earthquakes, mudslides, landslides, and floods often occur in the area. Subgrade slope protection projects must therefore be able to resist sudden natural disasters and prevent slope collapse from causing serious damage to the environment and railways.
From the above analysis, the more important factors in evaluating the green protection effect of a slope in the DHV area are summarized in Table 1.

3.2. Construction of an Evaluation Index System
Through the above analysis of the soil, vegetation, hydrology, climate, topography, landform, and other geographical features of the DHV area in southwest China and referring to Interim Provisions on Green Protection Techniques for Railway Subgrade Slopes (2003) and relevant research results [39], the three firstlevel indexes of the soil matrix quality, vegetation community quality, and protective performance of the slope are determined, and an evaluation index system is established through the further refinement and screening of 20 evaluation indexes as secondlevel indexes, as shown in Figure 2.
A fivelevel index standard, as described in Table 2, is established on the basis of the evaluation index system, the grading data of organic matter content of a soil survey conducted in China, relevant research results [40–42], and consultation with experts.

4. Combination Weighting Model
Methods of determining weights can be divided into subjective weighting methods and objective weighting methods. The disadvantage of a subjective weighting method is the excessive reliance on the opinions of experts while the disadvantage of an objective weighting method is the excessive reliance on quantitative data, but these two methods are complementary. Therefore, to ensure that the weight distribution reflects both the subjective information and objective information, this paper uses GT to combine the weights obtained from both the subjective weighting method (i.e., the IGAHP method) and the objective weighting method (i.e., the vectoranglecosine method), which make the weighting more reasonable and scientific.
4.1. IGAHP Method
The IGAHP method [43] is based on the traditional analytic hierarchy process (AHP) [44] method and considers the difference between different expert evaluations and assigns a weighting coefficient to each expert using the Euclidean distance to minimize the expert’s evaluation bias so that the weight of the index is more reasonable.
The weights of m indexes calculated by the tth expert using the AHP method are denoted aswhere is the AHP weight of the ith index given by the tth expert.
The difference between the evaluations of the th and hth experts is then expressed by the Euclidean distance :where and are, respectively, the weights of the ith index given by the th and hth experts, satisfying d_{tt} = 0 and d_{gh} = d_{hg} ≥ 0.
Equation (2) shows that closer values of and (i.e., smaller differences between the evaluations of the two experts) result in smaller d_{gh}; if d_{gh} = 0 and ≠ h, then the evaluations of the two experts are exactly the same.
The sum d_{t} of the values of the degree of similarity of the evaluations made by the tth expert and other experts is expressed bywhere T is the number of experts while d_{tj} is the degree of similarity of the evaluations of the tth and jth experts.
The evaluation weight coefficient λ_{t} of the tth expert is calculated as
In summary, the subjective weight W_{1} determined using the IGAHP method can be calculated as
4.2. VectorAngleCosine Method [45]
Step 1. Construct the optimal value vector and worst value vector of the evaluation index aswhere and arewhere is the optimal value vector of the evaluation index, is the optimal value of the ith evaluation index, is the worst evaluation vector of the evaluation index, is the worst evaluation value of the ith evaluation index, m is the number of evaluation indexes, n is the number of experts, I_{1} is a positive index set, I_{2} is a negative index set, and V_{ij} is the basic data of the ith index measured by the jth expert.
Step 2. Calculate the matrix R of the relative deviation between the evaluation object and the optimal value vector and calculate the matrix of the relative deviation between the evaluation object and the worstvalue vector:where r_{ij} is the optimal value of the relative deviation of the ith index calculated from the basic data given by the jth expert while δ_{ij}is the worst value of the relative deviation of the ith index calculated from the basic data given by the jth expert.
Step 3. Determine the value of the vector angle cosine c_{i} of the evaluation index (equation (12)) and then normalize it to obtain the objective weight W_{2} (equation (13)):where the equation is used in normalization and is the vectoranglecosine weight of the ith index.
4.3. GT Combination Weighting Model
The GT [46, 47] combination weighting model is constructed by finding the minimum deviation between subjective and objective weights to obtain a compromise value considering the conflict between subjective and objective weights. This achieves the effect of making an interactive decision between subjective and objective weights, so as to obtain the optimal combination of weights.
Supposing θ_{1} and θ_{2} are the weighting coefficients of subjective and objective weights and L is the number of methods used to determine the weights, the combination weight is calculated aswhere and are the transposed matrices of W_{1} and W_{2}.
Specifically, the combination weights are obtained by the following calculation.
The values of θ_{1} and θ_{2} are optimized according to the aggregation model of GT:
According to the differential characteristics of the matrix, the optimal first derivative is obtained from equation (15):
The weighted coefficients θ_{1} and θ_{2} are obtained by solving equation (16) and then normalizing the solution:where and are the normalized values of θ_{1} and θ_{2}.
The optimized combination weight is then calculated as
5. Evaluation Model
5.1. MCGC Model
The MCGC model [48] is a comprehensive evaluation method suitable for solving uncertainty and fuzziness problems. Through the conversion of qualitative concepts and quantitative data, the evaluation and measurement of the uncertainty problem are realized.
Suppose that U is a quantitative domain composed of precise values and C is a qualitative concept on the domain U. For any element x, there is a random number with stable tendency, which is called the membership degree of x to the qualitative concept C. The distribution of membership degree μ on the domain U forms the membership cloud , and constitutes a cloud drop [49]. The expectation Ex, entropy En, and hyperentropy He are introduced to represent the eigenvalues of the cloud. The expectation Ex is the point value on the domain U that can best represent the fuzzy qualitative concept C, indicating the central value of the fuzzy concept C in the domain, which is the central distribution of the cloud. Entropy En is the degree of uncertainty of qualitative concept C, which reflects the random probability of cloud droplet of qualitative concept C and the acceptable range of cloud droplet of qualitative concept C in domain U, and the greater the entropy value, the greater the randomness. Hyperentropy He is a measure of the uncertainty of entropy, which reflects the degree of cloud dispersion and cloud thickness, as shown in Figure 3.
In Figure 3, Ex is the average of the cloud, which represents the qualitative concept. En represents the value range of qualitative concepts in the domain space, which reflects the “wide” and “narrow” of the cloud. He represents the randomness of the sample values and provides a method to combine the ambiguity and randomness, which reflects the “thickness” and “thinness” of the cloud [50].
5.2. Determining the Eigenvalues of the Standard Evaluation Level
Experts in relevant fields, such as railway subgrade slope construction, slope green protection, and the natural environment of the DHV area, were invited to give numerical intervals for the five evaluation levels of excellent, good, general, poor, and very poor. The inverse cloud generator of MATLAB software was used to generate eigenvalues of the standard evaluation level. The relevant formulas arewhere is the numerical interval of the jth standard evaluation level given by the ith expert; are the values of expectation, entropy, and hyperentropy obtained from the numerical intervals of the standard evaluation levels given by each expert; and k is a constant that can be set according to needs.
Assuming that the weights of the evaluation experts are the same, the eigenvalues of the standard evaluation level arewhere are the eigenvalues obtained by integrating the expectation, entropy, and hyperentropy calculated from the numerical intervals of the jth standard evaluation level given by all experts.
5.3. Aggregation of Evaluation Indexes
Index aggregation refers to determining the corresponding standard evaluation level according to the data of secondlevel indexes collected by each expert to obtain the eigenvalues of secondlevel indexes and then performing parallel aggregation to obtain the weighted eigenvalues of secondlevel indexes. Finally, vertical upward aggregation is performed to obtain the eigenvalues of firstlevel indexes and the eigenvalues of the evaluation result.
5.3.1. Parallel Aggregation
where denotes the weight of the ith expert (i.e., the expert weight coefficient determined using the IGAHP method) while Ex, En, and He, respectively, refer to the expectation, entropy, and hyperentropy of the secondlevel indexes after parallel aggregation.
5.3.2. Vertical Upward Aggregation
where is the weight of the GT combination of the ith index while Ex, En, and He, respectively, refer to the expectation, entropy, and hyperentropy after vertical upward aggregation.
5.4. Generating a Cloud Map of Evaluation Results
The eigenvalues of the evaluation result are finally obtained through parallel aggregation and vertical upward polymerization. The eigenvalues of the evaluation result and the eigenvalues of the five evaluation levels are input to the forward cloud generator of MATLAB software to obtain the cloud map of evaluation results.
6. Case Study
6.1. Generating Eigenvalues of the Standard Evaluation Level
The standard evaluation levels of excellent, good, general, poor, and very poor do not have precise numerical intervals; i.e., there are no eigenvalues of the standard evaluation level. Therefore, two experts on the construction of railway subgrade slopes, two experts on the green protection of subgrade slopes, and two experts on the natural environment in the DHV area were invited to give numerical intervals for the five standard evaluation levels. The intervals are given in Table 3.

The numerical intervals in Table 3 and equations (19) and (20) are input into the reverse cloud generator of MATLAB software to obtain the eigenvalue of the standard evaluation level, as shown in Table 4.

6.2. Engineering Case
The SichuanTibet railway runs from Chengdu, the capital of Sichuan, to Lhasa, the capital of Xizang, which is divided into three sections: Chengdu to Ya’an (Chengya section), Ya’an to Nyingchi (Yalin section), and Nyingchi to Lhasa (Lalin section). The Chengya section was officially opened to traffic on December 28, 2018, the Lalin section has now entered the tracklaying phase, and the most difficult Yalin section is now under construction. Most areas the railway passes through are mountains, hills, and valleys, where the climate is extremely harsh and natural disasters such as floods, sandstorms, and earthquakes are common. In particular, Luding and surrounding areas of the Dadu River Basin, Batang and surrounding areas of the Jinsha River Basin, and parts of the Yapan River and other areas along the SichuanTibet railway have DHV areas.
After consulting the six experts as described above, a section of the railway from Ya’an to Kangding of the SichuanTibet railway, which is mainly located in Dadu River Basin in Luding County and passes through DHV area, is selected as the research object. The area has a subhumid climate on the QinghaiTibet Plateau and is influenced by the southeast and southwest monsoons and the cold air on the QinghaiTibet Plateau. The annual precipitation is 664.4 mm, and the evaporation is 1480.9 mm. November to April of the next year is the dry season, which is a typical DHV area [51]. The soil types mainly include red soil, mountain brown soil, and mostly leached soil, which is very likely to cause geological disasters. Vegetation is mostly alpine shrubs and meadows, as well as Opuntia sp., succulent shrubs, and Acacia farnesiana and Pistacia weinmannifolia communities. The six experts were invited to select different locations of the subgrade slope project to measure the index data through onsite investigation and sample testing. Specifically, the subordinate indexes of soil matrix quality, such as volumetric weight of soil and total soil porosity, are measured by the cutting ring method, and the soil infiltration rate and organic matter content are measured by the round block technique and the CODcr method, respectively. The subordinate indexes of vegetation community quality are determined by the plant community sample survey method, and the species richness index refers to the Shannon–Wiener index. The strength of qualitative indexes such as soil expansion and contraction, antiscouring ability, ability of slope to protect against wind, and sand is given by experts using the visual method according to the actual situation of the study area. The specific data are shown in Table 5.

The evaluation levels of each secondlevel index are obtained from the basic data in Tables 2–5, as shown in Table 6.

6.3. Determining the Combined Weight
6.3.1. IGAHP Weights
The six experts used the AHP method to analyze and compare the importance of the firstlevel indexes of the soil matrix quality C_{1}, vegetation community quality C_{2}, and protective performance of the slope C_{3} and then provided their judgment matrix, which passed a consistency test, as shown in Table 7.

According to the judgment matrix in Table 7, the firstlevel index weights determined using the AHP method are obtained, and these weights are then substituted into equations (1)–(4) to obtain the expert’s weight coefficient λ:
The weight coefficient λ of experts and the firstlevel index weights determined using the AHP method are substituted into equation (5) to obtain the IGAHP weight of the firstlevel indexes:
Similarly, the IGAHP weights of secondlevel indexes are obtained, as shown in Table 8. The detailed calculation process is not described here.

6.3.2. VectorAngleCosine Weights
For the qualitative indexes of the evaluation index system, supposing a score of 10 is the highest score, experts provided a score for the actual conditions. The weight of the vectoranglecosine method of the firstlevel indexes is obtained by substituting the scores of the firstlevel indexes given by experts into equations (6)–(13):
Similarly, the data of quantitative indexes in Table 5 and the scores of qualitative secondlevel indexes given by experts are substituted into equations (6)–(13), and the weights of the vectoranglecosine method for the secondlevel indexes are then obtained, as shown in Table 8. The scores of qualitative indexes and the detailed calculation process are not provided here.
6.3.3. Combined Weights
The GT combination coefficients of the firstlevel indexes are obtained by substituting the subjective and objective weights of the firstlevel indexes determined by the IGAHP method and the vectoranglecosine method into equations (14)–(17):
According to equation (18) and the subjective and objective weights calculated above, the combined weights of the firstlevel indicators are obtained:
Similarly, the GT combination coefficient of the secondlevel indexes is obtained and the combined weights of the secondlevel indexes are then obtained, as shown in Table 8. The detailed calculation process is not described here.
6.4. Index Aggregation
According to the calculation of the subjective weight adopting the IGAHP, the weight of each expert is .
Substituting the eigenvalues of the secondlevel indexes obtained according to Tables 4–6 and the weight of the experts into equation (21) gives the eigenvalues of the secondlevel indexes after parallel aggregation, as shown in Table 9.

According to the values in Tables 8 and 9, vertical upward aggregation is performed using equation (22), and the eigenvalues of the firstlevel indexes and the evaluation result are obtained, as shown in Table 10.

6.5. Cloud Map of Evaluation Results
The eigenvalues of the standard evaluation levels in Table 4 and the eigenvalues of the evaluation result in Table 10 are input into the forward cloud generator of MATLAB software, and a contrast cloud map of evaluation results and standard evaluation levels of the green protection effect of the subgrade slope in the DHV area of the SichuanTibet railway is obtained, as shown in Figure 4.
In Figure 4, the blue map is the cloud map of the standard evaluation level “excellent,” the pink map is the cloud map of the standard evaluation level “good,” the yellow map is the cloud map of the standard evaluation level “general,” the green map is the cloud map of the standard evaluation level “poor,” the red map is the cloud map of the standard evaluation level “very poor,” and the black map is the cloud map of the evaluation result. The figure shows that the cloud map of the evaluation result is basically contained in the cloud map of the standard evaluation level “good,” and the evaluation value corresponding to the vertex of the evaluation result cloud map (i.e., the expectation of the evaluation result) is close to the expectation of the standard evaluation level “good.” Therefore, the green protection effect of this railway subgrade slope in the DHV area is “good.” The evaluation result is consistent with the actual protection effect of the expert onsite investigation.
To further determine the accuracy of the evaluation results, an evaluation cloud map comparing the firstlevel indexes of the soil matrix quality C_{1}, vegetation community quality C_{2}, and protective performance of the slope C_{3} with the standard evaluation level is also generated, as shown in Figures 5–7.
Figures 5–7 show that the evaluation cloud maps of the soil matrix quality C_{1}, vegetation community quality C_{2}, and protective performance of the slope C_{3} (i.e., the firstlevel indexes) are all between the cloud maps of the standard evaluation levels “good” and “excellent” and very close to the cloud map of the standard evaluation level “good.” In particular, the evaluation cloud maps of C_{1} and C_{3} are basically completely contained in the cloud map of the standard evaluation level “good,” such that the evaluation level of the three firstlevel indexes is “good.” By determining the evaluation level of the firstlevel indexes, it is further confirmed that the green protection effect of the railway subgrade slope in this section is “good.”
7. Discussion
The evaluation index system basically contains all the factors that affect the green protection effect for the railway subgrade slope in the DHV area, and the evaluation results are thus convincing. However, the limitation of the evaluation index system is that half the evaluation indexes are qualitative and half are quantitative; i.e., there are too few quantitative indexes and too many qualitative indexes. In addition, the qualitative indexes are only giving their status levels according to the actual conditions observed by experts, which makes the assessment of qualitative indexes too subjective because of the influence of experts themselves. Therefore, how to measure and evaluate qualitative indexes needs further research to improve the accuracy of evaluation results. Moreover, the evaluation index system is established for DHV areas with fragile ecological environments, and some indexes used in the evaluation system may not be applicable to other regions. For instance, duration of soil PWP C_{14}, soil expansion and contraction C_{18}, ability of vegetation to resist drought and heat C_{26}, and slope resistance to leaching erosion C_{31} indexes are unique to DHV areas. The DHV areas are mostly torrid red soil, and soil water content at a depth of 20 cm for 7 to 8 months of the year is in a PWP state, which causes plants to die for lack of water. Although index C_{18} may be applicable to other areas, soil expansion and contraction is particularly prominent in DHV areas. There are distinct dry and wet seasons in the DHV area, which results in the large expansion and contraction of soil due to changes in dry and wet conditions, resulting in cracks in the soil. In severe cases, the plants will die due to root fracture. The DHV area is extremely dry and hot, and the slope protection vegetation must be able to resist drought and hot climate in order to play a better slope protection role; therefore, the index C_{26} is exclusively suitable for the study area. Although the DHV areas are very dry, the rainfall is very concentrated during the wet season, which makes the water infiltrated into the soil too large, dissolving the nutrients from the soil; therefore, index C_{31} is also unique to the areas. The four indexes above are selected based on the unique geological environment and climatic characteristics of the DHV areas, which is relatively applicable only to the DHV areas. The evaluation system therefore needs to be improved for the green protection effect of railway subgrade slopes in other areas according to the actual situation in those areas.
The advantage of the present study is that a subjective weighting method and objective weighting method are combined to determine the weighting, such that the weight distribution is based on not only the effects of subjective human factors but also the effects of objective sample data, which makes the weight more comprehensive and reasonable. In references [41, 44], only a subjective weight method (AHP) is used, which relies too much on a person’s information and fails to consider the importance of experts themselves. However, this paper uses the IGAHP method not only to consider the importance of the experts themselves but also to use the GT to combine the vectoranglecosine method with the IGAHP method so that the weights are more accurate. As for the evaluation model, reference [41] uses the grey system theory to calculate the grey level of the evaluation object by constructing the definite weighted function so as to evaluate the quality of slope protection engineering. Reference [44] uses fuzzy comprehensive evaluation method to calculate the comprehensive membership of the research object and determine the slope protection effect of the slope protection project by comparison with the standard membership. Both methods above solve the problem by considering the uncertainty and ambiguity of the research problem, and they have their own advantages. However, the MCGC model not only considers the fuzziness of the research problem but also evaluates it by calculating the sample data and drawing the contrast cloud map with MATLAB software. Furthermore, the model intuitively reflects the evaluation results through a combination of graphics and text, making the evaluation process more detailed and the evaluation results more accurate.
The evaluation results are in good agreement with the actual situation of the railway subgrade slope protection project in the studied section of the SichuanTibet railway. This shows that the construction of the subgrade slope protection project complies with the requirements of China to protect the ecological environment when constructing railways. However, only six experts were invited to collect data, and small data sample resulted in slight bias in the calculation results. Therefore, in future research, it will be necessary to take as many sample data as possible to improve the accuracy of the results. Moreover, it is worth mentioning that the present research provides a theoretical basis for government to evaluate the quality of the green protection engineering of the railway subgrade slope in the DHV area.
8. Conclusion
(1)Environmental and geographical characteristics were analyzed and combined with the characteristics and requirements of railway subgrade slope protection, and an evaluation index system for the green protection effect of railway subgrade slopes in a DHV area was then established.(2)Subjective and objective weights of the evaluation indexes were calculated using the IGAHP method and the vectoranglecosine method, respectively, and GT was then used to combine the subjective and objective weights such that the combined weight not only highlights the particularity of the DHV area but also eliminates the conflict and difference between subjective and objective factors. In addition, an evaluation model was established using the MCGC method, and an evaluation cloud map was drawn with MATLAB software to analyze the evaluation result.(3)The GT combination weighting model and MCGC evaluation model were used to evaluate the protection effect in a practical engineering case, which is a section of subgrade slope protection engineering passing through the DHV area of the SichuanTibet railway. The evaluation results are in good agreement with the actual engineering situation, demonstrating the applicability and effectiveness of the model.
Data Availability
The data used to support the findings of the study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (grant no. 51768034).
References
 K. Wu, Q. Sheng, and Y. H. Zhang, “Development of realtime remote monitoring and forecasting system for geological disasters at subgrade slopes of mountainous highways and its application,” Rock and Soil Mechanics, vol. 31, no. 11, pp. 3683–3687, 2010. View at: Google Scholar
 D. Wang and J. J. Wang, “Selection of optimal slope protection plants based on fuzzy neural network,” Journal of Northeast Forestry University, vol. 7, pp. 116–119, 2011. View at: Google Scholar
 K.S. Cha and T.H. Kim, “Evaluation of slope stability with topography and slope stability analysis method,” KSCE Journal of Civil Engineering, vol. 15, no. 2, pp. 251–256, 2011. View at: Publisher Site  Google Scholar
 T. Yang, S. C. Li, and H. L. Sun, “Key issues in research of ecoengineering for rock slope protection,” Research of Soil and Water Conservation, vol. 6, pp. 14–16, 2007. View at: Google Scholar
 Y.C. Hwang, “Change of slope stability due to slope inclination and surface conditions,” Journal of the Korean Geotechnical Society, vol. 31, no. 2, pp. 5–11, 2015. View at: Publisher Site  Google Scholar
 H. G. Jeong and Y. S. Seo, “Slope stability analysis and landslide hazard assessment in tunnel portal area,” Journal of Korean Tunnelling and Underground Space Association, vol. 15, no. 4, pp. 387–400, 2013. View at: Publisher Site  Google Scholar
 J. H. Shin, H. M. Kim, and K. J. Jang, “Effect of pore water pressure on slope stability by using coupled finite element analysis,” Journal of the Korean Geotechnical Society, vol. 25, no. 2, pp. 25–35, 2009. View at: Google Scholar
 Y. A. Makarova and A. Y. Manukovskiy, “The effect of floodwater on slope disturbance,” Bulletin of Higher Educational Institutions. Lesnoi Zhurnal (Forestry Journal), vol. 2, no. 2, pp. 70–76, 2018. View at: Publisher Site  Google Scholar
 J. F. Lian, Q. Luo, and T. Xie, “Mechanics effects of frame protection for shallow stability of soil slope and frame structure optimization,” Journal of the China Railway Society, vol. 39, no. 7, pp. 134–140, 2017. View at: Google Scholar
 W. X. Cui, “Research on the optimization of highspeed railway earthwork slope protections,” Journal of Railway Engineering Society, vol. 10, pp. 17–22, 2019. View at: Google Scholar
 S. Cheng, “Application of assembled skeleton slope protection in subgrade slope protection of QinghaiTibet railway,” Railway Engineering, vol. 10, pp. 69–71, 2003. View at: Google Scholar
 A. Wang, “Application of stabilizing piles with prestressed anchor cables retaining structure in subgrade slope protection,” Highway Engineering, vol. 43, no. 1, pp. 145–148, 2018. View at: Google Scholar
 G. Y. Wang, Y. Wang, and N. Xu, “Test and study on erosion characteristics of 3D geonet protection slop,” Research of Soil and Water Conservation, vol. 24, no. 1, pp. 79–83, 2017. View at: Google Scholar
 Z.X. Yan, Z.H. Ren, C.M. Yan, P. Jiang, and H.Y. Wang, “Study on original ecological tridimensional slope vegetation,” Journal of Mountain Science, vol. 10, no. 6, pp. 932–939, 2013. View at: Publisher Site  Google Scholar
 H. Y. Liu, L. Liu, and X. J. Li, “Comprehensive benefit evaluation on the protection technique of plant fiber blanket on the road side slope,” Journal of Soil and Water Conservation, vol. 33, no. 1, pp. 345–352, 2019. View at: Google Scholar
 Y. W. Luo, P. Han, and G. Y. Zheng, “Research on the environmentfriendly railway route selection of Ya’ankangding section of SichuanTibet railway,” Journal of Railway Engineering Society, vol. 8, pp. 1–5, 2017. View at: Google Scholar
 S. Q. Luo, F. D. Xu, and S. F. Li, “Ecological restoration measures for difficult sites in dry and hot mountain valleys in southwest China,” Journal of Green Science and Technology, vol. 14, pp. 11–13, 2019. View at: Google Scholar
 Y. Dong, D. Xiong, Z. Su et al., “The influences of mass failure on the erosion and hydraulic processes of gully headcuts based on an in situ scouring experiment in dryhot valley of China,” Catena, vol. 176, pp. 14–25, 2019. View at: Publisher Site  Google Scholar
 A. L. Khan, G. R. McMeeking, J. P. Schwarz et al., “Nearsurface refractory black carbon observations in the atmosphere and snow in the McMurdo dry valleys, Antarctica, and potential impacts of foehn winds,” Journal of Geophysical Research: Atmospheres, vol. 123, no. 5, pp. 2877–2887, 2018. View at: Publisher Site  Google Scholar
 C. R. Ou, Q. K. Zhu, and Y. Y. Sun, “Research progress in the landcape pattern of dryhot valley in southwest China,” Journal of West China Forestry Science, vol. 44, no. 6, pp. 137–142, 2015. View at: Google Scholar
 L. Rong, X. Duan, G. Zhang, Z. Gu, and D. Feng, “Impacts of tillage practices on ephemeral gully erosion in a dryhot valley region in southwestern China,” Soil and Tillage Research, vol. 187, pp. 72–84, 2019. View at: Publisher Site  Google Scholar
 T. Hubble, S. Clarke, A. Stokes, and C. Phillips, “4th international conference on soil bio and ecoengineering (SBEE2016) 'the use of vegetation to improve slope stability',” Ecological Engineering, vol. 109, pp. 141–144, 2017. View at: Publisher Site  Google Scholar
 G. Tardío and S. B. Mickovski, “Implementation of ecoengineering design into existing slope stability design practices,” Ecological Engineering, vol. 92, pp. 138–147, 2016. View at: Publisher Site  Google Scholar
 J. Geris, D. Tetzlaff, J. McDonnell, and C. Soulsby, “The relative role of soil type and tree cover on water storage and transmission in northern headwater catchments,” Hydrological Processes, vol. 29, no. 7, pp. 1844–1860, 2015. View at: Publisher Site  Google Scholar
 D.H. Xiong, H.Y. Zhou, Z. Yang, and X.B. Zhang, “Slope lithologic property, soil moisture condition and revegetation in dryhot valley of Jinsha river,” Chinese Geographical Science, vol. 15, no. 2, pp. 186–192, 2005. View at: Publisher Site  Google Scholar
 Q. B. Chen, K. Q. Wang, and Y. M. Li, “Effect of different vegetation types on soil amelioration in dryhot valley of Jinshajiang River basin,” Journal of Soil Water Conservation, vol. 19, no. 2, pp. 67–70, 2003. View at: Google Scholar
 A. Ghosh, A. Singh, R. Kumar et al., “Soil enzymes and microbial elemental stoichiometry as bioindicators of soil quality in diverse cropping systems and nutrient management practices of Indian Vertisols,” Applied Soil Ecology, vol. 145, 2020. View at: Publisher Site  Google Scholar
 Z. Q. Liu, N. J. Lang, and K. Q. Wang, “Research on soil infiltration law of different landues types in Jinsha River dryhot valley,” Advanced Materials Research, vol. 726731, pp. 3867–3871, 2013. View at: Publisher Site  Google Scholar
 N. Y. Musadji, L. Lemee, L. Caner, G. Porel, P. Poinot, and C. GeffroyRodier, “Spectral characteristics of soil dissolved organic matter: longterm effects of exogenous organic matter on soil organic matter and spatialtemporal changes,” Chemosphere, vol. 240, pp. 124808–124814, 2020. View at: Publisher Site  Google Scholar
 L. Xiao, P. Li, P. Shi, and Y. Liu, “Soil nutrient stoichiometries and enzymatic activities along an elevational gradient in the dryhot valley region of southwestern China,” Archives of Agronomy and Soil Science, vol. 65, no. 3, pp. 322–333, 2018. View at: Publisher Site  Google Scholar
 X.w. Duan, X. Han, J.m. Hu, D.t. Feng, and L. Rong, “A novel model to assess soil productivity in the dryhot valleys of China,” Journal of Mountain Science, vol. 14, no. 4, pp. 705–715, 2017. View at: Publisher Site  Google Scholar
 N. F. Zhao, W. M. Ye, B. Chen, Y. G. Chen, and Y. J. Cui, “Modeling of the swellingshrinkage behavior of expansive clays during wettingdrying cycles,” Acta Geotechnica, vol. 14, no. 5, pp. 1325–1335, 2019. View at: Publisher Site  Google Scholar
 Z. Xiang, “Degradation of ecosystem and ways of its rehabilitation and reconstruction in dryhot valley,” Resources and Environment in the Yangtze Basin, vol. 9, no. 3, pp. 376–383, 2000. View at: Google Scholar
 J. Pei, W. Yang, and Y. P. Cai, “Relationship between vegetation and environment in an aridhot valley in southwestern China,” Sustainability, vol. 10, no. 12, pp. 1–14, 2018. View at: Publisher Site  Google Scholar
 H. Wu, D.H. Xiong, L. Xiao et al., “Effects of vegetation coverage and seasonal change on soil microbial biomass and community structure in the dryhot valley region,” Journal of Mountain Science, vol. 15, no. 7, pp. 1546–1558, 2018. View at: Publisher Site  Google Scholar
 Y. F. Dong, D. H. Xiong, Z. G. Su et al., “The distribution of and factors influencing the vegetation in a gully in the dryhot Valley of southwest China,” Catena, vol. 116, pp. 60–67, 2014. View at: Publisher Site  Google Scholar
 S. Peng, A. Chen, H. Fang, J. Wu, and G. Liu, “Effects of vegetation restoration types on soil quality in Yuanmou dryhot valley, China,” Soil Science and Plant Nutrition, vol. 59, no. 3, pp. 347–360, 2013. View at: Publisher Site  Google Scholar
 Y. Cui, J. Li, A. Chen et al., “Fractal dimensions of trapped sediment particle size distribution can reveal sediment retention ability of common plants in a dryhot valley,” Catena, vol. 180, pp. 252–262, 2019. View at: Publisher Site  Google Scholar
 X. Hu, Z. Chen, and C. J. Li, “Grade evaluation on quality of slope protection engineering by vegetation,” Bulletin of Soil and Water Conservation, vol. 33, no. 3, pp. 180–185, 2013. View at: Google Scholar
 The First Editing Room of Standards Press of China, Compilation of National Standards for Soil and Water Conservation, China Standards Press, Beijing, China, 2011.
 B. T. Wan, X. Y. Bao, and A. C. Li, “Quality evaluation system and application of ecological slope protection project based on unascertainedset pair coupling,” Bulletin of Soil and Water Conservation, vol. 39, no. 2, pp. 108–114, 2019. View at: Google Scholar
 F. Long, S. C. Li, and H. L. Sun, “Evaluation indicator system for effect of ecoengineering for rock slope protection and its application,” Chinese Journal of Rock Mechanics and Engineering, vol. 28, no. supp.1, pp. 3095–3101, 2009. View at: Google Scholar
 J. An, T. X. Xu, and X. Zeng, “Equipment quality condition assessment under fusion information based on combination weighting,” Control and Decision, vol. 33, no. 9, pp. 1693–1698, 2018. View at: Google Scholar
 S. Chakraborty and S. Mukhopadhyay, “Assessing flood risk using analytical hierarchy process (AHP) and geographical information system (GIS): application in Coochbehar district of West Bengal, India,” Natural Hazards, vol. 99, no. 1, pp. 247–274, 2019. View at: Publisher Site  Google Scholar
 H. Shen, J. C. Xie, and J. X. Li, “Combined hydrological forecasting method based on vector cosine,” System Engineering Theory and Practice, vol. 32, no. 7, pp. 1591–1597, 2012. View at: Google Scholar
 S. Abapour, B. MohammadiIvatloo, and M. T. Hagh, “Robust bidding strategy for demand response aggregators in electricity market based on game theory,” Journal of Cleaner Production, vol. 243, p. 118393, 2020. View at: Publisher Site  Google Scholar
 K. X. Huang, C. J. Zhou, Y. Q. Qin, and W. Tu, “A gametheoretic approach to crosslayer security decisionmaking in industrial cyberphysical systems,” IEEE Transactions on Industrial Electronics, vol. 67, no. 3, pp. 2371–2379, 2020. View at: Publisher Site  Google Scholar
 Z. L. Wei, C. L. Chen, and Y. Q. Dai, “Research on quality evaluation of scheme for equipment system based on IGAHPIEWGTMCGC,” Fire Control & Command Control, vol. 42, no. 11, pp. 114–120, 2017. View at: Google Scholar
 Z. Sun and L. H. Zhao, “Reliability realtime evaluation of JTC based on cloud theory,” Journal of the China Railway Society, vol. 41, no. 2, pp. 89–96, 2019. View at: Google Scholar
 Z. Huang, H. L. Fu, and J. B. Zhang, “Comprehensive evaluation model of shield tunnel construction risk based on cloud theory,” Journal of Railway Science and Engineering, vol. 15, no. 11, pp. 3012–3020, 2018. View at: Google Scholar
 L. P. Liang, Y. G. Liu, and Z. H. Tang, “Geologic hazards susceptibility assessment based on weighted information value a case study in Ludi county, Sichuan province,” Bulletin of Soil and Water Conservation, vol. 39, no. 6, pp. 176–182, 2019. View at: Google Scholar
Copyright
Copyright © 2020 Xueying Bao and Haiwen Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.