基于碳排放流的综合能源系统碳排放监测方法.pdf

第 42 卷 第 10 期2023 年10 月Zhejiang Electric PowerVol.42，No.10Oct.25.2023基于碳排放流的综合能源系统碳排放监测方法刘子华1，曹瑞峰2，赵志扬1，李哲1，娄伟明1，陈婷1，邢海军3（1.国网浙江省电力有限公司台州供电公司，浙江 台州 318000；2.国网浙江省电力有限公司，杭州 310007；3.上海电力大学 电气工程学院，上海 200090）摘要：针对电-气耦合的综合能源系统，提出了一种基于碳排放流的碳排放监测方法，将发电侧碳排放分摊到各节点负荷与各支路功率上，从而实现碳排放的准确追踪与溯源。首先介绍了碳排放流原理；接着研究碳排放流在电网及气网中的应用，建立电网及气网碳排放监测模型，并基于碳排放流理论提出一种碳排放监测方法；最后通过仿真表明，该综合能源系统碳排放监测方法能够同时表明系统时间与空间上的碳足迹，克服传统碳计量方法统计过于宏观的局限性。关键词：综合能源系统；碳排放流；碳监测；电-气耦合DOI:10.19585/j.zjdl.202310008 开放科学（资源服务）标识码（OSID）：A carbon emission monitoring method for integrated energy systems based on carbon emission flowLIU Zihua1，CAO Ruifeng2，ZHAO Zhiyang1，LI Zhe1，LOU Weiming1，CHEN Ting1，XING Haijun3（1.State Grid Taizhou Power Supply Company，Taizhou，Zhejiang 318000，China；2.State Grid Zhejiang Electric Power Co.，Ltd，Hangzhou 310007，China；3.School of Electrical Engineering，Shanghai University of Electric Power，Shanghai 200090，China）Abstract:A carbon emission monitoring method based on carbon emission flow is proposed for the integrated electricity and gas systems（IEGS）.This method distributes carbon emissions on the generation side to various node loads and branch power，thereby achieving accurate tracking and traceability of carbon emissions.Firstly，the principle of carbon emission flow is introduced.Subsequently，the application of carbon emission flow in the power grid and gas network is investigated.A carbon emission monitoring model for the power grid and gas network is established，and a carbon emission monitoring method is proposed based on the theory of carbon emission flow.Finally，simulation results demonstrate that this carbon emission monitoring method for the integrated energy system can simultaneously elucidate the temporal and spatial carbon footprint of the system，overcoming the limitations of conventional carbon measurement methods that tend to be overly macroscopic.Keywords:integrated energy systems；carbon emission flow；carbon monitoring；integrated electricity and gas0引言全球气候变化已成为人类社会面临的一个重大难题，由温室气体排放导致的温升亟待控制。为减缓气候变化，降低温室气体排放，我国承诺“碳排放力争于2030年前达到峰值，争取2060年前实现碳中和”。电力行业作为碳排放的重点领域，也是实现“碳达峰、碳中和”目标的主力军，将承担更大的减排责任，以电力行业晚达峰支撑全社会尽早达峰。其中，综合能源系统以多种形式能源进行生产、输送和消费，实现多种能源之间的互动，是实现低碳的强有力工具。利用综合能源多能耦合机制促进协同互补运行，能够有效提升区域综合能源系统的经济性，此前已有较多研究，而区域能源系统的减碳潜力研究则较少1-5。碳计量作为“碳达峰、碳中和”的重要手段，需要精准地计量数据，并融入碳中和能源管理平台，实现碳足迹的可视化，进而有针对性地对电力系统进行节能降碳，实现低碳转型。从“汇”与“源”的角度来看，碳监测获取的基础信息包括碳汇状况、环境中温室气体浓度和温室气体基金项目：国网浙江省电力有限公司科技项目（5211TZ220002）第 42 卷排放强度三方面。“排放是源头，是 加 的过程，碳汇是消解，是 减 的过程，而环境中浓度可以理解为加减后的存量。”通过碳监测，能够起到主动适应气候变化需求，加强气候变暖对我国承受力脆弱地区影响的观测和评估，也可以促进社会从各个方面落实减污降碳、源头治理要求，因此碳监测具有重要意义6。目前，综合能源系统中的碳计量与监测方式主要为宏观统计法7-10。宏观统计法多用于国家层面电力系统宏观碳排放量的统计，可以得到系统在较长周期内的碳排放总量，直接且实用。然而，此方法未考虑电力系统具体的碳排放过程，难以实现电力系统碳排放的精细化分析，且易受总量统计以及碳排放因子测算误差影响。此外，部分学者从全生命周期法的时间角度对系统的碳排放量进行统计。一种产品从原材料开采生产开始，到加工、制造、运输，再到消费者使用、回收，最终进入循环或作废弃物处理，这一整个过程称为产品的全生命周期11。在电力系统碳排放计量中应用LCA（生命周期评价）进行考虑，可以有效规避仅关注燃料消耗产生的碳排放的局限性12-13。但是，电力系统碳排放具有动态转移特性，LCA仍然缺乏这方面的考虑，无法清晰描述时空转移特性14。碳排放流法从空间角度处理碳排放与系统潮流之间的关系，可以得到特定时间点或时间段的系统碳排放区域分布，进行系统各节点实时、准确、全面的监测15-16。用碳排放流法能够将发电侧碳排放分摊到各节点负荷、各支路功率以及网络损耗，从而实现碳排放具体流向的准确追踪与溯源17。因此，本文提出面向电-气耦合的综合能源系统的碳排放监测方法。首先，分别建立电网与气网碳排放流的碳监测模型，接着通过EH（能源枢纽）形成电-气耦合，进而研究电-气耦合下的综合能源系统的碳排放监测，获取动态实时的碳排放信息。1碳排放流原理1.1碳排放流概念碳排放流的概念最早用于分析地区间由贸易引起的碳排放转移。从空间角度处理碳排放与系统潮流之间的关系，可以得到特定时间点或时间段的系统碳排放区域分布。碳排放流被定义为与能源系统中的能量流相关的不加选择的CO2（二氧化碳）排放，其中碳排放在发电侧产生，碳排放流的目的是将能源生产和转换方的所有碳排放量累积到各节点及需求方。因此，每吨CO2排放量为综合能源系统从能源产生、输送、转化到终端消费全过程的能量流进行核算。在综合能源系统中，既有天然气、煤炭等一次能源，也有电力、热力、制冷等二次能源。因此，碳排放包括实际碳排放（如燃气燃烧）和虚拟碳排放（如与电相关的碳排放）。对于一次能源网络，碳排放流表示将在后期转换或消费过程中发生的碳排放流量。对于二次能源网络，碳排放流表示能量产生时系统的碳排放流量18。1.2碳排放流基本原则对电-气耦合的综合能源系统碳排放进行监测的关键是建立碳排放流理论及其计算模型，用以合理分析综合能源系统的碳排放。潮流追踪法的核心原则是比例共享原则，依托潮流追踪法建立碳排放流模型更加合适且准确19。下面用一个连接有4条支路的节点m进行说明。如图1所示，节点m的进线1和进线2的功率分别为P1和P2，出线3和出线4的功率分别为P3和P4。根据比例共享原则，输出支路的功率由各输入支路按功率比例提供，输入支路的功率也由各输出支路按功率比例分配。即出线功率P3由进线1提供的功率为P1P1+P2P3，由进线2提供的功率为P2P1+P2P3；而进线功率P2在支路3中分配的功 率 为P3P3+P4P2，在 支 路 4 中 分 配 的 功 率为P4P3+P4P2。同时，在该方法中需要考虑以下两个基本图1 比例共享原则Fig.1Proportional sharing principle66 第 10 期刘子华，等：基于碳排放流的综合能源系统碳排放监测方法原则。1）能量合并原理：在综合能源系统中，当来自不同分支的能量注入到一个节点时，流入的碳排放是混合的，这个节点的碳强度eNn被认为是所有流入支路碳流强度Bi的加权平均值。eNn=i+nBiQi/i+nQi（1）式中：+n代表节点n的流入分支集合；Qi表示流入分支i的能量流。2）能量调度原则：根据比例共享原则，从一个节点流出的能量按等比例分享该节点的能量注入。以节点n为例，第j个流出分支的能量流来自第i个流入分支的份额可以表示为：Qj，i=QjQi/s+nQs，j-n（2）式中：-n代表节点n的流出分支集合。结合式（1）和式（2），就可以得到节点n的每个流出分支的支路碳流强度等于该节点的碳势。需要注意的是，如果存在向节点n注入能量流的发电机，则每个发电机都可以看作是一个流入分支，将发电机的支路碳流强度替换为发电机的碳流强度，以上两个原则仍然成立。2碳排放监测模型2.1电网碳排放监测模型从电网角度来看，因为发电技术不同，所以会导致不同的碳排放特性与碳排放量。故依托潮流构建电网碳排放流程理论体系。从表面上看，电网的碳排放流是在给每条支路的潮流概念上再添加一个碳流的标记。从表层看，CO2是由发电厂排放到大气中的；但是从深层次来说，CO2排放是由用户侧消费所驱动产生的。因为碳流与潮流是紧密联系的，可以认为碳流随潮流一起流动，共同从发电厂流出并最终流入到负荷节点上。目前，随着碳排放流理论的逐步完善以及研究的不断深入，可以发现碳排放流模型基于比例分担原则能起到比较好的效果。以下结合图2简述碳排放流理论的概念与计算。碳流率指的是单位时间t内通过电网节点或者是电网支路潮流所对应的碳流排放，用R表示。R=dCdt（3）式中：C表示流入的CO2排放量；t表示时间。碳流密度指的是单位电量所对应的 CO2排放量。发电机组碳排放强度指的是发电厂依据发电特性实时发电所产生的CO2排放强度。支路碳流密度指的是在电网支路中传输单位电量导致发电厂侧所需排放的CO2量，即支路碳流率与有功潮流的比值，用表示。=RP（4）式中：P为支路对应的有功功率。支路碳流密度应与流入该支路首节点的碳势相等。节点碳势指的是将节点消耗单位电量所造成的等效于发电侧的碳排放量，等于所有流入节点的支路的碳流密度关于有功功率的加权平均，例如节点n的碳势为：En=k=1KPk，nk，n+PG，nEG，nk=1KPk，n+PG，n（5）式中：En为电网节点n的碳势；Pk，n为注入电网节点n的第k条支路的功率；PG，n为电网节点n所接发电机的输出功率；k，n为注入电网节点n的第k条支路的碳流密度；K为与此节点相连的注入功率的支路数量；EG，n为电网节点n所接发电机的碳排放强度。例如针对图2，节点6碳势为：E6=P1，61，6+PG，6EG，6P1，6+PG，6（6）式中：P1，6为注入电网节点6的第1条支路的功率；1，6为注入电网节点6的第1条支路的碳流密度；PG，6为电网节点6所接发电机的输出功率；EG，6为图2 电网碳排放流示意图Fig.2Schematic diagram of carbon emission flow of power grid67第 42 卷电网节点6所接发电机的碳排放强度。节点负荷的碳排放量可以通过节点碳势计算得到：Tn=EnLn（7）式中：Tn为节点n的碳排放量；Ln为节点n的负荷量。2.2气网碳排放监测模型在综合能源系统的气网中，因为天然气燃烧后会产生CO2，因此需要建立燃气网的碳排放流模型。气网的碳排放流模型与电网碳排放流类似，同理，根据能量分配原则，每条管道的支路碳流密度等于其流入节点的节点碳势。气网碳排放流如图3所示，满足如下关系：Egas，n=Q，n，n+k=1KgQgas，kgas，kQ，n+k=1KgQgas，k（8）式中：Egas，n为天然气网络节点n的碳势；Q，n为气井的气流量；，n为气井的碳排放强度；Qgas，k为管道k气流量；gas，k为管道k中气流量的碳势；Kg为节点n处流入气流量的支路数量。同理，气网节点n的碳势等于所有流入节点的支路碳流密度关于气网潮流的加权平均。3电-气耦合的综合能源系统碳排放监测方法综合能源系统碳足迹监测流程如图4所示，具体步骤如下：1）计算系统潮流，统计发电机组碳排放强度向量EG。2）生成支路潮流分布矩阵PB、机组注入分布矩阵 PG、负荷分布矩阵 PL和节点有功通量矩阵PN。3）判断PN的绝对值是否等于0。若PN等于0，将与PN中对角元为0的所在行对应节点相连的机组和线路从电网中去除，并返回到步骤2；若PN不等于0，计算节点碳势En、系统支路和各节点的碳流分布。4）根据碳排放流理论对所得区域能源系统中的碳流分布进行碳足迹监测，展示碳足迹数据及趋势。其中，源侧产生的碳排放传输到EH的输入端后，还需要将其通过EH内部能流网络分配至负荷端，因此需对EH进行碳排放流建模。EH结构如图5所示，每个EH的输入端为上级电网和上级气网输送过来的能源；内部元件主要包括光伏、电转气设备、热电联产、燃气锅炉、热泵和储能等；输出端由电、气、热负荷组成。EH中包含多个能量转换器，将EH中的能量转化过程建模为具有单个或多个输入和输出端口图4 系统碳排放监测流程Fig.4Flow chart of system carbon emission monitoring图3 气网碳排放流示意图Fig.3Schematic diagram of carbon emission flow of gas network68 第 10 期刘子华，等：基于碳排放流的综合能源系统碳排放监测方法的单元，与上级电网及上级气网中的节点碳势产生联系。其中，电转气设备、热泵、燃气锅炉为单输入单输出设备，热电联产为单输入多输出设备。对系统内各能量转换公式及碳排放平衡约束组合求解，可得EH负荷端的碳势以及碳排放量。电转气设备将输入的电能转换为天然气，热泵将输入的电能转换为热能，燃气锅炉将输入的天然气转换为热能，碳排放流模型如下：outptgQoutptg=inptgPinptg（9）outebHouteb=inebPineb（10）outgbHoutgb=ingbQingb（11）式中：inptg、ineb、ingb分别为电转气设备、热泵、燃气锅炉输入端的碳势；outptg、outeb、outgb分别为电转气设备、热泵、燃气锅炉输出端的碳势；Pinptg、Pineb、Qingb分别为电转气设备、热泵、燃气锅炉的输入能流功率；Qoutptg、Houteb、Houtgb分别为电转气设备、热泵、燃气锅炉的输出能流功率。上述转换设备的输入输出约束如下：Qoutptg=ptgPinptg（12）Houteb=ebPineb（13）Houtgb=gbQingb（14）式中：ptg、eb、gb分别为电转气设备、热泵、燃气锅炉的转换效率。热电联产作为单输入多输出设备，将输入的天然气转换为电能及热能，其碳排放流模型及设备转换约束如下：inchpQinchp=outchp，hHoutchp+outchp，pPoutchp（15）Houtchp=chp，hQinchp（16）Poutchp=chp，pQinchp（17）式中：inchp为热电联产输入端的碳势；Qinchp为热电联产输入端的能流功率；outchp，h和outchp，p分别为热电联产热、电输出端的碳势；Houtchp和Poutchp分别为热电联产热、电输出端的能流功率；chp，h和chp，p分别为热电联产的热、电转换效率。依据碳排放平衡原则，EH负荷端的碳排放与输入的碳排放应保持平衡，约束如下：EnPi+vPv+outchp，pPoutchp-inptgPinptg-inebPineb=Le，jle，j（18）Egas，nQi+outptgQoutptg-ingbQingb-inchpQinchp=Lg，jlg，j（19）outebHouteb+outchp，hHoutchp+outgbHoutgb=Lh，jlh，j（20）式中：v为光伏输入的碳排放强度，一般认为光伏为清洁能源，本文取值为0 kg/MWh；Pv为光伏的输出功率；Le，j、Lg，j、Lh，j分别为第j个EH的电、气、热负荷；le，j、lg，j、lh，j分别为第j个EH的电、气、热负荷对应的碳势。4算例分析对包含6节点电网、6节点天然气网和3个EH的综合能源系统进行碳排放监测。该系统结构如图6所示，G1G3为电厂燃煤机组，EH1EH3为能量枢纽，GW1、GW2为天然气气井。虚线部分代表通过EH中的转换设备形成电网与气网的耦合。EH1基础数据曲线如图7所示，EH2和EH3的负荷为EH1的一半。燃煤机组G1G3碳排放强度分别取800 kg/MWh、600 kg/MWh、400 kg/MWh。GW1和GW2的碳排放强度均取200 kg/MWh。其中，电网及其气网潮流计算模型参照文献 20建立。图5 EH结构Fig.5Structure of EH图6 综合能源系统结构Fig.6Diagram of the integrated energy system structure 69第 42 卷首先进行电网及气网的潮流计算，得到电网碳流计算相关参数PB、PG、PL、PN以及气网碳流计算相关参数QB、QG、QL、QN；进而监测综合能源系统中电网各节点、气网各节点的碳势及从EH角度产生的全天总碳排放量，监测结果见图 811。EH中的各设备将电网与气网进行耦合，结合碳排放流分析法，产生的碳排放传输到EH的输入端后，需要将其分配至各EH负荷端以进行荷侧的具体监测。通过各个能量转换公式的组合求解，可得出EH负荷端的碳势以及碳排放量。图1214分别为EH电、气、热负荷碳排放强度监测情况。可以看出，碳排放强度随潮流变化而改变，EH2与EH3所设定的电、气、热负荷及转换设备均一致，然而因所处的节点不同，碳势不一致，进而导致监测到的日碳排放强度曲线不同。图15为电-气耦合的综合能源系统碳排放量监测情况。可以看出，总体与负荷变化趋势一致，能真实反映负荷端的碳排放情况。5结语本文建立了电-气综合能源系统碳排放监测模型，采用碳排放流法分别对电网及气网进行建模，并对电-气耦合综合能源系统中的设备进行分析，提出一种电-气耦合的综合能源系统碳排放监测方图7 EH1基础数据曲线Fig.7Basic data curve of EH1图8 电网碳排放强度Fig.8Carbon emission intensity of power grid图9 电网碳排放量Fig.9Carbon emissions of power grid图10 气网碳排放强度Fig.10Carbon emission intensity of gas network图11 气网碳排放量Fig.11Carbon emissions of gas network70 第 10 期刘子华，等：基于碳排放流的综合能源系统碳排放监测方法法。通过算例验证分析，得出以下结论：根据系统本身的电网潮流及气网潮流，利用碳排放流法进行碳流分析，可以同时明晰系统时间与空间上的碳足迹，克服传统碳计量方法统计过于宏观的局限性；通过对电-气耦合的综合能源系统各个节点的碳排放强度及碳排放量进行实时监测，可为精准降碳及需求侧降碳提供信息基础。参考文献 1 王海洋，李珂，张承慧，等.基于主从博弈的社区综合能源系统分布式协同优化运行策略 J.中国电机工程学报，2020，40（17）：5435-5445.WANG Haiyang，LI Ke，ZHANG Chenghui，et al.Distributed coordinative optimal operation of community integrated energy system based on stackelberg game J.Proceedings of the CSEE，2020，40（17）：5435-5445.2 何宇斌.电-气综合能源系统的随机优化和分布式调度研究 D.杭州：浙江大学，2019.HE Yubin.Research on stochastic optimization and decentralized dispatch for integrated electricity-natural gas system D.Hangzhou：Zhejiang University，2019.3 王俐英，林嘉琳，宋美琴，等.考虑需求响应激励机制的园区综合能源系统博弈优化调度 J/OL.控制与决策：1-92022-11-20.DOI：10.13195/j.kzyjc.2022.0436.WANG Liying，LIN Jialin，SONG Meiqin，et al.Optimal dispatch of park integrated energy system considering demand response incentive mechanismJ/OL.Control and Decision：1-9 2022-11-20.DOI：10.1 3195/j.kzyjc.2022.0436.4 陈乾，张沈习，程浩忠，等.计及热网蓄热特性的多区域综合能源系统多元储能规划 J/OL.中国电机工程学报：1-142022-11-20.DOI：10.13334/j.0258-8013.pcsee.221515.CHEN Qian1，ZHANG Shenxi1，CHENGHaozhong，et al.Multiple energy storage planning of multi-district integrated energy system considering heat storage characteristics of heat network J/OL.Proceedings of the CSEE：1-142022-11-20.DOI：10.13334/j.0258-8013.pcsee.221515.5 朱建全，刘海欣，叶汉芳，等.园区综合能源系统优化运行研究综述 J.高电压技术，2022，48（7）：2469-2482.ZHU Jianquan，LIU Haixin，YE Hanfang，et al.Review on 图12 EH电负荷碳排放强度Fig.12Carbon emission intensity of EH electric load图13 EH气负荷碳排放强度Fig.13Carbon emission intensity of EH gas load图14 EH热负荷碳排放强度Fig.14Carbon emission intensity of EH thermal load图15 电-气耦合的综合能源系统碳排放量Fig.15Carbon emissions of IEGS71第 42 卷optimal operation of park-level integrated energy systemJ.High Voltage Engineering，2022，48（7）：2469-2482.6 黄雨涵，丁涛，李雨婷，等.碳中和背景下能源低碳化技术综述及对新型电力系统发展的启示 J.中国电机工程学报，2021，41（增刊1）：28-51.HUANG Yuhan，DING Tao，LI Yuting，et al.Decarbonization technologies and inspirations for the development of novel power systems in the context of carbon neutralityJ.Proceedings of the CSEE，2021，41（S1）：28-51.7 PAUSTIAN K，RAVINDRANATH N H，AMSTEL A V.2006 IPCC guidelines for national greenhouse gas inventories J.International Panel on Climate Change，2006.DOI：http：/dx.doi.org/10.1016/S1462-9011（99）00023-4.8 崔杨，邓贵波，赵钰婷，等.考虑源荷低碳特性互补的含风电电力系统经济调度 J.中国电机工程学报，2021，41（14）：4799-4815.CUI Yang，DENG Guibo，ZHAO Yuting，et al.Economic dispatch of power system with wind power considering the complementarity of low-carbon characteristics of source side and load side J.Proceedings of the CSEE，2021，41（14）：4799-4815.9 韩莹，于三川，李荦一，等.计及阶梯式碳交易的风光氢储微电网低碳经济配置方法 J.高电压技术，2022，48（7）：2523-2533.HAN Ying，YU Sanchuan，LI Luoyi，et al.Low-carbon and economic configuration method for solar hydrogen storage microgrid including stepped carbon tradingJ.High Voltage Engineering，2022，48（7）：2523-2533.10 孙惠娟，吴鹏，彭春华，等.基于源荷协同降碳的综合能源系统分布鲁棒经济调度 J.电网技术，2023，47（4）：1589-1603.SUN Huijuan，WU Peng，PENG Chunhua，et al.Distributionally robust economic scheduling for integrated energy systems based on source-load cooperative carbon reductionJ.Power System Technology，2023，47（4）：1589-1603.11 International Organization for Standardization.Environmental management：life cycle assessment；Principles and Framework M.ISO，2006.12 龙芸.燃煤电厂CO2排放计算模型与方法研究 D.重庆：重庆大学，2016.LONG Yun.Study on calculation model and method of carbon dioxide emissions for coal-fired power plantD.Chongqing：Chongqing University，2016.13 黄景光，熊华健，李振兴，等.基于生命周期法和碳权交易的综合能源系统低碳经济调度 J.电测与仪表，2022，59（3）：82-91.HUANG Jingguang，XIONG Huajian，LI Zhenxing，et al.Low-carbon economic dispatch of integrated energy system based on life cycle method and carbon tradingJ.Electrical Measurement&Instrumentation，2022，59（3）：82-91.14 CHAU C，LEUNG T，NG W.A review on Life Cycle Assessment，Life Cycle Energy Assessment and Life Cycle Carbon Emissions Assessment on buildingsJ.Applied energy，2015，143：395-413.15 周天睿，康重庆，徐乾耀，等.电力系统碳排放流分析理论初探 J.电力系统自动化，2012，36（7）：38-43.ZHOU Tianrui，KANG Chongqing，XU Qianyao，et al.Preliminary theoretical investigation on power system carbon emission flow J.Automation of Electric Power Systems，2012，36（7）：38-43.16 周天睿，康重庆，徐乾耀，等.电力系统碳排放流的计算方法初探 J.电力系统自动化，2012，36（11）：44-49.ZHOU Tianrui，KANG Chongqing，XU Qianyao，et al.Preliminary investigation on a method for carbon emission flow calculation of power system J.Automation of Electric Power Systems，2012，36（11）：44-49.17 康重庆，杜尔顺，李姚旺，等.新型电力系统的“碳视角”：科学问题与研究框架 J.电网技术，2022，46（3）：821-833.KANG Chongqing，DU Ershun，LI Yaowang，et al.Key scientific problems and research framework for carbon perspective research of new power systems J.Power System Technology，2022，46（3）：821-833.18 刘哲远，邢海军，程浩忠，等.考虑碳排放流及需求响应的综合能源系统双层优化调度 J.高电压技术，2023，49（1）：169-178.LIU Zheyuan，XING Haijun，CHENG Haozhong，et al.Bi-level optimal scheduling of integrated energy system considering carbon emission flow and demand responseJ.High Voltage Engineering，2023，49（1）：169-178.19 岳建坤.基于潮流追踪的输电费用分摊方法研究 D.北京：华北电力大学，2013.YUE Jiankun.Research on allocation method of transmission costs based on power flow tracing D.Beijing：North China Electric Power University，2013.20 CHENG Y H，ZHANG N，WANG Y，et al.Modeling carbon emission flow in multiple energy systemsJ.IEEE Transactions on Smart Grid，2018，10（4）：3562-3574.收稿日期：2022-12-06；修回日期：2023-01-15作者简介：刘子华（1988），男，工学学士，工程师，主要从事分布式电源、储能、非统调电厂并网等方面的研究工作。（本文编辑：方明霞）72