P2G与碳捕集协同的多能微网低碳经济调度建模实战

做了这么多年微网优化调度有个感受一直很深P2G和碳捕集这两项技术单拎出来都是热门方向但真正把它们塞进同一个调度框架里算经济账的却不多。这次的项目标题是考虑P2G与碳捕集机组的多能微网低碳经济调度探索听起来挺学院派实际落地就是回答三个问题微网里的富余可再生电力怎么变成可存储可运输的气燃气机组排出的CO2怎么变成P2G的原料而不是直接排到大气里这一套电-气-碳耦合下来到底省了多少钱、降了多少碳这篇博文就把建模思路、参数取值、求解工具和踩坑记录从头到尾捋一遍希望能给做园区微网、综合能源系统或碳管理方向的同行提供一份能直接参考的实操笔记。1. 立项思路为什么要把P2G和碳捕集放进同一个调度框里算1.1 多能微网的三笔账电、气、碳多能微网的调度表面上是电功率平衡实际上是三笔账一起算电是一笔气是一笔碳又是一笔。传统微网调度只盯前两笔把碳排放当成一个事后的统计数字最多加个碳约束就算完事。但引入P2G和碳捕集机组后碳不再是排放出来的废物而是参与了能量转化的物料。也就是说CO2从环境负担变成了系统里的中间产物能不能用得起、用得好直接决定了调度方案的优劣。我接手这个项目时园区微网的基本盘是1台燃气轮机、若干屋顶光伏、两台风机、一个电储能还有用户的电热负荷。原来调度目标函数里无非是购电成本、燃料成本和设备运维成本碳排放相关的内容就是每度电折算多少碳乘以碳价。这种粗放处理在碳价低的时候没什么问题但一旦涉及碳捕集就必须把CO2的捕捉、存储、消纳全过程纳进去否则系统里出现的CO2流量就无法闭环。项目标题里的探索两个字本质就是这三笔账联动后怎么做最优决策。1.2 传统微网调度方案漏掉的那块碳资产传统调度的目标函数里没有碳捕集也没有P2G所以它天然会导致两个结果。第一富余风电和光伏在负荷低谷期只能弃掉因为没有足够的储能容量去消纳第二燃气轮机产生的CO2全部直排不但要承担碳配额超标处罚还白白浪费了碳源。把P2G和碳捕集加进来之后逻辑就变了弃电时段用P2G电解水制氢再把氢和捕集来的CO2甲烷化产出天然气燃气轮机发电的同时烟气里的CO2被捕集装置回收一部分送往P2G一部分临时存储。用生活类比来说传统微网像是一间只记账不回收的厨房燃气灶燃气轮机做饭产生废气剩饭剩菜弃风弃光直接倒掉。引入P2G和碳捕集后厨房变成了循环系统废气被净化器吸收变成燃料原料剩饭被发酵罐变成新能源。这个项目的核心思路就是抛弃线性消耗的旧模型改用循环转化的新模型。只有把碳捕集机组的捕集能耗、P2G装置的电转气效率、CO2存储罐的容量约束全部耦合进调度模型才能真正回答减碳能不能省钱这个老板最关心的问题。2. 建模前的关键认知P2G与碳捕集的原理和耦合点2.1 P2G装置不是电制氢那么简单很多人一提P2G就想到电解水制氢实际上一个完整的P2G系统分两段第一段是电解槽把水和电变成氢气和氧气第二段是甲烷化反应器把氢气与CO2在催化剂作用下合成甲烷。为什么必须做第二段因为氢气直接注入天然气管网有限制目前管网掺氢比例一般在5%上下高了影响终端设备而甲烷与现有燃气基础设施完全兼容。这直接导致一个调度上的时间节奏问题电解槽要响应电价波动甲烷化要响应CO2的可获得性。两者不是简单的前后串联中间最好有一个储氢缓冲罐。我在这项目里把P2G建模成电解槽储氢甲烷化三个单元而不是一个大黑箱。为什么要拆这么细因为电解槽的启停响应很快分钟级甲烷化反应器启停比较慢至少30分钟到1小时如果不解耦建模调度优化器会给出一个每秒都在启停的不可行解工程上完全实现不了。关键参数取舍碱性电解槽的效率我用0.65-0.75甲烷化效率取0.75-0.80折算下来整个P2G系统从电到天然气的电转气效率大约0.48-0.60。凡是看到宣传里说电转气效率90%的基本都是只算了理想状态单段效率没有把电解槽散热、压缩机耗电、甲烷化冷却损失算进去。调度模型里的效率参数宁可保守也别激进否则优化出来的方案在真实运行中会发飘。2.2 碳捕集机组的净出力陷阱捕集能耗与净出力碳捕集机组和普通燃气机组最核心的区别是什么是净出力而非总出力。传统机组调度你说燃机输出多少电那就是多少电。碳捕集机组不是它有一个固定的能耗惩罚捕集CO2需要热能胺法再生塔耗蒸汽和电能压缩机、泵。我在DCS系统上看到的实际数据是捕集单位CO2大约需要0.3-0.5 kWh电能另加一定量蒸汽。所以模型里一定要把机组出力拆成总发电出力上网净出力捕集能耗分摊。这个拆法直接决定调度结果的经济性。有的文献把捕集能耗设成一个常数我实际跑下来发现不对——捕集能耗其实随捕集率变化捕集率越高单位能耗越大有个边际递增的关系。线性模型里可以分段处理比如捕集率85%时能耗系数取0.3590%时取0.4595%时取0.55。你在文献里也许看不到这个细节但在工程中这条边际曲线是CO2成本和发电效率之间的真实博弈。另一个容易踩坑的地方是捕碳但不存碳、不卖碳的工况。如果碳捕集系统捕集下来的CO2没有去处就只能临时存入CO2缓冲罐。缓冲罐容量有限满了以后要么停止捕集失去减排意义要么启动CO2压缩储存设备额外耗电。所以我在约束里给CO2存储罐设了容量边界并在目标函数里加了一个碳捕集后若无法消纳则产生惩罚的虚拟项这样模型才会在最优解里主动控制捕集节奏而不是盲目满负荷捕碳。2.3 碳闭环的工程可行性管道、储罐与时间节奏在园区微网里做碳闭环最容易被忽略的是CO2的供需时间错配。白天光伏大发P2G刚好有大量电解制氢能力但它需要的CO2来自燃气轮机的烟气捕集燃气轮机恰恰在白天负荷高、发电多、碳排放多。这个节奏天然匹配。晚上负荷低谷燃气轮机出力小捕集的CO2量也小P2G即使想多产甲烷也没有原料。所以P2G在晚间的产能受限不是电价不够便宜而是缺碳源。我在项目组里给的工程方案是CO2日存储罐容量按白天捕集高峰多余量夜间甲烷化需求折中设计。日存储罐的好处是不用搞大规模地质封存那套复杂工程一个中压球罐就能搞定投资可控。模型里CO2存储罐的物理方程就是三个状态进流量捕集输出、出流量P2G消耗、罐内储量。这看起来简单但在优化模型里必须作为连续状态变量参与约束否则CO2的物料守恒不成立可能出现今天捕了一吨碳明天P2G用掉了两吨碳这种数学上成立、物理上荒谬的结果。3. 低碳经济调度模型的搭建实操3.1 目标函数成本要算全碳价要递进模型的目标函数我最终定为四块成本之和购电购气成本从上级电网购电、从市政管网购天然气的费用注意P2G产气如果并网售卖要记为负成本。设备运维成本燃气轮机、P2G电解槽、甲烷化反应器、碳捕集装置的运维费率通常按启停次数和出力水平折算。碳交易成本采用阶梯式碳价机制不是一碳一价。免费配额用完之后超出部分价格递增。这个递进关系非常重要我从实际算例里观察到单碳价模型会让模型在高碳价环境下过度捕碳阶梯碳价则会让模型在刚好消化配额附近寻找更经济的平衡点。碳排放环境惩罚超配额排放部分在碳交易市场之外再加一个环境损害系数用来模拟政策监管带来的隐性成本。目标函数写成最小化各项之和但有一个陷阱运维成本和启停次数不是简单的线性关系。燃气轮机频繁启停会增加寿命损耗这个成本如果不进模型优化器会为了省燃料费而让机组频繁启停。我的解决方法是给燃气轮机加最小开停机时间约束再给P2G电解槽加启停次数上限约束。这种二次约束在MILP框架里很好表达但对普通经验试算场景你也能在启发式算法里用罚函数方式近似处理。3.2 约束条件六类约束一个都不能少具体约束我按六个模块整理电功率平衡约束所有电源出力购电储能放电电负荷P2G电解耗电捕集能耗储能充电。这里面P2G耗电和捕集能耗都是隐性电负荷很多初版模型容易把这两项漏掉导致优化结果里P2G在电价高峰时段疯狂耗电实际运行却无法执行。气平衡约束燃气轮机消耗天然气P2G产甲烷注入购气终端气负荷锅炉气耗。如果园区气负荷很小P2G产气没有消纳通道模型就会提示无解。这时候要么加储气罐要么允许P2G产气外售。热平衡约束燃气轮机余热和溴化锂机组的热输出满足用户热负荷注意热电比约束与机组出力水平的非线性关系我用分段线性拟合。碳捕集与CO2流量约束捕集量不超过烟气含碳总量CO2储存容量在上下限内P2G甲烷化的CO2消耗不超过可用CO2。储能约束电储能SOC连续充放电功率限制充放电切换最小间隔时间。爬坡与最小启停约束燃气轮机爬坡率±0.3 MW/min最小开机2小时最小停机2小时电解槽启动次数每日不超过5次。3.3 求解工具选型MILP是首选启发式做验证调度模型我做了线性化处理目标函数里的阶梯碳价用0-1变量加线性段模拟机组热电比关系用分段线性函数储能充放电效率用净效率折算。整体模型是标准的MILP规模大概3000多个连续变量和200多个0-1变量。我用的求解器是Gurobi也可以通过YALMIP在MATLAB里调用CPLEX收敛时间几秒到几十秒不等取决于场景数量。为什么不用粒子群或鲸鱼算法做主求解我在小规模算例里对比过启发式算法在离散变量多、约束耦合强的时候容易陷入局部最优而且没法给出最优性间隙。MILP求解器可以明确告诉你解的最优性gap是0.1%还是5%。对于工程调度5%的gap已经可以接受但如果你用PSO你永远不知道当前解离最优有多远。我的习惯是MILP求主调度方案启发式算法的结果当作破局思路——比如在MILP求解时间长到不可接受时用启发式出一个初值或热启动点。3.4 参数设置心得这几个数决定了模型能不能看参数这一块我觉得比公式本身更容易让模型失真。我列一个自己实际调参时的基准表参数取值说明电解槽效率0.68碱性电解槽实际效率含辅助能耗甲烷化效率0.78固定床反应器含压缩能耗P2G整体效率0.53电解和甲烷化连乘后乘0.9的整流损耗碳捕集率上限0.90单级胺法吸收工艺超过则能耗剧增捕集单位电耗0.38 kWh/kg CO2含压缩机功耗随捕集率递增碳交易基准价60元/吨按国内碳市场近两年区间中值阶梯碳价增量50元/吨超配额5%-15%后逐级递进光伏预测误差±8%典型日晴天用于场景法区间覆盖这些参数不是拍脑袋是我从设备厂家的技术手册和现场DCS报表里捞出来的。强烈建议你先做一次灵敏度分析再正式跑调度把P2G效率从0.5变到0.6看调度结果是否发生质变。如果某个参数的小幅波动会让优化方案翻盘那这个参数就是模型的痛点值得花更多精力去核准。我在这个项目中最大的参数教训是碳价不能拍脑袋定一个一定要结合最新的地方碳配额政策否则算出来的最优方案在真实碳交易市场里根本不经济。4. 典型算例从仿真结果里读出的三条规律4.1 典型日场景设置我拿夏季园区典型日来展示效果因为夏季光伏出力高、空调电负荷高弃光风险和碳排放量都比较有代表性。风电装机8 MW光伏装机5 MW燃气轮机6 MW电储能4 MW/8 MWhP2G容量1.5 MWCO2捕集装置按燃机排烟6 t/h设计CO2缓冲罐容量60 t。电网购电价用峰平谷三段曲线天然气价格波动0.5%/时段。这个场景在剥离P2G和碳捕集后初始碳排放约42 t/天弃光量约5.8 MWh。加上P2G和碳捕集后系统调度维度从电气扩到电气碳模型运算量有所上升但Gurobi求解时间仍能控制在1分钟内。4.2 三种方案对比协同优化的效果不是112我做了三个方案对比方案A是无碳捕集无P2G传统微网、方案B是只加碳捕集、方案C是碳捕集P2G协同。方案A成本购电6.4万元购气3.8万元碳交易成本1.2万元超配额总运行成本11.4万元碳排放42 t/天。方案B成本购电7.0万元捕集能耗使购电量增加购气3.5万元捕集运维0.6万元碳交易成本0.4万元未超配额总运行成本11.5万元。碳排放低于20 t/天但成本反而微增因为捕集能耗的电费比省下的碳费更多。方案C成本购电6.0万元P2G参与削峰填谷减少了高峰购电购气1.2万元P2G产气替代部分购气捕集运维P2G运维1.1万元碳交易成本0.1万元总运行成本8.4万元碳排放约12 t/天。方案C的单日成本比方案A低了约26%碳排放下降约70%。这个降碳幅度并非来自捕集率的堆高而是因为P2G把原本弃掉的光伏电力转化成了天然气替代了外部购气同时捕集下来的CO2直接消纳进P2G形成了闭环。方案B的教训就是典型的有减碳无经济只捕不用CO2捕了却滞销等于把碳费换成了电费。4.3 经济性拐点P2G何时会赔本赚吆喝最优调度结果里有一个值得警惕的规律P2G并不是全天开机而是集中在光伏大发时段和电价低谷时段。一旦P2G开机时段超过8小时它的边际收益就开始下降因为储氢罐和CO2缓冲罐的容量限制了连续生产。我算出当前参数组合下P2G的最佳日运行时间为6到8小时对应产甲烷量约1.2万立方米刚刚好够园区自身的燃气轮机补气消耗。这也解释了为什么很多示范项目里P2G长期处于演示多于投产的状态——如果电价低谷窗口不够宽、弃电不够多P2G的运行时间不足固定成本摊销不下去。我给园区算过一笔账如果P2G装置满负荷年运行时间低于1500小时单位产气成本就会超过管道天然气市场价经济性就完全站不住。所以调度模型不应该预设P2G必须用满而是要在目标函数里让模型自己决定什么时候开机。模型算出来的P2G运行窗口就是决策层判断投资是否合理的依据。4.4 灵敏度分析碳价和风光渗透率谁更影响调度我做了一组灵敏度扫描把碳价从40元/吨一路调到120元/吨P2G出力没有线性增长而是呈台阶式变化当碳价低于60元/吨P2G几乎不启动高于60元/吨后P2G出力迅速升到上限再往上增加的产出有限主要瓶颈在CO2捕集量和储罐容量。这个台阶效应告诉我们碳价政策低于某个阈值时光靠激励机制根本推不动P2G投资。风光渗透率的影响也很有意思。当风电光伏装机占比从30%提高到50%弃电量和P2G启停次数同步增加但系统的碳排放并不会线性下降。原因是渗透率提高后燃气轮机为保障供电可靠性仍然需要维持最小出力碳排放基数压不下来。真正压碳的关键是给燃机配上更高捕集率下的碳捕集能力而不是单纯扩大光伏板面积。如果你在项目汇报中只讲配了多少新能源领导会问碳降了多少这时候碳捕集的协同效果才是硬指标。5. 实际工程中的常见问题与排查实录5.1 问题1模型求解卡死提示不可行或无界搭建第一版MILP模型时最容易出现不可行的问题。我排查时最先检查六个约束模块里的物料平衡。一个典型错误是P2G的CO2消耗量大于捕集量加储罐可用量这时候约束无解。解决方法是检查CO2存储罐初值设置是否合理我吃过一个亏把CO2储罐初始储量设成与容量上限相同而捕集量又不够补充系统库存不足导致甲烷化出力被卡死在0。排查顺序建议先查电平衡再查气平衡最后查碳平衡。三套平衡里碳平衡的排查最费事因为碳捕集量、CO2存储变化率、甲烷化消耗量三个变量互相耦合。把这三个变量单独打印出来按24个时段画成折线图一眼就能看出哪个时段碳流出现断点。这种调试习惯和查电路板虚焊是一个道理先从最细的支路查起。5.2 问题2P2G在凌晨莫名开机经济性却很差第一次跑出完整调度结果我看到优化器让P2G在凌晨2点开机一度很兴奋以为它发现了什么深夜低谷电的套利空间。结果细看发现凌晨时段电网购电价确实低但P2G需要的CO2捕集量几乎为零因为燃气轮机在夜间接近停机。模型为了利用低谷电而强行产甲烷导致CO2储罐库存被一夜清空白天捕集的高峰时段反而出现CO2储罐容量不足、捕集系统被迫降负荷的情况。这不是模型逻辑错误而是缺少CO2储罐库存变化的阶梯惩罚仓库空的时候应该给储存成本打折鼓励多捕满的时候应该给储存压力加大惩罚鼓励多用。我在目标函数里加了一个储罐稀缺性惩罚项按罐内储量与容量上限的比值动态调整模型立刻收敛到一个更合理的方案P2G白天光伏大发时段开机晚间几乎不动作。5.3 问题3电解槽频繁启停导致设备寿命隐患优化模型本身不会考虑设备疲劳度它只会看成本和约束。所以不加约束时P2G电解槽可能在一个调度日内开关8次这在物理上很伤催化剂和隔膜。我给电解槽设了最小连续运行时间1小时和日启停次数≤5次的整数约束后模型照样能找到经济可行解只是开关次数少了总成本增加了不到2%。这2%相当于给设备买了保险非常值得。燃气轮机同样有此问题。某些场景下模型让燃机在午间光伏大发时段停机、晚间电价高峰再启动从购电成本看是最优的但一次启停的检修成本可以抹平大半天的燃料节省。建议在目标函数里给启停动作加一个小的固定成本项不要设为0因为现实中任何设备启停都有代价。5.4 问题4CO2缓冲罐容量设得太大投资成本爆炸有的同事在设计CO2储罐容量时习惯按照最大捕集量×最长储存时间来选型结果罐体从60 t一路加到200 t投资翻了三四倍。我算过一笔账储罐容量只要覆盖捕集高峰到甲烷化消耗高峰之间的净积累量即可典型值就是日捕集总量的1/3到1/2。容量继续加大边际效益下降得厉害因为捕集量和甲烷化消耗的节奏本身有同步性。另外一个工程细节是CO2储罐要考虑液态储存还是气态储存。微网场景下我用中压气态存储1.0-1.5 MPa单位设备成本比液态储存低很多虽然密度低一点但对日调度尺度够用。液态储罐的冷却能耗反而会增加系统的碳排放得不偿失。这些工程决策在纯仿真模型里体现不出来但实际落地阶段限制往往比模型约束更严格。5.5 排查技巧速查表单位统一是硬道理症状可能原因排查方法无解/不可行CO2物料不平衡或储罐初值异常单独输出碳平衡曲线检查CO2储罐日末储量是否回落到初值求解时间异常长0-1变量过多或分段线性化段数过多压缩碳价梯级数简化热电比分段段数减到误差可接受即可P2G出力与预期相反碳价参数过低或CO2可用量被约束卡死调碳价重跑灵敏度看P2G出力台阶变化碳排放量不降反升捕集能耗导致燃机多发电多排碳检查捕集率约束尝试围栏式降低燃机最小出力配合电储能结果数值波动大量纲不统一MW与kWh、吨与kg混用全模型统一基础单位我统一用MW和MWhCO2流量用t/h量纲问题其实排在我这里所有排查项里的第一位。第一次建模时我把光伏预测功率用了kW购电价格却用了元/MWh结果目标函数里出现了一个数量级106的隐性系数优化器给出的方案让储能疯狂充放电来赚差价。这个错误耽误了三天最后是逐项打印单位量纲才抓住。所以我的建议是模型做完后先跑一个全零输入的退化测试再把单个设备模型拿出来单独验证最后再耦合全系统这样可以省掉大量排错时间。6. 几点经验体会这套探索做下来我最大的感受是减碳技术不是堆出来的是算出来的。碳捕集与P2G的协同效应本质上是一道时间节奏题什么时候捕碳、什么时候用氢、什么时候产甲烷、什么时候储气、什么时候充放电每一步都需要调度模型去平衡电价格、气价格、碳价格三套信号。模型的价值不是替代人的判断而是帮你在成千上万种组合里快速排除掉那些看似合理、实际不经济的方案。如果你正打算在园区微网里引入P2G或碳捕集我的建议是先做一个不投资版本的调度敏感性分析把P2G效率和碳捕集率作为变量观察系统在不同技术成熟度下的最优运行成本曲线。如果曲线显示即便技术参数很理想经济性依然没有优势那说明这个项目在当前的园区负荷结构和碳价水平下不具备落地条件趁早放弃比勉强上马更明智。这一条我觉得比任何模型公式都实用。