计及电动汽车灵活性的微网多时间尺度协调调度模型

本文链接：https://blog.csdn.net/weixin_46671758/article/details/138448148

该博客介绍了考虑电动汽车灵活性的微网多时间尺度协调调度模型。通过建立数学模型，考虑了市场电价、柴油发电机、储能设备、需求响应参数以及不同品牌电动汽车的行驶和充电策略。通过Cplex求解器进行优化，以最小化购电费用、柴油机费用、需求响应负荷费用、电动汽车电池损耗费用和可转移负荷补偿费用。最后，展示了各种调度结果，包括负荷转移、分布式电源组合优化、电动汽车调度情况以及电池蓄电量随分时电价的变化情况。

摘要由CSDN通过智能技术生成

代码：

%% 程序初始化
clc
clear
close all

%% 定义变量
%%定义市场电价
%采用冀北电网分时电价工商业电价，18-21h为尖峰电价
xb=[277 277 277 277 277 277 659 978 978 978 659 659 659 659 659 659 659 1041 1041 1041 978 978 277 277];
xs=xb;
%% 定义柴油发电机参数
a=600;%固定开机费用
kcp=300;%分段线性化费用
sconv=600;%启停费用
gtmax=3.31;%出力上限
gtmin=1.3;%最小出力值
ramp=1.5;%爬坡率

%% 定义储能参数
gescmax=2;%充电功率上限
gesdmax=2;%放电功率上限
sessmax=8;%蓄电量上限
sessmin=0;%蓄电量最小值
uesc=0.95;%充电效率
uesd=0.95;%放电效率

%% 需求响应参数
kil=[100,600,800];%中断负荷补偿费用
kshift=[0 0 0 0 0 0 0 0 0 600 600 600 600 600 0 0 0 600 600 600 0 0 0 0];

pwt=0.5*[5.8,6.7,5.8,5.1,6.3,5,6.2,6,4.1,6,7,6.8,6.5,6.9,5,5.6,6,5.8,6.2,4.7,3.3,4.4,5.6,5];

%% 其他输入参数
pmgmax=15;%最大交易量

cil=[0.15,0.1,0.08];%最大中断负荷比例
%负荷值
pload=1.5*[2.2 2.4 2.8 2.5 2.9 3.2 3.5 5.6 6.8 6.2 3.6 3.2 3 2.8 3.3 3.5 3.8 6.3 8.4 7.6 6.3 5.8 3.2 2.6];
%光伏出力
ppv=[0 0 0 0 0 0 3.8 3.9 4.5 5.2 6.5 7.3 7.4 7.2 7.4 6.5 5.5 4.8 0 0 0 0 0 0];

%日产电动汽车行车距离
dvr=[0   0   0   0   0   0   0   11   0   0   0   0   0   0   0   0   11   0   0   0   0   0   0   0];
%比亚迪电动汽车行车距离
dvb=[0   0   0   0   0   0   0   0   22   0   0   0   0   0   0   0   0   0   0   0   22   0   0   0];
%三菱电动汽车行使距离
dvs=[0   0   0   0   0   0   8   0   0   0   0   8   0   8   0   0   0   8   0   0   0   0   0   0];
%宝马min电动汽车行使距离
dvm=[0   0   0   0   0   0   20   0   0   0   0   0   0   0   0   0   0   0   0   20   0   0   0   0];

%% 定义变量sdpvar/binvar
umob=binvar(1,24);%是否购电
umos=binvar(1,24);%是否售电
pmgb=sdpvar(1,24);%市场购电量
pmgs=sdpvar(1,24);%市场售电量
xconv=binvar(1,24);%柴油机工作状态变量，当处于工作状态，定义为1；否则取0
yconv=binvar(1,24);%柴油机启停状态变量
pmt=sdpvar(1,24);%柴油机出力
gesc=sdpvar(1,24);%储能充电功率
gesd=sdpvar(1,24);%储能放电功率
sess=sdpvar(1,24);%蓄电池蓄电量
pil=sdpvar(3,24);%中断负荷

svb=sdpvar(1,24);%比亚迪电动汽车蓄电量
gcvb=sdpvar(1,24);%充电功率
gdvb=sdpvar(1,24);%放电功率
ucvb=binvar(1,24);%充电状态变量
udvb=binvar(1,24);%放电状态变量

svr=sdpvar(