【随机潮流】基于半不变量的概率潮流计算【IEEE34节点】（Matlab代码实现）

荔枝科研社

已于 2024-04-30 09:16:25 修改

阅读量240

点赞数 4

文章标签： matlab 开发语言

于 2024-03-18 19:33:08 首次发布

本文链接：https://blog.csdn.net/weixin_46039719/article/details/136818805

版权

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据、文章

💥1 概述

参考文献：

随着分布式电源在电力系统中所占比例的不断扩大,研究分布式发电对系统稳态运行的影响势在必行。带分布式发电的潮流计算常常用来评估其并网后对系统的影响,同时它也是分析分布式发电对电网稳定性的影响等其他理论研究工作的基础。文献[3]针对辐射状配电系统,应用前推回推法求解潮流,将系统中分布式电源模拟成PV节点,应用PV节点敏感性矩阵来消除电压幅值的偏差。文献[4]则针对配电网中不同种类的分布式电源建立相应的稳态潮流计算模型:对于用同步电机接入电网的分布式电源﹐如微型燃气轮机等﹐用PV节点来模拟;而对用异步电机接入电网的分布式电源,如风力发电等,则用PQ节点模拟,考虑了节点电压对无功功率的影响。文献[5]则具体研究了风电机的潮流模型,提出了RX模型,此模型充分考虑到了风力发电机的输出功率特性,在常规潮流迭代基础上增加了异步发电机的滑差迭代计算。

本文旨在利用半不变量法进行概率潮流计算，并提供了三种不同的方法：蒙特卡洛模拟法、半不变量法结合Gram-Charlier级数展开，以及半不变量法结合Cornish-Fisher级数展开。通过对这些方法进行比较，可以评估它们的效果。在模型设计过程中，特别考虑了负荷不确定性（服从正态分布）和光伏系统的随机性。以IEEE34节点系统为例，该程序能够计算节点电压和支路潮流的概率密度，并生成累积概率曲线，从而提供更全面的数据分析和决策依据。

📚2 运行结果

部分代码：

%节点电压越限概率
%欧洲标准EN 50160规定电压越限的概率小于某一给定值：P(Ui《0.95 or Ui>1.05)<0.05。
%该m文件用于在传统半不变量法（PLF-CM）和Cornish-Fisher级数展开式求解各节点在标准规定下的越限概率。
function P=ProbCMCF(fwd,fws,n)
Vmin=0.95;
Vmax=1.05;
P1=zeros(n,1);
P2=zeros(n,1);
P=zeros(n,1);
for i=1:n
if sum(fwd(i,:)<=Vmin)>0
[m1,n1]=find(fwd(i,:)<Vmin);
ndown=max(n1);
P1(i,1)=fws(ndown);
else
disp('电压未越下限');
end
if sum(fwd(i,:)>=1.05)>0
[m2,n2]=find(fwd(i,:)>Vmax);
nup=min(n2);
P2(i,1)=1-fws(nup);
else
disp('电压未越上限');
end
P=P1+P2;
end

%负荷的八阶半不变量形成-------------------------------------------------------
%%pdfload(i,1)为负荷序号
%%pdfload(i,2)为负荷的节点号
%%pdfload(i,3)为负荷有功均值
%%pdfload(i,4)为负荷无功均值
%%pdfload(i,5)为负荷有功标准差
%%pdfload(i,6)为负荷无功标准差 %%标准差给定可以参照“3Sita原则”
%负荷的八阶半不变量---------------------------------------------------------
pdfload=textread('IEEE34load_30%.txt');%%负荷负荷正态分布
nload=length(pdfload(:,1));
PlPx=zeros(Nodes,8);
PlQx=zeros(Nodes,8);
PlPx(pdfload(:,2),:)=NcalPLCum(-pdfload(:,3),-pdfload(:,5));
PlQx(pdfload(:,2),:)=NcalPLCum(-pdfload(:,4),-pdfload(:,6));

%-------光伏随机特性建模-----------------------------
%选择上海31°8’N、121°35’E作为光照强度分布的考量位置，在HOMERE软件上获取光强分布的期望值和方差。
%miu=0.150314263;
%sita=0.049758487;