003文章解读与完整程序——电网技术EI\CSCD\北大核心《基于主从博弈理论的共享储能与综合能源微网优化运行研究》已提供下载资源

下载资源链接

摘要：综合能源微网与共享储能在提高系统能源利用率方面具有显著效益，逐渐成为目前研究热点，如何建立一套共享储能背景下综合能源微网优化运行模型是当下亟待解决的问题。文章首先介绍系统运行框架，分析系统内各利益体的功能。然后，分别针对微网运营商、共享储能服务商以及用户聚合商建立优化运行模型。分析微网运营商与用户聚合商间的博弈关系，提出共享储能背景下微网运营商与用户聚合商间的Stackelberg博弈模型，并证明Stackelberg均衡解的存在性与唯一性。最后，在MATLAB平台上进行算例仿真，通过Yalmip工具与CPLEX求解器进行建模与求解，利用启发式算法与求解器结合的方法优化微网运营商与用户聚合商的策略。结果表明，所提模型不仅能有效权衡微网运营商与用户聚合商的利益，也能够实现用户聚合商与共享储能运营商的收益双赢，采用的求解算法能够保护微网运营商与用户聚合商间的数据隐私。

部分代码展示：

%用户侧收益作为目标函数
function [P_MT,F_user,F_share,Eload,Hload,ES,P_h,Prl,P_buy,P_sell] = computeObj(x,load_e,load_h,P_PV,pe_grid_B)
P_MT=sdpvar(1,24,'full');%微燃轮机输出电功率
P_buy=sdpvar(1,24,'full');%用户向运营商买电电量
P_sell=sdpvar(1,24,'full');%用户向电网卖电电量
ES=sdpvar(1,24,'full');%储能余量
%% 需求侧定义变量
%电负荷：固定、可平移、可消减负荷、电替热
%热负荷：固定、可消减、热被电替
Pfl=sdpvar(1,24,'full');%可平移电负荷量
eload=0.8*(load_e); %消减之后的电负荷量
Pcl_h=sdpvar(1,24,'full');%可消减热负荷量
Prl=sdpvar(1,24,'full');%电制热设备供热量
P_h=sdpvar(1,24,'full');%微网运营商供热量
char=sdpvar(1,24,'full'); %充电功率
char_sign=binvar(1,24,'full');%充电标志 
dischar=sdpvar(1,24,'full'); %放电功率
dischar_sign=binvar(1,24,'full');%放电标志
%a 、b 、c为用户聚合商的用电效用函数的参数
a=-0.05;b=4;
%微燃电机系数
MT_e=0.4; %发电效率
MT_h=0.8;   %制热效率
MT_hh=0.05;%散热损失率

%约束条件
C =[];
%% 共享储能服务商
ESS_max=1350;ESS_char=0.95;ESS_dischar=0.95;%电储能容量/充电/放电
SOC0=0.5;
ES(1,1)=SOC0*ESS_max;%%初始电量
 for t=2:25  %在一个周期内的充放电功率为零
     C=[C,(ES(mod(t-1,24)+1)==(ES(mod(t-2,24)+1)+(ESS_char*char(mod(t-2,24)+1)-(1/ESS_dischar)*dischar(mod(t-2,24)+1))))];
     C=[C,ES(1,1)==SOC0*ESS_max];
 end
for i=1:24
     C=[C,300<=ES(1,i)<=1350]; %储能容量约束
     C=[C,0<=char(1,i)<=50*char_sign(1,i)]; %充放电约束
     C=[C,0<=dischar(1,i)<=50*dischar_sign(1,i)];  
end
%蓄电池充放电约束
 for i=1:24
     C=[C,char_sign(1,i)+dischar_sign(1,i)<=1];    
 end
  C=[C,10<=sum(char_sign(1,1:24)+dischar_sign(1,1:24))<=20];%考虑寿命

效果展示：