Next Article in Journal
Protection of the Liuzheng Water Source: A Karst Water System in Dawu, Zibo, China
Next Article in Special Issue
Robustness and Water Distribution System: State-of-the-Art Review
Previous Article in Journal
Looking Beyond Energy Efficiency: An Applied Review of Water Desalination Technologies and an Introduction to Capillary-Driven Desalination
Previous Article in Special Issue
The Impacts of Spatially Variable Demand Patterns on Water Distribution System Design and Operation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Comparative Study on a Hydraulic and Water-Quality Analysis Method for Determining Rechlorination Injection Points for a Water-Supply Network

1
School of Civil, Environmental and Architectural Engineering, Korea University, Seoul 02841, Korea
2
Research Center for Disaster Prevention Science and Technology, Korea University, Seoul 02841, Korea
3
Department of Civil Engineering, The University of Suwon, Gyeonggi-do 18323, Korea
*
Author to whom correspondence should be addressed.
Water 2019, 11(4), 697; https://doi.org/10.3390/w11040697
Submission received: 22 February 2019 / Revised: 21 March 2019 / Accepted: 2 April 2019 / Published: 4 April 2019

Abstract

:
Water supply facilities are vulnerable to extreme weather events, such as droughts and floods. To establish a sustainable solution that resists accidents and disasters, a distributed system is required. To supply high-quality tap water using the existing water-supply network, rechlorination facilities must be installed to secure residual chlorine at the pipe end. In this study, a process is developed to determine the injection points and dosages of rechlorination using the latest pressure-driven analysis. The method was compared to the results of demand driven analysis methods. The proposed model was applied to P City in Korea to draw results. A detailed evaluation was performed to study how water pressure head and demand-based hydraulic and water quality analysis results impact the injection points and dosages of rechlorination. Thus, the existing demand-based model shows significant spatial deviations in the pressure head in the presence of water pressure drops, which subsequently lead to over-estimation of chlorine injection dosages for maintaining the concentration of residual chlorine. However, the proposed model involves a numerically validated theory and draws more reasonable results for hydraulic, water quality, and rechlorination dosages. The proposed model can be used as a decision-making tool based on hydraulic analysis for the supply of water of a stable quality.

1. Introduction

Global climate change has resulted in frequent floods and extreme droughts, and accelerates urbanization. Additionally, aging infrastructure has led to a higher probability of man-made accidents. In this context, the maintenance of a city’s sustainability functions against natural and man-made disasters is an important issue. Water-supply networks (WSN) provide the infrastructure that supplies consumers with treated potable water, which must be provided with the correct volume, pressure, and quality. Thus, any occurrence of poor water quality directly impacts human life. A variety of water-quality issues in pipelines can thus arise. Direct damage to pipelines caused by man-made accidents, terrorist attacks, earthquakes, etc. can cause quality issues. Additionally, there can be shortages of raw water because of drought and residual chlorine concentrations (i.e., disinfection). This is related to abnormally high raw-water temperatures. In this aspect, maintaining an adequate residual chlorine concentration in a pipeline system is essential for secondary contamination control to inhibit microbial regrowth. Generally, highly concentrated residual chlorine occurs in the vicinity of water purification plants, where chlorine disinfection takes place. This is because of the pipeline system of a branched network and the characteristics of volatile chlorine disinfectant. Conversely, the chlorine level decreases at pipe ends. One of the methods to address this issue is to convert a branched network to a looped network or to install a rechlorination facility in the system.
In the former case, the scale of pipeline infrastructure that would require installation to achieve equalization effects is massive and would incur extreme costs. The latter case of installing a rechlorination facility incurs lower costs and is therefore widely used in the field (K-water [1]). Studies on the installation and operation of rechlorination facilities originally focused on minimizing the disinfectant mass at chlorine-injection facilities. Boccelli et al. [2] conducted a water-quality analysis in pipelines using the primary decay reaction of chlorine, aiming to minimize the disinfectant mass for predetermined rechlorination injection points. As a result, an adequate water quality was achieved with less disinfectant mass, compared to the existing method of injecting disinfectants only into water sources.
Munavalli and Mohan Kumar [3], Tryby et al. [4], and Propato and Uber [5] proposed a design and operational method to maximize the effects of rechlorination using a single-purpose optimization method. Munavalli and Mohan Kumar [3] used a genetic algorithm as an optimization method and set the minimization of the residual chlorine concentration as an objective function to determine the injection concentration at predetermined injection points. Tryby et al. [4] attempted to minimize the disinfectant mass for injection and applied mixed-integer linear programming to determine the rechlorination injection points. Propato and Uber [5] achieved optimization by setting an objective function as the minimal total of the squared deviations of residual concentrations to assess the effects of rechlorination, focusing on equalization. More recently, Yoo et al. [6] used a harmony search algorithm, a meta-heuristic optimization method, to study issues requiring the installation of a rechlorination facility to determine the rechlorination points for minimizing the disinfectant mass.
Prasad et al. [7], Ostfeld and Salomons [8], Islam et al. [9], Choi et al. [10], and other research teams evaluated residual chlorine by addressing not only residual chlorine concentrations, but also the hydraulic performance of the water-supply system or additional water-quality assessment factors. Prasad et al. [7] proposed a multi-purpose optimization method that minimized the total chlorine mass for injection (i.e., water-quality indicator) and simultaneously maximized the water-supply mass (i.e., hydraulic indicator) in an attempt to determine the location and operational methods (i.e., scheduling) of a rechlorination facility. Ostfeld and Salomons [8] optimized the operational method (i.e., hydraulic indicator) of pumps and the operational method (i.e., water quality indicator) of a rechlorination facility while minimizing the operational and construction costs of a water supply. Islam et al. [9] proposed a water-quality index (WQI) to provide a single value to reflect microbial, chemical, and aesthetic water quality standards. Using the WQI value as a basis, another study was conducted to determine an optimal rechlorination mass for injection points. Choi et al. [10] analyzed the effects of residual chlorine in relation to rechlorination for water-quality improvement in the supply process. They set three assessment indices of decreasing residual chlorine concentrations between a water purification plant and pipe ends. They also addressed the decrease in total chlorine mass for injection and the decrease in disinfection by-product generation.
When looking at past studies on rechlorination, most methodologies obtained residual chlorine concentration results through hydraulic and water-quality analyses and determined an optimal design or operational method based on the results. Most extant studies, however, were conducted using EPANET2 (Rossman [11]), a legacy model for a demand-driven analysis (DDA). With DDA, workers perform hydraulic analysis under the premise that a given water volume is supplied under all circumstances, and all water volumes supplied to the network are considered as known values. Therefore, it cannot directly simulate unplanned events, such as water-supply decreases caused by water pressure losses. This has generated distorted hydraulic statistics for water-quality analysis results. One method used to address this issue is a pressure-driven analysis (PDA). PDA provides a supply amount (i.e., total of the consumption amount and leak amount) as an unknown value, affected by residual heads and water pressures. When looking at a PDA from a hydraulic perspective in the water-supply pipeline system, studies have been conducted by Giustolisi et al. [12,13], Wu et al. [14], Lee et al. [15], Yoo et al. [16], etc. These studies suggested that DDA could not precisely mathematically simulate the event of pressure decreases and that the PDA could give more realistic results. From a water-quality analysis perspective, Seyoum and Tanyimboh [17] and Seyoum et al. [18,19,20] conducted studies on how a PDA could affect water-quality analysis. These studies showed that, for DDA, the assumption of supplying all demands, despite actual water-pressure decreases, could lead to an underestimation of water age. Thus, it was suggested that a decrease in the flow rate in the pipe caused by a pressure decrease could lead to relative increases in water age from a PDA perspective. However, these studies only suggested the general tendency of hydraulic and water-quality analysis from a PDA, implying that residual chlorine concentration levels could change according to actual flow rates in the pipes. Thus, there has yet to be any analysis on how hydraulic and water-quality analysis results impact the final chlorine injection amount and the selection of rechlorination injection points. Furthermore, a pressure-driven hydraulic analysis has been pre-installed in several commercial programs (WaterGEMs, Bentley [21]) and source codes (EPANET3 [22]). However, the development of a model for repeatedly conducting water-quality analysis and evaluating its result based on a pressure-driven hydraulic analysis has not been proposed.
An analysis of water quality changes in WSN is a direct extension of a hydraulic WSN analysis and can be used effectively. That is, a water-quality analysis involves a process of calculating chemical combinations, concentration changes, water ages, etc., in WSN, based on the results of hydraulic analysis. The flow direction, which indicates the flow rate in pipelines and flow routes of water within pipelines, is used to determine the combination and water age and other hydraulic characteristics that affect the transfer of pollution sources and concentration decrease. Most hydraulic models use one-dimensional flows and reaction transfers, which express the concentration in a particular pipe as a function of the distance and time to predict concentration changes in ingredients caused by transfers within WSN and calculate the formation of reactants and concentration decreases. A typical example of volume and wall reactions is the disinfection reaction, which is most commonly used in hydraulic analyses. When water leaves a reservoir and flows through pipelines, it generates various chemical reactions with the disinfectant. Chlorine reacts to natural organic matters within the volume, to the bio-membranes within the pipe walls, and to the pipe material, causing decay. Generally, chlorine decay models are classified into a first-order decay and power-law decay (nth order) of second or higher order, according to the form of the attenuation equation (Carrico and Singer [23]; Yoo et al. [6]). The primary chlorine attenuation model for the volume and wall reaction is commonly used in studies for its simplicity and availability (Courtis et al. [24]; Powell et al. [25]; Tamminen and Ramos [26]). Therefore, the first-order chlorine decay model, which accounts for bulk and wall effects, is used in this study.
In this paper, a model is developed for determining rechlorination injection points to maintain an adequate residual chlorine concentration level in the system using a PDA, including a water-quality analysis. Furthermore, results of the existing DDA on rechlorination injection points and dosages are obtained and compared with that of the PDA results. Lastly, the proposed model is applied to an actual WSN under normal and abnormal conditions, and a detailed evaluation is performed to study how hydraulic and water-quality analysis results affect the determination of rechlorination injection points.

2. Materials and Methods

2.1. Definition of Abnormal Condition in WSN

Abnormal situations that can occur in a water supply network system are very diverse. For instance, abnormal situations include pipeline breaks, shutdown of the pump, and reduction of the water level (reservoir, reservoir) due to drought from the hydraulic point of view. In addition, the exceptional variation of the effluent chlorine concentration at the water treatment plant can be considered as an abnormal condition of the water quality factor. The chlorine concentration is affected by the hydraulic factors (pipe flow, nodal water pressure head, etc.) in a water supply network. However, in the case of DDA, the hydraulic factors are calculated assuming that all the demands of the nodal point are satisfied regardless of the nodal water pressure head, so that the validity of the hydraulic analysis results cannot be secured in the event of abnormal conditions. Therefore, the reliability of the water quality analysis results also cannot be guaranteed. Therefore, in this study, we presented the results of an abnormal situation scenario when the water level of the reservoir decreases, in order to compare the DDA and the PDA.

2.2. K-NriskQ

K-NRisk (K-water [27]) is a pressure-driven demand simulation model that analyzes areas with water-supply suspensions and water-pressure shortages under abnormal conditions to quantify supply availability for demand. It typically performs the function of pressure-driven leak simulation for a precise hydraulic analysis simulation, according to changes in the leak amount. Its major difference from EPANET3, which recently became open source, is that it automatically performs a sensitivity analysis of input factors used for a PDA. It quantifies the uncertainty innate in simulation models and provides a function for novice users to easily operate typical analysis scenarios of abnormal situations. In particular, it aims to develop a user-friendly decision-making tool for analysis efficiency improvement, such as for local water-supply system operation efficiency boost projects and water-flow rate improvement projects, which require a WSN hydraulic analysis.
K-NRiskQ (K-water [28]) is an advanced model that focuses on the water-quality analysis function of K-NRisk. It was developed based on a parameter sensitivity analysis function and a simulation function for residual chlorine equalization. Parameter sensitivity analysis was designed to intuitively understand the changes in network water quality according to changes in the pipe wall and body coefficients and to simulate water-quality changes in relation to changes (e.g., minimum, average, maximum temperature) set in accordance with the water temperature of body coefficients in a bottle test. Residual chlorine equalization simulations were designed to provide an automatic comparative review of rechlorination injection, manually performed, to determine rechlorination injection points for boosting water-quality satisfaction. Our study is based on the residual chlorine equalization simulation function designed in K-NRiskQ. Table 1 shows the available functions of the major hydraulic and water-quality analysis models for WSN. It suggests that the water-quality analysis function, via PDA, is more valid in K-NRiskQ than other tools.

2.3. Methodology for Determining Rechlorination Points

2.3.1. Objective Function and Constraints

The objective function used in this model is as shown in Equation (1). It aims to minimize the total disinfectant mass additionally injected. The determinant variables of the model are the locations of the rechlorination points, the disinfectant concentration, and the mass injected into each rechlorination point:
minimize   f = i = 1 n n M i
where f indicates the total disinfectant mass (kg/d), Mi is the disinfectant mass at joint i (kg/d), and nn is the total number of joints. The constraints in this model can be mainly classified into two categories. The first is the residual chlorine concentration at each joint. As shown in Equation (2), the concentration at all joints must be between the proposed minimum and maximum values:
ClminCliClmax (i = 1, …, nn)
where Cli indicates the residual chlorine concentration at joint i (mg/L), and Clmin and Clmax are the minimum and maximum residual chlorine concentrations allowed, respectively (mg/L). The second category is the number of rechlorination facilities to be installed. The number of facilities to be installed is predetermined and set as constraints for this model. Therefore, this model supports different numbers of facilities in scenarios:
N = {1, …, nn}
where N denotes the number of possible designed boosters.

2.3.2. Methodology Procedure

In this study, two scenarios of normal and abnormal conditions were set and analyzed to study residual chlorine. A greedy algorithm was applied to explore rechlorination injection points. A greedy algorithm is one of the simplest methods of making the best optimal choice at an early stage. Whereas it does not always guarantee an optimal result, it has the advantage of promptly obtaining a value close to an optimal solution. Scenarios are selected from either a normal or an abnormal situation. Figure 1 shows the progress diagram of this study. The purpose of the greedy algorithm is to keep the residual chlorine standard above the minimum and below the maximum for all target joints. When the standard deviates from this acceptable range, it is configured to first select ones with the least error.
The procedure for determining rechlorination injection points is as follows:
(1).
Check the break points and water sources of the target SWN where chlorine injection points are set.
(2).
Distinguish between normal and abnormal scenarios, and detail the scenario.
(3).
Check the hydraulic and water-quality analysis results from the current chlorine injection for normal and abnormal scenarios via PDA and DDA analyses.
(4).
Set a candidate group for chlorine injection. The break points in the target area in Figure 2 were set as the candidate group for chlorine injection for this study. However, this was set as configurable by a model user.
(5).
Set quality constraints and the number of rechlorination injection points to determine the optimal conditions for residual chlorine equalization. Set the maximum and minimum values and the number of rechlorination injection points, which are the acceptable range for residual chlorine concentration.
(6).
Perform a greedy algorithm to satisfy the concentration conditions for residual chlorine set by a user and calculate the number of chlorine injection points within the acceptable range of the residual chlorine concentration and dosage. To confirm the constraints, a hydraulic analysis is required by trial calculation. After the calculation is complete, propose the rechlorination injection method as a result.

3. Application and Result

3.1. Subject WSN

The region used as a case study was the WSN of P City in Korea, which was first studied by Yoo et al. [6]. P City is a small-to-medium-sized city (826 km2) that supplies water to 119,750 residents. Its water-supply system comprises approximately 85-km-long water transmission pipes and 773-km-long water distribution pipes. Eight water distribution locations are used to supply 54,020-m3-per-day of water. Figure 2 shows the structure of the water transmission pipes of the water supply system in P City. P City is expected to have additional demand for 10,000 CMD of water at point GN-D (the circled area in Figure 2), owing to changes in water sources and additional urban development. According to the current Korean standards, the water quality of residual chlorine must range between 0.1 mg/L and 4.0 mg/L (KWWA [29]). Moreover, it is regulated to be 0.4 mg/L or higher when polluted by, or when it has the potential to be polluted by, pathogenic microbes. Therefore, a WSN operator in P City strives to maintain water quality from 0.4 mg/L to 4.0 mg/L, even when additional water is needed, by installing rechlorination injection facilities downstream of point GN-D.
The result of the maintenance and water-quality review for maintaining the current chlorine injection points and dosage, according to the additional water supply demand, is shown in Table 2.
The result of a DDA simulation for additional water-supply demand showed that 21 break points did not satisfy the lower required limit of 0.4 mg/L. In the case of YP-E and SBGI-J, both having long travel times, the expected values of chlorine concentration were 0.05 mg/L and 0.01 mg/L, respectively. They did not meet the legal minimum standard as well. Additionally, the average residual chlorine concentration of major break points was also observed to be 0.30 mg/L below the minimum standard. Therefore, to supply safe tap water through pipe ends, it is necessary to install rechlorination facilities to secure residual chlorine.
According to the PDA simulation result, it was confirmed that nine break points did not meet the lower limit of 0.4 mg/L requested by the operator. The break point with the lowest expected value was SBGI-J with 0.06 mg/L, which was below the legal minimum standard. The average residual chlorine concentration of breakpoints was 0.44 mg/L, satisfying the minimum standard, differing from the DDA analysis result.
The DDA and PDA results of the residual chlorine analysis showed disparity, and the number of breakpoints satisfying the acceptable standard showed a difference, because, during calculation, it was assumed that all basic demands were met with the corresponding supply, despite the fact that the available amounts should decrease, owing to water-pressure decreases in real situations. Such an assumption results in increased deviation in the water pressure head and residual chlorine concentration between break points. Additionally, water-pressure head decreases lead to decreases in the available supply in a DDA. Regarding break point YB-D, the simulation showed that the available amount of supply was 2692 CMD, compared to the basic demand of 3426 CMD, indicating that approximately 79% of demand could be met. Such a decrease in the amount of water supply can cause a relative increase in the residual chlorine concentration.

3.2. Setting Analysis Scenarios

This study aims to identify rechlorination points that can address the current below-standard residual chlorine concentrations under normal conditions and to determine rechlorination injection points when the head of the water source (i.e., main source) is operated at a 5-m-lower point because of abnormal conditions. Unlike the study of Yoo et al. [6], our reason for additionally accounting for an abnormal situation is that many studies have raised an issue regarding the adequacy of the DDA results under abnormal conditions. Thus, because water-quality analysis results are used as basic input data, it is also necessary to conduct a water-quality analysis for abnormal conditions. In this study, the model was operated with 2, 3, and 4 rechlorination injection points, and the results were compared and analyzed. A total of 40,000 repetitive trial calculations were conducted for each model.
The HOR used in the DDA for demand is an equation proposed by Wagner et al. [30]. The minimum water pressure for supplying break points was set to 0 m, and the ideal water pressure for 100% supply-to-demand was set to 15 m. Yoo et al. [6] suggested the bulk and wall decay coefficient of the P-city network based on a bottle test and the results of previous studies [31,32], separately. Reflecting on the results of the literature, the bulk and wall decay coefficients entered for modelling the residual chlorine concentration in pipes were set to −0.5760 d-1 and −0.1 d−1, respectively.

3.3. Application and Results

Table 3 shows the results of hydraulic and water quality analyses for abnormal conditions, assuming that the level of the water source decreased by 5 m. The DDA simulation result under abnormal conditions showed similar trends to those of normal conditions. However, it was confirmed that the number of points below the residual chlorine standard significantly decreased from 9 to 3, according to the PDA simulation result. This leads to a situation where the available amount of water supply in the system decreases as the amount of water source decreases, and the residual chlorine concentration in the system increases, as opposed to normal conditions.
Finally, when looking at the hydraulic and water-quality analysis results for both normal and abnormal conditions prior to additional rechlorination, it is shown that the change in water quality is insignificant, despite the decrease in the number of heads of water sources, because the available amount of supply is the same in the DDA. However, it was shown that the change in the available amount of supply greatly affected the water-quality analysis result in the case of PDA. Eventually, from a rechlorination injection perspective, it can be expected that the number of total rechlorination injections will decrease during abnormal conditions. The rechlorination injection points may appear different according to the system status, deviations by break points, etc.
Figure 3 diagrams an optimal model operation under normal and abnormal conditions where there is a need to install 2, 3, or 4 rechlorination injection facilities in the system. It shows the total rechlorination injection dosage, which is the objective function of the optimization methodology. When applying the PDA under normal conditions, the result shows that some rechlorination dosage increases for determining two points. However, it shows a decreasing tendency in the total injection dosage required, because the number of installation points increases under both normal and abnormal conditions and for both DDA and PDA analyses. Because the greedy algorithm does not guarantee an optimal solution for all issues, the optimization performance of the proposed methodology in this study was evaluated and compared to the result from the methodology proposed by Yoo et al. [6]. Figure 4 shows the result of the DDA under normal conditions using the harmony search (HS), a metaheuristic optimization methodology. Figure 3a and Figure 4 provide the results of the greedy algorithm and HS for the same scenario and the hydraulic analysis method. When comparing the results, the HS result is observed to be superior (lesser injection amount) in terms of the final total injection amount. However, the greedy algorithm also showed an objective function result at almost the same level. The purpose of this study is to compare the result of applying DDA and PDA and normal and abnormal conditions in determining the rechlorination injection amount, rather than to obtain the results of better injection amounts. Therefore, whereas the greedy algorithm does not show superior results compared to the HS optimization method, it was confirmed that it can provide a close-to-optimal result.
Table 4 and Table 5 provide the results of operating an optimal model under normal and abnormal conditions when needing to install 2, 3, or 4 rechlorination injection facilities in the system. These are the results for the expected residual chlorine concentrations at break points when applied to the result of the optimal injection points and concentrations. It was confirmed that the optimization results of all scenarios satisfied the constraints, which is the acceptable range of residual chlorine concentration proposed. It is notable that the total rechlorination injection amount indicated as the result of the DDA was relatively greater than that of the PDA. In terms of the proportion, it was between 11.2% and 41.8% under normal conditions and 9.8% and 25.5% under abnormal conditions. This occurs because of the decrease in the water supply and pressure head shown in the PDA results, leading to an increase in the residual chlorine concentration, which then causes a decrease in the additional rechlorination injection amounts. The standard deviation indicates the distribution of the concentration between break points and shows smaller deviations compared to the PDA result. This suggests that increases in the deviation between break points in the hydraulic analysis, caused by the DDA result, assuming an amount of supply for all demand, affect the water quality analysis.
Additionally, decreases in the supply are more evident under abnormal conditions, which causes additional decreases in the water pressure head than under normal conditions, leading to lower rechlorination injection amounts needed to satisfy the acceptable range of the residual chlorine concentration. This result can help intuitively understand the changes in the injection amount under various abnormal conditions. For instance, one can assume that the available amount of supply is infinite or sufficient and that there is a change in the system status, such as pipe damage or increases in demand. In this case, despite the same abnormal conditions, the demand in the system or leakage amount is expected to increase, contrary to the decrease in the amount of water source assumed in the scenario. Thus, there can be additional decreases in the residual chlorination concentration. Certainly, other external factors (e.g., system operation and distinct characteristics) can also affect the results. However, if hydraulic and water-quality analysis results can be obtained while reflecting the characteristics of the system based the method proposed in this study, it can help with decision-making to address increases in rechlorination injection amounts and water quality issues.
Figure 5 and Figure 6 provide network diagrams with locations of the optimal rechlorination injection points for normal and abnormal conditions. The injection points were selected focusing on injection points and end points of the system in the DDA result, whereas the center of the system or the points where flow rates split were selected in the PDA result. This result suggests that, because the deviations between break points for the residual chlorine concentration were relatively high in the DDA, it was necessary to select the injection points or end points as rechlorination points, because they had the minimum concentrations within the acceptable range of the concentration. However, in the case of PDA, the deviations were relatively low, and the amount of the additional chlorine injection was also low, producing an insignificant change. However, the flow-rate split, triggering a change of the uniform concentration as determined by the optimal point.
Table 6 and Table 7 provide the residual chlorine concentration levels, shown in each break point when rechlorination facilities are installed at two locations. As suggested, the deviation in the residual chlorine concentration in the DDA result is greater than that of the PDA result. However, it satisfies all constraints for an acceptable residual chlorine concentration for all analyses.

4. Conclusions

A methodology was proposed for determining the location of rechlorination facilities and dosages, using an in-depth comparative analysis, performed on the results of PDA and DDA. A simple greedy algorithm was used in the optimal solution calculation model, and the model was operated under four different scenarios, with the settings of normal and abnormal conditions. The model was operated for 2, 3, and 4 chlorine injection points, and the results showed that all break points satisfied the minimum and maximum residual chlorine concentration standard. Whereas the greedy algorithm did not provide a result superior to the complex metaheuristic optimization method, it was close enough and was confirmed to be an adequate model to serve the purpose of this study, which entailed the comparison of hydraulic analysis results. In terms of determining rechlorination injection points through DDA and PDA, we confirmed that the DDA model showed significant spatial deviations in the presence of water pressure head decreases, leading to an overestimation of injection amounts. However, the proposed model, which involved PDA, showed that it had less rechlorination injection amounts determined to meet the acceptable range, compared to the DDA model. Because the PDA model is applied as a numerically validated theory, it appears to offer a more rational hydraulic and water-quality analysis based on rechlorination injection amounts. Furthermore, the optimization of rechlorination injection was performed with a scenario of abnormal conditions. Thus, it was confirmed that an increase or decrease in the amount of supply of the overall system, caused by abnormal conditions, affected the amount of the rechlorination injection. Thus, the results could vary depending on whether the supply-flow rate of the system decreased or increased. When the model proposed in this study was used, more realistic hydraulic and water quality analysis results were obtained, and the tendency of the water-quality results under normal and abnormal conditions was confirmed. This model can be used as a hydraulic analysis-based decision-making tool to supply water of a stable quality.
In the future, an extended period simulation should be conducted to ensure an adequate level of residual chlorine throughout the entire period, rather than a single-period simulation based on average demand. Additionally, further studies will need to be applied to more diverse scales and complex network systems (branched and looped type systems) and their results should be analyzed. Furthermore, changes in water-quality analysis results should be examined in more variety, and there should be studies on the development of indices to be used in decision-making. Furthermore, the water age and pressure head can be used as an additional design and operational factor for the overall evaluation of WSN.

Author Contributions

S.M.L. and J.H.K. conceived and designed the original idea of proposed method; D.G.Y. and S.M.L. carried out survey of previous studies; S.M.L., D.G.Y., and H.M.L. analyzed the data; S.M.L. and D.G.Y. wrote the paper.

Funding

This work was supported by a grant from the National Research Foundation (NRF) of Korea, funded by the Korean government (MSIP) (No. 2016R1A2A1A05005306).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. K-water. Final Reports of Smart Water City Projects; K-water: Daejeon, Korea, 2016.
  2. Boccelli, D.L.; Tryby, M.E.; Uber, J.G.; Rossman, L.A.; Zierolf, M.L.; Polycarpou, M.M. Optimal scheduling of booster disinfection in water distribution systems. J. Water Resour. Plan. Manag. 1998, 124, 99–111. [Google Scholar] [CrossRef]
  3. Munavalli, G.R.; Kumar, M.M. Optimal scheduling of multiple chlorine sources in water distribution systems. J. Water Resour. Plan. Manag. 2003, 129, 493–504. [Google Scholar] [CrossRef]
  4. Tryby, M.E.; Boccelli, D.L.; Uber, J.G.; Rossman, L.A. Facility location model for booster disinfection of water supply networks. J. Water Resour. Plan. Manag. 2002, 128, 322–333. [Google Scholar] [CrossRef]
  5. Propato, M.; Uber, J.G. Linear least-squares formulation for operation of booster disinfection systems. J. Water Resour. Plan. Manag. 2003, 130, 53–62. [Google Scholar] [CrossRef]
  6. Yoo, D.G.; Lee, S.M.; Lee, H.M.; Choi, Y.H.; Kim, J.H. Optimizing Re-Chlorination Injection Points for Water Supply Networks Using Harmony Search Algorithm. Water 2018, 10, 547. [Google Scholar] [CrossRef]
  7. Prasad, T.D.; Walters, G.A.; Savic, D.A. Booster disinfection of water supply networks: Multiobjective approach. J. Water Resour. Plan. Manag. 2004, 130, 367–376. [Google Scholar] [CrossRef]
  8. Ostfeld, A.; Salomons, E. Conjunctive optimal scheduling of pumping and booster chlorine injections in water distribution systems. Eng. Optim. 2006, 38, 337–352. [Google Scholar] [CrossRef]
  9. Islam, N.; Sadiq, R.; Rodriguez, M.J. Optimizing booster chlorination in water distribution networks: A water quality index approach. Environ. Monit. Assess. 2013, 185, 8035–8050. [Google Scholar] [CrossRef]
  10. Choi, T.; Lee, D.; Bae, C.; Moon, J. Analysing the Effect of Residual Chlorine Equalization for Water Quality Improvement in Water Distribution System. J. Korean Soc. Water Wastewater 2016, 30, 587–596. [Google Scholar] [CrossRef]
  11. Rossman, L. EPANET User’s Manual; United States Environmental Protection Agency (EPA): Washington, DC, USA, 2000.
  12. Giustolisi, O.; Kapelan, Z.; Savic, D.A. Extended period simulation analysis considering valve shutdowns. J. Water Resour. Plan. Manag. 2008, 134, 527–537. [Google Scholar] [CrossRef]
  13. Giustolisi, O.; Savic, D.A.; Kapelan, Z. Pressure-driven demand and leakage simulation for water distribution networks. J. Hydraul. Eng. 2008, 134, 626–635. [Google Scholar] [CrossRef]
  14. Wu, Y.W.; Wang, R.H.; Walski, T.M.; Yang, S.Y.; Bowdler, D.; Baggett, C.C. Extended global-gradient algorithm for pressure-dependent water distribution analysis. J. Water Resour. Plan. Manag. 2009, 135, 13–22. [Google Scholar] [CrossRef]
  15. Lee, H.M.; Yoo, D.G.; Kang, D.; Jun, H.; Kim, J.H. Uncertainty quantification of pressure-driven analysis for water distribution network modeling. Water Sci. Technol. Water Supply 2016, 16, 599–610. [Google Scholar] [CrossRef]
  16. Yoo, D.G.; Chang, D.E.; Song, Y.H.; Lee, J.H. Optimal Placement of Pressure Gauges for Water Distribution Networks Using Entropy Theory Based on Pressure Dependent Hydraulic Simulation. Entropy 2018, 20, 576. [Google Scholar] [CrossRef]
  17. Seyoum, A.G.; Tanyimboh, T. Pressure dependent network water quality modelling. Proc. ICE-Water Manag. 2014, 167, 342–355. [Google Scholar] [CrossRef]
  18. Seyoum, A.G.; Tanyimboh, T.; Siew, C. Comparison of demand driven and pressure dependent hydraulic approaches for modelling water quality in distribution networks. In Proceedings of the 11th International Conference on Computing and Control for the Water Industry, Exeter, UK, 5–7 September 2011. [Google Scholar]
  19. Seyoum, A.G.; Tanyimboh, T.T.; Siew, C. Assessment of water quality modelling capabilities of EPANET multiple species and pressure-dependent extension models. Water Sci. Technol. Water Supply 2013, 13, 1161–1166. [Google Scholar] [CrossRef]
  20. Seyoum, A.G.; Tanyimboh, T.T. Integration of hydraulic and water quality modelling in distribution networks: EPANET-PMX. Water Resour. Manag. 2017, 31, 4485–4503. [Google Scholar] [CrossRef]
  21. WaterGEMs User’s Manual. Available online: https://www.bentley.com/en/products/product-line/hydraulics-and-hydrology-software/watergems (accessed on 4 April 2019).
  22. EPANET3. Available online: https://github.com/OpenWaterAnalytics/epanet-dev (accessed on 4 April 2019).
  23. Carrico, B.; Singer, P.C. Impact of booster chlorination on chlorine decay and THM production: Simulated analysis. J. Environ. Eng. 2009, 135, 928–935. [Google Scholar] [CrossRef]
  24. Courtis, B.J.; West, J.R.; Bridgeman, J. Chlorine demand-based predictive modeling of THM formation in water distribution networks. Urban Water J. 2009, 6, 407–415. [Google Scholar] [CrossRef]
  25. Powell, J.C.; Hallam, N.B.; West, J.R.; Forster, C.F.; Simms, J. Factors which control bulk chlorine decay rates. Water Res. 2000, 34, 117–126. [Google Scholar] [CrossRef]
  26. Tamminen, S.; Ramos, H.; Covas, D. Water supply system performance for different pipe materials Part I: Water quality analysis. Water Resour. Manag. 2008, 22, 1579–1607. [Google Scholar] [CrossRef]
  27. K-water. K-NRisk User’s Manual; K-water: Daejeon, Korea, 2017.
  28. K-water. K-NRiskQ User’s Manual; K-water: Daejeon, Korea, 2018.
  29. Korea Water Works Association (KWWA). Standards for Water Supply Facilities; Korea Water Works Association (KWWA): Seoul, Korea, 2010.
  30. Wagner, J.M.; Shamir, U.; Marks, D.H. Water distribution reliability: Simulation methods. J. Water Resour. Plan. Manag. 1988, 114, 276–294. [Google Scholar] [CrossRef]
  31. Al-Jasser, A.O. Chlorine decay in drinking-water transmission and distribution systems: Pipe service age effect. Water Res. 2007, 41, 387–396. [Google Scholar] [CrossRef]
  32. Hallam, N.B.; West, J.R.; Forster, C.F.; Powell, J.C.; Spencer, I. The decay of chlorine associated with the pipe wall in water distribution systems. Water Res. 2002, 36, 3479–3488. [Google Scholar] [CrossRef]
Figure 1. Proposed model operation procedure.
Figure 1. Proposed model operation procedure.
Water 11 00697 g001
Figure 2. Target WSN (P City, Republic of Korea) (Yoo et al. [6]).
Figure 2. Target WSN (P City, Republic of Korea) (Yoo et al. [6]).
Water 11 00697 g002
Figure 3. Changes in optimal injection amount by scenario.
Figure 3. Changes in optimal injection amount by scenario.
Water 11 00697 g003
Figure 4. Changes in optimal rechlorination injection amount for 2, 3, and 4 points, using harmony search (Yoo et al. [6]).
Figure 4. Changes in optimal rechlorination injection amount for 2, 3, and 4 points, using harmony search (Yoo et al. [6]).
Water 11 00697 g004
Figure 5. Optimal points under normal conditions.
Figure 5. Optimal points under normal conditions.
Water 11 00697 g005aWater 11 00697 g005b
Figure 6. Optimal points in abnormal conditions.
Figure 6. Optimal points in abnormal conditions.
Water 11 00697 g006aWater 11 00697 g006b
Table 1. Function comparison between major hydraulic/water-quality analysis models.
Table 1. Function comparison between major hydraulic/water-quality analysis models.
FunctionEPANET2EPANET3K-NRiskK-NRiskQ
Pressure driven analysis (head outflow relationship)X *O *OO
Quantification through leakage performance indexXXOO
Scenarios for multiple abnormal situations createdXXOO
Sensitivity analysis of main factors in hydraulic/water quality analysisXX△ * (hydraulic analysis)O
Analysis of representative tasks (e.g., iterative water quality simulation) in hydraulic/water quality analysisXX△ (hydraulic analysis)O
Residual chlorine equalizationXXXO
* O, △, X denoted “included”, “partially included”, and “not included” separately.
Table 2. DDA and PDA hydraulic/water quality analysis results, considering additional demand for water.
Table 2. DDA and PDA hydraulic/water quality analysis results, considering additional demand for water.
EPANET2 DDAK-NRiskQ PDA
Junction IDPressure Head (m)Available Demand (CMD)Residual Chlorine (mg/L)Pressure Head (m)Available Demand (CMD)Residual Chlorine (mg/L)
GN-D28.3321,723.000.4531.0521,723.000.73
GD5.900.000.585.930.000.58
GY40.790.000.5440.880.000.56
MH32.190.000.5232.290.000.55
HN-J59.010.000.4959.140.000.53
PC-J52.640.000.4752.770.000.53
SY57.340.000.4557.470.000.51
JJ-J38.541000.000.2238.671000.000.39
PC-D12.295846.000.4412.535342.000.51
GC-J57.28294.000.4757.46294.000.51
SB-J59.030.000.4159.370.000.44
SP-J62.430.000.3762.890.000.48
MS64.160.000.3664.660.000.47
ID-J65.850.000.3466.370.000.43
YJ-J59.510.000.2960.120.000.43
SI-J68.000.000.2568.880.000.40
YB-P50.050.000.2251.010.000.40
YM-J27.380.000.1730.190.000.36
SJ-J48.000.000.1550.890.000.37
YB-D6.073426.000.149.262692.260.35
SP1-J60.980.000.2561.320.000.45
SBGI-J36.31116.000.0136.65116.000.06
SP-P52.750.000.3053.210.000.45
SP-D59.60835.000.2560.06835.000.43
ID-P13.490.000.2214.010.000.42
ID-D36.362760.000.2136.882760.000.38
HH38.980.000.2639.500.000.43
YP75.390.000.1076.010.000.39
YP-E83.70200.000.0584.32200.000.39
SI-P22.460.000.2224.480.000.40
SI-D10.965000.000.2013.944820.850.39
Min.5.90 0.015.93 0.06
Max.83.700.5884.320.73
Average44.700.3045.550.44
SD21.100.1520.710.11
Table 3. Hydraulic and water quality results under abnormal conditions.
Table 3. Hydraulic and water quality results under abnormal conditions.
EPANET2 DDAK-NRiskQ PDA
Junction IDPressure Head (m)Available Demand (CMD)Residual Chlorine (mg/L)Pressure Head (m)Available Demand (CMD)Residual Chlorine (mg/L)
GN-D28.3321,723.000.4531.0521,723.000.73
GD0.900.000.580.960.000.58
GY35.790.000.5435.990.000.55
MH27.190.000.5227.420.000.55
HN-J54.010.000.4954.290.000.53
PC-J47.640.000.4747.930.000.52
SY52.340.000.4552.630.000.50
JJ-J33.541000.000.2233.831000.000.37
PC-D7.295846.000.447.894240.620.49
GC-J52.28294.000.4752.64294.000.51
SB-J54.030.000.4154.680.000.44
SP-J57.430.000.3758.290.000.49
MS59.160.000.3660.090.000.45
ID-J60.850.000.3461.820.000.42
YJ-J54.510.000.2955.650.000.43
SI-J63.000.000.2564.590.000.45
YB-P45.050.000.2246.770.000.43
YM-J22.380.000.1726.940.000.39
SJ-J43.000.000.1547.690.000.41
YB-D1.073426.000.146.222206.120.40
SP1-J55.980.000.2556.630.000.44
SBGI-J31.31116.000.0131.97116.000.01
SP-P47.750.000.3048.620.000.48
SP-D54.60835.000.2555.47835.000.43
ID-P8.490.000.229.460.000.41
ID-D31.362760.000.2132.332760.000.41
HH33.980.000.2634.950.000.41
YP70.390.000.1071.530.000.41
YP-E78.70200.000.0579.84200.000.37
SI-P17.460.000.2221.740.000.40
SI-D5.965000.000.2012.524567.840.42
Min.0.90 0.010.96 0.01
Max.78.700.5879.840.73
Average39.860.3041.370.45
SD20.990.1520.340.11
Table 4. Water-quality results of DDA and PDA under normal conditions.
Table 4. Water-quality results of DDA and PDA under normal conditions.
No. Injection PointsItemDDAPDA
2Total injection amount (kg/d)
(PDA/DDA × 100, %)
67.967.61 (11.2)
Min. (mg/L)0.400.56
Max. (mg/L)3.991.58
Avg. (mg/L)2.090.81
Standard deviation1.060.18
3Total injection amount (kg/d)24.5310.25 (41.8)
Min. (mg/L)0.450.48
Max. (mg/L)3.901.35
Avg. (mg/L)1.200.82
Standard deviation0.810.23
4Total injection amount (kg/d)23.928.68 (36.3)
Min. (mg/L)0.410.42
Max. (mg/L)3.531.80
Avg. (mg/L)1.080.81
Standard deviation0.640.39
Table 5. Water-quality results of DDA and PDA under abnormal condition.
Table 5. Water-quality results of DDA and PDA under abnormal condition.
No. Injection PointsItemDDAPDA
2Total injection amount (kg/d)
(PDA/DDA × 100, %)
66.716.53 (9.8)
Min. (mg/L)0.410.43
Max. (mg/L)3.901.56
Avg. (mg/L)2.050.79
Standard deviation1.100.22
3Total injection amount (kg/d)19.194.89 (25.5)
Min. (mg/L)0.410.41
Max. (mg/L)2.232.17
Avg. (mg/L)0.890.73
Standard deviation0.400.37
4Total injection amount (kg/d)14.833.38 (22.8)
Min. (mg/L)0.440.40
Max. (mg/L)3.863.17
Avg. (mg/L)1.000.77
Standard deviation0.760.74
Table 6. Two injection-point results under normal conditions.
Table 6. Two injection-point results under normal conditions.
EPANET2 DDAK-NRiskQ PDA
Junction IDInjection Rate (mg/L)Injection Mass (kg/d)Residual Chlorine (mg/L)Injection Rate (mg/L)Injection Mass (kg/d)Residual Chlorine (mg/L)
GN-D00.4500.73
GD00.5800.58
GY00.5400.56
MH3.4767.643.990.417.420.96
HN-J03.8300.92
PC-J03.6500.92
SY03.4600.90
JJ-J01.6900.71
PC-D03.4400.89
GC-J03.6500.89
SB-J03.1900.89
SP-J02.8800.89
MS02.7500.84
ID-J02.6300.84
YJ-J02.2600.76
SI-J01.9200.74
YB-P01.7000.71
YM-J01.3200.70
SJ-J01.1800.67
YB-D01.0700.68
SP1-J01.9600.87
SBGI-J2.700.312.801.580.181.58
SP-P02.3500.87
SP-D01.9100.84
ID-P01.7200.69
ID-D01.5900.77
HH01.9800.74
YP00.7700.75
YP-E00.4000.69
SI-P01.6900.75
SI-D01.5700.70
Table 7. Two injection-point results under abnormal conditions.
Table 7. Two injection-point results under abnormal conditions.
EPANET2 DDAK-NRiskQ PDA
Junction IDInjection Rate (mg/L)Injection Mass (kg/d)Residual Chlorine
(mg/L)
Injection Rate (mg/L)Injection Mass (kg/d)Residual Chlorine
(mg/L)
GN-D00.4500.73
GD00.5800.58
GY00.5400.55
MH00.5200.55
HN-J3.4066.283.9000.53
PC-J03.7200.52
SY03.5300.50
JJ-J01.721.181.181.56
PC-D03.5000.49
GC-J03.7200.51
SB-J03.250.505.350.96
SP-J02.94 00.98
MS02.8000.93
ID-J 02.6800.90
YJ-J02.3000.91
SI-J01.9500.93
YB-P01.74 00.87
YM-J01.3500.85
SJ-J01.21 00.86
YB-D01.0900.81
SP1-J02.0000.93
SBGI-J3.690.433.8000.43
SP-P02.3900.97
SP-D01.9400.90
ID-P01.7500.83
ID-D01.6200.86
HH02.0200.87
YP00.7800.84
YP-E00.4100.74
SI-P01.7200.81
SI-D 01.6000.86

Share and Cite

MDPI and ACS Style

Lee, S.M.; Lee, H.M.; Yoo, D.G.; Kim, J.H. A Comparative Study on a Hydraulic and Water-Quality Analysis Method for Determining Rechlorination Injection Points for a Water-Supply Network. Water 2019, 11, 697. https://doi.org/10.3390/w11040697

AMA Style

Lee SM, Lee HM, Yoo DG, Kim JH. A Comparative Study on a Hydraulic and Water-Quality Analysis Method for Determining Rechlorination Injection Points for a Water-Supply Network. Water. 2019; 11(4):697. https://doi.org/10.3390/w11040697

Chicago/Turabian Style

Lee, Sang Myoung, Ho Min Lee, Do Guen Yoo, and Joong Hoon Kim. 2019. "A Comparative Study on a Hydraulic and Water-Quality Analysis Method for Determining Rechlorination Injection Points for a Water-Supply Network" Water 11, no. 4: 697. https://doi.org/10.3390/w11040697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop