Prediction of Water Logging Using Analytical Solutions—A Case Study of Kalisindh Chambal River Linking Canal
Total Page:16
File Type:pdf, Size:1020Kb
Journal of Water Resource and Protection, 2013, 5, 624-632 http://dx.doi.org/10.4236/jwarp.2013.56063 Published Online June 2013 (http://www.scirp.org/journal/jwarp) Prediction of Water Logging Using Analytical Solutions—A Case Study of Kalisindh Chambal River Linking Canal Dipak N. Kongre, Rohit Goyal Civil Engineering Department, Malaviya National Institute of Technology, Jaipur, India Email: [email protected], [email protected] Received April 9, 2013; revised May 9, 2013; accepted May 31, 2013 Copyright © 2013 Dipak N. Kongre, Rohit Goyal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. ABSTRACT The canals are designed to transport water to meet irrigation and other water demands or to divert water from surplus basins to deficient basins to meet the ever increasing water demands. Though the positives of canal network are increase in agricultural output and improvement in quality of life, the negatives of canal introduction and irrigation, along its route, are inherent problems of water logging and salinity due to seepage from canals and the irrigation, when not man- aged properly. To plan strategies to prevent waterlogging and salinity, it is necessary to predict, in advance, the prob- able area which would be affected due to seepage. This paper presents a methodology to predict the area prone to water logging due to seepage from canal by using 2D seepage solutions to 3D field problem. The available analytical solu- tions for seepage from canals founded on pervious medium and asymmetrically placed drains, have been utilized. The area, prone to waterlogging, has been mapped using GIS. Keywords: Waterlogging; Seepage; Analytical Solutions; GIS; 3D Analyst 1. Introduction of the IGNP project [7]. The major reasons for water logging and salinization in this area are indiscriminate Growing need of food and fibre for increasing population use of irrigation water, canal seepage, sandy texture and needs more agricultural land to be irrigated, necessitating the transportation of water through canals from the res- absence of natural surface drainage [8]. The same factors ervoirs, wherever possible. In India there is a large varia- have led to rise in water table in the north-west region of tion in the rainfall in both time and space leading to sce- Haryana, India and water logging problems. About 500,000 nario of droughts and floods simultaneously. To over- ha area is waterlogged and unproductive [9,10]. It is es- come this challenge of droughts and floods and to pro- timated that in India nearly 8.4 million ha is affected by vide water to the water deficient regions, the Government soil salinity and alkalinity, of which about 5.5 million ha of India has muted river interlinking project. The Parbati is also waterlogged [11]. In the states of Bihar, Gujarat, Kalisindh Chambal river interlinking is a part of penin- Madhya Pradesh, Jammu & Kashmir, Karnataka, Kerala, sular river interlinking project. The surplus water of Par- Maharashtra, Odisha and Uttar Pradesh 117.808 thou- bati and Kalisindh basins is to be diverted to meet the sand hectare waterlogged area has been approved for re- demand in upper Chambal basin [1-4]. The construction clamation under command area development program by of this link would surely benefit the recipients of addi- Ministry of Water Resources [12]. The total waterlog- tional water but the problems of water logging and salin- ged area in canal command in India, in the year 1996, ity due to seepage from canals and irrigation need con- was around 2.189 million ha [13]. More than 33% of the sideration for planning preventive measures. In the north- world’s irrigated land is affected by secondary saliniza- west part of Rajasthan state, canal irrigation was intro- tion and/or waterlogging [14]. These statistics clearly show duced after commissioning Indira Gandhi NaharPariyo- the need to manage canal seepage and irrigation to re- jna (IGNP) to irrigate nearly 2.2 million ha of arid land. duce water logging and secondary salinization. It has increased the food production but it also introduced Canal linings are used to reduce seepage. Yao et al. waterlogging and secondary salinization problems [5,6]. [15] studied the effect of canal lining and multilayered The size of the waterlogged area in the year 1998 was soil system on canal seepage and found that the combi- 17,220 ha in stage I of the project and 800 ha in stage II nation of canal lining and a low-permeability layer below Copyright © 2013 SciRes. JWARP D. N. KONGRE, R. GOYAL 625 the canal is effective in reducing seepage. Though perfect tions for seepage problems and observed that though canal lining can prevent seepage loss, but cracks can de- numerical techniques have become more significant in velop in the lining, and the performance of the canal lin- solving practical problems of seepage theory but analyti- ing deteriorates with time [16]. Wachyan and Ruston [17] cal methods are necessary not only to develop and test studied several canals and concluded that significant see- the numerical algorithms but also to gain a deeper under- page losses occur even in a well maintained canal with standing of the underlying physics, as well as for the pa- good lining. rametric analysis of complex flow patterns and the opti- Seepage from canals can be calculated by physical, mization and estimation of the properties of seepage empirical and mathematical (analytical and numerical) fields. Swamee et al. [38] obtained an analytical solution techniques [18-21]. Empirical formulae and graphical so- for seepage from a rectangular canal in a soil layer of fi- lutions are generally used to estimate seepage losses from nite depth overlying a drainage layer using inversion of proposed canals whereas the direct measurements such as hodograph and conformal mapping technique. Bardet and inflow-outflow method, ponding method, seepage meter Tobita [39] presented finite difference approach for cal- method are used to evaluate seepage from existing canals culating unconfined seepage using spread-sheets. They [22,23]. Various empirical formulae for estimating seep- derived the finite difference equations using flux conser- age are discussed by Bakri and Awad [24]. The seepage vation in the general case of non-uniform and anisotropic losses from irrigation canals with different lining materi- permeability and boundary conditions. The flow lines als, subsurface flow and subsurface storage along the ca- and free surfaces can be obtained using their method but nal, channel longitudinal slope, have been estimated by the method cannot be adopted when systems of equations different researchers using direct measurements or elec- become large. trical resistivity [25-28]. Analytical and electrical anal- Sharma and Shakya [40] applied the Bousinessq equa- ogy solutions of seepage problems related to irrigation tion using the Laplace transform and Fourier cosine tran- canals have been presented by many authors. Zhukovsky sform in determining the phreatic surface elevation in was the first to introduce the method for solution of pro- horizontal unconfined aquifers along the canal sides. An blems involving unconfined seepage using a function exact analytical solution for the quantity of seepage from which is now well known as Zhukovsky’s function. Ve- a trapezoidal channel underlain by a drainage layer at a dernikov gave an exact mathematical solution to uncon- shallow depth obtained by using an inverse hodograph fined, steady-state seepage from a triangular and a trape- and a Schwarz-Christoffel transformation was presented zoidal canal in a homogeneous, isotropic, porous medium by Chahar [41]. It provides a set of parametric equations of large depth [29]. Vedernikov [30] solved the problem for the location of phreatic line. To reduce the convey- of seepage from a canal to the symmetrically placed col- ance losses due to seepage, many researchers have pro- lector drainages, neglecting the effect of the canal water vided methods to design the canals for minimum seepage depth and side slopes. Many Russian investigators ana- [20,38,42,43]. Ayvaz and Karahan [44] provided spread- lysed solutions for rectangular, triangular channel of zero sheet application of three-dimensional (3D) seepage mo- depth for different boundary conditions [31,32]. The so- deling with an unknown free surface. They derived gov- lution for seepage problem for a rectangular canal was erning equation using finite differences method in the ge- provided by Morel-Seytoux [33]. Sharma and Chawla neral case of anisotropic and non-uniform material prop- [34] presented solution of the problem of seepage from a erties and variable grid spacing and also modified it by canal to vertical and horizontal drainages, symmetrically the extended pressure method. Only one finite difference located at finite distances from the canal, in a homoge- equation was applied to the solution domain instead of nous medium extending up to a finite depth. The water derivation of additional finite difference equation to im- depth in the canal was assumed negligible in comparison pervious boundary conditions, inclined interfaces, etc. So- to the width. Exact solution of the problem of seepage lution to seepage under the dam has been provided using from a canal in a homogeneous medium to asymmetric this method. Ahmed and Bazaraa [45] investigated the drains located at finite distance from the canal was pro- problem of seepage under the floor of hydraulic struc- vided by Wolde-Kirkos and Chawla [35]. Goyal [36] pro- tures considering the compartment of flow that seeps vided solutions for seepage from canals founded on per- through the surrounding banks of the canal. A computer vious soil with asymmetric drainages. Algorithm to solve program, utilizing a finite-element method and capable the nonlinear integral equations to obtain the value of of handling (3D) saturated-unsaturated flow problems, seepage discharges and profile of the free surface was de- was used.