👨‍🎓个人主页：研学社的博客 

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

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

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

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

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码+数据

---

💥1 概述

文献来源：

综合以上研究现状,目前共享储能的研究仍处于初始阶段，现有的基础集中在共享储能服务运营商运行决策和商业模式上，用户负荷类型单一，对用户使用共享储能的充放电行为和用户与共享储能间费用结算分析较少,同时具体将共享储能应用于冷热电联供型多微网系统优化配置的相关研究较少。基于以上问题,本文提出基于储能电站服务的冷热电联供型多微网系统双层优化配置,建立同时考虑储能电站容量配置和冷热电联供型多微网系统优化运行的双层规划模型(bi-level programming,BLP)
首先提出储能电站服务模式;然后建立考虑规划和运行的双层优化配置模型，采用KKT条件和 Big-M法将双层模型转换为单层线性模型求解;最后以仿真算例进行分析，验证所提模型的经济性和有效性。本文的创新点总结如下:

1）提出储能电站服务模式，分析其运行方式和盈利机制;

2）将该模式应用于冷热电联供型多微网系统中,实现用户内部多种能源间互补共济;

3）采用双层模型建立同时综合考虑储能电站运营商和多微网系统利益的储能决策。

                  

📚2 运行结果

 

 

 

 

 

 

 

 

 

 部分代码：

%Big-M法处理割平面约束
for i=1:3
    for t=1:24
        Con=[Con,
            0<=u1itmin(i,t)<=M*v15(i,t),
            0<=PGTwi(i,t)-5<=M*(1-v15(i,t)),
            0<=u1itmax(i,t)<=M*v16(i,t),
            0<=3000-PGTwi(i,t)<=M*(1-v16(i,t)),
            0<=u2itmin(i,t)<=M*v17(i,t),
            0<=QACwi(i,t)-5<=M*(1-v17(i,t)),
            0<=u2itmax(i,t)<=M*v18(i,t),
            0<=4000-QACwi(i,t)<=M*(1-v18(i,t)),
            0<=u3itmin(i,t)<=M*v19(i,t),
            0<=PECwi(i,t)-5<=M*(1-v19(i,t)),
            0<=u3itmax(i,t)<=M*v20(i,t),
            0<=4000-PECwi(i,t)<=M*(1-v20(i,t)),
            0<=u4itmin(i,t)<=M*v21(i,t),
            0<=QGBwi(i,t)-5<=M*(1-v21(i,t)),
            0<=u4itmax(i,t)<=M*v22(i,t),
            0<=4000-QGBwi(i,t)<=M*(1-v22(i,t)),
            0<=u5itmin(i,t)<=M*v23(i,t),
            0<=PHXwi(i,t)-5<=M*(1-v23(i,t)),
            0<=u5itmax(i,t)<=M*v24(i,t),
            0<=4000-PHXwi(i,t)<=M*(1-v24(i,t)),
            0<=u6itmin(i,t)<=M*v25(i,t),
            0<=Pgridwi(i,t)<=M*(1-v25(i,t)),
            0<=u6itmax(i,t)<=M*v26(i,t),
            0<=4000-Pgridwi(i,t)<=M*(1-v26(i,t)),
            0<=u7itmin(i,t)<=M*v27(i,t),
            0<=Pessswi(i,t)<=M*(1-v27(i,t)),
            0<=u7itmin(i,t)<=M*v28(i,t),
            0<=(4000*Usalewi(i,t)-Pessswi(i,t))<=M*(1-v28(i,t)),
            0<=u8itmin(i,t)<=M*v29(i,t),
            0<=Pessbwi(i,t)<=M*(1-v29(i,t)),
            0<=u8itmin(i,t)<=M*v30(i,t),
            0<=(4000*Ubuywi(i,t)-Pessbwi(i,t))<=M*(1-v30(i,t)),
            0<=u9itmax(i,t)<=M*v31(i,t),
            0<=(1-Ubuywi(i,t)-Usalewi(i,t))<=M*v31(i,t),
            ];
    end
end

🎉3 参考文献

部分理论来源于网络，如有侵权请联系删除。

[1]吴盛军,李群,刘建坤,周前,汪成根.基于储能电站服务的冷热电多微网系统双层优化配置[J].电网技术,2021,45(10):3822-3832.DOI:10.13335/j.1000-3673.pst.2020.1838.

🌈4 Matlab代码+数据