Next Article in Journal
An Interface Extension for the New EPANET 2.X
Previous Article in Journal
Optimal Regulation of Variable Speed Pumps in Sewer Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Leakage Calibration in Water Distribution Networks with Pressure-Driven Analysis: A Real Case Study †

1
Departamento de Engenharia Civil, Instituto Politécnico de Coimbra, ISEC, Rua Pedro Nunes, 3030-199 Coimbra, Portugal
2
INESC Coimbra, DEEC, Universidade de Coimbra, Rua Sílvio Lima, Polo II, 3030-290 Coimbra, Portugal
3
Be Water-Águas de Ourém, Rua Dr. Carlos Vaz de Faria e Almeida, nº21 R/C, 2490-547 Ourém, Portugal
4
Departamento de Informática, Universidade da Beira Interior, Rua Marquês d’Ávila e Bolama, 6201-001 Covilhã, Portugal
5
Departamento de Engenharia Civil, Universidade de Coimbra, Rua Luís Reis Santos, Polo II, 3030-788 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Presented at the 4th EWaS International Conference: Valuing the Water, Carbon, Ecological Footprints of Human Activities, Online, 24–27 June 2020.
Published: 6 September 2020

Abstract

:
Leakage in water distribution networks (WDN) is still a major concern for water companies. In recent years, the scientific community has dedicated some effort to the leakage calibration issue to obtain accurate models. But leakage modelling implies the use of a pressure-driven approach as well as specific data to define the pressure/leakage relationship. This paper presents the calibration process of a real case study WDN model. The process started with pressure step tests, the model was built in WaterNetGen and the leakage calibration process was performed by a simulated annealing algorithm. As illustrated, after calibration the model was able to produce accurate results.

1. Introduction

Water loss in water distribution networks (WDN) is still a major concern for water companies. They are known to be the cause of many negative impacts, like economic losses, technical and social problems, environmental damages and difficulties in guaranteeing the water quality safety. Water losses include both apparent and real losses. While apparent losses are related to customer meter inaccuracies and unauthorized consumption, real losses correspond to the water lost due to leakage and overflows at storage tanks and leakage on transmission and distribution mains and in service connections up to the point of customer metering [1]. Unreported leaks (in mains and service connections) are a major cause of real losses because, even if the flow is relatively small, as they can last for a long time the total volume lost can be huge. Unreported leaks are difficult to locate, particularly in plastic pipes, large diameters and low-pressure conditions. In recent years, the scientific community has dedicated some effort to the leakage calibration issue in order to obtain accurate models for WDN, supported by flow and pressure measurements, and using modelling and optimization tools, considering both demand-driven analysis (DDA) [2,3] or pressure-driven analysis (PDA) [4,5]. However, the theoretically correct leakage modelling implies the use of PDA, which is far more complex than the more traditional DDA and needs specific data to define the pressure/leakage relationship.
This paper presents the calibration process of a real case study WDN model. Section 2 is dedicated to the theoretical background and presentation of the calibration methodology, Section 3 presents the case study and the results and finally Section 4 is dedicated to the conclusions.

2. Materials and Methods

2.1. Leakage Modelling

Leakage is assumed to be a function of the pressure and can be modelled with two components, bursts and background leakage [6]:
q k l e a k ( P k ) = { β k · L k · ( P k ) α k + C k · ( P k ) δ k P k > 0 0 P k 0
where qkleak is the total leakage along pipe k; Lk is the length of pipe k; αk and βk are parameters of the background leakage model; Ck and δk are parameters of the bursts leakage model; and Pk is the average pressure in pipe k computed as the mean of the pressure values of its end nodes.
The exponent αk depends on the network characteristics (pipe material and failure mode) and can take values between 0.5 and 2.5. The parameter βk must be set by calibration (initial values can be set around 10-7) [7].

2.2. Calibration Methodology

The first goal of this work is to calibrate leakage in WDN models, that is to estimate the pipe background leakage parameters (αk and βk) from Equation (1). The αk parameter is network specific and can be obtained by performing pressure step tests in the network being studied. The βk parameter is pipe specific and here is estimated by solving an optimization problem in which the objective function is the difference between measured and observed flows and pressures and the decision variables are the βk values for each pipe in the WDN model. This optimization problem is solved by a simulated annealing algorithm [8], with the help of a pressure-driven hydraulic simulation model, implemented in WaterNetGen [9,10]—an EPANET [11] extension, to assess the hydraulic constraints and estimate the pipe leakage.
The objective function (OF) can be written as:
Min .   OF = W P · Nodes | P m P e | + W F · Pipes | F m F e |
where WP is the weight for the pressure errors, Nodes is the set of nodes with pressure measurements, Pm and Pe are the pressures measured and estimated, respectively, WF is the weight for the flow errors, Pipes is the set of pipes with flow measurements, Fm and Fe are the flows measured and estimated, respectively.
The constraints that define the hydraulic behavior of the WDN, as well as the leakage estimates, are assessed by WaterNetGen.
The simulated annealing algorithm was implemented as described in Table 1. The procedure starts by assigning an initial value to the βk parameter for each pipe in the WDN model (in this case it was assigned βk = 0). A starting temperature is calculated and a certain number of candidate solutions are evaluated. Each candidate solution is obtained from the current solution by randomly choosing a pipe and changing its βk value (up or down) and its objective function is calculated with the results from the PDA simulation performed by WaterNetGen. The number of candidate solutions (Li) generated at each temperature (Ti) varies according to the percentage of solutions accepted at the last temperature (Pai−1). Each new solution is accepted or not, according to the Metropolis criterion. If it is accepted, this solution becomes the current solution and will be used to produce the next candidate solution. If not, the original current solution will be used. The algorithm ends if the stopping criteria are reached, that is, for two successive temperatures the number of solutions accepted remains lower than 5% and there is no improvement.
All this procedure was embedded in WaterNetGen and the user just has to build the WDN model (topology, tank properties, junction properties, pipe properties, base demand, demand patter, and other operational parameters or rules) and ask the software to solve the optimization problem (perform the calibration procedure).

3. Results and Discussion

The proposed methodology was applied to a real case study WDN model, namely the Loureira District Metered Area (DMA), Figure 1, which is a part of the Fatima WDN, managed by the water company Be Water—Águas de Ourém. This DMA is supplied by a pressure reduction valve (PRV) to manage the pressure, which was introduced in the model as a reservoir with a head pattern obtained from the pressure measurements obtained downstream the PRV. The pipes from this network are all made of PVC, with diameters between 63 and 315 mm, in a total of about 13 km. The WDN model was built from Geographical Information System (GIS) files containing information about the topology of the network and the characteristics of the network components.
The leakage calibration process started by performing pressure step tests to obtain the data for the pressure/leakage relationship (pressures were gradually reduced at the PRV while monitoring the effect on the minimum night flows). The data obtained (Figure 2) lead to a αk parameter equal to 1.8, which was introduced in the WDN model.
The flow measurements at the DMA Loureira entrance were used to estimate the leakage flow, by analyzing the minimum night flow. The difference between the flow measurements and the water demand is the water losses that should be assigned to the pipes to obtain a calibrated model, here estimated as about 5 m3/h (1.4 L/s). At the roundabout located in southeastern part of the network there is a garden irrigating during the night, with a flow of about 5 m3/h (1.4 L/s). As this consumption is known (value and location) it was decided to exclude it from the water demand and consider it as a leak to be located and, to a certain extent, to confirm the accuracy of the results.
The flow measurements, after excluding the water losses and the garden irrigation, resulted in the demand pattern used in the model (Figure 3).
The leakage calibration process was supported by the flow measurements obtained at the DMA Loureira entry point and the pressure measurements obtained in three points of the network (P3, P4 and P5), Figure 4.
All these data, together with the GIS files data, were used to build the WDN model (Figure 5) in WaterNetGen, an EPANET extension including a PDA engine.
The topology of this WDN, with a radial shape, does not ease the leakage calibration process since there are no long pipes with significant head losses making them quite sensitive to flow changes. On the other hand, this WDN has another particular characteristic that can make things even more difficult for the calibration. This network is in the vicinity of the Fatima sanctuary, a world-famous religious destination. In a few days of the year (particularly in the 13th of May) this place becomes full of pilgrims and the water demand increases a lot. Consequently, this WDN is designed to face those huge water demands but in most of the year it is oversized, presenting low sensitivity to pressure; even so, it is worth a try.
The leakage calibration process (assignment of the leakage parameters, βk, to the pipes) was performed by the simulated annealing algorithm used to solve the optimization problem aimed at minimizing the differences between the field measurements and the PDA simulation results. The process took about one hour to produce the results and the calibration reports (for flow and pressure) presented in Figure 6 seem to confirm the accuracy of the proposed methodology, showing a very good fitting between the observed (measurements) and the computed (results of the PDA simulation) values.
Although the calibration results seem to be quite accurate, many different combinations of flows in the model can fit the measurements and be quite different from the flow distribution occurring in the real network. The only way to check the accuracy of the calibration is to confirm the existence of leaks in the pipes pointed out by the methodology in the WDN model.
The leakage calibration process identified a few potential leaky pipes in the surroundings of the roundabout (probably related to the garden irrigation used as a known leak) and a concentration of potential leaky pipes at the western zone of the network (1.2 L/s), Figure 7. Before the closure of the paper submission period, it was still possible to execute a step-test (closing isolation valves during the night) and this confirmed the existence of a considerable night flow at the western zone of the network, probably due to the existence of leaks. Further work will include an acoustic survey in that part of the network to confirm/locate those leaks. There were also other pipes in the WDN that were pointed out as probable leaky pipes but without any considerable significance.
The first goal of this work was to estimate the pipe background leakage parameters (βk) in order to obtain a calibrated WDN model. The results obtained apparently lead to the conclusion that the goal was accomplished because the calibration reports show that the simulation results fit well with the field measurements. However, the final conclusion can only be drawn after carrying out field acoustic inspection works to check if there are really leaks in the pipes identified as probable leaky pipes.

4. Conclusions

This paper presents the calibration process of a real case study WDN model. The process started by performing pressure step tests to obtain the exponents of the pressure/leakage relationship (pressures were gradually reduced while monitoring the effect on the minimum night flows). The WDN model was built in WaterNetGen, an EPANET extension including a PDA engine, which is used to assess the hydraulic constraints and estimate the pipe leakage. The leakage calibration process (assignment of the pipe background leakage parameters, βk) was performed by a simulated annealing algorithm used to solve an optimization problem aimed at minimizing the differences between field measurements (flows and pressures) and simulation results. The calibration reports show a good fit between the simulation results and the field measurements, leading to the conclusion that the proposed leakage calibration methodology is accurate. However, the final conclusion can only be drawn after carrying out field acoustic inspection works to check if there are really leaks in the pipes identified as probable leaky pipes. By identifying the most probable leaky pipes, the results from this leakage calibration methodology will certainly be of great help in the context of active leakage control, optimizing the resources allocated to the leak location process.

Author Contributions

All authors have read and agree to the published version of the manuscript. The authors contributed equally to this work.

Funding

This research was partially funded by Portuguese Foundation for Science and Technology under project UIDB/00308/2020.

Acknowledgments

The authors thank to Be Water-Águas de Ourém for providing the data and the permission for its use.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lambert, A.O.; Brown, T.G.; Takizawa, M.; Weimer, D. A review of performance indicators for real losses from water supply systems. J. Water Supply Res. Tech. Aqua 1999, 48, 227–237. [Google Scholar] [CrossRef]
  2. Ribeiro, L.; Sousa, J.; Sá Marques, A.; Simões, N.E. Locating leaks with TrustRank algorithm support. Water 2015, 7, 1378–1401. [Google Scholar] [CrossRef]
  3. Wu, Z.Y. Water loss detection via genetic algorithm optimization-based model calibration. In Proceedings of the 8th Annual Water Distribution System Analysis Symposium (WDSA), Cincinnati, OH, USA, 27–30 August 2006. [Google Scholar] [CrossRef]
  4. Giustolisi, O.; Savic, D.; Kapelan, Z. Pressure-driven demand and leakage simulation for water distribution networks. J. Hydraul. Eng. 2008, 134, 626–635. [Google Scholar] [CrossRef]
  5. Wu, Z.Y.; Sage, P.; Turtle, D. Pressure-dependent leak detection model and its application to a district water system. J. Water Res. Plan. Manag. 2010, 136, 116–128. [Google Scholar] [CrossRef]
  6. Germanopoulos, G. A technical note on the inclusion of pressure dependent demand and leakage terms in water supply network models. Civ. Eng. Syst. 1985, 2, 171–179. [Google Scholar] [CrossRef]
  7. Lambert, A. What do we know about pressure: Leakage relationships in distribution systems? In Proceedings of the IWA System Approach to Leakage Control and Water Distribution Systems Management, Brno, Czech Republic, 16–18 May 2001. [Google Scholar]
  8. Sousa, J.; Muranho, J.; Sá Marques, A.; Gomes, R. Optimal Management of Water Distribution Networks with Simulated Annealing: The C-Town Problem. J. Water Res. Plan. Manag. 2016, 142, C4015010. [Google Scholar] [CrossRef]
  9. Muranho, J.; Ferreira, A.; Sousa, J.; Gomes, A.; Sá Marques, A. WaterNetGen: An EPANET extension for automatic water distribution network models generation and pipe sizing. Water Supply 2012, 12, 117–123. [Google Scholar] [CrossRef]
  10. Muranho, J.; Ferreira, A.; Sousa, J.; Gomes, A.; Sá Marques, A. Pressure dependent demand and leakage modeling with an EPANET extension—WaterNetGen. Procedia Eng. 2014, 89, 632–639. [Google Scholar] [CrossRef]
  11. Rossman, L. EPANET 2 Users’s Manual; Risk Reduction Eng. Lab., Office of Research and Development, U.S. Environmental Protection Agency: Cincinnati, OH, USA, 2000. [Google Scholar]
Figure 1. District Metered Area (DMA) Loureira water distribution networks (WDN) with the pressure reduction valve (PRV) and the pressure monitoring points (P3, P4 and P5).
Figure 1. District Metered Area (DMA) Loureira water distribution networks (WDN) with the pressure reduction valve (PRV) and the pressure monitoring points (P3, P4 and P5).
Environsciproc 02 00059 g001
Figure 2. Pressure/leakage relationship for DMA Loureira.
Figure 2. Pressure/leakage relationship for DMA Loureira.
Environsciproc 02 00059 g002
Figure 3. Demand pattern for DMA Loureira.
Figure 3. Demand pattern for DMA Loureira.
Environsciproc 02 00059 g003
Figure 4. Flow and pressure measurements for DMA Loureira.
Figure 4. Flow and pressure measurements for DMA Loureira.
Environsciproc 02 00059 g004
Figure 5. WDN model for DMA Loureira.
Figure 5. WDN model for DMA Loureira.
Environsciproc 02 00059 g005
Figure 6. Calibration reports (flow and pressure) for DMA Loureira.
Figure 6. Calibration reports (flow and pressure) for DMA Loureira.
Environsciproc 02 00059 g006
Figure 7. Probable leaky pipes in DMA Loureira.
Figure 7. Probable leaky pipes in DMA Loureira.
Environsciproc 02 00059 g007
Table 1. Simulated annealing algorithm implementation.
Table 1. Simulated annealing algorithm implementation.
Part of the AlgorithmImplementation
Initial solutionAssign a βk value to each pipe (Xo)
Initial temperature T i = 0.1 O F ( X 0 ) L o g ( 0.5 )
Control parametersIf Pa > 70% Ti+1 = 0.5 Ti and IAi+1 = 10
If Pa > 50% Ti+1 = 0.6 Ti and IAi+1 = 20
If Pa > 35% Ti+1 = 0.7 Ti and IAi+1 = 30
If Pa > 20% Ti+1 = 0.8 Ti and IAi+1 = 40
If Pa ≤ 20% Ti+1 = 0.9 Ti and IAi+1 = 50
Number of evaluations (Li) at
each temperature (Ti)
IAi × Number of pipes in the network
Stopping criteriaPa < 5% and 2 temperatures without improvement
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sousa, J.; Martinho, N.; Muranho, J.; Marques, A.S. Leakage Calibration in Water Distribution Networks with Pressure-Driven Analysis: A Real Case Study. Environ. Sci. Proc. 2020, 2, 59. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002059

AMA Style

Sousa J, Martinho N, Muranho J, Marques AS. Leakage Calibration in Water Distribution Networks with Pressure-Driven Analysis: A Real Case Study. Environmental Sciences Proceedings. 2020; 2(1):59. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002059

Chicago/Turabian Style

Sousa, Joaquim, Nuno Martinho, João Muranho, and Alfeu Sá Marques. 2020. "Leakage Calibration in Water Distribution Networks with Pressure-Driven Analysis: A Real Case Study" Environmental Sciences Proceedings 2, no. 1: 59. https://0-doi-org.brum.beds.ac.uk/10.3390/environsciproc2020002059

Article Metrics

Back to TopTop