Water Resource Engineering
Water Resource Engineering
Surface Subsurface
methods methods
Sprinkler Drip
irrigation irrigation
a) Surface Irrigation: Just flooding water. About 90% of
the irrigated areas in the world are by this method.
(e) The alignment should not pass through the thickly populated areas,
religious places and burial grounds.
(7) Final Survey
(8) Preparation of Drawing
(9) Office / Design Works
(10) Final Location Survey
(a) The centre line of the main and branch canals are
marked with concrete pillars with an interval of 30
m or 50 m
(b) The total land width required for the main and
branch canal should be marked with pillars at
suitable intervals.
IRRIGATION PROJECT REPORT
Important Irrigation Projects of India
(1) Bhakra Nangal Project
(2) Damodar Valley Project
(3) Farraka Barrage Project
(4) Gandak Project
(5) Godavari Barrage Project
(6) Hirakud Project
(7) Mahanadi Delta Project
(8) Mayurakshi Project
(9) Kansabati Project
(10) Kosi Project
(11) Tungavadra Project
(12) Ukai Project
Quality of Irrigation Water
Just as every water is not suitable for human
beings, in the same way, every water is not
suitable for plant life. Water containing
impurities, which are injurious to plant growth,
is not satisfactory for irrigation, and is called the
unsatisfactory water.
The various types of impurities which make the
water unfit for irrigation are classified as:
1) Sediment Concentration in water
2) Total concentration of soluble salts in water
3) Concentration of sodium ions compared to
other cations
4) Concentration of potentially toxic elements
present in water
5) Bicarbonate concentration as related to
concentration of Ca + Mg
6) Bacterial Concentration
(1) Sediment Concentration
The Sediment of fine silt improve the fertility of
land. Where as other types of sedimentation
decreases the fertility. Again excessive
sedimentation creates trouble in canal and
reservoirs. So water should not contain
excessive suspended sediment.
In general, ground water or surface water from
reservoirs etc does not have sufficient sediment
to cause any serious problems in irrigation
(2) Total concentration of soluble salts
Salts of calcium, magnesium, sodium and
potassium, present in the irrigation water may
prove injurious to plants. When present in
excessive quantities, they reduce the osmotic
activities of the plants, and may prevent
adequate aeration, causing injuries to plant
growth. The injurious effects of salts on the
plant growth depend upon the concentration of
slats left in the soil.
The salinity concentration of the soil solution
(Cs) after the Consumptive water (Cu) has been
extracted from the soil, is given by
The salt concentration is generally expressed by
ppm (parts per million or by mg/l, both units
being equal.
The critical salt concentration in the irrigation
water depends upon many factors, however,
amounts in excess of 700 ppm are harmful to
some plants and more than 2000 ppm are
injurious to all crops.
The salt concentration is also measured by
determining the electrical conductivity of water.
They are directly proportional to each other.
Based on the value, the irrigation water is
classified as given in the table:
(3) Concentration of sodium ions
Hence, it is S2.
As electrical conductivity is 200 micro mhos per
cm at 25 degree C, it is C1
Hence it is S2-C1 type of water.
(4) Concentration of potentially toxic elements
A large number of elements such as boron, selenium
etc may be toxic to plants. Traces of Boron are essential
to plant growth but its concentration above 0.3 ppm
may prove toxic to certain plants.
Even for the most tolerant crops, boron concentration
should not exceed 4 ppm. Boron is generally present in
various soaps. The waste water containing soap should
therefore be used with great care in irrigation.
Surface Subsurface
methods methods
Sprinkler Drip
irrigation irrigation
a) Surface Irrigation: Just flooding water. About 90% of
the irrigated areas in the world are by this method.
Crop period is slightly more than the base period, but for
all practical purposes, they are taken as one and the same
thing, and generally expressed in B days.
Gross Command Area (G.C.A.) : The gross command
area is the total which can be commanded or irrigated
by a canal system.
Wf = W s + R f + D f
Where Rf = Surface runoff from the field
Df = Deep percolation loss below the root
zone.
(c) Water Use Efficiency (ηu)
It is the ratio of water used beneficially to the amount
of water delivered.
ηc = (Wu / Wf) * 100
Wu= Water used beneficially
Wf = water delivered to the field
(d) Consumptive Use Efficiency (ηcu)
(e) Water Storage Efficiency (ηs)
❖ Penman Method
❖ Blaney-Criddle method
❖ Hargreaves class A pan evaporation method
Blaney-Criddle Method
Assumes that main factors are Temperature and Day-time hours.
Cu= k.f
Where, k = Crop factor
f = Monthly consumptive use factor
= p/40 (1.8t + 32)
p = percent of day time hours of the year
t = Temperature in degree centrigade
Penmann Method
ETo = C x [W x Rn + (1 - W) x f(U) x (Es – Ea)]
where
ETo = Reference crop evapotranspiration in mm/day
C = Factor to compensate for the effect of day and night weather
condition
W = Temperature related weighting factor
Rn = Net radiation in mm/day = (1 – α) x Rs - Rnl
α = Reflection coefficient, (albedo) the numerical values of which for
most crops is 0.25
Rs = Incoming short wave radiation
= (0.25 + 0.5 x (SH / SHM) x Ra in mm/day
SH = Mean value of actual sunshine hours for the period considered
SHM = Maximum possible sunshine hours for the period considered
Ra = Extra terrestrial radiation in mm/day
Rnl = Net long wave radiations in mm/day
= F1 x F2 x F3
F1 = (2.01 x 10-9) x (T + 273)4
F2 = 0.34 – 0.044 √Ea
F3 = 0.1 + 0.9 x (SH / SHM)
f(U) = Wind related function
= 0.27 x (1 + U2 / 100)
Ea = Actual vapor pressure in mbar
= Es x RH
Es = Saturation vapor pressure in mbar
= 6.11 x exp [17.27 T / (237.3 + T)]
RH = Relative humidity
Crop Evapotranspiration, Cu or ETc = Kc X ETo
Typical Kc for different types of full grown crops
Determination of Irrigation Requirement of
Crops
In order to determine the irrigation requirement of a
certain crop, during its base period, the following
terms are used:
• It is essential to
maintain readily
available water in the
soil if crops are to
make satisfactory
growth.
Depth and Frequency of Irrigation
Example 2
If wheat requires about 7.5 cm of water after every 28 days, and the base period for
wheat is 140 days, find out the value of delta for wheat.
Solution.
No. of required waterings = 140/28 = 5
The depth of water required each time = 7.5 cm.
:. Total depth of water reqd. in 140 days = 5 x 7.5 cm = 37.5 cm
Hence, ∆ for wheat = 37.5 cm. Ans.
Average Approximate Values of ∆ for Certain Important Crops
Duty of Water.
GATE Question
Total discharge (as both wheat and gram are rabi crops) =
0.46 +0.96 = 1.42 cumec
GATE 2015
Ans : (B)
GATE 2014
GATE 2013
Answer
Irrigation water requirement= 48 - 8 cm = 40 cm
=0.4 m
B=10 days
D=8.64 x 10/0.4
=216 hectares/cumec
GATE 2012
Answer :
• Nature of • Deforestation
drainage basin
• Coastal influences
Physical Causes:
• Climatological factors:
– E.g. Bangladesh
Human causes of
flooding
1. Deforestation
– Increase run-off
– Smooth surfaces
served by a dense
network of drains,
gutters and
underground sewers,
effectively increasing
drainage density
Other Human Causes
• Failure of hard engineering such as a dam
burst can produce catastrophic flooding
– Dam burst in an iron ore mine near Seville,
Spain 1999
Loss of industrial
Inconveniences of recovery
production, traffic
Indirect after flood, increased
disruption, loss of
vulnerability of survivors.
customers.
Benefits of flooding
• Replenishes groundwater reserves
▪ Ninety-five dead
▪ Over 2 million homeless.
▪ Half a million evacuees are living in relief
camps with disease prone conditions
▪ Damaging landslides hindered relief
operations
INDIA
UTTAR PRADESH
21.90%
OTHER
STATES
43.9%
Predicting To reduce flood damage, we need to
FLOODS know the magnitude and frequency of
the floods a watershed is likely to get.
Flood Hazard Mitigation
Measures
FLOOD MANAGEMENT
APPROACH
ORGANISATIONS
• STATE FLOOD CONTROL DEPARTMENTS
• BRAHMAPUTRA BOARD
STRUCTURAL MEASURES
• FLOOD FIGHTING
• FLOOD INSURANCE
CONCEPT
• AN IMPORTANT NON-STRUCTURAL MEASURE.
• REGULATES LAND USE IN FLOOD PLAINS TO
RESTRICT DAMAGE BY FLOODS.
• INVOLVES DEMARCATION OF ZONES IN FLOOD
PLAINS COMPATIBLE WITH FLOOD RISKS
INVOLVED.
FLOOD PLAIN ZONING
WATER SUPPLY
III III II I
II
I
FACTORY
HOSPITAL
GOVT. OFFICE POWER HOUSE
PARK
ZONE REGULATION
PRIORITY - I
ZONE REGULATION
PRIORITY - II
ZONE REGULATION
PRIORITY - III
RIVER CHANNEL
INDIA
FLOOD FORECASTING
BEGINNING
• STARTED BY CENTRAL WATER COMMISSION
• YEAR 1958
• RIVER YAMUNA
■ DATA COLLECTION
■ DATA TRANSMISSION
■ DISSEMINATION OF FORECAST
FLOOD FORECASTING
DATA COLLECTION
• HYDROLOGICAL
RIVER WATER LEVEL
RIVER DISCHARGE
• HYDROMETEOROLGICAL
RAINFALL
OTHER PRECIPITATION
eg. SNOW, HAIL ETC.
FLOOD FORECASTING
DATA TRANSMISSION
• WIRELESS ■ SATELLITE
• TELEPHONE ■ TELEGRAPH
• FAX ■ TELEX
FLOOD FORECASTING
DATA ANALYSIS & FORECAST FORMULATION
• CORRELATION
• MATHEMATICAL MODELS
FLOOD FORECASTING
DISSEMINATION
CIVIL AUTHORITY
PRESS, A I R
DOORDARSHAN
OTHER USERS
LOCAL
REPRESENTATIVES
FLOOD FORECASTING AND WARNING
PRESENT SYSTEM - SCHEMATIC DIAGRAM
RIVER STAGE & RIVER STAGE &
DISCHARGE AT DISCHARGE AT
WEATHER FORECASTS BASE STATION FORECASTING STATIONS RAINFALL
FORMULATION OF
C.E., S.E. CWC HQ MINISTRY OF WR
FORECAST
CIVIL AUTHORITIES
PRESS RADIO STATE GOVT.
ENGG. AUTHORITIES
P&T DOORDARSHAN CONTROL ROOM
DEFENCE INDUSTRIES
POLICE FLOOD COMMITTEE
RAILWAYS
FLOOD MANAGEMENT
scope
27
17 15
9 8
3 2
23
15 14
11 11
6
4 4 3 2 2 2 1
UP R AP B. A
A
A RA ND
A
AT
DE CT
AM SS
MP
W.
AN TTAR
NH
IH
AN
K
AL
LH
B S R I T A
N
TA
S JA R SH KH
RY
CH
A O D&
GU
A
A
U
’R R
HA
RN
A
M JH
KA
FLOOD MANAGEMENT ACTIVITIES
During flood
Infiltration Trenches are narrow ditches filled with gravel that intercept runoff
from upslope impervious areas. They provide storage volume and
additional time for captured runoff to infiltrate the native soil
below.
Continuous Porous Pavement systems are excavated areas filled with gravel
and paved over with a porous concrete or asphalt mix. Normally all
rainfall will immediately pass through the pavement into the gravel
storage layer below it where it can infiltrate at natural rates into the
site's native soil. Block Paver systems consist of impervious paver blocks
placed on a sand or pea gravel bed with a gravel storage layer
below. Rainfall is captured in the open spaces between the
blocks and conveyed to the storage zone and native soil below.
Rain Barrels (or Cisterns) are containers that collect roof
runoff during storm events and can
either release or re-use the rainwater during dry
periods.
CN (SCS) Rational
Indirect Methods
Strange’s Table
Estimation of Flood Magnitude
Deterministic
Probability
• The probability P(A) of an event A is the relative
number of occurrences of event A after a very
large number of trials.
• It will be convenient to define the sum of the
probabilities of all the possibilities as one.
• Σi (P(Xi) = 1
• This means that a probability can range from zero
to one
• 0≤P≤1
Probability
Probability may be described as the relative frequency with which an event occurs
or will occur in the long run.
It is measured normally on a scale ranging from 0 (meaning impossible) to 1 (meaning certain)
ABSOLUTE CERTAINTY 1
That you will die some day
0.9
0.8
0.7 That you ( a man) will be taller than 1.65m.
0.6
Probability
0.3
0.2
That you will win the lottery
0.1
ABSOLUTE IMPOSSIBILITY 0 That you can swim across the Atlantic
Probability Analysis
Histogram
0.24
0.2
0.16
Relative Frequency
0.12
0.08
0.04
0
450 500 550 600 650 700 750 800 850 900 950
0.954
0.997
80
Cumulative %
60
40
20
450 500 550 600 650 700 750 800 850 900 950
Rainfall classes (mm)
Return Period
• Random variable:
• Threshold level:
• Extreme event occurs if:
• Recurrence interval:
• Return Period:
Average recurrence interval between events equaling or
exceeding a threshold
• If p is the probability of occurrence of an extreme
event, then
or
21
Relation between Probability and Return period
• What is the probability ‘p’ that in a certain period of time (i.e. 1 year),
the intensity of one event ‘X’ (rainfall) of duration ‘d’ (hours or days)
is bigger than ‘xt’ (100 mm)?
• P(X>100) when X is rainfall of i.e. 24 hs duration expressed in mm.
• We assume that the events are independents
• So we start observing events. You wait for (1 year) and you note
down the rainfalls in this year. For each year observation, there are
two possible outcomes
– “success" X ≥ xt (probability p)
– "failure" X < xt (probability 1-p)
Since the observations are independent, the probability of a
recurrence interval of duration T is the product of the probabilities
of (T-1) failures followed by one success:
First year: failure -> (1-p)
Second year: failure -> (1-p)
Third year: failure -> (1-p)
Forth year: Success!! -> p
So the probability of success in general is (1-p)*(1-p)*(1-p)*p=
(1-p)(t-1) * p
Assuming that the series of data is infinite, the E(T) can be
expressed as:
29
EV-I (Gumbel) Distribution
The probability of occurrence of an event equal
to or larger than a value x0 is
(2)
In which y is a dimensionless variable given by
y= α (x-a) a =
Thus
(3)
Where = mean and σx = standard deviation
In Practice we require the value of X for a given P.
hence, we can write yp = -ln[-ln(1-P)]
And as P = 1/T
where
Theoretically Gumbel’s assumption is valid for an
infinite sample size i.e. N ∞
But practically N is finite. Hence modification
based on sample size is required for
practical use.
Frequency Factors
• If the T is outside the ranges of the observed data,
Theoretical probability distribution is used
• Chow proposed using:
(1)
• where f (x)
X
33
EV-I (Gumbel) Distribution
The probability of occurrence of an event equal
to or larger than a value x0 is
(2)
In which y is a dimensionless variable given by
y= α (x-a) a =
Thus
(3)
Where = mean and σx = standard deviation
In Practice we require the value of X for a given P.
hence, we can write yp = -ln[-ln(1-P)]
And as P = 1/T
where
Theoretically Gumbel’s assumption is valid for an
infinite sample size i.e. N ∞
But practically N is finite. Hence modification
based on sample size is required for
practical use.
Log-Pearson Type III Distribution
Steps:
1) Prepare the annual maximum flood data (x)
2) Take logarithm (base 10)
z= log(x)
3) The T year flood estimate is given by
4) From table calculate Kz based on T and Cs
T (years)
Cs 1 2 5 10 25 50 100 200
3 -0.667 -0.396 0.42 1.18 2.278 3.152 4.051 4.97
2.9 -0.69 -0.39 0.44 1.195 2.277 3.134 4.013 4.904
2.8 -0.714 -0.384 0.46 1.21 2.275 3.114 3.973 4.847
2.7 -0.74 -0.376 0.479 1.224 2.272 3.093 3.932 4.783
2.6 -0.769 -0.368 0.499 1.238 2.267 3.071 3.889 4.718
2.5 -0.799 -0.36 0.518 1.25 2.262 3.048 3.845 4.652
2.4 -0.832 -0.351 0.537 1.262 2.256 3.023 3.8 4.584
2.3 -0.867 -0.341 0.555 1.274 2.248 2.997 3.753 4.515
2.2 -0.905 -0.33 0.574 1.284 2.24 2.97 3.705 4.444
0.1 -2.252 -0.017 0.836 1.292 1.785 2.107 2.4 2.67
0 -2.326 0 0.842 1.282 1.751 2.054 2.326 2.576
5) After finding zT the corresponding value of xT
is obtained as xT =antilog(zT)
Log Normal Distribution
Special case of log-pearson type III distribution
where Cs=0
Procedures are same except the Kz values is
fixed and taken from the table as Cs=0
Confidence Interval
The Confidence interval indicates the limits
about the calculated value between which the
true value can be said to lie with a specific
probability based on sampling error
POPULATION
SAMPLE
For a confidence probability c, the confidence
interval of the flood is between values x1 and
x2. where x1 and x2 is given by :
Reliability = (1 - 1/T)n
Risk
Example
What is the probability of at least one 50 yr flood in a
30 year period, where the probability of success in
any year is 1/T = 1/50 = 0.02
6-h UH 92 79 64 52 40 30 20 14 10 6 0
ordinates
(m3/s)
Design Storm (Standard Project Storm)
Duration
4. Abstract Initial and infiltration loss to get the
design rainfall hyetograph
Time from 0 6 12 18 24 30 36 42 48 54 60
start (h)
6-h UH 92 79 64 52 40 30 20 14 10 6 0
ordinates
(m3/s)
The ordinates of a mass curve of rainfall from a
severe storm in a catchment is given.
Ordinates of a 12-h unit hydrograph applicable
to the catchment are also given. Using the
given mass curve, develop a design storm to
estimate the design flood for the catchment.
Take φ index as 0.15 cm/h, estimate the
resulting flood hydrograph. Assume the
base flow to be 50 m3/s.
Tim Cum Incr 12-h First Desig Infiltr Rain DRH DRH DRH DRH Final Base Flood
e Rain eme UH arra n ation fall due of due to due to DRH Flow Hydro
(h) (cm) ntal (m3/ nge Stor loss Exce to 1.7 10.4 16.5 0.2 cm (m3 (m3/s graph
Rain s) men m (cm) ss cm cm cm rain /s) ) (m3/s)
(cm) t(cm (cm) (cm) Rain Rain Rain
)
0 0 0 0 0 0 0 0 0 0 0 0 50 50
36 34 3.5 130 18.3 18.3 1.8 16.5 221 998.4 528 0 1747 50 1797
48 36 2 126 12.2 2 1.8 0.2 214.2 1352 1584 6.4 3157 50 3207
t
• Procedure to
determine the flood
Q
hydrograph at a
point on a
watershed from a
t
known hydrograph
upstream Q
• As the hydrograph
travels, it
– attenuates t
– gets delayed Q
t
Applications of routing techniques:
• Flood predictions
• Evaluation of flood control measures
• Assessment of effects of urbanization
• Flood warning
• Spillway design for dams
• Reservoir design
Routing methods can be broadly classified into
two categories :
i) Hydrologic routing
• Methods combine the continuity equation with
some relationship between storage and
outflow or possibly inflow.
• These relationships are usually assumed, or
empirical.
• An example of such a relationship might be a
stage-discharge relationship.
ii) Hydraulic Routing
= change in storage
Numerical Equivalent
I 1 + I 2 – Q 1 + Q2 S 2 – S1
=
2 2 Δt
The time interval ∆t should be sufficiently
short so that the inflow and outflow
hydrographs can be assumed to be
straight lines in that time interval.
∆t should be shorter than the time of transit
of the flood wave through the reach
Hydrologic Storage Routing or
Reservoir Routing or Level Pool Routing
Advancing
Flood
Wave
I>Q
Receding
Flood
Wave
Q>I
Combine and rearrange
Main Canal
Principal canal of a network of irrigation canals.
Takes off directly from a river or reservoir or from the tail end of a feeder canal.
Usually a large capacity canal which supplies water to branch canals and major distributaries.
Branch Canal
These irrigation canals take off from the main canal on either side.
Generally carry a discharge over 5 cumecs.
Main function is to supply water to major and minor distributaries.
From small ones, direct irrigation is done.
Major Distributary
These are the irrigation canals which take off from the branch canals and sometimes from
the main canal.
Carry a discharge varying from 0.25 to 5 cumec.
Generally used for direct irrigation and supply water through outlets to water courses.
Minor Distributary
Take off from major distributaries and branch canals
Carry discharge less than 0.25 cumec
Water Courses or Field Channels
Small channels which carry water from the outlets of a major or minor distributary or a
branch canal to the fields to be irrigated.
Owned, constructed, controlled and maintained by the cultivators
Based on Canal Alignment
Ridge canal
It is aligned along the ridge or the natural watershed line
When the canal runs on a watershed, it can irrigate areas on both sides and hence a large
area can be brought under cultivation
No drainage can intersect a ridge line, so cross drainage works are not required to be
provided.
Economical
Contour canal
Aligned nearly parallels to the contour of the country
Irrigate areas only on one side
No cross drainage
Side slope canal
Aligned at right angles to the contours of the country
Neither on the watershed nor on the valley but in between the two along the slope
Parallel to the natural drainge
Irrigate areas on only one side
Steep bed slope
Based on Financial Output
Productive canal
When fully developed yields enough revenue to cover up its running cost and a net saving at
the rate of more than six percent of the capital invested initially for its construction
Protective canal
It is constructed as a relief work during famine to provide employment to the inhabitants of
the famine affected area and to protect that area against famine in future
Revenue from this canal usually does not cover its running charges and also does not repay
its initial capital expenditure.
Based on Soil through which it is constructed
Alluvial canals
Constructed through alluvial soils
Readily scoured and deposited
Available in Indo-Gangetic Plains in Northern India
Non alluvial canals
Those canals are constructed through hard soils or disintegrated rocks
Usually found in Central and Southern India
Based on Lining being Provided or not
Unlined canals
An unlined canal is the one which has its bed and banks made up of natural soil through
which it is constructed and not provided with a lining of impervious material
Velocity of flow in these canals is kept low so that bed and banks of the canal may not be
scoured.
Lined canals
Provided with a lining of impervious material on its bed and banks to prevent the seepage
of water
Higher velocities of flow can be allowed in lined canals and hence their cross-sectional areas
can be reduced.
Figure:
Figure:
Figure:
Design of Lined Canal
Introduction
Conservation of water supplies is important as the demand continues to increase and new sources
of supply are becoming scarce. Reduction in the amount of water lost through seepage suring
conveyance in canals is the most important principle of conservation. Seepage losses from unlined
canal system not only entail sizeable wastage of investment made on the construction of dam to
impound water and canal ssytem to distribute it, but also calls for further investment to control
the waterlogging in the irrigated areas.In certain cases, lining may be required purely from
technical considerations. For example, a canal cosntructed partly on rock and partly on permeable
foundation may be unsafe unless it is lined.Sometimes a hard lined surface is required to withstand
the high velocities. Apart from this, the engineer is required to produce a good economic
justification for the capital that is likely to be invested on lining.To justify the lining, economics of
canal lining must be worked out.
Advantages of Lining
Lining of channels reduces the seepage loss thus resulting in a saving of water which can be
utilised for irrigating additional area. The cost of irrigating per unit area is therefore reduced.
It is an important anti waterlogging measure as it reduces the percolation of water to the
ground water reservoir and thus prevents rising of water table.
The lining provides a smooth surface having a low value of rugosity coefficient. Hence a lined
canal has a higher velocity of flow.
Due to higher velocity, greater discharge can pass through a smaller area of cross-section of the
channel.
The higher velocity of flow in a lined canal results in reducing the evaporation loss.
The problem of silting is reduced.
Due to smaller cross-section, there is saving in cost of earthwork, land width and various
structures.
Lining prevents weed growth.
Lining prevalents the water flowing in the channel to come in contact with harmful salts which
might be present in the soil of the bed and banks of the channel
Reduces maintenance cost and possibility of breaching due to increased stability of channel
section.
Provides a stable channel section which is easier to operate.
Disadvantages of Lining
The lining of channel will be economically justified if the extra cost of providing the lining is less
than or equal to the benefits resulting from it.Thus the maximum rate of expenditure on lining
that is economically justifiable can be found out as follows:
Let C = cost of lining in rupees per square meter including the additional cost of dresssing the
banks for lining and accounting for the saving, if any resulting from the smaller cross-sections and
hence smaller area of land, quantity of earthwork and structures required for the lined sections.
This saving will be available on new canals excavated to lined cross-sections from the beginning,
but not on lining of existing unlined canals.
s and S = seepage losses in unlined and lined canals respectively in cubic metres per square meter
of wetted surface per day of 24 hours
p and P = wetted perimeter in metres of unlined and lined canals respectively
T = Total perimeter of lining in metres
d = number of running days of the channel per year
W = value of water saved in rupees per cubic metre
L = Length of the channel in metres
Y = life of lining in years
M = annual saving in rupees in operation and maintenance due to lining
B = annual estimated value in rupees of other benefits for the length of canal under consideration
which include prevention of waterlogging, reduced cost of drainage for adjoining lands etc.
X = percent rate of interest
Economics of Canal Lining
The value of water lost by seepage from the unlined section = pLsdW rupees
The annual saving in value of water otherwise lost by seepage
= (pLsdW −PLSdW )rupees = [LdW (ps − PS)]rupees
Total annual benefits resulting from the lining of channel
(1 + X )Y − 1
NPW = a
X (1 + X )Y
For the lining to be economically feasible, the additional initial cost of the lined canal (TLC)
should be equal to or less than the net present worth (NPW) of the annual benefits, i.e.,
TLC ≤ NPW
Example
An unlined canal giving a seepage loss of 4.0 cumec per million square metres of wetted area is
proposed to be lined with 12 cm thick cement concrete lining which costs Rs. 20.00 per square
meter. Work out the economics of lining and show if the scheme is justified on the basis of the
following data:
Annual revenue from crops = Rs. 4,00,000 per cumec of water
Discharge of the channel = 100 cumec
Area of the channel = 50 m2
Wetted perimeter of the channel = 22.4 m
Wetted perimeter of the lining = 22 m
Annual maintenance cost of unlined channel = 12 paise per m2
Seepage loss in lined channel = 0.01 cumec per million square meter of wetted perimeter
Saving in annual maintenance as a consequence to lining = 40%
Life of lining = 40 years
Rate of interest = 6 percent
Assume additional suitable data, if required.
Solution
Figure:
Concrete lining
The concrete blocks, bricks or stones are laid flat on the compacted sides and
bed of the trapezoidal canal. The joints are filled with cement mortar, which
should have a cement-tosand ratio of 1:3 to 1:4 . A rectangular canal can be
constructed with a concrete or masonry bed and vertical masonry walls. Usually
the water side of the masonry structure is plastered, particularly if the bricks are
not of good quality.
Concrete block, brick or stone masonry lining
Design of Lined Canal
In case of channels lined with hard surfaced materials, two types of sections are used. A triangular
section with circular bottom is adopted for smaller channels and a trapezoidal section with
rounded corners is adopted for large channels. IS: 10430-1982 recommend only trapezoidal section
with rounded corners for all discharges. Typical cross-section of the lined channel is shown in the
figure. The various considerations for the design are given below.
Side Slope
It should be such that they are nearly equal to the angle of repose of the natural soil in the
subgrade so that no earth pressure is exerted over the back of the lining.
Freeboard
Freeboard is measured from the full supply level to the top of the lining.
Figure:
Steps for Design
P = B + 2D(θ + cotθ)
Example
Design a trapezidal shaped concrete lined channel to carry a discharge of 200 cumec at a slope of
30 cm/km. The side slopes of the channel are 1.5:1. The value of N may be taken as 0.017.
Assume limiting velocity in the channel as 2 m/s.
Solution
A = BD + D 2 (θ + cotθ)
or
P = B + 2D(θ + cotθ)
or
36.36 = B + 4.176D (2)
Solution(Contd.)
D = 3.43 m
Then
B = 36.36 − (4.176 × 3.43) = 22.04 m
Thus Bed width = 22.04 m and depth of flow = 3.43 m
Design of Unlined Alluvial Channels
Introduction
A river in alluvium carries huge quantity of sediments along with water. This
sediment comprises soil particles of various sizes ranging from fine silt to coarse
sand. Sediments are either carried in suspension or get deposited on the river bed.
A portion of this sediment is also received by the irrigation canals which draw
water directly from the river. In case of unlined alluvial canals, the quantity of
silt transported may vary from section to section due to scouring of bed and sides
of the canal, as well as silting (deposition of silt) at any section. If the velocity of
flow in a canal is high, then bed and sides of a canal will be scoured. On the
other hand, silting may take place on bed and sides if velocity of flow is low.
Both silting and scouring change the cross section of the canal. Unlined alluvial
canals should be designed for such a velocity that no silting and no scouring
would occur. This velocity is known as non-silting and non-scouring velocity.
Regime Method
Introduction
Cross section of a stable alluvial canal would depend on the flow rate, sediment transport rate and
the sediment size. One method of designing an alluvial canal section is based of the regime
approach in which a set of empirical equations is used. These equations have been obtained by
analyzing the data of stable field channels.
Regime methods for the design of stable channels were first developed by the British engineers
working for canal irrigation in India in the nineteenth century. At that time problem of sediment
deposition was one of the major problems of canal design in India. In order to find a solution to
this problem, some of the British engineers studied the behaviour of such stretches of the existing
canals where the bed was in a state of stable equilibrium.
Regime Method
These stable reaches had not required any sediment clearance for several years of canal operation.
Such canals were called regime canals. These canals generally carried a sediment load smaller that
500 ppm.
Suitable relationships for the velocity of flow in regime canals were evolved. These relationships
are known as regime equations which find acceptance in other parts of the world as well.
The regime relations do not account for the sediment load and hence, should be considered valid
when the sediment load is not large.
Two widely accepted regime method of canal design are Kennedy’s method and Lacey’s method.
Kennedy’s Method
R. G. Kennedy (Executive Engineer, Irrigation Department, Punjab) collected data from 22 canals
of Upper Bari Doab canal system in Punjab. His observations in this canal system led him to
conclude that the sediment in a canal is kept in suspension solely by the vertical component of the
eddies which are generated on the channel bed. In his opinion, the eddies generating on the sides
of the canal had horizontal movements for greater part and therefore did not have sediment
transport power.
This means the sediment supporting power of a canal is proportional to its width and not wetted
perimeter.
Kennedy’s Method
Kennedy’s Equation
On plotting the observed data, Kennedy obtained the following relation, known as Kennedy’s
equation.
V0 = 0.55D0.64
Kennedy termed V0 as the critical velocity (in m/s), defined as the mean velocity which do not
allow scouring or silting in a canal having depth of flow equal to D (in m). This equation is
obviously applicable to such channels which have the same type of sediment as was presented in
the Upper Bari Doab canal system.
On recognizing the effect of the sediment size on the critical velocity, Kennedy modified the above
equation, V = 0.55mD0.64 in which m is the critical velocity ratio, and is equal to V/V0 Here V is
the critical velocity for all sizes of sediment while V0 is the critical velocity for the Upper Bari
Doab sediment.
Kennedy did not try to establish any other relationship for the slope of regime canals in terms of
either the critical velocity or the depth of flow. He suggested the use of Kutter’s equation
alongwith Mannig’s roughness coefficient. The final result do not differ much if one uses the
Manning’s equation instead of Kutter’s equation. Thus the equations
Q = AV (2)
1
V = R 2/3S 1/2 (3)
n
enable one to determine the unknowns B, D and V for Given Q, n and m, if the longitudinal slope
Steps for design
Case 1
Given Discharge, critical velocity ratio, Manning’s n, Bed slope(1H:2V), Side slope
Assume a trial depth D.
Calculate velocity from equation 1.
V = 0.55mD 0.64
A = BD + D 2 /2
Design an irrigation canal carrying a discharge of 30 m3 /s with critical velocity ratio m and
Manning’s n equal to 1.0 and 0.0225 respectively. Assume bed slope as 1 in 5000 and side slopes
as 1H:2V.
Solution
Trial D V0 A B P R V
1 2.0000 0.8571 35.0025 16.5 20.97 1.669 0.8859
2 2.5 0.9887 30.3443 10.89 16.48 1.842 0.9462
3 2.25 0.9242 32.4610 13.3 18.33 1.771 0.9217
4 2.23 0.9189 32.6470 13.52 18.51 1.764 0.9192
5 2.22 0.9163 32.7411 13.64 18.6 1.76 0.9180
6 2.221 0.9165 32.7316 13.63 18.59 1.76 0.9181
7 2.222 0.9168 32.7222 13.62 18.58 1.761 0.9182
8 2.223 0.9171 32.7128 13.6 18.57 1.761 0.9184
9 2.224 0.9173 32.7034 13.59 18.57 1.761 0.9185
Figure:
Steps
Case 2
Given Discharge, critical velocity ratio, Manning’s n, B/D ratio.
Given B/D = c
Express A like this
A = BD + D 2 /2
A = cD 2 + D 2 /2 = D 2 (c + 0.5)
V = 0.55mD 0.64
Q = AV
or
Q = D 2 (c + 0.5)(0.55mD 0.64)
D = (Q/(0.55m(c + 0.5))1/2.64
Calculate B as B = cD
Calculate Hydraulic Radius
Example 2
Design an irrigation canal to carry a discharge of 5 m3/s. Assume Manning’s n = 0.0225, critical
velocity ratio m=1 and B/D = 3.24.
Solution
B/D = 3.24
B = 3.24D = 4.536m
√
P = B + D 5 = 7.6664m
R = A/P = 0.9561
n2 V 2
S = 4/3 = 1 in 4089
R
Lacey’s Theory
Gerald Lacey carried out a detailed study of the problem of designing stable channels in alluvium
and made significant contribution. He developed the regime theory and evolved a number of
equations on the basis of his own observations and the observations of the earlier investigators.
According to Lacey the width, depth and bed slope of a channel constructed in easily erodible
alluvium and carrying a fixed discharge and silt load are uniquely determined by nature. In other
words, for a channel constructed in alluvium to carry discharge, the width, depth and bed slope
will undergo modification by silting and scouring till equilibrium is attained. The channel is termed
as regime channel.
Lacey’s regime equations
Lacey considered hydraulic mean radius R as a more pertinent variable rather than depth of flow D
considered by Keneddy. Further as a measure of silt grade Lacey introduced a silt factor f with its
value equal to unity for the silt of channels of Upper Bari Doab Canal System. Lacey obtained a
relationship which is as follows
r
2
V = fR (4)
5
where V = mean velocity of flow in m/s
f = silt factor
R = hydraulic mean radius in m
This was the first regime equation given by Lacey.
For determining the channel dimension second equation would be necessary giving either
cross-sectional area A or perimeter P. Lacey on the basis of plotted data evolved a relationship
between A and V which is as follows
Af 2 = 140V 5 (5)
Lacey’s regime equation
According to Lacey, the silt factor f is dependent on the mean particle size of the boundary
material in the channel. The value of f may be determined by the following relation given by Lacey
√
f = 1.76 m (6)
where m = mean particle size in mm = d50in mm
Lacey’s regime equation
V-Q-f Relation
or
Qf 2 = 140V 6 (7)
or 1/6
Qf 2
V = (8)
140
where Q = discharge in cumec
Lacey’s regime equation
P-Q Relation
By combining equations 4 and equation 7, a relationship between P and Q can be obtained as follows.
Raising both sides of the equation 4 to the fourth power, we get
4 2 2
V4 = f R
25
Eliminating f2 between this equation and equation (7), we get
25V 4
Q = 140V 6
4R 2
or
25Q
= 140V 2
4R 2
Since
R = A/P
25Q
= 140V 2
4(A/P) 2
or
Q 560 2
= V
(A/P)2 25
or
560
P 2Q = (AV )2
25
or
P 2 = 22.4Q
or
p
P = 4.75 Q (9)
Regime flow equation
After determining the dimensions of a regime channel, the bed slope of the channel is required to
be determined. Thus by plotting a large number of data Lacey obtained the following flow
equation
S-f-R Relation
V 3 = 1260 R 2 S
From equation (4), we get
For the design of channel by Lacey’s method, the following equations are used:
√
f = 1.76 m
1/6
Qf 2
V =
140
Q = AV
√
P = 4.75 Q
r
2
V = fR
5
f 5/3
S=
3340Q 1/6
The items which must be known are Discharge Q, Silt factor f for mean particle size m
Steps for Design
Calculate the silt factor f from the 1st Equation if mean particle size is given
Calculate the velocity from the equation
1/6
Qf 2
V =
140
D2
A = BD +
2
√
P=B+ 5D
Both the values of R obtained in the above two steps should be same.
Calculate Bed Slope S from the following relation
f 5/3
S=
3340Q 1/6
Example
Design an irrigation canal in alluvial soil according to Lacey’s silt theory for the following data.
Full supply discharge = 10 cumec
Lacey’s silt factor = 0.9
Side slopes of the channel = 21 (H):1(V)
Solution h 2 i1/6 h i1/6
Qf 10×0.92
Velocity V = 140 = 140 = 0.62 m/s
Area A = VQ = 0.62
10
= 16.13m2
2 (0.62)2
Hydraulic mean radius R = 25 Vf = 5
2 × 0.9 = 1.07 m
Perimeter √P = RA = 16.13
1.07 = 15.07m
P = B +√D 5
orB + D 5 = 15.07
A = BD + D 2 /2
or
BD + D 2/2 = 16.13
From the two equations, we get
f 5/3 (0.9)5/3 1
Bed SlopeS = 3340(Q)1/6
= 3340(10)1/6
= 5844
Canal Irrigation
Planning of an irrigation canal project includes the determination of: (i) canal alignment,
and (ii) the water demand. The first step in the planning of an irrigation canal project is to
carry out a preliminary survey to establish the feasibility of a proposal. Once the feasibility
of the proposal has been established, a detailed survey of the area is carried out and,the
alignment of the canal is fixed. The water demand of the canal is, then, worked out.
The information on the following features of the area are to be collected: (i) Type of soil,
(ii) Topography of the area, (iii) Crops of the area, (iv) Rainfall in the area, (v) Water
table elevations in the area, (vi) Existing irrigation facilities, and (vii) General outlook of
the cultivators with respect to cultivation and irrigation
The preparation of plans for a large canal project is simplified in a developed area because
of the availability of settlement maps (also called shajra maps having scale of 16 inches
to a mile i.e., 1/3960 1/4000) and revenue records in respect of each of the villages of
the area. The settlement maps show the boundaries and assigned numbers of all the fields
of the area, location of residential areas, culturable and barren land, wells, ponds, and
other features of the area. Usually for every village there is one settlement (or shajra)
map. These maps and the revenue records together give information on total land area,
cultivated area, crop-wise cultivated area and the area irrigated by the existing ponds and
wells. With the help of settlement maps of all the villages in a doab, a drawing indicating
distinguishing features, such as courses of well-defined drainages of the area, is prepared.
On this drawing are then marked the contours and other topographical details not available
on the settlement maps but required for the planning of a canal irrigation project. Contours
are marked after carrying out levelling survey of the area. The details obtained from the
settlement maps should also be updated in respect of developments such as new roads,
additional cultivated area due to dried-up ponds, and so on.
Seepage water lost from the canals percolate deep into the ground and is the main cause
of the loss of water carried by a canal. In addition, some canal water is also lost due to
evaporation. The loss due to evaporation is about 10 per cent of the quantity lost due to
seepage. The seepage loss varies with the type of the material through which the canal
runs. Obviously, the loss is greater in coarse sand and gravel, less in loam, and still less
in clay soil. Between the headworks of a canal and the watercourses, the loss of water on
account of seepage and evaporation is also known as conveyance loss, and expressed as
m3 /s per million square metres of either wetted perimeter or the exposed water surface
area. This loss may vary from 20 to 50 per cent of water diverted at the headworks
depending upon the type of soil through which canal runs and the climatic conditions of
the region. The following empirical relation has been found to give comparable results
2
ql = (1/1200)(B + h) 3 (1)
In this relation, ql is the loss expressed in m3/s per kilometre length of canal and B and h
are, respectively, canal bed width and depth of flow in metres
The canal alignment does not remain straight all through the length of the canal, and
curves or bends have to be provided. There is a slight drop in the water surface at the
inner edge of the curve and a slight rise at the outer edge of the curve. This results
in slight increase in the velocity at the inner edge and slight decrease in the velocity at
the outer edge. As a result of this, the low-velocity fluid particles near the bed move to
the inner bank and the high-velocity fluid particles near the surface gradually cross to the
outer bank. The cross currents tend to cause erosion along the outer bank. Therefore,
wherever possible, curves in channels excavated through loose soil should be avoided. If it
is unavoidable, the curves should have a long radius of curvature.
The amount of water needed for the growth of a crop during its entire crop-growing period
is known as the water requirement of the crop, and is measured in terms of depth of water
spread over the irrigated area. This requirement varies at different stages of the growth of
the plant. The peak requirement must be obtained for the period of the keenest demand.
One of the methods to decide the water requirement is on the basis of kor watering.
When the plant is only a few centimetres high, it must be given its first watering, called
the kor watering, in a limited period of time which is known as the kor period. If the plants
do not receive water during the kor period, their growth is retarded and the crop yield
reduces considerably. The kor watering depth and the kor period vary depending upon the
crop and the climatic factors of the region. In UP, the kor watering depth for wheat is
13.5 cm and the kor period varies from 8 weeks in north-east UP (a relatively dry region)
to 3 weeks in the hilly region (which is relatively humid). For rice, the kor watering depth
is 19 cm and the kor period varies from 2 to 3 weeks.
∆ = 8.64b/D(meters) (2)
For the purpose of designing on the basis of the keenest demand (i.e., the kor period
requirement) the base period b and the water depth ∆ are replaced by the kor period and
kor water depth, respectively.
Since the water demands for wheat and rice are at different times, these are not cumulative.
Therefore, the distributary channel should be designed for the larger of the two discharges,
viz., 1.7 m3/s. The above calculations exclude channel losses and the water requirement
of other major crops during their kor period.
When the canal water has reached near the fields to be irrigated, it has to be transferred
to the watercourses. An outlet is a masonry structure through which water is admitted
from the distributary into a watercourse. It also acts as a discharge measuring device. The
discharge though an outlet is usually less than 0.085 m3 /s.If the total available supply is
insufficient, the outlets must be such that equitable distribution can be ensured.
Types of Outlet Canal outlets are of the following three types: (i) Non-modular outlets,
(ii) Semi-modular outlets, and (iii) Modular outlets.
Non-modular outlets The non-modular outlets are very suitable for low head conditions
as the loss of head in this outlet is less than that in a modular outlet. The discharge in
non-modular outlets capacity depends on the difference of water levels in the distributary
and the watercourse. The discharge through these outlets fluctuates over a wide range
with variations in the water levels of either the distributary or the watercourse. Such an
outlet is controlled by a shutter at its upstream end.
Modular outlets Modular outlets are those whose discharge is independent of the water
levels in the distributary and watercourse. A modular outlet supplies fixed discharge and,
therefore, enables the farmer to plan his irrigation accordingly. However, in case of excess
or deficient supplies in the distributary, the tail-end reach of the distributary may either
get flooded or be deprived of water. This is due to the reason that the modular outlet
would not adjust its discharge corresponding to the water level in the distributary. But, if
an outlet is to be provided in a branch canal which is likely to run with large fluctuations
in discharge, a modular outlet would be an ideal choice.
Q = C1 hn (4)
dQ dh
=n (5)
Q h
Similarly, the discharge Q0 through the outlet can be expressed as
Q0 = C1 H m (6)
dQ0 dH
=m (7)
Q0 H
Here, m and n are suitable indices and C1 and C2 are constants.
Canal Irrigation March 16, 2021 14 / 31
Flexibility
Thus,
m h dH
F = x x (8)
n H dh
For semi-modular outlets, the change in the head dH at an outlet would be equal to the
change in the depth of flow dh in the distributary. Therefore,
m h
F = x (9)
n H
If the value of F is unity, the rate of change of outlet discharge equals that of the distributary
discharge. For a modular outlet, the flexibility is equal to zero. Depending upon the value
of F, the outlets can be classified as: (i) proportional outlets (F = 1), (ii) hyperproportional
outlets (F > 1), and (iii) subproportional outlets (F < 1). When a certain change in the
distributary discharge causes a proportionate change in the outlet discharge, the outlet (or
semi-module) is said to be proportional. A proportional semi-module ensures proportionate
distribution of water when the distributary discharge cannot be kept constant. For a
proportional semi-modular outlet (F = 1),
m h
= (10)
n H
The ratio (H/h) is a measure of the location of the outlet and is termed setting.
Canal Irrigation March 16, 2021 15 / 31
Flexibility
The setting for a proportional outlet is equal to the ratio of the outlet and the channel
indices.
For hyper-proportional and sub-proportional outlets the setting must be, respec-
tively,less and more thanm/n.
For a wide trapezoidal (or rectangular) channel,n can be approximately taken as 5/3
and for an orifice type outlet, m can be taken as 1/2. Thus, an orifice-type module
will be proportional if the setting (H/h) is equal to (1/2)/(5/3), i.e., 0.3.
The module will be hyper-proportional if the setting is less than 0.3 and sub-proportional
if the setting is greater than 0.3.
Similarly, a free flow weir type outlet (m = 3/2) would be proportional when the
setting equals 0.9 which means that the outlet is fixed at 0.9 h below the water
surface in the distributary.
Sensitivity The ratio of the rate of change of discharge (dQ0 /Q0 ) of an outlet to the rate
of change in the water surface level of the distributary channel with respect to the depth
of flow in the channel is called the sensitivity of the outlet. Thus,
dQ0 /Q0
S= (11)
dG /h
dQ0 /Q0
F = (12)
dQ/Q
dh
= (dQ0 /Q0 )/n( ) (13)
h
1
= S (14)
n
S = nF (15)
The non-modular outlet is usually in the form of a submerged pipe outlet or a masonry
sluice‘which is fixed in the canal bank at right angles to the direction of flow in the
distributary.
The diameter of the pipe varies from 10 to 30 cm. The pipe is laid on a light concrete
foundation to avoid uneven settlement of the pipe and consequent leakage problems.
The pipe inlet is generally kept about 25 cm below the water level in the distributary.
When considerable fluctuation in the distributary water level is anticipated, the inlet
is so fixed that it is below the minimum water level in the distributary.
Pipe outlet
If H is the difference in water levels of the distributary and the watercourse then the
discharge Q through the outlet can be obtained from the equation,
V2 fL
H= [0.5 + + 1] (16)
2g D
where
Q p d
V = 2
= 2gH( )1/2 (17)
(π/4)d 1.5d + fL
where, d = diameter of pipe outlet L = length of pipe outlet and f = friction factor for
pipe. Alternatively, the discharge Q can be expressed as
Q = AV (18)
π 2p d
= d 2gH( )1/2 (19)
4 1.5d + fL
p
Q = CA 2gH (20)
d
C =( )1/2 (21)
1.5d + fL
The simplest type of semi-modular outlet is a pipe outlet discharging freely into the
atmosphere.
The pipe outlet, described as the non-modular outlet, works as semi-module when it
discharges freely into the watercourse.
The exit end of the pipe is placed higher than the water level in the watercourse.
In this case, the working head H is the difference between the water level in the
distributary and the centre of the pipe outlet.
The efficiency of the pipe outlet is high and its sediment conduction is also good.
The discharge through the pipe outlet cannot be increased by the cultivator by digging
the watercourse and thus lowering the water level of the watercourse.
Usually, a pipe outlet is set so that it behaves as subproportional outlet, i.e., its setting
is kept less than 0.3.
It mainly consists of an orifice with bellmouth entry, a long expanding delivery pipe,
and an intervening vertical air column above the throat
The air vent pipe permits free circulation of air around the jet. This arrangement
makes the discharge through the outlet independent of the water level in the water-
course.
The water jet enters the cast iron expanding pipe which is about 3 m long and at the
end of which a cement concrete pipe extension is generally provided. Water is then
discharged into the watercourse.
This outlet can be easily tampered with by the cultivator who blocks the air vent pipe
to increase the discharge through the outlet.
An open flume outlet is a weir with a sufficiently constricted throat to ensure super-
critical flow and long enough to ensure that the controlling section remains within the
throat at all discharges up to the maximum.
A gradual expansion is provided downstream of the throat. The entire structure is
built in brick masonry but the controlling section is generally provided with cast iron
or steel bed and check plates.
This arrangement ensures the formation of hydraulic jump and, hence, the outlet
discharge remains independent of the water level in the watercourse.
The discharge through the outlet is proportional to H3/2. The efficiency of the outlet
varies between 80 and 90 per cent.
For the range of outlet discharges normally used, the outlet is either deep and narrow,
or shallow and wide. While a narrow outlet gets easily blocked, a shallow outlet is not
able to draw its fair share of sediment.
Most of the modular outlets have moving parts which make them costly to install as well
as maintain. The following two types of modular outlets (also known as rigid modules),
however, do not have any moving part:
1 Gibbs rigid module, and
2 Khannas rigid module.
Gibbs rigid module This module has an inlet pipe under the distributary bank. This
pipe takes water from distributary to a rising spiral pipe which joins the eddy chamber
. This arrangement results in free vortex motion. Due to this free vortex motion,
there is heading up of water near the outer wall of the rising pipe. The water surface
thus slopes towards the inner wall. A number of baffle plates of suitable size are
suspended from the roof of the eddy chamber such that the lower ends of these plates
slope against the flow direction. With the increase in head, the water banks up at
the outer wall of the eddy chamber and impinges against the baffles and spins round
in the compartment between two successive baffle plates. This causes dissipation of
excess energy and results in constant discharge. The outlet is relatively more costly
and its sediment withdrawal is also not good.
An irrigation canal takes its supplies from a river or a stream. In order to divert water from the
river into the canal, it is necessary to construct certain works or structures across the river and at
the head of the offtaking canal. These are known as canal headworks. THe canal headworks may
be classified into the following two types.
Storage headworks
Consists of a dam constructed across the river to create a reservoir in which water is stored
during the period of excess flow in the river.
Diversion headworks
It serves to raise the water level in the river and divert the required quantity into the canal.
Command area is increased.
Regulates the supply of water into the canal.
Controls the entry of silt into the canal.
Provides some storage of water for a short period.
Reduces the fluctuations in the level of supply in the river.
Types of diversion headworks
River
Divide Wall
Head Regulator
Approach Channel
Canal
Weir
Scouring Sluices
Fish Ladder
The undersluices are openings provided in the weir wall with their crest at a low level. They are
located on the same as offtaking canal.
Functions of undersluices
They preserve a clear and well-defined channel towards the head regulator.
They scour the silt deposited on the river bed in the pocket upstream of the canal head
regulator.
They pass low floods without the necessity of dropping the weir crest shutters
Divide Wall
A divide wall is a long masonry wall which is constructed at right angles to the axis of the weir to
seperate the undersluices from the rest of the weir. The top width of divide wall is about 1.5 to
2.5m.
Fish ladder
Large rivers have many types of fishes. They move from one part to another according to season.
In our country generally fish move from upstream to downstream in search of warmth and return
upstream before monsoon for clear water. Due to the construction of a weir across the river, such
migration will be obstructed and large scale destruction of fish life may take place in the river.
Thus to enable the fish to migrate, fish ladder is provided. Since most fish can travel upstream
only if the velocity of flow does not exceed 3 to 3.5 m/s, the design of the fish ladder should be
such that it constantly releases water at a velocity not exceeding this value.
Canal Head regulator
A canal head regulator is a structure constructed at the head of the canal taking
off from the upstream of a weir or barrage. It consists of a number of spans
supported by piers which support the gates provided for regulation of flow into
the canal.
Functions
It regulates the supply of water into the canal.
It controls te entry of silt into the canal.
It excludes the high flood from entering the canal.
Weir
A weir is a structure constructed across a river to raise its water level and divert the water into the
canal. On the crest of the weirs usually shutters are provided so that part of the raising up of
water (or ponding) is carried out by shutters. During floods, the shutters may be dropped down to
allow water to flow over the crest of the weir. Weirs are usually aligned at right angles to the
direction of flow of the river.
Weirs may be classified according to the material of construction and certain design features into
the following three types
Masonry weirs with vertical drop weir
Rockfill weirs
Concrete weirs
Pond Level
Crest Shutter
Weir Wall
Block Protection
Inverted Filter
Impervious horizontal floor or apron and a masonry weir wall with with either both faces
inclined or both faces vertical or one face vertical and other face inclined
At the u/s and d/s ends, piles are provided
At the u/s end,a block protection and at the d/s end, a graded filter is provided
Launching aprons or pervious aprons are provided after the block protection or filter
Suitable for any type of foundation
Floor design based on Bligh’s theory
Barrage
The crest is kept at a low level and the raising up of water level is accomplished mainly by means
of gates. During floods, these gates can be raised clear off the high flood level and thus enable the
high flood to pass with mimimum of afflux. It provides a better control on the water level in the
river but is comparatively more costly.Design involves same procedure as a concrte weir.
High Flood Level
Pond Level
GATE
Inverted Filter
Launching Apron
U/S PILE
D/S PILE
Figure:
Causes of failure if Weirs on permeable foundation and their remedies
The various causes of failure of weirs on permeable foundations may be classified broadly into two
categories:
Due to seepage flow
Due to surface flow
The seepage flow may cause the failure of a weir in two ways
By piping
If the water percolating through the foundation has sufficient force when it emerges at the
downstream end of the impervious floor it may lift up the soil particles at the end of the floor
With the removal of surface soil, there is further concentration of flow into the resulting
depression and more soil is removed.
This process of erosion progressively extends backwards towards the u/s side and results in
the removal of soil and developing pipe like formation beneath the floor.
The floor may subside in the hollows so formed and fail which is known as failure due to
piping.
The following measures may be taken to prevent the failure due to piping
Providing sufficient length of the impervious floor so that the path of percolation is increased
and exit gradient(hydraulic gradient) is reduced.
Providing piles at the upstream and downstream ends of the impervious floor.
Causes of failure if Weirs on permeable foundation and their remedies
By Uplift Pressure
The water percolating through the foundation exerts an upward pressure on the imprevious
floor. This pressure is known as uplift pressure. If the uplift pressure is not counterbalanced by
the weight of the floor, it may fail by rupture.
The following measure may be taken to prevent the failure of the floor by rupture due to excessive
uplift pressure.
Providing sufficient thickness of the impervious floor.
Providing pile at the u/s end of the impervious floor so that the uplift pressure is reduced in
the d/s side.
Barrier
(H/L)2d1
(H/L)2d2
Subsoil hydraulic
gradient Line H
(H/L)2d3
d2
d1
d3
Bligh assumed that the percolating water creeps along the base profile of the structure which is in
contact with the subsoil. The length of path thus traversed by the percolating water is called the
length of creep. He further assumed that the heas loss per unit length of creep which is called
hydraulic gradient is constant throughout the percolating passage which means loss of head is
proportional to the length of creep. He made no distinction between vertical and horizontal creep.
In the figure, a barrier impounding water of depth H is provided with a horizontal floor of length b
with three vertical cutoffs (sheet piles). The percolating water then follows the path indicated by
arrows and the creep length L will be given by
To ensure the safety of the impervious floor against the two possible ways in which failure may be
caused by subsurface flow, following criteria are required to be satisfied.
Safety against piping
The length of creep should be sufficient to provide a safe hydraulic gradient according to the
type of soil. He has recommended certain values of C for different soils as given in Table.
According to Bligh, if the hydraulic gradient ≤ (1/C ), there will be no danger of piping.
The seepage head H is to be measured from the water level upstream to the corresponding
lowest water level downstream.
Type of soil Value of C Safe hydraulic gradient
Light sand and mud 18 1/18
Fine mica sand 15 1/15
Coarse grained sand 12 1/12
Sand mixed with boulder and gravel 5 to 9 1/9 to 1/5
Bligh’s Creep Theory
wh0 = (wG )t
h0
t= (3)
G
h0 = tG
∴ h0 − t = tG − t = t(G − 1)
h0 − t
∴t=
G −1
h
or t = (4)
G −1
The floor thickness given by the above equation is usually incrased by considering a factor of safety of (4/3) and hence
4 h
t= (5)
3 G −1
Example
Subsoil hydraulic
5m gradient Line
A B C
5m 3m
5m
8m 10 m
10 m 8m
25 m
Figure:
The figure shows the section of a hydraulic structure on permeable foundation. Calculate the
average hydraulic gradient according to Bligh’s creep theory. Also find the uplift pressures at
points A,B and C as shown in the figure and the floor thickness required at these points.
Example
Solution
According to Bligh’s creep theory, the total creep length is given by
L = (2 × 5) + (2 × 3) + (2 × 10) + 25 = 61 m
5 1
∴Hydraulic gradient = 61 = 12.2
The structure would be safe on coarse grained sand.
Uplift pressure at point A
Length of creep upto point A = (2×5)+5 = 15
Residual seepage head at point A
15
h1 = 5(1 − ) = 3.77 m
61
4 h 4 3.77
t1 = = = 4.05m
3 G − 1 3 2.24 − 1
4 h 4 2.87
t1 = = = 3.09m
3 G − 1 3 2.24 − 1
4 h 4 2.21
t1 = = = 2.38m
3 G − 1 3 2.24 − 1
Khosla’s Theory
A number of important structures were designed on the basis of this theory, someof which
remained stable while others gave trouble or failed. Khosla carried out investigation and led to the
following conclusions.
The outer faces of the end sheet piles were much effective than the inner ones and horizontal
length of the floor.
The intermediate sheet piles, if smaller in length than the outer ones were ineffective except for
local redistrbution of pressures.
Undermining of floors started from the tail end. If the hydraulic gradient at the exit was more
than the critical gradient for the particular soil, the soil particles would move with the flow of
water thus causing progressive degradation of the subsoil, resulting in cavities and ultimate
failure.
It was absolutely essential to have a reasonably deep vertical cutoff at the d/s end of the floor
to prevent undermining.
Khosla’s method of independent Variables
W.S. W.S.
H H
E1 C1 E C
d
d
b b
D1
D
Pile at Upstream end Pile at Downstream end
W.S.
E C
b1
d
Intermediate Pile
Figure:
The above cases have been analysed by Khosla for determining the residual seepage head or the uplift pressure head at
the key points and the exit gradient. The key points are the junction points of pile and floor, the bottom point of pile
and the bottom corners of depressed floor. The results of this analysis are presented in the form of curves from which the
values of φi.e., the ratio of the residual seepage head and the total seepage head can be determined at the key points.
Pile at Upstream end
100 λ−2
φC1 = 100 − cos−1
π λ
100 λ−1
φD1 = 100 − cos−1
π λ
where
√
1+ 1 + α2
λ=
2
b
α=
d
100 λ−2
φE = cos−1
π λ
100 λ−1
φD = cos−1
π λ
Intermediate Pile
100 λ1 − 1
φE = cos−1
π λ
100 λ1
φD = cos−1
π λ
100 λ1 + 1
φC = cos−1
π λ
b
α1 =
d
b − b1
α2 =
d
q q
λ1 = 1 + α12 − 1 + α22 /2
q q
λ1 = 1 + α12 + 1 + α22 /2
Corrections
Figure:
Ec E c t2 t3 E c
t1
c1 E1 c1 E1
Flow Flow
d3
d1 d2
Flow
D D D
In the standard forms with vertical cutoffs, the thickness of the floor is assumed to be negligible. So the pressures at
the junction points E and C pertain to the level at the top of the floor whereas the actual junction is with the bottom
of the floor. The pressures at the actual junction points E1and C1are interpolated by assuming a linear variation from
the hypothetical point E to D and also from D to C.
φD − φC
φC1 = φC + × t1
d1
φE − φD
φE1 = φE − × t2
d2
φD − φC
φC1 = φC + × t2
d2
For pile at downstream end
Correction for E1= φE d−φ
3
D
× t3 (subtractive)
∴Pressure at E1
φE − φD
φE1 = φE − × t3
d3
Correction for mutual interferecence of piles
EC
Pile 1 d EC
EC
bs
Pile 2
D D
b’
Pile 3
D
D
b
A suitable correction is to be applied for a sloping floor. The values of the correction in percent of
pressure are given in the following table.
Slope (Vertical to Horizontal) Correction % of pressure
1 in 1 11.2
1 in 2 6.5
1 in 3 4.5
1 in 4 3.3
1 in 5 2.8
1 in 6 2.5
1 in 7 2.3
1 in 8 2.0
The correction given by this table is to be multiplied by the horizontal length of the sloping floor
and divided by th distance between the two pile lines between which the sloping floor is located.
The correction is applicable only to those key points of the pile lines at the beginning or the end
of the sloping floor which are facing the sloping floor. So in the above figure, the correction for
slope is to be applied only to point E of pile no. 2 and the correction will be obtained multiplying
the appropriate value from the table by (bs /b’). The correction is additive for a downward sloping
floor and subtractive for an upward sloping floor in the direction of flow.
The uplift pressure distribution along the floor may be obtained by assuming a linear variation
between the key points after obtaining pressures at the key points.
Exit Gradient
Defined as the hydraulic or pressure gradient of subsoil flow at the downstream or the exit end
of the floor.
For a length of floor of length b with a vertical cutoff of depth d at its downstream end, Khosla
and his associates derived an expression for the exit gradient GE as follows:
H 1
GE = √
dπ λ
where H = total seepage head
√
1 + 1 + α2
λ=
2
and α = db
From the above equation, it can be seen that if no cutoff is provided at the downstream end of
the floor, i.e., d = 0, the exit gradient GE is infinite. But it is evident that if there is no cutoff at
the d/s end, a higher exit gradient will exist will exist which may lead to the failure of the floor
due to piping. It is essential to provide a vertical cutoff at the d/s end of the floor to reduce the
exit gradient. The exit gradient should always be less than critical hydraulic gradient which is
defined as the gradient at which the soil particles will be lifted up and lead to undermining. The
permissible exit gradients for different soils are given in the following table.
Using Khosla’s method, obtain the uplift pressures for the weir profile shown in the figure. Also
calculate the exit gradient and find whether it is safe against piping if it is founded on fine sand
with permissible exit gradient 1/6.
Solution
100 λ−1
φD1 = 100 − cos−1 = 77.418%
π λ
where
√
1+ 1 + α2
λ= = 4.145
2
b 65
α= = = 7.222
d 100 − 91
d = 98 − 91 = 7m
D = 98 − 89.5 = 8.5m
0
b = 16 + 32 = 48m
Solution(Contd.)
r
8.5 7 + 8.5
C = 19 = 1.9%(additive)
48 65
∴Corrected value of φC1 = 67.313+2.24+1.9 = 71.453%
Pile at Downstream end
100 λ−2
φE3 = cos−1 = 27.976%
π λ
100 λ−1
φ D3 = cos−1 = 19.452%
π λ
where
√
1+ 1 + α2
λ= = 5.525
2
b 65
α= = = 10
d 6.5
d = 94 − 89.5 = 4.5m
d = 94 − 89.5 = 4.5m
r
4.5 4.5 + 4.5
C = 19 = 1.395%(subtractive)
16 65
∴Corrected value of φE3 = 27.976 − 2.627 − 1.395 = 23.954%
Solution(Contd.)
Intermediate pile
b1 48.5
α1 = = = 7.461
d 6.5
b − b1 65 − 48.5
α2 = = = 2.538
d 6.5
q q
λ1 = 1+ α12 − 1+ α22 /2 = 2.401
q q
2 2
λ= 1 + α1 + 1 + α2 /2 = 5.128
100 λ1 − 1
φE2 = cos−1 = 41.19%
π λ
100 λ1
φ D2 = cos−1 = 34.49%
π λ
100 λ1 + 1
φC2 = cos−1 = 26.91%
π λ
Solution(Contd.)
D = 93 − 91 = 2m
d = 93 − 89.5 = 3.5m
Correction for interference of d/s pile at φC2
r r
D d +D 3.5 3.5 + 3.5
C = 19 = 19 = 1%(+)
b0 b 16 65
d = 93 − 89.5 = 3.5m
D = 93 − 89.5 = 3.5m
Correction for slope
Slope = 97−93
32
=1 in 8
From the table, for a slope of 1 in 8, correction in % of pressure = 2%
32
∴Correction = 2× bbs0 =2× 48 = 1.3%(+)
∴Corrected value of φE2 = 41.19 − 3.1 − 0.3 + 1.3 = 39.1%
∴Corrected value of φC2 = 26.91 + 3.5 + 1 = 31.41%
Solution(Contd.)
Exit Gradient
H 1
GE = √
dπ λ
where H = total seepage head
H = 103 -96 = 7m
d = 96 - 89.5 = 6.5 m
b 65
α= = = 10
d 6.5
√
1+ 1 + α2
λ= = 5.525
2
7 1 1
∴ GE = × √ =
6.5 π 5.525 6.856
Permissible gradient = 16
Hence safe against piping.
Land Drainage
I Layout - Aligned along the path of the natural drainage and should be provided with an outfall,
either into a bigger drain or natural stream
I Capacity - Drains are provided to carry storm runoffs, seepage water from subsoil and surface
flow resulting from excess irrigation of fields. Since storm runoffs are usually far in excess of
the other flows, the drains are designed for storm runoff only. A judicious estimate of surface
runoff should be made to obtain the design capacity of drains.
I Section - The drains have trapezoidal cross-section and are designed on the same principle of
irrigation channels. Since deep drains receive both seepage and storm water, a small drain in
the bed of larger drain is usually provided to carry the small seepage discharge. With this
arrangement, the full section of the drain would be operativeonly during te storm. This would
reduce the problem of maintenance because the seepage water which would otherwise have
flown at small depth and velocity over the entire width leading to weed growth over the entire
section would flow only thorugh the small drain and so only this sction would be required to
maintained. The depth and bottom width of deep drains generally varies from 0.6 to 1.5 m and
0.3 to 1.2m respectively.
I Maintenance - The main problem of maintenance is to keep from weeds. It may be retarded by
having larger velocities of flow in the drains.
Drawbacks of open drains
Closed drains are useful for prevention of waterlogging of land as well as for
relieving the land already waterlogged. These are tile drains usually of porous
earthen ware which are laid below the ground level butting each other with open
joints and covered up by earth. Thus these drains do not put any area out of
cultivation and do not obstruct any farming operations. The usual spacing of
these drains is from 15 to 45 m, closer spacing being used for soils of low
permeability. The usual diameter of tile drains is 100 mm. The drains are usually
located about 0.3m lower than the desired water table. These drains are laid at
at a gradient steeper than 1 in 500. As far as possible, the closed drains should
be placed in a permeable stratum to ensure effective drainge. The closed drains
have their outlets in natural or artificial channels.
Closed drains are designed to carry only the seepage water and their capacity is
determined by the rate of infiltration which can be estimated theoretically if the
soil permeabilities are known.
Layout of closed drains system
I Double main system - Two seperate mains/submains with seperate laterals are
provided. The system is provided when the bottom of the depression is wide.
This arrangement helps to reduce the length of laterals.
I Intercepting drain system - In this system, there are no laterals but only a
main/submain is provided at the toe of the slope. The arrangement is adopted
when the hilly land is to be drained.
Layout of closed drains system
Discharge and spacing of closed drains
Figure:
With reference to the figure, let L be the spacing of the drains and a be their height above the
impervious stratum. Let b be the maximum height of the drained water table above the
impervious stratum. At a distance x from the centre of drain, let the drained water table be at a
height y above the impervious stratum. The various assumptions are as follows:
i) The hydraulic gradient at a distance x from the centre of the drain is (dy/dx).
ii) Flow lines are parallel and area of flow section at a distance x from the centre of the drain per
unit length of drain is y×1=y
iii) Discharge towards the drain is inversely proportional to the distance from the drain.
Discharge and spacing of closed drains
dy
Qy = ky (1)
dx
where Qy is the discharge passing throguh the section at a distance x from the centre of drain per
unit length of the drain.
If QD is the total discharge per unit length carried by the drain, then (QD /2)enters the drain from
either side. Since discharge is inversely proportional to the distance from the drain, discharge = 0
when x=(L/2) and discharge = (QD /2)when x=0.
Thus
QD (L/2) − x
Qy =
2 (L/2)
or
QD
Qy = (L − 2x) (2)
2L
Discharge and spacing of closed drains
or
QD
(L − 2x)dx = ydy
2kL
Integrating we get,
QD y2
(Lx − x 2 ) = +C
2kL 2
Lk(y 2 − a2 )
QD =
Lx − x 2
4k(b2 − a2 )
∴ QD = (3)
L
and
4k(b2 − a2 )
L= (4)
QD
Discharge and spacing of closed drains
From equations 3 and 4, discharge and spacing of closed drains can be determined. For a given
drainage system, equation 3 can be used to predict the inflow. However, to determine the spacing
for new drainage system, by equation 4 we can get if QD is known which will depend on infiltration
discharge into the ground which should be removed by drain. Different values have been
suggested for QD out of which the commonly adopted value being 1% of the average annual rainfall
drained in 24 hours.
Example 1
In a system of closed drains, the drains are placed with their centres 7.7 m above the impervious
stratum and the maximum height of the drained water table above the impervious stratum is 8m.
If the spacing between the drains is 25m and discharge entering the drain per unit length is
4.5×10−6 cumec/m, find the coefficient of permeability of soil.
Solution
4k(b2 − a2 )
L=
QD
or
LQD
k=
4(b2 − a2 )
In a drainage system, closed drains are to be placed with their centres 2 m below the ground level
to keep the highest position of the water table 1.7m below the ground level. The impervious
stratum is at a depth of 9.6m below the ground level. If the average annual rainfall in the area is
850 mm, find the spacing of drains. Assume 1% of the average annual rainfall to be drained in 24
hours and coefficient of permability k=1×10−5 m/s
Solution
850×A×1 3
1% of average annual rainfall over an area A km2 catered by each drain = 1000×100 m = 85 × 10−4 A m3
This is to be drained in 24 hours and hence the volume of water to be drained per second
−4
= 85×10 A
24×3600 cumec
If L is the spacing beyween the drains, then the area catered by drain per unit length = L m2
85 × 10−4 × L
∴ QD = = 9.84 × 10−8 L cumec/m
24 × 3600
From the equation (4), we have,
4k(b2 − a2 )
L=
QD
b=(9.6-1.7) = 7.9m and a = (9.6-2) = 7.6m
Subtituing these values,
4 × 1 × 10−5 [(7.9)2 − (7.6)2 ]
L=
9.84 × 10−8 L
L2 = 1890.24
∴ L = 43.48m = 43.5 m
In a drainage system, closed drains are to be placed with their centres 10.5 m above the
impervious stratum and the maximum height of the drained water table above the centre of drains
is 0.3m. If the spacing of drains is 36m and the drrains carry 1% of the average annual rainfall in
24 hours, find the average annual rainfall. Take coefficient of permability k=5×10−6 m/s.
Solution
4k(b2 − a2 )
QD =
L
k= 5×10−6 m/s, b= (10.5+0.3) = 10.8m, a = 10.5 m and L= 36m
If the average annual rainfall is x mm, then the volume of water drained per second
x×A x×L
= 1000×100×24×3600 cumec = 1000×100×24×3600 cumec
Substituing the value of QD , we get
x × 36
3.55 × 10−6 =
1000 × 100 × 24 × 3600
∴
∴ x = 852mm