1. Introduction
In recent years, the construction volume of large underground projects in China has grown [
1]. During underground engineering excavation, a large amount of rubble waste material is temporarily piled up near the project area to form a large-scale waste rockfill body. With the development and utilization of rubble waste material, the temporary waste material dump needs to be excavated and transported. Excavation disturbance inevitably causes slope failure. Therefore, the safety of waste rockfill slopes during excavation and transportation has become an important factor [
2].
Slope deformation monitoring and landslide warning rank high among the most effective safety prevention and control measures [
3]. Displacement monitoring forms the basis of landslide deformation monitoring. Slope displacement and its related parameters are important reference indicators for landslide early warning research [
4]. In 1968, Saito proposed the “three-segment” change rule of slope displacement–time curves [
5]. Since then, scholars have conducted a vast amount of research on slope deformation monitoring and landslide warning from different perspectives and with different means. For instance, Hoek (1970) analyzed the displacement–time curve measured at the Chuquicanata mine and proposed an extension method to extrapolate the displacement–time curves and to predict the time of landslides according to the curve trend [
6,
7,
8]. Zhong proposed the Verhulst inverse function model according to the characteristics of displacement–time curves [
9]. Qiang et al. analyzed and studied many typical landslide displacement monitoring results and proposed a landslide early warning method using acceleration as the evaluation index [
10]. In addition, many scholars have tried to combine various mathematical models with slope displacement–time curves. Such examples include the golden section method for early warning of the time of landslides [
11,
12,
13], the markov early warning model [
14,
15], the fuzzy mathematics early warning model [
16,
17], the poisson cycle early warning model [
18,
19,
20,
21], and the gradient-sinusoidal early warning model [
20,
22,
23].
For natural landslides, the displacement vectors at various points on the landslide body can generally be kept consistent. Even under special circumstances, the deformation of each point in the initial stage of slope deformation is disordered, but, once entering deformation acceleration, the displacement of each point on the landslide body and the time curve also tend to be the same. Therefore, the displacement data of one or more monitoring points can be used as the basis for slope deformation analysis and landslide warning research [
24]. However, the waste rockfill are massive and unconsolidated and feature dynamic excavation and transportation processes. In slope deformation monitoring, it is difficult to apply contact monitoring methods such as embedded sensors and GPS. Because of the granular structure characteristics of a rockfill body, the rubble particles on the slope will rotate during slope excavation, which leads to inconsistent deformation among the rubble. It is irrational to describe the displacement and warn of the deformation of the slope by a single point. Therefore, traditional slope monitoring and landslide warning methods are not applicable for large abandoned rockfill slopes.
To effectively solve the problem of catastrophic events that may occur on the slopes of the rockfills, it is necessary to select appropriate monitoring methods, to describe the deformation process of slopes with scientific quantitative indicators, and to provide early warning for slope landslides. GB-InSAR (ground-based interferometric synthetic aperture radar) can provide full-field deformation information within the monitored area [
25], which makes it possible to use the change of deformation region on the slope as the characteristic value to study the slope deformation. Therefore, the GB-InSAR system (IBIS-L system from Innovative Interferometric Radar for Environmental and Civil Engineering Applications, Italy) was used to carry out on-site monitoring tests on the slope of China’s first temporary large-scale groundwater-sealed oil storage dumping site in Huangdao, Shandong Province. Through data processing of the full-field displacement results obtained from monitoring (data cleaning, data interpolation, generating equal displacement lines, etc.), we propose to describe the deformation process of the slope by using the curve of the cross-sectional area of the constant displacement body with time. By summarizing the results of the data, a landslide early warning method with an equal displacement cross-sectional area as an early warning index is proposed to achieve the real-time monitoring and early warning of rockfill landslides during construction.
2. Project Overview
The waste rockfill is located about 500 m southwest of the underground project entrance, and the floor area is approximately 9.6 × 10
4 m
2. The ground where the waste rockfill body is located is high in the north and low in the south, and the difference of elevation between the south and north is about 20 m. The east and west directions are relatively gentle, and there is a slightly protruding depression in the middle part. The difference of elevation is about 4 m. The waste rockfill is dumped by a dump truck [
26], and it is leveled by a bulldozer. The rockfill body slopes can be divided into two stages based on the formation process. The boundary between the two stages are shown in
Figure 1 (the yellow dotted line).
The waste is mainly rubble produced by the blasting and crushing of granite gneiss.
Table 1 shows the petrophysical and mechanical properties of granite gneiss in intact rock mass. The rockfill is highly porous and unconsolidated with high permeability of the rubble. The bottom of the rockfill body is composed of large rubbles with an average particle size of 1 m, and the rubble particles are independent of each other. The slope within the monitoring area (the red polygon area in
Figure 1) is the residual slope after the occurrence of the previous landslides, and many large boulders are already broken, having an average size of 10–50 cm (
Figure 2). Besides, due to the occurrence of landslide, the rock rolling and impact will be broken, resulting in much coarse-grained soil left on the slope.
Due to the stacking height limits and the economic value of rubble, rubble waste is transferred to a dumping site about 1 km south of the site of rubble transfer, as shown in
Figure 3.
To prevent uncontrollable large-scale slippage caused by an excavation-related disturbance, rubble excavation is carried out in two phases. In the first phase, a “C”-shaped platform is constructed above the platform of the original rockfill body. Excavation is carried out using construction equipment located on the platform (
Figure 4) to control the location and slope height of the area affected by the excavation, thereby controlling the scale of the landslide and the direction of the sliding body. After the removal of most of the rubble via the above the platform, the second phase of the operation is to transport the rubble under the platform.
4. Early Warning Method for Landslides Based on the Rate of Change of the Cross-Sectional Area of Equal Displacement Bodies
The uniqueness of waste rockfill bodies makes it unsuitable to use single-point displacement as an indicator for slope deformation analysis or empirical judgment. Therefore, in this study, the deformation characteristics of waste rockfill slopes were characterized by the change of the cross-sectional area of the equal displacement body. On this basis, a warning method for waste rockfill slope landslides during excavation is proposed.
4.1. Ideas on Early Warning
(1A) Based on the results in [
2] and the 26 landslides monitored in the previous period (the remaining five landslides were used to verify the early warning method), the surface characteristics of waste rockfill slopes can be identified, and a preliminary judgment can be made about the possible landslide patterns in the excavation process. During the accumulation of waste rockfill, the density of some accumulation layers and areas caused by rainfall, the rolling of construction vehicles, etc. is quite different, or there might be a certain amount of abnormally large-sized rubble on the slope surface and in the shallow layer. Both conditions play a great role in supporting the rubble particles above the dense zone. Therefore, combined with the surface characteristics of the rockfill slope, a possible landslide pattern can be determined through the analysis of the characteristics of the cross-sectional area of different displacement bodies over time.
(1B) For possible landslides being monitored, based on the real-time constant displacement cross-sectional area data obtained from slope deformation monitoring, a time series prediction model such as ARIMA-SVR [
28] (Autoregressive Integrated Moving Average model-Support Vector Regression: a hybrid ARIMA and SVR model of time series prediction mode) can be used. Based on the existing data and the prediction of future data, the future trend of constant displacement cross-sectional areas can be predicted.
(2) Real-time monitoring data can be extracted, the real-time changing process of the cross-sectional area of equal displacement can be studied and analyzed, and the mode of slope deformation can be determined according to the change trend of the cross-sectional area of the equal displacement body.
(3) According to the measured (1A) and predicted values of the cross-sectional area of the equal displacement body (1B), the relationship between the cross-sectional area change of the equal displacement body in different modes and landslides of the rockfill slope can be analyzed, and an early warning judgment of such landslides can be determined. Then, the change process of the cross-sectional area of different landslides and equal displacement bodies can be independently analyzed.
(4) On the basis of the determined early warning criteria of the rockfill slope, an early warning grade division standard of the rockfill slope can be established.
4.2. Early Warning Indicators and Grading Standard
Considering the needs of on-site early warning in terms of the number of data calculations and the data extraction speed, under the condition of ensuring the accuracy of radar displacement monitoring, the deformation of −0.5 mm on the slope is easier to identify in time. Therefore, the cross-sectional area of the −0.5 mm equal displacement body was used as an index to validate the early warning method of the rockfill slope during excavation and transportation [
2]. The on-site conditions of waste rockfill slope landslides are complicated, and both different and similar types of landslides often occur in the same area. Therefore, the cross-sectional area of the −0.5 mm equal displacement body was normalized, and then the rate of change was used as a warning indicator of landslides. The cross-sectional area of the −0.5 mm constant displacement body at each moment after normalization is
where
S−0.5(
t) is the measured or predicted value of the cross-sectional area of the −0.5 mm constant displacement body at time
t and S
max is the cumulative value of the cross-sectional area of the −0.5 mm equal displacement body obtained from the last measurement before the slope starts sliding.
The normalized −0.5 mm equal displacement body cross-sectional area change rate is
Based on the field measurement results and the analogy between landslide cases, the characteristics of the cross-sectional area of the −0.5 mm equal displacement body of the two types of 29 landslides that were previously monitored in the rockfill were statistically analyzed. Then, early warning standards were developed according to two different types of landslides. Since the data selected were actually measured landslide data, the lower limit of the 95% confidence band of the mean value of the cross-sectional area of the two types of landslides of the −0.5 mm equal displacement body was selected as the standard to formulate the landslide warning level. According to the characteristic values of 50%, 30%, and 10% of the lower limit of the 95% confidence band of the average displacement rate of the cross-sectional area of the −0.5 mm equal displacement body before sliding, the landslide warning was divided into five levels, as shown in
Figure 6. Note that the rate of change of the cross-sectional area of the “B1”-type landslide of the −0.5 mm equal displacement body shows a decreasing trend within a certain period of the initial stage of the development of the displacement field. Therefore, the maximum value of
before the start of deformation acceleration is defined as the upper limit of the yellow warning for class B landslides. To make the early warning method more convenient and to ensure the compatibility of the two types of landslides, the early warning standard was uniformly divided into four levels based on the statistical results of the early warning thresholds of the two types of landslides, as shown in
Figure 11. The grading results and construction regulations were then established, as shown in
Table 4.
.
5. Method Validation
To check the feasibility and reliability of the landslide warning method of rockfill slopes, the last five landslides of two different types of landslides recorded during field monitoring were used to verify the landslide warning method.
The waste rockfill slope was taken as the remaining broken surface after the previous landslide, and the slope angle was approximately 35–38°. The excavation and loading speed during construction was approximately a 10-t load muck truck loaded every 6–9 min. During the five landslide monitoring periods, there was no rainfall, the daytime temperature was about 20–25 °C, and the nighttime temperature was about 15–20 °C.
The landslides occurred at 4:55 on 5 October 2013, 2:45 on 10 October 2013, 3:00 on 15 October 2013, 21:00 on 15 October 2013, and 9:38 on 16 October 2013. The final landslide areas were 450, 1000, 1100, 1600, and 1100 m
2, respectively. The specific records are shown in
Table 5.
First, the slope was initially classified according to whether there was a dense belt formed by the rolling of the slope or whether there was abnormally large-sized rubble in the middle of the slope. Landslides 44 and 53 conformed to the B1-type failure mode, while Landslides 67, 71, and 72 conformed to the B2-type failure mode. Then, ARIMA-SVR model was used to predict the variation trend of −0.5 mm equal displacement cross-section area of these five landslides, and the results are shown in
Figure 12. The average deviations between the predicted value and the measured value of five test landslides are shown in
Table 6.
According to Equation (2), the prediction results of five −0.5 mm equal displacement bodies of five landslides were normalized. Then, the deformation rates of the cross-sectional area of the −0.5 mm equal displacement bodies were calculated after normalization using Equation (3). The results are shown in
Figure 13.
From the analysis results of the above examples, the landslide early warning method and the early warning classification established in this paper could provide a warning to avoid risks before the occurrence of all five tested landslides.
Compared with the traditional landslide warning method based on slope deformation, it mainly depends on the displacement of some points on the slope. For a loose rockfill slope, it is feasible to express the deformation of the slope by the displacement of a certain point, but it is difficult to find a suitable point. However, as described in
Section 3.3, the variation curve of a −0.5 mm equal displacement cross-sectional area over time conforms to the form of a “three-section” curve. Combined with previous research experience, this will make it easier for the equal displacement cross-sectional area to be used as an early warning indicator.
It should be noted that this verification was based on the last instance of landslide of each type on the site, and the actual working conditions were similar to the previous landslides. The establishment of the landslide warning standard also depends on the statistical results of these previous similar landslide cases. Therefore, these four successful landslide warnings are limited by certain working conditions, such as rock lithology, rubble particle gradation, rolling of vehicles during the accumulation of rock mass, rainfall, and so on. These conditions will directly affect the repose angle of the rockfill, and the repose angle is an important factor affecting rockfill landslides. Generally, when the angle of the slope is less than the repose angle, the slope is unstable. Rainfall has a direct impact on the water content of rockfill, which is directly proportional to the angle of internal friction. For loose rockfill, due to the structural characteristics of the large pores, leakage and evaporation are rapid, so water content is usually small. Particle gradation will also directly affect the repose angle of the rockfill body. Small-particle-size rubble plays a particular “lubrication” role in the rockfill body. When the content of small-particle-size rubble is relatively low, large-particle-size rubble relies on friction and the “aggregate interlock capacity of rubble” to keep the slope in a stable state. When the content of small-particle-size rubble is high, due to the “lubricating” effect of small-particle-size rubble particles, the relative movement between large-particle-size rubble particles is more likely to occur, resulting in slope instability. Therefore, for the slopes of other rockfill bodies under similar working conditions, the actual situation should be considered to make appropriate adjustments to the early warning standard.
In general, most landslide warning models are based on the displacement–time curve of a single point on the slope. According to the curve change that conforms to the three stages of “initial deformation–uniform deformation–accelerated deformation”, the landslide warning levels are divided into “caution–vigilance–danger”. For example, Wang [
29], He [
30], and Qiang [
31], among others, all studied the landslide early warning index and landslide early warning classification standard based on the “three stages” of the single point displacement–time curve. However, which point of displacement data should be selected as the basis in the process of landslide warning is a fundamental question. To solve this problem, Qiang [
24] proposed that the deformation of each point on the slope is disordered before the basic penetration of the sliding surface. However, after the basic penetration of the sliding surface, the slope begins to slide as a whole. As mentioned above, the rockfill body is loosely piled, and the deformation of adjacent positions on the slope is disordered before the occurrence of landslide. Moreover, under the condition of excavation disturbance, a shallow landslide is likely to occur in a short time, and the time between the connection of the sliding surface and the occurrence of landslide is generally quick. The landslide early warning method established in this paper describes the deformation process of the slope through the area of equal displacement cross-section, avoiding the problem of selecting appropriate points on the slope, which has significant importance for landslide warning of loose accumulation.
6. Conclusions
Due to the granular structure of waste rockfills, it is difficult to find a suitable single-point displacement as an early warning index. The advantage of GB-InSAR full-field scanning is that it helps to obtain slope displacement and other displacement information (the electromagnetic wave reflection intensity of the measured target, etc.) through data processing. In this paper, a landslide early warning system based on the rate of change of the cross-sectional area of equal displacement bodies was proposed for the landslides of unconsolidated coarse sedimentary rocks.
By summarizing the previous landslide monitoring and data processing results, using the normalized rate of change of the cross-sectional area of a −0.5 mm equal displacement body as the early warning criteria, two types of landslides and four levels of early warning classifications were established as follows: (a) Level 1 (green) warning area, < 0.01; (b) Level 2 (yellow) warning area, 0.01 ≤ < 0.02; (c) Level 3 (orange) warning area, 0.02 ≤ < 0.03; and (d) Level 4 (red) warning area, ≥ 0.03.
The validity of the early warning method was verified based on five landslides to compare the predicted and real values in this paper. Further research is needed to take into consideration comprehensive factors such as rock mechanics and rainfall for early warning of rockfill landslides.