Optimal Pump Scheduling for Urban Drainage under Variable Flow Conditions

The paper is focused on the optimal scheduling of a drainage pumping station, complying with variations in the pump rotational speed and a recurrent pattern for the inflow discharge. The paper is structured in several consecutive steps. In the first step, the experimental set-up is described and results of calibration tests on different pumping machines are presented to obtain equations linking significant variables (discharge, head, power, efficiency). Then, those equations are utilized to build a mixed-integer optimization model able to find the scheduling solution that minimizes required pumping energy. The model is solved with respect to a case study referred to a urban drainage system in Naples (Italy) and optimization results are analysed to provide insights on the algorithm computational performance and on the influence of pumping machine characteristics on the overall efficiency savings. With reference to the simulated scenarios, an average value of 32% energy can be saved with an optimized control. Its actual value depends on the hydraulic characteristics of the system.


Introduction
Nowadays, the issue of management efficiency is of great concern in the context of water infrastructures.Focus is usually given to the reduction of cost of production [1], energy recovery [2], reduction of water losses [3], reduction of pipe breaks and maintenance operations [4].Cost of production especially includes cost of the energy used for pumping, which, according to recent reports [5,6], constitutes 4% of the entire amount of national electricity consumed in U.S. [4,7] and 7% of the electrical energy worldwide [8,9].Reduction of pumping energy use within water networks is one of the most promising fields in the context of energy recovery and efficiency [10][11][12].
Pumping systems within water supply and drainage networks are equipped with multiple pumps, starting with a minimum number of two, one of them is kept for replacement purposes.Pumping systems has been commonly designed to work at fixed speed and constant hydraulic conditions (head and discharge) which are close to the Best Efficiency Point (BEP) of the pump, so to have the best possible performances.Given the presence of multiple pumps in the system, and possible variations in the operating conditions (variable discharges and variable tank level to make an example), a scheduling is needed to optimize system performances.A modern trend in the management of pumping systems is based on the use of variable speed drives to change the impeller rotational speed of one or more pumps.This approach works well in pipelines with large friction losses and the energy use will reduce if the pumps works for a longer time at lower speed and with smaller discharges [13] when compared to fixed speed operations.Several examples of scheduling optimization exist in literature, with different optimization algorithms, different variables and different objective functions; the evolution of research in the field of energy efficiency optimization complies with the development of more and more sophisticated optimization tools and algorithms.A first attempt of reducing operational costs in water networks concerns the use of linear programming [14], integer linear programming [15], non-linear programming [16,17] and dynamic programming [18], with limited possibilities of generalizing results to any water network different from those tested.More recently, heuristic algorithms, such as genetic algorithms, ant-colony or harmony search [19][20][21][22], were applied in the coupling with hydraulic simulators, that were often overcome by using Artificial Neural Networks to reproduce results of the hydraulic simulations [23].The use of ANNs improved computational times, making it possible to use these tools for real time control of pumping systems [24][25][26][27].Some authors [28][29][30] developed an energy efficiency model for pump scheduling with multiple objective functions.One last energy efficiency method is the data-mining approach, where neural network algorithms are applied to develop models for predicting the pump energy consumption and the flow rate after pumping, based on previously collected data [9,31].
In the cited literature, pump scheduling optimization output usually consists of a set of values describing whether each pump is working or not at a given time (ON/OFF scheduling) [4,32] and which its rotational speed is [33], usually with one-hour pace [9].For wastewater pumping systems (both in combined and separated drainage networks), the main control is given by the water level in a storage tank [34] upstream of the pumps; however, any efficient method for energy consumption reduction should reproduce a control system which is able to deal with the highly complex, non-linear behaviour of stormwater flow rates [11].This complexity increases if the pump speed is considered in scheduling optimization.

Aim of the Paper and Methodology
Variable speed control and pump scheduling has been applied to clear water pumping stations in order to reduce the energy use so far.This work focuses on the optimal scheduling of a drainage pumping station.This problem has not been deeply investigated in the technical literature yet, due to the complexity of the mathematical formulation.Generally, these kinds of utilities are equipped with one or more sewage pumps working at constant speed (CS), whose activation is controlled by the water level within the wet well.When the pumping station is equipped with only one pump, if the water level exceeds an assigned value, the pump turns on; conversely, as the water surface decreases below a minimum value, the pump turns off.Some pumping stations are equipped with more than one pump: in those cases there are usually different set points to control each of the pumps.
The main scopes of this study are: • Formulating a complete mathematical model of a drainage pumping system • Optimizing the pumping system aiming at reducing the required energy • Showing the benefit of such an optimization comparing the results with a classical CS plant.
To achieve such objectives, the water level set points are removed and a new control technique is presented, where both the instantaneous rotational speed of the pump and the ON/OFF sequence are optimized.Such a flexibility in the control system makes the mathematical model quite complex, since introducing integer variables and non linear equations.Thus, a mixed-integer non-linear optimization model has been developed and a deterministic algorithm, i.e., a branch-and-bound non-linear optimization algorithm, has been selected to solve it.To the author's knowledge, only a few studies on this topic have been presented so far [11,35], dealing with either simplified models or heuristic approaches.
A real flow pattern of a drainage system has been selected as a case study and different scenarios have been simulated to investigate the performance of the optimization model.The results show that the energy savings obtainable by the scheduling optimization depend on the characteristics of the hydraulic system: the comparison between the CS control and the optimized control demonstrate that an average value of 32% energy can be saved, with a peak value of more than 70% only in the best case.The presented mathematical model and the encouraging results achieved suggest the possibility to develop real-time algorithms to control drainage and sewer pumps to reduce their energetic impact.

Classical Design
A pumping station is designed to collect and transport water to a point of higher elevation.The inflow water is stored in a wet well, where the submersible pump is located.In larger plants, surface pumps could be alternatively installed into dry wells and connected to the tank with a suction pipe.The pump lifts the water through a pressurized pipe toward another collecting system (i.e., a tank or a channel).Generally, the pumps work at constant rotational speed and fairly constant discharge, while the inflow is time dependent.If the grid frequency is not modified (i.e., 50 Hz in EU and 60 Hz in USA) the motor velocity (N) attains its maximum value N max .The difference between the inflow and the outflow causes modifications to the water level within the wet well.The classical control of the constant speed (CS) pump is based on the water level: when it rises above a certain level the pump is started.The pump is selected so that its discharge is higher than the maximum inflow: in this way, as the pump is started, the water level decreases.As the water level reaches a minimum value, the pump is turned off and the cycle starts again.A sketch of the wet well is shown in Figure 1.In the larger pumping stations, two or more pumps are placed to work in parallel, and different water levels are set to start and stop each of the parallel pumps.Since each pump operates in constant conditions, it is selected so that, at the selected discharge (Q), it exhibits the maximum efficiency.In other words, the designed outflow corresponds to the Best Efficiency Point (BEP) of the pump (Q BEP N max ).The time interval between two consecutive starts of the pump, τ, with Q in inflow discharge, is then equal to: where W is the volume of water stored between the maximum and the minimum water levels within the well.The number of starts per hour, S h can be calculated as: where τ is expressed in seconds.If the inflow discharge Q in is considered as a θ portion of the Q BEP N max , then the number of starts per hour can be written as: where the volume W is expressed in m 3 , Q BEP N max in m 3 /s and S h in starts per hour.For a wet well with a single pump, in order to avoid that the number of starts per hour (S h ) exceeds the maximum value indicated by the pump manufacturer (S h max ), the value of θ must be set to 0.5, as it can be demonstrated by deriving Equation (3) with respect to θ.Then, the volume can be calculated as: This formula, that relates the volume of the wet well only to the characteristics of the chosen pump, ensures that the number of starts per hour will never exceed the maximum allowed value, whatever the inflow discharge.The maximum number of starts occurs when the inflow discharge is the half of Q BEP N max and reduces otherwise.The choice of the pump in the design is crucial: a larger pump allows dealing with larger inflow discharge, but implies an increase of the volume, with a subsequent increase of the building costs and of the retention time.The latter can be calculated as the ratio between the volume of the well and the average input discharge [36].The retention time should be kept short to avoid septicity in the tank if it collects sewage or wastewater.The effect of the pump choice is discussed hereafter.

Experimental Investigation of the Behaviour of Two Submersible Pumps
Two submersible sewage pumps have been tested in the HELAB, the HydroEnergy Laboratory of the University of Naples (CeSMA).The laboratory has been specifically designed for the experimental tests of turbomachines according to ISO 9906 regulation.A large underground storage reservoir (110 m 3 of water) is used to feed the pump.The pumps were located on the bottom of the tank and connected to a 250 mm pipe that conveyed the water through an electro-magnetic discharge meter (accuracy < 0.5%) and then back to the tank.The pressure was measured with a piezoelectric transducer (accuracy < 0.5%).A digital power meter (accuracy < 0.5%) was used to measure the input power.A SCADA system ensured the automatic and simultaneous acquisition of the data.A speed-driver was used to modify the electric frequency ( f ) and to regulate the speed of the asynchronous motor.For each of the six tested frequencies (50, 45, 40, 35, 30 and 25 Hz respectively), a gate valve was used to set different discharge values.For each discharge, all the measurements were averaged on 500 samples (with a 0.01 s sampling rate) to reduce the fluctuation of the signal.For each of the two pumps, the maximum number of starts per hour, namely S h max suggested by the manufacturer is 10.Thus, a complete set of experimental points of discharge (Q), head (H), power (P) has been obtained.The rotational speed (N, in rpm) of the pump can be calculated (if the slip of the motor is neglected) by: where p p is the number of pole pairs of the motor.The efficiency (η) has been calculated as the ratio between the hydraulic and the electric power: being γ the specific weight of water, equal to 9806 N/m 3 .For each frequency, the head curve (H N (Q)) and the power curve (P N (Q)) can be interpolated by a polynomial regression (second and third order respectively).Thus, the efficiency can be calculated by Equation ( 6).
Figure 2 displays the experimental points and the interpolation curves of the two pumps.The best efficiency of both pumps decreases as the frequency decreases.Turbomachinery affinity laws can be used to simulate the pump behaviour under variable speed [37,38] but, unfortunately, they do not get the efficiency reduction [39][40][41].Indeed, Figure 3 shows that the dispersion of experimental point on the chart (Q/N, H/N 2 ) is negligible, while it is considerable in the chart (Q/N, P/N 3 ).Thus, according to affinity laws, a unique interpolating polynomial curve (calculated on all the experimental points) has been used to calculate the head, namely: while, for the power, a polynomial best fit has been calculated only for f = 50 Hz, i.e., N = N max = 3000 rpm on the values of p N max In order to calculate the power for any N rotational speed, a correction to the affinity laws should be introduced.For each N rotational speed, the best efficiency (η N BEP ) can be calculated as: where N is the speed of the pump calculated by Equation ( 5), as shown by Figure 4. Thus, a new parameter can be introduced, namely the relative efficiency e [42]: As shown by Figure 5, the data (Q/N, e) lies on the same curve, with a little dispersion.Thus, for each value of Q and N, e can be calculated by a polynomial regression on (Q/N) while the efficiency can be calculated by Equation (10).Finally, the power P can be calculated by coupling Equations ( 6) and ( 8)-( 10):

Input Discharge Pattern
Data used in the present study describe the wastewater flow arriving at the pumping station of Coroglio (Naples, Italy) and collected along the "Arena S. Antonio" (ASA) urban basin.ASA is the largest drainage basin in Western Naples (Italy), having a total area of 1760 ha and a length of the main channel equal to 8.5 km.Originally, it was a stormwater sewer system, but after a strong urbanization of the area it now collects wastewater as well [43]; the catchment is mainly urbanized, with a total percentage of impervious area equal to 62%.The sewer starts at 159 m a.s.l., whereas its terminal section, located at Coroglio treatment plant, has an average altitude of 1 m a.s.l.Different cross-section configurations occur along the main sewer, with a terminal rectangular cross-section with a width of 9 m and a height of 3.7 m.
Runoff data is expressed in terms of water depth in the channel cross-section; data were recorded by means of an ultrasound level gauge at the terminal section of the sewer, 80 m before the pre-treatment plant in Coroglio.The time step of the recording is 90 s.The modelling of the whole system, described in [44], allowed for the estimation of the head-discharge relation for the measurement section and the subsequent estimation of flow rates corresponding to measured heads in dry periods of the year 2016.In the present paper, flow rates describing the daily wastewater pattern of the ASA basin for the day with maximum average daily discharge, which occurred on 18 June 2012 and were aggregated at the time scale of 15 min, were adopted.
The scaled daily pattern q(t) was obtained dividing the measured discharge Q exp (t) by its maximum value: Then, for each simulation, the input discharge was calculated setting the maximum input flow rate Q in max as a fraction of the BEP discharge Q N max BEP : where α is a coefficient that in this study has been set to 1, 1.5 and 2. The value of α affects two aspects of the design of the pumping station.For an assigned inflow pattern, higher values of α lead to an oversized plant, with a larger pump, and then a larger wet well (see ( 4)).On the contrary, the resilience of the plant is higher, because, even if the maximum inflow discharge increases more than the design value, the pump is able to deal with it.Once the value of α was set, then the input discharge pattern Q in (t) was obtained: The input pattern, and the effect of α on it, is shown in Figure 6 with reference to the BEP discharge of the first machine.

Plant Behaviour
The pumping head H man can be calculated from the pumping discharge, by using a quadratic head-loss formula, H 0 being the static head when the wet well is empty, H w the water level in the wet well, depending on time t and K a coefficient depending on the material and the diameter of the outlet pipe.H 0 can be set as a ratio of the BEP head at maximum speed H BEP N max , namely β: If the pump is correctly designed, the pumping head H man corresponding to Q BEP N max is close to H BEP N max .Then, the coefficient K can be calculated as follows: The effect of β on the plant head curve is showed in Figure 7, together with the pump head curves.The water level into the wet well H w can be calculated by the continuity equation: W w being the volume of water inside the wet well.If S is the cross section of the wet well, then W w = S • H w .Thus, from the continuity equation, the water level in the wet well can be written as:

Optimization Model
The optimization model has been written in order to minimize the energy requested to pump the water out of the wet well.Thus, the objective function E T has been set as follows: T being the time window and P(t) the instantaneous power requested by the pump at the instant t.This equation should be of course discretized, in order to be solved.Thus, the optimizing function becomes: i being the generic time interval whose length is ∆t.During ∆t, each variable is assumed constant.
For each i time interval, the power P i is a function of the rotational speed of the pump (N i ) and of the pumped discharge (Q i ), as expressed by Equation (12).Furthermore, the pump can be either switched on or off.Thus, a new binary variable has been added, namely the switch I i , which is equal to 1 if the pump is on at the i-th interval, and 0 otherwise.Thus, the requested power for each time interval can be calculated as follows: From Equation ( 23), the power at each time interval appears as a function of the discharge, the rotational speed and the switch.The rotational speed can be set between the value corresponding to the maximum grid frequency (namely 50 Hz, corresponding to 3000 rpm) and cannot be lower than 1500 rpm to avoid a danger heating of the machine [45].Once the rotational speed is selected, the discharge can be calculated: if the pump is switched off, the discharge is obviously zero; otherwise, it can be calculated by equating the head curve of the pump in Equation (7) and that of the plant in Equation ( 16), as follows: The water level in the wet well, H wi that appears in Equation ( 24) is related to the continuity equation, namely Equation (20), which can be written by applying a centered finite difference scheme: Actually, only the speed and the switch are independent variables, while the discharge and the water level in the wet well could be calculated by means of Equations ( 24) and (25).Nevertheless, in order to improve the coding, they can be set as decision variables, subjected to some constraints.In addition to Equation (25), which is a linear constraint involving H wi and Q i , the pumped discharge of each time step can be set between a maximum and a minimum value, as follows: The right side of such a constraint is written to force the discharge to zero if the pump is switched off.Q ad max is set to a high value, e.g., 3 times Q BEP .Furthermore, the momentum equation can be modified as follows, in order to include the switch (I i ): Such equality is satisfied if the pumping head equals the plant equation and collapses to an identity if the pump is switched off.Moreover, a minimum and a maximum bond are fixed also for the water level in the wet well: Finally, a last constraint that involves I i can be written, considering the maximum allowable starts per hour of the pump (S h max ).This last non-linear constraint can be written with reference to a new variable, that counts the number of starts for each one-hour window, S h i , that is a non linear function of I i .It can be calculated as the sum of the occurrence of the event (I j − I j−1 ) > 0), where j is the generic time interval of a one-hour window before the i instant.Hence, All the equations written above can be summarized in the following model, written for the T time window:

Resolution of the Optimization Model
The optimization problem involves linear and non linear constraints, continuous and integer variables and the objective is a non linear function of the selected variables.The problem is classified as mixed integer non linear problem (MINLP).The Basic Open-source Nonlinear Mixed INteger (BONMIN) [46] code has been chosen to solve the optimization.The selected algorithm is a branch-and-bound based algorithm [47] specifically designed for mixed integer non linear problems and has been successfully used in hydraulic problems in recent times [48].For each step, the algorithm solves the continuous relaxed problem through the Interior Point OPTimizer (IPOPT) [49] and the Coin-or branch and cut (Cbc) algorithm [50] to solve the mixed integer problem.Even if the algorithm is designed for convex problems, it can retrieve heuristic solutions in case of non-convex problems [51].The problem described by Equation ( 30) is quite complex and the convexity has not been proven herein.Thus, all the options for the resolution of non-convex problems have been selected.
The problem can be quite large, depending on the choice of the time window, T, and the length of the time interval, ∆t.Each of the four variables is a vector with n T elements, while the 11 constraints in Equation ( 30) can be written for any i interval.Thus the problem has 4 • n T variables and 11 • n T constraints.A Mixed Integer Non Linear Programming is classified as a NP-hard problem [52], where the computational time and resources exponentially increase with the number of variables.If the input pattern is considered periodic, i.e., the pattern keeps happening every day, the optimization should be performed on the whole daily pattern.The choice of the length of the time interval ∆t is crucial: it cannot be too short, in order to avoid a too large number of variables, while a too long interval produces unreliable results, due to the coarse discretization of the problem.Furthermore, the filling time (with the maximum input discharge) of a wet well is generally lower than thousands of seconds: then, ∆t should be lower than the filling time, in order to avoid the complete filling of the wet well when the pump is switched off.In this paper, ∆t has been chosen equal to 60 s.
Thus, in order to get reliable results in a reasonable time, once ∆t has been chosen, the whole day has been divided in n w time windows of T length (being n w • T = 1day), and the optimization problem has been solved for each of them.The optimal solution for the whole day, i.e., the optimal sequence of speeds and switches, has been obtained as the union of n w consecutive solutions.Hence, the optimal daily energy, E opt , can be calculated as the sum of the energy of each time window E T : In this way, only a nearly-optimal solution can be found, and in this paper the dependency of the results on the length of the window has been investigated.

Application and Results
The optimization model has been applied to several different situations.Two different wet wells have been designed for the two pumps, based on Equation (4), resulting in W = 26.6 m 3 and 15.19 m 3 respectively.A constant cross section wet well has been chosen, with a 10 m 2 base surface, then H max results 2.66 m and 1.52 m respectively.Different scenarios have been created by selecting the values of the α and the β coefficients.This means that, assigned the pump (machine 1 or 2), different input patterns have been calculated by setting the α value equal to 1, 1.5 and 2 respectively.The resulting values of Q in max are shown in Table 1.Different plants have been simulated by setting β equal to 0, 0.25, 0.5, 0.75 and 1 respectively, resulting in different values of H 0 and head loss.For each scenario, a reference value of daily requested energy, namely E re f has been calculated, as follows: Each scenario has been optimized, setting a time window of 30 min and a ∆t of 1 min, so that n T = 30.Thus, for each time window, 120 optimal variables have been found.For the first instant of the day, the water level within has been set as half of the maximum level.Table 2 shows the values of energy required by the optimal solution E opt , as well as the calculated values E cs , namely the daily energy required by CS pumping station, i.e., a pumping station where the pump operates only at its maximum speed and its starts and stops are controlled only by the water level has been also calculated.Two new parameters, namely η cs and η opt , calculated as the ratio between E cs or E opt respectively and E re f are also shown.The values of η CS demonstrate that, as the maximum inflow approaches the BEP discharge, the efficiency of the pump station increases.Furthermore, the values of η OPT are even larger than η CS , demonstrating that the optimal regulation allows a better energy efficiency of the system.
The benefit column, namely ε, shows the ratio between E cs and E opt .The higher ε values are larger than 2 and occur when β is zero.This means that, with such an optimization, more than half of the energy requested by the classical operation can be saved.The ε value decreases as β increases.Values of ε lower than 1 (resulting from values of E opt slightly higher than E cs ) occur when β is equal to 1.This can be considered a weakness of the model, since the constant speed operation is a feasible solution for the optimization model and should be selected by the algorithm (thus, a ε equal to 1 should be a lower bound).Values of ε slightly lower than 1 could be due either to the discretization or to the non convexity of the problem.Among the different scenarios, an average value of 32% energy can be saved with an optimized control.The last column of the Table shows the computational time.It is highly variable and, even if in the majority of the cases is reasonable, reaches a maximum value of 11.5 days.A pattern of the computational time is apparently not detectable.The unstable behaviour of the optimization algorithm is probably due to the non-convexity of the problem, together with the choice of dividing the entire day in multiple time windows.
Figure 8 shows the dependency of η cl and η opt on α and β.The efficiency of the optimal regulation is generally larger than the efficiency of the classical regulation, but, the difference reduces as β increases, i.e., the head loss reduces.The efficiency of both systems is higher when α is lower and when α is equal to 1 the efficiency of the optimal regulation is close to the maximum value, independently from the value of β.The efficiency of the optimal scheduling is fairly constant as β exceeds 0.25.This means that, if the head loss changes during the life of the plant, the efficiency of the system can be kept constant during the years through an optimal scheduling of the pump operations.Figure 9 shows the benefit, ε, that can be obtained when the optimal regulation is used to replace the classical system.Such a benefit decreases as β increases, i.e., as the head loss reduces, following the behaviour of η CS and η OPT .The reduction of the benefit as β increases is due to two reasons: when the head loss reduces, a reduction in the output discharge produces smaller reductions in the required head, as shown in Figure 7. Consequently, the total energy required to pump the entire volume of water slightly depends on the output discharge.Furthermore, as shown in Figure 7, as β approaches 1, the plant line intersects regions of the pump performance where the efficiency is lower.Thus, a speed regulation aimed at a reduction of the outflow discharge becomes less effective.
The highest value of ε is 3.44 and occurs when α is equal to 2.0 and β is equal to zero.This means that up to more than 70% energy can be saved.A similiar behaviour occurs for α equal to 1.5.Lower benefits can be obtained when α is equal to 1.
For the first five scenarios, the dependency of the best solution on the time window has been tested, repeating the optimization with a time window of 60 and 120 min respectively.In each of the three cases, the timestep is fixed and equal to 1 min.The results are shown in Table 3.The improvement in terms of saved energy is really low, demonstrating that a time window of 30 min is sufficient to get the optimal result.

Conclusions
In the present paper, an optimization problem is described concerning the pump scheduling of a drainage pumping station.Compared to current literature, the main difficulty lies in the pump rotational speed, that is able to vary in a given range altering the usual operational field of the machine.To comply with the peculiarities of the problem, the paper is structured into several consecutive steps: 1.An experimental campaign is undertaken to explore the effects of variable speed on the pumping efficiency.Specifically, in accordance with the affinity laws, an empirical equation is provided to compute the pumping head (Equation ( 7)), whereas a novel approach based on the concept of relative efficiency (Equation ( 10)) is provided to compute the pumping power under variable speed conditions (Equation ( 12)). 2. On the basis of the above-mentioned theoretical framework, a mixed-integer optimization problem (Equation ( 30)) is built that is made up of an objective function (the overall pumping energy) to be minimized and a set of constraints for the variables of interest.The model is also able to comply with the ON/OFF switch of the pump, and two parameters (α and β) are introduced to simulate different scenarios for the inflow discharge and the plant configuration.The influence of the time window and step for computations is also discussed.3. The model is solved for a case study (a literature sewage daily pattern provided for the City of Naples, Italy) relying on a literature algorithm, and some indicators are analyzed to test the computational performance of the algorithm and the overall energy savings given by the optimal solution for the different scenarios.The case study was developed assuming a known inflow pattern and a constant energy cost during the day.
The optimization results show that, if the pump scheduling, i.e., the pump starts and rotational speed, is optimized, meaningful energy saving can be pursued.The amount of saving depends on the plant characteristics, in terms of inflow discharge and head loss.Among the different simulated scenarios, the average ratio of saved energy is equal to 32%, with a peak value of more than 70% in the most convenient scenario.The efficiency of the optimal system increases as the difference between the pump maximum discharge and the inflow decreases, while the effect of the head losses is mild and the efficiency of the optimal scheduling is fairly constant if the head loss is not too high, i.e., the static head is higher than one fourth of the BEP head of the pump.This could mean that the efficiency of an optimal scheduled plant can be kept fairly constant over the years, even if the plant head loss changes due to the wear of the pipeline.Nevertheless, when the head loss of the pumping pipeline reduces to zero, the benefit due to an optimized control, when compared to a classical control based on the water level, becomes negligible.The optimization algorithm that has been implemented showed good performance in terms of computational time, even if in certain cases it extends to several days, with no detectable patterns.In order to keep the optimization problem small, the entire day is divided into multiple time windows.The final solution is the sum of the optimal scheduling of each time window.A final comparison shows that a time window of 30 min is sufficient to detect the optimal solution and a further extension of it does not give significant improvements in terms of energy saving.
This entire study is developed assuming hypothesis of knowing the inflow pattern and considering a constant cost of the energy.Future developments could include (i) the variability of the energy cost during the day, e.g., saving energy during the day and working more during the night and (ii) a real time optimized control to face the random variability of the inflow pattern.

Figure 1 .
Figure 1.Sketch of the wet well of the pumping station with a single submersible pump.

Figure 2 .
Figure 2. Performance curves of the two tested machines (Machine 1 on the left column, Machine 2 on the right column).

Figure 3 .
Figure 3. Experimental points of Q/N, H/N 2 and P/N 3 for the two machines (Machine 1 on the left column, Machine 2 on the right column).

Figure 4 .
Figure 4. Best efficiency variation with the rotational speed (Machine 1 on the left column, Machine 2 on the right column).

Figure 5 .
Figure 5. Relative efficiency of both machines (Machine 1 on the left column, Machine 2 on the right column).

Figure 6 .
Figure 6.BEP discharge of the first machine and input pattern depending on α.

Figure 7 .
Figure 7. Pump head curves and plant curve depending on β.

Figure 9 .
Figure 9. Variation of the benefit (ε) of the optimal system with α and β.

Table 1 .
Parameters of 30 the different studied scenarios, resulting from the choice of the machine and different values of α and β.

Table 2 .
Results of the optimization with a time window of 30 min (C.T.-computational time).

Table 3 .
Comparison of the optimization results with different values of time window (T).
Pumped head at N max rotational speed at the BEP of the pump H man Minimum allowable water level in the wet well w Volume of water inside the wet well