CN117371650B - 考虑负荷侧电能替代的配电网精准碳计量方法及系统 - Google Patents
考虑负荷侧电能替代的配电网精准碳计量方法及系统 Download PDFInfo
- Publication number
- CN117371650B CN117371650B CN202311306251.8A CN202311306251A CN117371650B CN 117371650 B CN117371650 B CN 117371650B CN 202311306251 A CN202311306251 A CN 202311306251A CN 117371650 B CN117371650 B CN 117371650B
- Authority
- CN
- China
- Prior art keywords
- carbon
- node
- branch
- load
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 321
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 314
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000006467 substitution reaction Methods 0.000 title claims abstract description 32
- 239000011159 matrix material Substances 0.000 claims abstract description 114
- 238000004364 calculation method Methods 0.000 claims abstract description 40
- 230000009467 reduction Effects 0.000 claims abstract description 25
- 230000004907 flux Effects 0.000 claims abstract description 19
- 238000002347 injection Methods 0.000 claims abstract description 17
- 239000007924 injection Substances 0.000 claims abstract description 17
- 230000008901 benefit Effects 0.000 claims abstract description 13
- 238000012937 correction Methods 0.000 claims description 16
- 239000002803 fossil fuel Substances 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 13
- 238000010438 heat treatment Methods 0.000 claims description 11
- 239000000446 fuel Substances 0.000 claims description 7
- 230000005611 electricity Effects 0.000 claims description 6
- 238000010257 thawing Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 230000036316 preload Effects 0.000 claims 1
- 238000004134 energy conservation Methods 0.000 abstract description 6
- 238000005265 energy consumption Methods 0.000 abstract description 6
- 239000003245 coal Substances 0.000 description 12
- 238000010248 power generation Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000006386 neutralization reaction Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- General Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Primary Health Care (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Development Economics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种考虑负荷侧电能替代的配电网精准碳计量方法及系统,包括:获取当前时刻系统相关参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量;基于碳排放流理论,计算节点碳势和支路碳流密度;计算各节点负荷、支路及损耗对应的碳流率和等效碳排放量;负荷侧电能替代后计算减碳效益并对碳势进行修正,根据修正值计算电能替代后下一时刻的碳势并提出针对负荷侧用能的节能降碳建议。本发明计算考虑用户侧电能替代的配电网各支路和负荷侧的碳排放量,有利于配电网的低碳运行和用户侧的节能降碳。
Description
技术领域
本发明属于配电网碳排放核算领域,具体涉及一种考虑负荷侧电能替代的配电网精准碳计量方法及系统。
背景技术
为实现节能降碳目标,能源是主战场,电力是主阵地。作为碳达峰、碳中和的重要手段,精确的碳计量对于高效地管理碳排放是至关重要的。
电力系统全环节碳排放的精准碳计量也成为了低碳发展背景下电力系统的关键科学问题之一。准确、全面的电力系统碳计量与分析方法是量化电力相关碳排放数据,理清电力行业碳排放现状与趋势,挖掘电力系统全环节碳减排潜力,引导电力用户互动减碳,是支撑电碳耦合市场建设的关键技术。
对电力行业而言,间接碳排放将电力系统碳排放责任从源侧传导至荷侧。电力系统碳排放流理论为电力系统的间接碳计量提供了新的研究思路。碳随二次能源流动而转移,虚拟的碳排放流依附于电力潮流存在且随潮流沿电网移动,直至用电终端。通过分析虚拟碳排放流,将发电过程中产生的碳排放传导至用户,从而研究用户用电行差异对碳排放的影响,以推进减碳政策的实施和引导用户的减碳行为。本发明在传统碳流理论的基础上充分考虑负荷侧电能替代的行为,提出一种考虑负荷侧电能替代的配电网精准碳计量方法,旨在实现配电网的低碳运行和用户侧的节能降碳。
发明内容
针对配电网碳计量缺乏用户侧参与的问题,本发明提供一种考虑负荷侧电能替代的配电网精准碳计量方法及系统,能够实现配电网的低碳运行和用户侧的节能降碳。
实现本发明目的的技术解决方案为:一种考虑负荷侧电能替代的配电网精准碳计量方法,包括以下步骤:
步骤1:获取当前时刻系统中各节点负荷(包括有功和无功负荷)、发电机有功和无功出力、网络拓扑结构以及电阻、电抗等相关技术参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;
步骤2:对步骤1中的结果进行处理,构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量;
步骤3:基于碳排放流理论,计算荷侧节点碳势和网测支路碳流密度;
步骤4:对步骤3中的结果进行处理,计算各节点负荷、支路及损耗对应的碳流率即等效碳排放量;
步骤5:负荷侧电能替代后计算减碳效益并对碳势进行修正,根据修正结果计算碳流率并提出针对负荷侧用能的节能降碳建议。
进一步地,步骤1)中电力系统的碳排放流与功率流直接相关,在碳流计算之前,需事先确定网络中的功率分布。对于一个含n个节点的电力系统,其节点有功和无功潮流平衡方程如下式:
式中,PGi、QGi分别为接入节点i的发电机有功和无功输出;PLi、QLi分别为接入节点i的有功和无功负荷;Gij、Bij分别为支路i-j电导和电纳的负值;Ui、Uj分别为节点i和节点j的电压幅值;δij为节点i和节点j的电压角度差,即δij=δi-δj,其中δi和δj分别为节点i和节点j的电压角度。
根据上述潮流计算结果,各支路传输有功和损耗计算公式如下:
Pij=ViVj(Gijcosδij+Bijsinδij)
式中,Pij为流经支路i-j的有功功率,其正方向为节点i到节点j;为支路i-j的有功损耗。
进一步地,步骤2)中所述物理量包括以下:
①支路潮流分布矩阵为N阶方阵,用PB=(PBij)N×N表示。定义该矩阵的目的是为了描述电力系统的有功潮流分布,从电力网络层面给定碳排放流分布的边界条件。该矩阵既包含电力网络的拓扑结构信息,又包含系统稳态有功潮流的分布信息。支路潮流分布矩阵中的元素具体定义如下:
若节点i与节点j间有支路相连,且经此支路从节点i到节点j流过的正向有功潮流为p,则PBij=p,PBji=0;若流经该支路的有功潮流p为反向潮流,则PBij=0,PBji=p;其他情况下PBij=PBji=0。特别地,对所有对角元素,有PBii=0。
②机组注入分布矩阵为K×N阶矩阵,用PG=(PGkj)K×N表示。定义该矩阵的目的是为了描述所有发电机组与电力系统的连接关系以及机组向系统中注入的有功功率,同时也是便于描述系统中发电机组产生碳排放流的边界条件。矩阵中的元素具体定义如下:
若第k台发电机组接入节点j,且从第k个含有发电机的节点注入节点j的有功潮流为p,则PGkj=p,否则PGkj=0。
③节点有功通量矩阵为N阶对角阵,用PN=(PNij)N×N表示。根据基尔霍夫电流定律,任何时刻对任意节点,所有流入、流出该节点的支路电流的绝对值相等,代数和恒等于0。由此,在潮流分析中,任意节点的净注入功率均为0。但在碳流计算中,节点碳势只受注入潮流的影响,从节点流出的潮流对节点碳势不产生影响。因此,相比流经节点电流和潮流的代数和,碳流计算更关注考虑潮流方向下流入节点有功潮流的“绝对量”,称之为节点有功通量。在碳流计算中,将利用此概念来描述系统中发电机组对节点以及节点对节点碳势的贡献。节点有功通量矩阵的元素具体定义如下:
对节点i,令I+表示有潮流流入节点i的支路集合,pBs为支路s的有功功率,则有:
式中,pGi为接入节点i的发电机组出力,若该节点无发电机组或发电机组出力为0,则pGi=0。该矩阵中所有非对角元素PNij=0。
④负荷分布矩阵为M×N阶矩阵,用PL=(PLmj)M×N表示。定义该矩阵的目的是描述所有用电负荷与电力系统的连接关系以及有功负荷量,以描述系统中电力用户消费碳排放流的边界条件,矩阵中的元素具体定义如下:
若节点j是第m个存在负荷的节点,且有功负荷为p,则则PLmj=p,否则PLmj=0。
⑤不同发电机组具有不同的碳排放特性,在碳流计算中为已知条件,可组成系统的发电机组碳排放强度向量。设第k台发电机组的碳排放强度为eGk,则发电机组碳排放强度向量可表示为:
EG=[eG1 eG2 eG3 … eGk]T
进一步地,步骤3)中碳势的物理意义都可以定义为在该节点消耗单位电量所造成的等效于发电侧的碳排放值,其量纲一般为kgCO2/(kWh)。由节点碳势的定义,可得系统中t时刻节点j的碳势eNj:
式中,J+表示有潮流流入节点j的支路集合;pBs为支路s的有功功率;pGj为接入节点j的发电机组出力,若该节点无发电机组或发电机组出力为0,则pGj=0;eGj为发电机组的碳排放强度;ρs为支路s的碳流密度。
根据碳排放流的性质,支路碳流密度ρs可由支路始端节点碳势替代,将上式改成以下矩阵形式:
式中:为N维单位行向量,其中第j个元素为1;
根据节点有功通量矩阵的定义,可得:
由以上两式可得:
由于PN矩阵为对角阵,将上式扩充至全系统维度,可得:
整理后可得系统所有节点的碳势计算公式为:
式中,为t时刻节点碳势向量,/>
进一步地,步骤4)中计算得到节点碳势向量后,可进一步得到系统各个支路的碳流率和负荷碳流率。
①支路碳流率分布矩阵为N阶方阵,用表示。支路碳流率分布矩阵元素定义与支路潮流分布矩阵相似。若节点i与节点间有支路相连,且经此支路从节点i到节点j流过正向碳流率为R,则/>若流经该支路的碳流率R为反向,则其他情况下/>特别地,对所有对角元素,有/>
式中,为支路潮流分布矩阵;/>为支路损耗分布矩阵即损耗碳流率,由支路i-j的有功损耗/>组成;/>为节点碳势向量。
②负荷碳流率为M维向量,用表示。物理意义为发电侧为供应节点负荷每单位时间产生的碳排放量。计算得到节点碳势向量后,节点负荷的用电碳排放强度与该节点碳势相等,结合负荷分布矩阵,可得所有负荷对应的碳流率。对第m个存在负荷的节点,与其负荷对应的碳流率为/>则负荷碳流率向量可表示为:
式中,为负荷分布矩阵;/>为节点碳势向量。
进一步地,步骤5)中进行负荷侧电能替代,以空气源热泵为例进行热负荷电能替代。
基于小时级的热负荷需求数据,计算传统化石能源供能所产生的碳排放。假定实施电能替代之前,节点热负荷主要由燃煤锅炉产生的热力满足,此时的碳排放量主要取决于化石燃料的消耗量,二者之间的关系见下式:
式中,NCVj为第j种化石燃料的平均低位发热量,FCj,t为第j种化石燃料t时段的消耗量,CFj为第j种化石燃料的碳排放因子,j为化石燃料类型;
计及管道传输的热损失,室外供热管网输送效率β,已知t时段节点的热负荷需求为Hl,t,则t时段所需供热量为:
根据能量守恒,由节点热力需求Ht推导出锅炉的燃煤量为:
式中,FCcoal,t为燃煤锅炉t时段的耗煤量,NCVcoal为标准煤的平均低位发热量,ηcoal为锅炉热效率,受锅炉型号、新旧程度的影响,ADfuel,t为t时段供热设备消耗的外购化石能源总量。
基于小时级的热负荷需求数据,计算实施电能替代后节点供热产生的碳排放量及其减碳效益。假定实施电能替代之后,节点热负荷主要由空气源热泵产生的热力满足。
根据能量守恒定律,由节点热力需求Ht推导出气源热泵的耗电量为:
式中,Pheat,t为t时段节点电能替代后用于热负荷的用电量,COPcor为机组实际制热性能系数,COPicy为忽略机组除霜影响后的机组制热性能系数,COPred为除霜过程中机组COP的衰减率。
此时的碳排放量主要取决于电能的消耗量,二者之间的关系见下式:
则节点的供热进行电能替代后的减碳效益为:
Crd,t=Cfuel,t-Cheat,t
电能替代之后,负荷侧用能效率提高,即消耗同样的电量用能效益更高,可等效于负荷侧相较于电能替代前碳排放量较少即减碳效益,则电能替代后负荷碳流率/>可表示为:
根据可得t时刻电能替代后的碳势/>表示为:
则t时刻电能替代前后碳势的修正值可表示为:
由于t+1时刻与t时刻潮流信息不相同,因而不能用修正后的碳势来表征下一时刻的状态,需要在t+1时刻的基础上,用电能替代后的修正值来表示t+1时刻修正后的碳势,按照碳势的计算公式可得修正前t+1时刻的碳势/>为:
则t+1时刻修正后的碳势为:
则电能替代后t+1时刻的碳流率为:
与现有技术相比,本发明的有益效果为:本发明负荷侧电能替代后计算减碳效益并对碳势进行修正,根据修正值计算电能替代后下一时刻的碳势并提出针对负荷侧用能的节能降碳建议,本发明能有效计算考虑用户侧电能替代的配电网各支路和负荷侧的碳排放量,有利于配电网的低碳运行和用户侧的节能降碳。
附图说明
图1是考虑电能替代的碳排放流计算流程图。
图2是电能替代后节点碳势对比曲线图。
图3是电能替代后节点碳流率对比曲线图。
图4是电能替代后节点碳排放量统计图。
具体实施方式
下面结合附图对本发明提出的考虑负荷侧电能替代的配电网精准碳计量方法做出详细说明。
本发明所涉及方法包括以下环节:获取当前时刻系统各节点信息并进行潮流计算、基于碳排放流理论对节点碳势和碳流率进行计算、对负荷侧进行电能替代并进行碳势修正等。
其中,获取当前时刻系统中各节点负荷(包括有功和无功负荷)、发电机有功和无功出力、网络拓扑结构以及电阻、电抗等相关技术参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;基于碳排放流理论,构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量,计算各节点负荷、支路及损耗对应的碳流率即等效碳排放量;负荷侧电能替代后计算减碳效益并对碳势进行修正,根据修正结果计算碳流率并提出针对负荷侧用能的节能降碳建议。
具体的,本发明提出考虑负荷侧电能替代的配电网精准碳计量方法,如图1所示,包括如下步骤:
1)获取当前时刻系统中各节点负荷(包括有功和无功负荷)、发电机有功和无功出力、网络拓扑结构以及电阻、电抗等相关技术参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;
电力系统的碳排放流与功率流直接相关,在碳流计算之前,需事先确定网络中的功率分布。对于一个含n个节点的电力系统,其节点有功和无功潮流平衡方程如下式:
式中,PGi、QGi分别为接入节点i的发电机有功和无功输出;PLi、QLi分别为接入节点i的有功和无功负荷;Gij、Bij分别为支路i-j电导和电纳的负值;Ui、Uj分别为节点i和节点j的电压幅值;δij为节点i和节点j的电压角度差,即δij=δi-δj,其中δi和δj分别为节点i和节点j的电压角度。
根据上述潮流计算结果,各支路传输有功和损耗计算公式如下:
Pij=ViVj(Gijcosδij+Bijsinδij)
式中,Pij为流经支路i-j的有功功率,其正方向为节点i到节点j;为支路i-j的有功损耗。
2)基于步骤1对步骤中的结果进行处理,构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量(即发电碳排放因子);
①支路潮流分布矩阵为N阶方阵,用PB=(PBij)N×N表示。定义该矩阵的目的是为了描述电力系统的有功潮流分布,从电力网络层面给定碳排放流分布的边界条件。该矩阵既包含电力网络的拓扑结构信息,又包含系统稳态有功潮流的分布信息。支路潮流分布矩阵中的元素具体定义如下:
若节点i与节点j间有支路相连,且经此支路从节点i到节点j流过的正向有功潮流为p,则PBij=p,PBji=0;若流经该支路的有功潮流p为反向潮流,则PBij=0,PBji=p;其他情况下PBij=PBji=0。特别地,对所有对角元素,有PBii=0。
②机组注入分布矩阵为K×N阶矩阵,用PG=(PGkj)K×N表示。定义该矩阵的目的是为了描述所有发电机组与电力系统的连接关系以及机组向系统中注入的有功功率,同时也是便于描述系统中发电机组产生碳排放流的边界条件。矩阵中的元素具体定义如下:
若第k台发电机组接入节点j,且从第k个含有发电机的节点注入节点j的有功潮流为p,则PGkj=p,否则PGkj=0。
③节点有功通量矩阵为N阶对角阵,用PN=(PNij)N×N表示。根据基尔霍夫电流定律,任何时刻对任意节点,所有流入、流出该节点的支路电流的绝对值相等,代数和恒等于0。由此,在潮流分析中,任意节点的净注入功率均为0。但在碳流计算中,节点碳势只受注入潮流的影响,从节点流出的潮流对节点碳势不产生影响。因此,相比流经节点电流和潮流的代数和,碳流计算更关注考虑潮流方向下流入节点有功潮流的“绝对量”,称之为节点有功通量。在碳流计算中,将利用此概念来描述系统中发电机组对节点以及节点对节点碳势的贡献。节点有功通量矩阵的元素具体定义如下:
对节点i,令I+表示有潮流流入节点i的支路集合,pBs为支路s的有功功率,则有:
式中,pGi为接入节点i的发电机组出力,若该节点无发电机组或发电机组出力为0,则pGi=0。该矩阵中所有非对角元素PNij=0。
④负荷分布矩阵为M×N阶矩阵,用PL=(PLmj)M×N表示。定义该矩阵的目的是描述所有用电负荷与电力系统的连接关系以及有功负荷量,以描述系统中电力用户消费碳排放流的边界条件,矩阵中的元素具体定义如下:
若节点j是第m个存在负荷的节点,且有功负荷为p,则则PLmj=p,否则PLmj=0。
⑤不同发电机组具有不同的碳排放特性,在碳流计算中为已知条件,可组成系统的发电机组碳排放强度向量。设第k台发电机组的碳排放强度为eGk,则发电机组碳排放强度向量可表示为:
EG=[eG1 eG2 eG3 … eGK]T
3)基于碳排放流理论,计算节点碳势和支路碳流密度;
碳势的物理意义都可以定义为在该节点消耗单位电量所造成的等效于发电侧的碳排放值,其量纲一般为kgCO2/(kWh)。由节点碳势的定义,可得系统中t时刻节点j的碳势eNj:
式中,J+表示有潮流流入节点j的支路集合;pBs为支路s的有功功率;pGj为接入节点j的发电机组出力,若该节点无发电机组或发电机组出力为0,则pGj=0;eGj为发电机组的碳排放强度;ρs为支路s的碳流密度。
根据碳排放流的性质,支路碳流密度ρs可由支路始端节点碳势替代,将上式改成以下矩阵形式:
式中:为N维单位行向量,其中第j个元素为1;
根据节点有功通量矩阵的定义,可得:
由以上两式可得:
由于PN矩阵为对角阵,将上式扩充至全系统维度,可得:
整理后可得系统所有节点的碳势计算公式为:
式中,为t时刻节点碳势向量,/>
4)对步骤3中的结果进行处理,计算各节点负荷、支路及损耗对应的碳流率即等效碳排放量;
计算得到节点碳势向量后,可进一步得到系统各个支路的碳流率和负荷碳流率。
①支路碳流率分布矩阵为N阶方阵,用表示。支路碳流率分布矩阵元素定义与支路潮流分布矩阵相似。若节点i与节点间有支路相连,且经此支路从节点i到节点j流过正向碳流率为R,则/>若流经该支路的碳流率R为反向,则其他情况下/>特别地,对所有对角元素,有/>
式中,为支路潮流分布矩阵;/>为支路损耗分布矩阵;/>为节点碳势向量。
②负荷碳流率为M维向量,用表示。物理意义为发电侧为供应节点负荷每单位时间产生的碳排放量。计算得到节点碳势向量后,节点负荷的用电碳排放强度与该节点碳势相等,结合负荷分布矩阵,可得所有负荷对应的碳流率。对第m个存在负荷的节点,与其负荷对应的碳流率为/>则负荷碳流率向量可表示为:
式中,为负荷分布矩阵;/>为节点碳势向量。
5)负荷侧电能替代后计算减碳效益并对碳势进行修正,根据修正结果计算碳流率;
进行负荷侧电能替代,以空气源热泵为例进行热负荷电能替代。
基于小时级的热负荷需求数据,计算传统化石能源供能所产生的碳排放。假定实施电能替代之前,节点热负荷主要由燃煤锅炉产生的热力满足,此时的碳排放量主要取决于化石燃料的消耗量,二者之间的关系见下式:
式中,NCVj为第j种化石燃料的平均低位发热量,FCj,t为第j种化石燃料t时段的消耗量,CFj为第j种化石燃料的碳排放因子,j为化石燃料类型;
计及管道传输的热损失,室外供热管网输送效率β,已知t时段节点的热负荷需求为Hl,t,则t时段所需供热量为:
根据能量守恒,由节点热力需求Ht推导出锅炉的燃煤量为:
式中,FCcoal,t为燃煤锅炉t时段的耗煤量,NCVcoal为标准煤的平均低位发热量,ηcoal为锅炉热效率,受锅炉型号、新旧程度的影响,ADfuel,t为t时段供热设备消耗的外购化石能源总量。
基于小时级的热负荷需求数据,计算实施电能替代后节点供热产生的碳排放量及其减碳效益。假定实施电能替代之后,节点热负荷主要由空气源热泵产生的热力满足。
根据能量守恒定律,由节点热力需求Ht推导出气源热泵的耗电量为:
式中,Pheat,t为t时段节点电能替代后用于热负荷的用电量,COPcor为机组实际制热性能系数,COPicy为忽略机组除霜影响后的机组制热性能系数,COPred为除霜过程中机组COP的衰减率。
此时的碳排放量主要取决于电能的消耗量,二者之间的关系见下式:
则节点的供热进行电能替代后的减碳效益为:
Crd,t=Cfuel,t-Cheat,t
电能替代之后,负荷侧用能效率提高,即消耗同样的电量用能效益更高,可等效于负荷侧相较于电能替代前碳排放量较少即减碳效益,则电能替代后负荷碳流率/>可表示为:
根据可得t时刻电能替代后的碳势/>表示为:
则t时刻电能替代前后碳势的修正值可表示为:
由于t+1时刻与t时刻潮流信息不相同,因而不能用修正后的碳势来表征下一时刻的状态,需要在t+1时刻的基础上,用电能替代后的修正值来表示t+1时刻修正后的碳势,按照碳势的计算公式可得修正前t+1时刻的碳势/>为:
则t+1时刻修正后的碳势为:
则电能替代后t+1时刻的碳流率为:
一种考虑负荷侧电能替代的配电网精准碳计量系统,包括各支路有功功率和损耗计算单元、碳排放因子确定单元、节点碳势和支路碳流密度计算单元、碳排放量计算单元和修正单元;其中:
所述各支路有功功率和损耗计算单元用于获取当前时刻系统参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;
所述碳排放因子确定单元用于构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量;
所述节点碳势和支路碳流密度计算单元用于通过通过碳排放流理论计算节点碳势和支路碳流密度;
所述碳排放量计算单元用于计算各节点负荷、支路及损耗对应的碳流率,即等效碳排放量;
所述修正单元通过负荷侧电能替代对碳势进行修正。
该系统包含上述方法的全部技术特征,在此不再累述。
为验证本方法的有效性,结合上述步骤,以IEEE14母线系统进行算例分析。执行算例仿真的软件MATLAB_R2021b。所使用的仿真平台处理器为AMD Ryzen 74800U,内存为16GB,操作系统为64位Windows 11。对节点4和13电能替代后对碳势进行修正,结果对比图如图2所示,并在此基础上对重新计算碳流率并以燃料排放因子法为对照,结果对比图如图3所示,电能替代后节点碳排放量统计如图4所示,由图2-图4可见,本方法对负荷侧进行电能替代以在传统碳排放流的基础上对碳势进行修正,实现了配网的低碳运行和用户侧的节能减碳。
Claims (10)
1.一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,包括以下步骤:
获取当前时刻系统参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;
基于各支路有功功率和损耗,构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量;
基于支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵、负荷分布矩阵和发电机组碳排放强度向量,通过碳排放流理论计算节点碳势和支路碳流密度;
基于节点碳势和支路碳流密度,计算各节点负荷、支路及损耗对应的碳流率,即等效碳排放量;
基于各节点负荷、支路及损耗对应的碳流率,通过负荷侧电能替代对碳势进行修正,进而获取配电网精准碳计量;
所述计算节点碳势和支路碳流密度具体包括:
确定系统中t时刻节点j的碳势为:
其中,J+表示有潮流流入节点j的支路集合;pBs为支路s的有功功率;pGj为接入节点j的发电机组出力,若该节点无发电机组或发电机组出力为0,则pGj=0;eGj为发电机组的碳排放强度;ρs为支路s的碳流密度;
所述支路碳流密度由支路始端节点碳势替代为:
式中:为N维单位行向量,其中第j个元素为1,K为发电机组,PBvj为路潮流分布矩阵中第v行第j列的元素,PGtj为机组注入分布矩阵中第t行第j列的元素;
则所述系统所有节点的碳势为:
式中,为t时刻节点碳势向量,/> 和EG分别为t时刻的节点有功通量矩阵、支路潮流分布矩阵、机组注入分布矩阵和发电机组碳排放强度向量;
所述计算各节点负荷碳流率包括:负荷碳流率为M维向量,用 表示,计算得到节点碳势向量后,节点负荷的用电碳排放强度与该节点碳势相等,结合负荷分布矩阵,可得所有负荷对应的碳流率,令第m个存在负荷的节点与其负荷对应的碳流率为则负荷碳流率向量可表示为:
式中,为负荷分布矩阵;/>为节点碳势向量;
所述计算支路碳流率包括:支路碳流率分布矩阵为N阶方阵,用 表示,支路碳流率分布矩阵元素定义为:若节点i与节点间有支路相连,且经此支路从节点i到节点j流过正向碳流率为R,则/>若流经该支路的碳流率R为反向,则其他情况下/> 对所有对角元素,有/>则有:
式中,为支路潮流分布矩阵;/>为支路损耗分布矩阵即损耗碳流率,由支路i-j的有功损耗/>组成;/>为节点碳势向量;
所述负荷侧电能替代为对空气源热泵进行热负荷电能替代,所述对碳势进行修正具体包括:
基于化石燃料的消耗量确定碳排放量Cfuel,t为:
式中,NCVj为第j种化石燃料的平均低位发热量,FCj,t为第j种化石燃料t时段的消耗量,CFj为第j种化石燃料的碳排放因子,j为化石燃料类型;
令室外供热管网输送效率为β,已知t时段节点的热负荷需求为Hl,t,则t时段所需供热量为:
根据能量守恒定律,由节点供热量需求Ht推导出气源热泵的耗电量为:
式中,Pheat,t为t时段节点电能替代后用于热负荷的用电量,COPcor为机组实际制热性能系数,COPicy为忽略机组除霜影响后的机组制热性能系数,COPred为除霜过程中机组COP的衰减率;
此时的碳排放量取决于电能的消耗量,二者之间的关系见下式:
则节点的供热进行电能替代后的减碳效益为:
Crd,t=Cfuel,t-Cheat,t
电能替代之后,可等效于负荷侧相较于电能替代前负荷碳流率较少/>即减碳效益,则电能替代后负荷碳流率/>表示为:
根据可得t时刻电能替代后的碳势/>表示为:
则t时刻电能替代前后碳势的修正值表示为:
确定修正前t+1时刻的碳势为:
则t+1时刻修正后的碳势为:
则电能替代后t+1时刻的碳流率为:
2.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述系统参数包括各节点有功和无功负荷、发电机有功和无功出力、网络拓扑结构以及电阻和电抗。
3.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗具体包括:
确定网络中的功率分布,对于一个含n个节点的电力系统,其节点有功和无功潮流平衡方程为:
式中,PGi、QGi分别为接入节点i的发电机有功和无功输出;PLi、QLi分别为接入节点i的有功和无功负荷;Gij、Bij分别为支路i-j电导和电纳的负值;Ui、Uj分别为节点i和节点j的电压幅值;δij为节点i和节点j的电压角度差,即δij=δi-δj,其中δi和δj分别为节点i和节点j的电压角度;
基于节点有功和无功潮流平衡方程,确定各支路有功功率和损耗。
4.根据权利要求3所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述各支路有功功率和损耗为:
Pij=ViVj(Gijcosδij+Bijsinδij)
式中,Pij为流经支路i-j的有功功率,其正方向为节点i到节点j;为支路i-j的有功损耗。
5.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述支路潮流分布矩阵为N阶方阵,用PB=(PBij)N×N表示,支路潮流分布矩阵中的元素具体定义为:若节点i与节点j间有支路相连,且经此支路从节点i到节点j流过的正向有功潮流为p,则PBij=p,PBji=0;若流经该支路的有功潮流p为反向潮流,则PBij=0,PBji=p;其他情况下PBij=PBji=0,对所有对角元素,有PBii=0。
6.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述节点有功通量矩阵为N阶对角阵,用PN=(PNij)N×N表示,节点有功通量矩阵的元素具体定义为:对节点i,令I+表示有潮流流入节点i的支路集合,pBs为支路s的有功功率,则有:
式中,pGi为接入节点i的发电机组出力,若该节点无发电机组或发电机组出力为0,则pGi=0,该矩阵中所有非对角元素PNij=0。
7.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述机组注入分布矩阵为K×N阶矩阵,用PG=(PGkj)K×N表示,元素具体定义为:若第k台发电机组接入节点j,且从第k个含有发电机的节点注入节点j的有功潮流为p,则PGkj=p,否则PGkJ=0。
8.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述负荷分布矩阵为M×N阶矩阵,用PL=(PLmj)M×N表示,该矩阵中的元素具体定义为:若节点j是第m个存在负荷的节点,且有功负荷为p,则PLmj=p,否则PLmj=0。
9.根据权利要求1所述的一种考虑负荷侧电能替代的配电网精准碳计量方法,其特征在于,所述发电机组碳排放强度向量为:
EG=[eG1,eG2,eG3,…,eGk]T
其中,eGk为k台发电机组的碳排放强度。
10.一种实现权利要求1-9任一项所述考虑负荷侧电能替代的配电网精准碳计量方法的考虑负荷侧电能替代的配电网精准碳计量系统,其特征在于,包括各支路有功功率和损耗计算单元、碳排放因子确定单元、节点碳势和支路碳流密度计算单元、碳排放量计算单元和修正单元;其中:
所述各支路有功功率和损耗计算单元用于获取当前时刻系统参数,并用牛顿-拉夫逊法进行潮流计算,得到各支路有功功率和损耗;
所述碳排放因子确定单元用于构造支路潮流分布矩阵、节点有功通量矩阵、机组注入分布矩阵和负荷分布矩阵并根据不同发电机组特性构造发电机组碳排放强度向量;
所述节点碳势和支路碳流密度计算单元用于通过通过碳排放流理论计算节点碳势和支路碳流密度;
所述碳排放量计算单元用于计算各节点负荷、支路及损耗对应的碳流率,即等效碳排放量;
所述修正单元通过负荷侧电能替代对碳势进行修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311306251.8A CN117371650B (zh) | 2023-10-09 | 2023-10-09 | 考虑负荷侧电能替代的配电网精准碳计量方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311306251.8A CN117371650B (zh) | 2023-10-09 | 2023-10-09 | 考虑负荷侧电能替代的配电网精准碳计量方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117371650A CN117371650A (zh) | 2024-01-09 |
CN117371650B true CN117371650B (zh) | 2024-06-07 |
Family
ID=89390387
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311306251.8A Active CN117371650B (zh) | 2023-10-09 | 2023-10-09 | 考虑负荷侧电能替代的配电网精准碳计量方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117371650B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117764798B (zh) * | 2024-02-22 | 2024-05-24 | 福建省计量科学研究院(福建省眼镜质量检验站) | 一种用户碳表计量数据的校核方法及系统 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105375472A (zh) * | 2015-11-23 | 2016-03-02 | 国家电网公司 | 一种智能配电网低碳效益评估方法 |
CN106251095A (zh) * | 2016-09-06 | 2016-12-21 | 清华大学 | 一种电力系统碳排放实时计量的方法及碳表系统 |
CN107769192A (zh) * | 2017-09-18 | 2018-03-06 | 华北电力大学 | 一种供用电网络节点/支路碳流计算方法 |
CN113886767A (zh) * | 2021-08-30 | 2022-01-04 | 国网江苏省电力有限公司 | 一种电力系统碳排放流实时计算方法 |
CN114186805A (zh) * | 2021-11-23 | 2022-03-15 | 中国矿业大学 | 考虑实时电源成分及电能替代减碳效益的园区碳计量方法 |
CN114943471A (zh) * | 2022-06-14 | 2022-08-26 | 浙江华云电力工程设计咨询有限公司 | 一种电力系统低碳指标体系及综合评价方法 |
CN115375183A (zh) * | 2022-09-19 | 2022-11-22 | 浙江大学 | 考虑能流-碳流耦合约束的电力系统配网优化调度方法 |
CN115392528A (zh) * | 2022-06-06 | 2022-11-25 | 安徽大学 | 基于碳排放流理论的用电侧碳排放精细计量方法 |
CN116361603A (zh) * | 2023-03-28 | 2023-06-30 | 吉林省电力科学研究院有限公司 | 一种电力系统碳排放流计算方法 |
CN116388239A (zh) * | 2023-03-31 | 2023-07-04 | 广东电网有限责任公司珠海供电局 | 一种基于碳排放流理论的含有储能设备的节点碳势计算方法 |
CN116663823A (zh) * | 2023-05-25 | 2023-08-29 | 国网江苏省电力有限公司连云港供电分公司 | 融合精准碳数据的配电能源网格碳排放最优规划方法及系统 |
CN116796966A (zh) * | 2023-06-15 | 2023-09-22 | 贵州电网有限责任公司 | 一种基于双碳目标的能源网架规划方法 |
CN116823020A (zh) * | 2023-05-25 | 2023-09-29 | 国网河南省电力公司平顶山供电公司 | 一种考虑负荷侧减碳潜力的台区低碳运行综合评估方法 |
-
2023
- 2023-10-09 CN CN202311306251.8A patent/CN117371650B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105375472A (zh) * | 2015-11-23 | 2016-03-02 | 国家电网公司 | 一种智能配电网低碳效益评估方法 |
CN106251095A (zh) * | 2016-09-06 | 2016-12-21 | 清华大学 | 一种电力系统碳排放实时计量的方法及碳表系统 |
CN107769192A (zh) * | 2017-09-18 | 2018-03-06 | 华北电力大学 | 一种供用电网络节点/支路碳流计算方法 |
CN113886767A (zh) * | 2021-08-30 | 2022-01-04 | 国网江苏省电力有限公司 | 一种电力系统碳排放流实时计算方法 |
CN114186805A (zh) * | 2021-11-23 | 2022-03-15 | 中国矿业大学 | 考虑实时电源成分及电能替代减碳效益的园区碳计量方法 |
CN115392528A (zh) * | 2022-06-06 | 2022-11-25 | 安徽大学 | 基于碳排放流理论的用电侧碳排放精细计量方法 |
CN114943471A (zh) * | 2022-06-14 | 2022-08-26 | 浙江华云电力工程设计咨询有限公司 | 一种电力系统低碳指标体系及综合评价方法 |
CN115375183A (zh) * | 2022-09-19 | 2022-11-22 | 浙江大学 | 考虑能流-碳流耦合约束的电力系统配网优化调度方法 |
CN116361603A (zh) * | 2023-03-28 | 2023-06-30 | 吉林省电力科学研究院有限公司 | 一种电力系统碳排放流计算方法 |
CN116388239A (zh) * | 2023-03-31 | 2023-07-04 | 广东电网有限责任公司珠海供电局 | 一种基于碳排放流理论的含有储能设备的节点碳势计算方法 |
CN116663823A (zh) * | 2023-05-25 | 2023-08-29 | 国网江苏省电力有限公司连云港供电分公司 | 融合精准碳数据的配电能源网格碳排放最优规划方法及系统 |
CN116823020A (zh) * | 2023-05-25 | 2023-09-29 | 国网河南省电力公司平顶山供电公司 | 一种考虑负荷侧减碳潜力的台区低碳运行综合评估方法 |
CN116796966A (zh) * | 2023-06-15 | 2023-09-22 | 贵州电网有限责任公司 | 一种基于双碳目标的能源网架规划方法 |
Non-Patent Citations (4)
Title |
---|
Accounting for Carbon Emission Factors Measured by Electricity Users with Proportional Power Flow Tracking Considering Industry Characteristics;Wang PX等;2023 IEEE/IAS INDUSTRIAL AND COMMERCIAL POWER SYSTEM ASIA, I&CPS ASIA;20231231;1115-1120 * |
Optimal allocation of resources for yard crane network management to minimize carbon dioxide emissions;Peng Y等;JOURNAL OF CLEANER PRODUCTION;20160910;第131卷;649-658 * |
智能电网规划指标体系构建的量化分析;王宇坤;林其友;;电气应用;20190515(05);73-80 * |
考虑负荷侧电能替代的区域配电网 碳溯源方法;岳付昌等;广 东 电 力;20240225;第37卷(第2期);47-55 * |
Also Published As
Publication number | Publication date |
---|---|
CN117371650A (zh) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Distributionally robust day-ahead scheduling of park-level integrated energy system considering generalized energy storages | |
CN117371650B (zh) | 考虑负荷侧电能替代的配电网精准碳计量方法及系统 | |
CN111428351B (zh) | 基于前推回代法的电-热综合能源系统潮流计算方法 | |
CN111313429A (zh) | 一种用于综合能源系统的可靠性评估方法及系统 | |
CN109919450B (zh) | 解决综合智慧能源系统调度的博弈优化方法 | |
CN112258021B (zh) | 用于家庭燃料电池热电联供建筑的能效评估方法及系统 | |
CN116542370A (zh) | 一种考虑碳捕集与碳交易的园区低碳经济运行方法 | |
CN113554296A (zh) | 园区综合能源系统规划的多指标评价方法 | |
CN103337039A (zh) | 一种多目标多约束的中长期购电决策方法 | |
CN114186805A (zh) | 考虑实时电源成分及电能替代减碳效益的园区碳计量方法 | |
CN113205208A (zh) | 计及一次能源渗透率的综合能源系统㶲效率能效评估方法 | |
CN112883571A (zh) | 一种考虑动态多能需求响应的综合能源系统规划方法 | |
CN113158547B (zh) | 计及经济性和可靠性的区域综合能源系统优化配置方法 | |
CN117455104A (zh) | 一种基于电网全要素的碳流分析与计算方法 | |
CN111310310A (zh) | 一种用于量调节的热力系统静态潮流快速解耦计算方法 | |
CN110277785A (zh) | 基于连续多能流的电气耦合系统负荷裕度计算方法及系统 | |
CN116308874A (zh) | 基于电力系统碳排放流理论的区域碳排放流分析方法 | |
CN113131513B (zh) | 考虑碳排放的电、热、气转换优化运行配置方法和存储介质 | |
CN115908047A (zh) | 一种综合能源系统多时间尺度运行方法及其应用 | |
CN112465228B (zh) | 一种用户侧综合能源系统优化配置方法 | |
CN112529251B (zh) | 一种天然气发电机组选址定容优化规划方法及系统 | |
CN115034453A (zh) | 考虑碳交易市场的综合能源系统灵活运行方法及计算平台 | |
Wei et al. | Low-carbon economic dispatch of integrated energy system considering carbon trading and integrated demand response with CCER | |
CN107465195A (zh) | 一种基于微电网结合潮流计算的最优潮流双层迭代方法 | |
Cao et al. | Low-Carbon Planning of Integrated Energy System Considering a Reward and Punishment Ladder-type Carbon Trading Mechanism |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |