做燃烧后碳捕集模拟这些年,我大部分时间都耗在“怎么让流程收敛”和“怎么让能耗下来”这两件事上。最近完成的一套百万吨级CO2捕集系统模拟,用的是ASPEN Plus V11,溶剂走的是复配胺路线,重点做了两个工艺优化——吸收液中间冷却、吸收段结构调整。整个流程跑下来,再沸器能耗从单一MEA路线的3.9 GJ/tCO2压到了3.2 GJ/tCO2,捕集率还能稳定在90%以上。如果你正在做电厂、水泥窑或者钢厂的烟气碳捕集工艺设计,或者正纠结于要不要上中间冷却、复配胺到底怎么配,这篇内容应该能帮你省掉几轮反复试算。
我先把结论放前面:这套流程里最值钱的改动不是换溶剂本身,而是把吸收塔中部的温度剖面压了下来。吸收液中间冷却加上吸收段重新布置以后,溶剂循环量降低了将近两成,再生负荷也明显下降。复配胺则把吸收容量和反应速率同时照顾到了,避免了单一MEA那种“反应快但容量小,再生还费汽”的尴尬。下面按整个模拟项目的推进顺序,把设计基准、物性配置、中间冷却搭法、吸收段调整和那些容易卡住的坑完整记录一下。
1. 百万吨碳捕集:先算清规模,再谈工艺
做流程模拟最大的忌讳是:流程还没搭完,就在那里调塔板数。其实规模都不算清楚,后面所有用吸收液循环量、塔径、再沸器热负荷全都站不住脚。百万吨级听起来很唬人,但它落到模拟输入上,其实就是几个数。
1.1 烟气量和捕集速率怎么定
我先按年捕集100万吨CO2、年运行8000小时来反推。100万吨除以8000小时,平均捕集速率就是125 t/h。燃烧后烟气里的CO2浓度一般按13%体积分数考虑,捕集率设定为90%,那么进吸收塔的CO2总量就是125除以90%,约139 t/h。
换算成摩尔流量:139 t/h除以CO2分子量44,得到3160 kmol/h。如果烟气中CO2占13%,总烟气摩尔流量就是3160除以13%,约24300 kmol/h。标准状况下乘以22.4,结果是54.4万Nm3/h,也就是说烟气处理量大约是55万标方每小时。
这个数字落在百万吨级项目里非常典型。烟气质量流量大概是55万乘以1.29,再折算一下,约70万kg/h。如果按常见的液气比4.2 kg吸收液每kg烟气,吸收液循环量就在3000 t/h上下。这个数决定了吸收塔规格、泵的扬程、贫富液换热器面积,也决定了模拟里循环流股的规模,别拍脑袋给定。
提示:如果你做的是部分烟气捕集或燃烧后烟气再循环,CO2浓度会变,反推方式不变。先把捕集速率算出来,再设捕集率,所有工艺参数才有锚点。
塔径也可以粗算一遍。55万标方每小时除以3600,约151 Nm3/s。吸收塔操作温度40℃附近,真实体积流量乘上313.15除以273.15的温度修正,约174 m3/s。塔内气体操作速度取1.6 m/s左右,截面积约108 m2,塔径接近11.8 m。百万吨级吸收塔直径10~14 m很正常,要是你模拟结果里塔径只有3米,那说明循环量和烟气量大概率对不上。
1.2 优化为什么重点放在吸收侧
碳捕集全流程能耗大头在溶剂再生,也就是解吸塔再沸器。但再沸器热负荷的来源,很大程度上被吸收端的温度和负载决定。如果吸收塔里温度很高,富液CO2负载提不上去,解吸塔就要处理更大的循环量,再生每吨CO2的能耗自然就高。
吸收塔内的CO2吸收反应是强放热过程。烟气从塔底进、吸收液从塔顶往下走,塔中下部的CO2分压最高、反应最剧烈,液相温度经常从入口的40℃一路升到60℃甚至70℃。问题是温度越高,胺液表面的CO2平衡分压越高,吸收推动力掉得很快。吸收反应进行到后半段,差不多就是在“边吸边放”的拉锯里浪费塔高。
中间冷却是直接针对这个温升来的。它在吸收塔温度最高的位置把液体抽出来、冷却后再回去,把反应热带走。吸收段调整则是优化塔内气液接触级数和进料位置,让温度剖面和浓度剖面分布更合理。这两个改动都是从吸收侧发力,效果比单纯换溶剂更直接,因为反应热还在,但热量被挪走了。
2. 复配胺溶液:配比选择和模拟物性配置
单一胺溶液做碳捕集,优势是简单,劣势也明显。MEA反应快但吸收容量有限,再生能耗高;AMP容量大但反应速率偏慢;PZ抗降解性好但溶解度低,加多了容易在低温段析出。复配的思路就是取长补短,用主胺保证反应活性,用辅胺提高容量和抗降解能力。
2.1 最终配比怎么选的
我这个项目最终选了25 wt% MEA + 10 wt% AMP + 余量水。没有追求三组分复配,原因是多一个组分就多一组电解质反应和交互参数,模拟收敛难度和工程补加控制的复杂度都会上升。MEA负责快速反应和氨基甲酸盐生成,AMP利用位阻效应主要在吸收中后段贡献碳酸氢盐路径,提高整个溶液的CO2负载能力。
做对比模拟的时候,我扫过几套配方,结果趋势可以作为参考:
| 溶剂体系 | 吸收速率相对值 | 循环量相对值 | 再生能耗参考(GJ/tCO2) | 备注 |
|---|---|---|---|---|
| 30% MEA | 基准 | 1.00 | 3.9左右 | 稳定但要防降解 |
| 30% MEA + 5% PZ | 提升明显 | 0.85左右 | 3.4左右 | 注意低温析出风险 |
| 25% MEA + 10% AMP | 略低于MEA | 0.78左右 | 3.1~3.3 | 兼顾容量与再生温度 |
AMP的优势是它的位阻结构让CO2主要生成碳酸氢盐而不是单一氨基甲酸盐,理论上每摩尔胺可以吸收更多CO2。代价是反应速率稍慢,需要更高的填料段或更大传质面积。好在这个项目吸收塔本来就大,填料高度充足,这个代价可以接受。
2.2 ASPEN Plus V11里的组分、物性和化学反应
物性方法这块,胺-水-CO2体系我直接用ENRTL-RK,也就是Electrolyte NRTL配Redlich-Kwong状态方程。这个物性包对胺水溶液吸收CO2的电解质反应回归得很成熟,是燃烧后胺法模拟的事实标准。
组分列表里要放进N2、O2、CO2、H2O、MEA、AMP。不要漏掉烟气的惰性组分,尤其是O2,它对后续核算胺氧化降解和尾气氧含量有用。然后用电解质向导自动生成离子反应和分子组分,这一步会生成MEA正离子、AMP正离子、碳酸氢根、碳酸根、氢氧根、水合氢离子等。
MEA和CO2的反应路径在模拟里通常写成MEA + HCO3- ⇌ MEACOO- + H2O,或者直接用MEA + CO2生成氨基甲酸盐的形式,关键是要把可逆平衡反应和速率反应分别设置好。AMP因为位阻效应基本不生成稳定的氨基甲酸盐,所以只考虑它的质子化反应,不要强行加上生成阴离子那一步,否则会高估AMP的吸收容量。
亨利组分别忘了指定CO2、O2、N2。特别是CO2,如果不加亨利组分,气相CO2在液相里的物理溶解度会算成零,导致吸收塔出口气相浓度严重失真。这个细节是我早期做类似模拟踩过的坑,很多新手的模拟结果异常,问题就出在这。
2.3 物性和反应数据的校核思路
ASPEN里的内置数据足够做趋势分析和方案比选,但如果你手里有小试或中试数据,一定要用来校正。哪怕只是验证“给定CO2分压下富液平衡负载”这一个点,也比直接跑流程靠谱得多。
我当时没有专门做实验,所以我做了两件事:一是把模拟出来的贫液负载、富液负载和解吸塔塔底温度对照工业经验值,贫液负载如果低于0.15 mol/mol,说明再生条件过于苛刻,MEA热降解风险很大;二是看再生温度是否落在110~125℃区间,如果解吸塔压力太高导致温度冲到130℃以上,这套配方的优势就没了。
3. 吸收液中间冷却:在模拟软件里如何落地
中间冷却是吸收塔内部的一个液体循环。原理不复杂,但在ASPEN Plus里搭起来有几个关键点容易出错,尤其是抽液级和返回级的相对位置、冷却器压降设置、以及流股类型的指定。
3.1 中间冷却为什么有效
吸收反应放热,烟气从塔底往上走,吸收液从塔顶往下走,塔内温度最高点通常在中下部,也就是CO2分压最高、反应最剧烈的区域。温度一高,胺液表面的CO2平衡分压也跟着高,吸收推动力降低,相当于塔的中段在“空转”。
中间冷却把最高温度段的液体抽出来,冷却到40℃再送回去,塔内温度剖面整体压低。模拟结果里,无中间冷却时塔内液相峰值温度接近67℃;加了中间冷却后,峰值降到52℃左右,塔底富液出口温度也回落到46℃上下。温度低了,平衡推动力大了,同一循环量下富液负载能高出不少,解吸塔的负担就减轻了。
一个容易被忽略的细节:中间冷却并不是在“节省吸收反应热”,它本质上是把低品位热量直接排给冷却水。所以冷却负荷本身是增加的,但换来的再沸器能耗下降远大于冷却侧代价,节能的是系统总能耗。算综合能耗时,冷端电耗要折算,不能只看蒸汽。
3.2 在ASPEN Plus V11里搭中间冷却回路的步骤
吸收塔我用的是RadFrac模块。粗调阶段用平衡级模型,理论级设20级左右;精调阶段再切到Rate-based。中间冷却回路的具体搭法如下:
- 在吸收塔的Streams页面设置液相信号侧线采出,就是Liquid Draw,选在中间某个理论级上。比如20级塔,我选第12级,这个位置对应塔内温度峰值。
- 采出流股进一台换热器,用Cooler模块模拟。设定出口温度40℃、压降0.2 bar左右。注意流股如果发生闪蒸出现气相,说明冷却后温度低于泡点或者压力下降太多,建议把压降控制在0.3 bar以内,并检查液相分率。
- 冷却后的液体压力可能不够回到塔内,需要在冷却器后加一个泵模块,升压约0.5~1 bar,抵消位差和管道阻力。
- 返回时不要回到抽液的那一级,而是回到下面第13级。这样更贴近工程实际,也避免同一级液相流量突然翻倍造成水力学过载。
- 返回流股类型选择Liquid,进料位置填塔编号和理论级编号,不要顺手填成Vapor。
如果想让中间冷却的效果更准确,冷却器出口富液的压力和吸收塔该级压力要能衔接。我给泵定了一个固定出口压力,让它等于吸收塔返回级的压力再加0.3 bar富余。
3.3 关键参数的敏感性分析
有两个变量对中间冷却效果影响最大:冷却器出口温度和抽出液流量占循环量的比例。我分别做了单变量扫描,结果比较直观。
冷却出口温度的影线趋势是这样:
| 冷却出口温度 | 捕集率(%) | 再沸器能耗(GJ/tCO2) | 中间冷却负荷(MW) |
|---|---|---|---|
| 45℃ | 88.5 | 3.40 | 6.8 |
| 40℃ | 90.2 | 3.20 | 9.5 |
| 35℃ | 90.8 | 3.10 | 13.6 |
| 30℃ | 91.0 | 3.05 | 18.2 |
温度越低,捕集率越高、再沸器能耗越低,但冷却负荷涨得很快。需要强调的是,40℃对应的冷却负荷9.5 MW,如果用循环水就能实现;35℃以下往往要上冷冻水或者凉水塔极限工况,工程投资和维护成本完全不一样。所以这个项目的设计点取在40℃,没有盲目追低温。
抽出流量的优化趋势也很有意思:刚开始增加抽出流量时,中间冷却效果提升明显;当抽出量达到循环量的70%左右后,再加大抽出量,捕集率几乎不动,冷却负荷反而继续线性上涨。最后取了个折中值,70%循环量作为设计点。
提示:不要只做温度灵敏度就完事。把抽出液流量和冷却温度放在一起做二维扫描,才会找到“冷却能力-再生能耗”的最佳平衡点。
4. 调整吸收段与全流程能耗平衡
吸收段调整这块,核心目标是在不大幅增加塔压降和塔高的前提下,让温度剖面和浓度剖面更有利。我在模拟里动了三处:理论级数、贫液进塔位置、水洗段布置。
4.1 吸收段结构调整的几个方向
理论级数是最直接的变量。我用平衡级先扫,16级、20级、24级、28级四组。16级时捕集率明显不够,20级够用,24级以后捕集率提升幅度变得很平缓。这说明塔的有效高度在这套体系和操作条件下已经基本到极限,再加级数只会增加压降和投资。
贫液进塔位置也值得调。有些流程为了让上段更“干净”,会把吸收液分成两股,大部分从塔顶进入,小部分从顶部往下第3级左右进入。这样上段保持高度贫液状态,可以把尾气里的CO2进一步往下压。模拟里试了两股进料后,塔顶出口CO2浓度从2.3%降到1.8%左右,不过贫液泵和分布器要稍微复杂一点。
水洗段必须加。吸收塔顶部预留2~3级,用冷凝水或脱盐水循环,把尾气夹带的胺蒸汽洗下来。在ASPEN里,水洗段就是吸收塔顶部几个理论级,加一股冷却水循环回塔。别看这事小,如果不加,MEA蒸汽损失会被模拟放大,溶剂补充量直接超标。
4.2 L/G比和贫液负载的权衡
吸收液循环量和贫液负载之间的关系,是整个工艺优化的灵魂。贫液负载低,吸收推动力大,循环量可以小;但贫液负载是靠解吸塔更彻底的再生换来的,再沸器热耗必然增加。贫液负载高,再生能耗低,但吸收段需要更大的液气比才能维持捕集率,泵和塔径都跟着涨。
我用手动调整配合设计规定做了几组工况:
| 贫液负载(mol/mol) | 液气比(kg/kg) | 再沸器能耗(GJ/tCO2) | 趋势判断 |
|---|---|---|---|
| 0.30 | 4.8 | 3.50 | 吸收推动力不足 |
| 0.25 | 4.2 | 3.20 | 综合较优 |
| 0.20 | 3.7 | 3.45 | 再生太狠,得不偿失 |
做了这个表以后,结论很清楚:这个体系的最优贫液负载在0.25附近。低于0.20,再沸器温度和蒸汽消耗都上去了,而且塔底温度超过125℃会显著加速胺液热降解,工程上不是不能做,是运行成本压不住。
4.3 全流程优化前后的数据对比
整轮调整完成后,我会习惯性把基准工况和优化工况并排放一次,确保每个改动都有据可查。这次项目最终对比数据如下:
| 对比项 | 基准工况 | 优化工况 |
|---|---|---|
| 溶剂配方 | 30% MEA | 25% MEA + 10% AMP |
| 吸收塔理论级 | 20级 | 24级(含中部抽出) |
| 中间冷却 | 无 | 40℃、70%循环量 |
| 液气比 | 5.5 kg/kg | 4.2 kg/kg |
| 捕集率 | 87% | 90.2% |
| CO2产品纯度 | 99.2% | 99.5% |
| 再沸器能耗 | 3.9 GJ/tCO2 | 3.2 GJ/tCO2 |
再沸器能耗从3.9压到3.2,降幅接近18%,捕集率不降反升。对一个百万吨级项目来说,这个能耗差距意味着每年蒸汽消耗差出十几万吨,换算成运行成本是非常可观的数字。
这里也补充一下解吸塔的条件,免得大家以为只有吸收塔在改。解吸塔压力我维持在4.5 bar左右,塔底温度约120℃,回流比经过优化控制在1.3附近。解吸塔本身没有做太激进的改动,因为吸收侧的改善已经让贫液负载更容易达标,解吸塔操作变得更从容了。
5. 实操中的问题排查与避坑记录
这部分我打算写成问题速查的形式,不说空话,只列我这次实际遇到过的情况和对应的处理办法。
5.1 收敛困难大部分出在循环回路
这种碳捕集流程本质上是个大循环:吸收塔底富液去解吸塔,解吸塔底贫液再回吸收塔。再加上中间冷却的内部小循环,流程撕裂流股少说有三四处。直接全流程一起算,十有八九是不收敛的。
我的做法是分两步走。第一步,先把吸收塔当成独立的吸收单元跑,贫液流量和贫液负载按输入值固定,解吸塔暂时断开,确认吸收塔单独收敛。第二步再接上解吸塔,但解吸塔再沸器热负荷先给一个固定值,等流股迭代稳定后再改成由塔底温度或设计规定控制。这个方法能规避大部分初值敏感问题。
电解质反应也是收敛杀手。第一次跑通时,我先把化学反应全关掉,物理吸收收敛了再逐条打开平衡反应,最后再开速率反应。如果你一上来就开着全部反应,ASPEN Plus会疯狂迭代,出现负流量或者物性计算失败都不奇怪。
5.2 中间冷却回路的水力学问题
中间冷却返回位置如果设置不当,会出现液相流量在某一级翻倍的情况。平衡级模型可能不报错,但切到Rate-based以后,RadFrac的水力学检查会提示该级液泛或者压降异常。解决办法就是前面说的:返回级放在抽出级下方1~2级,而不是同一级。
还有一个隐藏问题是冷却器压降设得太大。液体的饱和压力对温度很敏感,如果压降超过0.3 bar,冷却后的液体可能在管道里闪蒸,ASPEN返回流股就是气液混合相,这时候你一定要在塔的进料设置里分别指定气相和液相进料,不能图省事只填一股混合流。否则塔级能量平衡错得离谱,温度剖面全乱。
5.3 溶剂损失和烟气预处理的模拟边界
模拟里看不见的损失不代表工程上不存在。ASTM那套纯模拟默认烟气已经脱硫除尘,但如果你工艺包上游烟气SO2浓度没控制好,热稳定盐会在系统里越积越多,胺液再生效果直线下降。做物料平衡时,我建议单独设一股胺补充流股,把蒸发和氧化损失折算成补充量,否则整个溶剂平衡是假的。
吸收塔顶水洗段有没有,直接决定尾气里MEA含量。我在最后一次流程里给吸收塔顶部加了两级水洗,冷却水温度设25℃,模拟出的胺损失小了两个数量级。真实项目里这一步绝对不能省,并且水洗段循环液要定期排出一部分回吸收系统,避免杂质累积。
5.4 模拟精度和调试时间之间的取舍
做这类大系统,没有必要一开始就用Rate-based模型死磕。我建议先用平衡级模型把全流程的热负荷、循环量、理论级数和关键温度跑稳定,再换成Rate-based模型验证塔径和填料高度。Rate-based模型对填料参数很敏感,金属鲍尔环、规整填料、散堆填料的HETP和压降数据都要仔细填,填错了结果会偏得离谱。
我自己通常把流程模拟分成三个阶段:物性验证、平衡级调参、速率级校核。前两个阶段花的时间最多,因为需要理解系统趋势;Rate-based阶段反而只是确认,并不会改变前面定下来的最优操作条件。你要是时间紧,直接上Rate-based也不是不行,但调试时间至少要翻一倍。
最后再说个小技巧:每次改完胺配比、吸收塔段数或者中间冷却参数,先重新初始化一遍全局变量,再从已经收敛的工况继续迭代。不要贪快,拿上次结果硬算,电解质体系牵一发动全身,流程在边界附近的韧性比你想的差得多。做碳捕集模拟没有捷径,把每个变量的敏感性摸清,比任何高级收敛技巧都更管用。