Journal of Engineering Geology

---

In this study, seepage phenomena through the left abutment of Marun dam are investigated. The Marun dam is a 170 m high rock fill dam, which regulates the waters of the Marun River, serves power generation, and flood control and provides irrigation needs. The dam site lies in the Zagros Mountains of southwest Iran. This region presents continuous series of mainly karstic limestone, marl, shale and gypsum ranging in age from Cretaceous to Pliocene. The region has subsequently been folded and faulted. All underground excavations are sited in the left abutment. The spacing of the diversion tunnels and pressure tunnel is considered to be acceptable, meaning relatively short, thus requiring 2 row grouting curtain into both embankments. Prior the reservoir impoundment, the concrete plug was constructed into the middle section of second diversion tunnel. Upstream section of tunnel was not concreted. During the first reservoir impounding, the old karst channels along ‘Vuggy Zone’ cut by the second diversion tunnel were reactivated and leakage occurred. The total amount of water leakage through the left bank of Marun dam was about. The unlined second diversion tunnel had a key role in connecting reservoir with karst conduit system. On the basis of detailed engineering geological analysis, the concept of remedial works was carried out. The main points of this concept are one of row grout curtain extension up to the section with shaly interbeds declared as watertight Asmari sequence (close to the watertight Pabdeh formation) and plugging of accessible section of main karst channel by concrete. In order to determine the seepage direction and karstification pattern, solubility studies were done. Also pinhole, XRD and XRF tests were carried out. The major joint system and interbedding cracks have predominant role in karst evolution process. Hydrogeological role of joints, perpendicu-lar to geological structure, is not negligible. As a result of these studies, seepage paths have been identified in the karstic limestone in the left abutment of the dam.

---

River bank erosion is a mechanism initiated by soil particles movement due to subsurface flow. This may occur in a soil texture at a critical hydraulic gradient. With regard to the complexity of river bank erosion processes and limited research in this field, it is of significance to investigate and to identify the effective parameters. In the present study, a physical model of a river bank was developed to achieve in-depth understanding of the effects of bank material particle sizes on porous media under various hydraulic gradients. It is concluded that the length of scour hole caused by seepage erosion may depend on the hydraulic gradients as well as Reynolds number. Further, two empirical relationships are derived on the basis of observed experimental results to estimate scour hole length resulted from river bank erosion for laminar and turbulent conditions

---

Subsurface flow contributes considerably to river flow and plays an important role in river sediment loads. This research has been focused on investigation of soil properties and bankstream slope on seepage erosion. For this purpose a series of lysimeter experiments were performed for four different slopes of bankstream by varying the soil grain sizes. The obtained results indicate that Reynolds number in porous medium plays an effective role in depth of scour hole in noncohesive layer. It was observed that the time of beginning of sediment motion decreases with an increase in the soil grain size.

---

Introduction
Pressure tunnels in hydroelectric plants are used to convey water to powerhouses. These tunnels are the sources of seepage flow to the rock formation, thus, during the water filling, they will have a low resistance to seepage and, by increasing the internal water pressure of the tunnel, the inflow force will be transferred to the rock mass. In these conditions, the cracks, pores and all other elements of the rock mass are affected by the seepage forces in all directions. This hydro-mechanical interaction affects changing the stresses and displacements of the rock mass around the tunnel and causes modifications in the permeability of rock elements during the water filling. Therefore, changes in stress distribution lead to alterations in the permeability coefficient and redistribution of the seepage field. In these conditions, since the analytical solution of the problem is not possible, the numerical analysis based on the finite element method has been used in this study.
Material and methods
In this approach, the rock mass is considered as an equivalent continuum in which the effects of discontinuities are taken into account in its material behavior. High-pressure tunnels under internal water pressure requires reinforced concrete lining to prevent hydro-fracturing. The ABAQUS software is capable of analyzing such as seepage from the tunnel, modeling of the steel bars in concrete, and taking into account hydro-mechanical interaction. Thus, the software is used for numerical analysis.
The pressure tunnel of the Gotvand dam and hydroelectric power plant (HPP) scheme is taken as a case study for the numerical simulation. Pressure tunnel of the Gotvand dam located in the southwest of Iran is taken as a case study for the numerical simulation. Among behavioral models in the software, Mohr-Coulomb failure criterion is considered to describe the rock mass, but the principle of effective stress determines the rock mass behavior. Since the concrete lining of the pressure tunnel will undergo two mechanisms of the cracking due to tension and the crushing due to compression, concrete damaged plasticity model is used to predict the response of the concrete elements. The evolution of the yield surface of the concrete lining is also controlled with tensile and compressive equivalent plastic strains, correspondingly.
In this study, the hydro-mechanical interaction is implemented based on the analysis of the pore fluid/deformation analysis, and the direct-coupled method is used to solve the governing equations of the problem. To verify the proposed model, the elastic behavior of the media is simulated to compare the numerical and the analytical solutions and good agreement is obtained. The numerical analyses are carried out the hydro-mechanical interaction with constant permeability coefficient. When cracks develop in the concrete lining at high water pressure, the properties of the concrete lining change and as a result, the stress dependent permeability of the lining and surrounding rock mass in pressure tunnels should be considered. The coefficient of permeability controls the rate of seepage flow in porous and fractured media. Although permeability represents an original property of the porous media, it can be modified when subjected to the stress variations. Instead of changing aperture, the change in the void space or volume is the typical consequence that results to change the permeability coefficient. In order to bring the model closer to the real conditions and in the validation of the new model, the influence of the permeability coefficient variations of the concrete and rock mass on the deformations and stresses of the model has been added to nonlinear analysis by USDFLD code. Increasing the water head in the tunnel during water filling is also considered with the combination of DLOAD and DISP codes in the model.
Results and discussion
Since the lining and rock mass have nonlinear properties and complex behavior, for verification of the model in ABAQUS software, the model is simulated with homogeneous, isotropic and elastic behavior. The results of seepage flow on the interface of the concrete lining and rock mass obtained by analytical and numerical solutions indicate that there is a ±5 % difference between them. Then, the results of the elastic behavior of the model show a good agreement with the results of analytical solutions. Therefore, this numerical model has been employed for the nonlinear analyses.
Finally, the optimal thickness of the concrete lining with the appropriate arrangement of the reinforcement in the reinforced concrete linings is utilized to minimize water losses from the tunnel based on the new model. Thus, the results of the analysis with the aim of reducing the water losses from the tunnel indicate that the suitable arrangement of the steel bars in the concrete lining leads to the distribution of micro cracks in the lining, and the reinforcement stress stays at a lower value with high internal water pressure. Based on the new numerical model, it is suggested that the lining should be designed with the thickness of 40 cm and the reinforcement with the diameter of 16 mm and the spacing of 20 cm.
 Conclusion
The results of the numerical model indicate that to control the seepage outflow from concrete-lined pressure tunnels, the thickness of the lining and the suitable arrangement of the steel bars in the concrete lining play a significant role in preventing excessive seepage from the tunnel./files/site1/files/124/3dadashi%DA%86%DA%A9%DB%8C%D8%AF%D9%87.pdf