Hydroelectric Operation Optimization and Unexpected Spillage Indications

: It is widely known that hydroelectric power plants beneﬁt from optimized operation schedules, since the latter prevent water and, therefore, monetary wastes, contributing to signiﬁcant environmental and economic gains. The level of detail on the representation of such systems is related to how far ahead the planning horizon is extended. Aiming at the very short-term optimization of hydroelectric power plants, which usually requires the most detailed models, this paper addresses an undesired effect that, despite being already mentioned in the literature, has not been properly explored and explained yet. This effect is given by the indication of spillage by the optimizer, even when the reservoir does not reach its maximum capacity. Simulations implemented in Julia language using real power plant data expose this phenomenon. Possible ways to circumvent it are presented. Results showed that, in speciﬁc cases, spillage allows the achieving of more efﬁcient operating points by reducing the gross head and increasing the amount of water that ﬂows through turbines. Furthermore, it was veriﬁed that applying water outﬂow-based objective functions prevents undesired spillage indications, despite causing machines to operate at lower efﬁciency levels, compared with the utilization of power losses-based objective functions.


Introduction
Optimizing the short-term operation of power plants comes down to determining which generating units (GUs) must be online, and how much power each of them has to generate while meeting all system constraints. Regardless of the system being hydraulic, thermal, hydro-thermal, or even incorporating alternative energy sources in its formulation, optimal operation schedules contribute to economic and environmental gains, i.e., they help secure sustainable scenarios concerning energy supply. In Brazil, there are multiple possible planning horizons. The widest is the long-term operation planning, which takes several years into account. The narrowest is a very short-term one, in which only the day-ahead operation is considered. The former usually includes several energy sources simultaneously and does not implement detailed characteristics of these sources. As losses of accuracy cause an immediate decrease in performance, the day-ahead operation requires high levels of detail, as in [1][2][3][4][5][6][7][8][9][10][11]. These papers consider several hydroelectric operation aspects, such as head dependency when calculating efficiency, reservoir and (MILP) problem was solved by CPLEX. Yuan et al. [4] applied a modified genetic algorithm to the problem. Finardi and da Silva [5] wrote the formulation as a quadratic programming problem, which was solved by LR. Diniz and Maceira [6] added the spillage as an axis to the power generation function, thus defining a 4D model. Piece-wise linearization was applied and the problem was solved by dual dynamic programming. Yuan et al. [7] and Mahor and Rangnekar [8] proposed an enhanced differential evolution algorithm and a modified particle swarm optimization algorithm, respectively, to solve the problem. Finardi et al. [9] performed a comparison among an LR formulation, a mixed-integer nonlinear programming (MINLP) solver, and a MILP formulation obtained by linearization. The authors reported that overall better performance was obtained by the LR approach. Passos de Aragão et al. [10] proposed a model including water transfer among basins. Network Flow and Particle Swarm Optimization algorithms were applied to solve the problem, which considers a hydrothermal system. Brito et al. [11] exposed and analyzed the errors of four different linear approximation possibilities of the hydro production function. Furthermore, seven piece-wise linear models were optimized by Gurobi MILP solver and compared regarding both computational burden and quality of results.
In the Brazilian network, the National System Operator (ONS) is responsible for providing generation goals to power plants and modeling reservoirs mathematically, among other tasks. The latter allows the reservoir level to be described by a polynomial function of its volume (level-volume polynomial). In some HPPs, called run-of-the-river HPPs, the notably vast area of the reservoir causes level variation to be negligible during the horizon of one day (order of 10 −1 meters). Itaipu, located in the border of Brazil with Paraguay, is an example of run-of-the-river HPP. In such cases, ONS usually considers the level-volume polynomial as being of order zero, i.e., a constant value. Unfortunately, this makes it not possible to model spillage as a decision variable properly, as reservoir level variation cannot be accounted for.
Assuming that a reservoir level-volume function is available for the very short-term operation optimization (in this case, the day-ahead planning), this paper explains why the optimizer indicates spillage in certain situations, even if the reservoir has not reached its maximum capacity. It is worth emphasizing that the goal of this work is not to propose a new approach to solve the day-ahead operation planning of HPPs. It is, however, to raise awareness about the mentioned issue, so that researchers studying this subject can more rapidly identify the cause for the potentially undesirable results obtained. Real HPP data extracted from [11] were used to perform simulations. It is necessary to point out that some information about the HPP is not available in the mentioned paper, such as the name and location of the HPP, and turbines type (e.g., Kaplan or Pelton).
Concerning environmental aspects of power generation, papers [20][21][22] expose the importance of renewable energy sources to sustainable development. More specifically, these works study how to estimate the energy that HPPs can extract from rivers. Publications like the mentioned ones show how relevant HPPs are to Green Energy. Settling HPPs cause irreversible impacts, such as rivers course modifications and deforestation. However, once settled, their operation is cost-free regarding fuel, does not cause pollution, and emits low amounts of greenhouse gases, consequently contributing to climate change prevention. In addition, HPPs deliver social benefits, such as local economic growth due to tourism and job creation, and scientific incentives via local flora and fauna preservation projects. In this context, optimizing HPPs operation as well as avoiding unnecessary spillage help these systems become even more efficient and less harmful to the environment. Minimizing reservoir water usage while meeting power demands helps future loads to be attended and prevents severe increases in the downstream basin, which, depending on the region, can affect fishing and agriculture, among other activities. Positive and negative impacts of HPPs can be consulted in the details of [23][24][25][26][27][28].
Regarding this paper's organization, Section 2 describes the typical formulation, Section 3 exposes results obtained from simulations and discusses the aforementioned spillage effect, and Section 4 concludes this work.

Problem Formulation
The mathematical model of an HPP depends, among other aspects, on the simplifications adopted. For instance, Arce et al. [29] exploited the fact that the studied HPP is a run-of-the-river plant to model the reservoir level as a constant value during the daily planning horizon, which, in such case, imposes irrelevant impact on results. As a drawback, spillage cannot be duly used as a decision variable in this approach, as it becomes impossible to calculate how much the variation in the reservoir water will change its level. Furthermore, the fact that all GUs are identical allows the authors to assume that power demands are equally distributed among online GUs [30], which greatly simplifies the problem.
To fulfill this paper's purpose, i.e., exploring the cause of unexpected spillage indications by the optimizer during the day-ahead planning, this research aims at analyzing a more generic formulation that suits most HPPs and requires minor adaptations, if any. Therefore, the model here presented is the same from [11], which refers to very short-term operation planning and contemplates storage variations, spillage as decision variables, and GUs with different characteristics.
Expression (1) defines the objective function, which is given by the HPP's total release, i.e., the sum of total turbine water outflow (q) and spillage (s) during the whole planning horizon. Equation (2) considers incremental flow (Y), q, and s to calculate reservoir volume (v) variation. In this equation, H is equal to 3.6 · 10 −3 and is responsible for converting m 3 /s to hm 3 . Equation (3) secures demand (L) attendance according to each GU's power generation (g), with N representing the number of GUs in the HPP. Equations (4)-(7) calculate g, turbine efficiency (η), the gross head (h), and q, respectively, with each individual water flow being represented by w. In Equations (4) and (5), the net head is given by h minus the penstock power losses, which are given by the multiplication of the penstock losses coefficient (D) by w squared. In Equation (6), B represents the coefficients of the level-volume polynomial, whereas C represents the coefficients of the level-release polynomial, which grants the downstream basin level as a function of the total water release. Equation (8) establishes GUs' current ON/OFF status (u) according to previous u and to startup (i) and shutdown (o) status. Regarding u, i, and o, 1 represents being online, started-up, and shutdown, respectively, whereas 0 represents being offline, not started-up, and not shut down, respectively. Equation (9) limits how many times each GU is allowed to startup during the planning horizon based on predetermined limit values (I). Equations (10) and (11) impose minimum uptime (T up ) and minimum downtime (T down ) constraints, respectively. Equations (12)-(15) set the bounds of part of the decision variables. It is worth mentioning that Equations (14) and (15) refer to two-sided constraints that must be added to the problem since w's and g's bounds depend on u. Expression (16) represents the binary characteristic of u, i, and o. In all nomenclatures, j is an index referring to a particular GU, t is an index referring to a particular period of the planning horizon, and min and max denote minimum and maximum values, respectively, of a particular parameter or variable.
To expose and analyze the mentioned unexpected spillage indications, an alternative objective function based on power losses will be considered (Equation (17)). In general, the problem's objective function depends on the approach chosen. For instance, if the intention is to minimize operational costs in a scenario that takes into account start-up and shutdown costs, applying a power losses-based function is more intuitive since energy prices can be directly extracted from the energy market, as in [29,31].
As stated in Section 1, even though Garcia-Gonzalez and Castro [3] mentioned undesired spillage indications, they did not explain their cause. In [3], Equation (13), which represents the bounds for spillage variables, was replaced by Equations (18)- (20), in which is a tolerance that allows spillage indications only if v max − ≤ v t ≤ v max , and e t is a binary decision variable, whose unitary and null values mean "allowed" and "not allowed" spillage, respectively. In this approach, the optimization software user is responsible for choosing the value for , i.e., for determining which is the minimum reservoir volume that allows spillage. The drawbacks of this method will be discussed towards the end of the following section.

Simulations
All codes were written in Julia 1.3.1 and solved in a JuMP 0.21.1 [32] environment. The optimizations were performed by Juniper 0.6.3 [33], an open-source local solver for non-convex MINLP problems that implements a branch-and-bound algorithm. Juniper requires a separate nonlinear programming solver and recommends the usage of a separate mixed-integer programming solver, which, in this paper, are the open-source packages Ipopt 0.6.1 [34] and Cbc 0.6.6 [35], respectively. The former applies a primal-dual interior-point algorithm with a filter line-search method, whereas the latter implements a branch-and-cut algorithm together with primal and dual simplex algorithms.
As stated in [36], which performs research using Julia language, although open-source solvers usually cannot keep up with the speed of commercial solvers when studying large systems, they contribute to transparency, collaboration among developers, and savings due to licenses fees. Such affirmations are demonstrated by paper [37], which, despite using the Matlab commercial language, applies Ipopt as the solver to its studied problem.
Regarding the simulation data, all information was extracted from and can be consulted in [11]. It should be mentioned that the data concern one HPP, i.e., one reservoir. Interested readers are welcome to contact the corresponding author of this paper and request the Julia codes. This approach will not only provide the data but will also promptly enable the replication of the results yet to be presented. In [11], three distinct scenarios are considered regarding power demands, initial reservoir volume, and expected incremental water flow (water the will "feed" the reservoir throughout the planning horizon), which will be presented further ahead.
To reduce the computational burden, constraints concerning the maximum number of start-ups (Equation (9)) and minimum up and downtime (Equations (10) and (11), respectively) were discarded. It is emphasized that, by not applying these constraints, the GUs will have no limitations concerning start-ups and shutdowns, which has no impact on the intended analysis, i.e., on discussing undesired spillage indications. It is noteworthy that discarding the mentioned constraints automatically makes i and o variables unnecessary as they relate to starting up and shutting down, respectively. Consequently, Equation (8) can also be excluded from the implementation.
One relevant aspect of the HPP studied is that, among its six GUs, a first group (G 1 ) is composed of four identical GUs, whereas a second group (G 2 ) is composed of two identical GUs, though different from G 1 's GUs. For this reason, although the GUs are treated individually in the optimization process, results will be disposed in terms of the number of GUs operating in each group (N G1 and N G2 ) and their power generation and water flow.
Part of the data is crucial to perform the necessary analysis. Therefore, some information will be exhibited in this paper. Tables 1-3 show the power demands in the first, second, and third scenarios, respectively. In [11], the power demands were exposed graphically, so the data was extracted using the free software Graph Grabber, which is available at https://www.quintessa.org/graph-grabber. Turbine efficiency (η) functions' coefficients to be utilized in Equation (5) are given in Table 4, which also provides the bounds for power generation (g) and turbines water outflow (w), and the penstock losses coefficients (D). Reservoir's minimum and maximum capacities are equal to 721 and 1123.67 hm 3 . Constraint (6) is numerically calculated according to Equation (21). The upper bound for h t is equal to 75.2 m. The incremental flow in the first, second, and third scenarios are equal to 1380, 637.5, and 1300 m 3 /s, respectively.  Table 5 exhibits the results considering the power losses-based objective function described in Equation (17). It is noted that, although the reservoir has not reached its maximum capacity, there were spillage indications in periods 16 and 20. Similar cases have been reported in [3], which provides a straight forward way to change the reservoir constraints so that spillage can only be indicated by the optimizer if the reservoir level reaches a user-predefined value, as exposed in the previous section (Equations (18)-(20)). However, the mentioned paper did not explore nor explained the reason for this phenomenon to happen. To the best of this paper's authors' knowledge, this matter has not been addressed in detail in the literature yet.   Table 4. Turbine efficiency functions' coefficients, g's and w's bounds, and D values. Data from [11].  (17) as objective function was performed. However, the decision variable representing spillage had its upper bound set to zero, i.e., spillage was not allowed to be indicated. The results obtained in this approach are shown in Table 6. Researchers studying similar situations should be aware of the possibility of non-convergence of the optimization if any spillage is required to keep the reservoir volume below or at its maximum value.
By using Equation (5) together with the individual water flow (columns 4 and 7) and gross head (column 10) values provided in Table 5, one can verify that the mean efficiency values are equal to 92.98% and 92.54% in periods 16 and 20, respectively. When analyzing Table 6, these values are equal to 92.87% and 91.96%, respectively. Such results make it clear that the spillage values indicated in the first simulation increased GUs efficiency, hence reducing power losses, which is in accordance with the chosen objective function. In fact, the objective function values obtained with and without spillage were 1631.75 MW and 1636.04 MW, respectively. Table 5. Results with power losses-based objective function-scenario 1. From Equation (5), it is easy to verify that the gradient of the efficiency (which defines the direction of its steepest ascent), i.e., its partial derivatives regarding the gross head ∂η jt ∂h t , and the water flow ∂η jt ∂w jt are given by Equations (22) and (23), respectively. It is noteworthy that the gradient must be calculated for each GU. However, since this paper considers that each group is composed of identical units, it is enough to refer to each group's gradient (∇ G 1 and ∇ G 2 ), i.e., the gradients have the same value in a particular group.
It is observed that, at P 1 , G 1 's turbines efficiency will reach higher values if h decreases and w G1 increases; or if h's increase is lesser than 0.4397 times w G1 's increase; or if h's decrease is greater than 0.4397 times w G1 's decrease. In fact, when observing G 1 's GUs in period 16 of Table 5 (which concerns the simulation in which spillage was allowed to happen), it is noted that h decreased 1.04 m, and w G1 increased 2.39 cubic meters per second, compared with period 16 of Table 6, which means that such alterations fit the criteria for the efficiency of G 1 's turbines to be increased.
At P 2 , G 1 's turbines efficiency will reach higher values if h decreases and w G2 increases; or if h's increase is lesser than 6.7854 times w G2 's increase; or if h's decrease is greater than 6.7854 times w G2 's decrease. By analyzing G 2 's GUs in period 16 of Table 5, it is possible to note that h decreased 1.04 m, and w G2 increased 5.22 cubic meters per second, compared with period 16 of Table 6. Once again, the variation in variables' values fit the criteria for the efficiency to be increased. In this case, for G 2 's turbines.
Although the analysis presented in the last two paragraphs will not be carried out for every period where unexpected spillage occurred, it demonstrates in detail that the observed increases in efficiency due to spillage indications are in accordance with the mathematical model of the HPP.
Regarding the second and third scenarios simulated with the power losses-based objective function (Equation (17)), Table 7 exposes the results for scenario 2. In this scenario, convergence is impossible if the spillage decision variable has its upper bound set to zero, regardless of the objective function chosen. This fact is due to the incremental flow and initial reservoir volume values in this scenario, which requires that water is spilled during the operation so that the reservoir can be kept within its bounds. Tables 8 and 9 show the results for scenario 3. The former refers to a simulation in which spillage was allowed to occur, whereas the latter refers to a simulation in which spillage variables were blocked, i.e., had their upper bounds set to zero.
Tables 7 and 8 show that the indicated spillage values are more than the necessary to keep the reservoir below or at its maximum capacity, as it did not reach its maximum volume of 1123.67 hm 3 . This means that hydro resources were wasted. The study presented so far shows that these undesired spillage indications are connected to the characteristics of the turbine efficiency curves, which are typically hill-shaped curves, regardless of the HPP. Therefore, if the spillage is considered in the formulation as a decision variable, it is very likely that the optimizer will indicate spillage when using a power losses-based or efficiency-based objective function, even if the reservoir does not reach its maximum capacity. It is important to emphasize that the cause for unexpected spillage indications is not the objective function itself, given the fact that the optimizer simply mathematically searches optimal values while not violating the constraints. However, due to how turbine efficiency curves behave, there are situations where the optimizer obtains these optimal values by indicating spillage from the reservoir to the downstream basin. As previously exposed, such excessive spillage allows the achieving of gross head and turbines water outflow values that imply higher efficiency. Table 7. Results with power losses-based objective function-scenario 2.  Table 8. Results with power losses-based objective function-scenario 3.  Table 9. Results with power losses-based objective function and null spillage upper bound-scenario 3. A direct and simple way to circumvent such a problem is to use a water outflow-based objective function, as described in Equation (1). As the goal in this approach is to use the least possible amount of water, the optimizer will indicate spillage only in situations where the reservoir would exceed its maximum capacity. Furthermore, this approach is aligned with environmental concerns, which is a relevant gain. To illustrate this situation, Tables 10-12 show the results of simulations using the same data of previous simulations, however applying Equation (1) as the objective function. These tables correspond to the first, second, and third scenarios, respectively.

Period N G1 g G1 (MW) w G1 (m 3 /s) N G2 g G2 (MW) w G2 (m 3 /s) s (m 3 /s) v (hm 3 ) h (m)
The following observations can be done regarding the water outflow-based approach: • In Table 10, which refers to scenario 1, it is noted that no spillage was indicated as the reservoir is always within its bounds. It is worth mentioning that this approach resulted in the usage of 111.22 hm 3 of water in total, whereas the simulations exposed in Table 5 (power losses with spillage) and Table 6 (power losses without spillage) resulted in 111.51 hm 3 and 111.26 hm 3 , respectively; • Table 11, which refers to scenario 2, clearly resulted in less spillage than Table 7 (power losses with spillage). The former used 51.96 hm 3 of water in total, whereas the latter used 54.08 hm 3 . An important aspect to be highlighted is the fact that the reservoir reached its maximum volume in period 18 of Table 11, i.e., some water must be spilled to keep the reservoir within its bounds.

•
Regarding scenario 3, which, as scenario 1, does not require any spillage to keep the reservoir volume within its bounds when utilizing a water outflow minimization, applying the water outflow-based approach resulted in total usage of 133.84 hm 3 of water, whereas the power losses-based approach with spillage (Table 8) would require 134.49 hm 3 . The power losses-based approach without spillage (Table 9) would also result in total water usage of 133.84 hm 3 . Table 10. Results with water outflow-based objective function-scenario 1.  Table 11. Results with water outflow-based objective function-scenario 2.   Table 12. Results with water outflow-based objective function-scenario 3. As a summary of the research here presented, Table 13 exhibits GUs mean efficiency values in each period of each scenario in which unexpected spillage was observed. In this table, PL, PLBS, and WO refer to the power losses approach, the power losses with blocked spillage approach, and the water outflow approach, respectively. This table makes it clear that, compared with the simulations in which spillage variables had their upper bounds set to zero, the occurrences of spillage actually allows the optimizer to obtain operating points that present higher efficiency when applying a power losses-based objective function. As stated before, such a phenomenon is due to the nature of hydraulic turbine efficiency curves. Spillage can decrease the gross head and increase turbine water outflow in a way that GU efficiency values are increased, which leads to fewer power losses.

Period N G1 g G1 (MW) w G1 (m 3 /s) N G2 g G2 (MW) w G2 (m 3 /s) s (m 3 /s) v (hm 3 ) h (m)
As stated before, the second scenario cannot be solved unless spillage decision variables are "free". This fact is due to the initial incremental flow and reservoir volume values. Nonetheless, Table 13 shows that minimizing water usage (WO column) implies lower efficiency levels in comparison with the power losses minimization approach (PL column). However, the former strategy does not indicate excessive spillage.
It is important to highlight that undesired spillage indications are consequences of the hill-shaped form intrinsic to hydraulic turbine efficiency curves, regardless of the turbine type. Figure 1 depicts the efficiency function for G 2 's turbines. In this curve, one can notice that if the water flow (w) is equal to 195 m 3 /s, the efficiency reaches higher values as the net head (h net ) decreases from 75 to 73 m. For any h net value, the efficiency reaches higher values as w increases from 195 to 225 m 3 /s. Intuitively, spilling water from the reservoir causes a decrease in its level and an increase in the downstream basin level, hence decreasing the net head and requiring more water to flow through the turbines to meet the given power demand. In other words, no matter the turbine type and as long as the function is given by a hill-shaped curve, spillage can lead to higher efficiency depending on the operating point. The affirmations presented in the previous paragraph, in conjunction with the development exposed in this section, make it clear that the HPP operation optimizer can indicate spillage even when the reservoir is not at its maximum capacity, leading to a waste of resources as a consequence. As the cause of these indications relates to the hill-shaped efficiency curve, they can occur in any HPP in which this type of function rules the turbines. This affirmative is true regardless of any other operation characteristic, such as the number of generating units (either being identical or not) in the power plant, the length of the planning horizon, and the study context concerning single or multiple HPPs. In summary, this study contributes to the development of coherent operation optimization codes, which can enable more reliable results, help save resources, and prevent operators from making questionable decisions regarding reservoir water spillage.
Regarding the approach presented in [3] (Equations (18)- (20)), since it only allows water to be spilled when the reservoir volume (v) is greater than or equal to v max − , the method tends to indicate less undesired spillage in comparison with the results presented in Tables 5, 7 and 8. However, the optimization can still result in excessive water usage if non-null spillage values are feasible, as exposed in this paper. In other words, this approach does not prevent the optimizer from indicating wasteful spillage. Furthermore, it is possible that spilling water only when v ≥ v max − is not enough to keep the reservoir below or at its maximum capacity if 's value is too low, hence making it impossible for the optimizer to reach convergence. Therefore, it is more appropriate to adapt the formulation and apply a water outflow-based objective function.
The limitation of this work relates exclusively to the simple modeling of turbine efficiency curves (Equation (5)), which was performed by second-order polynomial fitting, as in [11]. In general, these curves are provided by turbines manufacturers as a data set that contains numerous three-dimensional samples (net head, water flow, and efficiency), which are then used to obtain a "black box" that grants the efficiency (given net head and water flow values). For instance, these samples can be fitted by one polynomial (higher polynomial order implies higher accuracy), or by multiple polynomials obtained via segmentation of the original data set. More complex and precise approaches to model turbine efficiency may impact the exposed undesired spillage indications regarding their magnitude, which defines a possibility of future investigations concerning this subject. However, no matter the model, these indications are always prone to happen when applying power losses-based objective functions due to the hill-shaped efficiency curve, as explained in this paper. For the readers who intend to apply a monetary cost-based objective function that includes start-up and shutdown costs, though considering spillage as a decision variable, writing part of the objective function as costs related to power losses can be troublesome, as shown in this paper. To engage in such an approach and avoid undesired spillage indications, one could set monetary values to the storage water, thus making it possible to use water outflow values instead of power losses in the objective function, as in [2].
The "BenchmarkTools" Julia package was used to carry out a computational effort analysis of the simulations, which were performed in a computer with the following characteristics: AMD Ryzen 5 2600 six-core processor with 3.4 GHz, 8 GB of RAM, and a 64 bits Windows 10 operational system. Table 14 exposes the mean and standard deviation values obtained by performing four simulations of each scenario. Both power losses (PL) and water outflow (WO) objective functions were considered. In Table 14, the low standard deviation values (in comparison to mean values) show the consistency of the optimization lengths given a particular scenario and objective function. In contrast to scenarios 1 and 2, scenario 3 presents power demands that vary lightly throughout the day (see Tables 1-3). Such characteristic simplifies the process and is responsible for scenario 3's much lower optimization duration. Although both scenarios 1 and 2 present significant demand variation, the latter is the only one in which water must be spilled from the reservoir so that its volume remains within its bounds. The assignment of non-null values to spillage decision variables also impacts computational effort, which justifies scenario 2's higher optimization duration. Regarding the objective functions, WO required less time in scenarios 1 and 3, whereas PL prevailed in scenario 2. However, the differences were considerably small in all cases. In general, the required times to optimize the daily operation of the power plant were all viable, with the longest one being equivalent to less than six minutes.

Conclusions
This paper explained the cause of undesired spillage indications when optimizing the operation of hydroelectric power plants. Analysis of mathematical aspects showed that such a phenomenon is related to the nature of hydraulic turbine efficiency curves. In specific cases, power losses-based objective functions lead the optimizer to reach higher efficiency levels when spilling reservoir water and, consequently, reducing the gross head and increasing turbines water flow. Results have shown that reservoir resources can be wasted according to the objective function applied in the problem formulation when spillage is considered a decision variable. It is emphasized that renewable energy sources greatly contribute to environmental balance. Therefore, properly optimizing their operation helps secure resources for future demands. By exposing the studied matter, the authors of this paper intend to raise awareness on its cause so that fellow researchers more quickly identify similar situations. For future studies, it is intended to assess how much different turbine efficiency models impact the magnitude of undesired spillage indications.