基于粒子群算法的电力系统无功优化研究(IEEE14节点)（Matlab代码实现）

💥💥💥💞💞💞欢迎来到本博客❤️❤️❤️💥💥💥
🎉作者研究：***🏅🏅🏅主要研究方向是电力系统和智能算法、机器学习和深度学习。目前熟悉python网页爬虫、机器学习、群智能算法、深度学习的相关内容。希望将计算机和电网有效结合！*⭐️⭐️⭐️

📝目前更新：🌟🌟🌟电力系统相关知识，期刊论文，算法，机器学习和人工智能学习。
🚀****支持：🎁🎁🎁如果觉得博主的文章还不错或者您用得到的话，可以关注一下博主，如果三连收藏支持就更好啦！这就是给予我最大的支持！

                                🎉🎉**欢迎您的到来**🎉🎉

                 ⛅⛅⛅ 📃**个人主页：**科研室**🌈🌈🌈**

                📚📚📚📋**所有代码目录：**电气工程科研社👨‍💻👨‍💻👨‍💻

                       ![](https://img-blog.csdnimg.cn/ec532a82b3c84f0f91482f00dd8e1116.png?x-oss-process=image/watermark,type_d3F5LXplbmhlaQ,shadow_50,text_Q1NETiBA5piv5qKm5ZCn77yM5piv5L2g5ZCn77yB,size_20,color_FFFFFF,t_70,g_se,x_16)

【现在公众号名字改为：荔枝科研社】

📋📋📋本文目录如下：⛳️⛳️⛳️

1 概述

近年来，多种启发式优化算法被应用于电力系 统无功优化，并取得了较好的效果，但容易出现早熟 现象和陷入局部收敛等。粒子 群 优 化 算 法 ( Particle Swarm optimization，PSO) 是一种源于对鸟 群捕食行为的研究而发明的进化计算方法。粒 子群优化算法，对优化问题无可微性与连续性要求， 具有全局收敛性、通用性及鲁棒性强等优点，其 具有全局寻优能力，编程简单。本文采用粒子群改进算法对电力系统的无功优化进行研究，以达到对实际电力传输中的无功优化。

2 无功优化数学模型

无功优化的数学模型由目标函数、功率约束条件、变量约束条件组成，通常采用适当地调整发电机

机端电压、投切无功补偿容量和调节变压器分接头，在保证电压质量的前提下降低网损。从经济性考虑，通常把有功网损最小作为目标函数：

式中: N 为参与有功网损计算的系统支路条数; Gij为节点 i，j 之间的支路电导; Vi、Vj 分别为节点 i 和 j的电压。无功优化的控制变量包括发电机的机端电压、补偿点的补偿量和可调变压器变比等; 状态变量包括发电机无功出力和各节点电压值。各变量的约束条件包括等式约束和不等式约束，其中等式约束即有功、无功的潮流方程［7］，如式( 2) :

           ![\left\{\begin{array}{l} P_{i}=V_{i} \sum_{j=1}^{N_{B}} V_{j}\left(G_{i j} \cos \delta_{i j}+B_{i j} \sin \delta_{i j}\right) \\ Q_{i}=V_{i} \sum_{j}^{N_{B}} V_{j}\left(G_{i j} \sin \delta_{i j}+B_{i j} \cos \delta_{i j}\right) \end{array}\right.](https://latex.codecogs.com/gif.latex?%5Cleft%5C%7B%5Cbegin%7Barray%7D%7Bl%7D%20P_%7Bi%7D%3DV_%7Bi%7D%20%5Csum_%7Bj%3D1%7D%5E%7BN_%7BB%7D%7D%20V_%7Bj%7D%5Cleft%28G_%7Bi%20j%7D%20%5Ccos%20%5Cdelta_%7Bi%20j%7D+B_%7Bi%20j%7D%20%5Csin%20%5Cdelta_%7Bi%20j%7D%5Cright%29%20%5C%5C%20Q_%7Bi%7D%3DV_%7Bi%7D%20%5Csum_%7Bj%7D%5E%7BN_%7BB%7D%7D%20V_%7Bj%7D%5Cleft%28G_%7Bi%20j%7D%20%5Csin%20%5Cdelta_%7Bi%20j%7D+B_%7Bi%20j%7D%20%5Ccos%20%5Cdelta_%7Bi%20j%7D%5Cright%29%20%5Cend%7Barray%7D%5Cright.)

其他约束条件见第4部分。

3 运行结果

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者，博主专门做了一个专栏目录，整个专栏只放了一篇文章，足见我对其重视程度：博主专栏目录。做到极度细致，方便大家进行学习！亲民！！！还有我开了一个专栏给女朋友的，很浪漫的喔，代码学累的时候去瞧一瞧，看一看：女朋友的浪漫邂逅。有问题可以私密博主，博主看到会在第一时间回复。

4 Matlab代码及文章详细阅读

本文仅展现部分代码，全部代码及更多详细文章见：🍞正在为您运送作品详情

本文以IEEE14节点作为算例：，这里仅展现部分代码，全部代码点击上面链接即可获得。

function [TransFormer_Branch, Normal_Branch, PQ_Node, PV_Node, Swing_Node, Node_Num] = RE_IEEE14_data()
%输入IEEE-14节点系统数据
%===============定义支路数据======================

%================变压器支路=========================
%        首节点   末节点     电抗      非标准变比（首节点：末节点）
TransFormer_Branch = [ ...
    6     5       0.252020      0.932
    7     4       0.209120      0.978
    9     4       0.556180      0.969];

%==================输电线路===============================
%         首节点   末节点     支路电阻    支路电抗  接地电纳(－B/2)
Normal_Branch = [ ...
1     2     0.019380       0.059170      -0.026400      
2     3     0.046990       0.197970      -0.021900      
2     4     0.058110       0.176320      -0.018700      
1     5     0.054030       0.223040      -0.024600      
2     5     0.056950       0.173880      -0.017000      
3     4     0.067010       0.171030      -0.017300      
4     5     0.013350       0.042110      -0.006400      
7     8     0.000000       0.176150      0.000000      
7     9     0.000000       0.110010      0.000000      
9     10    0.031810       0.084500      0.000000      
6     11    0.094980       0.198900      0.000000      
6     12    0.122910       0.255810      0.000000      
6     13    0.066150       0.130270      0.000000      
9     14    0.127110       0.270380      0.000000      
10    11    0.082050       0.192070      0.000000      
12    13    0.220920       0.199880      0.000000      
13    14    0.170930       0.348020      0.000000];
% PQ节点定义：（已知值）1节点号  2有功负荷   3无功负荷  | （需求解值）4节点电压幅值    5节点电压相角    6并联元件（电容或电抗）
PQ_Node = [ ...
5      0.076000       0.016000       1.0        0.0         0.0 
7      0.000000       0.000000       1.0        0.0         0.0 
9      0.295000       0.166000       1.0        0.0         -0.19000 
10     0.090000       0.058000       1.0        0.0         0.0 
11     0.035000       0.018000       1.0        0.0         0.0 
12     0.061000       0.016000       1.0        0.0         0.0 
13     0.135000       0.058000       1.0        0.0         0.0 
14     0.149000       0.050000       1.0        0.0         0.0 ];
% PV节点定义：（已知值）1节点号    2注入有功    3节点电压幅值 4节点电压相角 | （需求解值）5注入无功      6有功负荷   7无功负荷  8并联元件（电容或电抗）
PV_Node = [...
1      0.0      1.06            0.0     0.0     0.000000       0.000000     0.0000
2      0.4      1.045000        0.0     0.0     0.217000       0.127000     0.0000
3      0.0      1.010000        0.0     0.0     0.942000       0.190000     0.0000
6      0.0      1.070000        0.0     0.0     0.112000       0.075000     0.0000
8      0.0      1.090000        0.0     0.0     0.000000       0.000000     0.0000];

% 平衡节点定义：（已知值）1节点号   2节点电压幅值   3节点电压相角 | （需求解值）4注入有功      5无功    
 Swing_Node = [ ...
4       1.02        0.0       0.0          0.0     0.478000       -0.039000]; 

Node_Num = 14;

5 结论

从图 2 和表 4 结果可以看出，电力系统的网损由 11． 505 MVA 下降到 7． 91 MVA，各节点电压也比没有用粒子群优化算法前有了相应的提高。这说明，粒子群算法用于无功优化是完全可行性的，其全局寻优能力强，收敛速度快，优化效果好，容易实现，是一种比较理想的无功优化算法。可以得出，基于粒子群算法的无功优化，收敛速度快，能精确地找到全局最优解，系统的网络损耗明显降低，电压值稳定且有了相应的提升，达到了改善电压质量、减少网络损耗和提高电压稳定性的无功优化目的。

👨‍🎓博主课外兴趣：中西方哲学，送予读者：

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。在我这个专栏记录我有空时的一些哲学思考和科研笔记：科研和哲思。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“真理”上的尘埃吧。

 或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

6 参考文献

[1]陆景,庹先国,彭桂力.基于粒子群算法的电力系统无功优化[J].微型机与应用,2017,36(13):19-21+25.DOI:10.19358/j.issn.1674-7720.2017.13.007.

[2]孙华,申方.改进粒子群算法的电力系统无功优化应用[J].黑龙江科技信息,2016(29):80-81.

[3]徐雷. 基于改进粒子群算法的电力系统无功优化[D].西华大学,2016.

[4]赵迪迪,王梦迪.粒子群算法在电力系统无功优化的应用[J].信息技术与信息化,2021(09):246-248.

7 写在最后

部分理论引用网络文献，若有侵权请联系博主删除。