Next Article in Journal
A Review on Rainfall Data Resolution and Its Role in the Hydrological Practice
Previous Article in Journal
The Influence Mechanism of Freeze-Thaw on Soil Erosion: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Estimation of Energy Balance Components in Water Networks for Top-Down Approach

Department of Water Resources Engineering, Faculty of Engineering, Kasetsart University, Bangkok 10900, Thailand
*
Author to whom correspondence should be addressed.
Water 2021, 13(8), 1011; https://doi.org/10.3390/w13081011
Submission received: 19 February 2021 / Revised: 22 March 2021 / Accepted: 6 April 2021 / Published: 7 April 2021
(This article belongs to the Section Urban Water Management)

Abstract

:
The energy balance calculation for pressurized water networks is an important step in assessing the energy efficiency of water distribution systems. However, the calculation generally requires mathematical modelling of the water networks to estimate three important energy components: outgoing energy through water loss ( E l ), friction energy loss ( E f ) and energy associated with water loss ( E W L ). Based on a theoretical energy balance analysis of simplified pipe networks, a simple method is proposed to estimate E l , E f and E W L with minimum data requirements: input energy, water loss (WL) and head loss between the source and the minimum energy point (ΔH). By inclusion of the head loss in water networks into the estimation, the percentages of E l and E W L are lower and higher, respectively, than using only the percentage of WL. The percentage of E f can be a function of the percentage of ΔH. By demonstrating our analysis with the simulation results from the mathematical models of 20 real water networks, the proposed method can be used to effectively estimate E l , E f and E W L as a top-down energy balance approach.

1. Introduction

The demand for energy is intensifying in the water sector. The International Energy Agency (IEA) [1] estimated that 4% of total global electricity was consumed by the water sector in 2014 and it may rise by 80% by 2040. In Portugal, urban water and wastewater systems use 3%–4% of the total national electricity [2]. In Brazil, over 1.9% of total electricity energy is consumed by water supply systems [3]. Pelli and Hitz [4] estimated that supplying drinking water and industrial water worldwide can consume as much as 2%–10% of a country’s total electricity usage. The World Bank [5] reported that electricity costs are 5%–30% of the total operating cost of water and wastewater utilities worldwide. In Latin America and the Caribbean, energy is generally 30%–40% of the operational costs of water supply services [6]. Nevertheless, IEA [1] suggested that there is a great opportunity to save energy in the water sector by 15% by 2040 if economically available energy efficiency and energy recovery potentials are exploited. Achieving this will require accurate and comprehensive energy assessment in the water sector.
In the past, energy assessment in water supply systems was focused on pump inefficiency. Later, the interest extended to assess the energy losses due to friction and leakage in water networks [4,7,8]. Bylka and Mroz [9] reviewed energy assessment methods for water supply systems and categorized them into two groups. While the first group is based solely on data collected in water utilities, the second one is based on both the data and modelling of physical processes. Thus, more detailed results from the second group can support a strategy for increasing the energy efficiency of each process. Since mathematical modelling requires collecting a large amount of data and is time consuming to analyze, it is usually applied on a relatively small network. Cabrera et al. [10] presented one of the first energy audit methods for water networks based on the modelling of physical processes. Over the past decade, more complex and comprehensive methods have been developed for energy balance and assessment based on physical processes to evaluate energy transformation and efficiency in water supply systems [11,12,13,14,15,16,17,18,19,20,21,22]. In addition, there are successful energy assessments of real networks based on physical processes in many countries [19,20,21,22,23,24,25,26,27].
Figure 1 shows two simplified versions of energy balances in water networks considering only leakage (as water loss) and pipe friction proposed by Cabrera et al. [10] and Mamade et al. [19,20,21]. The energy balance in Cabrera et al. [10] (Figure 1a) shows the input energy ( E i n ) divided into three components: the energy delivered to users ( E u ), the outgoing energy through water loss ( E l ) and the friction energy loss ( E f ). To be consistent with the International Water Association (IWA) water balance [28], on the other hand, Mamade et al. [19,20,21] divided E i n into two components: the energy associated with authorized consumption ( E A C ) and the energy associated with water loss ( E W L ), as shown in Figure 1b. Theoretically, these components should be calculated through mathematical models, which are based on physical processes. For the balance in Cabrera et al. [10], a calibrated mathematical model is sufficient to evaluate E u , E l and E f . However, for the balance in Mamade et al. [19,20,21], an additional model without water loss is needed to estimate the friction energy loss for a water loss-free network ( E f o ). Mamade et al. [21] proposed the top-down and bottom-up approaches for their energy balance. In the top-down approach, the ratio between E W L and E i n is assumed to be the ratio between water loss (WL) and system input volume (SIV). To prove the assumption, the comparison between E W L / E i n and WL/SIV was carried out through simulation using the mathematical models of real networks. The linear trendline shows that the value E W L / E i n is slightly higher than WL/SIV. In this study, the evaluation of E W L / E i n is reanalyzed by using both theoretical analysis and mathematical models and the estimations of E l / E i n and E f / E i n are also investigated. In Section 2, the theoretical analysis of energy balance describes the behaviors of each energy balance component and establishes the basic relationships between each variable. In Section 3, the mathematical models of 20 real district metering area (DMA) networks are introduced, and the relationships between the model results and basic parameters are investigated. Section 4 combines the results from the theory and the models, and the coefficient in the theory is calibrated using the model simulation results. The accuracy of the estimation of E l , E f and E W L is evaluated in order to apply as the top-down energy balance approach. Finally, the conclusion is provided in Section 5.

2. Theoretical Analysis of Energy Balance

In this section, the energy balance principle is introduced and analyzed theoretically on two simple pressurized water networks. The energy balance components are calculated and then normalized to the non-dimensional version. They are used for the top-down energy balance estimation.

2.1. Single Pipe Network

As the first case, Figure 2 shows a single pipe network consisting of a source and a demand node at the pipe end. The input energy head is defined as H, and the system inflow is Q. Using the water balance concept, Q can be divided into
Q = Q u + Q l
where Q u is the flow to supply authorized consumption and Q l is the flow due to water loss.
Then, we introduce the ratio of water loss (p) as
p = Q l Q
At the demand node at the pipe end, Q is out of the network and thus, the friction head loss along the pipe ( Δ H ) can be expressed as
Δ H = K Q n
where K is the loss coefficient and n is the flow exponent.
According to the energy balance concepts in Figure 1, the energy balance components can be calculated as
  • Input energy ( E i n )
    E i n = γ Q H
  • Outgoing energy through water loss ( E l )
    E l = γ Q l ( H K Q n )
  • Friction energy loss ( E f )
    E f = γ Q ( K Q n )
  • Friction energy loss for a water loss-free network ( E f o )
    E f o = γ Q u ( K Q u n )
  • Energy associated with water loss ( E W L )
    E W L = E l + E f E f o
    where γ is specific weight of water.
The nondimensionalization is introduced here. According to Mamade et al. [21], E W L divided by E i n can be approximated by p. Therefore, the energy balance components in Equations (4)–(8) are normalized by E i n as
( E i n , E l , E f , E f o , E W L ) = ( E i n , E l , E f , E f o , E W L ) E i n
and Δ H is normalized by H as
Δ H = K Q n H = Δ H H
where the variables with superscript ‘ are the normalized versions of the variables.
Thus, the normalized energy balance components for the first case in Figure 2 can be written as
  • Normalized input energy ( E i n )
    E i n = 1
  • Normalized outgoing energy through water loss ( E l )
    E l = ( Q l Q ) ( H K Q n H ) = p p Δ H
  • Normalized friction energy loss ( E f )
    E f = Δ H
  • Normalized friction energy loss for a water loss-free network ( E f o )
    E f o = ( Q u Q ) ( K Q u n H ) = ( 1 p ) n + 1 Δ H
  • Normalized energy associated with water loss ( E W L )
    E W L = p + p n Δ H
where
p n = ( 1 p ) [ 1 ( 1 p ) n ]
The energy components relating to water loss, E l and E W L in Equations (11) and (14), have a linear relationship with p only if Δ H is null. Thus, the hypothesis of Mamade et al. [21] is valid under the no energy loss condition. Both E l and E W L depend on Δ H in Equation (9b). The value of Δ H should be between 0 and 1. If the head loss ( Δ H ) is sufficiently smaller than the input head ( H ), Δ H will be close to 0. As a result, E l is slightly smaller than p, while E W L is slightly larger than p. On the other hand, if Δ H is comparable to H , Δ H will approach unity, and the discrepancies between E l , E W L and p grow with an increase in Δ H . Figure 3 shows pn as a function of p from Equation (15) for n = 1, 1.852 and 2. If n = 1, it implies a laminar flow, and the curve is a negative parabola with the maximum pn of 0.25 at p = 0.5 and the minimum pn of 0 at p = 0 and 1. For the case of complete turbulence, rough pipes (n = 2), the curve is skewed to the right with the maximum pn of around 0.385 at p ≈ 0.423. For the Hazen–Williams formula (n = 1.852), the curve differs from the n = 2 curve slightly. Thus, E W L approaches p when p is very low or extremely high, and the maximum deviation of E W L from p occurs when the percentage of water loss is around 40%–50%. In this study, n is assumed to be 2 for simplicity.

2.2. Branched Pipe Network with Uniformly Distributed Demand Nodes

The effects of branching and demand distribution on energy balance are investigated theoretically in this section. Here, the second network is introduced. It consists of branched pipes, and each branched pipe has uniformly distributed demand nodes as shown in Figure 4. A parameter m denotes the number of branches, and j denotes the number of demand nodes in each branch. Each branch has the same hydraulic properties, where the demand nodes have the same outflow of Q/m/j and they distribute uniformly along the branch. In addition, the constant friction slope along each branch is assumed, and the head loss equals K(Q/m)n at the end of each branch. These assumptions are used to simplify the problem and to make it possible to be solved theoretically. They will be discussed in the later section using the results from the real network models.
For the network in Figure 4, the energy balance components can be calculated as follows.
  • Input energy ( E i n )
    E i n = γ Q H
  • Outgoing energy through water loss ( E l )
    E l = m i = 1 j γ Q l i m [ H i j K ( Q m ) n ] = γ Q l [ H ( 1 2 m n ) ( 1 + 1 j ) K Q n ]
  • Friction energy loss ( E f )
    E f = m i = 1 j γ Q i m K ( Q m ) n j = ( 1 + 1 j ) γ Q ( K Q n ) 2 m n
  • Friction energy loss for a water loss-free network ( E f o )
    E f o = ( 1 + 1 j ) γ Q u ( K Q u n ) 2 m n
  • Energy associated with water loss ( E W L )
    E W L = E l + E f E f o
Using the nondimensionalization in Equation (9a,b) the normalized energy balance components can be written as
  • Normalized input energy ( E i n )
    E i n = 1
  • Normalized outgoing energy through water loss ( E l )
    E l = p C m j p Δ H
  • Normalized friction energy loss ( E f )
    E f = C m j Δ H
  • Normalized friction energy loss for a water loss-free network ( E f o )
    E f o = C m j ( 1 p ) n + 1 Δ H
  • Normalized energy associated with water loss ( E W L )
    E W L = p + C m j p n Δ H
    where
    C m j = ( 1 2 m n ) ( 1 + 1 j )
The normalized energy balance components in Equations (22)–(25) of the second case are the generalized versions of the results from the first case in Equations (11)–(14), respectively, considering the effects of branching and demand distribution on the energy balance. The parameter C m j in Equation (26) is a multiplier in the Δ H term and has a value between 0 and 1. For the network with one branch (m = 1) and one demand node (j = 1), C m j equals to 1 and the components in Equations (22)–(25) are simplified into Equations (11)–(14), respectively. An increase in branches (m > 1) and demand nodes (j > 1) produces a smaller value of C m j and subsequently, an increase in E l but a decrease in E f , E f o and E W L . Thus, C m j is the coefficient representing different water distribution patterns. Since C m j in Equation (26) is derived theoretically under many assumptions, we are going to find appropriate values of C m j for each energy balance component using the results of our network models later.

2.3. Utilization of Theory to Real Networks

The estimations of three normalized energy balance components, E l , E f and E W L , for real networks are discussed here. Up to this point, our theoretical analysis considers the energy balance of a single time period. To apply the theory to a real application, the extended period analysis is performed to evaluate the time-averaged values of E l , E f and E W L . However, the time-averaged Δ H terms may be too complicated for a top-down approach. Thus, for simplicity, we proposed the following theoretical estimations for real networks.
E l , t h e o = p ¯ C m j p Δ H ¯     p C m j p Δ H *
E f , t h e o = C m j Δ H ¯     C m j Δ H *
E W L , t h e o = p ¯ + C m j p n Δ H ¯     p + C m j p n Δ H *
where the subscript t h e o denotes the theoretical estimation, the bar ¯ expresses time average and Δ H * is the renormalized head loss, written as
Δ H * = H ¯ max H ¯ min H ¯ max
where H ¯ max is the maximum time-averaged energy head and H ¯ min is the minimum time-averaged energy head. The source with the maximum head can be chosen as the H ¯ max point, and the location where the maximum head loss occurs is the H ¯ min point. Thus, the H ¯ min location depends on network topology as well as demand distribution. The highest possibility is the furthest dead end from the source.

3. Application to Real Water Networks

3.1. Characteristics of Water Networks

In this section, 20 real water networks (Table 1) are used to test our theoretical analysis of energy balance. They are real district metering areas (DMAs) in the service area of the Samut Prakan branch office of the Metropolitan Waterworks Authority (MWA), Thailand. Samut Prakan is a province located at the mouth of the Chao Phraya River. Thus, the elevations refer to mean sea level (MSL). While the first 10 DMAs have 1 inlet (district meter), the remaining 10 DMAs have 2 inlets. The number of customers ranges from 739 to 11,545 connections with an average of 3508 connections. The average length and distribution pipe diameter are 40.8 km and 159 mm, respectively. While the average friction slope (Sf) from our network models ranges from 0.07 m/km to 0.70 m/km with an average of 0.23 m/km, Sf from the study of Mamade et al. [21] ranges from 0.1 m/km to 1.2 m/km with an average of 0.19 m/km. Based on the average value, the networks in our study have larger values of Sf. Mamade et al. [21] explained that their networks are overdesigned; thus, the impact on head loss by adding water loss is small. This implies that the estimation of the normalized energy associated with water loss ( E W L ) is not sensitive to head loss in their study. However, using our theory, E W L , t h e o in Equation (29) is a function of the normalized water loss (p) and head loss ( Δ H * ), not Sf. Our case study covers values of p between 2.8% and 54.9%, and values of Δ H * between 7.6% and 65.3%. Thus, the impact of head loss on the estimation of energy balance components can be investigated in our study.
Our pipe networks are modeled using the EPANET version 2.00.12.01 software (United States Environmental Protection Agency (US EPA), Columbus, OH, USA) [29]. Each customer is connected to distribution pipes by creating imaginary valves with no friction loss (not considering the friction loss on the service connection). Water loss is assumed to be pressure-dependent and simulated using the emitter function in the EPANET software. The emitter coefficients are distributed to all nodes of the networks. These extended period simulation (EPS) models are calibrated and provided by MWA and simulated on the average of the hourly base using the flow and pressure measurements and monthly water sales in March 2019. Thus, the models provide high levels of detail and accuracy. As examples, Figure 5 shows four DMAs (ID1, ID2, ID11 and ID13) used in this study, where the water distribution patterns of our networks vary from branching to semi-gridiron.

3.2. Basic Relationship for Energy Balance Components

The energy balance components extracted from the simulation results using EPANET (United States Environmental Protection Agency (US EPA), Columbus, OH, USA) models of 20 water networks are compared with the basic network parameters (p and Δ H * ) in Figure 6. The normalized outgoing energy through water loss by the models ( E l , m o d ) shows a good relationship with p (Figure 6a). The values of E l , m o d are slightly to moderately smaller than p. This tendency corresponds with the estimation of E l , t h e o in Equation (27) by our theoretical analysis, showing that E l , t h e o is smaller than p due to the second term that relates to Δ H * . Figure 6b shows a good relationship between the normalized friction energy loss by the models ( E f , m o d ) and Δ H * . While E f , m o d for 1 inlet is comparable to Δ H * , E f , m o d for 2 inlets seems to be smaller than Δ H * . According to Equation (28), E f , t h e o depends on only Δ H * with C m j as the slope. We hypothesize that the number of inlets affects the value of C m j . Thus, the calibration process of C m j in the next section will divide into two groups as 1 inlet and 2 inlets. In Figure 6c, the normalized energy loss associated with water loss by the models ( E W L , m o d ) shows a good relationship with p, similar to the case of E l , m o d in Figure 6a. However, the values of E W L , m o d are slightly to moderately larger than p. The result agrees with the theoretical E W L , t h e o in Equation (29) that E W L , t h e o is higher than p because of the additional friction energy due to an increased flow by water loss. The effect increases as Δ H * increases.

4. Estimation of Energy Balance Components

In this section, the theoretical estimations of the energy balance components ( E l , t h e o , E f , t h e o and E W L , t h e o in Equations (27)–(29)) are compared with the model results. Using the method of least squares, the values of the coefficient C m j are evaluated separately for each energy balance component and the number of inlets. Figure 7 shows the comparisons between the three energy balance components using the EPANET models and the theoretical estimations. Table 2 shows the values of the calibrated C m j used for Figure 7. In the table, “Before” means the basic relationship case that E l and E W L are estimated to be p and E f equals to Δ H * . This case corresponds to the results in Figure 6. “After” refers to the theoretical estimations after considering the calibrated C m j in Figure 7. The correlation (r) and the root mean square error (RMSE) for the “Before” and “After” cases are evaluated. The agreements improve substantially as the terms of Δ H * with the calibrated C m j are considered. E l and E W L can be estimated accurately as their RMSEs are around 2%. It is possible to estimate E f roughly as its RMSE is around 7%. An increase in the number of inlets from 1 to 2 causes a decrease in C m j , implying that the effects of branching and demand distribution become stronger. Although a number of assumptions have been used, the statistical evaluation in Table 2 shows a good performance of our theory over a wide range of p and Δ H * . Thus, our theoretical method can be used as an effective tool to estimate the energy balance components as a top-down energy balance approach.
Figure 8 shows the boxplots of the distribution of water demand in the network subarea corresponding to the head loss distribution. We cluster demand nodes of the network into 10 subareas according to the ratio of the head loss in each subarea to the total head loss ( Δ H i / Δ H ) as 0.1, 0.2, …, 1. For example, the subarea of Δ H i / Δ H = 0.1 means the subarea covering the demand nodes having head loss between 0%–10% of the total head loss, and the subarea of Δ H i / Δ H = 0.2 means the one with head loss between 10%–20% of the total head loss. Then, the water demands including water loss of all nodes in each subarea are summarized and normalized by the total water demand ( Q i / Q i n ). In our theory, the uniform distribution of water demand is assumed; thus, Q i / Q i n must be constant. Under this assumption, Q i / Q i n should equal to 0.1 for every subarea in this example. However, for the networks with 1 inlet in Figure 8a, more Q i / Q i n is distributed in the subareas with higher Δ H i / Δ H , implying that larger water demand is located further from the source in the energy sense. This nonuniform distribution of the demand amplifies the effect of head loss on the energy balance components. For the networks with 2 inlets, Q i / Q i n is distributed in the subareas with smaller Δ H i / Δ H to a greater extent than for the case of 1 inlet. Thus, the effect of head loss on the energy balance components is less, as the values of C m j for 2 inlets are smaller than the ones for 1 inlet in Table 2.

5. Conclusions

To perform the energy audit based on physical processes without mathematical modelling (top-down approach), the normalized energy loss associated with water loss by the input energy ( E W L ) was simply assumed to equal the water loss ratio (p) as in previous studies [19,20,21]. Using the theoretical analysis of energy balance, the renormalized head loss by the input head ( Δ H * ) was another important parameter to estimate E W L accurately. For a practical application, Δ H * can be calculated by measuring maximum and minimum energy heads in a network. In addition, the normalized outgoing energy through water loss ( E l ) and the normalized friction energy loss ( E f ) were first derived theoretically. These energy components are fundamental and crucial for an accurate energy assessment. The effects of demand distribution in the networks with 1 and 2 inlets on the estimation of energy balance components showed that the greater number of inlets caused less head loss effects on the energy balance components. These effects can be evaluated effectively through the coefficient C m j in our theory.

Author Contributions

Conceptualization, S.L. and A.P.; methodology, A.P.; software, S.L., S.K., N.C., J.C. and A.P.; validation, S.L., S.K. and A.P.; formal analysis and investigation, S.K., N.C. and A.P.; writing—original draft preparation, A.P.; writing—review and editing, S.L., J.C. and A.P.; visualization, N.C. and A.P.; supervision, S.L., J.C. and A.P.; funding acquisition, A.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the MWA Waterworks Institute of Thailand (MWAIT), grant number 109/2561. S.K. was supported by the Ph.D. scholarship by the Faculty of Engineering, Kasetsart University, grant number 61/01/WE/D.ENG.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data was provided from Metropolitan Waterworks Authority, Thailand (MWA). Direct requests for these materials may be made to the provider, as indicated in the acknowledgments.

Acknowledgments

The authors would like to thank three anonymous reviewers for their constructive reviews that helped to improve the manuscript greatly. In addition, the authors gratefully acknowledge Metropolitan Waterworks Authority (MWA) for providing data sets using in this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following symbols are used in this paper:
DMAdistrict metering area
IWAInternational Water Association
MWAMetropolitan Waterworks Authority, Thailand
C m j parameter in Equation (26)
Dpipe diameter
E A C energy associated with authorized consumption
E f friction energy loss
E f normalized friction energy loss
E f , m o d normalized friction energy loss evaluated by mathematical model
E f , t h e o normalized friction energy loss estimated by theory
E f o friction energy loss for a water loss-free network
E f o normalized friction energy loss for a water loss-free network
E i n input energy
E i n normalized input energy
E l outgoing energy through water loss
E l normalized outgoing energy through water loss
E l , m o d normalized outgoing energy through water loss by mathematical model
E l , t h e o normalized outgoing energy through water loss by theory
E W L energy associated with water loss
E W L normalized energy associated with water loss
E W L , m o d normalized energy associated with water loss by mathematical model
E W L , t h e o normalized energy associated with water loss by theory
Hinput energy head
jnumber of demand nodes in each branch
Kloss coefficient
mnumber of branches
nflow exponent in head loss formula
p ratio of water loss
p n coefficient as a function of p
Q inflow
Q i flow in subarea i
Q l flow due to water loss
Q u flow to supply authorized consumption
Sffriction slope
SIVsystem input volume
WLwater loss
γ specific gravity
Δ H head loss between the source and the minimum energy point
Δ H normalized head loss between the source and the minimum energy point
Δ H * renormalized head loss in Equation (30)
Δ H i head loss in subarea i

References

  1. International Energy Agency (IEA). Water Energy Nexus-Excerpt from the World Energy Outlook 2016; World Energy Outlook; IEA: Paris, France, 2016; p. 60. [Google Scholar]
  2. Loureiro, D.; Silva, C.; Cardoso, M.A.; Mamade, A.; Alegre, H.; Rosa, M.J. The Development of a Framework for Assessing the Energy Efficiency in Urban Water Systems and Its Demonstration in the Portuguese Water Sector. Water 2020, 12, 134. [Google Scholar] [CrossRef] [Green Version]
  3. Vilanova, M.R.N.; Balestieri, J.A.P. Exploring the water-energy nexus in Brazil: The electricity use for water supply. Energy 2015, 85, 415–432. [Google Scholar] [CrossRef]
  4. Pelli, T.; Hitz, H. Energy indicators and savings in water supply. J. Am. Water Work. Assoc. 2000, 92, 55–62. [Google Scholar] [CrossRef]
  5. World Bank. A Primer on Energy Efficiency for Municipal Water and Wastewater Utilities; Technical Report; 001.12. 2012b; ESMAP: Washington, DC, USA, 2012. [Google Scholar]
  6. WWAP. The United Nations World Water Development Report 2014: Water and Energy; UNESCO: Paris, France, 2014; Volume 1. [Google Scholar]
  7. Colombo, A.F.; Karney, B.W. Energy and Costs of Leaky Pipes: Toward Comprehensive Picture. J. Water Resour. Plan. Manag. 2002, 128, 441–450. [Google Scholar] [CrossRef] [Green Version]
  8. Colombo, A.F.; Karney, B.W. Impacts of Leaks on Energy Consumption in Pumped Systems with Storage. J. Water Resour. Plan. Manag. 2005, 131, 146–155. [Google Scholar] [CrossRef]
  9. Bylka, J.; Mroz, T. A Review of Energy Assessment Methodology for Water Supply Systems. Energies 2019, 12, 4599. [Google Scholar] [CrossRef] [Green Version]
  10. Cabrera, E.; Pardo, M.A.; Cobacho, R. Energy Audit of Water Networks. J. Water Resour. Plan. Manag. 2010, 136, 669–677. [Google Scholar] [CrossRef]
  11. Cabrera, E.; Cobacho, R.; Soriano, J. Towards an Energy Labelling of Pressurized Water Networks. Procedia Eng. 2014, 70, 209–217. [Google Scholar] [CrossRef] [Green Version]
  12. Cabrera, E.; Gomez, E.; Cabrera, E., Jr.; Soriano, J. Calculating the Economic Level of Friction in Pressurized Water Systems. Water 2018, 10, 763. [Google Scholar] [CrossRef] [Green Version]
  13. Cabrera, E.; Gómez, E.; Soriano, J.; Del Teso, R. Eco-Layouts in Water Distribution Systems. J. Water Resour. Plan. Manag. 2019, 145, 04018088. [Google Scholar] [CrossRef]
  14. Cabrera, E.; Gomez, E.; Soriano, J.; Espert, V. Energy Assessment of Pressurized Water Systems. J. Water Resour. Plan. Manag. 2015, 141, 04014095. [Google Scholar] [CrossRef]
  15. Gomez, E.; Cabrera, E.; Balaguer, M.; Soriano, J. Direct and Indirect Water Supply: An Energy Assessment. Procedia Eng. 2015, 119, 1088–1097. [Google Scholar] [CrossRef] [Green Version]
  16. Hashemi, S.; Filion, Y.R.; Speight, V.L. Pipe-level Energy Metrics for Energy Assessment in Water Distribution Networks. Procedia Eng. 2015, 119, 139–147. [Google Scholar] [CrossRef] [Green Version]
  17. Sarbu, I. A Study of Energy Optimisation of Urban Water Distribution Systems Using Potential Elements. Water 2016, 8, 593. [Google Scholar] [CrossRef] [Green Version]
  18. Gómez, E.; Del Teso, R.; Cabrera, E.; Cabrera, J.E.; Soriano, J. Labeling Water Transport Efficiencies. Water 2018, 10, 935. [Google Scholar] [CrossRef] [Green Version]
  19. Mamade, A.; Sousa, C.; Marques, A.; Loureiro, D.; Alegre, H.; Covas, D. Energy Auditing as a Tool for Outlining Major Inefficiencies: Results from a Real Water Supply System. Procedia Eng. 2015, 119, 1098–1108. [Google Scholar] [CrossRef] [Green Version]
  20. Mamade, A.; Loureiro, D.; Alegre, H.; Covas, D. A comprehensive and well tested energy balance for water supply systems. Urban Water J. 2017, 14, 853–861. [Google Scholar] [CrossRef]
  21. Mamade, A.; Loureiro, D.; Alegre, H.; Covas, D. Top-Down and Bottom-Up Approaches for Water-Energy Balance in Portuguese Supply Systems. Water 2018, 10, 577. [Google Scholar] [CrossRef] [Green Version]
  22. Cabrera, E.; Del Teso, R.; Gómez, E.; Cabrera, E.; Estruch-Juan, E. Deterministic Model to Estimate the Energy Requirements of Pressurized Water Transport Systems. Water 2021, 13, 345. [Google Scholar] [CrossRef]
  23. Lenzi, C.; Bragalli, C.; Bolognesi, A.; Artina, S. From energy balance to energy efficiency indicators including water losses. Water Supply 2013, 13, 889–895. [Google Scholar] [CrossRef]
  24. Dziedzic, R.; Karney, B.W. Energy Metrics for Water Distribution System Assessment: Case Study of the Toronto Network. J. Water Resour. Plan. Manag. 2015, 141, 04015032. [Google Scholar] [CrossRef] [Green Version]
  25. Wong, H.G.; Speight, V.L.; Filion, Y.R.; Wong, H.G. Impact of Urban Development on Energy Use in a Distribution System. J. Am. Water Work. Assoc. 2017, 109, E10–E18. [Google Scholar] [CrossRef]
  26. Lapprasert, S.; Pornprommin, A.; Lipiwattanakarn, S.; Chittaladakorn, S. Energy Balance of a Trunk Main Network in Bangkok, Thailand. J. Am. Water Work. Assoc. 2018, 110, E18–E27. [Google Scholar] [CrossRef]
  27. Lipiwattanakarn, S.; Kaewsang, S.; Pornprommin, A.; Wongwiset, T. Real benefits of leak repair and increasing the number of inlets to energy. Water Pract. Technol. 2019, 14, 714–725. [Google Scholar] [CrossRef]
  28. Alegre, H.; Baptista, J.M.; Cabrera, E.; Cubillo, F.; Duarte, P.; Hirner, W.; Merkel, W.; Parena, R. Performance Indicators for Water Supply Services; IWA Publishing: London, UK, 2006. [Google Scholar]
  29. Rossman, L.A. EPANET 2 Users Manual, Water Supply and Water Resources Division; USEPA: Cincinnati, OH, USA, 2000.
Figure 1. Two simplified energy balance concepts in water networks considering only leakage (as water loss) and pipe friction, where (a) is proposed by Cabrera et al. [10] and (b) is proposed by Mamade et al. [19,20,21].
Figure 1. Two simplified energy balance concepts in water networks considering only leakage (as water loss) and pipe friction, where (a) is proposed by Cabrera et al. [10] and (b) is proposed by Mamade et al. [19,20,21].
Water 13 01011 g001
Figure 2. Energy balance for single pipe network.
Figure 2. Energy balance for single pipe network.
Water 13 01011 g002
Figure 3. Relationship between parameter pn and ratio of water loss, p.
Figure 3. Relationship between parameter pn and ratio of water loss, p.
Water 13 01011 g003
Figure 4. Energy balance for simplified branched pipe network with uniformly distributed demand nodes, where m is number of branches and j is number of demand nodes in each branch.
Figure 4. Energy balance for simplified branched pipe network with uniformly distributed demand nodes, where m is number of branches and j is number of demand nodes in each branch.
Water 13 01011 g004
Figure 5. Examples of water networks, where (a,b) are networks with 1 inlet, and (c,d) are networks with 2 inlets. Red dashed circles show the inlets connecting to the networks by imaginary valves with no friction loss.
Figure 5. Examples of water networks, where (a,b) are networks with 1 inlet, and (c,d) are networks with 2 inlets. Red dashed circles show the inlets connecting to the networks by imaginary valves with no friction loss.
Water 13 01011 g005
Figure 6. Relationships between energy balance components calculated by EPANET models and basic network parameters, where (a) E l , m o d vs. p, (b) E f , m o d vs. Δ H and (c) E W L , m o d vs. p.
Figure 6. Relationships between energy balance components calculated by EPANET models and basic network parameters, where (a) E l , m o d vs. p, (b) E f , m o d vs. Δ H and (c) E W L , m o d vs. p.
Water 13 01011 g006
Figure 7. Comparison between energy balance components calculated by EPANET models and ones estimated by theory, where (a) E l , m o d vs. E l , t h e o , (b) E f , m o d vs. E f , t h e o and (c) E W L , m o d vs. E W L , t h e o .
Figure 7. Comparison between energy balance components calculated by EPANET models and ones estimated by theory, where (a) E l , m o d vs. E l , t h e o , (b) E f , m o d vs. E f , t h e o and (c) E W L , m o d vs. E W L , t h e o .
Water 13 01011 g007
Figure 8. Boxplots of normalized demand distribution ( Q i / Q i n ) vs. normalized head loss distribution ( Δ H i / Δ H ), where (a) is for networks with 1 inlet and (b) is for networks with 2 inlets. Circles show the outliers.
Figure 8. Boxplots of normalized demand distribution ( Q i / Q i n ) vs. normalized head loss distribution ( Δ H i / Δ H ), where (a) is for networks with 1 inlet and (b) is for networks with 2 inlets. Circles show the outliers.
Water 13 01011 g008
Table 1. Characteristics of district metering areas.
Table 1. Characteristics of district metering areas.
IDNo. of
Inlets
No. of CustomersLengthAvg. DAvg. SfWater Loss, pΔH*
(km)(mm)(m/km)(%)(%)
11266924.51610.1737.116.6
21265726.41470.1728.635.7
31439952.31480.0844.614.1
41262646.41740.2038.544.2
51359454.71390.1144.218.4
61481251.01430.3654.940.7
71460743.21300.1732.444.8
81169528.82080.0912.913.7
91363418.11830.1629.714.4
101182022.51320.142.814.8
112192122.21660.5030.035.0
122215119.01540.1650.97.6
132229724.91540.2231.925.9
14273917.31910.3333.932.0
152146815.91780.707.734.8
162420447.41530.4836.337.0
17211,545129.61500.1430.745.3
182446073.71800.1530.065.3
192495751.51430.0731.218.4
202389747.41540.2847.227.8
Avg.1.5350840.81590.2332.829.3
Table 2. Performance of proposed theoretical methods to evaluate energy balance components before and after considering normalized head loss Δ H with calibrated coefficient C m j .
Table 2. Performance of proposed theoretical methods to evaluate energy balance components before and after considering normalized head loss Δ H with calibrated coefficient C m j .
ComponentEquationNo. of Inlets Value   of   C m j rRMSE (%)
BeforeAfterBeforeAfter
E l (27)10.74660.9390.9908.571.65
20.50470.9570.9856.021.91
E f (28)11.08330.9050.9057.326.92
20.75380.8340.83411.306.99
E W L (29)10.42190.9920.9944.751.83
20.30950.9780.9844.462.17
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lipiwattanakarn, S.; Kaewsang, S.; Charuwimolkul, N.; Changklom, J.; Pornprommin, A. Theoretical Estimation of Energy Balance Components in Water Networks for Top-Down Approach. Water 2021, 13, 1011. https://doi.org/10.3390/w13081011

AMA Style

Lipiwattanakarn S, Kaewsang S, Charuwimolkul N, Changklom J, Pornprommin A. Theoretical Estimation of Energy Balance Components in Water Networks for Top-Down Approach. Water. 2021; 13(8):1011. https://doi.org/10.3390/w13081011

Chicago/Turabian Style

Lipiwattanakarn, Surachai, Suparak Kaewsang, Natchapol Charuwimolkul, Jiramate Changklom, and Adichai Pornprommin. 2021. "Theoretical Estimation of Energy Balance Components in Water Networks for Top-Down Approach" Water 13, no. 8: 1011. https://doi.org/10.3390/w13081011

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