目錄
主要內容???
??部分代碼???
??結果一覽???
下載鏈接
主要內容???
該程序實現的模型為綜合能源系統利益分配優化調度,采用合作博弈方法,模型針對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:24C = [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:23C=[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:23C=[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];