孙轶恺,漆淘懿,张利军,胡怡霜,俞楚天
(1. 国网浙江省电力有限公司经济技术研究院,杭州310008;2. 浙江大学电气工程学院,杭州310027)
为实现“碳达峰,碳中和”的目标,推进能源改革,提高能源利用效率势在必行[1]。综合能源系统由新能源发电、储能、燃气轮机等设备组成,通过电、热、冷各种能源间的协调互补,在满足系统多能需求的基础上,能够提高能源效率,促进可再生能源的消纳[2]。基于综合能源系统的多种优势,世界上许多国家都在大力推进综合能源系统的研究和落地工程[3 - 8]。以我国为例,北京大兴机场建设了包含分布式地源热泵,光伏发电和储能等设施的综合能源系统,能够对机场运营中所需的电力、供热和制冷进行有效补充[9]。上海迪士尼园区基于能量梯级利用的冷、热、电和压缩空气四联供技术,将一次能源利用率提升到80%以上[10]。
然而,综合能源系统的发展一直被高昂的建设成本和复杂的运营优化所困扰[11 - 12]。随着通信技术的飞速发展,综合能源系统对内部负荷资源的控制能力不断增强[13]。通过对灵活负荷更加精准地进行调控,综合能源系统不仅能够响应现货市场中实时电价,还能利用灵活负荷的需求响应能力参与辅助服务市场。山东、江苏、浙江等地针对负荷侧资源参与辅助服务市场已出台相关政策,制定了负荷侧主体参与辅助服务市场的准入条件、结算方法、考核机制等规则[14 - 16]。因此,综合能源系统通过提供调频备用(regulation reserve, RR)、同步备用(synchronized reserve, SR)等辅助服务,能够获取市场补贴从而进一步降低系统运行成本。
由于燃气轮机、蓄电池等传统设备的响应潜力有限,亟需挖掘新的灵活资源为综合能源系统提供响应容量。冰蓄冷空调包括制冷机、制冰机和蓄冰槽等组件,可以在夜间低电价时段制冰并储存在蓄冰槽内,白天高电价时段融冰制冷从而降低制冷机功率,是综合能源系统重要的制冷和储冷设备[17]。文献[18]对深圳地区的电冰蓄冷空调的典型工程应用进行了经济型研究,结果证明了冰蓄冷空调的削峰填谷作用有利于电网的健康运行和有序发展。文献[19]针对冷热电联供型微网提出了含冰蓄冷空调的多时间尺度优化模型,有效平抑了日内需求侧的随机波动,实现微网的经济稳定运行。文献[20]在含冰蓄冷空调的工厂综合能源系统中,利用冰蓄冷空调夜间蓄冰,白天融冰制冷,降低工厂运营成本。
然而,以上研究只关注了冰蓄冷空调的蓄冰时段优化,忽视了冰蓄冷空调在白天的负荷调节能力。白天冰蓄冷空调的制冷量来源于制冷机和融冰两部分,可以根据需求优化其各自的制冷量。在尖峰电价时段,冰蓄冷空调可以增加融冰制冷量,从而降低制冷机的电功率以减少用电成本。此外,不同负荷率下制冷机的能效比不同,优化过程中能效比也是必须要考虑的因素。在高电价时期,大幅降低制冷机功率而增加融冰制冷功率,不仅会使制冷机能效比显著下降,还会加速蓄冰量消耗。因此,需要综合考虑冰蓄冷空调的能效比、蓄冰量、电价等因素,实现其经济最优运行。
因此,本文充分挖掘了综合能源系统内冰蓄冷空调、蓄电池、燃气轮机等可调节资源的需求响应潜力,并在现货市场和辅助服务市场中,优化了综合能源系统的多能需求和响应容量,从而降低系统的运营成本。首先,本文针对综合能源系统,提出了考虑需求响应的设备通用运行模型。随后,针对综合能源系统中的冰蓄冷空调、燃气轮机、蓄电池,分别建立了不同的精细化需求响应模型,并拟合了冰蓄冷空调在不同负载率下的能效比曲线。进而以经济最优为目标,基于设备的需求响应模型,采用大M法对优化系统运行中的购电量和需求响应容量。最后通过算例证明了所提优化方法在市场环境下能够有效提高综合能源系统经济性。
图1展示了含冰蓄冷空调综合能源系统的结构和能量转移过程,主要包括能量输入、能量生产、能量转换、能量储存和能量消耗5个环节。综合能源系统内能量的主要形式为电、热和冷,因此系统内包含电流,热流和冷流3种能量流动。
图1 含冰蓄冷空调的综合能源系统结构Fig.1 Integrated energy system structure including ice-storage air-conditioning
综合能源系统中能量的来源主要有两种形式,分布式发电、燃气轮机和余热锅炉等系统自身设备的生产能量和外部购买的电力、天然气等能量。为满足系统实时变化的各种能量需求,不同类型的能量可通过能量转换设备彼此转换。同时,能量储存设备可以为综合能源系统储存一定的能量,提高系统的灵活性。以冰蓄冷空调为例,其不仅可以通过制冷机将电量转换为冷量,还可以利用自身的蓄冰槽储存冷量。根据需求优化调整能量生产、能量转换和能量储存3个环节,以满足不同类型负荷的用能需求。
综合能源系统中设备具有灵活的运行方式和能量储存能力,因此其可以调整内部设备的运行状态,改变其对外界能源输入需求,从而为电力系统提供辅助服务。因此,针对综合能源系统中的冰蓄冷空调、蓄电池等设备,需要评估其能够提供的需求响应能量,通过对设备的优化调节以提供辅助服务获取补贴,从而降低综合能源系统的运行成本。
综合能源系统可以提供的辅助服务主要有调频备用和同步备用,在提供需求响应服务时,响应资源需要符合系统运营商设定的响应标准。因此,首先需要对设备的需求响应潜力进行评估,主要包括备用容量、响应时间和持续时间。美国PJM辅助市场中对需求侧资源最小备用容量、响应时间和持续时间有要求[21]:需求侧资源提供调频备用时,最小备用容量为0.1 MW,在5 min内响应并持续1 h;提供同步备用时,最小备用容量0.1 MW,在15 min内响应并持续1 h。最小备用容量是单个参与者需要提供的最小响应容量;响应时间是指从接收需求响应信号到满足需求的时间间隔;持续时间是指需求响应资源可以保持备用容量或跟随需求响应信号的持续时间。尽管现阶段我国包含负荷参与的辅助服务市场建设还不够完善,但是随着电力市场改革的进一步深化,具有负荷调节能力的综合能源系统会逐渐参与到市场中来,响应电力系统的调度指令或者市场价格信号,为系统提供调峰、备用等辅助服务。
设备在实时运行过程中,只有满足响应要求时,才能提供辅助服务。为评估各设备的需求响应潜力,提出了考虑需求响应的设备运行模型,如图2所示。当设备提供调频备用时,主要用于平抑可再生能源的出力波动,它需要在持续时间内提供足够的上下容量,如图2(a)所示。当设备提供同步备用时,需要提供下调容量以解决需要削减负荷的能源短缺问题,如图2(b)所示。此外,设备的工作功率可能会在调频备用或同步备用的响应持续时间内发生变化。因此,需要在整个响应期间确保有足够的调频备用或同步备用容量。
图2 设备提供辅助服务时的运行模型Fig.2 Operation model when equipment provides auxiliary service
综合能源系统中包含多种设备,不同设备具有显著的响应特性差异。为准确描述设备的运行状态及需求响应潜力,本节针对不同设备,建立了考虑需求响应的精细化模型。
冰蓄冷空调主要包括制冷机、制冰机和蓄冰槽3个组件,其运行可以分为制冰蓄冷和融冰制冷两个阶段。夜间气温低于白天,冷负荷需求较少,制冷机可以独立满足冷负荷需求。同时,夜间的电价较低,在夜间制冰蓄冷然后在白天融冰制冷,能够降低运行成本。
2.1.1 蓄冰过程
冰蓄冷空调在夜间蓄冰过程中,蓄冰槽不能同时融冰制冷,因此融冰制冷功率为0,由制冷机满足所有的系统冷功率需求,如式(1)所示。
(1)
蓄冰过程中,冰蓄冷空调的功率为制冷机和制冰机的功率之和,制冷电功率、制冰电功率与冷功率之间的关系可以表示为:
(2)
2.1.2 融冰过程
白天冰蓄冷空调的制冷量由制冷机和蓄冰槽融冰共同提供。根据实时的供冷需求和冰蓄冷空调的运行状态,在可调区间内,制冷功率区间可以灵活分配融冰制冷量和制冷机制冷量,如式(3)所示。
(3)
(4)
2.1.3 冰蓄冷空调能效比
能效比是描述制冷机或者制冰机效率的指标,能效比越高,单位电量就能转换为更多的冷量。冰蓄冷空调制冰机的工作功率是恒定的,其制冰能效比可看作常量。而制冷机的功率会根据冷量需求实时调节,其能效比ηEER和负载水平密切相关,不同负载下冰蓄冷空调的效率具有显著差异。本文对冰蓄冷空调进行了部分负荷性能曲线拟合,结果如下:
(5)
式中:a、b、c和d分别为冰蓄冷空调能效比的拟合系数;RPL为负载比;Preal和Prated分别为制冷机的实际功率和额定功率。
从图3冰蓄冷空调能效比与负载率的拟合结果可以看出,在负载率小于80%的时候,制冷机的能效比与负载率大致呈线性关系,随着负载率的增加,能效比也随之提高。制冷机的最高能效比能够达到5.5,但是最低能效比只有2.23。因此,在冰蓄冷空调制冷机制冷量和融冰制冷量的优化分配过程中,能效比是不可忽视的一个因素,应尽量保证冰蓄冷空调运行在高效区间。
图3 冰蓄冷空调能效比与负载率之间的关系Fig.3 Relationship between EER and PLR of ice-storage air-conditioner
2.1.4 备用约束
(6)
(7)
此外,在提供调频备用或者同步备用过程中,从响应开始到响应结束期间,冰蓄冷空调提供的调频备用容量和同步备用容量都必须保持不变,如式(8)所示。
(8)
式中DSR为提供辅助服务的持续时间。
燃气轮机是综合能源系统中的典型能源生产设备,大多数燃气轮机都配备有余热锅炉。燃气轮机以天然气作为能源,输出电能和热能。电能和热能的生产可以表示为:
(9)
燃气轮机参与需求响应时的约束和冰蓄冷空调类似,可参照式(6)—(8)得到,此处不再赘述。
燃气锅炉是综合能源系统中重要的能源转换设备,燃气锅炉以天然气为能量来源,产生热蒸汽从而满足综合能源系统的热负荷需求,其产热功率可以表示为:
(10)
燃气锅炉既不消耗也不产生电能,不能提供辅助服务。通过优化其运行时间,可以补充其他设备由于提供辅助服务而减少的热量供应。
吸收式制冷机将热量转化为冷量,可以有效利用能源生产过程中的余热,补充综合能源系统的冷负荷需求。吸收式制冷机的制冷功率和吸热量成正比,如式(11)所示[22]。
(11)
吸收式制冷机既可以作为综合能源系统供冷的补充,也可以为电力系统提供辅助服务,其参与需求响应的约束可以参照式(6)—(8)得到。
能量存储设备可以存储能量,并在需要时释放能量,打破能量供需同时性的限制。以蓄电池为例,蓄电池的放电和充电过程可以表示为:
(12)
在蓄电池充放电过程中,充放电速率不能超过最大值,其约束可表示为:
(13)
为保证在提供需求响应服务的持续时间内都能够提供足够的同步备用和调频备用容量,避免在响应未结束时电池就无法提供响应容量,因此电池的实时电量不能超过最大和最小电量,其约束可表示为:
(14)
(15)
站在2018年的岁尾回望一路走来的历程,我相信,自豪感会充斥在每一个人的心中。而面对未来,我们的任务将更加艰巨,工业互联网、人工智能和智慧工厂等全新的词汇出现在大众视野中,学习创新是必不可少的,这一代人肩负着由大到强的光荣使命,站在新的转折点上,很多企业都积极进行技术改革创新,抓住时代机遇。
基于以上考虑需求响应的设备模型,本节提出综合能源系统的最优运行模型。然后,对部分非线性约束进行线性化处理,以降低计算难度。最后,通过求解线性化后形成的整数线性规划问题,获得最优运行结果。
(16)
(17)
综合能源系统的运行除了受电力系统约束,供热系统约束,制冷系统约束,系统中各组件还需要满足文章第2节“综合能源系统设备模型”中所述的运行约束。
如式(18)所示,综合能源系统产生的电力和购买的电力之和应始终等于电力消耗。同时,综合能源系统中通常同时存在交流系统和直流系统。式(18)为交流系统约束,式(19)为直流系统约束。式(20)表示这两个系统之间的转换。
(18)
(19)
(20)
系统的热量应满足热水、蒸汽、空间加热以及吸收式制冷机的需求。冷却能量主要满足部分设备的冷却需求和空间冷却。因此,供热和冷却系统的约束可以分别表示为:
(21)
(22)
由于设备参与需求响应的过程中存在非线性约束,求解难度较大,因此,本文采用大M法对约束进行线性化处理。大M法能够将包含0-1变量非线性约束线性化,通过引入人为设定的M值,在0-1变量的取值为0或1时,分别松弛掉两个互补条件中不起作用的一个条件[23]。以式(6)为例,使用大M法可以将其转换为式(23)。
(23)
式中M为一个非常大的正整数。其余各式的非线性约束,可以通过类似的方法转换成线性约束,此处不再赘述。
式(20)为分段函数,可以用以下线性表达式替换。
(24)
基于第3节给出的设备需求响应模型和求解方法,本文以中国广州某工厂为研究对象,建立了含冰蓄冷空调的综合能源系统模型,并利用美国PJM市场中的实际数据进行了仿真分析[24]。综合能源系统中各设备的相关参数如表1所示,每日的交流负荷和直流负荷如图4所示,每日的冷负荷和热负荷的需求如图5所示。
表1 综合能源系统设备参数Tab.1 Equipment parameters of integrated energy system
图4 综合能源系统的日交直流负荷曲线Fig.4 AC/DC load curves of integrated energy system
图5 综合能源系统的日冷热负荷曲线Fig.5 Daily cooling and heating load curves of integrated energy system
表2中列出了综合能源系统内部可以提供不同类型需求响应的设备,其中Y和N分别表示设备可以提供或不能提供相应的辅助服务。
表2 综合能源系统设备及需求响应能力Tab.2 Components of the IES and their DR ability
图6展示了2020年6月2日PJM数据中的电价,调频备用的价格和同步备用的价格[25]。并将天然气价格设定为5美元/千立方英尺,蒸汽价格设置为5美元/千立方英尺(1立方英尺约为0.028 317 m3)。仿真的优化步长设置为15 min。
图6 实时电价、调频备用价格和同步备用价格Fig.6 Real time electricity prices, RR prices and SR prices
图7为综合能源系统内部的交流侧电量、直流侧电量、冷量和热量的生产及消耗情况。从图7(a)和(b)中可以看出,在49—80时段内,电价相对较高。为降低综合能源系统的购电费用,交流侧的燃气轮机开始工作,同时直流侧的电池开始放电以满足系统用电需求。从图7(c)中可以看出,冰蓄冷空调也降低了用电功率,减少的制冷量由融冰制冷和吸收式制冷机共同补充以满足系统需求。图7(d)中余热锅炉随着燃气轮机的启动也开始产生热量,因此系统的购热量也随之降低。
图7 综合能源系统优化运行结果Fig.7 Optimized operation results of integrated energy system
图8为冰蓄冷空调的运行情况,从图中可以看出,在1—40时段,此时电价最低,为冰蓄冷空调的蓄冰过程,因此融冰制冷量小于0,此时冰蓄冷空调的制冷量全部由制冷机提供。在40—84时段,冷负荷需求较高,同时电价也高于夜间,因此蓄冰槽开始融冰制冷,和制冷机联合满足冷负荷需求。在考虑不同负载率下制冷机的能效比差异显著的情况下,制冷机一直维持在相对稳定的功率,因此能效比也一直维持在较高水平。
图8 冰蓄冷空调运行情况Fig.8 Operation status of ice-storage air-conditioning
图9为冰蓄冷空调和电池的辅助服务提供情况。冰蓄冷空调可以同时提供调频备用和同步备用,空调会根据实时的调频备用和同步备用价格调整制冷机制冷量和融冰制冷量。在50—56时段,调频备用价格较高,因此冰蓄冷空调降低电功率,增加融冰量满足系统需求,从而增加调频备用。例如在64—67时段,调频备用价格大幅降低,因此制冷机增加功率减少调频备用,减少融冰量以用在随后需要降低制冷机功率的时期。
图9 冰蓄冷空调和电池的辅助服务提供情况Fig.9 Ancillary services provided by ice storage air-conditioning and battery system
表3为冰蓄冷空调优化前后的运行费用对比。从表中可以看出,冰蓄冷空调运行优化后,购电费用从原本的278美元上升至282美元,增加了4美元。这是由于在提供辅助服务时,冰蓄冷空调需要实时调整电制冷量和融冰制冷量,导致效率稍有降低。但是通过提供调频备用和同步备用服务,冰蓄冷空调分别获得了115美元和24美元的收益,使得最终运行成本从278美元下降为153美元,大幅降低了冰蓄冷空调的运行成本。
表3 冰蓄冷空调优化前后运行费用Tab.3 Operation cost of ice-storage air-conditioner before and after optimization美元
电池相较于冰蓄冷空调,主要提供调频备用。在11—12时段,由于调频备用价格突然增加,电池的充电功率迅速从1.14 MW下降到0,将调频备用从270 kW迅速提升到了1 500 kW。在65—69时段,调频备用价格大幅下降,电池也随之将调频备用功率降为0。但是在提供辅助服务期间,电池不会完全充满电或者用完电,需要一直为调频备用保留一定的余量。
参与辅助服务可以为电力系统提供更多可调节的容量,并减少综合能源系统的运营成本。表4为综合能源系统优化前后的运行费用对比,综合能源系统运行成本从3 762美元降低至3 160美元,降低了16%。和冰蓄冷空调类似,优化前后系统的购电和购热费用几乎没有差别。系统运行总费用的降低主要来自于提供同步备用和调频备用获得的收益。
表4 综合能源系统优化前后运行费用Tab.4 Operation cost of IES before and after optimization 美元
综合能源系统的最优运行在消纳可再生能源和提高能源效率方面发挥着重要作用。为降低综合能源系统运营成本并促进其发展,本文充分发掘了冰蓄冷空调的调度潜力,通过优化冰蓄冷空调的电制冷量和融冰制冷量,提高了综合能源系统的需求响应容量。仿真结果表明,市场环境下综合能源系统能够根据实时的电价、调频备用价格和同步备用价格,调度系统内各设备以实现系统最优运行。经过优化,综合能源系统中冰蓄冷空调的运行成本降低了45%,系统运行成本降幅也达到了16%,具有显著的经济效益。