Next Article in Journal
Convolutional Neural Network and Bidirectional Long Short-Term Memory (CNN-BiLSTM)-Attention-Based Prediction of the Amount of Silica Powder Moving in and out of a Warehouse
Previous Article in Journal
Strategic Siting of Direct Air Capture Facilities in the United States
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Monthly Hydropower Scheduling of Cascaded Reservoirs Using a Genetic Algorithm with a Simulation Procedure

by
Deji Baima
1,
Guoyuan Qian
1,
Jingzhen Luo
1,
Pengcheng Wang
1,
Hao Zheng
2,* and
Jinwen Wang
3,4
1
School of Civil and Hydraulic Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, China
2
Water Resources Department, Changjiang River Scientific Research Institute, 23 Huangpu Road, Wuhan 430010, China
3
Institute of Water Resources and Hydropower, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430074, China
4
Hubei Key Laboratory of Digital River Basin Science and Technology, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 650041, China
*
Author to whom correspondence should be addressed.
Energies 2024, 17(15), 3756; https://doi.org/10.3390/en17153756
Submission received: 17 May 2024 / Revised: 15 July 2024 / Accepted: 26 July 2024 / Published: 30 July 2024

Abstract

:
This study integrates genetic algorithms with simulation programs, applying the genetic algorithm’s (GA) fitness calculation within the simulation to reduce complexity and significantly improve the efficiency of the optimization process. Additionally, the simulation introduces the concept of “Field Leveling” (FL), utilizing a push–pull strategy to explore more space for absorbing and utilizing unnecessary spillage for energy generation, thereby maximizing electricity production and ensuring optimal reservoir scheduling. Two methods are provided, namely the field-leveling genetic algorithms GAFL1 and GAFL2. GAFL1 involves only pushing and does not include a push–pull process; thus, it cannot optimize spillage. On the other hand, GAFL2 implements a complete push–pull strategy, continuously exploring additional space to absorb and utilize unnecessary spillage. Both GAFL1 and GAFL2 achieved reasonable results; specifically, compared to SQP, GAFL1 improved firm yield by 8.3%, spillage increased by 2.2 times, and total energy decreased by 1.2%. GAFL2, building on the basis of GAFL1, effectively reduces spillage under all hydrological conditions without affecting the highest priority of stable output. However, the impact of reducing spillage on energy generation is not consistent; in wet and dry years, reducing spillage increases energy generation. However, in normal years, a reduction in spillage corresponds with decreased energy generation.

1. Introduction

The continuous growth in global energy demand and the pressures of environmental conservation compel us to explore sustainable energy solutions. In this context, hydropower systems are considered important contributors to advancing sustainable energy development, despite the necessity to minimize their impact on ecosystems during implementation [1]. Specifically, the functions of reservoirs extend beyond power generation to include the regulation and management of water resources, which are crucial for sustaining human development. However, the effective management of water resources faces increasing challenges, such as water scarcity and environmental degradation, necessitating urgent improvements [2,3]. In basin development, the strategic construction of cascaded reservoirs and hydroelectric stations maximizes the utilization of water resources and hydraulic potential, becoming a vital strategy for comprehensive river resource development [4]. The optimized operation of these facilities, especially hydraulic scheduling, not only enhances the efficiency of reservoirs but also promotes the simultaneous conservation of water resources and energy while effectively reducing operational costs and environmental impacts. Therefore, developing optimal operational strategies to address water scarcity and meet the demands for responsible resource management in water resource engineering is critical to achieving sustainable economic, social, and environmental goals [5].
Scheduling cascaded hydroelectric stations in a complex interconnected reservoir system is a daunting challenge due to intricate nonlinear relationships, numerous constraints, input uncertainties, and dynamic correlations among key variables such as reservoir storage and head. This issue has garnered widespread attention in both academic and engineering circles. Two main categories of optimization methods are employed to tackle this challenge. Mathematical programming methods—such as linear programming (LP) proposed by Yoo et al. [6], dynamic programming (DP) by Bellman [7], and mixed integer programming (MIP) [8]—excel in solving nonlinear and non-convex problems, making them well suited for long-term hydro-scheduling. On the other hand, intelligent algorithms, including particle swarm optimization (PSO) generated from the foraging behavior of birds [9] and genetic algorithms (GAs) simulating the survival of the fittest principle in nature [10,11], offer practical solutions by effectively handling nonlinear constraints and data mining algorithms [12]. However, complex scenarios, especially in multidimensional, nonlinear, and large-scale reservoir optimization, may pose challenges for traditional optimization algorithms, leading to issues like local optima and feasibility problems within the solution space [13]. The significant challenge often encountered with intelligent algorithms lies in their ability to efficiently generate a diverse set of solutions evenly distributed across the entire feasible region, which is crucial for optimizing complex problems. This challenge becomes particularly pronounced when dealing with intricate constraints. A comprehensive review by Nicklow et al. [14] highlights the importance of intelligent algorithms in water resource planning and management but underscores their limitations in handling constraints effectively. To tackle this issue, intelligent algorithms employed in optimal reservoir operation frequently resort to penalty functions to convert constrained problems into unconstrained ones [15,16]. Therefore, the integration of intelligent algorithms with mathematical programming methods has become a growing trend in solving complex hydro-scheduling problems. Saadat et al. [17] introduced an enhanced stochastic dynamic programming method, RISDP, which overcomes infeasible conditions and dimensionality issues, improving reservoir operation efficiency and boosting performance by 15% in real-case studies. Rani et al. [18] introduced the DP-PSO algorithm, combining dynamic programming with particle swarm optimization. Li et al. [19] presented IGA-SA, a hybrid approach of genetic algorithms and simulated annealing. Through parallel computing, efficiency in solving is enhanced, making them highly suitable for large-scale parallel processing.
The genetic algorithm (GA), proposed by American scholar John Holland in the 1970s, is an optimization algorithm based on natural selection and genetic principles used to solve search and optimization problems. However, solving the optimization scheduling of cascaded reservoirs with complex constraints and multiple objectives is a formidable challenge for the GA [10,11]. In practice, Zhang et al. [20] addressed issues such as a reduced convergence capability and low search efficiency observed in traditional algorithms during the evolutionary process by proposing an improved genetic algorithm combined with parallel computing technology, effectively enhancing model convergence. Wang et al. [21] introduced a parallel genetic simulated annealing (GSA) algorithm, effectively combining genetic algorithm and simulated annealing to avoid premature convergence. Reddy et al. [22] successfully proposed a multi-objective genetic algorithm (MOGA), applied to the Badra reservoir system. Li et al. [23] improved the simulated annealing genetic algorithm by introducing niche technology to avoid local optima. Li et al. [24] introduced IGA-SA, a hybrid approach of genetic algorithm and simulated annealing, superior to traditional genetic algorithms. Arunkumar et al. [11] combined the chaos technique with a genetic algorithm to achieve better results. Fayaed et al. [25] proposed an integrated random dynamic genetic algorithm utilizing artificial neural networks. Kan Yang et al. [26] proposed effective optimization strategies for cascaded reservoir systems. Liu, S et al. [27] addressed the HUC problem by combining DP with the GA, minimizing the overflow volume and unit operation time in inefficient areas. Xu, B. et al. [28] proposed an improved genetic algorithm utilizing random operators for optimization within the feasible range of variables by gradually converting reservoir constraints into singular variable constraints and considering crossover points or mutation points in adjacent periods.
This study addresses the common challenges encountered when using intelligent algorithms to manage complex system scheduling, particularly in generating diverse and uniformly distributed feasible solutions for reservoir management. Our main objective is to optimize reservoir energy management and scheduling by integrating genetic algorithms (GA) with simulation programs, thereby improving overall efficiency and output. This paper introduces an improved method for genetic algorithms, where the fitness function calculations are performed directly within the simulation program, significantly reducing the computational load on the genetic algorithm. This integration minimizes the risks of local optima and dimensionality disasters while significantly enhancing the efficiency of the optimization process. This study incorporates the concept of field leveling (FL) [29] into the simulation program, which mimics the practice of soil leveling using a push–pull strategy to explore additional spaces for absorbing and utilizing excess spillage for power generation. We offer two strategies: GAFL1 and GAFL2. GAFL1 involves only forward pushing and does not engage in a push–pull process; thus, it does not optimize spillage. In contrast, GAFL2 implements a complete push–pull strategy, continuously exploring additional spaces to absorb and utilize unnecessary spillage for optimized energy use. This research optimizes reservoir scheduling through innovative strategies, achieving a maximized power output and enhanced operational efficiency. The practical application of these strategies can significantly benefit the system, including increased energy generation and enhanced adaptability to various hydrological conditions, providing a new optimization tool for reservoir management.

2. Problem Formulation

Monthly hydropower scheduling is formulated to maximize the firm power output and energy production sequentially during a planning horizon, expressed as
m a x W 1 · F + W 2 i = 1 N t = 0 T 1 P i t
where W1 and W2 are weights represented as W 1 W 2 , prioritizing firm power output (F) over energy production; i and t are subscripts for the reservoir and time-step, respectively; Pit is the power output in MW in time-step t.
Constraints include the following:
(1)
Water balance is
V i , t + 1 = V i t + j Ω i Q j t + I i t Q i t Δ t 24 × 3600 100,000,000
with
V i , 0 = V i i n i V i , T = V i e n d
where Vit is the storage in hm3 at the beginning of time-step t of reservoir i; Qit is the release in m3/s in time-step t from reservoir i; Iit is the local inflow in m3/s in time-step t into reservoir i; Ω ( i ) is the set of reservoirs immediately upstream of reservoir i; t is the number of days in time-step t; Viini and Viend are initial and target storages in hm3 at the beginning and end of the planning horizon, respectively. In this study, both are substituted with the normal water level.
(2)
Upper and lower bounds on storage or release are expressed as
V i dead V i t V i t max Q i min Q i t Q i max
where V i dead is dead storage in hm3 of reservoir i, representing the lower storage volume. V i t max is upper bound on the storage at the beginning of t of reservoir i, equal to the flood control limited storage during flooding seasons and the normal storage during dry seasons; Qimin and Qimax are lower and upper bounds on the release from reservoir i in time-step t.
(3)
Firm hydropower output is expressed as
i = 1 N   P i t F
where N = number of hydroplants/reservoirs.
(4)
The capacity of hydropower output due to the release available is
P i t A i · Q i t · h i t
and water heads are expressed as
  P i t G i m a x h i t = m i n c i 0 h i t + d i 0 , c i 1 h i t + d i 1
which, as illustrated in Figure 1, is equivalent to
P i t c i 0 h i t + d i 0 P i t c i 1 h i t + d i 1
with
h i t = Z i u V ¯ i t Z i d Q i t = α i · V ¯ i t V ^ i ( 0 ) β i + Z ^ i u χ i Q i t Q ^ i 0 δ i + Z ^ i d
V ¯ i t = V i t + V i , t + 1 2
where Ai is the power generating efficiency in MW·s/m4; hit is the water head in time-step t of hydroplant i; Gimax (·) is the capacity of hydropower output of i, a function of water head; c i 0 , c i 1 , d i 0 and d i 1 are coefficients to be estimated to fit the power output capacity with piecewise linearization; Ziu (·) and Zid (·) are forebay and tailwater elevations, dependent of the water storage and release, respectively, of reservoir i; α i , β i , V ^ i ( 0 ) and Z ^ i u are coefficients/parameters to be estimated for the relationship curve between storage and forebay water level.

3. Solution Techniques

Instead of solving the original problem directly, this work applies a genetic algorithm (GA) to
m i n S i t f = S i m u S 1 , S 2 , , S N type = 1
which also determines the solution
V , Q , P = S i m u S 1 , S 2 , , S N type = 2
where Sit represents the decision variables, interpreted as the storage targeted in hm3 at the beginning of time-step t for reservoir i; Simu (···) is a simulation procedure that returns the fitness of an individual solution represented by a decision matrix
S = S 1 , S 2 , , S N
when type = 1 or the hydro-scheduling solution, then
V = V 1 , V 2 , , V N Q = Q 1 , Q 2 , , Q N P = P 1 , P 2 , , P N
when type = 2, then
S i = S i 0 , S i 1 , , S i , T 1 , S i , T V i = V i 0 , V i 1 , , S i , T 1 , V i , T Q i = Q i 0 , Q i 1 , , Q i , T 1 P i = P i 0 , P i 1 , , P i , T 1
where the initial storage is always
S i 0 = V i 0 = V i ini
All the constraints will be satisfied after the simulation, and the fitness of an individual child/solution in the GA will be determined as
S i m u S t y p e = 1 = W 1 F + W 2 i = 1 N t = 0 T 1 P i t
With
F = m i n 0 t T 1 i = 1 N P i t
An individual solution in the GA is represented by the decision matrix S, which, serving as the target storages at each time-step t, may not be the same as the storages at the scheduling solution determined after the simulation. However, it is essential that the S, when generated randomly in the GA and happening to be optimum, should be equal to the storage V at the scheduling solution to theoretically ensure the capability of the GA to find the optimum solution to the problem.

3.1. Methodology

3.1.1. Genetic Algorithm

The genetic algorithm (GA) is an optimization technique inspired by the principles of biological evolution. It is a powerful tool for addressing complex problems by iteratively exploring and enhancing potential solutions. Emulating the essence of Darwin’s theory of evolution, the GA progressively refines solutions through a process that mirrors natural selection, genetic crossover, and mutation. This approach allows the GA to navigate solution spaces efficiently, identifying optimal outcomes in a manner akin to the evolutionary mechanisms observed in nature.
As illustrated in Figure 2, an improved description of the genetic algorithm encompasses the following pivotal stages:
(1)
Randomization generates an initial population of individuals randomly representing potential solutions to the problem under consideration.
(2)
Evaluation to assess an individual’s fitness in the population, involving the execution of a simulation program to gauge its suitability or quality.
(3)
Termination criterion to make a crucial decision to determine whether the algorithm should proceed with its iterative process, relying on whether the predefined number of iterations or other criteria have been satisfied.
(4)
Selection to choose parent solutions from the present population at probabilities determined by the fitness values, serving as a fundamental aspect of the genetic algorithm’s underlying principles.
(5)
Crossover between parent individuals by exchanging genetic information to generate child solutions as a fundamental genetic mechanism to create novel potential solutions.
(6)
Mutation to introduce diversity and explore a more comprehensive solution space, enhancing the algorithm’s capacity to break free from local optima and uncover improved solutions.
(7)
Iteration to steadily approach improved solutions until the predefined termination criteria are satisfied, guaranteeing that the algorithm iteratively hones its solutions across multiple generations.

3.1.2. Simulation at a Time-Step

The simulation will travel from the upstream to downstream reservoirs one by one, moving forward in time from the initial storage at the beginning of the scheduling horizon or backward from the target storage at the end of the scheduling horizon, either guided with or without minimizing spillages during the simulation procedure.
Figure 3 illustrates how the simulation works at a time-step when moving forward to determine the storage (Vi,t+1) at the end of this time-step (t), given the storage (Vit) at the beginning of this time-step and the target storage (Si,t+1) at the end of this time-step. The simulation looks as follows.
(1)
The lower and upper bounds on storage and outflow (4) are enforced with
V i , t + 1 m a x V i dead , m i n V i , t + 1 m a x , V i , t + 1
Q i t max Q i m i n , m i n Q i m a x , Q i t
(2)
The water is balanced when updating the storage and release with the following:
Q i t V i t V i , t + 1 Δ t 0.24 × 0.36 + j Ω i Q j t + I i t
V i , t + 1 V i t + j Ω i Q j t + I i t Q i t Δ t 0.24 × 0.36
The hydropower output is determined by considering its capacity and available release, as follows:
P i t m i n A i Q i t 0 h i t , c i 0 h i t + d i 0 , c i 1 h i t + d i 1
where Q i t 0 is the available release originally determined by the target storage (Si,t+1) and should also be targeted when updating the water head to accurately estimate the hydropower output capacity.
(3)
The iteration of updating water heads, indexed with k, continues until convergence, as follows:
V i , t + 1 k V i , t + 1 k 1 ε
where ε is a small converging accuracy value.
(4)
The release will be updated with
Q i t = P i t A i h i t
to reduce spillages if the simulation is guided by minimizing spillages.
This one-step simulation forward in time, shown in Figure 3, is readily adaptable to the one-step backward simulation, only by changing all the subscripts “t + 1” to “t” and keeping all the others unchanged, except for the storage update, which should be replaced with
  V i t : = V i , t + 1 ( j Ω ( i ) Q j t + I i t Q i t ) ( Δ t 0.24 × 0.36 )
The one-step backward simulation determines the storage (Vit) at the beginning of this time-step (t), given the storage (Vi,t+1) at the end of this time-step and the target storage (Sit) at the beginning of this time-step.

3.1.3. Simulation Procedure during a Scheduling Horizon

The simulation procedure will take different strategies dependent on whether or not guided with minimizing spillages, which in many cases will help yield a satisfactory solution very fast but may not be able to theoretically ensure the optimum be derived since the optimum, when randomly generated, is very likely to be modified when updating the outflow with Equation (25) to reduce spillages. The simulation procedures are summarized as follows:
(1)
Simulation procedure without minimizing spillages:
For hydropower reservoirs taking turns one by one from the upstream to downstream, let
V i 0 = V i ini S i T = V i end
Then, do the one-step forward simulation for t = 0,1, …, T − 1. This procedure ensures that a solution randomly generated in the GA, if it is also a feasible solution to the problem, will be the same as the storages at the scheduling solution since the simulation procedure can achieve the target storages without violating any constraints.
(2)
Simulation procedure with minimizing spillages:
Here, the FL algorithm is used to minimize spillages. The FL algorithm simulates the process of land leveling and iteratively adjusts the storage capacity through the “one-step forward,” “one-step backward,” and “one-more-step forward” steps to optimize the storage capacity for the following calculation. The following are the specific steps used, and each step is performed sequentially for upstream to downstream reservoirs:
(a)
Set boundary condition with Equation (27), then move forward to absorb spillages but keep the storage at the end of the scheduling horizon to be decided later by performing a one-step forward simulation for t = 0,1, …, T − 2.
(b)
Take Equation (3), then move backward to absorb spillages but keep the storage at the beginning of the scheduling horizon unchanged by performing a one-step backward simulation for t = T − 1, T − 2, …, 1.
(c)
Carry out step (a) but for t = 0,1, …, T − 1 to allow the storage to be changeable at the end of the scheduling horizon.
Depending on whether spillages are minimized or not, this study will investigate two strategies, GAFL1 and GAFL2, representing genetic algorithms without and with minimizing spillages in the field-leveling procedure, respectively.

4. Case Studies

4.1. Engineering Background

The Jinsha River, situated in the upper reaches of the Yangtze River in China, boasts abundant hydropower potential owing to its substantial elevation drop and ample water resources. The river hosts several large-scale hydropower facilities, including the Wudongde plant, with an installed capacity of 10.2 gigawatts (GW), and the Xiangjiaba plant in Yunnan Province, which generates an impressive annual output of 30 billion kilowatt-hours (kWh). The Jinsha River’s cascaded reservoir system comprises four critical reservoirs arranged from upstream to downstream: Baihetan, Wudongde, Xiluodu, and Xiangjiaba. Table 1 provides an overview of some fundamental characteristics of these four reservoirs.

4.2. Results of Curve Fitting

Table 2 provides the coefficients defining the two linear equations to establish the connection between the maximum hydropower output and water head, as a segmented linearization adopted.
Table 3 and Table 4 presents the parameters estimated using the nonlinear least squares method in PyCharm, which are used to define two exponential functions representing the relationship curves of the forebay water level vs. reservoir storage and the tailwater level vs. release. The nonlinear least squares method is a powerful optimization technique that fits nonlinear models to data by minimizing the sum of squared residuals. The coefficient of determination (R²) measures the fitting accuracy, which has an outstanding value of 99% for all four reservoirs. Figure 4 and Figure 5 illustrate the relationship curves for forebay water levels and reservoir storage volumes, as well as tailwater levels and discharge volumes, respectively, providing a visual representation of the fitting outcomes.

4.3. Data Sources

The runoff series generally follows a Pearson’s Type III probability distribution; statistical parameters are determined using frequency analysis to establish criteria for categorizing runoff into wet, normal, and dry years. For instance, the hydrological years 2000 (wet year), 1978 (normal year), and 1976 (dry year) are selected as examples, with corresponding probability percentages of 7.5%, 50.9%, and 86.0%, respectively. To illustrate hydrological trends, Figure 6 depicts the inflow volume into the headwater power station for these representative years, along with the inflow volume data from 1983 used for a comparison with SQP [30].

4.4. Results Analysis

4.4.1. Comparison GAFL1 vs. SQP

This study explores various methods employing traditional genetic algorithms, which, however, prove to be insufficient for a comparison with the current GAFL1 due to challenges in effectively handling the complexity of constraints. Thus, a comparison has was out between the outcomes obtained from the present method and those derived from a validated Sequential Quadratic Programming (SQP) approach [30]. The SQP approach formulates quadratic programming by linearizing the hydropower output, capacity, and constraints for cascaded reservoirs through a first-order Taylor expansion.
To ensure a scientifically rigorous and impartial comparison, both algorithms are configured with the same objective function, two weight factors (in a 1000:1 ratio), and water head, outflow, and generating discharge constraints. To attain essential parity and maintain fairness and accuracy in the comparison, the original water head formula (Equation (9)) is converted into Formula (28), ensuring that both models adhere to identical constraints.
h i t = Z i u V ¯ i t Z i d Q i t         = α i V ¯ i t V ^ i 0 + Z ^ i u β i Q i t Q ^ i 0 + Z ^ i d
Table 5 presents the firm yield, total energy, and spillage results obtained with the GAFL1 (genetic algorithm with minimizing spillages in field-leveling procedure) and SQP (Sequential Quadratic Programming) methodologies. Despite the GAFL1 causing a spillage that is more than double and subsequently reducing total energy by approximately 1.2%, it outperforms SQP with a notable 8.3% improvement in firm yield at the highest priority level.
Table 6 elucidates the scheduling outcomes for total energy and spillage in various reservoirs using the GAFL1 and SQP methods. A discernible observation reveals that the increase in spillage from the GAFL1 approach is primarily concentrated in the Baihetan and Xiluodu Reservoirs. This augmentation leads to a less favorable total energy outcome when compared to SQP for these two reservoirs. Notably, the Xiangjiaba Reservoir under the GAFL1 methodology registers a spillage result of zero, a remarkable deviation from SQP, entailing a substantial reduction of 1884.4 hm3 in spillage and, consequently, a conspicuously prominent enhancement in total energy.

4.4.2. Comparison between GAFL1 and GAFL2

The field-leveling procedure [29] utilizes an iterative push and pull strategy to assist the simulation in reducing spillages, thereby creating the potential to generate more hydropower energy. Figure 7 illustrates the convergence processes of the objective function values in GAFL1 and GAFL2, with both showing a monotonically increasing trend, ensuring convergence toward the limit. The most notable improvements for both methods were observed in the second iteration, and convergence was achieved in the 15th and 10th iterations, respectively.
Table 7 outlines the CPU times for both methods and the genetic algorithm’s parameters. It is noteworthy that the CPU time for GAFL2 exceeds that of GAFL1, primarily due to a significant enhancement in the push–pull process within GAFL2. The parameters for these methods are based on prior research and have undergone subsequent refinement [31]. The stopping criterion for the algorithm is termination upon the completion of the designated number of iterations.
Table 8 compares the performance of GAFL1 and GAFL2 in terms of firm yield, total energy, and spillage for four reservoirs on the upper Jinsha River, considering three typical hydrologic scenarios: a wet year, a normal year, and a dry year. The results indicate that in all hydrological conditions, the FL technique successfully minimizes spillage without compromising the priority of firm yield (Increases or decreases are marked with ↑ and ↓ in the table). Based on these findings, it is evident that reducing spillage enhances energy in both dry and wet years. However, the impact of this reduction is not consistent; during normal years, while spillage decreases, the total energy also increases. This suggests a variable influence of spillage control on energy, depending on the hydrological conditions of the year.
To elucidate the anomaly observed in normal years, as shown in Table 8, Figure 8 compares the scheduling processes of the GAFL1 and GAFL2 algorithms across three typical hydrological years, examining monthly energy (P), water head (h), and spillage (spl) in cascade reservoirs. It is observed that during both wet and dry years, the differences in h under the two methods are minimal, thus allowing for an increase in energy generation concurrent with reduced spillage. However, in normal years, despite a reduction in spillage, there is a significant decrease in head, which does not substantially aid in increasing energy generation. This leads to the phenomenon reported in Table 7, where a decrease in spillage during normal years does not increase energy generation, but rather causes a decrease. This discrepancy highlights the complex interdependencies between head, spillage, and energy generation under varying hydrological conditions.
Figure 9 shows the monthly scheduling processes of (a) GAFL1 and (b) GAFL2, detailing the operation of the four reservoirs along the Jinsha River during a normal year. In this context, Zmin represents the minimum water level for the year, typically set as the dead water level, while Zmax corresponds to the upper water level limit, enforced as the normal water level during the dry season and the flood control level during the flood season. By comparing the reservoir outflow (area charts of spillage and generation discharge) with the inflow (red line), the storage and release dynamics at different times can be analyzed. When the inflow exceeds the outflow, the reservoir stores water, causing the water level to rise; conversely, the reservoir releases water, causing the water level to drop. Both GAFL1 and GAFL2 methods yield reasonable results, storing water during the dry season and releasing it during the flood season, thereby minimizing spillage, maximizing hydropower generation, and ensuring optimal utilization of water resources.

5. Conclusions

This study integrates genetic algorithms (GAs) with a simulation program, specifically transferring the computation of the fitness function to the simulation in order to address issues of “dimensionality disaster” and local optima commonly encountered in the scheduling of cascade reservoirs. The objective is to maximize energy generation and ensure optimal output. Additionally, this study incorporates the concept of FL into the simulation program. This approach utilizes a push–pull strategy to explore additional spaces for harnessing and utilizing unnecessary spillage for energy generation, thus achieving optimal reservoir scheduling. It achieves this using two specific strategies: GAFL I, which includes only a forward push step in the simulation, and GAFL2, which features a sequence of forward push, backward pull, and another forward push to enhance the utilization of resources. Thus, this reveals the following:
(1)
This study verifies the feasibility and scientific basis of GAFL by comparing the data from the results of GAFL1 and Sequential Quadratic Programming (SQP). The results show that although GAFL1’s spillage increased by 2.2 times compared to SQP and total energy slightly decreased, GAFL1 significantly improved firm yield by 8.3% at the highest priority level. This highlights GAFL1’s advantage in prioritizing firm yield despite increased spillage, demonstrating its potential for optimizing reservoir operation under specific operating priorities.
(2)
Because the solution of the fitness function part is implemented within the simulation program, the workload of the genetic algorithm is reduced, thereby significantly improving the efficiency of the optimization process. Both GAFL1 and GAFL2 exhibit good convergence effects, achieving convergence of the objective function values in the 15th and 10th iterations, respectively. This indicates that both versions of the algorithm can quickly reach optimal or near-optimal solutions within a limited number of iterations, making them highly effective in addressing complex optimization problems in reservoir scheduling.
(3)
The comparative results between GAFL1 and GAFL2 show that GAFL2’s strategy comprising a push forward, a pull back, and a push forward again can minimize spillage to the greatest extent. This proves that it is effective in reducing spillage under all hydrological conditions without impacting the highest priority stable output.
(4)
The results indicate that across all hydrological years, FL can maximize the reduction in spillage without affecting the highest priority stable output. Based on this, it is observed that in both drought and wet years, reducing spillage increases energy generation. However, the impact of this reduction is variable; in normal years, although spillage decreases, total energy generation also reduces.
In summary, although the comparison between GAFL1 and SQP shows an increase in spillage with a significant improvement in firm yield, it highlights the potential of GAFL1 in prioritizing outputs; GAFL2, on the other hand, demonstrates the efficiency of the optimization algorithm in quickly achieving optimal solutions. Both strategies effectively address complex reservoir scheduling optimization issues. However, despite this, this study has limitations due to its reliance on specific simulation environments and parameter settings, which may affect the generalizability of the results and their applicability to other types of reservoir systems. Furthermore, although the research indicates that spillage can be effectively reduced under various hydrological conditions, the relationship between the increase in total energy and the reduction in spillage is not evident in normal hydrological conditions, suggesting that further research and strategy optimizations are needed for adaptations to take place regarding various hydrological changes.

Author Contributions

J.W. and H.Z. contributed to this study’s conception and design. Modeling, data collection, simulation calculation, and analysis were performed by D.B., P.W., J.L. and G.Q. The first draft of this manuscript was written by D.B. and all authors commented on the previous versions of this manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The authors declare that no funds, grants, or other support were received during the preparation of this manuscript.

Informed Consent Statement

All authors agreed with the content of this study and all of them provided explicit consent to submit this study for publication.

Data Availability Statement

Some or all data, models, or codes that support the findings of this study are available from the corresponding author upon reasonable request. Availability of data and materials is contained within the article.

Conflicts of Interest

The authors declare no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

References

  1. Shen, J.; Zhang, X.; Wang, J.; Cao, R.; Wang, S.; Zhang, J. Optimal Operation of Interprovincial Hydropower System Including Xiluodu and Local Plants in Multiple Recipient Regions. Energies 2019, 12, 144. [Google Scholar] [CrossRef]
  2. Yu, X.; Sun, H.; Wang, H.; Liu, Z.; Zhao, J.; Zhou, T.; Qin, H. Multi-Objective Sustainable Operation of the Three Gorges Cascaded Hydropower System Using Multi-Swarm Comprehensive Learning Particle Swarm Optimization. Energies 2016, 9, 438. [Google Scholar] [CrossRef]
  3. Chen, L.; McPhee, J.; Yeh, W.W.-G. A diversified multiobjective GA for optimizing reservoir rule curves. Adv. Water Resour. 2007, 30, 1082–1093. [Google Scholar] [CrossRef]
  4. Huang, D.; Shu, Y.; Ruan, J.; Hu, Y. Ultra high voltage transmission in China: Developments, current status and future prospects. Proc. IEEE 2009, 97, 555–583. [Google Scholar] [CrossRef]
  5. Cheng, C.; Yan, L.; Mirchi, A.; Madani, K. China’s booming hydropower: Systems modeling challenges and opportunities. J. Water Resour. Plan. Manag. 2017, 143, 02516002. [Google Scholar] [CrossRef]
  6. Yoo, J.-H. Maximization of hydropower generation through the application of a linear programming model. J. Hydrol. 2009, 376, 182–187. [Google Scholar] [CrossRef]
  7. Bellman, R. Dynamic programming. Science 1966, 153, 34–37. [Google Scholar] [CrossRef] [PubMed]
  8. Feng, S.; Zheng, H.; Qiao, Y.; Yang, Z.; Wang, J.; Liu, S. Weekly hydropower scheduling of cascaded reservoirs with hourly power and capacity balances. Appl. Energy 2022, 311, 118620. [Google Scholar] [CrossRef]
  9. Coello, C.A.C.; Toscano-Pulido, G.T.; Lechuga, M.S. Handling multiple objectives with particle swarm optimization. IEEE Trans. Evol. Comput. 2004, 8, 256–279. [Google Scholar] [CrossRef]
  10. Wardlaw, R.; Sharif, M. Evaluation of genetic algorithms for optimal reservoir system operation. J. Water Resour. Plan. Manag. 1999, 125, 25–33. [Google Scholar] [CrossRef]
  11. Arunkumar, R.; Jothiprakash, V. Chaotic evolutionary algorithms for multi-reservoir optimization. Water Resour. Manag. 2013, 27, 5207–5222. [Google Scholar] [CrossRef]
  12. Parvez, I.; Shen, J.; Hassan, I.; Zhang, N. Generation of Hydro Energy by Using Data Mining Algorithm for Cascaded Hydropower Plant. Energies 2021, 14, 298. [Google Scholar] [CrossRef]
  13. El-Ela, A.A.; Fetouh, T.; Bishr, M.; Saleh, R. Power systems operation using particle swarm optimization technique. Electr. Power Syst. Res. 2008, 78, 1906–1913. [Google Scholar] [CrossRef]
  14. Nicklow, J.; Reed, P.; Savic, D.; Dessalegne, T.; Harrell, L.; Chan-Hilton, A.; Karamouz, M.; Minsker, B.; Ostfeld, A.; Singh, A.; et al. State of the art for genetic algorithms and beyond in water resources planning and management. J. Water Resour. Plan. Manag. 2010, 136, 412–432. [Google Scholar] [CrossRef]
  15. Vasan, A.; Raju, K.S. Comparative analysis of simulated annealing, simulated quenching and genetic algorithms for optimal reservoir operation. Appl. Soft Comput. 2009, 9, 274–281. [Google Scholar] [CrossRef]
  16. Hossain, M.S.; El-Shafie, A. Intelligent systems in optimizing reservoir operation policy: A review. Water Resour. Manag. 2013, 27, 3387–3407. [Google Scholar] [CrossRef]
  17. Saadat, M.; Asghari, K. Reliability improved stochastic dynamic programming for reservoir operation optimization. Water Resour. Manag. 2017, 31, 1795–1807. [Google Scholar] [CrossRef]
  18. Bilal; Rani, D.; Pant, M.; Jain, S.K. Dynamic programming integrated particle swarm optimization algorithm for reservoir operation. Int. J. Syst. Assur. Eng. Manag. 2020, 11, 515–529. [Google Scholar] [CrossRef]
  19. Li, H.; Qin, T.; Wang, W.; Lei, X.; Wu, W. The improvement of particle swarm optimization: A case study of optimal operation in Goupitan Reservoir. In Proceedings of the IOP Conference Series: Earth and Environmental Science, Harbin, China, 8–10 December 2017; IOP Publishing: Bristol, UK, 2018. [Google Scholar]
  20. Zhang, R.; Zhang, S.; Wen, X.; Yue, Z.; Zhou, Y. Optimization of short-term hydropower scheduling with dynamic reservoir capacity based on improved genetic algorithm and parallel computing. J. Hydrol. 2024, 636, 131238. [Google Scholar] [CrossRef]
  21. Wang, Z.; Wong, Y.; Rahman, M. Development of a parallel optimization method based on genetic simulated annealing algorithm. Parallel Comput. 2005, 31, 839–857. [Google Scholar] [CrossRef]
  22. Reddy, M.J.; Kumar, D.N. Optimal reservoir operation using multi-objective evolutionary algorithm. Water Resour. Manag. 2006, 20, 861–878. [Google Scholar] [CrossRef]
  23. Li, C.; Xu, B.; Gao, H.; Yin, X.; Xu, L. Method of Reservoir Optimal Operation Based on Improved Simulated Annealing Genetic Algorithm. Sens. Transducers 2013, 159, 160. [Google Scholar]
  24. Li, X.-G.; Wei, X. An Improved Genetic Algorithm-Simulated Annealing Hybrid Algorithm for the Optimization of Multiple Reservoirs. Water Resour. Manag. 2008, 22, 1031–1049. [Google Scholar] [CrossRef]
  25. Fayaed, S.S.; El-Shafie, A.; Jaafar, O. Integrated artificial neural network (ANN) and stochastic dynamic programming (SDP) model for optimal release policy. Water Resour. Manag. 2013, 27, 3679–3696. [Google Scholar] [CrossRef]
  26. Yang, K.; Zheng, J.; Yang, M.; Zhou, R.; Liu, G. Adaptive Genetic Algorithm for Daily Optimal Operation of Cascade Reservoirs and its Improvement Strategy. Water Resour. Manag. 2013, 27, 4209–4235. [Google Scholar] [CrossRef]
  27. Liu, S.; Wang, P.; Xu, Z.; Feng, Z.; Zhang, C.; Wang, J.; Chen, C. Hydropower Unit Commitment Using a Genetic Algorithm with Dynamic Programming. Energies 2023, 16, 5842. [Google Scholar] [CrossRef]
  28. Xu, B.; Zhong, P.-A.; Wan, X.; Zhang, W.; Chen, X. Dynamic Feasible Region Genetic Algorithm for Optimal Operation of a Multi-Reservoir System. Energies 2012, 5, 2894–2910. [Google Scholar] [CrossRef]
  29. Wang, J.; Chen, C.; Liu, S. A new field-levelling procedure to minimize spillages in hydropower reservoir operation. Energy 2018, 160, 979–985. [Google Scholar] [CrossRef]
  30. Liu, S.; Luo, J.; Chen, H.; Wang, Y.; Li, X.; Zhang, J.; Wang, J. Third-Monthly Hydropower Scheduling of Cascaded Reservoirs Using Successive Quadratic Programming in Trust Corridor. Water 2023, 15, 716. [Google Scholar] [CrossRef]
  31. Rashid, M.U.; Abid, I.; Latif, A. Optimization of hydropower and related benefits through Cascade Reservoirs for sustainable economic growth. Renew. Energy 2022, 185, 241–254. [Google Scholar] [CrossRef]
Figure 1. Capacity of hydropower output due to water head.
Figure 1. Capacity of hydropower output due to water head.
Energies 17 03756 g001
Figure 2. Flow chart of the genetic algorithm.
Figure 2. Flow chart of the genetic algorithm.
Energies 17 03756 g002
Figure 3. Forward simulation in a time-step.
Figure 3. Forward simulation in a time-step.
Energies 17 03756 g003
Figure 4. Forebay relationship fitting curves.
Figure 4. Forebay relationship fitting curves.
Energies 17 03756 g004
Figure 5. Tailwater relationship fitting curves.
Figure 5. Tailwater relationship fitting curves.
Energies 17 03756 g005
Figure 6. Inflow trend at Wudongde during a hydrologic year.
Figure 6. Inflow trend at Wudongde during a hydrologic year.
Energies 17 03756 g006
Figure 7. Iteration process of the GAFL2 optimization.
Figure 7. Iteration process of the GAFL2 optimization.
Energies 17 03756 g007
Figure 8. Comparing spillage and water heads between GAFL1 and GAFL2.
Figure 8. Comparing spillage and water heads between GAFL1 and GAFL2.
Energies 17 03756 g008
Figure 9. The scheduling process for the four reservoirs derived using the two methods.
Figure 9. The scheduling process for the four reservoirs derived using the two methods.
Energies 17 03756 g009aEnergies 17 03756 g009b
Table 1. Basic characteristics of cascaded hydropower reservoirs.
Table 1. Basic characteristics of cascaded hydropower reservoirs.
NameAnnual
Inflow
(m3/s)
Installed
Capacity
(MW)
Storage
Capacity
(108 m3)
Minimum
Storage
(108 m3)
Minimum
Release
(m3/s)
Maximum
Release
(m3/s)
Operability
Wudongde344210,20074.0833.0490037,000Seasonal
Baihetan364716,00020685.6960037,000Annual
Xiluodu464012,600126.751.1280037,000Annual
Xiangjiaba4255640051.6340.74120037,000Seasonal
Table 2. Capacity of hydropower output of each hydropower plant.
Table 2. Capacity of hydropower output of each hydropower plant.
Name c i ( 0 ) d i ( 0 ) c i ( 1 ) d i ( 1 )
Wudongde147.3−10,723010,200
Baihetan179.1−19,518016,000
Xiluodu83.523−3746.1012,600
Xiangjiaba99.7154498.506400
Table 3. Parameters estimated for the forebay relationship curve.
Table 3. Parameters estimated for the forebay relationship curve.
Name α i V ^ i 0 β i Z ^ i u R2
Wudongde2.27947825.604850.789139939.114830.99995
Baihetan23,475.4654.039110.00446923,234.350.99919
Xiluodu16.755110.607250.486809438.434231
Xiangjiaba65.04197.6207280.287402171.703110.99993
Table 4. Parameters estimated for the tailwater relationship curve.
Table 4. Parameters estimated for the tailwater relationship curve.
Name χ i Q ^ i 0 δ i Z ^ i d R2
Wudongde0.126124456.6640.565465808.66720.99961
Baihetan0.087444432.38180.598788582.32390.99978
Xiluodu0.0930965025.770.583048365.23060.99977
Xiangjiaba0.0781171051.570.566308261.41790.9998
Table 5. Comparison of results between GAFL1 and SQP (1983).
Table 5. Comparison of results between GAFL1 and SQP (1983).
ModelFirm Yield
(MW)
Total Energy
(106 kWh)
Spillage
(hm3)
GAFL19769.25152,951.623550.03
SQP9020.0154,778.41607.7
Improvement8.3%−1.2%2.2 times
Table 6. Comparison of the result of each hydropower plant.
Table 6. Comparison of the result of each hydropower plant.
NAMEReservoir TypeGAFL1SQP
Total Energy
(106 kWh)
Spillage
(hm3)
Total Energy
(106 kWh)
Spillage
(hm3)
WudongdeSeasonal27,962.70.030,458.00.0
BaihetanAnnual45,131.22426.347,798.30.0
XiluoduAnnual47,101.61123.848,784.8323.4
XiangjiabaSeasonal32,756.10.027,737.31284.4
Table 7. GAFL function setting attributes.
Table 7. GAFL function setting attributes.
PerformancesPopulation SizeIteration NumberCPU TimeCrossover FractionMigration Fraction
GAFL1500100218.077 s0.890.89
GAFL2500100639.010 s0.890.89
Table 8. Comparison without vs. with minimizing spillages.
Table 8. Comparison without vs. with minimizing spillages.
PerformancesWet (2000)Normal (1978)Dry (1976)
GAFL1GAFL2GAFL1GAFL2GAFL1GAFL2
Firm yield
(MW)
15,743.815,743.811,596.111,596.110,388.610,388.6
Total energy (106 kWh)306,753.7307,081↑188,385.1186,975.2↓163,464.2163,464.6↑
Spillage
(hm3)
29,674.529,298.2↓7743.76899.5↓1417.41414.5↓
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Baima, D.; Qian, G.; Luo, J.; Wang, P.; Zheng, H.; Wang, J. Monthly Hydropower Scheduling of Cascaded Reservoirs Using a Genetic Algorithm with a Simulation Procedure. Energies 2024, 17, 3756. https://doi.org/10.3390/en17153756

AMA Style

Baima D, Qian G, Luo J, Wang P, Zheng H, Wang J. Monthly Hydropower Scheduling of Cascaded Reservoirs Using a Genetic Algorithm with a Simulation Procedure. Energies. 2024; 17(15):3756. https://doi.org/10.3390/en17153756

Chicago/Turabian Style

Baima, Deji, Guoyuan Qian, Jingzhen Luo, Pengcheng Wang, Hao Zheng, and Jinwen Wang. 2024. "Monthly Hydropower Scheduling of Cascaded Reservoirs Using a Genetic Algorithm with a Simulation Procedure" Energies 17, no. 15: 3756. https://doi.org/10.3390/en17153756

APA Style

Baima, D., Qian, G., Luo, J., Wang, P., Zheng, H., & Wang, J. (2024). Monthly Hydropower Scheduling of Cascaded Reservoirs Using a Genetic Algorithm with a Simulation Procedure. Energies, 17(15), 3756. https://doi.org/10.3390/en17153756

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop