Next Article in Journal
Application of Multi-Source Data Fusion Method in Updating Topography and Estimating Sedimentation of the Reservoir
Previous Article in Journal
Exploring the Dominant Runoff Processes in Two Typical Basins of the Yellow River, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deformation and Stability Characteristics of Layered Rock Slope Affected by Rainfall Based on Anisotropy of Strength and Hydraulic Conductivity

1
School of Geosciences and Info-Physics, Central South University, Changsha 410083, China
2
Key Laboratory of Metallogenic Prediction of Nonferrous Metals and Geological Environment Monitoring Ministry of Education, School of Geoscience and Info-physics, Central South University, Changsha 410083, China
3
Hunan Key Laboratory of Nonferrous Resources and Geological Hazards Exploration, Changsha 410083, China
4
Yunnan Diqing Nonferrous Metals Co., Ltd., Diqing 474400, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(11), 3056; https://doi.org/10.3390/w12113056
Submission received: 9 October 2020 / Revised: 25 October 2020 / Accepted: 27 October 2020 / Published: 30 October 2020
(This article belongs to the Section Hydrology)

Abstract

:
The strength and hydraulic conductivity anisotropy of rock slopes have a great impact on the slope stability. This study took a layered rock slope in Pulang, Southwestern China as a case study. The strength conversion equations of the seriously weathered rock mass were proposed. Then, considering the anisotropy ratio and anisotropy angle (dip angle of bedding plane) of strength and hydraulic conductivity, the deformation and stability characteristics of rock slope were calculated and compared with field monitoring data. The results showed that the sensitivity analysis of strength and hydraulic conductivity anisotropy could successfully predict the occurrence time, horizontal displacement (HD), and the scope of the rock landslide. When the anisotropy ratio was 0.01 and the dip angle was 30°, the calculated HD and scope of the landslide were consistent with the field monitoring data, which verified the feasibility of the strength conversion equations. The maximum horizontal displacement (MHD) reached the maximum value at the dip angle of 30°, and the MHD reached the minimum value at the dip angle of 60°. When the dip angle was 30°, the overall factor of safety (FS) and the minimum factor of safety (MFS) of the rock slope were the smallest. By assuming that the layered rock slope was homogeneous, the HD and MHD would be underestimated and FS and MFS would be overestimated. The obtained results are likely to provide a theoretical basis for the prediction and monitoring of layered rock landslides.

1. Introduction

Rainfall is one of the main factors that trigger landslides [1]. In the past two decades, research on rainfall-induced soil landslides has been gaining popularity [2,3,4,5,6], but research on rock landslides is less reported [1,7,8,9,10,11,12]. According to available statistical data, more than 90% of rock slope failures are related to fluid hydraulic conduction [13]. The strength parameters of rock slopes such as Poisson’s ratio and elastic modulus also influence the stability of rock slopes [14]. Collins and Znidarcic [14] pointed out that the depths of slope failures were not only controlled by hydraulic characteristics, but also by strength characteristics; thus, both factors should be considered to assess the stability of slopes.
Anisotropy is everywhere but isotropy is rare [15]. Barton and Quadros [15] questioned whether the a priori assumption of homogeneous-isotropic-elastic behavior had any significant place in the scientific practice of realistic rock mechanics. Hence, the physical and mechanical behaviors of anisotropic rocks, such as sandstone [16,17] shale [18,19,20], slate [21,22], and phyllite [23,24], have become a hot topic over the last decades [25]. Seriously weathered carbonaceous slate, as a kind of layered rock, has an obvious bedding structure [26]. The layered structure of slate rock induces strength anisotropy. The rock strength in the direction parallel to the layers is considerably smaller than in other directions [27]. Based on the triaxial compression test of shale, Geng et al. [28] found out that the brittleness anisotropy of Longmaxi shale samples was dependent on the bedding angle. Ding et al. [25] conducted laboratory tests of anisotropic slate with various dip angles and found out that the uniaxial compression strength exhibited a typical U-type trend with the variance of dip angle, the elastic modulus first decreased and then increased, and the variation of Poisson’s ratio showed an opposite trend with the variation in elastic modulus. Based on uniaxial and triaxial compression tests, Gholami and Rasouli [22] found out that the strength of slate exhibited obvious U-shaped anisotropy with an increasing dip angle. Chen et al. [21] found out that the apparent elastic modulus increased remarkably with the dip angle, and the compressive strength at any confining pressure varied in a typical U-shaped trend. The above research mainly focuses on the anisotropy of layered rock mass, but the strength anisotropy of actual rock slope is less studied.
The deformation and hydraulic conduction of layered rock masses display anisotropic characteristics [29]. The ratio of the hydraulic conductivity in the parallel direction to the vertical direction of the bedding plane in the western foothills of Taiwan was 10 to 100 [1]. Dong et al. [11] found out that the anisotropy ratio and dip angle of the bedding plane had significant effects on the pore pressure distributions and slope stability; however, the above research is restricted to dip angles of 30, 0, and -30°, and the influence of anisotropy ratio was ignored. Yeh et al. [1] also found the heterogeneity in hydraulic conduction and strength for layered rock masses. It was considered that the effects of anisotropy of hydraulic conductivity and strength on groundwater flow and slope stability were mainly reflected in the direction of the bedding plane. However, the rock slope was seriously weathered, while the Poisson’s ratio and elastic modulus used for numerical simulation were from intact rocks, which would lead to inaccurate results. Yu et al. [3] and Xia et al. [9] studied the seepage, deformation, and stability of soil slopes based on the anisotropy of hydraulic conductivity; however, the influence of the anisotropy of elastic modulus and Poisson’s ratio on the slope deformation and stability was ignored.
In view of the shortcomings of previous studies, the theory of fluid–solid coupling and the unsaturated seepage are introduced. Based on rock uniaxial compression tests, the strength conversion equations of weathered rock are proposed. A seriously weathered layered rock slope in the Pulang region of Southwestern China was selected as the study objective. By considering different anisotropy ratios and anisotropy angles of hydraulic conductivity and strength, numerical simulations were carried out to analyze the stability and deformation characteristics of rock slope. Major parameters, including the horizontal displacement (HD), maximum horizontal displacement (MHD), the factors of safety (FS), the minimum factors of safety (MFS), and factors of safety of landslide (FSL) were proposed to characterize the slope deformation and stability pattern. Finally, by comparing the actual detection data of the rock landslide with the numerical simulation results, the feasibility of the strength conversion formula and the accuracy of the landslide prediction were verified. The research results could provide some references for monitoring and predicting rock landslides.

2. Materials and Methods

2.1. Control Differential Equation

To simulate the seepage processes of rock slopes, rainfall infiltration was regarded as the variable flow boundary of the unsaturated zone, and a saturated–unsaturated seepage theory was applied to perform the simulations. Under the principle of mass conservation, the saturated–unsaturated seepage control differential equation is as follows [3,9,30]:
x ( k x H x ) + y ( k x H y ) + Q = m w γ w H t
where H is the total head, the unit is m; x and y are the coordinates in the direction of x and y; kx is the hydraulic conductivity in the x-direction, the unit is m/s; ky is the hydraulic conductivity in the y-direction, the unit is m/s; Q is the applied boundary flux, mw is the slope of the storage curve, t is the time, the unit is s; γw is the unit weight of water, the unit is N/m3.
Applying the Galerkin method of weighted residual to the governing differential equation, the finite element seepage equation in two-dimension can be derived as [3]:
τ A ( [ B ] T [ C ] [ B ] ) d A { H } + A ( λ N T N ) d A { H } , t = q τ L ( N T ) d L
where [B] is the gradient matrix; [C] is the element hydraulic conductivity matrix; [H] is the vector of nodal heads; <N> is the vector of interpolating function; q is the unit flux across the edge of an element; τ is the thickness of an element; λ is the storage term for a transient seepage equal to mwγw; A is a designation for summation over the area of an element; and L is a designation for summation over the edge of an element. In an abbreviated form, the finite element seepage equation can be expressed as:
[ K ] { H } + [ M ] { H } , t = { Q }
where [K] is the element characteristic matrix; [M] is the element mass matrix; and {Q} is the element applied flux vector.

2.2. Theory of Fluid–Solid Coupling

For a two-dimensional space, the incremental stress–strain relationship can be written as [3]:
{ Δ σ } = [ D ] { Δ ε } [ D ] { m H } { u a u w } + { Δ u a }
where ε is the normal strain; γ is the engineering shear strain; σ is the normal stress, the unit is kN; τ is the shear stress, the unit is kN; ua is the pore-air pressure, the unit is kPa; uw is the pore water pressure, the unit is kPa; [D] is the drained constitutive matrix; and {mH} can be expressed as:
{ m H } T = 1 H 1 H 1 H 0
where H’ is the unsaturated rock modulus for rock structure with respect to matric suction (uauw).

2.3. Establishment of the Factor of Safety for Unsaturated Rock

This study adopted a Morgenstern–Price (M–P) method based on the limit equilibrium theory to calculate the factors of safety. The M–P method [31] for the analysis of rock slope stability offered various user-specified inter-slice force functions and considered both the shear and normal inter-slice forces, which made the analysis more consistent with the actual engineering situation [32,33,34]. It was found that the method strictly satisfied the force balance and torque balance, and the calculation accuracy was found to be high.
Taking the rock slice as the research object, as shown in Figure 1, when the boundary conditions E0 = 0 and En = 0, according to the balance of the force in the normal and tangential directions of the slip surface of the rock slice, it can be obtained:
S i = ( N i tan φ + c b sec α ) / F s
N i = ( W i + λ f i 1 E i 1 λ f i E i + Q i cos ω i + J i sin β i ) cos α + E i E i 1 + Q i sin ω i J i cos β i ) sin α i
where Ni is the effective normal force on the slip surface of the rock slice and the unit is kN; Si is the shear force on the slip surface of the rock slice and the unit is kN, and the safety factor Fs can be derived as:
F s = i = 1 n 1 ( R i j = i n 1 ψ j ) + R n i = 1 n 1 ( T i j = i n 1 ψ j ) + T n
where Ri is the resistance force, R i = [ W i cos α i + Q i cos ( ω i α i ) + J i sin ( β i α i ) ] tan φ i + c i b i sec α i and the unit is kPa; Ti is the sliding force, T i = W i sin α Q i sin ( ω i α i ) + J i cos ( β i α i ) ; ψi is the transfer coefficient, ψ i 1 = ( sin α i λ f i 1 cos α i ) tan φ i + ( cos α i + λ f i 1 sin α i ) F s ( sin α i 1 λ f i 1 cos α i 1 ) tan φ i 1 + ( cos α i 1 + λ f i 1 sin α i 1 ) F s ; ci is the effective cohesion for every rock slice and the unit is kPa; φi represents the shear strength angle for every rock slice, the unit is the degree. Ji is the seepage force of the rock slip which adopts a simplified seepage field treatment method, assuming that its position of action is located at the center of gravity of the rock below the underwater line, and the distance from the center of the slip surface of the rock strip is hwi/2, and hwi is the distance from the underwater line to the slip surface. Ei and Ei−1 are the horizontal effective forces between the two sides of the rock slice and the unit is kN; λfiEi and λfi−1Ei−1 are the shear forces between the two sides of the rock strip, fi is the function of the inter-strip force, and λ is the proportional coefficient. The distances between the action positions on both sides and the center of the slip surface of the rock strip are Zi and Zi−1, respectively; Wi is divided into two parts based on the underwater line: Wi = Wi1 + Wi2, where Wi is the gravity of the rock above the underwater line and Wi2 is floating weight of the rock below the underwater line. Qi is an external force on the rock surface; ωi is the angle between the direction of the external force and the normal direction; αi is the angle of the slip surface.

2.4. Establishment of the Numerical Calculation Model

The layered rock slope in the Pulang region of Southwestern China is shown in Figure 2. In Figure 3, the stratum lithology of the slope is seriously weathered carbonaceous slate with a layered structure. This finite element model of the slope body is 260 m high and 540 m long. The left and right borders extend horizontally for 50 m, and the lower border extends downward for 50 m. The grid model is divided into 18,233 nodes and 18,196 elements. To study the characteristics of the deformation at different positions, three monitoring lines were set at x = 150 m (top of the slope), x = 250 m (middle of the slope), and x = 330 m (bottom of the slope). This rock slope experienced a landslide on 4 July 2019 and the scope of the landslide is shown in Figure 2a. Displacement monitoring stations were set between the bottom and the middle of the slope: They were located 1 m below the slope and were set at x = 345 m (BPK2-1), x = 370 m (BPK2-2), and x = 395 m (BPK2-3).
In Figure 3, the boundary conditions are as follows: AB and FD are the fixed water level boundaries of 175 and 90 m, respectively; CD indicates the boundary of the rainfall infiltration; BC is the free boundary with a limited flow; AF is an impervious boundary. Since the rock landslide occurred on 4 July 2019, to restore the rainfall boundary, the average daily rainfall from 1 March to 2 August 2019 (155 days) in the Pulang region was adopted as the simulated rainfall conditions, as detailed in Figure 4.

2.5. Determination of the Maximum Initial Matric Suction

Since the studied rock slope had a strong degree of weathered and the rock layers were so intensely fractured that the properties of a porous medium could be assumed, it had a certain degree of matric suction. The extent of the initial maximum matric suction may have varied, and the seepage states and stability levels of slopes may have also been quite different. Therefore, to make the initial matric suction more in line with the actual situation (refer to previous studies [9,35]), this study set the different maximum initial matric suctions as −50, −75, and −100 kPa, respectively, and the annual average rainfall of the Pulang area as 0.06 mm/h in the steady-state comparative analysis process. As detailed in Figure 5, the changes in the matric suction of the average annual rainfall were similar to the setting of the maximum initial matric suction of −75 kPa. Therefore, the initial maximum matric suction of all this study’s simulation experiments was set as −75 kPa, which was determined to be in line with the real situation in the study area.

2.6. Establishment of Unsaturated Permeability Coefficient

Due to the difficulties of measuring the unsaturated characteristic curves of rock masses, the rock mass was assumed as porous media or unsaturated soils, and the water retention curves of classic unsaturated porous media [36,37,38,39], and relative hydraulic conduction models [40] were applied to describe the unsaturated flow in fissures and weakly permeable rocks. Therefore, this study adopted the van Genuchten model to describe unsaturated features of rock masses. The van Genuchten model is as follows [37]:
θ = θ r + θ s θ r [ 1 + ( p / a ) n ] m ,
k = k w [ 1 | p / a | n 1 ( 1 + | p / a | n ) m ] 2 ( 1 + | p / a | n ) m 2 .
where p is the suction and the unit is kPa; θ indicates the adjusted volumetric water content, the unit is m−1; θr is the residual volumetric water content, the unit is m−1; θs indicates the saturated volumetric water content, the unit is m−1; a is the fitting parameter closely related to the air-entry value of the unsaturated rock mass and the unit is kPa; n and m (m = 1 − 1/n and n > 1) are fitting parameters that control the slope at the inflection point in the volumetric water content function [18]. kw is the saturated hydraulic conductivity, unit is m/s; k is the adjusted hydraulic conductivity, the unit is m/s.
Parameters such as a, n, m, and θr are known to have major influences on the unsaturated flow of rock slope [9]. The slope in the present study is located on the west side of the slope in a previous study [9], and the difference of unsaturated parameters between the two rock slopes is negligible. Therefore, the parameters of the water retention curve of the present slope were referred from a previous study [9]. The unsaturated parameter values are shown in Table 1 [9]. Accordingly, the water retention curve of the rock mass is shown in Figure 6.

2.7. Definition of Anisotropy and the Calculation Conditions

The different dip angles of the bedding plane of the rock slopes may result in different hydraulic conductivity [9]. In this study, the hydraulic conductivity matrix [C] is expressed as follows:
[ C ] = [ C 11 C 12 C 21 C 22 ]
In Equation (11), C11 = kxcos2β + kysin2β, C22 = kxsin2β + kycos2β, and C12 = C21 = kxsinβcosβ + kysinβcosβ. The kx and ky, as well as the dip angle β, were defined according to Figure 3. In the present study, kx represents the horizontal hydraulic conductivity; ky is the vertical hydraulic conductivity; β is the direction between kx and x-axis. Therefore, if β = 0°, then [C] will be reduced to:
[ C ] = [ k x 0 0 k y ]
Equation (12) was adopted in this study [4] with only the anisotropy ratio kr = ky/kx considered. However, for layered rock masses, due to the existence of dip angles of the bedding plane, the conditions that the anisotropy angle β was not equal to 0 were also considered. Therefore, in this study, to restore a reasonable rock formation state, anisotropy angles ranging between 0 and 90° in the layered rock slope in the Pulang area were taken into consideration. As shown in Figure 3, the anisotropic angle of hydraulic conductivity was equivalent to the dip angle of the layered rock slope.
In addition, the dip angle of the rock slope also affected the strength parameters of the rock mass, such as elastic modulus and Poisson’s ratio. The rock was processed into rock specimens with dip angles of 0, 15, 30, 45, 60, 75, and 90°, as shown in Figure 7. The WAW-1000kN microcomputer-controlled electro-hydraulic Servo universal testing machine was used to adopt the uniaxial compression test.
In Figure 8, through the uniaxial compression experiment, the elastic modulus of the rock specimens with different dip angles is in the shape of “U”, reaching the maximum when the dip angle is 90°, and reaching the minimum when the dip angle is 30° [21]. The influence of dip angle variance on Poisson’s ratio was minimal [21].
The strength parameters of the rock mass were not only affected by the dip angle, but also by the integrity of the rock mass [42]. Since the stratum lithology of the rock slope in the present study belonged to seriously weathered carbonaceous slate, the strength parameters of the intact rock could not be directly applied to the numerical simulation calculation. According to damage mechanics [42,43], the rock masses with different weathered states had different elastic modulus. There was an assumption that with the change of the dip angle, the change trends of the elastic modulus of different weathered states were similar. Therefore, the conversion equation of elastic modulus between intact rocks and seriously weathered rock masses was expressed as follows:
E d E min d E max d E min d = E E min E max E min
From Equation (13), the following can be derived:
E d = E E min E max E min ( E max d E min d ) + E min d
The same assumption for Poisson’s ratio is expressed as follows:
υ d = υ υ min υ max υ min ( υ max d υ min d ) + υ min d
E and v are the elastic modulus and Poisson’s ratio of the intact rock; Ed and vd are the elastic modulus and Poisson’s ratio of the seriously weathered rock mass; Emax, vmax, Emin, and vmin are the maximum and minimum elastic modulus and Poisson’s ratio of the intact rock; Edmax, vdmax, Edmin, and vdmin are the maximum and minimum elastic modulus and Poisson’s ratio of the seriously weathered rock mass, respectively. Edmax, vdmax, Edmin, and vdmin were measured by the uniaxial compression experiment. According to Equations (13)–(15), the strength parameters of the seriously weathered rock mass for different dip angles in the present study are shown in Figure 8.
Based on previous research [1,3,9,12,44], the anisotropy ratios of hydraulic conductivity were set as 0.002, 0.01, 0.02, 0.1, and 1, and the dip angles (anisotropy angle of hydraulic conductivity) were set as 0, 15, 30, 45, 60, 75, and 90°. The failure criterion of the rock layered slope simulation adopted the Mohr–Coulomb criterion. Correspondingly, the elastic modulus and Poisson’s ratio of the rock slope were also set according to Figure 8. The average values of the cohesion and friction angle for the whole rock slope were measured by the direct shear experiment on site, as detailed in Table 2.

3. Results

3.1. Effects of the Hydraulic Conductivity and Strength Anisotropy on the Deformation Characteristics

3.1.1. Analysis of the Horizontal Displacement of Monitoring Line

Based on the fluid–solid coupling theory, the layered rock slope produced the corresponding displacement under the action of the fluid. The seepage field exerted normal seepage pressure and tangential tension on the rock mass skeleton and affected the rock mass stress distribution. Under such stress conditions, the rock mass deformed accordingly. At the same time, the stress field affected the pores of the rock mass, and the hydraulic conductivity changed with the change of the pore size; thus, the stress field also affected the seepage field of the rock mass. According to the calculation conditions in Table 2, a total number of 70 numerical simulations was carried out, and 105 sections of horizontal displacement (HD) variation were obtained. The variations of the HD of layered rock slope with different β values and kr = 0.01 and kr = 0.1 are shown in Figure 9 to illustrate the influence of anisotropy direction β on deformation characteristics. Different kr values with β = 0, β = 30, and β = 90° are also displayed in Figure 10 to show the impacts of anisotropy ratio kr.
In Figure 9, the initial state is the distribution of the horizontal displacement (HD) of the rock slope before rainfall, and the remainder shows the distribution of the HD of day 155. The HDs of day 155 changed differently at different positions. At the top and middle of the slope, the HDs first increased and then decreased with the elevation. The maximum HDs at the top and middle of the slope were approximately 50 and 15 m from the surface, respectively. At the slope bottom, the HDs increased with elevation. The HDs at the bottom and middle of the slope were larger than the top of the slope. According to previous studies [3,9,30], the hydrodynamic pressure generated by rainfall had a large impact on the bottom and middle parts of the slope, and the landslide between the bottom and middle of the slope, which verified the results of the numerical simulation. Generally speaking, the hydraulic conductivity and strength anisotropy characteristics also had great impacts on the variation of HD. When the dip angle β reached 30°, the HD reached the maximum. When the HDs at different positions of the rock slope reached the minimum, the corresponding dip angles of the rock formation were different.
The trend of HD distribution in Figure 10 was similar in Figure 9. It is worth noting that when the dip angle β reached 30°, as the anisotropy ratio kr increased, the HD decreased. When the kr was smaller than 0.01, the maximum HD of the slope top appeared at 10 m from the surface. When the anisotropy angle was greater than 30° and kr = 1, the HD was the smallest. Therefore, considering the layered rock slope as a homogeneous medium may underestimate its horizontal displacement.
In Figure 11 and Figure 12, the initial state is the distribution of the pore water pressure (PWP) of the rock slope before rainfall and the remainder shows the distribution of the PWP of day 155. With the elevation, the initial PWP at the slope top and slope middle first decreased and then remained stable, and the initial PWP at the slope bottom kept decreasing. When the rock slope was regarded as an isotropic medium, the PWP distribution trend of day 155 was similar to the initial PWP distribution trend. The corresponding HDs meanwhile reached the minimum. When the rock slope was regarded as a heterogeneous medium, the PWP distribution trend of day 155 was different from the initial PWP distribution trend, and this difference was mainly manifested in the shallow layer. At the top and middle of the slope, as the kr decreased, the shallow PWP gradually increased, and at the slope bottom, as the kr decreased, the shallow PWP gradually decreased. The distribution of shallow HD was also similar to shallow PWP. At the top, middle, and bottom of the slope, as the β increased, the shallow PWP gradually increased. However, the shallow HD first increased, then decreased, and reached the maximum dip angle of 30°. Since the HD was not only affected by water pressure but also by its strength parameters (elastic modulus and Poisson’s ratio), the variations of HD were not synchronous with PWP.

3.1.2. Analysis of the Horizontal Displacement of Monitoring Point

The variations of horizontal displacement for the monitoring point of the layered rock slope during rainfall are shown in Figure 13. BPK2-1, BPK 2-2, and BPK 2-3 were the field monitoring data obtained by the displacement monitoring station. As a landslide occurred just near the field station during the monitoring activity, the deformation behavior of the rock slope could be detected [45]. In correspondence to the rainfall event of 4 July 2019 (day 91), the monitoring point of the rock slope experienced a sudden increase in horizontal displacement. Through numerical simulations, it was found that the HD of the monitoring point was underestimated when the rock slope was considered isotropic, compared to anisotropy during the rainfall period. Under varied kr and β, the HDs calculated by numerical simulations all increased abruptly on 4 July 2019, indicating that the rock slope model could accurately predict the occurrence time of the landslide. It is worth noting that when kr = 0.01 and β = 30°, the calculated HD was consistent with the field monitoring data.

3.1.3. Analysis of the Maximum Horizontal Displacement

In order to study the deformation at different slope positions under varied kr and β conditions more prominently, the maximum horizontal displacement (MHD) of each monitoring line was defined. Figure 14a–c shows that MHD at the different positions of the slope is greatly affected by kr and β. When β = 30°, the MHD reached the maximum, and the elastic modulus and Poisson’s ratio of the rock slope were the smallest in this case. When the anisotropy angle was close to 60°, and the MHD reached its minimum. When the rock slope was considered to be isotropic (kr = 1), the simulated MHD was almost negligible. Previous studies only took kr or β into consideration, but ignored their collective effect. Table 3 shows the MHD of the differences between only considering kr or β and considering both kr and β compared with the isotropy condition.
As can be seen from Table 3, the values of MHD for the rock slope are quite different under the conditions when only considering kr or β and when considering both kr and β. Therefore, the kr and β must be considered in the deformation analysis of rock slope, and MHD is more affected by β than kr.

3.2. Effects of the Hydraulic Conductivity and Strength Anisotropy on the Stability of the Rock Slope

3.2.1. Analysis of the Factor of Safety

The variations of the factor of safety (FS) for the rock slope under varied anisotropy ratios and dip angles are shown in Figure 15. When β = 30°, the overall FS of the rock slope was the smallest. During the period of heavy rainfall (from 21 June to 26 July 2019), the overall FSs obviously decreased and the FS of the rock slope decreased with the increase of β. Xia et al. [9] conducted relevant simulations and the results were similar to current research results. This is because a larger dip angle of the bedding plane would lead to a smaller hydraulic conductivity in the horizontal direction and a larger hydraulic conductivity in the vertical direction, which would promote the rainwater infiltrating into the deep part of the slope body [3,9,30], increase the height of the groundwater level and pore water pressure, and finally reduce the effective stress and shear strength. When the rainfall was heavy and the rainfall time was short, the FS increased sharply. This is because the shallow layer of the rock slope belongs to the unsaturated zone, and the hydraulic conductivity is small. A large amount of rainwater rushing into the rock slope for a short time caused blockage, which led to the accumulation of rainwater on the surface of the rock slope. The growth rate of pore water pressure inside the slope was less than the growth rate of the rainfall water pressure on the slope surface, which formed a reverse pressure on the rock slope. This would make it more difficult for the slope to transfer from an unsaturated state to a saturated state. The negative pressure zone gradually increased, leading to a corresponding increase in the FS of the rock slope.
It can be seen from Figure 16 that during the period of heavy rainfall (from 21 June to 26 July 2019), when β = 0°, the overall FSs of the rock slope decreased with the increase of kr. In the horizontal rock bedding plane, the hydraulic conductivity in the horizontal direction reached the maximum, but as kr increased, the hydraulic conductivity of the vertical direction increased, and the wet front of rain could reach the deep part of the slope. During the period of heavy rainfall, when β = 90°, the smaller the anisotropy ratio kr, the more notable the decrease of the FS was. In the vertical rock bedding plane, the hydraulic conductivity of the vertical direction reached the maximum, and as kr decreased, the hydraulic conductivity of the horizontal direction decreased, and the rainwater infiltration capacity was much greater than the horizontal drainage capacity. When the anisotropy angle β was close to 30°, with the change of the anisotropy ratio kr, the change law of FS was not notable. However, when rock slope was treated as an isotropic medium, the FS was overestimated.

3.2.2. Analysis of the Minimum Factor of Safety

In order to study the stability of rock slopes under varied kr and β conditions more prominently, the minimum factor of safety (MFS) and the factor of safety of landslide (FSL) on 4 July 2019 were defined. Variations in the factor of safety under varied kr and β are shown in Figure 17. When the kr is between 0.002 and 0.02, the MFS and FSL reached the minimum and maximum at β = 30 and β = 0°, respectively. It is worth noting that when the β was greater than 75°, the MFS reached the minimum, and the reason is discussed in Section 3.2.1. Table 4 shows the MFS and FSL of the differences between only considering kr or β and considering both kr and β compared with the isotropy condition.
As can be seen from Table 4, the values of MFS and FSL for rock slope are similar under the conditions when only considering β to considering both kr and β. Therefore, we can consider only the effect of β but ignore kr to simplify the calculation.

4. Discussion

The strength parameters of seriously weathered rock slopes like the elastic modulus and Poisson’s ratio were obtained by the conversion Equations (13)–(15). Through numerical simulations, the anisotropic rock slope model could predict the time of landslide occurrence. Figure 18 shows the cloud map of the horizontal displacement when the landslide occurred on 4 July 2019 (day 91). It could be found that when kr = 0.01 and β = 30°, the simulated landslide scope almost coincided with the field landslide scope. When the rock slope was regarded as a homogeneous medium (kr = 1), the distance between the simulated landslide scope and the field landslide scope was about 45 m, and the horizontal displacement was significantly small, which made it difficult to predict the occurrence of the landslide. According to the field survey, the field dip angle of the rock slope was 32°, as shown in Figure 19, which was similar to the dip angle of 30° in the simulation. Thus, in general, the landslide scope, horizontal displacement, dip angle, and landslide occurrence time of the field investigation were consistent with the results of the simulation. Therefore, the feasibility of the conversion Equations (11)–(14) of elastic modulus and Poisson’s ratio was also verified.

5. Conclusions

A layered rock slope in the Pulang Region, southwestern China, is selected as a case study. A numerical model of the landslide is constructed by considering different anisotropy ratios and anisotropy angles of hydraulic conductivity and strength. The research results were as follows:
(1)
The strength conversion equations of elastic modulus and Poisson’s ratio were feasible, and the rock slope model could accurately predict the occurrence time, horizontal displacement, and scope of the landslide;
(2)
The different anisotropy ratios and dip angles of the bedding plane were found to have major impacts on the deformation and stability of the layered rock slope;
(3)
The horizontal displacement (HD) and maximum horizontal displacement (MHD) were determined to characterize the deformation characteristics of the rock slope. Considering the layered rock slope as a homogeneous medium could underestimate its HD and MHD. When the dip angle was 30°, the MHD reached the maximum. When the anisotropy angle was close to 60°, and the MHD reached its minimum;
(4)
The factor of safety (FS), the minimum factor of safety (MFS), and the factor of safety of landslide (FSL) were determined to characterize the stability characteristics of the rock slope. When the dip angle was 30°, the FS, MFS, and FSL of the rock slope reached the minimum. However, when the rock slope was treated as an isotropic medium, the FS, MFS, and FSL were overestimated;
(5)
The changing law of the anisotropy ratio was not obvious, and it was difficult to verify by field data. This could be a focus of future research work. The obtained results are likely to provide a theoretical basis for the prediction and monitoring of layered rock landslide.

Author Contributions

C.X. wrote the manuscript, completed the experiment and method design, analyzed the data, and performed the experiments. G.L. and S.L. helped provide analysis of raw data. L.W. and J.D. helped polish the manuscript. Z.Z. provided research funding. L.Z. changed the format of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 41974148, and Hunan province 2019 special fund for safety production prevention and emergency, grant number 201907.

Acknowledgments

The authors gratefully acknowledge the finances support provided by The National Natural Science Fund (grant no. 41974148).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yeh, P.T.; Lee, K.Z.Z.; Chang, K.T. 3D Effects of permeability and strength anisotropy on the stability of weakly cemented rock slopes subjected to rainfall infiltration. Eng. Geol. 2020, 266. [Google Scholar] [CrossRef]
  2. Wei, X.S.; Fan, W.; Chai, X.Q.; Cao, Y.B.; Nan, Y.L. Field and numerical investigations on triggering mechanism in typical rainfall-induced shallow landslides: A case study in the Ren River catchment, China. Nat. Hazards 2020. [Google Scholar] [CrossRef]
  3. Yu, S.; Ren, X.; Zhang, J.; Wang, H.; Wang, J.; Zhu, W. Seepage, Deformation, and Stability Analysis of Sandy and Clay Slopes with Different Permeability Anisotropy Characteristics Affected by Reservoir Water Level Fluctuations. Water 2020, 12, 201. [Google Scholar] [CrossRef] [Green Version]
  4. Yeh, H.F.; Tsai, Y.J. Analyzing the Effect of Soil Hydraulic Conductivity Anisotropy on Slope Stability Using a Coupled Hydromechanical Framework. Water 2018, 10, 905. [Google Scholar] [CrossRef] [Green Version]
  5. Bai, D.; Tang, J.; Lu, G.; Zhu, Z.; Liu, T.; Fang, J. The design and application of landslide monitoring and early warning system based on microservice architecture. Geomat. Nat. Hazards Risk 2020, 11, 928–948. [Google Scholar] [CrossRef]
  6. Greco, R.; Marino, P.; Santonastaso, G.F.; Damiano, E. Interaction between Perched Epikarst Aquifer and Unsaturated Soil Cover in the Initiation of Shallow Landslides in Pyroclastic Soils. Water 2018, 10, 17. [Google Scholar] [CrossRef] [Green Version]
  7. Fan, X.Y.; Tang, J.J.; Tian, S.J.; Jiang, Y.J. Rainfall-induced rapid and long-runout catastrophic landslide on July 23, 2019 in Shuicheng, Guizhou, China. Landslides 2020, 17, 2161–2171. [Google Scholar] [CrossRef]
  8. Jordanova, G.; Gariano, S.L.; Melillo, M.; Peruccacci, S.; Brunetti, M.T.; Auflic, M.J. Determination of Empirical Rainfall Thresholds for Shallow Landslides in Slovenia Using an Automatic Tool. Water 2020, 12, 15. [Google Scholar] [CrossRef]
  9. Xia, C.; Lu, G.; Bai, D.; Zhu, Z.; Luo, S.; Zhang, G. Sensitivity Analyses of the Seepage and Stability of Layered Rock Slope Based on the Anisotropy of Hydraulic Conductivity: A Case Study in the Pulang Region of Southwestern China. Water 2020, 12, 2314. [Google Scholar] [CrossRef]
  10. Singh, P.K.; Singh, K.K.; Singh, T.N. Slope failure in stratified rocks: A case from NE Himalaya, India. Landslides 2016, 14, 1319–1331. [Google Scholar] [CrossRef]
  11. Dong, J.J.; Tu, C.H.; Lee, W.R.; Jheng, Y.J. Effects of hydraulic conductivity/strength anisotropy on the stability of stratified, poorly cemented rock slopes. Comput. Geotech. 2012, 40, 147–159. [Google Scholar] [CrossRef]
  12. Chen, Y.F.; Yu, H.; Ma, H.Z.; Li, X.; Hu, R.; Yang, Z. Inverse modeling of saturated-unsaturated flow in site-scale fractured rocks using the continuum approach: A case study at Baihetan dam site, Southwest China. J. Hydrol. 2020, 584, 124693. [Google Scholar] [CrossRef]
  13. Xu, Y.; Zhang, Z. Introduction to rock hydraulics. Geol. Hazards Environ. Prot. 1997, 6, 56. [Google Scholar]
  14. Collins, B.D.; Znidarcic, D. Stability Analyses of Rainfall Induced Landslides. J. Geotech. Geoenviron. Eng. 2004, 130, 362–372. [Google Scholar] [CrossRef]
  15. Barton, N.; Quadros, E. Anisotropy is Everywhere, to See, to Measure, and to Model. Rock Mech. Rock Eng. 2015, 48, 1323–1339. [Google Scholar] [CrossRef]
  16. Hu, S.C.; Tan, Y.L.; Zhou, H.; Guo, W.Y.; Hu, D.W.; Meng, F.Z.; Liu, Z.G. Impact of Bedding Planes on Mechanical Properties of Sandstone. Rock Mech. Rock Eng. 2017, 50, 2243–2251. [Google Scholar] [CrossRef]
  17. Yin, P.F.; Yang, S.Q. Experimental investigation of the strength and failure behavior of layered sandstone under uniaxial compression and Brazilian testing. Acta Geophys. 2018, 66, 585–605. [Google Scholar] [CrossRef]
  18. Jin, Z.; Li, W.; Jin, C.; Hambleton, J.; Cusatis, G. Anisotropic elastic, strength, and fracture properties of Marcellus shale. Int. J. Rock Mech. Min. Sci. 2018, 109, 124–137. [Google Scholar] [CrossRef] [Green Version]
  19. Fjær, E.; Nes, O.M. The impact of heterogeneity on the anisotropic strength of an outcrop shale. Rock Mech. Rock Eng. 2014, 47, 1603–1611. [Google Scholar] [CrossRef]
  20. Wu, Y.; Li, X.; He, J.; Zheng, B. Mechanical Properties of Longmaxi Black Organic-Rich Shale Samples from South China under Uniaxial and Triaxial Compression States. Energies 2016, 9, 1088. [Google Scholar] [CrossRef] [Green Version]
  21. Chen, Y.F.; Wei, K.; Liu, W.; Hu, S.H.; Hu, R.; Zhou, C.B. Experimental Characterization and Micromechanical Modelling of Anisotropic Slates. Rock Mech. Rock Eng. 2016, 49, 3541–3557. [Google Scholar] [CrossRef]
  22. Gholami, R.; Rasouli, V. Mechanical and elastic properties of transversely isotropic slate. Rock Mech. Rock Eng. 2014, 47, 1763–1773. [Google Scholar] [CrossRef]
  23. Xu, G.; He, C.; Su, A.; Chen, Z. Experimental investigation of the anisotropic mechanical behavior of phyllite under triaxial compression. Int. J. Rock Mech. Min. Sci. 2018, 104, 100–112. [Google Scholar] [CrossRef]
  24. Ramamurthy, T.; Rao, G.V.; Singh, J. Engineering behaviour of phyllites. Eng. Geol. 1993, 33, 209–225. [Google Scholar] [CrossRef]
  25. Ding, C.D.; Hu, D.W.; Zhou, H.; Lu, J.J.; Lv, T. Investigations of P-Wave velocity, mechanical behavior and thermal properties of anisotropic slate. Int. J. Rock Mech. Min. Sci. 2020, 127, 13. [Google Scholar] [CrossRef]
  26. Nasseri, M.H.B.; Rao, K.S.; Ramamurthy, T. Anisotropic strength and deformational behavior of Himalayan schists. Int. J. Rock Mech. Min. Sci. 2003, 40, 3–23. [Google Scholar] [CrossRef]
  27. Debecker, B.; Vervoort, A. Experimental observation of fracture patterns in layered slate. Int. J. Fract. 2009, 159, 51–62. [Google Scholar] [CrossRef]
  28. Geng, Z.; Chen, M.; Jin, Y.; Yang, S.; Yi, Z.C.; Fang, X.; Du, X.Y. Experimental study of brittleness anisotropy of shale in triaxial compression. J. Nat. Gas Sci. Eng. 2016, 36, 510–518. [Google Scholar] [CrossRef]
  29. Gonzaga, G.G.; Leite, M.H.; Corthésy, R. Determination of anisotropic deformability parameters from a single standard rock specimen. Int. J. Rock Mech. Min. Sci. 2008, 45, 1420–1438. [Google Scholar] [CrossRef]
  30. Yu, S.; Ren, X.; Zhang, J.; Wang, H.; Zhang, Z. Sensibility Analysis of the Hydraulic Conductivity Anisotropy on Seepage and Stability of Sandy and Clayey Slope. Water 2020, 12, 277. [Google Scholar] [CrossRef] [Green Version]
  31. Morgenstern, N.R.; Price, V.E. The Analysis of the Stability of General Slip Surfaces. Géotechnique 1965, 15, 79–93. [Google Scholar] [CrossRef]
  32. Oo, H.Z.; Ai, L.Z.; Qiu, Z. Numerical Analysis of River Bank Slope Stability during Rapid Drawdown of Water Level. Study Civ. Eng. Archit. 2013, 2, 98–103. [Google Scholar]
  33. Johansson, J.M.A.; Edeskr, T. Effects of external water-level fluctuations on slope stability. Electron. J. Geotech. Eng. 2014, 19, 2437–2463. [Google Scholar]
  34. Mandal, A.K.; Li, X.; Shrestha, R. Influence of Water Level Rise on the Bank of Reservoir on Slope Stability: A Case Study of Dagangshan Hydropower Project. Geotech. Geol. Eng. 2019, 37, 5187–5198. [Google Scholar] [CrossRef]
  35. Tang, D.; Li, D.; Zhou, C. Slope stability analysis considering the antecedent rainfall process. Rock Soil Mech. 2013, 34, 3239–3248. [Google Scholar]
  36. Fredlund, D.G.; Morgenstern, N.R.; Widger, R.A. The shear strength of unsaturated soils. Can. Geotech. J. 1978, 15, 313–321. [Google Scholar] [CrossRef]
  37. Genuchten, V.T.M. A Closed-form Equation for Predicting the Hydraulic Conductivity of Unsaturated Soils1. Soil Sci. Soc. Am. J. 1980, 44, 892–898. [Google Scholar] [CrossRef] [Green Version]
  38. Fredlund, D.G.; Xing, A. Equations for the soil-water characteristic curve. Can. Geotech. J. 1994, 31, 521–532. [Google Scholar] [CrossRef]
  39. Brooks, R.H.; Corey, A.T. Hydraulic Properties of Porous Media; Colorado State University, Hydrology and Water Resources Program: Fort Collins, CO, USA, 1964; pp. 1–27. [Google Scholar]
  40. Mualem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef] [Green Version]
  41. Huang, C.; Zuo, S.; Wang, S.; Qu, C.; Zhao, Y. Laboratory uniaxial compression test analysis of layered anisotropic rock mass. J. Yangtze River Sci. Res. Inst. 2016, 5, 58–62. [Google Scholar]
  42. Lemaitre, J. How to use damage mechanics. Nucl. Eng. Des. 1984, 80, 233–245. [Google Scholar] [CrossRef]
  43. Guo, Q.G.; Che, A. Research on relationship of rock mass integrity index and rock mass elastic modulus. Chin. J. Rock Mech. Eng. 2011, 9, 3914–3919. [Google Scholar]
  44. Jiang, Z.; Xiong, X.; Zeng, L. Unsaturated seepage analysis of slope under rainfall condition based on FLAC3D. Geotech. Mech. 2014, 35, 855–861. [Google Scholar]
  45. Montrasio, L.; Valentino, R.; Meisina, C. Soil Saturation and Stability Analysis of a Test Site Slope Using the Shallow Landslide Instability Prediction (SLIP) Model. Geotech. Geol. Eng. 2018. [Google Scholar] [CrossRef]
Figure 1. Morgenstern–Price (M–P) method calculation principle and model. (modified after Morgenstern et al. [31]).
Figure 1. Morgenstern–Price (M–P) method calculation principle and model. (modified after Morgenstern et al. [31]).
Water 12 03056 g001
Figure 2. Geographic location of the study area: (a) Study area location; (b) contour map of rock slope; (c) topography of rock slope; (d) right elevation of rock slope; (e) left elevation of rock slope.
Figure 2. Geographic location of the study area: (a) Study area location; (b) contour map of rock slope; (c) topography of rock slope; (d) right elevation of rock slope; (e) left elevation of rock slope.
Water 12 03056 g002
Figure 3. Grid model: (a) Rock landslides of rock slopes; (b) displacement monitoring station (BPK2-1, BPK2-2, and BPK2-3).
Figure 3. Grid model: (a) Rock landslides of rock slopes; (b) displacement monitoring station (BPK2-1, BPK2-2, and BPK2-3).
Water 12 03056 g003
Figure 4. Daily rainfall from 1 March to 2 August 2019 (155 days in total).
Figure 4. Daily rainfall from 1 March to 2 August 2019 (155 days in total).
Water 12 03056 g004
Figure 5. Initial pore pressure distribution.
Figure 5. Initial pore pressure distribution.
Water 12 03056 g005
Figure 6. (a) Water retention curve; (b) hydraulic conductivity curve.
Figure 6. (a) Water retention curve; (b) hydraulic conductivity curve.
Water 12 03056 g006
Figure 7. Rock specimens with different dip angles, β (modified after Huang et al. [41]).
Figure 7. Rock specimens with different dip angles, β (modified after Huang et al. [41]).
Water 12 03056 g007
Figure 8. The anisotropic strength parameters for different dip angles: (a) Elastic modulus; (b) Poisson’s ratio.
Figure 8. The anisotropic strength parameters for different dip angles: (a) Elastic modulus; (b) Poisson’s ratio.
Water 12 03056 g008aWater 12 03056 g008b
Figure 9. Horizontal displacement levels at different positions on the layered rock slope under different β values with kr = 0.01 and kr = 0.1: (a) Top of the slope with kr = 0.01; (b) middle of the slope with kr = 0.01; (c) middle of the slope with kr = 0.01; (d) top of the slope with kr = 0.1; (e) middle of the slope with kr = 0.1; (f) bottom of the slope with kr = 0.1.
Figure 9. Horizontal displacement levels at different positions on the layered rock slope under different β values with kr = 0.01 and kr = 0.1: (a) Top of the slope with kr = 0.01; (b) middle of the slope with kr = 0.01; (c) middle of the slope with kr = 0.01; (d) top of the slope with kr = 0.1; (e) middle of the slope with kr = 0.1; (f) bottom of the slope with kr = 0.1.
Water 12 03056 g009aWater 12 03056 g009b
Figure 10. Horizontal displacement levels at different positions on the layered rock slope under different kr values with β = 0, β = 45, and β = 90°: (a) Top of the slope with β = 0°; (b) middle of the slope with β = 0°; (c) bottom of the slope with β = 0°; (d) top of the slope with β = 30°; (e) middle of the slope with β = 30°; (f) bottom of the slope with β = 30°; (g) top of the slope with β = 90°; (h) middle of the slope with β = 90°; (i) bottom of the slope with β = 90°.
Figure 10. Horizontal displacement levels at different positions on the layered rock slope under different kr values with β = 0, β = 45, and β = 90°: (a) Top of the slope with β = 0°; (b) middle of the slope with β = 0°; (c) bottom of the slope with β = 0°; (d) top of the slope with β = 30°; (e) middle of the slope with β = 30°; (f) bottom of the slope with β = 30°; (g) top of the slope with β = 90°; (h) middle of the slope with β = 90°; (i) bottom of the slope with β = 90°.
Water 12 03056 g010aWater 12 03056 g010b
Figure 11. Pore water pressure levels at different positions on the layered rock slope under different β values with kr = 0.01 and kr = 0.1: (a) Top of the slope with kr = 0.01; (b) middle of the slope with kr = 0.01; (c) middle of the slope with kr = 0.01; (d) top of the slope with kr = 0.1; (e) middle of the slope with kr = 0.1; (f) bottom of the slope with kr = 0.1.
Figure 11. Pore water pressure levels at different positions on the layered rock slope under different β values with kr = 0.01 and kr = 0.1: (a) Top of the slope with kr = 0.01; (b) middle of the slope with kr = 0.01; (c) middle of the slope with kr = 0.01; (d) top of the slope with kr = 0.1; (e) middle of the slope with kr = 0.1; (f) bottom of the slope with kr = 0.1.
Water 12 03056 g011
Figure 12. Pore water pressure levels at different positions on the layered rock slope under different kr values with β = 0, β = 30, and β = 90°: (a) Top of the slope with β = 0°; (b) middle of the slope with β = 0°; (c) bottom of the slope with β = 0°; (d) top of the slope with β = 30°; (e) middle of the slope with β = 30°; (f) bottom of the slope with β = 30°; (g) top of the slope with β = 90°; (h) middle of the slope with β = 90°; (i) bottom of the slope with β = 90°.
Figure 12. Pore water pressure levels at different positions on the layered rock slope under different kr values with β = 0, β = 30, and β = 90°: (a) Top of the slope with β = 0°; (b) middle of the slope with β = 0°; (c) bottom of the slope with β = 0°; (d) top of the slope with β = 30°; (e) middle of the slope with β = 30°; (f) bottom of the slope with β = 30°; (g) top of the slope with β = 90°; (h) middle of the slope with β = 90°; (i) bottom of the slope with β = 90°.
Water 12 03056 g012aWater 12 03056 g012b
Figure 13. Variations in the horizon displacement of numerical simulation and field monitoring at different monitoring points during rainfall: (a) BPK2-1; (b) BPK2-2; (c) BPK2-3.
Figure 13. Variations in the horizon displacement of numerical simulation and field monitoring at different monitoring points during rainfall: (a) BPK2-1; (b) BPK2-2; (c) BPK2-3.
Water 12 03056 g013aWater 12 03056 g013b
Figure 14. Variations in the maximum horizon displacement (MHD) at the different positions of the slope: (a) MHD of the top of the slope; (b) MHD of the middle of the slope; (c) MHD of the bottom of the slope.
Figure 14. Variations in the maximum horizon displacement (MHD) at the different positions of the slope: (a) MHD of the top of the slope; (b) MHD of the middle of the slope; (c) MHD of the bottom of the slope.
Water 12 03056 g014
Figure 15. Variations in the factor of safety during rainfall under different β values with kr = 0.01 and kr = 0.1: (a) kr = 0.01; (b) kr = 0.1.
Figure 15. Variations in the factor of safety during rainfall under different β values with kr = 0.01 and kr = 0.1: (a) kr = 0.01; (b) kr = 0.1.
Water 12 03056 g015
Figure 16. Variations in the safety factors during rainfall under different kr values with β = 0, β = 30, and β = 90°: (a) β = 0; (b) β = 30; (c) β = 90°.
Figure 16. Variations in the safety factors during rainfall under different kr values with β = 0, β = 30, and β = 90°: (a) β = 0; (b) β = 30; (c) β = 90°.
Water 12 03056 g016
Figure 17. Variations in the factor of safety: (a) Factor of safety of landslide (FSL); (b) minimum factor of safety (MFS).
Figure 17. Variations in the factor of safety: (a) Factor of safety of landslide (FSL); (b) minimum factor of safety (MFS).
Water 12 03056 g017
Figure 18. Cloud map of horizontal displacement of rock slope when landslide occurs (day 91): (a) Isotropic (kr = 1); (b) kr = 0.01, β = 30°.
Figure 18. Cloud map of horizontal displacement of rock slope when landslide occurs (day 91): (a) Isotropic (kr = 1); (b) kr = 0.01, β = 30°.
Water 12 03056 g018
Figure 19. Stereographic projections of the rock slope.
Figure 19. Stereographic projections of the rock slope.
Water 12 03056 g019
Table 1. Parameters of the water retention curve.
Table 1. Parameters of the water retention curve.
LayerMaterialsFitting ParametersHydraulic Conduction Coefficient
A (kPa)mnθsθrk (m/s)
IStrongly weathered carbonaceous slate100.331.50.2420.0018.08 × 10−5
Table 2. Different calculation conditions.
Table 2. Different calculation conditions.
Rock TypesAnisotropy Ratio
kr = ky/kx
Anisotropic Angle/ Dip Angle β (°)Elastic Modulus Ed (GPa)Poisson Ratio vdUnit Weight (kN/m3)Cohesion (kPa)Friction Angle (°)
seriously weathered carbonaceous slate [ 0.002 0.01 0.02 0.1 1 ] [ 0 15 30 45 60 75 90 ] [ 2.902 3.014 2.645 2.728 2.953 4.428 5.561 ] [ 0.3569 0.3546 0.3500 0.3731 0.3638 0.3800 0.3615 ] 22.493.633.3
Table 3. Differences between only considering kr or β and considering both kr and β.
Table 3. Differences between only considering kr or β and considering both kr and β.
MHDOnly Considering krOnly Considering βConsidering Both kr and β
Slope top23%118%2.25 × 105%
Slope middle86%127%3.38 × 105%
Slope bottom81%123%1.17 × 105%
Table 4. Differences between only considering kr or β and considering both kr and β.
Table 4. Differences between only considering kr or β and considering both kr and β.
MHDOnly Considering krOnly Considering βConsidering Both kr and β
FSL7%16%16.4%
MFS6%17%17.6%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xia, C.; Lu, G.; Zhu, Z.; Wu, L.; Zhang, L.; Luo, S.; Dong, J. Deformation and Stability Characteristics of Layered Rock Slope Affected by Rainfall Based on Anisotropy of Strength and Hydraulic Conductivity. Water 2020, 12, 3056. https://doi.org/10.3390/w12113056

AMA Style

Xia C, Lu G, Zhu Z, Wu L, Zhang L, Luo S, Dong J. Deformation and Stability Characteristics of Layered Rock Slope Affected by Rainfall Based on Anisotropy of Strength and Hydraulic Conductivity. Water. 2020; 12(11):3056. https://doi.org/10.3390/w12113056

Chicago/Turabian Style

Xia, Chengzhi, Guangyin Lu, Ziqiang Zhu, Lianrong Wu, Liang Zhang, Shuai Luo, and Jie Dong. 2020. "Deformation and Stability Characteristics of Layered Rock Slope Affected by Rainfall Based on Anisotropy of Strength and Hydraulic Conductivity" Water 12, no. 11: 3056. https://doi.org/10.3390/w12113056

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