Abstract

A series of settlement, cracking, and collapse of the Zhengjiaping slope has been increasingly developing since October 2015 during impoundment at the Dagangshan Hydropower Station. It is a dangerous signal, which means thatlandslide potential will be existed for the Zhengjiaping slope and poses greatthreat to the operation of hydropower station and traffic safety. Related slope protection measures and research studies have been implemented accordingly in the Zhengjiaping slope to prevent and control adverse influences on property security and human life in the reservoir area. In order to understand the geological and hydrological settings of the Zhengjiaping slope and its surroundings, a field investigation was carried out to determine the lithological composition and toppling/sliding deformation characteristics using exploratory adit and site boreholes. The large deformation process in the apparently deformed area was analyzed using borehole inclinometers and global navigation satellite systems (GNSS). It was found that the apparently deformed area zone was characterized by crushed rock masses, with only a small amount of slope deposits and the sliding deformation occurring in Zone I. The deformation process of the reservoir landslide was considered to be a complex integration of the geological effects of various adverse factors. Impoundment and heavy rainfall are the direct causes of sliding deformation. During the preparation of the basic conditions for sliding, lithology, tectonic activity, and artificial disturbances play an important role, including the sliding mass and the sliding surface zone.

1. Introduction

Reservoir landslides are various types of gravitational mass movements of the earth’s surface that occur on the banks of the reservoir area [1, 2]. According to a field investigation, the origin of reservoir landslides is complicated and multiple: (i) slope deposits by alluviation, proluvial action, colluviation, ancient landslide action, or hybrid origin; (ii) broken rock mass by weathering, runoff erosion, tectonic activities, seismic load, artificial disturbances; and (iii) rock mass with the involvement of soft interlayer. Furthermore, differences also exist in proportion of the shallow earth with deformation occurring in the mass above the sliding surface. Integral failure or disintegrated failure has been ever presented in previous occurrences of reservoir landslides because of the differences in the amount and development degree of the sliding surface zone [3, 4].

The preparation and occurrence of reservoir landslides is a very complex and dynamic process that is a long-playing integration of multiple factors. It is generally accepted that impoundment, reservoir level fluctuations, rainfall, and artificial disturbances are the main contributing factors [510]. There have been many studies focusing on the origin of reservoir landslides from slope deposits [1114]. For slope deposits, impoundment and the coupling of heavy rainfall and reservoir level drawdown during the operational period weakens the physical and mechanical properties of the slope material, forms transient seepage, and changes the mechanical behavior, which are highly conducive to the occurrence of reservoir landslides [15]. Understanding how adverse factors contribute to the formation of sliding conditions for the earth's surface and the causative factors that cause sliding deformation is particularly important for landslide prevention and mitigation. A comprehensive understanding of reservoir landslides is important for predicting or identifying precursor phenomena to large reservoir landslides and further generating a series of effective prevention and control measures. Few articles have made specific reference to reservoir rock landslides, but previous studies have made important contributions to the understanding of different aspects related to general rock slide failures. From these studies, there are implications for the understanding of reservoir landslides. Field investigation of geological conditions and characterization of the velocity, distribution, and evolution of movement are fundamental tasks for understanding unstable slope behavior and potential failure [1618]. According to many cases, topography, lithology, geological structure, and tectonic activity play an important role in the development of slope instability, which should usually be a long-term evolutionary process from its inception to catastrophic failure [19]. For rock slopes with potential geologic hazards, the temporal evolution of geomorphology and structure, as well as the influence of lithologic differences and tectonic activity on slope deformation, are of great interest [2025]. For rock slopes that have undergone landslides, mathematical description of the dynamic mechanisms and kinematic processes of disaster triggers is a key point [2632]. To prevent the occurrence of rock slides, methods for evaluating the stability of slopes in reservoir areas have been investigated, such as developmental characterization of cracked rock masses, reliability analysis of rock slopes, or predictive models based on monitoring or statistical analysis [3339].

Reservoir landslides, as a notable geological hazard in reservoir areas, are characterized by their large scale and widespread impact. Reservoir landslides pose a great threat to human beings, environmental stability, property safety, and reservoir operation [2, 15]. The reservoir area in southwest China has steep topography and fragile geology, with broken rock masses in shallow earth and quaternary accumulations covered. The complex geological conditions of the reservoir area are not conducive to slope stability. With the gradual development of cascading hydropower stations, the frequent occurrence of reservoir landslides has seriously threatened the safety of power station buildings and the stability of the reservoir area slopes. At present, the study of reservoir landslides has become one of the important tasks of hydropower development in southwest China.

This paper takes the Zhengjiaping slope as an example to study the deformation characteristics and mechanism of the rock slope deformation occurring in the reservoir area. The geological conditions, lithologic composition, and toppling deformation characteristics revealed by the field survey are the basis for understanding the deformation behavior of the Zhengjiaping slope. The monitoring mainly includes lateral displacement and surface point displacement. The results of lateral displacement and geologic cognition can confirm the existence of basic conditions for sliding deformation, including sliding mass and sliding surface zone. The analysis of surface point displacement is used to identify the factors involved in the deformation of the Zhengjiaping slope, including the adverse factors in the preparation of deformation and causative factors. Then, the influence of unfavorable factors on the process of formation and evolution of sliding deformation is discussed in order to refine some peculiarities of rock slides in the reservoir area.

2. Study Area

The Zhengjiaping slope is located in Xingfu Village, Tianwan Township, Sichuan Province, in the central alpine region between the Tibetan Plateau and the Sichuan Basin in southwestern Sichuan Province, as shown in Figure 1(a). The Zhengjiaping slope is located upstream of Dagangshan Power Station on the right bank of the Dadu River, about 15 km from the dam site of Dagangshan Power Station,as shown in Figure 1(b). According to the relationship between rock inclination and slope ([40], Figure 2), the Zhengjiaping slope belongs to the anaclinal slope. The slope of the Zhengjiaping side slope is 30°–40°, and the surface is mostly covered by accumulated deposits.

Figure 3(a) shows the lithology of the Zhengjiaping slope and the surrounding area. The evident deformation zone in the Zhengjiaping slope was developed in the stratified sandstone and shale of the Baiguowan Formation. The Baiguowan Group is distributed between the East Branch Fault (F1) and the West Branch Fault (F2) of the Dadu River Fracture Zone. The lithology of the Baiguowan Group is mostly siltstone mudstone, with gray sandstone and gray-black carbonaceous shale in the interlayer. The Chengjiang stage, characterized by granite, is distributed on the left bank of the Dadu River and the footwall of the East Branch Fault of the Dadu River Fracture Zone on the right bank. The Zhongdian period is characterized by tonalite gneiss and diorite gneiss, which are distributed in the footwall of the West Branch Fault (F2) of the Dadu River Fracture Zone. The area below 1130 m elevation of the Zhengjiaping slope is locally covered by the accumulation of quaternary system. As shown in Figure 4, this study investigated the deformation quality of the Zhengjiaping slope by exploratory holes and eight site boreholes. The horizontal length of one exploratory borehole was 100 m and the investigation depth of the site borehole was 80 m. The Zhengjiaping slope has strong and weak weathering zones, as well as noncorresponding strong and weak unloading zones, with intact fresh rock at depth. Figure 5 shows the attitude of slope surface, bedding plane, and major joints. Table 1 shows the volume content of the different lithologies, and Figure 6 shows the degree of weathering and unloading of the rock masses at different depths.

The impoundment in Dagangshan Hydropower Station began in May 2015. In October 2015, during the process of rising to normal water level, cracks were found in the middle and leading edge of the Zhengjiaping slope. By the end of February 2016, the width, number, and distribution range of cracks increased. In addition, a localized area of the Zhengjiaping slope has a more pronounced deformation in the 1100–1260 m elevation range. A number of cracks have been intermittently distributed along the back edge, extending about 100 m in length, 10–60 cm in width, and 20–50 cm in vertical offset. On both sides, cracks have a maximum crack aperture of 35 cm and a vertical offset of about 20 cm. On 22 April 2016, a large number of bubbles were found emerging from the water surface of the reservoir near the upstream of the study area.

3. Methods

The distinct deformation zone of the Zhengjiaping slope extends from 1,125 m from the front elevation to 1,330 m from the back elevation along the direction orthogonal to the flow direction of the Dadu River. According to the different deformation characteristics, the study area can be divided into Zone I and Zone II. The volume of Zone I is about 5 × 105 m3 and the volume of Zone II is about 1.5 × 106 m3. According to the geological background of the Zhengjiaping slope, the soil surface process and lithology of the obvious deformation zone are analyzed by field investigation and basic mechanical test. The deep lateral deformation monitoring operation was carried out using three borehole inclinometers (IN01-IN03). The Zhengjiaping slope point deformation monitoring operation was conducted using GNSS (Global Navigation Satellite System), including a continuously operating satellite positioning reference station and several deformation monitoring points. For better analysis of the monitoring data, a three-dimensional coordinate system was used, as shown in Figure 4.

4. Results

4.1. Earth Surface Processes

For the anaclinal slope, the soft plastic features of the rock masses and the thin stratified structure in the Baiguowan Formation favor the development of toppling deformation under long-term maximum principal stresses parallel to the slope surface [41, 42]. Different processes of flexural toppling deformation of bedrock occur on the Zhengjiaping slope, which can be roughly divided into four categories: (1) the disintegration and relaxing deformation zone; (2) the strong toppling deformation zone; (3) the weak toppling deformation zone; and (4) the normal rock stratum zone (Table 2).

Disintegration and the relaxing deformation zone are mainly distributed in the shallow part of the Zhengjiaping slope at a depth of 6.6–46.25 m. As shown in Figure 7, rock stratum has a dip angle 0°–20°, partly gently slope out. The deformation is characterized by strong tensile fracture. The fragmented relaxed rock mass has evidently internal hollow parts, which are filled with large pieces of block broken stones and breccia rock cuttings.

The distribution depths of the strong toppling deformation zone range from 6.6–46.25 m. As shown in Figure 8(a), the stratified rock mass is characterized by strong tensile deformation between the bedding planes, and the extension-shear fracture surfaces of the cutting beds are commonly developed. The deformation behavior of the rocks in the strong toppling deformation zone is complex and can be summarized into three basic failure modes: intercalated shear sliding, discontinuous tensile fracture, and cut-bed tensile sliding (Figure 9). Intercalated shear sliding is a common phenomenon in schistosity planes and weak rock belts. Stratified rock masses with large dip angles have a tendency to dip out of the slope under bending moments from ground stresses, similar to cantilever beams. The flexural toppling deformations have a process of gradual development from the shallow part to the deep part, and the higher elevation causes the more intense deformation. In the local stratigraphy, tensile deformation occurs between the bedding planes under the high flexural toppling deformation. As flexural toppling deformation develops, the level of toppling of the stratified rock mass is increasingly enhanced and the tensile stresses applied to the stratified rock mass are increasingly expanded. If the increased tensile stress reaches or exceeds the tensile strength, tensile fractures occur in the stratified rock mass accompanied by intercalated shear dislocations, resulting in discontinuous tensile fractures. When the development of toppling deformation is sufficient, the inclination of the deformed section changes considerably, the ground stress bending moment increases greatly, and the shear effect is greatly amplified. Under the joint action of shear and tension, cutting-bed tension sliding occurred on the deformed section.

Wedged cracks can be observed where there is evident toppling deformation, including the strong toppling deformation zone and the weak toppling deformation zone. As shown in Figures 8(b) and 10, the existence of intercalated shear sliding and discontinuous tensile rupture rebuild the morphology of the original intact bedrock and form wedged cracks. The weak toppling deformation zone is majorly distributed in the deep part of the Zhengjiaping slope, with a depth more than 38.4 m. This zone is characterized by the weak toppled level and slightly intercalated shear sliding with a dip direction of 40°–60°. Because of the weak extension between bedding planes and the rare development of extension-shear fracture surfaces cutting the rock stratum, the stratified rock mass has a few tensile ruptures in the weak toppling deformation zone.

4.2. Deep Deformation

As shown in Figure 11, the data obtained by borehole inclinometers show lateral displacements of monitoring points for approximately half a year. As shown in Figure 4, the borehole inclinometers IN01 and IN03 were installed in Zone II, and IN02 was installed in Zone I. The evident deformation zone has a sliding surface of development with a depth of approximately 50 m on IN02 and several shearing dislocation zones of development on both borehole inclinometers. There are differences in the timing, number, and development of shear deformation zones at the three borehole locations. The records of the borehole inclinometers indicate that the rock mass of Zone I actually participates in the sliding deformation. However, as can be seen in Table 3, the rock mass below 80 m depth is identified as a weak toppling deformation zone or a normal rock stratum zone, which makes it difficult to produce a sliding surface zone. Therefore, it is certain that the displacement of Zone I is characterized by a combination of sliding and toppling deformation processes that need to be focused on.

4.3. Slope Deformation

As shown in Figure 12, the monitoring data obtained by the GNSS have a record of accumulative displacements of monitoring points for approximately one and a half years. The preconstruction geological survey overlooked the potential of reservoir landslides in the Zhengjiaping slope. When the monitoring displacement began on 21 April 2016, it was approximately one year from the beginning of impoundment, and creep/sliding deformations had already been occurring in the Zhengjiaping slope. Due to the late recording of the initial process of large deformation by the monitoring project, it is difficult to obtain a correlation between deformation and the contributing factors. At the same time, the Zhengjiaping slope was reinforced as it posed a great threat to traffic and residents’ lives. Reinforcement was carried out from mid-May to late August 2016, mainly targeting the excavated area in Zone I. Although the reinforcement had a restraining effect on the deformation of the Zhengjiaping slope, a continuous increase in cumulative displacement was still observed through GNSS.

Firstly, the temporal evolution of displacement and rainfall are compared. This is an oversight of the monitoring program that no rainfall measuring device is arranged in the range within 5 km of the targeted slope. There are two meteorological stations near the Dagangshan Hydropower Station. The Luding Meteorological Station has a distance of approximately 72 km from the upstream of the dam site and the Shimian Meteorological Station has a distance of approximately 40 km from the downstream of the dam site. As shown in Figure 13, the rainfall is mainly focused on a period from June to September, which accounts for 70–80% of the year’s precipitation. It can be used as a reference for the rainfall characters for this case. The two mutations of displacement growths both happened in the flood season, that is, from June to September.

Then, the temporal evolutions of displacement and water level are compared. This two-year record can be superficially divided into two phases. The first phase is approximately from the beginning of deformation monitoring to 31 October 2016, which was characterized by the gradually weakened growth of displacement. The drawdown of the reservoir level was presented in the early monitoring period, which presented the decrease in the deformation velocity, and the accumulative displacements still increased. In late August 2016, the reservoir level began to rise gradually, but the deformation velocities maintained a reducing trend. After that, the deformation velocities tended to be stable and kept a small value, which had little correlation to the changes of the reservoir level. After entering the second phase, the growth rate of displacements tended to be stable; only in two moments, there was an obvious mutation. There was a rise of water level in both displacement mutations. However, when the water level raised at other time, the growth rate of displacement had not been affected. However, the increase of accumulative displacements shows deceleration immediately because of the following drawdown of the reservoir level. As shown in Figure 12, the occurrence time of the increases in the two small-scale velocities is presented in the flood season with the frequent appearance of strong rainfall and increasing of water level. In conclusion, by using a monitoring analysis, impoundment and strong rainfall are considered to be the major disaster-inducing factors for sliding deformation on the Zhengjiaping slope.

5. Discussion and Conclusions

Taking the Zhengjiaping slope as a research example, this paper analyzed the influence of geological conditions and water level on the anaclinal rock slope through field investigation and monitoring analysis. The study area of the Zhengjiaping slope is characterized by a fragmentized rock mass and a small amount of slope deposits. According to the difference in deformation characteristics, the study area can be divided into Zone I and Zone II. Since the impoundment, the deformation characteristics of the Zhengjiaping slope are toppling deformation, sliding deformation, and partial shallow collapse. Based on the degree of deformation and the identification of the sliding surface area, it was found that the sliding deformation only occurred in Zone I. The accumulative displacements on TP01, TP02, and TP03 had been 109 mm, 427.1 mm, and 388.8 mm, respectively, by the end of 6 January 2018.

As shown in Figure 14, the sliding deformation on the Zhengjiaping slope is an integration of many adverse factors. The hydrological activities, lithology, tectonic activities, and artificial disturbances, etc., make great contributions to the formation of the sliding mass and the sliding surface zone, for which a major performance assists in generating fragmentized rock mass in the shallow part. The Zhengjiaping slope has geological settings that are prone to toppling deformation under long-term maximum principal stress parallel to the slope surface. The main contents of the lithology in the sliding mass are sandstone and shale, which have low strength and lamellar structures with a large density of bedding surfaces, which aggravates the development of toppling deformation. Lithology is not the direct reason for the observed deformation after impoundment but has a specific lithologic distribution that ignores the potential for reservoir landslides in the early investigation. As shown in Figure 3(a), the region over the F1 Fault is characterized by argillaceous siltstone with intercalations of sandstone and carbonaceous shale, and the region below the F1 Fault is characterized by granite. The argillaceous siltstone, sandstone, and carbonaceous shale are different types of cementitious sedimentary rock. As shown in Table 4, the argillaceous siltstone and sandstone undergo a decrease in strength when they are saturated. When the water level rises significantly, the large volume of rock is immersed in the reservoir water, which in turn leads to a decrease in the shear capacity of the slope and a faster increase in displacement.

The effects of impoundment and strong rainfall directly cause the sliding deformation, which can be represented as both the physical and chemical actions of water on fragmentized rock mass. Monitoring analysis showed a correlation between the time series of displacement and reservoir water level, but it is not always observed. The phenomenon that the value of cumulative displacement was suffering the change of water level is only existed in the flood season. The evidence of these appearances suggests that both rainfall and water level changes are involved in the deformation. The infiltration of strong rainfall is aimed at a shallow part of sliding mass and sliding surface zone by a rear scarp. The influencing range of impoundment is referred to as the part of the sliding mass between the normal water level (the highest level of water that a reservoir can store under normal operating conditions) and the original water level (the level of the reservoir before the impoundment).

Artificial disturbance has both good and bad effects on slope stability. The good effects referred to the reinforcement on Zone I, which caused the cumulative displacements tend to be converged. The bad effects referred to the slope excavation in the reconstruction of the S211 road. First, the impact loads, which have been produced by the excavation work, cause the generation of new fractures and the extension of existing fractures. The development of fractures leads to a decrease in the rock strength around the excavation face. Second, excavation leads to a change in the distribution of stresses within the slope, which is detrimental to slope stability.

Data Availability

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

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Key R&D Program of China (2017YFC1501102), the National Natural Science Foundation of China (41977229), and the Sichuan Youth Science and Technology Innovation Research Team Project (2020JDTD0006).