CN115640894A - 一种基于双时间尺度的综合能源系统优化调度方法 - Google Patents

一种基于双时间尺度的综合能源系统优化调度方法 Download PDF

Info

Publication number
CN115640894A
CN115640894A CN202211294624.XA CN202211294624A CN115640894A CN 115640894 A CN115640894 A CN 115640894A CN 202211294624 A CN202211294624 A CN 202211294624A CN 115640894 A CN115640894 A CN 115640894A
Authority
CN
China
Prior art keywords
power
hydrogen
time scale
time period
photovoltaic
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
Application number
CN202211294624.XA
Other languages
English (en)
Other versions
CN115640894B (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.)
Shandong University
Original Assignee
Shandong University
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 Shandong University filed Critical Shandong University
Priority to CN202211294624.XA priority Critical patent/CN115640894B/zh
Publication of CN115640894A publication Critical patent/CN115640894A/zh
Application granted granted Critical
Publication of CN115640894B publication Critical patent/CN115640894B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开一种基于双时间尺度的综合能源系统优化调度方法,该方法的具体步骤如下:基于综合能源系统的光伏、风电历史运行数据,建立光伏、风电长时间尺度下的功率波动区间模型;在长时间尺度下,考虑多类型电‑热‑气调节设备的调节能力,以综合用能成本最低为目标,建立长时间尺度下的鲁棒优化调度模型,求解得到长时间尺度下各调控设备的调度方案;在短时间尺度下,考虑光伏、风电功率与其在长时间尺度下的预测偏差,利用氢储能的快速响应调节能力,建立短时间尺度修正调度模型,求解得到氢储能的修正调度方案,平抑光伏、风电实时功率波动。本发明提供的调度方法,能够提高综合能源系统的用能经济性。

Description

一种基于双时间尺度的综合能源系统优化调度方法
技术领域
本发明涉及电力低压配电网技术领域,尤其涉及一种基于电力电子技术和通信技术的双时间尺度综合能源系统优化调度方法。
背景技术
近年来,我国风光装机容量增长迅猛,在电网中发电比例逐年提升,而综合能源系统作为一种能够实现多能互补、不同种类能源协同运行的有机载体,它能够有效促进可再生能源就地消纳,但因可再生能源固有的较大随机性与波动性等特征,其大规模的接入带来的不确定性不断给电网的可靠运行带来挑战,并且影响了综合能源系统的经济效益。现阶段提出了相关的综合能源协调调度方法,但缺乏对多类型调控设备在不同时间尺度上的协同考虑。
基于此,本专利根据综合能源系统的电-热-气多类型调节设备特性,通过长、短双时间尺度的协同调度,实现对光伏、风电等间歇性可再生能源实时功率波动平抑,减少综合能源系统实际用能与计划用能的偏差,提高综合能源系统的用能经济性。
发明内容
本发明的目的在于通过双时间尺度弥补可再生能源发电不确定性,有助于提高综合能源系统的供电可靠性和降低系统运行成本,具体技术方案如下:
一种基于双时间尺度的综合能源系统优化调度方法,共包括以下步骤:
步骤1,基于综合能源系统的光伏、风电历史运行数据,建立光伏、风电长时间尺度下的功率波动区间模型;
步骤2,在长时间尺度下,考虑多类型电-热-气调节设备的调节能力,以综合用能成本最低为目标,建立长时间尺度下的鲁棒优化调度模型,求解得到长时间尺度下各调控设备的调度方案;
步骤3,在短时间尺度下,考虑光伏、风电功率与其在长时间尺度下的预测偏差,利用氢储能的快速响应调节能力,建立短时间尺度修正调度模型,求解得到氢储能的修正调度方案。
进一步地,所述步骤1中所述的可再生能源功率波动区间,具体表示如下:
Figure BDA0003902020000000021
Figure BDA0003902020000000022
Figure BDA0003902020000000023
Figure BDA0003902020000000024
Figure BDA0003902020000000025
Figure BDA0003902020000000026
式中,
Figure BDA0003902020000000027
分别对应长时间尺度下光伏、风电有功功率实际值;
Figure BDA0003902020000000028
分别对应长时间尺度下光伏、风电有功功率预测值;
Figure BDA0003902020000000029
为因高估光伏、风电出力后需额外购买电能的惩罚成本,
Figure BDA00039020200000000210
为高估后购买电能惩罚价格,
Figure BDA00039020200000000211
为高估光伏、风电出力;
Figure BDA00039020200000000212
为因低估光伏、风电出力出售额外电能收益,
Figure BDA00039020200000000213
为低估后售电惩罚价格,
Figure BDA00039020200000000214
为低估光伏、风电出力;
Figure BDA00039020200000000215
分别为光伏、风电有功功率预测偏差;ΓPV、ΓWIN是鲁棒系数,为了避免决策过于保守。
进一步地,所述步骤2中所述的考虑多类型电-热-气调节设备的调节能力,以综合用能成本最低为目标,建立长时间尺度下的鲁棒优化调度模型,具体为:
多类型电-热-气调节设备包括电动汽车、燃气轮机、氢储能、燃气锅炉;
长时间尺度下的鲁棒优化调度模型的目标函数与约束条件具体为:
(1)目标函数:
Figure BDA0003902020000000031
Figure BDA0003902020000000032
Figure BDA0003902020000000033
Figure BDA0003902020000000034
Figure BDA0003902020000000035
Figure BDA0003902020000000036
式中:
Figure BDA0003902020000000037
分别为t时间段内购买天然气成本、向电动汽车用户提供补偿成本、燃气轮机环境成本、综合能源系统向上级电网售电收益以及氢储能运维成本;
Figure BDA0003902020000000038
为天然气价格;
Figure BDA0003902020000000039
为t时间段天然气负荷;
Figure BDA00039020200000000310
为燃气轮机t时间段有功出力;ηgt为燃气轮机气电效率系数;
Figure BDA00039020200000000311
为t时间段锅炉有功出力;ηgb为锅炉气电效率系数;eev为综合能源体向提供放电服务的电动汽车用户提供的补偿成本;
Figure BDA00039020200000000312
为电动汽车蓄电池t时间段内的放电功率;N 为电动汽车总数量;Ql为燃气轮机第l项污染气体的排放强度;Yl、Bl分别为第l项污染气体的罚款成本系数和环境成本系数;
Figure BDA00039020200000000313
为t时间段综合能源体向上级电网售、购电功率大小;
Figure BDA00039020200000000314
分别为综合能源体向上级电网售、购电电价;
Figure BDA00039020200000000315
分别对应t时间段氢储能系统电制氢和氢燃料电池发电功率;ηOMH2为氢储能系统运维成本系数;
(2)约束条件
1)电动汽车调节约束:
Figure BDA0003902020000000041
Figure BDA0003902020000000042
Figure BDA0003902020000000043
Figure BDA0003902020000000044
Figure BDA0003902020000000045
式中:SOCev,t,i为t时间段第i组电动汽车蓄电池荷电状态;
Figure BDA0003902020000000046
分别对应第i组电动汽车蓄电池荷电状态下限值和上限值;
Figure BDA0003902020000000047
分别对应 t时间段第i组电动汽车充放电功率;
Figure BDA0003902020000000048
分别对应电动汽车蓄电池的充放电效率系数;
Figure BDA0003902020000000049
分别对应第i组电动汽车蓄电池最大充放电功率;0-1变量
Figure BDA00039020200000000410
分别对应第i组电动汽车t时间段充放电标志位,充电则取1,否则取0,两者不能同时为1;
2)燃气轮机调节约束:
Figure BDA00039020200000000411
Figure BDA00039020200000000412
式中:
Figure BDA00039020200000000413
分别对应燃气轮机最小和最大输出功率;
Figure BDA00039020200000000414
分别对应燃气轮机向下、向上爬坡功率;
4)氢储能调节约束:
Figure BDA00039020200000000415
式中:
Figure BDA00039020200000000416
为t时间段电制氢功率;
Figure BDA00039020200000000417
为t时间段氢燃料电池发电功率;
Figure BDA00039020200000000418
Figure BDA00039020200000000419
分别对应电制氢和氢燃料电池发电的0/1变量,目的为了避免同时制氢和发电;
5)燃气锅炉调节约束:
Figure BDA0003902020000000051
Figure BDA0003902020000000052
式中:
Figure BDA0003902020000000053
为锅炉最小和最大输出功率;
Figure BDA0003902020000000054
为燃气锅炉向下、向上爬坡功率;
6)综合能源系统电能供需平衡约束:
Figure BDA0003902020000000055
式中:
Figure BDA0003902020000000056
分别对应t时间段风力发电和光伏发电功率大小;
Figure BDA0003902020000000057
为t时间段电负荷功率;
7)综合能源系统热能供需平衡约束:
Figure BDA0003902020000000058
式中:
Figure BDA0003902020000000059
为t时间段热负荷功率大小;
Figure BDA00039020200000000510
为t时间段弃热功率大小;
8)综合能源系统天然气供需平衡约束:
Figure BDA00039020200000000511
式中:
Figure BDA00039020200000000512
为t时间段综合能源主体天然气总消耗量;
Figure BDA00039020200000000513
分别对应t时间段燃气轮机、燃气锅炉耗气量;
Figure BDA00039020200000000514
为t时间段天然气负荷;
9)综合能源系统氢气供需平衡约束:
Figure BDA00039020200000000515
式中:
Figure BDA00039020200000000516
为t时间段储氢罐容量;
Figure BDA00039020200000000517
为t时间段氢气负荷;
Figure BDA00039020200000000518
为每千瓦时电能通过电解槽电解水产生氢气质量;η为每千克氢气通过燃料电池燃烧所产生的电能;θel、θfc分别对应电解槽和燃料电池的对应转换效率;
Figure BDA00039020200000000519
代表氢气的摩尔质量。
进一步地,所述步骤2中所述的求解得到长时间尺度下各调控设备的调度方案,采用C&CG列约束生成法进行求解,具体步骤如下:
1)设置最大迭代次数kmax与收敛精度ε;设置鲁棒优化模型目标函数值上界UB与下界LB;设置初始迭代次数k;输入综合能源系统电-热-气和氢气负荷数据、各设备机组参数和风、光和负荷的预测数据;
2)求解主问题:根据输入的风、光、负荷预测值求解长时间尺度主问题最优解,寻求综合能源体长时间尺度下最小成本调度策略,并作为鲁棒优化模型下界;
3)求解子问题:将主问题优化结果带入子问题,寻求当前调度方案下风、光发电最恶劣申报场景,并再次求解该场景下的调度方案,得到方案后反馈给主问题,将子问题最优目标函数值和主问题最优目标函数值之和作为鲁棒优化模型新的上界;
4)判断收敛:判断鲁棒优化调度模型子问题与主问题的优化目标函数之差是否小于收敛精度,若满足收敛条件,则停止迭代,输出综合能源系统各设备长时间尺度的调度方案;否则,将子问题求得的恶劣场景代入主问题中,并令k=k+1,返回步骤2)继续迭代求解,直到满足收敛条件。
进一步地,所述步骤3中所述的建立短时间尺度修正调度模型,其目标函数为:
Figure BDA0003902020000000061
式中:
Figure BDA0003902020000000062
为短时间尺度下因光伏、风电预测偏差带来的购售电量申报偏差惩罚成本;
Figure BDA0003902020000000063
为氢储能运维成本增量;
短时间尺度修正调度模型的约束条件同式(20)、(23)~(26)。
本发明的有益效果在于:充分利用了综合能源系统中电-热-气不同类型设备在不同时间尺度下的调节特性,通过双时间尺度下的协同调度,有效平抑光伏、风电等间歇性可再生能源的实时功率波动,减少综合能源系统实际用能与计划用能的偏差,提高综合能源系统的用能经济性。
附图说明
图1是基于双时间尺度的综合能源系统优化调度方法流程图。
图2是长时间尺度下具有鲁棒特性的综合能源系统电平衡调度方案。
图3是长时间尺度下综合能源系统热平衡调度方案。
图4是长时间尺度下综合能源系统天然气平衡调度方案。
图5是短时间尺度下氢储能弥补可再生能源功率波动调度方案。
图6即表1是对比仅考虑长时间尺度调度方案(场景1)和同时考虑双时间尺度协同调度(场景2)的综合能源系统运行成本和各设备运行成本。
具体实施方式
以下结合附图和实例对本发明的具体实施做进一步说明。
如图1为本发明实施例提供的一种基于双时间尺度的综合能源系统优化调度方法,包括以下步骤:
步骤1,基于综合能源系统的光伏、风电历史运行数据,建立光伏、风电长时间尺度下的功率波动区间模型,具体表示如下:
Figure BDA0003902020000000071
Figure BDA0003902020000000072
Figure BDA0003902020000000073
Figure BDA0003902020000000074
Figure BDA0003902020000000081
Figure BDA0003902020000000082
式中,
Figure BDA0003902020000000083
分别对应长时间尺度下光伏、风电有功功率实际值;
Figure BDA0003902020000000084
分别对应长时间尺度下光伏、风电有功功率预测值;
Figure BDA0003902020000000085
为因高估光伏、风电出力后需额外购买电能的惩罚成本,
Figure BDA0003902020000000086
为高估后购买电能惩罚价格,取值为购电电价的1.5倍;
Figure BDA0003902020000000087
为高估光伏、风电出力;
Figure BDA00039020200000000812
为因低估光伏、风电出力出售额外电能收益,
Figure BDA0003902020000000088
为低估后售电惩罚价格,取值为售电电价的0.5 倍;
Figure BDA00039020200000000813
为低估光伏、风电出力;
Figure BDA0003902020000000089
分别为光伏、风电有功功率预测偏差,取值为±15%;ΓPV、ΓWIN是鲁棒系数,为了避免决策过于保守。
步骤1建立功率波动区间模型的作用是对长时间尺度下光伏、风电有功功率进行预估,以用作步骤2中综合能源系统电能供需平衡约束;及步骤3中对于短时间尺度下同等光伏、风电有功功率预测偏差的情况下,对进行短时间尺度修正调度。
步骤2,在长时间尺度下,考虑多类型电-热-气调节设备的调节能力,以综合用能成本最低为目标,建立长时间尺度下的鲁棒优化调度模型,求解得到长时间尺度下各调控设备的调度方案,具体地,
多类型电-热-气调节设备包括电动汽车、燃气轮机、氢储能、燃气锅炉;
长时间尺度下的鲁棒优化调度模型的目标函数与约束条件具体为:
(1)目标函数:
Figure BDA00039020200000000810
Figure BDA00039020200000000811
Figure BDA0003902020000000091
Figure BDA0003902020000000092
Figure BDA0003902020000000093
Figure BDA0003902020000000094
式中:
Figure BDA0003902020000000095
分别为t时间段内购买天然气成本、向电动汽车用户提供补偿成本、燃气轮机环境成本、综合能源系统向上级电网售电收益以及氢储能运维成本;
Figure BDA0003902020000000096
为天然气价格;
Figure BDA0003902020000000097
为t时间段天然气负荷;
Figure BDA0003902020000000098
为燃气轮机t时间段有功出力;ηgt为燃气轮机气电效率系数;
Figure BDA0003902020000000099
为t时间段锅炉有功出力;ηgb为锅炉气电效率系数;eev为综合能源体向提供放电服务的电动汽车用户提供的补偿成本;Pdisc,t为电动汽车蓄电池t时间段内的放电功率;N 为电动汽车总数量;Ql为燃气轮机第l项污染气体的排放强度;Yl、Bl分别为第l项污染气体的罚款成本系数和环境成本系数;
Figure BDA00039020200000000910
为t时间段综合能源体向上级电网售、购电功率大小;
Figure BDA00039020200000000911
分别为综合能源体向上级电网售、购电电价;
Figure BDA00039020200000000912
分别对应t时间段氢储能系统电制氢和氢燃料电池发电功率;ηOMH2为氢储能系统运维成本系数;
(2)约束条件
1)电动汽车调节约束:
Figure BDA00039020200000000913
Figure BDA00039020200000000914
Figure BDA00039020200000000915
Figure BDA00039020200000000916
Figure BDA00039020200000000917
式中:SOCev,t,i为t时间段第i组电动汽车蓄电池荷电状态;
Figure BDA0003902020000000101
分别对应第i组电动汽车蓄电池荷电状态下限值和上限值;
Figure BDA0003902020000000102
分别对应 t时间段第i组电动汽车充放电功率;
Figure BDA0003902020000000103
分别对应电动汽车蓄电池的充放电效率系数;
Figure BDA0003902020000000104
分别对应第i组电动汽车蓄电池最大充放电功率;0-1变量
Figure BDA0003902020000000105
分别对应第i组电动汽车t时间段充放电标志位,充电则取1,否则取0,两者不能同时为1;
2)燃气轮机调节约束:
Figure BDA0003902020000000106
Figure BDA0003902020000000107
式中:
Figure BDA0003902020000000108
分别对应燃气轮机最小和最大输出功率;
Figure BDA0003902020000000109
分别对应燃气轮机向下、向上爬坡功率;
4)氢储能调节约束:
Figure BDA00039020200000001010
式中:
Figure BDA00039020200000001011
为t时间段电制氢功率;
Figure BDA00039020200000001012
为t时间段氢燃料电池发电功率;
Figure BDA00039020200000001013
Figure BDA00039020200000001014
分别对应电制氢和氢燃料电池发电的0/1变量,目的为了避免同时制氢和发电;
5)燃气锅炉调节约束:
Figure BDA00039020200000001015
Figure BDA00039020200000001016
式中:
Figure BDA00039020200000001017
为锅炉最小和最大输出功率;
Figure BDA00039020200000001018
为燃气锅炉向下、向上爬坡功率;
6)综合能源系统电能供需平衡约束:
Figure BDA0003902020000000111
式中:
Figure BDA0003902020000000112
分别对应t时间段风力发电和光伏发电功率大小;
Figure BDA0003902020000000113
分别对应t时间段电储能系统充放电功率;
Figure BDA0003902020000000114
为t时间段电负荷功率;
7)综合能源系统热能供需平衡约束:
Figure BDA0003902020000000115
式中:
Figure BDA0003902020000000116
分别对应t时间段热储能储热、放热功率;
Figure BDA0003902020000000117
为t时间段热负荷功率大小;
Figure BDA0003902020000000118
为t时间段弃热功率大小;
8)综合能源系统天然气供需平衡约束:
Figure BDA0003902020000000119
式中:
Figure BDA00039020200000001110
为t时间段综合能源主体天然气总消耗量;
Figure BDA00039020200000001111
分别对应t时间段燃气轮机、燃气锅炉耗气量;
Figure BDA00039020200000001112
为t时间段天然气负荷;
9)综合能源系统氢气供需平衡约束:
Figure BDA00039020200000001113
式中:
Figure BDA00039020200000001114
为t时间段储氢罐容量;
Figure BDA00039020200000001115
为t时间段氢气负荷;
Figure BDA00039020200000001116
为每千瓦时电能通过电解槽电解水产生氢气质量;η为每千克氢气通过燃料电池燃烧所产生的电能;θel、θfc分别对应电解槽和燃料电池的对应转换效率;
Figure BDA00039020200000001117
代表氢气的摩尔质量。
进一步地,所述步骤2中所述的求解得到长时间尺度下各调控设备的调度方案,采用C&CG列约束生成法进行求解,具体步骤如下:
1)设置最大迭代次数kmax=50,收敛精度ε=0.01;设置鲁棒优化模型目标函数值上界UB=+∞与下界LB=-∞;设置初始迭代次数k=1;输入综合能源系统电-热-气和氢气负荷数据、各设备机组参数和风、光和负荷的预测数据;风、光的预测数据由步骤1可得,负荷的预测数据通过步骤1光伏、风电有功功率预测值,将光伏、风电有功功率预测值进行结合得到负荷的预测值;
2)求解主问题:根据输入的风、光、负荷预测值求解长时间尺度主问题最优解,寻求综合能源体长时间尺度最小成本调度策略,并作为鲁棒优化模型下界;
3)求解子问题:将主问题优化结果带入子问题,寻求当前调度方案下风、光发电最恶劣申报场景,并再次求解该场景下的调度方案,得到方案后反馈给主问题,将子问题最优目标函数值和主问题最优目标函数值之和作为鲁棒优化模型新的上界;
其中,主问题是综合能源体长时间尺度最小成本调度策略;
子问题是当前调度方案下风、光发电最恶劣申报场景,并再次求解该场景下的调度方案;
当前调度方案下风、光发电最恶劣申报场景为:风、光发电上下限快速变化的场景;
得到方案后反馈给主问题:即在得到该场景下的调度方案后,根据该调度方案对鲁棒优化模型目标函数值上界进行修改,以子问题最优目标函数值和主问题最优目标函数值之和作为鲁棒优化模型新的上界;
4)判断收敛:判断鲁棒优化调度模型子问题与主问题的优化目标函数之差是否小于收敛精度,若满足收敛条件,则停止迭代,输出综合能源系统各设备长时间尺度调度方案;否则,将子问题求得的恶劣场景代入主问题中,并令k=k+1,返回步骤2)继续迭代求解,直到满足收敛条件。
具体来说,子问题求得的恶劣场景为风、光发电上下限快速变化的场景,该场景中风、光预测值和负荷预测值范围缩小,故将缩小后的风、光预测值和负荷预测值作为新的输入,并再次求解;
一般不存在迭代次数达到最大时仍不满足收敛条件;若发生,则返回步骤1重新调整参数。
步骤3,在短时间尺度下,考虑光伏、风电功率与其在长时间尺度下的预测偏差,利用氢储能的快速响应调节能力,建立短时间尺度修正调度模型,求解得到氢储能的修正调度方案,具体地,
短时间尺度修正调度模型的目标函数为:
Figure BDA0003902020000000131
式中:
Figure BDA0003902020000000132
为短时间尺度下因光伏、风电预测偏差带来的购售电量申报偏差惩罚成本;
Figure BDA0003902020000000133
为氢储能运维成本增量;
短时间尺度修正调度模型的约束条件同式(20)、(23)~(26),短时间尺度修正调度模型本质为混合整数线性规划问题,采用基于matlab平台yalmip工具进行建模,采用CPLEX求解器进行求解。
以下是本发明方法的一个算例,长时间尺度下电能、热能与天然气的供需调度结果分别如图2-4所示,短时间尺度下的修正调度结果如图5所示,并对比了仅采用长时间尺度优化调度与所提双时间尺度协同调度下综合能源系统的用能成本差异,如表1所示,可知,所提双时间尺度协同调度的运行成本比单一时间尺度的减少850.88元,验证了本发明提供方法的有效性与良好的经济性。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他任何未背离本发明的精神实质和原理下所作的修改、修饰、替代、组合、简化,均应为等效的置换方式,都应包含在本发明的保护范围之内。

Claims (4)

1.一种基于双时间尺度的综合能源系统优化调度方法,其特征在于包括以下步骤:
步骤1,基于综合能源系统的光伏、风电历史运行数据,建立光伏、风电长时间尺度下的功率波动区间模型;
步骤2,在长时间尺度下,考虑多类型电-热-气调节设备的调节能力,以综合用能成本最低为目标,建立长时间尺度下的鲁棒优化调度模型,求解得到长时间尺度下各调控设备的调度方案;
步骤3,在短时间尺度下,考虑光伏、风电功率与其在长时间尺度下的预测偏差,利用氢储能的快速响应调节能力,建立短时间尺度修正调度模型,求解得到氢储能的修正调度方案。
2.根据权利要求1所述的一种基于双时间尺度的综合能源系统优化调度方法,其特征在于,步骤1所述基于光伏、风电历史数据建立长时间尺度下的光伏、风电功率波动区间模型,具体表示如下:
Figure FDA0003902019990000011
Figure FDA0003902019990000012
Figure FDA0003902019990000013
Figure FDA0003902019990000014
Figure FDA0003902019990000015
Figure FDA0003902019990000016
式中,
Figure FDA0003902019990000017
分别对应长时间尺度下光伏、风电有功功率实际值;
Figure FDA0003902019990000018
分别对应长时间尺度下光伏、风电有功功率预测值;
Figure FDA0003902019990000019
为因高估光伏、风电出力后需额外购买电能的惩罚成本,
Figure FDA0003902019990000021
为高估后购买电能惩罚价格,
Figure FDA0003902019990000022
为高估光伏、风电出力;
Figure FDA0003902019990000023
为因低估光伏、风电出力出售额外电能收益,
Figure FDA0003902019990000024
为低估后售电惩罚价格,
Figure FDA0003902019990000025
为低估光伏、风电出力;
Figure FDA0003902019990000026
分别为光伏、风电有功功率预测偏差;ΓPV、ΓWIN是鲁棒系数,为了避免决策过于保守。
3.根据权利要求1所述的一种基于双时间尺度的综合能源系统优化调度方法,其特征在于,步骤2所述的考虑多类型电-热-气调节设备的调节能力,以综合用能成本最低为目标,建立长时间尺度下的鲁棒优化调度模型,其中:
多类型电-热-气调节设备包括电动汽车、燃气轮机、氢储能、燃气锅炉;
长时间尺度下的鲁棒优化调度模型的目标函数与约束条件具体为:
(1)目标函数:
Figure FDA0003902019990000027
Figure FDA0003902019990000028
Figure FDA0003902019990000029
Figure FDA00039020199900000210
Figure FDA00039020199900000211
Figure FDA00039020199900000212
式中:
Figure FDA00039020199900000213
分别为t时间段内购买天然气成本、向电动汽车用户提供补偿成本、燃气轮机环境成本、综合能源系统向上级电网售电收益以及氢储能运维成本;
Figure FDA0003902019990000031
为天然气价格;
Figure FDA0003902019990000032
为t时间段天然气负荷;
Figure FDA0003902019990000033
为燃气轮机t时间段有功出力;ηgt为燃气轮机气电效率系数;
Figure FDA0003902019990000034
为t时间段锅炉有功出力;ηgb为锅炉气电效率系数;eev为综合能源体向提供放电服务的电动汽车用户提供的补偿成本;Pdisc,t为电动汽车蓄电池t时间段内的放电功率;N为电动汽车总数量;Ql为燃气轮机第l项污染气体的排放强度;Yl、Bl分别为第l项污染气体的罚款成本系数和环境成本系数;
Figure FDA0003902019990000035
为t时间段综合能源体向上级电网售、购电功率大小;
Figure FDA0003902019990000036
分别为综合能源体向上级电网售、购电电价;
Figure FDA0003902019990000037
分别对应t时间段氢储能系统电制氢和氢燃料电池发电功率;ηOMH2为氢储能系统运维成本系数;
(2)约束条件
1)电动汽车调节约束:
Figure FDA0003902019990000038
Figure FDA0003902019990000039
Figure FDA00039020199900000310
Figure FDA00039020199900000311
Figure FDA00039020199900000312
式中:SOCev,t,i为t时间段第i组电动汽车蓄电池荷电状态;
Figure FDA00039020199900000313
Figure FDA00039020199900000314
分别对应第i组电动汽车蓄电池荷电状态下限值和上限值;
Figure FDA00039020199900000315
分别对应t时间段第i组电动汽车充放电功率;
Figure FDA00039020199900000316
分别对应电动汽车蓄电池的充放电效率系数;
Figure FDA00039020199900000317
分别对应第i组电动汽车蓄电池最大充放电功率;0-1变量
Figure FDA0003902019990000041
分别对应第i组电动汽车t时间段充放电标志位,充电则取1,否则取0,两者不能同时为1;
2)燃气轮机调节约束:
Figure FDA0003902019990000042
Figure FDA0003902019990000043
式中:
Figure FDA0003902019990000044
分别对应燃气轮机最小和最大输出功率;
Figure FDA0003902019990000045
分别对应燃气轮机向下、向上爬坡功率;
4)氢储能调节约束:
Figure FDA0003902019990000046
式中:
Figure FDA0003902019990000047
为t时间段电制氢功率;
Figure FDA0003902019990000048
为t时间段氢燃料电池发电功率;
Figure FDA0003902019990000049
分别对应电制氢和氢燃料电池发电的0/1变量,目的为了避免同时制氢和发电;
5)燃气锅炉调节约束:
Figure FDA00039020199900000410
Figure FDA00039020199900000411
式中:
Figure FDA00039020199900000412
为锅炉最小和最大输出功率;
Figure FDA00039020199900000413
为燃气锅炉向下、向上爬坡功率;
6)综合能源系统电能供需平衡约束:
Figure FDA00039020199900000414
式中:
Figure FDA0003902019990000051
分别对应t时间段风力发电和光伏发电功率大小;
Figure FDA0003902019990000052
为t时间段电负荷功率;
7)综合能源系统热能供需平衡约束:
Figure FDA0003902019990000053
式中:
Figure FDA0003902019990000054
为t时间段燃气锅炉输出的热功率;
Figure FDA0003902019990000055
为t时间段热负荷功率大小;
Figure FDA0003902019990000056
为t时间段弃热功率大小;
8)综合能源系统天然气供需平衡约束:
Figure FDA0003902019990000057
式中:
Figure FDA0003902019990000058
为t时间段综合能源主体天然气总消耗量;
Figure FDA0003902019990000059
分别对应t时间段燃气轮机、燃气锅炉耗气量;
Figure FDA00039020199900000510
为t时间段天然气负荷;
9)综合能源系统氢气供需平衡约束:
Figure FDA00039020199900000511
式中:
Figure FDA00039020199900000512
为t时间段储氢罐容量;
Figure FDA00039020199900000513
为t时间段氢气负荷;
Figure FDA00039020199900000514
为每千瓦时电能通过电解槽电解水产生氢气质量;η为每千克氢气通过燃料电池燃烧所产生的电能;θel、θfc分别对应电解槽和燃料电池的对应转换效率;
Figure FDA00039020199900000515
代表氢气的摩尔质量。
4.根据权利要求1所述的一种基于双时间尺度的综合能源系统优化调度方法,其特征在于,步骤3所述的建立短时间尺度修正调度模型,考虑短时间尺度与长时间尺度存在的风电、光伏功率预测偏差,在其他设备调度方案不变的前提下,修正氢储能的调度方案,利用其快速响应能力平抑光伏、风电实时功率波动;
短时间尺度修正调度模型的目标函数为:
Figure FDA0003902019990000061
式中:
Figure FDA0003902019990000062
为短时间尺度下因光伏、风电预测偏差带来的购售电量申报偏差惩罚成本;
Figure FDA0003902019990000063
为氢储能运维成本增量;
短时间尺度修正调度模型的约束条件同式(20)、(23)~(26)。
CN202211294624.XA 2022-10-21 2022-10-21 一种基于双时间尺度的综合能源系统优化调度方法 Active CN115640894B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211294624.XA CN115640894B (zh) 2022-10-21 2022-10-21 一种基于双时间尺度的综合能源系统优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211294624.XA CN115640894B (zh) 2022-10-21 2022-10-21 一种基于双时间尺度的综合能源系统优化调度方法

Publications (2)

Publication Number Publication Date
CN115640894A true CN115640894A (zh) 2023-01-24
CN115640894B CN115640894B (zh) 2024-03-26

Family

ID=84944906

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211294624.XA Active CN115640894B (zh) 2022-10-21 2022-10-21 一种基于双时间尺度的综合能源系统优化调度方法

Country Status (1)

Country Link
CN (1) CN115640894B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090632A (zh) * 2018-01-23 2018-05-29 南方电网科学研究院有限责任公司 基于鲁棒优化的新能源并网电力系统多时间尺度调度方法
CN108521132A (zh) * 2018-04-16 2018-09-11 广西大学 多能互补支撑电网频率调节的多时间尺度优化控制方法
CN110009152A (zh) * 2019-04-03 2019-07-12 东南大学 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法
CN110197312A (zh) * 2019-06-28 2019-09-03 东南大学 一种基于多时间尺度的用户级综合能源系统优化调度方法
WO2019165701A1 (zh) * 2018-02-28 2019-09-06 东南大学 一种交直流混联微网的随机鲁棒耦合型优化调度方法
CN110417006A (zh) * 2019-07-24 2019-11-05 三峡大学 考虑多能协同优化的综合能源系统多时间尺度能量调度方法
WO2022100091A1 (zh) * 2020-11-10 2022-05-19 中国能源建设集团湖南省电力设计院有限公司 一种广义源储系统调度的集中控制方法
CN114781817A (zh) * 2022-03-30 2022-07-22 全球能源互联网集团有限公司 电-气互联系统多时间尺度联合优化调度方法及系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090632A (zh) * 2018-01-23 2018-05-29 南方电网科学研究院有限责任公司 基于鲁棒优化的新能源并网电力系统多时间尺度调度方法
WO2019165701A1 (zh) * 2018-02-28 2019-09-06 东南大学 一种交直流混联微网的随机鲁棒耦合型优化调度方法
CN108521132A (zh) * 2018-04-16 2018-09-11 广西大学 多能互补支撑电网频率调节的多时间尺度优化控制方法
CN110009152A (zh) * 2019-04-03 2019-07-12 东南大学 一种考虑电转气和不确定性的区域综合能源系统运行鲁棒优化方法
CN110197312A (zh) * 2019-06-28 2019-09-03 东南大学 一种基于多时间尺度的用户级综合能源系统优化调度方法
CN110417006A (zh) * 2019-07-24 2019-11-05 三峡大学 考虑多能协同优化的综合能源系统多时间尺度能量调度方法
WO2022100091A1 (zh) * 2020-11-10 2022-05-19 中国能源建设集团湖南省电力设计院有限公司 一种广义源储系统调度的集中控制方法
CN114781817A (zh) * 2022-03-30 2022-07-22 全球能源互联网集团有限公司 电-气互联系统多时间尺度联合优化调度方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何国鑫等: "园区集中式混合储能系统的响应控制策略", 《可再生能源》, no. 09 *
翟晶晶等: "考虑需求响应与光伏不确定性的综合能源系统鲁棒优化", 《中国电力》, vol. 53, no. 08, 5 August 2020 (2020-08-05) *

Also Published As

Publication number Publication date
CN115640894B (zh) 2024-03-26

Similar Documents

Publication Publication Date Title
CN111738497B (zh) 计及需求侧响应的虚拟电厂双层优化调度方法
Moradi et al. Operational strategy optimization in an optimal sized smart microgrid
CN108009748B (zh) 一种电能替代项目节能减排效果评价方法
CN111950808B (zh) 基于综合需求响应的综合能源系统随机鲁棒优化运行方法
CN113792953B (zh) 一种虚拟电厂优化调度方法及系统
CN111861783A (zh) 考虑负荷转移的综合能源系统多目标优化配置方法
CN111900725B (zh) 一种电动汽车充电站与光伏发电站的协同互动方法和系统
CN112990523A (zh) 基于多目标模型预测控制的区域综合能源系统分层优化运行方法
CN115117940A (zh) 考虑风电、光伏及负荷不确定的风光水火储系统低碳调度模型建模方法
CN115207972A (zh) 一种容量电价与风光火配比相协调的电源规划方法
CN113410854B (zh) 一种多类型储能系统优化运行方法
CN115115130A (zh) 一种基于模拟退火算法的风光储制氢系统日前调度方法
CN115115087A (zh) 一种考虑氢燃料汽车与氢储能的虚拟电厂协调调度方法
Geng et al. Optimal allocation model of virtual power plant capacity considering Electric vehicles
CN117254502A (zh) 基于电-氢混合储能综合能源系统多目标优化调度方法
CN115640894A (zh) 一种基于双时间尺度的综合能源系统优化调度方法
CN115995832A (zh) 一种计及柔性负荷的户用氢能系统优化运行方法及装置
CN114784896A (zh) 面向虚拟电厂的规模化充电桩能量优化管理方法和系统
CN114971154A (zh) 一种包含碳交易机制的可再生能源消纳方法
CN113742944A (zh) 一种考虑电制氢系统的虚拟电厂建模方法
CN113449987A (zh) 可再生能源配额制下的园区能源系统优化调度方法
CN116061742B (zh) 一种分时电价光伏园区内电动汽车的充电控制方法和系统
CN112149339A (zh) 风电-光伏-光热-电加热器互补发电系统容量优化模型
CA2511632A1 (en) Energy network
CN117977644A (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