Abstract

Predicting landslide displacement is of great significance in geotechnical engineering. An iteration-based combined prediction method was proposed for predicting the landslide displacement in this paper. Firstly, the landslide displacement was predicted by 10 latest multivariable grey models, and then the final landslide displacement prediction value was obtained through an iteration-based combined strategy. Concurrently, the reliability of the quadratic programming-based combined prediction method (QPCPM) and the iteration-based combined prediction method (ICPM) was rigorously proved in this paper. In addition, the inapplicability conditions of the optimal weight-based combined prediction method (OWCPM) were pointed out. ICPM could be regarded as a simplified version of QPCPM. The Bazimen and Baishuihe landslides in the Three Gorges Reservoir area of China were used as numerical examples to elaborate the performance of ICPM. This paper also demonstrated the reliability of ICPM by considering the effects of rainfall and reservoir water level on landslide displacement. Overall, ICPM features in simple and easy calculation and has rosy engineering application prospects.

1. Introduction

The landslide time prediction is a common concern of geological engineers [1]. However, it is beset with difficulties to solve this problem at present. Although the inverse velocity method [2] has achieved some success in the landslide time prediction, its application effect is affected by the operator’s operant level, and the theory also has some implausible aspects [3]. Given this, people turn to predict the landslide displacement to realize the warming of landslides indirectly. Prediction of landslide displacement plays a critical role in reducing the damage of landslide disaster. In recent years, prediction methods based on mathematical models have been popular topics in studying landslide displacement prediction [4, 5].

Many observations indicate that the landslide displacement is nonlinear and nonstationary and is related to groundwater and rainfall [6]. A support vector machine-based prediction method considering nonlinear characteristics of landslide displacement was developed by Feng et al. [7], and the method has been used to analyze the high rock slope of the permanent ship-lock of the Three Gorges project and the area deformation at depth in the Bachimen landslide in Fujian Province, China. In order to consider the instability of landslide displacement, a prediction method based on modified ensemble empirical mode decomposition and extreme learning machine was proposed [8]. The method showed promising results in predicting the landslide displacement of the Baishuihe landslide in the Three Gorges Reservoir area of China. A landslide displacement prediction method based on a multivariate chaotic extreme learning machine was developed for the sake of considering the impact of rainfall and groundwater on the displacement of the landslide [9]. These methods have successfully predicted the displacement of the Bazimen and Baishuihe landslides in the Three Gorges Reservoir area of China. Besides, a variety of landslide displacement prediction methods based on mathematical models have been proposed. As a case in point, a landslide displacement prediction method based on the neural network was proposed by Du and Yin [10] and obtained good prediction accuracy. Simultaneously, the prediction model based on the extreme learning adaptive neuro-fuzzy inference system [11] and the novel kernel extreme learning machine [12] can also be applied to landslide displacement. Undoubtedly, these prediction methods could achieve great results on specific issues. The essence of these methods is to extrapolate the mathematical model obtained by approximating the simulated value of the landslide displacement and the measured value to realize the prediction [13, 14]. Provided that people use different methods to predict landslide displacements simultaneously, they tend to trust the predicted values with the smallest simulated error. However, this approach also faces a dilemma that the better the simulated effect does not necessarily mean the better the predicted effect. It is incredibly challenging to select the best landslide displacement prediction model based on simulated results.

Undoubtedly, an accurate prediction method should have smaller simulated errors. The prediction methods with similar simulated errors have different performance on different samples, and there is no perfect prediction method that applies to all samples. In order to ensure the reliability of the landslide displacement prediction method, it may be a feasible way to predict the landslide displacement from the perspective of combined prediction [15, 16]. Unfortunately, there are relatively few studies on landslide displacement prediction based on the combined prediction strategy. A landslide displacement prediction method based on the combined method was developed by Li et al. [17]. The method was used to predict the landslide displacement of the Huanglongxicun landslide and the Saleshan landslide. However, due to the defects of the prediction method itself, the weights will appear negative, so this is one of the issues that this paper needs to solve.

This paper developed an iteration-based combined prediction method (ICPM) for landslide displacement prediction based on the above discussions. The feasibility and advantages of this approach have been demonstrated with rigorous mathematical formulas. The structure of this paper is as follows: in Section 2, ICPM developed in this paper was introduced; in Section 3, the Bazimen landslide and the Baishuihe landslide in the Three Gorges Reservoir area of China were used as numerical examples to present the effectiveness of the proposed method in this paper; finally, the conclusions were drawn.

2. Methods

In this section, ICPM for landslide displacement prediction was proposed. Firstly, two conditions for constructing a combined prediction method were given. Condition 1: the performance of the mathematical models used to establish the combined prediction method should not be too far apart because it is of little significance to establish a combined method by using models with significant performance differences. Condition 2: due to the fact that complex combined strategies tend to lose practicality, combined strategies should be as simple as possible. Then, according to condition 1, this paper selected 10 different multivariable grey models to predict landslide displacement. The multivariate grey model is an essential type of the prediction method compared with neural networks and other prediction methods. The multivariate grey model has low complexity in the calculation and stable predicted results; thus, it can be applied to engineering practice. According to condition 2, a new combined method with relative low complexity in the calculation is put forward in this paper.

2.1. 10 Multivariable Grey Models

The grey prediction model is fundamental initially put forward by Deng and evolved from the original unvaried prediction model to the multivariate prediction model. Here, we introduced the original multivariate grey model GM (1, N).

Consider the observed data, , where is the dependent variable, are the independent variables, and

For the first-order accumulation of the observation set, equation (2) is obtained:where

X1 can be described by the following differential equation [18]:where bi(i = 1, 2, …, N) is the coefficient, which can be estimated by least squares, i.e.,where

The initial value conditions are considered:

The solution of equation (4) can be obtained, and the fitted value X1 can be obtained by substituting the parameters bi(i = 1, 2, …, N), and the fitted value of the original sequence can be obtained by subtractive damage reduction:

The GM (1, N) has a clear principle and simple calculation, but it has limited accuracy in predicting complex time series in reality [19]. For this reason, some improved multivariate grey models have been proposed. The 10 latest multivariate grey prediction models were used in this paper to predict landslide displacement. Detailed information can be found in Table 1. The 10 multivariate grey models in Table 1 were improvements to GM (1, N).

To ensure that these 10 multivariate grey models (DGM, FDGM, GMCG, GMCN, GMCPQ, GMCT, NGBMC, NGM, OBGM, and RDGM) meet condition 1, the performance of these 10 prediction models is firstly compared through published data in literature references. The source of the reference data is given in Table 2. The data for all examples are divided into two parts: the simulated samples and the predicted samples. Model performance is measured by mean absolute percent error (MAPE), i.e.,

The smaller the MAPE is, the better the performance of the model will be.

The programs of the 10 multivariate grey models have been compiled on MATLAB and applied to the 9 cases in Table 2 and then combined equation (9) to measure the performance of each model. Figure 1 shows MAPE of the simulated samples and the predicted samples of the 10 models in 9 cases. As can be seen from Figure 1(a), the MAPE of the simulated samples of the 10 models is all within 6%. The difference is not large. As can be seen from Figure 1(b), the MAPE of the predicted samples is significantly larger than that of the simulated samples, but the differences between the models are still small. Figure 2 shows a box plot of the 10 models. It can be seen from Figures 1 and 2 that although the performance of these 10 models has a certain gap, the differences are generally small.

2.2. The Combined Strategy

A combined strategy for predicting landslide displacement was proposed [17], which is called the optimal weight strategy. The strategy is derived as follows:

Using m different models to predict landslide displacement and weighting with , the residual of the combined method can be written as follows:

From equation (10), equation (11) can be obtained:

Further, the residual sum of squares is as follows:

According to the Lagrange multiplier method, can be solved by equation (12), i.e.,

Equation (13) is the combined strategy of Li et al. [17]. However, the equation does not restrict weight’s positive and negative properties. The calculated value of equation (13) will appear negative, which does not meet the definition of the weight. A negative weight will also result in a large predicted error (refer to Section 2.3 or Section 3).

If the condition that “weight” cannot be negative is set in equation (12), the following mathematical problem can be obtained:

Equation (14) is a quadratic programming problem that often needs to be solved by numerical methods. This paper will give a simple combined strategy to avoid the complicated quadratic programming process. Firstly, a definition and a theorem were given.

Definition 1. When  = 1/m, the method is called the equal weight-based combined prediction method.

Theorem 1. Ja is set to the sum of squares due to error (SSE) of the equal weight combined prediction method, and then Ja < Jmax must be true, where Jmax is the maximum SSE in all single prediction models.

Proof. According to the definition of the equal weight-based combined prediction method, equation (15) can be obtained:From the Cauchy formula [28],That is,According to Theorem 1, we can design an iteration-based combined prediction strategy. The basic idea is to replace the value of the prediction model corresponding to Jmax with the predicted value corresponding to Ja until the algorithm converges (SSE is less than 0.001), i.e.,In this paper, the combined prediction methods established by equations (13)–(15) are, respectively, called the optimal weight-based combined prediction method (OWCPM), quadratic programming-based combined prediction method (QPCPM), and iteration-based combined prediction method (ICPM).

2.3. Mathematical Properties of the Method

Property 1. The simulated samples’ SSE is least better than the worst single model.

Proof. Property 1 is obvious, so it will no longer be proved.

Property 2. It is impossible to estimate the upper bound of SSE in OWCPM’s predicted samples. However, QPCPM and ICPM are at least better than the worst single model.

Proof. Set the weights of m models to be [1, 2, …, m].OWCPM could not guarantee that the inequality sign is established. However, QPCPM and ICPM could ensure that the inequality sign is established.
Property 1 and Property 2, respectively, answer the two questions of “why the combined prediction method is used” and “what advantage using combined prediction method has.” From Property 1, it can be further obtained that SSE of OWCPM for simulated samples is the smallest. Unfortunately, Property 2 directly indicates that OWCPM is unreliable. Although QPCPM and ICPM are reliable, the latter is much simpler. ICPM could be regarded as a simplified version of QPCPM. In summary, engineers have the confidence to use ICPM in actual engineering.

2.4. Calculation Flow Chart

The calculation flow chart of ICPM is shown in Figure 3. The third step (that is, “check for abnormal conditions”) is mainly to check whether there are different models with the same simulated values in the 10 multivariate grey models. If different models have the same simulated values, then using only one of these models can reduce computation (see Section 3.2 for a detailed example).

3. Numerical Examples

The Bazimen and Baishuihe landslides in the Three Gorges Reservoir area in China are two examples with significant focus in the study of landslide displacement prediction. Figure 4 shows the geographic location of the Bazimen landslide and the Baishuihe landslide in the Three Gorges Reservoir area. This paper takes these two examples to further illustrate the effectiveness of ICPM.

3.1. Bazimen Landslide

Located on the right bank of Xiangxi River, Xiangxi village, Zigui County, Hubei province, the landslide is 0.8 km away from the river mouth and 38 km away from the Three Gorges Dam. The landslide is 350 m long, 350∼500 m wide, 30 m in average thickness, and 4 × 107 m3 in volume. It is an accumulation landslide, and the bank slope is the counter-inclined slope.

The Bazimen landslide was formed in earlier days and belonged to the old landslide. According to the drilling, there are two layers of slip zones in the front, and the upper part is the secondary slip surface. The main slip zone is continuously distributed between the bottom of the slip body and the landslide bed interface, and the inclination of the slip surface is steeper in the middle and rear, about 20° to 30°. The front part is gently inclined and slightly reversed, and its thickness varies from 0.9 to 3.6 m. The landslide slides along the soil-rock contact surface or the top surface of the flood deposit. The thickness of the slip zone’s overall trend is that the front edge and upper part are thin, and the lower part is thick. The slip material is breccia mass, with partial surface being polished and scratched, and the compositions are mainly sandstone and mudstone, filled with silty clay.

The monitoring data of ZG110 and ZG111 monitoring points of the Bazimen landslide from November 2012 to November 2014 were selected (the topographic map of the landslide is shown in Figure 5, and the geological profile is shown in Figure 6). The data unit is on a monthly basis with a total of 25 sets of monitoring data. The first 18 sets of data were used as simulated samples and considered the influence of rainfall and reservoir water level (the grey correlation degree between the landslide displacement data and the rainfall data at ZG111 is 0.60, and the grey correlation degree between the landslide displacement data and the reservoir water level data is 0.94; the grey correlation degree between the landslide displacement data and the rainfall data at ZG110 is 0.61, and the grey correlation degree between the landslide displacement data and the reservoir water level data is 0.93, so the correlation degree calculation results show that it is feasible to use rainfall and reservoir water level as the influencing factors of landslide displacement). The detailed data are shown in Figure 7. It can be seen from Figure 7 that the flood season is concentrated from May to September each year. In July 2014, the rainfall reached the maximum in the two years, and the reservoir water level has a large increase in July and August each year. The annual cumulative displacements measured by GPS monitoring points ZG110 and ZG111 in 2013 were 14.003 mm and 15.075 mm, respectively, and the annual cumulative displacements measured in 2014 were 11.852 mm and 14.544 mm, and their average velocity was 1.077 mm/month and 1.234 mm/month, respectively. The cumulative displacement of ZG110 on the frontal part of the landslide grew more quickly than that of ZG111 on the upper part. Moreover, it can be seen from Figure 7 that the cumulative displacement of the monitoring point is affected by the reservoir water level, and there is a sudden increase in the months when the reservoir water level and rainfall change suddenly, with a time interval of 1a.

Figure 8 shows the calculation results of 10 multivariate grey models. It can be seen from the figure that for most prediction models, the prediction results are quite different from the actual values. Among them, the predicted result of FDGM far exceeds the actual value. Only RDGM model’s predicted values at the two GPS monitoring points are closer to the actual values. Tables 3 and 4 give the MAPE values of simulated samples and predicted samples of the 10 multivariate grey models, respectively. Tables 5 and 6 give SSE of the 10 multivariate grey models. The FDGM model shows better performance in the Bazimen landslide, while the NGBMC model shows better performance in the Baishuihe landslide.

The calculation results of the 3 combined prediction methods are shown in Figure 9. MAPE and SSE of the 3 combined methods are shown in Tables 7 and 8 according to Figure 8. From the comparison of the simulated samples’ results in Tables 3 and 7, it can be seen that the 3 combined models obtain better MAPE than most single models. Only MAPE of FDGM’s ZG111 point is 0.488 better than QPCPM. Comparing Tables 4 and 7, MAPE of OWCPM is larger than that of all the single models (the reason for this problem is that the weight calculation is unreasonable, as shown in Figure 10, and the theoretical description is shown in Section 2). QPCPM and ICPM can produce better results than a single model. It can be seen from Table 8 that 2 methods have obtained relatively reliable prediction results.

3.2. Baishuihe Landslide

The Baishuihe landslide is located in Zigui County of Three Gorges Reservoir area and on the south bank of the Yangtze River, 56 km away from the Three Gorges Dam. The landslide’s front edge submerges into the Yangtze River, the left and right sides are bounded by the bedrock ridge, and the back edge is bounded by the boundary of the rock and soil. The landslide’s main slip direction is about NE20°, with the NS of 600 m long, the EW of about 700 m wide, the average thickness of 30 m, and the volume is about 1.26 × 107 m3. The landslide’s front and back edges are relatively steep, and the middle part is relatively gentle. The overall slope is in a “staircase” shape toward the Yangtze River, with an overall slope of 30°.

The Baishuihe landslide is an old landslide and has witnessed many large deformations in history. In June 2003, a lateral crack of more than 300°m emerged in the eastern front of landslide. In July 2004, the landslide deformation destroyed 21 houses. From August 2005 to August 2006, sinking and cracking occurred in many places on the landslide surface. On June 30, 2007, a large-scale collapse occurred at the back edge of the early warning zone. In August 2009, the cracks in the western boundary of the landslide spread in a feather-like manner. At present, the landslide presents the characteristics of traction deformation and failure as a whole. The deformation is mainly concentrated in the early warning zone at the front of the landslide, and the displacement of the monitoring point at the rear is relatively stable.

The monitoring data of the Baishuihe landslide from January 2010 to December 2011 are selected (the landslide topographic map is shown in Figure 11, the geological profile is shown in Figure 12, and the monitoring data are shown in Figure 13). The monitoring point’s cumulative displacement continues to increase, and the maximum displacement increase rate occurs during the reservoir water level’s falling to the bottom level. The annual increase value of cumulative displacement of ZG118 is 156.1 mm and 146.34 mm, respectively, and the annual increase value of cumulative displacement of XD01 is 263.42 mm and 185.37 mm, respectively. It can be seen that the cumulative displacement of XD01 increases faster than that of ZG118. With the long-term adjustment of the reservoir water level, the landslide has undergone long-term stress adjustment. The landslide gradually adapts to the new reservoir scheduling mode, and its displacement growth tends to be stable.

The former 18 sets of data are still used to establish the model (the grey correlation degree between the landslide displacement data at ZG118 and the rainfall data is 0.71, and the grey correlation degree between ZG118 and the reservoir water level data is 0.96; the grey correlation degree between the landslide displacement data and the rainfall data at ZG110 is 0.70, and the grey correlation degree between ZG110 and the reservoir water level data is 0.96. The correlation degree calculation indicates that it is feasible to use rainfall and reservoir water levels as factors influencing landslide displacement). Figure 14 provides the results of 10 multivariate grey prediction models. Figure 15 provides the prediction results of 3 combined methods. The simulated and predicted samples’ MAPE and SSE in Tables 36 can be calculated from Figure 14. MAPE and SSE of the three methods derived from Figure 15 are shown in Tables 7 and 8, respectively. It is worth noting that the 10 multivariate grey prediction models obtained very different results at 2 GPS monitoring points. In the prediction results of ZG118, the FDGM model showed promising results, and the NGM model showed good performance in the prediction results of XD01. However, the overall effect of using a single method is not satisfactory. It can be seen from Table 7 that MAPE value of OWCPM is the smallest at ZG118, and the minimum MAPE at XD01 is QPCPM. The difference is also confirmed in Figure 15.

Based on the above two examples, it can be seen that it is beset with difficulties to only rely on a single model to predict landslide displacement. For example, in the Bazimen landslide, FDGM achieves the best performance on the simulated samples. According to the general standard, the prediction results of FDGM are considered reliable. Unfortunately, the calculation of prediction samples shows that the prediction performance of FDGM is not good. OWCPM is not recommended in practice because of its theoretical flaws. In terms of calculation results, although QPCPM is better than ICPM, the relative level of the prediction performance of the two methods is at random. In the Bazimen landslide simulation samples, the two combined methods’ MAPE and SSE are very close. However, MAPE of QPCPM is slightly lower than that of ICPM in prediction samples (MAPE of QPCPM is 2.555% and 4.412%; MAPE of ICPM is 2.580% and 5.028%; SSE of QPCPM is 28.051 mm2 and 58.271 mm2; and SSE of ICPM is 28.272 mm2 and 76.464 mm2). In the Baishuihe landslide, the prediction samples based on QPCPM have a lower MAPE in ZG118 than ICPM, while the prediction samples’ results at XD01 are the opposite. Therefore, the iteration-based combined prediction method is still recommended in terms of computational complexity.

4. Discussion

Prediction of landslide displacement plays an essential role in reducing landslide disaster damage. A large number of landslide displacement prediction models have been proposed. However, as Intrieri et al. [29] pointed out: “Nonetheless, this field has experienced a great pulse of research development in the very recent past, even to the point that the literature seems to be quite redundant and not conclusive.” Many papers have proposed the prediction methods of landslide displacement prediction, but their effectiveness is obtained only by one example and lacks explicit elaboration for the method’s mathematical properties. The combined strategies are elaborated in this paper. Since the optimal weight strategy does not restrict the weights’ positive and negative properties, this paper adds the restriction conditions in equation (14) to the optimal weight strategy to obtain the quadratic programming problem. Simultaneously, to avoid a complicated solution process, the value of the prediction model corresponding to Jmax is replaced with the predicted value corresponding to Ja, and ICPM can be obtained.

Artificial neural networks are often used for landslide displacement prediction. A landslide displacement prediction method based on the ANNs-switched method was proposed by Lian et al. [30]. This method considers rainfall and reservoir level elevation as parameters as in this paper. MAPE of predicted samples of the ANNs-switched prediction method in the Bazimen landslide is 1.8940%. Its predicted performance is higher than 7 of the 10 latest grey multivariate prediction models. Only the performance of GMCN, NGM, and OBGM is slightly better than the ANNs-switched prediction method.

It can be concluded from Tables 7 and 8 that regardless of whether MAPE or SSE is considered, the performance of the prediction method using the combined strategy is better than that of the ANNs-switched prediction method. However, due to OWCPM’s theoretical flaws, the results in the predicted samples of the Bazimen landslide are not satisfactory. At the same time, it was also found that the MAPE and SSE of QPCPM and ICPM are also very close. Therefore, in the actual use of these two prediction methods, very similar results can be obtained. Nevertheless, QPCPM has too much calculation and is not a convenient method. As a simplified QPCPM, ICPM greatly simplifies the calculation and reduces the complexity. Its mathematical reliability is proven in Section 2.3.

A suitable landslide displacement prediction method must have a relatively minor simulated error and predicted error. However, the prediction method applicable to all landslides does not exist. The same prediction method has different performances on different landslide samples. To be precise, ICPM cannot be regarded as an accurate method but a method with high reliability. It is worth noting that the proposed method cannot be applied for regions where soils may exhibit considerable strain softening, which may trigger sudden rapid dynamic landslide movement such as that occurred during the Vaiont landslide [3133].

5. Conclusions

In this paper, an iteration-based combined strategy for landslide displacement is proposed. The conclusions are as follows:(1)ICPM is based on the 10 latest multivariate grey prediction models’ calculation. The proposed method obtains the predicted value of landslide displacement by simple iterative calculation, which has the advantages of simplicity and high reliability.(2)It proved that OWCPM is challenged to guarantee the prediction performance, and it is not recommended to use it in practice. Compared with QPCPM, ICPM has low complexity and is convenient for engineers to adopt.(3)The landslide displacement prediction cases indicate that ICPM is useful, and the method is expected to be applied in engineering.

Data Availability

The data used to support the findings of this study are included within the article.

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 General Program of Chongqing Natural Science Foundation (grant no. cstc2019jcyj-mscm1865).