目录
主要内容
该程序实现的模型为综合能源系统利益分配优化调度,采用合作博弈方法,模型针对IES系统的P2G、电解槽、甲烷反应器、储氢罐、CHP和燃气锅炉等设备进行建模,实现基于合作博弈的IES协同优化调度模型,利用Shapley值法对合作剩余按贡献进行分配,程序采用matlab+cplex求解,程序注释清楚,方便学习!
部分代码
%% 决策变量定义 P_wt=sdpvar(1,T); %风机消纳功率 P_pv=sdpvar(1,T); %光伏消纳功率 V_gas=sdpvar(1,T); %P2G产生天然气体积 Q_co2=sdpvar(1,T); %制气厂消耗二氧化碳的量 P_H2=sdpvar(1,T); %电解池产生的氢气量 P_P2G=sdpvar(1,T); %P2G消耗电功率量 P_h2cha=sdpvar(1,T); %储氢罐充气功率 P_h2dis=sdpvar(1,T); %储氢罐放气功率 H_h2cha=binvar(1,T); %储氢罐充气变量 H_h2dis=binvar(1,T); %储氢罐放气变量 P_H2total=sdpvar(1,T); E_h=sdpvar(1,T); %储氢罐放气功率 P_CHPe=sdpvar(1,T); %CHP发电功率 P_CHPh=sdpvar(1,T); %CHP发热功率 P_CHPg=sdpvar(1,T); %CHP耗气功率 Q_co2out=sdpvar(1,T); %CHP燃烧产生的co2 P_be=sdpvar(1,T); %购电功率 P_bg=sdpvar(1,T); %购气功率 V_CB=sdpvar(1,T); %燃气锅炉进气量 Q_CB=sdpvar(1,T); %燃气锅炉发热量 P_MRh=sdpvar(1,T); %甲烷反应器输入氢功率 P_MRg=sdpvar(1,T); %甲烷反应器输出天然气功率 %% 风光模型 C=[C,0<=P_wt<=S_w]; C=[C,0<=P_pv<=S_p]; %% P2G设备运行模型 C=[C,Q_co2==alpha_co2*nita_ch4*load_H]; C=[C,P_H2==nita_H2*P_P2G]; %电解槽产生氢气功率 C=[C,P_H2total==P_H2+P_h2dis-P_h2cha]; C=[C,0<=P_H2<=500]; %% 甲烷反应器运行模型 C=[C,P_MRg==nita_ch4*P_MRh]; %电解槽产生氢气功率 C=[C,0<=P_MRg<=250]; %% 储氢罐模型 C = [C, H_h2dis + H_h2cha <= 1]; C = [C, 0 <= P_h2cha <=H_h2cha*50]; C = [C, 0 <= P_h2dis<= H_h2dis*50]; for t = 2:24 C = [C, E_h(t) == E_h(t-1) + 0.9*P_h2cha - P_h2dis/0.9]; end C = [C, 0.1 <= E_h(:)/200 <= 0.9]; C = [C, E_h(1)==E_h(24)]; %% CHP运行约束 C=[C,P_CHPe==P_CHPg*CHP_e]; C=[C,P_CHPh==P_CHPg*CHP_h]; C=[C,Q_co2out==P_CHPg*alpha_co2]; C=[C,0<=P_CHPg<=600]; for t=1:1:23 C=[C,-50<=P_CHPg(t+1)-P_CHPg(t)<=50]; end %% 燃气锅炉约束 C=[C,Q_CB==K_gas*nita_CB*V_CB]; C=[C,0<=Q_CB<=800]; for t=1:1:23 C=[C,-100<=Q_CB(t+1)-Q_CB(t)<=100]; end %% 功率平衡约束 C=[C,0<=P_be<=1000]; C=[C,0<=P_bg<=1000]; %电功率平衡 C=[C,load_e+P_P2G==P_wt+P_pv+P_CHPe+P_be]; %热功率平衡 C=[C,load_h==P_CHPh+Q_CB]; %气功率平衡 C=[C,load_g+P_CHPg+V_CB==P_MRg+P_bg]; %氢功率平衡 C=[C,load_H+P_MRh==P_H2+P_h2dis-P_h2cha];