Microsoft Word - ETASR_V12_N1_pp8086-8089 Engineering, Technology & Applied Science Research Vol. 12, No. 1, 2022, 8085-8089 8085 www.etasr.com Bashir et al.: Understanding the Seepage Behavior of Nai Gaj Dam through Numerical Analysis Understanding the Seepage Behavior of Nai Gaj Dam through Numerical Analysis Afsheen Bashir US-Pakistan Center for Advanced Studies in Water Mehran University of Engineering and Technology Jamshoro, Sindh, Pakistan afsheenbasheer@yahoo.com Rabia Chaudhry PhD Researcher/Expert in Geotechnical Engineering Water & Power Development Authority Punjab, Pakistan rabiach100@gmail.com Abdul Latif Qureshi US-Pakistan Center for Advanced Studies in Water Mehran University of Engineering and Technology Jamshoro, Sindh, Pakistan alqureshi.uspcasw@faculty.muet.edu.pk Uroosa Memon US-Pakistan Center for Advanced Studies in Water Mehran University of Engineering and Technology Jamshoro, Sindh, Pakistan uroosamemon11@gmail.com Naraindas Bheel Department of Civil and Environmental Engineering Universiti Teknologi PETRONAS Perak, Malaysia naraindas04@gmail.com Abstract—This study presents the seepage patterns of earth-fill dams, using critical situations by employing the finite element approach. The Nai Gaj dam is 65km northwest of Dadu city in Sindh Province, Pakistan. In this study, the seepage through the main dam body and foundation is computed and simulated for different scenarios, i.e. maximum, normal, and minimum reservoir level. Seepage analysis was conducted by using the SEEP/W sub-program of GEO-SLOPE software. Dam design parameters and dam geometry data were used as input data to compute the unknown seepage. The seepage behavior of the Nai Gaj dam is shown in terms of net flow which consists of equipotential lines, streamlines, velocity vectors, and phreatic lines. The results show the seepage flux, maximum seepage, and exit gradient at different reservoir levels. The results show that the average flow rate at normal, maximum, and minimum reservoir levels are 1.49×10 -7 cumec/m, 3.38×10 -7 cumec/m, and 2.108×10 -8 cumec/m respectively. In addition, the overall stability of the side slope of the dam is discussed. Keywords-seepage; earth fill dam; phreatic line; SEEP/W; finite element modelling I. INTRODUCTION Pakistan is basically an agricultural country. The agriculture in Pakistan depends mainly on water irrigation. In this regard, many dams have been constructed in the past and some others are planned to be constructed in the future. Embankment dams are important and costly civil engineering structures that provide an essential infrastructure for the management of water [1, 2]. One of the critical aspects of dam design is the analysis of stability and safety of the earth structure under various operating conditions. Dams are built for specific functions, i.e. storage of water for domestic use and irrigation, flood control, entertainment, fishing, and hydroelectric power generation. Design and construction of dams have to meet safety criteria and satisfy specific requirements. Stability of the dam against slope failure is an essential component for the design, operation level under all flood and drought circumstances [3, 4]. The estimated number of the earth fill dams around the world is 11192 [5, 6]. A review of historical dam failures indicates that nearly half of all failures of embankment dams have occurred due to seepage induced internal erosion. Detailed studies of seepage analysis through and beneath the dam body need to be parts of design, construction, and maintenance of any dam component. Most geotechnical studies related to dams involved seepage and stability analysis considering the pore pressure variations in unsaturated soil conditions. Seepage failures occur by rainfall or water infiltration that changes the distribution of pore water pressure. This effect leads to changes in the stress and affects the seepage process of the soil. An analyst must have a good understanding of seepage, its solutions, and the preventive measures for its monitoring. Finite element analysis is extensively used in the geotechnical studies for seepage and stability analysis considering the pore pressure variations. In case of dam seepage, failures occur by rainfall or water infiltration that changes the distribution of pore water pressure. This effect leads to changes in the stress and affects the seepage process of the existing soil. To make the dam safe, it is important to take great care on stability and seepage control Corresponding author: Naraindas Bheel Engineering, Technology & Applied Science Research Vol. 12, No. 1, 2022, 8085-8089 8086 www.etasr.com Bashir et al.: Understanding the Seepage Behavior of Nai Gaj Dam through Numerical Analysis while designing embankment dams. Embankment dams are generally constructed with various zones of materials such as two permeable shells and impermeable core in the middle [7]. The main reasons behind embankment dams’ failures are seepage, hydraulic failures, piping through the dam body, and structural damage due to earthquakes [8]. Permeable materials such as gravel or coarse sand when used in dam foundation lead to heavy seepage [5, 9-11]. In the past, embankment dams mostly faced seepage problems. As water in the reservoir seeks the paths of least resistance through the dam body and its foundation, it may pose some hazards on the dam safety. Excessive seepage may lead to dam failure if it is not treated and controlled properly. If seepage flow can continue unhindered in appreciable quantities, then the seepage force may erode fine soil particles and wash them out causing piping failure of the dam or create uplift problems. Seepage, therefore, may be considered as one of the most common safety hazards for embankment dams. The basic problem the designers and operators face are to discern how seepage is affecting a particular dam and what measures, if any, must be taken to ensure that the seepage does not unfavourably affect the safety of the dam. This paper presents a case study of the Nai Gaj Dam and highlights the probability of pre and post seepage issues. The SEEP/W GEO-SLOPE software was utilized. It can work quite well on seepage analysis in homogenous and non- homogenous embankment dams [12]. Slope stability analysis confirms the safety of the embankment dam. Different loading conditions can be studied, for example end of construction, steady state, rapid draw down, and the earthquake loading scenario. The aim of this study is to analyze the steady state seepage behavior of the dam foundation of the Nai Gaj Dam for various scenarios, namely normal, maximum, and minimum reservoir level conditions. In addition, the stability of the dam was analyzed. II. DESCRIPTION OF THE STUDY AREA Nia Gaj Dam is located on the Gaj river in the Kirther mountain range. It is located 65km north-west of the Dadu city, Sindh province. The purpose of the Dam Project is to provide irrigation to agriculture lands downstream of the dam. The foundation rocks of this dam project are generally weak and fragile. Sandstone and silty sandstone beds encountered at the core trench are weak and friable with jointings and open bedding planes [13, 14]. These rock characteristics may indicate poor rock conditions, low shear strength, and probability of increased seepage through the foundation [15]. III. MATERIALS AND METHODS Initially, the cross section of the dam was drawn in CAD format and imported in SEEP/W as shown in Figure 1. The material zones in the dam are sandy gravel, random fill, central impervious clay core, sand filter, coarse filter, and drainage blanket. The foundation of the dam mainly consists of sandy siltstone. The dam is 59m high, 1137m long, having a maximum design water level of 56.6m with 10m crest width. The core has an upstream side slope of 1H:3V, and a downstream slope of 1H:4V. The upstream and downstream shoulders have slopes 1V:2.5H. Drainage features, including a chimney drain and horizontal drainage blanket below the downstream shell, have been incorporated in the design to relieve the pore water pressure within the dam body and to dispose the seepage water in a safe manner to prevent erosion and boiling. Numerical analysis was carried out on a typical section of the dam. A non- homogenous section is selected for the analysis of the response of complex materials by varying the water levels in order to understand and analyze the seepage phenomenon. The physical, mechanical and index characteristics of the materials were taken from the references mentioned in Table I. These material properties were assigned to the finite element model of the dam. The generated FEM mesh is shown in Figure 2. Dirichlet and Neumann boundary nodes were assigned at the inlet and outlet boundary conditions [16, 17]. It was accepted that the discharge is gravity driven flow while the capacity changes due to water compressibility, soil structure compressibility, and changes in matric suction (i.e. drainage), in spite of the fact that soil structure compressibility is related to pore water pressure changes. Subsequently, all the stresses inside the domain are accepted as constant. No loading and unloading of soil mass is included. Once the required data were imported to the SEEP/W, it verified the mesh development (7047 nodes and 6913 elements). Therefore, now the model can be used for simulations and analysis of the results. Fig. 1. X-section of th eNai Gaj dam. CL CREST EL=184.00 1 2.50 2.50 1 D / S U / S 1 4 1 3 H.F.L EL=181.700 MCL=177.50 1 2 4a 5R2 5R1 3 4b 2 6 7 CONSOLIDATION GROUTING CURTAIN GROUTING Engineering, Technology & Applied Science Research Vol. 12, No. 1, 2022, 8085-8089 8087 www.etasr.com Bashir et al.: Understanding the Seepage Behavior of Nai Gaj Dam through Numerical Analysis Fig. 2. Mesh formation for the selected non-homogeneous section. Results are generated for the following scenarios: (i) normal reservoir level (150m), (ii) maximum reservoir level (177.5m), and (iii) minimum reservoir level (130m). The boundary conditions which are used in the analysis are: • Dirichlet and Neumann boundary nodes are assigned at the inlet and outlet boundary conditions [17]. • Zero flux (Neuman) condition is considered as the boundary node at the bottom of foundation for the all the given conditions. The results are analyzed in terms of development of net flows for different scenarios, along with velocity vectors, quantity of seepage value through the core and foundation of the dam for different scenarios, and exit gradient. In addition, the stability of the dam was determined based on slope stability analysis. A. Fundamental Equations In this study, finite element method is used for the modelling of the groundwater flow in porous media by utilizing the following differential equations to analyze the seepage through the body of the dam at its foundation. Equations are derived on the principle of mass conservation law [17, 18]. The partial differential equation (1) is used in SEEP/W: � �� (Kx �� �� � + � �� � �� �� � � � �� �� (1) where H is the hydraulic head, Kx and Ky are the horizontal and vertical hydraulic conductivity respectively, Q is the discharge, t is the time θ-domain and volumetric water content. The above equation is a two-dimensional and nonlinear 2 nd order PDE and is used for transient flow. It is derived from Darcy’s Law: v = − K ∇H (2) where: ∇H � � �� ���� �� � (3) where v is the Darcian velocity, K is the hydraulic conductivity of the soil material, and ∇H is the hydraulic gradient. Equation (1) is dependent on time, and tells that the flow that enters an elemental volume minus the flow leaving the elemental volume at a point is equal to the volumetric water content at that point. If the volume of water that enters is equal to the volume of water that leaves at steady state situations, then the right side of (1) becomes equal to zero. � �� ( Kx �� �� � + � �� � K� �� �� � � Q � 0 (4) IV. RESULTS AND DISCUSSION The material properties which are used in the selected dam sections were taken from past studies and from a feasibility report on the Nai Gaj dam [19]. The inherent material parameters used in the dam are presented in Table I. The mesh was generated by SEEP/W and got a no error report for further computations. Also, the given up- and down-stream boundary conditions were assigned. The physical and mechanical parameters of the materials during the time of construction were adopted. SLOPE/W then verified the development model and found no error in the generation of the numerical model. TABLE I. INHERENT MATERIAL PARAMETERS OF THE NAI GAJ DAM [1, 19-23] Dam component Material Unit weight (kN/m 3 ) Cohesion (kPa) Friction angle (deg) Hydraulic conductivity (m/day) C ᴓ Foundation Sandy siltstone 20.4* 12 29* 0.00063* Upstream Sandy gravel 21.5** 0 37* 86.4*** Downstream Random fill 18.85* 0 34* 0.0263* Downstream shell Riprap 19.5* 0 34* 8640*** Core Clay 18.85* 9.57 30* 0.000263* Filter blanket Sand filter 18.85* 0 36* 26.33* The seepage analysis of the Nai Gaj dam and its foundation at different reservoir levels was conducted with the use of the SEEP/W. It is assumed that the downstream drain will drain out the seepage and no seepage will develop downstream. For this purpose, flow net has been drawn as shown in Figures 3-5 for the selected dam sections. It was found that if the reservoir is kept at 150m and 130m there is no pronounced seepage observed. Engineering, Technology & Applied Science Research Vol. 12, No. 1, 2022, 8085-8089 8088 www.etasr.com Bashir et al.: Understanding the Seepage Behavior of Nai Gaj Dam through Numerical Analysis Fig. 3. Flow net of the selected section at maximum reservoir level (177.5m). Fig. 4. Flow net of the selected section at normal reservoir level (150m). Fig. 5. Flow net of the selected section at minimum reservoir level (130m). TABLE II. CALCULATED SEEPAGE FLUX AND EXIT GRADIENT Upstream reservoir level (m) Water head (m) Seepage flux (cumec/m) F o u n d a ti o n 177.5 (max) 177.5 1.29 x 10 -7 150 (normal) 150 4.53 x 10 -8 130 (min) 130 6.92 x 10 -9 C o r e 177.5 (max) 177.5 3.38 x 10 -7 150 (normal) 150 1.49 x 10 -7 130 (min) 130 2.106 x 10 -8 The exit gradient was calculated at the downstream toe of the dam for 3 reservoir levels, i.e. 177.5, 150, and 130 m (Table III). The factor of safety against piping is described as follows: FS = critical hydraulic gradient/exit gradient (5) The value of the critical hydraulic gradient is taken as 1. A dam is considered as safe against piping, if the value of the factor of safety is between 4 and 5. The results indicate that the values of factor of safety for normal, average, and maximum water level are 7, 6, and 5 respectively. This shows that the dam is safe against piping failure even at the maximum reservoir level of 177.5m. Table III shows that the exit gradient of the selected section was mostly found less than 1 for all the water levels which conforms the safety criteria of the dam. It further support the current design of the dam. From Tables II and III, it can be observed that both the exit gradient and seepage flux are functions of water level. Generally, the stability of a slope at the downstream embankment is critical at steady state seepage condition. It is assumed that the steady seepage condition may develop when the reservoir is full at maximum retention level, i.e. at EL 177.5m. The phreatic surface is estimated using the flow net development by the SEEP/W program. Furthermore, to support the outcomes, deterministic slope stability analysis was Engineering, Technology & Applied Science Research Vol. 12, No. 1, 2022, 8085-8089 8089 www.etasr.com Bashir et al.: Understanding the Seepage Behavior of Nai Gaj Dam through Numerical Analysis performed for different loading conditions, such as the end of construction and steady state seepage. The material properties used are mentioned in Table I and the results were discussed in terms of their factor of safety as shown in Table III. The computed FOSs at different loading conditions were further compared with the USACE standards [24]. It was noted that in slope stability analysis, the computed FOSs for different loading conditions are greater than the standard values which indicates safety of the side slopes of the dam as given in Table III. TABLE III. EXIT GRADIENT AND FACTOR OF SAFETY AGAINST PIPING Upstream pond level (m) Exit gradient (m) Factor of safety against piping 177.5 (max) 0.19 5 150 (average) 0.16 6 130 (normal) 0.14 7 TABLE IV. FACTOR OF SAFETY Load Slope Static equilibrium condition Required* Computed End of construction Downstream 1.30 1.7 Steady state condition Downstream 1.50 1.8 V. CONCLUSION In this paper, seepage analysis of Nai Gaj dam was carried out with the sub program SEEP/W of the GEO-SLOPE software. The flow rates of the seepage computed at the core of the dam at maximum (177.5m), normal (150m), and minimum (130m) reservoir level were found to be 1.49×10 -7 cumec/m, 3.38×10 -7 cumec/m, and 2.108× 10 -8 cumec/m respectively. It was observed that as the reservoir level increased, the seepage quantity also increased. For the three tested scenarios of the reservoir level, the value of the exit gradient was found to be less than 1 which confirms that the dam is safe against piping due to seepage. In addition, slope stability analysis shows that the dam is stable at the end of construction and at steady state seepage condition. ACKNOWLEDGEMENT The authors are thankful to the US-Pakistan Center for Advanced Studies in Water, Mehran University of Engineering and Technology, Jamshoro, Sindh, Pakistan for providing the facilities to carry out this research study. REFERENCES [1] A. H. Bhutto et al., "Mohr-Coulomb and Hardening Soil Model Comparison of the Settlement of an Embankment Dam," Engineering, Technology & Applied Science Research, vol. 9, no. 5, pp. 4654–4658, Oct. 2019, https://doi.org/10.48084/etasr.3034. [2] A. H. Bhutto, S. Zardari, G. S. Bhurgri, M. A. Zardari, R. Bhanbhro, and B. A. Memon, "Post Construction and Long Term Settlement of an Embankment Dam Computed with Two Constitutive Models," Engineering, Technology & Applied Science Research, vol. 9, no. 5, pp. 4750–4754, Oct. 2019, https://doi.org/10.48084/etasr.3070. [3] A. Liaghat, A. Adib, and H. R. Gafouri, "Evaluating the Effects of Dam Construction on the Morphological Changes of Downstream Meandering Rivers (Case Study: Karkheh River)," Engineering, Technology & Applied Science Research, vol. 7, no. 2, pp. 1515–1522, Apr. 2017, https://doi.org/10.48084/etasr.969. [4] M. A. M. Ismail, S. M. Ng, and E. K. Gey, "Stability analysis of kelau earth-fill dam design under main critical conditions," Electronic Journal of Geotechnical Engineering, vol. 17 W, pp. 3209–3219, 2012. [5] M. Foster, R. Fell, and M. Spannagle, "A method for assessing the relative likelihood of failure of embankment dams by piping," Canadian Geotechnical Journal, vol. 37, no. 5, pp. 1025–1061, Oct. 2000, https://doi.org/10.1139/t00-029. [6] M. Foster, R. Fell, and M. Spannagle, "The statistics of embankment dam failures and accidents," Canadian Geotechnical Journal, vol. 37, no. 5, pp. 1000–1024, Oct. 2000, https://doi.org/10.1139/t00-030. [7] S. S. Athani, Shivamanth, C. H. Solanki, and G. R. Dodagoudar, "Seepage and Stability Analyses of Earth Dam Using Finite Element Method," Aquatic Procedia, vol. 4, pp. 876–883, Jan. 2015, https://doi.org/10.1016/j.aqpro.2015.02.110. [8] S. O. Osuji and S. A. Adegbemileke, "Phreatic Line and Pore Pressure Stresses in Zoned Rockfill Dam," Asian Journal of Science and Technology, vol. 6, no. 5, pp. 1447–1454, 2015. [9] R. Fell, C. F. Wan, J. Cyganiewicz, and M. Foster, "Time for Development of Internal Erosion and Piping in Embankment Dams," Journal of Geotechnical and Geoenvironmental Engineering, vol. 129, no. 4, pp. 307–314, Apr. 2003, https://doi.org/10.1061/(ASCE)1090- 0241(2003)129:4(307). [10] K. S. Richards and K. R. Reddy, "Critical appraisal of piping phenomena in earth dams," Bulletin of Engineering Geology and the Environment, vol. 66, no. 4, pp. 381–402, Nov. 2007, https://doi.org/10.1007/s10064- 007-0095-0. [11] C. F. Wan and R. Fell, "Investigation of Rate of Erosion of Soils in Embankment Dams," Journal of Geotechnical and Geoenvironmental Engineering, vol. 130, no. 4, pp. 373–380, Apr. 2004, https://doi.org/10.1061/(ASCE)1090-0241(2004)130:4(373). [12] T. Mohammed et al., "Seepage through homogenous and non- homogenous earth dams: Comparison between observation and simulation," Electronic Journal of Geotechnical Engineering, vol. 11, Jan. 2006. [13] N. Adamo, N. Al-Ansari, S. Knutsson, J. Laue, and V. Sissakian, "Mosul Dam : A Catastrophe yet to unfold," Engineering, vol. 9, no. 3, pp. 263– 278, 2017. [14] Report on Foundation Geology, Nai Gaj Dam Project. Islamabad, Pakistan: WAPDA, 2013. [15] V. K. Sissakian, N. Adamo, N. Al-Ansari, S. Knutsson, and J. Laue, "Defects in Foundation Design Due to Miss-Interpretation of the Geological Data : A Case Study of Mosul Dam," Engineering, vol. 9, no. 7, pp. 683–702, 2017. [16] I. Arshad and M. M. Babar, "Comparison of SEEP/W Simulations with Field Observations for Seepage Analysis through an Earthen Dam (Case Study: Hub Dam - Pakistan)," International Journal of Research, vol. 1, no. 7, pp. 57–70, Sep. 2014. [17] Seepage Modelling with SEEP/W. Calgary, Canada: GEO-SLOPE International Ltd, 2012. [18] I. Arshad and M. M. Babar, "Finite Element Analysis of Seepage through an Earthen Dam by using Geo-Slope (SEEP/W) software," International Journal of Research, vol. 1, no. 8, pp. 619–634, Sep. 2014. [19] Project Planning Report of Nai Gaj Dam. Islamabad, Pakistan: WAPDA, 2006. [20] A. J. Schleiss and R. M. Boes, Dams and Reservoirs under Changing Challenges. Florida, USA: CRC Press, 2011. [21] J. Charles, C. Abbiss, and E. Gosschalk, An engineering guide to seismic risk to dams in the United Kingdom. London, UK: Building Research Establishment, 1991. [22] R. F. Craig, Soil Mechanics. Boston, MA, USA: Springer, 2013. [23] T. L. Bergman, F. P. Incropera, D. P. DeWitt, and A. S. Lavine, Fundamentals of Heat and Mass Transfer. New York, NY, USA: John Wiley & Sons, 2011. [24] Slope Stability. Washington, DC, USA: USACE, 2003.