阶梯碳下考虑 P2G-CCS 与供需灵活响应的 IES 优化调度(Matlab代码实现）

💥💥💞💞欢迎来到本博客❤️❤️💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。

⛳️座右铭：行百里者，半于九十。

📋📋📋本文目录如下：🎁🎁🎁

目录

💥1 概述

📚2 运行结果

🎉3 参考文献

🌈4 Matlab代码实现

---

💥1 概述

自“十二五”规划以来，我国在北京、天津、上海等 7 个省市启动了碳排放权交易试点工作。2021 年7 月 16 日，全国碳排放权交易市场启动上线交易，交易中心设在上海，登记中心设在武汉。对碳交易机制从理论研究逐步走向实际应用，说明碳交易市场对 IES 的降碳减排起到了良好的引导。文献［3］引入碳交易机制，以 IES 发电能源成本与碳交易成本之和最小为目标，建立了电-气互联 IES 低碳经济运行模型，提升了 IES 的低碳性；文献［4］引入光热电站，结合 P2G 等装置在碳交易机制下建立了一种IES 系统低碳优化运行方法；上述文献仅以固定碳价进行交易，对碳市场的引导作用有限；为此文献［5］引入阶梯型碳交易机制并建立了包含电转气和燃气轮机的电-热-气 IES 低碳经济调度模型；文献［6］在阶梯碳交易机制下，构建了园区综合能源系

统多阶段规划方法；上述文献验证了阶梯型碳交易机制能进一步约束 IES 碳排放。电转气（Power-to-Gas，P2G）技术可以将富余风电转化为天然气，较大程度地解决风电反调峰特性带来的弃风问题。文献［7］构建了含 P2G 设备的电-气-热多能源系统优化调度模型，分析了 P2G 消纳风电的经济效益；文献［8］在 IES 中细化电转气两阶段运行过程，提高了能源利用率；文献［9］考虑P2G 消纳多余水电，建立了水-电-气系统低碳鲁棒优化调度模型；然而，上述文献均未考虑 P2G 技术甲烷化过程中 CO2的来源。

在“双碳”背景下，为了进一步提高综合能源系统(Integrated energy system, IES)的低碳性和经济性，研究人员提出了一种园区级IES低碳经济优化调度模型。该模型的核心思想在于，通过引入氢能参与IES实现降碳减排，并采用阶梯式碳交易机制进一步约束碳排放，从而实现低碳经济。具体而言，该模型采用P2G-CCS(power to gas and carbon capture system，P2G-CCS)耦合，建立电转气和碳捕集系统的模型，将氢燃料电池与燃气轮机共同实现热电联产，并引入有机朗肯循环实现供应侧热电联产机组的灵活响应。同时，该模型还考虑了需求侧的多元负荷的可转移和可替代性，建立了需求侧灵活响应模型。最终，该模型以购能成本、碳成本、需求响应补偿成本、弃风成本之和最小为目标，建立IES优化调度模型，并采用MATLAB+YALMIP调用Cplex商业求解器求解。这一模型的提出，为园区级IES低碳经济的实现提供了新的思路和方法。

📚2 运行结果

部分代码：

%%平方线性化
T=24;
Num = 5 ; %分成5段折现
%设定P的上下限
Pup_BUY  =   1500;
Plow_BUY =   0;
%用变量名P2 当作 P^2
P2_BUY  =  sdpvar(1,T ); 
%确定出力区间约束
Constraints = [Constraints ,   Plow_BUY <= E_BUY <= Pup_BUY  ] ;
%那么端点数目就是 Num + 1
%找出所有端点的坐标（x,y）
x_Duandian_BUY = linspace(Plow_BUY , Pup_BUY , Num + 1 );
y_Duandian_BUY =  x_Duandian_BUY.^2 ;
%下面利用Num + 1个点，确定 Num 个分段的斜率， k = dita_y./dita_x
dita_y_BUY = y_Duandian_BUY(1,2:end) - y_Duandian_BUY(1,1:end-1 );
dita_x_BUY = x_Duandian_BUY(1,2:end) - x_Duandian_BUY(1,1:end-1 );
k_BUY =  dita_y_BUY./dita_x_BUY;
%发现，当 x_Duandian_i <=  P  <= x_Duandian_(i+1)时，
% P^2 可以用  y_Duandian_i + k_i*( P - x_Duandian_i  ) 进行近似替代
%难点在于两个，找出P所在分段的i序号
%为此，需要引入一个2进制变量标识位置
LuodianQuduan_BUY  =  binvar(T,Num );
%引入约束
%落点唯一，且在相邻两点之间
%相邻两点之间,意思是属于同一个分段，所以最好是将Num+1个点，做成Num个组合
Fenduan_low_BUY = x_Duandian_BUY(1,1:Num );
Fenduan_up_BUY  = x_Duandian_BUY(1,2:Num+1 ); 

🎉3 参考文献

文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。

[1]喻鑫,胡志坚,陈锦鹏等.阶梯碳下考虑P2G-CCS与供需灵活响应的IES优化调度[J/OL].武汉大学学报(工学版):1-14[2023-12-23].http://kns.cnki.net/kcms/detail/42.1675.T.20230518.0913.002.html.

🌈4 Matlab代码实现