目录
1 主要内容
该程序复现博士文章《互动环境下分布式电源与电动汽车充电站的优化配置方法研究》第五章《单向/双向V2G环境下分布式电源与电动汽车充电站联合配置方法》，本章构建了单向/双向V2G环境下分布式电源与电动汽车充电站联合配置模型，以确定配电系统中分布式电源与电动汽车充电站的最佳安装位置及安装容量。单向V2G环境下，电动汽车的充电时间段被认为是可调节的对象，在充分满足电动汽车充电需求的前提下合理部署相应的充电方案可以有效提高配电系统运行的经济性和安全性。而在双向V2G环境下，电动汽车不仅可以从配电系统获取电能，还能在必要的时候向配电系统释放自身储存的电能，从而为配电系统的运行提供帮助。模型构建的联合配置模型以年化社会总成本最小为目标，采用二阶锥技术求解，最终呈现为一个在多项式时间内可解的混合整数二阶锥规划模型，并以江苏地区一个实际的电气- 地理信息耦合系统进行了测试。
目标函数
电动汽车负荷建模
算例系统图
程序亮点
本程序不仅将31节点系统和实际地理图有机结合，而且考虑了不同情形下（无序充电/单向V2G/双向V2G）配置方案，方便对照学习！
2 部分代码
T = 96;%时段数为1小时
nb = 31;%节点数,根节点为31
nl = 30;%支路数
P1=1e-5.*[500 200 250 180 180 120 240 240 240 240 240 160 160 400 320 300 300 300 300 800 0 80 200 100 60 40 500 200 80 30 0]';
Q1=1e-5.*[300 90 160 90 60 60 80 80 80 120 120 90 100 75 100 125 100 100 100 320 0 40 100 40 25 15 160 120 30 10 0]';
load pc_jm_w;
load pc_jm_wd;
load px_jm_w;
load px_jm_wd;
load pq_jm_w;
load pq_jm_wd;
load pd_jm_w;
load pd_jm_wd;
load pc_sc_w;
load pc_sc_wd;
load px_sc_w;
load px_sc_wd;
load pq_sc_w;
load pq_sc_wd;
load pd_sc_w;
load pd_sc_wd;
load pc_bg_w;
load pc_bg_wd;
load px_bg_w;
load px_bg_wd;
load pq_bg_w;
load pq_bg_wd;
load pd_bg_w;
load pd_bg_wd;
%电动汽车
prl=100;%容量
pev=30*0.25;
cdz=[7 30 60;400 3250 5600;40 325 560;10 10 10];%充电桩数据
cr=106.5;%单位容量增容成本
cl=80;%单位网损费用
plmax=600*1e-4;%线路最大电流
d=0.03;%折现率
num_peak=[1 2 5 1.*ones(1,28)];%测试数据，下为正式数据，运行速度太慢
%num_peak=[30 12 15 11 11 7 14 14 14 14 14 10 10 24 19 18 18 18 18 48 0 5 12 6 4 2 30 12 5 2 0]';%峰值电动汽车停车数量
load arr_jm_w;
load arr_jm_wd;
load arr_sc_w;
load arr_sc_wd;
load arr_bg_w;
load arr_bg_wd;
load stay_jm_w;
load stay_jm_wd;
load stay_sc_w;
load stay_sc_wd;
load stay_bg_w;
load stay_bg_wd;
sty_jd=[2 2 1 3 3 2 1 1 1 1 1 1 1 2 2 2 2 1 1 2 1 3 2 2 3 3 3 1 2 3 1];%节点区域类型
%构建负荷矩阵
for i=1:nb
if sty_jd(i)==1
pload(i,1:T)=P1(i).*pc_jm_w;%春 工作日
pload(i,T+1:2*T)=P1(i).*pc_jm_wd;%春 周末
pload(i,2*T+1:3*T)=P1(i).*px_jm_w;%夏 工作日
pload(i,3*T+1:4*T)=P1(i).*px_jm_wd;%夏 周末
pload(i,4*T+1:5*T)=P1(i).*pq_jm_w;%秋 工作日
pload(i,5*T+1:6*T)=P1(i).*pq_jm_wd;%秋 周末
pload(i,6*T+1:7*T)=P1(i).*pd_jm_w;%冬 工作日
pload(i,7*T+1:8*T)=P1(i).*pd_jm_wd;%冬 周末
elseif sty_jd(i)==2
pload(i,1:T)=P1(i).*pc_sc_w;%春 工作日
pload(i,T+1:2*T)=P1(i).*pc_sc_wd;%春 周末
pload(i,2*T+1:3*T)=P1(i).*px_sc_w;%夏 工作日
pload(i,3*T+1:4*T)=P1(i).*px_sc_wd;%夏 周末
pload(i,4*T+1:5*T)=P1(i).*pq_sc_w;%秋 工作日
pload(i,5*T+1:6*T)=P1(i).*pq_sc_wd;%秋 周末
pload(i,6*T+1:7*T)=P1(i).*pd_sc_w;%冬 工作日
pload(i,7*T+1:8*T)=P1(i).*pd_sc_wd;%冬 周末
else
pload(i,1:T)=P1(i).*pc_bg_w;%春 工作日
pload(i,T+1:2*T)=P1(i).*pc_bg_wd;%春 周末
pload(i,2*T+1:3*T)=P1(i).*px_bg_w;%夏 工作日
pload(i,3*T+1:4*T)=P1(i).*px_bg_wd;%夏 周末
pload(i,4*T+1:5*T)=P1(i).*pq_bg_w;%秋 工作日
pload(i,5*T+1:6*T)=P1(i).*pq_bg_wd;%秋 周末
pload(i,6*T+1:7*T)=P1(i).*pd_bg_w;%冬 工作日
pload(i,7*T+1:8*T)=P1(i).*pd_bg_wd;%冬 周末
end
end
qload=repmat(Q1,1,8*T);
ndjzd=[42,14,0.0492131803564864,1,5;4,27,0.407730811355982,27,3;72,4,0.822394006759020,19,2;95,6,0.827209084946655,19,2;92,17,0.0302701634197687,16,5];
[num_w,~]=size(djz);
[num_wd,~]=size(djzd);
%% 1.设参
mpc = real31;
%pload = mpc.Pload(:,t)*(pload1(t)/sum(mpc.Pload(:,t)))/100;%节点有功负荷
%qload = mpc.Qload(:,1);%节点无功负荷
branch = mpc.branch;
branch(:,3) = branch(:,3)*100/(12.66^2);%求阻抗标幺值
r=real(branch(:,3));
x=imag(branch(:,3));
r=r(1:30);
x=x(1:30);
upstream=zeros(nb,nl);%代表流入节点支路
dnstream=zeros(nb,nl);%代表流出节点支路
for i=1:nl
upstream(i,i)=1;
end
for i=[1:5,7:8,21:25]
dnstream(i,i+1)=1;
end
dnstream(1,[7 10 12 14 15])=1;
dnstream(10,11)=1;
dnstream(12,13)=1;
dnstream(2,[16 18 20])=1;
dnstream(16,17)=1;
dnstream(18,19)=1;
dnstream(22,[27 28])=1;
dnstream(23,[29 30])=1;
dnstream(31,1)=1;
dnstream(31,21)=1;
Vmax=[1.05*1.05*ones(30,1);1.05*1.05*ones(1,1)];
Vmin=[0.95*0.95*ones(30,1);1.05*1.05*ones(1,1)];
Pgmax=[zeros(30,1);ones(1,1)];
Qgmax=[zeros(30,1);ones(1,1)];
%定义变量
V = sdpvar(nb,8*T);%电压的平方
I = sdpvar(nl,8*T);%电流的平方
P = sdpvar(nl,8*T);%线路有功
Q = sdpvar(nl,8*T);%线路无功
Pg = sdpvar(nb,8*T);%发电机有功
Qg = sdpvar(nb,8*T);%发电机无功
Ns = intvar(nb,1);%充电桩数量
Npv = intvar(8,1);%光伏节点安装数量
Ng = intvar(6,1);%燃气轮机安装数量
evc = binvar(num_w,T);%工作日电动汽车充电状态变量
evd = binvar(num_w,T);%工作日电动汽车放电状态变量
evdc = binvar(num_wd,T);%周末电动汽车充电状态变量
evdd = binvar(num_wd,T);%周末电动汽车放电状态变量
3 程序结果
