考虑来流频率差异的多年调节水库年末消落水位确定方法技术

本发明专利技术公开了一种考虑来流频率差异的多年调节水库年末消落水位确定方法，属于水资源管理领域。包括：将梯级系统中多年调节水库的死水位至正常蓄水位范围进行离散；以每个离散水位值作为多年调节水库的年末消落水位约束，以总发电量最大为目标建立梯级系统联合优化调度模型；以多年来流频率作为模型输入，以多维动态规划算法对模型进行求解；将发电量最大对应的离散水位作为每个来流频率的最优消落水位；基于最小二乘原理拟合来流频率与最优消落水位散点图，得到不同来流频率下的最优消落水位。本发明专利技术能够充分发挥多年调节水库的调节性能，提高梯级系统水能源转换效率和梯级总发电量，对指导含多年调节水库的梯级水库群实际调度运行具有重要意义。

【技术实现步骤摘要】
考虑来流频率差异的多年调节水库年末消落水位确定方法
本专利技术属于水资源管理领域，更具体地，涉及一种考虑来流频率差异的多年调节水库年末消落水位确定方法。
技术介绍
水电能源是目前全世界正在大力发展的一种优质高效能源，其具有清洁无污染、可再生性、运行成本低、出力响应迅速等特点。尤其是其对电力负荷的快速响应特性，使其在电力系统的安全稳定运行中扮演了极其重要的角色，发挥了极其重要的作用。目前水电能源的开发利用方式主要是通过在河道中筑坝建库，将分散的河道水能资源集中并抬高水头，再通过水电站-水轮发电机组使得水能高效地转化为电能。对于水库的类型，按照其调节性能的强弱，可以分为日调节，季调节、年调节及多年调节类型。一般而言，河流来水一般具有年内丰枯变化规律、年际间水量不同的特征。考虑到水文年的周期性，实际调度一般以一个水文年作为一个调度周期，此时年调节水库就能很好的对年内来水不均性进行调节。但是，实际水文过程中，年际间水量差异也存在，因此需要在年际间进行水量分配与调节，而多年调节水库就可以做到年际间的水量分配。因此对于多年调节水库，如何确定其在一个调度周期内的运行方式，对多年调节水库最终发电效益的影响较大。许多学者也在这方面做过相关的研究，取得了一些成果，但是已有研究大多集中在多年调节水库单独运行时的情形。随着众多水库的建设和投运，逐渐形成了具有上下游水量和水头联系的梯级水库群。梯级水库形成后，原有的单库调度方法已经不能适应库群联合调度的需求，需要针对性的对梯级水库系统中的多年调节水库调度方式进行研究,即确定梯级系统中多年调节水库的最优年末消落水位。此外，已有针对多年调节水库运行方式的研究仅从多年平均的角度考虑其年末消落水位的大小，没有考虑年际间来水差异，即忽略了不同年份来水频率差异的影响，存在一定的不足和尚待改进的地方。
技术实现思路
针对现有技术的以上缺陷或改进需求，本专利技术提供了一种考虑来流频率差异的多年调节水库年末消落水位确定方法，其目的在于在考虑不同年份来流频率差异的情况下，确定梯级水库系统中多年调节水库的最优年末最佳消落水位。为实现上述目的，本专利技术的提供了一种一种考虑来流频率差异的多年调节水库年末消落水位确定方法，包括：S1.根据设定的水位离散精度，将梯级系统中多年调节水库的死水位至正常蓄水位范围离散为一系列离散水位值；S2.以每一个离散水位值作为多年调节水库的年末消落水位约束，以总发电量最大为目标建立梯级系统联合优化调度模型；S3.以多年来流频率作为所述联合优化调度模型的输入数据，采用多维动态规划算法对上述模型进行求解，得出对应不同离散水位及不同来流年份下梯级系统的总发电量；S4.对每一个来流频率遍历所有离散水位值，找出其中发电量最大对应的离散水位值，作为当前来流频率下的最优消落水位；S5.以来流频率为横坐标，以最优消落水位为纵坐标，点绘各来流频率与其对应的最优消落水位散点图，并基于最小二乘原理，以多项式拟合散点的最佳趋势线，最终得到不同来流频率下的最优消落水位。进一步地，步骤S1中水位离散精度根据计算精度要求设定。进一步地，上述梯级系统联合优化调度模型目标函数为：其中，E为梯级系统在整个调度期的总发电量，T为整个调度期的调度时段数，梯级水库从上游到下游依次编号为1,2,...,N，Ks为第s个电站的出力系数，qts为第s个水库在第t个时段的发电引用流量，Hts为第s个水库在第t个时段的平均水头，Δt为一个调度时段长度。进一步地，上述梯级系统联合优化调度模型调度期初、末水位约束条件为：其中Z0s为第s个水库在第1个时段初的水位，Zbs为第s个水库在整个调度期初的水位，ZTs为第s个水库在最后一个时段末的水位，Zes为第s个水库在整个调度期末的水位，即为步骤S1得到的离散水位值。进一步地，上述联合优化调度模型约束条件还包括水量平衡约束、水库库容约束、下泄流量约束和出力约束。进一步地，步骤S4包括：S4.1.对第1个来流频率P1，遍历所有离散水位值(Z1,Z2,…,Zn)，找出其中发电量最大值Ei,1，将其所对应的离散水位值Zi作为当前来流频率P1下的最优消落水位Z1*；S4.2.对第2个来流频率P2，遍历所有离散水位值(Z1,Z2,…,Zn)，找出其中发电量最大值Ei,2，将其所对应的离散水位值Zi作为当前来水频率P2下的最优消落水位Z2*。S4.3.重复上述过程直至找出所有来流频率Pj(j＝1,2,…,m)发电量最大所对应的最优优落水位Zj*。进一步地，步骤S5包括：S5.1.在可行范围内等步长产生一系列多项式的项数和次数；S5.2.根据每一个多项式的项数和次数，确定一个多项式结构，根据散点数据和最小二乘原理，找出当前结构下的最佳多项式系数及其相关性系数R2；S5.3.以相关性系数R2最大为目标，遍历所有多项式结构，确定多项式最优结构及对应的系数；S5.4.根据所得最优多项式结构，得到不同来流频率下的最优消落水位值。总体而言，通过本专利技术所构思的以上技术方案与现有技术相比，能够取得下列有益效果。(1)本专利技术基于多维动态规划的最优计算结果，以最小二乘原理提取了考虑来流频率的多年调节水库年末消落水位确定规则，相比于多年固定水位的消落方式，所提取的消落规则考虑了来水的年际差异性，能够充分发挥多年调节水库的调节性能，提高梯级系统水能资源转换效率和梯级发电量，对指导含多年调节水库的梯级水库群实际调度运行具有重要意义。(2)本专利技术提取得到的考虑不同来流频率的消落规则，直接通过来流频率确定年末消落水位，不需要其他额外信息，操作性强。且来流频率直接反应了来水量的大小，消落规则的确定以其作为决策指标，能够很好的协调多年调节水库来流频率(水量)与消落水位(水头)的耦合关系，充分发挥梯级水库的整体效益，使梯级系统总发电量与最优情形下的总发电量非常接近。附图说明图1是本专利技术提供的一种考虑来流频率差异的多年调节水库年末消落水位确定方法流程图；图2是本专利技术提供的雅砻江流域中、下游梯级水库地理位置图；图3是本专利技术提供的不同来流频率下的最优消落水位图；图4是本专利技术提供的枯水情形下不同来流频率的最优消落水位图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本专利技术，并不用于限定本专利技术。此外，下面所描述的本专利技术各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。参考图1，本专利技术实施例提供了一种考虑来流频率差异的多年调节水库年末消落水位确定方法，包括：S1.根据设定的水位离散精度，将梯级系统中多年调节水库的死水位至正常蓄水位范围离散为一系列离散水位值；具体地，本专利技术实施例以如图2所示的雅砻江流域中下游梯级水库群为例，对本本文档来自技高网...

【技术保护点】
1.一种考虑来流频率差异的多年调节水库年末消落水位确定方法，其特征在于，包括：/nS1.根据设定的水位离散精度，将梯级系统中多年调节水库的死水位至正常蓄水位范围离散为一系列离散水位值；/nS2.以每一个离散水位值作为多年调节水库的年末消落水位约束，以总发电量最大为目标建立梯级系统联合优化调度模型；/nS3.以多年来流频率作为所述联合优化调度模型的输入数据，采用多维动态规划算法对上述模型进行求解，得出对应不同离散水位及不同来流年份下梯级系统的总发电量；/nS4.对每一个来流频率遍历所有离散水位值，找出其中发电量最大对应的离散水位值，作为当前来流频率下的最优消落水位；/nS5.以来流频率为横坐标，以最优消落水位为纵坐标，点绘各来流频率与其对应的最优消落水位散点图，并基于最小二乘原理，以多项式拟合散点的最佳趋势线，最终得到不同来流频率下的最优消落水位。/n

【技术特征摘要】
1.一种考虑来流频率差异的多年调节水库年末消落水位确定方法，其特征在于，包括：
S1.根据设定的水位离散精度，将梯级系统中多年调节水库的死水位至正常蓄水位范围离散为一系列离散水位值；
S2.以每一个离散水位值作为多年调节水库的年末消落水位约束，以总发电量最大为目标建立梯级系统联合优化调度模型；
S3.以多年来流频率作为所述联合优化调度模型的输入数据，采用多维动态规划算法对上述模型进行求解，得出对应不同离散水位及不同来流年份下梯级系统的总发电量；
S4.对每一个来流频率遍历所有离散水位值，找出其中发电量最大对应的离散水位值，作为当前来流频率下的最优消落水位；
S5.以来流频率为横坐标，以最优消落水位为纵坐标，点绘各来流频率与其对应的最优消落水位散点图，并基于最小二乘原理，以多项式拟合散点的最佳趋势线，最终得到不同来流频率下的最优消落水位。


2.根据权利要求1所述的一种考虑来流频率差异的多年调节水库年末消落水位确定方法，其特征在于，步骤S1中水位离散精度根据计算精度要求设定。


3.根据权利要求1或2所述的一种考虑来流频率差异的多年调节水库年末消落水位确定方法，其特征在于，上述梯级系统联合优化调度模型目标函数为：



其中，E为梯级系统在整个调度期的总发电量，T为整个调度期的调度时段数，梯级水库从上游到下游依次编号为1,2,...,N，Ks为第s个电站的出力系数，qts为第s个水库在第t个时段的发电引用流量，Hts为第s个水库在第t个时段的平均水头，Δt为一个调度时段长度。


4.根据权利要求3所述的一种考虑来流频率差异的多年调节水库年末消落水位确定方法，其特征在于，所述梯级系统联合优化调度模型调度期初、末水位约束条件为：

【专利技术属性】
技术研发人员：蒋志强，胡德超，覃晖，刘懿，冯仲恺，陈璐，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：湖北;42