%% 决策变量定义 P_wt=sdpvar(1,T); %风机消纳功率P_pv=sdpvar(1,T); %光伏消纳功率V_gas=sdpvar(1,T); %P2*生天然气体积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燃烧产生的co2P_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];
