👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html
附带附录,灾害前后程序
针对极端灾害导致配电网大面积停电场景,利用多种分布式资源协同恢复重要负荷。在电 网与交通网融合背景下,提出了一种考虑移动储能预布局与动态调度的两阶段配电网韧性提升策 略。在灾前预防阶段,考虑光伏出力不确定性与网络重构建立两阶段鲁棒优化模型,采用Big-M法 和列约束生成算法迭代求解移动储能的配置数量与位置方案。在灾后恢复阶段,考虑灾害对交通 网通行时间的影响,基于移动储能、电动汽车与柴油发电机的动态调度与时序出力特性,建立了多 源协同恢复的混合整数二阶锥规划模型,求解最优的负荷削减与电源出力方案。最后,通过IEEE 33节点配电系统算例验证了所提策略提升配电网韧性的有效性。
部分代码展示:
%% 清除内存空间
clc
clear
close all
warning off
%% 系统参数
mpc = IEEE33; % 配电网参数
SB = mpc.baseMVA; % 基准功率,MW
VB = 12.66; % 基准电压,kV
IB = SB/VB; % 基准电流,kA
T0 = 5; % 故障发生时刻
load curve % 读取负荷,光伏曲线
load index
O_B = mpc.bus(:,1); % 节点集合
NB = length(O_B); % 节点数
O_L = mpc.branch(:,1:2); % 支路集合
NL = length(O_L); % 支路数
O_DG = [2 9 29 6 25 13 18 20 24 33]; % DG节点集合,1-5为DEG,6-10为PV
NDG = length(O_DG); % DG节点数
O_PV = O_DG(6:10); % PV节点集合
NPV = length(O_PV); % PV节点数
O_T = 4:14; % 故障时段集合
NT = 11; % 故障时段总数
O_ME = [12 24]; % 灾前预配置的移动储能集合
alpha_ME0 = zeros(NB,1); % 初始连接的节点
alpha_ME0(O_ME) = 1;
NME = 2; % 灾前预配置的移动储能总数
O_E = [5 15 21]; % 电动汽车节点集合
NE = 3; % 电动汽车总数
D_jk0 = cal_distance(O_L,2); % 相邻节点间电气距离为2km
v0 = 25; % 零交通流量条件下的理想车速为25km/h
c = [0.7470 0.8155 0.8045 0.7297 ... % 灾害场景下交通网络的拥堵程度
0.6100 0.4991 0.4093 0.3559 ...
0.3441 0.3362 0.3415];
v_ME = v0*exp(c); % 实际车速
D_jk = zeros(NB,NB,NT); % t时刻节点j与节点k间的道路距离
T_ME = zeros(NB,NB,NT); % 移动储能在节点j与节点k间的交通通行时间
for t = 1:NT
D_jk(:,:,t) = D_jk0*(1+1/v_ME(t));
T_ME(:,:,t) = D_jk(:,:,t)/v_ME(t);
end
T_ME0 = 0.5; % 安装配置时间
P_M_max = 0.2; % 移动储能有功充放电功率上限值,0.2MW
Q_M_max = 0.17; % 移动储能无功充放电功率上限值,0.17MW
eta_M = 0.98; % 移动储能充放电效率
E_ME_max = 0.6; % 移动储能容量上限,0.6MWh
E_ME_min = E_ME_max*0.15; % 移动储能容量下限,0.09MWh
P_E_max = 0.06; % 电动汽车有功充放电功率上限值,0.06MW
eta_E = 0.98; % 电动汽车充放电效率
E_E_max = 0.12; % 电动汽车容量上限,0.12MWh
E_E_min = E_ME_max*0.15; % 电动汽车容量下限,0.018MWh
essential_user = [3 4 6 10 11 15 ... % 重要用户
17 19 24 26 28 33];
ordinary_user = setdiff(O_B,essential_user);
C_ME = 500; % 预配置单位移动储能设备的成本,500元
w_i = ones(NB,1); % 单位负荷削减成本,元/kW
w_i(essential_user) = 10; % 重要用户负荷削减成本为10元/kW
N_ME_max = 2; % 预布局阶段的移动储能数量上限
tan_max = 0.8; % 功率因数正切的最大值
tan_min = 0.4; % 功率因数正切的最小值
M0 = 99999; % 大M法的中间变量
P_L_max = ones(1,NB); % 负荷有功功率最大值,MW
Q_L_max = zeros(1,NB); % 负荷无功功率最大值,MW
P_L_max(essential_user) = essential_load(T0)*mpc.bus(essential_user,3);
P_L_max(ordinary_user) = ordinary_load(T0)*mpc.bus(ordinary_user,3);
Q_L_max(essential_user) = essential_load(T0)*mpc.bus(essential_user,4);
Q_L_max(ordinary_user) = ordinary_load(T0)*mpc.bus(ordinary_user,4);
R_ij=mpc.branch(:,3); % 线路电阻,pu
X_ij=mpc.branch(:,4); % 线路电抗,pu
P_DG_max = [[0.12 0.12 0.12 ... % DG有功出力上限,MW
0.08 0.08]'*ones(1,NT);pv_curve'/1000];
Q_DG_max = [[0.1 0.1 0.1 0.06... % DG无功出力上限,MW
0.06]'*ones(1,NT);pv_curve'*tan_max/1000];
P_MESS_max = 0.2; % 移动储能有功功率最大值,MW
Q_MESS_max = 0.17; % 移动储能有功功率最大值,Mvar
V_max = [1;1.1*1.1*ones(32,1)]; % 电压平方的上限,pu
V_min = [1;0.9*0.9*ones(32,1)]; % 电压平方的下限,pu
I_max = 1; % 电流平方的下限,pu
alpha_DG = zeros(1,NB); % 各类分布式电源与节点 i 的连接状态
alpha_DG(O_DG) = 1;
delta_i = zeros(NB,NL); % 节点i作为末端节点的支路集合
gamma_i = zeros(NB,NL); % 节点i作为首端节点的支路集合
fault_branch = [1 3 7 28 35]; % 故障支路
for k=1:NL
delta_i(O_L(k,2),k)=1;
gamma_i(O_L(k,1),k)=1;
end
%% 决策变量
alpha_ME1 = binvar(NB,NT); % 0-1变量,表示移动储能1与节点i的连接状态
alpha_ME2 = binvar(NB,NT); % 0-1变量,表示移动储能2与节点i的连接状态
alpha_ij = binvar(NL,NT); % 支路 ij 的开断状态
P_Lsu = sdpvar(NB,NT); % 负荷有功削减功率
Q_Lsu = sdpvar(NB,NT); % 负荷无功削减功率
P_DG = sdpvar(NDG,NT); % DG有功出力
Q_DG = sdpvar(NDG,NT); % DG无功出力
V_sqr = sdpvar(NB,NT); % 节点电压的平方
I_sqr = sdpvar(NL,NT); % 支路电流的平方
F_ij = sdpvar(NL,NT); % 支路虚拟功率流
F_vs = sdpvar(NDG,NT); % 虚拟电源出力
S_vs = binvar(NDG,NT); % 虚拟电源出力标志
P_ij = sdpvar(NL,NT); % 线路有功功率
Q_ij = sdpvar(NL,NT); % 线路无功功率
E_ME = sdpvar(NME,NT); % 移动储能容量
U_Mch = binvar(NME,NT); % 移动储能充电标志
U_Mdch = binvar(NME,NT); % 移动储能放电标志
P_Mch = sdpvar(NME,NT); % 移动储能有功充电功率
P_Mdch = sdpvar(NME,NT); % 移动储能有功放电功率
Q_Mch = sdpvar(NME,NT); % 移动储能无功充电功率
Q_Mdch = sdpvar(NME,NT); % 移动储能无功放电功率
alpha_MCS = binvar(NME,NT); % 移动储能中间变量
E_E = sdpvar(NE,NT); % 电动汽车容量
U_Ech = binvar(NE,NT); % 电动汽车充电标志
U_Edch = binvar(NE,NT); % 电动汽车放电标志
P_Ech = sdpvar(NE,NT); % 电动汽车有功充电功率
P_Edch = sdpvar(NE,NT); % 电动汽车有功放电功率
alpha_ME = binvar(NME,NB,NT);
%% 目标函数
objective = sum(sum((w_i*ones(1,NT)).*P_Lsu*1000*SB));
效果展示:
111号资源-源程序:论文可在知网下载《面向配电网韧性提升的移动储能预布局与动态调度策略》本人博客有解读资源-CSDN文库https://download.csdn.net/download/LIANG674027206/89319950👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆
《《《《《《《《更多资源还请持续关注本专栏》》》》》》》
论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html