配电网重构知识及matlab实现
网重构中,很重要的一个约束条件为配网应随时保持开环、辐射的状态: 配电网系统是属于闭环设计但是开环运行的系统,因此,在开关的开闭过程中,随时保持配电网的开环状态时很重要。Mendoza等利用图论,尤其是基本环向量的方法,不仅仅有效地保持了网络的辐射性,还大大地节省了搜寻时间和空间。还有文献对如何利用基本环向量去创造可行解的规则进行了详细阐述,文 中定义了基本节点,外节点,内节点,环向量,基本支路向量以及禁止的组向量 避免不可行解。本文中,应用了一个更加简单但同样有效的基于基本环向量的 方法保持配电网的开环,下面将进行详细阐述:
目录
一、理论分析
二、程序介绍
1.基本环矩阵M的matlab代码
2.智能算法重构代码
完整粒子群算法的重构程序代码链接
一、理论分析
本文中算例采用的是IEEE 33节点系统,因此以IEEE 33节点系统为例说明,配网中保持配网辐射状态的方法,IEEE 33节点系统有33个节点,32个分段开关和5个联络开关,为了使得网络保持辐射运行状态,具体的步骤如下:
1.将系统中的各个支路进行标号,0代表支路上的开关是断开的,1代表支 路上的开关是闭合的; 2.利用图论决定基本的环路,环路的个数与联络开关的个数是一样的,对33 节点系统而言即5个,如下图所示。
3.然后,根据基本的环路列出基本环矩阵M。 M的行数即基本环的个数, 列数为支路数最多的那个环的支路数,其他小一些的环中没有这么多支路的空余 的地方值都填充为0。本文中,为了更简便,每条支路,包括环共有的支路,都只会包含在一个支路环里,这样就省略了节点分类以及还需创建其他支路的麻烦。 编码的规则是从图的最上面的支路编码到最下面的支路,从最左边的支路编码到最右边的支路。公共的支路只会被包含在首次将其编码进去的环内。通过上面的 编码即得到IEEE 33节点的基本环矩阵M:
4.从矩阵M的每一行中选取一个非零元素,在每一次的网络重构中,每个环路中的都有且只有一条支路(非零元素)能被选为断开,每个环中被选中的断开的支路开关状态改为0,其余支路的开关状态都保持为1。
通过上面的步骤,在重构过程中没有非可行解的产生,因此整个重构过程的步骤大大简化。
二、程序介绍
上述理论非常适合应用于智能算法的配电网重构中,首先需要通过配电网参数形成M矩阵,具体代码如下:
1.基本环矩阵M的matlab代码
%产生回路编号矩阵M
b=32; %支路数,开断开关数
n=33; %节点数
LL=5; %联络开关数
Sb=10; %MW
Vb=12.66; %KV
Zb=Vb^2/Sb; %ohm
%第一列存支路号,第二列存首节点号,第三列存尾节点号,第四列存支路自阻抗(ohm),第五列存尾节点给定功率(MW\Mvar)
%其中:33-37支路为5个联络开关
Z=[1 0 1 0.0922+i*0.047 0.1000+i*0.06002 1 2 0.4930+i*0.2511 0.0900+i*0.04003 2 3 0.3660+i*0.1864 0.1200+i*0.08004 3 4 0.3811+i*0.1941 0.0600+i*0.03005 4 5 0.8190+i*0.7070 0.0600+i*0.02006 5 6 0.1872+i*0.6188 0.2000+i*0.10007 6 7 0.7114+i*0.2351 0.2000+i*0.10008 7 8 1.0300+i*0.7400 0.0600+i*0.02009 8 9 1.0440+i*0.7400 0.0600+i*0.020010 9 10 0.1966+i*0.0650 0.0450+i*0.030011 10 11 0.3744+i*0.1238 0.0600+i*0.035012 11 12 1.4680+i*1.1550 0.0600+i*0.035013 12 13 0.5416+i*0.7129 0.1200+i*0.080014 13 14 0.5910+i*0.5260 0.0600+i*0.010015 14 15 0.7463+i*0.5450 0.0600+i*0.020016 15 16 1.2890+i*1.7210 0.0600+i*0.020017 16 17 0.3720+i*0.5740 0.0900+i*0.040018 1 18 0.1640+i*0.1565 0.0900+i*0.040019 18 19 1.5042+i*1.3554 0.0900+i*0.040020 19 20 0.4095+i*0.4784 0.0900+i*0.040021 20 21 0.7089+i*0.9373 0.0900+i*0.040022 2 22 0.4512+i*0.3083 0.0900+i*0.050023 22 23 0.8980+i*0.7091 0.4200+i*0.200024 23 24 0.8960+i*0.7011 0.4200+i*0.200025 5 25 0.2030+i*0.1034 0.0600+i*0.025026 25 26 0.2842+i*0.1447 0.0600+i*0.025027 26 27 1.0590+i*0.9337 0.0600+i*0.020028 27 28 0.8042+i*0.7006 0.1200+i*0.070029 28 29 0.5075+i*0.2585 0.2000+i*0.600030 29 30 0.9744+i*0.9630 0.1500+i*0.070031 30 31 0.3105+i*0.3619 0.2100+i*0.100032 31 32 0.3410+i*0.5362 0.0600+i*0.0400];C=[33 7 20 2.0+i*2.0 0 %联络开关矩阵C34 8 14 2.0+i*2.0 035 11 21 2.0+i*2.0 036 17 32 0.5+i*0.5 037 24 28 0.5+i*0.5 0 ];
% 定义节点到母节点路径矩阵A
A=zeros(b); %对矩阵A进行初始化
for i=1:bA(i,1)=i; %将路径矩阵A的第一列储存为支路号
end for j=1:b %利用Z矩阵中的首末节点列向量进行路径搜索n1=j;k=2;while(j~=0) %一直搜寻到母节点0结束for m=1:bif(Z(m,3)==j) %在末节点定位jA(n1,k)=Z(m,2); j=Z(m,2); %将首节点赋给j,然后再在末节点定位jk=k+1;endendendend%遍历A矩阵形成单回路矩阵Hfor i=1:LLjd=max(intersect(A(C(i,2),:),A(C(i,3),:))); %取回路交叉点jd,保证是单回路h=1;while(A(C(i,2),h)~=jd)M(i,h)=A(C(i,2),h);h=h+1;endh1=1;while(A(C(i,3),h1)~=jd)M(i,h)=A(C(i,3),h1);h=h+1;h1=h1+1;endM(i,h)=C(i,1);end
2.智能算法重构代码
如粒子群算法中,变量可以选取M矩阵每行中的非零元素,通过迭代得到最佳的重构后网络,重构子函数如下所示。
%该为智能算法程序计算子函数,加入不可行解判断程序(在指定节点处加入DG)
b=32;
k=1;
n=33;
LL=5; %联络开关数 z
Sb=10; %MW
Vb=12.66; %KV
Zb=Vb^2/Sb; %ohm
check=1;
checkhl=1;
checkgd=1;
H=[ 2 3 4 5 6 7 33 20 19 18 0 0 0 0 0 0 0 0 0 0 09 10 11 12 13 14 34 0 0 0 0 0 0 0 0 0 0 0 0 0 02 3 4 5 6 7 8 9 10 11 35 21 20 19 18 0 0 0 0 0 06 7 8 9 10 11 12 13 14 15 16 17 36 32 31 30 29 28 27 26 253 4 5 25 26 27 28 37 24 23 22 0 0 0 0 0 0 0 0 0 0];%由函数matrixH生成for i1=1:LLa(1,i1)=H(i1,Swarm1(1,i1));end
%1、判断是否形成环路,F为支路环路关联矩阵(行表示回路,列表示断开开关,若任意两行相同,则表示形成了环路)
F=zeros(5);
for i1=1:LL %回路for i2=1:LL %断开开关if max(a(1,i2)==H(i1,:))F(i1,i2)=1;endend
end
for i1=1:LL-1for i2=i1+1:LLif min(F(:,i1)==F(:,i2))checkhl=0; %disp('出现环路')%出现环路时guan=10000;endend
end
for i1=1:LL %按照断开开关矩阵,剔除Z矩阵中的断开支路j=i1-1;for i2=1:b+LL-jif BranchM(i2,1)==a(1,i1)BranchM(i2,:)=[];breakendend
end
NodeN=zeros(n); %节点-节点关联矩阵A
for i1=1:bNodeN(BranchM(i1,2),BranchM(i1,3))=1;NodeN(BranchM(i1,3),BranchM(i1,2))=1;
endLayerM=[1]; %节点分层矩阵,电源节点号记“1”
NU=zeros(1,n); %上层节点矩阵(有33列的行矩阵)while(checkhl==1)%以下用循环求取矩阵LayerM和NUwhile(checkgd==1)h=1;
while(min(NU(2:33)~=0)==0) %NU矩阵的2-最后都有上层节点了,表示循环结束了
m=max(find(LayerM(:,h))); %m为矩阵LayerM第h列非零元素的个数k=1;
for i1=1:m g=LayerM(i1,h); %LayerM的第i1行第h列元素ss=find(NodeN(g,:)==1);for i2=1:length(ss)if LayerM~=ss(1,i2) %排除相同节点LayerM(k,h+1)=ss(1,i2);NU(1,ss(1,i2))=g;k=k+1; %k表示第h层含有的节点数endend
end
h=h+1; %h表示网络分层的层数
if length(LayerM(1,:))==h-1 %如果网络分层矩阵没有搜索到下层节点,说明形成了断点,后边网络形成了孤岛,与电源节点没有连通回路checkgd=0; %disp('形成孤岛')guan=10000;check=0;break %结束循环end
end
if min(NU(2:33)~=0) %若解可行,已经计算完LayerM,则让其跳出最外层while循环checkgd=0;%disp('可行解')end
endwhile(check==1)
BranchM(:,4)=BranchM(:,4)/Zb; %阻抗标幺化
NodeM(:,2)=NodeM(:,2)/Sb; %功率标幺化
%下面进行分层前推回代法潮流计算
V=ones(n,1); %节点电压
J=zeros(n,1); %支路电流
k=1; %记录迭代次数
V0=zeros(n,1);
t=0;
while(max(abs(V-V0))>1e-3) %判断收敛性(收敛精度设为1e-6)
V0=V; %记录上一次迭代的电压值
%1、回代求支路电流矩阵J
I=conj(NodeM(:,2)./V); %节点注入电流
for i1=h:-1:2ss=find(LayerM(:,i1)~=0);for i2=1:length(ss)if min(LayerM(i2,i1)==NU)J(LayerM(i2,i1))=I(LayerM(i2,i1));elsesumJ=0;gx=find(NU==LayerM(i2,i1)); %求得节点LayerM(i2,i1)的下层节点矩阵(大小不定)for i3=1:length(gx)sumJ=sumJ+J(gx(1,i3));endJ(LayerM(i2,i1))=I(LayerM(i2,i1))+sumJ;endend
end
%2、前推求节点电压矩阵V
Z=zeros(n,1); %支路阻抗
for i1=2:hss=find(LayerM(:,i1)~=0);for i2=1:length(ss)m1=NU(1,LayerM(i2,i1)); %首节点 n1=LayerM(i2,i1); %尾节点for i3=1:b %从BranchM中搜寻得到各支路阻抗矩阵if (BranchM(i3,2)==m1&&BranchM(i3,3)==n1)||(BranchM(i3,2)==n1&&BranchM(i3,3)==m1)Z(n1,1)=BranchM(i3,4);endendV(n1,1)=V(m1,1)-Z(n1,1)*J(n1,1);end
end
k=k+1; %迭代次数k
if k>20guan=10000;break
end
end
%%
Vby=abs(V);
V=V*Vb; %反标幺
Vm=abs(V); %节点电压幅值
%————————————————————————————————
%节点电压上下限约束(排除a=[7 9 2 12 3];这种状况)
%if min(Vm)<Vb*0.9||min(Vm)>Vb%guan=10000;
%end
%————————————————————————————————
Va=angle(V); %节点电压弧度
%下面计算有功网损fploss
Sloss=zeros(n,1);
Z=zeros(n,1); %支路阻抗for i1=2:h %!!!此for循环程序段是为了搜寻网络各支路的首末节点与支路阻抗,应对网络拓扑变化,造成的潮流流向变化ss=find(LayerM(:,i1)~=0);for i2=1:length(ss)m1=NU(1,LayerM(i2,i1)); %首节点 n1=LayerM(i2,i1); %尾节点for i3=1:b %从BranchM中搜寻得到各支路阻抗矩阵if (BranchM(i3,2)==m1&&BranchM(i3,3)==n1)||(BranchM(i3,2)==n1&&BranchM(i3,3)==m1)Z(n1,1)=BranchM(i3,4);endendSloss(n1)=(abs(V(m1)-V(n1)))^2/conj(Z(n1,1)*Zb)*1000; %各支路的功率损耗(单位为kW)end
end
Ploss=real(Sloss); %有功损耗
Qloss=imag(Sloss); %无功损耗
fPloss=sum(Ploss); %系统的总有功损耗
guan=fPloss;
check=0;
end
checkhl=0;end
if k>20guan=10000;
end
完整粒子群算法的重构程序代码链接:
考虑分布式电源的配电网重构优化程序matlab_基于分布式电源的配电网重构优化程序-Matlab文档类资源-CSDN下载配网重构中,很重要的一个约束条件为配网应随时保持开环、辐射的状态:配电网系统是属于闭环设计但是开环基于分布式电源的配电网重构优化程序更多下载资源、学习资料请访问CSDN下载频道.https://download.csdn.net/download/zhangxd212489/73526211
配电网重构知识及matlab实现相关推荐
- 配电网重构知识及matlab实现(2)
上一篇配电网重构文章解决了如何利用基本环方法实现配电网重构,解决了智能算法求解配电网重构的断开线路组合难题,大大提高了重构运算效率,但是对于环网检测和孤岛判断并没有深入分析,这篇文章主要通过配电网分层 ...
- 【配电网重构】基于matlab负荷平衡的配电网重构【含Matlab源码 2180期】
⛄一. 负荷平衡的配电网重构 1 引言 配电网网络重构是指在正常或非正常运行条件下通过改变开关的开关状态来改变网络的拓扑结构.配电网故障恢复供电是指配电网发生故障后,在故障已被定位和隔离的基础上,研究 ...
- 【配电网重构】基于matlab粒子群算法33节点配电网重构【含Matlab源码 1829期】
⛄一.配电网重构简介 1 配电网重构的数学模型 1.1 配电网重构的目标函数 配电网重构的数学模型实质就是多目标多约束的最优解问题, 配网重构的目标一般由以下几类: (1) 配电网最小网损为目标; ( ...
- 基于遗传算法的配电网重构研究(Matlab代码实现)
- 智能软开关 配电网重构matlab 二阶锥 以33节点为研究对象,编制配电网故障重构模型
智能软开关 配电网重构matlab 二阶锥 编程方法:matlab+yalmip(cplex为求解器) 基本内容:以33节点为研究对象,编制配电网故障重构模型,采用图论知识保证配电网的连通性和辐射性, ...
- 基于改进二进制粒子群算法的配电网重构(matlab实现)
目录 一.引言 1.问题背景 2.二进制粒子群算法 2.1简介 2.2 S i g m o i d Sigmoid Sigmoid函数 2.3二进制粒子群算法 2.4算法的改进 二.配电网重构模型 1 ...
- 基于混合整数二阶锥(MISOCP)的配电网重构(附matlab代码)
参考资料:主动配电网网络分析与运行调控 (sciencereading.cn) 配电网重构是指在满足配电网运行基本约束的前提下,通过改变配电网中一个或多个开关的状态对配电网中一个或多个指标进行优化.通 ...
- 【配电网重构】基于粒子群算法实现配电网重构含Matlab源码
1 简介 随着大规模,跨区域的配电网不断发展,对配电网运行的经济性和可靠性要求越来越高,在配电网发生大范围停电事故后,需要对配电网的拓扑结构进行重新组合,从而达到恢复供电的目的,这个重新组合配电网拓扑 ...
- 配电网故障重构(含matlab代码)
之前分享了配电网重构的相关知识和代码,有智能算法(粒子群)和二阶锥算法的,具体内容可翻看往期文章,最近有很多同学在私信咨询粒子群算法的故障重构的程序,本期就这个专题进行分享. 一.故障重构和优化重构的 ...
- 基于粒子群算法的配电网重构研究matlab程序
基于粒子群算法的配电网重构研究matlab程序 参考文献:基于改进灰狼算法的含分布式电源配电网重构研究 (本文未考虑分布式电源) 摘要:使用基本环矩阵编码的智能优化算法在处理配电网重构问题中,通常使用 ...
最新文章
- 14PS中的切图基本操作
- 支持向量机(SVM)的约束和无约束优化、理论和实现
- 工作 10 年后转行程序员,如今 35+ 岁并不焦虑!
- SQL Server数据库中、获得刚插入新记录的自动ID号
- webpack [记录]
- 买房一定要知道的购房误区 买涨不买跌的心态可能得改
- java ant 详解
- 【CCCC】L2-023 图着色问题 (25分),,图的染色判定,遍历
- java chatat delete,StringBuffer deleteCharAt(int index)
- 鼠标屏幕取词技术的原理和实现 (转)
- 千兆路由器怎么设置网速最快_千兆路由器
- ug如何导入excel数据点_ug怎么导入 点数据库
- excel如何去重统计户数_如何用好excel统计函数
- levene ttest
- vi使用手册(zt)
- Excel如何统计数据个数
- AssertionError: ‘.../best.pt‘ acceptable suffix is [‘.pt‘]和ValueError: loaded state dict contains...
- 十一个Python在线学习网站
- 快速查询每个快递单号物流详情,支持国内常用的快递公司
- 首批!工信部下达2021年国家工业专项节能监察任务,涉及270 个数据中心(附名单)...
热门文章
- js中鼠标事件mouseover、mouseenter和mouseleave、mouseout的区别
- GUI 自动测试工具[2021清单]
- 基于Java毕业设计/疫情展示平台源码+系统+mysql+lw文档+部署软件
- 【优秀的截图软件】jietu 2.2.2 for mac
- 基于JAVA的旅游企业财务管理系统源码【包调试】
- 游程编码详解(C语言)
- 2、idea热部署插件JRebel+2020年Jrebel激活码+Springboot web开发+Springboot配置文件详解+thymeleaf模板引擎的使用【Springboot】
- win10禁止易升最有效的办法
- H5网站接入微信支付(H5支付+JSAPI支付)
- SPSS入门教程—对数据进行去重操作的相关方法