南方电网技术SouthernPowerSystemTechnologyISSN1674-0629,CN44-1643/TK《南方电网技术》网络首发论文题目：市场环境下含冰蓄冷空调的综合能源系统优化运行作者：孙轶恺，漆淘懿，张利军，胡怡霜，俞楚天收稿日期：2021-10-21网络首发日期：2022-02-23引用格式：孙轶恺，漆淘懿，张利军，胡怡霜，俞楚天．市场环境下含冰蓄冷空调的综合能源系统优化运行[J/OL]．南方电网技术.https://kns.cnki.net/kcms/detail/44.1643.TK.20220223.1129.024.html网络首发：在编辑部工作流程中，稿件从录用到出版要经历录用定稿、排版定稿、整期汇编定稿等阶段。录用定稿指内容已经确定，且通过同行评议、主编终审同意刊用的稿件。排版定稿指录用定稿按照期刊特定版式（包括网络呈现版式）排版后的稿件，可暂不确定出版年、卷、期和页码。整期汇编定稿指出版年、卷、期、页码均已确定的印刷或数字出版的整期汇编稿件。录用定稿网络首发稿件内容必须符合《出版管理条例》和《期刊出版管理规定》的有关规定；学术研究成果具有创新性、科学性和先进性，符合编辑部对刊文的录用要求，不存在学术不端行为及其他侵权行为；稿件内容应基本符合国家有关书刊编辑、出版的技术标准，正确使用和统一规范语言文字、符号、数字、外文字母、法定计量单位及地图标注等。为确保录用定稿网络首发的严肃性，录用定稿一经发布，不得修改论文题目、作者、机构名称和学术内容，只可基于编辑规范进行少量文字的修改。出版确认：纸质期刊编辑部通过与《中国学术期刊（光盘版）》电子杂志社有限公司签约，在《中国学术期刊（网络版）》出版传播平台上创办与纸质期刊内容一致的网络版，以单篇或整期出版形式，在印刷出版之前刊发论文的录用定稿、排版定稿、整期汇编定稿。因为《中国学术期刊（网络版）》是国家新闻出版广电总局批准的网络连续型出版物（ISSN2096-4188，CN11-6037/Z），所以签约期刊的网络版上网络首发论文视为正式出版。南方电网技术第16卷，2022年SOUTHERNPOWERSYSTEMTECHNOLOGY网络首发文章中图分类号：TK01+9文献标志码：A市场环境下含冰蓄冷空调的综合能源系统优化运行孙轶恺1，漆淘懿2，张利军1，胡怡霜2，俞楚天1（1.国网浙江省电力有限公司经济技术研究院，杭州310008；2.浙江大学电气工程学院，杭州310027）摘要：综合能源系统能够提高能源效率，促进可再生能源消纳。综合能源系统的进一步发展需要优化运行技术的支撑，通过协调其内部各种设备在市场环境下的能源消耗，降低整体的运营成本。因此，本文建立的优化运行技术不仅能够挖掘传统的燃气轮机、蓄电池等设备的响应潜力，还重点优化冰蓄冷空调的电制冷量和融冰制冷量，提高系统的需求响应容量，帮助综合能源系统响应实时电价和辅助服务市场。首先建立了考虑需求响应的设备通用运行模型，并对冰蓄冷空调、燃气轮机、蓄电池等设备进行了精细化建模。随后提出了市场环境下综合能源系统调度运行和需求响应容量的优化方法。然后，将提出的非凸模型转换为可解决的混合整数线性规划问题进行求解。最后，基于实际的综合能源系统运营数据，通过算例仿真验证了所提方法的有效性。关键词：电力市场；综合能源系统；冰蓄冷空调；需求响应；优化运行OptimalOperationofIntegratedEnergySystemIncludingIce-StorageAir-ConditioninginPowerMarketSUNYikai1,QITaoyi2,ZHANGLijun1,HUYishuang2,YUChutian1(1.EconomicandTechnicalResearchInstituteofStateGridZhejiangElectricPowerCo.,Ltd.,Hangzhou310008,China;2.CollegeofElectricalEngineering,ZhejiangUniversity,Hangzhou310027,China)Abstract:Theintegratedenergysystemcanimproveenergyefficiencyandpromotetheconsumptionofrenewableenergy.Thefurtherdevelopmentoftheintegratedenergysystemrequiresthesupportofoptimizedoperationtechnology,andreducestheoveralloperatingcostbycoordinatingtheenergyconsumptionofvariousinternalequipmentintheelectricitymarket.Therefore,theoptimizedoperationtechnologyestablishedinthisarticlecannotonlytaptheresponsepotentialoftraditionalgasturbines,batteriesandotherequipment,butalsofocusonoptimizingtheelectriccoolingcapacityandmeltingicecoolingcapacityoficestorageairconditioners,increasingthedemandresponsecapacityofthesystem,andhelpingtheintegratedenergysystemtorespond.Real-timeelectricitypriceandancillaryservicemarket.First,ageneralequipmentoperationmodelconsideringdemandresponseisestablished,andrefinedmodelingoficestorageairconditioners,gasturbines,batteriesandotherequipmentiscarriedout.Subsequently,anoptimizationmethodforthedispatchoperationanddemandresponsecapacityoftheintegratedenergysystemunderthemarketenvironmentisproposed.Then,theproposednon-convexmodelisconvertedintoasolvablemixedintegerlinearprogrammingproblemforsolution.Finally,basedontheactualoperatingdataoftheintegratedenergysystem,theeffectivenessoftheproposedmethodisverifiedbysimulationexamples.Keywords:electricitymarket;integratedenergysystem;icestorageairconditioning;demandresponse;optimizedoperation0引言为实现“碳达峰，碳中和”的目标，推进能源改革，提高能源利用效率势在必行[1]。综合能源系统由新能源发电、储能、燃气轮机等设备组成，通基金项目：国网浙江省电力有限公司科技项目（5211JY19000V）。Foundationitem:SupportedbytheScienceandTechnologyProjectofStateGridZhejiangElectricPowerCo.,Ltd.(5211JY19000V).过电、热、冷各种能源间的协调互补，在满足系统多能需求的基础上，能够提高能源效率，促进可再生能源的消纳[2]。基于综合能源系统的多种优势，世界上许多国家都在大力推进综合能源系统的研究和落地工程[3-8]。中国北京的大兴机场建设了包含分布式地源热泵，光伏发电和储能等设施的综合能源系统，能够对机场运营中所需的电力、供热和制冷进行有效补充[9]。上海迪士尼园区基于能量梯级利网络首发时间：2022-02-2318:13:53网络首发地址：https://kns.cnki.net/kcms/detail/44.1643.TK.20220223.1129.024.html2南方电网技术第16卷用的冷、热、电和压缩空气四联供技术，将一次能源利用率提升到80%以上[10]。然而，综合能源系统的发展一直被高昂的建设成本和复杂的运营优化所困扰[11-12]。随着通信技术的飞速发展，综合能源系统对内部负荷资源的控制能力不断增强[13]。通过对灵活负荷更加精准的调控，综合能源系统不仅能够响应现货市场中实时电价，还能利用灵活负荷的需求响应能力参与辅助服务市场。山东、江苏、浙江等地针对负荷侧资源参与辅助服务市场已出台相关政策，制定了负荷侧主体参与辅助服务市场的准入条件、结算方法、考核机制等规则[14-16]。因此，综合能源系统通过提供调频备用（regulationreserve,RR）、同步备用（synchronizedreserve,SR）等辅助服务，能够获取市场补贴从而进一步降低系统运行成本。由于燃气轮机、蓄电池等传统设备的响应潜力有限，亟需挖掘新的灵活资源为综合能源系统提供响应容量。冰蓄冷空调包括制冷机、制冰机和蓄冰槽等组件，可以在夜间低电价时期制冰并储存在蓄冰槽内，白天高电价时期融冰制冷从而降低制冷机功率，是综合能源系统重要的制冷和储冷设备[17]。文献[18]对深圳地区的电冰蓄冷空调的典型工程应用进行了经济型研究，结果证明了冰蓄冷空调的削峰填谷作用有利于电网的健康运行和有序发展。文献[19]针对冷热电联供型微网提出了含冰蓄冷空调的多时间尺度优化模型，有效平抑了日内需求侧的随机波动，实现微网的经济稳定运行。文献[20]在含冰蓄冷空调的工厂综合能源系统中，利用冰蓄冷空调夜间蓄冰，白天融冰制冷，降低工厂运营成本。然而，以上研究只关注了冰蓄冷空调的蓄冰时段优化，忽视了冰蓄冷空调在白天的负荷调节能力。白天冰蓄冷空调的制冷量来源于制冷机和融冰两部分，可以根据需求优化其各自的制冷量。在尖峰电价时期，冰蓄冷空调可以增加融冰制冷量，从而降低制冷机的电功率以减少用电成本。此外，不同负荷率下制冷机的能效比不同，优化过程中能效比也是必须要考虑的因素。在高电价时期，大幅降低制冷机功率而增加融冰制冷功率，不仅会使制冷机能效比显著下降，还会加速蓄冰量消耗。因此，需要综合考虑冰蓄冷空调的能效比、蓄冰量、电价等因素，实现其经济最优运行。因此，本文充分挖掘了综合能源系统内冰蓄冷空调、蓄电池、燃气轮机等可调节资源的需求响应潜力，并在现货市场和辅助服务市场中，优化了综合能源系统的多能需求和响应容量，从而降低系统的运营成本。首先，本文针对综合能源系统，提出了考虑需求响应的设备通用运行模型。随后，针对综合能源系统中的冰蓄冷空调、燃气轮机、蓄电池，分别建立了不同的精细化需求响应模型，并拟合了冰蓄冷空调在不同负载率下的能效比曲线。进而以经济最优为目标，基于设备的需求响应模型，采用大M法对优化系统运行中的购电量和需求响应容量。最后通过算例证明了所提优化方法在市场环境下能够有效提高综合能源系统经济性。1综合能源系统结构及其参与市场的条件图1展示了含冰蓄冷空调综合能源系统的结构和能量转移过程，主要包括能量输入、能量生产、能量转换、能量储存和能量消耗五个环节。综合能源系统内能量的主要形式为电、热和冷，因此系统内包含电流，热流和冷流3种能量流动。配电网可再生能源天然气管网蒸汽网络余热锅炉燃气轮机分布式发电吸收式制冷机电热设备制冷机储热罐蓄冰槽蓄电池电负荷热负荷冷负荷电流热流冷流能量输入能量生产能量转换能量储存能量消耗冰蓄冷空调图1含冰蓄冷空调的综合能源系统结构Fig.1Integratedenergysystemstructureincludingice-storageair-conditioning综合能源系统中能量的来源主要有两种形式，分布式发电、燃气轮机和余热锅炉等系统自身设备的生产能量和外部购买的电力、天然气等能量。为满足系统实时变化的各种能量需求，不同类型的能量可通过能量转换设备彼此转换。同时，能量储存设备可以为综合能源系统储存一定的能量，提高系统的灵活性。以冰蓄冷空调为例，其不仅可以通过制冷机将电量转换为冷量，还可以利用自身的蓄冰孙轶恺，等：市场环境下含冰蓄冷空调的综合能源系统优化运行3槽储存冷量。根据需求优化调整能量生产、能量转换和能量储存3个环节，以满足不同类型负荷的用能需求。综合能源系统中设备具有灵活的运行方式和能量储存能力，因此其可以调整内部设备的运行状态，改变其对外界能源输入需求，从而为电力系统提供辅助服务。因此，针对综合能源系统中的冰蓄冷空调、蓄电池等设备，需要评估其能够提供的需求响应能量，通过对设备的优化调节以提供辅助服务获取补贴，从而降低综合能源系统的运行成本。综合能源系统可以提供的辅助服务主要有调频备用和同步备用，在提供需求响应服务时，响应资源需要符合系统运营商设定的响应标准。因此，首先需要对设备的需求响应潜力进行评估，主要包括备用容量，响应时间和持续时间。美国PJM辅助市场中对需求侧资源最小备用容量、响应时间和持续时间有要求[21]：需求侧资源提供调频备用时，最小备用容量为0.1MW，在5min内响应并持续1h；提供同步备用时，最小备用容量0.1MW，在15min内响应并持续1h。最小备用容量是单个参与者需要提供的最小响应容量；响应时间是指从接收需求响应信号到满足需求的时间间隔；持续时间是指需求响应资源可以保持备用容量或跟随需求响应信号的持续时间。尽管现阶段我国包含负荷参与的辅助服务市场建设还不够完善，但是随着电力市场改革的进一步深化，具有负荷调节能力的综合能源系统会逐渐参与到市场中来，响应电力系统的调度指令或者市场价格信号，为系统提供调峰、备用等辅助服务。设备在实时运行过程中，只有满足响应要求时，才能提供辅助服务。为评估各设备的需求响应潜力，提出了考虑需求响应的设备运行模型，如图2所示。当设备提供调频备用时，主要用于平抑可再生能源的出力波动，它需要在持续时间内提供足够的上下容量，如图2(a)所示。当设备提供同步备用时，需要提供下调容量以解决需要削减负荷的能源短缺问题，如图2(b)所示。此外，设备的工作功率可能会在调频备用或同步备用的响应持续时间内发生变化。因此，需要在整个响应期间确保有足够的调频备用或同步备用容量。miniPmaxiPt时刻第i台设备的调频容量第i台设备的最小运行功率t时刻第i台设备的功率,itRR,itRRRRDt调频备用的持续时间要求,itP第i台设备的最大运行功率maxiPminiP(a)调频备用miniPmaxiPt时刻第i台设备的同步容量t时刻第i台设备的功率,itSR,itRRSRDt同步备用的持续时间要求,itP第i台设备的最小运行功率第i台设备的最大运行功率maxiPminiP(b)同步备用图2设备提供辅助服务时的运行模型Fig.2Operationmodelwhenequipmentprovidesauxiliaryservice2综合能源系统设备模型综合能源系统中包含多种设备，不同设备具有显著的响应特性差异。为准确描述设备的运行状态及需求响应潜力，本节针对不同设备，建立了考虑需求响应的精细化模型。2.1冰蓄冷空调冰蓄冷空调主要包括制冷机、制冰机和蓄冰槽3个组件，其运行可以分为制冰蓄冷和融冰制冷两个阶段。夜间气温低于白天，冷负荷需求较少，制冷机可以独立满足冷负荷需求。同时，夜间的电价较低，在夜间制冰蓄冷然后在白天融冰制冷，能够降低运行成本。2.1.1蓄冰过程冰蓄冷空调在夜间蓄冰过程中，蓄冰槽不能同时融冰制冷，因此融冰制冷功率为0，由制冷机满足所有的系统冷功率需求，如式（1）所示：e,cooling,ice,1ice,iceice,icem,0(1)ttiititttiiiQQQSSQT(1)式中：e,tiQ和ice,tiQ分别为第i台冰蓄冷空调在t时刻的制冷机制冷量和融冰制冷量；cooling,tiQ为冰蓄冷空4南方电网技术第16卷调的总制冷量；1ice,tiS和ice,tiS分别为t和t+1时刻的蓄冰量；ice为蓄冰槽的自损耗系数；icem,tiQ为蓄冰制冷量；T为优化周期。蓄冰过程中，冰蓄冷空调的功率为制冷机和制冰机的功率之和，制冷电功率、制冰电功率与冷功率之间的关系可以表示为：e,ice,cooling,e,e,EER,ice,icem,icemtttiiitttiiittiiPPPPQPQI(2)式中：e,tiP和ice,tiP分别为制冷机和制冰机的功率；cooling,tiP为冰蓄冷空调的总功率；EER,ti为制冷机的能效比；icemI为制冰机的能效比。2.1.2融冰过程白天冰蓄冷空调的制冷量由制冷机和蓄冰槽融冰共同提供。根据实时的供冷需求和冰蓄冷空调的运行状态，在可调区间内，制冷功率区间可以灵活分配融冰制冷量和制冷机制冷量，如下式所示：e,ice,cooling,maxe,e,maxe,icee,ice00tttiiitiitQQQQQQQ≤≤≤≤(3)式中：e,tiQ和ice,tiQ分别为第i台冰蓄冷空调在t时刻的制冷机制冷量和融冰制冷量；maxe,iQ和maxe,iceQ分别为制冷机的最大制冷量和融冰的最大制冷量。e,cooling,ice,e,e,EER,0ttiititttiiiPPPPQ(4)式中：,teiP和,ticeiP分别为制冷机功率和制冰机功率；,tcoolingiP为总功率；,tEERi为制冷机的能效比。2.1.3冰蓄冷空调能效比能效比是描述制冷机或者制冰机效率的指标，能效比越高，单位电量就能转换为更多的冷量。冰蓄冷空调制冰机的工作功率是恒定的，其制冰能效比可看作常量。而制冷机的功率会根据冷量需求实时调节，其能效比和负载水平密切相关，不同负载下冰蓄冷空调的效率具有显著差异。本文对冰蓄冷空调进行了部分负荷性能曲线拟合，结果如下：32EERPLPLPLrealPLratedaRbRcRdPRP(5)式中：式中a、b、c和d分别为冰蓄冷空调能效比的拟合系数；PLR为负载比；realP和ratedP分别为制冷机的实际功率和额定功率。从图3中的拟合结果可以看出，在负载率小于80%的时候，制冷机的能效比与负载率大致呈线性关系，随着负载率的增加能效比也随之提高。制冷机的最高能效比能够达到5.5，但是最低能效比只有2.23。因此，在冰蓄冷空调制冷机制冷量和融冰制冷量的优化分配过程中，能效比是不可忽视的一个因素，应尽量保证冰蓄冷空调运行在高效区间。012345600.20.40.60.81能效比负载比图3冰蓄冷空调能效比与负载率之间的关系Fig.3RelationshipbetweenEERandPLRoficestorageairconditioner2.1.4备用约束参与辅助服务市场时，冰蓄冷空调提供调频备用和同步备用的约束如式(6)所示。t时刻冰蓄冷空调的功率e,tiP加上同步备用容量和向上调频备用容量不能超过空调的最大功率；同时，t时刻冰蓄冷空调的功率e,tiP减去向下调频备用容量必须大于等于最低功率。并且在冰蓄冷空调提供的调频备用容量和同步备用容量都不为负值，如式(7)所示。maxe,ice,SRSRIice,RRII,ice,RRRRI0tttttittteiPxSxRRPPxS≤≤(6)RRISRI00ttSS≥≥(7)式中：ice,RRtx和ice,SRtx是0-1变量，当冰蓄冷空调能够提供调频备用或者同步备用的时候值为1；RRItS和SRItS分别为t时刻的同步备用容量和调频备用容量；maxIP为冰蓄冷空调地最大功率。此外，在提供调频备用或者同步备用过程中，从响应开始到响应结束期间，冰蓄冷空调提供的调孙轶恺，等：市场环境下含冰蓄冷空调的综合能源系统优化运行5频备用容量和同步备用容量都必须保持不变，如式(8)所示。SRRR11SRISRISRI11RRIRRIRRItDtttDttSSSSSS(8)式中SRD为提供辅助服务的持续时间。2.2燃气轮机燃气轮机是综合能源系统中的典型能源生产设备，大多数燃气轮机都配备有余热锅炉。燃气轮机以天然气作为能源，输出电能和热能。电能和热能的生产可以表示为：eGGgasGheGGGgasG(1)ttttPFHF(9)式中：GtP,GtH和GtF分别为t时刻的发电量，供热量和天然气消耗量；eG和hG分别表示发电率和发热率；gas为天然气的热值。燃气轮机参与需求响应时的约束和冰蓄冷空调类似，可参照式(6)—(8)得到，此处不再赘述。2.3燃气锅炉燃气锅炉是综合能源系统中重要的能源转换设备，燃气锅炉以天然气为能量来源，产生热蒸汽从而满足综合能源系统的热负荷需求，其产热功率可以表示为：GB,GBgasGB,ttiiHF(10)式中：GB,tiH为燃气锅炉的产热量，GB为燃气锅炉的效率，GB,tiF为天然气消耗速率。燃气锅炉既不消耗和也不产生电能，不能提供辅助服务。通过优化其运行时间，可以补充其他设备由于提供辅助服务而减少的热量供应。2.4吸收式制冷机吸收式制冷机将热量转化为冷量，可以有效利用能源生产过程中的余热，补充综合能源系统的冷负荷需求。吸收式制冷机的制冷功率和吸热量成正比，如式(11)所示[22]。AR,ARAR,ttiiQH(11)式中：AR,tiQ和AR,tiH分别为时刻t第i台吸收式制冷机的制冷功率和耗热功率；AR为吸收式制冷机的转换效率。吸收式制冷机既可以作为综合能源系统供冷的补充，也可以为电力系统提供辅助服务，其参与需求响应的约束可参照式(6)—(8)得到。2.5蓄电池能量存储设备可以存储能量，并在需要时释放能量，打破能量供需同时性的限制。以蓄电池为例，蓄电池的放电和充电过程可以表示为：1ssrrminmaxstartend(1)(/)tttttSSPPTSSSSS≤≤(12)式中：tS为t时刻的存储容量；T为仿真步长；stP和rtP分别为t时刻的放电和充电速率；为电能损耗率；s和r分别为充电和放电速率；maxS和minS分别为最大和最小存储容量。在蓄电池充放电过程中，充放电速率不能超过最大值，其约束可表示为：maxsssmaxrrrsr0001ttttttPPPP≤≤≤≤≤≤(13)式中：st和rt为0-1变量，当电池在t时刻处于充电或放电过程，其值为1；maxsP和maxrP为最大充电和放电速率。为保证在提供需求响应服务的持续时间内都能够提供足够的同步备用和调频备用容量，避免在响应未结束时电池就无法提供响应容量，因此电池的实时电量不能超过最大和最小电量，其约束可表示为：down11rRRsSRsrup11ssRRs(1)(1)(1)(1)()tDtjjjjtDtDjjttDtjtDtDjjjjtSSPSSTSSPST(14)SRRRmindownmaxupmax(,)ttDDDSSSS≤≤(15)式中：D为参与需求响应要求的持续时间；SRD和RRD分别为同步备用和调频备用的持续时间；downtS和uptS分别为蓄电池在tD时刻的最大和最小电量。3市场环境下综合能源系统优化运行模型基于以上考虑需求响应的设备模型，本节提出6南方电网技术第16卷综合能源系统的最优运行模型。然后，对部分非线性约束进行线性化处理，以降低计算难度。最后，通过求解线性化后形成的整数线性规划问题，获得最优运行结果。3.1优化目标模型的目标函数是使综合能源系统的运行成本最低，包括能源购买成本entC（电力购买成本，天然气购买成本和蒸汽购买成本）、设备运行和维护成本omtC、需求响应服务收入drtI和需求响应服务成本drtC。enomdrdrmin()tttttfCCICT(16)enbuyegasbuysteamomomdrSRe,srRRe,rrdrSRc,srRRc,rr()()tttttttttttttttCPFHCPISSCSS(17)式中：为单位价格或单位成本；T为仿真步长；buytP和buytH分别为从外部电网传输的电力和从外部加热网购买的蒸汽。3.2约束条件综合能源系统的运行除了受电力系统约束，供热系统约束，制冷系统约束，系统中各组件还需要满足文章前一节“综合能源系统设备模型”中所述的运行约束。如式（18）所示，综合能源系统产生的电力和购买的电力之和应始终等于电力消耗。同时，综合能源系统中通常有两个交流和直流系统。式（18）是交流系统约束，式（19）是直流系统约束。式（20）表示这两个系统之间的转换。buyGcooling,ACacac-dcttttttiPPPPLP(18)dcG,dcdcsr()tttttPPLPP(19)ac-dcdcdcac-dcdc-acdcdc,0,0tttttPPPPP≥(20)式中：actL和dctL分别是交流和直流负载；ac-dc和dc-ac分别为交流到直流和直流到交流的转换比；G,dctP是直流系统中的分布式光伏等发电设备。系统的热量应满足热水、蒸汽、空间加热以及吸收式制冷机的需求。冷却能量主要满足部分设备的冷却需求和空间冷却。因此，供热和冷却系统的约束可以分别表示为：buyGsrsteamspacewaterAR()()ttttttttHHHHHHHH(21),,ttttARicoolingiloadspaceQQQQ(22)式中：stH和rtH为储热装置t时刻储存和释放热量的速率；steamtH、spacetH和watertH分别表示ｔ时刻系统的蒸汽需求、热量需求和热水需求；loadtQ和spacetQ分别表示系统的设备冷却能量需求和制冷需求。3.3求解方法由于设备参与需求响应的过程中存在非线性约束，求解难度较大，因此，本文采用大M法对约束进行线性化处理。大M法能够将包含0-1变量非线性约束线性化，通过引入人为设定的M值，在0-1变量的取值为0或1时，分别松弛掉两互补条件中不起作用的一个条件[23]。以式(6)为例，使用大M法可以将其转换为式（23）。maxeSRIRRIIeRRISRIice,SRRRIice,RR0tttttttttPSSPPSSxMSxM≤≤≤≤(23)式中，M是一个非常大的正整数。其余各式的非线性约束，可以通过类似的方法转换成线性约束，此处不再赘述。式(20)是一个分段函数，可以用以下线性表达式替换。dcdcdcac-dcdcdcac-dcac-dcac-dcdcdcac-dcdcdcac-dcac-dcac-dcdcdc(1)(1)(1)tttttttttxMPxMPxMPPPxMPxMPPPxM≤≤≤≤≤≤(24)式中，dcx是一个0-1变量，当dctP大于0时为1，反之其值为0。4算例分析基于前述的设备需求响应模型和求解方法，本文以中国广州某工厂为研究对象，建立了含冰蓄冷空调的综合能源系统模型，并利用美国PJM市场中的实际数据进行了仿真分析[24]。综合能源系统中各设备的相关参数如表1所示，每日的交流负荷和直流负荷如图4所示，每日的冷负荷和热负荷的需求如图5所示。孙轶恺，等：市场环境下含冰蓄冷空调的综合能源系统优化运行7表1综合能源系统设备参数Tab.1Equipmentparametersofintegratedenergysystem设备名称参数名称取值冰蓄冷空调制冷机最大功率/kW600蓄冰槽制冷能效比3蓄冰槽融冰效率0.9燃气轮机最大发电功率/MW1最小发电功率/kW50发电效率0.33燃气锅炉最大输出功率/MW1供热效率0.9吸收式制冷机最大耗热功率/MW1制冷能效比0.8蓄电池最大充/放电功率/MW2充/放电效率0.9容量/MWh1初始荷能状态/kWh60001020304050012345670163248648096直流负荷功率(kW)交流负荷功率(MW)时间(15min)交流负荷直流负荷图4工厂日交直流负荷曲线Fig.4AC/DCloadcurvesofintegratedenergysystem01234560163248648096交流负荷功率(MW)时间(15min)冷负荷热负荷图5工厂的日冷热负荷曲线Fig.5Dailycoolingandheatingloadcurvesofintegratedenergysystem表2中列出了综合能源系统内部可以提供不同类型需求响应的设备，Y和N分别表示设备可以提供或不能提供相应的辅助服务。图6展示了2020年6月2日PJM数据中的电价，调频备用的价格和同步备用的价格[25]。并将天然气价格设定为5美元/千立方英尺，蒸汽价格设置为5美元/千立方英尺。仿真的优化步长设置为15min。表2综合能源系统设备及需求响应能力Tab.2ComponentsoftheIESandtheirDRability设备名称调频备用同步备用燃气轮机及余热锅炉NY燃气锅炉NN分布式光伏NN吸收式制冷机YY电池系统YY冰蓄冷空调YY0102030405060708001224364860728496价格($/MWh)时间（15min）实时电价调频备用价格同步备用价格图6实时电价、调频备用价格和同步备用价格Fig.6Electricityprices,RRpricesandSRprices图7为综合能源系统内部的交流侧电量、直流侧电量、冷量和热量的生产及消耗情况。从图7(a)和(b)中可以看出，在49到80的时间段内，电价相对较高。为降低综合能源系统的购电费用，交流侧的燃气轮机开始工作，同时直流侧的电池开始放电以满足系统用电需求。从图7(c)中可以看出，冰蓄冷空调也降低了用电功率，减少的制冷量由融冰制冷和吸收式制冷机共同补充以满足系统需求。图7(d)中余热锅炉随着燃气轮机的启动也开始产生热量，因此系统的购热量也随之降低。图8为冰蓄冷空调的运行情况，从图中可以看出，在1-40时段，此时电价最低，为冰蓄冷空调的蓄冰过程，因此融冰制冷量小于0，此时冰蓄冷空调的制冷量全部由制冷机提供。在40-84时段，冷负荷需求较高，同时电价也高于夜间，因此蓄冰槽开始融冰制冷，和制冷机联合满足冷负荷需求。在考虑不同负载率下制冷机的能效比差异显著的情况下，制冷机一直维持在相对稳定的功率，因此能效比也一直维持在较高水平。8南方电网技术第16卷00.40.81.21.62024681001224364860728496燃气轮机及冰蓄冷空调功率(MW)购电及负荷功率(MW)时间（15min）购电负荷燃气轮机冰蓄冷空调(a)交流侧电力生产及消耗0100200300400-4-202401224364860728496光伏及直流负荷功率(kW)电池功率(MW)时间（15min）电池光伏直流负荷(b)直流侧电力生产及消耗00.511.522.5301224364860728496冷功率(MW)时间(15min)冷负荷需求冰蓄冷空调吸收式制冷机(c)冷量生产及消耗012345601224364860728496热功率(MW)时间(15min)热负荷购热余热锅炉燃气锅炉(d)热量生产及消耗图7综合能源系统优化运行结果Fig.7Optimizedoperationresultsofintegratedenergysystem图9为冰蓄冷空调和电池的辅助服务提供情况。冰蓄冷空调可以同时提供调频备用和同步备用，空调会根据实时的调频备用和同步备用价格调整制234567-101234501224364860728496能效比冷功率(MW)时间(15min)冷负荷需求电制冷量融冰制冷量EER图8冰蓄冷空调运行情况Fig.8Operationstatusoficestorageair-conditioning冷机制冷量和融冰制冷量。在50-56阶段，调频备用价格较高，因此冰蓄冷空调降低电功率，增加融冰量满足系统需求。例如在64-67时段，调频备用价格大幅降低，因此制冷机增加功率，减少融冰量以用在后期需要降低制冷机功率的时期。040080012001600200001224364860728496功率（kW）时间（15min）空调RR电池RR空调SR电池SR图9冰蓄冷空调和电池的辅助服务提供情况Fig.9Ancillaryservicesprovidedbyicestorageair-conditioningandbatterysystem表3为冰蓄冷空调优化前后的运行费用对比。从表中可以看出，冰蓄冷空调运行优化后，购电费用稍有增加。这是由于在提供辅助服务时，冰蓄冷空调需要实时调整电制冷量和融冰制冷量，导致效率稍有降低。但是通过提供调频备用和同步备用服务，冰蓄冷空调分别获得了115和24美元的收益，最终运行成本从278美元下降为153美元，大幅降低了冰蓄冷空调的运行成本。表3冰蓄冷空调优化前后运行费用Tab.3Operationcostofice-storageairconditionerbeforeandafteroptimization美元费用类别购电费用调频备用收益同步备用收益总费用优化前27800278孙轶恺，等：市场环境下含冰蓄冷空调的综合能源系统优化运行9优化后28211524153电池相较于冰蓄冷空调，主要提供调频备用。在11-12时段，由于调频备用价格突然增加，电池的充电功率迅速从1.14MW下降到0，将调频备用从270kW迅速提升到了1500kW。在65-69时段，调频备用价格大幅下降，电池也随之将调频备用功率降为0。但是在提供辅助服务期间，电池不会完全充满电或者用完电，需要一直为调频备用保留一定的余量。参与辅助服务可以为电力系统提供更多可调节的容量，并减少综合能源系统的运营成本。表4为综合能源系统优化前后的运行费用对比，综合能源系统运行成本从3762美元降低至3160美元，降低了16％。和冰蓄冷空调类似，优化前后系统的购电和购热费用几乎没有差别。系统运行总费用的降低主要来自于提供同步备用和调频备用获得的收益。表4综合能源系统优化前后运行费用Tab.4OperationcostofIESbeforeandafteroptimization美元费用类别优化前优化后购电总费用22322236购热总费用15301534同步备用收益0577调频备用收益033运行总费用376231605结语综合能源系统的最优运行在消纳可再生能源和提高能源效率方面发挥着重要作用。为降低综合能源系统运营成本并促进其发展，本文充分发掘了冰蓄冷空调的调度潜力，通过优化冰蓄冷空调的电制冷量和融冰制冷量，提高了综合能源系统的需求响应容量。仿真结果表明，市场环境下综合能源系统能够根据实时的电价、调频备用价格和同步备用价格，调度系统内各设备以实现系统最优运行。经过优化后的综合能源系统，冰蓄冷空调的运行成本降低了45%，系统运行成本降幅也达到了16%，具有显著的经济效益。参考文献[1]朱继忠,何晨可,陈婧韵,等.综合能源系统环境下电动汽车充换电设施规划的综述[J/OL].南方电网技术:1-19.[2022-01-18].http://kns.cnki.net/kcms/detail/44.1643.TK.20211208.1010.002.html.[2]汤翔鹰,胡炎,耿琪,等.考虑多能灵活性的综合能源系统多时间尺度优化调度[J].电力系统自动化,2021,45(4):81-90.TANGXiangying,HUYan,GENGQi,etal.Multi-time-scaleopti-malschedulingofintegratedenergysystemconsideringmulti-energyflexibility[J].AutomationofElectricPowerSystems,2021,45(4):81-90.[3]陈皓勇,陈思敏,陈锦彬,等.面向综合能源系统建模与分析的能量网络理论[J].南方电网技术,2020,14(2):62-74.CHENHaoyong,CHENSimin,CHENJinbin,etal.Energynetworktheoryformodelingandanalysisofintegratedenergysystems[J].SouthernPowerSystemTechnology,2020,14(2):62-74.[4]郇嘉嘉,张小辉,黄学劲,等.考虑负荷不确定性与运行可靠性的综合能源系统资产利用效率分析[J].南方电网技术,2021,15(2):56-63.XUNJiajia,ZHANGXiaohui,HUANGXuejin,etal.Assetutiliza-tionefficiencyanalysisofintegratedenergysystemconsideringloaduncertaintyandoperationreliability[J].SouthernPowerSystemTechnology,2021,15(2):56-63.[5]郑涛,戴则梅,姚家豪,等.综合能源系统控制自由度指标及其对经济调度的影响[J].中国电力,2021,54(4):94-106,118.ZHENGTao,DAIZemei,YAOJiahao,etal.Controlfreedomindexofintegratedenergysystemanditsimpactoneconomicdispatch[J].ElectricPower,2021,54(4):94-106,118.[6]廖宗毅,万文略,陈曦.考虑价格引导机制的园区综合能源系统日前优化调度方法[J].南方电网技术,2021,15(9):53-60.LIAOZongyi,WANWenluo,CHENXi.Day-aheadoptimalsched-ulingmethodofparkintegratedenergysystemconsideringpriceguidancemechanism[J].SouthernPowerSystemTechnology,2021,15(9):53-60.[7]王昀,谢海鹏,孙啸天,等.计及激励型综合需求响应的电-热综合能源系统日前经济调度[J].电工技术学报,2021,36(9):1926-1934.WANGYun,XIEHaipeng,SUNXiaotian,etal.Day-aheadeconom-icdispatchforelectricity-heatingintegratedenergysystemconsider-ingincentiveintegrateddemandresponse[J].TransactionsofChinaElectrotechnicalSociety,2021,36(9):1926-1934.[8]李鹏,吴迪凡,李雨薇,等.基于综合需求响应和主从博弈的多微网综合能源系统优化调度策略[J].中国电机工程学报,2021,41(4):1307-1321.LIPeng,WUDifan,LIYuwei,etal.Optimaldispatchofmul-ti-microgridsintegratedenergysystembasedonintegrateddemandresponseandstackelberggame[J].ProceedingsoftheCSEE,2021,41(4):1307-1321.[9]湖南睿日光伏.100%绿电供应大兴机场成为全球绿色能源新地标[EB/OL].搜狐网.(2019-10-14).https://www.sohu.com/a/346905110_100038373.[10]彭克,张聪,徐丙垠,等.多能协同综合能源系统示范工程现状与展望[J].电力自动化设备,2017,37(6):3-10.PENGKe,ZHANGCong,XUBinggen,etal.Statusandprospectofpilotprojectsofintegratedenergysystemwithmulti-energycollabo-ration[J].ElectricPowerAutomationEquipment,2017,37(6):3-10.[11]BAOM,DINGY,SINGHC,etal.Amulti-statemodelforreliabilityassessmentofintegratedgasandpowersystemsutilizinguniversalgeneratingfunctiontechniques[J].IEEETransactionsonSmartGrid,2019,10(6):6271-6283.10南方电网技术第16卷[12]陈晚晴,穆云飞,贾宏杰,等.考虑设备变工况特性的区域综合能源系统优化调度方法[J].电网技术,2021,45(3):951-958.CHENWanqing,MUYunfei,JIAHongjie,etal.Operationoptimiza-tionmethodforregionalintegratedenergysystemconsideringpart-loadperformancesofdevices[J].PowerSystemTechnology,2021,45(3):951-958.[13]HUIH,DINGY,SHIQ,etal.5Gnetwork-basedInternetofThingsfordemandresponseinsmartgrid:Asurveyonapplicationpotential[J].AppliedEnergy,2020(257):113972.[14]国家能源局山东监管办公室.山东能源监管办关于《山东电力辅助服务市场运营规则(试行)(2021年修订版)(征求意见稿)》公开征求意见的公告[EB/OL].(2021-07-21).http://www.cnea.gov.cn/tzgg/content_6873.[15]国家能源局江苏监管办公室.关于印发《江苏电力市场用户可调负荷参与辅助服务市场交易规则（试行）》的通知[EB/OL].(2020-11-10).http://jsb.nea.gov.cn/news/2020-11/20201110100742.htm.[16]新华网.浙江开展负荷响应专项市场试点交易[EB/OL].http://www.zj.xinhuanet.com/2021-09/06/c_1127833292.htm,2021.[17]徐笑强,陈振乾.基于负荷预测的冰蓄冷空调系统优化运行[J].建筑热能通风空调,2018,37(1):1-5,10.XUXiaoqiang,CHENZhenqian.Studyoneconomicoptimalcontrolstrategyoficestorageair-conditioningsystembasedonneuralnet-workprediction[J].BuildingEnergy&Environment,2018,37(1):1-5,10.[18]丁庆,段绍辉,王执中,等.冰蓄冷空调在高峰谷负荷差地区应用的经济性[J].电力系统及其自动化学报,2014,28(1):72-75,80.DINGQing,DUANShaohui,WANGZhizhong,etal.Economyofice-storageair-conditionusedintheareaofhighpeak-valleyloaddifference[J].ProceedingsoftheCSU-EPSA,2014,28(1):72-75,80.[19]程杉,黄天力,魏荣宗.含冰蓄冷空调的冷热电联供型微网多时间尺度优化调度[J].电力系统自动化,2019,43(5):30-38.CHENGShan,HUANGTianli,WEIRongzong.Multi-time-scaleoptimalschedulingofcchpmicrogridwithice-storageair-conditioning[J].AutomationofElectricPowerSystems,2019,43(5):30-38.[20]孙可,何德,李春筱,等.考虑冰蓄冷空调多模式的工厂综合能源系统多能协同优化模型[J].电力建设,2017,38(12):12-19.SUNKe,HEDe,LIChunxiao,et.al.Multi-energycooperativeop-timizationmodeloffactoryiesconsideringmulti-modeloficestor-age[J].ElectricPowerConstruction,2017,38(12):12-19.[21]PJM.Day-AheadandReal-TimeMarketOperations[Z]//PJM.PJMManual11:Energy&AncillaryServiceMarketOperations,Penn-sylvania:PJM,2019.[22]徐航,董树锋,何仲潇,等.考虑能量梯级利用的工厂综合能源系统多能协同优化[J].电力系统自动化,2018,42(14):123-130.XUHang,DONGShufeng,HEZhongxiao,et.al.Multi-energyco-operativeoptimizationofintegratedenergysysteminplantconsid-eringsteppedutilizationofenergy[J].AutomationofElectricPowerSystems,2018,42(14):123-130.[23]丁涛,孙宏斌,柏瑞,等.考虑最大风电容量接入的带禁止区间实时经济调度模型[J].中国电机工程学报,2015,35(4):759-765.DINGTao,SUNHongbin,BORui,etal.Real-timeeconomicdis-patchmodelwithprohibitedzonesconsideringmaximumwindpow-erinjection[J].ProceedingsoftheCSEE,2015,35(4):759-765.[24]何仲潇.多能协同的综合能源系统协调调度方法研究[D].杭州：浙江大学,2018.[25]PJM.Markets&operations[EB/OL].(2020-05-13).https://pjm.com/markets-andoperations.aspx.收稿日期：2021-10-21作者简介：孙轶恺（1988），男，工程师，硕士，研究方向为储能技术，综合能源系统，sunyikai@zj.sgcc.com.cn；漆淘懿（1998），男，通信作者，硕士研究生，研究方向为需求管理，储能技术，qty_sgool@163.com；张利军（1973），男，高级工程师（教授级），硕士，研究方向为电力系统运行与控制，zhanglijun@zj.sgcc.com.cn。