Watercad Eraq
Watercad Eraq
net/publication/305498922
CITATIONS READS
2 1,154
1 author:
Aqeel AL-Adili
University of Technology, Iraq
96 PUBLICATIONS 92 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Drainage investigation of surface runoff for highway pavement ( experimental study). View project
Management And Feasibility Study For Shut Al-Arab - Basarh Region -Analysis, Design and Decision For Water Resources System. Ph.D. project. View project
All content following this page was uploaded by Aqeel AL-Adili on 22 July 2016.
ABSTRACT:
In this research, several software models have been used to design and systemize management
of water resources with adopting meteorological data and measurements at field of the study area.
Soil-Plant-Air-Water (SPAW) computer model was used to analyze the type and properties of soil
of the study area, and then using the results of this model as input data for the CropWAT, 2008
model, to calculate the wheat crop water requirements, and hence irrigation requirements for a field
area of 54 donums with wheat cultivation, and a system of scheduling irrigation works was
designed, after proposing and designing a periodic hand move sprinkler irrigation system, to reduce
the capital cost for the project to enable land owners benefit from it by using (WaterCAD) model.
The results of the designed models were compared in terms of discharges with actually applied
discharge on the ground. This research has indicated that there is a big difference between the
applied discharges and water consumption by the plant. The final design of the system capacity is
76.56 m3/hr to irrigate area of 54 donums, within 9 days of 15 hours per day, while the average of
filed measurements of six fields of wheat in the same area showed that the existing applied
discharge is 327 m3/hr, with spending 10 days to irrigate 85 donums with not less ten hours per day,
so that will contribute to vivification and irrigate about 984400 donums of wheat lands.
1 INTRODUCTION
The term "water resources systems" is comprised of water sources, means of their control and
transportation to water users. However, these sources are limited , and the competition for it is
increasing day by day. The Middle East is the scarcestwater region in the world. Worldwide, the
average water availability per person is close to 7,000m3/person/year, whereas in the Middle East
region, only around 1,200m3/person/year is available. One half of Middle East’s population lives
under conditions of water stress. Moreover, with the population expected to grow, per capita
availability of water is expected to halve by 2050, (FAO, 2012).Agriculture accounts for the vast
majority ofwater resources consumption in Iraq, withdrawing 92%of total freshwater for irrigation
and foodproduction, (Iraqi Ministry of Water Resources, 2012).
Adiwaniyah governorate suffers from network deteriorations of water resources and irregular
distribution of water, despite the multiplicity of sources, which led to the deterioration of
agricultural lands and poor in good soil in the province. Because of poor management of its water
resources and low efficiency of irrigation to less than 27%, thus, low productivity of donum that
significantly affected level of living for the people of the province. For the purpose of addressing
this situation and due to the fact that this province from the provincial agricultural and large core in
the country with a total of arable lands where up to 1.85 million donums, the data showed that the
largest proportion of planted area was 47% in 2008, so it became necessary to find a proper
management for the its water resources (Adiwaniyah Water Resources Directorate, 2012).
New methods to design water resources, irrigation systems, and management for decision
making are being developed which incorporate computer based capabilities with data management,
163
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
analysis, and making decisions on operation and maintenance for better productions. Such tools is a
Water CAD as a Simulation Technique for Decision Support System (DSS), which is defined by
Sprague and Carlson, (1982) as “an interactive computer based system that helps decision makers
utilizes data and models to solve unstructured or under-specified problems.”
The aim of the research is to design and establish rational irrigation planning systems by
employing advanced scientific simulation techniques, and suggest a design management model to
estimate multiple irrigation water networks and scheduling systems at field level for Adiwaniyah
governorate agriculture lands.Thus, the hope is to lead to improving the water use efficiency in
agriculture and spreading of irrigation benefits to tail end areas.
Prior to design any irrigation project, it is necessary to calculate the project water requirements.
This means that exact (correct) amounts of water and correct timing of application should be
adopted. In addition, it will need more wide spread adoption of deficit irrigation, especially in arid
and semi-arid regions. Recent advances in new irrigation technologies will help to identify
irrigation scheduling strategies that minimize water demand with minimal impacts on yields and
yield quality, leading to improved food security. So, details about computed reference
evapotranspiration, crop water requirement, irrigation requirement, and soil investigation will be
discussed in detail.
Evapotranspiration is one of the most important factors to be known prior to making and
implementing any decisions for the design and management of water resources systems.
Evapotranspiration varies day-to-day on the weather patterns and specific crop; therefore it is not
suitable to consider a fixed value. There are many methods proposed by different scientists and
organizations for estimating the evapotranspiration of a reference crop.
Land and Water Development Division of Food Agriculture Organization, (FAO) has developed
a software named CropWAT in 1992 to evaluate the irrigation water requirement, CropWAT (8.0)
(2009) which will be adopted in this research.
The study have been used Water CAD to model a sprinkler network as well as solving a
problem related to an irrigation system. The models consist of an optimization technique from a
network solver Water CAD used to reach the proposed scheme.
Water CAD (8.0) model was chosen, because it handles both steady state and extended period
simulation of water distribution network. This section presents discussions on Water CAD
simulation model using a sample model of a sprinkler irrigation system and a real irrigation system.
Water CAD from Haestad method is a software application for construction of simple and
complex pipe networks. The Water distribution system program is developed by the U.S.A
Environmental Protection Agency’s Water Supply and Water Resources Division. It is compatible
with Geographical Information System (GIS) software, (Maksimovicand Prodanovic, 1996).
164
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
Water CAD uses several equations to calculate the head loss along the pipes and through the
fittings and display it on each element of the model with color coding within the model, the head
loss equation frequently used is the Hazen-Williams formula. It is frequently used in the analysis of
pressure pipe systems (such as water distribution networks and sewer force mains). The formula is
as follows:
The Hazen-Williams C-Factor is given by WaterCAD when the user selects a respective type of
pipe material, the lower the C-Factor, the rougher the inside of the pipe. In this research study, the
default value in WaterCAD was used for the water distribution model. For example, the default C-
Factors for PVC pipes and Cast Iron pipes of 150 and 130 were used, respectively.
The WaterCAD computer model used for water distribution network analysis is composed of
two parts: the input data file and the WaterCAD computer program.The data file defines the
characteristics of the pipes, the nodes (ends of the pipe), and the control components (such as
pumps and valves) in the pipe network. The computer program solves the nonlinear energy
equations and linear mass equations for pressures at nodes and flow rates in pipes. Simulation
carried out is presented to the user in a readable format on the GUI, (Figure-1).
165
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
The total area of Adiwaniyah governorate is 8366.443 km2, which represents 1.9% of the total
area of Iraq. The gross arable area of Adiwaniyah governorate does not exceed 4643.43 km2
(approximately 55% of the total area of Adiwaniyah governorate), (Figure -3).
However, the percentage of actually cultivated land was about 26% of the total area of
Adiwaniyah governorate (Adiwaniyah Agriculture Directorate, 2012), which represents 47% of the
arable land, due to limited and bad management for water resources and inequity of distribution of
water resource of the area, (Figure-3).
166
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
167
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
Basic
Farm Data
Map of Design
Area
Amount Per
Application
Application
Rate
Irrigation
Interval
Irrigation Period
Determination of Pipe
Sprinkler Selection
Size
Total Dynamic
System Capacity
Head
Selection of Pump
Operation Instruction
Where: ETO = Reference crop evapotranspiration (mm/day).RN =Net radiation at crop surface
(MJ-2 /day).G = Soil heat flux (MJ-2/day).T = Average temperature (oC).
The available data of wind speed at the meteorological station of the study area measured at 10
m height, so it must be converted toCropWAT(8.0) models format by using the following
formula(eq.3) (Shariati 1997),
𝟒. 𝟖𝟔𝟖
𝑼𝟐 = ∗ 𝑼𝒁 … … … … … … … … … … … … … … … . … (𝟑)
𝒍𝒏(𝟔𝟕. 𝟕𝟓𝒁 − 𝟓. 𝟒𝟐)
Where:U2 = Wind speed at two meters.Uz= Wind speed in Z meter elevation, and Z= altitude of
measured wind speed.
The result of ET0 of CropWAT (8.0) has been validated with Blaney-Criddle which is modified
by N. Kharofa (1985) Eq. (4).
ETo = C P Tc1.3 ……………………………………… (4)
168
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
Where:C = coefficient calculated for the local site, Kharofa (1985) fixed it by 0.34 for the middle of
Iraq.P= the percentage of the number of daylight hours in the month relative to the number in the
year.Tc = average monthly temperature (Co).
This variation in ETo values is attributed to combined effects of temperature, sun shine hours,
radiation, wind speed and humidity. The increase in ET o during May to Sep. in the study area can
be related to the change in temperature and wind speed.
In specific case, the computation of the net depth ofwater application requires additional inputs
corporate CropWAT (8.0) model as the following:
a. The available soil moisture [FC - Permanent Wilting Point (PWP)].
b. The allowable soil moisture depletion (P) as a percentage from the available soil moisture.
c. The effective depth of root zone of the crop (RZD).
Soil survey and tests should be done to determine the field capacity (FC) and permanent wilting
point (PWP) of the soil. Reports and data from Huriyah-Daghara Project's field laboratory were
used to estimate the soil texture, FC, PWP and P. by using SPAW computer model simulate. The
difference between the field capacity and permanent wilting point will give the available soil
moisture (water holding capacity), which is the total amount of water that the crop can use.
Depending on the crop sensitivity to stress, the soil moisture should be allowed to be depleted only
partially. For most field crops, a depletion of 50% of the available moisture is acceptable.
The effective root zone depth of the crop under consideration can be established from specified
tables. Some tests and measurements for root depth of wheat crop at the field of the study were
taken and compared with measurements of Agriculture Directorate of Adiwaniyah Governorate,
these tests and measurements showed that the RZD of wheat is (0.8-1m), so that it was taken as 1m
as the worst case for this research.
169
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
The maximum net depth to be applied per irrigation will be at the maximum root zone depth at
the mid-season of the crop depending on the crop cycle, can be calculated, using the following
equation, (FAO, 1992):
Where:dnet = Readily available moisture or net depth of water application per irrigation for the
selected crop (mm), FC = Soil moisture at field capacity (mm/m), PWP = Soil moisture at the
permanent wilting point (mm/m), RZD = the depth of soil that the roots exploit effectively (m), P=
the allowable portion of available moisture permitted for depletion by the crop before the next
irrigation.
Depending on Huriyah-Daghara’s field laboratory data and the SPAW model results with
applying equation-4in CropWAT model, the maximum dnet for the area of the case study is 60.4
mm.
In order to express the depth of water in terms of the volume, eq.5 calculated this volume, the
area proposed for irrigation must be multiplied by the depth:
Where:A = Area proposed for irrigation (ha), d = Depth of water application (mm), 10 = Conversion
factor.
For the field of the study area of an area of 54 donums (13.5 hectare), depending on the results of
CropWAT (8.0) model and using equation -5,a net application of 8154 m3of water will be required
per irrigation to bring the root zone depth of soil from the 50% allowable depletion level to the field
capacity at worst case within the period between 16 November and 30 of April.
Crop evapotranspiration per decade was calculated by multiplication of the number of effective
crop days. To convert monthly rainfall data to decade values, a linear interpolation was carried out.
Values for first and third decades of each month were respectively calculated by interpolation with
the preceding and successive months.
Crop water requirements were then calculated as the difference between the crop
evapotranspiration and effective rainfall, (Table-3).
170
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
Where:IF= irrigation frequency (days), dnet= net depth of water application (mm), ETC= peak daily
water use (mm/day)
According to CropWAT (8.0) model result, the peak water demand of wheat at Adiwaniyah
meteorological station was estimated to be 6.38 mm/day. Therefore, using equation 4 will give
irrigation frequency IF equal to 60.4/6.38=9.47 days.
The system should be designed to provide 60.4 mm every 9.47days (113.64 hours). For practical
purposes, fractions of days are not used for irrigation frequency purposes. Hence, the irrigation
frequency in this case should be 9 days, with a corresponding dnet of 57.42 mm (6.38 x 9), and by
using equation-4,the moisture depletion should be 34%.
The gross depth can be calculated from equation (7), (Irrigation Manual, 2001):
𝒅𝒏𝒆𝒕
𝒅𝒈𝒓𝒐𝒔𝒔 = ……………………………………. 𝟕
𝑬
The farm irrigation efficiency of sprinkler systems varies from climate to other as in table (4),
(FAO, 1992).
𝟏𝟎 ∗ 𝑨 ∗ 𝒅𝒈𝒓𝒐𝒔𝒔
𝑸= … … … … … … … … … … … … … … … … (𝟖)
𝑰 ∗ 𝑵𝒔 ∗ 𝑻
Where:Q = System capacity (m3/hr.), A = Design area (ha), d = Gross depth of water
application (mm), I = Irrigation cycle (days), Ns = Number of shifts (movement) per day,T =
Irrigation time per shift (hr.).
For the area of 13.5 ha, in order to achieve the maximum degree of equipment utilization, it is
desirable, but not always necessary, that the irrigation system should operate for 15 hours per shift
at one shift per day during peak demand and take an irrigation cycle of 9 days to complete irrigating
the 13.5 ha. Substituting the values in equation-8gives a system capacity of 76.56 m3/hr.
8 DESIGN CRITERIA
8.1 Design Steps for Periodic-Move System
Once the preliminary design parameters are obtained, the design adjustment can commence. The
adjustment allows for the revision of the preliminary design parameters, in order to suit the
physical, human, financial, and equipment performance limitations or impositions. The next design
step is to select the Sprinkler system design and the spacing.
172
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
manufacturer usually states the value of the discharge coefficient as the flow rate in CMH through
the device at a 1 psi pressure drop.
Emitters are used to model flow through sprinkler systems and irrigation networks. Users can
also be used to simulate the leakage in a pipe connected to the junction (if a discharge coefficient
and pressure exponent for the leaking crack or joint can be estimated) and compute a flow at the
junction. In the latter case, a very high value of the discharge coefficient e.g., 100 times the
maximum flow expected can be used.
When both an emitter and a normal demand are specified for a junction, the demand that Bentley
Water CAD reports in its output results includes both the normal demand and the flow through the
emitter.
The flow through an emitter is calculated as; (WaterCAD user’s manual, 2008):
𝑸 = 𝒌𝒑𝒏 … … … … … … … … … … … … … . . … (𝟗)
Where:Q = Flow (m3/hr), k = emitter coefficient (property of the node), P = Pressure (n/m2), n = is
the emitter exponent and is set globally in the calculation
Options for the run; it is dimensionless but affects the units of k. The default value for n is 0.5
which is a typical value for an orifice, (WaterCAD user’s manual, 2008).
Figure-5, shows that, the designed sprinkler network consist of the main line of 445 m in length,
three pairs of laterals are laid perpendicularly on the main line; each lateral pipe is 135m in length,
the distance between two laterals is 135 m as well, and 66 Sprinklers carried by the laterals with 11
sprinklers per lateral.
173
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
The selection of the correct sprinkler depends on how the best fit spacing with a certain pressure
and nozzle size can provide the water at an application rate that does neither cause runoff nor
damage the crop and at the best possible uniformity under the prevailing wind conditions.
In this study, where a precipitation rate of 6.44 mm/hr of chosen sprinkler is compatible with the
soil and wheat crop, there are several nozzle size, pressure and sprinkler spacing combinations to
choose from (Table-5), so there is another aspect to consider in selecting a sprinkler which is the
energy cost.
Lower pressures are preferable as long as the uniformity of application is not compromised. The
coefficient of uniformity (CU) is a measure of the uniformity of water application. A value of 100%
indicates perfect uniformity, which means that the water is applied to the same depth at each point
in the field. As a rule, the selected sprinkler should have a CU not less than 85% or more. It is
advisable to avoid using the lowest pressure since usually this is the pressure that corresponds to
low CU values,(Naser, 2003).
Referring to the sprinklers manufacture tables, there area number of spacing fits on the land of
the case study, such as 12*12m, 12*15m, 12*18m, and 15*15m. The next step is to find out how
the winds will affect the spacing between two sprinklers and laterals. For this purpose, the mean
wind velocity of the windiest month of the year is considered. Most designers set the maximum
spacing of sprinklers based on the information of table-6. It should be noted also that in the
rectangular pattern, better distribution is obtained when the lateral is placed across the prevailing
wind direction.
174
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
In the case of the research, where the average wind velocity in June is 4.1 m/sec (14.7 km/hr)
and in March is 3.4m/sec (12.2 km/hr), the sprinkler spacing should be based on 45% of sprinkler
wetted diameter (D) for square pattern, and 40% of D∗ 60% of D for rectangular pattern.
From table-5, the sprinkler of the 4 mm nozzle at 350 kPa with 30.50 m of wetted diameter and
6.44 mm/hr precipitation rate was chosen for this design. From table-5,the spacing for a rectangular
pattern for 14.7 km/hr wind speed (10-15 km/hr), 40% of D and 60% of D for the 12 m x 15 m
spacing are 12.2 m (> than 12 m sprinkler spacing) and 18.3 m (> than 15 m lateral spacing),
respectively. Therefore, the wind requirements are satisfied the 12 m x 15 m spacing.
Therefore, 6.44mm/hr sprinkler precipitation rate is less than infiltration rate of the soil of the
study area (7.26mm/hr), so it is compatible with the infiltration rate, and the 4 mm nozzle can be
considered.
It is advisable to determine whether the sprinkler with a 4.0 mm nozzle would satisfy the wind
requirements at the 12 m x 18 m spacing at 300 kPa. At this pressure, the wetted diameter is 26.60
m. 40% of D and 60% of D are 10.64 m (< than 12 m sprinkler spacing) and 15.96 m (< than 18 m
lateral spacing), respectively. For the 15 m x 15 m spacing, 45% of D is 11.9 m (0.45 * 26.60),
which is less than the sprinkler and lateral spacing of 15 m each. Therefore, the 4.0 mm nozzle
operating at 300 kPa pressure does not meet the wind requirements either under 12 m x 18 m
spacing or 15 m x 15 m spacing as the wetted diameter is too small compared to the desired spacing
requirement.
The next step is to determine the set time (Ts), which is the time each set of sprinklers should
operate at the same position in order to deliver the gross irrigation depth, and establish whether it is
acceptable. So, it can be found out according to equation 10,
𝒅𝒈𝒓𝒐𝒔𝒔
𝑻𝒔 = … … … … … … … … … … … … . … … … (𝟏𝟎)(FAO,1992)
𝑷𝒓
Where:Ts= Set time (hr), Pr= Sprinkler precipitation rate (mm/hr), Substituting the values in
Equation (4.6) gives:-
Hence, each set of sprinklers should operate at the same position for 12 hours in order to deliver
the 76.56 mm gross application per irrigation for nine days. If it is assumed that the design is a
permanent system, as in (Figure-6), this would have been ideal, because it will have full utilization
of the equipment by having two sets per day. However, if the design is a semi-portable system, as in
(Figure-7), where the laterals have to be moved from one position to the next, there would be time
available to move the laterals between each of the two shifts for the next day during the peak water
demand period. In this case, the following choices will be adopted:
1. To increase twice the number of operating laterals so that extra laterals are moved while the
other laterals are operating, or
2. To re-assess the moisture depletion level, or
175
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
3. To use a different sprinkler with the same or different spacing, nozzle, pressure and
precipitation rate.
176
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
the irrigation cycle would be 9 days. In order to apply the 57.42 mm net per irrigation, the gross
application at 75% efficiency should be 76.56 mm (57.42/0.75).
Laterals in a semi-portable system are PVC pipes with multi-outlets (sprinklers) along their
length. The friction losses were simulated and calculated by WaterCAD model with Hazen-
Williams equation, the model corrects the pipe diameter, since the flow reduces along the lateral.
While the class 6 pipe is used for surface irrigation, the most commonly used classes for
pressurized irrigation systems are the class 12.
𝐯𝟐
……………………….. (11)
𝟐𝒈
Keller &Bliesner, (1990) recommended that for centrifugal pumps the diameter of the suction
pipe should be selected such that the water velocity v<3.3 m/s in order to ensure good pump
performance. Assuming this maximum velocity for the flow and applying the above formula, then
the velocity head corresponding to the minimum diameter of the suction pipe that can be selected to
satisfy this condition is 0.56 m (3.32/ (2 x 9.81)).
Considering that the water level is at 18.00 m, so there will be no difference in elevation between
the water level and the eye of the impeller (located at 18 m elevation). Since the maximum velocity
head is 0.56 m, the suction lift is 0.56 m (0 + 0.56). Assuming minor losses in fittings and a short
suction pipe, the suction lift is rounded up to 2 m.
The elevation difference is obtainable from the contour map and is approximately 0.5 m.The
length of the supply line is 25 m. The friction losses for the supply line are computed as follows:
177
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
The friction losses in the riser could be assumed about 0.25 m per m of riser, and for fittings, it
is usually taken 10% of the total head losses. Additionally, the height of the riser should be included
in the calculations. This was assumed to be 2 meters in this design, to cater for tall crops, such as
maize/corn. So that the total head requirement in the system must be not less the minimum
operation sprinkler pressure.
The main line of the system consists of three parts according to diameter, carried discharge and
unit of the head losses along each part of the pipe, part one (main line -1) includes pipes 1,2 and 3
with diameter of 150 mm and operation discharge of 76.56, 76.56, 63.8 m3/hr,respectively. Main
Line -2 includes pipes 4 and 5, which are 125 mm in diameter with operation discharges of 51.04
and 38.2m3/hr, respectively. Pipe 6 and pipe 7 form main line 3 which has diameter of 100 mm and
discharge of 25.52 and 12.76m3/hr. respectively, (Figure-8),(table-7).
178
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
From WaterCAD model calculation, the maximum pressure in the system is 395 kPaoccurring at
the first sprinkler on the field(Figure-10), whilst the minimum pressure is 351 kPa occurring at the
last sprinkler (Figure-11). Thus, the pressure variation is within the constraints limits.
179
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
CONCLUSION
The main concern in the design water resources system of pressurized irrigation networks
is to achieve an optimal solution which satisfies the constraints of the network, especially if
the network is large.The following conclusions were derived from the present study as follows;
1- The designed models employed and applied to irrigate about 368200, 304370 donums of
wheat and barley respectively at winter seasons, and to irrigate about 17485, 46543 donums
of maize/corn and sorghum respectively, at summer season, to save more water quantities
to vivification and irrigate about 984400 donums of heathy land.
2- The final design of the system capacity is 76.56 m3/hr to irrigate area of 54 donums, within
9 days of 15 hours per day. While the average of existing filed measurements of six fields
of wheat in the same area showed that the existing applied discharge is 327 m3/hr, and
spending 10 days to irrigate 85 donums with not less ten hours per day. This means each
donum gates 385 m3 per irrigation.
3- This study showed that each donum requires only 150 m3per irrigation, so that the
application of sprinkler irrigation will assist to doublicate the area of cultivation by about
2.5 times, in case of wheat crop.
4- This study revealed the ability of applying these models on other regions inAdiwaniyah and
other Iraq governorates.Itis recommended due to the efficiency and compatibility of such
simulation especially if the network is large.
5- The model was able to simulate water management, the calibration, and validation results
with Kharoofa formula indicated that the results were compatible.
6- However, this study recommend the possibility of employing this simulation model on
other crops such as barley, maize, andsorghum can be achieved to assess the irrigation
system.
REFERENCES
Addink, H. and Bytat, 1989. “Design and operation of farm irrigation systems”, ASAE,
USA, 31(3):pp. 821-829
180
International Water Technology Journal, IWTJ Vol. 5 –No.3, September 2015
181