简介
本程序参考EI论文《计及碳捕集电厂低碳特性的含风电电力系统源-荷多时间尺度调度方法》，文章综合考虑碳捕集电厂综合灵活运行方式的能量时移特性与快速调节特性、DR资源的多时间尺度特性与零碳排放特性，提出考虑源荷配合的电力系统日前-日内-实时三阶段低碳经济调度方法。。程序中算例丰富、注释清晰、干货满满，创新性很高！下面对文章和程序做简要介绍！
创新点及主要工作
1. 挖掘源荷两侧低碳资源并分析其低碳特性，源侧在碳捕集电厂中装设烟气旁路系统与溶液存储器，形成碳捕集电厂综合灵活运行方式进而与风电协调配合；
2. 荷侧调用不同响应速度的价格型、激励型需求响应资源克服多时间尺度下碳捕集电厂综合灵活运行方式的局限，通过源荷资源协调优化，从而提高系统的低碳性能。
3. 构建源荷协调的日前-日内-实时三阶段低碳经济调度模型，优化系统的负荷及旋转备用分配计划，并改善失负荷与弃风问题。
多时间尺度滚动框架
文中结果
程序结果
部分程序
%% 数据导入load A_IDRA_IDR_P_tran=A_IDR(1,:); % 日前调度需求响应--电力转移A_IDR_P_cut=A_IDR(2,:); % 日前调度需求响应--电力削减A_IDR_H_tran=A_IDR(3,:); % 日前调度需求响应--热力转移A_IDR_H_cut=A_IDR(4,:); % 日前调度需求响应--热力削减% 调度周期T=96;a=[0.00048,0.00031,0.0002]; % 成本系数ab=[16.2,17.3,16.6]; % 成本系数bc=[1000,970,700]; % 成本系数ceg=[0.91,0.95,0.98]; % 碳排放系数GT_P_max=[150,80];E_bata=0.81; % 碳捕集效率P_yita=[1.05,1.05,1.05]; % 再生塔和压缩机最大工作状态系数P_lamda=[0.269,0.269,0.269]; % 单位捕碳量能耗P_G_max=[400,455,200]; % 机组最大出力P_G_min=[200,120,100]; % 机组最小出力P_D1=0.25*0.5*[20,20,15]; % 碳捕集固定能耗endR_u_intra=[50,50,25]; % 机组上爬坡M_MEA=61.08; % M_MEA的摩尔质量M_co2=44; % 二氧化碳的摩尔质量sita=0.4; % 再生塔解析量CR=30; % 醇胺溶液浓度(%)rou_R=1.01; % 醇胺溶液密V_CR=[2000,2000,2000]; % 溶液储液装置容量K_R=1.17; % 乙醇胺溶剂成本系数fai=1.5; % 溶剂运行损耗系数 h_min=[200,120,100];h_max=[401,456,201];hl1=(h_max-h_min)./h_gn;hl2=zeros(3,h_gn+1);k_min=[-400,-455,-200];k_max=[-199,-119,-99];kl1=(k_max-k_min)./k_gn;kl2=zeros(3,k_gn+1);%% 决策变量% 电力源出力GT_P = sdpvar(2,96,'full'); % 燃气轮机电出力P_w = sdpvar(1,96,'full'); % 风电机组出力P_G = sdpvar(3,96,'full'); % 火电机组出力EB=sdpvar(2,96,'full'); % 电锅炉出力P_s_intra=sdpvar(1,96,'full'); % 失电负荷功率% 热力源出力GT_H = sdpvar(2,96,'full'); % 燃气轮机热出力EB_H=sdpvar(2,96,'full'); % 电锅炉热出力H_s_intra=sdpvar(1,96,'full'); % 失热负荷功率% 天然气P_gas=sdpvar(2,96,'full'); % 天然气需求% 碳捕集相关E_G=sdpvar(3,96,'full'); % 碳捕集机组产生的总碳排放E_total_co2=sdpvar(3,96,'full'); % 机组捕获的总碳排放% 电力需求响应P_tran_intra=sdpvar(1,96,'full'); % 系统可转移电负荷
部分内容源自网络，侵权联系删除！
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！
更多创新性预测及优化程序请点击公众号首页：《预测优化》|《综合能源》链接！
