Journal cover Journal topic
Proceedings of the International Association of Hydrological Sciences An open-access publication for refereed proceedings in hydrology
Journal topic
Proc. IAHS, 379, 199-203, 2018
https://doi.org/10.5194/piahs-379-199-2018
Proc. IAHS, 379, 199-203, 2018
https://doi.org/10.5194/piahs-379-199-2018

Pre-conference publication 05 Jun 2018

Pre-conference publication | 05 Jun 2018

# Mid and long-term optimize scheduling of cascade hydro-power stations based on modified GA-POA method

Mid and long-term optimize scheduling of cascade hydro-power stations
Jiqing Li and Xiong Yang Jiqing Li and Xiong Yang
• School of Renewable Energy, North China Electric Power University, Beijing 102206, China
Abstract

In this paper, to explore the efficiency and rationality of the cascade combined generation, a cascade combined optimal model with the maximum generating capacity is established, and solving the model by the modified GA-POA method. It provides a useful reference for the joint development of cascade hydro-power stations in large river basins. The typical annual runoff data are selected to calculate the difference between the calculated results under different representative years. The results show that the cascade operation of cascaded hydro-power stations can significantly increase the overall power generation of cascade and ease the flood risk caused by concentration of flood season.

1 Introduction

As more and more large and medium hydro-power stations are put into use in our country, the joint dispatching of cascaded hydro-power plants has been getting more and more widespread attention. In order to encourage all hydro-power stations in cascade hydro-power stations to actively participate in the joint dispatch in order to maximize the benefits, the hydro-power industry advocates the implementation of “Cascade, rolling, integrated” mode of operation (Xu, 2015; He, 2014; Zhou, 2016). In the joint operation of cascade hydro-power stations, the leading hydro-power station in the upper reaches has good regulating performance, which can greatly increase the power generation capacity of each downstream hydro-power station and has objective power generation compensation benefits (Li, 2016; Wang, 2014). So, reflecting the compensating function of power plants with good performance in cascade joint scheduling and solving the joint operation process by highly efficient and feasible algorithms mean very much.

2 Modified GA-POA method

According to the idea of successive optimization, the T-stage optimization problem can be divided into several two-stage optimization sub-problems, while the sub-problem decision-making, using real-coded genetic algorithm to solve. Traditional genetic algorithms generate initial populations in a random way. However, randomized methods often fail to meet the requirements of “global distribution” of initial populations, which leads to their low representation and inefficiency. Based on this, the idea of uniform design is introduced into the generation of initial population of genetic algorithm to improve the initial population quality by increasing the uniform distribution of the initial population.

## 2.1 Coding method and initial population generation

Taking the water level of the library as a gene and coding with an integer, the population size is set as Size, $i=\mathrm{1}\sim$ Size. The initial population is generated as follows:

$\begin{array}{}\text{(1)}& & {P}_{i,t}={Z}_{t,min}+\mathrm{0.1}n\text{(2)}& & n=\mathrm{10}×\left({Z}_{t,max}-{Z}_{t,min}\right)=\mathrm{1},\mathrm{2},\phantom{\rule{0.125em}{0ex}}\mathrm{\dots }\phantom{\rule{0.125em}{0ex}},\mathrm{size},\end{array}$

where Pi,t is the integer coding of the i individual in the parent population at time t, $t=\mathrm{1},\mathrm{2},\mathrm{\dots },T+\mathrm{1}$; T is the number of scheduling period; ${Z}_{t,min},{Z}_{t,max}$ is upper and lower water level.

## 2.2 Fitness function setting

The fitness function is used to evaluate the individual's advantages and disadvantages. In this paper, the size of power generation is mainly considered. When the system initial solution does not meet the minimum system output requirements, the penalty function is used to adjust. The objective function value:

$\begin{array}{}\text{(3)}& \mathit{\varphi }\left(Z\right)=\sum _{t=\mathrm{1}}^{T}{N}_{t}\cdot \mathrm{\Delta }t+\mathit{\alpha }\cdot min\left(\frac{{N}_{t}-{N}_{min}}{{N}_{min}},\mathrm{0}\right),\end{array}$

where T is Each time period; Nt is the output of the cascade at time t; Nmin is system minimum output.

## 2.3 Cross operation

Using the single point crossover method, it is assumed that i1 and i2 cross at the time point:

$\begin{array}{}\text{(4)}& {P}_{k,t}^{\prime }=\left\{\begin{array}{l}{p}_{i\mathrm{1},t}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}t<\mathrm{time}\\ {p}_{i\mathrm{2},t}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}t\ge \mathrm{time}\end{array}\right\{P}_{k+\mathrm{1},t}^{\prime }=\left\{\begin{array}{l}{p}_{i\mathrm{2},t}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}t<\mathrm{time}\\ {p}_{i\mathrm{1},t}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}t\ge \mathrm{time},\end{array}\right\\end{array}$

where time $\in \left[\mathrm{1},\mathrm{2},\phantom{\rule{0.125em}{0ex}}\mathrm{\dots }\phantom{\rule{0.125em}{0ex}},T+\mathrm{1}\right]$, after crossing, a population P of the same size as Size is formed.

## 2.4 Variation operation

Using random mutation, the gene Py is controlled by the probability Py, a new gene is generated at the mutation point:

$\begin{array}{}\text{(5)}& {P}_{i,t}^{\prime \prime }=\left\{\begin{array}{l}{Z}_{t,min}+\left({Z}_{t,max}-{Z}_{t,min}\right)\cdot {P}_{y}^{\prime }{P}_{y}^{\prime }<{P}_{y}\\ {P}_{i,t}{P}_{y}^{\prime }\ge {P}_{y},\end{array}\right\\end{array}$

where Py, ${P}_{y}^{\prime }$ are all random numbers distributed between [0, 1], mutation to population ${P}^{\prime \prime }$ with the same population size of size.

## 2.5 Choose operation

The league method (Howson, 1975) was used to score a total of 3Size individuals in the parental group P, crossover groupP, and mutation group P. The scoring rules of individual i were as follows: size random individuals were randomly selected without duplication, and the adaptability of i exceeded that of the competition individuals counts Score as i score:

$\begin{array}{}\text{(6)}& {\mathrm{Score}}_{i}=\mathrm{Count}\mathit{\left\{}{\mathit{\varphi }}_{i}>{\mathit{\varphi }}_{j}\left|i,j\in {\mathrm{\Omega }}_{\mathrm{size}}\right\mathit{\right\}}.\end{array}$

3Size individuals ranked according to the score, the highest score of Size individuals as a parent group into the next generation of evolution.

## 2.6 Evolution termination conditions

The evolution termination condition is to reach the maximum evolutionary algebra set in advance.

The multi-stage optimization problem in the progressive optimality algorithm is split into multiple two-stage sub-problems, and then the modified genetic algorithm is applied to the decision of the sub-problems. The modified GA-POA method solving flowchart is shown as Fig. 1.

Figure 1The modified GA-POA method solving flowchart.

3 Mid and Long-Term Optimize Scheduling of Cascade Hydro-power Stations

## 3.1 The objective function

Model to maximize the power generation as the optimization goal, the mathematical description of:

$\begin{array}{}\text{(7)}& maxF=\sum _{i=\mathrm{1}}^{n}\sum _{j=\mathrm{1}}^{T}{N}_{ij}\mathrm{\Delta }t,\end{array}$

where F is the goal of optimal scheduling is the total power generation of all stations in the cascade throughout the entire scheduling period; n is the number of reservoirs; T is scheduling period; Nij is the i reservoir's output of j period; Δt is Time interval.

## 3.2 Restrictions

### 3.2.1 Water balance constraints

$\begin{array}{}\text{(8)}& {V}_{t+\mathrm{1}}-{V}_{t}=\left({Q}_{t}-{q}_{t}\right)\cdot \mathrm{\Delta }t,\end{array}$

where ${V}_{t},{V}_{t+\mathrm{1}}$ is the beginning and the end of the reservoir's storage capacity of t period;

Qt, qt is the average inflow and the average outflow.

### 3.2.2 Water level constraints

$\begin{array}{}\text{(9)}& {Z}_{min}\phantom{\rule{0.125em}{0ex}}\mathit{<}\phantom{\rule{0.125em}{0ex}}{Z}_{t}\phantom{\rule{0.125em}{0ex}}\mathit{<}\phantom{\rule{0.125em}{0ex}}{Z}_{max},\end{array}$

where Zmax⁡, Zmin⁡ is the reservoir's upper and lower water level of t period.

### 3.2.3 Output constraint

$\begin{array}{}\text{(10)}& {N}_{min}\phantom{\rule{0.125em}{0ex}}\mathit{<}\phantom{\rule{0.125em}{0ex}}{N}_{t}\phantom{\rule{0.125em}{0ex}}\mathit{<}\phantom{\rule{0.125em}{0ex}}{N}_{max},\end{array}$

where Nmax⁡, Nmin⁡ is the reservoir's upper and lower output.

### 3.2.4 Discharge flow constraints

$\begin{array}{}\text{(11)}& {q}_{min}\phantom{\rule{0.125em}{0ex}}\mathit{<}\phantom{\rule{0.125em}{0ex}}{q}_{t}\phantom{\rule{0.125em}{0ex}}\mathit{<}\phantom{\rule{0.125em}{0ex}}{q}_{max},\end{array}$

where qmax⁡, qmin is the reservoir's upper and lower outflow.

### 3.2.5 Hydraulic equation

$\begin{array}{}\text{(12)}& {I}_{i,l}={q}_{i-\mathrm{1},l}+{Q}_{qj,l}-{Q}_{ks,l},\end{array}$

where Ii,l is the reservoir is inflow of the l period; ${q}_{i-\mathrm{1},l}$ is the reservoir i−1's outflow of the l period; Qqj,l, Qks,l is the interval flow and the loss flow of the l period.

### 3.2.6 Non-negative conditional constraints

All of the above variables satisfy the non-negative condition.

## 3.3 Output and power generation

$\begin{array}{}\text{(13)}& & N=KQ\cdot \left(H-\mathrm{\Delta }H\right)\text{(14)}& & \mathrm{\Delta }H=\mathit{\theta }\cdot {Q}^{\mathrm{2}}\text{(15)}& & E=N\cdot t,\end{array}$

where K is output coefficient; Q is used flow; H is net head; ΔH is head loss; θ is coefficient of head loss; t is time; N,E is the output and power generation of hydro-power station.

4 Case study

In this study, an example is given to four hydro-power stations in a river basin cascade in China, which are represented by A, B, C and D. Among them, A is the leading power plant with huge storage capacity and significant power generation benefit which has not been completed and put into production. The adjustment capacity of B, C and D is small, and the runoff data used in this study is the monthly typical runoff which are regarded as run-of-run power plants. In this study, population size = 200, maximum iteration generation = 300, crossover time = 3 and mutation probability = 0.03. Solving the model by the established modified GA-POA method with C$++$ programming language, the cascade power generation under different representative years are shown as Table 1 and Fig. 2.

Table 1Cascade power generation 108 kW h.

When cascade hydro-power station is operated in combination, the leading hydro-power station can effectively increase the power generation of the downstream hydro-power plants during the dry period, and the total generating capacity of the hydro-power station increases significantly. The power generation of joint operation growth by 4.46 % in flood year, 2.25 % in mid-year and 1.76 % in dry year than alone operation. The compensatory benefits meet the requirements of Flood Year > Mid-Year > Dry Year.

The joint operation can significantly improve the overall power generation of the cascade. When the water inflow is small, the joint operation has a greater significant. The upstream power station compensates the power plants of the downstream power stations and improves the output of the dry season and the seasonal power generation. It can also ease the flood risk caused by concentration of flood season. This will help mobilize the enthusiasm of the leading hydro-power station and downstream hydro-power stations to participate in the joint cascade scheduling and achieve the goal of joint operation, which bring the maximum benefit in the basin.

The monthly output process of each station in flood year when it is operated independently and jointly is shown as Figs. 3 to 6. When each station of cascades operates independently, the output of the power plant fluctuates greatly. The specific performance is that when in the flood season, the output is large; when in the dry season, the output is small.

Figure 3The output of plant A.

Figure 4The output of plant B.

Figure 5The output of plant C.

Figure 6The output of plant D.

When the water supply is reduced rapidly in some months, there are efforts to make a sharp turn; when the cascade power plants are jointly operated, the power plant output process is relatively gradual. Due to the regulation and storage of the leading power station, the output of the power plant during the flood season is smaller than that of the power station alone. In the dry season, the power plant output is greater than when the power station is operated alone. Output, the overall change is smaller than that of each power station when it is running alone, and the leading power station has a significant effect on the adjustment of the steps.

5 Conclusion

In this paper, the long-term optimal scheduling model is established, and the modified GA-POA method is proposed to solve the model by C$++$ programming. This model is applied to cascade reservoirs in a river basin in China, which achieved good results. The study provides a beneficial reference for the joint development of large-scale cascaded hydro-power stations.

Data availability
Data availability.

The data used in this paper are provided by the reservoir management departments, which can not be provided externally because the data are confidential.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

This article is part of the special issue “Innovative water resources management – understanding and balancing interactions between humankind and nature”. It is a result of the 8th International Water Resources Management Conference of ICWRS, Beijing, China, 13–15 June 2018.

Acknowledgements
Acknowledgements.

This study was financially supported by National Key Projects of China “Water resources efficient development and utilization” (2017YFC0405900, 2016YFC0402208, 2016YFC0401903) and National Natural Science Foundation of China (No.51641901). The author would like to give special thanks to anonymous reviewers.

Edited by: Yangbo Chen
Reviewed by: two anonymous referees

References

He, X., Wu, Y., and Jiang, G.: Profit allocation model of supply network based on revised Shapely value, Soft Science, 28, 70–73, 2014.

Howson, H. R. and Sancho, N. G. F.: A new algorithm for the solution of multi-state dynamic programming problems, Math. Program., 8, 104–106, 1975.

Li, Y., Zhou, J., Zhang, S., Huang, C., and Yan, D.: Combined dispatching chart of cascade hydropower stations on Yalong river, Jinsha river and Yangtze river, Journal of Hydroelectric Engineering, 35, 32–40, 2016.

Wang, Y., Guo, S., Zhou, Y., Hong, X., Hu, T., and Liu, D.: Uncertain analysis of long-term generation optimal operation for reservoirs, Water Resources and Power, 32, 61–65, 2014.

Xu, W.: Research on some issues in the joint optimal operation of cascaded reservoirs between Long yang xia and Liu jia xia, Tsinghua University, 2015.

Zhou, B.: Research on compensative operation of two-reservoirs and multi-cascaded hydropower plants, Da Lian University of Technology, 2016.