Abstract

During the reservoir operation process, the long-term security and stability of the bank slope is affected by dynamic response characteristics of its seismic action directly. Aimed at the typical bank slope existing in the actual reservoir environment, an experiment considering reservoir water level fluctuation and soaking-air-drying cyclic water-rock interaction has been designed and conducted while the cyclic loading test was performed in different water-rock cycles. Research results indicate the following: Firstly, in the process of water-rock interaction, the dynamic characteristics of sandstone show evident degradation trend, with the increase of the damping ratio and Poisson’s ratio and decrease of dynamic elastic modulus, and the former six water-rock cycle degradation effects are particularly obvious. Secondly, the numerical analog computation analysis of dynamic response in typical bank slope shows that as the water-rock interaction period is increased, the dynamic response of the slope hydro-fluctuation belt zone increases gradually, while the other parts weaken. Thirdly, under the long-term water-rock interaction process, the hydro-fluctuation belt zone gradually becomes a “soft layer” which is sensitive to the earthquake effect and dynamic response, resulting in a direct influence on long-term seismic performance of the bank slope. Therefore, it is necessary to make better protection for the bank slope hydro-fluctuation belt zone.

1. Introduction

In the process of hydraulic and hydroelectric engineering construction, reservoir impoundment will change geological and mechanical conditions naturally formed over the years [1]. This has triggered series of new geological disasters, one of which is inevitable reservoir earthquake. The reservoir earthquakes can be classified as karst type, ore collapse type, and superficial stress adjustment type [2]. And, some of the reservoir earthquake mechanisms are related to their hydrogeological structure and the nearby fault coulomb stress [3], while the others correspond with karst effect [4]. Recent research results show that small earthquakes in the Three Gorges Reservoir Area appeared to be linear or clumpy distribution [5].

For example, a 6.1 magnitude earthquake that occurred near Xinfengjiang reservoir in 1962 caused significant damage to its dam. Since the impoundment of the Three Gorges Reservoir on June 2003, a series of reservoir earthquakes ensued along the river section in Zigui and Badong County, which has reached over ten thousand times so far. Most of them are small earthquakes below the level of 3-4, and the magnitude of 4.0 or above has occurred six times. On December 16, 2013, a 5.1 magnitude earthquake occurred in Badong County of the Three Gorges reservoir area, with a depth of 5.0 km and about 100 kilometers from the Three Gorges Dam. It was the largest earthquake since the impoundment of the Three Gorges Reservoir. Analysis indicates that its recorded seismic waveform has long period with more low-frequency components and less high-frequency ones. Except for deformation and damage of hydraulic structures, like the rock collapse in the higher water pressure [6, 7], reservoir-induced earthquake also generates significant impacts on the deformation and safety of the bank slopes. Thus, the dynamic response to the seismic action of the bank slope plays an important role in the long-term deformation and stability analysis. In the dynamic response analysis of the bank slope, it is necessary to consider the influence of reservoir water on the dynamic characteristics.

Previous studies on water-rock interaction show that the saturated condition or dry-wet circulation influenced obviously the static characteristics of rock mass. Liu et al. [8], Fu et al. [9], Yao et al. [10], Jiang et al. [11], Hale and Shakoor [12], and Zhu et al. [13] successively carried out a series of dry-wet cycling tests on sandstone and shale, whose results indicated that distinct degradation generated on the rock’s strength parameters during the dry-wet circulation. Jeng et al. [14] and Lin et al. [15] studied the microscopic mechanism of sandstone degradation under the dry-wet circulation. They found that the sandstone strength decreased obviously under the dry-wet circulation and the porosity increased in the nonlinear form. Combined with self-designed rock soak device YRK-1 to simulate repeated fluctuation of water level from 145 m to 175 m, Deng et al. [1619] carried out the water-rock test considering changes in water pressure and soaking dry cycle, which indicated that water-rock damage is greatly affected by the varied soaking water pressure.

Research results above indicate that the water-rock interaction of the bank slope has obvious weakening effect on its static characteristics, and this effect shows distinct accumulation and timeliness. So, does the rock dynamic response under water-rock interaction fit with the similar variation law?

From the current research of the water-rock interaction, most of the works mainly focus on the rock static characteristics degradation law, with the rare study of its dynamic characteristics. The rock dynamic parameter deterioration influenced by water-rock interaction is also less considered in the long-term evaluation of the bank slope stability.

In this paper, the typical bank slope was selected as the object from the Three Gorges Reservoir Area. Considering the repeated fluctuation of the reservoir water level and the process of soaking-air-drying circulation, the water-rock interaction test is designed and carried out. The cyclic loading and unloading tests were performed in different water-rock interaction cycles, for the analysis of sandstone dynamic characteristics degradation rule. By the numerical simulation calculation of FLAC3D software on the dynamic module, the research on the slope dynamic response under long-term water-rock interaction was carried out in this paper.

2. Experimental Programs

2.1. Sample Preparation

Originating from the bank slope of Baishuihe in the Three Gorges Reservoir Area, known as weak weathering sandstone with medium-grained quartz inside, a sample was selected as the rock sample in this paper. According to the standard requirements [20] in the Specifications in the Testing Method of Engineering Rock Mass(GB/T 50266-2013), the standard rock samples with a diameter of 50 mm and a height of 100 mm are prepared. After the longitudinal wave velocity and mass measurement, the rock samples with relative concentration of wave velocity and density are selected as test samples. Typical test rock samples are shown in Figure 1.

2.2. The Water-Rock Interaction Test Scheme

According to the operation requirements of the Three Gorges Reservoir, the water level of the reservoir fluctuates between 145 m and 175 m each year after its completion, which generates the hydro-fluctuation belt with a height of 30 m between the two sides of the 600 km long channel. It is a typical zone affected greatly by water-rock interaction. Under substantial fluctuation of the reservoir water, the reservoir water supplies groundwater occasionally and sometimes supplies the others in the reverse order, which generate alternate state of soaking and dewatering in the rock mass of the hydro-fluctuation belt. This alternate effect is a kind of “fatigue” [21, 22] for rock mass resulting in its property degradation; each of the processes may not be significant enough, but cumulative damage effect will develop through numerous repeated processes of soaking and dewatering. In other words, the hydro-fluctuation belt is a sensitive deformation zone of the bank slope. Therefore, mainly considering the water-rock interaction of the hydro-fluctuation belt, it is crucial to simulate the change process of the soaking water pressure, including rising process, consistent process, declining process, and the natural air drying process after that. Hence, as is shown in Figure 2, the self-developed YRK-2 rock immersion-air-dry circulation test instrument was used to simulate the water pressure fluctuation and soak-air-dry cycle process. This instrument realizes the water-rock interaction process such as pressure soaking and temperature controlled air-drying, with automatic collection of parameters like water pressure and temperature.

According to the previous experimental experience, the experiment scheme of soaking-air-dry circulation to simulate water-rock interaction process was designed, and the process is shown in Figure 3. A single water-rock cycle lasts for 40 days and is divided into two stages: the first stage is the soaking period—firstly, the samples are soaked in a container filled with water for 10 days; during the time, water pressure increases to 0.3 MPa (to simulate the water level rise from 145 m to 175 m), then the water pressure remains unchanged for 10 days (to simulate the constant water level of 175 m), and in the end, the uniform water pressure reduces uniformly to 0 within the last 10 days (to simulate the water level decrease from 175 m to 145 m); the second stage is the air-drying stage, during which rock samples are soaked in a container with the constant temperature of 35°C (simulation process of dry during low initial water level) and then air dried for 10 days. Then, the soak-air-dry procedure is repeated, and the total design cycle is 10 times. For the study of the dynamic characteristics of the sample degradation rule at different water-rock cycles, cyclic loading and unloading experiments and the uniaxial compression test were performed on a set of samples selected separately at the end of 1st, 2nd, 4th, 6th, 8th, and 10th cycle of soaking.

2.3. The Cyclic Loading and Unloading Test Scheme

As for the seismic simulation test, cyclic loading and unloading tests are the most common ones in the laboratory. Previous research shows that the accumulation of irreversible deformation of rock, growth trend, and total fatigue under cyclic loading is directly related to fatigue damage [23]. According to the standard requirements [24, 25], saturated rock uniaxial compressive strength is about 50 MPa, and the constant loop unloading test was carried out considering the frequent emergence of reservoir earthquake with low frequency and magnitude [26]. The lower and upper limits of cyclic load stress are, respectively, defined as 10 MPa and 20 MPa, which refer to 20% and 50% of rock uniaxial compressive strength [27]. The sine wave repeated 30 times was used for cyclic loading and unloading tests, with a frequency of 0.1 Hz. Typical cyclic loading and unloading stress-strain curve is shown in Figure 4.

As is shown in Figure 5, the RMT-150CROCk Mass Zesting system is used to carry out the cyclic test. It can switch from different load models including displacement, time, and their combination, which is suitable for uniaxial or triaxial compression unloading circulation tests, and so on.

3. Analysis of the Sandstone Dynamic Characteristics Deterioration under Water-Rock Interaction

3.1. Calculation Principle of Rock Dynamic Parameters under Cyclic Loading and Unloading Test

In the loading and unloading circulation process, the stress-strain curve often developed hysteresis loop due to the nonideal elastic medium of the rock, as is shown in Figure 6. Damping ratio and dynamic elastic modulus are defined [28] in the following equations:

3.2. Analysis of the Sandstone Dynamic Characteristics Deterioration Rule under the Water-Rock Interaction

To facilitate the comparison analysis of sample’s characteristics degradation effect in different water-rock circulations, the 16th stress-strain hysteresis loop was voted as calculation for the related rock dynamic parameters. The variation curves of damping ratio, dynamic modulus of elasticity, and passion ratio are as given in Figures 79, respectively.

As is depicted in Figure 7, the damping ratio of the rock sample increased gradually as water-rock interaction process increased. After the 1st, 2nd, 4th, 6th, 8th, and 10th cycles of immersion air-dry, damping ratio increased to 7.71%, 14.49%, 20.55%, 22.81%, 24.35%, and 25.45%, respectively. The curve of the damping ratio grew fast during the first six immersion air-dry cycles and then gradually tended to be slow, which shares the similar law with the static strength deterioration and deformation tendencies of sandstone [1619].

As is depicted in Figure 8, the dynamic elasticity modulus of the rock sample decreased gradually as water-rock interaction process increased. After the 1st, 2nd, 4th, 6th, 8th, and 10th cycles of immersion air-dry, the dynamic elasticity modulus decreased to 8.23%, 15.80%, 22.10%, 26.52%, 28.35%, and 29.42%, respectively. The curve of the dynamic elasticity modulus declined fast during the first six immersion air-dry cycles and then gradually tended to be slow.

As is depicted in Figure 9, Poisson’s ratio of the rock sample increased gradually as water-rock interaction process increased. After the 1st, 2nd, 4th, 6th, 8th, and 10th cycles of immersion air-dry, Poisson’s ratio increased to 4.52%, 8.14%, 13.12%, 16.29%, 17.65%, and 18.55%, respectively. The curve of Poisson’s ratio grew fast during the first six immersion air-dry cycles and then gradually tended to be slow.

According to energy theory, the rock dynamic elastic modulus is the embodiment of dynamic elastic parameters under dynamic load, and its value reflects the quality of rock’s elastic bearing performance; the rock damping ratio represents the ratio of the total energy consumed by the rock to its elastic strain energy during the single loading period. With the increase of water-rock interaction cycle, the dynamic elastic modulus gradually reduces while damping ratio increases, which indicates that under the water-rock interaction, sandstone gradually becomes soft, with internal pores and fractures being developed fast at the same time, and a single load produced more rock-dissipated energy than before. The related microstructure deterioration mechanism has been elaborated in the literature [1619].

4. Analysis of the Dynamic Response of the Bank Slope under the Water-Rock Interaction

4.1. Slope Model

In order to analyze the dynamic response change law of the rock mass slope under the long-term water-rock interaction, a numerical simulation analysis of the typical rock slope is established, as shown in Figure 10. The area with the length, width, and height of 1000 m, 400 m, and 10 m, respectively, is taken for the calculation of range of slope, and the slope trailing edge length is 400 m with slope angle of 45°. The model was totally composed of 4493 nodes and 12353 units. Based on the previous study by Liu Xinrong et al [29], with the consideration of reservoir water level of 145 m and 175 m, it is of great importance to analyze the dynamic response influenced by dynamic parameters degradation in the hydro-fluctuation belt zone of the bank slope. The rock mass under the reservoir was regarded as a saturated material in the numerical simulation analysis.

According to the test results above, the dynamic parameters of slope rock mass were shown in Table 1. In the calculation process, the corresponding values in Table 1 were used as the dynamic parameters in the hydro-fluctuation belt at different water-rock cycles, but other areas of the slope rock mass use the parameters at the initial state for all the cycles.

4.2. Boundary Conditions and Record Point Settings

FLAC3D uses static boundary and free field boundary in dynamic analysis. In this paper, the rock slope was selected as the study object. Due to its large value of bedrock modulus, the bottom of the slope is set as a rigid foundation for simulation processing, and free field boundary around the slope model, to make side border of the main body grid couple with the free field through the damper. In the simulation process, the rock mass in the model is regarded as an elastoplastic material with the use of local damping and the Mohr–Coulomb strength criterion. The calculation model and boundary condition settings are shown in Figure 11.

There are 12 record points set on the slope, respectively, tracking the acceleration time of each point, and their elevation and point numbers are shown in Table 2. The P01 point is at top of the slope, and P07, P08, and P09 are from the hydro-fluctuation belt area.

5. The Seismic Dynamic Response of the Slope under the Parameter Deterioration Effect

On the basis of previous research studies, the earthquake max magnitudes of the Three Gorges Reservoir induced by water storage ranged from 5.5 to 6; thus, the peak value of the horizontal acceleration is set as 0.05 g in this paper. After filtering and baseline correction through seismo signal, the famous El Centro seismic wave was selected as the input seismic wave, with its 30 s acceleration time input to FLAC3D, as is shown in Figure 12.

And the earthquake waves were put in from the bottom of the slope model, of which horizontal waves take the prime. From recent research reports, it can be found that the actual dynamic response in bank slope fits with the results of test and numerical simulation process sometimes, and the model set and parameters standard for evaluation are of great importance [30, 31].

The acceleration response and its distribution law are the basic information to evaluate the slope dynamic response. Therefore, the seismic dynamic response of the slope under the degradation condition is analyzed mainly from its acceleration response.

5.1. Dynamic Response Rules of the Bank Slope

In order to describe the slope acceleration response rule under the earthquake action, the ratio of dynamic response peak acceleration of any point in the slope body and that of point in the slope toe was defined as peak value of acceleration amplification factor (PGA). The variation rule of the horizontal direction amplification coefficient (HPGA) of each record point on the slope surface in different water-rock cycles is shown in Figure 13.

As seen in Figure 13, without considering the water-rock interaction, the HPGA of each record point of the slope increases gradually as its height rises and reaches the peak at the top of the slope. And the initial line (T0) shows the same rule as the sandstone static parameters in the slope surface. However, its regular curve is increasing in the form of fold line, like the line of T4, which is consistent with the conclusion summarized by Zhixin Yan et al. [32].

And the variation trend of the HPGA of the record point at the top of the slope under the water-rock interaction cycle is shown in Figure 14, with that of the hydro-fluctuation belt shown in Figure 15.

When the water-rock interaction is considered, the HPGA of the record point in the hydro-fluctuation belt area increases rapidly, while that of other places in the slope decreases slowly. When the water level is 145 m, for example, after the 1st cycle of water-rock interaction, the HPGA of the record point P07 increases from 1.23 to 1.37, while that of the record point P01 goes down from 1.39 to 1.28. It can be found that the HPGA of the hydro-fluctuation belt shows greater magnification than that of the slope top after the 10th cycle of water-rock interaction and the HPGA of record point P07 increases to 1.60, while that of record point P01 goes down to 1.14.

With the increase of the water-rock interaction period, the HGPA of the record points in the hydro-fluctuation belt area increases gradually with the total increase of 29.03% to 34.17% at the end of the 10th water-rock interaction cycle; the HGPA of the record points in the other area decreases gradually with the total decrease of 17.93% to 18.65% at the same time. During the first six water-rock interaction processes, the HGPA of each record point shows obvious variation trend depicted in the figures above.

With the reservoir water level of 145 m and 175 m, dynamic response of bank slope share the similar law under earthquake action and so does the HPGA change law of each record point in the process of water-rock interaction. In comparison, with the water level of 145 m, the HPGA of the record point on the slope is relatively less than that of 175 m, and the difference between the two is about 0.05.

5.2. Dynamic Response Mechanism Analysis of the Slope under Water-Rock Interaction

As seen from the results above, the main influence of water-rock interaction in the bank slope focuses on the enhanced dynamic response in the hydro-fluctuation belt and the relatively weaken response of its other parts. There are two main reasons as follows:

First of all, without considering the water-rock interaction, the hydro-fluctuation belt shares the same rock properties with the other parts of the bank slope rock mass; thus, the HPGA of the slope surface presents the similar rule to general seismic response of the slope, appeared as the top acceleration amplification effect. However, during the process of water-rock interaction, the dynamic characteristics of rock mass are gradually deteriorating under long-term immersion-air-drying alternation, with the increase of damping ratio and the damping coefficient and decrease of the dynamic elasticity modulus. Similar to the near-surface low-velocity layer of slope [33], the hydro-fluctuation belt has a greater magnifying effect on the dynamic response of seismic waves than the normal slope surface, resulting in the more obvious seismic amplification effect.

Secondly, the slope shows vertical exaggeration and near-surface magnifying effect on the input seismic wave. Based on the previous studies of rock mass under higher water compression, the rock of the bank slope appeared with varied degrees of damage and destruction. After the long-term rock mass accumulation damage [34], the bank slope body dynamical response gradually strengthened, leading to slope collapse occasionally. As for the reason, it mainly concerns about the energy storage and release of rock mass. Due to the strong amplification effect of relatively low-intensity rock mass on seismic wave, according to the law of energy conservation [35, 36], the strong dynamic response in the hydro-fluctuation belt released much energy, leading to the relatively less energy released by the other parts, resulting in weaken dynamic response corresponding. With the dynamic parameters continued deterioration of rock mass as the water-rock interaction cycle increased, the increasing consumption and release of energy in the hydro-fluctuation belt gradually weakened the acceleration response of the other parts of the slope.

6. Conclusions

Based on the results and analysis obtained above, the following conclusions can be made:

In the process of water-rock interaction, the damping ratio and Poisson's ratio of rock sample gradually increased, and the dynamic elastic modulus decreased. At the end of the 6th water-rock interaction, the variation trend of each dynamic parameter is tended to slow down, and the variable quantity of the damping ratio, dynamic elastic modulus, and Poisson ratio, under the corresponding single water-rock interaction cycle, is gradually reduced.

Without considering the water-rock interaction, the HPGA of each record point of the slope increases gradually as its height rises. When the water-rock interaction is considered, the HPGA of the record point in the hydro-fluctuation belt area increases rapidly, while that of other places in the slope decreases on the whole.

With the increase of the water-rock interaction period, the HGPA of the record points in the hydro-fluctuation belt area increases gradually, while that of the other area decreases. During the first six water-rock interaction processes, the HGPA of each record point shows obvious variation trend, and then tends to slow down in the later period.

During the process of long-term water-rock interaction, the rock mass in the hydro-fluctuation belt area of the bank slope generates accumulative damage. With the continuous degradation of its static and dynamic parameters, seismic dynamic response gradually strengthened with the increasing energy consumption, resulting in the weakened dynamic responses of the other parts of the slope.

Therefore, the hydro-fluctuation belt is the sensitive zone of the seismic dynamic response of the bank slope, which is closely related to the long-term seismic performance of the bank slope, and its protection should be strengthened in the slope prevention.

Nomenclature

:Maximum dynamic stress of dynamic stress-dynamic strain hysteresis curve
:Minimum dynamic stress of dynamic stress-dynamic strain hysteresis curve
:Maximum dynamic strain of dynamic stress-dynamic strain hysteresis curve
:Minimum dynamic strain of dynamic stress-dynamic strain hysteresis curve
:The area of ABCD in the hysteresis loop
:The area of the triangle AEF
ς:Damping ratio of rock sample
Ed:Dynamic elastic modulus of rock sample
Poisson’s ratio of rock sample
PGA:Peak value of acceleration amplification factor
HPGA:Horizontal peak value of acceleration amplification factor.

Data Availability

The data in this article are obtained from original laboratory tests with all the authors’ help. Anyone who wants to use the data should acknowledge the provenance of the article and apply for permission from us if convenient.

Conflicts of Interest

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

Acknowledgments

The research was funded by National Natural Science Foundation of China (Grant no. 51439003), Natural Science Foundation of Hubei Province (Grant no. 2015CFA140), and National Natural Science Foundation of China (Grant no. 51679127), and it was also sponsored by Research Fund for Excellent Dissertation of China Three Gorges University(Grant no. 2018SSPY023).