能源系统的优化
12个视图(30天)
显示旧的评论
你好,
硕士论文我已经开发出下面的脚本,找到最优的技术组合和规模,减少酒店的能源系统的运营成本,同时满足电力和供暖需求(一年半)。
问题是我现在想要添加一个电池随着另一个变量(和优化其能力)允许系统来储存多余的一代在任何时候放电时代低于需求,同时保护目标函数和约束。知道怎么做这个吗?
提前谢谢,
赫克托耳
清晰的
%的负载需求概要
负载(“DemandProfile.mat”)
负载(“WeatherDataFayoum.mat”,“GTilt”,“vwind”,“Tamb”)
%的参数
%光伏
nompv = 0.3;%千瓦
areapv = 1.6;%平方米
co2pv = 6 e - 3;% g二氧化碳/ Wh
effpv = 0.19;%选择效率
% PVT
nompvt = 0.2;%千瓦
areapvt = 1.2;%平方米
co2pvt = 17 e - 3;% g二氧化碳/ Wh
effpvt = 0.15;%选择效率
effpvth = 0.48;%热效率
%圣
nomst = 1.4;%千瓦
areast = 2.14;%平方米
co2st = 17 e - 3;% g二氧化碳/ Wh
effst = 0.64;%热效率
% WT
nomwt = 2.4;%千瓦
areawt = 10.87;%平方米
rhoair = 1.225;%公斤/立方米
co2wt = 4 e - 3;% g二氧化碳/ Wh
effwt = 0.35;
% CHP
co2CHP = 32 e - 3;% g二氧化碳/假名
effCHP = 0.26;%选择效率
effCHPh = 0.62;%热效率
%生物量
massbio = 0.2036;%公斤/千瓦时
nombio = 15;%千瓦
co2bio = 32;% g二氧化碳/千瓦时
effbio = 0.93;%热效率
为i = 1:17520
%的需求数据不同的决议万博 尤文图斯
demandelec (i) =((我)+冷却加热器(我))/ 0.5;% W
demandheat (i) =加热(i);% W
%计算天气条件在不同的决议万博 尤文图斯
辐照度(i) = GTilt(我);
风(i) = vwind(我);
临时(i) = Tamb(我);
%的约束:满足电力和供暖需求
(我:)= -[发光(i) * effpv * 0.5,辐照度(i) * effpvt * 0.5, 0, 0.5 * rhoair * 0.5 * effwt * areawt *(风(i))。^ 3, effCHP * 0.5, 0];
(我+ 17520,)= -[0,辐照度(i) * effpvth * 0.5,辐照度(i) * effst * 0.5, 0, effCHPh * 0.5, 0.5 * nombio * 1000 * effbio];
(我:)= -demandelec(我);
b(我+ 17520,)= -demandheat(我);
结束
Aeq = [];
说真的= [];
%变量
x0 = [0, 0, 0, 0, 0, 0);
磅= [0,0,0,0,0,0);
乌兰巴托=[正无穷,正无穷,正无穷,正无穷,正无穷,正);%待定区域局限等。
%的功能
[x, V] = fmincon (Aeq @fObjFunc x0, A, b,说真的,磅,乌兰巴托);
PV = x (1);
PVT = x (2);
圣= x (3);
WT = x (4);
CHP = x (5);
生物= x (6);
运营成本= V;
%优化
函数f = fObjFunc (x)
%打开输入向量
x1 = x (1);x2 = x (2);x3 = x (3);x4 = x (4);x5 = x (5);x6 = x (6);
%重返参数
%光伏
opexpv_fixed = 3.96;%磅/平方米/年(资本支出的2%)
opexpv_var = 0;%磅/千瓦时
effpv = 0.19;%选择效率
% PVT
opexpvt_fixed = 9.58;%磅/平方米/年(资本支出的1.8%)
opexpvt_var = 0.00125;%磅/千瓦时
effpvt = 0.15;%选择效率
effpvth = 0.48;%热效率
%圣
opexst_fixed = 12.74;%磅/平方米/年(资本支出的1.8%)
opexst_var = 0.00125;%磅/千瓦时(1250 L /兆瓦时,0.1便士/升)
effst = 0.64;%热效率
% WT
opexwt_fixed = 138.84;%磅/单位/年(资本支出的5.2%)
opexwt_var = 0.002;%磅/千瓦时
areawt = 10.87;%平方米
rhoair = 1.225;%公斤/立方米
effwt = 0.35;
% CHP
opexCHP_fixed = 0.073;%磅/ W /年
opexCHP_var = 0.005;%磅/千瓦时
effCHP = 0.26;%选择效率
effCHPh = 0.62;%热效率
%生物量
opexbio_fixed = 243.75;%磅/单位/年(资本支出的5%)
opexbio_var = 0.003;%磅/千瓦时
effbio = 0.93;%热效率
nombio = 15;%千瓦
%指标计算
负载(“WeatherDataFayoum.mat”,“GTilt”,“vwind”,“Tamb”)
%计算每年的一代
为i = 1:17520
%计算天气条件在不同的决议万博 尤文图斯
辐照度(i) = GTilt(我);
风(i) = vwind(我);
%代(千瓦时)
GenPV (i) = x1 *辐照度(我)* effpv * 0.5/1000;
GenPVTe (i) = x2 *辐照度(我)* effpvt * 0.5/1000;
GenPVTh (i) = x2 *辐照度(我)* effpvth * 0.5/1000;
GenST (i) = x3 *辐照度(我)* effst * 0.5/1000;
GenWT (i) = (x4) * 0.5 * rhoair * effwt * areawt *(风(i)。^ 3) * 0.5/1000;
GenCHPe (i) = x5 * effCHP * 0.5/1000;
GenCHPh (i) = x5 * effCHPh * 0.5/1000;
GenBIO (i) = x6 * nombio * effbio * 0.5;
%计算运营成本
OpexPV = opexpv_fixed * x1 + opexpv_var *总和(GenPV);
OpexPVTe (i) = opexpvt_fixed * x2 + opexpvt_var *总和(GenPVTe);
OpexST (i) = opexst_fixed * x3 + opexst_var *总和(GenST);
OpexWT (i) = opexwt_fixed * x4 + opexwt_var *总和(GenWT);
OpexCHPe (i) = opexCHP_fixed * x5 + opexCHP_var *总和(GenCHPe);
OpexBio (i) = opexbio_fixed * x6 + opexbio_var *总和(GenBIO);
结束
OPEXe =总和(OpexPV) +总和(OpexPVTe) + (OpexWT) +金额总和(OpexCHPe);
OPEXt =和(OpexST) +和(OpexBio);
%目标函数
f = OPEXe + OPEXt;
结束
