程序名称：考虑共享储能接入的工业园区多类型负荷需求响应经济运行研究
实现平台：matlab-yalmip-cplex/gurobi
代码简介：提出一种含共享储能的园区多类型负荷需求响应经济运行模型。首先综合考虑3个园区的不同用户侧柔性负荷的可平移、可转移、可削减的负荷特性，进行多类型差异性负荷需求响应标准建模。然后将共享储能电站应用到园区负荷的经济优化调度中，通过协调各用户使用共享储能电站进行充电和放电的功率，实现园区群日运行成本最优。研究了各园区调度结果，以及不同共享储能服务费下的储能使用率和运行成本变化趋势。原创代码！附带参考文献，注释详细。代码非常极品，可拓展性高！适合电力系统优化调度需求响应和共享储能方向。
参考文献：《考虑用户侧柔性负荷的社区综合能源系统日前优化调度_刘蓉晖》《考虑需求响应的社区综合能源系统两阶段优化调度_刘蓉晖》《考虑需求侧协同响应的社区综合能源系统低碳经济调度_刘蓉晖》《考虑需求侧响应的含储能区域综合能源系统运行优化_李子林》《考虑含多种可控负荷的社区冷热电联供系统协同优化调度_刘蓉晖》《计及综合需求响应的社区综合能源系统优化调度_刘蓉晖》《基于共享储能电站的工业用户日前优化经济调度》
需求侧资源按照负荷参与需求响应的方式，将用户负荷分为4类：1）基础负荷：属于不可控负荷，对用户的需求完全响应，系统不能改变其用能方式和用能时间；2）可平移负荷：负荷供电时间可按计划变动，负荷需整体平移，用电时间跨越多个调度时段；3）可削减负荷：可承受一定中断或降功率、减少时间运行的负荷，根据供需情况对其进行部分或全部削减；4）可转移负荷：各时间段用电量可灵活调节，但要满足转移后整个周期负荷总量与转移前保持不变。
可平移负荷和可转移负荷都具有负荷供电时间按计划变动的特性，但二者也有区别：可平移负荷需整体平移，其用电时间不能中断且持续时间固定，用电时段所需的功率不可改变，如洗衣机、消毒柜等；可转移负荷比可平移负荷更灵活，用电时段的用电量可灵活调节，用电时段允许中断且持续时间不固定，只需满足转移前后的负荷需求总量不变，电动汽车是典型的可转移负荷，有序充电模式下电动汽车的充电时间、充电功率均可调整，但所需的总充电量不变。
共享储能模式下用户共用集中式储能，储能接受运营商统一运维管理，为同一配电网区域内的多个负荷聚合商提供相应服务。集中式储能无须用户独立投资建设储能设备，节省了投资建设成本。用户通过缴纳服务费换取共享储能使用权，各个负荷聚合商均可利用共享储能进行电能存取，提高电能利用率，保证电力供应稳定性，减少电网购电费用。各个负荷聚合商可以根据历史用电数据预测用电负荷，定制符合自身实际情况的共享储能使用计划，通过合理利用储能装置减小运行成本。
该代码基于MATLAB-YALMIP-Cplex/Gurobi平台，构建了考虑共享储能接入的工业园区多类型负荷需求响应经济运行模型，其核心实现原理围绕多类型柔性负荷精准建模、共享储能协同调度机制及园区群经济运行目标优化展开，具体如下：首先针对3个园区的差异化用户侧柔性负荷特性（可平移、可转移、可削减）进行标准建模——对于可转移负荷，定义转移功率变量（如P_tr）及二进制标识（U_tr），通过约束功率范围（-P_tr_max*U_tr ≤ P_tr ≤ P_tr_max*U_tr）、总转移量平衡（sum(P_tr)=0，确保转入转出相等）及最短连续转移时间（sum(U_tr(t:t+T_tr_min)) ≥ T_tr_min*(U_tr(t)-U_tr(t-1))，如园区1转移负荷需至少持续3小时）保障其调节合理性；对于可削减负荷，定义削减功率变量（P_cut）及二进制标识（U_cut），通过功率上下限约束（0.5*P_cut_max*U_cut ≤ P_cut ≤ P_cut_max*U_cut）、最短/最长连续削减时间及最大削减次数（如园区1最多削减7次、每次最长6小时）限制其调节边界；对于可平移负荷，通过二进制平移标志（Temp_Pshift1）和平移电量变量（PPshift1）建模，约束平移总量（如sum(Temp_Pshift1(5:21))=2表示需平移2个连续时段）、连续性（如sum(Temp_Pshift1(i:i+1)) ≥ 2*(标志变化)确保平移时段连续）并关联平移量（如PPshift1=25*Temp_Pshift1表示每个平移时段转移25单位电量）。其次，将共享储能作为园区群公共调节资源，建立其充放电功率（P_ESS_char/P_ESS_dis）、状态标识（U_ESS，避免同时充放电）、储能电量（E_ESS）及容量约束模型：通过初始与末期电量相等（E_ESS(1)=0.5*E_ESS_max、E_ESS(24)=E_ESS(1)）保证日内循环，结合效率因子（0.9）建立时段间电量动态方程（E_ESS(t) = 0.9*E_ESS(t-1) + P_ESS_char(t-1) - P_ESS_dis(t-1)），并限制充放电功率不超过最大值（0≤P_ESS_char≤U_ESS*P_ESS_max、0≤P_ESS_dis≤(1-U_ESS)*P_ESS_max），且总充放电量等于三园区交互量之和（P_ESS_char-P_ESS_dis=P_ESS3+P_ESS2+P_ESS1）以实现协同调度。最后，基于功率平衡关系（如园区1：P_bs1 = PL0_1 + P_tr + P_tr_2 - P_cut + P_ESS1，即电网交互量等于固定负荷加转移负荷、减削减负荷再加储能交互量）构建经济运行目标，以园区群日运行成本最优为目标函数——第一部分为分时电价下的购售电成本（sum((P_bs1+P_bs2+P_bs3).*C1)，C1为峰谷平分时电价），第二部分为共享储能服务费（随服务费系数k变化，即0.02*k*sum(P_ESS_char+P_ESS_dis)），通过求解混合整数线性规划问题得到各园区柔性负荷调节策略、共享储能充放电计划及电网交互功率的最优解，并进一步分析不同共享储能服务费下储能使用率与运行成本的变化趋势，最终实现多园区柔性负荷与共享储能的协同优化，降低园区群整体运行成本。
专题推荐：代码分享，典藏级代码，原创代码改进，论文思路，免责声明（点击即可跳转）【New Power System预测和优化理论】
代码获取方式：【原创代码改进】考虑共享储能接入的工业园区多类型负荷需求响应经济运行研究
运行结果展示
