周传炽1. 项目背景与核心价值
电力系统碳排放流计算是当前能源低碳转型中的关键技术之一。这个项目基于IEEE 14节点系统，实现了碳排放流的精确计算与可视化分析。在实际电网运行中，我们经常需要回答这样的问题：某个工厂消耗的1度电，究竟产生了多少碳排放？传统的"平均排放因子"法无法反映电网潮流的时空差异，而碳排放流理论则能提供更精确的溯源分析。
我在参与多个省级电网碳核算项目时发现，现有开源工具在碳排放流计算方面存在两个主要痛点：一是缺乏完整的IEEE标准节点系统案例，二是计算过程的可视化程度不足。这个项目正好填补了这两个空白，特别适合电力系统专业的研究人员和工程师快速上手。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 碳排放流计算原理详解
2.1 理论基础与数学模型
碳排放流计算的核心是建立发电机组碳排放与负荷用电之间的映射关系。其理论基础可以概括为以下三个关键方程：
-
节点碳势方程：
matlab复制
φ_i = (Σ_j∈G_i E_j + Σ_k∈N φ_kP_ki)/P_i其中φ_i表示节点i的碳势（gCO2/kWh），E_j是发电机j的碳排放率，P_ki是支路k→i的功率。
-
支路碳流方程：
matlab复制
C_ij = φ_i × P_ij表示从节点i到j的碳流量（gCO2/h）
-
负荷碳责任分配：
matlab复制
Carbon_load = Σ(φ_i × P_load_i)
2.2 IEEE 14节点系统特点
IEEE 14节点系统包含：
- 5台发电机（节点1,2,3,6,8）
- 11个负荷节点
- 20条输电线路
- 基准功率100MVA
在这个系统中，节点1作为平衡节点，其碳势计算需要特殊处理。我们采用反向潮流追踪法，从末端负荷节点开始逐步回溯到发电节点。
3. Matlab实现关键代码解析
3.1 数据准备与初始化
matlab复制% 发电机碳排放强度 (gCO2/kWh)
gen_emission = [820 650 700 550 450];
% 构建节点导纳矩阵
Ybus = makeYbus(baseMVA, bus, branch);
% 初始化碳势矩阵
phi = zeros(