Abstract

Reservoir filling can affect the landslide stability and cause the landslide deformation within the reservoir area. In this paper, a physical model test and a series of numerical analyses were combined to investigate the landslide deformation and stability under reservoir filling. The surface deformation, deep displacement, and pore water pressure were recorded during the physical model test. In the model test, the increasing period of the pore water pressure inside the landslide is proposed to be a critical period for the landslide subjected to the reservoir filling. During this period, large landslide deformation occurred. The numerical analyses show that a greater factor of safety (FOS) appeared under a higher water level rising rate or a lower permeability coefficient during the water level rising stage when other variables are fixed, due to the domination of the reservoir buttressing effect which can increase the landslide resistance. During the reservoir maintaining stage, the reservoir water infiltrated into the landslide continuously to cause the matric suction dissipation and pore water pressure increase, which reduced the landslide shear strength and then decreased the landslide stability.

1. Introduction

The Three Gorges Dam is the largest hydropower project in the world. The reservoir impoundment affects the landslide stability within the reservoir area, and more than 5300 landslides have been identified in the Three Gorges Reservoir area [1, 2]. The drawdown of reservoir water is a critical condition for destabilizing landslides [3]. Observations that landslide deformation is associated with the reservoir drawdown have been widely reported [48]. The dynamic seepage pressure pointing outwards the slope is a crucial factor for deforming landslides after water level drawdown [911], especially when the landslide material has a low hydraulic conductivity [12]. Landslides were often found to show a retrogressive failure characteristic under drawdown conditions [10, 1315]. The reservoir filling can also trigger or reactivate landslides. For example, the Qianjiangping landslide in the Three Gorges Reservoir area occurred shortly after the water level reached 135 m, which claimed 24 lives and made 1200 people homeless [16]. Jones [17] observed landslides from 1941 to 1953 in the area around Roosevelt Lake in northeast Washington State, USA. It was found that 49% of the landslides occurred during the reservoir filling stage. Riemer [18] examined 60 publications regarding landslides in reservoir areas and found that 85% of landslides were induced during the reservoir filling or 2 years after the filling.

The field survey and monitoring, physical model tests, and numerical simulation tests are three effective methods to investigate the deformation and stability of landslides influenced by the reservoir impoundment [1922]. Comparatively, the field survey and monitoring are more accurate than the other two methods. However, considerable amount of labor, material resources, and time are needed to install the monitoring equipment and collect data in the field monitoring [23]. Less labor and expense are required in the numerical simulation or physical model tests than in the field monitoring. Moreover, different test conditions and landslide geometries can be implemented for comparative analyses both in the numerical simulation and physical model tests [2426]. Luo and Zhang [26] indicated that the slip surface became deeper inside the slope as the slope gradient increases under drawdown conditions using model tests.

However, there are still very few studies related to landslide stability and deformation under reservoir filling conditions [27], and few researches have analyzed the relation between the seepage field and landslide deformation during the reservoir filling. The mechanism for the landslide stability variation and deformation under reservoir filling also needs exploration. The accuracy of the landslide deformation analysis can be enhanced when considering the landslide deep displacement. The existing studies conducted with model tests mainly focused on the investigation of landslide surface displacement [13, 2830], and the analysis involving the landslide deep displacement in model tests is rare.

In this paper, the physical model test and numerical simulation were combined to study the deformation and stability of a landslide subjected to the reservoir filling. The landslide surface and deep displacement are analyzed in the model test. The correlation between the landslide deformation and the pore water pressure is examined using the grey relational analysis. In the numerical tests, the effect of hydraulic conductivity and water level rising rate on the landslide stability is studied, and the mechanism for the landslide stability variation under reservoir filing is also discussed. The results provide improved insight into the deformation and stability of a landslide subjected to the reservoir filling.

2. Geology of the Majiagou Landslide

The Majiagou landslide was selected as an example to develop the physical and numerical landslide model. The landslide is located on the left bank of the Zhaxi River in the Three Gorges Reservoir (Figure 1). The landslide moved after the reservoir was first impounded in 2003. The toe of this landslide is 135 m above the sea level, and the crown is at an elevation of 280 m. The average slope of the Majiagou landslide is 15°, and the sliding direction is almost perpendicular to the Zhaxi River [31] (Figure 2).

Surficial deposits and sedimentary bedrock are the main material of the Majiagou landslide based on the field investigation. The surficial deposit is composed of gravel mixed with silty clay. The bedrock consists of weathered interbedded gray sandstone and purple-red mudstone of the Jurassic Suining Formation (Figure 3). The purple-red mudstone can be softened by water and has low strength, which causes the silty mudstone beneath the surficial deposits to weather and become the sliding zone eventually [2, 31].

3. Physical Model Test

3.1. Landslide Model
3.1.1. Model Geometry

The physical model is comprised of a model frame, a landslide model, and a monitoring system as shown in Figure 4. A testing frame with a length of 2.8 m, a width of 1.1 m, and a height of 1.5 m was used in the test. The interior of the model frame is waterproof. The front and one side of the frame were made of transparent tempered glass for observing the water level and landslide deformation inside the frame. The other side and back of the frame were made of steel.

3.1.2. Landslide Materials

The landslide model consists of a sliding mass, sliding zone, and bedrock. The landslide model is 2.7 m long and 1 m wide. The sliding mass and sliding zone have an average thickness of 0.32 m and 40 mm, respectively (Figure 4).

The material properties used for constructing the landslide model were determined by the law of similitude [32]. The landslide model was required to be similar to the prototype in terms of the geometry, material properties, and boundary conditions. The ratio of the prototype to model in the certain parameters (e.g., length, density, and cohesion) is defined as the similar ratio. The parameters between the Majiagou landslide and the physical model were converted based on the Buckingham π theorem [33]. The similar ratio of gravity and density was set as and , respectively, and the similar ratio of geometry was set as in this study. The similar ratio of other parameters was derived using the dimensional method. The results are shown in Table 1.

The bedrock of the landslide model was constructed with masonry and cement mortar. The sliding mass was made from a mixture of sand, clay, water, and bentonite with a ratio of 28.5 : 62.5 : 8 : 1, and the sliding zone was composed of glass beads, clay, and water in proportions of 60 : 32 : 8. These ratios were obtained after a series of repeated matching tests. The clay both used in the sliding mass and sliding zone was obtained from the Majiagou landslide. The material properties of the Majiagou landslide and the model are listed in Table 2.

The clay used in the tests was dried, crushed, and screened by a 2 mm mesh before the mixture of the materials. To obtain homogeneous and mixed materials, a concrete mixer was employed during the mixing process and the mixing time was set at 20 min each time. In the construction of the landslide model, the layered construction and compaction were applied. The sliding zone was constructed in a 40 mm layer. The sliding mass was constructed in 50 mm thick layers. Each layer was compacted by hand tamping 160 times with a rubber hammer. The landslide model was hand trimmed to obtain the designed shape before finishing the construction.

To simulate the landslide subjected to the reservoir filling, the reservoir water level was gradually raised from 0.4 m (the landslide toe) to an elevation of 0.89 m over 30 min and maintained at that elevation for 120 min.

3.2. Monitoring

A RIEGL VZ-400 3D laser scanner was used to measure the surface deformation of the model (Figure 4). White spherical pushpins with a diameter of 8.5 mm were placed as surface monitoring points on the landslide model surface (Figure 5). The 3D laser scanner was set about 0.5 m in front of the model frame, and it scanned the surface every 5 min to obtain point cloud data of the surface. The movement of the model surface was obtained after the surface monitoring point locations were further analyzed.

Seven pore water pressure gauges (27 mm in diameter and 15 mm thick) were placed in the sliding zone along the center axis of the landslide model, as shown in Figures 4 and 5.

Two flexible inclinometer probes C1 and C2 were set in a vertical orientation, perpendicular to the center axis of the model to monitor the landslide deep displacement (Figures 4 and 5). The measurement device for the deep displacement monitoring includes a flexible probe, controller, and PC software (Figure 6). The measurement principle of the probe and its verification in physical model tests were presented by Zhang et al. [34]. The flexible probe contains the gravitational acceleration measurement units (Figures 5 and 6). The horizontal displacement of each measurement unit in the probe was obtained based on the angle measurement using the gravitational acceleration measurement units. The collected displacement data were sent to the computer by the controller through the wireless communication [34].

3.3. Results of the Model Test
3.3.1. Pore Water Pressure

Figure 7 shows the variation of the pore water pressure and water level versus time in the model test. Lagged responses at all pore water pressure gauges are observed compared with the water level. The pore water pressure of P1 responds earliest and increases 13 min after the water level has risen. P2 has the largest value of pore water pressure (2.92 kPa). The gauge of P7 has no change through the whole test which indicates the water does not reach the elevation of P7. The pore water pressure at P1, which has the lowest elevation among the seven gauges, does not rise as expected to exceed the value of P2. The rising rate of P1 decline after 20 min, which is likely due to gauge failure. After 65 min, the pore water pressures at all gauges become stable and keep almost constant until the test end.

3.3.2. Surface Displacement

The 3D laser scanner cannot capture the landslide surface movement under the water during the test. To better describe the feature of surface deformation, the landslide surface above the water is divided into three zones (Figure 5). The average displacements of monitoring points in the three zones versus time are shown in Figure 8. Zone 1, zone 2, and zone 3 start to deform 20 min, 31 min and 35 min after the water rises, respectively. Zone 1 which is closest to the reservoir has the largest average displacement (11 mm), and the deformation of zone 3 is the smallest (1 mm). The displacement rates of these three zones become smaller after 65 min.

3.3.3. Deep Displacement

Figure 9 presents the variation of the deep displacement in the flexible inclinometer probes C1 and C2. The displacement of C1 (located in zone 1) increases obviously during 15 to 60 min, which is almost consistent with the surface displacement increase of zone 1. C2 has less total displacement compared with C1. Similar to C1, the rapid increase period (30~60 min) of the displacement in C2 (located in zone 2) is also in accord with the surface displacement increase of zone 2. Both in C1 and C2, the deformation at the deeper section (140 mm~300 mm) is a bit larger than that at the shallow section (0~140 mm) after 15 min and 30 min, respectively. The reason may be that the sliding zone containing the glass beads has smaller shear strength than the sliding mass under the softening effect of reservoir water. Then, the deeper soil which is closer to the sliding zone deforms greater than the soil in the shallow area.

3.3.4. Overall Deformation Distribution

Based on the deep and surface displacement, the overall deformation contour of the three zones can be obtained using the triangulation and linear interpolation [35, 36]. The displacements of C1 and C2 (located in zone 1 and zone 2, respectively) can be used as the deep displacements of zone 1 and zone 2 in the contour plotting of the overall deformation. However, the deep displacement of zone 3 is unknown, as few surface deformations appeared in zone 3 (less than 1 mm), and the surface displacement is close to the deep displacement according to zone 1 and zone 2. Therefore, the deep displacement was set to be the same with the surface displacement in the overall deformation contour plotting of zone 3. The overall deformation contour of zone 1, zone 2, and zone 3 and saturation lines at different times are shown in Figure 10.

The displacement beneath the water (grey portion in Figure 10) was not presented due to the limitation of the 3D laser scanner. From Figure 10(a), few deformations (less than 0.8 mm) are found in zone 1 at 20 min, and the deformation of zone 2 and zone 3 is almost 0 mm. Then, the displacement of zone 1 increases to about 4 mm at 30 min when the water level reaches the highest elevation (0.89 m) (Figure 10(b)). At 65 min, the largest deformation (occurs in zone 1) rises to about 13 mm as the saturation line becomes almost parallel to the water level line. However, from 65 min to 150 min, the saturation line almost keeps constant and few displacements (about 2 mm) appear. The result shows that most landslide deformations occurred during the first 35 min (30~65 min) of the water level maintaining stage.

4. Numerical Tests

4.1. Numerical Model

The seepage field of the landslide subjected to the first reservoir impoundment was simulated by the two-dimensional finite-element seepage module SEEP/W of the GeoStudio 2007 [37]. The Morgenstern-Price method is a vigorous limit equilibrium method of slices for slope stability calculation. It completely satisfies the equilibrium condition and involves the least numerical difficulties [38]. The factor of safety (FOS) of the Majiagou landslide was then calculated using the Morgenstern-Price method in the SLOPE/W module, which employed the transient seepage results obtained from SEEP/W. The vertical profile of the Majiagou landslide (Figure 3) was selected for establishing the numerical model.

The established finite-element model composed of the sliding mass, sliding surface, and bedrock is partitioned into 786 elements and 843 nodes (Figure 11). In the transient seepage analysis, the reservoir water is raised from 145 m to 175 m and then the water level is maintained at 175 m. The boundary conditions of the model are set as follows: The boundary for the varying water head is at the front edge of the sliding mass from the elevation of 145 m to 175 m. The groundwater level at the trailing edge of the landslide is set as the constant water head with an elevation of about 232 m based on the field groundwater monitoring data. The bottom line of the model is set as the impermeable boundary.

4.2. Test Scheme

To analyze the effect of the reservoir filling on the seepage field and stability of the landslide, four distinct permeability coefficients (, , , and ) for the sliding mass and seven different rising rates (0.167 m/day, 0.2 m/day, 0.25 m/day, 0.33 m/day, 0.5 m/day, 1 m/day, and 2 m/day) were considered to conduct the comparative study. Before the transient seepage simulation, the steady-state seepage analysis was conducted to obtain the seepage field of the landslide under the water level of 145 m. Then, the result was used as the initial condition for the transient seepage analysis under the reservoir impoundment (water level rising from 145 m to 175 m and maintaining at 175 m).

The estimated soil-water characteristic curve and the hydraulic conductivity function of the sliding mass are shown in Figure 12. The hydraulic permeability functions of the sliding mass were evaluated using the Van Genuchten method [39]. The parameters for estimating the soil-water characteristic curve and the hydraulic conductivity function presented by Paronuzzi et al. [22] were applied in this study (Table 3). In the landslide stability evaluation, the parameters of the landslide material for the stability calculation are shown in Table 4.

4.3. Results of the Numerical Simulation
4.3.1. Transient Seepage Field

Figure 13 shows the saturation lines of the landslide under different permeability coefficients as the reservoir water rises from 145 m to 175 m over 90 days and is then maintained. From Figure 13, the permeability of the landslide material affects the transient seepage field inside the landslide. The saturation lines at the slope toe are flatter under greater permeability coefficients during the water level rising (Figures 13(a) and 13(b)), and the front part of the saturation lines tends to show a concave shape under a smaller permeability coefficient ( and ) (Figures 13(c) and 13(d)). The concave saturation lines indicate that the groundwater level lags behind the reservoir water level. During the water level maintaining stage, less time is needed for the groundwater equilibrium under a high permeability coefficient. For example, the groundwater almost follows the change of the reservoir water level under a permeability coefficient of (Figure 13(a)). However, under a permeability coefficient of , the groundwater cannot reach equilibrium though the reservoir is maintained for 3910 days (Figure 13(d)).

4.3.2. Stability of the Landslide

The FOS variation of the landslide during the water level rising and maintaining periods is presented in Figures 14 and 15. During the water level rising period, a greater FOS appears under a higher water level rising rate or a smaller permeability coefficient when other conditions are constant. In the water level maintaining stage, all FOSs decrease as the reservoir water infiltrates into the landslide and the FOSs drop to the same value under the same hydraulic conductivity. The FOS needs more time to decline to be stable under a smaller permeability coefficient due to a smaller water infiltration rate. From the deformation analysis in the physical model test, most landslide deformations occurred in the water level maintaining stage, which is consistent with the result here.

5. Discussion

5.1. Relation between Landslide Deformation and Seepage Field

The correlation between the landslide deformation and its impact factors (such as the reservoir fluctuation or rainfall) can be evaluated using the grey relational analysis [10]. The grey relational degree which ranges from 0 to 1 can be obtained by the grey relational analysis. The closer the grey relational degree to 1, the higher the correlation between the landslide deformation and impact factors. In this study, the pore water pressure is the only impact factor in the relational analysis. To find the correlation between the landslide deformation and the pore water pressure, four steps are described as follows: (1)Determination of Calculation Sequences. The landslide displacement rate is set as the primary sequence. The only one subsequence is the variation rate of pore water pressure. (2)Normalization of the Value in Each Sequence. The data of each sequence are equalized using where ; ; , the number of influencing factors ( in this study); and , the number of data points(3)Calculation of Correlation Coefficients between Primary-Sequence () and Subsequence (). The calculation method is shown as follows: where is the correlation coefficients between sequences and sequences and is the resolution coefficient, which is set as 0.5 normally(4)Calculation of the Grey Relational Degree. The grey relational degree can be calculated from

P2, P4, and P5 are located in the front part of zone 1, zone 2, and zone 3, respectively. The pore water pressures of P2, P4, and P5 (Figure 7) are selected to conduct the grey relational analysis with the landslide surface deformation of zone 1, zone 2, and zone 3 (Figure 8), respectively. After the calculations, the grey relational degrees of the three pairs of analyses are 0.8, 0.8, and 0.7. The result indicates that the landslide surface deformation is related to the pore water pressure variation.

To further study the relation between the landslide deformation and the pore water pressure, the average displacement of the landslide surface (zone 1 to zone 3) was calculated and presented in Figure 16(c), which is compared with the pore water pressure variation of P2 to P6 (Figure 16(a)) and the average variation trend of these five pore water pressure curves (Figure 16(b)).

The grey relational degree between the average pore water pressure and the surface average displacement is 0.92 by the calculation. In addition, the main increasing stage (from 15 min to 65 min) of the landslide surface displacement (Figure 16(c)) is consistent with the rising period of the pore water pressure (P2 to P6) (Figure 16(a)). It indicates that the deformation is closely related to the increase of the pore water pressure. The increasing period of the pore water pressure inside the landslide is proposed to be a critical period for the landslide subjected to the reservoir filling. During this critical period, large landslide deformation appears, especially in the area close to the reservoir (zone 1).

5.2. Mechanism for the Landslide Stability Variation

Shear strength of the unsaturated soil can be expressed as [40] where is the shear stress, is the effective cohesion, is the normal total stress, is the pore air pressure, is the effective friction angle, is the pore water pressure, is the matric suction of the soil on the plane of failure, is the volumetric water content, is the residual volumetric water content, and is the volumetric water content at a saturation of 100%.

This equation is applied at the landslide stability calculation in the SLOPE/W module. The first part of equation (1) is the saturated shear strength when the pore air pressure is equal to the pore water pressure [40]. During the reservoir water infiltration, the matric suction dissipation and pore water pressure increase reduced the shear strength and then caused the deformation or reduced the FOS (Figures 10, 14, and 15). The buttressing effect of the reservoir can explain the increase of landslide stability during the water level rising [22]. The buttressing force caused by the reservoir impoundment acting on the landslide surface increased the landslide resistance and thus raised the FOS during the water level rising period (Figures 14 and 15).

5.3. Interpretation of the Landslide Stability Variation and Deformation

The monitoring points P1~P6 were set in numerical models for recording the pore water pressure during the numerical tests (Figure 11). Figure 17 presents the pore water pressure variation of P1~P6 under the water level rising rate of 2 m/day and permeability coefficients of and . The evolution trend of the pore water pressure in numerical tests is similar to that of the model test in Figure 7. Both in the model and in the numerical tests, the pore water pressure in the landslide front part (P1, P2, P3, and P1 and P2) increases quickly to the maximum value and then keeps constant and the pore water pressure in the rear part (P6, P7, and P5 and P6) almost has the same variation characteristic (Figures 7 and 17). The pore water pressure of P1 and P2 under a higher permeability coefficient () changes synchronously with the water level (Figure 17(a)). However, under a lower permeability coefficient (), the pore water pressure needs much more time to rise and get stable (Figure 17(b)).

Buttressing effect of the reservoir is the same in all numerical tests. A larger permeability coefficient or a slower water level rising rate allows more water to flow into the landslide, which is negative to the landslide stability during the reservoir filling. Therefore, under a greater permeability coefficient or a lower water level rising rate, a smaller FOS occurs in the reservoir filling stage (Figure 14). In the maintaining stage, the buttressing effect keeps constant but the reservoir water flows into the landslide continuously. Therefore, all FOSs decline due to the matric suction dissipation and pore water pressure increase. In the physical model test (), less deformation occurred in the water level rising stage and most deformation appeared during the early water level maintaining period. It agrees well with the numerical results and also can be explained by the reservoir buttressing effect and shear strength reduction caused by the reservoir water infiltration.

6. Conclusions

In this study, the physical model test and numerical tests were combined to study the deformation and stability of a landslide subjected to the reservoir filling. The Majiagou landslide in the Three Gorges Reservoir area was selected as an example to develop the physical and numerical models. During the physical model test, the surface deformation, deep displacement, and pore water pressure were recorded and the relationship between the landslide deformation and the seepage field was discussed. In the numerical tests, the influence of the permeability coefficients and water level rising rates on the seepage field and landslide stability was analyzed. The main findings are summarized as follows:

The landslide deformation was found to be associated with the pore water pressure variation during the reservoir filling using the grey relational analysis. In an area, when the pore water pressure started to rise, the deformation increased at the same time. Then, the deformation rate declined following the equilibrium of the pore water pressure. The main increasing stage of the landslide surface displacement was consistent with the general rising period of the pore water pressure inside the landslide. The increasing period of the pore water pressure within the landslide under reservoir filling was proposed to be a critical period, during which large landslide deformation occurred.

The mechanism for the landslide deformation and stability variation subjected to the reservoir filling was analyzed. The impounded reservoir provided the buttressing force acting on the landslide surface, which increased the landslide resistance and thus raised the landslide stability. However, the reservoir water infiltration caused the matric suction dissipation and pore water pressure increase within the landslide soil, which reduced the landslide shear strength and then reduced the landslide stability. The landslide stability subjected to the reservoir filling was governed by these two effects. When the former effect prevailed over the latter, the landslide stability increased and vice versa.

The landslide stability under different hydraulic conductivities and water level rising rates subjected to the reservoir filling was discussed. The buttressing effect of the reservoir enlarged gradually over the water level rising period. During this stage, less water flowed into the landslide under a higher water level rising rate or a lower permeability coefficient. Therefore, the buttressing effect of the reservoir was dominated, and thus, a greater FOS appeared. However, during the maintaining stage, the reservoir water flowed into the landslide continuously to cause the matric suction dissipation and pore water pressure increase, which reduced the landslide shear strength and then decreased the landslide stability.

Data Availability

The data used in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

This study is funded by the Key Program of the National Natural Science Foundation of China (41630643), the National Key Research and Development Program of China (2017YFC1501302), and the Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan) (CUGCJ1701).