电力系统直流潮流计算研究【IEEE9节点】（Matlab代码实现）

 👨‍🎓个人主页：研学社的博客    

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

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

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

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

目录

💥1 概述

📚2 运行结果

2.1 IEEE9节点

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

直流潮流发的特点是用电力系统的交流潮流（有功功率和无功功率）用等值的直流电流来代替。甚至只用直流电路的解析法来分析电力系统的有功潮流，而根不如你不考虑无功分布对有功的影响。这样一来计算速度加快，但计算的准确度有所降低，本方法适用于对潮流计算准确度要求不高的计算场景。
下面先对直流潮流法的原理进行简单介绍：

上图为直流法的等值图，在上图所示的输电线路中，有功潮流为：

 为了快速计算的需要，将上式进行了三项简化：
（1）考虑一般高压电网中线路的电阻远小于电抗，对地电导也可以忽略即 Gii=0 Gij =0
（2）按照标幺值计算时，节点电压与其额定电压相差不大，故有：Ui≈Uj≈1.0；
（3）线路两端的电压相角差（θi-θj）较小，所以有：

这样，上式前两项均为零，只剩第三项

这就相当于线路两端的直流电位分别为θi和θj。线路的直流电阻是Xij。则用矩阵表示为如下式所示。

式中：B0为正常运行时网络的节点电纳矩阵；
θ为网络中各节点的电压相位角的向量；
P为节点注入的有功功率向量；
 

📚2 运行结果

2.1 IEEE9节点

clc
close all
clear 
%% 算例
mpc = case9;

%% 潮流计算
[theta1,P_branch,M,Z,slackbus] = DCpowerflow(mpc);

%% 输出结果
disp('=============================');
disp('支路潮流矩阵')
disp('=============================');
disp('')
disp([num2str(P_branch)]);



disp('=============================');
disp('节点相位矩阵')
disp('=============================');
disp('')
disp([num2str(theta1)]);

clc
close all
clear 
%% 算例
mpc = case9;

%% 潮流计算
[theta1,P_branch,M,Z,slackbus] = DCpowerflow(mpc);

%% 输出结果
disp('=============================');
disp('支路潮流矩阵')
disp('=============================');
disp('')
disp([num2str(P_branch)]);

disp('=============================');
disp('节点相位矩阵')
disp('=============================');
disp('')
disp([num2str(theta1)]);

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]连宏汇,潘欢.直流潮流计算下的孤岛处理模型[J].电工电气,2018(12):8-11+16.

[2]雷婧婷,安婷,杜正春,袁峥.含直流配电网的交直流潮流计算[J].中国电机工程学报,2016,36(04):911-918.DOI:10.13334/j.0258-8013.pcsee.2016.04.003.

🌈4 Matlab代码实现