【复现】考虑泊位优化和多能协同的港口综合能源系统运行优化（Matlab代码实现）

砌墙_2301

于 2024-09-16 10:35:42 发布

阅读量475

点赞数 8

文章标签：能源 matlab 人工智能

本文链接：https://blog.csdn.net/2301_79155639/article/details/142299617

版权

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码、数据

💥1 概述

参考文献：

随着清洁能源的发展，世界各国都在利用清洁能源改善能源结构，提高能源利用率，追求节能减排的效果。区别于传统能源系统，综合能源系统充分结合可再生能源，优化调度内部多能源的使用，为能源需求增加，能源成本上升，环境污染加重等相关问题提供了有效的解决方案。
其中，CCHP系统是同时具备发电、供热和制冷能力的系统，因为其具有能源利用率高的特性，受到了广泛关注。文献[29]建立了一种混合整数非线性规划(Mixed Integer Non-Linear Programming, MINLP)模型，解决CCHP系统的多目标优化问题，在现实约束下，确定满足地区建筑物能源需求的最佳组合和最佳操作策略。求解分析了四种不同场景下的结果，即标准场景，低能耗成本场景，低CO2排放场景和综合场景。文献[30]提出了一个适用于建筑物冷热电联产系统的线性规划优化模型，以优化具有制冷功能的建筑物CCHP系统的规模，相较传统方案，建筑物CCHP系统表现出更高的能源利用效率。文献[31]研究了一种基于混合式制冷系统的CCHP系统的新策略，在优化电力生产的同时实现指定建筑物零多余发电量的目标,结果表明该策略在冬季可以减少高热负荷区域的能源需求和CO2排放量。文献[32]提出了CCHP微电网的两阶段能源管理策略，第一阶段为日前调度阶段，决定了包含机组的输出功率等在内的运行计划，以最大程度地降低运营成本;第二阶段则为实时调整阶段，根据房屋特性带来的功率波动来调整输出功率。文献[33]使用适当的分段线性逼近方法将MINLP问题转换为MILP问题，将其应用于规划CCHP系统的短期调度运行，以最大程度地降低CCHP系统的运行成本。

随着海港全面电气化和能源互联渐趋紧密,未来的海港将成为一种多能源的运输系统。在岸上电力深入推进的背景下，船舶在停泊期间完全由港口供应能源成为响应绿色港口发展目标的思路之一。由于海运行业的船舶类型众多，具有不同电、热、冷负荷需求的船舶进港停泊，需要建立多能协同的港口综合能源系统向在泊船舶供应能源。
如图3-1，论文将PV机组，WT机组，CCHP系统，EES设备及TES设备结合考虑为港口综合能源系统的能源供应侧,该系统在需要时可以从主电网购电。系统内部的机组通过消耗天然气或从主电网购电,发出电功率,热功率或冷功率，向在泊船舶提供相关服务。

📚2 运行结果

其他结果就不一一展示。

部分代码：

B=value(B); %船舶停泊决策变量
Pgrid=value(Pgrid); %主网购电功率
Fpgu=value(Fpgu); %PGU机组燃料消耗量
Ppgu=value(Ppgu); %PGU机组电出力
upgu=value(upgu); %PGU机组出力状态
Qre=value(Qre); %余热回收装置回收量
QHX=value(QHX); %HX机组出力
QEB=value(QEB); %EB输出热功率
PEB=value(PEB); %EB输出电功率
uEB=value(uEB); %EB机组出力状态
Qacex=value(Qacex); %AC输出制冷功率
Qac=value(Qac); %AC消耗热功率
uAC=value(uAC); %AC机组出力状态
QEC=value(QEC); %EC输出冷功率
PEC=value(PEC); %EC消耗电功率
uEC=value(uEC); %EC机组出力状态
p_dch = value(p_dch);%ESS放电功率
p_ch = value(p_ch );%ESS充电功率
u_dch = value(u_dch);%ESS放电状态
u_ch = value(u_ch );%ESS充电状态
E_ess = value(E_ess);%ESS的电量
p_tdch = value(p_tdch);%TES放热功率
p_tch = value(p_tch);%TES充热功率
u_tdch = value(u_tdch);%TES放热状态
u_tch = value(u_tch);%TES充热状态
E_tes = value(E_tes);%TES的储热量
p_wt =value(p_wt);%wt功率
p_pv = value(p_pv);%pv功率