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

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

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

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

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

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

目录

💥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代码+数据