Energy Saving in Water Distribution Network through Pump as Turbine Generators: Economic and Environmental Analysis

Complex systems of water distribution networks (WDS) are used to supply water to users. WDSs are systems where a lot of distributed energy is available. Historically, this energy is artificially dissipated by pressure reduction valves (PRVs), thanks to which water utilities manage the pressure level in selected nodes of the network. The present study explores the use of economic hydraulic machines, pumps as turbines (PATs) to produce energy in a small network located in a town close to Palermo (Italy). The main idea is to avoid dissipation in favor of renewable energy production. The proposed study is applied to a WDN typical of the Mediterranean countries, where the users, to collect water during the period of water scarcity conditions, install private tanks. The presence of private tanks deeply modifies the network from its designed condition. In the proposed analysis, the economic benefit of PATs application in water distribution networks has been investigated, accounting for the presence of users' private tanks. The analysis, carried out by mean of a mathematical model able to dynamically simulate the water distribution network with PATs, shows the advantage of their installation in terms of renewable energy recovery, even though the energy production of PATs is strictly conditioned by their installation position.


Introduction
Water and energy savings are major concerns in the management of water distribution networks (WDNs).Nowadays, energy for pumping represents the main cost for water companies which have to guarantee the delivery of enough water with good quality to population: indeed, water losses reflect directly in energy wasting.This strong connection between water and energy consumption gives rise to programs, such as Watergy [1], aimed at investigating technical and managerial strategies to realize energy, water, and monetary savings through cost-effective actions.
A large effort has been made in recent years to investigate management strategies aimed at reducing net energy consumption in WDNs.Traditional solutions focus on pump and pipe selection, leak detection and repair, and pressure control policies.In fact, water and energy thrifts in WDNs are strictly linked to pipeline water leakages (pipe breaking, joint failures, etc.), which are related to water pressure.
Traditionally, in existing WDNs, head reduction has been achieved by inserting pressure reducing valves (PRVs) [2][3][4][5] that are able to maintain pressure variation in a controlled manner.Recently, there has been growing interest in recovering energy instead of dissipating excess pressure through valves, analyzing the potential of turbines to replace PRV or break-pressure tanks, achieving two main objectives: leakage reduction and renewable energy production.Unfortunately, the hydraulic turbine installation is a suitable solution just where the available hydraulic power is fairly constant.In fact, where flow rates and heads are highly variable, turbines are not cost effective, so other strategies have been investigated in recent years.
One of the easiest technological solutions, which can combine efficiency and economical convenience, is the use of a centrifugal pump in reverse mode, with the engine acting as a generator [6,7].Pumps as turbines (PATs) require low cost hydraulic and electric equipment and are often claimed as the cheapest and most sustainable solutions [8] for energy production [9][10][11], especially for small size power plants (<40-50 kW).In fact, pumps have various advantages compared to turbines, which are linked to ease of availability, proven technology, low initial installation and maintenance cost, availability for a wide range of heads, and flows [12,13].
The concept of electricity generation through reverse running centrifugal pumps is not new but dates back to around 80 years ago.Recently, Jain and Patel [12] analyzed the historical development and state-of-the-art of PATs, which found application in micro-hydropower plants as well as in water supply piping systems.
Pumps, which work in reverse mode, do not show the same performance when acting in direct (pump) and reverse (turbine) modes, and pumps manufacturers usually provide the characteristic curves only for the direct mode.The evaluation of the characteristic curves of PATs has been one of the most important limitations of their applicability [14][15][16][17][18].There have been many performance prediction techniques of the turbine operation of pumps based on best efficiency or on specific speed, which have tried to establish a correlation enabling the passage from the direct to the reverse characteristics (see [13,16] for a review).Derakhshan and Nourbakhsh [16] tested several centrifugal pumps and used experimental data to derive some relations to predict the best efficiency point characteristic of PAT.More recently, Nautiyal et al. [8] experimentally tested a centrifugal pump of specific speed to relate pump and turbine mode characteristics.They also used experimental data from other researchers to develop relations to obtain PATs efficiency from pump mode characteristics.These relations gave lower deviation from experimental results respect to those obtained by previous researchers [15][16][17].
A large number of theoretical and experimental studies have been done for prediction of PAT performance, but there is still a need to explore this area more deeply since the results predicted by these methods are not reliable for pumps with different specific speeds and capacities.
Most recent attempts to predict performance of PATs have been made using computational fluid dynamics (CFD).Unfortunately, in many cases, large deviations in results have been found between experimental and CFD results [13,15,18,19].
Recently, PATs have been investigated in terms of economic benefits.Arriaga [20] carried out the cost analysis for a pico-hydro plant development in the Lao People's Democratic Republic, showing that the implementation of PAT reduces the cost 53% as compared with cross-flow turbine.Carravetta et al. [19] proposed a method which allows PAT's geometry selection for a given flow-head distribution pattern to find the solution which maximizes the produced energy in WDNs.Carravetta et al. [21] extended the methodology to investigate on best economic efficiency between a hydraulic regulation via a series-parallel hydraulic circuit or an electrical regulation via inverter, finding that hydraulic regulation is more flexible and efficient than electrical regulation.Recently, Puleo et al. [22] investigated energy recovery by PATs in WDNs characterized by the presence of private tanks and intermittent service, highlighting that further analysis was needed for these operating conditions to investigate the possibility and the efficiency of the PATs to recover energy from the WDN.Fontana et al. [23] analyzed a water distribution network in the city of Naples to assess the potential revenues from energy recovery instead of dissipating excess pressure by replacing some of PRVs with PATs, which turned out to guarantee adequate pressure service.
They found that the optimal location of PRVs to reduce water losses did not maximize energy production.Maximum energy recovery needed a multi-objective approach.Nevertheless, results showed attractive profits and capital payback period (CPP), which further justify efforts in the analysis of PAT application in WDNs [21,24].
The aim of the present study is to investigate energy recovery in WDNs by PATs to guarantee adequate pressure service and recovery energy from the fluid instead of classical energy dissipation devices (i.e., PRVs).The analysis has been carried out in terms of energy and capital cost recovery.
CFD analysis has been used to analyze energy recovery in WDNs characterized by the presence of private tanks where users who experience water resources rationing store water resources (see [25,26]).Such local water storage volumes are actually very common in those countries (e.g., in the Mediterranean area) where water supply is not reliable and users are brought to collect water to cover their needs when the service in unavailable, due to shortage of supply or intermittent distribution [27].The local private storage interposed between the user and the network is subject to filling/emptying processes, which modify user's demand profile and pressure distribution in the network, causing modifications on PATs energy production capability in the network.Water heads and discharges are calculated through the method of characteristics (MOC).The model, already presented in De Marchis and Freni [28], was used to simulate PATs in small water distribution network.
The above-mentioned research focused the attention of the effect of PATs in water distribution network without actually taking into account the dynamic variation of the energy produced in light of the WDN users' behavior.The proposed mathematical model is able to modify the regulation point of the PAT with the discharge flowing in the WDN, simulating different scenarios and taking into account the presence of private tanks, if installed.The model can thus be considered a management tool for the water utilities to investigate on the possibility of recover energy in WDNs and to optimize the position of the PATs.With respect to other approaches, the proposed mathematical model which is able to simulate a lot of possible work combinations in WDNs (i.e., steady state condition, intermittent distribution, presence of water tanks, etc.).
The model has been applied to the Misilmeri network, a small town located close to Palermo (Italy), characterized by the presence of several private tanks on the roofs, and the results have been analyzed in terms of a management strategy based on the economic analysis and the payback period, taking into account users' satisfaction in terms of water volume supplied in each node for over 24 h.

Mathematical Model and Numerical Procedure
The mathematical tool has a main part given by the hydrodynamic model, able to simulate complex water distribution networks [25] and a second module able to reproduce the PAT operation mode [28].
The hydrodynamic model solves the momentum and continuity equations using MOC which transform the system of partial differential equations into ordinary differential equations, known as compatibility equations: where t is the time, V is the velocity averaged over the pipe cross-section, h is the water head, g is the acceleration due to gravity, c is the celerity of pressure waves and θ is the slope of the pipeline.In the present analysis, a steady state friction model has been used, thus J s , that is the head loss per unit length, has been calculated according to the classical formulation by Darcy-Weisbach.The formula is strictly dependent on the roughness of the pipe, which plays a fundamental role in the discharge estimation [29,30].The model is developed to take into account the unsteady friction model of Vítkovský et al. [31] (details can be found in De Marchis and Freni [28]), even though in the present analysis a steady state friction model was used.
Energies 2016, 9, 877 4 of 15 The compatibility equations are valid along the well known characteristic lines (positive and negative): Equations ( 1) and ( 2) can be solved through the finite difference technique and read: where V i,n+1 j and h i,n+1 j are the velocity and the water head in the jth section (of abscissa (j -1)L i /N i ) of the ith pipe at the time step t n + ∆t; θ i is the slope of the ith pipe; j m and j v are the sections upstream and downstream to the jth section, respectively.The time step advancement ∆t n i , function of the length and the celerity of the ith pipe, is calculated for each pipe and the minimum value is chosen as the unique time step integration: The flow rate into users' tanks is calculated as: where Q j,up is the discharge at the jth node entering the tank connected to the node, C v is the non-dimensional float valve emitter coefficient, a is the valve effective discharge area, g is the gravity acceleration, h i j is the water head at the jth node, and h i j,tank is the height of the private tank.Further details on the numerical model can be found in De Marchis et al. [25,27].
The accuracy of the hydrodynamic model was verified comparing numerical results to measured water heads in a real WDN of the city of Palermo (Italy).The comparison, shown in De Marchis et al. [32], clearly showed the ability of the numerical model to predict the variation of the pressure in 24 h.The qualitative agreement between measured and simulated pressure profiles were confirmed by the quantitative analysis carried out calculating the root square mean error (RSME), obtaining an accuracy of the model of about 5% of the average measured pressure.
To investigate on the energy recovery in water distribution networks, the hydrodynamic model was developed to simulate PAT behavior.Even though several analyses have been carried out to derive the characteristic curve of a water turbine [33], here the characteristic curves achieved by Derakhshan and Nourbakhsh [15,16], through experimental analysis, have been used: where H (m), P (W), and Q (m 3 /s) are head, power, and flow rate, respectively.Subscripts t and b are related to turbine and best efficiency point (BEF).Following, among others, Derakhshan and Nourbakhsh [15,16] PATs parameters are related to the specific rotational speed per second (rps).
In the present case of study, described in the next paragraph, a pipeline system is investigate, thus justifying the use of 1D model (MOC method and PAT curve module) to predict the energy recovery Energies 2016, 9, 877 5 of 15 in urban water distributions.The model is developed thus to work even for flow rate or the pressure conditions which are not compatible with the energy production by means of PATs.In these cases the PATs are by-passed and PRVs are simulated, according to the scheme shown in Figure 1, where turbines are coupled to valves [23].

The Case Study
The proposed numerical model has been applied to a small town close to Palermo city (Sicily), called Misilmeri, designed to supply about 22,000 inhabitants (11,300 users).The network is made of HDPE and cast iron pipes whose diameters range between 80 and 300 mm.Two main pipelines connected to two water reservoirs feed the network.The first tank, called Castle, is located at 190 m referred to the mean sea level (MSL) and has a volume of about 500 m 3 .The second reservoir, known as Cave, is located at 155 m MSL, having a total water volume of about 1,000 m 3 .
In Figure 2a, schematic representation of the water network is depicted.As shown, the pipeline system can be subdivided into two sub-network, supplied by the two water tanks.70 pipes and 57 nodes compose the network under investigation.In Appendix A (Tables A1 and A2), details on the distribution network (diameter, pipe length, water supplied, etc.) are reported.In Appendix A it is also reported the value of the roughness, clearly a key parameter, even though the calibration process (not shown here) pointed out that all the pipes fall into a hydraulically smooth regime, with values of equivalent sand grain roughness ks [34] in the range: 0.001 mm < ks < 0.005 mm.For this reason, the roughness ks was set equal to the average value (ks = 0.003 mm).
The water distribution network configuration is characterized by significant differences between the geodetic quote of water tanks and the users, suggesting high levels of energy recovery.The network orography thus seems to have an excellent configuration in order to test the application of micro energy recovery, introducing PATs in the network pipes.
One of the main problems in PATs positioning is to find the optimal installation point, thus to achieve both maximum energy production and desired water supply.Due to the small number of users, in the present case the choice of the optimal installation point has been carried out following the hydraulic regulation of the network.In small WDNs, the hydraulic regulation can be easily applied, achieving high performance levels, as demonstrated by De Marchis and Freni [28].
Taking advantage of the high variation of the geodetic elevation between the water tanks and the network, the main objective of the present analysis is to recover energy, in terms of water head, from water distribution networks.
At the moment, in WDNs the energy is dissipated through float valves.As specified in the following, the cost of electrically regulated valves is similar to the cost of a pump working in turbine mode, with the fundamental difference that PATs have the capability of producing small amounts of renewable energy.
Due to the specific case of study, three different scenarios have been explored.In the first PAT

The Case Study
The proposed numerical model has been applied to a small town close to Palermo city (Sicily), called Misilmeri, designed to supply about 22,000 inhabitants (11,300 users).The network is made of HDPE and cast iron pipes whose diameters range between 80 and 300 mm.Two main pipelines connected to two water reservoirs feed the network.The first tank, called Castle, is located at 190 m referred to the mean sea level (MSL) and has a volume of about 500 m 3 .The second reservoir, known as Cave, is located at 155 m MSL, having a total water volume of about 1,000 m 3 .
In Figure 2a, schematic representation of the water network is depicted.As shown, the pipeline system can be subdivided into two sub-network, supplied by the two water tanks.70 pipes and 57 nodes compose the network under investigation.In Appendix A (Tables A1 and A2), details on the distribution network (diameter, pipe length, water supplied, etc.) are reported.In Appendix A it is also reported the value of the roughness, clearly a key parameter, even though the calibration process (not shown here) pointed out that all the pipes fall into a hydraulically smooth regime, with values of equivalent sand grain roughness k s [34] in the range: 0.001 mm < k s < 0.005 mm.For this reason, the roughness k s was set equal to the average value (k s = 0.003 mm).
The water distribution network configuration is characterized by significant differences between the geodetic quote of water tanks and the users, suggesting high levels of energy recovery.The network orography thus seems to have an excellent configuration in order to test the application of micro energy recovery, introducing PATs in the network pipes.
One of the main problems in PATs positioning is to find the optimal installation point, thus to achieve both maximum energy production and desired water supply.Due to the small number of users, in the present case the choice of the optimal installation point has been carried out following the hydraulic regulation of the network.In small WDNs, the hydraulic regulation can be easily applied, achieving high performance levels, as demonstrated by De Marchis and Freni [28].
referred to TC2) analyzed the case of a PAT installed in the connection pipe, number 01, between the Castle tank and the network, having a diameter of 300 mm.In Table 1 some details on the BEPs of the scenarios, used in Equations ( 9) and ( 10) of the PAT model, are reported.Taking advantage of the high variation of the geodetic elevation between the water tanks and the network, the main objective of the present analysis is to recover energy, in terms of water head, from water distribution networks.
At the moment, in WDNs the energy is dissipated through float valves.As specified in the following, the cost of electrically regulated valves is similar to the cost of a pump working in turbine mode, with the fundamental difference that PATs have the capability of producing small amounts of renewable energy.Due to the specific case of study, three different scenarios have been explored.In the first PAT installation scenario (hereafter referred to TC1), the PAT has been installed in the pipe 23 having a DN equal to 200 mm and connecting the Cave tank with the WDN.The second scenario (hereafter referred to TC2) analyzed the case of a PAT installed in the connection pipe, number 01, between the Castle tank and the network, having a diameter of 300 mm.In Table 1 some details on the BEPs of the scenarios, used in Equations ( 9) and ( 10) of the PAT model, are reported.Finally, a third scenario (hereafter referred to TC3) has been simulated considering the simultaneous installation of both the turbines in the main pipes.The analyzed network, as briefly described above, is characterized by the presence of private tanks that strongly modify the network behavior.
The local reservoirs, installed by the users to cover the water scarcity condition periods, can affect the energy production and the efficiency of the PAT installation.In this framework, the proposed dynamic numerical model was applied to investigate on the effect of the user water tanks in the energy recovery from the WDS.
In the following, N-will indicate scenarios without the private tanks, whereas with the capital letter T-will specify the cases where user's tanks are considered.
In light of the year to year economic analysis, in the present case a single demand condition was simulated (see Appendix A for details), representative of an average condition of water supply.Even though the demand varies throughout a single day as well as throughout the seasons, the aims of the present research is to give indication about the energy recovery during several years, thus the use of an averaged condition of water demand is clearly justified.A daily demand pattern is request for diurnal analysis of energy recovery.
Considering that average demands are used and a year-to-year analysis is performed, the use of complex mathematical model able to simulate the transient condition could be not necessary.Nevertheless, the presence of private tanks could create small changes in water distribution networks.De Marchis et al. [25], in fact, demonstrated that small water hammer events occur in WDN due to the opening and closing of the ball valves supplying the private tanks.Even though, in the proposed analysis a steady state version of the mathematical model is applied, the presence of user's tanks suggest the use of a model, based on the MOC method, able to take into account for these small dynamic phenomenon, interacting with the PAT module.

Economic Analysis
The results of a generalized model to simulate both WDN behavior and power production are discussed with the aim of analyzing the energy recovery in water network of small towns.The analysis is based on the average daily energy production (ADEP (kWh/day)) obtained by integrating PAT power profile along 24 h, through the equation: where T = 86,400 s, η t is the PAT efficiency, γ is the specific weight of water, q is the flow rate through the supply pipe, and H t is the head drop due to the PAT, calculated according to equations of the characteristic curves achieved by Derakhshan and Nourbakhsh [15] (Equations ( 9) and ( 10)).Average yearly energy production (AYEP (kWh/year)) is calculated assuming that ADEP is generated daily, Energies 2016, 9, 877 8 of 15 also in consideration that the network supplies a mainly resident population and that an averaged demand pattern was considered as input for the hydraulic model.In Table 2 details on the installed power and related costs are reported.Future NCF are actualized to present time by applying a discount rate calculated as the weighted average cost of capital (WACC) of a potential investor.
In the present case, the investor would be a public utility company so that WACC is effectively equal to the current bank interest rate.Here, coherently with the period and the Italian National Bank, it is assumed to be 3%.
The Net Present Value (NPV) is defined by the equation: where the (NCF) j is the cash flow at year j related to the conditions under analysis, I 0 is the initial capital investment which is assumed to financially occur all in year 0, and i is the WACC.The PbT and IRR are respectively the time t and the interest rate i at which the following equation is true: whereas the PI is given by Equation ( 14): Finally, the ROI at year j, is defined as: The economic analysis is based on the use of the MOC model to solve the equation of water distribution networks, presented above.We analyzed two different scenarios related to two different operational condition of the network.In the first one, private tanks are simulated in the hydrodynamic model whereas, in the second, the users are directly connected with the WDN and the presence of private reservoirs is neglected.
The values obtained for NPV, PI, IRR, and ROI are reported on Table 3 calculated assuming the unit energy revenue, equal to 0.22 €/kWh (the current feed in tariff).

Indicator N-TC1 N-TC2 N-TC3 T-TC1 T-TC2 T-TC3
I 0 (€) 4900 9200 Table 3 reports the principal economic indicators for investments both in case of no-tanks (N-) and with tanks (T-).Results show that in case of No tanks, Scenario 1 offers the best option in terms of PbT and IRR, whilst, Scenario 3 offers a better option with respect to NPV.
However, among the No-Tank cases, Scenario 1 is to be considered the preferred one due to a favorable combination between low investment and high NPV.Despite the fact that NPV is higher in Scenario 3, the lower investment of Scenario 1, implies a significant lower risk and significant lower yearly working capital (operation and maintenance expenses).
Among the cases with tanks installed, Scenarios 1 and 3 offer largely similar returns both in terms of NPV and other indicators.
Figures 3 and 4 show results in graphic form for both the No tank and Tank cases.
Energies 2016, 9, 877 9 of 15 Figure 4 shows, in particular, how the N-TC3 case becomes more attractive than N-TC1 case only after the seventh year of operation, at which, the cumulated actualized cash flows becomes higher than the N-TC1 case.In addition, Figure 3 shows that the Cases T-TC1 and T-TC3 are very similar in terms of cumulated cash flow, throughout all the period and both considerably more attractive than T-TC2.
The proposed model can be considered an innovative tool, able to guide the water utilities to the optimal positioning of hydraulic machines for the energy saving in WDN.Similarly to the analysis of De Marchis and Freni [28], all the simulations were carried out ensuring the same level of user satisfaction in terms of water supply.Figure 4 shows, in particular, how the N-TC3 case becomes more attractive than N-TC1 case only after the seventh year of operation, at which, the cumulated actualized cash flows becomes higher than the N-TC1 case.In addition, Figure 3 shows that the Cases T-TC1 and T-TC3 are very similar in terms of cumulated cash flow, throughout all the period and both considerably more attractive than T-TC2.
The proposed model can be considered an innovative tool, able to guide the water utilities to the optimal positioning of hydraulic machines for the energy saving in WDN.Similarly to the analysis of De Marchis and Freni [28], all the simulations were carried out ensuring the same level of user satisfaction in terms of water supply.

Environmental Analysis
In addition to the economic indicators, the installation of PATs implies a significant savings of greenhouse gas (GHG) emissions.
According to the National Energy Authority, the current national energy mix entails a consumption of 0.49 kg equivalents of CO2 per each kWh produced [35][36][37][38]

Environmental Analysis
In addition to the economic indicators, the installation of PATs implies a significant savings of greenhouse gas (GHG) emissions.
According to the National Energy Authority, the current national energy mix entails a consumption of 0.49 kg equivalents of CO 2 per each kWh produced [35][36][37][38].The renewable energy produced yearly by PATs in each of the scenarios directly translates into savings of energy which would otherwise be produced by traditional sources, thus in a direct saving of CO 2 emission.Figure 5 shows the GHG savings for all the scenarios.The no-tank scenarios induce higher CO 2 saving: N-TC1 and N-TC3 are the most beneficial.
Energies 2016, 9, 877 10 of 15 produced yearly by PATs in each of the scenarios directly translates into savings of energy which would otherwise be produced by traditional sources, thus in a direct saving of CO2 emission.Figure 5 shows the GHG savings for all the scenarios.The no-tank scenarios induce higher CO2 saving: N-TC1 and N-TC3 are the most beneficial.
If we now couple the economic results described above with the benefits induced by savings of GHG, we may conclude that N-TC3 is the best scenario.In general terms, should the No-Tank case apply, TC1 is the preferred option from the economic point of view, and TC3 becomes the preferred one if we consider also the environmental benefits.Should instead the Tank case apply, TC3 proves the best option both in economic and environmental terms.Another environmental consideration worthy of attention is leak reduction.It is known that even though several efforts have been made to reduce hydraulic leakages in WDN, almost all the networks are characterized by water losses.The leakages are strictly connected with the pressure in the networks (see De Marchis et al. [39]) thus the use of PATs, able to produce renewable energy reducing the pressure, represent a double environmental challenge, worthy of extensive future investigations.

Conclusions
In the proposed study the capability of recover energy in small water distribution network is If we now couple the economic results described above with the benefits induced by savings of GHG, we may conclude that N-TC3 is the best scenario.In general terms, should the No-Tank case apply, TC1 is the preferred option from the economic point of view, and TC3 becomes the preferred one if we consider also the environmental benefits.Should instead the Tank case apply, TC3 proves the best option both in economic and environmental terms.
Another environmental consideration worthy of attention is leak reduction.It is known that even though several efforts have been made to reduce hydraulic leakages in WDN, almost all the networks are characterized by water losses.The leakages are strictly connected with the pressure in the networks (see De Marchis et al. [39]) thus the use of PATs, able to produce renewable energy reducing the pressure, represent a double environmental challenge, worthy of extensive future investigations.

Conclusions
In the proposed study the capability of recover energy in small water distribution network is explored.Basically, a mathematical model was applied to investigate the installation of hydraulic machines known as PATs in small districts of WDNs.The economic and environmental analysis was performed including, in the hydraulic simulation, the presence of private tanks, installed by the users.This way, the results can be considered representative of almost all the networks of the Mediterranean countries, where the practice to use small single reservoirs is well established.
The analysis showed that the use of PATs can be attractive in economic and energy recovery terms, considering the lower capital cost, when compared to traditional turbines.Nevertheless, the PAT installation point affects energy saving.The computed compound Payback Time achieved values of about 2 years or 12 years as a function of the PAT position, clearly showing that, in some cases, the investment is economically unfeasible.
The proposed model can be thus considered a tool, for the water utilities, to maximize the energy saving in WDN.It is useful to specify that all the simulations were carried out ensuring the same level of user satisfaction in terms of water supply.
Overall, the results confirmed that private tanks cause a high decrease of energy production.Considering that, in Mediterranean countries, water utilities have no control on these tanks, to give a more realistic analysis on the performance of energy recovery in WDS, private tanks must be considered.
Finally, the model indicates the possibility of a significant reduction of atmosphere emission of CO 2 equivalent greenhouse gases, thus confirming the utility to apply distributed micro energy generators along the several and several WDNs.

Figure 1 .
Figure 1.Installation scheme of a coupled system of pumps as turbine (PAT) and pressure reduction valve (PRV).

Figure 3 .
Figure 3. Net actualized cash flows generated in all scenarios.

Figure 4 .
Figure 4. Actualized cumulated cash flows generated by all scenarios.

Figure 4 .
Figure 4. Actualized cumulated cash flows generated by all scenarios.

Figure 5 .
Figure 5. Amount of CO2 saved (kg) yearly in all scenarios analyzed.

Figure 5 .
Figure 5. Amount of CO 2 saved (kg) yearly in all scenarios analyzed.

Table 1 .
PATs characteristic in turbine mode.

Table 1 .
PATs characteristic in turbine mode.

Table 2 .
Economic feasibility of PATs installation in a WDN: Production and Cost analysis for cases with no tank (N-) and with tank (T-).EGE: energy generation equipment; CWC: civil work cost; and MC: maintain cost.
. The renewable energy