【精彩论文】考虑天然气季节性存储的综合能源系统年度运行方式研究
考虑天然气季节性存储的综合能源系统年度运行方式研究
刘赫川1, 周孝信1, 杨小煜1, 李亚楼1, 李雄2
(1. 中国电力科学研究院有限公司 国家电网仿真中心,北京 100192; 2. 华北电力大学 电气与电子工程学院,北京 102206)
引文信息
刘赫川, 周孝信, 杨小煜, 等. 考虑天然气季节性存储的综合能源系统年度运行方式研究[J]. 中国电力, 2022, 55(4): 145-155.
LIU Hechuan, ZHOU Xiaoxin, YANG Xiaoyu, et al. Annual operation mode study of integrated energy system considering seasonal natural gas storage[J]. Electric Power, 2022, 55(4): 145-155.
引言
本文采用一种多时间尺度分层优化的综合能源系统生产模拟方法,研究电气耦合系统的年度优化运行,通过结合短期与长期的综合能源系统生产模拟,将冷热电气负荷、天然气大规模存储等季节性因素考虑在优化之中。同时应用此方法对CCHP装机大幅增加的综合能源系统年度运行方式进行评估。
1 电气耦合的综合能源系统模型
1.1 燃气轮机
燃气轮机以天然气作为燃料,是电气耦合综合能源系统中的核心设备。在生产模拟计算中燃气轮机的模型重点关注能量转换效率、出力、经济性指标等约束条件。燃气轮机发电功率式中:热率 RHi 反映天然气发电输入输出的效率;αi 、 βi 、 γi 为拟合成二次函数形式的效率函数中的参数。
天然气的消耗量
式中:VLH 为天然气的低热值。
为便于生产模拟的计算,将综合能源系统中各种能源形式的能量单位统一转换成电力系统单位,功率单位为kW,能量单位为kW·h。将式(1)~(2)转化为
式中:
燃气轮机的输出功率上、下限约束为
式中:
机组的启停约束为
式中:
燃气轮机的爬坡速率约束条件表示为
式中:
燃气轮机的运行总费用
式中:
本文中燃气锅炉采用固定效率模型,其天然气消耗量
式中:
与燃气轮机相同,将燃气锅炉中涉及气、热的单位统一转换成电力系统单位。燃气锅炉的输出功率
式中:
燃气轮机的总运行费用
式中:
近年来冷热电联产机组受到广泛关注,通过消耗天然气满足冷热电负荷,提供了一条减少燃煤、缓解污染的途径。图1为常见的CCHP装置结构,燃气轮机发电后将余热收集,通过余热锅炉以及溴化锂吸收式制冷装置提供热和冷的输出功率。
图1 冷热电联产装置结构
Fig.1 Structure of CCHP
CCHP中的燃气轮机采用固定效率模型以简化计算,其电输出功率
式中:
CCHP的燃气轮机可回收热量
式中:
将上述公式中的单位统一转换成电力系统单位,则有
式中:
CCHP输出功率的上、下限约束为
式中:
CCHP的运行费用
式中:
电气耦合的综合能源系统中包含了两种天然气存储装置,一种是满足天然气季节性供需平衡的大型天然气存储设备,另一种是城市天然气管网中配置的一定量的中小型天然气存储设备,用于满足短时(日或周)内的燃气供需平衡。在生产模拟计算中两者可以用同样一套模型表示其中的能量平衡,其约束条件如下。
天然气存储装置的存储量为
式中:
储气装置储气量的约束条件为
式中:
式(34)和(35)分别表示储气装置储气、放气的功率限制,式(36)确保储气装置的存储与释放两种工作状态不同时出现。
式中:
储气装置的运行总费用
式中:
生产模拟计算中除了上述关于能量转换、存储设备的约束条件外,还有表示能量平衡的约束条件。式(38)表示系统中电能供需的功率平衡,其中电负荷等于燃气轮机、CCHP发电以及与外电网交互功率
式中:I、B、C、S、Le、Lc、Lh和Lg分别为燃气轮机、锅炉、CCHP、储气装置以及电、冷、热、气负荷的数量;St为时段t天然气供应量;Lei,t 、 Lni,t 、 Lci,t 、 Lgi,t 分别为电、热、冷、气负荷。
2 综合能源系统多时间尺度生产模拟方法
本文提出一种多时间尺度的综合能源系统生产模拟方法,选取周与年两个时间尺度,采用分层优化的方法处理生产模拟计算中的MILP问题,同时将天然气年度存储因素纳入综合能源系统的研究之中,流程如图2所示。
图2 多时间尺度分层优化计算流程
Fig.2 Calculation flow of multiple time scale layered optimization
首先以周为仿真时长、小时为计算步长,根据上一章节的约束条件对全年52周逐个进行生产模拟计算,得到综合能源系统能量转换设备每周、每小时的工况。在此环节不考虑系统中年度约束条件,如天然气年度总量约束等。这一层级优化中的约束条件在上一章节进行了介绍,包括燃气轮机(式(1)~(8))、燃气锅炉(式(13)~(15))、CCHP(式(19)~(28))、小型周平衡天然气存储装置(式(32)~(36))以及能量平衡约束(式(38)~(41))。优化的目标函数为整个系统的运行费用最小,由上一章节中各设备运行费用累加得到,如(42)所示,在这一层级仿真时长为168小时,即T为168。
式中:
3 算例介绍
本节搭建电气耦合的综合能源系统算例,并在其中应用多时间尺度的综合能源系统生产模拟方法,验证该方法的有效性,同时对系统中关键设备进行分析。
3.1 算例系统为验证方法有效性,本节在已有文献中电气耦合系统基础上[24]搭建了包含6节点电力系统、7节点天然气系统、3节点热系统以及3节点冷系统的综合能源系统算例,如图3所示。6节点电力系统通过4台燃气轮机(GT1、GT2、GT3、GT5)、2台CCHP机组(CCHP2、CCHP5)与天然气系统相连接,电力系统中包含3处电负荷(EL2、EL3、EL6)以及电制热(EH3)、电制冷负荷(EC6)。天然气系统中包含2处气源(GS1、GS2)、3处气负荷(GL1、GL2、GL3),电系统通过EB1母线与外部电力系统相连起到后备作用。气系统中还配有两套储能设备,分别为城市中用于满足供气短时平衡的小型天然气存储装置GTS2和大型天然气季节性存储装置GST1。供热系统中包含3处热负荷(HL1、HL2、HL3),供冷系统中有3处冷负荷(CL1、CL2、CL3)。
本文利用分层优化方法对天然气季节性存储进行年度生产模拟计算,在算例中对于年度生产模拟计算影响不大的模型进行了简化。电、气、热、冷系统中的网络约束仅考虑能量传输容量约束,不涉及电压、气压等约束条件。对上述模型的简化主要是基于仿真场景以及计算难度上的考虑。首先生产模拟问题研究在仿真时段内的运行优化问题,仿真时间步长一般为时段(周、小时、分钟),而电压、气压等约束条件多应用于系统安全性、稳定性分析,与生产模拟研究多能互补、能量分配的内容不同。另一方面,电压、气压等潮流计算方面的约束条件将引入大量微分方程,给生产模拟中应用商业求解器求解MILP问题带来线性化等额外困难,降低计算效率;且实际生产中天然气、热力等流体管网系统压强允许波动范围很大,基本不会触及安全约束,故在本文的生产模拟中可以对上述网络约束进行简化。
算例关于电气网络的相关参数详见文献[23],冷热网络中设备的相关参数详见文献[24],算例中的单位统一转换成电力系统常用单位,功率单位采用MW,能量单位采用MW·h。仿真采用便携式计算机,处理器为intel core i7,主频2.5 GHz,优化计算调用商业求解器。
3.2 对比方案描述本文在上述耦合算例系统中对所提算法进行验证。首先将所提多时间尺度分层优化方法与典型周(不计年度约束条件)累加方法进行对比,验证所提算法的正确性和必要性。其次应用所提算法对天然气季节性存储装置存储量占比不同的几个场景进行生产模拟计算,对季节性储气装置的运行状态进行分析。最后对系统加装CCHP前后的年度运行方式进行对比分析,研究CCHP对储气装置运行方式带来的影响。4 结果与分析
4.1 典型周累加方法
采用典型周的方法在不考虑天然气相关年度约束条件(季节性存储、总量约束)的情况下对上述综合能源系统进行生产模拟计算。文中算例选取春夏秋冬4个典型周,其冷、热、电、气负荷如图4中a)~d)所示。
图4 四季典型周负荷功率
Fig.4 Load power in four typical weeks of four seasons
根据上述数据进行生产模拟计算,其年度总用气量折算为3.07×107MW·h。此处计算得到的天然气消耗量没有考虑天然气的年度供应量上限的约束条件,为理想化结果,在实际系统中可能出现天然气供应量不足的情况。
4.2 多时间尺度分层优化方法
采用本文所提方法的算例,考虑天然气年度存储及相关约束。算例中参考典型周方法中4个典型周数据,确保生成的全年52周数据中负荷的周平均值与4.1中的负荷平均值相同。本节算例中设置天然气年度供应上限为3.00×107MW·h,季节性天然气存储装置的存储量参设定为年度用气量的8%、10%、13%三个等级,观察季节性天然气存储装置对综合能源系统的年度运行方式影响。上述各等级依据中国的天然气存储现状进行选取,中国现阶段天然气季节性可调峰储气量与全年用气总量相比占比不高,仅为6%左右,与世界平均水平的差距颇大(12%~14%)。根据现有规划,预计2025年储气库储量达到8%,而未来要求的全年的储气规模占比要达到15%以上才能保障中国能源安全。4.2.1 与典型周累加方法对比
与4.1中不计及天然气季节性存储的典型周方法生产模拟进行比较,多时间尺度分层优化的生产模拟方法将天然气系统年度的约束条件考虑进去。以储气量站年度消费总量13%的场景为例,进行生产模拟计算,得出的全年燃气消耗量折算后由3.07×107MW·h降至2.90×107MW·h,小于设定的3.00×107MW·h的年供应上限。与典型周累加方法相比耗气量有所不同的原因是在第二层优化计算中,上一层计算得到的天然气消耗量不满足年度总量约束条件,因此返回上一层对天然气消耗进行了修正,减小了燃气发电、供热的年度消耗量。对于因天然气总量限制而产生的供热、供电不足,本算例由外接电网购电、电供暖满足。在冬季取暖负荷较大时,通过外电网购电,利用电制热装置满足一部分热负荷,降低燃气锅炉分担的热负荷,减少天然气消耗。4.2.2 天然气储气量占比不同场景对比
图5a)、b)、c)分别为季节性天然气存储量占年度总用气量的8%、10%、13%三种情况下储气装置一年52周的运行方式。由于天然气源年度总量以及周供应量的约束,季节性存储装置的存储量越大,其储满所需时间越长,因此全年52周中的工作时间越长。
图5 季节性天然气存储装置年度运行方式
Fig.5 Annual operation mode of seasonal natural gas storage device
以季节性天然气存储装置存储量占比为13%的情况为例,图6中a)、b)分别展示了天然气用气高峰周、用气低谷周的168小时季节性存储装置的工作情况。每一周的存储、释放天然气的总量由第二层的优化计算得出,经过168小时完成该周天然气的消耗或存储目标。
图6 季节性天然气存储装置用气高峰、低谷周工作方式
Fig.6 Operation mode of seasonal natural gas storage device in weeks of gas consumption peaks and troughs
4.2.3 CCHP投运前后场景对比
在天然气季节性存储量容量固定为年度总用气量13%的情况下配置一定量的CCHP,研究其投入运行后对电气耦合系统年度运行方式的影响。算例中的CCHP装机容量选取实际工程应用中的常用数据,由2台78 MW燃气轮机和溴化锂制冷设备组成CCHP[17]。图7为初始无CCHP算例与配置CCHP后的天然气存储装置运行状态对比。
图7 天然气存储装置年度运行方式对比
Fig.7 Annual operation mode comparison of natural gas storage device
结果表明:配置CCHP后,天然气存储装置的燃气存储量峰值小于没有配置CCHP的情况,同时在存储装置全年的工作(存储、释放)时间也少于初始状态。造成这种结果的主要原因是CCHP的运行在年度上减少了天然气的消耗量。这部分减少的原因是供冷季CCHP以天然气为原料通过热泵制冷的效率要高于燃气发电电制冷的效率,在相同制冷负荷的情况下所消耗的天然气比较少。
5 结论
本文首先介绍了因天然气消耗增加而形成了电气耦合逐渐加深的综合能源系统,并说明了综合能源系统年度仿真计及季节性天然气存储装置的必要性。进而提出了一种多层级多时间尺度的综合能源系统生产模拟方法,考虑有关负荷波动、天然气季节性存储装置等年度约束。通过搭建冷热电气耦合系统算例将传统的典型周生产模拟方法与多时间尺度分层优化的生产模拟方法进行对比。结果表明,上述分层优化的方法能有效评估天然气资源总量、大型天然气存储装置等年度约束条件对综合能源系统年度运行方式的影响。算例中典型周方法的天然气消耗量大于年度天然气总量,分层优化方法中则通过对燃气消耗进行修正,将天然气年度消耗降至天然气总量限制之下,证明了该方法的必要性。此外,算例对CCHP装机日益增加的情况进行分析,评估了CCHP装机对电气耦合系统的年度运行方式影响。这种影响通过大型季节性天然气存储装置全年工况的先后变化对比来体现,天然气年度消耗的减少导致季节性存储装置的存储量峰值降低以及全年中工作时间的减少。
(责任编辑 许晓艳)作者介绍
刘赫川(1990—),男,通信作者,博士研究生,从事电力系统仿真、综合能源系统建模与优化等研究,E-mail:liuhechuan90@163.com;★
周孝信(1940—),男,中国科学院院士,博士生导师,从事电力系统分析、稳定和控制等研究;
★
杨小煜(1971—),男,博士,高级工程师,从事电力系统数字仿真、综合能源分析计算研究;
★
李亚楼(1974—),男,博士生导师,高级工程师(教授级),从事电力系统建模与仿真研究;
★
李雄(1992—),男,博士研究生,从事电力系统分析、高电压系统仿真建模等研究.
往期回顾
审核:方彤
根据国家版权局最新规定,纸媒、网站、微博、微信公众号转载、摘编《中国电力》编辑部的作品,转载时要包含本微信号名称、二维码等关键信息,在文首注明《中国电力》原创。个人请按本微信原文转发、分享。欢迎大家转载分享。