Free Residual Chlorine Calibration by WaterCAD at El-Nozha Water Network in Alexandria Governorate , Egypt

Most of developing countries suffer from decreasing and poor quality of drinking water which led to emergence of many dangerous diseases. In addition, there isn’t any methodology followed to predict and disinfect drinking water by using advanced software program. Maintaining water quality in the water distribution system has become a prominent issue in the study of the water network. Residual chlorine concentration is the indicator to ensure the quality of water in the water network because it eliminates contaminants in the distribution network beginning of the treatment plant down to the consumer. In collaboration with Alexandria (Egypt) Water Company, samples were taken from El-NOZHA water plant station and El-HADARA water distribution network to know the free residual chlorine. In this paper, WaterCAD software has been used to make hydraulic analysis and calibration of residual chlorine in water distribution network to know the ideal chlorine dose that should be added at the water treatment plant and to know the areas of strength and weakness in the concentration of free residual chlorine in the water distribution network. In addition different scenarios have been found to know the free residual chlorine at the weakness areas after injecting chlorine in some junctions and the impact of a fire case or breaking in the water pipe distribution network on the residual chlorine. Results showed ensuring in the water quality in the distribution network by adding chlorine dose in water less than the existing dose which has been added in the El-NOZHA water treatment plant. It is possible to maintain the percentage of free residual chlorine concentration at different locations without relying on adding chlorine only in water treatment plant by injecting low percentage of chlorine dose in the junctions.


Introduction
Water is essential for all forms of life on the earth.A human can survive for weeks without food but only a few days without water.The water is essential for growth and preserves our bodies.The treated water which does not contain any contaminates revives and heals our body from many diseases.On the contrary if the water is polluted that's a case setback in human health and cause many diseases that are known today [1].Disinfection is an important step in ensuring that water is safe to drink.Water treatment plants add disinfectants to destroy microorganisms that can cause disease in humans.There are several methods of disinfection\n such as chlorination, chloramines, ozone, and ultraviolet light.Other disinfection methods include chlorine dioxide, potassium permanganate, and nanofiltration [2].The most important methods of disinfection are chlorine dioxide (the most common disinfection method).The chlorine was discovered in 1774 by Carl Wilhelm Scheele [3].The first use was as a weapon against humans in World War I by Germany.After that, scientists knew its advantages and its effectiveness in water disinfection [4].Water quality management has become very important, especially with the spread of diseases transmitted by water.Computer software lets sanitary engineering make more controls at water disinfection and ensure water quality for the customer.The Alexandria Water Company setup plans to develop control methods at water quality management which includes:  Application of computer software to develop the water plant station control. Using GIS and SKADA system to know strengths and weaknesses at distribution system. Encouragement scientific research and analysis of the results by using advanced software. Obtaining international certification at water quality control.
This research focuses on the water quality control by measuring the residual chlorine concentration in water distribution system and calibrates the doses system to reach the optimal chlorine dose.

Study Area
Alexandria city consists of six districts (Table 1) in addition to Burj Al Arab city.These districts are divided into sixteen administrative sections.The Alexandria total area is 723,271 acres of which 73,392 acres is a construction area and represents 10.15% from the total area of the province.The other areas are deserts lands, agricultural land and surface water (Lake Mareotis) [5].
In 2007 the total population was 194.206 which represent about 32.3% of the total population of the middle district.The section consists of six subsections; it's called Sheyakha.The total population for each subsection is shown in (Table 2) [6].
As shown in (Figure 1) the study area is located in Hadara Kebili and Ezbet El-Gamei with a total population about 87,000 people.These two areas are in the service area of El-Nozha water plant station; after investigation, El-Nozha water plant station was found to be the suitable location for this research.El-Nozha zone is chosen because it is almost mediates the Alexandria city.Furthermore, it is located between three water stations:  In the northeast, Seyouf water plant station. In the west, Furn El Jaraya water plant station. In the southwest, Mansheya water plant station.

Description of Selected Water Distribution System
Alexandria water networks are mixed between two systems (Grid-Iron & Dead-End system); this is due to the urban extension there are many zones with dead ends or tree system.This system causes problems in hydraulic analysis and residual chlorine investigation because there are varieties of pipe diameter, material and dead pressure zones.At the GIS Central Laboratory, all existed water network drawn by using GIS ARCMAP software   version 10.This facility provides the engineers and the technicians to do surveys for the water network to know detailed information about network type, diameter, length and age.Also, they can record any pipe replacement and renewal, to enable them to develop future plan, and identify the strength and weakness parts in Alexandria water company (AWCO) service area [5] [7].

Locations & Dates of Sampling Points
This section covers the selection criteria of sampling points, which called lab samples in this search.The consideration about sampling is detailed in these following lines:  The lab samples should be from multiple areas (near, centre & far) from the water plant station. The lab samples are from closed zone areas to avoid the effect of residual chlorine dose in water network from another water station. Selected areas should have varieties of pipe material, to know the effective of the each pipe martial (wall reaction) on the residual chlorine. Samples also should be from the dead end of pipe to ensure the residual chlorine percentage. Samples also should be taken from the water treatment plant to know the residual chlorine dose before its interaction in the water distribution network.
In order to maintain the above criteria, thirty three locations were selected from various parts of the studied area.Samples were collected from the actual water consumption tap that may be at restaurant, shops or other activity place, in addition other locations were selected in a specific area which is nearly from El-Nozha water plant station as shown in (Figure 2).
Four months throughout the year were chosen (December, January, April & July), and samples were taken along the month to know the actual residual chlorine concentration in water network.These months were chosen to cover all probable seasonal variations.Total numbers of the samples were about two hundred and sixty samples per month.Samples were collected two times daily; the first time at eleven o'clock a.m. and the second time at four o'clock p.m.These two times were chosen to know the chlorine concentration at two different times during the day.Samples were also taken from the water plant, which is the main source of chlorine.

Analysis Tool
The WaterCAD versions started from (V4) to latest version (V8i).Each version has new options and more possibilities to be useful to any sanitary engineer to do complex hydraulic analysis and ensure water quality in any distribution system.In this work Bentely WaterCAD (V8i) has been used.This software is used to perform steady-state and extended hydraulic analysis of water distribution network.
The steady state analysis determines the operating behavior of the system at a specific point in time or under steady-state conditions (flow rates and hydraulic grades remain constant over time).This type of analysis is useful for determining pressures and flow rates under minimum, average, peak, or short term effects on the system due to fire flows.For this type of analysis, the network equations are determined and solved with tanks being treated as fixed grade boundaries.The results that are obtained from this type of analysis are instantaneous values and may or may not be representative of the values of the system a few hours, or even few minute.
The extended state type of analysis allows engineer to model tanks filling and draining, regulating valves opening and closing, and pressures and flow rates changing throughout the system in response to varying demand conditions and automatic control.
WaterCAD also perform water quality simulations to determine the water source and age, or track the growth or decay of a chemical constituent throughout the network, perform fire flow analysis on the system to determine how it will behave under extreme conditions.
It uses the powerful Scenario Management features to mix and match a variety of "What If?" alternatives to the system.Create multiple sets of hydraulic, physical property, operations, initial setting, fire flow, cost, and water quality alternatives.This software is also running any number of scenarios by mixing and matching alternatives, then view and compares the results quickly and easily.After drawing the water distribution system it appears on the WaterCAD as shown in Figure 3.
Hydraulic analysis in WaterCAD is based on one of three equations, Hazen-Williams or Darcy-Weisbach or Manning Formula [8].

Water Quality Theory
The governing equations for WaterCAD water quality solver are based on the principles of conservation of mass coupled with reaction kinetics.When chlorinated water is introduced in the distribution system, the residual chlorine gradually dissipates.The exhaustion of the residual chlorine may then result in undesirable effects including biological decontamination in the distribution system.There are three factors influencing at chlorine consumption: 1) The consumption of chlorine as it reacts with organic and inorganic chemicals (e.g., ammonia, sulfides, ferrous iron, manganous ion, humic material).
2) The consumption of chlorine due to reactions with biofilms at the distribution pipe wall.
3) The consumption of chlorine in the corrosion process.The consumption of chlorine by reaction with organic and inorganic chemicals in the bulk aqueous phase is reasonably well defined.Certain of the chemicals react rather quickly and the reaction rate of the remainder is reasonably well characterized as a simple first-order decay process.
The first order kinetic model for the disappearance of residual chlorine due to reactions with materials in the aqueous phase at different residence times in the network expressed as:  where c is the chlorine concentration (mg/L) and k is the first-order decay constant (min −1 ).The residence time is defined as the pipe length divided by the mean flow velocity in the pipe.Integrating equation gives: where C(t) is the chlorine concentration (mg/L) at time t, c 0 is the initial chlorine concentration (mg/L), and it is the residence time in the pipe (min).The chlorine decay constant k is site specific and must be verified by field measurements.The k values vary with:  The water quality.
 The water temperature.
 The flow velocity.
 The pipe material.
 The area of contact with the pipe.
 The second-order rate process expressed as: ( ) ( ) where K' is a second-order rate constant.The use of this model was illustrated by Murphy (1985) [9].

Calibration Methodology
To ensure a good investment return and correct use of the models, the model should be capable of correctly simulating flow conditions encountered at the site.This is achieved by calibrating the models.A calibration involves the process of adjusting model characteristics and parameters so that the model's predicted flows and pressures match actual observed field data.Calibration of a water distribution model is a complicated task.There are many uncertain parameters that need to be adjusted to reduce the discrepancy between the model predictions and field observations of junction HGL and pipe discharges.Pipe roughness coefficients are often considered for calibration.However, there are many other parameters that are uncertain and affect junction HGL and pipe flow rate such as junction demand, operation status of pipes and valves.
Calibrating water distribution network models relies upon field measurement data, such as junction pressures, pipe flows, and water levels in storage facilities, valve settings, pump operating status (on/off), and pump speeds.Among all the possible field observation data, junction HGL and pipe flows are most often used to evaluate the goodness-of-fit of the model calibration.Other parameters, such as tank levels, valve settings, and pump operating status/speed are used as boundary conditions that are recorded when collecting a set of calibration observations of junction pressures and pipe flow rates.
Traditional calibration of a water distribution model is based on a trial-and-error procedure by which an engineer or modeller first estimates the values of model parameters, run the model to obtain a predicted pressure and flow, and finally compares the simulated values to the observed data.If the predicted data does not compare closely with the observed data, the engineer returns to the model, makes some adjustments to the model parameters, and calculates it again to produce a new set of simulation results.This may have to be repeated many times to make sure that the model produces a calibrated prediction of the water distribution network in the real world [10].

Scenario Management
Scenario management feature can dramatically increase the productivity in the "What If?" areas of modelling, including calibration, operations analysis, and planning.The main advantages of scenario management are:  A single project file makes it possible to generate an unlimited number of "What If?" conditions without becoming overwhelmed with numerous modelling files and separate results. The software maintains the data for all the scenarios in a single project, it can provide with powerful automated tools for directly comparing scenario results.Any set of results is immediately available at any time. The scenario/alternative relationship empowers to mix and match groups of data from existing scenarios without having to re-declare any data. The designer doesn't have to re-enter data if it remains unchanged in a new alternative or scenario, avoiding redundant copies of the same data.Inheritance also enables to correct a data input error in a parent scenario and automatically update the corrected attribute in all child scenarios [10].

Contour Map Manager
The Contour Map Manager contains the information required to generate contours for a calculated network, organized as follows [11]:  Minimum: lowest value to be included in the contour map.
 Maximum: highest value for which contours will be generated. Increment: step by which the contours increase.
The created contours will be evenly divisible by the increment, and are not directly related to the minimum and maximum values.For example, a contour set with 10 minimum, 20 maximum, and an increment of 3 would result in the following set: [12, 15, 18] not [10,13,16,19].

Hydraulic Analysis at Steady State Case
The study area is divided in two zones, zone one is located nearly from the NOZHA WTP, so it's in the circle of influence.Zone two is located nearly the EASTERN WTP.At average consumption rate = 3.04 m 3 /hr, hydraulic grade & pressure have been calculated.Figure 3 shows the pressures at the two zones by contour lines.
At maximum consumption rate = 3.87 m 3 /hr the hydraulic grade has been changed, Figure 4 shows the pressures at two zones by contour lines.

Hydraulic Analysis at Variable Consumption Rate (Extended State)
At this case the consumption is variable throughout the day, the multiplier multiplied in an average consumption rate.To adjust the pressure, two pump stations have been added to the hydraulic modelling analysis system to be as close as possible to the reality.Figure 5 and Figure 6 show the contour lines of pressure at two o'clock and 14 o'clock.There are some differences between existing hydraulic grade & software analysis because of the effect of other WTP in water distribution network.

Water Analysis
Residual chlorine concentration in the last stage of water treatment plant has been measured, as shown in Table 3 to know the residual chlorine after spreading in the network.
Unfortunately numbers shows the high percentage of residual chlorine which can lead to problems at human health especially the near from treatment plant.network in the study area at 2 p.m.
All chlorine measurements in water distribution network showed an increment and decrement in residual chlorine.The important reasons which cause these variables are:  The water distribution system is not separated so each WTP doesn't have its separate service area. The chlorine dose in WTP is not constant, it has been changed according to the quality of the water. Changing in consumption rate from winter to summer season, so the concentration changes with the change of consumption.For example if the chlorine dose is 2 mg/l = 2 g/m 3 and the flow rate during winter is 100.000m 3 /day, and in summer is 200.000m 3 /day, by compensation in this equation.

( )
Flow rate chlorine dose Chlorine dose at WTP 1000 convert to kg × =  The high rate of ammonia in the water at some areas leads to loss a percentage of free residual chlorine. There aren't stable systems for adding a chlorine dose in water treatment plants, all percentages of chlorine doses has been determined only by the laboratory results of water quality and this is not enough to make control on chlorine dose. There aren't specific link between water rate and chlorine dose in water treatment plants, such as SCADA system. The near areas from water treatment plants have a high percentage of residual chlorine and the contrast in remote areas.

Results of Residual Chlorine at WaterCAD after Adding the Actual Residual Chlorine in WTP
This part shows the results of residual chlorine in water distribution system before calibration by adding the same dose of residual chlorine which was added in the water treatment plant, to know the difference between the existing and software results.Table 4 shows the average of residual chlorine in water treatment plant, that's will spread in water distribution network (WDN).Figures 11-14 show the results of residual chlorine at WaterCAD software in the four months December, January, April and July.There are significant differences between the existing and software results.
To know the variation, scatter diagrams has been founded between the average existing residual chlorine and WaterCAD results before calibration shown in Figure 15.

Concentration of Residual Chlorine after Calibration at WaterCAD
The residual chlorine was recalculated in WaterCAD to calibrate the chlorine concentration to be as close as reality.The chlorine dose which has been added in the software model shown in the Table 5.As observed the percentage of chlorine dose is less than the existing and the residual chlorine concentration after calibration.Scatter diagrams have been found between the average existing residual chlorine and WaterCAD results after calibration shown in Figure 16.
The average residual chlorine concentration after calibration in WaterCAD software was as close as possible to the field results at the study area.The differences between the existing and modelling results are resulted from: -The contamination in the water network.
-Dead end & tree type of water network.
-The difference between software conditions and the other variables in fact.
-Temperature difference between existing and software condition.
-Differences between existing and calculated consumption rate by equations.

Injecting Chlorine Dose in Some Nodes at Weakness Areas
The World Health Organization (WHO) recommends that the residual chlorine concentration must not be less than 0.2 mg/l.There are some locations have residual chlorine concentration less than 0.2 mg/l, so operator can increase chlorine concentration at these nodes by inject chlorine at weakness residual chlorine area as shown in Figure 17.
Figure 18 shows that after injecting of chlorine dose = 0.5 mg/l at July month in five junctions at different lo-  cations, the minimum residual chlorine dose will increase to be at safety level.

Fire Flow Scenario
During firefighting time the flow rate in the water network will increase to be about 1189.5 m 3 /hr.Which will affect on residual chlorine concentration, this case called emergency case.By increasing the rate of water pumping at the treatment plant to avoid the problems of water consumption in water distribution system without increasing  of the chlorine dose, that's will lead to decreasing at residual chlorine concentration.Figure 19 shows the difference between the minimum residual chlorine at before and after fire case scenario.

Break in Pipe of Water Supply
This scenario could happen when some old pipes are broken and have a leakage, especially at areas which have old  pipes.Figure 20 shows the difference between the maximum and minimum residual chlorine after leakage without chlorine dose in critical junction and after inject 1 mg/l in nodes which are near from leakage points.These different scenarios systems and their solutions are more practical, effective and more economic because operator can control the concentration of chlorine at any junction in the water network, by addition chlorine at critical areas to ensure water quality.

Conclusions
1) To be more economical and useful, disinfection by low chlorine concentration can be used in water network after using software, because the operator can know the perfect doses of chlorine that will be added.That will lead to saving cost and use of a difference in cost in using other disinfection methods in the water treatment plants.
2) Decline in the concentration of free residual chlorine at some locations in water network is related to: * The pipes material and the lifespan of pipes.* Water treatment plants are only the main source of the chlorine dose in water network; other sources can be founded such as pumping station to ensure water quality at weakness areas.4) After calibration has been made by WaterCAD, chlorine dose is decreased while at the same time residual chlorine is at a safety level, and this leads to saving costs.
5) The sudden rise of water consumption such as fire flow case without increasing chlorine dose, leads to decreasing in residual chlorine concentration at many different locations.
6) The lack of separation of each distribution network leads to inequality of residual chlorine concentration in different water network locations.

Figure 1 .
Figure 1.Map showing the location of the study area.

Figure 2 .
Figure 2. Lab samples locations and service area.

Figure 3 .
Figure 3. Contour lines of pressure at study area (average consumption).

Figures 7 -
10 show the residual chlorine in the water distribution

Figure 4 .
Figure 4. Contour lines of pressure at study area (maximum consumption rate).

Figure 7 .
Figure 7.The residual chlorine at water distribution network during (Dec).

Figure 8 .
Figure 8.The residual chlorine at water distribution network during (Jan).

Figure 9 .
Figure 9.The residual chlorine at water distribution network during (Apr).

Figure 10 .
Figure 10.The residual chlorine at water distribution network during (Jul).

Figure 11 .
Figure 11.The residual chlorine by WaterCAD before calibration at Dec.

Figure 12 .
Figure 12.The residual chlorine by WaterCAD before calibration at Jan.

Figure 13 .
Figure 13.The residual chlorine by WaterCAD before calibration at Apr.

Figure 14 .
Figure 14.The residual chlorine by WaterCAD before calibration at Jul.

Figure 15 .
Figure 15.Standard deviation of residual chlorine before calibration.

Figure 16 .
Figure 16.Standard deviation of residual chlorine after calibration.

Figure 17 .
Figure 17.Injecting of chlorine dose in some nodes.

Table 2 .
Sheyakhat population at eastern section.

Table 3 .
Residual chlorine concentration in water treatment plant.