CN110009152A - 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 - Google Patents
一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 Download PDFInfo
- Publication number
- CN110009152A CN110009152A CN201910264889.7A CN201910264889A CN110009152A CN 110009152 A CN110009152 A CN 110009152A CN 201910264889 A CN201910264889 A CN 201910264889A CN 110009152 A CN110009152 A CN 110009152A
- Authority
- CN
- China
- Prior art keywords
- power
- energy
- hydrogen
- gas
- indicate
- 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
- 230000005611 electricity Effects 0.000 title claims abstract description 86
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000005457 optimization Methods 0.000 claims abstract description 39
- 230000008569 process Effects 0.000 claims abstract description 23
- 238000012913 prioritisation Methods 0.000 claims abstract description 15
- 230000000295 complement effect Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims abstract description 7
- 238000003487 electrochemical reaction Methods 0.000 claims abstract description 6
- 239000001257 hydrogen Substances 0.000 claims description 159
- 229910052739 hydrogen Inorganic materials 0.000 claims description 159
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 141
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 106
- 239000007789 gas Substances 0.000 claims description 101
- 239000003345 natural gas Substances 0.000 claims description 48
- 238000006243 chemical reaction Methods 0.000 claims description 36
- 230000006870 function Effects 0.000 claims description 28
- 238000007906 compression Methods 0.000 claims description 20
- 238000004146 energy storage Methods 0.000 claims description 15
- 230000006835 compression Effects 0.000 claims description 13
- 238000003860 storage Methods 0.000 claims description 13
- 238000005265 energy consumption Methods 0.000 claims description 12
- 150000002431 hydrogen Chemical class 0.000 claims description 12
- 238000004519 manufacturing process Methods 0.000 claims description 10
- 238000010438 heat treatment Methods 0.000 claims description 8
- 238000005868 electrolysis reaction Methods 0.000 claims description 7
- 240000002853 Nelumbo nucifera Species 0.000 claims description 6
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 6
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims description 6
- 238000005338 heat storage Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 6
- 230000002596 correlated effect Effects 0.000 claims description 5
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 5
- 238000011112 process operation Methods 0.000 claims description 5
- 238000010977 unit operation Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 claims description 4
- 230000000875 corresponding effect Effects 0.000 claims description 4
- 239000000567 combustion gas Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 3
- 230000021715 photosynthesis, light harvesting Effects 0.000 claims description 3
- 238000012887 quadratic function Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 244000131316 Panax pseudoginseng Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 6
- 238000005516 engineering process Methods 0.000 description 8
- 238000009412 basement excavation Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013439 planning Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 206010033307 Overweight Diseases 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000005804 alkylation reaction Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 239000003054 catalyst Substances 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000036425 denaturation Effects 0.000 description 1
- 238000004925 denaturation Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000003446 ligand Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 235000020825 overweight Nutrition 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000013519 translation Methods 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- 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)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Electrolytic Production Of Non-Metals, Compounds, Apparatuses Therefor (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,包括根据电化学反应的原理对电转气过程详细建模,基于能源中心模型用矩阵的形式分析不同能量形式之间的耦合关系;建立考虑多种负荷预测功率不确定性的区域综合能源***运行优化模型,并根据多面体形式的不确定区间建立两阶段鲁棒优化模型;将两阶段鲁棒优化模型分解为主问题与子问题,将子问题转化为单一目标的优化问题,并进行迭代求解得到综合能源***鲁棒优化方案。本发明充分挖掘多种能源形式之间互补互济、柔性调度的潜力,能够减小弃风,提高综合能源***运行的经济性和灵活性;能够根据区域综合能源***实际情况协调运行风险和成本间的矛盾。
Description
技术领域
本发明属于综合能源***运行优化领域,具体涉及一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法。
背景技术
综合能源***(RIES)被认为是未来人类社会能源的主要承载形式。其通过实现多种异质能源子***之间的协调规划和互补互济,在满足***内多元化用能需求的同时,有效地促进可再生能源就地消纳。由于RIES含有多样化的负荷类型,功率预测不确定性和时变性强,综合能源***能量优化调度方式与传统微电网的调度方法相比面临更大的挑战。
电转气(P2G)技术作为RIES的重要支撑技术之一,通过电-气网络的深度融合,为解决弃风消纳问题提供了有效途径。已有的关于电转气技术在综合能源***中的应用,侧重于利用P2G技术的能量转换和时空平移特性解决可再生能源过剩的问题;或侧重于电-气混联综合能源***的协同规划。当前研究通常将电-气-氢耦合考虑为固定效率,对氢储能建模采用广义储能模型。实际中,作为P2G技术核心能源转换装置的电解槽在运行过程中受输入功率、额定容量等因素影响转换效率并不固定;氢储能的充放过程分为压缩-存储-再压缩3个阶段,压缩机的工作过会带来耗能问题,这与蓄电池的充放电过程有所区别。
区域综合能源***还需考虑供需侧不确定性因素给***调度带来的风险。处理不确定量的优化方法包括:基于概率分布的随机规划、机会约束优化、鲁棒优化等。其中,RIES中不确定变量众多,随机规划、机会约束优化很难获得较为精确的概率分布函数,影响方法的精确性。鲁棒优化只需知道不确定参数的波动范围,通过一个不确定集合表征参数的波动,当在不确定集合中取值时,鲁棒优化模型就一定能得到确定的可行解。由于无需获得其精确的概率分布,计算量明显降低。在经济调度问题中,鲁棒优化方法按照建模思路的不同可分为两阶段鲁棒优化模型和工程博弈模型。目前有关不确定性因素对电力网络的影响及分析方法已有较多研究,相比而言,不确定性因素对综合能源***的影响分析尚处于起步阶段。因此,一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法急需提出。
发明内容
发明目的:为了解决现有区域综合能源***运行技术中对电转气环节和不确定性因素考虑不足的问题,本发明提供一种考虑电转气详细建模和源荷不确定性的区域综合能源***运行鲁棒优化方法。
技术方案:为实现上述发明目的,本发明采用以下技术方案:
一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,包括以下步骤:
(1)根据电化学反应的原理对电转气过程详细建模,基于能源中心模型用矩阵的形式分析不同能量形式之间的耦合关系;
(2)建立考虑风电出力和电、热、天然气、氢气多种负荷预测功率不确定性的区域综合能源***运行优化模型,并根据多面体形式的不确定区间建立两阶段鲁棒优化模型;
(3)将步骤(2)得到的两阶段鲁棒优化模型分解为主问题与子问题,使用KKT条件和bigM方法将子问题转化为单一目标的优化问题,并进行迭代求解得到综合能源***鲁棒优化方案。
进一步的,步骤(1)包括:
(11)根据电化学反应的原理对电转气过程详细建模,其中电解槽制氢气是一切电转气工艺的基础,电解输出氢气的量由如下的数学模型表示:
uEL,e,tPEL,e,min≤PEL,e,t≤uEL,e,tPEL,e,max;
其中,表示电解槽效率函数,是关于电解槽t时刻消耗电功率PEL,e,t与额定电功率PEL,rated比值的二次函数;aEL、bEL、cEL分别表示与电解槽效率相关的常数;EEL,H2,t表示电解槽装置在t时段制得氢气的量;πEL,rated表示电解槽模块的额定容量;uEL,e,t是表征电解槽在t时段启停状态的0-1变量,uEL,e,t=1表示启动,uEL,e,t=0表示停机;PEL,e,max和PEL,e,min表示用于电解制氢的电功率上下限约束;
(12)能源中心(EH)是具有多种输入和输出的能源传输和转化中心,其内部包括能量的传输、转换、存储三种元素,假设输入端有ξ种能源载体,输出端有种能源载体,用EH表示的多能源平衡方程表示如下:
式中,L表示能源中心的负荷向量;P表示能源中心的输入向量;耦合系数Cab表示输入的b能源经过EH内部各种设备转换到a能源的系数,该耦合系数由转换单元的的运行效率和EH的调度系数共同决定。
更进一步的,步骤(11)中电转气详细建模包括:
1)储氢罐与压缩机模型
储氢罐存储模型表达式为:
EHT,t=EHT,t-1+EHT,in,t-EHT,out,t;
EHT,min≤EHT,t≤EHT,max;
其中,EHT,t表示储氢罐在t时段存储氢气的量;EHT,min、EHT,max表示储氢罐存储容量的最小值和最大值;EHT,in,t和EHT,out,t分别为t时段存入、输出储氢罐的氢气的量;
前压缩过程能耗模型与当前时段输入氢气量EHT,in,t相关,关系式为:
Ppre,e,t=λpre,eEHT,in,t;
0≤EHT,in,t≤EHT,pre,max;
其中,Ppre,e,t表示前压缩机需要的电能;λpre,e是前压缩设备的能耗系数;EHT,pre,max表示前压缩机允许的氢气量上限;
储氢罐内部压强数学关系式为:
其中,Ptank,t表示t时段高压容器内部压强;EET,t为t时段储氢罐内氢气的量;R为理想气体常数;Zh表示输入端的压缩系数;Ttank是储氢罐内的热力学温度;V表示储氢罐的容积;
后压缩过程的能耗与高压储氢罐内部压强和输出的氢气流量有关,模型如下:
0≤EHT,out,t≤EHT,post,max;
式中,Wpost,t表示压缩机t时段每压缩1kmol氢气的能耗;Ppost,e,t为后压缩机在t时段消耗的电能;表示输出端压缩系数;Rcomp表示输出端气体常数;k是氢气的热容比常数;Pout是储氢罐输出端压强;ηpost表示后压缩机输出效率;EHT,post,max表示经由后压缩装置输出的氢流量上限;
2)甲烷化反应模型
为简化计算,采用如下的氢气制甲烷模型:
其中,PM,g,t表示t时段甲烷化过程制得的天然气功率;ηm是反应转化效率;表示参与甲烷化反应的氢气输入量;HL表示天然气低热值;4是氢气的摩尔质量,可用来折算出生成天然气的量;κ表示天然气管网内的气体密度;和分别表示参与甲烷化反应的氢气输入量最大、最小值;是表示甲烷化反应在t时段是否启动的0-1变量,为1表示启动,为0表示停机。
进一步的,步骤(2)包括:
(21)建立考虑风电和负荷功率不确定性的区域综合能源***日前优化模型;
日前经济调度的目标函数为:
其中,T为调度的时段数,对于日前优化调度T=24;Pe,t为综合能源***t时段从电网购电的功率;ce,t表示从电网购买的实时电价;Pg,t为综合能源***t时段从天然气网购气的功率;cg,t表示天然气单位售价;为t时段风电出力的预测值;Pwind,e,t表示t时段RIES实际利用的风电功率;β表示弃风惩罚费用系数;
各时段约束条件包括不确定性功率区间约束、***功率平衡约束、外网传输功率限制、可再生能源出力约束、储能单元约束、热电联供机组运行约束、电转气过程运行约束以及锅炉运行约束;
(22)约束条件中建立的不确定集是一个多面体有界集合,不确定变量在该有界集合内任意变化时,需要始终保证优化方案的可行性;选择鲁棒优化方法处理上述考虑不确定性的区域综合能源***日前调度问题,将上述运行优化模型写成两阶段鲁棒模型的形式,如下所示:
式中,F是原目标函数;x、y为优化变量;u为定义的源、荷不确定量,U为不确定集;外层问题的决策变量x是RIES各可控单元的调度方案,即等效EH的耦合矩阵系数;内层问题的决策变量是风电和各类型负荷的波动程度;对于某一个固定的外层问题解x和不确定量u,内层问题的可行域为Ω(x,u);内层问题为max-min形式,其中min部分表示对于某一个固定的x和u求解最小化的内层问题目标,而max则表示对于整个不确定集合中所有的不确定量可能确定最差情况。
更进一步的,约束条件具体为:
1)不确定性功率区间约束
式中,分别为考虑不确定性的风电功率以及电、热、气、氢负荷功率值;为基于日前功率预测的值; 为功率波动系数,即每种功率相对预测值的最大偏差;Γwind、Γe,L、Γh,L、Γg,L、为负荷功率不确定性预算控制参数,Γe,L,Γh,L,Γg,L,具体的取值根据实际情况而定;
2)***功率平衡约束
综合能源***涉及的电能、热能、天然气、氢气四种能源形式,平衡关系式如下:
其中,分别为t时段***电负荷、热负荷、天然气负荷和氢能负荷;PGT,e,t、Pbt,ch,t、Pbt,dis,t、PEB,e,t、PEL,e,t、Ppre,e,t、Ppost,e,t是与电能平衡相关的量,分别表示t时段燃气轮机发电功率、蓄电池充电功率、蓄电池放电功率、电锅炉功率、电解制氢耗电功率、前压缩机耗电功率、后压缩机耗电功率;Ptst,ch,t、Ptst,dis,t、PGB,h,t、PEB,h,t、PGT,h,t是与热能平衡相关的量,分别表示t时段储热单元充热功率、放热功率、燃气锅炉发热功率、电锅炉发热功率、燃气轮机发热功率;PGB,g,t、PGT,g,t、PM,g,t是与***内天然气平衡相关的量,分别表示t时段燃气锅炉用气功率、燃气轮机用气功率、甲烷化反应制得天然气的功率;EHT,in,t、EHT,out,t分别为t时段电解水制得的氢量、参与甲烷化反应用到的氢量、注入储氢罐的氢量、从储氢罐释放的氢量;氢能平衡式中的4表示氢气摩尔质量折算系数;
3)外网传输功率限制
综合能源***与电网、气网的最大交互功率受到网络传输能力的限制:
Pe,min≤Pe,t≤Pe,max;
Pg,min≤Pg,t≤Pg,max;
其中,Pe,max和Pe,min表示RIES与电网联络线功率的最大、最小值;Pg,max和Pg,min表示RIES与天然气***传输功率的最大、最小值;
4)可再生能源出力约束
其中,Pwind,e,t表示RIES实际利用的风电出力,为风力发电机组在t时刻考虑不确定性区间的发电预测值,表征可利用的最大风电出力;
5)储能单元约束
蓄电池和蓄热槽的建模采用储能通用模型;
6)热电联供机组运行约束
燃气轮机的电、热输出功率与燃气输入采用线性数学模型,考虑爬坡约束;
7)电转气过程运行约束
8)锅炉运行约束
电锅炉和燃气锅炉输入、输出采用线性模型,考虑爬坡约束。
进一步的,步骤(3)包括:
(31)采用计算效率和求解性能都较优的列约束生产算法将考虑不确定性的区域综合能源***日前调度模型分解;分解为主问题(MP)和子问题(SP);
主问题:
MP:minη
式中,η为求解目标和s.t.为约束条件,分别表示原问题约束中电能、热能、天然气、氢气能量平衡等式除了不确定变量,其余相关变量之和;区别在于,原问题约束条件中考虑不确定性的变量替换为由子问题求解获得的最恶劣波动情况下的值Pe,L'、Ph,L'、Ph,L'、
主问题是在外网传输功率、储能单元、热电联供机组、锅炉和电转气过程约束条件下以综合能源***日前调度成本最小为目标进行优化,决策变量是各可控单元的机组调度方案;通过子问题不确定性集合中有限个可能的波动场景,将主问题中涉及不确定量的约束条件用部分枚举场景替代,由此可知主问题即为单目标的混合整数线性规划问题;
子问题:
s.t.潮流方程、电压方程、运行上下限和功率区间约束
上式中,燃气锅炉、燃气轮机、电解槽、甲烷化反应元件的机组启停计划已在主问题求出,为已知变量;
子问题目标函数包含max-min双层目标,较难直接求解;由于子问题为线性问题,满足强对偶条件,因此选择增加KKT条件将内层min问题转化为max问题,从而获得单一max问题,增加相应的对偶约束和互补松弛约束并使用big-M方法进行线性处理;转化后的子问题目标函数如下所示:
(32)根据分解原两阶段鲁棒优化模型得到的主问题和子问题,设计迭代方法求解鲁棒优化方案。
更进一步的,迭代方法求解鲁棒优化方案具体步骤如下:
1)设置目标函数值的上界UB=+∞,下界LB=-∞,风力发电和负荷的功率波动初始值为0,当前迭代次数k=0;
2)设置最小公差ε;
3)求解主问题,得到的优化结果包括各个可控单元的调度系数,得到主问题目标函数最优解η(k),并更新下界为LB=max{LB,η(k)};
4)基于当前主问题求解得到的可控单元的调度系数,求解子问题,得到的优化结果包括考虑不确定性情况下风电利用功率、与电网、气网交互功率和最优子问题目标函数值f(k),并更新上界为UB=min{UB,f(k)};
5)判断UB-LB<ε是否成立。若UB-LB<ε,则输出鲁棒优化的结果;若UB-LB≥ε,则令k=k+1,并根据第k次迭代得到的子问题考虑不确定量波动程度最优解获得该场景下风电最大值和负荷功率情况,添加关于风电功率利用和***功率平衡的列约束至第k+1次迭代的主问题,并转到2);具体列约束如下:
式中,表示原问题约束中电能、热能、天然气、氢气能量平衡等式除了不确定变量,其余相关变量之和; 为第k次迭代由子问题确定的不确定量的值,经过上述基于C&CG方法的迭代求解,最终可以得到区域综合能源***运行鲁棒优化方案。
有益效果:与现有技术相比,本发明的方法具有以下优点:
(1)本发明结合电化学原理建立了详细的电转气模型,弥补了当前区域综合能源***运行优化中对于电转气过程考虑过于简单这一不足,充分挖掘电能、热能、天然气、氢能多种能源形式之间互补互济、柔性调度的潜力,提高了综合能源***运行的经济性和灵活性。
(2)本发明所研究***融合了热电联产、电转气技术,并把弃风成本纳入目标函数,有助于实现节能减排、实现能量的梯级利用,丰富了现有的可再生能源消纳技术。
(3)本发明提出的鲁棒优化,同时考虑了源侧的风电不确定性,负荷侧电、热、天然气、氢能多种负荷的预测波动影响,有助于应对RIES多能源特点带来的不确定性增强,能根据***实际情况协调运行风险和成本之间的矛盾。
附图说明
图1为本发明方法的流程图;
图2为区域综合能源***能量耦合关系图;
图3为基于C&CG方法求解的流程图。
具体实施方式
下面对本发明技术方案进行详细说明,但是本发明的保护范围不局限于所述实施例。
如图1所示,本发明是一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,包括以下步骤:
步骤1:根据电化学反应的原理对电转气过程详细建模,基于能源中心模型,分析不同区域综合能源***中不同能量形式之间的耦合关系。具体实现过程如下:
(1)电解槽制氢气是一切电转气工艺的基础,本文选择基于质子交换膜技术的电解槽模型。在实际电解过程中,电解槽的转换效率由电解设备性能参数、输入电能功率等因素共同决定,而不是固定的。为模拟这一过程,电解槽输出氢气的量由如下的数学模型表示:
uEL,e,tPEL,e,min≤PEL,e,t≤uEL,e,tPEL,e,max (3);
其中,表示电解槽效率函数,它是关于电解槽t时刻消耗电功率PEL,e,t(kW)与额定电功率PEL,rated(kW)比值的二次函数;aEL、bEL、cEL(kW)表示与电解槽效率相关的常数;表示电解槽装置在t时段制得氢气的量,为了计算的便利,采用物质的量单位,kmol;πEL,rated表示电解槽模块的额定容量,kmol/h;uEL,e,t是表征电解槽在t时段启停状态的0-1变量,uEL,e,t=1表示启动,uEL,e,t=0表示停机;PEL,e,max和PEL,e,min表示用于电解制氢的电功率上下限约束。
从电转氢的利用方式出发,考虑氢气直接供应氢燃料电池汽车抵消一部分氢负荷,以及氢气进一步制甲烷抵消一部分天然气负荷,这两类用途。详细的电转气模型包括:
1)储氢罐与压缩机模型
本发明采用高压储气罐作为氢气的存储设备,以供应氢能负荷。当氢能应用在燃料电池汽车领域时,对氢气能量密度有一定要求。氢气在高压储气罐内部仍以气态存在,储气量与内部压力成正比,虽然增大压强可以提高氢气的储量,但同时增加了压缩机的耗能。电解产生的氢气进入储氢罐之前需压缩以提升气压,这一过程称为前压缩;输气端压力根据负荷需求也有固定的范围,氢气从储氢罐再经过加压供应氢负荷,称为后压缩过程。因此,储氢罐模型分为前压缩机、存储、后压缩机三部分。
储氢罐存储模型有如下表达式:
EHT,t=EHT,t-1+EHT,in,t-EHT,out,t (4);
EHT,min≤EHT,t≤EHT,max (5);
式中,EHT,t表示储氢罐在t时段存储氢气的量,kmol;EHT,min、EHT,max表示储氢罐存储容量的最小值和最大值。EHT,in,t和EHT,out,t分别为t时段存入、输出储氢罐的氢气的量,kmol。
前压缩过程能耗模型与当前时段输入氢气量EHT,in,t相关,关系式为:
Ppre,e,t=λpre,eEHT,in,t (6);
0≤EHT,in,t≤EHT,pre,max (7);
其中,Ppre,e,t表示前压缩机需要的电能,kW;λpre,e是前压缩设备的能耗系数;EHT,pre,max表示前压缩机允许的氢气量上限,kmol。
储氢罐内部压强有如下数学关系式:
其中,Ptank,t表示t时段高压容器内部压强,bar;EET,t为t时段储氢罐内氢气的量,kmol;R为理想气体常数;Zh表示输入端的压缩系数;Ttank是储氢罐内的热力学温度,K;V表示储氢罐的容积,m3。
后压缩过程的能耗与高压储氢罐内部压强和输出的氢气流量有关,模型如下:
0≤EHT,out,t≤EHT,post,max (11);
式中,Wpost,t表示压缩机t时段每压缩1kmol氢气的能耗,kW/kmol;Ppost,e,t为后压缩机在t时段消耗的电能,kW;表示输出端压缩系数;Rcomp表示输出端气体常数;k是氢气的热容比常数;Pout是储氢罐输出端压强,bar;ηpost表示后压缩机输出效率;EHT,post,max表示经由后压缩装置输出的氢流量上限,kmol。
2)甲烷化反应模型
由氢气制备甲烷时,甲烷化反应的效率受到反应温度、压力、碳氢配比、催化剂等因素的影响。为简化计算,采用如下的氢气制甲烷模型:
其中,PM,g,t表示t时段甲烷化过程制得的天然气功率,kW;ηm是反应转化效率;表示参与甲烷化反应的氢气输入量,kmol;HL表示天然气低热值,kW/m3;4是氢气的摩尔质量,kg/kmol,可用来折算出生成天然气的量;κ表示天然气管网内的气体密度,kg/m3;和分别表示参与甲烷化反应的氢气输入量最大、最小值;是表示甲烷化反应在t时段是否启动的0-1变量,为1表示启动,为0表示停机。
(2)本发明所研究的综合能源***的能量耦合关系如图2所示。负荷形式有电、氢、天然气、热负荷四类;供能来自风电机组、电网和天然气网输入的能量。RIES内有多种能量转化设备:燃气轮机消耗天然气输出电能和热能,燃气锅炉消耗天然气输出热能,电锅炉消耗电能输出热能,P2G消耗电能输出氢能和人工合成天然气。储能单元包括储电、储热和储氢设备。
为便于分析不同能量形式的耦合关系,引入能源中心(Energy Hub,EH)的概念。EH是具有多种输入和输出的能源传输和转化中心,其内部包括能量的传输、转换、存储三种元素。假设输入端有ξ种能源载体,输出端有种能源载体,用EH表示的多能源平衡方程可以表示如下:
式中,L表示能源中心的负荷向量(电,热,气,氢);P表示能源中心的输入向量(电网、气网输入功率,风力机组输入功率);耦合系数Cab表示输入的b能源经过EH内部各种设备转换到a能源的系数,该耦合系数由转换单元的的运行效率和EH的调度系数共同决定。
由此,对于给定的EH负荷向量,可以把综合能源***的调度问题等价为对EH内部调度系数以及能源输入向量的优化。
步骤2:根据步骤1,建立考虑风电和负荷功率不确定性的区域综合能源***日前优化模型,并根据多面体形式的不确定集建立RIES运行两阶段鲁棒优化模型。具体实现过程如下:
(1)为鼓励可再生能源的消纳,减少弃风,RIES的日运行成本包括电能、天然气的购买成本,以及弃风惩罚成本。
日前经济调度的目标函数为:
其中,T为调度的时段数,对于日前优化调度T=24;Pe,t为综合能源***t时段从电网购电的功率;ce,t表示从电网购买的实时电价;Pg,t为综合能源***t时段从天然气网购气的功率;cg,t表示天然气单位售价;为考虑不确定性的t时段风电出力预测值;Pwind,e,t表示t时段RIES实际利用的风电功率;β表示弃风惩罚费用系数。
各时段具体的约束条件如下:
1)不确定性功率区间约束
式中,分别为考虑不确定性的风电功率以及电、热、气、氢负荷功率值;为基于日前功率预测的值; 为功率波动系数,即每种功率相对预测值的最大偏差;Γwind、Γe,L、Γh,L、Γg,L、为负荷功率不确定性控制参数,Γe,L,Γh,L,Γg,L,具体的取值根据实际情况而定。
2)***功率平衡约束
本发明研究的综合能源***涉及电能、热能、天然气、氢气四种能源形式,平衡关系式如下:
其中,分别为t时段***电负荷、热负荷、天然气负荷和氢能负荷;PGT,e,t、Pbt,ch,t、Pbt,dis,t、PEB,e,t、PEL,e,t、Ppre,e,t、Ppost,e,t是与电能平衡相关的量,分别表示t时段燃气轮机发电功率、蓄电池充电功率、蓄电池放电功率、电锅炉功率、电解制氢耗电功率、前压缩机耗电功率、后压缩机耗电功率;Ptst,ch,t、Ptst,dis,t、PGB,h,t、PEB,h,t、PGT,h,t是与热能平衡相关的量,分别表示t时段储热单元充热功率、放热功率、燃气锅炉发热功率、电锅炉发热功率、燃气轮机发热功率;PGB,g,t、PGT,g,t、PM,g,t是与***内天然气平衡相关的量,分别表示t时段燃气锅炉用气功率、燃气轮机用气功率、甲烷化反应制得天然气的功率;EHT,in,t、EHT,out,t分别为t时段电解水制得的氢量、参与甲烷化反应用到的氢量、注入储氢罐的氢量、从储氢罐释放的氢量;氢能平衡式中的4表示氢气摩尔质量折算系数,kg/kmol。
3)外网传输功率限制
综合能源***与电网、气网的最大交互功率受到网络传输能力的限制:
Pe,min≤Pe,t≤Pe,max
Pg,min≤Pg,t≤Pg,max
其中,Pe,max和Pe,min表示RIES与电网联络线功率的最大、最小值;Pg,max和Pg,min表示RIES与天然气***传输功率的最大、最小值。
4)可再生能源出力约束
其中,Pwind,e,t表示RIES实际利用的风电出力,为风力发电机组在t时刻考虑不确定性区间的发电预测值,表征可利用的最大风电出力。
5)储能单元约束
蓄电池和蓄热槽的建模采用储能通用模型。
6)热电联供机组运行约束
燃气轮机的电、热输出功率与燃气输入采用线性数学模型,考虑爬坡约束。
7)电转气过程运行约束:
电解制氢、甲烷化反应、储氢罐模型和约束见步骤1。
8)锅炉运行约束
电锅炉和燃气锅炉输入、输出采用线性模型,考虑爬坡约束
(2)由于约束条件中建立的不确定集是一个多面体有界集合,不确定变量在该有界集合内任意变化时,需要始终保证优化方案的可行性。选择鲁棒优化方法处理上述考虑不确定性的区域综合能源***日前调度问题,将上述运行优化模型写成两阶段鲁棒模型的形式,如下所示:
式中,F是原目标函数;x、y为优化变量;u为定义的源、荷不确定量,U为不确定集。外层问题的决策变量x是RIES各可控单元的调度方案,即等效EH的耦合矩阵系数;内层问题的决策变量是风电和各类型负荷的波动程度。对于某一个固定的外层问题解x和不确定量u,内层问题的可行域为Ω(x,u)。内层问题为max-min形式,其中min部分表示对于某一个固定的x和u求解最小化的内层问题目标,而max则表示对于整个不确定集合中所有的不确定量可能确定最差情况。
步骤3:将步骤2得到的两阶段鲁棒优化模型分解为主问题与子问题,使用KKT条件和bigM方法将子问题转化为单一目标的优化问题,并进行迭代求解得到鲁棒优化方案,具体实现过程如下:
(1)采用计算效率和求解性能都较优的列约束生产算法(C&CG)将考虑不确定性的区域综合能源***日前调度模型分解。根据步骤2建立的两阶段鲁棒优化模型可知,外层问题的决策变量各可控单元的机组调度方案,含有离散变量;内层问题的决策变量是风电和各类型负荷波动情况,还有从电网、气网购进的功率,是连续变量。可将该问题分解为主问题(MP)和子问题(SP)。
主问题为:
MP:minη
式中,求解目标η和约束条件s.t.与步骤2(1)建立的原问题相同,变量定义同原问题。T为调度的时段数,对于日前优化调度T=24;Pe,t为综合能源***t时段从电网购电的功率;ce,t表示从电网购买的实时电价;Pg,t为综合能源***t时段从天然气网购气的功率;cg,t表示天然气单位售价;为考虑不确定性的t时段风电出力预测值;Pwind,e,t表示t时段RIES实际利用的风电功率;β表示弃风惩罚费用系数。 表示原问题约束中电能、热能、天然气、氢气能量平衡等式除了不确定变量,其余相关变量之和。区别在于,原问题约束条件中考虑不确定性的变量替换为由子问题求解获得的最恶劣波动情况下的值Pe,L'、Ph,L'、Ph,L'、
主问题是在外网传输功率、储能单元、热电联供机组、锅炉和电转气过程约束条件下以综合能源***日前调度成本最小为目标进行优化,决策变量是各可控单元的机组调度方案。通过子问题不确定性集合中有限个可能的波动场景,将主问题中涉及不确定量的约束条件用部分枚举场景替代,由此可知主问题即为单目标的混合整数线性规划问题。
子问题为:
s.t.潮流方程、电压方程、运行上下限和功率区间约束
上式中,变量定义与步骤2(1)建立的原问题相同,各可控单元的机组调度方案已在主问题求出,为已知变量。
子问题用于生成新的枚举场景,并以约束条件的形式添加到主问题中。对于任何给定的机组调度方案,子问题都能生成最恶劣的波动场景,求解出相应的最优解,因此子问题为获取目标函数最大值的混合整数线性规划问题。
子问题目标函数包含max-min双层目标,较难直接求解。由于子问题为线性问题,满足强对偶条件,因此选择增加KKT条件将内层min问题转化为max问题,从而获得单一max问题,增加相应的对偶约束和互补松弛约束并使用big-M方法进行线性处理。转化后的子问题目标函数如下所示:
(2)根据分解原两阶段鲁棒优化模型得到的主问题和子问题,设计迭代方法求解鲁棒优化方案,算法流程图如图3所示。具体步骤如下:
1)设置目标函数值的上界UB=+∞,下界LB=-∞,风力发电和负荷的功率波动初始值为0,当前迭代次数k=0。
2)设置最小公差ε。
3)求解主问题,得到的优化结果包括各个可控单元的调度系数,得到主问题目标函数最优解η(k),并更新下界为LB=max{LB,η(k)}。
4)基于当前主问题求解得到的可控单元的调度系数,求解子问题,得到的优化结果包括考虑不确定性情况下风电利用功率、与电网、气网交互功率和最优子问题目标函数值f(k),并更新上界为UB=min{UB,f(k)}。
5)判断UB-LB<ε是否成立。若UB-LB<ε,则输出鲁棒优化的结果;若UB-LB≥ε,则令k=k+1,并根据第k次迭代得到的子问题考虑不确定量波动程度最优解获得该场景下风电最大值和负荷功率情况,添加关于风电功率利用和***功率平衡的列约束至第k+1次迭代的主问题,并转到2)。具体列约束如下:
式中,表示原问题约束中电能、热能、天然气、氢气能量平衡等式除了不确定变量,其余相关变量之和。 为第k次迭代由子问题确定的不确定量的值。经过上述基于C&CG方法的迭代求解,最终可以得到区域综合能源***运行鲁棒优化方案。
本发明的一种考虑电转气和不确定性的区域综合能源***(RIES)运行鲁棒优化方法。包括对电转气(P2G)过程详细建模,基于能源中心模型用矩阵的形式分析不同能量形式之间的耦合关系;建立考虑风电和负荷功率不确定性的区域综合能源***日前调度模型,并根据多面体形式的不确定区间建立RIES两阶段鲁棒优化模型;利用列约束生产算法(C&CG)将两阶段鲁棒优化模型分解为主问题与子问题,迭代求解得到区域综合能源***鲁棒优化方案。本发明弥补了当前区域综合能源***运行研究中对于电转气过程考虑过于简单这一不足,充分挖掘电能、热能、天然气、氢能多种能源形式之间互补互济、柔性调度的潜力,能够减小弃风,提高综合能源***运行的经济性和灵活性;本发明考虑了源侧的风电不确定性,负荷侧电、热、天然气、氢能多种负荷的预测波动影响,能够根据区域综合能源***实际情况协调运行风险和成本间的矛盾。
Claims (7)
1.一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,包括以下步骤:
(1)根据电化学反应的原理对电转气过程详细建模,基于能源中心模型用矩阵的形式分析不同能量形式之间的耦合关系;
(2)建立考虑风电出力和电、热、天然气、氢气多种负荷预测功率不确定性的区域综合能源***运行优化模型,并根据多面体形式的不确定区间建立两阶段鲁棒优化模型;
(3)将步骤(2)得到的两阶段鲁棒优化模型分解为主问题与子问题,使用KKT条件和bigM方法将子问题转化为单一目标的优化问题,并进行迭代求解得到综合能源***鲁棒优化方案。
2.根据权利要求1所述的一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,步骤(1)包括:
(11)根据电化学反应的原理对电转气过程详细建模,其中电解槽制氢气是一切电转气工艺的基础,电解输出氢气的量由如下的数学模型表示:
uEL,e,tPEL,e,min≤PEL,e,t≤uEL,e,tPEL,e,max;
其中,表示电解槽效率函数,是关于电解槽t时刻消耗电功率PEL,e,t与额定电功率PEL,rated比值的二次函数;aEL、bEL、cEL分别表示与电解槽效率相关的常数;EEL,H2,t表示电解槽装置在t时段制得氢气的量;πEL,rated表示电解槽模块的额定容量;uEL,e,t是表征电解槽在t时段启停状态的0-1变量,uEL,e,t=1表示启动,uEL,e,t=0表示停机;PEL,e,max和PEL,e,min表示用于电解制氢的电功率上下限约束;
(12)能源中心(EH)是具有多种输入和输出的能源传输和转化中心,其内部包括能量的传输、转换、存储三种元素,假设输入端有ξ种能源载体,输出端有种能源载体,用EH表示的多能源平衡方程表示如下:
式中,L表示能源中心的负荷向量;P表示能源中心的输入向量;耦合系数Cab表示输入的b能源经过EH内部各种设备转换到a能源的系数,该耦合系数由转换单元的的运行效率和EH的调度系数共同决定。
3.根据权利要求2所述的一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,步骤(11)中电转气详细建模包括:
1)储氢罐与压缩机模型
储氢罐存储模型表达式为:
EHT,t=EHT,t-1+EHT,in,t-EHT,out,t;
EHT,min≤EHT,t≤EHT,max;
其中,EHT,t表示储氢罐在t时段存储氢气的量;EHT,min、EHT,max表示储氢罐存储容量的最小值和最大值;EHT,in,t和EHT,out,t分别为t时段存入、输出储氢罐的氢气的量;
前压缩过程能耗模型与当前时段输入氢气量EHT,in,t相关,关系式为:
Ppre,e,t=λpre,eEHT,in,t;
0≤EHT,in,t≤EHT,pre,max;
其中,Ppre,e,t表示前压缩机需要的电能;λpre,e是前压缩设备的能耗系数;EHT,pre,max表示前压缩机允许的氢气量上限;
储氢罐内部压强数学关系式为:
其中,Ptank,t表示t时段高压容器内部压强;EET,t为t时段储氢罐内氢气的量;R为理想气体常数;Zh表示输入端的压缩系数;Ttank是储氢罐内的热力学温度;V表示储氢罐的容积;
后压缩过程的能耗与高压储氢罐内部压强和输出的氢气流量有关,模型如下:
0≤EHT,out,t≤EHT,post,max;
式中,Wpost,t表示压缩机t时段每压缩1kmol氢气的能耗;Ppost,e,t为后压缩机在t时段消耗的电能;表示输出端压缩系数;Rcomp表示输出端气体常数;k是氢气的热容比常数;Pout是储氢罐输出端压强;ηpost表示后压缩机输出效率;EHT,post,max表示经由后压缩装置输出的氢流量上限;
2)甲烷化反应模型
为简化计算,采用如下的氢气制甲烷模型:
其中,PM,g,t表示t时段甲烷化过程制得的天然气功率;ηm是反应转化效率;表示参与甲烷化反应的氢气输入量;HL表示天然气低热值;4是氢气的摩尔质量,用来折算出生成天然气的量;κ表示天然气管网内的气体密度;和分别表示参与甲烷化反应的氢气输入量最大、最小值;是表示甲烷化反应在t时段是否启动的0-1变量,为1表示启动,为0表示停机。
4.根据权利要求1所述的一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,步骤(2)包括:
(21)建立考虑风电和负荷功率不确定性的区域综合能源***日前优化模型;
日前经济调度的目标函数为:
其中,T为调度的时段数,对于日前优化调度T=24;Pe,t为综合能源***t时段从电网购电的功率;ce,t表示从电网购买的实时电价;Pg,t为综合能源***t时段从天然气网购气的功率;cg,t表示天然气单位售价;为t时段风电出力的预测值;Pwind,e,t表示t时段RIES实际利用的风电功率;β表示弃风惩罚费用系数;
各时段约束条件包括不确定性功率区间约束、***功率平衡约束、外网传输功率限制、可再生能源出力约束、储能单元约束、热电联供机组运行约束、电转气过程运行约束以及锅炉运行约束;
(22)约束条件中建立的不确定集是一个多面体有界集合,不确定变量在该有界集合内任意变化时,需要始终保证优化方案的可行性;选择鲁棒优化方法处理上述考虑不确定性的区域综合能源***日前调度问题,将上述运行优化模型写成两阶段鲁棒模型的形式,如下所示:
式中,F是原目标函数;x、y为优化变量;u为定义的源、荷不确定量,U为不确定集;外层问题的决策变量x是RIES各可控单元的调度方案,即等效EH的耦合矩阵系数;内层问题的决策变量是风电和各类型负荷的波动程度;对于某一个固定的外层问题解x和不确定量u,内层问题的可行域为Ω(x,u);内层问题为max-min形式,其中min部分表示对于某一个固定的x和u求解最小化的内层问题目标,而max则表示对于整个不确定集合中所有的不确定量可能确定最差情况。
5.根据权利要求4所述的一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,约束条件具体为:
1)不确定性功率区间约束
式中,分别为考虑不确定性的风电功率以及电、热、气、氢负荷功率值;为基于日前功率预测的值; 为功率波动系数,即每种功率相对预测值的最大偏差;Γwind、Γe,L、Γh,L、Γg,L、为负荷功率不确定性预算控制参数,Γe,L,Γh,L,Γg,L,具体的取值根据实际情况而定;
2)***功率平衡约束
综合能源***涉及的电能、热能、天然气、氢气四种能源形式,平衡关系式如下:
其中,分别为t时段***电负荷、热负荷、天然气负荷和氢能负荷;PGT,e,t、Pbt,ch,t、Pbt,dis,t、PEB,e,t、PEL,e,t、Ppre,e,t、Ppost,e,t是与电能平衡相关的量,分别表示t时段燃气轮机发电功率、蓄电池充电功率、蓄电池放电功率、电锅炉功率、电解制氢耗电功率、前压缩机耗电功率、后压缩机耗电功率;Ptst,ch,t、Ptst,dis,t、PGB,h,t、PEB,h,t、PGT,h,t是与热能平衡相关的量,分别表示t时段储热单元充热功率、放热功率、燃气锅炉发热功率、电锅炉发热功率、燃气轮机发热功率;PGB,g,t、PGT,g,t、PM,g,t是与***内天然气平衡相关的量,分别表示t时段燃气锅炉用气功率、燃气轮机用气功率、甲烷化反应制得天然气的功率;EHT,in,t、EHT,out,t分别为t时段电解水制得的氢量、参与甲烷化反应用到的氢量、注入储氢罐的氢量、从储氢罐释放的氢量;氢能平衡式中的4表示氢气摩尔质量折算系数;
3)外网传输功率限制
综合能源***与电网、气网的最大交互功率受到网络传输能力的限制:
Pe,min≤Pe,t≤Pe,max;
Pg,min≤Pg,t≤Pg,max;
其中,Pe,max和Pe,min表示RIES与电网联络线功率的最大、最小值;Pg,max和Pg,min表示RIES与天然气***传输功率的最大、最小值;
4)可再生能源出力约束
其中,Pwind,e,t表示RIES实际利用的风电出力,为风力发电机组在t时刻考虑不确定性区间的发电预测值,表征可利用的最大风电出力;
5)储能单元约束
蓄电池和蓄热槽的建模采用储能通用模型;
6)热电联供机组运行约束
燃气轮机的电、热输出功率与燃气输入采用线性数学模型,考虑爬坡约束;
7)电转气过程运行约束
8)锅炉运行约束
电锅炉和燃气锅炉输入、输出采用线性模型,考虑爬坡约束。
6.根据权利要求1所述的考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,步骤(3)包括:
(31)采用计算效率和求解性能都较优的列约束生产算法将考虑不确定性的区域综合能源***日前调度模型分解;分解为主问题(MP)和子问题(SP);
主问题:
式中,η为求解目标和s.t.为约束条件,分别表示原问题约束中电能、热能、天然气、氢气能量平衡等式除了不确定变量,其余相关变量之和;区别在于,原问题约束条件中考虑不确定性的变量替换为由子问题求解获得的最恶劣波动情况下的值Pe,L'、Ph,L'、Ph,L'、主问题是在外网传输功率、储能单元、热电联供机组、锅炉和电转气过程约束条件下以综合能源***日前调度成本最小为目标进行优化,决策变量是各可控单元的机组调度方案;通过子问题不确定性集合中有限个可能的波动场景,将主问题中涉及不确定量的约束条件用部分枚举场景替代,由此可知主问题即为单目标的混合整数线性规划问题;
子问题:
SP:
s.t.潮流方程、电压方程、运行上下限和功率区间约束
上式中,燃气锅炉、燃气轮机、电解槽、甲烷化反应元件的机组启停计划已在主问题求出,为已知变量;
子问题目标函数包含max-min双层目标,较难直接求解;由于子问题为线性问题,满足强对偶条件,因此选择增加KKT条件将内层min问题转化为max问题,从而获得单一max问题,增加相应的对偶约束和互补松弛约束并使用big-M方法进行线性处理;转化后的子问题目标函数如下所示:
SP:
(32)根据分解原两阶段鲁棒优化模型得到的主问题和子问题,设计迭代方法求解鲁棒优化方案。
7.根据权利要求6所述的考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法,其特征在于,迭代方法求解鲁棒优化方案具体步骤如下:
1)设置目标函数值的上界UB=+∞,下界LB=-∞,风力发电和负荷的功率波动初始值为0,当前迭代次数k=0;
2)设置最小公差ε;
3)求解主问题,得到的优化结果包括各个可控单元的调度系数,得到主问题目标函数最优解η(k),并更新下界为LB=max{LB,η(k)};
4)基于当前主问题求解得到的可控单元的调度系数,求解子问题,得到的优化结果包括考虑不确定性情况下风电利用功率、与电网、气网交互功率和最优子问题目标函数值f(k),并更新上界为UB=min{UB,f(k)};
5)判断UB-LB<ε是否成立;若UB-LB<ε,则输出鲁棒优化的结果;若UB-LB≥ε,则令k=k+1,并根据第k次迭代得到的子问题考虑不确定量波动程度最优解获得该场景下风电最大值和负荷功率情况,添加关于风电功率利用和***功率平衡的列约束至第k+1次迭代的主问题,并转到2);具体列约束如下:
式中,表示原问题约束中电能、热能、天然气、氢气能量平衡等式除了不确定变量,其余相关变量之和; 为第k次迭代由子问题确定的不确定量的值,经过上述基于C&CG方法的迭代求解,最终得到区域综合能源***运行鲁棒优化方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910264889.7A CN110009152B (zh) | 2019-04-03 | 2019-04-03 | 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910264889.7A CN110009152B (zh) | 2019-04-03 | 2019-04-03 | 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110009152A true CN110009152A (zh) | 2019-07-12 |
CN110009152B CN110009152B (zh) | 2021-12-07 |
Family
ID=67169619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910264889.7A Active CN110009152B (zh) | 2019-04-03 | 2019-04-03 | 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110009152B (zh) |
Cited By (59)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110490386A (zh) * | 2019-08-26 | 2019-11-22 | 苏州树森信息科技有限公司 | 一种综合能源调度方法和综合能源调度*** |
CN110503250A (zh) * | 2019-08-08 | 2019-11-26 | 燕山大学 | 考虑电-热转移负荷转移量不确定性的综合能源规划方法 |
CN110555618A (zh) * | 2019-09-06 | 2019-12-10 | 东北大学 | 基于改进樽海鞘算法的网络化综合能源***优化调度方法 |
CN110598904A (zh) * | 2019-08-08 | 2019-12-20 | 华北电力大学 | 市场环境下考虑可再生能源消纳的车网能量交互优化方法 |
CN110729722A (zh) * | 2019-10-22 | 2020-01-24 | 东北大学 | 一种基于能源细胞架构的区域能源网优化调度方法 |
CN110729766A (zh) * | 2019-09-27 | 2020-01-24 | 国网河北省电力有限公司 | 一种多区域电-气综合能源***分布式鲁棒优化方法 |
CN110807560A (zh) * | 2019-11-12 | 2020-02-18 | 国网重庆市电力公司电力科学研究院 | 氢混天然气能源***调度方法、装置及可读存储介质 |
CN111008723A (zh) * | 2019-09-24 | 2020-04-14 | 华北电力大学 | 一种用于分布式能源pev充电站设计的优化方法 |
CN111027846A (zh) * | 2019-12-06 | 2020-04-17 | 东南大学 | 一种考虑热氢联产的电氢综合能源***及其容量配置方法 |
CN111030110A (zh) * | 2019-12-28 | 2020-04-17 | 福州大学 | 一种考虑电转气消纳风电的电力-天然气耦合***鲁棒协同调度方法 |
CN111049135A (zh) * | 2019-12-30 | 2020-04-21 | 国网吉林省电力有限公司 | 一种多区域电气耦合***的分布式两阶段合作运行方法 |
CN111049134A (zh) * | 2019-12-26 | 2020-04-21 | 清华大学 | 一种多能互补园区响应电力***削峰需求的计算方法 |
CN111144620A (zh) * | 2019-12-06 | 2020-05-12 | 东南大学 | 一种考虑季节储氢的电氢综合能源***及其鲁棒规划方法 |
CN111159842A (zh) * | 2019-11-29 | 2020-05-15 | 国网江苏省电力有限公司经济技术研究院 | 一种提高综合能源***运行灵活性的校验方法 |
CN111161809A (zh) * | 2019-12-10 | 2020-05-15 | 南京华盾电力信息安全测评有限公司 | 基于电转氢的能源***操作灵活性提升方法 |
CN111259309A (zh) * | 2020-01-14 | 2020-06-09 | 沈阳工业大学 | 一种基于电转气和燃料电池的多能源***稳定性优化方法 |
CN111313442A (zh) * | 2020-03-06 | 2020-06-19 | 东南大学 | 考虑富氧气体利用的综合能源***经济性提升方法 |
CN111476394A (zh) * | 2020-02-28 | 2020-07-31 | 浙江工业大学 | 一种适用于电热气等多能源***的鲁棒运行优化方法 |
CN111509713A (zh) * | 2020-05-18 | 2020-08-07 | 无锡隆玛科技股份有限公司 | 区域综合能源***模型优化配置方法及*** |
CN111523249A (zh) * | 2020-05-27 | 2020-08-11 | 福州大学 | 一种考虑源荷不确定性的数据驱动鲁棒优化调度方法 |
CN111626464A (zh) * | 2020-04-03 | 2020-09-04 | 国网江苏省电力有限公司常州供电分公司 | 考虑多余氢气和氧气售卖的综合能源***经济性提升方法 |
CN111639824A (zh) * | 2020-06-20 | 2020-09-08 | 东北电力大学 | 一种含电转气的区域综合能源***热电优化调度方法 |
CN111695793A (zh) * | 2020-05-29 | 2020-09-22 | 浙江大学 | 一种综合能源***用能灵活性评价方法及*** |
CN111737873A (zh) * | 2020-06-24 | 2020-10-02 | 广东电网有限责任公司电力调度控制中心 | 基于电解铝负荷提供辅助服务的机组组合方法及相关装置 |
CN111898795A (zh) * | 2020-06-09 | 2020-11-06 | 东南大学 | 考虑多运行模式储氢罐的电制氢/甲烷***优化调度方法 |
CN111950807A (zh) * | 2020-08-26 | 2020-11-17 | 华北电力大学(保定) | 计及不确定性与需求响应的综合能源***优化运行方法 |
CN111950808A (zh) * | 2020-08-26 | 2020-11-17 | 华北电力大学(保定) | 基于综合需求响应的综合能源***随机鲁棒优化运行方法 |
CN112101756A (zh) * | 2020-09-02 | 2020-12-18 | 合肥工业大学 | 数据驱动的多能互补能量枢纽负荷优化调度方法和*** |
CN112234632A (zh) * | 2020-09-30 | 2021-01-15 | 国网新疆电力有限公司经济技术研究院 | 一种季节性氢储能规划方法 |
CN112257229A (zh) * | 2020-09-18 | 2021-01-22 | 西安理工大学 | 一种微网两阶段鲁棒调度方法 |
CN112263766A (zh) * | 2020-10-23 | 2021-01-26 | 子量量子十星医疗仪器(上海)有限公司 | 一种大容量吸氢氧机 |
CN112290533A (zh) * | 2020-09-17 | 2021-01-29 | 国家电网公司西南分部 | 一种氢能-天然气混合储能的综合能源微电网调度方法 |
CN112465307A (zh) * | 2020-11-09 | 2021-03-09 | 青海大学 | 工业园区综合能源配置***及方法 |
CN112633575A (zh) * | 2020-12-22 | 2021-04-09 | 山东大学 | 基于鲁棒优化的多能互补综合能源***容量配置方法及*** |
CN112668188A (zh) * | 2020-12-31 | 2021-04-16 | 福州大学 | 多园区综合能源***的分布式鲁棒协同优化调度方法 |
CN112818559A (zh) * | 2021-02-23 | 2021-05-18 | 上海交通大学 | 基于随机微分方程的区域综合能源连续调度方法及其*** |
CN112966855A (zh) * | 2021-02-09 | 2021-06-15 | 西安理工大学 | 考虑风电不确定性的电-气耦合能源配网协调优化方法 |
CN112989594A (zh) * | 2021-03-09 | 2021-06-18 | 南方电网科学研究院有限责任公司 | 一种考虑氢能的综合能源***运行优化方法 |
CN113013872A (zh) * | 2021-01-29 | 2021-06-22 | 国网江苏省电力有限公司 | 一种考虑社区供暖负荷灵活性的配电网协同优化运行方法 |
CN113052498A (zh) * | 2021-04-23 | 2021-06-29 | 国核电力规划设计研究院有限公司 | 基于综合能源***的电转气两阶段转化装置调度方法 |
CN113065682A (zh) * | 2020-01-02 | 2021-07-02 | 中国电力科学研究院有限公司 | 一种气电综合能源***调度方法及装置 |
CN113065707A (zh) * | 2021-04-08 | 2021-07-02 | 国网甘肃省电力公司平凉供电公司 | 能源调度方法及装置 |
CN113095791A (zh) * | 2021-04-29 | 2021-07-09 | 长沙理工大学 | 一种综合能源***运行方法及*** |
CN113256045A (zh) * | 2020-08-04 | 2021-08-13 | 四川大学 | 考虑风光不确定性的园区综合能源***日前经济调度方法 |
CN113381400A (zh) * | 2021-05-28 | 2021-09-10 | 国网青海省电力公司 | 储能消纳新能源能力评估方法和装置 |
CN113393054A (zh) * | 2021-07-05 | 2021-09-14 | 华北电力大学 | 一种风储联合***的优化调度方法及优化调度*** |
CN113452020A (zh) * | 2021-02-08 | 2021-09-28 | 国网江苏省电力有限公司电力科学研究院 | 一种考虑柔性氢需求的电氢能源***调度方法 |
CN113659580A (zh) * | 2021-08-12 | 2021-11-16 | 国网福建省电力有限公司 | 一种区域电量传输可行域确定方法及*** |
CN113783226A (zh) * | 2021-08-27 | 2021-12-10 | 中国科学院广州能源研究所 | 一种海上风电制氢并网发电***的分层预测能量管理方法 |
CN113837589A (zh) * | 2021-09-18 | 2021-12-24 | 清华大学 | 一种综合能源***n-1调度方法及装置 |
CN113850474A (zh) * | 2021-08-26 | 2021-12-28 | 西南交通大学 | 一种热电氢多能流综合能源***及其优化调度方法 |
CN113902182A (zh) * | 2021-09-29 | 2022-01-07 | 大连理工大学 | 考虑能耗不确定性的电动公交车队鲁棒充电优化方法 |
CN113988714A (zh) * | 2021-11-26 | 2022-01-28 | 湖南大学 | 计及多重不确定性的园区综合能源***动态规划方法、设备及介质 |
CN114358383A (zh) * | 2021-12-03 | 2022-04-15 | 湖南大学 | 一种复杂不确定场景下多能源微网鲁棒优化方法和*** |
CN114819700A (zh) * | 2022-05-12 | 2022-07-29 | 国网能源研究院有限公司 | 区域能源互联网分布式鲁棒经济调度方法及装置 |
CN115173444A (zh) * | 2022-08-11 | 2022-10-11 | 内蒙古工业大学 | 一种风光可再生能源耦合氢储综合能源*** |
CN115640894A (zh) * | 2022-10-21 | 2023-01-24 | 山东大学 | 一种基于双时间尺度的综合能源***优化调度方法 |
CN116027757A (zh) * | 2023-03-27 | 2023-04-28 | 中国科学院电工研究所 | 光催化-热化学耦合制燃料的综合能源***运行调控方法 |
CN116316717A (zh) * | 2023-02-10 | 2023-06-23 | 河海大学 | 一种电-氢综合能源***机会约束调度方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130346040A1 (en) * | 2012-06-20 | 2013-12-26 | Schlumberger Technology Corporation | Shale gas production forecasting |
CN106447152A (zh) * | 2016-06-21 | 2017-02-22 | 中国南方电网有限责任公司电网技术研究中心 | 基于能源中心的电‑气‑热***协同调度方法和*** |
CN107994571A (zh) * | 2017-12-04 | 2018-05-04 | 国网山东省电力公司德州供电公司 | 一种含多模糊参数的区域能源网多时间尺度管理方法 |
CN108599154A (zh) * | 2018-05-14 | 2018-09-28 | 东南大学 | 一种考虑不确定性预算的三相不平衡配电网鲁棒动态重构方法 |
CN108832665A (zh) * | 2018-07-04 | 2018-11-16 | 四川大学 | 一种考虑风电不确定性的电热综合***分布式鲁棒协调优化调度模型 |
CN109193668A (zh) * | 2018-10-31 | 2019-01-11 | 四川大学 | 一种基于分布鲁棒优化的合同电量分解方法 |
-
2019
- 2019-04-03 CN CN201910264889.7A patent/CN110009152B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130346040A1 (en) * | 2012-06-20 | 2013-12-26 | Schlumberger Technology Corporation | Shale gas production forecasting |
CN106447152A (zh) * | 2016-06-21 | 2017-02-22 | 中国南方电网有限责任公司电网技术研究中心 | 基于能源中心的电‑气‑热***协同调度方法和*** |
CN107994571A (zh) * | 2017-12-04 | 2018-05-04 | 国网山东省电力公司德州供电公司 | 一种含多模糊参数的区域能源网多时间尺度管理方法 |
CN108599154A (zh) * | 2018-05-14 | 2018-09-28 | 东南大学 | 一种考虑不确定性预算的三相不平衡配电网鲁棒动态重构方法 |
CN108832665A (zh) * | 2018-07-04 | 2018-11-16 | 四川大学 | 一种考虑风电不确定性的电热综合***分布式鲁棒协调优化调度模型 |
CN109193668A (zh) * | 2018-10-31 | 2019-01-11 | 四川大学 | 一种基于分布鲁棒优化的合同电量分解方法 |
Non-Patent Citations (4)
Title |
---|
XINGQUANXIAO 等: "A Two-Stage Distributionally Robust Coordinated Dispatch for Integrated Electricity and Natural-gas Energy Systems Considering Uncertainty of Wind Power", 《2018 IOP CONFERENCE SERIES: MATERIALS SCIENCE AND ENGINEERING》 * |
徐俊俊等: "考虑分布式电源不确定性的配电网鲁棒动态重构", 《中国电机工程学报》 * |
朱兰等: "计及电转气精细化模型的综合能源***鲁棒随机优化调度", 《电网技术》 * |
李杨等: "含电转气的电-气-热***协同调度与消纳风电效益分析", 《电网技术》 * |
Cited By (101)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110598904B (zh) * | 2019-08-08 | 2023-06-27 | 华北电力大学 | 市场环境下考虑可再生能源消纳的车网能量交互优化方法 |
CN110503250A (zh) * | 2019-08-08 | 2019-11-26 | 燕山大学 | 考虑电-热转移负荷转移量不确定性的综合能源规划方法 |
CN110598904A (zh) * | 2019-08-08 | 2019-12-20 | 华北电力大学 | 市场环境下考虑可再生能源消纳的车网能量交互优化方法 |
CN110490386B (zh) * | 2019-08-26 | 2022-12-06 | 苏州树森信息科技有限公司 | 一种综合能源调度方法和综合能源调度*** |
CN110490386A (zh) * | 2019-08-26 | 2019-11-22 | 苏州树森信息科技有限公司 | 一种综合能源调度方法和综合能源调度*** |
CN110555618A (zh) * | 2019-09-06 | 2019-12-10 | 东北大学 | 基于改进樽海鞘算法的网络化综合能源***优化调度方法 |
CN110555618B (zh) * | 2019-09-06 | 2023-01-13 | 东北大学 | 基于改进樽海鞘算法的网络化综合能源***优化调度方法 |
CN111008723A (zh) * | 2019-09-24 | 2020-04-14 | 华北电力大学 | 一种用于分布式能源pev充电站设计的优化方法 |
CN111008723B (zh) * | 2019-09-24 | 2023-12-22 | 华北电力大学 | 一种用于分布式能源pev充电站设计的优化方法 |
CN110729766A (zh) * | 2019-09-27 | 2020-01-24 | 国网河北省电力有限公司 | 一种多区域电-气综合能源***分布式鲁棒优化方法 |
CN110729722A (zh) * | 2019-10-22 | 2020-01-24 | 东北大学 | 一种基于能源细胞架构的区域能源网优化调度方法 |
CN110729722B (zh) * | 2019-10-22 | 2023-11-07 | 东北大学 | 一种基于能源细胞架构的区域能源网优化调度方法 |
CN110807560A (zh) * | 2019-11-12 | 2020-02-18 | 国网重庆市电力公司电力科学研究院 | 氢混天然气能源***调度方法、装置及可读存储介质 |
CN110807560B (zh) * | 2019-11-12 | 2022-10-04 | 国网重庆市电力公司电力科学研究院 | 氢混天然气能源***调度方法、装置及可读存储介质 |
CN111159842A (zh) * | 2019-11-29 | 2020-05-15 | 国网江苏省电力有限公司经济技术研究院 | 一种提高综合能源***运行灵活性的校验方法 |
CN111159842B (zh) * | 2019-11-29 | 2022-08-02 | 国网江苏省电力有限公司经济技术研究院 | 一种提高综合能源***运行灵活性的校验方法 |
CN111144620A (zh) * | 2019-12-06 | 2020-05-12 | 东南大学 | 一种考虑季节储氢的电氢综合能源***及其鲁棒规划方法 |
CN111027846A (zh) * | 2019-12-06 | 2020-04-17 | 东南大学 | 一种考虑热氢联产的电氢综合能源***及其容量配置方法 |
CN111161809A (zh) * | 2019-12-10 | 2020-05-15 | 南京华盾电力信息安全测评有限公司 | 基于电转氢的能源***操作灵活性提升方法 |
CN111161809B (zh) * | 2019-12-10 | 2021-05-25 | 南京华盾电力信息安全测评有限公司 | 基于电转氢的能源***操作灵活性提升方法 |
CN111049134A (zh) * | 2019-12-26 | 2020-04-21 | 清华大学 | 一种多能互补园区响应电力***削峰需求的计算方法 |
CN111030110B (zh) * | 2019-12-28 | 2021-05-18 | 福州大学 | 一种考虑电转气消纳风电的电力-天然气耦合***鲁棒协同调度方法 |
CN111030110A (zh) * | 2019-12-28 | 2020-04-17 | 福州大学 | 一种考虑电转气消纳风电的电力-天然气耦合***鲁棒协同调度方法 |
CN111049135B (zh) * | 2019-12-30 | 2023-02-03 | 国网吉林省电力有限公司 | 一种多区域电气耦合***的分布式两阶段合作运行方法 |
CN111049135A (zh) * | 2019-12-30 | 2020-04-21 | 国网吉林省电力有限公司 | 一种多区域电气耦合***的分布式两阶段合作运行方法 |
CN113065682B (zh) * | 2020-01-02 | 2024-06-07 | 中国电力科学研究院有限公司 | 一种气电综合能源***调度方法及装置 |
CN113065682A (zh) * | 2020-01-02 | 2021-07-02 | 中国电力科学研究院有限公司 | 一种气电综合能源***调度方法及装置 |
CN111259309B (zh) * | 2020-01-14 | 2023-05-26 | 沈阳工业大学 | 一种基于电转气和燃料电池的多能源***稳定性优化方法 |
CN111259309A (zh) * | 2020-01-14 | 2020-06-09 | 沈阳工业大学 | 一种基于电转气和燃料电池的多能源***稳定性优化方法 |
CN111476394B (zh) * | 2020-02-28 | 2022-07-22 | 浙江工业大学 | 一种适用于电热气等多能源***的鲁棒运行优化方法 |
CN111476394A (zh) * | 2020-02-28 | 2020-07-31 | 浙江工业大学 | 一种适用于电热气等多能源***的鲁棒运行优化方法 |
CN111313442B (zh) * | 2020-03-06 | 2022-06-07 | 东南大学 | 考虑富氧气体利用的综合能源***经济性提升方法 |
CN111313442A (zh) * | 2020-03-06 | 2020-06-19 | 东南大学 | 考虑富氧气体利用的综合能源***经济性提升方法 |
CN111626464A (zh) * | 2020-04-03 | 2020-09-04 | 国网江苏省电力有限公司常州供电分公司 | 考虑多余氢气和氧气售卖的综合能源***经济性提升方法 |
CN111509713A (zh) * | 2020-05-18 | 2020-08-07 | 无锡隆玛科技股份有限公司 | 区域综合能源***模型优化配置方法及*** |
CN111509713B (zh) * | 2020-05-18 | 2023-05-23 | 无锡隆玛科技股份有限公司 | 区域综合能源***模型优化配置方法及*** |
CN111523249B (zh) * | 2020-05-27 | 2022-06-03 | 福州大学 | 一种考虑源荷不确定性的数据驱动鲁棒优化调度方法 |
CN111523249A (zh) * | 2020-05-27 | 2020-08-11 | 福州大学 | 一种考虑源荷不确定性的数据驱动鲁棒优化调度方法 |
CN111695793B (zh) * | 2020-05-29 | 2022-11-29 | 浙江大学 | 一种综合能源***用能灵活性评价方法及*** |
CN111695793A (zh) * | 2020-05-29 | 2020-09-22 | 浙江大学 | 一种综合能源***用能灵活性评价方法及*** |
CN111898795B (zh) * | 2020-06-09 | 2024-06-04 | 东南大学 | 考虑多运行模式储氢罐的电制氢/甲烷***优化调度方法 |
CN111898795A (zh) * | 2020-06-09 | 2020-11-06 | 东南大学 | 考虑多运行模式储氢罐的电制氢/甲烷***优化调度方法 |
CN111639824A (zh) * | 2020-06-20 | 2020-09-08 | 东北电力大学 | 一种含电转气的区域综合能源***热电优化调度方法 |
CN111639824B (zh) * | 2020-06-20 | 2022-05-24 | 东北电力大学 | 一种含电转气的区域综合能源***热电优化调度方法 |
CN111737873A (zh) * | 2020-06-24 | 2020-10-02 | 广东电网有限责任公司电力调度控制中心 | 基于电解铝负荷提供辅助服务的机组组合方法及相关装置 |
CN113256045A (zh) * | 2020-08-04 | 2021-08-13 | 四川大学 | 考虑风光不确定性的园区综合能源***日前经济调度方法 |
CN113256045B (zh) * | 2020-08-04 | 2023-04-07 | 四川大学 | 考虑风光不确定性的园区综合能源***日前经济调度方法 |
CN111950807A (zh) * | 2020-08-26 | 2020-11-17 | 华北电力大学(保定) | 计及不确定性与需求响应的综合能源***优化运行方法 |
CN111950808A (zh) * | 2020-08-26 | 2020-11-17 | 华北电力大学(保定) | 基于综合需求响应的综合能源***随机鲁棒优化运行方法 |
CN111950807B (zh) * | 2020-08-26 | 2022-03-25 | 华北电力大学(保定) | 计及不确定性与需求响应的综合能源***优化运行方法 |
CN112101756A (zh) * | 2020-09-02 | 2020-12-18 | 合肥工业大学 | 数据驱动的多能互补能量枢纽负荷优化调度方法和*** |
CN112101756B (zh) * | 2020-09-02 | 2022-09-30 | 合肥工业大学 | 数据驱动的多能互补能量枢纽负荷优化调度方法和*** |
CN112290533A (zh) * | 2020-09-17 | 2021-01-29 | 国家电网公司西南分部 | 一种氢能-天然气混合储能的综合能源微电网调度方法 |
CN112257229B (zh) * | 2020-09-18 | 2024-04-16 | 西安理工大学 | 一种微网两阶段鲁棒调度方法 |
CN112257229A (zh) * | 2020-09-18 | 2021-01-22 | 西安理工大学 | 一种微网两阶段鲁棒调度方法 |
CN112234632B (zh) * | 2020-09-30 | 2023-12-22 | 国网新疆电力有限公司经济技术研究院 | 一种季节性氢储能规划方法 |
CN112234632A (zh) * | 2020-09-30 | 2021-01-15 | 国网新疆电力有限公司经济技术研究院 | 一种季节性氢储能规划方法 |
CN112263766A (zh) * | 2020-10-23 | 2021-01-26 | 子量量子十星医疗仪器(上海)有限公司 | 一种大容量吸氢氧机 |
CN112263766B (zh) * | 2020-10-23 | 2022-12-06 | 子量量子十星医疗仪器(上海)有限公司 | 一种大容量吸氢氧机 |
CN112465307A (zh) * | 2020-11-09 | 2021-03-09 | 青海大学 | 工业园区综合能源配置***及方法 |
CN112633575B (zh) * | 2020-12-22 | 2023-01-31 | 山东大学 | 基于鲁棒优化的多能互补综合能源***容量配置方法及*** |
CN112633575A (zh) * | 2020-12-22 | 2021-04-09 | 山东大学 | 基于鲁棒优化的多能互补综合能源***容量配置方法及*** |
CN112668188A (zh) * | 2020-12-31 | 2021-04-16 | 福州大学 | 多园区综合能源***的分布式鲁棒协同优化调度方法 |
CN112668188B (zh) * | 2020-12-31 | 2023-02-24 | 福州大学 | 多园区综合能源***的分布式鲁棒协同优化调度方法 |
CN113013872B (zh) * | 2021-01-29 | 2022-09-27 | 国网江苏省电力有限公司 | 一种考虑社区供暖负荷灵活性的配电网协同优化运行方法 |
CN113013872A (zh) * | 2021-01-29 | 2021-06-22 | 国网江苏省电力有限公司 | 一种考虑社区供暖负荷灵活性的配电网协同优化运行方法 |
CN113452020A (zh) * | 2021-02-08 | 2021-09-28 | 国网江苏省电力有限公司电力科学研究院 | 一种考虑柔性氢需求的电氢能源***调度方法 |
CN112966855A (zh) * | 2021-02-09 | 2021-06-15 | 西安理工大学 | 考虑风电不确定性的电-气耦合能源配网协调优化方法 |
CN112966855B (zh) * | 2021-02-09 | 2023-10-24 | 西安理工大学 | 考虑风电不确定性的电-气耦合能源配网协调优化方法 |
CN112818559A (zh) * | 2021-02-23 | 2021-05-18 | 上海交通大学 | 基于随机微分方程的区域综合能源连续调度方法及其*** |
CN112818559B (zh) * | 2021-02-23 | 2022-06-28 | 上海交通大学 | 基于随机微分方程的区域综合能源连续调度方法及其*** |
CN112989594A (zh) * | 2021-03-09 | 2021-06-18 | 南方电网科学研究院有限责任公司 | 一种考虑氢能的综合能源***运行优化方法 |
CN112989594B (zh) * | 2021-03-09 | 2023-05-19 | 南方电网科学研究院有限责任公司 | 一种考虑氢能的综合能源***运行优化方法 |
CN113065707B (zh) * | 2021-04-08 | 2024-02-02 | 国网甘肃省电力公司平凉供电公司 | 能源调度方法及装置 |
CN113065707A (zh) * | 2021-04-08 | 2021-07-02 | 国网甘肃省电力公司平凉供电公司 | 能源调度方法及装置 |
CN113052498A (zh) * | 2021-04-23 | 2021-06-29 | 国核电力规划设计研究院有限公司 | 基于综合能源***的电转气两阶段转化装置调度方法 |
CN113052498B (zh) * | 2021-04-23 | 2024-04-05 | 国核电力规划设计研究院有限公司 | 基于综合能源***的电转气两阶段转化装置调度方法 |
CN113095791B (zh) * | 2021-04-29 | 2024-02-27 | 长沙理工大学 | 一种综合能源***运行方法及*** |
CN113095791A (zh) * | 2021-04-29 | 2021-07-09 | 长沙理工大学 | 一种综合能源***运行方法及*** |
CN113381400A (zh) * | 2021-05-28 | 2021-09-10 | 国网青海省电力公司 | 储能消纳新能源能力评估方法和装置 |
CN113393054A (zh) * | 2021-07-05 | 2021-09-14 | 华北电力大学 | 一种风储联合***的优化调度方法及优化调度*** |
CN113393054B (zh) * | 2021-07-05 | 2023-11-24 | 华北电力大学 | 一种风储联合***的优化调度方法及优化调度*** |
CN113659580A (zh) * | 2021-08-12 | 2021-11-16 | 国网福建省电力有限公司 | 一种区域电量传输可行域确定方法及*** |
CN113659580B (zh) * | 2021-08-12 | 2022-08-05 | 国网福建省电力有限公司 | 一种区域电量传输可行域确定方法及*** |
CN113850474A (zh) * | 2021-08-26 | 2021-12-28 | 西南交通大学 | 一种热电氢多能流综合能源***及其优化调度方法 |
CN113850474B (zh) * | 2021-08-26 | 2023-04-07 | 西南交通大学 | 一种热电氢多能流综合能源***及其优化调度方法 |
CN113783226A (zh) * | 2021-08-27 | 2021-12-10 | 中国科学院广州能源研究所 | 一种海上风电制氢并网发电***的分层预测能量管理方法 |
CN113837589B (zh) * | 2021-09-18 | 2024-06-21 | 清华大学 | 一种综合能源***n-1调度方法及装置 |
CN113837589A (zh) * | 2021-09-18 | 2021-12-24 | 清华大学 | 一种综合能源***n-1调度方法及装置 |
CN113902182B (zh) * | 2021-09-29 | 2022-09-20 | 大连理工大学 | 考虑能耗不确定性的电动公交车队鲁棒充电优化方法 |
CN113902182A (zh) * | 2021-09-29 | 2022-01-07 | 大连理工大学 | 考虑能耗不确定性的电动公交车队鲁棒充电优化方法 |
CN113988714B (zh) * | 2021-11-26 | 2024-05-31 | 湖南大学 | 计及多重不确定性的园区综合能源***动态规划方法、设备及介质 |
CN113988714A (zh) * | 2021-11-26 | 2022-01-28 | 湖南大学 | 计及多重不确定性的园区综合能源***动态规划方法、设备及介质 |
CN114358383A (zh) * | 2021-12-03 | 2022-04-15 | 湖南大学 | 一种复杂不确定场景下多能源微网鲁棒优化方法和*** |
CN114819700A (zh) * | 2022-05-12 | 2022-07-29 | 国网能源研究院有限公司 | 区域能源互联网分布式鲁棒经济调度方法及装置 |
CN115173444A (zh) * | 2022-08-11 | 2022-10-11 | 内蒙古工业大学 | 一种风光可再生能源耦合氢储综合能源*** |
CN115640894B (zh) * | 2022-10-21 | 2024-03-26 | 山东大学 | 一种基于双时间尺度的综合能源***优化调度方法 |
CN115640894A (zh) * | 2022-10-21 | 2023-01-24 | 山东大学 | 一种基于双时间尺度的综合能源***优化调度方法 |
CN116316717B (zh) * | 2023-02-10 | 2024-05-28 | 河海大学 | 一种电-氢综合能源***机会约束调度方法 |
CN116316717A (zh) * | 2023-02-10 | 2023-06-23 | 河海大学 | 一种电-氢综合能源***机会约束调度方法 |
CN116027757A (zh) * | 2023-03-27 | 2023-04-28 | 中国科学院电工研究所 | 光催化-热化学耦合制燃料的综合能源***运行调控方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110009152B (zh) | 2021-12-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110009152A (zh) | 一种考虑电转气和不确定性的区域综合能源***运行鲁棒优化方法 | |
Wu et al. | Optimal scheduling for microgrids with hydrogen fueling stations considering uncertainty using data-driven approach | |
CN105703369B (zh) | 一种多能耦合输配网最优能量流建模及求解方法 | |
CN113095791B (zh) | 一种综合能源***运行方法及*** | |
Zhang et al. | A multi-agent deep reinforcement learning approach enabled distributed energy management schedule for the coordinate control of multi-energy hub with gas, electricity, and freshwater | |
CN106773704A (zh) | 多***联合优化调度方法及装置 | |
CN108009693A (zh) | 基于两级需求响应的并网微电网双层优化方法 | |
CN108206543A (zh) | 一种基于能源梯级利用的能源路由器及其运行优化方法 | |
CN113890023A (zh) | 一种综合能源微网分布式经济调度优化方法及*** | |
CN112068436B (zh) | 工业园区的综合能源***分层分布式控制方法及*** | |
CN114358431A (zh) | 计及供需双向需求响应的多能源***优化调度方法及装置 | |
CN114374232A (zh) | 一种考虑需求响应的综合能源***两阶段优化调度方法 | |
CN114077934A (zh) | 一种综合能源微网互联***及其调度方法 | |
Li et al. | Intraday multi-objective hierarchical coordinated operation of a multi-energy system | |
CN115455709A (zh) | 计及碳捕集设备加装的低碳综合能源***仿真及配置方法 | |
Si et al. | Cost-effective operation of the urban energy system with variable supply and demand via coordination of multi-energy flows | |
Lu et al. | Two-stage robust scheduling and real-time load control of community microgrid with multiple uncertainties | |
CN109976155A (zh) | 参与气电市场的虚拟电厂内部随机优化控制方法及*** | |
CN113298407A (zh) | 一种工业园区电-气综合能源***优化调度模型建立方法 | |
CN116502921A (zh) | 一种园区综合能源***优化管理***及其协调调度方法 | |
CN109149630A (zh) | 一种新能源电量消纳能力分析规划方法 | |
CN115495906A (zh) | 高比例可再生能源接入的配网多能互补协同优化方法 | |
CN107273619B (zh) | 一种钢铁企业静态能量流网络优化设计方法 | |
CN112116131A (zh) | 计及碳排放的综合能源***多层级优化方法 | |
Lust et al. | Case Study of a Hydrogen-Based District Heating in a Rural Area: Modeling and Evaluation of Prediction and Optimization Methodologies |
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 |