技術共有

한어Русский языкEnglishFrançaisIndonesianSanskrit日本語DeutschPortuguêsΕλληνικάespañolItalianoSuomalainenLatina

目次

1 主な内容

2部構成の手順

3 成果を上げる

4つのダウンロードリンク

---

1 メインコンテンツ

この手順は、「炭素回収型発電所の低炭素特性を考慮した風力発電所」のレビューです。パワーシステム「ソースロードマルチタイムスケールスケジューリング方式」方式を再現するが、完全に再現するわけではない。前日と当日の部分のみを行い、これをベースに改良し、電力・火力一体型電力にアップグレードする。ノードシステムを考慮せずにマイクログリッドシステムを提供します。

このプログラムは、上記の文献モデルに基づいて仮想発電所/マイクログリッドを実装します。時間スケール電気と熱の統合エネルギーシステムの低炭素経済的供給モデルでは、供給側は炭素回収発電所に排ガスバイパスシステムと溶液貯蔵所を設置し、炭素回収発電所の包括的かつ柔軟な運用モードを形成して、風力発電と調整し、負荷側は異なる応答速度を要求します。 価格ベースおよびインセンティブベースのデマンド レスポンス リソースは、複数の時間スケールにおける炭素回収発電所の包括的かつ柔軟な運用モードの制限を克服し、低炭素パフォーマンスを向上させます。ソース負荷リソースの調整と最適化を通じてシステムの最適化を実現します。第二に,システム負荷と配電計画を最適化するために,ソース負荷調整のための前日2段階低炭素経済配電モデルを構築した。

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:24
    for i=1:3
    C = [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つのダウンロードリンク