【顶级EI完美复现】电力系统碳排放流的计算方法【IEEE 14节点】（Matlab代码实现）-CSDN博客

本文链接：https://blog.csdn.net/weixin_46039719/article/details/147318974

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

⛳️赠与读者

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章

⛳️赠与读者

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能解答你胸中升起的一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。

或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

💥1 概述

文献来源：

摘要：电力系统碳排放流分析的理念为低碳电力带来了新的研究方向。根据潮流计算的结果准确而系统地求解电力系统中碳排放流的分布成为亟待解决的问题。在电力系统碳排放流分析理论的基础上,进一步分析了碳排放流和电力系统潮流计算之间的异同,以及电力系统碳排放流的影响因素、计算体系和计算思路;根据碳排放流求解的需要,定义了一些关键矩阵和向量;结合电力系统潮流计算方法与高等电力网络分析方法,在忽略网损的情况下,建立了电力系统碳排放流的基本计算方法,并通过算例系统验证了该方法的正确性。

关键词：

碳排放流;低碳电力技术;潮流计算;节点碳势;

面对问题：

电力系统碳排放流分析的理念为低碳电力带来了新的研究方向。根据潮流计算的结果准确而系统地求解电力系统中碳排放流的分布成为亟待解决的问题。

创新点（解决办法）：

在电力系统碳排放流分析理论的基础上，进一步分析了碳排放流和电力系统潮流计算之间的异同，以及电力系统碳排放流的影响因素、计算体系和计算思路；根据碳排放流求解的需要，定义了一些关键矩阵和向量；结合电力系统潮流计算方法与高等电力网络分析方法，在忽略网损的情况下，建立了电力系统碳排放流的基本计算方法，并通过算例系统验证了该方法的正确性。

文章目录：

📚2 运行结果

论文结果：

复现结果：

论文结果：

代码复现结果：

部分代码：

%% 求系统潮流
power_flow;
%% 论文复现
K=5;%发电机数
N=14;%系统节点数
M=11;%负荷数
%% 支路潮流分布矩阵PB
PB=Pij*baseMVA;
for k=1:length(branch(:,1))
branch(k,14)=PB(branch(k,1),branch(k,2));
branch(k,15)=PB(branch(k,2),branch(k,1));
end
for k=1:14
for kk=1:14
if PB(k,kk)<0
PB(k,kk)=0;
end
end
end
%% 机组注入分布矩阵
PG=zeros(K,N);
for k=1:length(gen(:,1))
PG(k,gen(k,1))=gen(k,2);
end
%% 辅助变量PZ
xigama=ones(1,K+N);
PZ=[PB;PG];
%% 节点有功通量矩阵
PN=diag(xigama*PZ);
disp('节点有功通量：')
disp((xigama*PZ)')
%% 发电机组碳排放强度向量
EG=[875;525;0;520;0];
%% 节点碳势分布向量
EN=(PN-PB')^(-1)*PG'*EG;
disp('节点碳势分布向量：')
disp(EN)
%% 负荷分布矩阵
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第三列为有功功率
bus_copy(kk,3)=0;
break
end
end
end
%% 支路碳硫率分布矩阵
RB=PB*diag(EN);
RB=RB/1000;
disp('支路碳流率分布矩阵：')
disp(RB)
%% 负荷碳流率向量
RL=PL*EN;
RL=RL/1000;%转换单位
disp('负荷碳流率向量：')
disp(RL)