CN105427052A - 一种基于参考线的并联水库确定性优化调度方法 - Google Patents
一种基于参考线的并联水库确定性优化调度方法 Download PDFInfo
- Publication number
- CN105427052A CN105427052A CN201510895532.0A CN201510895532A CN105427052A CN 105427052 A CN105427052 A CN 105427052A CN 201510895532 A CN201510895532 A CN 201510895532A CN 105427052 A CN105427052 A CN 105427052A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- water
- period
- water storage
- output
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 91
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 201
- 230000007774 longterm Effects 0.000 claims abstract description 4
- 230000008569 process Effects 0.000 claims description 25
- 230000002265 prevention Effects 0.000 claims description 24
- 238000010248 power generation Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 14
- 230000009467 reduction Effects 0.000 claims description 7
- 238000005457 optimization Methods 0.000 description 8
- 239000002351 wastewater Substances 0.000 description 3
- 238000004146 energy storage Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于参考线的并联水库确定性优化调度方法,属于水利电力技术领域;本发明包括以下步骤:(1)确定调度期;(2)建立水库群确定性长期优化调度发电量最大模型;(3)确定水库群初始参考线和初始调度方案;(4)判断初始调度方案是否可行。本发明将水库最小决策限制线和防弃水线结合起来使用,可准确确定水库增发单位出力引起的不蓄电能损失值(K值)计算的起始和终止时刻,不要求各水库使用统一的调度起始和终止时刻,避免了各水库蓄、供水期划分的不同对调度结果的影响;将专家经验和水库调度理论相结合,避免了水库群调度“维数障碍”问题,适用于大规模并联水库群优化调度。
Description
技术领域
本发明属于水利电力技术领域,更具体地,涉及一种基于参考线的并联水库确定性优化调度方法。
背景技术
水库群调度是一个典型的多目标、多阶段、多约束、非线性优化问题,广泛应用的水库供蓄水判别系数法的最主要的缺陷是需人为地将调度期划分为供水期、蓄水期、不蓄不供期等时期,由于蓄、供水期的划分与水库调度周期内的来水过程、水库发电运行方式等因素密切相关,其起止时刻具有不确定性,难以准确确定;传统的优化方法为动态规划及其改进算法,其中改进的算法主要包括逐步优化法、动态规划逐次逼近法、状态逐密动态规划法和增量动态规划法等,与全状态枚举法相比,这些方法虽然大大减小了计算工作量,但却未从根本上克服水库群调度“维数障碍”。
发明内容
针对现有技术中存在的问题,本申请提供的是一种基于参考线的并联水库确定性优化调度方法,其中通过对并联水库的参考线进行研究和设计,在保证调度方案最优性的前提下,极大地缩短了计算时间,因而尤其适用于大规模并联水库群的调度应用场合。
为实现上述目的,本发明提供了1、一种基于参考线的并联水库确定性优化调度方法,其特征在于,所述方法包括以下步骤:
步骤(1),设置调度期,将所述调度期均匀划分为T个等间距的时段,其中T为正整数;设置所述并联水库调度期内的初始条件和约束条件;所述初始条件包括水库调度期期初蓄水状态、期末蓄水状态和调度期内各时段入库流量;所述约束条件包括水库蓄水量约束、出力约束、出库流量约束以及水量平衡方程;
步骤(2),根据所述并联水库中各水库的库容特性曲线、尾水关系曲线、水库优化调度的目标函数和所述约束条件,建立水库群长期优化调度发电量最大模型:
Nm,t=KmQm,tHm,t
其中,E表示水库群全时期最优发电量;Nm,t、Qm,t和Hm,t分别表示t时段水库m的出力、发电流量和发电水头;Km表示水库m的综合出力系数;
步骤(3),确定所述各水库的初始参考线,包括确定所述各水库的防弃水线和最小决策限制线;
步骤(4),以所述各水库的初始参考线为其初始调度线,根据入库流量信息和水量平衡方程从前往后逐时段递推至最后一个时段,计算所述各水库各时段的出库流量和出力值,直至所有水库计算完毕,统计所述并联水库各时段的总出力值;
步骤(5),从前往后逐时段检查所述并联水库各时段的总出力值是否大于保证出力,是,则转下一时段;否,则计算不足出力并按所述并联水库不蓄电能损失最小原则调整所述各水库的时段出力值,调整完毕后根据调整后的所述各水库的蓄水状态按所述步骤(3)所述方法调整后续时段所述各水库的最小决策限制线,调整完毕后转下一时段;上述过程一直持续至调度期末;根据等微增率原理,调整后所述各水库增发单位出力引起的不蓄电能损失值,即K值,应大致相当;以待调整时段t0所述各水库调整出力之后的K值为搜索变量,以t0时段各水库调整之后的出力总和等于保证出力为搜索条件采用一维搜索方法调整所述各水库的初始参考线;
步骤(6),输出所述各水库调整后的参考线和各时段出力值。
优选地,所述步骤(3)包括以下子步骤:
(3-1)以减少弃水能为原则确定所述各水库的防弃水线,其包括以下子步骤:
(3-1-1)已知所述并联水库中的某一水库一个时段的时段末蓄水状态和入库流量信息,令所述时段出库流量等于水轮机最大过流,由水量平衡方程求得所述某一水库所述时段初的蓄水状态;
(3-1-2)根据所述某一水库所述时段初的蓄水状态、时段末的蓄水状态、入库流量和最大容许出力,计算求得所述某一水库的弃水流量,若所述弃水流量大于零,则减小所述某一水库的出库流量使得所述弃水流量为零,并根据所述水量平衡方程重新计算所述时段初的蓄水状态,若所述弃水流量等于零,则所述时段初的蓄水状态取所述步骤(3-1-1)所求得的值,将计算得到的所述时段初的蓄水状态作为上一时段末的蓄水状态;
(3-1-3)按所述步骤(3-1-1)和(3-1-2)从调度期末递推计算至调度期初,得到所述某一水库的一条蓄水过程线;
(3-1-4)对(3-1-3)所述蓄水过程线从调度期初至调度期末逐时段检查是否有弃水发生,若各时段均无弃水,则所述蓄水过程线即为所述某一水库的防弃水线;若某时段有弃水发生,则进行相应调整,使弃水集中在最后一个或若干个时段,调整后的蓄水过程线即为所述某一水库的防弃水线;
(3-1-5)按所述步骤(3-1-1)~(3-1-4)的方法求得所述各水库的防弃水线;
(3-2)根据保持高水位运行以减少不蓄电能损失的原则,确定所述各水库的最小决策限制线,其包括以下子步骤:
(3-2-1)根据所述并联水库中的一个水库某一时段初的蓄水状态和入库流量信息,计算求得所述一个水库所述某一时段的最小出力对应的出库流量,并与最小下泄流量比较,取两者中较大者作为所述一个水库所述某一时段的出库流量,根据水量平衡方程获得所述一个水库所述某一时段末的蓄水状态,并将其作为所述一个水库所述某一时段的下一时段初的蓄水状态;
(3-2-2)按所述步骤(3-2-1)从调度期初递推计算至调度期末,得到一条蓄水过程线,即为所述一个水库的最小决策限制线;
(3-2-3)按所述步骤(3-2-1)和(3-2-2)的方法求得所述各水库的最小决策限制线;
(3-3)对所述一个水库相同时刻的所述防弃水线和所述最小决策限制线的蓄水状态值取小,获得一条蓄水过程线,作为所述一个水库的初始参考线;
(3-4)按步骤(3-3)的方法求得所述各水库的初始参考线。
优选地,所述步骤(5)中所述各水库K值计算包括以下子步骤:
(5-1)确定所述各水库K值计算的起始时段和终止时段;所述各水库K值计算的起始时段取待调整时段t0;所述各水库K值计算的终止时段的确定包括如下步骤:
(5-1-1)计算所述各水库t0时段增发发单位出力后的时段末蓄水状态;根据所述各水库t0时段初的蓄水状态、入库流量和出力值以及水量平衡方程计算所述各水库t0时段末的蓄水状态。
(5-1-2)根据步骤(5-1-1)所得所述各水库t0时段末的蓄水状态,按步骤(3-2)的方法重新计算所述各水库从t0+1时段末至调度期末的最小决策限制线,设其与所述各水库的防弃水线相交于时段t1,则时段t1即为所述各水库K值计算的终止时段;
(5-2)计算所述各水库t0时段增发单位出力后在t0时段至t1时段内总的不蓄电能损失值,即为所述各水库的K值;所述各水库中任一水库某一时段的不蓄电能损失值为所述任一水库所述某一时段的水头减少量与入库流量的乘积。
总体而言,按照本发明的上述技术构思与现有技术相比,主要具备以下的技术优点:
1、本发明将水库最小决策限制线和防弃水线结合使用,可达到在尽量不弃水的情况下尽可能抬高水电站发电水头的效果,提高水电站能量转换效率;
2、本发明所涉及的方法体现了专家调度经验,在保证跳读方案最优性的同时可大大减小计算工作量,且所得结果可理解性强,便于调度人员在实际运用中动态决策;
3、本发明与水库蓄供水判别系数法相比,其优势在于可准确确定水库增发单位出力引起的不蓄电能损失值(K值)计算的终止时刻;
4、本发明与传统的优化方法(动态规划及其改进算法)相比,其优势在于不存在水库群调度“维数障碍”问题,适用于大规模并联水库群优化调度。
附图说明
图1是本发明基于参考线的并联水库确定性优化调度的总体技术路线;
图2是本发明基于参考线的并联水库确定性优化调度流程图;
图3是采用0.618法以不蓄电能损失最小原则将不足出力分配至各水库的流程图;
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
水库发电过程中涉及到的能量有水库初始蓄能、输入能、损耗能和输出能,根据能量守恒定律,在水库初始蓄能和输入能一定的情况下,损耗能越少,则输出的电能越多。基于上述原理,可获得水库发电优化调度须遵循的两条基本原则:(1)尽量避免弃水,以减少弃水能损耗;(2)尽量抬高水库运行水位,以减少不蓄电能损失。根据水库初始蓄水位和计划期内入库流量过程,运用上述原则可确定水库最优运行策略。
图2所示为基于参考线的并联水库确定性优化调度方法的流程图,具体包括以下步骤:
步骤1设置调度期,将调度期(1年)均匀划分为T(T为正整数,在本发明实施例中,T=12或36)个等间距的时段;设置水库调度期内的约束条件,包括蓄水量约束:出力约束:出库流量约束:水量平衡方程:Vm,t+1=Vm,t+(Im,t-Rm,t)·Δt和水库期初、期末水位限制:Zm,1=Zms,Zm,T+1=Zms。其中,V m,t、Vm,t、 N m,t、Nm,t、 R m,t、Rm,t、分别表示t(t=1,2,…,T)时段水库m的最小蓄水量、蓄水量、最大蓄水量、最小出力、出力、最大出力、最小下泄流量、下泄流量和最大下泄流量;Im,t表示m水库t时段的入库流量;Δt表示时段步长;Zms、Zmend分别为水库m期初和期末水位,为给定值。
步骤2根据水库的库容特性曲线、尾水关系曲线、以及水库优化调度的目标函数和上述约束条件,以水库各时段入库流量作为输入,建立水库群长期优化调度发电量最大模型如下:
Nm,t=KmQm,tHm,t
其中,E表示水库群全时期最优发电量;Nm,t、Qm,t和Hm,t分别表示t时段水库m的出力、发电流量和发电水头;Km表示水库m的综合出力系数。
步骤3确定水库群初始参考线,包括确定各水库防弃水线和最小决策限制线,优选地,步骤3包括以下子步骤:
(3-1)确定水库防弃水线。水库运行过程中,减少弃水有利于提高发电量。水库用于发电的下泄水量在以下两种情况下会形成弃水:一是下泄流量大于水库水轮机最大过流量PlantsMaxQm;二是水库按PlantsMaxQm下泄,其出力值大于时段最大出力限制具体包括以下子步骤:
(3-1-1)已知水库m计划期末水位Zm,T+1和时段入库流量Im,T,由库容特性线性插值得其时段末库容Vm,T+1,令所述时段其出库流量Rm,T等于水库水轮机最大过流PlantsMaxQm,由Vm,T+1、Rm,T、Im,T和水量平衡方程,可求得水库时段初库容Vm,T,若则取若Vm,T<V m,T,则取Vm,T=V m,T,求得时段初库容Vm,T后由库容特性线性插值获得时段初水位Zm,T;
(3-1-2)由Vm,T、Vm,T+1、Im,T、PlantsMaxQm、和水量平衡方程可求得弃水流量若则调整出库流量为根据水量平衡方程重新计算Vm,T,若则取若Vm,T<V m,T,则取Vm,T=V m,T,求得Vm,T后由库容特性线性插值获得时段初水位Zm,T;若则所述时段初水位Zm,T取上述步骤(3-1-1)所获得的值;将Zm,T作为上一时段的时段末水位;
(3-1-3)按照步骤(3-1-1)和步骤(3-1-2)从调度期末递推计算至调度期初,得到一条水位过程线(t=1,2,…,T,T+1),相应的库容为(t=1,2,…,T,T+1);
(3-1-4)对步骤(3-1-3)中所述水位过程线从调度期初至调度期末逐时段检查是否有弃水发生。若各时段均无弃水发生,则所述水位过程线即为水库防弃水线;若时段t1有弃水发生,则进行相应调整,调整原则是使弃水集中在最后一个或若干个时段,以保持水库高水位发电,减小不蓄电能损失。具体调整方法如下:调整水库t1时段出库流量等于水库水轮机最大过流PlantsMaxQm,由和水量平衡方程,可求得水库时段末库容若则取 若 则取 由PlantsMaxQm、和水量平衡方程可求得弃水流量若则调整出库流量为根据水量平衡方程重新计算若 则取 若 则取 求得后由库容特性线性插值获得时段末水位并计算所述时段较调整前增加的蓄水量后续时段水库末库容均增加若遇t2(t1<t2<T)时段末库容大于最大库容的情况,则令计算所述时段增加的蓄水量后续时段水库末库容在的基础上均增加调整后的水位过程线即为防弃水线;
(3-1-5)按所述步骤(3-1-1)~(3-1-4)的方法求得所有水库的防弃水线;
(3-2)确定最小决策限制线。水电站运行过程应尽量抬高水位发电,以减小不蓄电能损失,增发电量,具体包括以下子步骤:
(3-2-1)已知水库调度期初水位Zm,1和时段入库流量Im,1,计算所述时段水库按最小出力N m,1运行对应的出库流量并与所述时段容许的最小下泄流量R m,1比较,取两者中取较大者作为所述时段出库流量Rm,1,即所述时段出库流量根据水量平衡方程,计算所述时段末库容Vm,2,若则取若则取Vm,2=V m,2,再根据库容特性线性插值得所述时段末水位Zm,2,将其作为下一时段的时段初水位;
(3-2-2)按步骤(3-2-1)从调度期初递推计算至调度期末,得到一条水位过程线(t=1,2,…,T,T+1),即为水库m最小决策限制线。
(3-2-3)按所述步骤(3-2-1)和(3-2-2)的方法求得所有水库的最小决策限制线;
(3-3)对水库相同时刻的防弃水线和最小决策限制线的水位值取小,得到水库初始参考线(t=1,2,…,T,T+1),即
(3-4)按步骤(3-3)的方法求得所有水库初始参考线;
步骤4以各水库初始参考线为其初始调度线,根据各水库各时段入库流量Im,t和水量平衡方程从前往后逐时段计算各水库各时段出库流量Rm,t和出力Nm,t,直至计算完调度期最后一个时段为止。若时段t遇则取所有水库计算完毕后统计水库群各时段总出力值
步骤5从前往后逐时段检查水库群时段总出力值是否大于保证出力Nb,是,则转下一时段;否,则计算不足出力Nbz并按水库群不蓄电能损失最小原则将不足出力Nbz分配至各水库,调整完毕后根据调整后的水库时段末蓄水状态按步骤(3-2)所述方法调整后续时段各水库的最小决策限制线,调整完毕后转下一时段。上述过程一直持续至调度期末。根据等微增率原理,调整后各水库增发单位出力引起的不蓄电能损失值(K值)应大致相当。采用0.618法调整水库群初始参考线,以待调整时段t0各水库调整出力之后的K值为搜索变量,以待调整时段t0各水库调整后的出力总和等于保证出力为搜索条件,优选地,包括以下步骤:
(5-1)确定K值的下限值。包括以下步骤:
(5-1-1)计算各水库增发单位出力引起的不蓄电能损失值Km,包括以下步骤:
(5-1-1-1)将水库待调整时段t0的出力增加一个单位值,计算出时段末水位,作为下一时段的时段初水位,按所述步骤(3-2)的方法计算待调整时段t0至调度期最后一个时段T的水库最小决策限制线设其与所述水库防弃水线相交于时刻t1,记t1时刻所在时段为t2。
(5-1-1-2)计算所述t0~t2时期内水库各时段水头减少量ΔHm,t和不蓄电能损失值Km。水库m时段t的水位减少量应为调整前后水库平均水位之差,即 水库m时段t内不蓄电能损失值Km,t应为所述时段入库流量Im,t与水头减少量ΔHm,t之乘积,即Km,t=Im,t·ΔHm,t。
(5-1-1-3)按所述步骤(5-1-1-1)和(5-1-1-2)的方法计算所有水库增发单位出力引起的不蓄电能损失值Km。
(5-1-2)K值的下限值取各水库增发单位出力引起的不蓄电能损失值中的最小者,即Kmin=min(Km),其中(m=1,2,…,M)。
(5-2)确定K值的上限值。包括以下步骤:
(5-2-1)计算各水库减发单位出力引起的不蓄电能增加值K'm,包括以下步骤:
(5-2-1-1)取水库待调整时段t0末的水位为水库死水位根据水量平衡方程计算所述水库待调整时段t0的出库流量Rm,t、出力Nm,t和弃水流量若则减小出库流量Rm,t,使弃水流量并计算时段末水库水位值Zm,t+1和时段出力值Nm,t;否则,保持出库流量Rm,t、时段末水位Zm,t+1和时段出力值Nm,t不变。按所述步骤(3-2)的方法计算待调整时段t0至调度期最后一个时段T的最小决策限制线设其与所述水库防弃水线相交于时刻t3,记t3时刻所在时段为t4。
(5-2-1-2)将水库待调整时段t0的出力减小一个单位值,计算出时段末水位Zm,t+1,作为下一时段初水位,按所述步骤(3-2)的方法计算所述t0~T时期内的水库最小决策限制线与所述水库防弃水线相交于时刻t5,记t5时刻所在时段为t6。
(5-2-1-3)计算所述t0~t6时期内水库各时段水头增加量ΔH'm,t和不蓄电能增加值K'm。水库m时段t的水位增加量应为调整后与调整前水库平均水位之差,即 水库m时段t内不蓄电能增加值K'm,t应为所述时段入库流量Im,t与水头增加量ΔH'm,t之乘积,即K'm,t=Im,t·ΔH'm,t。
(5-2-1-4)按所述步骤(5-2-1-1)和(5-2-1-3)的方法计算所有水库增发单位出力引起的不蓄电能增加值K'm。
(5-2-2)取K值的上限值为各水库减发单位出力引起的不蓄电能增加值K'm中的最大者,即Kmax=max(K'm),其中(m=1,2,…,M)。
(5-3)采用0.618法在闭区间[Kmin,Kmax]内搜索得到一个K*,使待调整时段t0各水库调整出力之后的K值与K*之差的绝对值最小,即min(K-K*),且使水库群出力总和等于保证出力Nb,即包括以下步骤:
(5-3-1)在闭区间[Kmin,Kmax]内取K1=Kmin+0.382(Kmax-Kmin),各水库以待调整时段t0调整出力之后的K值与K1之差的绝对值最小为目标,即以min(K-K1)为目标,采用0.618法搜索得到其出力值,所有水库搜索完毕后计算待调整时段t0水库群的总出力值,记为N1,并保存各水库调整之后的参考线取K2=Kmin+0.618(Kmax-Kmin),各水库以待调整时段t0调整出力之后的K值与K2之差的绝对值最小为目标,即以min(K-K2)为目标,采用0.618法搜索得到其出力值,所有水库搜索完毕后计算待调整时段t0水库群的总出力值,记为N2,并保存各水库调整之后的参考线若|N1-Nb|<|N2-Nb|,则令Kmax=K2,否则令Kmin=K1。
所述步骤(5-3-1)中采用0.618法搜索水库出力值,包括以下步骤(以水库m和K=K1为例):
(5-3-1-1)确定水库m出力上限和下限下限值取初始方案中水库待调整时段t0的出力,上限值取所述水库所述时段最大容许出力
(5-3-1-2)在闭区间内在取并令水库待调整时段t0的出力等于N3,即计算时段末水库水位值并按(5-1-1)所述方法计算所述水位值下水库增发单位出力引起的不蓄电能损失值K3;取并令水库待调整时段t0的出力等于N4,即计算时段末水库水位值并按(5-1-1)所述方法计算所述水位值下水库增发单位出力引起的不蓄电能损失值K4。若|K3-K1|<|K4-K1|,则令否则令
(5-3-1-3)重复步骤(5-3-1-2),直至满足精度要求为止,即输出水库m调整之后的出力值N'm,t;
(5-4)重复步骤(5-3),直至满足精度要求为止,即|Kmax-Kmin|<ε',输出K*;
步骤(6)输出各水库调整后的出力值N'm,t和调整后参考线
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于参考线的并联水库确定性优化调度方法,其特征在于,所述方法包括以下步骤:
步骤(1),设置调度期,将所述调度期均匀划分为T个等间距的时段,其中T为正整数;设置所述并联水库调度期内的初始条件和约束条件;所述初始条件包括水库调度期期初蓄水状态、期末蓄水状态和调度期内各时段入库流量;所述约束条件包括水库蓄水量约束、出力约束、出库流量约束以及水量平衡方程;
步骤(2),根据所述并联水库中各水库的库容特性曲线、尾水关系曲线、水库优化调度的目标函数和所述约束条件,建立水库群长期优化调度发电量最大模型:
Nm,t=KmQm,tHm,t
其中,E表示水库群全时期最优发电量;Nm,t、Qm,t和Hm,t分别表示t时段水库m的出力、发电流量和发电水头;Km表示水库m的综合出力系数;
步骤(3),确定所述各水库的初始参考线,包括确定所述各水库的防弃水线和最小决策限制线;
步骤(4),以所述各水库的初始参考线为其初始调度线,根据入库流量信息和水量平衡方程从前往后逐时段递推至最后一个时段,计算所述各水库各时段的出库流量和出力值,直至所有水库计算完毕,统计所述并联水库各时段的总出力值;
步骤(5),从前往后逐时段检查所述并联水库各时段的总出力值是否大于保证出力,是,则转下一时段;否,则计算不足出力并按所述并联水库不蓄电能损失最小原则调整所述各水库的时段出力值,调整完毕后根据调整后的所述各水库的蓄水状态按所述步骤(3)所述方法调整后续时段所述各水库的最小决策限制线,调整完毕后转下一时段;上述过程一直持续至调度期末;根据等微增率原理,调整后所述各水库增发单位出力引起的不蓄电能损失值,即K值,应大致相当;以待调整时段t0所述各水库调整出力之后的K值为搜索变量,以t0时段各水库调整之后的出力总和等于保证出力为搜索条件采用一维搜索方法调整所述各水库的初始参考线;
步骤(6),输出所述各水库调整后的参考线和各时段出力值。
2.如权利要求1所述的方法,其特征在于,所述步骤(3)包括以下子步骤:
(3-1)以减少弃水能为原则确定所述各水库的防弃水线,其包括以下子步骤:
(3-1-1)已知所述并联水库中的某一水库一个时段的时段末蓄水状态和入库流量信息,令所述时段出库流量等于水轮机最大过流,由水量平衡方程求得所述某一水库所述时段初的蓄水状态;
(3-1-2)根据所述某一水库所述时段初的蓄水状态、时段末的蓄水状态、入库流量和最大容许出力,计算求得所述某一水库的弃水流量,若所述弃水流量大于零,则减小所述某一水库的出库流量使得所述弃水流量为零,并根据所述水量平衡方程重新计算所述时段初的蓄水状态,若所述弃水流量等于零,则所述时段初的蓄水状态取所述步骤(3-1-1)所求得的值,将计算得到的所述时段初的蓄水状态作为上一时段末的蓄水状态;
(3-1-3)按所述步骤(3-1-1)和(3-1-2)从调度期末递推计算至调度期初,得到所述某一水库的一条蓄水过程线;
(3-1-4)对(3-1-3)所述蓄水过程线从调度期初至调度期末逐时段检查是否有弃水发生,若各时段均无弃水,则所述蓄水过程线即为所述某一水库的防弃水线;若某时段有弃水发生,则进行相应调整,使弃水集中在最后一个或若干个时段,调整后的蓄水过程线即为所述某一水库的防弃水线;
(3-1-5)按所述步骤(3-1-1)~(3-1-4)的方法求得所述各水库的防弃水线;
(3-2)根据保持高水位运行以减少不蓄电能损失的原则,确定所述各水库的最小决策限制线,其包括以下子步骤:
(3-2-1)根据所述并联水库中的一个水库某一时段初的蓄水状态和入库流量信息,计算求得所述一个水库所述某一时段的最小出力对应的出库流量,并与最小下泄流量比较,取两者中较大者作为所述一个水库所述某一时段的出库流量,根据水量平衡方程获得所述一个水库所述某一时段末的蓄水状态,并将其作为所述一个水库所述某一时段的下一时段初的蓄水状态;
(3-2-2)按所述步骤(3-2-1)从调度期初递推计算至调度期末,得到一条蓄水过程线,即为所述一个水库的最小决策限制线;
(3-2-3)按所述步骤(3-2-1)和(3-2-2)的方法求得所述各水库的最小决策限制线;
(3-3)对所述一个水库相同时刻的所述防弃水线和所述最小决策限制线的蓄水状态值取小,获得一条蓄水过程线,作为所述一个水库的初始参考线;
(3-4)按步骤(3-3)的方法求得所述各水库的初始参考线。
3.如权利要求1所述的方法,其特征在于,所述步骤(5)中所述各水库K值计算包括以下子步骤:
(5-1)确定所述各水库K值计算的起始时段和终止时段;所述各水库K值计算的起始时段取待调整时段t0;所述各水库K值计算的终止时段的确定包括如下步骤:
(5-1-1)计算所述各水库t0时段增发发单位出力后的时段末蓄水状态;根据所述各水库t0时段初的蓄水状态、入库流量和出力值以及水量平衡方程计算所述各水库t0时段末的蓄水状态。
(5-1-2)根据步骤(5-1-1)所得所述各水库t0时段末的蓄水状态,按步骤(3-2)的方法重新计算所述各水库从t0+1时段末至调度期末的最小决策限制线,设其与所述各水库的防弃水线相交于时段t1,则时段t1即为所述各水库K值计算的终止时段;
(5-2)计算所述各水库t0时段增发单位出力后在t0时段至t1时段内总的不蓄电能损失值,即为所述各水库的K值;所述各水库中任一水库某一时段的不蓄电能损失值为所述任一水库所述某一时段的水头减少量与入库流量的乘积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510895532.0A CN105427052B (zh) | 2015-12-08 | 2015-12-08 | 一种基于参考线的并联水库确定性优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510895532.0A CN105427052B (zh) | 2015-12-08 | 2015-12-08 | 一种基于参考线的并联水库确定性优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105427052A true CN105427052A (zh) | 2016-03-23 |
CN105427052B CN105427052B (zh) | 2019-01-22 |
Family
ID=55505245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510895532.0A Expired - Fee Related CN105427052B (zh) | 2015-12-08 | 2015-12-08 | 一种基于参考线的并联水库确定性优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105427052B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825437A (zh) * | 2016-03-24 | 2016-08-03 | 大连理工大学 | 一种复杂水库群共同供水任务分配方法 |
CN106320257A (zh) * | 2016-08-24 | 2017-01-11 | 武汉大学 | 基于水文观测的湖库槽蓄曲线确定方法 |
CN106873372A (zh) * | 2017-03-22 | 2017-06-20 | 中国水利水电科学研究院 | 基于防洪调度数据自适应控制的水库防洪调度优化方法 |
CN107958306A (zh) * | 2017-11-22 | 2018-04-24 | 华中科技大学 | 一种基于参考线的水电站随机优化调度方法 |
CN108564213A (zh) * | 2018-04-10 | 2018-09-21 | 中国水利水电科学研究院 | 一种基于gpu加速的并联水库群防洪优化调度方法 |
CN109886543A (zh) * | 2019-01-16 | 2019-06-14 | 河南大学 | 基于均匀设计的水库优化调度方法 |
CN112232659A (zh) * | 2020-10-15 | 2021-01-15 | 华中科技大学 | 一种梯级水库发电调度方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537577A (zh) * | 2014-12-25 | 2015-04-22 | 华中科技大学 | 一种基于参考线的水电站确定性优化调度方法 |
-
2015
- 2015-12-08 CN CN201510895532.0A patent/CN105427052B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104537577A (zh) * | 2014-12-25 | 2015-04-22 | 华中科技大学 | 一种基于参考线的水电站确定性优化调度方法 |
Non-Patent Citations (1)
Title |
---|
谢维: "水电站(群)优化调度与运行规则研究", 《中国博士学位论文全文数据库工程科技II辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105825437A (zh) * | 2016-03-24 | 2016-08-03 | 大连理工大学 | 一种复杂水库群共同供水任务分配方法 |
CN105825437B (zh) * | 2016-03-24 | 2019-04-19 | 大连理工大学 | 一种复杂水库群共同供水任务分配方法 |
CN106320257A (zh) * | 2016-08-24 | 2017-01-11 | 武汉大学 | 基于水文观测的湖库槽蓄曲线确定方法 |
CN106873372A (zh) * | 2017-03-22 | 2017-06-20 | 中国水利水电科学研究院 | 基于防洪调度数据自适应控制的水库防洪调度优化方法 |
CN106873372B (zh) * | 2017-03-22 | 2018-05-11 | 中国水利水电科学研究院 | 基于防洪调度数据自适应控制的水库防洪调度优化方法 |
CN107958306A (zh) * | 2017-11-22 | 2018-04-24 | 华中科技大学 | 一种基于参考线的水电站随机优化调度方法 |
CN108564213A (zh) * | 2018-04-10 | 2018-09-21 | 中国水利水电科学研究院 | 一种基于gpu加速的并联水库群防洪优化调度方法 |
CN108564213B (zh) * | 2018-04-10 | 2022-05-13 | 中国水利水电科学研究院 | 一种基于gpu加速的并联水库群防洪优化调度方法 |
CN109886543A (zh) * | 2019-01-16 | 2019-06-14 | 河南大学 | 基于均匀设计的水库优化调度方法 |
CN109886543B (zh) * | 2019-01-16 | 2022-09-16 | 河南大学 | 基于均匀设计的水库优化调度方法 |
CN112232659A (zh) * | 2020-10-15 | 2021-01-15 | 华中科技大学 | 一种梯级水库发电调度方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105427052B (zh) | 2019-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105427052B (zh) | 一种基于参考线的并联水库确定性优化调度方法 | |
CN103942728B (zh) | 梯级水电站群日发电计划编制方法 | |
CN104537577B (zh) | 一种基于参考线的水电站确定性优化调度方法 | |
CN109345010B (zh) | 一种梯级泵站的多目标优化调度方法 | |
CN107248015B (zh) | 基于集合预报的水库实时蓄水调度方法 | |
CN109447405B (zh) | 一种承担调峰任务的一库多级式梯级库群短期计划制定方法 | |
CN110059345B (zh) | 一种抽水蓄能机组相继甩负荷关机规律优化方法与系统 | |
CN102182159A (zh) | 一种梯级水库汛限水位联合运用调度方法 | |
CN108710972B (zh) | 一种基于poa算法的水库防洪调度方法 | |
CN104965954B (zh) | 一种基于全特性空间曲面建模的梯级负荷经济分配方法 | |
CN110598983A (zh) | 一种自适应改进粒子群算法的梯级水库优化调度方法 | |
CN116227863B (zh) | 一种基于哈里斯鹰优化算法的梯级水电站优化调度方法 | |
CN114358492B (zh) | 水电站水库调度的确定方法 | |
CN107506894A (zh) | 一种考虑非恒定耦合约束的水电群调度方法 | |
CN111864811B (zh) | 一种考虑电网送出需求的水光互补系统光伏容量配置方法 | |
CN107958306B (zh) | 一种基于参考线的水电站随机优化调度方法 | |
CN106300441A (zh) | 一种面向电网多级输电断面的水电站机组负荷分配方法 | |
CN113363976B (zh) | 一种基于场景图的风光水互补发电系统中期优化调度方法 | |
CN115409234A (zh) | 一种基于混合算法的梯级水电站优化调度模型求解方法 | |
CN112766565B (zh) | 梯级水库群防洪优化调度分数阶逐步优化方法 | |
CN109635999A (zh) | 一种基于粒子群-细菌觅食的水电站调度方法及系统 | |
CN105976058B (zh) | 一种基于灰色系统预测的离散微分动态规划方法 | |
CN108181938B (zh) | 一种逆时消纳水电站弃水的水位调整方法 | |
CN114362242A (zh) | 基于双层优化的风水火电力系统多目标动态调度方法 | |
CN104635769B (zh) | 一种梯级水电站群水电站汛期最优水位范围控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190122 |
|
CF01 | Termination of patent right due to non-payment of annual fee |