企业官网建设流程全解析

燃烧后CO2捕集MEA吸收过程模拟教程（01）：燃烧后CO₂捕集为什么难模拟——MEA吸收体系的物理与化学全景

版本声明块

• 工具/软件：Aspen Plus（官方产品页未公示数字版本号，本系列一律以官方文档为准）；Python 3.10 及以上
• 语言/环境：Python，仅使用标准库dataclasses，无需安装任何第三方包
• 本文目标：建立 MEA 吸收体系的化学、流程与工况量级认知，并说明平衡级模型的能力边界

一句话结论：MEA 吸收塔之所以难模拟，是因为 CO₂ 与单乙醇胺的反应总级数为二级（对 CO₂ 一级、对胺一级，Penny 1983），并同时走氨基甲酸盐与 HCO₃⁻ 两条路径，因此必须先用 30 wt% MEA、富液 0.275–0.476 mol CO₂/mol MEA 这类带工况的负载量级约束模型，再决定用 RadFrac 平衡级还是速率基塔。

〇、本篇要解决的认知问题

1. 燃烧后 CO₂ 捕集（PCC）到底让电厂付出多大代价，为什么"溶剂再生"是主要矛盾？
2. MEA 与 CO₂ 之间究竟发生了什么化学反应，所谓"两条路径"指什么？
3. 一套 MEA 吸收—再生流程的最低限度骨架由哪些设备构成？
4. 30 wt% MEA、捕集率 90%、富液 0.275–0.476 mol CO₂/mol MEA 这些数字该怎么用，绝不能怎么用？
5. 为什么"平衡级"假设在 MEA 吸收塔上会明显失真，必须走向速率基模型？

一、机制解析

1.1 先看代价：PCC 为什么吃掉十来个效率百分点

底账 [C5]（经其转引 Kothandaraman et al. 2009）给出两条可直接引用的结论：其一，溶剂再生约占整个 PCC 能耗惩罚的 60%——注意这是二次引用，不得当作普适常数；其二，PCC 使电厂净效率下降约 10–12 个百分点，主因是溶剂再生，其次是 CO₂ 压缩与辅机。

"再生贵"的物理原因很直白：再沸器要同时付出三笔热量——升温显热、解吸反应热、水汽化潜热（stripping steam）[C5]。这三项折算成工程 KPI 就是SRD（Specific Reboiler Duty，比再沸器负荷），单位 GJ/t CO₂。底账已核实、可直接引用的基线如下（口径为 MJ/kg CO₂ 与 GJ/t 等价，4.10 MJ/kg = 4.10 GJ/t）：

出处：TCM 系列见 [S28][S29]；Sherman 与 CCSI 见 [C14]；江油见 [C20]；综述区间见 [C5]。

这里必须建立第一条写作纪律：不许把 3.5–4.0 GJ/t 拍平成一个"通用值"。同一台装置，捕集率从 83.5% 提到 90%，TCM 的 SRD 从 3.90 变到 4.10；捕集率再提到 99%，CCSI 模型给到 3.77；而燃气 CCGT 机组的整体水平显著高于煤电。能耗口径与拆分方法会在第 10 篇展开，本篇只建立量级与"数值必须带工况"的引用习惯。

1.2 化学：MEA 为什么"难算"

MEA（monoethanolamine，单乙醇胺）是伯胺。底账 [C5] 明确指出 CO₂ 在 MEA 水溶液中走两条并存的亲核路径：

1. 氨基甲酸盐路径（主）：MEA 与 CO₂ 结合生成氨基甲酸盐（carbamate），在溶液中主要以 MEACOO⁻ 形式存在；
2. 水合路径：CO₂ 与 OH⁻ 结合生成 HCO₃⁻（碳酸氢根）。

两条路径的平衡常数不同、对温度的敏感度不同，而溶液里同时存在 MEA、MEAH⁺、MEACOO⁻、HCO₃⁻、CO₃²⁻、H⁺、OH⁻、H₂O、CO₂ 等多个物种——这就是"MFA 体系属于电解质体系"的真实含义。这些物种的活度系数怎么算，正是第 03 篇三类热力学模型（eNRTL / Kent-Eisenberg / Deshmukh-Mather）的分水岭所在。底账给出的平衡反应式编号（式 2–6）请回 [C5] 原文核对。

反应级数方面，有一条可以放心引用的硬事实来自 Penny (1983)：CO₂ 与伯胺在水溶液中的反应总级数为二级，对 CO₂ 一级、对胺一级[C15]。经典实验文献是 Pinsent, Pearson & Roughton (1956) 与 Hikita, Asai 等 (1977) [C16][C17]。二级速率常数普遍写成 Arrhenius 形式k₂ = A·exp(−Ea/RT)，但A 与 Ea 的具体数值本系列一律不写——底账 §7-② 明确其未获一手核对，需要时请回原文补。可用的旁证只有一条：LLNL/CCSI 报告引用 Hikita 等 (1977) 回归得到的反应速率系数平均偏差 3.5%[C13]，说明"速率参数本身带不确定度"应当作为建模前提（第 07 篇专门做不确定性传播）。

1.3 流程骨架

贫液 (lean, 进塔温度实测 49 ℃ [C14]) 净化气 (vent) | ^ v | +---------------------+ +------------------+ | 吸收塔 Absorber |<--- 烟气 (DCC 冷却至 20–25 ℃ 并增湿 [C19]) | 规整填料 15–24 m | +------------------+ | Sherman 15 m / TCM 24 m (12+6+6 m) +---------------------+ ^ | 富液 (rich) | 贫液回流 (经 trim cooler) v | +---------------------+ +----------------------+---------+ | 贫富液换热器 Lean- |---->| 再生塔 / 汽提塔 Stripper | | Rich Cross Exchanger|<----| 再沸器 SRD 3.62–4.10 GJ/t CO₂ | +---------------------+ | 再生压力 1.5–2.4 bar | +----------------------+-----------+ | ^ v | CO₂ 产品气 (压缩) 再沸器蒸汽 (低压蒸汽)

这套拓扑的信息量集中在四处：吸收塔进液温度、富/贫液负载差、再生压力、填料高度。贫富液换热器（lean-rich cross exchanger）负责用再生塔出来的热贫液预热冷富液，必要时后面再接 trim cooler 把进吸收塔的贫液温度压回目标值 [C19]；这是"吸收塔性能"与"再生塔能耗"之间最直接的耦合通道，第 09 篇会专门拆开。

1.4 关键工况量级（全部带出处与条件）

关于表中"吸收塔温度 / 再生塔温度"这一行：吸收塔 40–60 ℃、再生塔 100–120 ℃ 属于行业常见区间，但底账明确未锁定单一权威出处，因此正文只写"典型设计区间（经验值）"，不得声称其出自 ISO 或厂商标准。真正可锚定的是：贫液进塔 49 ℃、DCC 冷却至 20–25 ℃、再生压力 1.5–2.4 bar。

还须强调负载的单位基：一律写成mol CO₂/mol MEA，且富/贫液分列。把用户口径的"负载 0.4–0.5"直接当全区间使用是错的——只有富液上限成立（约 0.45–0.48），贫液通常落在 0.2–0.3；实测区间分别是富液 0.275–0.476、贫液 0.075–0.239 [C14]。凡是省略基、把富贫混合成一句的写法，本系列一律不接受。

1.5 为什么平衡级不够

平衡级（理论级）模型的核心假设是"离开每一级的汽液两相达到相平衡"，塔高被抽象成"理论级数"。但 MEA 吸收 CO₂ 是快速反应：反应主要发生在液膜内，而不是在塔板上等待整体达到平衡。速率基模型用速率方程与传质系数计算界面上的质量与能量传递，替代"理论板"假设 [C10]；代价是它显式要求黏度、表面张力、导热系数、二元扩散系数等传递物性模型 [C1]，并且求解变量远多于平衡级。

工程上的正确路径被写进了 AspenTech 官方课程：先把 RadFrac 平衡级 CO₂ 吸收塔收敛，再转换成使用速率基模型的塔[S15]。这条"先平衡级、后速率基"的路径是本系列第 08 篇的主线（第 05 篇讲原理，第 15 篇讲不收敛时的诊断）。本篇后面第三段代码会用最朴素的平衡级模型算出"理论级数"，让你亲眼看到：这个数既不含反应速率，也不含有效传质面积，因此不能直接换算成填料高度。

二、完整代码与逐行剖析

三段代码全部无第三方依赖，可直接复制运行。所有物性与热效应均为示意参数，用于建立数量级直觉；真实的定量结果必须交给 Aspen Plus 这类带电解质物性包的过程模拟器。

代码 1：PCC 物流量与 CO₂/溶剂摩尔衡算 + 流程骨架打印

# 01_code1_pcc_mass_balance.py# 运行环境：Python 3.10+（仅标准库，无第三方依赖）# 用途：PCC-MEA 吸收体系的量级衡算（CO2 / 溶剂循环 / L/G）+ 流程骨架 ASCII 打印# 【重要】下列物性与工况均为“示意参数”，仅用于建立数量级直觉，# 不能替代 Aspen Plus 等过程模拟器的电解质物性计算。fromdataclassesimportdataclass# ---------------- 示意参数区（非任何装置的设计值） ----------------M_MEA=61.08# 单乙醇胺摩尔质量, g/molM_FG_AVG=28.5# 烟气平均摩尔质量, g/mol（示意）MEA_WT=0.30# MEA 质量分数：底账基准为 30 wt%（[C19][C14][C5]）CAPTURE=0.90# 捕集率：底账经典设计基准 90%（[C19][C14]）RICH=0.45# 富液负载, mol CO2/mol MEA（落在底账实测 0.275-0.476 内, [C14]）LEAN=0.23# 贫液负载, mol CO2/mol MEA（落在底账实测 0.075-0.239 内, [C14]）@dataclassclassFlueGas:# 烟气总量与 CO2 含量都是示意值；底账核实的烟气 CO2 量级为 3.5-15%（[C14][C19]）total_kmol_per_h:float=10_000.0y_co2:float=0.035defpcc_mass_balance(fg:FlueGas)->dict:# 先把 CO2 摩尔流固定下来：它是全流程的被计量物，负载、循环量、SRD 都以它为基准，# 先钉住分母，后面才不会出现“按烟气算”和“按 CO2 算”两套口径串味。n_co2_in=fg.total_kmol_per_h*fg.y_co2 n_co2_cap=n_co2_in*CAPTURE# 由“目标负载差”反推溶剂循环量，而不是先拍一个循环量再回算负载：# 吸收塔的贫液负载 / L/G / 填料高度三者强耦合（第 12 篇展开），# 先定负载差再反推循环量，才能保证输出结果是一致性可检验的。d_loading=RICH-LEAN n_mea=n_co2_cap/d_loading# kmol/hm_mea=n_mea*M_MEA/1000.0# t/hm_sol=m_mea/MEA_WT# t/h（按 30 wt% 折算成溶液）m_fg=fg.total_kmol_per_h*M_FG_AVG/1000.0return{"n_co2_in_kmol_h":n_co2_in,"n_co2_cap_kmol_h":n_co2_cap,"d_loading":d_loading,"n_mea_kmol_h":n_mea,"m_sol_t_h":m_sol,"m_fg_t_h":m_fg,"lg_kg_per_kg":m_sol/m_fg,}PROCESS_ASCII=r""" 贫液 (lean, 进塔温度实测 49 ℃ [C14]) 净化气 (vent) | ^ v | +---------------------+ +------------------+ | 吸收塔 Absorber |<--- 烟气 (DCC 冷却至 20-25 ℃ 并增湿 [C19]) | 规整填料 15-24 m | +------------------+ | Sherman 15 m / TCM 24 m (12+6+6 m) +---------------------+ ^ | 富液 (rich) | 贫液回流 (经 trim cooler) v | +---------------------+ +----------------------+---------+ | 贫富液换热器 Lean- |---->| 再生塔 / 汽提塔 Stripper | | Rich Cross Exchanger|<----| 再沸器 SRD 3.62-4.10 GJ/t CO2 | +---------------------+ | 再生压力 1.5-2.4 bar | +----------------------+-----------+ | ^ v | CO2 产品气 (压缩) 再沸器蒸汽 (低压蒸汽) """defmain()->None:fg=FlueGas()r=pcc_mass_balance(fg)print("===== PCC-MEA 吸收体系量级衡算（示意参数） =====")print(f"烟气总摩尔流率 :{fg.total_kmol_per_h:,.0f}kmol/h")print(f"烟气 CO2 摩尔流率 :{r['n_co2_in_kmol_h']:.1f}kmol/h")print(f"被吸收 CO2（捕集率{CAPTURE:.0%}）:{r['n_co2_cap_kmol_h']:.1f}kmol/h")print(f"富/贫液负载（示意） :{RICH:.3f}/{LEAN:.3f}mol CO2/mol MEA")print(f"负载差 :{r['d_loading']:.3f}mol CO2/mol MEA")print(f"所需 MEA 摩尔流率 :{r['n_mea_kmol_h']:.1f}kmol/h")print(f"所需溶液质量流率 :{r['m_sol_t_h']:.1f}t/h（30 wt% MEA）")print(f"烟气质量流率 :{r['m_fg_t_h']:.1f}t/h")print(f"L/G（液气比） :{r['lg_kg_per_kg']:.3f}kg/kg")print("参照：Sherman FEED L/G = 1.069 kg/kg [C14]")print(PROCESS_ASCII)if__name__=="__main__":main()

预期输出与判读：L/G会落在 1.02 kg/kg 左右，与 Sherman FEED 的 1.069 kg/kg [C14] 相差约 4%，说明本脚本的"烟气量—CO₂ 量—负载差—循环量"链条在量级上自洽。但这不代表对上了中试：示意参数是人为设定的，Sherman 的真实烟气组成、MEA 浓度（35 wt%）与贫富液负载都与这里不同，任何偏差都必须回到工况层面解释，而不是调参数去"凑"一个漂亮数字。

代码 2：再沸器负荷三项拆分的量级演示

# 01_code2_srd_split.py# 运行环境：Python 3.10+（仅标准库）# 用途：把 SRD 按三项热量拆分，并与底账已核实的中试基线并排对照# 【重要】三项热效应的绝对量是“示意参数”，真实占比须用实际物性包与工况计算；# 对照用的 SRD 数值才是底账已核实的引用值，不能反过来把示意值当结论。# 示意参数：解吸反应热 / 水汽化潜热 / 升温显热，单位 GJ/t CO2Q_DESORB=1.90# 示意Q_WATER=1.10# 示意Q_SENS=0.55# 示意# 底账已核实的 SRD 基线：每条都带工况条件，禁止拍平BASELINES=[("TCM MEA 基线 2014",4.10,"30 wt% MEA, 3.5-3.6% CO2 烟气, 捕集 90.0%"),("TCM MEA 2015",3.90,"捕集 83.5%"),("TCM MEA + 消泡剂 2015",3.62,"捕集 87.4%（现行基线）"),("Sherman FEED (Bechtel/DOE)",3.65,"35 wt% MEA, 捕集约 95%"),("CCSI / Aspen 模型",3.50,"35 wt% MEA, 捕集 95%"),("CCSI / Aspen 模型",3.77,"35 wt% MEA, 捕集 99%"),]defsplit(total:float)->dict:# 三项热效率的比例是示意值：把它当作“固定的拆分模型”，# 再按比例映射到已核实的 SRD 基线上，只是为了让读者看清口径差异，# 绝不是声称中试装置的真实占比就等于这个比例。base=Q_DESORB+Q_WATER+Q_SENSreturn{"ratio":(Q_DESORB/base,Q_WATER/base,Q_SENS/base),"rescale":(total*Q_DESORB/base,total*Q_WATER/base,total*Q_SENS/base),}defmain()->None:base=Q_DESORB+Q_WATER+Q_SENS rd,rw,rs=split(base)["ratio"]print("===== 再沸器负荷三项拆分（示意参数） =====")print(f"示意合计 ={base:.2f}GJ/t CO2 (解吸{rd:.1%}/ 汽化{rw:.1%}/ 显热{rs:.1%})")print("\n----- 底账已核实 SRD 基线（按同一示意拆分比例映射） -----")forname,srd,condinBASELINES:qd,qw,qs=split(srd)["rescale"]print(f"{name:<28s}SRD={srd:.2f}GJ/t | "f"解吸{qd:.2f}+ 汽化{qw:.2f}+ 显热{qs:.2f}|{cond}")print("\n注意：解吸热与汽化潜热的占比随贫液负载与再生压力显著变化，")print(" 第 10 篇将用带出处的区间做完整拆分，本篇仅演示拆分方法。")if__name__=="__main__":main()

反直觉点：解吸反应热在直觉上"最贵"，但水汽化潜热往往占到相当份额，而降汽化潜热的直接手段是提高贫液负载、缩小循环量——这又会让吸收塔更难达标。这正是"负载差—循环量—再生能耗"三角耦合的根源，也是第 10 篇的主线。

代码 3：平衡级理论级数——看清理想化假设的代价

# 01_code3_equilibrium_stages.py# 运行环境：Python 3.10+（仅标准库）# 用途：用最朴素的平衡级模型（线性平衡关系 + Kremser 公式）反推达到目标# 捕集率所需的理论级数，用来说明平衡级假设遗漏了什么。# 【重要】线性平衡关系 m 是示意参数；真实的 MEA-H2O-CO2 体系必须用# 电解质活度模型（第 03 篇）与速率基传质（第 05 篇）才能定量。M_EQ=0.60# 示意表观分配系数 y = m * x（气相摩尔分数 / 液相摩尔分数）L_OVER_V=1.4# 示意摩尔液气比（真实值须由溶剂循环量与烟气量算出）defstages_needed(A:float,e_target:float,n_max:int=200)->int:"""Kremser 吸收公式 E = (A^(N+1) - A) / (A^(N+1) - 1)，数值扫描求最小 N。"""forninrange(1,n_max):den=A**(n+1)-1.0ifden==0:continue# A=1 是公式的奇点，教学模型直接跳过if(A**(n+1)-A)/den>=e_target:returnnreturn-1# -1 表示在 n_max 级内无法达标defmain()->None:a=L_OVER_V/M_EQ# 吸收因子 A = L / (m * V)print(f"示意吸收因子 A = L/(m*V) ={a:.3f}（L/V={L_OVER_V}, m={M_EQ}）")fortargetin(0.80,0.90,0.95,0.99):print(f"目标捕集率{target:.0%}-> 理论级数 N ={stages_needed(a,target)}")print("\n判读：理论级数只描述“需要多少次平衡接触”，")print("它既不含反应速率，也不含填料的有效传质面积，")print("因此不能直接换算成填料高度——这就是必须走向速率基的原因。")if__name__=="__main__":main()

预期输出：目标捕集率 80%、90%、95%、99% 对应的理论级数约为 2、3、3、5。请注意这个模型里"传质"是隐形的：它假设每一级都能瞬间达到平衡，于是塔只要"够多级"就行。真实填料塔里，快速反应发生在液膜内、有效传质面积随液体负荷变化、黏度随负载升高而上升，这些全部被平衡级假设抹掉了 [C10][C1]。

三、常见报错与排查

本节按"现象 → 根因 → 解法"组织，覆盖四个概念性误用与一个脚本级问题。

现象 1：把捕集率 90% 与 SRD 直接绑定，写成"MEA 捕集能耗 = 3.62 GJ/t"。
根因：SRD 与捕集率、溶剂浓度、机组类型强耦合，TCM 现行基线 3.62 GJ/t 对应捕集 87.4%、30 wt% MEA [C19]。解法：引用 SRD 时必须同时写出捕集率与溶剂浓度；燃气 CCGT 要另取 4–5 GJ/t 区间 [C5]。

现象 2：文中出现"MEA 负载 0.4–0.5"这样的无基写法，且富贫液混为一谈。
根因：把用户口径的富液上限当成了全工况区间。解法：富液 0.275–0.476、贫液 0.075–0.239 mol CO₂/mol MEA 分列并带工况（35 wt% MEA、CCGT 烟气）[C14]。

现象 3：把 eNRTL 的出处写成 Chen & Song (2004)。
根因：这是流传较广的错误引用。解法：eNRTL 的规范引用是Chen & Evans (1986)（AIChE J. 32, 444–454）及 1982 年单电解质原始文献 [C3][C4][C5]；Chen & Song (2004) 是 NRTL-SAC 溶解度模型，两者不是同一模型 [C8]。第 03 篇会专门纠偏。

现象 4：用平衡级 RadFrac 直接预测出口净化气 CO₂ 浓度并当作设计依据。
根因：快速反应发生在液膜内，平衡级假设不成立。解法：把平衡级解当作初值，再切换到速率基塔（AspenTech 官方课程给出的教学路径 [S15]）；原理见第 05 篇，实操见第 08 篇。

现象 5：stages_needed返回 -1，误以为程序出错。
根因：在n_max=200级内无法达到目标捕集率，通常是吸收因子 A 太小（液气比不足或平衡关系太不利）。解法：先确认L_OVER_V / M_EQ是否明显小于 1；A < 1 时单级吸收率极低，理论上需要极多级，这不是 bug 而是物理结论。

四、动手练习

练习 1（校验量级）：把代码 1 的CAPTURE依次设为 0.90 与 0.99，记录L/G的变化。
成功判定：0.99 时L/G应比 0.90 时高出约 10%（因为要吸收的 CO₂ 增多、负载差不变）。若你的结果不满足该比例，请检查是否误把d_loading写成了常量以外的量。

练习 2（校验单位基）：在代码 1 中把RICH改成 0.48、LEAN改成 0.09，重算L/G。
成功判定：负载差从 0.22 增大到 0.39，因此所需溶剂循环量与L/G应下降到约 0.58 kg/kg（低于原先的 1.023 kg/kg）。记录该值并在笔记中标注"负载单位为 mol CO₂/mol MEA，且富/贫分列"；若你的结果反而变大，说明把负载差用反了。

练习 3（理解平衡级边界）：在代码 3 中把M_EQ从 0.60 改为 1.20，观察 90% 与 99% 对应的理论级数。
成功判定：理论级数应明显增加（平衡关系变差、吸收因子 A 变小）。请写下结论：“平衡关系越不利，平衡级模型要求越多理论级；但真实塔高还受传质系数控制，二者不是一回事。”

五、小结与下一篇预告

本篇建立了四件事：PCC 的代价结构（净效率下降约 10–12 个百分点、溶剂再生是主要矛盾，后者为二次引用）、MEA 吸收的双路径化学（氨基甲酸盐 + HCO₃⁻）、流程骨架与带工况的工况量级、以及平衡级模型的能力边界。三段代码分别演示了物流衡算、SRD 三项拆分与理论级数计算，并全部标注为示意参数。

下一篇（02）转向工具层：把 Aspen Plus、Aspen HYSYS、DWSIM、COCO/COFE、ProSimPlus、gPROMS、CHEMCAD 的能力边界逐项摆开，说明"速率基 + 电解质双支持"目前最成熟的是哪一类工具，开源路线缺什么部件，并给出一个探测本机 COM 环境自检脚本。再往后，第 03 篇进入热力学模型的第一道分水岭，第 04 篇就动手搭第一个能跑起来的平衡级吸收塔。

本篇认知问题回显（FAQ）

Q1：燃烧后 CO₂ 捕集（PCC）让电厂净效率下降多少个百分点？

A：PCC 使电厂净效率下降约 10–12 个百分点，主因是溶剂再生，其次是 CO₂ 压缩与辅机。

Q2：MEA 吸收 CO₂ 有哪两条反应路径，生成什么物种？

A：主路径是 MEA 与 CO₂ 生成氨基甲酸盐（溶液中为 MEACOO⁻），并存水合路径 CO₂ 与 OH⁻ 生成 HCO₃⁻。

Q3：MEA 捕集流程由哪些主要设备组成？

A：主要由吸收塔、贫富液换热器、再生塔（汽提塔）、再沸器和 CO₂ 压缩组成，贫液经 trim cooler 回流吸收塔。

Q4：MEA 吸收体系的富液与贫液负载实测区间是多少？

A：35 wt% MEA、CCGT 烟气条件下，富液负载实测 0.275–0.476、贫液负载实测 0.075–0.239 mol CO₂/mol MEA，两者必须分开列出并标注单位基。

Q5：MEA 吸收塔为什么不能用平衡级模型直接预测？

A：MEA 吸收 CO₂ 是快速反应，主要发生在液膜内，平衡级的理论板假设不成立，需要速率基模型计算界面传质。