CN108388954B - 一种基于随机安全域的梯级水电鲁棒优化调度方法 - Google Patents

一种基于随机安全域的梯级水电鲁棒优化调度方法 Download PDF

Info

Publication number
CN108388954B
CN108388954B CN201810012187.5A CN201810012187A CN108388954B CN 108388954 B CN108388954 B CN 108388954B CN 201810012187 A CN201810012187 A CN 201810012187A CN 108388954 B CN108388954 B CN 108388954B
Authority
CN
China
Prior art keywords
hydropower
hydropower station
robust
water
model
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.)
Active
Application number
CN201810012187.5A
Other languages
English (en)
Other versions
CN108388954A (zh
Inventor
葛晓琳
钟俊玲
金言
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai University of Electric Power
Original Assignee
Shanghai University of Electric Power
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai University of Electric Power filed Critical Shanghai University of Electric Power
Priority to CN201810012187.5A priority Critical patent/CN108388954B/zh
Publication of CN108388954A publication Critical patent/CN108388954A/zh
Application granted granted Critical
Publication of CN108388954B publication Critical patent/CN108388954B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Educational Administration (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)将鲁棒两阶段模型中第二阶段的再调度模型转化为max‑min的数学优化问题,并对预调度方案的鲁棒可行性进行判定,通过反馈修正协调优化最终获得具有鲁棒性的调度方案。与现有技术相比,本发明具有精确可靠、考虑来水时序相关性、鲁棒性和经济性相协调等优点。

Description

一种基于随机安全域的梯级水电鲁棒优化调度方法
技术领域
本发明涉及梯级水电调度领域,尤其是涉及一种基于随机安全域的梯级水电鲁棒优化调度方法。
背景技术
跨年梯级水电优化调度的主要任务是以两年为调度周期,月为调度时段,结合两年所有可能来水情况,在满足各种约束条件下通过优化水力发电和水库调度计划。相对于短期优化调度,跨年的梯级水电调度周期更长,可优化空间更大,一直以来都是研究的热点问题。
受自然因素的影响,跨年水电调度周期内水库来水是随机变化的,为此假定来水为固定值的确定性调度方法适用性有限,而考虑来水不确定性的随机优化调度方法更具实用性。现有方法中大多根据假定的概率分布函数或者多场景技术来模拟来水的不确定性。然而中长期来水的概率分布函数通常都是难以获得的,此外为了保证调度方案的可靠性需要考虑大量场景,计算效率偏低,实际应用效果受限。当难以获得不确定参数的概率分布函数时,鲁棒优化法提供了更加可靠的调度策略。不同于随机规划法,鲁棒优化通过不确定集来描述不确定性,即将不确定因素的所有可能的实现情况划定在集合中。鲁棒优化理论目前已初步应用于随机优化问题的求解。然而现有文献中鲁棒优化理论多用于处理风电的不确定性,针对梯级水电来水不确定性的应用较少涉及。梯级水电之间的水电耦合关系会进一步增加鲁棒优化理论应用的难度。
鲁棒优化是通过应对不确定集中的最恶劣情况,来保证所有可能情况的可行性。当不确定集仅用一个波段描述时,其最恶劣情况为所有参数均处于最极端的情况,然而由于不确定参数之间的内在联系,所有参数同时处于极端情况的发生概率极小甚至不存在,考虑到梯级水电站水库来水是连续,非平稳的随机过程,表现出或强或弱的时序相关性,直接采用传统的不确定集构建方法会导致不确定集的不准确性。
因此,急需一种新的梯级水电鲁棒调度方法,更加准确和严谨的描述水库来水的内在特性和时序相关性,以获取安全性和经济性相协调的鲁棒调度方案。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于随机安全域的梯级水电鲁棒优化调度方法。
本发明的目的可以通过以下技术方案来实现:
一种基于随机安全域的梯级水电鲁棒优化调度方法,包括以下步骤:
1)根据中长期梯级水电来水的不确定性及其时序相关性建立随机安全域,即考虑来水时序相关性的多波段不确定集合;
2)根据多波段不确定集合构建经济性和安全性协调的鲁棒两阶段模型;
3)求解鲁棒两阶段模型中第一阶段的跨年梯级水电预调度模型,得到水电站的预调度方案和预测场景下的水库库容;
4)将鲁棒两阶段模型中第二阶段的再调度模型转化为max-min的数学优化问题,并对预调度方案的鲁棒可行性进行判定,通过反馈修正协调优化最终获得具有鲁棒性的调度方案。
所述的步骤1)具体包括以下步骤:
11)根据中长期梯级水电来水的随机特性,建立来水不确定集,并将来水不确定集划分为B段,得到来水多波段不确定集;
12)基于来水的时序分布特性以及流域相关性,建立来水时序相关性约束;
13)获取考虑来水时序相关性的多波段不确定集合Rf,即:
Figure BDA0001540688820000021
Figure BDA0001540688820000022
Figure BDA0001540688820000023
Figure BDA0001540688820000024
式中:h为水电站编号,t为时段编号,b为波段编号,b∈[1,B],B为波段总数。Rh,t为水电站h在时段t的预测均值,
Figure BDA0001540688820000025
分别为水电站h在时段t来水不确定第b波段子集正负偏离端点值,
Figure BDA0001540688820000026
均为二进制变量,
Figure BDA0001540688820000027
取值为1时为不确定变量落入第b正偏离波段,取值为0为未落入波段,
Figure BDA0001540688820000031
取值为1时为不确定变量落入第b负偏离波段,取值为0为未落入波段,
Figure BDA0001540688820000032
为水电站h在时段t的来水不确定值,πb为第b波段的比重系数,即来水落入第b波段的占比,其比重系数之和为1,ΔH1为不确定预算,表示所有
Figure BDA0001540688820000033
求和的上限,取值范围为[0,NT*NH],其中NT为调度时段数,NH为水电站数,
Figure BDA0001540688820000034
为向上取整,
Figure BDA0001540688820000035
为比重系数,ΔH2为不确定预算。
所述的步骤2)中,鲁棒两阶段模型的第一阶段为跨年梯级水电预调度模型,其表达式为:
Figure BDA0001540688820000036
ph,j,t=Ah,jqh,j,tHh,t
其中,f为调度周期内梯级水电系统的总发电效益,NH为水电站总数,T为一年的时段总数,Jh为水电站h中的机组总数目,ζh,t为水电站h在t时段的市场电价,ph,j,t为水电站h机组j在t时段的发电功率,Δt为时间尺度,Ah,j为水电站h机组j的出力系数,qh,j,t为水电站h机组j在t时段的发电流量,Hh,t为水电站h时段t的平均水头。
所述的跨年梯级水电预调度模型的约束条件包括:
1、水量平衡约束:
Figure BDA0001540688820000037
Figure BDA0001540688820000038
其中,vh,t为水电站h在t时段的库容,vh,t-1为水电站h在t-1时段的库容,qh,t为水电站h在t时段的发电流量,sh,t为水电站h在t时段的弃水流量,Rh,t为水电站h在t时段的来水的预测值,m为水电站编号,Uh为水电站h的上游水电站集合,qm,t和sm,t分别为上游电站的发电流量和弃水流量,Δt′为时段t所包含的秒数;
2、库容约束:
Figure BDA0001540688820000039
其中,
Figure BDA00015406888200000310
为水电站h的库容最小值和最大值;
3、弃水约束:
Figure BDA00015406888200000311
其中,
Figure BDA0001540688820000041
分别为水电站h弃水流量的最小值和最大值;
4、发电流量约束:
Figure BDA0001540688820000042
其中,
Figure BDA0001540688820000043
分别为水电站h发电流量的最小值和最大值;
5、机组出力约束:
Figure BDA0001540688820000044
其中,
Figure BDA0001540688820000045
为水电站h中机组j的出力最小值和最大值;
6、初始期末库容约束:
vh,0=vh,ini
vh,2T=vh,term
其中,vh,ini为水电站h第一年年初库容,vh,term为水电站h第二年年末库容,vh,0为水电站h期初库容,vh,2T为水电站h期末库容。
所述的步骤3)中具体包括以下步骤:
设定来水为预测值,采用混合整数线性规划法求解跨年梯级水电预调度模型,获取水电站的预调度方案以及在预测场景下的水库库容
Figure BDA0001540688820000046
所述的步骤4)中,第二阶段再调度模型的表达式为:
Figure BDA0001540688820000047
其中,z为博弈均衡值,λ1,h,t、λ2,h,t、λ3,h,t、λ4,h,t、λ5,h,t、λ6,h,t为相应约束条件中引入的松弛变量。
所述的第二阶段再调度模型的约束条件包括:
Figure BDA0001540688820000048
Figure BDA0001540688820000049
Figure BDA00015406888200000410
Figure BDA00015406888200000411
Figure BDA00015406888200000412
Figure BDA00015406888200000413
Figure BDA00015406888200000414
Figure BDA00015406888200000415
Figure BDA00015406888200000416
vh,0=vh,ini
vh,2T=vh,term
其中,
Figure BDA0001540688820000051
为不确定来水下的库容,
Figure BDA0001540688820000052
分别为t时段水库调节的上下限值,
Figure BDA0001540688820000053
为预测场景下的水库库容,m表示水电站h的上游水电站编号,qm,t和sm,t分别为上游电站的发电流量和弃水流量。
所述的步骤4)中,将max-min的数学优化问题转换为max问题,并采用混合整数二次线性规划法求解,若博弈均衡值z≤ε,ε为取值为0.001,则判定预调度方案能够满足所有可能来水情景下的安全运行约束,若博弈均衡值z>ε,则返回预调度模型中修正预调度方案,通过迭代直至满足收敛条件。
所述的返回预调度模型中修正预调度方案,按照下式进行:
Figure BDA0001540688820000054
其中,u为再调度模型的对偶变量,z为博弈均衡值,x为预调度变量,
Figure BDA0001540688820000056
为预调度中传递到下层的变量,即预测场景下的水库库容
Figure BDA0001540688820000055
本梯级水电鲁棒优化调度方法的调度周期为2年。
与现有技术相比,本发明具有以下优点:
一、精确可靠:本发明公开的方法通过多波段不确定集能够精确可靠地描述中长期梯级水电来水特性。
二、考虑来水时序相关性:本发明公开的方法对来水的时序相关性和流域相关性进行了建模,这使得模型更加符合实际。
三、鲁棒性和经济性相协调:本发明公开的方法保证了两年内所有可能来水场景都能满足安全经济运行的需求,建立了两阶段的鲁棒调度模型,第一阶段为跨年梯级水电预调度模型,第二阶段为安全校核的再调度模型,两阶段通过反馈修正协调优化最终获得具有鲁棒性的调度方案。
附图说明
图1为水电站预测信息,其中,图(1a)为1#水电站的预测信息,图(1b)为2#水电站的预测信息。
图2为两种模型的发电效益与不确定预算关系曲线图。
图3两种模型的跨年水库调度计划。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
本发明提出了一种基于随机安全域的梯级水电鲁棒优化调度方法,本发明首先建立了来水多波段不确定集约束,以往的模型中,针对中长期梯级水电来水的不确定性,大多根据假定的概率分布函数或者多场景技术来模拟来水的不确定性。然而中长期来水的概率分布函数通常都是难以获得的,此外为了保证调度方案的可靠性需要考虑大量场景,计算效率偏低,实际应用效果受限。故提出了多波段不确定集描述来水的随机特性和时序相关性。
接着,构建安全性和经济性协调的鲁棒两阶段模型,具体步骤如下:
步骤1:以跨年整体效益最大化为目标构建预调度模型;
步骤2:在预调度中假定来水为预测值,进而利用混合整数线性规划法求解模型,得到水电站的预调度方案以及预测场景下的水库库容
Figure BDA0001540688820000061
步骤3:结合来水多波段不确定集构建安全校核的再调度max-min模型。
步骤4:基于对偶理论将max-min模型转化为max问题;
步骤5:利用混合整数二次线性规划法求解上述模型得到博弈均衡值z。若z≤ε(ε为一个很小的正数)则预调度方案能够满足所有可能来水情景下的安全运行约束,否则,将返回修正约束到预调度模型中修正预调度方案,通过迭代直至满足收敛条件,得到相应的调度方案;
实施例1:
本发明提出一种基于随机安全域的梯级水电鲁棒优化调度方法,包括以下步骤:
1)针对来水随机性,采用多波段不确定集描述来水的不确定性和时序相关性,针对中长期梯级水电来水的随机特性和时序相关性建立来水多波段不确定集,具体表达式如下:
Figure BDA0001540688820000062
Figure BDA0001540688820000063
Figure BDA0001540688820000064
Figure BDA0001540688820000065
式中:h为水电站编号,t为时段编号,b为波段编号,b∈[1,B],B为波段总数。Rh,t为水电站h在时段t的预测均值。
Figure BDA0001540688820000066
分别表示水电站h在时段t来水不确定第b波段子集正负偏离端点值,
Figure BDA0001540688820000071
为二进制变量,取值为1时表示不确定变量落入第b正/负偏离波段。
Figure BDA0001540688820000072
为水电站h在时段t的来水不确定值。πb为第b波段的比重系数,即来水落入第b波段的占比,比重系数之和为1;ΔH1为不确定预算,表示所有
Figure BDA0001540688820000073
求和的上限,取值范围为[0,NT*NH],其中NT表示调度时段数,NH为水电站数。第三式表示水电站在所有时段上的总的偏离程度是有一定范围的,被πb和不确定预算ΔH1所限制。
Figure BDA0001540688820000074
表示向上取整,
Figure BDA0001540688820000075
取值越大,表示来水偏离其预测值的场景越多。第四式表示在波段b上,整个调度周期内水电站h的偏离程度。由于来水具有时序相关性,同时在整个调度周期内偏离程度是在一定范围内的,为此引入波段的比重系数
Figure BDA0001540688820000076
和不确定预算ΔH2来具体描述来水时序相关性对模型带来的影响;
2)以梯级水电的整体效益最大化为目标建立梯级水电跨年预调度模型,结合来水的预测信息利用混合整数二次规划法进行求解,得到相应预调度方案和预测场景下的水库库容,跨年梯级水电预调度模型是以梯级水电的整体发电效益最大化为目标:
Figure BDA0001540688820000077
其中,f为调度周期内梯级水电系统的总发电效益;NH为水电站总数;T为一年的时段总数;Jh为水电站h中的机组总数目;ζh,t为水电站h在t时段的市场电价,由于我国目前是垂直一体化的电力体制,上网电价通常由政府制定,这里按照标杆电价给定,这里暂不考虑电价的随机变化。ph,j,t为水电站h机组j在t时段的发电功率;Δt为时间尺度,本文以月为调度时段,因此Δt为各月所包含的小时数;ph,j,t为水电站h机组j在t时段的发电功率。
约束条件包括:
水头约束
ph,j,t=Ah,jqh,j,tHh,t
其中Ah,j为水电站h机组j的出力系数,qh,j,t为水电站h机组j在t时段的发电流量;Hh,t为水电站h时段t的平均水头。上式为机组出力和发电流量、水头之间的非线性约束条件。
水量平衡约束
Figure BDA0001540688820000078
Figure BDA0001540688820000081
式中:vh,t表示水电站h在t时段的库容;qh,t为水电站h在t时段的发电流量;sh,t为水电站h在t时段的弃水流量;Rh,t为水电站h在t时段的来水的预测值;m为水电站编号;Uh为水电站h的上游水电站集合;qm,t和sm,t分别为上游电站的发电流量和弃水流量;Δt′为时段t所包含的秒数。
库容约束
Figure BDA0001540688820000082
式中:
Figure BDA0001540688820000083
为水电站h的库容最小、最大值。
弃水约束
Figure BDA0001540688820000084
式中:
Figure BDA0001540688820000085
分别为水电站h弃水流量的最小、最大值。
发电流量约束
Figure BDA0001540688820000086
式中:
Figure BDA0001540688820000087
分别为水电站h发电流量的最小、最大值。
机组出力约束
Figure BDA0001540688820000088
式中:
Figure BDA0001540688820000089
为水电站h中机组j的出力最小、最大值。
初始期末库容约束
vh,0=vh,ini
vh,2T=vh,term
式中:vh,ini为水电站h第一年年初库容,vh,term为水电站h第二年年末库容,水库的期末库容通常由更长期的优化调度所决定;
3)构建相应的安全校核的再调度max-min模型,基于对偶理论将其转化为max问题进行求解,max-min的数学优化问题,即:
Figure BDA00015406888200000810
式中:z表示博弈均衡值;λ1,h,t、λ2,h,t、λ3,h,t、λ4,h,t、λ5,h,t、λ6,h,t为相应约束条件中引入的松弛变量,具体如下:
再调度过程约束条件:
Figure BDA00015406888200000811
Figure BDA0001540688820000091
Figure BDA0001540688820000092
Figure BDA0001540688820000093
Figure BDA0001540688820000094
Figure BDA0001540688820000095
Figure BDA0001540688820000096
Figure BDA0001540688820000097
Figure BDA0001540688820000098
vh,0=vh,ini u10,h
vh,2T=vh,term u11,h,t
式中:
Figure BDA0001540688820000099
为不确定来水下的库容。
Figure BDA00015406888200000910
分别为t时段水库调节的上下限值。u1,h,t、u2,h,t、u3,h,t、u4,h,t、u5,h,t、u6,h,t、u7,h,t、u8,h,t、u9,h,t、u10,h、u11,h,t为约束对应式的对偶变量。
结合多波段来水不确定集构成了跨年梯级水电再调度博弈模型。
运用对偶原理,将再调度max-min模型转换为max模型,具体如下:
Figure BDA00015406888200000911
-u1,h,t+u2,h,t+u3,h,t+u4,h,t-u5,h,t≤0(t=1)
-u1,h,t+u2,h,t-u3,h,t-1+u3,h,t+u4,h,t-u5,h,t≤0(1<t<2T)
-u1,h,t+u2,h,t-u3,h,t-1+u3,h,t+u4,h,t-u5,h,t+u11,h,t≤0(t=2T)
u10,h-u3,h,t≤0(t=1)
Figure BDA0001540688820000101
Figure BDA0001540688820000102
u1,h,t≥-1
u2,h,t≥-1
-1≤u3,h,t≤1
u4,h,t≥-1
u5h,t≥-1
Figure BDA0001540688820000103
其中u1,h,t、u2,h,t、u4,h,t、u5,h,t、u6,h,t、u7,h,t、u8,h,t、u9,h,t≤0,u3,h,t、u10,h、u11,h,t无符号限制,
Figure BDA0001540688820000104
为0/1整数变量;
4)对预调度方案的鲁棒可行性进行了判定;针对不可行的情况,反馈鲁棒可行割集,对预调度决策进行校正,通过迭代获得能够适应不确定集中所有取值的安全经济调度方案,具体为:
利用混合整数二次线性规划法求解上述模型得到目标函数z。若z≤ε(ε为一个很小的正数)则预调度方案能够满足所有可能来水情景下的安全运行约束,否则,将下式返回到预调度模型中修正预调度方案,通过迭代直至满足收敛条件。
Figure BDA0001540688820000105
式中,
Figure BDA0001540688820000106
为预调度中传递到下层的变量,这里为预测场景下的水库库容
Figure BDA0001540688820000107
实施例1
下面具体结合一个包含2个水电站梯级水电系统进行详细分析。为了验证本发明的合理性和有效性与场景法进行对比,来水的预测信息如图1所示,两种模型对应优化结果比较情况如表1所示。从表中可以看出,场景数20和场景数100时的发电效益略大于鲁棒调度策略,这是因为鲁棒优化调度模型考虑了未来一切可能来水情况,保证最差来水情况下最大化水电站发电效益,因此其他场景下来水的发电效益不会小于这一最差情况下的发电效益。当场景数为500时,由于比较充分考虑来水的不确定性,所以场景法和鲁棒优化法得到的发电效益较为接近。但是随着考虑场景数的增大,场景法的求解计算时间明显增大,达到了分钟级别,而鲁棒优化调度求解时间仅为17秒。因此在准确表征来水不确定性的基础上,鲁棒优化调度方法的求解效率更高。
表1两种模型对应优化结果比较情况
Figure BDA0001540688820000111
为了验证鲁棒调度法和场景法处理不确定性的能力,将随机生成的10000组来水数据代入到相应模型中,观察这些场景的约束条件越限的情况,并进行统计,如表1所示。从表中可以看出,场景法均有一定程度的越限,而鲁棒优化调度的约束越限概率为0,这说明优化所得的调度方案可以应对第2年所有可能来水情况,具有鲁棒性。
实施例2:
为了进一步分析本发明的适用性,下面具体结合一个我国南方地区某大型梯级水电系统中进行仿真分析。该系统包含10个水电站。首先考虑来水的不确定波动波段为其预测值的10%,则传统的单波段不确定集为[1-10%,1+10%]。多波段的不确定集是以2%为间隔,将整个波段均分为5段,各波段对应比重系数分别为0.25,0.25,0.25,0.125,0.125。两种跨年梯级水电鲁棒调度模型的发电效益与不同预算约束值的曲线关系如图2所示。
从图2中可以看出,随着不确定预算值的增大,两种调度模型的发电效益都是呈现单调递减的趋势,为此通过改变不确定预算可以调整鲁棒调度模型的经济性。此外,相比较于单波段模型,所建立的多波段模型的发电效益对不确定预算的数值大小变化更加敏感,并且在不确定预算取值小于28之前多波段跨年梯级水电模型的发电效益均比单波段模型高,这表明多波段模型有效改善了跨年梯级水电鲁棒调度决策的保守性。
为了分析来水时序相关性约束对调度方案的影响,分别采用两种调度模型进行仿真分析。模型1考虑了来水时间相关性约束,模型2不考虑相关性约束,其他约束条件均相同。两种调度模型下水电站#6的跨年水库调度方案如图3所示。从图3中可以看出,两种调度模型所得到的水库库容调度方案变化趋势大致相同,因此跨年梯级水电多波段鲁棒调度模型具有实用性。然而考虑相关性约束的模型在枯水期水库库容明显大于不考虑相关性约束的模型,因此考虑相关性约束的模型更加细致表达了来水的季节性,得到了优化效果更好的调度方案。

Claims (9)

1.一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,包括以下步骤:
1)根据中长期梯级水电来水的不确定性及其时序相关性建立随机安全域,即考虑来水时序相关性的多波段不确定集合,具体包括以下步骤:
11)根据中长期梯级水电来水的随机特性,建立来水不确定集,并将来水不确定集划分为B段,得到来水多波段不确定集;
12)基于来水的时序分布特性以及流域相关性,建立来水时序相关性约束;
13)获取考虑来水时序相关性的多波段不确定集合Rf,即:
Figure FDA0003008038340000011
Figure FDA0003008038340000012
Figure FDA0003008038340000013
Figure FDA0003008038340000014
式中:h为水电站编号,t为时段编号,b为波段编号,b∈[1,B],B为波段总数,Rh,t为水电站h在时段t的预测均值,
Figure FDA0003008038340000015
分别为水电站h在时段t来水不确定第b波段子集正负偏离端点值,
Figure FDA0003008038340000016
均为二进制变量,
Figure FDA0003008038340000017
取值为1时为不确定变量落入第b正偏离波段,取值为0为未落入波段,
Figure FDA0003008038340000018
取值为1时为不确定变量落入第b负偏离波段,取值为0为未落入波段,
Figure FDA0003008038340000019
为水电站h在时段t的来水不确定值,πb为第b波段的比重系数,即来水落入第b波段的占比,其比重系数之和为1,ΔH1为不确定预算,表示所有
Figure FDA00030080383400000110
求和的上限,取值范围为[0,NT*NH],其中NT为调度时段数,NH为水电站数,
Figure FDA00030080383400000111
为向上取整,
Figure FDA00030080383400000112
为比重系数,ΔH2为不确定预算;
2)根据多波段不确定集合构建经济性和安全性协调的鲁棒两阶段模型;
3)求解鲁棒两阶段模型中第一阶段的跨年梯级水电预调度模型,得到水电站的预调度方案和预测场景下的水库库容;
4)将鲁棒两阶段模型中第二阶段的再调度模型转化为max-min的数学优化问题,并对预调度方案的鲁棒可行性进行判定,通过反馈修正协调优化最终获得具有鲁棒性的调度方案。
2.根据权利要求1所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的步骤2)中,鲁棒两阶段模型的第一阶段为跨年梯级水电预调度模型,其表达式为:
Figure FDA0003008038340000021
ph,j,t=Ah,jqh,j,tHh,t
其中,f为调度周期内梯级水电系统的总发电效益,NH为水电站总数,T为一年的时段总数,Jh为水电站h中的机组总数目,ζh,t为水电站h在t时段的市场电价,ph,j,t为水电站h机组j在t时段的发电功率,Δt为时间尺度,Ah,j为水电站h机组j的出力系数,qh,j,t为水电站h机组j在t时段的发电流量,Hh,t为水电站h时段t的平均水头。
3.根据权利要求2所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的跨年梯级水电预调度模型的约束条件包括:
1、水量平衡约束:
Figure FDA0003008038340000022
Figure FDA0003008038340000023
其中,vh,t为水电站h在t时段的库容,vh,t-1为水电站h在t-1时段的库容,qh,t为水电站h在t时段的发电流量,sh,t为水电站h在t时段的弃水流量,Rh,t为水电站h在t时段的来水的预测值,m为水电站编号,Uh为水电站h的上游水电站集合,qm,t和sm,t分别为上游电站的发电流量和弃水流量,Δt′为时段t所包含的秒数;
2、库容约束:
Figure FDA0003008038340000024
其中,
Figure FDA0003008038340000025
为水电站h的库容最小值和最大值;
3、弃水约束:
Figure FDA0003008038340000026
其中,
Figure FDA0003008038340000027
分别为水电站h弃水流量的最小值和最大值;
4、发电流量约束:
Figure FDA0003008038340000028
其中,
Figure FDA0003008038340000031
分别为水电站h发电流量的最小值和最大值;
5、机组出力约束:
Figure FDA0003008038340000032
其中,
Figure FDA0003008038340000033
为水电站h中机组j的出力最小值和最大值;
6、初始期末库容约束:
vh,0=vh,ini
vh,2T=vh,term
其中,vh,ini为水电站h第一年年初库容,vh,term为水电站h第二年年末库容,vh,0为水电站h期初库容,vh,2T为水电站h期末库容。
4.根据权利要求1所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的步骤3)中具体包括以下步骤:
设定来水为预测值,采用混合整数线性规划法求解跨年梯级水电预调度模型,获取水电站的预调度方案以及在预测场景下的水库库容
Figure FDA0003008038340000034
5.根据权利要求3所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的步骤4)中,第二阶段再调度模型的表达式为:
Figure FDA0003008038340000035
其中,z为博弈均衡值,λ1,h,t、λ2,h,t、λ3,h,t、λ4,h,t、λ5,h,t、λ6,h,t为相应约束条件中引入的松弛变量。
6.根据权利要求5所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的第二阶段再调度模型的约束条件包括:
Figure FDA0003008038340000036
Figure FDA0003008038340000037
Figure FDA0003008038340000038
Figure FDA0003008038340000039
Figure FDA00030080383400000310
Figure FDA00030080383400000311
Figure FDA00030080383400000312
Figure FDA00030080383400000313
Figure FDA00030080383400000314
vh,0=vh,ini
vh,2T=vh,term
其中,
Figure FDA0003008038340000041
为不确定来水下的库容,
Figure FDA0003008038340000042
分别为t时段水库调节的上下限值,
Figure FDA0003008038340000043
为预测场景下的水库库容,m表示水电站h的上游水电站编号,qm,t和sm,t分别为上游电站的发电流量和弃水流量。
7.根据权利要求6所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的步骤4)中,将max-min的数学优化问题转换为max问题,并采用混合整数二次线性规划法求解,若博弈均衡值z≤ε,ε为取值为0.001,则判定预调度方案能够满足所有可能来水情景下的安全运行约束,若博弈均衡值z>ε,则返回预调度模型中修正预调度方案,通过迭代直至满足收敛条件。
8.根据权利要求7所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,所述的返回预调度模型中修正预调度方案,按照下式进行:
Figure FDA0003008038340000044
其中,u为再调度模型的对偶变量,z为博弈均衡值,x为预调度变量,
Figure FDA0003008038340000045
为预调度中传递到下层的变量,即预测场景下的水库库容
Figure FDA0003008038340000046
9.根据权利要求1所述的一种基于随机安全域的梯级水电鲁棒优化调度方法,其特征在于,本梯级水电鲁棒优化调度方法的调度周期为2年。
CN201810012187.5A 2018-01-05 2018-01-05 一种基于随机安全域的梯级水电鲁棒优化调度方法 Active CN108388954B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810012187.5A CN108388954B (zh) 2018-01-05 2018-01-05 一种基于随机安全域的梯级水电鲁棒优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810012187.5A CN108388954B (zh) 2018-01-05 2018-01-05 一种基于随机安全域的梯级水电鲁棒优化调度方法

Publications (2)

Publication Number Publication Date
CN108388954A CN108388954A (zh) 2018-08-10
CN108388954B true CN108388954B (zh) 2021-07-20

Family

ID=63077003

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810012187.5A Active CN108388954B (zh) 2018-01-05 2018-01-05 一种基于随机安全域的梯级水电鲁棒优化调度方法

Country Status (1)

Country Link
CN (1) CN108388954B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109787821B (zh) * 2019-01-04 2020-06-19 华南理工大学 一种大规模移动客户流量消费智能预测方法
CN112465323B (zh) * 2020-11-19 2024-02-23 大连理工大学 一种耦合日电量分解及日前市场竞价的梯级水电站短期鲁棒调度方法
CN113312747B (zh) * 2021-04-19 2023-07-18 广西大学 考虑不确定性的水电动态库容经济调度方法
CN116388301B (zh) * 2023-06-06 2023-08-08 河海大学 梯级水光互补日前鲁棒调峰调度方法、系统、设备及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102832614A (zh) * 2012-08-03 2012-12-19 国电南瑞科技股份有限公司 不确定性环境下发电计划的鲁棒优化方法
CN103606967A (zh) * 2013-11-26 2014-02-26 华中科技大学 一种实现电力系统鲁棒运行的调度方法
EP3033826A1 (en) * 2013-08-12 2016-06-22 Zentrum Mikroelektronik Dresden AG Adaptive dual stage identification control method for a power stage of a power converter
CN107069716A (zh) * 2017-05-05 2017-08-18 山东大学 一种计及注入转移分布因子估计误差的鲁棒实时调度方法
CN107516909A (zh) * 2017-08-31 2017-12-26 华北电力大学(保定) 一种网架重构过程中可接入风电出力的优化方法和装置
CN107528348A (zh) * 2017-09-14 2017-12-29 华北电力大学 一种基于来水不确定性的梯级水电站负荷调整方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102832614A (zh) * 2012-08-03 2012-12-19 国电南瑞科技股份有限公司 不确定性环境下发电计划的鲁棒优化方法
EP3033826A1 (en) * 2013-08-12 2016-06-22 Zentrum Mikroelektronik Dresden AG Adaptive dual stage identification control method for a power stage of a power converter
CN103606967A (zh) * 2013-11-26 2014-02-26 华中科技大学 一种实现电力系统鲁棒运行的调度方法
CN107069716A (zh) * 2017-05-05 2017-08-18 山东大学 一种计及注入转移分布因子估计误差的鲁棒实时调度方法
CN107516909A (zh) * 2017-08-31 2017-12-26 华北电力大学(保定) 一种网架重构过程中可接入风电出力的优化方法和装置
CN107528348A (zh) * 2017-09-14 2017-12-29 华北电力大学 一种基于来水不确定性的梯级水电站负荷调整方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
一类电力系统鲁棒优化问题的博弈模型及应用实例;梅生伟;《中国电机工程学报》;20130705;第33卷(第19期);47-56 *
含多风电场的电力系统储能鲁棒优化配置方法;韩杏宁 等;《中国电机工程学报》;20141224;第35卷(第9期);2120-2127 *
基于分布鲁棒优化的光伏电站并网极限容量研究;向加佳 等;《电力建设》;20160601;第37卷(第6期);151-156 *
基于可行性检测的考虑风电和需求响应的机组组合鲁棒优化方法;仉梦林 等;《中国电机工程学报》;20171120;3184-3194 *
考虑风电时空相关性的仿射可调鲁棒机组组合;吴巍 等;《中国电机工程学报》;20170216;第37卷(第14期);4089-4097 *

Also Published As

Publication number Publication date
CN108388954A (zh) 2018-08-10

Similar Documents

Publication Publication Date Title
CN108388954B (zh) 一种基于随机安全域的梯级水电鲁棒优化调度方法
Apostolopoulou et al. Robust optimization for hydroelectric system operation under uncertainty
CN107240932B (zh) 一种水光互补系统中光伏电站容量优化方法
Liu et al. Deriving optimal refill rules for multi-purpose reservoir operation
WO2015154318A1 (zh) 梯级水电站群日发电计划编制方法
CN102930155A (zh) 获取电力需求的预警参数的方法及装置
CN107453408B (zh) 一种考虑不确定性的微电网能量优化调度方法
Zimmermann et al. Cross-border effects of capacity remuneration mechanisms: The Swiss case
Mi et al. Low-carbon generation expansion planning considering uncertainty of renewable energy at multi-time scales
CN111104750A (zh) 一种耦合集合预报信息的水库调度图绘制方法
CN114529100A (zh) 一种面向区域电网风光消纳的储能优化配置方法及系统
CN103714400A (zh) 一种任意输入参数的风能预测系统和方法
CN112803491B (zh) 一种耦合弃电风险的风光水多能互补短期优化调度方法
Henden et al. Economic analysis of large-scale pumped storage plants in Norway
Zhanga et al. Optimized operation of cascade reservoirs considering complementary characteristics between wind and photovoltaic based on variational auto-encoder
CN106529698B (zh) 一种日内梯级水电两阶段优化调度方法
Mirzapour-Kamanaj et al. Optimal look-ahead strategic bidding/offering of integrated renewable power plants and CAES with stochastic-robust approach
CN105095999A (zh) 一种基于改进轻鲁棒模型的分布式发电站规划方法
CN114254946A (zh) 新能源发电等值年费用比较方法、系统、设备及存储介质
Scarcelli et al. Aggregated inflows on stochastic dynamic programming for long term hydropower scheduling
Xu et al. Day-ahead battery scheduling in microgrid considering wind power uncertainty using ordinal optimization
Grønvik et al. Feasibility of linear decision rules for hydropower scheduling
Abdulla et al. Accounting for forecast uncertainty in the optimized operation of energy storage
Helseth et al. Approximating Ramping Constraints in Hydropower Scheduling
CN111709605A (zh) 一种基于多重反调节作用的水库电站调峰能力评估方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant