MATLB|分布式能源的选址与定容IEEE30节点实现

news2024/11/25 8:21:05

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

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者
📝目前更新：🌟🌟🌟电力系统相关知识，期刊论文，算法，机器学习和人工智能学习。
🚀支持：🎁🎁🎁如果觉得博主的文章还不错或者您用得到的话，可以关注一下博主，如果三连收藏支持就更好啦！这就是给予我最大的支持！

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

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

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

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

目录

1 分布式能源的选址与定容

2 运行结果

3 Matlab代码实现

1 分布式能源的选址与定容

分布式电源（DG）作为发展智能电网的重要组成部分，具有减轻环境污染、降低网络损耗以及提

高电能质量等优点。按 DG 所使用的一次能源可以主要分为两类：可再生能源 DG 和不可再生能源

DG。可再生能源 DG 主要包括小型风力发电及太阳能光伏发电，其特点是出力受自然条件因素制约，不受调度人员控制；不可再生能源 DG 主要包括燃气轮机发电、燃料电池发电及天然气发电等，其特点是具有很好的调节能力，调度人员可对其发出的有功 P 和无功 Q 进行控制。研究表明，DG 并网对配电网有很大影响，包括电压水平、网损、可靠性等方面，且其影响程度与 DG 的安装位置及容量密切相关。合理的安装位置及容量可以有效改善配电网电压质量、减小有功损耗、提高系统负荷率。配置不合理将适得其反，甚至威胁电网的安全稳定运行。因此，为充分发挥DG 对配电网的积极影响，在规划阶段应对接入配电网的 DG 在最佳安装位置与装机容量两个方面进行重点研究。

2 运行结果

部分代码：

 function ybus=ybuspg_ds33
line_dat=Ldat33;
b_dat=Bdat33;
fb=line_dat(:,2);tb=line_dat(:,3);
Res=line_dat(:,4);
Xe=line_dat(:,5);
Bs=line_dat(:,6);   % 接地导纳, B/2...     
a = line_dat(:,7);              % 抽头设定值..
z = Res + 1i*Xe;                    % 阻抗矩阵...
y = 1./z;                       % 节点导纳
b = 1i*Bs;                        % 使 B 为虚部...
nb = max(max(fb),max(tb));      % 节点个数...
ybus = zeros(nb,nb);               % 初始化节点导纳矩阵..
 br=max(line_dat(:,1));  % br= 支路数
 bus=max(line_dat(:,3)); % bus= 节点个数
 baseMva=max(b_dat(:,7));   % 基准功率
 baseKv=12.62; % 基准电压
 z=zeros(bus,bus);
 zpu=zeros(bus,bus);
 for ii=1:br
      x=line_dat(ii,2);
      y=line_dat(ii,3);
 z(x,y)=complex(line_dat(ii,4),line_dat(ii,5)); % 在 pu 中转换 a 和 y
 zpu(x,y)=z(x,y)*baseMva/(((baseKv)^2)*1000);
 end
 ybus=zeros(bus,bus);
for ii=1:br
    x=line_dat(ii,2);
      y=line_dat(ii,3);
    ybus(x,y)= -(1/zpu(x,y))/a(x);
    ybus(y,x)= -(1/zpu(x,y))/a(x);
    ybus(x,x)= ((ybus(x,x)-ybus(x,y))/a(x))+(1i*(line_dat(ii,6)/2));
    ybus(y,y)= (ybus(y,y)-ybus(y,x)/a(x));
end
ybus;

function ybus=ybuspg_ds33
line_dat=Ldat33;
b_dat=Bdat33;
fb=line_dat(:,2);tb=line_dat(:,3);
Res=line_dat(:,4);
Xe=line_dat(:,5);
Bs=line_dat(:,6); % 接地导纳, B/2...
a = line_dat(:,7); % 抽头设定值..
z = Res + 1i*Xe; % 阻抗矩阵...
y = 1./z; % 节点导纳
b = 1i*Bs; % 使 B 为虚部...
nb = max(max(fb),max(tb)); % 节点个数...
ybus = zeros(nb,nb); % 初始化节点导纳矩阵..
br=max(line_dat(:,1)); % br= 支路数
bus=max(line_dat(:,3)); % bus= 节点个数
baseMva=max(b_dat(:,7)); % 基准功率
baseKv=12.62; % 基准电压
z=zeros(bus,bus);
zpu=zeros(bus,bus);
for ii=1:br
x=line_dat(ii,2);
y=line_dat(ii,3);
z(x,y)=complex(line_dat(ii,4),line_dat(ii,5)); % 在 pu 中转换 a 和 y
zpu(x,y)=z(x,y)*baseMva/(((baseKv)^2)*1000);
end
ybus=zeros(bus,bus);
for ii=1:br
x=line_dat(ii,2);
y=line_dat(ii,3);
ybus(x,y)= -(1/zpu(x,y))/a(x);
ybus(y,x)= -(1/zpu(x,y))/a(x);
ybus(x,x)= ((ybus(x,x)-ybus(x,y))/a(x))+(1i*(line_dat(ii,6)/2));
ybus(y,y)= (ybus(y,y)-ybus(y,x)/a(x));
end
ybus;