简介
该程序参考王成山老师的中国电机工程学报文章《冷热电联供微网优化调度通用建模方法》,文章在分析典型冷热电联供(combined cooling,heat and power,CCHP)系统的基础上,提出描述其组成和结构的母线式结构,并围绕该系统结构设计了微网调度优化模型构架。在结构中,选取电气、烟气、蒸汽、热水、空气作为基本母线,与源、负荷、储能和转换装置联接形成微网。使用该结构对各设备进行独立建模,有助于CCHP系统的灵活配置和通用建模。
分析:文章虽然提及了电气、烟气、蒸汽、热水、空气等空气母线,但实际上只用了电、热、空气(冷)母线,也就是3个母线平衡方程,尝试复现此文章的小伙伴需要注意!
文中提到的母线结构
实际用到的母线结构
等式/不等式约束
母线平衡
电母线
热母线
冷母线
需要注意的是,冷母线包含了两个公式(5)和(6),复现时只需考虑公式(5)即可,公式(5)和(6)放在一起,是为了显示热电比的作用!
热电比:有的取1.7,有的取2.3,小编参考了大量文献,结果表明热电比取值在1.7~2.3中均可,因微燃机额定功率和型号改变。
文中结果:
编程步骤:important!
参数输入—变量设置(sdpvar+binvar)—约束条件—目标函数—求解器设置!
程序结果
小编使用了两种复现方法,关注公众号并截图发送,小编会将代码打包给您,感谢您!
部分代码
%系统参数
%费用的单位都采用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; %热交换器输出下限
%电网参数
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; %爬坡下限
部分内容源自网络,侵权联系删除!
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码,小编会不定期更新高质量的学习资料、文章和程序代码,为您的科研加油助力!
更多创新性预测及优化程序请点击公众号首页:《预测优化》|《综合能源》链接!