一种针对多年调节水库的智能调度方法

1.本发明涉及水库调度领域，尤其涉及一种针对多年调节水库的智能调度方法。


背景技术：

2.自上世纪五十年代以来，由于动态规划原理和马尔科夫决策方法的提出，同时随着电子计算机的兴起，水库中长期优化调度研究得到了迅速发展，并取得丰硕成果。
3.水库优化调度问题主要通过建立水库调度的数学模型进行，根据输入径流的描述方法不同，水库优化调度可分为随机优化调度和确定性优化调度。随机优化调度又可分为显随机优化和隐随机优化调度。显随机优化调度是指将径流过程描述成符合一定概率分布的不确定性条件，进而运用确定性优化原理实现长系列优化，典型代表方法有随机动态规划；隐随机优化调度是指以确定性优化调度为样本，从样本分析中提取具有指导意义的调度规则，典型代表方法为智能算法、多元线性回归方法等。
4.本技术发明人在实施本发明的过程中，发现现有技术中存在如下技术问题：
5.随机优化调度输入资料即径流量是一个不完全确定量，也为最优策略以及最优决策选取带来很大的不确定性。确定性优化调度研究起始于单库优化调度问题，模型结构比较简单，它将长系列径流过程视为已知条件，所求最优效益是最大可能的极限值，典型代表方法有运筹学理论里的线性规划、动态规划类方法以及智能算法等。
6.由此可知，现有技术中的方法存在因预报不确定性而导致的调度效果不佳的技术问题。


技术实现要素：

7.本发明提出一种针对多年调节水库的智能调度方法，用于解决或者至少部分解决现有技术中的方法存在因预报不确定性而导致的调度效果不佳的技术问题。
8.为了解决上述技术问题，本发明提供了一种针对多年调节水库的智能调度方法，包括：
9.s1：将决策过程分为预见期末、预见期末所在月份的月末和年末三个时间节点，选取发电量最大作为目标函数，根据水库运行实际建立约束条件，其中约束条件包括水量平衡约束、水电站物理特性约束；
10.s2：根据历史资料采用变动流量方式计算各离散点位年末余留效益，其中，变动流量过程计算方式为：以末水位等于初水位为前提，根据不同年份的历史月径流资料进行计算最后取平均，得到的结果作为每一年年末各点的余留效益；
11.s3：根据预见期末所在月份的月末至年末的预报月径流资料，进行以月为时段的优化调度计算，获得该月末各离散点位的余留效益；
12.s4：再根据预见期末至该月末预报日流量，进行以日为时段的优化调度计算，得到预见期末各离散点位的余留效益；
13.s5：以当前时刻水位为初始值，以预见期末各离散点位的余留效益为边界条件，在
预见期内进行以小时为时段的水电站短期经济运行，得到最优决策过程以及相应的最优效益，基于所述最优决策过程对水电站完成一日的调度决策。
14.在一种实施方式中，所述方法还包括：
15.s6：完成一日的调度决策后，再进行下一日开始的调度，首先比较预见期末时间是否超出本月，若超出，则以下一个月作为预见期末执行步骤s3～s5，若未超出，则执行步骤s4～s5。
16.在一种实施方式中，步骤s1中的目标函数为：
[0017][0018]
式中：为计划期t内水电站总发电收入，元；p(t)为t时段水电站总出力，kw；c(t)为t时段电价，元/(kwh)，可反映峰谷电价差异，若c(t)＝1，目标函数就表示发电量最大；
[0019]
水量平衡约束包括水量平衡方程：
[0020]
v(t+1)＝v(t)+[q
rk
(t)
‑
q
fd
(t)
‑
q
qs
(t)]
×
δt(t)
[0021]
v(t+1)、v(t)代表t+1时刻与t时刻的库容，δt(t)代表时段长，q
rk
(t)、q
fd
(t)、q
qs
(t)分别代表水库t时刻的入库、发电、弃水流量；
[0022]
水电站物理特性约束包括库水位或库容约束、水位库容曲线约束、尾水位下泄曲线约束、闸门泄流曲线约束和水头出力限制曲线限制；
[0023]
(1)库水位或库容约束：
[0024][0025]
z
sy
(t+1)为t+1时刻的上游水位，v(t+1)为t+1时刻的库容，为t+1时刻的上游水位下限，通常为死水位，为t+1时刻的上游水位上限，通常在汛期为汛限水位，非汛期为正常蓄水位，v
min
(t+1)与v
max
(t+1)分别为上游水位下限和上游水位上限所对应的库容；
[0026]
(2)水位库容曲线约束：
[0027]
z
sy
(t+1)＝f
zv
[v(t+1)]
[0028]
z
sy
(t+1)为t+1时刻的上游水位，v(t+1)为t+1时刻的库容，f
zv
为水位库容曲线关系；
[0029]
(3)尾水位下泄曲线约束：
[0030]
z
xy
(t)＝f
zq
[q
ck
(t)]
[0031]
z
xy
(t)为t时刻的下游水位，q
ck
(t)为t时刻的出库流量，f
zq
为尾水位下泄流量曲线关系；
[0032]
(4)闸门泄流曲线约束：
[0033]
q
qs
(t)≤f
qz
[z
xy
(t)]
[0034]
q
qs
(t)为t时刻的弃水流量，z
xy
(t)为t时刻的下游水位，f
qz
为水位与闸门泄流能力曲线关系；
[0035]
(5)水头出力限制曲线限制；
[0036]
p(t)≤f
ph
[h0(t)]
[0037]
p(t)为t时刻的出力，h0(t)为t时刻的水头，f
ph
为水头与出力曲线关系。
[0038]
在一种实施方式中，步骤s2包括：
[0039]
获取不同年份的历史月径流资料，具体包括n年历史月径流数据；
[0040]
根据每年的历史月径流资料进行一年内的优化调度计算，则年末余留效益计算公式为：
[0041][0042]
其中，表示第n年所得的优化调度最优发电量，为计算期初末水位均为xm时所对应的余留效益，表示采用第t年径流资料进行优化调度计算所得到的最优发电量。
[0043]
在一种实施方式中，步骤s3包括：
[0044]
在预见期末到年末这一时段内以月为单位采用离散微分动态规划算法，计算到年末离散出的每一末水位点所对应的发电量，得到的视作阶段效益其中n指对应的第n个水位点，则该初水位到年末各个水位离散点总效益为：
[0045][0046]
其中代表初水位为xm时到年末第n个离散水位的总效益，表示年末第n个点对应的余留效益，代表初水位为xm时到年末第n个离散水位的阶段效益；
[0047]
从初水位为xm时到年末各个离散点对应的总效益中，筛选总效益的最大值作为初水位为xm时对应的余留效益：
[0048][0049]
代表水位为xm时的月末余留效益，代表初水位为xm时到年末各个离散水位点的总效益。
[0050]
本技术实施例中的上述一个或多个技术方案，至少具有如下一种或多种技术效果：
[0051]
1.可在计算过程中得到预见期末、月末、年末的各离散点位余留效益，并将其与预见期阶段效益一起确定最优策略，使所得的最优策略更为科学合理；
[0052]
2.考虑历史长系列径流、中长期径流预报和短期径流预报信息，得到水电站的最优决策过程以及对应的最优效益更符合实际情况；
[0053]
3.调度决策可随径流预报的滚动更新计算，可极大减少人工干预，实时性强。
附图说明
[0054]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现
有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0055]
图1为本发明实施例中变动流量过程计算示意图；
[0056]
图2为本发明实施例中水库蓄水状态(水位)离散示意图；
[0057]
图3为本发明实施例中离散点位示意图；
[0058]
图4为本发明实施例中智能优化调度模型三段计算方式；
[0059]
图5为本发明实施例中循环计算框图。
具体实施方式
[0060]
针对现有方法调度效果不佳的问题，本发明提出了一种既考虑不同预见期入库流量预报精度，又考虑末水位不定的情况下，通过对历史径流过程的模拟计算确定年末不同末水位时的余留效益，再结合入库流量的滚动预报及不同预见期的预报精度情况，滚动寻优确定各时段最优决策的新型智能调度模式，本发明的目的在于提供一种实用合理的多年调节水库智能优化调度方式。
[0061]
本发明的主要构思为：
[0062]
一、按照精度需要对水库正常运行情况下的蓄水位进行离散；二、根据水库历史资料计算离散各水位点所对应的年末余留效益；三、根据预见期末所在月份至年末预报月径流资料计算离散各水位点所对应的月末余留效益；四、根据预见期末至月末预报日径流资料计算离散各水位点所对应的预见期末余留效益；五、在预见期内进行短期经济运行得到逐小时的水电站最优决策过程以及相应的最优效益；六、完成一日的调度决策后，再进行下一日开始的调度，首先比较预见期末时间是否超出该月，若超出，则计算下一月各离散水位点所对应余留效益，接着重复上述步骤进行下一次计算，若未超出，则直接根据本月余留效益进行下一次计算。
[0063]
本发明所提供的针对多年调节水库的智能调度方法，既具有实时性又极大减少了人工干预，且该方法概念清晰，考虑全面，对充分发挥水电站水库的经济效益、提高水电站水库的运行调度水平具有非常重要的现实意义。
[0064]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0065]
本发明实施例提供了一种针对多年调节水库的智能调度方法，包括：
[0066]
s1：将决策过程分为预见期末、预见期末所在月份的月末和年末三个时间节点，选取发电量最大作为目标函数，根据水库运行实际建立约束条件，其中约束条件包括水量平衡约束、水电站物理特性约束；
[0067]
s2：根据历史资料采用变动流量方式计算各离散点位年末余留效益，其中，变动流量过程计算方式为：以末水位等于初水位为前提，根据不同年份的历史月径流资料进行计算最后取平均，得到的结果作为每一年年末各点的余留效益；
[0068]
s3：根据预见期末所在月份的月末至年末的预报月径流资料，进行以月为时段的
优化调度计算，获得该月末各离散点位的余留效益；
[0069]
s4：再根据预见期末至该月末预报日流量，进行以日为时段的优化调度计算，得到预见期末各离散点位的余留效益；
[0070]
s5：以当前时刻水位为初始值，以预见期末各离散点位的余留效益为边界条件，在预见期内进行以小时为时段的水电站短期经济运行，得到最优决策过程以及相应的最优效益，基于所述最优决策过程对水电站完成一日的调度决策。
[0071]
在一种实施方式中，所述方法还包括：
[0072]
s6：完成一日的调度决策后，再进行下一日开始的调度，首先比较预见期末时间是否超出本月，若超出，则以下一个月作为预见期末执行步骤s3～s5，若未超出，则执行步骤s4～s5。
[0073]
具体地，本发明的技术方案主要包括四个部分：一(s1)、将决策过程分为预见期(当前日期开始的若干天，如3天)末、预见期末所在月份的月末和年末三个重要时间节点，选取发电量最大作为目标函数，根据水库运行实际建立约束条件，根据水位范围离散出符合实际需要的水位点；二(s2
‑
s4)、根据历史资料采用变动流量方式计算各离散点位年末平均余留效益，根据预见期末所在月份的月末至年末预测月径流资料计算月末各离散点位余留效益，再根据预见期末至该月末预报日流量计算预见期末各离散点位余留效益；三(s5)、从当前时刻到预见期末以短期预报入库流量为依据，进行以小时为时段的短期经济运行，得到水电站最优决策过程以及相应的最优效益；四(s6)、完成一日的调度决策后，再进行下一日开始的调度，首先比较预见期末时间是否超出该月，若超出，则计算下一月末各离散水位点所对应余留效益，再根据预见期末至该月末预报日流量计算预见期末各离散点位余留效益；若未超出，则直接根据本月余留效益进行下一次计算。
[0074]
本发明所提出的针对多年调节水库的智能优化调度方法，创造了一种结合历史径流资料与预报径流资料获得水电站最优发电量的调度方式，获得预见期内的最优决策以及最优效益，且该方法概念清晰，考虑全面。
[0075]
本发明中涉及的计算主要包括：
[0076]
1.水库蓄水状态离散
[0077]
在死水位与正常蓄水位(汛期为汛限水位)之间离散出不同的水位点，参见附图2和图3，分别表示特征水位范围以及离散点位，图2中，zs表示死水位，z
z
表示正常蓄水位。
[0078]
2.各重要时间节点余留效益计算
[0079]
在一年的总调度期内，假设运行一年后水位不变，即末水位等于初水位，采用不同年份的历史月径流资料进行多次计算最后取平均，得到的结果作为年末各点的余留效益；之后根据预见期末所在月的月末至年末的中长期月径流预报资料采取优化调度方式计算当月月末各水位点的余留效益；最后根据预见期末至月末预报日流量资料采取优化调度方式计算预见期末余留效益。
[0080]
3.预见期内短期优化调度
[0081]
采用dddp算法(离散微分动态规划)，根据已知的初水位、短期水文预报成果和预见期末每个末水位点的余留效益，以小时为时段进行优化计算，获得预见期内水电站的最优决策过程并得到最优效益。
[0082]
4.循环计算
[0083]
根据附图5过程，完成一日的调度决策后，再进行下一日开始的调度，首先比较预见期末时间是否超出该月，若超出，则计算下一月末各离散水位点所对应余留效益，接着重复上述步骤进行下一次计算，若未超出，则直接根据本月余留效益进行下一次计算。
[0084]
本发明在水库正常运行情况下的蓄水状态(水位或库容)进行离散时，考虑了蓄水状态下限和上限，并利用了多年历史资料，计算余留效益时考虑了来水的影响。采用的预报流量可以随着时间推移进行更新，操作具有实时性。在预见期末水位未知的情况下，最优决策可以提供一个合理的末水位用于调度管理者参考。并且可根据时间推移而逐日计算，极大减少了人工干预，体现了智能的要求。
[0085]
在一种实施方式中，步骤s1中的目标函数为：
[0086][0087]
式中：为计划期t内水电站总发电收入，元；p(t)为t时段水电站总出力，kw；c(t)为t时段电价，元/(kwh)，可反映峰谷电价差异，若c(t)＝1，目标函数就表示发电量最大；
[0088]
水量平衡约束包括水量平衡方程：
[0089]
v(t+1)＝v(t)+[q
rk
(t)
‑
q
fd
(t)
‑
q
qs
(t)]
×
δt(t)
[0090]
v(t+1)、v(t)代表t+1时刻与t时刻的库容，δt(t)代表时段长，q
rk
(t)、q
fd
(t)、q
qs
(t)分别代表水库t时刻的入库、发电、弃水流量；
[0091]
水电站物理特性约束包括库水位或库容约束、水位库容曲线约束、尾水位下泄曲线约束、闸门泄流曲线约束和水头出力限制曲线限制；
[0092]
(1)库水位或库容约束：
[0093][0094]
z
sy
(t+1)为t+1时刻的上游水位，v(t+1)为t+1时刻的库容，为t+1时刻的上游水位下限，通常为死水位，为t+1时刻的上游水位上限，通常在汛期为汛限水位，非汛期为正常蓄水位，v
min
(t+1)与v
max
(t+1)分别为上游水位下限和上游水位上限所对应的库容；
[0095]
(2)水位库容曲线约束：
[0096]
z
sy
(t+1)＝f
zv
[v(t+1)]
[0097]
z
sy
(t+1)为t+1时刻的上游水位，v(t+1)为t+1时刻的库容，f
zv
为水位库容曲线关系；
[0098]
(3)尾水位下泄曲线约束：
[0099]
z
xy
(t)＝f
zq
[q
ck
(t)]
[0100]
z
xy
(t)为t时刻的下游水位，q
ck
(t)为t时刻的出库流量，f
zq
为尾水位下泄流量曲线关系；
[0101]
(4)闸门泄流曲线约束：
[0102]
q
qs
(t)≤f
qz
[z
xy
(t)]
[0103]
q
qs
(t)为t时刻的弃水流量，z
xy
(t)为t时刻的下游水位，f
qz
为水位与闸门泄流能力曲线关系；
[0104]
(5)水头出力限制曲线限制；
[0105]
p(t)≤f
ph
[h0(t)]
[0106]
p(t)为t时刻的出力，h0(t)为t时刻的水头，f
ph
为水头与出力曲线关系。
[0107]
在约束条件限制的基础上进行优化调度计算，可以获得余留效益和阶段效益。
[0108]
在一种实施方式中，步骤s2包括：
[0109]
获取不同年份的历史月径流资料，具体包括n年历史月径流数据；
[0110]
根据每年的历史月径流资料进行一年内的优化调度计算，则年末余留效益计算公式为：
[0111][0112]
其中，表示第n年所得的优化调度最优发电量，为计算期初末水位均为xm时所对应的余留效益，表示采用第t年径流资料进行优化调度计算所得到的最优发电量。
[0113]
具体实施过程中，在一年的总调度期内，假设一年运行后水位不变，即末水位等于初水位，计算过程可参见图1，z
tn
表示第t个时刻离散出的第n个水位点，图中指年末。采用不同年份的历史月径流资料进行多次计算最后取平均，得到的结果作为年末各点的余留效益，例如当初水位为100m时对应余留效益为例根据上述公式可以计算期初末水位均为100m时所对应的余留效益。
[0114]
在一种实施方式中，步骤s3包括：
[0115]
在预见期末到年末这一时段内以月为单位采用离散微分动态规划算法，计算到年末离散出的每一末水位点所对应的发电量，得到的视作阶段效益其中n指对应的第n个水位点，则该初水位到年末各个水位离散点总效益为：
[0116][0117]
其中代表初水位为xm时到年末第n个离散水位的总效益，表示年末第n个点对应的余留效益，代表初水位为xm时到年末第n个离散水位的阶段效益；
[0118]
从初水位为xm时到年末各个离散点对应的总效益中，筛选总效益的最大值作为初水位为xm时对应的余留效益：
[0119][0120]
代表水位为xm时的月末余留效益，代表初水位为xm时到年末各个离散水位点的总效益。
[0121]
具体来说，当初水位为100m时，代表初水位为100m
时到年末各个离散水位点(共n个)的总效益，然后从这n个总效益中筛选出最大值作为初水位为100m时对应的余留效益。
[0122]
然后执行步骤s4,根据预见期末至该月末预报日流量，进行以日为时段的优化调度计算，得到预见期末各离散点位的余留效益，这里计算方式与第三步骤s3相同，唯一的区别在于优化调度的计算期长度以及计算期尺度(由月尺度变为日尺度)，这里不予赘述。
[0123]
其中，图4示出了本发明实施例中智能优化调度模型三段计算方式。
[0124]
步骤s5中以当前时刻水位为初始值，根据预见期内的径流预报，以预见期末各离散点位的余留效益为边界条件，在预见期内进行以小时为时段的水电站短期经济运行，得到最优决策过程，即在预见期内以小时为单位的水位过程，以及相应的最优效益(发电量)，计算公式参见目标函数，该决策过程可以作为水电站的最优决策过程，选取第一日的决策过程进行第一日的调度决策。
[0125]
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。