C类核心论文复现：基于电力系统碳排放流理论的碳排放分摊模型研究程序代码！

本程序参考中文核心论文《基于电力系统碳排放流理论的碳排放分摊模型研究》，文中运用潮流追踪对系统网损引起的碳排放进行分配定量，最后得到系统实时碳排放的产权分配结果，程序较为简单和基础，具有较强的可扩展性和适用性，注释清晰，干货满满，下面对文章和程序作简要介绍。

创新点：

提出一种电力系统碳排放产权界定分配的模型。分摊模型以电力系统碳排放结构分解为切入点，借鉴负荷分析理论将电力系统碳排放以其所对应的电能用途分成了综合用电负荷碳排放、网损碳排放、厂用电碳排放等３部分。确定结构成分并与电力系统碳排放流的计算指标对应后，计算得到系统中各碳排放流成分分布；运用潮流追踪对系统网损引起的碳排放进行 分配定量，最后得到系统实时碳排放的产权分配结果。在对ＩＥＥＥ１４标准节点进行算例分析后，其结果证明了本文所述模型的可行性与正确性。

文中结果：

程序结果：

部分程序

%% 论文复现
K=5;%发电机数
N=14;%系统节点数
M=11;%负荷数
%% 式(2)的计算
PE=zeros(N);%支路潮流分布矩阵
for k=1:length(branch(:,1))
    if branch(k,14)>0
        PE(branch(k,1),branch(k,2))=branch(k,14);
    else
        PE(branch(k,2),branch(k,1))=-branch(k,14);
    end
end
PG=zeros(K,N);%机组注入分布矩阵
for k=1:length(gen(:,1))
    PG(k,gen(k,1))=gen(k,2);
end
xigama=ones(1,K+N);
PZ=[PE;PG];
PN=diag(xigama*PZ);%节点有功通量矩阵
EG=[875;525;0;520;0];%碳排放强度
EN=(PN-PE')^(-1)*PG'*EG;%节点碳势分布向量
%% 式(3)的计算
PL=zeros(M,N);%负荷分布矩阵
bus_copy=bus;
for k=1:M
    for kk=1:N
        if bus_copy(kk,3)>0
            PL(k,kk)=bus_copy(kk,3);
            bus_copy(kk,3)=0;
            break
        end
    end
end
RL=PL*EN;%负荷碳流率向量
RL=RL/1000;%转换单位
%% 式(4)的计算
Pl=zeros(N);%支路有功损耗分布矩阵
for k=1:length(branch(:,1))
    if sum(branch(k,[14 16]))>0
        Pl(branch(k,1),branch(k,2))=sum(branch(k,[14 16]));
    else
        Pl(branch(k,2),branch(k,1))=-sum(branch(k,[14 16]));
    end
end
Rl=diag(EN)*Pl;%支路功率损耗引起的碳排放分布
Rl=Rl/1000;%转换单位
%% 式(5)的计算
Pa=[124.5;40;60;19;20];%厂用电功率？？？？
Ra=diag(EG)*Pa;%厂用电碳流率向量
Ra=Ra/1000;%转换单位
%% 式(6-10)的计算
Ad=eye(N)-PE*PN^(-1);%上游分布矩阵
PG_l=zeros(K,1);%机组对系统有功网损的贡献
RG_l=zeros(K,1);%各个机组应承担的碳排放
J=[1 1;2 2;3 3;4 6;5 8];%各个发电机与其连接的节点
entaN=zeros(K,N);%第j个元素为1，其余分量为0的N维向量
for k=1:K
    entaN(k,J(k,2))=1;
end
for k=1:K
    PG_l(k)=((PG(k,J(k,2))/PN(J(k,2),J(k,2))*entaN(k,:)*Ad^(-1))*Pl)*ones(N,1);
    RG_l(k)=PG_l(k)*EG(k);
end
RG_l=RG_l/1000;%转换单位

欢迎感兴趣的小伙伴关注下方公众号获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！