EI论文复现：基于线性化方法的交直流混合配电系统网架规划程序代码！

适用平台：Matlab+Yalmip+Cplex

程序以系统的投资运行成本最小为目标，考虑了系统配置方案的所有可能性，对系统中每个节点之间的线路连接进行规划，为配电网提供最优的网架配置方案。程序算例丰富、注释清晰、干货满满，可扩展性和创新性很高！下面对文章和程序做简要介绍！

程序创新点：

1）在配电网系统规划时考虑了新能源的接入。

2）考虑了线路的输电容量，使用二阶锥约束对模型进行求解。

主要工作：

随着配电系统的快速发展，含高比例新能源的配电系统得到越来越广泛的应用。为适应配电系统中不同类型的新能源的接入，提出一种新的基于线性化方法的配电系统网架结构规划方法，该方法考虑了网架配置的所有可能性。首先以二进制网络矩阵描述配电系统的网架结构，其次建立网架单层规划模型，将规划变量和运行变量同时进行优化以提升全局寻优能力，最后将单层规划模型进行线性化处理，转化为混合整数线性规划问题以提升求解效率。仿真结果表明，所提方法能够为配电系统提供最佳的网架配置方案，相比于传统的规划方案具有更高的经济性和供电能力，在算法上与经典遗传算法求解的双层规划方法相比具有更好的全局寻优能力和计算效率。

Tip：本程序在复现时没有考虑配电网中交直流输电网架！

1、网架规划的必要性

近年来，全国大部分地区负荷逐年攀升，电网在为地区经济社会发展做出巨大贡献的同时，也暴露出供电能力不足、网架结构薄弱、可靠性有待提高、电网建设难度大等突出问题，为尽快解决现状电网存在的问题，迅速提高供电能力与供电质量，深入推进电网与其他基础设施的协调发展，不断满足社会发展的需求，支撑并促进经济发展规划的顺利实施，全面建设坚强智能电网，特明确规划编制目的：

1)指导高中低压配电网中、长期建设与改造，保证配电网与市政建设改造及大电网的协调发展。

2)促进地区配电网的适度超前发展，优化电网结构，保证向社会和居民的优质、可靠供电

3)促进配电网的优化发展，提高供电可靠率，降低网络损耗，提高电力企业的投入产出效益。

4)提出安全、可靠、经济和灵活的远景目标网架和中期网架，做到近中远期目标网架的平滑、合理过渡。

5)预留变电站用地、电缆通道和架空走廊，将配电网规划纳入市政规划，为电网的持续发展赢得时间和空间，避免无序发展造成的征地和建设费用高、难度大等问题。

6)确定合理的配电网建设改造工程量和投资规模;进行投入产出经济评价，为领导审议、决策提供科学依据;同时提高各地区的经济效益，保证电力企业的可持续发展。

2、网架规划遵循的原则

1)经济：对建设方案进行经济分析，尽可能的减小投资成本，使电力建设项目保持合理的投资回报率，充分发挥投资效益。

2)环保：淘汰可能造成环境影响的设备，采用新技术、新设备，切实做好项目环评工作。

3)优质：合理优化配网结构，完善电力生产、供应和销售的全过程，提高供电效率和供电可靠性。

4)高效：切合实际地规划高压的网架结构，电网中各电压层容量之间，按一定的变电容量比配置，容载比选取应符合相关导则规定。

程序结果：

1）线路投资成本：6680000元；

2）线路损耗：2915.9454元；

3）链接主网成本：0元；

4）光伏建设成本：13753.812元

部分程序：

%% 建线路 、光电  数量约束
C_1=[C_1,sum(x_E1_line)==34];
% 各光伏接入节点均可正常运行
C_1=[C_1,sum(x_pvg)==Num_PVG];
%% sd：首节点，md：末节点
SJD=zeros(Node,5);MJD=zeros(Node,5);
%% 有功无功集合 RIGHT1 有功 RIGHT2 无功
LINE1=zeros(1,Node);      % 以该节点为起始节点线路有几条
LINE2=zeros(1,Node);      % 以该节点为终止节点线路有几条
% 以该节点为起始节点线路有几条LINE1(j)=sum(SJD(j,:)~=0,2);
% 以该节点为终止节点线路有几条LINE2(j)=sum(MJD(j,:)~=0,2);
P_ALL=P_ALL+P_PVG(find(PVG_jr==j),t);Q_ALL=Q_ALL+P_PVG(PVG_jr==j,t)/tan(o2);
P_ALL=P_ALL+P_TR(find(TR_jr==j),t);Q_ALL=Q_ALL+Q_TR(find(TR_jr==j),t);
  % 没有建设线路
C_2=[C_2,P_ALL==0];C_2=[C_2,Q_ALL==0];
  % 首节点功率分布
C_2=[C_2,-sum(P_ij(JD_gl(MJD(j,1:LINE2(j)),j),t))+sum(I_ij(JD_gl(MJD(j,1:LINE2(j)),j),t).*rline(MJD(j,1:LINE2(j)),j))/1000==P_ALL];  % KVA为单位进行换算
C_2=[C_2,-sum(Q_ij(JD_gl(MJD(j,1:LINE2(j)),j),t))+sum(I_ij(JD_gl(MJD(j,1:LINE2(j)),j),t).*xline(MJD(j,1:LINE2(j)),j))/1000==Q_ALL];  % KVA为单位进行换算            
% 末节点功率分布
C_2=[C_2,sum(P_ij(JD_gl(j,SJD(j,1:LINE1(j))),t))==P_ALL];    %  KVA为单位进行换算
C_2=[C_2,sum(Q_ij(JD_gl(j,SJD(j,1:LINE1(j))),t))==Q_ALL];    %  KVA为单位进行换算
% 正常线路功率分布
C_2=[C_2,-sum(P_ij(JD_gl(MJD(j,1:LINE2(j)),j),t))+sum(I_ij(JD_gl(MJD(j,1:LINE2(j)),j),t).*rline(MJD(j,1:LINE2(j)),j))/1000+sum(P_ij(JD_gl(j,SJD(j,1:LINE1(j))),t))==P_ALL];  %  KVA为单位进行换算
C_2=[C_2,-sum(Q_ij(JD_gl(MJD(j,1:LINE2(j)),j),t))+sum(I_ij(JD_gl(MJD(j,1:LINE2(j)),j),t).*xline(MJD(j,1:LINE2(j)),j))/1000+sum(Q_ij(JD_gl(j,SJD(j,1:LINE1(j))),t))==Q_ALL];  %  KVA为单位进行换算           
%% 线路电压平衡约束
C_3=[C_3,U_j(j,t)*1000==U_j(i,t)*1000-2*(P_ij(JD_gl(i,j),t)*rline(i,j)+Q_ij(JD_gl(i,j),t)*xline(i,j))+I_ij(JD_gl(i,j),t)*(rline(i,j)^2+xline(i,j)^2)/1000];  
C_3=[C_3,U_j(i,t)*1000-U_j(j,t)*1000-2*(P_ij(JD_gl(i,j),t)*rline(i,j)+Q_ij(JD_gl(i,j),t)*xline(i,j))+I_ij(JD_gl(i,j),t)*(rline(i,j)^2+xline(i,j)^2)/1000<= M*(1-x_E1_line(find(E1==JD_gl(i,j))))]
C_3=[C_3,U_j(i,t)*1000-U_j(j,t)*1000-2*(P_ij(JD_gl(i,j),t)*rline(i,j)+Q_ij(JD_gl(i,j),t)*xline(i,j))+I_ij(JD_gl(i,j),t)*(rline(i,j)^2+xline(i,j)^2)/1000>=-M*(1-x_E1_line(find(E1==JD_gl(i,j))))]
%% 二阶锥约束
C_4=[C_4,norm([2*P_ij(JD_gl(i,j),t),2*Q_ij(JD_gl(i,j),t),I_ij(JD_gl(i,j),t)-U_j(j,t)],2)<=I_ij(JD_gl(i,j),t)+U_j(j,t)]
%% 节点电压、线路电流约束
 if ismember(JD_gl(i,j),E0)==1
C_5=[C_5,0<=I_ij(JD_gl(i,j),:)<=(SLmax/12.66)^2]
C_5=[C_5,0<=I_ij(JD_gl(i,j),:)<=(SLmax/12.66)^2*x_E1_line(find(E1==JD_gl(i,j)))]
%% 变压器节点功率约束 
C_6=[C_6,P_TRmin<=P_TR<=P_TRmax];C_6=[C_6,P_TRmin<=Q_TR<=P_TRmax];
%% 风光功率约束 考虑是否建设
 % 如果不需要，注释掉
%C_7=[C_7,(1-a_max_PVG)*x_pvg(find(PVG_jr==j))*P_PVG_yc(find(PVG_jr==j),t)<=P_PVG(find(PVG_jr==j),t)<=x_pvg(find(PVG_jr==j))*P_PVG_yc(find(PVG_jr==j),t)];

部分内容源自网络，侵权联系删除！

欢迎感兴趣的小伙伴关注并私信获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！