Next Article in Journal
Photochemically Induced Changes of Dissolved Organic Matter in a Humic-Rich and Forested Stream
Next Article in Special Issue
Nonlinear Dynamic Modeling of Urban Water Consumption Using Chaotic Approach (Case Study: City of Kelowna)
Previous Article in Journal
Impacts of Tile Drainage on Phosphorus Losses from Edge-of-Field Plots in the Lake Champlain Basin of New York
Previous Article in Special Issue
Sediment Transport in Sewage Pressure Pipes, Part II: 1 D Numerical Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Bi-Objective Approach for Optimizing the Installation of PATs in Systems of Transmission Mains

1
Dipartimento di Ingegneria Civile e Architettura, Università degli Studi di Pavia, Via Ferrata 3, 27100 Pavia, Italy
2
Dipartimento di Ingegneria Civile e Architettura, Università degli Studi di Catania, Via Santa Sofia 64, 95123 Catania, Italy
*
Authors to whom correspondence should be addressed.
Water 2020, 12(2), 330; https://doi.org/10.3390/w12020330
Submission received: 11 December 2019 / Revised: 20 January 2020 / Accepted: 21 January 2020 / Published: 23 January 2020
(This article belongs to the Special Issue Advances in Modeling and Management of Urban Water Networks)

Abstract

:
This paper proposes the bi-objective optimization for the installation of pumps operating as turbines (PATs) in systems of transmission mains, which typically operate at steady flow conditions to cater to tanks in the service of water distribution networks. The methodology aims to find optimal solutions in the trade-off between installation costs and generated hydropower, which are to be minimized and maximized, respectively. While the bi-objective optimization is carried out by means of a genetic algorithm, an inner optimization sub-algorithm provides for the regulation of PAT settings. The applications concerned a real Italian case study, made up of nine systems of transmission mains. The methodology proved able to thoroughly explore the trade-off between the two objective functions, offering solutions able to recover hydropower up to 83 KW. In each system considered, the optimal solutions obtained were postprocessed in terms of long-life net profit. Due to the large geodesic elevation variations available in the case study, this analysis showed that, in all systems, the optimal solution with the highest value of generated hydropower was the most profitable under usual economic scenarios, with payback periods always lower than 3 years.

1. Introduction

The contribution of the water industry is not negligible in the total demand of energy in the world [1]. Within this contribution, a major role is played by the extraction and treatment of source water [1], mainly due to the operation of pumping stations. To reduce this component of energy demand, water utilities are more and more often implementing practices to reduce wastes of water in distribution networks [2]. On the other hand, attempts have recently been made to recover energy from water networks and thus to compensate part of the intensive energy consumption associated with water supply. Energy recover can be accomplished by installing devices such as hydro-turbines [3,4,5,6,7] or pumps operating in reverse mode, also called pumps as turbines (PATs) [8,9,10,11,12], within pressurized pipes. Though featuring lower efficiencies than hydro-turbines, PATs feature much lower costs, thus offering themselves as a very cost-effective solution for energy recovery from water distribution networks.
The analysis of the scientific and grey literature concerning the use of hydro-turbines and PATs for energy recovery shows that research efforts have been made mainly in three different directions. The first direction is concerned with the analysis of the performance of the pump operating in reverse-mode, in terms of head-discharge and efficiency-discharge relationships. In fact, although acknowledging the possibility of using their devices in turbine mode, manufacturers of pumps do not usually provide the performance relationships. Therefore, several works in recent years were dedicated to analyzing, experimentally (e.g., [13,14]) or numerically through computational fluid dynamics (e.g., [15,16]), how the head-discharge and efficiency-discharge relationships of pumps relate to the relationships of the same devices operated in reverse mode, as well as to setting up procedures to derive the latter from the former (e.g., [17,18]).
The second research direction concerns the modeling, operation, and regulation of energy recovery devices in pressurized water networks in an attempt to both generate hydropower and reduce service pressure. In fact, the hydraulic regulation uses a bypass line and control valves to regulate flow, such that the PAT operates at or near the best efficiency point (BEP) [10]. Furthermore, the electric regulation of the PAT may be implemented by connecting the PAT to an inverter [11,12] to adjust the impeller rotational speed as a function of inflow and available head drop. Real time control algorithms can be applied to dynamically adjust the settings of PAT, control valves, and inverter under time varying operational conditions [19,20].
The third main research direction, to which the present paper belongs, concerns the optimization of sites and settings of energy recovery devices in pressurized pipes. In this context, energy production can also be optimized in combination with service pressure to help mitigating leakage in pressurized water networks. In this context, Giugni et al. [21] showed the application of a genetic algorithm to optimize location and settings. These authors considered the following two simplifying assumptions:
  • the number of energy recovery devices is a surrogate for their cost;
  • these devices have unitary efficiency (that is, their performance curves are neglected).
Additionally, the subsequent works by Corcoran et al. [22], Coelho and Andrade-Campos [23], and Fernández García and McNabola [24] use the number of devices as a surrogate for their total cost in the optimization phase, although removing the simplifying assumption of unitary efficiency made in [21]. Notably, Corcoran et al. [22] compared the effectiveness of three different kinds of algorithms, namely nonlinear programming, mixed integer nonlinear programming, and evolutionary optimization, in the maximization of produced hydropower. Coelho and Andrade-Campos [23] proposed a two-step optimization, which was carried out in a decoupled way. The first step concerned identification of site locations, which was performed using an optimization algorithm tailored to identifying the best locations of pressure reducing valves. A second manual step followed, aiming to identify the turbine models that could be installed at each site. Finally, a financial analysis was carried out to look at installation cost and at income associated with the sale of produced hydropower. Fernández García and Mc Nabola [24] proposed a nonlinear programming algorithm based on the sequential addition of energy recovery devices along the water paths present in a water distribution network. Device installation was iterated until there was no more significant power to be recovered from the network. A financial analysis was also carried out by the two authors to evaluate installation cost, annual income associated with PAT operation, and investment payback period. Unlike the works in [21,22,23,24], the works of Fecarotta and McNabola [25] and Tricarico et al. [26] proposed considering the installation cost in the context of the optimization, rather than the number of energy recovery devices. In fact, Fecarotta and McNabola [25] proposed an objective function combining installation costs and annual net income associated with device operation to be used inside a nonlinear mixed integer programming algorithm. The more recent work by Tricarico et al. [26] used a similar objective function into the context of a multi-objective optimization, leading to the identification of optimal sites for energy recovery devices. However, the device settings were not optimized in the methodology proposed in [26].
In this paper, a new multi-objective optimization algorithm is proposed for the installation of energy recovery devices in pressurized water networks, with the main focus on the installation of PATs in systems of transmission mains. As compared to work in [26], the novel algorithm enables optimization of both installation sites and settings. The algorithm is applied with two objective functions, namely the total installation cost of the devices and the generated hydropower, to be minimized and maximized, respectively. The Pareto fronts obtained with the algorithm lend themselves to being postprocessed in various economic scenarios to help decision makers in the choice of the ultimate solution.
In the following sections, first the methodology is described, followed by the applications to a real case study.

2. Methodology

2.1. PAT Installation

The methodology presented in this paper can be applied to systems of transmission mains that connect water sources to tanks feeding water distribution networks (Figure 1a). These systems are typically branched and operate under steady flow conditions over the whole day. While water distribution tanks are fed by transmission mains with a daily constant water discharge, their storage enables balancing the peaks of network consumption. In existing systems of transmission mains, the maximum water discharges depend only on the available heads at sources and tanks. In some cases, the values are larger than required or limited by source water availability. Therefore, the desired water discharges are enforced by valve regulation, which dissipates energy surpluses.
The redundant energy can be recovered through installation of rotary machines for electric power generation. The choice of number, sites, and models of machines represents an optimization problem aimed at minimizing installation costs and at maximizing the produced hydropower, under the hydraulic constraint of desired water discharges at mains, of fixed heads at source and tank nodes, and of desired heads at internal nodes. It is worth noting that the optimal solutions may not use the whole surplus of energy. In this case, the remaining surplus is dissipated through regulation valves.
In this work, the layout shown in Figure 1 was considered for PAT installation at the generic system main. The PAT was installed in-line with the system main inside a vault. Actually, the asynchronous PAT motor was connected to an inverter, which was in turn connected to the electrical grid (Figure 1b). The presence of the inverter enabled modification of the speed setting of the PAT and, therefore, the maximization of its production of hydropower for each daily value of water discharge in the main. Along each main of the system, a regulation valve was assumed to be present outside the vault to regulate flows as a function of fixed head values at system tanks and of desired head values at inner system nodes. Indeed, while the PAT was regulated to maximize the produced hydropower, some residual head surplus could be dissipated in the regulation valve. At a generic value of the PAT rotation speed, let Q (m3/s) be the water discharge in the PAT. As a result of the PAT curves of performance, the PAT produces a head drop YPAT (m) and features a certain efficiency η (-), based on its performance curves. The hydropower PPAT (KW) generated by the PAT is
P P A T = ρ   g   Q   Y P A T   η 1000 ,
where ρ = 1000 kg/m3 is water density, and g = 9.81 m/s2 is the gravity acceleration.
Finally, let Yv (m) be the head drop across the regulation valve assumed present along the main line outside of the PAT vault, and let J (-) be the energy friction slope in the main, evaluated by means of the Strickler formula. The friction-related head loss in the main is equal to J L, where L (m) is the length of the main.

2.2. Optimization Algorithm

The multi-objective genetic algorithm NSGAII [27] was used in this work to search for solutions in the trade-off between two objective functions, namely the total installation cost Cinst,tot (€) of the PATs (to be minimized) and the total hydropower PPAT,tot (KW) generated (to be maximized) in each system of transmission mains. While the presence of the inverter enables PAT operation to adjust to yearly variable values of the daily constant water discharge, a yearly average value of the water discharge was considered in the optimizations of the present work.
The first objective function Cinst,tot is calculated as the sum of the installation costs of the PAT at the various mains, i.e.,
C i n s t , t o t = i = 1 n m C i n s t , i
where Cinst,i (€) is the cost of the PAT installed in the i-th of the nm mains. In the case of no PAT installation suggested by the optimizer at the generic main, Cinst,i = 0.
Indeed, the total installation cost should include other issues such as the costs associated with the approaching of the electrical grid to the PAT site in the case when the grid is not available at close sites.
The second objective function PPAT,tot is calculated as the sum of the hydropower contributions at each main, i.e.,
P P A T , t o t = i = 1 n m P P A T , i
where PPAT,i is the hydropower generated in the i-th of the nm mains. In the case of no PAT installation suggested by the optimizer at the generic main, PPAT,i = 0.
Inside NSGAII, each individual is encoded in a set of nm genes, i.e., decisional variables. Specifically, the i-th of the nm genes represents the ID of the PAT installed in the i-th of the nm system mains.
For each solution (set of PATs installed in the mains) proposed by NSGAII, an inner optimization is performed through an embedded single objective genetic algorithm (sub-algorithm) to search for the PAT revolution speeds that might maximize the total hydropower generation in the system. Starting from the curves H(Q) and P(Q) provided by the manufacturer, the corresponding curves H(Q) and P(Q) curves for the PAT were estimated using the methodology proposed in [13]. In this methodology, first the BEP of a PAT is estimated starting from the BEP of the pump operated in direct mode. Then, the performance of the PAT in terms of H(Q) and P(Q) away from the BEP is obtained through polynomial equations.
For each solution proposed by NSGAII and for each set of PAT speeds proposed by the sub-algorithm, a hydraulic simulator based on the global gradient algorithm [28] is used to search for feasible values of Yv (i.e, Yv ≥ 0) for the regulation valves present at the system mains, while respecting the minimum desired heads at the inner nodes and the fixed heads at the system tanks.
The flowchart of instructions carried out for each individual in NSGAII is shown in Figure 2. First, the PATs are located in the system based on individual genes. Then, the total cost Cinst,tot of the PATs installed in the system, which is the first objective function of the optimization, is calculated. Then, a sub-algorithm is applied to the PATs for estimating PAT speeds. Then, the total power PPAT,tot generated by the system is calculated, with a penalization applied to infeasible solutions (i.e., solutions with negative value(s) of Yv).

2.3. Economic Postprocessing

The optimization described in the previous subsection produces a Pareto front of optimal solutions, each of which features its own distribution of PATs in the system mains, in the trade-off between Cinst,tot and PPAT,tot. For the choice of the ultimate solution in the Pareto front, an economic criterion can be applied. In this context, the net profit Cp (€) of each solution can be calculated as follows:
C P = C s a l e C i n s t , t o t C a d d
where Csale (€) is the gross profit associated with the sale of hydropower, converted into value at the initial time. This variable can be evaluated as
C s a l e = i = 1 T 8760   c   P P A T , t o t ( 1 + r ) i
where c (€/KW h), T (years), and r (-) are the unit selling price of the electrical energy produced, the useful life of the PAT, and the discount rate for price conversion, respectively. In Equation (5), the numerical value 8760 is the number of hours in a year. Cadd (€) represents potential additional costs associated with the PAT (e.g., maintenance costs). Finally, the solution with the highest Cp value in the Pareto front can be selected as the ultimate solution.
Another meaningful postprocessing variable is the payback period, i.e., the amount of time Tb (years) it takes to recover the cost of the investment. While authors in [23,24] used a simplified calculation of Tb, which did not consider the decrease in the net present value as a function of the selling year, Tb was calculated here iteratively, by searching for the lowest value of T for which Cp is larger than 0 in Equation (4).

3. Applications

3.1. Case Study

The described methodology was applied to the system of networks operated by Azienda Consorziale Servizi Etnei (ACOSET), the major water distribution manager in Catania, Italy. The system was already analyzed in [3] and is characterized by transmission mains that cater to twenty urban centers situated at the foot of the volcano Etna. The main sources for these mains are the Ciapparazzo extraction tunnel, the Maniace source, and numerous ACOSET and private wells (Figure 3).
For the calculations in this paper, nine separate systems of water transmission mains (features reported in Table 1) were considered, in which no PATs are currently present, namely, the systems feeding Ciapparazzo, Adrano, Gravina, Mascalucia, S. Gregorio, S. M. La Stella, Maniace, Camporotondo, and S. Giovanni Galermo. As an example, the qualitative elevation layout of the Ciapparazzo system is shown in Figure 4.
Due to the high geodesic gradients, the considered systems of transmission mains always show a surplus of energy that is currently dissipated through already installed regulation valves (see Figure 4 as an example). As a virtuous alternative option to the full dissipation of the available surplus of energy, the installation of PATs in this system is proposed in this work.
For the pumps to be installed as turbines in the systems, the 84 pumps of the NM4 series in the catalogue of the manufacturer Calpeda were considered. In their normal operation as pumps, they operate at 1450 rpm in the 50-Hz electrical grid. In their operation as PATs, the device revolution speed was assumed to be variable within the range 1450 to 3000 rpm. Furthermore, the electrical grid was assumed to be available at all sites, thus incurring no extra installation costs.
The cost of the electro-mechanical equipment for the generic PAT was obtained starting from the Calpeda price list and by applying some steps of increase, including the replacement of the asynchronous motor, the installation of an inverter, and a standard amount of +10% for additional hydraulic and electric equipment (e.g., gates and by-pass pipe, and electric wiring and meter). As for the replacement of the asynchronous motor, this must be done because of the different performance curves of the PAT from the same device used in direct mode [29]. Since the cost of the electro-mechanical equipment typical adds up to about 40% of the total installation cost [30], the latter was obtained by multiplying the former by 100/40 = 2.5.

3.2. Results

3.2.1. Optimization

The methodology described above yielded the Pareto fronts shown in the graphs in Figure 5. As expected, each front reports increasing values of the produced hydropower PPAT,tot as Cinst,tot grows. Overall, the maximum values of PPAT,tot, at the right end of the Pareto fronts, ranged from about 6 KW to 83 KW in the systems of Camporotondo and Ciapparazzo, respectively. As expected, the larger maximum values of PPAT,tot were associated with larger geodesic height variations and larger water discharges available in the system.
The analysis of the results revealed irregularities in the growth of PPAT,tot as a function of Cinst,tot, with changes in slope in the Pareto front for all systems. To shed light on this aspect, the Pareto front in the system of Camporotondo, made up of nine solutions (see Table 2), was taken as an example. The first solution was made up of no PAT in the system. As Cinst,tot grew, the optimizer first suggested installing only one PAT in main three, which was the last main of the system with the lowest water discharge. For a larger value of Cinst,tot = 6468 €, the installation of the PAT at main two was proposed. When the cost further increased, combinations of two PATs, involving first mains one and three and then mains two and three, were proposed. Finally, the two most expensive solutions featured PATs installed at all sites. Globally, the only PAT model suggested for main one was NM4 40/16C, which was operated at different speeds in the various solutions. The same consideration could be made for main two, in which the only PAT model installed was the NM4 32/20A. Various PAT models, operated at different speeds, were proposed, instead, for main three. For growing values of Cinst,tot, the change in the PAT model and PAT speed, as well as the introduction of the PAT in new sites, was responsible for slope changes in the Pareto front PPAT,tot(Cinst,tot).
The solutions obtained in the various systems for the highest values of Cinst,tot (and therefore of PPAT,tot) are presented in Table 3, which reports the features in terms of installed PAT models and adopted speed settings, hydraulic variables, and produced hydropower. As an example of these solutions, the system of S. M. La Stella is shown in Figure 6, also including the qualitative pattern of the head grade line along the system mains. This pattern points out that the head decreased due to friction losses in the mains as well as to local head dissipations at the two PATs and at the two regulation valves.

3.2.2. Economic Postprocessing

The Pareto fronts of optimal solutions shown in Figure 5 were postprocessed to obtain curves of net profit Cp as a function of the total installation cost Cinst,tot (see Figure 7). In these calculations, a useful life T = 20 years was considered for the PAT. Cadd in Equation (4) was assumed to include the annual maintenance costs, which were set at 3% of Cinst,tot and converted into present value at year 0. As an example, the Pareto front for the system of Camporotondo (reported in Figure 5c) was postprocessed by considering four alternative scenarios.
The main features of the four scenarios are:
  • scenario 1, with c = 0.1561 €/KW h, r = 0.04;
  • scenario 2, with c = 0.1561 €/KW h, r = 0.02;
  • scenario 3, with c = 0.1561 €/KW h, r = 0.06;
  • scenario 4, with c = 0.08 €/KW h, r = 0.04.
In scenario 1, the values of c and r represent currently typical conditions in Italy for the sale of hydropower and for the discount rate for cost conversion, respectively. Scenario 2 differs from scenario 1 in the lower value of r, thus resulting in higher present worth values for sold hydropower and for maintenance costs. Compared to scenario 1, scenario 3 has a larger value of r, resulting in lower present worth values for sold hydropower and for maintenance costs. Finally, scenario 4 differs from scenario 1 because of its lower value of c.
Overall, Figure 7 shows that, in scenario 1, Cp tends to grow, over almost the whole range of Cinst,tot values, up to a value of about 81,000 €. Therefore, the most profitable solution in terms of Cp is that featuring the highest value of Cinst,tot = 18,293 € (solution 9 in Table 2). In comparison with scenario 1, due to the different values of r, the pattern of Cp(Cinst,tot) in scenarios 2 and 3 is slightly higher and lower, respectively. In both scenarios, the superiority of solution 9 persists because of its large generated hydropower. In fact, for all the considered values of r, the sale of hydropower pays back the large investment and maintenance costs associated with solution 9. Under scenario 4, Cp rapidly grows up to a quite stable value around 28,000 €, and the maximum of Cp is at 30,964 €, obtained for Cinst,tot = 10,388 € (solution 6 in Table 2). In this case, due to the low value of c, the large investment and maintenance costs of solution 9 are less paid back in the long run, in comparison with other solutions.
Overall, the economic postprocessing in the system of Camporotondo shows that the most expensive solution in terms of Cinst,tot yielded by the optimization tends to be the most profitable solution unless the economic scenario used for the postprocessing differs largely from the typical ones in terms of c.
Scenario 1 for the economic postprocessing was applied for all the considered systems. In all systems, the most profitable solution was the right end solution in the Pareto front in Figure 5, thus confirming that the largest optimal investment is well paid back by the sale of electrical energy. Overall, the net profit values Cp ranged from about 8.1 × 104 € to about 1.5 × 106 €, for the systems of Camporotondo and Ciapparazzo, respectively. Higher values of net profit were associated with larger values of the produced hydropower PPAT,tot. The payback period Tb was lower than 2 years in all cases except for the system of Camporotondo, for which Tb = 3 years. Overall, these values of Tb were even lower than those observed in [24], due to the more favorable conditions of geodesic elevation variations in the present case study.

4. Discussion and Conclusions

In this paper, the bi-objective genetic optimization was proposed for tackling the problem of optimal installation of PATs in existing systems of transmission mains. The proposed methodology is able to search for optimal locations and speed settings of PATs in each system in the trade-off between the total installation cost and the produced hydropower, to be minimized and maximized, respectively. The applications concerned various topologically simple systems of transmission mains managed by ACOSET in Italy, while including in the analysis the estimated installation costs and the expected reverse-mode performance of the pumps produced by a manufacturer.
The application of the bi-objective optimization enabled Pareto fronts of optimal PAT installation solutions in the trade-off between installation costs and hydropower produced to be obtained. Even if the optimization of PAT locations and settings was performed in this case on yearly average water discharges at the mains, the presence of the inverter in the installation enables adjustment of the operation of each PAT to variable values of the daily water discharge in the year. The analysis of these Pareto fronts revealed that, as the cost grows, the power produced grows irregularly. In fact, as the installation cost varies in the generic system, optimal solutions are found with variable numbers, models, and speeds of PATs. Based on this result, it is expected that algorithms based on the sequential addition of PATs and keeping PAT speeds unchanged during the iterations, such as that proposed in [24], may lead to sub-optimal solutions. However, it must be remarked that the application of the bi-objective optimization took computational time with order of magnitude of hours/days for the simple case systems considered in this work. Therefore, computational times would become prohibitive for more topologically complex systems.
Starting from the trade-off between installation costs and generated hydropower, which is explored in a single optimization on each system, a quick economic postprocessing enables analysis of how optimal solutions perform in terms of long-cycle net profit and payback periods in various economic scenarios, thus supporting decision makers in the choice of the ultimate design solution. This marks a significative improvement compared to the single objective optimization, in which one optimization must be carried out for each considered economic scenario.
The following conclusions can be drawn from this work:
  • the installation of PATs confirms itself as a valuable solution for energy recovery from pressurized water networks;
  • in the case-study considered, made up of transmission mains featuring high geodesic elevation variations, the payback period for the installation investment ranges from 1 to 3 years, lower than those previously remarked in the scientific literature;
  • in the case-study considered, the most profitable solutions at all systems in typical economic scenarios are those associated with the maximum hydropower recovered, which ranges from about 6 KW to 83 KW.
Future work will be dedicated to making the methodology feasibly applicable to more topologically complex systems under time varying operating conditions.

Author Contributions

Conceptualization, E.C. and G.P.; methodology, E.C. and G.G.; writing—original draft preparation, E.C. and G.G.; writing—review and editing, A.C., C.C., and G.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IEA (International Energy Agency). Water Energy Nexus: Excerpt from the World Energy Outlook. 2016. Available online: https://www.iea.org/publications/freepublications/publication/WorldEnergyOutlook2016ExcerptWaterEnergyNexus.pdf (accessed on 22 May 2019).
  2. Xu, Q.; Chen, Q.; Ma, J.; Blanckaert, K.; Wan, Z. Water saving and energy reduction through pressure management in urban water distribution networks. Water Resour. Manag. 2014, 28, 3715–3726. [Google Scholar] [CrossRef]
  3. Pezzinga, G.; Tosto, G. Adeguamento energetico di reti di adduzione idrica in pressione. L’Acqua 2001, 1, 27–34. (In Italian) [Google Scholar]
  4. Sammartano, V.; Sinagra, M.; Filianoti, P.; Tucciarelli, T. A Banki-Michell turbine for in-line water supply systems. J. Hydraul. Res. 2017, 55, 686–694. [Google Scholar] [CrossRef]
  5. Samora, I.; Manso, P.; Franca, M.; Schleiss, A.; Ramos, H. Energy recovery using micro-hydropower technology in water supply systems: The case study of the city of Fribourg. Water 2016, 8, 344. [Google Scholar] [CrossRef] [Green Version]
  6. Sinagra, M.; Sammartano, V.; Morreale, G.; Tucciarelli, T. A new device for pressure control and energy recovery in water distribution networks. Water 2017, 9, 309. [Google Scholar] [CrossRef] [Green Version]
  7. Bonthuys, G.J.; van Dijk, M.; Cavazzini, G. Water distribution system energy recovery and leakage reduction optimisation through hydro turbines. J. Water Resour. Plan. Manag. 2016, 142, 04015045. [Google Scholar]
  8. Ramos, H.; Borga, A. Pumps as turbines: An unconventional solution to energy production. Urban Water 1999, 1, 261–263. [Google Scholar] [CrossRef]
  9. Agarwal, T. Review of pump as turbine (PAT) for microhydropower. Int. J. Emerg. Technol. Adv. Eng. 2012, 2, 163–169. [Google Scholar]
  10. Carravetta, A.; Del Giudice, G.; Fecarotta, O.; Ramos, H. Energy production in water distribution networks: A PAT design strategy. Water Resour. Manag. 2012, 26, 3947–3959. [Google Scholar] [CrossRef]
  11. Carravetta, A.; Del Giudice, G.; Fecarotta, O.; Ramos, H. PAT design strategy for energy recovery in water distribution networks by electrical regulation. Energies 2013, 6, 411–424. [Google Scholar] [CrossRef] [Green Version]
  12. Alberizzi, J.C.; Renzi, M.; Nigro, A.; Rossi, M. Study of a Pump-as-Turbine (PaT) speed control for a Water Distribution Network (WDN) in South-Tyrol subjected to high variable water flow rates. Energy Procedia 2018, 148, 226–233. [Google Scholar] [CrossRef]
  13. Derakhshan, S.; Nourbakhsh, A. Experimental study of characteristic curves of centrifugal pumps working as turbines in different specific speeds. Exp. Therm. Fluid Sci. 2008, 32, 800–807. [Google Scholar] [CrossRef]
  14. Pugliese, F.; De Paola, F.; Fontana, N.; Giugni, M.; Marini, G. Performance of vertical-axis pumps as turbines. J. Hydraul. Res. 2018, 56, 482–493. [Google Scholar] [CrossRef]
  15. Nautiyal, H.; Kumar, V.; Thakur, A. CFD Analysis on Pumps working as Turbines. Hydro Nepal 2010, 6, 35–37. [Google Scholar] [CrossRef]
  16. Fecarotta, O.; Carravetta, A.; Ramos, H. CFD and comparisons for a pump as turbine: Mesh reliability and performance concerns. Int. J. Energy Environ. 2011, 2, 39–48. [Google Scholar]
  17. Marchiori, I.N.; Lima, G.M.; Brentan, B.M.; Junior, E.L. Effectiveness of methods for selecting pumps as turbines to operate in water distribution networks. Water Supply 2018, 19, 417–423. [Google Scholar] [CrossRef]
  18. Novara, D.; McNabola, A. A model for the extrapolation of the characteristic curves of pumps as turbines from a datum best efficiency point. Energy Convers. Manag. 2018, 174, 1–7. [Google Scholar] [CrossRef]
  19. Fontana, N.; Giugni, M.; Glielmo, L.; Marini, G. Real time control of a prototype for pressure regulation and energy production in water distribution networks. J. Water Resour. Plan. Manag. 2016, 142, 04016015. [Google Scholar] [CrossRef]
  20. Fontana, N.; Giugni, M.; Glielmo, L.; Marini, G.; Zollo, R. Hydraulic and Electric Regulation of a Prototype for Real-Time Control of Pressure and Hydropower Generation in a Water Distribution Network. J. Water Resour. Plan. Manag. 2018, 144, 04018072. [Google Scholar] [CrossRef]
  21. Giugni, M.; Fontana, N.; Ranucci, A. Optimal location of PRVs and turbines in water distribution systems. J. Water Resour. Plan. Manag. 2014, 140, 06014004. [Google Scholar] [CrossRef]
  22. Corcoran, L.; McNabola, A.; Coughlan, P. Optimization of water distribution networks for combined hydropower energy recovery and leakage reduction. J. Water Resour. Plan. Manag. 2015, 142, 04015045. [Google Scholar] [CrossRef]
  23. Coelho, B.; Andrade-Campos, A. Energy recovery in water networks: Numerical decision support tool for optimal site and selection of micro turbines. J. Water Resour. Plan. Manag. 2018, 144, 04018004. [Google Scholar] [CrossRef]
  24. Fernández García, I.; Mc Nabola, A. Maximizing Hydropower Generation in Gravity Water Distribution Networks: Determining the Optimal Location and Number of Pumps as Turbines. J. Water Resour. Plan. Manag. 2020, 146, 04019066. [Google Scholar] [CrossRef]
  25. Fecarotta, O.; McNabola, A. Optimal location of pump as turbines (PATs) in water distribution networks to recover energy and reduce leakage. Water Resour. Manag. 2017, 31, 5043–5059. [Google Scholar] [CrossRef]
  26. Tricarico, C.; Morley, M.S.; Gargano, R.; Kapelan, Z.; Savic, D.; Santopietro, S.; Granata, F.; De Marinis, G. Optimal energy recovery by means of pumps as turbines (PATs) for improved WDS management. Water Sci. Technol. Water Supply 2018, 18, 1365–1374. [Google Scholar] [CrossRef]
  27. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II. IEEE Trans. Evolut. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
  28. Todini, E.; Pilati, S. A gradient algorithm for the analysis of pipe networks. Computer Application in Water Supply. In System Analysis and Simulation; Coulbeck, B., Choun-Hou, O., Eds.; John Wiley & Sons: London, UK, 1988; Volume I, pp. 1–20. [Google Scholar]
  29. Novara, D.; Carravetta, A.; McNabola, A.; Ramos, H. Cost Model for Pumps as Turbines in Run-of-River and In-Pipe Microhydropower Applications. J. Water Resour. Plan. Manag. 2019, 145, 04019012. [Google Scholar] [CrossRef]
  30. Chapallaz, J.M.; Eichenberger, P.; Fischer, G. Manual on Pumps Used as Turbines; VIEWEG, 1992. Available online: http://skat.ch/book/manual-on-pumps-used-as-turbines-volume-11/ (accessed on 22 January 2020).
Figure 1. Example of transmission main (a). Layout of pump as turbine (PAT) installation, adapted from [11] (b).
Figure 1. Example of transmission main (a). Layout of pump as turbine (PAT) installation, adapted from [11] (b).
Water 12 00330 g001
Figure 2. Flow chart of the instructions to carry out for each individual of NSGAII.
Figure 2. Flow chart of the instructions to carry out for each individual of NSGAII.
Water 12 00330 g002
Figure 3. Planview of the water mains; adapted from Pezzinga and Tosto [3].
Figure 3. Planview of the water mains; adapted from Pezzinga and Tosto [3].
Water 12 00330 g003
Figure 4. Qualitative elevation layout of the system of transmission mains of Ciapparazzo. The numbers indicate the nodes in Table 1, which also reports the geometrical data of the system.
Figure 4. Qualitative elevation layout of the system of transmission mains of Ciapparazzo. The numbers indicate the nodes in Table 1, which also reports the geometrical data of the system.
Water 12 00330 g004
Figure 5. Pareto front of optimal solutions in the trade-off between Cinst,tot and PPAT,tot in the systems of (a) Ciapparazzo, Adrano, Mascalucia, (b) Gravina, S. Gregorio, S. G. Galermo, (c) S. M. La Stella, Manciace, and Camporotondo.
Figure 5. Pareto front of optimal solutions in the trade-off between Cinst,tot and PPAT,tot in the systems of (a) Ciapparazzo, Adrano, Mascalucia, (b) Gravina, S. Gregorio, S. G. Galermo, (c) S. M. La Stella, Manciace, and Camporotondo.
Water 12 00330 g005
Figure 6. Layout of the system of transmission mains of S. M. La Stella. Qualitative pattern for the head grade line for the solution with Cinst,tot = 12,021 € and PPAT,tot = 8.74 KW. See Table 3 for numerical values.
Figure 6. Layout of the system of transmission mains of S. M. La Stella. Qualitative pattern for the head grade line for the solution with Cinst,tot = 12,021 € and PPAT,tot = 8.74 KW. See Table 3 for numerical values.
Water 12 00330 g006
Figure 7. Curves of total net profit Cp as a function of Cinst,tot, obtained in scenarios 1–4, by postprocessing the Pareto front in Figure 5c for the system of Camporotondo.
Figure 7. Curves of total net profit Cp as a function of Cinst,tot, obtained in scenarios 1–4, by postprocessing the Pareto front in Figure 5c for the system of Camporotondo.
Water 12 00330 g007
Table 1. Features of the systems of transmission mains.
Table 1. Features of the systems of transmission mains.
SystemMainUpstream Nodeits,un (0/1)Downstream Nodeits,dn (0/1)Q (m3/s)Hun (m)Hdn (m)L (m)D (m)k (m1/3/s)
Ciapparazzo111210.475773310,3000.80075
221300.473371940900.80075
330400.32571971617000.80075
440500.377167155900.80075
550610.00871559510000.20075
661710.00859547816000.20075
750800.36271569537300.70075
880910.04695667300.20075
9801000.32269569547100.70075
101001110.0926956661700.40075
111001200.2369566523000.70075
121201300.276656055500.70075
131301410.0260549325000.15075
141301500.2560567013800.60075
151501610.17167068337900.40075
Adrano111210.0757337052500.60075
221310.06170561911500.60075
Gravina111200.0244354109000.15075
220310.0054103826000.08075
320410.01941032219000.15075
Mascalucia111200.04966558319000.25075
220310.0075835773500.15075
320400.04258348321000.25075
440510.0184834849000.15075
S Gregorio111200.03847639814000.20075
220300.01939830316000.15075
330410.013033896000.12575
430510.0130334716000.12575
S. M. La Stella111200.0168356936000.12575
220310.00556942631000.12575
Maniace111210.05565064516780.45090
221310.05564564136900.45090
331410.05564163917370.45090
441510.05563963817100.45090
551610.05563863516500.45090
661710.09563563127960.45090
771810.09563162828360.45090
881910.09562862520430.45090
9911010.09562562319280.45090
101011100.09562362014560.45090
111101210.05762061914560.40090
121211310.05761961719450.40090
131311410.05761761331160.40090
141411510.05761361023650.40090
151511610.00861047619000.25090
161511710.04961060725670.40090
171711810.04960760428730.40090
181811910.04960460228.80.40090
191912000.04960260014120.35090
202002110.02600599500.20090
212002210.02960059714120.35090
222212310.02959759419870.30090
232312410.02959459126960.30090
Camporotondo111200.00859956217500.15075
220300.00656249616500.12575
330410.0044964618000.10075
S. G. Galermo111200.03642035113000.20075
220310.0183512492200.15075
320410.01835125319000.15075
Footnote: the columns represent the system name, the ID of the main, the upstream node of the main, the absence (its,un = 0) or presence (its,un = 1) of a tank/source at the upstream node, the downstream node, the absence (its,dn = 0) or presence (its,dn = 1) of a tank/source at the downstream node, the flowing water discharge Q, the head at the upstream and downstream nodes (Hun and Hdn) (prescribed for a tank/source and minimum desired value for an intermediate node), the main length L, diameter D, and Strickler roughness coefficient k.
Table 2. Solutions in the Pareto front for the system of Camporotondo. Features in terms of PATs installed in the system, PAT speeds, installation cost Cinst,tot, and hydropower PPAT,tot.
Table 2. Solutions in the Pareto front for the system of Camporotondo. Features in terms of PATs installed in the system, PAT speeds, installation cost Cinst,tot, and hydropower PPAT,tot.
SolutionPAT in Main 1PAT in Main 2PAT in Main 3PAT Speed in Main 1PAT Speed in Main 2PAT Speed in Main 3Cinst,tot (€)PPAT,tot (KW)
1 00
2 NM4 25/12A 300039201.85
3 NM4 25/160B 300043913.43
4 NM4 32/20A 3000 64683.72
5NM4 40/16C NM4 25/160B2577 300097484.68
6 NM4 32/20ANM4 25/12A 2499198810,3884.79
7 NM4 32/20ANM4 32/20A 3000208312,9364.84
8NM4 40/16CNM4 32/20ANM4 32/16A25372510213217,0945.40
9NM4 40/16CNM4 32/20ANM4 32/20A23162510167118,2935.73
Table 3. Solutions obtained in the various systems for the highest values of Cinst,tot.
Table 3. Solutions obtained in the various systems for the highest values of Cinst,tot.
SystemMainPATCinst,i (€)Hact,un (m)Hact,dn (m)Yv (m)J L (m)YPAT (m)Speed (rpm)PPAT (kW)
Ciapparazzo1 756.50733.0013.589.92
2 733.00729.060.003.94
3 729.06727.980.001.08
4 727.98727.490.000.49
5NM4 32/16B5107727.49595.0072.210.6359.6530003.83
6NM4 32/16A5269595.00477.5030.261.0086.2430004.90
7 727.49721.500.006.00
8NM4 65/16S10,908721.50667.300.000.4753.72293216.71
9 721.50715.500.005.99
10NM4 125/25E19,697715.50665.502.560.3547.09300041.58
11 715.50714.010.001.49
12 714.01713.520.000.49
13NM4 40/16A6262.5713.52493.4015.7745.38158.96300015.62
14 713.52711.110.002.41
15 711.11683.001.3726.74
Total 47,243 82.65
Adrano1NM4 125/25A42,428733.00705.000.000.0427.96151116.05
2NM4 80/20C12,206705.00619.004.890.1280.99300041.24
Total 54,634 57.29
Gravina1 434.00410.470.0023.53
2NM4 40/16C5357410.47382.000.8519.458.1816180.30
3NM4 50/20C8784410.47322.000.9331.1356.4229767.27
Total 14,141 7.58
Mascalucia1NM4 80/20B13,675665.00588.850.0013.3062.85300024.76
2 588.85577.0011.180.67
3NM4 65/20A13,858588.85508.320.0011.0369.50300021.01
4 508.32484.0011.0813.23
Total 27,533 45.77
S Gregorio1NM4 80/20C12,206476.00421.710.0019.7834.51267210.97
2 421.71395.500.0026.21
3 395.50389.000.006.50
4NM4 40/16C5357395.50347.002.0917.3329.0830002.07
Total 17,563 13.04
S. M. La Stella1NM4 40/20B7630683.00569.3014.2038.9960.5128204.29
2NM4 25/160B4391569.30425.7012.617.53123.4530004.45
Total 12,021 8.74
Maniace1 650.00645.404.140.46
2 645.40641.403.001.00
3 641.40638.802.130.47
4 638.80637.700.630.47
5 637.70634.852.400.45
6 634.85631.001.582.27
7 631.00627.750.952.30
8 627.75625.300.791.66
9 625.30622.850.891.56
10 622.85620.001.671.18
11 620.00619.100.100.80
12 619.10616.901.131.07
13 616.90613.002.191.71
14 613.00610.201.501.30
15NM4 32/20A6468610.20475.505.770.25128.6830008.21
16 610.20607.002.161.04
17 607.00603.502.341.16
18 603.50602.300.041.16
19 602.30601.140.001.16
20 601.14599.002.000.14
21 601.14596.853.880.41
22 596.85593.851.691.31
23 593.85591.001.081.77
Total 6468 8.21
Camporotondo1NM4 40/16C5357599.00574.430.005.0819.4923161.20
2NM4 32/20A6468574.43497.190.007.1370.1225103.49
3NM4 32/20A6468497.19461.000.005.0531.1416711.03
Total 18,293 5.73
S. G. Galermo1NM4 65/16S10,908420.10362.760.0016.4940.85241011.79
2NM4 50/25B11,437362.76248.550.003.23110.98261612.91
3NM4 40/16C5357362.76253.000.0027.9481.8324247.74
Total 27,702 32.44
Footnote: the columns represent the system name, the main ID, the PAT model installed, the installation cost Cinst,i, the calculated upstream (Hact,un) and downstream (Hact,dn) head, the head drop Yv in the valve, the head loss J L in the main, the head drop YPAT in the PAT, the PAT speed setting, and produced hydropower PPAT.

Share and Cite

MDPI and ACS Style

Creaco, E.; Galuppini, G.; Campisano, A.; Ciaponi, C.; Pezzinga, G. A Bi-Objective Approach for Optimizing the Installation of PATs in Systems of Transmission Mains. Water 2020, 12, 330. https://doi.org/10.3390/w12020330

AMA Style

Creaco E, Galuppini G, Campisano A, Ciaponi C, Pezzinga G. A Bi-Objective Approach for Optimizing the Installation of PATs in Systems of Transmission Mains. Water. 2020; 12(2):330. https://doi.org/10.3390/w12020330

Chicago/Turabian Style

Creaco, Enrico, Giacomo Galuppini, Alberto Campisano, Carlo Ciaponi, and Giuseppe Pezzinga. 2020. "A Bi-Objective Approach for Optimizing the Installation of PATs in Systems of Transmission Mains" Water 12, no. 2: 330. https://doi.org/10.3390/w12020330

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