并联水库防洪优化调度方法、系统、介质及设备与流程

1.本发明涉及水利工程技术领域，具体为并联水库防洪优化调度方法、系统、介质及设备。


背景技术：

2.本部分的陈述仅仅是提供了与本发明相关的背景技术信息，不必然构成在先技术。
3.并联水库相互之间影响较小，能够为下游干流河道防洪提供了更广泛的调度空间。同时，由于并联水库之间水力联系较小，对两水库之间下泄流量过程的优化往往面临巨大的不确定性。并联水库一方面需要考虑上游来水和区间来水情况，需要根据水库水位库容关系及防汛水位的限制要求，保障自身大坝安全；另一方面需要根据下游河道或者防洪对象的安全等级，控制并联水库各自的下泄流量过程，实现错峰下泄，减小下游控制断面洪峰的要求。
4.根据不同水库以及下游控制断面等防洪对象的重要程度，实现协调调度，最大限度的流域内防洪对象的整体安全性，目前尚无成熟的方法。


技术实现要素：

5.为了解决上述背景技术中存在的技术问题，本发明提供并联水库防洪优化调度方法、系统、介质及设备，基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
6.为了实现上述目的，本发明采用如下技术方案：
7.本发明的第一个方面提供并联水库防洪优化调度方法，包括以下步骤：
8.基于并联的水库在t时刻的水位和下游控制断面在t时刻的流量，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，构建最优值函数；
9.基于最优值函数依次经水位离散、末阶段寻优、中间阶段寻优和时段递推后，得到并联水库中其中一个水库的最优泄流策略，作为单库防洪优化策略；
10.将得到的单库防洪优化策略作为初始解，基于最优值函数求取并联水库中下一水库的最优泄流策略，重复执行该步骤，得到并联水库中各水库的下泄流量和对应的水位。
11.最优值函数为：
[0012][0013]
式中：表示第i个水库在t时刻且状态为k时，对应策略的最优值；hi′
(t)为第i个水库在t时刻时的水位经归一化处理后的值；q
′
(t)为下游控制断面在t时刻的流量经归一化处理后的值；αi为第i个水库的重要性权重，β为下游控制断面的重要性权重，αi+β＝1。
[0014]
得到并联水库中其中一个水库的最优泄流策略，包括：
[0015]
水位离散；具体为：第i个水库设定调度周期内的初始水位和末期水位分别为z
i,1
和z
i,t+1
，在水位允许范围内离散中间时段t的初水位分别为其中，t＝2～t；
[0016]
末阶段寻优；具体为：令时段t初水位若由时段平均来水、初末水位依据水量平衡方程确定的时段平均泄流量不满足对应的泄流约束，则在调度期末水位的限制下，该初水位不可行；当初水位可行时，将该初水位对应的时段最优泄流策略带入递推方程得到该初水位对应的最优余留期效益令k＝k+1，重复上述末阶段寻优的过程直至k＝k，得到t时段各初水位对应的时段最优泄流策略和余留期效益
[0017]
中间阶段寻优；具体为：令t＝t-1，时段t初水位末水位若由该初末水位、时段平均来水依据水量平衡方程确定的时段平均泄流量不满足对应的泄流约束或者经末时段寻优已判定该末水位不可行，则在初水位时，该末水位不可行，此时令k2＝k2+1继续判断，若直至k2＝k时，所有末水位都不可行，则该初水位不可行；若在初水位时，末水位可行，则计算此时时段平均泄流结合末水位的最优余留期效益推算当前初水位的余留期效益；令k2＝k2+1，重复上述中间阶段寻优的过程，直至k2＝k；按照递推方程优选出初水位时，时段最优泄流决策和最优余留效益令k1＝k1+1，重复上述中间阶段寻优的过程，直至k1＝k。
[0018]
时段递推；具体为：令t＝t-1，重复中间阶段寻优的过程，直至t＝1；
[0019]
推求最优决策；具体为：令t＝1，时段初水位为调度期初水位z
i,1
，由z
i,1
对应的最优泄流决策依据水量平衡确定本末水位z
i,t+1
；令t＝2，时段初水位z
i,1
＝z
i,t-1
，由z
i,1
对应的最优泄流决策依据水量平衡确定本末水位z
i,t+1
；令t＝t+1，重复推求最优决策的过程，直至t＝t，则该调度周期的最优泄流决策
[0020]
本发明的第二个方面提供实现上述方法所需的系统，包括：
[0021]
调度目标模块，被配置为：基于并联的水库在t时刻的水位和下游控制断面在t时刻的流量，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，构建最优值函数；
[0022]
单库求解模块，被配置为：基于最优值函数依次经水位离散、末阶段寻优、中间阶段寻优和时段递推后，得到并联水库中其中一个水库的最优泄流策略，作为单库防洪优化策略；
[0023]
多库求解模块，被配置为：将得到的单库防洪优化策略作为初始解，基于最优值函数求取并联水库中下一水库的最优泄流策略，重复执行该步骤，得到并联水库中各水库的下泄流量和对应的水位。
[0024]
本发明的第三个方面提供一种计算机可读存储介质。
[0025]
一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现
如上述所述的并联水库防洪优化调度方法中的步骤。
[0026]
本发明的第四个方面提供一种计算机设备。
[0027]
一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述所述的并联水库防洪优化调度方法中的步骤。
[0028]
与现有技术相比，以上一个或多个技术方案存在以下有益效果：
[0029]
基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
附图说明
[0030]
构成本发明的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。
[0031]
图1是本发明一个或多个实施例提供的并联水库示意图；
[0032]
图2是本发明一个或多个实施例提供的水库a下泄流量随时间变化的曲线示意图；
[0033]
图3是本发明一个或多个实施例提供的水库b下泄流量随时间变化的曲线示意图；
[0034]
图4是本发明一个或多个实施例提供的控制断面总流量随时间变化的曲线示意图。
具体实施方式
[0035]
下面结合附图与实施例对本发明作进一步说明。
[0036]
应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
[0037]
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0038]
正如背景技术中所描述的，根据不同水库以及下游控制断面等防洪对象的重要程度，实现协调调度，最大限度的流域内防洪对象的整体安全性，目前尚无成熟的方法。
[0039]
因此以下实施例给出并联水库防洪优化调度方法、系统、介质及设备，基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
[0040]
实施例一：
[0041]
并联水库防洪优化调度方法，包括以下步骤：
[0042]
1、假定各库的初始策略和状态序列。
[0043]
确定并联水库及下游控制断面防洪优化调度的目标函数、分配各自权重，分配权重的目的在于各水库和防洪控制断面执行设定的防汛水位(流量)的程度。
[0044]
确定各水库和控制断面的约束条件和决策变量，设定计算参数，包括最大迭代次数m、各水库初始离散步长、初始离散数目k、收敛精度ε、水库数量n和时段总数t；
[0045]
采用人工经验决策或常规方法生成满足各约束条件的各个水库的初始调度过程。
[0046]
2、单库求解策略。
[0047]
先对第一个水库进行优化，其余n-1个水库的运行策略和状态变量序列暂时保持不变，相当于对n维决策向量加n-1个约束条件。此时，该子问题只含有一个状态变量和一个决策变量，就可以用一维的优化算法(如poa)求解，从而可得第一个水库的运行策略。
[0048]
3、多库求解策略。
[0049]
dpsa(dynamic programming successive approximation，动态规划逐次逼近)的优点是将n维动态规划问题转化成n个一维子问题，从而使得计算工作量只是随维数n成线性增长而不是成指数增长，间接的降低计算维数。因此，当面临待解决问题的维数越高，该方法所节省的计算工作量越显著。在解决工程实际问题时常将该方法与其他单库优化算法相结合，首先用dpsa将n维问题分解成只有一维的n个子问题，然后用优化算法实现每个子问题的最优化，降低并联水库调度寻优的维数。
[0050]
获取第一个水库优化调度过程的基础上，再对第二个水库进行优化，除第一个水库保持新的运行策略外，其余水库仍保持初始策略，用优化算法求解第二个水库的最优运行策略。
[0051]
用同样的方法对剩下的水库分别进行优化，得到各个水库的新的运行策略。
[0052]
再从第一个水库开始，重复上述步骤，进行第2轮、第3轮
……
迭代计算，直到目标函数不再有改善或前后两轮迭代目标函数值满足收敛要求为止。
[0053]
关于dpsa-poa算法
[0054]
1、dp-poa求解单库优化
[0055]
本实施例采取动态规划算法(dynamic programming，dp)嵌套逐步优化算法(progressive optimality algorithm，poa)的方法推求水库的下泄流量的最优策略，其中，dp采用逆序递推。
[0056]
对于单库防洪调度问题，以调度时段为阶段，时段号t代表阶段变量(t＝1～t)，以i(i＝1，2)水库各时段初水位z
ik
(t)为i水库第t阶段的状态变量，时段平均泄流量qi(t)为决策变量。
[0057]
目标函数同时考虑水库水位和下游控制断面的流量，如式(1)所示：
[0058][0059]
式中：表示i水库t阶段状态为z
ik
(t)(k＝1～k)时，对应最优子策略u
i,t,t
＝{qi(t),qi(t-1),
…
,qi(t)}的最优值；αi为第i个水库的重要性权重，β为下游控制断面的重要性权重，αi+β＝1。
[0060]
状态转移方程为水量平衡方程，如式(2)所示，
[0061][0062]
式中：表示i水库t阶段状态为z
ik
(t)(k＝1～k)时，对应最优子策略u
i,t,t
＝
{qi(t),qi(t-1),
…
,qi(t)}为第t～t阶段目标函数的最优值。
[0063]
运用dp求解水库防洪优化调度问题步骤如下。
[0064]
(1)水位离散。对于水库i(i＝1～m)，给定调度周期初末水位分别为z
i,1
和z
i,t+1
，在水位允许范围内离散中间时段t(t＝2～t)的初水位为
[0065]
(2)末阶段寻优。令时段t初水位若由时段平均来水、初末水位依据水量平衡方程确定的时段平均泄流量不满足泄流相应约束，则说明在调度期末水位的限制下该初水位不可行。当初水位可行时，即为该初水位对应的时段最优泄流策略，并将带入递推方程求出该初水位对应的最优余留期效益令k＝k+1，重复以上过程直至k＝k，得到t时段各初水位对应的时段最优泄流策略和余留期效益
[0066]
(3)中间阶段寻优。令t＝t-1，时段t初水位末水位若由该初末水位、时段平均来水依据水量平衡方程确定的时段平均泄流量不满足泄流相应约束或者经末时段寻优已判定该末水位不可行，则说明在初水位时，该末水位不可行，此时令k2＝k2+1，继续以上判断，若直至k2＝k，所有末水位都不可行，说明该初水位不可行。若在初水位时，末水位可行，则计算此时时段平均泄流结合末水位的最优余留期效益推算当前初水位的余留期效益。令k2＝k2+1，重复以上过程，直至k2＝k。按照递推方程优选出初水位时，时段最优泄流决策和最优余留效益令k1＝k1+1，重复以上过程，直至k1＝k。
[0067]
(4)时段递推。令t＝t-1，重复步骤2，直至t＝1。
[0068]
(5)推求最优决策。令t＝1，时段初水位为调度期初水位z
i,1
，由z
i,1
对应的最优泄流决策依据水量平衡确定本末水位z
i,t+1
。令t＝2，时段初水位z
i,1
＝z
i,t-1
，由z
i,1
对应的最优泄流决策依据水量平衡确定本末水位z
i,t+1
。令t＝t+1，重复以上过程，直至t＝t，则该调度周期的最优泄流决策
[0069]
2、流域内多个水库优化求解
[0070]
dpsa算法求解流域内并联水库防洪优化调度问题计算步骤如下：
[0071]
(1)依据dp求解得到的单库防洪优化策略，将其作为初始值求解各库防洪优化调度问题。
[0072]
各单库防洪调度最优策略对应的水位过程组成流域水库群防洪调度初始策略序列
[0073]
(2)以各水库的各时段初水位z
i,t
作为决策变量。在初始策略基础上，调整z
1,t
(t＝2～t)依次为第1个水库水位离散值，固定其他变量，运用流域水库群联合防洪优化调度目标函数对z
1,t
寻优，舍弃不满足约束的解，求得令调整z
2,t
(t＝2～t)依次为第2个水库水位离散值，固定其他变量，运用流域水库群联合防洪优化调度目标函数对z2,t
寻优，舍弃不满足约束的解，求得令继续调整z
3,t
,z
4,t
,
…zm,t
的值；以此类推直至t＝t，得到新的决策序列
[0074]
(3)以新的决策序列ui′
,1,t
作为初始轨迹，重复步骤(2)的过程，得到决策序列ui″
,
″
1,t
，比较ui″
,1,t
和ui″
,
″
1,t
，若差值小于给定的精度要求值，则输出ui″
,
″
1,t
，作为最优决策序列。
[0075]
关于河道洪水演进计算
[0076]
采用分段马斯京根法对河流洪水过程进行分河段连续演算。假定n为上下游两断面间的河段数，k为两个断面间河段稳定流的传播时间(h)，l为河段总长度(km)，l1为单元河段长(km)，x
l
为单元河段流量比重系数，l为特征河长(km)。在单元河段稳定流传播时间k
l
＝δt的条件下，马斯京根法汇流参数可用下式进行计算：
[0077][0078]
计算出演算河段的单元河段数n及单元河段流量比重系数x
l
；根据单元河段流量比重系数x
l
，按下面公式计算出单元河段演算系数c0,c1和c2：
[0079][0080]
根据单元河段数及c0,c1和c2值，按下面公式可以计算出河段汇流系数p
mn
：
[0081][0082]
式中，m为汇流系数所对应的时段数；p
0n
单元河段数为n的零时段汇流系数；p
mn
单元河段数为n的第m时段汇流系数。
[0083]
在得到p
mn
的基础上，按照倍比、叠加原理使用卷积公式将上游大型水库的泄洪流量过程qi汇流演算到下游控制断面，按照下式可以计算出下游断面时段末的出流量：
[0084][0085]
上述方法基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
[0086]
实例演示
[0087]
如图1-4所示，有水库a和水库b两座并联水库，假设已知水库a的入流过程为q
r,a
(t)，水库b的入流过程为q
r,b
(t)，以及水库a和水库b坝址至下游控制断面间的区间汇流过程qq(t)。水库a和水库b的汛限水位、起调水位、各级控泄水位及对应最大泄流量、下游控制断面的最大流量等均作为已知条件。
[0088]
定义水库a的下泄过程为q1(t)，水库b的下泄过程为q2(t)，利用马斯京根法演进至流域下游控制断面处的流量过程分别为q
’1(t)和q
’2(t)；区间汇流过程q3(t)；则下游控制断面处的洪水过程可表示为q(t)＝q
’1(t)+q
’2(t)+q3(t)。
[0089]
考虑目标函数：
[0090]
minmaxq(t)＝min{maxq’(t)+maxh
’1(t)+maxh
’2(t)}(7)
[0091]
式中，h1(t)为水库a在t时刻的水位，h2(t)为水库b在t时刻的水位，q(t)为下游控制断面在t时刻的流量。对应的h
’1(t)为经归一化后，水库a在t时刻的水位；h
’2(t)为经归一化后，水库b在t时刻的水位；q’(t)为经归一化后，下游控制断面在t时刻的流量。
[0092]
约束条件如下：
[0093]
水库水量平衡约束：
[0094][0095]
式中：ii(t)、ii(t+1)分别为i水库t和t+1时刻的入库流量；qi(t)、qi(t+1)分别为i水库t时刻、t+1时刻出库流量；vi(t)和vi(t+1)分别为i水库t时刻和t+1时刻的库容值。
[0096]
泄流约束：
[0097]
水库a：
[0098]
水库b：
[0099]
末水位约束：
[0100]
zi(t+1)＜z
i,js
[0101]
式中：zi(t+1)为i水库调度期计算末水位；z
i,js
为i水库调度期末的预期水位。水库a的预期末水位为179.95m，水库b的预期末水位为177.51m。
[0102]
流域下游控制断面水位约束：
[0103]hk
(t)≤h
k,js
[0104]
式中，hk(t)为第k个控制断面t时刻的水位，单位为m；h
k,js
为第k个控制断面的最高允许水位，单位为m。本例中h
k,js
＝91.71m。
[0105]
所有变量均为非负数值。
[0106]
实施例二：
[0107]
实现上述方法的系统，包括：
[0108]
调度目标模块，被配置为：基于并联的水库在t时刻的水位和下游控制断面在t时刻的流量，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，构建最优值函数；
[0109]
单库求解模块，被配置为：基于最优值函数依次经水位离散、末阶段寻优、中间阶段寻优和时段递推后，得到并联水库中其中一个水库的最优泄流策略，作为单库防洪优化策略；
[0110]
多库求解模块，被配置为：将得到的单库防洪优化策略作为初始解，基于最优值函数求取并联水库中下一水库的最优泄流策略，重复执行该步骤，得到并联水库中各水库的下泄流量和对应的水位。
[0111]
上述系统基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
[0112]
实施例三：
[0113]
本实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如上述实施例一所述的并联水库防洪优化调度方法中的步骤。
[0114]
上述方法基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
[0115]
实施例四：
[0116]
本实施例提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述实施例一所述的并联水库防洪优化调度方法中的步骤。
[0117]
基于单库动态规划寻优和多库逐步优化结合的方法，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，对并联水库各自的下泄流量进行优化调度，实现水库和下游控制断面的防洪安全。
[0118]
以上实施例二至四中涉及的各步骤或模块与实施例一相对应，具体实施方式可参见实施例一的相关说明部分。术语“计算机可读存储介质”应该理解为包括一个或多个指令集的单个介质或多个介质；还应当被理解为包括任何介质，所述任何介质能够存储、编码或承载用于由处理器执行的指令集并使处理器执行本发明中的任一方法。
[0119]
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技
术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

技术特征：
1.并联水库防洪优化调度方法，其特征在于，包括以下步骤：基于并联的水库在t时刻的水位和下游控制断面在t时刻的流量，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，构建最优值函数；基于最优值函数依次经水位离散、末阶段寻优、中间阶段寻优和时段递推后，得到并联水库中其中一个水库的最优泄流策略，作为单库防洪优化策略；将得到的单库防洪优化策略作为初始解，基于最优值函数求取并联水库中下一水库的最优泄流策略，重复执行该步骤，得到并联水库中各水库的下泄流量和对应的水位。2.如权利要求1所述的并联水库防洪优化调度方法，其特征在于，所述最优值函数为：式中：表示第i个水库在t时刻且状态为k时，对应策略的最优值；h
i
′
(t)为第i个水库在t时刻时的水位经归一化处理后的值；q
′
(t)为下游控制断面在t时刻的流量经归一化处理后的值；α
i
为第i个水库的重要性权重，β为下游控制断面的重要性权重，α
i
+β＝1。3.如权利要求1所述的并联水库防洪优化调度方法，其特征在于，所述水位离散，具体为：第i个水库设定调度周期内的初始水位和末期水位分别为z
i,1
和z
i,t+1
，在水位允许范围内离散中间时段t的初水位分别为分别为其中，t＝2～t。4.如权利要求1所述的并联水库防洪优化调度方法，其特征在于，所述末阶段寻优，具体为：令时段t初水位若由时段平均来水、初末水位依据水量平衡方程确定的时段平均泄流量不满足对应的泄流约束，则在调度期末水位的限制下，该初水位不可行；当初水位可行时，将该初水位对应的时段最优泄流策略带入递推方程得到该初水位对应的最优余留期效益令k＝k+1，重复上述末阶段寻优的过程直至k＝k，得到t时段各初水位对应的时段最优泄流策略和余留期效益5.如权利要求1所述的并联水库防洪优化调度方法，其特征在于，所述中间阶段寻优，具体为：令t＝t-1，时段t初水位末水位若由该初末水位、时段平均来水依据水量平衡方程确定的时段平均泄流量不满足对应的泄流约束或者经末时段寻优已判定该末水位不可行，则在初水位时，该末水位不可行，此时令k2＝k2+1继续判断，若直至k2＝k时，所有末水位都不可行，则该初水位不可行；若在初水位时，末水位可行，则计算此时时段平均泄流结合末水位的最优余留期效益推算当前初水位的余留期效益。
6.如权利要求1所述的并联水库防洪优化调度方法，其特征在于，所述时段递推，具体为：令t＝t-1，重复中间阶段寻优的过程，直至t＝1。7.如权利要求1所述的并联水库防洪优化调度方法，其特征在于，得到并联水库中其中一个水库的最优泄流策略，还包括：推求最优决策，具体为：令t＝1，时段初水位为调度期初水位z
i,1
，由z
i,1
对应的最优泄流决策依据水量平衡确定本末水位z
i,t+1
；令t＝2，时段初水位z
i,1
＝z
i,t-1
，由z
i,1
对应的最优泄流决策依据水量平衡确定本末水位z
i,t+1
；令t＝t+1，重复推求最优决策的过程，直至t＝t，则该调度周期的最优泄流决策8.并联水库防洪优化调度系统，其特征在于，包括：调度目标模块，被配置为：基于并联的水库在t时刻的水位和下游控制断面在t时刻的流量，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，构建最优值函数；单库求解模块，被配置为：基于最优值函数依次经水位离散、末阶段寻优、中间阶段寻优和时段递推后，得到并联水库中其中一个水库的最优泄流策略，作为单库防洪优化策略；多库求解模块，被配置为：将得到的单库防洪优化策略作为初始解，基于最优值函数求取并联水库中下一水库的最优泄流策略，重复执行该步骤，得到并联水库中各水库的下泄流量和对应的水位。9.一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如上述权利要求1-7任一项所述的并联水库防洪优化调度方法中的步骤。10.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如权利要求1-7任一项所述的并联水库防洪优化调度方法中的步骤。

技术总结
本发明涉及并联水库防洪优化调度方法、系统、介质及设备，其中的并联水库防洪优化调度方法包括以下步骤：基于并联的水库在t时刻的水位和下游控制断面在t时刻的流量，以下游控制断面的洪峰流量最小和水库蓄洪达到的最高水位最低为调度目标，构建最优值函数；基于最优值函数依次经水位离散、末阶段寻优、中间阶段寻优和时段递推后，得到并联水库中其中一个水库的最优泄流策略，作为单库防洪优化策略；将得到的单库防洪优化策略作为初始解，基于最优值函数求取并联水库中下一水库的最优泄流策略，重复执行该步骤，得到并联水库中各水库的下泄流量和对应的水位。的下泄流量和对应的水位。的下泄流量和对应的水位。


技术研发人员：王俊 刘阳 贾菁菁 王海军 刘薇 陈翠英 曹升乐 王家昌 郭成
受保护的技术使用者：山东省水文中心
技术研发日：2023.04.14
技术公布日：2023/7/13