💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势:🌞🌞🌞博客内容尽量做到思维缜密,逻辑清晰,为了方便读者。
⛳️座右铭:行百里者,半于九十。
📋📋📋本文目录如下:🎁🎁🎁
目录
⛳️赠与读者
👨💻做科研,涉及到一个深在的思想系统,需要科研者逻辑缜密,踏实认真,但是不能只是努力,很多时候借力比努力更重要,然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学,什么是电的时候,不要觉得这些问题搞笑。哲学是科学之母,哲学就是追究终极问题,寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。建议读者按目录次序逐一浏览,免得骤然跌入幽暗的迷宫找不到来时的路,它不足为你揭示全部问题的答案,但若能让人胸中升起一朵朵疑云,也未尝不会酿成晚霞斑斓的别一番景致,万一它居然给你带来了一场精神世界的苦雨,那就借机洗刷一下原来存放在那儿的“躺平”上的尘埃吧。
或许,雨过云收,神驰的天地更清朗.......🔎🔎🔎
💥1 概述
"含氢气氨气综合能源系统优化调度研究"这一主题涉及了能源系统工程、化学工程以及系统优化等多个领域的交叉知识,旨在探讨如何高效、经济、环保地利用氢气和氨气作为能源载体,在综合能源系统中进行优化调度。这类研究对于促进能源结构的转型、提高能源使用效率以及减少碳排放具有重要意义。下面是一个概要性的框架,帮助理解这一领域研究的基本内容和方向:
1. 研究背景与意义
- 背景:随着全球对可持续发展和应对气候变化的关注加深,寻找清洁、高效的能源存储和转换技术成为迫切需求。氢气作为一种清洁能源,具有高能量密度、零排放的特点;氨气作为一种含氢分子,易于液化和储存运输,被视为有潜力的氢能载体。
- 意义:探索含氢气氨气的综合能源系统,不仅能够促进氢能经济的发展,还能为可再生能源的储存和远距离输送提供解决方案,有助于实现能源结构的多元化和低碳化。
2. 系统组成与工作原理
- 系统组成:通常包括氢气和氨气的生产(如水电解制氢、氮气与氢气合成氨)、储存、输运、转换(如燃料电池、燃烧发电)及最终的能源应用环节。
- 工作原理:详细阐述各组成部分如何协同工作,以实现从原料到最终能源输出的全过程,并探讨在不同应用场景下的适应性。
3. 优化调度模型
- 目标函数:建立以经济效益最大化、能源效率最优化或碳排放最小化为目标的数学模型。
- 约束条件:考虑资源限制(如原料供应、储能能力)、运行限制(如设备容量、安全操作范围)、环境政策(如碳排放标准)等因素。
- 算法设计:采用线性规划、非线性规划、混合整数规划、遗传算法、粒子群优化等方法,解决多目标、多约束条件下的调度问题。
4. 关键技术与挑战
- 关键技术:高效制氢/制氨技术、低温液化技术、氢/氨储运技术、燃料电池及氨燃料发动机技术等。
- 挑战:安全性问题(如氢气泄漏、氨气腐蚀)、成本控制(尤其是氢气生产成本)、能效提升、系统集成与管理复杂性等。
5. 案例分析与实证研究
- 通过具体案例,展示含氢气氨气综合能源系统在不同场景(如工业区、偏远地区供电、交通领域)的应用效果,评估其经济性、环境影响和社会接受度。
6. 结论与展望
- 总结当前研究的主要发现,提出未来研究方向,可能包括更高效的能源转化技术、低成本的储运方案、系统运行策略的智能化提升,以及政策与市场机制的设计等。
进行此类研究时,需跨学科合作,整合能源科学、化学工程、计算机科学、经济学等多方面的知识,以推动含氢气氨气综合能源系统的实际应用和技术进步。
📚2 运行结果
部分代码:
%燃气轮机建模约束
C=[C, ECHPmin<=ECHP,ECHP<=ECHPmax,
HCHPmin<=HCHP,HCHP<=HCHPmax,
ECHP(2:24)+HCHP(2:24)-ECHP(1:23)-HCHP(1:23)>= ditaEHCHPmin,
ECHP(2:24)+HCHP(2:24)-ECHP(1:23)-HCHP(1:23)<= ditaEHCHPmax,
(m2qCH4*mCH4CHP)/q2e*nHCHP == HCHP,
(m2qCH4*mCH4CHP)/q2e*nECHP == ECHP,
mCH4CHP == v2mCH4*vch4CHP,
mco2CHP == vch42mco2*vch4CHP,
mco2CHP == v2mco2*vco2CHP,
];
%火电机组
C=[C,
EMmin*YEM<=EM,EM<=EMmax*YEM,
ditaEMmin<=EM(2:24)-EM(1:23),EM(2:24)-EM(1:23)<=ditaEMmax,
mEM==0.320*EM-mNH3*18720/23022,%一度电约等于320g标准煤
% mEM==0.2705489*EM+0.011537-mNH3*18720/23022,%
mco2EM == 2.57*mEM, %1 kW·h电折算标准煤0.32 kg,co2排放0.7976kg
mco2EM == v2mco2*vco2EM,
YEMqi(2:24)-YEMting(2:24)==YEM(2:24)-YEM(1:23),
YEMqi(1)==0,
YEMting(1)==0,
mNH3*18720<=0.2*mEM*23022;
mNH3*18720>=0*mEM*23022;
];
%% 平衡约束
C=[C,ECHP+EM+Ewind ==Eload+ENH3 ];
C=[C,HCHP+HNH3 ==Hload ];
C=[C,VGgrid== vch4CHP ];
F = sdpvar(1,1); %总成本
Fth = sdpvar(1,1);
C=[C, Fth== 15e4*sum(YEMqi+YEMting)+PriM*sum(mEM/1000)];
Fcur = sdpvar(1,1); %弃风成本
C=[C, Fcur == 0.6*sum(Ewindcur) ];
Fbuy = sdpvar(1,1); %购气成本
C=[C, Fbuy == PriG*sum(VGgrid) ];
C=[C, F == Fth+Fcur+Fbuy ];
%% 调用cplex求解器求解
% ops=sdpsettings('verbose', 2, 'solver', 'cplex');
ops = sdpsettings('solver', 'cplex+', 'verbose', 2, 'debug', 1);
ops.cplex.MIPGap =0.0001;
ops.cplex.TimeLimit = 200;
ops.cplex.MIPFocus=3;
ops.cplex.Presolve=2;
ops.cplex.DualReductions = 0;
%进行求解计算
result=optimize(C,F,ops);
if result.problem==0
fprintf('求解成功 !!!\n');
else
error('求解出错!请查找错误来源');
end
%% 绘图
%% 电功率平衡图
ECHP=value(ECHP);EM=value(EM);Ewind=value(Ewind);
Eload=value(Eload);ENH3=value(ENH3);
figure
Plot_E_yuan=zeros(3,24);
for t=1:24
Plot_E_yuan(1,t)=ECHP(t)/1000;
🎉3 参考文献
文章中一些内容引自网络,会注明出处或引用为参考文献,难免有未尽之处,如有不妥,请随时联系删除。
[1]袁文腾,陈亮,王春波,等.基于氨储能技术的电转氨耦合风–光–火综合能源系统双层优化调度[J].中国电机工程学报,2023,43(18):6992-7003.DOI:10.13334/j.0258-8013.pcsee.223152.
🌈4 Matlab代码实现
资料获取,更多粉丝福利,MATLAB|Simulink|Python资源获取