Abstract

Risk assessment of debris flow is conducted by multicriteria decisions. Based on the shortcomings of the existing methods in determining the weight of assessment factors, this paper proposes a new approach to conduct a risk assessment of debris flow. This new approach regards the weight of factors as a uniform random variable, whose bounds could be determined by the equal weight method, maximal deviation method, and entropy method. The results of this new approach are obtained by Monte Carlo simulation. According to the risk of 72 debris flows collected in Beichuan, Sichuan, China, this new approach proves convergent. It is suggested that the minimum sample amount of Monte Carlo simulation should be 63095. The result also demonstrates that sorted results with different weights of factors vary a lot, so it is not convincing to sort samples with a specific weight.

1. Introduction

The debris flow is a common geological hazard in rural areas [1], and debris flows are usually caused by intense rainfall [2]. The risk assessment of debris flow could provide pivotal references for policymakers to make prevention measures that achieve unity and economic and safety [3].

There are so many factors that can determine the risk of economic efficiency and safety. Therefore, the risk assessment should be regarded as a multicriteria problem [4]. People have conducted systematic researches on risk assessment of debris flow for the long term. A neural network model for the risk assessment of debris flow was established by Chong and Chao [5], which regarded the length of the creek, the average slope, the effective watershed area, the shape coefficient, the median size of soil grain, the effective cumulative rainfall, and the effective rainfall intensity as significant factors. Zhang et al. [6] established an extension model considering the impact of loose material volume per square kilometer, loose material supply length ratio, average gradient of the main channel, average hill slope, drainage density, curvature of the main channel, and poor vegetation area ratio for the risk assessment of debris flow. A professional debris flows early warning system considered the effective watershed, the length of the effective channel, the slope of the effective channel, the rocks in the effective watershed, the collapsed area within the effective watershed, the effective accumulated precipitation, the effective rainfall intensity, and the impact of the vegetation index proposed by Kung et al. [7]. The genetic algorithm was adopted to predict debris flow due to the excellent capacity of the global search [8]. The risk assessment of debris flow method based on a cloud model was developed by Cao et al. [9], and it was applied to the risk assessment of 29 debris flow gullies in Zhirui (Inner Mongolia, China). The Bayesian network model for debris flow prediction was proposed by Banihabib et al. [10], and this research indicated that average basin elevation, watershed area, current rainfall, and discharge of 1 day ago are the effective predictors in forecasting debris flows. The application of the fuzzy C-means algorithm in the risk assessment of debris flow has also achieved more reliable results [11]. The aforementioned studies provided many available approaches to conduct a risk assessment of debris flow. However, they failed to discuss the weight definition of assessment factors further.

As for multicriteria decision problems, weight calculation is one of the most critical problems. When the weight is clear in a continuous number multicriteria decision problem (assessment factors contain no interval number), decision schemes can quickly be ranked by the matrix multiplication formula. There are standard weight determination methods, including equal weight method, analytic hierarchy process [12], maximum deviation method [13], entropy method [14], and projection pursuit method [15]. There are subjective and objective methods, the equal weight method and analytic hierarchy process are subjective methods, while the maximum deviation method, entropy method, and projection pursuit method are objective methods. Different weight definition methods could lead to weights of comparatively large variance. As a result, it has long been focused on multicriteria decisions on adopting different methods to get the weight reasonably [1618]. Some methods combining subjective and objective approaches (e.g., [17, 19]) are prompted and applied in real circumstances. However, these methods lack sufficient underpinnings.

This is in contrast to other engineering problems, such as the global sensitivity analysis of explicit equations [20] which usually have analytical expressions. There are no theoretically precise answers to multicriteria decision problems. This is where the complexity of the multicriteria decision problem lies. As a result, weights calculated through any weight definition methods could not be deemed precise. Therefore, it only can be a reference for policymakers to take measures. To settle the multicriteria decision, the weight was regarded as a variable, which is uniformly distributed and whose bounds could be determined by getting results based on different weight definition methods in this paper.

This paper proposed a new method to conduct the risk assessment of debris flow, which is a probabilistic method. In Section 2, 6 influencing factors selected in this paper and the calculation steps of the proposed method are introduced. In Section 3, the detailed justification of the proposed method’s feasibility and rationality is demonstrated. In Section 4, the engineering geological conditions of the study area are introduced, and the risk assessment results are analyzed. The final section is the conclusions.

2. Study Area

Beichuan County lies in Sichuan province, China. After the 5.12 Wenchuan earthquake, the geological environment in Beichuan County is very fragile. There are many regional debris flows that were induced after a long period of rainfall. Therefore, this paper takes 72 debris flows in Beichuan County as the cases (Figure 1).

Mountains dominate the terrain in Beichuan County, and the maximum relative elevation difference can reach 1000 m. The overall trend is higher in the northwest and lowers in the southeast. The debris flow slope is generally greater than 25°, and the steeper gully slope even exceeds 50°. The high-mountain area accounts for 46.5% of Beichuan County’s area, due to its high altitude, harsh climate, sparsely populated areas, and being mostly covered by virgin forests. The area of low mountain area accounts for 16.0%, but geological disasters are less due to the slow terrain. The middle mountain, which accounts for 37.5% of the area, is an area prone to geological disasters. Because residents are concentrated on the banks of the Baicaohe River, the Sushuihe River, and the Bainigou River, it is also the area where most of the debris flows studied in this paper are located.

In Beichuan County, there are Paleozoic Sinian, Cambrian, Silurian, Devonian, Carboniferous strata, and Quaternary loose source materials of Mesozoic and Cenozoic. Debris flows mainly occur in the loose accumulation layers of the Quaternary system in the east and south. Alluvial deposits and flood alluvial deposits are widely distributed on both sides of the river and consist of gravel-clay and clay layers. Residual overburden is distributed in high, deep, and gentle areas on slopes and gullies and is mainly composed of clayey sand and gravel soil.

Beichuan County is located in the subtropical humid monsoon climatic zone. Simultaneously, because the rainfall is mostly concentrated in summer, July–September are also months with a high incidence of debris flow disasters. The rainfall is affected by the topography of the mountain area and shows a decreasing tendency in space from the southeast to the northwest.

3. Methods

3.1. Material

The debris flow is a major geological hazard in Beichuan County. Sorting levels of severity of debris flow could provide policymakers references when they tackle disasters. Many factors should be taken into consideration to assess the risk of debris flow. Referring to current researches, the assessment factors are simplified to 6. This paper takes the following aspects as assessment factors based on the development characteristics of debris flow in Beichuan County after the earthquake: the loose source material reserves, the basin area, the drainage density, the basin relative relief, the shifting bed proportion, and the main channel length. The term named “the loose source material reserves” is defined as the volume of the source of loose solids that form debris flows, which is composed of colluvium deposits, eluvium deposits, and channel deposits. It is one of the most direct indicators that affect the risk of debris flow. The term named “basin area” is defined as the catchment area surrounded by the drainage divide. The term named “drainage density” is defined as the ratio of the total length of a gully to the basin area. The term named “basin relative relief” is defined as the difference between the highest and the lowest points in the basin area. The term named “main channel length” is defined as the plane projection length of the main channel head to the channel end, and the term named “shifting bed proportion” is defined as the ratio of unstable channel length to the main channel length. Table 1 lists the survey data of the loose source material reserves, the basin area, the drainage density, the basin relative relief, the shifting bed proportion, and the main channel length of 72 debris flows located in Beichuan County. These data are from Wang [21].

3.2. Assessment Method

In this section, a method of sorting the risk of debris flow is given. Prior to this, two definitions are introduced.

Decision matrix: the matrix in the form of Table 1 is called the decision matrix Z, where Z = {Zij}i=1, 2, …,m;j=1, 2,n, m represents the number of samples (in this paper, the number of debris flows, that is 72), and n is the number of assessment factors (in this paper, n = 6).

Weights: the vector  = [, , …, ]T that satisfies is called the weight.

In order to eliminate the influence between the dimensions of the assessment factors, the decision matrix should be normalized first, i.e.,where max( ) means taking the maximum function and min( ) means taking the minimum function. Because the larger the value of the six assessment factors in this paper, the higher the risk grade, it is only necessary to use equation (1) to normalize the decision matrix. A positive relationship exists between the normalized value and the risk grade.

Multiply the normalized decision matrix by the weight to get the score vector u.where u = [u1, u2, …, um]T. The larger the score value ui, the higher the risk grade of the i-th debris flow. Determining the weight is the core problem of risk assessment of debris flow.

Several more commonly used methods of weight determination are described hereinafter.

3.2.1. Equal Weight Method

The equal weight method is the most widely used in practice. For some problems where it is often difficult for people to judge the importance of assessment factors, the equal weight method is needed. The equal weight method can be expressed as

The equal weight method can be regarded as a special case of the analytic hierarchy process. If the importance of the assessment factors is considered the same in the analytic hierarchy process, the equal weight method results can be obtained.

3.2.2. Maximum Deviation Method

The maximum deviation method determines weights based on a decision matrix [13].

Calculate the deviation according to the following equation:

The total deviation is defined as

In the maximum deviation method, it is considered that the weights should maximize the deviation of all assessment factors. That is,

Using the Lagrangian method, the solution yields weights of

3.2.3. Entropy Method

The entropy method is also a method to determine the weight based on the decision matrix [14].

The entropy vector H of the evaluation factor is defined as

According to equation (8), the entropy weight of the assessment factor is

It can be seen from equations (8) and (9) that there is a negative relationship between the entropy value and the entropy weight, which indicates that the more effective the information of the assessment factor, the more significant the assessment factor. Conversely, the larger the entropy value is, the smaller the entropy weight is, and the less important the assessment factor is. The entropy weight reflects that the function of the assessment factor in objective information is the objective weight.

The weights calculated by different methods are entirely different. How to combine the weights obtained by different methods is a classic problem. A common approach is to average the weights from different methods. However, this approach is too simple and lacking insufficient evidence. Hereinafter, we will propose a new method to synthesize the weights obtained by different methods.

3.2.4. Proposed Method

According to the definition of weight, any vector  = [, , …, ]T satisfying can be used as a weight. The weights of assessment factors are obtained from different data analysis perspectives, including the equal weight method, the maximum deviation method, and the entropy method. Through the three methods, the upper and lower limits of weights of assessment factors can be obtained. The proposed method drew on the probability analysis idea to obtain weights, which regarded the weight as a random variable with uniform distribution within the boundary. To sum up, a method combining the weights of different methods through probability analysis and calculating the score vector is proposed in this paper. Its implementation steps (Figure 2) are as follows.Step 1: The weights of assessment factors are calculated using the equal weight method, the maximum deviation method, and the entropy method.Step 2: Calculate the upper and lower limits of weights of the assessment factor.Step 3: Calculate the normalized decision matrix.Step 4: Generate the random variable sample W = {}i = 1, 2, …, l; j = 1, 2, …, n that satisfies the uniform distribution of weight within the boundary, where l represents the number of samples.Step 5: Calculate the scoring matrix of each random weight sample according to equation (2). It should be noted that since the random sample of each weight does not meet the definition of weight, it is necessary to normalize each random sample.Step 6: Take the average value of l scoring matrices to obtain the final scoring matrix.

Steps 3–6 above are often referred to as the Monte Carlo simulation strategy. It is worth noting that the Monte Carlo simulation strategy is only used to generate random weights within the boundary, while the debris flow decision matrix uses field statistics. There are two essential questions that need to be answered about the proposed approach. One is whether the method converges, and the other is how much l should be taken. These two issues will be discussed below.

4. Results

4.1. Calculation Results of the Three Methods

Figure 3 is the calculation results of the equal weight method, the maximum deviation method, and the entropy method. The maximum deviation method and the entropy method are calculated by Matlab [22]. As shown in Figure 2, there is a large difference between the calculations of the three methods. The entropy method considers that loose source material reserves are much more important than other assessment factors, while the maximum deviation method considers that the importance of the basin area is more significant than other assessment factors. The results obtained by several methods are consistent with only one evaluation factor (i.e., the main channel length).

Because the essence of the entropy method is to reflect the degree of disorganization, the greater the randomness and degree of disorganization, the more critical the assessment factor. The loose source material reserves have a large difference among the 72 debris flows in this paper. The maximum is 1967.9 × 103 m3, and the minimum is only 0.04 × 103 m3. Therefore, the weight of the loose source material reserves in the entropy method is as high as 0.44. However, this result overemphasizes the importance of loose source material reserves.

The maximum deviation method considers the deviation between the basic data and its average value. Although the sample range of loose source material reserves is huge, the deviation is smaller than that of the basin area. Moreover, the weight of the basin area is assigned 0.21 in the maximum deviation method. However, of the 72 debris flows in Beichuan County, the main ones are small- and medium-sized debris flows. The deviation of each debris flow assessment factor is small, and the calculated weights have little difference, which fails to reflect the characteristics of assessment factors. Because of the advantages and disadvantages of the two methods, the proposed method within the boundary conditions considers the characteristics of assessment factors and avoids overemphasizing the importance of a specific factor.

4.2. Convergence Analysis

The value range of 6 assessment factors, respectively, obtained from Figure 3 is

Random samples of different l(31, 63, 125, 251, 501, 1000, 1995, 3981, 7943, 15848, 31622, 63095, 125892, 251188, 501187, 1000000, 1995262) were generated in Matlab 2018a. Figure 4 shows the change curve of the score values of some debris flows with the number of samples l. As shown in Figure 4, with the increase of l, the scoring value gradually tends to be stable, indicating that the proposed method is convergent.

There was some fluctuation in the sorting of the 5th and 63rd debris flow samples in the calculations (Figure 5). However, with the increase of l to 63095, the sorting of the risk level has been entirely stable. The limit for l seems to be reduced to a certain range between 31622 and 63095. However, even if 63095 is taken, the calculation time is only 0.3510 seconds, and 63095 is completely feasible.

4.3. Comparison of Results

The weights calculated by equal weight method, maximum deviation method, entropy method, and average method are substituted into equation (2) to sort the debris flow risk. Table 2 shows the detailed sorting results of the 4 methods. It can be seen that the sorting results of these methods are vastly divergent. Figure 6 more intuitively shows the differences in the sorting of debris flow risk by different methods, which is calculated from Table 2. It reflects the compliance degree of the five methods on the same debris flow’s sorting result (for example, #30 debris flow sorts first among all methods in terms of risk, so its compliance degree is 100%; #25 debris flow only ranks 11th in this method, so its compliance degree is 0%). The results of several methods on the top three dangerous debris flows (i.e., #30, #55, #59) are exactly the same. From the field survey data (Table 1), these three debris flows are highly developed and potentially dangerous.

The correlation matrix is defined as the ratio of the same number of sorting to the total debris flows. Table 3 is calculated from Table 2 and shows the correlation matrix of the 5 methods. As shown in Table 2, the correlation coefficients of the equal weight method, the maximum deviation method, and the entropy method are all very small due to the difference in weights. Although the average method of the weight obtained by the 3 methods also considers the differences between different methods, it is only a special case of the proposed method. Take the highest dangerous debris flow (i.e., #30) as an example. The loose source material reserves are 1967.9 × 103 m3, and the loose source materials in the channel are abundant. Simultaneously, the basin area is 15.7 km2, and the water source in the basin is sufficient. The considerable basin relative relief and shifting bed proportion reflect the high erosion degree of the channel, and the sediment and loose source materials along the way are adequately recharged. Therefore, from the perspective of hydrological and topographical conditions, #30 is more dangerous.

5. Discussions

5.1. Risk Sorting Results

As can be seen from Table 2, the equal weight method and entropy method consider #11 to be more dangerous than #69. However, the 50-year debris flow scale in #11 is only 16.2 × 103 m3, and the debris flow scale in #69 is 184.61 × 103 m3. This result is not credible. Although the loose source material reserves of #11 have reached 1500 × 103m3, there are not enough water sources in the basin to provide conditions for the development of debris flows, and low basin relative relief is not conducive to the occurrence of large-scale debris flows. This error is due to the excessive consideration of #11 loose source reserves when using the single weight method. In the proposed method, this error is avoided.

Since there are 72 debris flows in this paper, typical debris flows #30 and #42 are taken as examples in discussing the assessment results. Figure 7 shows the activity of the #30 and #42 debris flows. The formation of debris flow is determined by topographic, geomorphic, and water source sources and loose source materials. The #30 is the debris flow with the highest risk score in the proposed method. It can be seen from Table 1 that the loose source material reserves of #30 are 1967.9 × 103 m3, while the average value is only 191.3 × 103 m3. Water source conditions are the driving factors that induce debris flow determined by the basin area. The basin area of #30 reaches 15.7 km2, while the average basin area is only 6.2 km2. Therefore, for #30, the loose source materials are extensive, the water source in the basin is sufficient, and the necessary conditions for debris flow are available. The basin relative relief exceeds 1.8 km, and the larger drainage density and shifting bed proportion reflect the high erosion degree of the channel and sufficient sediment replenishment along the way. Therefore, the proposed method which regards #30 as the most dangerous debris flow is reasonable.

The loose source material reserves of #42 are 106.5 × 103 m3, and the basin area is only 0.3 km2, which are both much smaller than the average. The terrain of #42 is relatively flat, with the basin relative relief of 0.52 km. The drainage density and shifting bed proportion are far below the average. Consequently, #42 is relatively safe in terms of terrain, hydrological, and provenance conditions.

5.2. Further Discussion of the Proposed Method

The proposed method is to sort the risk grade of debris flow from the perspective of probability. The purpose of this method is to avoid the problem of determining weights. Weight determination has always been an essential issue in decision science. Although many weight determination methods have been put forward, there is no absolute criterion to judge the rationality of the obtained weight. It is because there will always be some deviation existing when using a certain weight to sort the risk grade of debris flow. In contrast, the sorting method in this paper can be regarded as a reasonable way considering all possible weights.

It is worth noticing that the proposed method is different from the study of Zhang et al. [6]. The proposed method just sorted the risk grade of debris flow and did not categorize the debris flow. In the study of Zhang et al. [6], a classification standard was artificially given, and a series of calculations were made to determine the risk grade of debris flow, which should be more appropriately called a classification issue. Such an approach is reasonable, but the classification standard needs to be subjectively given which often will be unconvincing. The study of classification standards is beyond the scope of this paper.

It should also be noted that this paper only uses the equal weight method, the maximum deviation method, and the entropy method to determine weights. Using more methods will be more beneficial in determining the value range of weight in practical application.

In a recent study, Alireza [23] reviewed the development of fuzzy multicriteria decision problems in the past decade. The sorting of fuzzy numbers is still a big challenge, but the sorting of real numbers has been well solved. In other words, it is obvious that the difficulty of debris flow risk assessment lies only on the determination of weight. However, different sorting models, such as the grey target model [24] or the fuzzy assessment model [25], have little difference in terms of their essence with equations (1) and (2). For classification problems, after the classification standard is given, the classification boundary can be substituted into equations (1) and (2) to obtain the boundary score value. And the sample grade can be determined by comparing the size of the sample score value and the boundary score value.

Although some studies (such as Wang et al., [17]) claim that classification standards also have fuzziness, the fuzziness of the assessment factor should be the main factor. Therefore, in future research, we will use fuzzy number operation rules to consider debris flow risk assessment under a fuzzy environment.

6. Conclusion

A new method for sorting the risk grade of debris flow is proposed. The new method regards weight as a variable, which is uniformly distributed and whose boundaries could be determined by the results calculated by the equal weight strategy, the maximal deviation method, and the entropy method.(1)This new method is adopted to sort the risk grade of 72 debris flow gullies in Beichuan. Furthermore, the result proves that the proposed method is convergent and finds out the minimum sample volume needed by the proposed method.(2)The new method considers the characteristics of the assessment factor and avoids the defects of the traditional weight method. Therefore, this method is more objective and accurate.(3)The comparison between the assessment results of the new method and the actual results verifies the reliability of the new method.

Our suggestions for further study are to apply the proposed method to sort the risk grade of debris flow in fuzzy environments.

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

The authors thank the Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant nos. KJQN201901221, KJQN201901239, KJQN201901241, KJQN202001218, and KJQN202001219) and the General Program of Chongqing Natural Science Foundation (Grant no. cstc2019jcyj-mscm1865).