考虑季节性氢储和燃氢轮机利用的综合能源系统优化方法与流程

1.本发明涉及综合能源系统规划运行技术领域，尤其涉及一种考虑季节性氢储和燃氢轮机利用的综合能源系统优化方法。


背景技术：

2.随着新型能源结构的升级改造，以及为了促进低碳目标的实现，传统火电机组的装机容量占比不断被风电与光伏发电稀释，一定程度上削弱了系统的安全稳定性。为应对电力系统的调节能力不足的问题，国家不断出台针对储能的相关政策指导意见，电化学储能、抽水蓄能等储能形式的研究实践日渐深入。但风电、光伏、水电以及负荷用电特性受自然环境影响具有明显季节性，能源系统在冬夏季节供不应求，而在春秋季节供大于求，在特殊极端天气情况下，新能源可能连续多日处于低出力状态，进一步扩大供能缺口。抽水蓄能以及锂电池储能受制于建设地理要求高、存储容量与经济性制约、储存能量耗散率高、后续回收利用难等多种因素，难以参与周、月等更长时间尺度的运行优化中，因此需要进一步研究大规模、长时间、宽领域的季节性储能(seasonal energy storage，ses)。
3.季节性储能主要关键特征为长时间尺度储能、跨能源形式、空间可运输3个方面，随着电解水制氢技术的不断进步，利用可再生能源电解水制氢成本不断降低，氢气作为能量转化以及二氧化碳处理和资源化利用过程中关键因素之一，不仅可以满足可再生能源规模化、长周期储能需求，而且作为无污染的绿色燃料可以有效与电能深度耦合，降低化石能源的使用，减少碳排量，因此氢储是具有应用潜力的季节性储能。曹军文、覃祥富、耿嘎等在《石油学报(石油加工)》第37卷第6期第1461-1478页(2021年)发表的题为《氢气储运技术的发展现状与展望》一文中总结归纳了氢气目前的储运方式及技术研究现状，uchman w,skorek-osikowska a,jurczyk m等在《能源》(energy)第213卷(2020年)发表的题为《the analysis of dynamic operation of power-to-sng system with hydrogen generator powered with renewable energy,hydrogen storage and methanation unit》一文中评估了氢储对电力系统时域分析影响，张红、袁铁江、谭捷在《中国电机工程学报》第40卷第10期第3364-3372页(2021年)发表的题为《统一能源系统氢负荷中长期预测》一文中分析新能源电制-储氢耦合的能源系统结构形态、利用模式的改变影响；pan g,gu w,lu y等在《ieee可持续能源交易》(ieee transactions on sustainable energy)第11卷第4期(2020年)发表的题为《the analysis of dynamic operation of power-to-sng system with hydrogen generator powered with renewable energy,hydrogen storage and methanation unit》一文中提出了一种考虑发电负荷不确定性、机组n-1安全性和电氢柔性转换过程的电氢综合能源系统鲁棒规划模型，分析了一年4个季节的中长时间尺度季节性氢储的影响；侯慧、刘鹏、黄亮等在《电工技术学报》第36卷增刊1第133-144页(2021年)发表的题为《考虑不确定性的电-热-氢综合能源系统规划》一文中结合chp运行模式切换策略，提出的电-热-氢综合能源系统的多目标规划模型；丁剑、方晓松、宋云亭等在《电力系统自动化》第45卷第24期第1-9页(2021年)年发表的题为《碳中和背景下西部新能源传输的电氢
综合能源网构想》一文中针对未来中国西部高比例新能源电源对外输送及消纳问题，提出了近期和远期西部电氢综合能源电力网的初步构建方案。目前大多数含氢的综合能源运行优化策略研究普遍聚焦在较小的园区内，部分研究内容使用氢能燃料电池进行氢转电，从而减小系统峰谷差，但考虑到电-氢-电的能量转换过程中能量损失率太高且燃料电池本身造价成本高昂等弊端，难以实际运用在大容量的综合能源系统中，需要寻求更好的氢电耦合方式。
4.燃气轮机是促进多能耦合的重要能源设备，李海波、潘志明、黄耀文等在《电力设备管理》2020年第8期第94-96页发表的题为《浅析氢燃料燃气轮机发电的应用前景》一文中分析了氢燃料燃气轮机发电的技术特点，秦锋、秦亚迪、单彤文等在《广东电力》第34卷第10期第10-16页(2021年)发表的题为《碳中和背景下氢燃料燃气轮机技术现状及发展前景》一文中指出混合氢的天然气燃烧更充分，产生的能量更多且化学产物副作用更小。目前各大燃机厂商纷纷展开富氢燃料甚至是纯氢燃料燃机的研发投产，通用电气发电公司已有多个混合氢燃料燃机项目投入商业运行中。由此可见，混合燃料的燃气机组具有承担多能耦合且高效使用氢燃料的巨大潜力。氢气替代天然气出力是能源转型过程中重要的减碳排方式之一，结合碳捕集、利用与封存(carbon capture,utilization and storage，ccus)技术，能够实现从产生源头到最终使用全过程的碳管控，因此需要在综合能源系统构架增加考虑碳流的耦合的研究。
5.因此，本领域的技术人员致力于开发一种考虑季节性氢储和燃氢轮机利用的综合能源系统优化方法，以克服上述存在的问题。


技术实现要素：

6.有鉴于现有技术的上述缺陷，本发明所要解决的技术问题是如何解决可再生能源出力和负荷需求在中长时间尺度上具有季节性倒挂的问题，以及提供一种更好的氢电耦合方式。
7.为实现上述目的，本发明提供了一种考虑季节性氢储和燃氢轮机利用的综合能源系统优化方法，所述方法包括以下步骤：
8.步骤1、提供含碳流的多能耦合综合能源系统，所述综合能源系统包括风电、光伏、水电、火电机组、传统燃气机组、混合燃料燃气机组、电解水装置、储氢装置、氢气甲烷化装置、碳捕集-封存装置、吸收式制冷机、电制冷设备、燃气锅炉；
9.步骤2、建立以系统改造升级成本、运行成本以及惩罚成本综合最小为目标函数的综合能源系统双层规划-运行模型；
10.步骤3、利用改进差分进化算法进行求解，达到配置最优设备容量的同时优化不同机组运行策略。
11.进一步地，所述步骤1中所述电解水装置能量转换模型为：
[0012][0013][0014][0015]
[0016]
式中，分别是电解水装置在t时刻的耗电功率、产氢容量、产热功率，为电解水装置电制氢效率和余热利用效率，为电-氢气折算单位换算系数，为氢气燃烧热值142500kj/m3，表示电解水装置最大制氢容量，α为千瓦与千焦每小时的转换系数取3600。
[0017]
进一步地，所述步骤1中所述储氢装置采用季节性储氢方式，其约束模型为：
[0018][0019][0020][0021]sshs
(0)＝0.5q
shs
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0022][0023]
0≤s
shs
(t)≤q
shs
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0024]
式中，分别表示在t时刻储氢装置储存和释放的功率，分别表示在t时刻充入和放出的0-1状态量，v
shs-max
表示储氢装置最大功率，q
shs
表示储氢装置最大容量，s
shs
(0)、s
shs
(t)、s
shs
(t-1)分别为储氢装置存储能量初始值、t时刻剩余能量、t-1时刻剩余能量，分别表示储氢装置充入和放出的效率，
△
t表示储氢装置充放的单位时间。
[0025]
进一步地，所述步骤1中所述氢气甲烷化能量转换模型为：
[0026][0027][0028]
ξ
e-gas
＝q
gas
/α
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0029][0030][0031]
式中，表示在t时刻氢气甲烷化容量以及氢气甲烷化反应过程余热的功率，分别表示在t时刻氢气甲烷化消耗的氢气量和二氧化碳量，ω1表示反应过程中二氧化碳气体混入比例，表示氢气制甲烷的甲烷和热能转换效率，为氢气甲烷化装置的最大容量，ξ
e-gas
为电-天然气折算单位换算系数，q
gas
为天然气燃烧热值33486.8kj/m3。
[0032]
进一步地，所述步骤2中建立的综合能源系统双层规划-运行模型包括上层规划投资成本目标函数和下层运行优化成本目标函数；
[0033]
所述上层规划投资成本为：
[0034][0035]
λ
crf
＝r
·
(1 r)y/((1 r)
y-1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0036]
式中，c
inv
表示上层规划投资成本，分别为混合燃料
的燃气机组改造升级、电解水装置、储氢装置功率、储氢装置容量、氢气甲烷化装置的单位容量投资成本；p
gt-max
、v
shs-max
、q
shs
、表示上层规划模型需求得最优的系统燃气机组改造升级最大容量、电解水装置最大制氢容量、储氢装置最大功率、储氢装置最大容量、氢气甲烷化装置的最大容量，提供给下层模型优化约束条件；λ
crf
为资金回收系数，r为年利率，y为系统设计平均寿命期限；
[0037]
所述下层运行优化成本包括运行成本c
op
和惩罚成本c
pw
，
[0038]
所述运行成本c
op
包括燃料成本c
opf
、其他装置运行成本c
opa
和机组启动成本c
opu
，其中火电机组的燃料成本进行分段线性化处理：
[0039]cop
＝c
opf
c
opa
c
opu
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0040][0041][0042][0043][0044]
式中：c
gas
为天然气价格成本，为在t时刻系统购买的天然气量；a1、a2、a3、b1、b2、b3为火电机组发电成本分段线性化的系数，cg′
en,i
(t)为在t时刻火电机组i的燃料成本变量；p
gen,i
(t)为火电机组i在t时刻的出力功率；为电解水装置、储氢装置、氢气甲烷化装置、碳捕集-封存装置、吸收式制冷机的单位容量运行成本；为火电机组、混和燃料燃气机组、传统燃气机组的启动成本，为在t时刻火电机组、混和燃料燃气机组、传统燃气机组的启动状态变量，n
gen
为火电机组的数量，为在t时刻碳捕集后封存的二氧化碳气体量、吸收式制冷机发出冷功率，
[0045]
所述惩罚成本c
pw
为弃风光水的惩罚：
[0046][0047]
式中，c
cut
为弃风光水惩罚单位成本，为在t时刻弃风、弃光、弃水功率。
[0048]
进一步地，所述综合能源系统双层规划-运行模型总的目标函数为：
[0049]
minc
total
＝c
inv
c
op
c
pw
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(24)
[0050]
式中，c
total
为总成本。
[0051]
进一步地，所述综合能源系统双层规划-运行模型的约束条件包括：
[0052]
系统电、热、冷、氢功率平衡方程：
[0053][0054][0055][0056][0057][0058]
式中，le(t)、lh(t)、l
co
(t)、分别表示电、热、冷、氢常规负荷在t时刻的需求量，
[0059]
式(25)表示综合能源系统的电功率平衡，为火电机组i在t时刻净出力，p
pv
(t)、p
wt
(t)、p
hp
(t)为风电、光伏、水电机组在t时刻出力，为传统燃气机组n在t时刻净出力，为改造后混合燃料燃气机组在t时刻净出力，为电解水装置、电制冷设备在t时刻消耗的电功率，
[0060]
式(26)表示综合能源系统的天然气量平衡，为传统燃气机组n在t时刻燃烧的天然气量，为改造后混合燃料燃气机组t时刻燃烧的天然气量，为燃气锅炉t时刻燃烧的天然气量，为氢气甲烷化装置在t时刻制得的天然气量，
[0061]
式(27)表示综合能源系统的热功率供需平衡，为传统燃气机组n在t时刻燃烧余热功率的收集转换功率，为改造后混合燃料燃气机组在t时刻燃烧余热功率的收集转换功率，为在t时刻氢气甲烷化和电解水反应过程中热能的收集利用功率，为燃气锅炉在t时刻补充热功率，为在t时刻吸收式制冷机设备吸收的热功率，为在t时刻系统热能无法有效利用而弃掉的功率，
[0062]
式(28)表示综合能源系统的冷功率供需平衡，为在t时刻电制冷设备以及吸收式制冷机设备的制冷功率，
[0063]
式(29)表示综合能源系统的制-储-用氢量的平衡，表示在t时刻电解水反应产生的氢气量，表示在t时刻氢气甲烷化和改造后混合燃料燃气机组消耗的氢气量，表示在t时刻季节性氢储释放和储存的氢气量；
[0064]
可再生能源出力约束：
[0065][0066][0067]
[0068]
为在t时刻风电、光伏、水电预测最大出力值；
[0069]
火电机组出力约束：
[0070]ugen,i
(t)p
gen-min,i
≤p
gen,i
(t)≤u
gen,i
(t)p
gen-max,i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(33)
[0071][0072][0073][0074]
式中，p
gen-max,i
、p
gen-min,i
分别表示火电机组i出力的最大值、最小值，u
gen,i
(t)表示在t时刻火电机组i的状态变量，为在t时刻火电机组i启动状态变量，为在t时刻火电机组i停机状态变量；
[0075]
燃气机组出力约束，包括混合燃料燃气机组出力约束：
[0076][0077][0078][0079][0080][0081][0082]
式中，分别表示混合燃料燃气机组在t时刻输出的电功率和热功率，分别表示混合燃料燃气机组在t时刻天然气、氢气输入体积流量，ξ
e-gas
、分别表示混合燃料燃气机组电-天然气、电-氢气折算单位换算系数，ω2表示混合燃料燃气机组氢气和天然气的混合系数；分别表示混合燃料燃气机组转换电能和热能的效率；u
gtc
(t)、为0-1的状态量，分别表示在t时刻混合燃料燃气机组的状态变量、启动状态变量、停机状态变量，p
gtc-max
、p
gtc-min
为混合燃料燃气机组的最大、最小功率；
[0083]
碳捕集-封存装置约束：
[0084][0085][0086][0087][0088]
式中，为燃气机组在t时刻产生的电力功率，为燃气机组在t时刻净输出电功率，为碳捕集-封存装置在t时刻运行能耗，为碳捕集-封存装置在t时刻基础固定能耗，表示碳捕集-封存装置在t时刻捕集二氧化碳气体量，
表示碳捕集-封存装置在t时刻捕集后封存的二氧化碳气体量，表示碳捕集-封存装置在t时刻捕集后加入甲烷化反应的二氧化碳气体量，为捕集单位碳的运行能耗系数，为捕集效率，为燃气机组单位能量碳排放强度；
[0089]
冷热能量补充机组约束，所述冷热能量补充机组约束包括吸收式制冷机约束、电制冷设备约束、燃气锅炉约束。
[0090]
进一步地，所述吸收式制冷机约束为：
[0091][0092][0093]
式中，表示吸收式制冷机在t时刻发出冷功率，表示吸收式制冷机在t时刻吸收的热功率，表示吸收式制冷机能量转换效率，表示吸收式制冷机能转化发出最大冷功率，u
ac
(t)为表示吸收式制冷机在t时刻启停状态的0-1变量；
[0094]
所述电制冷设备约束为：
[0095][0096][0097]
式中，表示电制冷设备在t时刻发出冷功率，表示电制冷设备在t时刻消耗的电功率，表示电制冷设备能量转换效率，表示电制冷设备转化发出最大冷功率，u
fr
(t)为表示电制冷设备在t时刻启停状态的0-1变量；
[0098]
所述燃气锅炉约束为：
[0099][0100][0101]
式中，表示燃气锅炉在t时刻发出热功率，表示燃气锅炉在t时刻消耗的电功率，表示燃气锅炉能量转换效率，表示燃气锅炉转化发出最大热功率，u
gb
(t)为表示燃气锅炉在t时刻启停状态的0-1变量。
[0102]
进一步地，所述步骤3中所述改进差分进化算法具体包括以下步骤：
[0103]
步骤3.1、初始化：确定规划容量的边界范围，设置种群数量n
p
，随机生成初始种群，其中，种群个体为：
[0104][0105]
步骤3.2、变异：
[0106][0107]
f＝2
λ
f0ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(55)
[0108][0109]
式中，为第g代中随机选取的3个不同个体，为变异种群中的个体，f0为初始设定的变异参数；g表示当前进化代数，gm表示最大进化代数；
[0110]
步骤3.3、交叉：
[0111][0112]
式中，为交叉后所得种群中第i个体的第n维变量；cr为交叉因子，取值[0,1]之间；
[0113]
步骤3.4、竞争：
[0114][0115]
式中，为对应个体的适应度函数，即模型中综合成本目标函数，采用最小化问题的选择方式；
[0116]
步骤3.5、当g》gm时，算法终止，得到最优解，否则，g＝g 1，返回步骤3.2，进行下一次优化。
[0117]
进一步地，所述步骤3中所述改进差分进化算法采用gurobi求解器配合求解。
[0118]
本发明的有益效果为：本发明在提出含碳流的多能耦合综合能源系统框架的基础上，建立了季节性氢储关键设备数学模型，富余可再生能源制氢不仅满足常规氢负荷需求之外，还可以输送到混合燃料的燃气轮机机组进行发电，以及配合捕集的二氧化碳转换为天然气；其次，在传统综合能源系统基础上，建立了以升级改造设备的规划投资以及有效反映季节特性的年度运行成本最小的目标函数模型，结合改进的差分进化算法，进行最优升级改造的容量配置，并分析了中长时间尺度上多种机组运行策略、季节性制-储-用氢能和二氧化碳捕集利用的过程。本文提出的综合能源系统可以有效应对净负荷季节性的波动，促进可再生能源的消纳，减少整体系统碳排量。
[0119]
以下将结合附图对本发明的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本发明的目的、特征和效果。
附图说明
[0120]
图1是本发明的一个较佳实施例的含碳流的多能耦合综合能源系统框架图；
[0121]
图2是本发明的一个较佳实施例的系统优化的改进差分进化算法的求解流程图。
具体实施方式
[0122]
以下参考说明书附图介绍本发明的优选实施例，使其技术内容更加清楚和便于理解。本发明可以通过许多不同形式的实施例来得以体现，本发明的保护范围并非仅限于文中提到的实施例。
[0123]
本发明在考虑季节性氢储和混合氢燃料的燃气轮机原理的基础上建立了含碳流的多能耦合综合能源系统优化方法。
[0124]
为了促进双碳目标的实现，新能源渗透率不断提高，加剧了季节性供需不平衡的矛盾，因此需对新型综合能源系统进行中长期时间尺度规划运行研究。首先，提出了含碳流的多能耦合综合能源系统框架，研究了氢气作为季节性储能的制-储-用过程，建立混合氢
燃料的燃气机组作为氢-电耦合的模型；其次，建立以系统改造升级成本、运行成本以及惩罚成本综合最小为目标函数的综合能源系统双层规划-运行模型，利用改进的差分进化算法进行求解，达到配置最优设备容量的同时优化不同机组运行策略；最后通过具体算例，证明了季节性氢储的引入可有效促进可再生能源的消纳，平抑净负荷曲线季节性峰谷差，提高综合能源系统的经济性，减少系统的碳排量。
[0125]
本发明在传统综合能源系统上增加了可用于中长期储能的电解水制氢和储氢模型，对传统燃气轮机模型进行改造升级，使其可高效利用氢气作为燃料进行发电出力，促进氢-电紧密耦合，充分考虑了碳捕集-封存和氢气甲烷化反应的利用过程，并且将电解水反应和氢气甲烷化反应过程余热再利用也增加耦合到热能网中。
[0126]
一、含碳流的综合能源系统框架及季节性氢储关键设备模型
[0127]
为促进可再生能源的消纳利用，平抑系统净负荷季节性波动，在传统综合能源系统上增加了可用于中长期储能的电解水制氢和储氢设备，同时对燃气轮机进行部分改造升级，使其可高效利用氢气作为燃料进行发电出力，促进氢-电紧密耦合。增加考虑了碳捕集、封存以及和氢气甲烷化反应的利用过程。考虑电解水反应和氢气甲烷化反应都是高温反应，将反应过程余热再利用也增加耦合到热能网中。冷热负荷在中长时间尺度的季节特性可以通过吸收式制冷机设备互补，燃气锅炉、电制冷设备作为备用补充。本文构建了一种能适用于中长时间尺度的含碳流的多能耦合综合能源系统框架如图1所示。
[0128]
下面介绍季节性氢储关键设备模型：
[0129]
1、电解水制氢模型
[0130]
目前电解制氢主要分为碱性液体电解(alkaline electrolysis，alk)、质子交换膜电解(polymer electrolyte membrane electrolysis，pem)、碱性固体阴离子交换膜(anion exchange membrane，aem)，以及高温固体氧化物电解(solid oxide electrolyzer cell，soec)，其中soec方法能量转化效率可达到100％，是未来高效制氢的重要途经。
[0131]
电解水装置(electrolysis device,ed)能量转换模型：
[0132][0133][0134][0135][0136]
式中，分别是电解水装置在t时刻的耗电功率、产氢容量、产热功率，为电解水装置电制氢效率和余热利用效率，为电-氢气折算单位换算系数，为氢气燃烧热值142500kj/m3，表示电解水装置最大制氢容量，α为千瓦与千焦每小时的转换系数取3600。
[0137]
2、储氢模型
[0138]
季节性氢储(储氢装置)区别于普通的储能装置，可以在日内进行多次充放，季节性氢储在每个周期内由供需关系决定，只有一个充电或者放电的状态，且考虑到大容量的储存和充放电的过程细化，类比锂电池从2个变量角度考虑，季节性储氢方式约束模型为：
[0139]
[0140][0141][0142]sshs
(0)＝0.5q
shs
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(8)
[0143][0144]
0≤s
shs
(t)≤q
shs
ꢀꢀꢀꢀꢀꢀ
(10)
[0145]
式中，分别表示在t时刻储氢装置充入(储存)和放出(释放)的功率，分别表示在t时刻充入和放出的0-1状态量，v
shs-max
表示储氢装置最大功率，q
shs
表示储氢装置最大容量，s
shs
(0)、s
shs
(t)、s
shs
(t-1)分别为储氢装置存储能量初始值、t时刻剩余能量、t-1时刻剩余能量，分别表示储氢装置充入和放出的效率，
△
t表示储氢装置充放的单位时间。
[0146]
3、氢气甲烷化
[0147]
可以协同碳捕集、利用与封存装置，利用常规机组产生的二氧化碳和氢气进行甲烷化，进一步推动双碳目标的实现。目前co2加氢合成甲烷的过程主要通过气固多相催化反应或者生物法来实现，受温度、压力、具体催化剂类型和原料气碳氢比等因素有关。co2甲烷化反应是强放热反应，同时主要产物之一是纯水，为电解水制氢的重要原料，因此过程产物和热量可参与综合能源框架里充分利用。
[0148][0149][0150]
ξ
e-gas
＝q
gas
/α
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0151][0152][0153]
式中，表示在t时刻氢气甲烷化容量以及氢气甲烷化反应过程余热的功率，分别表示在t时刻氢气甲烷化消耗的氢气量和二氧化碳量，ω1表示反应过程中二氧化碳气体混入比例，表示氢气制甲烷的甲烷和热能转换效率，为氢气甲烷化装置的最大容量，ξ
e-gas
为电-天然气折算单位换算系数，q
gas
为天然气燃烧热值33486.8kj/m3。
[0154]
二、综合能源系统双层规划-运行模型
[0155]
1、上层规划投资成本目标函数
[0156]
主要考虑在现有能源系统基础上进一步发展的规划，因此常规风电机组、火电机组、水电机组、燃气机组、以及冷热设备的投资成本暂忽略不计，只考虑燃气机组使用混合燃料的改造升级成本、氢气制-储-用过程设备成本投资，碳封存的成本受捕集的二氧化碳容量影响，因此合并在后面运行成本当中。
[0157][0158]
λ
crf
＝r
·
(1 r)y/((1 r)
y-1)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0159]
式中，c
inv
表示上层规划投资成本，分别为混合燃料的燃气机组改造升级、电解水装置、储氢装置功率、储氢装置容量、氢气甲烷化装置的单位容量投资成本；p
gt-max
、v
shs-max
、q
shs
、表示上层规划模型需求得最优的系统燃气机组改造升级最大容量、电解水装置最大制氢容量、储氢装置最大功率、储氢装置最大容量、氢气甲烷化装置的最大容量，提供给下层模型优化约束条件；λ
crf
为资金回收系数，r为年利率，取4％，y为系统设计平均寿命期限，取20年。
[0160]
2、下层运行优化成本目标函数
[0161]
运行成本c
op
主要分为燃料成本c
opf
、其他装置运行成本c
opa
和机组启动成本c
opu
，其中火电机组的燃料成本进行分段线性化处理：
[0162]cop
＝c
opf
c
opa
c
opu
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0163][0164][0165][0166][0167]
式中：c
gas
为天然气价格成本，为在t时刻系统购买的天然气量；a1、a2、a3、b1、b2、b3为火电机组发电成本分段线性化的系数，cg′
eni
(t)为在t时刻火电机组i的燃料成本变量；p
gen,i
(t)为火电机组i在t时刻的出力功率；为电解水、储氢、氢气甲烷化、碳捕集-封存、热制冷(吸收式制冷机)等装置的单位容量运行成本；为火电机组、混和燃料燃气机组、传统燃气机组的启动成本，为在t时刻火电机组、混和燃料燃气机组、传统燃气机组的启动状态变量，n
gen
为火电机组的数量，为在t时刻碳捕集后封存的二氧化碳气体量、吸收式制冷机发出冷功率。
[0168]
为了促进新能源的消纳，引入惩罚成本c
pw
，主要为弃风光水的惩罚：
[0169][0170]
式中，c
cut
为弃风光水惩罚单位成本，为在t时刻弃风、弃光、弃水功率。
[0171]
因此本节提出的考虑季节性氢储的综合能源系统双层规划-运行模型总的目标函数为：
[0172]
minc
total
＝c
inv
c
op
c
pw
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(24)
[0173]
式中，c
total
为总成本。
[0174]
3、约束条件
[0175]
1、系统电、热、冷、氢功率平衡方程
[0176][0177][0178][0179][0180][0181]
式中，le(t)、lh(t)、l
co
(t)、分别表示电、热、冷、氢常规负荷在t时刻的需求量。
[0182]
式(25)表示综合能源系统的电功率平衡，为火电机组i在t时刻净出力，p
pv
(t)、p
wt
(t)、p
hp
(t)为风电、光伏、水电机组在t时刻出力，为传统燃气机组n在t时刻净出力，为改造后混合燃料燃气机组在t时刻净出力，为电解水装置、电制冷设备在t时刻消耗的电功率。
[0183]
式(26)表示综合能源系统的天然气量平衡，为传统燃气机组n在t时刻燃烧的天然气量，为改造后混合燃料燃气机组t时刻燃烧的天然气量，为燃气锅炉t时刻燃烧的天然气量，为氢气甲烷化装置在t时刻制得的天然气量。
[0184]
式(27)表示综合能源系统的热功率供需平衡，为传统燃气机组n在t时刻燃烧余热功率的收集转换功率，为改造后混合燃料燃气机组在t时刻燃烧余热功率的收集转换功率，为在t时刻氢气甲烷化和电解水反应过程中热能的收集利用功率，为燃气锅炉在t时刻补充热功率，为在t时刻吸收式制冷机设备吸收的热功率，为在t时刻系统热能无法有效利用而弃掉的功率。
[0185]
式(28)表示综合能源系统的冷功率供需平衡，为在t时刻电制冷设备以及吸收式制冷机设备的制冷功率。
[0186]
式(29)表示综合能源系统的制-储-用氢量的平衡，表示在t时刻电解水反应产生的氢气量，表示在t时刻氢气甲烷化和改造后混合燃料燃气机组消耗的氢气量，表示在t时刻季节性氢储释放和储存的氢气量。
[0187]
2、可再生能源出力约束
[0188][0189][0190]
[0191]
为在t时刻风电、光伏、水电预测最大出力值。
[0192]
3、火电机组出力约束
[0193]
由于本节研究的是综合能源系统中长时间尺度的运行优化，以一周为具体的优化时间颗粒，因此机组具体出力变动的爬坡速率等约束忽略，为保证短时间尺度运行有充分的灵活资源调度，因此系统开机的火电机组至少有2台以上。
[0194]ugen,i
(t)p
gen-min,i
≤p
gen,i
(t)≤u
gen,i
(t)p
gen-max,i
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(33)
[0195][0196][0197][0198]
式中，p
gen-maxi
、p
gen-mini
分别表示火电机组i出力的最大值、最小值，u
geni
(t)表示在t时刻火电机组i的状态变量，为在t时刻火电机组i启动状态变量，为在t时刻火电机组i停机状态变量。
[0199]
4、燃气机组出力约束
[0200]
本节设置2种类型燃气机组，一种是传统只以天然气为原料，一种是天然气和氢气混合为原料，以混合燃料的燃气机组模型为例表示为：
[0201][0202][0203][0204][0205][0206][0207]
式中，分别表示混合燃料燃气机组在t时刻输出的电功率和热功率，分别表示混合燃料燃气机组在t时刻天然气、氢气输入体积流量，ξ
e-gas
、分别表示混合燃料燃气机组电-天然气、电-氢气折算单位换算系数，ω2表示混合燃料燃气机组氢气和天然气的混合系数；分别表示混合燃料燃气机组转换电能和热能的效率；u
gtc
(t)、为0-1的状态量，分别表示在t时刻混合燃料燃气机组的状态变量、启动状态变量、停机状态变量，p
gtc-max
、p
gtc-min
为混合燃料燃气机组的最大、最小功率。
[0208]
5、碳捕集-封存装置约束
[0209]
碳捕集-封存装置通常直接安装在火电机组和燃气机组附近，便于就地消纳碳排。
[0210]
以燃气机组为例，模型公式为：
[0211]
[0212][0213][0214][0215]
式中，为燃气机组在t时刻产生的电力功率，为燃气机组在t时刻净输出电功率，为碳捕集-封存装置在t时刻运行能耗，为碳捕集-封存装置在t时刻基础固定能耗，表示碳捕集-封存装置在t时刻捕集二氧化碳气体量，表示碳捕集-封存装置在t时刻捕集后封存的二氧化碳气体量，表示碳捕集-封存装置在t时刻捕集后加入甲烷化反应的二氧化碳气体量或者其他的碳利用量，为捕集单位碳的运行能耗系数，为捕集效率，为燃气机组单位能量碳排放强度。
[0216]
6、冷热能量补充机组约束
[0217]
吸收式制冷机模型：
[0218][0219][0220]
式中，表示吸收式制冷机在t时刻发出冷功率，表示吸收式制冷机在t时刻吸收的热功率，表示吸收式制冷机能量转换效率，表示吸收式制冷机能转化发出最大冷功率，u
ac
(t)为表示吸收式制冷机在t时刻启停状态的0-1变量。
[0221]
电制冷设备：
[0222][0223][0224]
式中，表示电制冷设备在t时刻发出冷功率，表示电制冷设备在t时刻消耗的电功率，表示电制冷设备能量转换效率，表示电制冷设备转化发出最大冷功率，u
fr
(t)为表示电制冷设备在t时刻启停状态的0-1变量。
[0225]
燃气锅炉：
[0226][0227][0228]
式中，表示燃气锅炉在t时刻发出热功率，表示燃气锅炉在t时刻消耗的电功率，表示燃气锅炉能量转换效率，表示燃气锅炉转化发出最大热功率，u
gb
(t)为表示燃气锅炉在t时刻启停状态的0-1变量。
[0229]
三、改进差分进化算法(modified differential evolution，mde)求解
[0230]
差分进化算法(differential evolution，de)是一种基于群体的启发式搜索算法，与遗传算法类似都包括变异、杂交和选择操作，为了提高寻优效率，现已有较多研究根据算法中变异过程原理改进变异算子f，本文基本流程如图2所示。
[0231]
1)初始化：确定规划容量的边界范围，设置种群数量n
p
，取50，随机生成初始种群，
种群个体：
[0232][0233]
2)变异：
[0234][0235]
f＝2
λ
f0ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(55)
[0236][0237]
式中，为第g代中随机选取的3个不同个体，为变异种群中的个体，f0为初始设定的变异参数，取0.5；g表示当前进化代数，gm表示最大进化代数，取80。
[0238]
3)交叉：
[0239][0240]
式中，为交叉后所得种群中第i个体的第n维变量；cr为交叉因子，取值[0,1]之间。
[0241]
4)竞争：
[0242][0243]
式中，为对应个体的适应度函数，即模型中综合成本目标函数，采用最小化问题的选择方式。
[0244]
5)当g》gm时，算法终止，得到最优解，否则，g＝g 1，流程返回2)，进行下一次优化。
[0245]
优选地，采用gurobi求解器与改进差分进化算法配合求解。
[0246]
计算示例：
[0247]
本发明用于送端地区负荷量较大的综合能源系统年度规划运行，通常这种地区已经有部分传统综合能源系统的能源设备，以及风电、光伏、水电、火电机组的规划建设，本发明在这些现有能源设备的基础上，优化新增季节性储能以及二氧化碳的捕集利用等设备的容量配置，并以周为优化周期分析了地区年度系统各设备出力情况，具有实际参与电网规划、运行参考的意义，具体实施方式为：
[0248]
1、根据前5年历史数据以及未来天气情况预测，确定规划年度内风电、光伏、水电等可再生能源的出力，以及根据负荷需求增长确定多能负荷年度需求曲线。
[0249]
2、统筹地区内常规火电机组、燃气机组、燃气锅炉等常规能源设备参数，确定新增电解水装置、储氢装置、碳捕集装置等参数及购置成本等。
[0250]
3、根据前面提出的系统架构和双层模型，增加各能源设备的运行和能源转换的约束，通过改进差分优化算法，针对目标函数进行模型优化从处理，确定最优的混合燃料燃气机组改进最大出力、电解水装置最大制氢容量、储氢装置最大容量、储氢装置最大功率以及氢气甲烷化装置的最大容量。
[0251]
4、在确定规划设备容量新增的基础上，将各设备的最大最小约束系数代入模型，
进行运行优化结果分析，确定常规火电机组和燃气机组的年度出力特性，可以以此为根据安排检修计划；确定制氢和储氢设备充放的季节性特性，分析可再生能源的季节性波动平抑效果；确定燃气机组和多能耦合反应过程热出力以及冷能机组的年度曲线，对规划年的多能需求进行掌握及调度优化；分析年度碳排量的各项指标，分析氢气甲烷化反应的可行性，为未来低碳目标的实现探索多样途径。
[0252]
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。