CN103093282B - 一种水电站群期末蓄能最大短期优化调度方法 - Google Patents
一种水电站群期末蓄能最大短期优化调度方法 Download PDFInfo
- Publication number
- CN103093282B CN103093282B CN201210579060.4A CN201210579060A CN103093282B CN 103093282 B CN103093282 B CN 103093282B CN 201210579060 A CN201210579060 A CN 201210579060A CN 103093282 B CN103093282 B CN 103093282B
- Authority
- CN
- China
- Prior art keywords
- power station
- constraint
- period
- sigma
- reservoir
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000009825 accumulation Methods 0.000 title claims abstract description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 29
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 238000003860 storage Methods 0.000 claims description 14
- 238000011144 upstream manufacturing Methods 0.000 claims description 14
- 230000009194 climbing Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 5
- 238000003892 spreading Methods 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 238000007689 inspection Methods 0.000 claims description 2
- 238000002372 labelling Methods 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000002688 persistence Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种水电站群期末蓄能最大短期优化调度方法,属于水电调度领域。首先将原模型中各时段总出力与总负荷需求的一致性约束,代替为各时段总出力与总负荷需求最大差值为零的约束;而后借助凝聚函数法将简化后的最大差值转换为便于处理的形式。在此基础上,采用拉格朗日松弛法建立对偶优化模型,并以次梯度法求解。本发明方法将存在规模为96的拉格朗日乘子向量的期末蓄能最大优化调度问题简化为只有一个拉格朗日乘子的问题,降低了问题的复杂性,采用次梯度法和逐次逼近搜索算法相结合的方法可有效求解,具有很强的实用性,应用广泛。<!--1-->
Description
技术领域
本发明属于水电调度领域,是一种涉及复杂约束条件下,水电站群期末蓄能最大短期优化调度方法。
背景技术
水电站短期优化调度一般指以15分钟为时段,未来几天内的优化调度问题,可以大致分为两类,第一类是给定各水库末水位,在调度期内按照调度目标调节水量时空分布,此类模型偏重于在满足各电站的中长期控制方式前提下,提高水能利用效率,是广义的以水定电方法,目标通常为发电量最大、发电效益最大等;第二类是给定各时段水电***负荷需求,在各电站间分配出力,是广义的以电定水方法,主要采用蓄能最大等模型,适于梯级和跨流域水电站群负荷分配问题。
蓄能最大模型求解受到负荷需求约束条件限制,较困难。目前流行的拉格朗日松弛算法求解效率较低,一方面由于短期调度每天包括96个时段,拉格朗日乘子向量规模较大,需进行很多次向量更新;另一方面,由于对偶问题包括出力持续性、爬坡速度等复杂约束,因而求解耗时较多。
发明内容
针对以拉格朗日松弛法求解水电站群短期蓄能最大模型时,乘子向量规模较大对计算效率的限制问题,提出:首先以各时段总出力与其目标间最大差值为零代替各时段总出力目标约束;而后借助凝聚函数法将简化后的最大差值进行转换。在此基础上,采用拉格朗日松弛法建立对偶模型,并采用次梯度法和逐次逼近搜索算法相结合的方法。
本发明的技术方案如下:
一种水电站群期末蓄能最大优化调度方法,以15分钟为时段,进行未来1天内的梯级和跨流域水电站群负荷分配;首先将原模型中各时段总出力与总负荷需求的一致性约束,代替为各时段总出力与总负荷需求最大差值为零的约束;而后借助凝聚函数法将简化后的最大差值进行转换,问题复杂度较原问题大大降低。最后模型采用拉格朗日松弛和逐次逼近搜索算法求解。
本发明首先将原模型中各时段总出力与总负荷需求的一致性约束,代替为各时段总出力与总负荷需求最大差值为零的约束,而后借助凝聚函数法将简化后的最大差值转换为便于处理的形式。该方法大大降低了问题的复杂性,采用次梯度法和逐次逼近搜索算法相结合的方法可有效求解,具有很强的实用性,应用广泛。
附图说明
图1是初始发电流量过程;
图2是扩展调整范围过程示意;
图3是流量变化速度修正示意;
图4是流量波动修正示意;
图5是出力升降时段修正示意;
图6是各电站出力过程;
图7是***总出力过程
具体实施方式
以下结合技术方案详细叙述本发明的具体实施例。
(一)建立短期蓄能最大优化调度模型
1.1目标函数
其中:M为跨流域水电站群参与计算电站的总数(1≤m≤M);t为时段标号,T为计算时段总数,控制期内1≤t≤T;ESm表示m号电站及其全部上游电站死水位以上水量在m号电站可产生的电量,采用式(2)计算;Elm表示m号电站上游电站在m号电站产生的滞时电量,采用式(3)计算。
其中:表示水库m第t时段初的库容;ηm表示m号电站平均耗水率;WT(m)表示m号电站水库及其全部上游电站水库调度期末死水位以上蓄水量;Um表示m号电站直接上游电站标号数组;Km为m号水电站水库的直接上游电站数目。
其中:,WL(m)表示m号电站全部上游电站滞时水量计算函数;,表示有直接上下游关系的上游Um[k]号水电站水库和下游m号水电站水库间的最大最小滞时时段数;Δt表示第t时段的小时数;表示Um[k]号电站在n时段出库流量为时,在下游m号电站T-1时段以后产生的入库流量,按式(4)计算:
为Um[k]号电站水库出库流量为时的滞时时段数。
1.2约束条件
除常规水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束等条件外,该模型尚需满足如下约束:
●水电站出力爬坡限制
(MW)表示m号水电站相邻时段最大出力升降限制。
●水电站出力波动限制
tvm是m号水电站出力升降最小间隔时段数,即在一轮出力升降过程最高最低点需持续最少tvm个时段。
●最小出力升降时段数限制
m号水电站出力上升开始至下降开始,或由下降开始至上升开始的时间间隔不少于tpm个时段。
●水电总负荷约束
其中:Nt表示第t时段的水电负荷。
(二)短期蓄能最大优化调度方法
2.1模型简化
将约束,t=1,2,…,T,转换为:。
而后采用凝聚函数法将转换为,。P为参数,取为0.01。
2.2建立对偶问题
针对水电***总负荷约束采用拉格朗日松弛法建立对偶优化问题。在不采用2.1模型时,建立以式(8)为目标函数,以水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束、水电站出力爬坡限制、最小出力升降时段数限制为约束条件的优化问题。
λt为t时段拉格朗日乘子向量分量, ,B={λ1,λ2,…,λT}。
而在采用简化方式时,对偶优化问题约束条件相同,目标函数为:
γ为拉格朗日乘子。
2.3优化求解
上述两种模型,均可采用次梯度法通过交替进行对偶优化问题求解和乘子向量的更新实现原问题求解,式(8)含有96个乘子,而式(9)仅包含一个乘子,显然后者复杂性要小得多。计算流程如下:
(1)初设i=1和拉格朗日乘子γi;
(2)采用逐次逼近的搜索算法求解以式(9)为目标函数,以水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束、水电站出力爬坡限制、最小出力升降时段数限制为约束条件的优化问题,此时;
(3)采用次梯度法更新乘子,即γi+1=γi+β*F2,β为乘子更新系数,设i=i+1;
(4)若|γi+1-γi|<ε,ε为收敛精度,转步骤(5),否则返回步骤(2)。
(5)在优化结果的基础上,以为目标消除对偶间隙,仍采用逐次逼近的算法求解。
流程中采用的逐次逼近的搜索算法需要使得搜索模式满足出力爬坡、波动和升降时段约束,本发明采用如下方法:
(1)在t时段(如图1)发电流量由变为后,依次向前后方向扩展,以满足出力波动约束(如图2)。计原具体方法是依次向前调整t-1,t-2,t-tvm/2时段的出力,直到。向后搜索的方法类似。
(2)在扩展范围边缘修正出力变化幅度,以满足出力爬坡约束。如图3,设扩展范围的范围为[t1,t2],则首先搜索最大的满足,(t0<t1)的时段t0,由时段依次按照爬坡出力值调整t1-1,t1-2,…t0+1的发电流量;而后搜索最小的满足,(t3>t2)的时段t3,而后由时段依次按照爬坡出力值调整t2+1,t2+2,…t3-1的发电流量。
(3)在前两步发电流量变化范围的两侧,检验出力升降时段约束是否满足,必要时进行修正。如图4,(2)步后发电流量的变化范围为[t0,t3],则由t0时段向前搜索最大满足的t4,若|t0-t4|<tvm,则搜索t5,t6使得|t5-t6|≥tvm,t5≥t0,t6≤t4,,否则t6=t4。由t3点向后修整的方法类似。
(4)如图5,第(3)步后发电流量的变化范围为[t6,t7],t7为由t3时段向后修正范围的最大时段号。则由t5时段向前搜索最大的满足,t8≤t5的t8,若|t5-t8|<tpm,则搜索距离最大的t9,使得,t9<t8,修正t9和t5间的流量过程满足出力升降时段要求。由t7点向后修整的方法类似。
3应用实例
采用本发明的方法求解乌江流域某日短期优化调度问题,图6为各电站优化出力过程,图7为给定和计算总负荷过程。
本发明的具体实施方式在各方面应被视为例示性而非限制性实施例,所有的改变只要合乎本发明权利要求书所定义的范围或为其技术实施方式等效者,均应包含在本发明的保护范畴中。
Claims (1)
1.一种水电站群期末蓄能最大优化调度方法,其特征包括以下步骤,
(一)建立短期蓄能最大优化调度模型
1.1目标函数
其中:M为跨流域水电站群参与计算电站的总数,1≤m≤M;t为时段标号,T为计算时段总数,控制期内1≤t≤T;ESm表示m号电站及其全部上游电站死水位以上水量在m号电站可产生的电量,采用式(2)计算;Elm表示m号电站上游电站在m号电站产生的滞时电量,采用式(3)计算;
其中:表示水库m第t时段初的库容;ηm表示m号电站平均耗水率;WT(m)表示m号电站水库及其全部上游电站水库调度期末死水位以上蓄水量;Um表示m号电站直接上游电站标号数组;Km为m号水电站水库的直接上游电站数目;
Elm=WL(m)/ηm(3)
其中: WL(m)表示m号电站全部上游电站滞时水量计算函数;表示有直接上下游关系的上游Um[k]号水电站水库和下游m号水电站水库间的最大最小滞时时段数;Δt表示第t时段的小时数;表示Um[k]号电站在n时段出库流量为时,在下游m号电站T-1时段以后产生的入库流量,按式(4)计算:
为Um[k]号电站水库出库流量为时的滞时时段数;
1.2约束条件
除常规水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束条件外,该模型尚需满足如下约束:
水电站出力爬坡限制:
(MW)表示m号水电站相邻时段最大出力升降限制;
水电站出力波动限制:
tvm是m号水电站出力升降最小间隔时段数,即在一轮出力升降过程最高最低点需持续最少tvm个时段;
最小出力升降时段数限制:m号水电站出力上升开始至下降开始,或由下降开始至上升开始的时间间隔不少于tpm个时段;
水电总负荷约束:
其中:Nt表示第t时段的水电负荷;
(二)短期蓄能最大优化调度方法
2.1模型简化
将约束 转换为:
而后采用凝聚函数法将 转换为 P为参数,取为0.01;
2.2建立对偶问题
针对水电***总负荷约束采用拉格朗日松弛法建立对偶优化问题;在不采用2.1模型时,建立以式(8)为目标函数,以水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束、水电站出力爬坡限制、最小出力升降时段数限制为约束条件的优化问题;
λt为t时段拉格朗日乘子向量分量, B={λ1,λ2,...,λT};
而在采用简化方式时,对偶优化问题约束条件相同,目标函数为:
γ为拉格朗日乘子;
2.3优化求解
计算流程如下:
(1)初设i=1和拉格朗日乘子γi;
(2)采用逐次逼近的搜索算法求解以式(9)为目标函数,以水量平衡、发电流量约束、电站出力约束、库水位约束、出库流量约束、水电站出力爬坡限制、最小出力升降时段数限制为约束条件的优化问题,此时
(3)采用次梯度法更新乘子,即γi+1=γi+β*F2,β为乘子更新系数,设i=i+1;
(4)若|γi+1-γi|<ε,ε为收敛精度,转步骤(5),否则返回步骤(2);
(5)在优化结果的基础上,以为目标消除对偶间隙,仍采用逐次逼近的算法求解;
流程中采用的逐次逼近的搜索算法需要使得搜索模式满足出力爬坡、波动和升降时段约束,具体方法如下:
1)在t时段发电流量由变为后,依次向前后方向扩展,以满足出力波动约束;计算的具体方法是依次向前调整t-1,t-2,t-tvm/2时段的出力,直到向后搜索的方法类似;
2)在扩展范围边缘修正出力变化幅度,以满足出力爬坡约束,设扩展范围的范围为[t1,t2],则首先搜索最大的满足(t0<t1)的时段t0,由时段依次按照爬坡出力值调整t1-1,t1-2,...t0+1的发电流量;而后搜索最小的满足(t3>t2)的时段t3,而后由时段依次按照爬坡出力值调整t2+1,t2+2,...t3-1的发电流量;
3)在前两步发电流量变化范围的两侧,检验出力升降时段约束是否满足,必要时进行修正;2)步后发电流量的变化范围为[t0,t3],则由t0时段向前搜索最大满足的t4,若|t0-t4|<tvm,则搜索t5,t6使得|t5-t6|≥tvm,t5≥t0,t6≤t4,否则t6=t4;由t3点向后修整的方法类似;
4)第3)步后发电流量的变化范围为[t6,t7],t7为由t3时段向后修正范围的最大时段号;则由t5时段向前搜索最大的满足t8≤t5的t8,若|t5-t8|<tpm,则搜索距离最大的t9,使得t9<t8,修正t9和t5间的流量过程满足出力升降时段要求;由t7点向后修整的方法类似。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210579060.4A CN103093282B (zh) | 2012-12-27 | 2012-12-27 | 一种水电站群期末蓄能最大短期优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210579060.4A CN103093282B (zh) | 2012-12-27 | 2012-12-27 | 一种水电站群期末蓄能最大短期优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103093282A CN103093282A (zh) | 2013-05-08 |
CN103093282B true CN103093282B (zh) | 2016-06-22 |
Family
ID=48205828
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210579060.4A Active CN103093282B (zh) | 2012-12-27 | 2012-12-27 | 一种水电站群期末蓄能最大短期优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103093282B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104063808B (zh) * | 2014-06-27 | 2015-06-24 | 大连理工大学 | 一种跨省送电梯级水电站群调峰调度两阶段搜索方法 |
CN104182804B (zh) * | 2014-08-21 | 2017-03-29 | 大连理工大学 | 一种预测出力不确定小水电与大中型水电站协调的日前发电方法 |
CN105225017B (zh) * | 2015-10-30 | 2019-02-26 | 南京南瑞集团公司 | 一种多Agent的水电站群短期优化调度方法 |
CN106099960B (zh) * | 2016-07-18 | 2018-12-28 | 华中科技大学 | 一种小水电机群组成分布式储能***的方法 |
CN106327022B (zh) * | 2016-08-31 | 2019-11-12 | 华北电力大学 | 一种稳定梯级水电站发电流量波动的方法 |
CN109902401A (zh) * | 2019-03-05 | 2019-06-18 | 东北电力大学 | 一种基于拉格朗日松弛法、正割法和差分演化法的发电机组组合求解方法 |
CN110350537B (zh) * | 2019-07-10 | 2021-06-01 | 佳源科技有限公司 | 基于可再生能源整合的微电网优化调度方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430358A (zh) * | 2008-12-10 | 2009-05-13 | 河南国网宝泉抽水蓄能有限公司 | 抽水蓄能电站首机首次水泵工况整组启动试验方法 |
-
2012
- 2012-12-27 CN CN201210579060.4A patent/CN103093282B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101430358A (zh) * | 2008-12-10 | 2009-05-13 | 河南国网宝泉抽水蓄能有限公司 | 抽水蓄能电站首机首次水泵工况整组启动试验方法 |
Non-Patent Citations (2)
Title |
---|
梯级水电站短期周优化调度规律探讨;李亮,周云,黄强;《水力发电学报》;20090825(第04期);38-42 * |
贵州电网梯级水电站群联合调度运行方式研究;苏华英,孙斌,林成,程春田;《水电能源科学》;20100625(第06期);115-117 * |
Also Published As
Publication number | Publication date |
---|---|
CN103093282A (zh) | 2013-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103093282B (zh) | 一种水电站群期末蓄能最大短期优化调度方法 | |
CN102855591B (zh) | 梯级水库群短期联合发电优化调度方法及*** | |
CN104167730B (zh) | 一种复杂约束限制下的梯级水电站群实时优化调度方法 | |
CN101705671B (zh) | 黄河上游梯级水电站运行设计与优化调度方法及设备 | |
CN103971198B (zh) | 梯级水电站厂间负荷实时分配***及其实现方法 | |
CN102182159B (zh) | 一种梯级水库汛限水位联合运用调度方法 | |
CN103745023A (zh) | 水电站发电出力方案制作与最优负荷分配耦合建模方法 | |
CN107681695B (zh) | 一种储能辅助火电机组调频的容量配置方法 | |
CN105006843A (zh) | 一种应对风电不确定性的多时间尺度柔性负荷调度方法 | |
CN109447405A (zh) | 一种承担调峰任务的一库多级式梯级库群短期计划制定方法 | |
CN105574605A (zh) | 一种梯级泵站输水***日优化调控方法 | |
CN104239980B (zh) | 一种梯级水电站群短期优化调度方法 | |
CN104077635A (zh) | 一种基于光伏发电***的电动汽车充电站充电优化方法 | |
CN102080366A (zh) | 一种梯级水库联合调度图制作方法 | |
CN103020742A (zh) | 具有多限制运行区的梯级水电站群短期优化调度方法 | |
CN103065206A (zh) | 复杂时段耦合型约束下的水电站群变尺度优化调度方法 | |
CN103023073A (zh) | 一种水电站群电站与机组混合优化调度方法 | |
CN105184426B (zh) | 一种基于随机连续寻优策略的梯级电站调峰方法 | |
CN104701888A (zh) | 一种电网水、火电工况自匹配模式调整实时负荷的方法 | |
CN104504455B (zh) | 一种梯级蓄能控制下水电站群长期优化调度方法 | |
CN109726907A (zh) | 一种水电站群中长期出力计算方法及*** | |
CN112186802A (zh) | 一种动态经济调度的多时间尺度滚动调度方法及其*** | |
CN110942212A (zh) | 一种基于梯级水库运行系数的梯级水库优化运行方法 | |
CN104674748B (zh) | 一种梯级水电站群弃水概率量化方法 | |
CN116544978A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |