[종합에너지] 탄소포집발전소의 저탄소 특성과 수요반응을 고려한 집단에너지 시스템의 다시간 규모 급전 모델

2024-07-12

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),]; % 旋转备用约束

기술나눔