基于储能电站服务的冷热电多微网系统双层优化配置附Matlab代码

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。​

⛄ 内容介绍

随着储能技术的进步和共享经济的发展，共享储能电站服务模式将成为未来用户侧储能应用的新形态。提出基于储能电站服务的冷热电联供型多微网系统双层优化配置方法。首先，提出储能电站服务这种新型的共享储能模式，分析共享储能电站的运行方式和盈利机制。其次，将储能电站服务应用到冷热电联供型多微网系统中，建立考虑两个不同时间尺度问题的双层规划模型，上层模型负责求解长时间尺度的

储能电站配置问题，下层模型负责求解短时间尺度的多微网系 统 优 化 运 行 问 题 。再 次 ， 根 据 下 层 优 化 模型的Karush-Kuhn-Tucher(KKT)条件将下层模型转换为上层模型的约束条件，采用 Big-M 法对非线性问题线性化。最后，通过 3 个场景的算例分析验证所提双层规划模型的合理性和有效性，并证明所提出的共享储能服务能够有效降低用户成本，节约储能资源，实现用户与储能电站运营商的互利共赢。

⛄ 部分代码

%========================================================================

%电网技术论文复现

%基于储能电站服务的冷热电多微网系统双层配置

%双层规划,KKT条件

%========================================================================

figure

bar(-Pessswi(1,:));

hold on

bar(Pessbwi(1,:));

hold on 

legend('MG1向储能电站售电功率','MG1向储能电站购电功率');

xlabel('时间/h');

ylabel('功率/kW');

title('微网向储能电站购售电功率结果图');

box off

figure

bar(-Pessswi(3,:));

hold on

bar(Pessbwi(3,:));

hold on 

legend('MG3向储能电站售电功率','MG3向储能电站购电功率');

xlabel('时间/h');

ylabel('功率/kW');

title('微网向储能电站购售电功率结果图');

box off

figure

bar(-Pessswi(2,:));

hold on

bar(Pessbwi(2,:));

hold on 

legend('MG2向储能电站售电功率','MG2向储能电站购电功率');

xlabel('时间/h');

ylabel('功率/kW');

title('微网向储能电站购售电功率结果图');

box off

figure

bar(PGTwi(1,:),'stacked');%第i个微网的燃气轮机输出功率

% hold on

% bar(-Pgridwi(1,:),'stacked');%第i个微网从电网的购电功率

hold on 

plot(PECwi(1,:),'r-o','LineWidth',1.5)%第i个微网电制冷机消耗的电功率

hold on

plot(QACwi(1,:),'b-*','LineWidth',1.5)%第i个微网制冷机的输出制冷功率

hold on

plot(QGBwi(1,:),'g-*','LineWidth',1.5)%第i个电网的燃气轮机输出热功率

hold on

plot(PHXwi(1,:),'b-d','LineWidth',1.5)%第i个电网的换热装置输出制热功率

legend('燃气轮机输出功率','制冷机消耗的电功率','制冷机的输出制冷功率','燃气轮机输出热功率','换热装置输出制热功率');

xlabel('时间/h');

ylabel('功率/kW');

title('MG1调度结果');

box off

figure

bar(PGTwi(2,:),'stacked');%第i个微网的燃气轮机输出功率

% hold on

% bar(-Pgridwi(1,:),'stacked');%第i个微网从电网的购电功率

hold on 

plot(PECwi(2,:),'r-o','LineWidth',1.5)%第i个微网电制冷机消耗的电功率

hold on

plot(QACwi(2,:),'b-*','LineWidth',1.5)%第i个微网制冷机的输出制冷功率

hold on

plot(QGBwi(2,:),'g-*','LineWidth',1.5)%第i个电网的燃气轮机输出热功率

hold on

plot(PHXwi(2,:),'b-d','LineWidth',1.5)%第i个电网的换热装置输出制热功率

legend('燃气轮机输出功率','制冷机消耗的电功率','制冷机的输出制冷功率','燃气轮机输出热功率','换热装置输出制热功率');

xlabel('时间/h');

ylabel('功率/kW');

title('MG2调度结果');

box off

figure

bar(PGTwi(3,:),'stacked');%第i个微网的燃气轮机输出功率

% hold on

% bar(-Pgridwi(1,:),'stacked');%第i个微网从电网的购电功率

hold on 

plot(PECwi(3,:),'r-o','LineWidth',1.5)%第i个微网电制冷机消耗的电功率

hold on

plot(QACwi(3,:),'b-*','LineWidth',1.5)%第i个微网制冷机的输出制冷功率

hold on

plot(QGBwi(3,:),'g-*','LineWidth',1.5)%第i个电网的燃气轮机输出热功率

hold on

plot(PHXwi(3,:),'b-d','LineWidth',1.5)%第i个电网的换热装置输出制热功率

legend('燃气轮机输出功率','制冷机消耗的电功率','制冷机的输出制冷功率','燃气轮机输出热功率','换热装置输出制热功率');

xlabel('时间/h');

ylabel('功率/kW');

title('MG3调度结果');

box off

⛄ 运行结果

⛄ 参考文献

[1]吴盛军, 李群, 刘建坤,等. 基于储能电站服务的冷热电多微网系统双层优化配置[J]. 电网技术, 2021(10):3822-3829.

⛄ Matlab代码关注

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料