电、冷、热、光伏CCHP的调度附matlab代码.zip

clc; clear; Pel=[0,0,0,0,0,0,0,0,106,106,108,106,106,106,106,106,107,111,56,32,20,14,0,0]; Cooling=[0,0,0,0,0,0,0,28,60,57,66,74,62,57,61,50,50,30,15,11,13,0,0,0]; Heating=[1,1,1,1,1,1,1,88,136,146,145,146,145,150,161,174,171,166,108,95,87,2,1,1]; Pv=[0,0,0,0,0,2.5,9,15.5,21.5,26,28.5,32,31,27,21.5,15,9,3,0,0,0,0,0,0]; %上述代码实现的功能为：电、冷、热、光伏的实时出力 %% %系统参数 %费用的单位都采用RMB Cgas=0.349; %天然气价格 %燃汽轮机参数 Cmt_onoff=3; %启停费用 % Cmt_om=0.015; %燃汽轮机运维费 Pmt_max=60; %燃气轮机运行上限 Pmt_min=0; %燃气轮机运行下限 % Xmc=0.3; %燃气轮机发电效率 % Xmcloss=0.02; %燃气轮机的能量损失 % Pmt_up=30; %爬坡上限 % Pmt_down=-30; %爬坡下限 %燃气锅炉参数 % CBoiler_om=0.0162; %燃气锅炉运维费用 Xboiler=0.9; %燃气锅炉发电效率 Pboiler_max=100; %燃气锅炉运行上限 Pboiler_min=0; %燃气锅炉运行下限 %光伏参数 CPV_om=0.0123; %光伏运维费 Cap_PV=50; %光伏容量 %吸附式制冷机、电制冷机参数 % Cac_om=0.0144; %吸收式制冷机运维费 COPac=1.2; %吸附式制冷机 Pac_max=300; %吸附式制冷机输出上限 Pac_min=0; %吸附式制冷机输出下限 Cec_om=0.0096; %电制冷机运维费 COPec=4; %电制冷效率 Pec_max=300; %电制冷机输出上限 Pec_min=0; %电制冷机输出下限 %热回收设备、热交换器参数 Xhr=0.73; %余热回收效率 Che_om=0.006; %热交换器的运行费用 Xhe=0.9; %热交换器的效率 Phe_max=120; %热交换器输出上限 Phe_min=0; %热交换器输出下限 %电网参数 % 向大电网购电的电价 Toprice=[0.47110,0.47110,0.47110,0.47110,0.47110,0.47110,0.47110,0.47110,0.87590,0.87590,0.87590,0.87590,0.87590,0.87590,1.09470,1.09470,1.09470,0.87590,0.87590,1.09470,1.09470,1.09470,0.87590,0.87590]; % 向大电网买售电的电价 % 没有售电项 Goprice=[0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12,0.12]; Pgrid_max=60; %电网交互功率上限 Pgrid_min=-60; %电网交互功率下限 %蓄电池参数 % CBT_om=0.00636; %蓄电池运维费 Cap_BT=200; %蓄电池容量 BT_ch=0.95; %蓄电池充电效率 BT_dis=0.95; %蓄电池放电效率 BT_Pch=0.2; %蓄电池最大充电率 BT_Pdis=0.4; %蓄电池最大放电率 BT_loss=0.04; %蓄电池自放电率 WBT_max=0.9*Cap_BT; %蓄电池最大储能量 WBT_min=0.2*Cap_BT; %蓄电池最小储能量 % Pch_up=30; %爬坡上限 % Pch_down=-30; %爬坡下限 % Pdisch_up=30; %爬坡上限 % Pdisch_down=-30; %爬坡下限 U0=0; P0=0; WBT0=0.5*Cap_BT; PBT0=0; PTST0=0; k=1; %% %变量定义 P=sdpvar(6,24); %燃气轮机产生的电(四个值，第一个是总的功率表征值，三个是三段的对应值)、 燃气轮机产生的热（一个）、余热锅炉 吸附式制冷机（输入） 电制冷机（输入） %蓄电池（充、放功率） 蓄热槽（充、放功率） 电网（买电为正，卖电为负） 共14个变量 %后面9个是对偶变量（冷热电） W=sdpvar(1,25); %蓄电池电量 蓄热槽热量λc、Πc1、Πc2、λh、Πh1、Πh2、λnl、Πl1、Πl2、Πpv1、Πpv2,光伏的各时段出力，多出一个表示初始能量 %DP=sdpvar(9,25); %对偶变量 U=binvar(2,24); %蓄电池充放状态、蓄电池充放标志位 %% %约束条件 Gmama_c=0; Gmama_h=0; Gmama_nl=0; BM=100; A=zeros(3,23); %燃机 燃气锅炉 充电 放电 电制冷 电网购电 A=[0 0 0 0 4 0; 2.3 0.9 0 0 0 0; 1 0 -1 1 -1 1]; %B=[冷 热/Xhe 电-光伏]; B=zeros(3,24); %B是一个3行24列的矩阵 B(1,:)=Cooling(1:24); %B的第一行是制冷机数据 B(2,:)=Heating(1:24)/Xhe; %B的第二行是热数据 B(3,:)=Pel(1:24)-Pv(1:24); %B的第三行是燃气轮机与光伏出力之差,意味着光伏完全被消纳 lb=[]; ub=[]; F=[]; %B=[冷 热/Xhe 电-光伏]; for t=1:24 %等式约束 F=F+[A*P(:,t)==B(:,t)]; %冷热电三个平衡方程 %上下限约束 lb=[ 0 0 0 0 0 0 ]'; ub=[ 60 100 40*U(1,t) 80*U(2,t) 300 120 ]'; F = F + [ lb <= P(:,t) <= ub ]; % 燃气轮机的相关约束 %蓄电池 if t==1 F = F+ [ W(1,1)==WBT0 ]; end F = F+ [ W(1,t+1)==W(1,t)*(1-BT_loss)+(P(3,t)*BT_ch-P(4,t)/BT_dis) ]; F = F + [ WBT_min<=W(1,t+1)<=WBT_max ]; F = F + [ W(1,25)==WBT0]; % 蓄电池充放电状态互斥 F = F +[ U(1,t)+U(2,t) <= 1 ]; % 蓄电池充放电爬坡约束 if t==1 F = F + [ -30 <= P(3,t)-P0 <= 30]; % else % F = F + [ -30 <= P(11,t)-P(11,t-1) <= 30 ]; end % 蓄电池充放电爬坡约束 if t==1 F = F + [ -30 <= P(4,t)-P0 <=30 ]; % else % F = F + [ -30 <= P(12,t)-P(12,t-1) <= 30 ]; end % 蓄电池充放电状态计数 % if t==1 % F = F +[ U(1,t)-U0 <= U(2,t) ]; % else % F = F +[ U(4,t)-U(4,t-1) <= U(6,t) ]; % end % if t==1 % F = F +[ U(5,t)-U0 <= U(7,t) ]; % else % F = F +[ U(5,t)-U(5,t-1) <= U(7,t) ]; % end % CCHP微电网与大电网交互功率状态互斥 % F = F +[ U(1,t)+U(2,t) <= 1 ]; end %% %费用 f=0; for t=1:24 f=f+(P(1,t)/0.3)*Cgas+(P(2,t)/0.9)*Cgas+Toprice(t)*P(6,t) end ops = sdpsettings('solver', 'cplex'); %启动cplex sol = optimize(F,f,ops); %solvesdp(F,f); m=double(f') aa=double(P); bb=double(W); cc=double(U);