【优化选址】基于遗传算法求解分布式电源的选址定容问题附matlab代码

本文链接：https://blog.csdn.net/matlab_dingdang/article/details/126237716

1 内容介绍

随着我国经济持续高速发展,能源、特别是电能的消耗量越来越大；为满足电能需求,今后一个时期,我国电力行业仍需大规模建设。在化石能源逐渐枯竭、环境压力逐年增大的背景下,在现有配电网上引入分布式发电技术将是电力系统未来发展方向之一。分布式电源既能用于电力调峰,也可为边远地区独立供电,还能在节能环保同时给供电企业带来可观的经济效益。但是,分布式电源的引入会对配电网的节点电压、线路电流、功率分配和供电可靠性造成影响,这种影响与分布式电源安装的位置和容量密切相关。因此,研究配电网中分布式电源的选址定容问题是非常重要的工作。为了描述配电网的物理结构、计算配电网的运行状态,本文建立了配电网物理结构的等效模型,并用牛顿-拉夫逊方法计算配电网的潮流。在此基础上,对普通配电网的供电成本进行了计算。本文建立了一种在已有配电网中解决分布式电源选址定容问题的模型,并应用遗传算法对其进行求解。该模型考虑了四项经济指标：发电总成本、分布式电源安装费、停电罚款、碳排放罚款及环境效益损失。以IEEE 33节点配电网系统为算例,对模型适用性和算法的有效性进行验证。

2 仿真代码

function orgloss(B1,B2,isb,H,pr,Vb,Sbz)

n=33;%n=39 ;

n1=32;%n1=38;

Zb=(Vb^2/Sbz)*100;

for i=1:n1

B1(i,3)=B1(i,3)./Zb;

B1(i,4)=B1(i,4)./Zb;

end

B1;

for i=1:n

B2(i,3)=B2(i,3)/Sbz;

B2(i,4)=B2(i,4)/Sbz;

B2(i,5)=B2(i,5)/Vb;

end

B2;

Y=zeros(n); %zeros就是生成一个全0的矩阵

%创建节点导纳矩阵电导统一为零

for i=1:n1

p=B1(i,1);

q=B1(i,2);

Y(p,q)=Y(p,q)-1/((B1(i,3)+B1(i,4))*B1(i,5));

Y(q,p)=Y(p,q);

Y(p,p)=Y(p,p)+1/(B1(i,3)+B1(i,4))+0.5*B1(i,6);%低压侧不变

Y(q,q)=Y(q,q)+1/((B1(i,3)+B1(i,4))*B1(i,5)^2)+0.5*B1(i,6);%高压侧阻抗乘以变比平方输入时注意低压侧在前

end

%disp('节点导纳矩阵：') ;

G=real(Y);

B=imag(Y);

[B1,B2,Times]=powerflow(B1,B2,n,H,pr,isb,Y);

disp('初始潮流计算结果如下：');

disp('迭代次数：');

Times

disp('节点电压幅值如下（按节点由小到大的顺序）：');

B2(:,5)

%创建Sb，用于存储平衡节点功率

Sb=zeros(1);

for i=1:n

if i==isb

for j=1:n

Sb=Sb+(B2(i,5)+sqrt(-1)*B2(i,6))*conj(Y(i,j))*conj(B2(j,5)+sqrt(-1)*B2(j,6));

end

disp('初始平衡节点功率：')

%求解各条支路的功率

Sij=zeros(n);

for i=1:n

for j=1:n

Sij(i,j)=Sij(i,j)+(B2(i,5)+sqrt(-1)*B2(i,6))^2*conj(B1(1,4))+(B2(i,5)+sqrt(-1)*B2(i,6))*conj(Y(i,j))*(conj((B2(i,5)+sqrt(-1)*B2(i,6)))-conj((B2(j,5)+sqrt(-1)*B2(j,6))));

end

%disp('线路功率分布：')

Sij;

%%%求解系统网损

Ploss=0;

for i=1:n

for j=1:n

Ploss=Ploss+B2(i,5)*B2(j,5)*(G(i,j)*cos(B2(i,6)-B2(j,6))+B(i,j)*sin(B2(i,6)-B2(j,6)));

end

disp('初始网损：')

Ploss*10000000