Next Article in Journal
Geometric Analysis of Sun-Assisted Lunar Transfer Trajectories in the Planar Bicircular Four-Body Model
Next Article in Special Issue
Study on the Change in Vegetation Coverage in Desert Oasis and Its Driving Factors from 1990 to 2020 Based on Google Earth Engine
Previous Article in Journal
Adaptive Visual Saliency Feature Enhancement of CBCT for Image-Guided Radiotherapy
Previous Article in Special Issue
Lithological Mapping of Kohat Basin in Pakistan Using Multispectral Remote Sensing Data: A Comparison of Support Vector Machine (SVM) and Artificial Neural Network (ANN)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning Approaches for Slope Deformation Prediction Based on Monitored Time-Series Displacement Data: A Comparative Investigation

1
School of Engineering and Technology, China University of Geosciences (Beijing), Beijing 100083, China
2
China Institute of Geo-Environment Monitoring, Beijing 100081, China
3
Center for Hydrogeology and Environmental Geology, China Geological Survey, Baoding 071000, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2023, 13(8), 4677; https://doi.org/10.3390/app13084677
Submission received: 7 March 2023 / Revised: 26 March 2023 / Accepted: 3 April 2023 / Published: 7 April 2023
(This article belongs to the Special Issue Applications of Machine Learning on Earth Sciences)

Abstract

:
Slope deformation prediction is one of the critical factors in the early warning of slope failure. Establishing an accurate slope deformation prediction model is important. Time-series displacement data of slopes directly reflect the deformation characteristics and stability properties of slopes. The use of existing data analysis approaches, such as statistical methods and machine learning algorithms, to establish a reasonable and accurate prediction model based on the monitored time-series displacement data is a common solution to slope deformation prediction. In this paper, we conduct a comparative investigation of machine learning approaches for slope deformation prediction based on monitored time-series displacement data. First, we established eleven slope deformation prediction models based on the time-series displacement data obtained from seven in situ monitoring points of the Huanglianshu landslide using machine learning approaches. Second, four evaluation metrics were used to comparatively analyze the prediction performance of all models at each monitoring point. The experimental results of the Huanglianshu landslide indicated that the long-short-term memory (LSTM) model with an attention mechanism and the transformer model achieved the highest prediction accuracy. The comparative analysis of model characteristics suggested that the Transformer model is better adapted to predict nonlinear landslide displacements that are affected by multiple factors. The drawn conclusion could help select a suitable slope deformation model for early landslide warnings.

1. Introduction

Landslides are one of the most serious geological disasters, which frequently occur worldwide and seriously affect people’s property and environmental safety. Therefore, slope deformation monitoring and early warning are of great significance [1,2]. Slope deformation prediction is one of the key factors in the early warning of slope failure. Establishing an accurate prediction model for forecasting slope deformation is critical. Slope deformation prediction models are mainly divided into physics-driven models and data-driven models [3]. Data-driven models, which are widely used and possess higher prediction accuracy, fully consider the complexity and nonlinearity of the deformation characteristics of slopes. [4,5].
The time-series forecasting method [6,7] is a typical data-driven approach for slope deformation, which predicts the future by learning the potential laws in the historical time-series displacement data of slopes. The time-series displacement data of slopes, which can directly reflect the deformation characteristics of slopes, are a series of displacement observation values obtained by the slope field monitoring equipment at a certain time interval. In practical engineering applications, field monitoring devices of the slopes are expensive, and the monitoring process is influenced by many factors [8,9]. Data analysis approaches such as statistical methods and machine learning approaches have been widely applied to various engineering problems with outstanding achievements in recent years. Therefore, exploiting the existing time-series forecasting approaches to establish a reasonable and accurate prediction model based on the monitored displacement data is a common solution to slope deformation prediction.
Currently, time-series forecasting methods mainly include three types: statistical approaches, traditional machine learning approaches, and deep learning approaches. The most popular statistical model is the autoregressive integrated moving average (ARIMA). The ARIMA model is mainly used to predict linear time-series data with accurate forecasting performance on short-term time series [10,11]. However, although the ARIMA can be adopted to achieve a highly accurate linear prediction, it cannot accurately forecast nonlinear time series. The ARIMA model needs to be combined with machine learning methods such as neural networks to improve the prediction performance on nonlinear time series [12,13]. For example, Shi et al. [14] pointed out that prediction models combining machine learning approaches and ARIMA outperformed a single ARIMA model.
The slope system is a complex nonlinear dynamic system, and the stability and deformation of slopes are affected by many uncertain factors, such as rock and soil properties, geomorphological features, rainfall, human engineering activities, and temperature changes. The interaction of these features leads to the complexity, ambiguity and randomness of slope deformation. Therefore, simple statistical methods are not suitable for complex nonlinear time-series problems of slope deformation, while machine learning is a promising alternative to traditional statistical methods due to its nonlinear learning capabilities.
The most applied traditional machine learning approaches for time-series forecasting are support vector regression (SVR) [15,16] and extreme gradient boosting (XGBoost). Compared with statistical approaches, machine learning approaches more easily capture the complexity, dynamics, and nonlinear characteristics of the nonlinear time-series displacement data of slopes. Currently, machine learning approaches have made some achievements in slope displacement prediction research [17,18]. For example, Miao et al. [19] proposed a coupled model of GA-SVR based on time series and a genetic algorithm and achieved an accurate landslide displacement prediction.
In recent years, deep learning has developed rapidly and performed better than traditional machine learning models in many studies [20]. The common deep learning approaches for time-series forecasting mainly include deep neural networks (DNNs) [21,22], convolutional neural networks (CNNs) [23], recurrent neural networks (RNNs) [24] and their variant models, as well as the transformer. Among them, RNNs have a strong model fitting ability for serialized data; and the most widely used RNN variants are the long-short-term memory (LSTM) [25] and gate recurrent unit (GRU) [26], which improve the shortcoming of the gradient explosion problem in conventional RNNs.
The LSTM and GRU are more effective on long sequence data than conventional RNNs and have achieved satisfactory prediction results in many engineering applications [27,28]. For example, Yang et al. [25] proposed an LSTM-based dynamic prediction model for landslide displacement and verified the performance of the model with observed data from two stepped landslides. Ma et al. [29] presented a novel deep learning method, T-GCN, which combined a GCN and a GRU to predict the slope deformation. However, when the input sequence is long, the performance of the LSTM model is still poor. LSTM models that include an attention mechanism are effective to improve the performance of the LSTM and have been applied to slope deformation in recent years [24,30,31].
The transformer is a new deep learning model proposed in 2017. The transformer model replaces the traditional sequential structure of RNNs with a self-attention mechanism, which allows the model to be trained in parallel and have access to global information. Currently, transformer models have been successfully applied to time-series forecasting problems [32,33]. However, few studies have applied the transformer to slope deformation prediction.
Currently, most studies for slope deformation focus on individual machine learning approaches compared to the ARIMA to show the nonlinear prediction performance of machine learning approaches. However, there are significant performance differences among different machine learning approaches, and few studies have conducted a comprehensive comparative analysis of the common machine learning models. In addition, some of the latest deep learning methods, such as the transformer, although outstanding in time series prediction problems in other fields, have hardly been studied for and applied to slope deformation prediction. Each model has different characteristics and application scenarios; thus, choosing an accurate model for slope deformation prediction is critical.
To address the above problem, in this paper, we conducted a comparative investigation of machine learning approaches for slope deformation prediction based on monitored time-series displacement data. The contributions in this paper can be summarized as follows:
(1)
We established two traditional machine learning models of the SVR and XGBoost and eight deep learning models, including the conventional DNN, CNN, RNN, LSTM, BiLSTM, Attention-LSTM, GRU, and transformer, while employing the linear ARIMA model as a baseline for comparison. The displacement data of seven in situ monitoring points of the Huanglianshu landslide located in eastern Chongqing were used as an example for verification.
(2)
We comparatively analyzed the prediction performance of the employed models based on four error metrics and provided a suggestion for selecting an appropriate prediction model for slope deformation.
The rest of this paper is organized as follows. Section 2 introduces the overview of the study area. Section 3 describes the applied prediction models in detail. Section 4 describes the predicted results of the models. Section 5 comparatively analyses the results. Section 6 concludes the paper.

2. Study Area

2.1. Overview of the Huanglianshu Landslide

The Huanglianshu landslide is located in Dabao village, Anping Town, Fengjie County, Eastern Chongqing, on the south bank of the Changjiang River in the upper reaches of the Qutang gorge. It is adjacent to the Yangtze River in the north and approximately 30 km upstream of the Qutang gorge. The geographical location is illustrated in Figure 1.
The Huanglianshu landslide is an old landslide that has been deformed many times in history. Local deformation of the landslide occurred after 836 mm rainfall on 30 May 2012. The landslide as a whole slid to form a sliding area with an area of approximately 6.25 × 10 4 m2 and a volume of approximately 80 × 10 4 m3. In the early morning of 31 May, a large-scale overall sliding occurred in the lower deformation area of the landslide, and the landslide deformation area disintegrated seriously. After the landslide occurred, its deformation continued.
The Huanglianshu landslide is a representative landslide in the Three Gorges reservoir area. Since the storage operation of the Three Gorges Project, the original equilibrium of the landslides in the reservoir area has been seriously damaged, and the Three Gorges reservoir area has become a landslide-disaster-prone area. This paper analyzes the Huanglianshu landslide as an example, which is of great practical significance for the early warning of such landslides.

2.2. Monitoring Data

Fourteen GPS displacement monitoring points (FJ001–FJ014) were established in the Huanglianshu landslide working area, and the monitoring equipment recorded landslide displacements from May 2012 to June 2013. Since the large-scale local sliding resulted in the destruction of several GPS monitoring devices, we selected the monitoring points FJ003, FJ004, FJ005, FJ006, FJ009, FJ010, and FJ011 with complete records for research to reveal the deformation characteristics of the slope after local deformation and establish an appropriate prediction model. The employed data in this paper was provided by the Center for Hydrogeology and Environmental Geology of the China Geological Survey. The geographical location relationship of the monitoring points studied in this paper is illustrated in Figure 2.
Figure 3 shows the displacement change curve of the monitoring points of the Huanglianshu landslide. There is a complicated nonlinear deformation characteristic of the Huanglianshu landslide. After the local large deformation of the Huanglianshu landslide occurred on 31 May 2012, the deformation characteristics of each monitoring point were different in the following years. Among the points, the monitored average displacements on FJ003, FJ004, FJ005, and FJ006 were small, while the monitored average displacements on FJ009, FJ010, and FJ011 were large.
The deformation curves in Figure 3 illustrate that the displacement changes of monitoring points were closely related to the relative positions of the monitoring points and the deformation area. As is shown in Figure 2, monitoring points FJ010 and FJ011 were located in the local deformation area; monitoring point FJ009 was located on the right side of the deformation area; and the monitoring points FJ003, FJ004, FJ005, and FJ006 were located outside the local deformation area. Within one year after the occurrence of the local deformation, the sliding velocity of the monitoring points within the deformation was approximately 10 times higher than that of monitoring points above the deformation area.

3. Machine Learning Approaches for Slope Deformation Prediction

3.1. Overview

The prediction process for slope deformation is shown in Figure 4. First, the displacement data at seven monitoring points of the Huanglianshu landslide are divided into a training set and a test set after pre-processing. Second, the training set and test data are divided into multiple datasets according to the time-series forecasting strategy. Third, various statistical and machine learning models are established for training datasets. Finally, the prediction accuracy of each model is verified on the test set data. In this paper, we established a statistical model and 10 machine learning models for slope deformation prediction. The structures of the employed models are illustrated in Figure 5.

3.2. Data Preprocessing

3.2.1. Linear Interpolation

The collected data at the Huanglianshu landslide include the daily recorded data of seven monitoring points from June 2012 to June 2013, with an original data size of 7 × 158. However, since the recording interval of each month was different, some data included more than one day. Since the recording conditions at each monitoring point were the same, we took the FJ003 monitoring point as an example, and the recording interval is shown in Figure 6. On this basis, we adopted the linear interpolation method to transform the original data into daily data from 7 June 2012 to 31 May 2013. The total amount of data after interpolation was 7 × 359.

3.2.2. Data Standardization

Before establishing the time-series prediction model, it is necessary to standardize the series data to eliminate the influence of dimensions. In this paper, the maximum and minimum normalization method was used to normalize the displacement data to [0, 1]. The calculation method is as shown in Equation (1).
x = x x m i n x m a x x m i n
where x represents the calculated value, x m a x is maximum value in the sequence data, and x m i n is minimum in the sequence data.

3.3. Time-Series Forecasting Strategy

The time-series problem can generally be expressed as follows. For a time-series containing N historical observation values [ y 1 , y 2 , , y N ], the recent d historical observation values [ y N d + 1 , , y N ] are used to predict the values [ y N + h , , y N + h ] at the next H time points. When predicting according to the direct strategy [34], first, the sliding window method [35] is used to obtain H datasets that contain input data and target data, i.e., the d historical data as the input and next h data as the output. Second, H independent regression models are trained according to the functional relationship shown in Equation (2). Finally, H historical datasets are input into each model to predict the value at the corresponding time.
y ^ N + h = f ^ h ( y N d + 1 , , y N ) , h 1 , , H
where [ y N d + 1 , , y N ] is the input data, y N + h is the target data, H is the length of the prediction time, and f ^ h is the regression model.
In this paper, first, the time-series displacement data at each monitoring point were divided into a training set and a test set. Second, each model was trained to obtain optimal parameters based on the training data. Third, the test data were predicted based on the trained model, as shown in Figure 4.
Specifically, in the training process, first, the training data were divided into an array of size 283 × 7 based on a sliding window with a window size of seven and a time step of one. Each array has a size of seven, where the first six data are used as inputs to the deep learning model, and the seventh data is used as the output. Second, the model learned the features of 283 displacement arrays and established a correspondence between the input data and output data. In the predicting process, first, the test data were divided into an array of size 66 × 7 based on the sliding window method. Second, the first six data in each array were fed into the trained model to predict the seventh displacement data. Finally, the predicted displacements were compared with the true displacements to evaluate the model performance.

3.4. Employed Statistical and Machine Learning Approaches

3.4.1. Statistical Approach: ARIMA

The most popular statistical model for time-series prediction is the ARIMA. The ARIMA is a traditional linear prediction model that has three undetermined parameters: autoregressive term order p, sequence difference order d, and moving average term order q. The process of determining p, d, and q is called the ARIMA model order determination.
In the entire order determination and optimation process of the ARIMA model, three tests need to be conducted: (1) the sequence stationarity test; (2) the residual white noise test; and (3) the parameter significance test. The workflow of the ARIMA model is shown in Figure 5a. In this paper, the Akaike information criterion (AIC) and the Bayesian information criterion (BIC) were used to automatically determine the optimal parameters (p, d, and q). The optimal parameters of the ARIMA model for the time-series displacement data of the Huanglianshu landslide are shown in Table 1.

3.4.2. Traditional Machine Learning Approaches

(1)
SVR
The SVR involves the application of a support vector machine (SVM) [36] to a regression problem. As is shown in Figure 5b, the principle of the SVR is to find a hyperplane so that the distance between the sample point farthest from the hyperplane and the plane is the shortest. In the nonlinear regression problem, the SVR can achieve a better regression effect by adding a kernel function. The choice of the kernel function is critical for performance. The main optional kernel functions include the linear function, polynomial function, and radial basis function (RBF).
(2)
XGBoost
As shown in Figure 5c, XGBoost is a tree-based model. In this paper, the sliding window representation was used to transform the time-series dataset of Huanglianshu landslide into supervised learning, and the k-fold cross validation method was exploited to evaluate the performance of the model. When the input is a time series and the output is a predicted value, each leaf node in the model is no longer a category but instead is a value, i.e., the prediction score. The calculation principle is shown in Equation (3).
y i ^ = k = 1 K f k ( x i )
where y i ^ is the predicted value, x i is the feature vector, f k ( x i ) is the calculated value for each tree, and K is the total number of trees; in this paper, K is 1000.

3.4.3. Deep Learning Approaches

(1)
Conventional DNN
The DNN is a neural network that contains multiple hidden layers. DNNs are divided into three types of neural network layers: the input layer, hidden layers, and the output layer. Generally, the first layer is the input layer, the middle layers are hidden layers, and the last layer is the output layer. In this paper, we established a conventional DNN model as a comparison baseline, and the layers of the model were fully connected to each other.
(2)
CNN
The CNN is type of feedforward neural network. Common CNNs consist of a convolution layer, a pooling layer, and a fully connected layer. CNNs play an important role in the field of computer vision; furthermore, CNNs can also be adapted to time-series problems. The difference between them is that a two-dimensional convolution kernel is applied to an image, while a one-dimensional convolution kernel is applied to a time series.
(3)
Conventional RNN, LSTM, GRU, BiLSTM and Attention-LSTM
The RNN is a recurrent neural network that is mainly used to model sequence data. Common RNNs consist of an input layer, hidden layers, and an output layer. In contrast to the DNN and CNN, the RNN is quite effective for data with sequential characteristics and possesses a strong model fitting ability, which can mine information between time-series data.
The LSTM and GRU are two variants that improve the limitation of RNNs for long-term dependence. The LSTM is a special network structure with three “gates” that can be used to remove or increase the information of the “cell state”. The gate used to increase information is called the “update gate”, the gate used for forgetting is called the “forget gate”, and the “output gate” is used for the output. The GRU is a variant of the LSTM model that combines the input gate and forget gate into an update gate. The internal structures of the RNN, LSTM, and GRU are shown in Figure 7.
BiLSTM [37] is a combination of the forward LSTM and backward LSTM. Modeling time series with an LSTM cannot encode back-to-back information, but a BiLSTM can better capture the relationships between sequence data in both directions.
The Attention-LSTM model is an LSTM model with an attention mechanism. In recent years, the attention mechanism, which can focus on the key part of time-series data to reduce the impact of nonkey data for prediction and is considered a fully connected layer and softmax function, has been widely used.
(4)
Transformer
The transformer is a sequence-to-sequence (seq2seq) structure [38] with encoder and decoder networks that use a self-attention mechanism. Compared with the traditional RNN and CNN, the transformer model proposes a new position encoding mechanism to capture the time-series information between input data; thus, the transformer model can fully capture the time-series information with richer dimensions. In the time-series prediction problem, the encoder takes the history of the time-series as the input, while the decoder predicts the future value via an autoregression. The employed transformer model in this paper included the encoder part of the transformer and finally connected a fully connected layer.

4. Results

4.1. Results of Predicted Displacements

In this paper, a total of eight deep learning models were established. Hyperparameters [39] are essential for optimizing the training performance of deep learning models. The optimal hyperparameters for the eight models were finally determined based on multiple training experiments, as are shown in Table 2.
The predicted displacements of the eleven models and the measured displacements on the test data at seven monitoring points of the Huanglianshu landslide are shown in Figure 8b. As is shown in Figure 8b, the eleven models generally showed that there was a lag in the predicted displacements, i.e., the predicted value at time t was usually close to the measured value at time t 1. The lag in the predicted displacements suggests the existence of autocorrelation in the time-series displacement data of the Huanglianshu landslide.
We conducted an autocorrelation function on the seven time-series displacement data to verify the autocorrelation of the landslide displacements. The autocorrelation function describes the correlation between the data at the current moment and the subsequent moments. As are shown in Figure 8a, the autocorrelation results indicate the varying degrees of autocorrelation in the seven sets of displacement data, which were consistent with the predicted results in Figure 8b.

4.2. Prediction Metrics and Accuracy

In this paper, we used four metrics to evaluate the prediction performance of the eleven models for slope deformation: the mean absolute error (MAE), mean absolute percentage error (MAPE) [40], root mean square error (RMSE), and R 2 , which are defined in Table 3. The results of the prediction accuracy when using the above four metrics are shown in Table 4. The smaller the MAE or the RMSE, or the MAPE value, the higher the prediction accuracy of the model, and MAPE = 0 indicates a perfect model; the closer the value of R 2 is to 1, the higher the prediction accuracy of the model.

5. Comparative Analysis

5.1. Stability Test of the Displacement Data and Evaluation of Model Performance

5.1.1. Stability Test of the Mointored Displacement Data

In order to investigate the deformation pattern of the time-series displacement data within and outside the deformation area of the Huanglianshu landslide after local sliding, we conducted a stability test on all monitoring points using the augmented Dickey Fuller (ADF) test [41]. The ADF test determines whether there is a unit root in the time series. If the time series is stationary, there is no unit root. Therefore, the ADF test assumes the existence of a unit root; and if the obtained test statistic p-value is less than significance levels 10, 5, or 1%, it is equivalent to having 90, 95, or 99% certainty against the original hypothesis, respectively. In this paper, we adopted 5% as the significance level, i.e., the situation of the p-value > 5% indicated that the time-series was nonstationary.
As are listed in Table 5, the displacement data at monitoring points FJ003, FJ004, FJ005, FJ006, and FJ009 were stationary, and the displacement data at monitoring points FJ010 and FJ011 within the deformation area were not stationary. The above results indicate that the displacement data within the deformation area had more complex nonlinear relationships.

5.1.2. Evaluation of Model Prediction Performance

According to Table 4, we evaluated the prediction performance of the eleven models when using the four different evaluation metrics in Figure 9. The evaluation results were as follows:
(1)
R 2 reflected the overall fit of the prediction models to the real displacement data. As shown in Figure 9b, when using R 2 as the evaluation index, all models had the highest prediction accuracies at the monitoring points FJ003 and FJ009, while the prediction accuracies at monitoring points FJ004, FJ005, FJ006, and FJ011 were relatively low. Specifically, monitoring point FJ006 had the lowest prediction accuracy. As a whole, except for the monitoring point FJ006, the R 2 of all other monitoring points reached more than 0.5, which showed that, after large-scale local sliding of the slope, the displacement data both within the deformation area and its surrounding area had a certain linear or nonlinear change characteristic, i.e., the slope deformation could be learned to achieve accurate prediction through a specific model.
(2)
Upon combing through the results of the four evaluation metrics in Figure 9, it can be concluded that the optimal prediction models at each monitoring point were the same when using the MAE and MAPE as evaluation metrics, and the optimal prediction models at each monitoring point were the same when using the R 2 and RMSE as evaluation metrics. In summary, the ARIMA, the Attention-LSTM, and the transformer model had higher prediction accuracies. In order to compare the accuracy differences of the models for the displacement data within and outside the deformation area, we list the best and worst prediction models at the seven monitoring points of the Huanglianshu landslide in Table 6. Table 6 indicated that, although the ARIMA model performed with outstanding prediction accuracy for the displacement data outside the deformation area, it did not perform as well as the transformer and Attention-LSTM models for the displacement data inside the deformation area. Second, the traditional machine learning models XGBoost and SVR had the worst prediction accuracies.

5.2. Comparative Analysis of the Model Characteristics Based on Prediction Results

Typically, the prediction performance of the model is mainly affected by the model characteristics and the original displacement data. The analysis of the influence of the characteristics of the models on the prediction results is as follows.
The ARIMA statistical method was effective for predicting short-term stationary time-series data. In this paper, we employed the displacement data of the next two months for prediction, which belonged to the range of short-term prediction; and the prediction results proved the ability of the ARIMA model to predict short-term problems. However, the ARIMA utilized a traditional linear method, and it is difficult to capture the nonlinear relationships between the displacement data of slopes.
As is shown in Table 6, the prediction accuracy of the ARIMA on the nonstationary displacement data at the monitoring points of FJ010 and FJ011 within the deformation area was much lower than that of most machine learning models, which suggests that ARIMA is better suited to predicting more stable displacement data at monitoring points outside the deformarion area. However, for the actual time-series displacement data of slopes, because of the various uncertain factors, such as weather changes and rainfall, the original data were normally nonstationary. Therefore, the ARIMA model was actually not the best choice for slope deformation prediction in most cases.
The prediction performance of the SVR model was closely related to the selection of the model kernel function. As mentioned in Section 3.4.2, we adopted three kernel functions to predict the displacement data at each monitoring point separately. The calculation results show that the displacement data at monitoring points FJ006, FJ010, and FJ011 were suitable for the RBF kernel function, and other monitoring points were suitable for the linear kernel function. As is shown in Table 3, the SVR model achieved more accurate prediction results for the displacement data of monitoring points FJ006, FJ010, and FJ011, which indicates that the SVR was better at handling nonstationary displacement data within the deformation area than the stationary displacement data outside the deformation area.
Moreover, the traditional machine learning models SVR and XGBoost have outlier prediction values. As shown in Figure 10, the SVR model had outlier prediction values that deviated seriously from the real situation at monitoring point FJ011 (Figure 10a), and the XGBoost model had outlier prediction values at monitoring point FJ003 (Figure 10b). Although the SVR model and XGBoost model can roughly predict the changing trend of displacement data, they have difficulty capturing unusual fluctuations in the data. In other words, the SVR and XGBoost models have more stable performance fitting the data in the normal range, while the ARIMA model is better than the SVR and XGBoost for outlier fitting. Therefore, when the displacement data have noise, the traditional machine models SVR and XGBoost would be more sensitive to overfitting.
For the deep learning models, overall, the conventional RNN and its variant LSTM, GRU, BiLSTM, and Attention-LSTM models achieved more accurate predictions compared with the conventional DNN model and CNN model, and the RNN variant model had better prediction performance than the conventional RNN at the monitoring points within and outside the deformation area. In a conventional RNN model and its variant models, the network structure can store the information on each neuron; whereas in a conventional DNN, information only passes through the network in the forward direction. As mentioned in Section 4.2, there was a first-order autocorrelation between the before and after data in the time-series displacement data of the Huanglianshu landslide; thus, the conventional RNN model and its variant models were more suitable for slope deformation prediction.
Among all deep learning models, the CNN model employed in this paper accomplished the lowest prediction accuracy at all monitoring points. Actually, the CNN method was placed before the input of the LSTM model for the first step to reduce the features [24,42] in most applications. The CNN method was used to capture the deep features in the data, and then the LSTM method was adopted to extract the time-series features in the data. Therefore, using a single CNN model is not suitable for predicting slope deformation.
Among the conventional RNN and its variant models, the BiLSTM and the Attention-LSTM models were more prominent, especially the Attention-LSTM model. The better prediction performance of the LSTM and GRU models was due to their improvement of the long-term dependence problem of the RNN models. As shown in Figure 10, when using the R 2 as the metric, the prediction accuracy of the Attention-LSTM model at five monitoring points of FJ004, FJ006, FJ009, FJ010, and FJ011 was much higher than that of the ordinary LSTM. As listed in Table 5, the values of the ADF test at FJ006, FJ010, and FJ011 were greater than 0.01, and the value at FJ004 was close to 0.01. Therefore, when the displacement data are more unsteady, the fitting performance of the Attention-LSTM model is more outstanding.
The above-predicted results at seven monitoring points of the Huanglianshu landslide indicate that the introduction of the attention mechanism has an important influence on the nonlinear displacement data. The deformation of a slope is influenced by many factors, such as rainfall and weather changes, and the influence weight of these potential factors on slope deformation is reflected in the time-series displacement data. The conventional LSTM model cannot focus on the important parts of the displacement data. However, applying the attention mechanism to the LSTM model allows the neural network to adaptively focus on the local displacement data of input time-series that are more important to the current output. Therefore, the LSTM model with an attention mechanism was more effective for slope deformation prediction.
As listed in Table 4, the values of the four prediction metrics of the transformer model on monitoring points FJ010 and FJ011 within the deformation area were significantly higher than other models, while the displacement data of monitoring points FJ003 and FJ005 outside the deformation area also had high prediction accuracy. This indicates that the transformer model has strong nonlinear processing ability, and the displacement within and outside the deformation area of the Huanglianshu landslide can be accurately predicted by the transformer model. The network structure of the transformer model is entirely composed of the self-attention mechanism. The self-attention mechanism is a variant of the attention mechanism, is less dependent on external information, and is better at capturing internal correlations between displacement data. The self-attention mechanism focuses on the important parts of the time-series data by assigning different weights, which reduces the degree of influence of abnormal points on the prediction results to improve the prediction accuracy.
In addition, for some special slopes, such as the rainfall-induced landslides [43,44], the landslide displacement is mainly influenced by the amount of rainfall; furthermore, the rainfall in various regions of the world shows obvious annual periodicity, which means that the displacement data of the slope have year-long dependency. Compared with the LSTM model, the transformer model had stronger long-term dependency modeling capabilities and has a better effect on long sequences. Therefore, the transformer model is better adapted to predict nonlinear landslide displacements affected by multiple factors.

5.3. Discussion

Based on the above analysis, we first conclude that the displacement data at the monitoring points within the deformation area of the Huanglianshu landslide have more complex nonlinear relationshisps than the displacement data at the monitoring points outside the deformation area. Second, the ARIMA statistical method is effective for predicting short-term stationary time-series displacement data within the deformation area, and the machine learning models are better at handling the nonstationary displacement data outside the deformation area. Third, the traditional machine learning models of the SVR and XGBoost were less effective than the ARIMA model at fitting outlier points in displacement data. Finally, deep learning approaches have better application prospects in slope deformation prediction problems, especially LSTM models with an attention mechanism, as well as transformer models.
At present, there is almost no research on applying the transformer model to slope deformation predictions. For time-series forecasting problems, the transformer model far surpassed the traditional long- and short-term memory networks in a variety of application scenarios due to its unique network structure. In recent years, with the rapid development of data science, deep learning has made outstanding achievements in various disciplines, which are not limited to the natural language processing (NLP) field. The application of deep learning to solve slope deformation prediction problems will become a powerful solution. In this paper, we established different prediction models based on the monitoring data of the Huanglianshu landslide and the development history of the time-series prediction method. We hope to apply the latest and most applicable computer technologies to slope prevention and early warning work, and to make advanced computer technologies an effective means of solving various geological disasters to compensate for the limitations of traditional physical means.
In the future, we will first comprehensively consider multiple influencing factors of slope deformation and exploit the existing deep learning approaches to establish a multivariable slope displacement prediction model. Second, more data interpolation methods will be added to the analysis. Finally, at present, spatiotemporal prediction [45,46] has been widely supported in traffic flow forecasting, that is, forecasting from two aspects, time and space, which is a great inspiration for multimonitoring landslides. Therefore, we will develop a spatiotemporal model for slope deformation prediction according to the development of the spatiotemporal prediction in the transportation field.

6. Conclusions

In this paper, to obtain an accurate slope deformation prediction, we conducted a comparative investigation between the ARIMA statistical method and 10 machine learning approaches based on monitored time-series displacement data of the Huanglianshu landslide located in eastern Chongqing. Four error metrics were used to comparatively analyze the prediction performance of all prediction models. The prediction results and comparative analysis of the model characteristics indicate the following: (1) the displacement data within the deformation area of the Huanglianshu landslide have more complex nonlinear relationships than the displacment data outside the deformation area; (2) machine learning approaches have a higher capacity to capture the nonlinear relationships between the time-series displacement data of slopes than the ARIMA model; and (3) among all the models, the Attention-LSTM model and the transformer model achieved the highest prediction accuracies for slope deformation, and the transformer model was better adapted to predict nonlinear landslide displacements affected by multiple factors. In the future, we will consider more factors that affect slope deformation and employ state-of-the-art deep learning approaches to address multivariable slope deformation prediction problems.

Author Contributions

Conceptualization, N.X. and G.M.; methodology, N.X., Q.Y., Y.S. and G.M.; validation, N.X., Q.Y., Y.S. and G.M.; formal analysis, N.X. and G.M.; investigation, N.X., Q.Y., Y.S. and G.M.; resources, Q.Y. and Y.S.; data curation, N.X., Q.Y. and Y.S.; writing—original draft preparation, N.X., Q.Y., Y.S. and G.M.; writing—review and editing, N.X. and G.M.; funding acquisition, G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of China (Grant No. 11602235 and 42277161).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the editor and the reviewers for their contributions.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ARIMAAutoregressive Integrated Moving Average
ADFAugmented Dickey Fuller
AICAkaike Information Criterion
BICBayesian Information Criterion
CNNConvolutional Neural Network
DNNDeep Neural Network
GRUGate Recurrent Unit
LSTMLong Short Term Memory
MAEMean Absolute Error
MAPEMean Absolute Percentage Error
MSEMean Square Error
NLPNatural Language Processing
RMSERoot Mean Square Error
RNNRecurrent Neural Network
RBFRadial Basis Function
SVMSupport Vector Machine
SVRSupport Vector Regression
XGBoosteXtreme Gradient Boost

References

  1. Intrieri, E.; Gigli, G.; Mugnai, F.; Fanti, R.; Casagli, N. Design and implementation of a landslide early warning system. Eng. Geol. 2012, 147–148, 124–136. [Google Scholar] [CrossRef] [Green Version]
  2. Dai, F.; Lee, C.; Ngai, Y. Landslide risk assessment and management: An overview. Eng. Geol. 2002, 64, 65–87. [Google Scholar] [CrossRef]
  3. Yao, W.; Zeng, Z.; Lian, C.; Tang, H. Training enhanced reservoir computing predictor for landslide displacement. Eng. Geol. 2015, 188, 101–109. [Google Scholar] [CrossRef]
  4. Li, Z.; Cheng, P.; Zheng, J. Prediction of time to slope failure based on a new model. Bull. Eng. Geol. Environ. 2021, 80, 5279–5291. [Google Scholar] [CrossRef]
  5. Feng, X.T.; Zhao, H.; Li, S. Modeling non-linear displacement time series of geo-materials using evolutionary support vector machines. Int. J. Rock Mech. Min. Sci. 2004, 41, 1087–1107. [Google Scholar] [CrossRef]
  6. De Gooijer, J.; Hyndman, R. 25 years of time series forecasting. Int. J. Forecast. 2006, 22, 443–473. [Google Scholar] [CrossRef] [Green Version]
  7. Deb, C.; Zhang, F.; Yang, J.; Lee, S.; Shah, K. A review on time series forecasting techniques for building energy consumption. Renew. Sustain. Energy Rev. 2017, 74, 902–924. [Google Scholar] [CrossRef]
  8. Wasowski, J.; Bovenga, F. Investigating landslides and unstable slopes with satellite Multi Temporal Interferometry: Current issues and future perspectives. Eng. Geol. 2014, 174, 103–138. [Google Scholar] [CrossRef]
  9. Gariano, S.; Guzzetti, F. Landslides in a changing climate. Earth-Sci. Rev. 2016, 162, 227–252. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, H.; Tian, H.Q.; Li, Y.F. Comparison of two new ARIMA-ANN and ARIMA-Kalman hybrid methods for wind speed prediction. Appl. Energy 2012, 98, 415–424. [Google Scholar] [CrossRef]
  11. Reichenbach, P.; Rossi, M.; Malamud, B.; Mihir, M.; Guzzetti, F. A review of statistically-based landslide susceptibility models. Earth-Sci. Rev. 2018, 180, 60–91. [Google Scholar] [CrossRef]
  12. Chen, W.; Shahabi, H.; Shirzadi, A.; Hong, H.; Akgun, A.; Tian, Y.; Liu, J.; Zhu, A.X.; Li, S. Novel hybrid artificial intelligence approach of bivariate statistical-methods-based kernel logistic regression classifier for landslide susceptibility modeling. Bull. Eng. Geol. Environ. 2019, 78, 4397–4419. [Google Scholar] [CrossRef]
  13. Li, H.; Xu, Q.; He, Y.; Deng, J. Prediction of landslide displacement with an ensemble-based extreme learning machine and copula models. Landslides 2018, 15, 2047–2059. [Google Scholar] [CrossRef]
  14. Shi, J.; Guo, J.; Zheng, S. Evaluation of hybrid forecasting approaches for wind speed and power generation time series. Renew. Sustain. Energy Rev. 2012, 16, 3471–3480. [Google Scholar] [CrossRef]
  15. Huang, Y.; Zhao, L. Review on landslide susceptibility mapping using support vector machines. Catena 2018, 165, 520–529. [Google Scholar] [CrossRef]
  16. Pradhan, B. A comparative study on the predictive ability of the decision tree, support vector machine and neuro-fuzzy models in landslide susceptibility mapping using GIS. Comput. Geosci. 2013, 51, 350–365. [Google Scholar] [CrossRef]
  17. Ma, Z.; Mei, G.; Piccialli, F. Machine learning for landslides prevention: A survey. Neural Comput. Appl. 2021, 33, 10881–10907. [Google Scholar] [CrossRef]
  18. Deng, L.; Smith, A.; Dixon, N.; Yuan, H. Machine learning prediction of landslide deformation behaviour using acoustic emission and rainfall measurements. Eng. Geol. 2021, 293. [Google Scholar] [CrossRef]
  19. Miao, F.; Wu, Y.; Xie, Y.; Li, Y. Prediction of landslide displacement with step-like behavior based on multialgorithm optimization and a support vector regression model. Landslides 2018, 15, 475–488. [Google Scholar] [CrossRef]
  20. Ma, Z.; Mei, G. Deep learning for geological hazards analysis: Data, models, applications, and opportunities. Earth-Sci. Rev. 2021, 223. [Google Scholar] [CrossRef]
  21. Xiao, L.; Zhang, Y.; Peng, G. Landslide susceptibility assessment using integrated deep learning algorithm along the china-nepal highway. Sensors 2018, 18, 4436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lian, C.; Zeng, Z.; Yao, W.; Tang, H. Multiple neural networks switched prediction for landslide displacement. Eng. Geol. 2015, 186, 91–99. [Google Scholar] [CrossRef]
  23. Pei, H.; Meng, F.; Zhu, H. Landslide displacement prediction based on a novel hybrid model and convolutional neural network considering time-varying factors. Bull. Eng. Geol. Environ. 2021, 80, 7403–7422. [Google Scholar] [CrossRef]
  24. Wang, H.; Zhang, L.; Luo, H.; He, J.; Cheung, R. AI-powered landslide susceptibility assessment in Hong Kong. Eng. Geol. 2021, 288. [Google Scholar] [CrossRef]
  25. Yang, B.; Yin, K.; Lacasse, S.; Liu, Z. Time series analysis and long short-term memory neural network to predict landslide displacement. Landslides 2019, 16, 677–694. [Google Scholar] [CrossRef]
  26. Zhang, Y.G.; Tang, J.; He, Z.Y.; Tan, J.; Li, C. A novel displacement prediction method using gated recurrent unit model with time series analysis in the Erdaohe landslide. Nat. Hazards 2021, 105, 783–813. [Google Scholar] [CrossRef]
  27. Tien Bui, D.; Tuan, T.; Klempe, H.; Pradhan, B.; Revhaug, I. Spatial prediction models for shallow landslide hazards: A comparative assessment of the efficacy of support vector machines, artificial neural networks, kernel logistic regression, and logistic model tree. Landslides 2016, 13, 361–378. [Google Scholar] [CrossRef]
  28. Liu, Z.; Gilbert, G.; Cepeda, J.; Lysdahl, A.; Piciullo, L.; Hefre, H.; Lacasse, S. Modelling of shallow landslides with machine learning algorithms. Geosci. Front. 2021, 12, 385–393. [Google Scholar] [CrossRef]
  29. Ma, Z.; Mei, G.; Prezioso, E.; Zhang, Z.; Xu, N. A deep learning approach using graph convolutional networks for slope deformation prediction based on time-series displacement data. Neural Comput. Appl. 2021. [Google Scholar] [CrossRef]
  30. Yan, L.; Chen, C.; Hang, T.; Hu, Y. A stream prediction model based on attention-LSTM. Earth Sci. Inform. 2021, 14, 723–733. [Google Scholar] [CrossRef]
  31. Cornia, M.; Baraldi, L.; Serra, G.; Cucchiara, R. Predicting human eye fixations via an LSTM-Based saliency attentive model. IEEE Trans. Image Process. 2018, 27, 5142–5154. [Google Scholar] [CrossRef] [Green Version]
  32. Schwaller, P.; Petraglia, R.; Zullo, V.; Nair, V.; Haeuselmann, R.; Pisoni, R.; Bekas, C.; Iuliano, A.; Laino, T. Predicting retrosynthetic pathways using transformer-based models and a hyper-graph exploration strategy. Chem. Sci. 2020, 11, 3316–3325. [Google Scholar] [CrossRef] [Green Version]
  33. Schwaller, P.; Laino, T.; Gaudin, T.; Bolgar, P.; Hunter, C.; Bekas, C.; Lee, A. Molecular Transformer: A Model for Uncertainty-Calibrated Chemical Reaction Prediction. ACS Cent. Sci. 2019, 5, 1572–1583. [Google Scholar] [CrossRef] [Green Version]
  34. Sorjamaa, A.; Hao, J.; Reyhani, N.; Ji, Y.; Lendasse, A. Methodology for long-term prediction of time series. Neurocomputing 2007, 70, 2861–2869. [Google Scholar] [CrossRef] [Green Version]
  35. Hindriks, R.; Adhikari, M.; Murayama, Y.; Ganzetti, M.; Mantini, D.; Logothetis, N.; Deco, G. Can sliding-window correlations reveal dynamic functional connectivity in resting-state fMRI? NeuroImage 2016, 127, 242–256. [Google Scholar] [CrossRef] [Green Version]
  36. Cai, Z.; Xu, W.; Meng, Y.; Shi, C.; Wang, R. Prediction of landslide displacement based on GA-LSSVM with multiple factors. Bull. Eng. Geol. Environ. 2016, 75, 637–646. [Google Scholar] [CrossRef]
  37. Zhang, K.; Zhang, K.; Cai, C.; Liu, W.; Xie, J. Displacement prediction of step-like landslides based on feature optimization and VMD-Bi-LSTM: A case study of the Bazimen and Baishuihe landslides in the Three Gorges, China. Bull. Eng. Geol. Environ. 2021. [Google Scholar] [CrossRef]
  38. Guo, Y.; Wu, T.; Liu, Y.; Jiang, Y. Univariate Forecasting of time series data Based on Seq2Seq Model. In Proceedings of the EITCE 2020: 2020 4th International Conference on Electronic Information Technology and Computer Engineering, Xiamen, China, 6–8 November 2020; pp. 642–646. [Google Scholar] [CrossRef]
  39. Dong, X.; Shen, J.; Wang, W.; Shao, L.; Ling, H.; Porikli, F. Dynamical hyperparameter optimization via deep reinforcement learning in tracking. IEEE Trans. Pattern Anal. Mach. Intell. 2021, 43, 1515–1529. [Google Scholar] [CrossRef]
  40. de Myttenaere, A.; Golden, B.; Le Grand, B.; Rossi, F. Mean Absolute Percentage Error for regression models. Neurocomputing 2016, 192, 38–48. [Google Scholar] [CrossRef] [Green Version]
  41. Caspi, I. Rtadf: Testing for bubbles with EViews. J. Stat. Softw. 2017, 81, 1–16. [Google Scholar] [CrossRef] [Green Version]
  42. Xie, H.; Zhang, L.; Lim, C. Evolving CNN-LSTM Models for Time Series Prediction Using Enhanced Grey Wolf Optimizer. IEEE Access 2020, 8, 161519–161541. [Google Scholar] [CrossRef]
  43. Liu, Y.; Liu, D.; Qin, Z.; Liu, F.; Liu, L. Rainfall data feature extraction and its verification in displacement prediction of Baishuihe landslide in China. Bull. Eng. Geol. Environ. 2016, 75, 897–907. [Google Scholar] [CrossRef]
  44. Zhu, X.; Xu, Q.; Tang, M.; Nie, W.; Ma, S.; Xu, Z. Comparison of two optimized machine learning models for predicting displacement of rainfall-induced landslide: A case study in Sichuan Province, China. Eng. Geol. 2017, 218, 213–222. [Google Scholar] [CrossRef]
  45. Tien Bui, D.; Tuan, T.; Hoang, N.D.; Thanh, N.; Nguyen, D.; Van Liem, N.; Pradhan, B. Spatial prediction of rainfall-induced landslides for the Lao Cai area (Vietnam) using a hybrid intelligent approach of least squares support vector machines inference model and artificial bee colony optimization. Landslides 2017, 14, 447–458. [Google Scholar] [CrossRef]
  46. Rahimi, M.; Shafieezadeh, A.; Wood, D.; Kubatko, E.; Dormady, N. Bayesian calibration of multi-response systems via multivariate Kriging: Methodology and geological and geotechnical case studies. Eng. Geol. 2019, 260. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the Huanglianshu landslide.
Figure 1. Geographical location of the Huanglianshu landslide.
Applsci 13 04677 g001
Figure 2. Schematic diagram of monitoring arrangement of the Huanglianshu landslide.
Figure 2. Schematic diagram of monitoring arrangement of the Huanglianshu landslide.
Applsci 13 04677 g002
Figure 3. Displacement change curve of monitoring points of the Huanglianshu landslide.
Figure 3. Displacement change curve of monitoring points of the Huanglianshu landslide.
Applsci 13 04677 g003
Figure 4. Workflow of slope deformation prediction.
Figure 4. Workflow of slope deformation prediction.
Applsci 13 04677 g004
Figure 5. Employed statistical and machine learning models. (a) Workflow of ARIMA model, (b) Workflow of SVR, (c) Workflow of XGBoost, (d) Workflow of Conventional DNN, (e) Workflow of CNN, (f) Workflow of RNN/LSTM/GRU, (g) Workflow of BiLSTM, (h) Worflow of Attention-LSTM, (i) Workflow of Transformer.
Figure 5. Employed statistical and machine learning models. (a) Workflow of ARIMA model, (b) Workflow of SVR, (c) Workflow of XGBoost, (d) Workflow of Conventional DNN, (e) Workflow of CNN, (f) Workflow of RNN/LSTM/GRU, (g) Workflow of BiLSTM, (h) Worflow of Attention-LSTM, (i) Workflow of Transformer.
Applsci 13 04677 g005
Figure 6. Data recording interval and linear interpolation representation of the displacement data of the Huanglianshu landslide.
Figure 6. Data recording interval and linear interpolation representation of the displacement data of the Huanglianshu landslide.
Applsci 13 04677 g006
Figure 7. Internal structures of the conventional RNN, LSTM, and GRU. (a) Conventional RNN, (b) LSTM, (c) GRU.
Figure 7. Internal structures of the conventional RNN, LSTM, and GRU. (a) Conventional RNN, (b) LSTM, (c) GRU.
Applsci 13 04677 g007
Figure 8. Predicted displacements and the autocorrelation functions of the seven monitoring points. (a) The autocorrelation functions; (b) the predicted and measured displacements.
Figure 8. Predicted displacements and the autocorrelation functions of the seven monitoring points. (a) The autocorrelation functions; (b) the predicted and measured displacements.
Applsci 13 04677 g008
Figure 9. Prediction performance of the models when using different evaluation metrics. (a) MAE, (b) R 2 , (c) MAPE, (d) RMSE.
Figure 9. Prediction performance of the models when using different evaluation metrics. (a) MAE, (b) R 2 , (c) MAPE, (d) RMSE.
Applsci 13 04677 g009
Figure 10. Outlier prediction values of SVR and XGBoost model. (a) Outlier prediction values of SVR. (b) Outlier prediction values of XGBoost.
Figure 10. Outlier prediction values of SVR and XGBoost model. (a) Outlier prediction values of SVR. (b) Outlier prediction values of XGBoost.
Applsci 13 04677 g010
Table 1. The optimal parameters of the ARIMA model for the displacement data of the Huanglianshu landslide.
Table 1. The optimal parameters of the ARIMA model for the displacement data of the Huanglianshu landslide.
ParameterFJ003FJ004FJ005FJ006FJ009FJ010FJ011
p1311011
d0101122
q0021000
Table 2. Hyperparameters of the deep learning models.
Table 2. Hyperparameters of the deep learning models.
Learning RateBatch SizeNumber of EpochsOptimizerLoss Function
0.0021001000AdamMSE
Table 3. Metrics of the prediction performance.
Table 3. Metrics of the prediction performance.
MetricDefinitionDescription
Mean absolute errorMAE = 1 n t = 1 n | y t y t | y t is the predicted value, y t is the measured
value, and n is the number of test data.
MAE represents the average absolute error
between measured and predicted values.
Mean absolute percent errorMAPE = 100 % n t = 1 n | y t y t | y t MAPE represents the percentage of errors
between measured and predicted values of
all samples.
Root mean square errorRMSE = 1 n t = 1 n ( y t y t ) 2 RMSE represents the root mean square error
between measured and predicted values.
Coefficient of determination R 2 = 1 t = 1 n ( y t y t ) 2 t = 1 n ( y ¯ y t ) 2 The closer the value of R 2 is to 1, the better
the prediction performance of the model.
Table 4. Prediction accuracies of all models.
Table 4. Prediction accuracies of all models.
MetricModelFJ003FJ004FJ005FJ006FJ009FJ010FJ011
MAEARIMA2.75142.69871.71342.07784.51724.21542.5601
SVR3.11673.05271.96692.02364.52564.57683.2575
XGBoost3.21072.76151.95412.22265.66544.24742.4652
DNN3.42442.94131.82092.19795.64164.44443.0304
CNN2.95973.07871.76002.23046.15004.93732.6765
RNN2.80052.85381.77162.14414.74174.63312.6581
LSTM2.72783.03491.77502.05825.37514.68903.0843
BiLSTM2.75483.08431.82522.07335.14704.51052.5058
Attention-LSTM2.82363.03411.76012.00815.27494.62902.3208
GRU2.69422.79171.79282.05494.89124.28502.5024
Transformer2.41173.00251.71172.02225.96673.99412.6732
MAPEARIMA0.17610.39370.13100.45460.01520.01130.0229
SVR0.19010.48820.15710.46320.01530.01220.0289
XGBoost0.20180.43640.14330.49970.01880.01140.0220
DNN0.22020.53300.13970.49140.01940.01200.0271
CNN0.18790.52940.13980.50620.02110.01330.0241
RNN0.17310.44480.13520.50380.01590.01230.0236
LSTM0.16740.52020.13400.46490.01800.01240.0275
BiLSTM0.16890.54080.13900.47330.01720.01200.0224
Attention-LSTM0.17350.52180.13050.47790.01810.01250.0208
GRU0.16980.43220.13750.47390.01650.01140.0223
Transformer0.16140.49460.13670.44440.02030.01060.0238
RMSEARIMA4.26553.88682.59592.93877.23596.90434.4374
SVR4.65444.05962.73872.82947.41127.37955.7565
XGBoost5.50684.10253.02213.23359.04566.88724.2473
DNN4.85584.05412.59113.02238.17836.44104.8310
CNN4.65764.14452.60573.00688.54566.73134.5103
RNN4.39233.98462.57072.96047.42216.85914.4002
LSTM4.31404.03512.52162.92648.24477.05404.5134
BiLSTM4.34293.95162.57712.92417.87106.66524.3341
Attention-LSTM4.37213.95832.54332.75597.48046.55114.2251
GRU4.37243.99432.58242.85767.61676.53044.2470
Transformer3.96803.94502.59142.82478.46816.31294.3167
R 2 ARIMA0.84230.50880.51700.32530.87590.56720.4649
SVR0.81220.46420.46250.37460.86980.50560.0995
XGBoost0.73710.45280.34540.18320.80610.56930.5098
DNN0.79560.46560.51880.28640.84150.62330.4057
CNN0.81190.44150.51340.29370.82690.58860.4472
RNN0.83270.48380.52640.31530.86950.57280.4738
LSTM0.83860.47060.54430.33100.83890.54820.4464
BiLSTM0.83650.49230.52400.33200.85320.59670.4895
Attention-LSTM0.83430.49060.53640.40670.86740.61030.5149
GRU0.83420.48130.52200.36210.86250.61280.5036
Transformer0.86350.49400.51870.37680.83010.63820.4936
Bold of numbers denotes the best performance value for each data set.
Table 5. Results of the ADF test for time-series displacement data at 7 monitoring points.
Table 5. Results of the ADF test for time-series displacement data at 7 monitoring points.
Monitoring PointFJ003FJ004FJ005FJ006FJ009FJ010FJ011
p-value0.0000710.0099870.0064900.0133650.0024630.0740620.797659
Table 6. The best and worst prediction models for displacement data within and outside the local deformation area.
Table 6. The best and worst prediction models for displacement data within and outside the local deformation area.
Monitoring PointOutside the Deformation AreaInside the Deformation Area
FJ003FJ004FJ005FJ006FJ009FJ010FJ011
Best modelTransformerARIMAAttention-LSTMAttention-LSTMARIMATransformerAttention-LSTM
Worst modelXGBoostCNNXGBoostXGBoostXGBoostSVRSVR
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xi, N.; Yang, Q.; Sun, Y.; Mei, G. Machine Learning Approaches for Slope Deformation Prediction Based on Monitored Time-Series Displacement Data: A Comparative Investigation. Appl. Sci. 2023, 13, 4677. https://doi.org/10.3390/app13084677

AMA Style

Xi N, Yang Q, Sun Y, Mei G. Machine Learning Approaches for Slope Deformation Prediction Based on Monitored Time-Series Displacement Data: A Comparative Investigation. Applied Sciences. 2023; 13(8):4677. https://doi.org/10.3390/app13084677

Chicago/Turabian Style

Xi, Ning, Qiang Yang, Yingjie Sun, and Gang Mei. 2023. "Machine Learning Approaches for Slope Deformation Prediction Based on Monitored Time-Series Displacement Data: A Comparative Investigation" Applied Sciences 13, no. 8: 4677. https://doi.org/10.3390/app13084677

APA Style

Xi, N., Yang, Q., Sun, Y., & Mei, G. (2023). Machine Learning Approaches for Slope Deformation Prediction Based on Monitored Time-Series Displacement Data: A Comparative Investigation. Applied Sciences, 13(8), 4677. https://doi.org/10.3390/app13084677

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop