技术特征：
1.一种梯级水库汛期运行水位协同浮动的运用方法，其特征在于，包括如下步骤：步骤一，剖析梯级水库防洪库容的补偿机理，厘清梯级水库汛期运行水位协同浮动关系；步骤二，对多个梯级水库，考虑梯级水库间水力联系，依从上游至下游顺序采用预报预泄法推求各水库汛期运行水位上浮的浮动量；步骤三，建立以有效预泄时间内防洪风险最小化和发电量最大化为优化目标函数、水库出库流量为决策变量的协同优化模块，并根据该协同优化模块优化防洪库容等效置换系数和梯级水库汛期运行水位浮动量。2.根据权利要求1所述的梯级水库汛期运行水位协同浮动的运用方法，其特征在于，步骤一中采用库容置换方法剖析梯级水库防洪库容补偿机理，其中，库容置换过程如下：若梯级水库按上游至下游依次编号为1，2，...，n，梯级水库相应的防洪库容分别为随机生成一个防洪库容置换系数矩阵随机生成一个防洪库容置换系数矩阵则梯级水库的汛期运行水位浮动为则梯级水库的汛期运行水位浮动为其中，f1(
·
)，f2(
·
)，...，f
n
(
·
)为第1，2，...，n个水库的库容转水位的函数；梯级水库的防洪库容置换需满足如下限定性约束条件：式中：为第i个水库的防洪库容置换系数；当时，代表水库的防洪库容减少且相应的汛期运行水位向上浮动；当时，代表水库的防洪库容不变且相应的汛期运行水位维持分期汛限水位；当时，代表水库的防洪库容增加且相应的汛期运行水位向下浮动。3.根据权利要求1所述的梯级水库汛期运行水位协同浮动的运用方法，其特征在于，所述步骤二中梯级间水库间水力联系模拟如下：先利用水文学法和水动力学法模拟第2，3，...，n个水库的入库流量，得到水库入库流量样本，然后利用这些样本训练机器学习模型，得到梯级间水力联系的代理模型，相邻水库水力联系计算式如下：式中，为第i个水库t-1时刻的出库流量；为第i+1个水库t-1时刻的出库流量；为第i+1个水库t时刻的入库流量；为第i+1个水库t时刻的区间入流；g(
·
)为机器学习模型的映射函数。4.根据权利要求3所述的梯级水库汛期运行水位协同浮动的运用方法，其特征在于，所述步骤二中基于预报预泄法推求水库汛期运行水位上浮的浮动量如下：综合考虑降雨预报及洪水预报预见期长度或精度、预泄能力约束和下游安全泄量，水库在有效预泄时间内按预泄能力对应的提前泄流水量将水库汛期运行水位向上浮动，浮动
量δz计算如下：量δz计算如下：量δz计算如下：量δz计算如下：量δz计算如下：量δz计算如下：式中：δz
i
为第i个水库汛期运行水位上浮的浮动量；t
i
为第i个水库有效预泄时间；为有效预泄时间内第i个水库的平均预报入库流量，为有效预泄时间内第i个水库的平均出库流量；为第i个水库t时刻的入库流量，为第i个水库t时刻的出库流量；ε
i
(t)为第i个水库t时刻入库流量的预报误差；为第i个水库下游防护目标的堤防过流能力；为第i个水库的水文预报预见期，含降雨预报和洪水预报预见期；为第i个水库作业预报时间、决策时间和闸门操作时间之和；δt为计算时段；w
i
为有效预泄时间内第i个水库的预泄水量。5.根据权利要求1所述的梯级水库汛期运行水位协同浮动的运用方法，其特征在于，步骤三中，以有效预泄时间内防洪风险最小化和发电量最大化为优化目标函数分别为：最大化为优化目标函数分别为：式中：fr为防洪风险最小化目标函数，hg为发电量最大化目标函数；v
i
(t)为第i个水库t时刻的库容；为第i个水库分期汛限水位对应的库容；m为年数或场次洪水数；p
i
(t)为第i个水库t时刻的出力；min{
·
}为最小化函数，max{
·
}为最大化函数；在建立上述的优化目标函数后，以水库出库流量为决策变量，综合考虑防洪库容置换、预报预泄过程、防洪调度规则、水力联系和水库水位、流量、出力约束，采用元启发式优化算法协同优化梯级水库调度过程，推求防洪库容等效置换系数和梯级水库汛期运行水位协同浮动量。6.根据权利要求5所述的梯级水库汛期运行水位协同浮动的运用方法，其特征在于，步骤三中的优化目标函数及各水库/电站需满足以下约束条件：a.水量平衡约束：
式中，v
i
(t)和v
i
(t+1)分别为第i个水库t时刻和t+1时刻的库容；为第i个水库t时刻的入库流量，为第i个水库t时刻的出库流量，分别为第i个水库t时刻的损失流量；b.相邻水库水力联系：式中，为第i个水库t-1时刻的出库流量，分别为第i个水库t时刻的出库流量；为第i+1个水库t-1时刻的出库流量，分别为第i+1个水库t时刻的入库流量；为第i+1个水库t时刻的区间入流；g(
·
)为机器学习模型的映射函数；c.水库水位约束：式中，z(t)为第i个水库t时刻的库水位；为第i个水库运行水位下限，汛期取分期汛限水位；为第i个水库运行水位上限，汛期取步骤二中推求的汛期运行水位：d.水库出库流量约束：式中，为第i个水库出库流量下限，综合考虑灌溉、航运和生态环境需求来定；为第i个水库出库流量上限，由下游堤防过流能力、水电站过流能力和水库泄流能力决定；e.水电站出力约束：式中，为第i个水库水电站的出力下限，分别为第i个水库水电站的出力上限。7.根据权利要求6所述的梯级水库汛期运行水位协同浮动的运用方法，其特征在于，步骤三中以元启发式优化算法根据协同优化模块求解防洪库容等效置换系数和梯级水库汛期运行水位浮动量的过程如下：
①
初始化蝗虫算法参数和编码决策变量；设置蝗虫种群规模n
pop
、最大迭代次数i
max
、衰减系数的最小值c
min
和最大值c
max
；采用实数编码方式编码并随机生成一定规模的决策变量，即水库出库流量
②
评价适应度；对种群进行快速非支配排序，将蝗虫种群划分等级；计算蝗虫种群的拥挤距离；执行轮盘操作选择适应度最大个体，以生成规模仍为n
pop
的子代蝗虫种群；
③
更新衰减系数和蝗虫个体在集群中位置；模拟梯级水库防洪库容置换、预报预泄过程、防洪调度规则和水力联系，且需满足水库水位/流量/出力约束，依据蝗虫集群行为，将蝗虫个体活动区域划分为吸引区(0，a1]、舒适区(a1，a2]、排斥区(a2，a3]，其中0＜a1＜a2＜
a3；对第k次迭代，衰减系数c
k
更新计算式和蝗虫个体在集群中位置的计算式分别如下：c
k
＝c
max-k(c
max-c
min
)/i
max
；；式中，c
k
为第k次迭代的衰减系数，用于调整吸引区、舒适区和排斥区范围；为第j个蝗虫个体第d个决策变量的位置向量，其中，d＝1，2，...，d，d为决策变量数目；为第d个决策变量的最佳位置向量；u
d
为第d个决策变量约束的上限，l
d
为第d个决策变量约束的下限；第n个蝗虫个体的第d个决策变量的位置，为第j个蝗虫个体的第d个决策变量的位置；为从第j个蝗虫个体位置到第n个蝗虫个体位置的单位向量；s(x)为自变量x的集群行为函数，其中参数f和l分别为吸引区的吸引强度和吸引距离尺度；
④
判断计算终止条件；若当前迭代次数小于最大迭代次数i
max
，则重复步骤
②
至
④
；否则终止计算并输出优化结果。
技术总结
本发明提供一种梯级水库汛期运行水位协同浮动的运用方法，包括：剖析梯级水库防洪库容的补偿机理，厘清梯级水库汛期运行水位协同浮动关系；对多个梯级水库，考虑梯级水库间水力联系，依从上游至下游顺序采用预报预泄法推求各水库汛期运行水位上浮的浮动量；建立以有效预泄时间内防洪风险最小化和发电量最大化为优化目标函数、水库出库流量为决策变量的协同优化模块，并根据该协同优化模块优化防洪库容等效置换系数和梯级水库汛期运行水位浮动量。本发明可大幅提升仿生进化算法的寻优效率，对梯级水库汛期运行水位协同浮动运用的预报调度一体化程度更高且可实现梯级水库防洪库容的等效置换。库容的等效置换。库容的等效置换。
技术研发人员：周研来 郭生练 许崇育 王俊 尹家波
受保护的技术使用者：中国长江三峡集团有限公司
技术研发日：2021.11.17
技术公布日：2022/2/28