Lecture 2 Irrigation Engineeing II
Lecture 2 Irrigation Engineeing II
The field irrigation requirement represents the rate (in mm/d) at which water must be delivered to
the field to prevent the crop suffering a shortage of water. Use data from the feasibility study.
fb
1
D
d
14
To minimize losses, the canal banks should be built from clayey soil and constructed in
layers, with each layer compacted using heavy rammers.
The required size of the canal can be decided using Manning’s formula:
0.5
A x R 0.67 x S o
Q
n
Where: Q = discharge (m3/s). Note: 1 m3/s = 1000 l/s
A = wetted area (m2)
R = hydraulic radius (m) (= wetted area/wetted perimeter)
So = slope (fraction) for example 0.001 (or 1 m/km)
n = Manning’s roughness coefficient (commonly taken as 0.03 for small irrigation
canals)
Refer to a design chart for other assumptions.
Example 2
During a baseline study for a 10ha farm in Loitoktok, the following information was gathered from a
farmer. Average crop factor (Kc) = 1.2 and reference crop evapo-transpiration (ETO) = 7.6 mm/day.
Table 1 presents details of crops grown
16
Table 1: Commonly grown crops
Crop Rooting depth Depletion factor (MAD)
D(m) n
Sorghum 0.6-1.5 0.6
Maize 0.6-1.0 0.5
Table 2: Suitable side slopes for channels built in various types of materials (chow, 1959)
The most prevalent soil type is sandy loam, with moisture storage capacity of 150mm/m.
A gravity fed surface irrigation scheme is to be established. The irrigation system will be
composed of the following design parameters: conveyance, field canal and application
efficiencies as 90%, 85% and 55% respectively.
(a) Crop water demand, ETC, is obtained using Equation 1
ETC = K c ETO , Reference evapo-transpiration ETo = 7.6 mm/day
ETC = 1.2 7.6 mm/day = 9.12 mm/day
(b) Taking conveyance along the earth canal, distribution and application efficiencies as
90%, 85% and 50% respectively. The overall efficiency is determined by:
Ea = Ec Ed E f ,
Where, Ea is overall efficiency, Ec is conveyance efficiency, Ed is field canal efficiency and Ef is
application efficiency all in decimals.
Ea = 0.90 0.85 0.55 0.421 or 42.1%
(c)Take overall efficiency of 42 %, the unit gross water demand (UGWD), which is total
amount needed including losses is obtained by:
NWR ET C
UGWD =
Ea Efficiency E a
9.12
UGWD = = 21.71 mm/day
0.42
But, 1 mm/day is equivalent to 0.116 l/s/ha which is also equivalent to 10m3/ha/day,
Since 1 mm/day ≡ 0.116 l/s/ha, Then,
17
UGWD A
(d) SWR = (l/s)
h d
*
24 7
Where, ”SWR” is scheme water requirement, “h” is number of irrigation hours per day, “d” is
number of irrigation days per week; “A” is area to be irrigated in ha.
Therefore,
2.94 100
SWR =
15 6
*
24 7
2.52 60
SWR = = 235.2 l/s Therefore, take Q design = 0.235 m3/s
18 * 6
24 7
Longitudinal slope of the main canal was determined using the following Equation
E E H *100
S O = WL
L
Where “So” is longitudinal slope as a percentage (%), “EWL” is river water level in “m”, “EH” is
elevation at the head of scheme and “L” is the length of canal in “m”. For example,
SO =
97.589 97.21*100 0.026%
1456
The main and sub-mains will be earth canals and thus not lined. Therefore, the canals will be
trapezoidal in shaped and the section is presented in the following figure
18
B
fb
1
D
d
Z=2
Canal design
The design is based on characteristics of sediment load i.e silt in canal water.
The irrigation canals are of a particular shape say trapezoidal in shape with side slopes of 2:1,
for sandy loam soils, select minimum allowed canal slope of say 0.05% based on guidelines.
Q=AV; V = (1/ n) R 2/3 so½
Permissible Velocities: Minimum and Maximum
canals carrying water with higher velocities may scour the bed and the sides of the channel
leading to the collapse of the canal. On the other hand the weeds and plants grow in the
channel when the nutrients are available in the water and the flow velocity is low.
Therefore, the minimum permissible velocity should not allow the growth of vegetation
such as weed, hycinth as well you should not be permitting the settlement of suspended
material (non-silting velocity).
"Minimum permissible velocity" refers to the smallest velocity which will prevent both
sedimentation and vegetative growth in general.
Maximum permissible velocities entirely depend on the material that is used and the bed
slope of the channel.
For example: in case of chutes, spillways the velocity may reach as high as 25 m/s. As the
dam heights are increasing the expected velocities of the flows are also increasing and it
can reach as high as 70 m/s in exceptional cases. Thus, when one refers to maximum
permissible velocity, it is for the normal canals built for irrigation purposes and Power
canals in which the energy loss must be minimised. Hence, following table gives the
maximum permissible velocity for some selected materials.
19
Table 4: Maximum permissible velocities for selected materials
Freeboard
The term freeboard refers to the vertical distance between either the top of the channel or
the top of the channel is carrying the design flow at normal depth. The purpose of
freeboard is to prevent the overtopping of either the lining or the top of the channel
fluctuations in the water surface caused by:
(1) Wind - driven waves,
(2) Tidal action,
(3) Hydraulic jumps,
(4) Super elevation of the water surface as the flow goes round curves at high velocities,
(5) The interception of storm runoff by the channel,
(6) The occurrence of greater than design depths of flow caused by canal sedimentation or
an increased coefficient of friction, or
(7) Temporary mis-operation of the canal system.
There is no universally accepted rule for the determination of free board since, waves,
unsteady flow condition, curves etc., and influence the free board.
Q= (A/ n) R2/3 so½, Where, Q is water flow rates in m3/s, n is roughness factor, A is cross
sectional area of canal in m2, R is the hydraulic radius of canal in m and So is the canal slope
in m/m.
20
A= bd + zd ; R=A/P ; P= b + 2 d √(1 + z 2) , Z for sandy loam should be above 2, to
prevent canal banks from collapsing
Take the normal water level in the river as 0.5m, slope is say 0.06%. Using trial and error
method, the main canal characteristics can be computed and presented in Tabular form as
shown below.
In summary canal characteristics are: Z=2, d=…m, fb =….m, D=…m, b=……, B=……m and Q=
………m3/s, V=…..m/s
Q(m3/s) z(m-1) S (%) B(m) D(m) Fb(m) H(m) V(m/s) Fr T(m)
Take overall efficiency of 36 %, the unit gross water demand (UGWD), which is total amount
needed including losses is obtained using Equation 3
NWR ET C
UGWD = [3]
Ea Efficiency E a
9.12
UGWD = = 25.33 mm/day
0.36
But, 1 mm/day is equivalent to 0.116 l/s/ha which is also equivalent to 10m 3/ha/day,
therefore
UGWD =0.116 * 25.33 = 2.94 l/s/ha
Scheme water requirement (SWR) is water needed to irrigate all the farms in the scheme.
Given that the proposed Irrigation scheme has 100ha to be irrigate 6 days a week and for 15
hours a day, SWR needed is obtained from Equation 4.
UGWD A
SWR = (L/s) [4]
h d
*
24 7
Where, h is number of irrigation hours per day, d is number of irrigation days per week; A is
area to be irrigated in ha.
Therefore,
2.94 100
SWR = = 548.8 l/s
15 6
*
24 7
Therefore, take Q design = 0.549 m3/s
22
Conveyance system design
Longitudinal slope of the main canal was determined using Equation 5
E E H *100
S O = WL [5]
L
Where So is longitudinal slope as a percentage (%), EWL is river water level in m, EH is
elevation at the head of scheme and L is the length of canal in m.
SO =
97.589 97.21*100 0.026%
1456
The main and sub-mains will be earth canals and thus not lined. Therefore, the canals will be
trapezoidal in shaped and the section is presented by Figure 5.1
fb
1
D
d
The flow velocity in the main canal was determined through Equation 6.
2 1
R3S 2
V [6]
n
Where V = Flow Velocity in the main canal (m/s)
R = Hydraulic radius of channel (m)
= Channel-cross sectional area/wetted perimeter
S = Canal (channel slope /channel bed longitudinal slope (m/m)
n = Manning’s roughness coefficient/Manning’s resistance coefficient for open channel
flow.
Hydraulic radius (R) was obtained through Equation 7
R= A [7]
P
Where, “A” is channel cross section area in m2 and “P” is the wetted perimeter in m.
“A” was obtained through Equation 8
A = bd 2d 2 = d (b+2d) [8]
23
Where, “b” and “d” are internal width and water depth (level) in the canal, respectively.
“P” was determined through Equation 9
P = b 2d 5 [9]
The canal sectional factor determines the size of the canal through Equation 10b, given by
Manning’s Equation(10a) derived by combining Equations 6 and 8, and then rearranging
which gives Equation 10b.
1 2
S 2 AR3
Q [10a]
n
2
Qn
Canal sec tion factor A R 3 [10b]
S 1/ 2
Taking S= 0.05% and n= 0.023 (FAO, 2002) for an earth channels, then canal section factor
was computed as:
2
Q n 0.549 * 0.023
Canal sec tion factor A R 0.5647
3
S 1 / 2 0.05 *10 2 1 / 2
Then specifications are obtained through trial and error using excel spread sheet as
indicated in Table 5.3.
24
From Table 5.3, the main canal bed width was determined as 2.0m, water depth of 0.43m
and Free board is determined by Equation 11
fb = 0.3D or 0.2m which ever value of the two that is higher [11]
Where fb is free board in m, 0.3 is discharge coefficient and “d” is depth of water in canal in
“m”.
Therefore, fb = 0.3*0.43 m = 0.129m.
Since 0.2m is higher, take fb = 0.2m,
The total depth of the canal was determined by Equation 12
D = d + f b , Where D is total depth in m [12]
D 0.43 0.2 0.63 m
Velocity of flow in the canal was determined by previous equation (Equation 6) and is given
as:
1
2/3 2
R S
V= [13]
n
1
2
0.05 2
0.31 3
V= 100
= 0.445m/s
0.023
Therefore the flow is within safe velocity since it is less than maximum allowable velocity of
0.6-0.7 m/s (For earth canal, allowable velocity should not exceed 0.7m/s and the flow
should be laminar (Fr<1). Higher velocity more than 0.7 m/s for sandy loam soils
encourages soil erosion in the canal (FAO, 2002).
Froude number, Fr determines whether the flow is lamina or tranquil. Value beyond 1 gives
supercritical energy status, hence turbulent water flow. Whereas the value less than 1 gives
sub- critic energy level and flow is lamina. The following Equation 14 adopted from Michael,
(1978) was used:
V
Fr [14]
gy
Where (Fr), is the Froude number, V is velocity of flow in m/s, g is acceleration due to gravity
in m/s 2 and y is the depth of flow of water in m. Therefore,
0.445
Fr 0.217
9.81 * 0.43
The flow velocity is okay since Froude number is less than 1, hence sub- critical (lamina flow).
Erosion during peak flow
When canal fills during flood flows, peak flow occurs. To estimate erosion, peak values are
computed using Equation 15 as:
1 2/3
VP RP S 1 / 2 [15]
n
25
Where: V p is peak flow velocity and R p is hydraulic radius (peak flow)
D = d + fb = (0.43 + 0.2) m = 0.63m
Given that R2/3 = 0.565, S = 0.05%, R = A/P= 0.31m
Pp = b + 2D (1 +z2)0.5 = 2.0+ 2*0.63(1+22)0.5= 4.82 m
Ap = D (b+zD) = 0.63(2+2*0.63) m2 = 2.05m2
B = b + 2zD = 2+ 2*2*0.63= 4.52m
Rp =Ap/Pp= 2.05m2/4.82 m = 0.425 m
Rp 2/3= 0.565
1 2/3
VP RP S 1 / 2
n
Vp = (1/0.023) * 0.4252/3 * (0.05/100)1/2 = 0.550 m/s
Velocity of peak flow (Vp) will not be erosive since it is still within the safe margins (less than
0.7 m/s)
The peak discharge (Qp) was estimated from Equation 16 as:
Qp = Ap*Vp [16]
3
Qp = 2.05*0.550 m /s = 1.127 m /s 3
Specifications for canal are bottom width (b) is 2.0 m, flow depth (d) is 0.43m, D is 0.63 m, Q
is 0.549m3/s, V is 0.445 m/s, QP is 1.127m3/s, VP is 0.550 m/s and Fr is 0.217.
The velocity of water along the canal was computed as 0.445 m/s, hence the width of intake
orifice was computed as 0.8m. This width is too wide for a single gate, hence 2 gates of
dimensions 0.4m are preferred and with a height of 0.63m.
Intake Design
Through the information gather at the intake site (Figure below) using participatory
approaches (Focus group discussions and Key informant interviews) on traditional
knowledge in water resources, the following data was obtained:
River bed level = 97.072m
Normal water level = 97.822m
Flood water level = 99.387m
These values were corresponding to river stage of:
Normal and flood flows of 0.75m and 2.315m respectively.
The river is very wide about 100m and with very loose bed materials that are sandy in nature.
Hence, the most appropriate head works is side weir, which depends on water level in the
26
river. Hence height of wing wall to be 2.4m high. During normal water flow, the orifice
functions as rectangular open channel. The sketch of the intake section is presented in the
Figure below.
99.387m
97.572m
H
0.5m
97.072m
The average altitude is 549m asl from the topographical map. It should be noted that survey
data elevations are different from those of the national grid and have been based on
elevations of temporary bench marks established during the survey.
During such flood flows, the minimum orifice opening is computed using the orifice
Equation as:
QOrifice = CA 2 gh
Where Qorifice = intake design discharge / flow rates through orifice (m3/s);
C = coefficient of discharge =0.65
h = head loss over the orifice/head of water measured from the center of orifice
(m);
27
However, there is need to contain all the water in the event of flooding and wave action.
Hence, the height of the steel gates shall be (99.387-97.287) m = 2.1 m. Ideally a beam 0.5
m breadth will be constructed to reduce the gates to a height of 2.1m-0.5m = 1.6 m.
Dimensions of the gate are: Depth = 1.6 m and width = 0.40m
Velocity of flow across gate (Vg) was calculated using Equation 18 as:
Q Q 0.549
Vg 1.373 m / s
A B * D 0.8 * 0.5
To reduce silt entering the intake structure, a chamber will be constructed in front of the
intake. Wooden planks will be used to control the level of water that can freely flow into the
intake. The height of the wing walls will be of minimum of 2.1m, therefore constructing to
the height of river bank of 2.4m. The intake consists of 2 No. rectangular orifices with sluice
gates as shown in the drawing. The orifice opens an outlet into a trapezoidal earth canal
section as shown in the drawing. Orifice flow occurs when the upstream flow is sufficient to
seal the opening/gate. In other words, the bottom of the gate is lower than the upstream
water elevation.
Sedimentation tank or Settling Basin
28
Settling Model
Where: t = V/Q
A = surface area of the basin
Furrow systems
Shape and size of furrows
• Must be sufficiently high to take discharge without overtopping
• Depth typically 15 to 20 cm
• Top width typically 25 to 30 cm for most situations but wide, flat ridges needed for
some (clay) soils or crops
Furrow slopes
Discharge is function of S1/2
Upper limit is about 2% slope but if there is high intensity rainfall and erodible soils,
upper limit may be as little as 0.3%
Allowable slopes for furrows:-
Duration of irrigation
Evaluation of required duration of irrigation, T is required to be calculated unless
tensiometers or neutron probes are used to indicate when sufficient water has been
applied. Phillips equation is not easy to solve for t. For this situation, the use Kostiakov’s
equation is suggested:
I = Kta
Where:
I = Cumulative Infiltrated depth (I)
t = Elapsed time t since the start
log I = log K + a log t
Plot log I against log t to obtain K (from the intercept) and a (from the slope).
29
T ( IKr )
1
Ta
Tf
4
Tf = Time to reach the opposite end of the field (hr)
Ta = Application time (hr)
Distribution problem and control of inflow
Advance & recession wave
Need water to infiltrate to bottom of root zone at end of field (this changes as the
crop is established)
Apply maximum possible inflow to get wetting front to the bottom of the field then
cut back to reduce wastage - ideally water should continue to reach bottom of field
but only a small amount should drain off
30