CN117151926A - 一种电力系统负荷侧碳责任核算系统和方法及介质 - Google Patents
一种电力系统负荷侧碳责任核算系统和方法及介质 Download PDFInfo
- Publication number
- CN117151926A CN117151926A CN202311102357.6A CN202311102357A CN117151926A CN 117151926 A CN117151926 A CN 117151926A CN 202311102357 A CN202311102357 A CN 202311102357A CN 117151926 A CN117151926 A CN 117151926A
- Authority
- CN
- China
- Prior art keywords
- load
- power
- carbon
- responsibility
- node
- 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.)
- Pending
Links
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 237
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 225
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 104
- 230000005611 electricity Effects 0.000 claims abstract description 46
- 238000009826 distribution Methods 0.000 claims abstract description 34
- 238000012937 correction Methods 0.000 claims description 21
- 230000002457 bidirectional effect Effects 0.000 claims description 15
- 230000014509 gene expression Effects 0.000 claims description 12
- 238000004458 analytical method Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 11
- 238000010248 power generation Methods 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 3
- 239000007924 injection Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 238000011144 upstream manufacturing Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 7
- 230000009467 reduction Effects 0.000 description 7
- 238000011161 development Methods 0.000 description 4
- 230000006399 behavior Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000002146 bilateral effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009472 formulation Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 230000001737 promoting effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000010792 warming Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- 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/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- Power Engineering (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种电力系统负荷侧碳责任核算系统和方法及介质,该系统的系统潮流计算模块用于得到电力系统的有功潮流分布及各支路网损;网损公平分摊模型模块用于根据有功潮流分布,计算系统各净负荷及各机组的净出力;二次潮流计算模块用于构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,确定系统各负荷成员的逐时碳排放情况;总碳排放量计算模块用于得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;最终碳责任计算模块用于对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。本发明考虑系统网损和用户的峰荷特性对系统碳排放量的影响,解决以往核算结果公平性不足的问题,提高了责任计算的准确性。
Description
技术领域
本发明涉及电力系统低碳优化技术领域,尤其涉及一种电力系统负荷侧碳责任核算系统和方法及介质。
背景技术
全球气候变暖问题日益突出,深刻影响着社会的发展与进步,电力行业作为重要的能源部门,其碳排放量一直高居不下,是气候变暖的重要原因之一。为了确保电力系统碳减排目标的顺利实现,缓解全球温室效应,其碳责任的划分显得尤为重要。碳责任分摊是指在对电力系统发电侧、电网侧以及负荷侧的碳排放量进行统计和分析的基础上,根据碳流理论,采用适宜的分摊方法划分各节点应承担的排放责任。因此,准确、合理、公平的碳责任计算结果,有助于精准锁定重点减排区域,促进生产端和运输端提升电力生产和传输调度能力,增强消费端的节约意识和合理用电能力,对促进各区域之间的碳减排合作和经济协调发展具有重要意义。
近年来,越来越多的研究者将责任划分的视角从发电侧转移到了负荷侧,依照碳排放流理论,将潮流与碳流相结合,摒弃了原先只计算发电侧碳责任的思想,有效克服了责任分摊公平性和准确性缺失的问题。但是,传统方法大多没有明确系统网损产生的碳排放对责任分摊结果的影响,无法精准核算该部分的责任量,不利于减排方案的制定。同时,忽略了各个负荷节点的峰谷用电行为对系统碳排放量的影响,使得分摊结果仅仅聚焦于系统潮流分布,无法建立用户的用电特征和碳责任的有机联系,严重影响责任分摊结果的合理性,因此,亟需新的责任核算方法助力行业低碳发展。
在当前新型电力系统大量接入可再生能源机组的背景下,系统的低碳经济性得到了大幅提升,安全稳定性却无法维持,导致夏、冬季节区域性断电的事件频发,同时,新能源大量替代常规机组,系统频率调节能力和电压调节能力显著下降。为确保系统平稳运行,提高系统运行鲁棒性,需留存部分燃煤机组作为备用电源,以满足用电高峰期仍然能够快速响应,保证系统运行稳定性的需求。在尽可能促进可再生能源消纳,保证系统平稳运行的前提下,负荷侧用电量剧增将导致发电侧部分常规机组发电量增加,进而增加系统的碳排放量,使系统的低碳属性进一步降低,不利于我国双碳目标的实现。
发明内容
本发明提供了一种电力系统负荷侧碳责任核算系统和方法及介质,该系统综合考虑各支路的网损归属及各用户的峰荷特征对碳责任分摊结果的影响,充分考虑“源随荷动”的因素,将负荷量激增特性引起的电源侧变化进而引起的碳排放变化考虑在内,很好地解决了以往责任分摊仅关注潮流断面分布情况的责任核算结果,分摊过程颗粒度较粗糙的问题。
本发明提供一种电力系统负荷侧碳责任核算系统,包括系统潮流计算模块、网损公平分摊模型模块、二次潮流计算模块、总碳排放量计算模块和最终碳责任计算模块:
所述系统潮流计算模块用于基于电力系统的网络拓扑结构、元件参数和源荷参量,完成系统潮流计算,得到电力系统的有功潮流分布及各支路网损;
所述网损公平分摊模型模块用于根据所述有功潮流分布,建立系统有功潮流顺逆流双向追踪模型,再将所述各支路网损等效为虚拟电源或负荷、各新能源机组等效为负的负荷,结合潮流追踪算法,建立电力系统的网损公平分摊模型,将系统网损公平分摊给源荷两侧,计算系统各净负荷及各机组的净出力;
所述二次潮流计算模块用于在完成网损公平分摊的基础上,构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,重新确定其功率分布,并结合机组的碳排放特性,计算得到系统碳流排放的流动状态和进行追踪分析,确定系统各负荷成员的逐时碳排放情况;
所述总碳排放量计算模块用于基于各负荷成员的逐时碳排放情况,得到各负荷成员典型日的总碳排放量,以此计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;
所述最终碳责任计算模块用于综合考虑不同用户的峰荷特点,采用综合负荷法,建立典型用电量修正模型,并对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。
进一步地,在所述系统潮流计算模块中,按照以下方式进行系统潮流计算:
解析给定的电力系统的网络拓扑结构,收集系统潮流计算所需的基础参数;
结合所述基础参数对电力系统进行建模仿真,并进行系统的潮流计算,得到电力系统的有功潮流分布。
进一步地,在所述网损公平分摊模型模块中,按照以下方式建立系统有功潮流顺逆流双向追踪模型:
在确定系统潮流分布的基础上,建立系统有功潮流顺逆流双向追踪模型,其表达式如下所示:
上式中[Ad]ij为系统功率的顺流跟踪矩阵;[Au]ij为系统功率的逆流跟踪矩阵;Si为节点i的流过功率;Sij为以i-j为正方向,线路ij流经过的功率;Di为节点i的下游节点集合;Ui为节点i的上游节点集合。
进一步地,在所述网损公平分摊模型模块中,按照以下方式计算系统各净负荷及各机组的净出力:
根据比例共享原则,计算得到各个发电机组和支路对各个负荷节点的贡献,各表达式如下所示:
式中:k为支路中增设的虚拟节点;节点i负荷、支路i-t流出功率可表示为各发电机的功率分量之和;SLk表示新增的虚拟k节点的负荷需求;SGi表示节点i电源出力;Sfi表示支路f-i的流入功率;n为原始网络中的节点数,m为存在负荷的节点个数;SGk是k节点虚拟成电源的网损,SLi是k节点虚拟成负荷的网损,Sit是线路i-t上的功率;式中所有包含字母X的参量为各机组、支路对节点i的贡献量,其中XGi为节点i电源出力对各负荷节点SLk的贡献量,Xfi为支路f-i的流入功率对各负荷节点SLk的贡献量,XLi为节点i负荷表示为各发电机的功率分量,Xit为支路i-t流出功率表示为各发电机的功率分量;
结合上述计算结果,构建网损的公平分摊模型,其中发电侧和用户侧各分摊50%的支路损耗,各负荷及电源承担的网损量的表达式分别如下所示:
再根据各支路网损和各负荷及电源承担的网损量,计算系统各净负荷及各机组的净出力,其表达式分别如下所示:
S'Lk=SLk+△SLk (9)
S'Gk=SGk-△SGk (10)
式中:S'Lk和S'Gk分别表示净负荷和各机组的净出力。
进一步地,在所述二次潮流计算模块中,按照以下方式确定系统各负荷成员的逐时碳排放情况:
在将网络损耗公平分摊给负荷侧和电源侧的基础上,通过直流潮流方程重新确定虚拟无损网络中各支路的有功分布及机组和支路对各负荷节点的贡献,考虑到“源随荷动的情况”,因此再采用潮流逆流追踪的方式,完成网损双向分摊和虚拟网络构建:
P'=Bθ' (11)
式中:θ'∈Rn×1为虚拟网络中各节点电压相角构成的矢量;P'∈Rn×1为各节点注入功率构成的矢量;节点k的注入功率为P'k=P'Gk-P'Lk,P'Gk和P'Lk分别为节点k的电源有功净出力和有功净负荷;B∈Rn×n是以支路电抗形成的节点导纳矩阵;X'Li和X'it分别为节点负荷以及支路功率中各发电机的有功贡献量;P'Li和P'it分别是虚拟无损网络的i节点负荷量和支路i-t的有功通量;
结合机组的碳排放强度,建立系统的碳流追踪与解析模型:
RGk=Re[△SGk]EGk (16)
式中:RLi为负荷节点i的碳流率,其含义为用户每小时用电所产生的等效于发电侧的碳排放量;Rit为支路i-t的碳流率,表示单位时间内随潮流通过的碳流量;RGk为机组k承担的网损碳流率;Re表示取复数ΔSGk的实部;EGk为各机组的碳排放强度向量。
进一步地,在所述总碳排放量计算模块中,按照以下方式计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度:
基于所述系统的碳流追踪与解析模型,计算得到各负荷成员典型日内的总碳排放量,进而得到月内的碳责任量,即:
yLi=30xLi (18)
式中:xLi为负荷i在典型日内的总碳排放量;yLi为负荷i的月内碳责任量;
基于负荷典型日内的总碳排放量结果,构建典型用电量及碳强度计算模型:
式中:Pi为负荷i的典型用电量;Ci为负荷i的典型碳强度。
进一步地,在所述最终碳责任计算模块中,采用如下方式得到得到最终碳责任计算结果:
基于用户的月内用电量数据,采用综合负荷法,构建修正比例计算模型:
式中:Wla为负荷成员a的修正比例;Lat、Lbt、…、Lit分别为各负荷成员典型月内的t时刻最大用电量;εt为用户的典型月最大用电负荷占各时刻用户的典型月最大用电负荷之和的比例;
综合考虑用户的峰荷特性对系统碳排放产生的影响,构建典型用电量修正模型:
式中:为系统负荷侧的修正比例均值;/>为系统负荷侧的典型用电量均值;Pci为典型用电量的修正值;
基于所述典型碳强度和修正后的典型用电量,生成反映负荷碳属性的最终碳责任计算结果:
YLi=30CiPci (23)
式中:YLi为负荷i的最终碳责任。
本发明还提供一种电力系统负荷侧碳责任核算方法,包括以下步骤:
基于电力系统的网络拓扑结构、元件参数和源荷参量,完成系统潮流计算,得到电力系统的有功潮流分布及各支路网损;
根据所述有功潮流分布,建立系统有功潮流顺逆流双向追踪模型,再将所述各支路网损等效为虚拟电源或负荷、各新能源机组等效为负的负荷,结合潮流追踪算法,建立电力系统的网损公平分摊模型,将系统网损公平分摊给源荷两侧,计算系统各净负荷及各机组的净出力;
在完成网损公平分摊的基础上,构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,重新确定其功率分布,并结合机组的碳排放特性,计算得到系统碳流排放的流动状态和进行追踪分析,确定系统各负荷成员的逐时碳排放情况;
基于各负荷成员的逐时碳排放情况,得到各负荷成员典型日的总碳排放量,以此计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;
综合考虑不同用户的峰荷特点,采用综合负荷法,建立典型用电量修正模型,并对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现所述电力系统负荷侧碳责任核算方法的步骤。
相比与现有技术,本发明的有益效果为:
1、本发明考虑系统网损和用户的峰荷特性对系统碳排放量的影响,在实现网损的“源-荷”双侧公平分摊的基础上,结合综合负荷法,对责任分摊结果进行修正,解决以往核算结果公平性不足的问题,提高了责任计算的准确性,为电力行业早日实现双碳目标提供参考,利于社会低碳水平的提高;
2、本发明建立的碳责任核算方法在公平分摊各支路网损给源荷两侧的基础上,充分考虑了不同用户的峰荷特性,对原始的碳责任进行了合理修正,所得结果可为电力系统低碳转型策略的制定指明方向,具有很好的应用价值。
附图说明
图1为本发明的一种电力系统负荷侧碳责任核算系统的结构框图;
图2为本发明的一种电力系统负荷侧碳责任核算方法的流程示意图;
图3为本发明实施例中具体的拓扑网络结构示意图;
图4为本发明实施例中修正前后各负荷成员责任量结果对比汇总图。
具体实施方式
以下结合附图1~4、表1~4和具体实施例对本发明作进一步的详细说明。
如图1所示,本发明提供一种电力系统负荷侧碳责任核算系统,包括系统潮流计算模块、网损公平分摊模型模块、二次潮流计算模块、总碳排放量计算模块和最终碳责任计算模块:
所述系统潮流计算模块用于基于电力系统的网络拓扑结构、元件参数和源荷参量,完成系统潮流计算,得到电力系统的有功潮流分布及各支路网损;
所述网损公平分摊模型模块用于根据所述有功潮流分布,建立系统有功潮流顺逆流双向追踪模型,再将所述各支路网损等效为虚拟电源或负荷、各新能源机组等效为负的负荷,结合潮流追踪算法,建立电力系统的网损公平分摊模型,将系统网损公平分摊给源荷两侧,计算系统各净负荷及各机组的净出力;
所述二次潮流计算模块用于在完成网损公平分摊的基础上,构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,重新确定其功率分布,并结合机组的碳排放特性,计算得到系统碳流排放的流动状态和进行追踪分析,确定系统各负荷成员的逐时碳排放情况;
所述总碳排放量计算模块用于基于各负荷成员的逐时碳排放情况,得到各负荷成员典型日的总碳排放量,以此计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;
所述最终碳责任计算模块用于综合考虑不同用户的峰荷特点,采用综合负荷法,建立典型用电量修正模型,并对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。
本发明考虑系统网损和用户的峰荷特性对系统碳排放量的影响,在实现网损的“源-荷”双侧公平分摊的基础上,结合综合负荷法,对责任分摊结果进行修正,解决以往核算结果公平性不足的问题,提高了责任计算的准确性,为电力行业早日实现双碳目标提供参考,利于社会低碳水平的提高。同时本发明建立的碳责任核算方法在公平分摊各支路网损给源荷两侧的基础上,充分考虑了不同用户的峰荷特性,对原始的碳责任进行了合理修正,所得结果可为电力系统低碳转型策略的制定指明方向,具有很好的应用价值。
本实施例中,在所述系统潮流计算模块中,按照以下方式进行系统潮流计算:
解析给定的电力系统的网络拓扑结构,收集系统潮流计算所需的基础参数;
结合所述基础参数对电力系统进行建模仿真,并进行系统的潮流计算,得到电力系统的有功潮流分布。
本实施例中,在所述网损公平分摊模型模块中,按照以下方式建立系统有功潮流顺逆流双向追踪模型:
在确定系统潮流分布的基础上,建立系统有功潮流顺逆流双向追踪模型,其表达式如下所示:
上式中[Ad]ij为系统功率的顺流跟踪矩阵;[Au]ij为系统功率的逆流跟踪矩阵;Si为节点i的流过功率;Sij为以i-j为正方向,线路ij流经过的功率;Di为节点i的下游节点集合;Ui为节点i的上游节点集合。
本实施例中,在所述网损公平分摊模型模块中,按照以下方式计算系统各净负荷及各机组的净出力:
根据比例共享原则,计算得到各个发电机组和支路对各个负荷节点的贡献,各表达式如下所示:
式中:k为支路中增设的虚拟节点,并将支路损耗等效为虚拟负荷或电源;将新能源机组等效为负的负荷,可以量化评估可再生能源消纳对系统的减碳贡献;节点i负荷、支路i-t流出功率可表示为各发电机的功率分量之和;SLk表示新增的虚拟k节点的负荷需求;SGi表示节点i电源出力;Sfi表示支路f-i的流入功率;n为原始网络中的节点数,m为存在负荷的节点个数;由于支路的损耗已经等效为虚拟电源或虚拟负荷,将功率SGi(i∈[n+1,n+m])相加即为网络总损耗的相反数,SLi(i∈[n+1,n+m]])相加即为网络总损耗,SGk是k节点虚拟成电源的网损,SLi是k节点虚拟成负荷的网损,Sit是线路i-t上的功率;式中所有包含字母X的参量为各机组、支路对节点i的贡献量,其中XGi为节点i电源出力对各负荷节点SLk的贡献量,Xfi为支路f-i的流入功率对各负荷节点SLk的贡献量,XLi为节点i负荷表示为各发电机的功率分量,Xit为支路i-t流出功率表示为各发电机的功率分量;
结合上述计算结果,构建网损的公平分摊模型,其中发电侧和用户侧各分摊50%的支路损耗,各负荷及电源承担的网损量的表达式分别如下所示:
再根据各支路网损和各负荷及电源承担的网损量,计算系统各净负荷及各机组的净出力,其表达式分别如下所示:
S'Lk=SLk+△SLk (9)
S'Gk=SGk-△SGk (10)
式中:S'Lk和S'Gk分别表示净负荷和各机组的净出力。
本实施例中,如图3所示,以某电力14节点系统为例进行计算和分析。基于图3给定电力系统的网络拓扑结构,收集的元件参数和源荷参量,根据上述公式(1)~(10),完成系统潮流计算。在本实施例中以IEEE14系统拓扑为基准,结合某地区11个负荷节点某典型日的用电情况,完成潮流计算并输出结果。其中,G1-G5为发电机组的编号,G1、G2和G8为常规机组,其碳排放强度分别为0.875tCO2/WM、0.525tCO2/WM和0.520tCO2/WM,G3和G4均为新能源机组;L2、L3至L14等为负荷成员编号;BUS1-BUS14等为各节点编号,各节点的基准电压均为110KV。
本实施例中,建立电力系统的网损公平分摊模型,得到典型日内虚拟无损网络的各负荷需求,汇总情况详见下表1。
表1.电力14节点系统的负荷信息汇总表
如上表1所示,各负荷成员均在不同程度上承担了系统网损。其中,L9所承担的网损量最大,为0.23MW;相反,负荷L11和L13承担的网损最小,均为0.02MW;系统的用电需求总量为194.30MW,网损公平分摊后的系统用电需求总量为195.34MW。由方法定义可知,负荷在网损公平分摊后的典型日内负荷需求总量为该负荷节点的典型用电量。
本实施例中,在所述二次潮流计算模块中,按照以下方式确定系统各负荷成员的逐时碳排放情况:
在将网络损耗公平分摊给负荷侧和电源侧的基础上,通过直流潮流方程重新确定虚拟无损网络中各支路的有功分布及机组和支路对各负荷节点的贡献,考虑到“源随荷动的情况”,因此再采用潮流逆流追踪的方式,完成网损双向分摊和虚拟网络构建:
P'=Bθ' (11)
式中:θ'∈Rn×1为虚拟网络中各节点电压相角构成的矢量;P'∈Rn×1为各节点注入功率构成的矢量;节点k的注入功率为P'k=P'Gk-P'Lk,P'Gk和P'Lk分别为节点k的电源有功净出力和有功净负荷;B∈Rn×n是以支路电抗形成的节点导纳矩阵;X'Li和X'it分别为节点负荷以及支路功率中各发电机的有功贡献量;P'Li和P'it分别是虚拟无损网络的i节点负荷量和支路i-t的有功通量。
结合机组的碳排放强度,建立系统的碳流追踪与解析模型:
RGk=Re[△SGk]EGk (16)
式中:RLi为负荷节点i的碳流率,其含义为用户每小时用电所产生的等效于发电侧的碳排放量;Rit为支路i-t的碳流率,表示单位时间内随潮流通过的碳流量;RGk为机组k承担的网损碳流率;Re表示取复数ΔSGk的实部;EGk为各机组的碳排放强度向量。
本实施例中,在所述总碳排放量计算模块中,按照以下方式计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度:
基于所述系统的碳流追踪与解析模型,计算得到各负荷成员典型日内的总碳排放量,进而得到月内的碳责任量,即:
yLi=30xLi (18)
式中:xLi为负荷i在典型日内的总碳排放量;yLi为负荷i的月内碳责任量。
基于负荷典型日内的总碳排放量结果,构建典型用电量及碳强度计算模型:
式中:Pi为负荷i的典型用电量;Ci为负荷i的典型碳强度。
表2.系统各负荷节点典型日内的碳排放量和碳强度汇总表
由上表信息可知,负荷L4的日内碳排放总量最大,为24.63tCO2;L11的日内碳排放总量最低,仅为2.38tCO2;各负荷的典型碳强度与碳排放量和典型用电量有关,负荷L6的典型碳强度最大,高达0.7892tCO2/MW,其原因在于虽然L6的日内总碳排放量较低,仅为8.88tCO2,但其典型用电量同样较低,为11.20MW,因此,从单位用电量产生的碳排放量角度分析可知,L6的碳强度高居不下。同理,从总量分析的角度出发可知,负荷L9的典型碳强度最低的原因为其单位发电量产生的碳排放量较少,为0.5910tCO2/MW。同时,经计算可知,G3与G4机组对系统的减碳贡献分别为-24.37tCO2和-15.76tCO2,其典型碳强度为-0.6100tCO2/MW和-0.7892tCO2/MW,其含义为新能源机组每发出1MW绿色电能,可为系统减少0.6100tCO2和0.7892tCO2,充分表明该核算方法对可再生能源的推广和应用有着良好的促进作用。
本实施例中,在所述最终碳责任计算模块中,采用如下方式得到得到最终碳责任计算结果:
基于用户的月内用电量数据,采用综合负荷法,构建修正比例计算模型:
式中:Wla为负荷成员a的修正比例,此处假设负荷成员为从成员a-i,以此类推可计算各负荷节点的修正比例;Lat、Lbt、…、Lit分别为各负荷成员典型月内的t时刻最大用电量;εt为用户的典型月最大用电负荷占各时刻用户的典型月最大用电负荷之和的比例;
综合考虑用户的峰荷特性对系统碳排放产生的影响,构建典型用电量修正模型:
式中:为系统负荷侧的修正比例均值,由公式(21)计算的加和平均值得到;Wli为负荷成员i的修正比例,/>为系统负荷侧的典型用电量均值,由公式(19)计算的加和平均值得到;Pci为典型用电量的修正值;
基于所述典型碳强度和修正后的典型用电量,生成反映负荷碳属性的最终碳责任计算结果:
YLi=30CiPci (23)
式中:YLi为负荷i的最终碳责任。
表3.系统各负荷节点的初始碳责任计算结果汇总表
表4.基于用户峰荷特性的碳责任计算结果汇总表
如上表3~4和图4所示,系统负荷侧的总碳责任为3920.63tCO2,其中,负荷L4的碳责任最大,原因在于该节点典型日内的碳排放总量最高;L11的碳责任最低,从总量分析的角度出发,虽然该用户的典型碳强度最高,但是,典型日内的用电量较少,碳排放总量最低,因此无需承担过多的碳责任。之后,充分考虑不同负荷节点在月内的峰荷特性,对典型用电量进行修正,结果表明,节点L2、L3及L13高于修正比例的基准线,意味着这些用户在月内存在过度的高峰用电行为,导致机组的电力供应负担加剧,系统的碳排放量随着增加,不利于低碳电力的发展,因此,需要增加相应比例的碳责任以限制该用户的用电行为。至于其它用户,以用户L4为例,因为其月内的高峰用电水平较低,其碳责任有所减少。依此类推,参考月内的峰荷情况,其余负荷节点的碳责任也做了相应的修正,确保最终结果的公平性和合理性,充分表明了本发明合理性和可行性,为后续系统减碳方案的制定指明了方向。
如图2所示,本发明还提供一种电力系统负荷侧碳责任核算方法,包括以下步骤:
基于电力系统的网络拓扑结构、元件参数和源荷参量,完成系统潮流计算,得到电力系统的有功潮流分布及各支路网损;
根据所述有功潮流分布,建立系统有功潮流顺逆流双向追踪模型,再将所述各支路网损等效为虚拟电源或负荷、各新能源机组等效为负的负荷,结合潮流追踪算法,建立电力系统的网损公平分摊模型,将系统网损公平分摊给源荷两侧,计算系统各净负荷及各机组的净出力;
在完成网损公平分摊的基础上,构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,重新确定其功率分布,并结合机组的碳排放特性,计算得到系统碳流排放的流动状态和进行追踪分析,确定系统各负荷成员的逐时碳排放情况;
基于各负荷成员的逐时碳排放情况,得到各负荷成员典型日的总碳排放量,以此计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;
综合考虑不同用户的峰荷特点,采用综合负荷法,建立典型用电量修正模型,并对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现所述电力系统负荷侧碳责任核算方法的步骤。
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现所述调峰需求响应方法的步骤。
以上所述发明仅表达了本发明实施例的实施方式,并不能因此理解为对发明专利范围的限制,也并非对本发明实施例的结构作任何形式上的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明实施例构思的前提下,还可以做出若干变化和改进,这些都属于本发明实施例的保护范围。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本发明实施例中的方案可以采用各种计算机语言实现,例如,面向对象的程序设计语言Java和直译式脚本语言JavaScript等。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
Claims (9)
1.一种电力系统负荷侧碳责任核算系统,其特征在于,包括系统潮流计算模块、网损公平分摊模型模块、二次潮流计算模块、总碳排放量计算模块和最终碳责任计算模块:
所述系统潮流计算模块用于基于电力系统的网络拓扑结构、元件参数和源荷参量,完成系统潮流计算,得到电力系统的有功潮流分布及各支路网损;
所述网损公平分摊模型模块用于根据所述有功潮流分布,建立系统有功潮流顺逆流双向追踪模型,再将所述各支路网损等效为虚拟电源或负荷、各新能源机组等效为负的负荷,结合潮流追踪算法,建立电力系统的网损公平分摊模型,将系统网损公平分摊给源荷两侧,计算系统各净负荷及各机组的净出力;
所述二次潮流计算模块用于在完成网损公平分摊的基础上,构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,重新确定其功率分布,并结合机组的碳排放特性,计算得到系统碳流排放的流动状态和进行追踪分析,确定系统各负荷成员的逐时碳排放情况;
所述总碳排放量计算模块用于基于各负荷成员的逐时碳排放情况,得到各负荷成员典型日的总碳排放量,以此计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;
所述最终碳责任计算模块用于综合考虑不同用户的峰荷特点,采用综合负荷法,建立典型用电量修正模型,并对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。
2.根据权利要求1要求的一种电力系统负荷侧碳责任核算系统,其特征在于,在所述系统潮流计算模块中,按照以下方式进行系统潮流计算:
解析给定的电力系统的网络拓扑结构,收集系统潮流计算所需的基础参数;
结合所述基础参数对电力系统进行建模仿真,并进行系统的潮流计算,得到电力系统的有功潮流分布。
3.根据权利要求2要求的一种电力系统负荷侧碳责任核算系统,其特征在于,在所述网损公平分摊模型模块中,按照以下方式建立系统有功潮流顺逆流双向追踪模型:
在确定系统潮流分布的基础上,建立系统有功潮流顺逆流双向追踪模型,其表达式如下所示:
上式中[Ad]ij为系统功率的顺流跟踪矩阵;[Au]ij为系统功率的逆流跟踪矩阵;Si为节点i的流过功率;Sij为以i-j为正方向,线路ij流经过的功率;Di为节点i的下游节点集合;Ui为节点i的上游节点集合。
4.根据权利要求3要求的一种电力系统负荷侧碳责任核算系统,其特征在于,在所述网损公平分摊模型模块中,按照以下方式计算系统各净负荷及各机组的净出力:
根据比例共享原则,计算得到各个发电机组和支路对各个负荷节点的贡献,各表达式如下所示:
式中:k为支路中增设的虚拟节点;节点i负荷、支路i-t流出功率可表示为各发电机的功率分量之和;SLk表示新增的虚拟k节点的负荷需求;SGi表示节点i电源出力;Sfi表示支路f-i的流入功率;n为原始网络中的节点数,m为存在负荷的节点个数;SGk是k节点虚拟成电源的网损,SLi是k节点虚拟成负荷的网损,Sit是线路i-t上的功率;式中所有包含字母X的参量为各机组、支路对节点i的贡献量,其中XGi为节点i电源出力对各负荷节点SLk的贡献量,Xfi为支路f-i的流入功率对各负荷节点SLk的贡献量,XLi为节点i负荷表示为各发电机的功率分量,Xit为支路i-t流出功率表示为各发电机的功率分量;
结合上述计算结果,构建网损的公平分摊模型,其中发电侧和用户侧各分摊50%的支路损耗,各负荷及电源承担的网损量的表达式分别如下所示:
再根据各支路网损和各负荷及电源承担的网损量,计算系统各净负荷及各机组的净出力,其表达式分别如下所示:
S'Lk=SLk+△SLk (9)
S'Gk=SGk-△SGk (10)
式中:S'Lk和S'Gk分别表示净负荷和各机组的净出力。
5.根据权利要求4要求的一种电力系统负荷侧碳责任核算系统,其特征在于,在所述二次潮流计算模块中,按照以下方式确定系统各负荷成员的逐时碳排放情况:
在将网络损耗公平分摊给负荷侧和电源侧的基础上,通过直流潮流方程重新确定虚拟无损网络中各支路的有功分布及机组和支路对各负荷节点的贡献,考虑到“源随荷动的情况”,因此再采用潮流逆流追踪的方式,完成网损双向分摊和虚拟网络构建:
P'=Bθ' (11)
式中:θ'∈Rn×1为虚拟网络中各节点电压相角构成的矢量;P'∈Rn×1为各节点注入功率构成的矢量;节点k的注入功率为P'k=P'Gk-P'Lk,P'Gk和P'Lk分别为节点k的电源有功净出力和有功净负荷;B∈Rn×n是以支路电抗形成的节点导纳矩阵;X'Li和X'it分别为节点负荷以及支路功率中各发电机的有功贡献量;P'Li和P'it分别是虚拟无损网络的i节点负荷量和支路i-t的有功通量;
结合机组的碳排放强度,建立系统的碳流追踪与解析模型:
RGk=Re[△SGk]EGk (16)
式中:RLi为负荷节点i的碳流率,其含义为用户每小时用电所产生的等效于发电侧的碳排放量;Rit为支路i-t的碳流率,表示单位时间内随潮流通过的碳流量;RGk为机组k承担的网损碳流率;Re表示取复数ΔSGk的实部;EGk为各机组的碳排放强度向量。
6.根据权利要求5要求的一种电力系统负荷侧碳责任核算系统,其特征在于,在所述总碳排放量计算模块中,按照以下方式计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度:
基于所述系统的碳流追踪与解析模型,计算得到各负荷成员典型日内的总碳排放量,进而得到月内的碳责任量,即:
yLi=30xLi (18)
式中:xLi为负荷i在典型日内的总碳排放量;yLi为负荷i的月内碳责任量;
基于负荷典型日内的总碳排放量结果,构建典型用电量及碳强度计算模型:
式中:Pi为负荷i的典型用电量;Ci为负荷i的典型碳强度。
7.根据权利要求6要求的一种电力系统负荷侧碳责任核算系统,其特征在于,在所述最终碳责任计算模块中,采用如下方式得到得到最终碳责任计算结果:
基于用户的月内用电量数据,采用综合负荷法,构建修正比例计算模型:
式中:Wla为负荷成员a的修正比例;Lat、Lbt、…、Lit分别为各负荷成员典型月内的t时刻最大用电量;εt为用户的典型月最大用电负荷占各时刻用户的典型月最大用电负荷之和的比例;
综合考虑用户的峰荷特性对系统碳排放产生的影响,构建典型用电量修正模型:
式中:为系统负荷侧的修正比例均值;/>为系统负荷侧的典型用电量均值;Pci为典型用电量的修正值;
基于所述典型碳强度和修正后的典型用电量,生成反映负荷碳属性的最终碳责任计算结果:
YLi=30CiPci (23)
式中:YLi为负荷i的最终碳责任。
8.一种电力系统负荷侧碳责任核算方法,其特征在于,包括以下步骤:
基于电力系统的网络拓扑结构、元件参数和源荷参量,完成系统潮流计算,得到电力系统的有功潮流分布及各支路网损;
根据所述有功潮流分布,建立系统有功潮流顺逆流双向追踪模型,再将所述各支路网损等效为虚拟电源或负荷、各新能源机组等效为负的负荷,结合潮流追踪算法,建立电力系统的网损公平分摊模型,将系统网损公平分摊给源荷两侧,计算系统各净负荷及各机组的净出力;
在完成网损公平分摊的基础上,构建虚拟无损网络,采用直流潮流计算方法对虚拟无损网络进行二次潮流计算,重新确定其功率分布,并结合机组的碳排放特性,计算得到系统碳流排放的流动状态和进行追踪分析,确定系统各负荷成员的逐时碳排放情况;
基于各负荷成员的逐时碳排放情况,得到各负荷成员典型日的总碳排放量,以此计算得到各负荷成员的月内碳责任量、典型用电量和典型碳强度;
综合考虑不同用户的峰荷特点,采用综合负荷法,建立典型用电量修正模型,并对各负荷成员的典型用电量进行修正,得到最终碳责任计算结果。
9.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现如权利要求8所述电力系统负荷侧碳责任核算方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311102357.6A CN117151926A (zh) | 2023-08-29 | 2023-08-29 | 一种电力系统负荷侧碳责任核算系统和方法及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311102357.6A CN117151926A (zh) | 2023-08-29 | 2023-08-29 | 一种电力系统负荷侧碳责任核算系统和方法及介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117151926A true CN117151926A (zh) | 2023-12-01 |
Family
ID=88911274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311102357.6A Pending CN117151926A (zh) | 2023-08-29 | 2023-08-29 | 一种电力系统负荷侧碳责任核算系统和方法及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117151926A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117745310A (zh) * | 2024-02-21 | 2024-03-22 | 中国电力科学研究院有限公司 | 一种基于孪生碳网络空间的碳排放量核算方法及系统 |
-
2023
- 2023-08-29 CN CN202311102357.6A patent/CN117151926A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117745310A (zh) * | 2024-02-21 | 2024-03-22 | 中国电力科学研究院有限公司 | 一种基于孪生碳网络空间的碳排放量核算方法及系统 |
CN117745310B (zh) * | 2024-02-21 | 2024-05-10 | 中国电力科学研究院有限公司 | 一种基于孪生碳网络空间的碳排放量核算方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106877397B (zh) | 一种考虑需求侧响应的基于博弈论的主动配电网孤岛恢复方法 | |
CN107634518A (zh) | 一种“源‑网‑荷”相协调的主动配电网经济调度方法 | |
CN113241757B (zh) | 一种考虑柔性负荷与ess-sop的多时间尺度优化调度方法 | |
CN110570327B (zh) | 一种计及源荷互动响应的主动配电网双层规划方法 | |
CN105322535A (zh) | 含统一潮流控制器的电力系统两阶段最优潮流计算方法 | |
CN117151926A (zh) | 一种电力系统负荷侧碳责任核算系统和方法及介质 | |
Suo et al. | New energy wide area complementary planning method for multi-energy power system | |
Deshun et al. | Research on planning and configuration of multi-objective energy storage system solved by improved ant colony algorithm | |
CN109787231A (zh) | 一种综合能源系统分布式能量优化方法及系统 | |
Lei et al. | Energy management considering energy storage and demand response for smart energy hub in Internet of Things | |
Zhou et al. | Energy management for smart energy hub considering gas dispatch factor and demand response | |
CN115293645A (zh) | 一种考虑碳排放累积熵的区域源荷协调调度方法 | |
CN116070754A (zh) | 考虑能量共享的多主体综合能源系统优化运行方法及系统 | |
CN115796447A (zh) | 一种运行虚拟电厂的温室气体减排量计算方法及系统 | |
CN113759705A (zh) | 一种基于合作博弈论的多能源优化配置方法 | |
CN111952977A (zh) | 一种基于潮流的电力市场耦合模型优化方法、装置及介质 | |
Li et al. | Research on Interactive Trading Mode between Thermal Storage Electric Boiler and Wind Power Based on Dynamic Time Sharing | |
CN117689184B (zh) | 考虑荷侧与低碳-经济协同的电力系统规划方法及系统 | |
CN117713176B (zh) | 源网荷储低碳运行方法、装置、电子设备及存储介质 | |
CN113794244B (zh) | 含多微网主动配电系统的定价与最优能量调度方法和系统 | |
CN115983435B (zh) | 一种跨省市电能量市场可行域优化方法及系统 | |
Qing et al. | Determining User's Indirect Carbon Emissions Considering Green Electricity Trading Decomposition | |
WO2021000484A1 (zh) | 跨区域电网的调峰量获取方法和装置 | |
Zhang et al. | Research on deep peaking cost allocation mechanism considering peaking demand subject and thermal power unit | |
Ding et al. | A novel approach to transmission loss rate calculation for electricity transactions |
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 |