目錄

1?主要內容

2?部分程序

3?實現效果

4?下載鏈接

---

1?主要內容

本程序是對《計及碳捕集電廠低碳特性的含風電電力系統源-荷多時間尺度調度方法》方法復現，非完全復現，只做了日前日內部分，并在上述基礎上改進升級為電熱綜合電源微網系統，未考慮節點系統。

本程序依據上述文獻模型實現一個虛擬電廠/微網多時間尺度電熱綜合能源系統低碳經濟調度模型，源側在碳捕集電廠中裝設煙氣旁路系統與溶液存儲器，形成碳捕集電廠綜合靈活運行方式進而與風電協調配合；荷側調用不同響應速度的價格型、激勵型需求響應資源克服多時間尺度下碳捕集電廠綜合靈活運行方式的局限，通過源荷資源協調優化，從而提高系統的低碳性能。其次，構建源荷協調的日前-日內兩階段低碳經濟調度模型，優化系統的負荷及分配計劃。

2?部分程序

%% 決策變量
% 電力源出力
GT_P = sdpvar(2,24,'full'); % 燃氣輪機電出力
P_w = sdpvar(1,24,'full'); % 風電機組出力
P_G = sdpvar(3,24,'full'); % 火電機組出力
EB=sdpvar(2,24,'full'); % 電鍋爐出力
% 熱力源出力
GT_H = sdpvar(2,24,'full'); % 燃氣輪機熱出力
EB_H=sdpvar(2,24,'full'); % 電鍋爐熱出力
% 天然氣
P_gas=sdpvar(2,24,'full'); % 天然氣需求
% 碳捕集相關
E_G=sdpvar(3,24,'full'); % 碳捕集機組產生的總碳排放
E_total_co2=sdpvar(3,24,'full'); % 機組捕獲的總碳排放
E_CG=sdpvar(3,24,'full'); % 儲液裝置提供的待捕集二氧化碳量
P_B=sdpvar(3,24,'full'); % 機組運行能耗
P_J=sdpvar(3,24,'full'); % 機組凈出力
V_CA=sdpvar(3,24,'full'); % 機組凈出力
V_FY=sdpvar(3,24,'full'); % 富液體積
V_PY=sdpvar(3,24,'full'); % 貧液體積
P_tran=sdpvar(1,24,'full'); % 系統可轉移電負荷
P_cut=sdpvar(1,24,'full'); % 系統可削減電負荷
P_DE=sdpvar(1,24,'full'); % 系統經過過需求響應后的電負荷
H_tran=sdpvar(1,24,'full'); % 系統可轉移熱負荷
H_cut=sdpvar(1,24,'full'); % 系統可削減熱負荷
H_DE=sdpvar(1,24,'full'); % 系統經過過需求響應后的熱負荷
gn=5;  
P_G_line= sdpvar(3,24,'full'); % 火電機組出力
%%?約束條件
C = [];  %約束條件初始
for t=1:24for i=1:3C = [C,         0<=E_CG(i,t),0<=P_B(i,t),E_G(i,t)==eg(i)*P_G(i,t), % 碳捕集機組產生的總碳排放   E_total_co2(i,t)==E_CG(i,t)+0.25*E_bata*eg(i)*(y1(i,t)-y2(i,t)), % 機組捕獲的二氧化碳總量0<=E_total_co2(i,t)<=P_yita*E_bata*eg(i)*P_G_max(i),   P_B(i,t)==P_lamda(i)*E_total_co2(i,t), % 機組運行能耗P_G(i,t)==P_J(i,t)+P_D(i,t)+P_B(i,t), % 機組輸出總功率P_G_min(i)-P_lamda(i)*P_yita*E_bata*eg(i)*P_G_max(i)-P_D(i)<=P_J(i,t)<=P_G_max(i)-P_D(i), % 碳捕集電廠凈出力范圍0<= P_w(t)<= P_prew(t), % 風電出力區間約束 sum(EB(:,t))+P_w(t)<=P_prew(t);P_G_min(i)<= P_G(i,t)<=P_G_max(i), % 火電機組出力約束];end
end
C=[C,min(sum(R_u),sum(P_G_max)-sum(P_G))>=0.01*max(P_DE),]; % 旋轉備用約束

3?實現效果

4?下載鏈接