Assessment Slope Stability Based on Deformation of Rock Joints and Soil with Simulation Method ()
1. Introduction
Today, fast developments in civil projects due to the fast population growth as well as providing the human needs, development of cities, and villages, and their development toward the slope areas, placement of transportation roots such as roads and railroads, and focused structures such as cities and villages, occurrence of landslides and consequent hazards are highly considered. Unfortunately, due to lack of attention to slope stability in the civil projects, every year especially after heavy precipitations and earthquakes, catastrophic happenings related to slope instability, in particular rock fall are observed in Iran. Since Iran has a lot of mountainous areas, the road construction projects face with passing through the mountains, which necessitates the stability studies in terms of natural and man-made slopes to ensure successful achievement and safety of the projects. A voluminous number of studies have been done related to slope stability including rock block and soil analysis [1] - [8] . Landslide occurs in mountainous areas of Iran every year, which results in significant damages. There are no accurate statistical data and related damages of landslides in Iran. Although no damage surveying has been done yet, the hazards both direct and indirect ones, occur highly every year, including destruction of residential areas, roads, farmlands, gardens, jungles, soil erosion and production of sediments.
2. Study Area
The study area is situated in Baluchestan province, which is located in the south east of Iran, between Khash and Iran shah cities: N60047, E27˚36 (area A) and N27˚26, E27, 26 (Area B), as shown in Figure 1.
3. General Geology
The study area belongs to Makran and East Iran zone (Nahbandan-Khash zone) based on the Iran’s structural classification [9] . Stöcklin [9] has interpreted this area as separated zone from Zagros zone along with Oman uplift, which is the continuous part of Baluchestan and India. Camp and Griffis [10] believed that occurrence of rifting in Middle Cretaceous, the Lut zone is separated from Afgan Block and an oceanic basin occurred in which flysch sediments are deposited. Bagheri [11] said that Sistan suture
Figure 1. Geographic location of the study area.
zone in east south of Iran, called flysch belt occurred post Cretaceous, which joins to Makran zone in east south of Iran. Flysch is deposited on the Upper Cretaceous ophiolites and related higher pressure metamorphic rocks, similar to accretionary wedges sequences. Based on the stratigraphic studies of Iranshahr map (1:250,000), there is no evidence of existence of older units prior to Cretaceous. The youngest sediments belonged to Quaternary and the oldest sediments occurred in the Middle Eocene.
Khash-Iranshahr road in area a passes from the column of sediments with 13 m thickness composed of conglomerate of Eocene (PLQC). In this area large scale rock slides are seen sometimes having over 10 m in dimension. In area B, the road is located on soil slopes, composed of sand gravel, shale and sandstone (EK2) (Figure 2).
4. Methodology
Based on the field studies, the potential unstable slopes including rock and soil slopes were defined. Slope morphology, discontinuity features such as dip and strike direction as well as groundwater condition were investigated. The soil and rock sampling were done to define the engineering properties of the components. The lab work were conducted and the engineering properties of soil and rock samples, including especial weight, saturated especial weight, cohesion, angle of internal friction, and grain size were determined (Table 1).
Figure 2. Geological map of the study area.
Table 1. The Geomechanic features of rock and soil materials in study area.
In this study, to measure and interpret the rock slope stability, at first, geometric of slope dip, elevation of dip, especial weight, and angle of internal friction were entered to Rock plane software. In this analysis, the slope stability was investigated in dry condition. Then, after importing the groundwater condition, the rock slope stability was interpreted and safety factor for every condition was considered. The analysis of the soil slopes was done via shear method, by which slope was divided to several sections and force conditions for every section was calculated. According to the effect of forces, several methods can be used, among them; simplified Bishop Method can provide a more realistic safety factor. The analysis points investigated in this study are included the slope dip condition, soil layering, and engineering properties of soil such as cohesion, especial weight, saturated especial weight, angle of internal friction, and survey soil slope stability.
5. Results and Discussion
5.1. Area A, Rock Slope
5.1.1. Fractures Conditions and Failures Analysis
In this area, three groups of fractures are recognized, consisting near horizontal fractures and the other two groups are perpendicular to each other, which created large blocks of rocks. These fractures were measured (Figure 3).
The measured fractures have been interpreted stereographically, showing that the majority of the fractures place in three main groups, one of them near horizontal fractures and two others are perpendicular to each other; making large rock blocks (Figure 4).
The failures analysis of the area was done by analyzing the morphological properties of the slope such as slope dip and strike, failures dip and strike, applying the angle of internal friction, using Dips program. It showed that the failures are block toppling type. Other conclusions are:
・ 32% of fractures that are almost horizontal belong to critical zone, which cause the block toppling.
・ 35% of fractures are developed in N-S strike.
・ 33% of fractures are developed in near E-W strike (Figure 5).
5.1.2. Rock Slope Stability of Area A under Dry and Saturated Conditions
At first, the study area and the inter-space of fractures were modeled in dry condition
Figure 3. Measurement of joints in the study area.
Figure 4. Rose diagram of the fractures of study area.
(a)(b)
Figure 5. Fractures location of the study area.
and the stability level was measured. To do so, the block landslide model was created by applying failures of the area and the engineering properties of the rocks such as cohesion and angle of internal friction and then the safety factor was calculated (Figure 6).
Fractures in this study area play a role as water channels and lead the precipitations to different part of this block. So, the slope condition was modeled under penetration of water condition to the fractures. The percentage of filling of the slope was interpreted and the safety factor was calculated (Figure 7).
Since the safety factor of stable and unstable boundary of a slope is 1.5, the conditions presented herein display the stability condition passing the critical line. Also, conditions in which 100% of the failure space is filled is studied, which can provide the worst condition, resulting in preventive solutions.
The results of this investigation include:
Figure 6. Wall condition under penetration of water. A = Area 1; B = When 100% of fractures spaces is filled by water; C = When 79% of fractures spaces is filled by water; D = When 80% of fractures spaces is filled by water; E = Area 2; F = When 100% of fractures spaces is filled by water; G = When 79% of fractures spaces is filled by water; H = When 80% of fractures spaces is filled by water.
Figure 7. Filling percentage of failures spaces and the safety factor.
1) The stability behavior of the block is different when water is added and the failures spaces are filled. This behavior depends on the percentage of the filled open spaces and by increasing the water penetration to the slope, the shear strength of slope decrease, which decrease the safety factor (Figure 8).
2) When 79% of the fractured spaces are filled by water, the slope stability index is 1.51.
Figure 8. Effect of water penetration on shear strength.
3) When 80% of the fractured slope is filled by water, the slope stability index is 1.49, which is the yielding point of starting instability in the block, causing the block toppling.
4) Eventually, if 100% of the fractured spaces is filled by water, the slope stability index becomes 1.18 and cause the landslide.
5.1.3. Stability Condition under Seismic Phase
In this part, the stability condition and safety factor is acquired by applying seismic index in the model. According to the seismic studies of crust in Iran, based on the standard of bylaw 2800 and based on the same acceleration seismic hazard zoning map of the study area in 75 year back period is defined, having maximum probable earthquake of 0.30. So, the aforementioned seismic acceleration is added to the slope and its stability in the seismic condition is calculated. Table 2 shows the slope stability condition in different situations. According to the slope stability, safety factor presented by US as 1.5, the results of this study are presented in Table 2.
5.2. Interpretation of Soil Slope Stability in Area B
After defining the engineering properties and morphological properties, slope dip of mentioned slope was modeled and the slope stability was interpreted. Based on the results, a part of dip of the slope shows safety factor of 0.89, displaying the unstable condition and potential to slide occurrence (Figure 9).
5.2.1. Strategies of the Stabilizing of Slope
1) Changing the dip and elevation
After defining the location of slide, different approaches for stabilizing the slope were investigated. One of the strategies of stabilizing slope is changing the dip appearance and elevation. To stabilize the aforementioned slope, the dip of sliding part has been changed based on trial and error approach and in all changes, slide was happened with different safety factors in the slope. The optimal condition of stabilizing of slope is seen in Figure 4, in this condition, slope safety factor becomes 4.77 (Figure 10).
Table 2. Results of slope stability in different situations.
Figure 9. Slope condition and unstable part.
If the dip and elevation changes applied on the soil slope stabilization, this slope will be stabilized, but these changes will form a slope with new morphology, causing formation of new sliding areas (Figure 5) and the safety factor of the new slope condition becomes 1.42 which shows the slope slide potential (Figure 11).
2) Stabilizing by construction of retaining wall
One of the other common methods of landslide hazards are retaining walls typically place in a vertical manner. To construct these walls, the vertical soil digging of the potential sliding is necessary. Figure 12 shows the modeled of sliding slope and then retaining wall is builds up by using the local materials (limestone blocks). By applying the retaining wall, the slope safety factor becomes 11.5 and the slope in this condition is completely stable (Figure 12).
3) Slope stabilizing by anchor construction
After slope modeling, the best anchor applied to the slope and the results as well as anchor features itself as shown in Figure 13, which shows the slope safety factor of 1.92 when anchor was used.
Figure 10. Slope stabilizing by changing of dip and elevation.
Figure 11. Occurrence of landslide after changing the dip and elevation.
Figure 12. Slope condition by applying retaining wall.
Figure 13. Slope condition when anchor is applied.
6. Conclusions
The existing faults in this slope act as the water channel and pass the precipitations to different part of this block. By increasing the amount of water penetration, the shear strength decreases, and therefore the slope safety factor decreases. According to the studies, if over 79% of the fracture spaces fill out with water, the slope will be unstable and landslide will occur. The slope instability to seismic activities shows that activating the seismic phase results in slope instability due to presence of fracture group perpendicular to each other, causing the slope slide in block sliding manner.
In area B, soil slope is 2 km in length and a road is constructed on. The slope stability investigation of this slope shows soil slide in this area. By changing the dip and stepped- style strategy, the slope stability was achieved. However, the new condition forms new morphology for the slope, causing spoon-like slide in the upper most part of the stair part. To solve this issue, the retaining walls were applied, which were formed by local material to ensure the slope stability. Also, the anchor approach can be used to stabilize the soil slope.