配电网前推后带法求电力系统潮流(Python&Matlab实现)

目录

 1 概述

2 数学模型

3 节点分层前推回代潮流计算及步骤

3.1 计算方法

3.2 计算步骤

 4 算例及数据

5 Matlab&Python代码实现


 1 概述

配电网通常是单电源全网连接、开环运行,即呈树状。针对配电系统分析,其根本就是进行潮流计算。潮流计算的重要性在配电网规划与设计、配电网重构、配电网故障分析等方面尤为突出。由于配电网的运行特点,其馈线上的分段开关、联络开关会根据不同的情况打开或闭合,使其网络结构发生变化,其潮流的流向也随之改变。因此,进行有效合理的拓扑分析,就大大有利于潮流计算。

提出了一种新的基于节点分层的配网潮流前推回代方法。这种方法采用节点分层来描述配电网,只需得到两个关键的辅助矩阵进行前推回代就能方便快速地进行潮流计算,适合于拓扑结构变化的配网进行快速的潮流计算。

2 数学模型

结合配电网呈辐射状的特点,用节点分层方法来描述,配电网的原始数据采用下面的格式:

支路参数矩阵BranchM

{支路所连节点1 支路所连节点2 支路阻抗参数}

节点参数矩阵NodeM

{节点号 有功负荷 无功负荷}

根据以上的参数矩阵容易求得配网此种结构 运行时的节点-节点关联矩阵 NodeN ,再利用下述

方法求得两个关键矩阵:节点分层矩阵 LayerM 及其对应的上层节点矩阵 NU 。其中 LayerM 矩阵大小根据具体分层来定, NU 为1× N N 为节点数,其中每列的元素是以此列数为节点号对应的上层节点号。

3 节点分层前推回代潮流计算及步骤

3.1 计算方法

因配电网呈辐射状的特点,使其网络上任一节点到源点具有唯一路径,所以可采用结合本文的节点分层前推回代法进行潮流计算。它的基本原理是已知末端负荷和根节点的电压,一般给定配电网络的初始电压,以馈线为计算基本单位,设全网电压为额定电压,根据负荷功率由最后一层向上一层逐层推导,经过一个回代过程可求得各支路电流。再根据给定的始端电压和求得的支路电流,由第一层向下一层逐段计算各段的电压降,求得各点电压,如此反复,直至各个节点的电压率偏差满足容许条件为止。LayerM和NU其具体的算法流程如下。

3.2 计算步骤

第一步:各节点电压初始值为 1,\dot{U}_{i}(0)=1.0,i=1,2,3,.....n,n为配网总节点数。

第二步:利用式(1)和 NU 矩阵,从节点分层矩阵 LayerM 的最后一列向前一列前推计算各支路

电流。

                                   

第三步:利用式(2)和 NU 矩阵,从节点分层矩阵 LayerM 的第一列向下一列回代计算各节点电

压。

                                    

第四步:

                                        

如果不满足条件则跳到第二步进行下一次迭代,直到满足。如果满足则输出最后一次迭代电压的结果。

第五步:根据所求电压求解支路潮流及网损。

详细步骤在附录全部上传啦。

 4 算例及数据

变电站电压=12.66kv,功率基准值=10MVA
 

5 Matlab&Python代码实现

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
基于回代配电网潮流计算程序代码如下: ```matlab function [V, I] = forward_backward_power_flow(A, P, Q, V0, tol, max_iter) % 回代配电网潮流计算程序 % 输入参数: % A: 配电网拓扑结构的节点导纳矩阵 % P: 负载有功功率 % Q: 负载无功功率 % V0: 初始电压 % tol: 收敛精度 % max_iter: 最大迭代次数 % 输出参数: % V: 节点电压 % I: 支路电流 % 获取配电网节点数 n = size(A, 1); % 初始化电压和电流 V = V0; I = zeros(size(A, 2), 1); % 迭代计算 for iter = 1:max_iter % 计算支路电流 for i = 1:n I(A(i, i+1:end) ~= 0) = (conj(P(i)) - 1j*conj(Q(i))) ./ conj(V(i)); end % 回代计算节点电压 for i = n:-1:1 V(i) = (conj(P(i)) + 1j*conj(Q(i))) / (A(i, i)*conj(V(i)) - I(i)); end % 判断收敛 if norm(A*V - I) < tol break; end end % 输出结果 if iter == max_iter warning('未收敛!'); end disp('节点电压:'); disp(V); disp('支路电流:'); disp(I); end ``` 这个程序使用了基于回代的算,通过输入配电网拓扑结构的节点导纳矩阵、负载有功和无功功率、初始电压等参数,计算出每个节点的电压和支路电流等参数。其中,`A` 表示配电网拓扑结构的节点导纳矩阵,`P` 和 `Q` 分别表示负载有功和无功功率,`V0` 表示初始电压,`tol` 表示收敛精度,`max_iter` 表示最大迭代次数。 通过回代的迭代计算,我们可以得到每个节点的电压和支路电流等参数,最终输出结果。需要注意的是,这个代码仅仅是一个简单的示例,实际的配电网潮流计算程序可能会涉及到更复杂的电力系统数学模型和算,需要根据具体情况进行开发和优化。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值