电网优化系列知识整理(一)——直流潮流计算

直流潮流发的特点是用电力系统的交流潮流(有功功率和无功功率)用等值的直流电流来代替。甚至只用直流电路的解析法来分析电力系统的有功潮流,而根不如你不考虑无功分布对有功的影响。这样一来计算速度加快,但计算的准确度有所降低,本方法适用于对潮流计算准确度要求不高的计算场景。
下面先对直流潮流法的原理进行简单介绍:
在这里插入图片描述
上图为直流法的等值图,在上图所示的输电线路中,有功潮流为:
在这里插入图片描述
为了快速计算的需要,将上式进行了三项简化:
(1)考虑一般高压电网中线路的电阻远小于电抗,对地电导也可以忽略即 Gii=0 Gij =0
(2)按照标幺值计算时,节点电压与其额定电压相差不大,故有:Ui≈Uj≈1.0;
(3)线路两端的电压相角差(θi-θj)较小,所以有:
在这里插入图片描述
这样,上式前两项均为零,只剩第三项
在这里插入图片描述
这就相当于线路两端的直流电位分别为θi和θj。线路的直流电阻是Xij。则用矩阵表示为如下式所示。

  • 29
    点赞
  • 132
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个简单的直流电网潮流计算MATLAB程序: ``` % 定义微电网拓扑 branches = [1 2; 1 3; 2 4; 2 5; 3 6; 3 7; 4 8; 5 8; 6 8; 7 8]; % 定义节点功率注入和负载 P_inj = [-1.8; -1.6; -1.2; -2.5; 1.2; 1.6; 0.9; 0.7]; P_load = [1.8; 1.6; 1.2; 2.5; -1.2; -1.6; -0.9; -0.7]; % 定义节点电压 V = [1.05; 1.05; 1.05; 1.05; 1.05; 1.05; 1.05; 1.05]; % 定义节点电阻和导纳 R = [0.1; 0.1; 0.2; 0.2; 0.15; 0.15; 0.1; 0.1]; Y = 1./R; % 构造导纳矩阵 G = zeros(length(V), length(V)); B = zeros(length(V), length(V)); for i = 1:length(branches) G(branches(i, 1), branches(i, 2)) = -Y(i); G(branches(i, 2), branches(i, 1)) = -Y(i); G(branches(i, 1), branches(i, 1)) = G(branches(i, 1), branches(i, 1)) + Y(i); G(branches(i, 2), branches(i, 2)) = G(branches(i, 2), branches(i, 2)) + Y(i); B(branches(i, 1), branches(i, 2)) = -imag(Y(i)); B(branches(i, 2), branches(i, 1)) = -imag(Y(i)); B(branches(i, 1), branches(i, 1)) = B(branches(i, 1), branches(i, 1)) + imag(Y(i)); B(branches(i, 2), branches(i, 2)) = B(branches(i, 2), branches(i, 2)) + imag(Y(i)); end % 定义节点注入功率和负载功率矩阵 P = P_inj - P_load; % 定义节点电流矩阵 I = (G + 1j*B)*V; % 定义节点电压矩阵 V_new = zeros(length(V), 1); % 迭代计算节点电压 while norm(V_new - V) > 1e-6 V = V_new; P_net = P + I.*conj(V); V_new = (G + 1j*B)\P_net; end % 输出结果 disp(V_new); ``` 这个程序使用节点注入功率和负载功率计算潮流,并通过迭代计算节点电压。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值