基于前推回代法的连续潮流计算研究【IEEE33节点】（Matlab代码实现）

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

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

目前 的配电网 潮 流 计算 方法 主要 有 改 进 的 牛 顿 － 拉 夫 逊 法 ［1］ 、 改 进 PQ 解 耦 法 ［2］ 、 回 路 阻 抗法［ 3 － 4 ］ 、 前 推 回 代 法 ［ 5 － 7 ］ 等 。 牛 顿 － 拉夫 逊 法 及 改进 PQ 解 耦 法 由 于 其 自身 原 理的 缺 点， 直 接 应用 于配电网潮 流 计算 时 ， 收 敛 性 较 差且 计算 误 差 大 。 回路阻抗 法 主要 应用 于 处 理 环 网 结构 ， 收 敛 性 较 好 ， 但其节点 和 支 路 编 号 处 理 复 杂 且 计算 量 大 。 针对 配电网辐 射 状结构 的 特 点， 前 推 回 代 法是 配电网 潮 流 计算最优 算 法 之一 ， 收 敛 性 好 且 计算 速 度 快 。 传 统的前推 回 代 法 在 计算 前 需要 对 网 络 的 节 点 、 支 路进行预编 号 ， 可 能 导致 计算 错 误 ， 因此文献 ［ 5 － 6 ］ 提出了一 种 自动 编 号 法 ， 但该方法 未考虑 日 益 增加 的分布式电源对 配电网的影响 。 文献 ［7］ 指 出 传 统 前 推 回代 法 的 无 功功 率 取值 可 能与实 际值 出入 较大， 将导致计算 误 差 增 大影响 迭 代 结 果 ， 因此 需要 将 PV

节 点 转 化为 PQ 节 点， 然后对 无 功 和 电 压方 程 进行处理; 然 而 ， 将 所 有 PV 节 点 转 化为 PQ 节 点 又 会导致潮 流 计算 收 敛 性变 差 。

各 类 配电网 潮 流 算 法性 能通 常 从 以 下 几 个 方面进行分析 :

1 ) 算 法 的 收 敛 速 度 。 改 进 的 牛 顿 － 拉夫 逊 法将非 线性 的 方 程 通过 一 定变 换 转 化为 对 相 应线性方

程 的 反 复 求解 ， 且 是 二 阶 算 法 ， 具有 平 方 收 敛 的 特点，能在个 位 次 数 的 迭 代 中 迅 速 求解 。

2 ) 稳定性 。 在配电网中，电力系统网 络 结构 、线路参 数 及 各 种扰 动 因 素 对 计算 结 果 影响的 程 度 即

为 稳定性 。 由 于 牛 顿 拉夫 逊 算 法是 二 阶 算 法 ， 且 受到配电网 Ｒ /X 比值 较 高 的影响 而 无 法 输 出 准 确 的计算结 果 。 而前 推 回 代 法 受 到 的影响 则 很 小 。

3 ) 算 法 的 复 杂 程 度 。 用 简 单 原 理的 算 法 通 常更可 靠 ， 由 于 前 推 回 代 法 不 需要 计算 节 点 导 纳 矩 阵 ，

其 计算 效 率 高 。综上 所述 ， 前 推 回 代 法 的 原 理 更 符 合 实 际 配电网结构 的要 求 ， 其 处 理 环 网 结构 方面 的能力较 弱 ; 但配电网的主要特 点 就是运行 时 呈 辐 射 状 ， 因此并 不会受 到 太 大的影响 。 此外 ， 该方法 计算 误 差 小 、 计算时间 少 ，在系统 异 常 时 仍 能 保 证 输 入有效 结 果 ， 且 收敛性 能 不 会 被 配电网 高 Ｒ /X 数值 影响， 目前 已 经 广泛被 用作 计算 配电网 潮 流 的主要 算 法 。

📚2 运行结果

 

%% 系统支路参数  节点  节点  电阻  电抗
Branch=[1 ,1 ,2 ,0.0922,0.0407;
2 ,2 ,3 ,0.4930,0.2511;
3 ,3 ,4 ,0.3660,0.1864;
4 ,4 ,5 ,0.3811,0.1941;
5 ,5 ,6 ,0.8190,0.7070;
6 ,6 ,7 ,0.1872,0.6188;
7 ,7 ,8 ,0.7144,0.2351;
8 ,8 ,9 ,1.0300,0.7400;
9 ,9 ,10,1.0440,0.7400;
10,10,11,0.1966,0.065;
11,11,12,0.3744,0.1238;
12,12,13,1.4680,1.1550;
13,13,14,0.5416,0.7129;
14,14,15,0.5910,0.5260;
15,15,16,0.7463,0.5450;
16,16,17,1.2890,1.7210;
17,17,18,0.7320,0.5740;
18,2, 19,0.1640,0.1565;
19,19,20,1.5042,1.3554;
20,20,21,0.4095,0.4784;
21,21,22,0.7089,0.9373;
22,3, 23,0.4512,0.3083;
23,23,24,0.8980,0.7091;
24,24,25,0.8960,0.7011;
25,6, 26,0.2030,0.1034;
26,26,27,0.2842,0.1447;
27,27,28,1.0590,0.9337;
28,28,29,0.8042,0.7006;
29,29,30,0.5075,0.2585;
30,30,31,0.9744,0.9630;
31,31,32,0.3105,0.3619;
32,32,33,0.3410,0.5302;]; 

🎉3 参考文献

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

[1]张捷,徐焰,张学飞.一种改进前推回代法的配电网潮流计算[J].四川电力技术,2020,43(03):85-90.DOI:10.16527/j.cnki.cn51-1315/tm.2020.03.017.

[2]张宾.基于改进Matpower连续潮流算法的静态电压稳定性分析[J].西安文理学院学报(自然科学版),2020,23(02):25-29.

[3]张小敏,林群,吴旭光,林振兴,杨惠源,唐爱红.适用于低压配电网潮流计算的前推回代法[J].水电能源科学,2016,34(09):183-186+70.

🌈4 Matlab代码实现