Applied Strategy to Characterize the Energy Improvement Using PATs in a Water Supply System

Sustainable development has been an idea raised in recent years. The results are related to the improvement and the use of new technologies to maximize efficiency in water management. However, energy consumption has been increasing as a consequence of new management and uses of water. Especially in pressurized water distribution systems, the use of pressure reduction valves (PRVs) increases the water usage efficiency but it decreases the energy consumption efficiency, since the valves dissipate energy that could be recovered. This research presents a proposal of a recovery system based on the installation of pumps used as turbines (PATs). These machines are located in different points of the high-pressure water distribution system in the Valencia Metropolitan System (Spain). An annual estimate of the theoretical recoverable energy as well as the “ideal” pump for each point were proposed. The theoretical recovered energy value was 847,301 kWh/year for a specific analyzed point. Besides, the characteristic curves of the PATs from a selected point were determined, estimating an improvement in the sustainable indexes. The calculus of these green parameters showed that the implementation of this solution caused a reduction in consumed energy of 1.50 kWh/m3.


Introduction
Water supply distribution systems have a fundamental objective: assuring the availability of water in the desired quantity and the required quality to satisfy the health parameter. To fulfill this intent, urban distribution systems depend on a great amount of energy to supply the water at a minimum pressure for it to be rightly used.
Energy has always been considered an important factor in society. Its significance has grown pressingly in recent years. This is due to the uncontrolled increase in the world population, the concentrations in urban cities and the rise in the demand of resources such as water and energy which threaten to force a future global deficit in natural resources [1].
In the water sector, the water-energy nexus has constantly been an important issue, since one of the most energy-consuming activities in a city is the cycle of water consumption. However, water supply and distribution are services which must be provided to the population as an obligation, since it is a right that concerns everyone. Therefore, water managers have continuously tended to find ways to enhance their systems. These strategies are focused on looking for approaches to have economic and sustainable development [2].
Water distribution systems consume different amounts energy around the world. This consumption depends on the topography and topology of the water network mainly. As an example, the water cycle represented around 7% of the total energy consumption in Spain in 2010. A total of 20% of the energy used for the water cycle was employed in the distribution of water to the consumption points in the different sectors [2]. Particularly in the city of Valencia (Spain), the energy consumption was 53. 21 GWh/year when the water cycle was considered in 2009. This value can be normalized to 1124 kWh/m 3 when the water consumption was considered. The catchment, potabilization and distribution of water represented 30% of the energy used for the water cycle in Valencia [3].
Water distribution serves as the one phase of the water cycle, which consumes the most energy. Therefore, this phase is the broadest to improve the efficiency indexes. Undoubtedly, one of the easiest methods to reduce the energy required for water distribution is pressure management. If the pressure is low in the system, it has less leaks, thus the flow decreases as well as the energy input. There are three basic methods to achieve pressure management. These operational strategies are described in Table 1. Table 1. Advantages and disadvantages of different pressure management methods.

Method Advantages Disadvantages References
Pressure reducing valves (PRVs) • Cheapest initial investment.

•
Decreases water leaks in the network.

•
Manages the pressure downstream of the valve.

•
Relatively simple design and operation.

•
The energy dissipates but is not recovered.

•
Flawed use if the energy is introduced into the system. [4,5] Pressure management zones (PMZs) • Ease of pressure management due to isolation of topographically different zones. • Better use of energy introduced into the system. •

•
Decreases water leaks in the network.

•
Manages the pressure downstream of the PAT.

•
Relatively complex design and operation. • Need to install PRVs to regulate remaining overpressure. [9][10][11][12] Water distribution systems mostly depend on energy [13]. This connection between water and energy creates a nexus, being part of the near-future sustainable development main consideration. In this framework, the use of small and micro-hydropower systems (specially using pumps working as turbines-PAT) has been increasing in recent years. PATs have demonstrated to be a feasible solution to recover the potential energy in water systems, both urban [14,15] and irrigation [16]. This pump's feasibility is linked to economic, availability and construction criteria. Theoretically, all types of pumps could be used in water systems [17,18]. However, when the water managers choose a PAT, it is advisable to avoid choosing axial pumps since these are less efficient than radial machines [19]. In addition, the axial pumps are the machines that have had fewer experiments in this regard and therefore, the efficiency prediction is most difficult [13].
This research establishes a strategy to choose PATs according to the methodology proposed by [20]. The novelty of this manuscript is that it was applied using data from pairs (i.e., flow and pressure) of the real water systems. Therefore, the proposed energy analysis, the selection of the machines as well as the comparison between different scenarios were applied to reality. This research aims to estimate the potentially recoverable energy of the supply network at the chosen study point as well as the selection of the machines to be used in those points. This selection was based on a manufacturer catalogue of pumps. For a specific case study point, the complete hydraulic design was made in order to approximate, in a more precise way, the recoverable energy in addition to the calculation of the proposed sustainability indicators to compare the impact the installation of PATs has on the system.

Materials
This manuscript aims to analyze several chosen points from the high-pressure water distribution system, which is called the Valencia Metropolitan System (VMS) in Spain. Those points have PRVs operating at the moment to try and maintain a steady pressure throughout the network. The pressure values upstream and downstream of the valves at each point are registered with an interval of 10 min, alongside with the date and the flow. Using these data available, one can only calculate the theoretical recoverable energy (ETR) and the theoretical unrecoverable energy (ENTR). The available data consist of more than 50,000 entries for variables of date, flow (Q), upstream (Pu) and downstream pressure (Pd) for each of the chosen study points. They are the results of the data registered from a PRV installed at every case, with a time interval between each entry of 10 min. The used data considered one year from January to December. Really, the operational conditions are similar each year because the water system case-studied is the main distribution ring. This network distributes the different origin resources to the different municipalities where there are other tanks and pipes, which supply to users. Table 2 exemplifies the data available for each study point. The meaning of every column of data is: • Time (∆t = 10min) is the date when the data were registered, they are in dd-mm-yy hh:mm format with 10 min of interval; • Q is the instant flow registered at the date, measured in liters per second (L/s); • is the pressure upstream of the valve registered at the date, measured in meters of water column (m w.c.); • is the pressure downstream of the valve registered at each time (m w.c.). This is the setting value, which should be established to guarantee the pressure conditions in the water supply system. The valves are chosen to keep the pressure head close to the setting.  Figure 1 describes the proposed methodology, which was followed in this research. It is focused on the selection of the PAT for each analyzed point, highlighting two important analyses to be made: the energy balance and the estimation of the theoretical recoverable energy. Besides, the plan establishes the operation strategy when the system operates considering hydraulic regulation. The procedure is integrated on three main blocks: A. energy balance, B. PATs selection and C. theoretical recoverable energy estimation and analysis of energy indexes. These blocks are described in the next subsections.

Energy Balance (A) for the Supply Network at the Chosen Study Points
The first block ( Figure 1) is an energy balance (Step I). It aims to determine the total energy supplied to the system while differentiating the required energy for the proper water consumption from the energy loss, which could be either operational losses or structural losses [21].
The proposed energy balance to quantify the theoretical recoverable energy is based on the energy distribution plan shown in Figure 2. The difference between piezometric lines can be observed in Figure 2a  The energies that define the scheme presented below are described as follows [16]:  The values for each of the previous energies described are defined by the expression presented in Table 3. These terms were defined by [5]. Table 3. Calculation expression for the different types of energy in a water distribution system.

Energy Calculation Expression
Total energy (ET) E T (kWh) = 9.81 3600 Friction energy (EFR) E FRi (kWh) = 9.81 3600 Q i (z 0 − (z i + P i )) ∆t (2) Theoretical energy necessary (ETN) E TNi (kWh) = 9.81 3600 Q i P min i ∆t Energy required for consumption (ERI) E RIi (kWh) = 9.81 3600 Theoretical available energy (ETA) E TAi (kWh) = 9.81 3600 Theoretical recoverable energy (ETR) E TRi (kWh) = 9.81 3600 Theoretical unrecoverable energy (ENTR) where subscript "i" indicates each analyzed point; Qi is the flow circulating by a line (m 3 /s); and zi is the geometry level above the reference plane of the analyzed point. In this case, the reference is sea level (m); z0 is the geometry level above the reference plane of the free water surface of the reservoir. In this case as well, the reference is sea level (m); ∆t is the time interval (s); Pi is the service pressure in the analyzed point when consumption exists (m w.c.); Pmin i is the minimum pressure of service of an analyzed point to ensure the minimum in the most disadvantageous consumption node (m w.c.); Pmin C is the minimum pressure of service of an analyzed point required by the local authorities, this value is set to be 25 m w.c. in VMS; and Hi is the head in the analyzed point (m w.c.), obtained as = − max( min ; min ). This information is used to calculate the total instant power supplied to the analyzed point defined by the following equation: where Qi is the flow circulating by the valve in L/s; and η is the efficiency of the turbine. It is a value between 0 and 1. For estimation purposes, an efficiency of 1 is selected; Pu is the pressure upstream of the PRVs. It is measured in meters of water column (m w.c.). In a similar calculus, the theoretical recoverable instant power is estimated for the analyzed point with the next equation: where η is the efficiency of the turbine. It is a value between 0 and 1. For estimation purposes and comparison, efficiencies of 1 and 0.6 have been selected; Hi is the dissipated pressure by the PRVs. It is obtained by the difference between and measured in m w.c. When is known, the total energy supplied to the analyzed point is calculated using the next equation: where ∆t is the time interval (s). Subsequently, the theoretical recoverable energy is estimated by the following equation: Having determined the previous values, the theoretical unrecoverable energy is estimated. The value is defined by the following equation: To ease the analysis, and for better understanding of the data, the results are grouped in time intervals of 1 h, 24 h, 1 month and 1 year.

Selection of PATs for Each Study Case (B)
The implementation of PATs in a water distribution system is widely justified by many advantages such as the ability to recover energy from the network, the decrease in the carbon footprint, the lower investment cost (compared to the installation of a traditional turbine), the readily available pumps and spare parts and the simplicity of the system [13].
However, the use of this technology has been depleted by the lack of characteristic curves for the pumps working as turbines. The absence of these curves creates an uncertainty in the design of the energy recovery systems. Although the right methodology is to ask the manufacturers for the curves of the machine working in reverse or to carry out experimentation in a laboratory to obtain these curves, this is not common because knowledge of curves is difficult to achieve. Usually, the design is based on empirical methods developed through the year. The following figure illustrates the methodology used in this paper: The block "B" is divided on four main course procedures derived from Figure 1. These are: (B-1) determinate the theoretical best efficiency point (QBEP,T, HBEP,T) for the ideal PAT at an analyzed point (steps II and III); (B-2) select a standard water pump by the catalog data that matches the requirement for energy recovery at the analyzed point (steps IV to VII); (B-3) estimate the operation point of the machine in turbine mode, based on the manufacture's data in pump mode (steps VIII to X); and (B-4) approximate the characteristic curves of the selected PAT by the proposed empirical curves (steps XI and XII).

B-1. Determination of the theoretical best efficiency point
The analysis of the operation range depends on a detailed statistical study. For each of the eight different study points, there are more than 50,000 entries for every variable (date, flow, upstream pressure and downstream pressure). A high variability for the values is expected, due to the great amount of data available and the time length of the series. As a result, the best efficiency point is the most representative operation point of the data series. To determine the most typical flow and the dissipated pressure in the analyzed point, a statistical study was based on the use of histograms for each variable (flow and head).
A fruitful technique, which is based on programming, was developed in this manuscript in order to be used to enhance the precision of the selection of the most characteristic variable. This was to delimitate the data to the desired work zone and repeat the process of the histogram and scatter plots.

B-2. Selection of a standard water pump by catalog data
The PAT selection is based on the methodology, which is proposed in Figure 1. The steps are the following: (IV) calculate the specific speed for the "ideal" turbine for the analyzed point; (V) determination of the coefficients to adapt the previous operational point as a turbine to a best efficiency point as a pump; (VI) estimate the best efficiency point (BEP) for the required machine in pump mode, using the previous coefficients; and (VII) select the standard pump from the available machines in the market based on the estimated BEP.
The first three steps mentioned above are to estimate the operational point for the standard water pump that matches the requirement previously established for a specific analyzed point. The head to be recovered by a machine in turbine mode does not match with the head for the same machine in pump mode [13]. The estimation of the behavior of pumps in their turbine mode has been an important issue and has led to numerous studies. Unfortunately, most of the studies did not have enough pumps for experimentation and therefore their results lack a certain precision [20]. The proposed strategy uses both Pérez-Sánchez et al.'s method and other proposed methods for comparison and validation.

B-3. Estimate the operation point of the machine in turbine mode
Three steps are considered in this sub-phase according to Figure 1: (VIII) determine the BEP from the catalog of the selected pump, estimating the conversion factors to transform the BEP in pump mode to a BEP in turbine mode; (IX) estimate the best efficiency point for the machine as a turbine by the previous coefficients; (X) calculate the approximation of C value. A valid selection must have a value of C ≤ 1.
The characteristics curves of a machine can be defined by polynomial equations depending on the specific speed (ns) of the machine. When there is only a suction pipe, this value is used to classify either pumps ( ) or pumps in turbine mode ( ). These values are determined by the following expressions [20]: where QBEP,P is the flow in the best efficiency point in m 3 /s when the machine operates as a pump, HBEP,P is the head in the best efficiency point in m w.c. when the machine operates as a pump, QBEP,T is the flow in the best efficiency point in m 3 /s when the machine operates as a turbine, HBEP,T is the head in the best efficiency point in m w.c. when the machine operates as a turbine and n is the rotational speed in rpm. Previous research [20] proposed the following empirical relationship to define the best efficiency point of the machine as a function of the specific speed (Table 4): Table 4. Empirical relationship to define the best efficiency point of a machine as a function of the specific speed proposed by [20]. In Table 4, the variable , is the best efficiency obtained from the machine working as a pump. These values are used to estimate a flow, head or efficiency of the machine working as a pump or as a turbine, based on the same variables when the machine is working in the opposite way. They are applied in the following equations: where βQ is the conversion coefficient for the flow, βH is the conversion coefficient for the head and βη is the conversion coefficient for efficiency. This research also includes the empirical equations proposed in previous studies to compare errors. These formulas are summarized in Table 5: The criteria to accept the pump selected is based on calculating an approximation value between the estimated best efficiency point of the machine in turbine mode and the theoretical best efficiency point determined for the analyzed point. The concept is based in an ellipse form which enables an error varying from ±30% at the major axis and ±10 at the minor axis for both the flow and head. The approximation value must be C ≤ 1, and the value is defined by the following equation [20]:

B-4. Estimation of characteristic curve of PATs from empirical curves
This last sub-phase related to PAT selection considers the steps XI and XII from Figure 1. The goal of this sub-phase is to estimate the characteristic curves of PATs, [20] proposed curves as function of the nondimensional numbers of discharge (φ) and the head (ψ) defined by the following equations.
where D is the impeller diameter in m, g is the gravity constant in m/s 2 , Pt is the shaft power in the PAT, ηt is the efficiency of the PAT and n is the rotational speed in rps. The empirical curves mentioned as a function of the nondimensional numbers of discharge and head, as well as the specific speed of the test pumps, returns estimate curves of Qt-Ht and Qt-ηt. They are presented in Figures 3a and 3b.
where HRunaway is the head in m w.c. at which the machine runs with a power output of zero; kRunaway is the coefficient proposed by Pérez-Sánchez et al.; Hzero-speed is the head in m w.c. while the rotational speed is zero; and kzero-speed is the coefficient proposed by Pérez-Sánchez et al. (2020).

Estimation of Theoretical Recoverable Energy and Calculation of Sustainability Indicators
Once the characteristic curves from the chosen PATs were known for an analyzed point, a more precise theoretical recoverable energy estimation was made. The characteristics curves were meant to be used to evaluate the hydraulic efficiency of the recovery system. However, this research considers the electrical efficiency of the generator for a more realistic approximation since the electric generator efficiency has a great significance on the global efficiency of the machine when hydraulic and electrical machines are considered [31]. The efficiency used for the electric generator represents the lowest permitted efficiencies in the marked motors associated to the IE code. They are defined by IEC/EN 60034-30-1 [32]. It has four different classes of motor: IE4 Super-Premium efficiency, IE3 Premium efficiency, IE2 High efficiency and IE1 Standard efficiency. This manuscript chose the IE1 motor, which varies as a function of the power. A common practice to increase the recovered energy is to install a variable frequency drive for the machines. However, for this estimation, a variable frequency drive was not considered for the installed PATs as the worst-case scenario is being analyzed.
An energy balance is repeated after the implementation of the energy recovery system to compare the existing conditions to the post-solution condition. The sustainable indicators shown in Table 6 are meant to reveal the changes connected to the installation of the PATs and it helps to assess the performance of the water distribution system against a sustainable focus [33]. This research proposes the sustainability indicators for pressurized water systems recommended by [34]. These indicators emphasize the focus on the energetic and environmental development for these systems. Two other indicators were proposed. They are based on the previously mentioned indicators. They were the absolute annual consumed energy (IAAE) and the absolute consumed energy per unit volume (IAEFW). These indicators present better the reduction in needed energy in a system where the recovered energy is being consumed.

Case Study Description
This research featured eight different network points from the high-pressure water distribution system of VMS in Spain. The available data were an annual series related to flow, upstream pressure and downstream pressure from an installed pressure reducing valve at each point, referenced to a piece of date information with an interval of 10 min.
The water distribution network of VMS is schemed in Figure 4, highlighting the study points of the system. The points between 1 and 5 are presented in Figure 4a and the remaining points are presented in Figure 4b. The legend from the figure specifies the color code to differentiate each one of them through the continuing length of the study.  Figure 5 summarizes the irregularity of the data. A high variability can be observed. This is because there are more than 50,000 entries for each variable and the water system operated under demand:

Theoretical Recoverable Energy Estimation
Based on the mentioned methods, the theoretical recoverable energy is estimated for two different scenarios: i.
Ideal efficiency (η = 100%): an ideal and practically unreachable scenario, it is used as a reference to compare the ETR and the approximation made after the PATs are selected; ii.
Conservative efficiency (η = 60%): a conservative scenario, using the right PATs selection and the right hydraulic design, this value is likely to be obtained or surpassed.
To ease the comparison, the results of the theoretical recoverable energy estimation were grouped at bigger intervals such as 1 h, 24 h, 1 month and the complete data series of a year length. The results for point #7 are summarized in Tables 7 and 8:  As shown in the previous tables, a theoretical recoverable energy of 847,301.47 kWh/year was available at point #7, representing 57.35% of the total energy consumed at this point. For a more conservative approach, considering an efficiency of 60%, the theoretical recoverable energy was 508,380.88 kWh/year which represented 34.41% of the total energy consumed at this point. From the difference between the average theoretical recoverable energy monthly and both, the minimum and maximum theoretical recoverable energies were −10.79% and 16.40%, respectively. Therefore, the conditions were fairly stable for the system, facilitating the selection for a PAT.
The same estimations were made for each study. Figure 6 compares the percentage of theoretical recoverable and unrecoverable energy for each point:  Figure 6 highlights that points such as #1, #2 and #5 consumed more than 1,000,000 kWh/year but their theoretical recoverable energy was at maximum equal to 50%. However, point #2 had the largest amount of theoretical recoverable energy. In contrast, point #3 presented the highest percentage of ETR (above 70%) for a total of 1,055,420 kWh/year.

Most Representative Operational Point Estimation
For each point, a preliminary statistical analysis was made. The results of the analysis for point #7 are presented in Figure 7d. Following the proposed methodology, Figure 7a,b shows the flow and head histograms, respectively. When point #7 was analyzed, the registered flow was concentrated in two main ranges: (i) 45-60 L/s and (ii) 390-405 L/s (Figure 7a). Both ranges tended to operate for different heads, which are shown in Figure 7c. Therefore, if the recovery systems want to be defined for this point, two PATs groups were selected with one of the groups consisting of two PATs working in parallel.
A representative operational point was determined for the PATs to be installed in the system for each analyzed point. Some points would have various operational points as the hydraulic design contemplates different PATs groups to optimize the recoverable energy available for the specific point. The used criterion to choose the PAT considered the frequency histogram of the flow and head for each machine. Different iterations were checked using programming to choose the best point. The results are summarized in Table 9.  Figure 8 highlights the graphical classification for the PATs of each point. This was based on the proposed general range of application of different turbine types. The standard pump was selected from the catalog data available in the market, looking for the type of machine specified in the previously mentioned image.

PATs Selection
Point #7 had two PATs groups with different characteristics, the conversion for both operational points will be presented in the same table. Table 10 offers a summary of the operational point for either of the groups. Based on those points, the best efficiency point was estimated comparing different methods. The results are presented in Table 11. Using the prior estimated BEP, the selection of standard pumps that matched these requirements was made. The selected pumps, which were chosen using manufacturer catalogs [33,34], and their best efficiency points are shown in Table 12.  Having selected the pumps for the groups of point #7, the transformation from BEP as a standard pump to a BEP in turbine mode was made. The results of the coefficients calculation and the BEP for each group are presented in the Tables 13 and 14. Both PATs meet the selection criteria for every method of estimation, except for the Alatorre-Frenk one, where the C value was 1.38. However, the main method used for this manuscript was the proposed method by [20]. This approximation presented good C values, which were 0.14 and 0.07 for group 1 and group 2, respectively. The same methodology was performed for each case of study. Table 15 presents the selected pumps for each analyzed point along with their best efficiency point as a standard pump. In addition, it shows the BEP estimation in turbine mode once the pump was selected. The pumps were chosen using manufacturer catalogs [34][35][36][37].

Estimation of PATs Characteristic Curves by Empirical Method
As previously mentioned, the following sections depend solely on the data from point #7, although they can be applied to the rest of the points. The empirical methodology proposed by Pérez-Sánchez et al. is based on three variables: (1) non-dimensional discharge number (φ); (2) nondimensional head number (ψ); and (3) specific speed for the pump used as a turbine ( ). These variables for the selected machine are presented in Table 16. Using the earlier values, the empirical curves were selected to estimate the values for the characteristic curves of each PAT (Figure 3). The curves, which were defined by the specific speed of nsT = 38.49 and nsT = 40.85, were used to interpolate the required curve defined by = 39.47. Based on these curves, the characteristic curves for the PATs from point #7 were obtained. Figure  9 shows the estimation of the curve for both selected machines at point #7. The operational curve is defined considering the machines operate a uniform rotational speed. The minimum and maximum operational flows are defined in the figure 9. The characteristic curves for group 1 are defined by the following equations: • The head and efficiency for flows lower than 60 L/s are obtained by these equations: For group 2, the characteristic curves shown in Figure 9 are defined by the following equations:

Estimation of the Theoretical Recoverable Energy with the Installed PATs
Once the characteristic curves for each PAT installed in the analyzed point were defined, the next step was to estimate the theoretical recoverable energy with the newly selected PATs installed. The theoretical recoverable energy was estimated considering the PAT efficiency curve as a function of the flow (Figure 9). In Figure 10a, a schematic arrangement is presented for the installed pump's important elements such as: the installed PATs, the PRVs for the bypass and the final regulation. The pressure reduction valves are that important for the system because they are responsible for the dissipation of the remaining energy in the system. The installation of these valves allows water managers to guarantee the flow and the pressure at the consumption points. The flow variation over time does not have significance if the machine is regulated correctly (when there is variation in the rotational speed) and the pressure reduction valves are defined, installed and controlled by pressure sensors. Another important factor for the estimation is the electric generator. Two scenarios were considered for this matter: (i) variable efficiency as a function of the outpower and (ii) constant efficiency of 85%. The results for the estimation in each scenario are presented in Figure 10b. Figure  10c shows the energy distribution for every month of the analyzed data series assuming an electric generator with variable efficiency. Lastly, Figure 11b presents the energy recovered by each machine alongside the energy dissipated by the valves.
To complete the hydraulic regulation for this energy recovery system, a suitable pressure reducing valve was to be selected by using a manufacturer catalog [38] and defining its open degree. The variable to choose the valve was the flow coefficient ( ) defined by the following equation: where Qt is the flow of the PRV in m 3 /h, ∆Pt is the dissipated head at the valve in kp/cm 2 and ρ is the relative density of the fluid (1 for the water). The selected valves are presented in Table 17. For each valve, the risk of cavitation was calculated and, based on the information supplied by the manufactures, this was not an issue.  The last step is to calculate the sustainability indicator proposed in Table 6. The results are presented in the Table 18: The installation of PATs in the water distribution systems represented a decrease of 64.16% in the dissipated energy, obtaining an IED equal to 0.21. As Table 18 states, the annual consumed energy (IAE) and the consumed energy per unit volume (IEFW) stayed the same after the PATs installation as these values depend on the total energy consumed in the analyzed point, a value that it not going to change. This manuscript proposed two new indicators, the absolute annual consumed energy (IAAE) and the absolute consumed energy per unit volume (IAEFW). Both of them presented a reduction of 36.91%, resulting in values of 932,028 kWh/year and 2.56 kWh/m 3 , respectively. Due to the increase in recovered energy in the system, the IER value was 545,335.26 kWh/year, representing an ERP of 36.91%. These changes are included in the IEC, presenting a reduction of 35.01% for a new value of 0.023 EUR/m 3 . The feasibility of these installations was demonstrated in different research in which the payback periods were lower than six years [9,16,18], including this feasibility index in the optimization procedure. Related to investment at a preliminary level, an average value around 545 EUR/kW is necessary to establish the recovery systems, although the final cost should be considered in a specific project [9,16].

Conclusions
This research presented a deep analysis of the available data from eight different points from the high-pressure water distribution system from the Valencia Metropolitan System in Spain. Three main focuses were considered: the theoretical recoverable energy available in each analyzed point; the selection of a suitable standard pump to be used as a turbine in each case; and an estimation of the characteristic curves for the selected PATs for a specific point alongside the analysis of sustainability indicators to evaluate the impact of the installation of the machines on the system. The manuscript proposed a strategy joining different published methods. This strategy shows as a novelty the possibility to define the hydraulic regulation in a real system using real pumps and valves in the selection. Therefore, the analysis includes both theoretical as well as real approximations using real machines in the selection. The approximation of the "C" values for the selected pumps for each of the analyzed pumps were below 1, with an average value of 0.297. This supports the theory that there are available pumps to be used as turbines in the market. It is a step ahead in the proposal of the use of PATs in real systems when the consumed volume is high in the network.
Finally, assuming an electric generator with an efficiency defined by the motor IE1, the energy recovery system for point #7, theoretically, recuperates a total of 545,335 kWh/year or 36.91% of the total consumed energy at the point. This value was around 64% of the theoretical recoverable energy ( Figure 6) at point #7. This recovery results in a 64.16% decrease in the dissipated energy in the system and a 36.91% decrease in the IAEE and the IAEFW. Based on the Spaniard electrical rate for an installation of more than 15 kW of power capacity, the system has an IEC reduction of 35.01%, improving the water cost from 0.036 EUR/m 3 to 0.023 EUR/m 3 . In this sense, the efforts of the water supply network managers of the VMS to analyze the recovery potential in the water system are remarkable. In the future, cities must be resilient and adapt to meet sustainable development goals, which imply maximum efficiency in all aspects related to the water-energy nexus.