无锡哪里有做网站,中国建盏logo,wordpress可视化建站,河北智慧团建官网登录入口目录 1 主要内容
模型示意图
目标函数
程序亮点
2 部分程序
3 程序结果
4 下载链接 1 主要内容
本程序参考文献《考虑可再生能源消纳的建筑综合能源系统日前经济调度模型》模型#xff0c;建立了电热综合能源系统优化调度模型#xff0c;包括燃气轮机、燃气锅炉、余热…目录 1 主要内容
模型示意图
目标函数
程序亮点
2 部分程序
3 程序结果
4 下载链接 1 主要内容
本程序参考文献《考虑可再生能源消纳的建筑综合能源系统日前经济调度模型》模型建立了电热综合能源系统优化调度模型包括燃气轮机、燃气锅炉、余热锅炉、蓄电池、储热、风电和光伏等设备模型注意不含文献中的地源热泵和空调设备可以作为比较经典的电热综合能源优化调度学习资料。 模型示意图 注意本程序不含地源热泵模型。 目标函数
目标函数:经济性燃料成本设备运维成本环境保护 目标函数代码 for t1:24 zz0.025*P_GT(t)0.0097*Q_GB(t)0.025*Q_RB(t)...%设备运维费用 0.03*(abs(P_ESc(t))abs(P_ESd(t)))0.02*(abs(P_HSc(t))abs(P_HSd(t)))...%储能设备 3.15*(F_GT(t)F_GB(t))...%买天然气的价格 0.088*(P_GT(t)*0.242Q_GB(t)*0.085)0.624*(P_GT(t)*0.0036Q_GB(t)*0.764)1.95*(P_GT(t)*0.2Q_GB(t)*0.54);%污染物处理 end 程序亮点
程序对于蓄电池和储热设备的约束形式较常规方式有所区别见部分程序59-66行虽然这种约束方式略显笨重但提供了另外一种程序编写方式顺便了解implies命令应用方式。程序对于蓄电池和储热采用自放电/热系数方式因此可以合理解释结果图2中储热满足24小时充放电功率一致但是储热量却一直下降的原因。
2 部分程序
%% 程序初始化
clc;
clear;
close all;
%% 定义变量
%电设备
P_PVsdpvar(1,24,full);%光伏发电
P_WTsdpvar(1,24,full);%风力发电
P_GTsdpvar(1,24,full);%燃气轮机电功率出力
Ebsdpvar(1,24,full);% 电储能的容量
P_ESsdpvar(1,24,full);%电储能输出有正有负
P_EScsdpvar(1,24,full);% 电储能充电功率
P_ESdsdpvar(1,24,full);% 电储能放电功率
%热设备
Q_GTsdpvar(1,24,full);%燃气轮机输出的热功率
Q_RBsdpvar(1,24,full);%余热锅炉输出的热功率
Q_GBsdpvar(1,24,full);%燃气锅炉输出的热功率
Ehsdpvar(1,24,full);%热储能的容量
P_HSsdpvar(1,24,full);%热储能输出有正有负
P_HScsdpvar(1,24,full);%热储能储热功率
P_HSdsdpvar(1,24,full);%热储能放热功率
%辅助变量
Ubcbinvar(1,24,full);%电储能充电状态1表示充电
Ubdbinvar(1,24,full);%电储能放电状态1表示放电
Uhcbinvar(1,24,full);%热储能储热状态1表示储热
Uhdbinvar(1,24,full);%热储能放热状态1表示放热
F_GTsdpvar(1,24,full);%燃气轮机消耗的天然气量
F_GBsdpvar(1,24,full);%燃气锅炉消耗的天然气量
%% 常量数据,包括电热冷负荷和风光出力
Pwt[66.9 68.2 71.9 72 78.8 94.8 114.3 145.1 155.5 142.1 115.9 127.1 141.8 145.6 145.3 150 206.9 225.5 236.1 210.8 198.6 177.9 147.2 58.7];
Ppv[0,0,0, 0, 0.06, 6.54, 20.19, 39.61, 49.64, 88.62, 101.59, 66.78, 110.46, 67.41, 31.53, 50.76,20.6,22.08,2.07,0,0,0,0,0];
load_e[88.24 83.01 80.15 79.01 76.07 78.39 89.95 128.85 155.45 176.35 193.71 182.57 179.64 166.31 164.61 164.61 174.48 203.93 218.99 238.11 216.14 173.87 131.07 94.04];
load_h[122 125 128 125 128 119 113 109 100,90 88,79 66,55 63,71 98 109 139 150 142 139 136 134];
%% 目标函数:经济性燃料成本设备运维成本环境保护
z0;
for t1:24zz0.025*P_GT(t)0.0097*Q_GB(t)0.025*Q_RB(t)...%设备投资运行费用0.043*P_WT(t)0.0096*P_PV(t)0.03*(abs(P_ESc(t))abs(P_ESd(t)))0.02*(abs(P_HSc(t))abs(P_HSd(t)))...%储能设备3.15*(F_GT(t)F_GB(t))...%买天然气的价格0.088*(P_GT(t)*0.242Q_GB(t)*0.085)0.624*(P_GT(t)*0.0036Q_GB(t)*0.764)1.95*(P_GT(t)*0.2Q_GB(t)*0.54);%污染物处理
endFz%% 约束条件
C[];
%机组约束,决策变量
C[C,0.8*PpvP_PVPpv];%光伏发电上下限约束
C[C,0.8*PpvP_WTPwt];%风力发电上下限约束
for t1:24C[C,0P_GT(t)550];%燃气轮机上下限约束√C[C,0Q_RB(t)750];%余热锅炉上下限约束√C[C,0Q_GB(t)700];%燃气锅炉上下限约束√C[C,0Q_GT(t)200];%燃气轮机输出的热功率
%设备模型
C[C,Ubc(t)Ubd(t)1];C[C,Uhc(t)Uhd(t)1];%确定充放电状态不可能同时发生C[C,-80P_ES(1,t)80,0P_ESc(1,t)80,-80P_ESd(1,t)0]; C[C,-80P_HS(1,t)80,0P_HSc(1,t)80,-80P_HSd(1,t)0]; C[C,implies(Ubc(1,t),[P_ES(1,t)0,P_ESc(1,t)P_ES(1,t),P_ESd(1,t)0])];%储电负荷C[C,implies(Ubd(1,t),[P_ES(1,t)0,P_ESd(1,t)P_ES(1,t),P_ESc(1,t)0])];C[C,implies(Uhc(1,t),[P_HS(1,t)0,P_HSc(1,t)P_HS(1,t),P_HSd(1,t)0])];%储热负荷C[C,implies(Uhd(1,t),[P_HS(1,t)0,P_HSd(1,t)P_HS(1,t),P_HSc(1,t)0])];3 程序结果 4 下载链接
