【水库调度】01 差分进化算法（DE）求解单库发电优化问题

置顶卢家波

已于 2024-04-27 17:48:27 修改

阅读量2.9k

点赞数 4

分类专栏： # 水库调度文章标签：水库调度差分进化算法 DE 发电

于 2019-11-20 16:08:46 首次发布

本文链接：https://blog.csdn.net/weixin_43012724/article/details/103152533

版权

水库调度专栏收录该内容

2 篇文章 5 订阅

订阅专栏

算法简介

差分进化算法是进化算法的一种，与遗传算法（GA）、粒子群算法（PSO）、蚁群算法、狼群算法等群体智能算法一样，常用来求解复杂优化问题的全局最优解。
具体算法参看差分进化算法
采用MATLAB编程，代码参考算法代码

博客作者简介

很高兴认识您！
我叫卢家波，河海大学水文学及水资源博士研究生，研究兴趣为高效洪水淹没预测、洪水灾害预警、机器学习、替代模型和降阶模型。
变化环境下，极端洪水事件多发，我希望能通过研究为水灾害防御做出贡献，为人民服务。
欢迎交流讨论和研究合作，vx Jiabo_Lu。
主页 https://lujiabo98.github.io
简历 https://lujiabo98.github.io/file/CV_JiaboLu_zh.pdf
博客 https://blog.csdn.net/weixin_43012724?type=blog
来信请说明博客标题及链接，谢谢。

测试函数

$f(x,y)=x^2+y^2$

$max f(x,y)=x^2+y^2$
$x,y\in[-1,1],x^2+y^2\le1$
在这里插入图片描述

$f(x,y,z)=x^2+y^2+z^2$

$maxf(x,y,z)=x^2+y^2+z^2$
$x,y,z\in[-1,1],x^2+y^2+z^2\le1$
在这里插入图片描述

$R a s t r i g in 函数$

$minf(x)=10d+\sum_{i=1}^d[x_i^2-10cos(2\pi x_{i})]$
$x_i\in[-5.12,5.12],i=1,...,d$
在这里插入图片描述

$S c ha ff er 函数$

$f(x)=0.5+\frac{sin^2(x_1^2+x_2^2)-0.5}{[1+0.001(x_1^2+x_2^2)]^2}$
$x_i\in[-100,100],i=1,2$
在这里插入图片描述

$A c k l ey 函数$

$f(x)=-aexp(-b\sqrt{\frac{1}{d}\sum_{i=1}^dx_i^2})-exp(\frac{1}{d}\sum_{i=1}^dcos(cx_i))+a+e$
$a=20,b=0.2,c=2\pi,x_i\in[-32.768,32.768],i=1,...,d$
在这里插入图片描述

$G r i e w ank 函数$

$f(x)=\sum_{i=1}^d\frac{x_i^2}{4000}-\prod_{i=1}^{d}cos(\frac{x_i}{\sqrt{i}})+1$
$x_i \in [-600, 600], i = 1, …, d$
在这里插入图片描述

单库发电优化建模

目标函数

以调度期内发电量最大为优化调度目标，其目标函数为：
$maxE=\sum_{t=1}^{T}N_t(qe_t,H_t)\Delta t$
式中， $N_t(·)$ 为 $t$ 时段的电站平均出力， $H_t$ 为水库 $t$ 时段平均水头， $qe_t$ 为 $t$ 时段平均发电流量， $T$ 为调度期时段数， $\Delta t$ 为时段长。

约束条件

水量平衡约束
$V_t=V_{t-1}+(Q_t-q_t)\Delta t$
式中， $V_t、V_{t +1}$ 分别为水库 $t$ 时段始末库蓄水量， $Q_t$ 为 $t$ 时段平均入库流量， $q_t$ 为 $t$ 时段平均出库流量。
水位约束
$\underline{Z_t}\le Z_t \le \overline{Z_t}$
式中， $\overline{Z_t}、\underline {Z_t}$ 分别为 $t$ 时段初水库水位上下限。
出库流量约束
$\underline{q_t}\le q_t \le \overline{q_t}$
式中， $\overline{q_t}、\underline {q_t}$ 分别为 $t$ 时段初水库出库流量上下限。
出力约束
$\underline{N_t}\le N_t \le \overline{N_t}$
式中， $\overline{N_t}、\underline {N_t}$ 分别为 $t$ 时段初水库出力上下限。
调度期初、末水位约束
$Z_0=Z_s$
$Z_T=Z_e$
式中： $Z_s$ 为调度期初水库水位， $m$ ； $Z_e$ 为调度期末水库控制水位， $m$ 。

计算实例

基本资料

以三峡水库的发电优化调度问题为例，其兴利库容165 亿 $m^3$ ，库容系数0.04，具有不完全年调节能力。水库正常蓄水位175 m，防洪限制水位145 m，枯季消落低水位155 m。电站32台机组，总装机容量2250万kW，保证出力499万kW，下游航运及生态基流要求5 600 $m^3/s$ 。选择调度期为1年，计算时段为旬，调度期初水位为174 m，调度期末水位也控制为173 m，调度期内的逐旬入流过程为模拟入库过程。

模型结构

以调度期内发电量最大为目标函数，最小出库流量和保证出力作为罚函数，最大出力和水轮机组最大过水流量作为阈值，约束条件如上。
$E=\sum_{t=1}^{T}{N_t(qe_t,H_t)\Delta t + Inf*min(\frac{N_t-\underline {N_t}}{\underline {N_t}},0)+Inf*min(\frac{qe_t-\underline {q_t}}{\underline {q_t}},0)}$
式中， $N_t(·)$ 为 $t$ 时段的电站平均出力， $H_t$ 为水库 $t$ 时段平均水头， $qe_t$ 为 $t$ 时段平均发电流量， $T$ 为调度期时段数， $\Delta t$ 为时段长， $I n f$ 为罚因子，这里取一个很大的正数。

计算结果

在这里插入图片描述

MATLAB代码