燃烧后CO2捕集MEA吸收过程模拟教程03热力学模型的第一道分水岭——eNRTL、Kent-Eisenberg 与 Deshmukh-Mather版本声明块工具/软件Aspen Plus官方产品页未公示数字版本号以官方文档为准Python 3.10 及以上语言/环境Python仅使用标准库math、re、itertools无需第三方包本文目标掌握三类电解质热力学模型的假设、适用区间与选型判据并纠正一处流传甚广的错误引用一句话结论MEA-H₂O-CO₂ 体系的热力学模型分三道eNRTL基于同号离子互斥与局部电中性假设规范引用为 Chen Evans (1986), AIChE J. 32, 444–454、Kent-Eisenberg半经验表观平衡常数高负载精度下降、Deshmukh-MatherDebye-Hückel 型加相互作用参数低 CO₂ 分压时比 Kent-Eisenberg 更准选错哪一道都会让吸收塔的 CO₂ 平衡分压算错。〇、本篇要解决的认知问题eNRTL 的物理假设是什么它的规范引用到底是哪一篇文献为什么eNRTL Chen Song 2004这个说法是错的正确的对应关系是什么Kent-Eisenberg 与 Deshmukh-Mather 分别在哪一段工况更准能不能互相替代Aspen 中平衡常数的温度关联式ln K(T) a b/T c·lnT d·T意味着什么风险给定一个 MEA 吸收模拟任务怎么用一条决策路径选出热力学模型一、机制解析1.1 为什么热力学模型是第一道分水岭MEA 吸收塔的推动力来自气相 CO₂ 分压与液相 CO₂ 平衡分压之差。气相好算难的是液相溶液里同时存在 MEA、MEAH⁺、MEACOO⁻、HCO₃⁻、CO₃²⁻、H⁺、OH⁻、H₂O、CO₂ 等物种见第 01 篇。给定负载与温度液相 CO₂ 的平衡分压要由所有物种的活度共同决定。这就是分水岭同一个实验点换一套活度模型算出的平衡分压可以差出一大截进而把塔高、液气比、SRD 全部带偏。1.2 eNRTL假设、来源与规范引用eNRTLelectrolyte NRTL的核心假设是同号离子互斥与局部电中性它先由 Chen、Britt、Boston、Evans (1982) 针对单溶剂完全解离体系提出随后由Chen Evans (1986)推广到多组分水溶液体系 [C3][C4][C5]。它的能力边界是明确的可计算所有物种离子与分子、溶质与溶剂的活度系数——这正是 MEA 这种离子物种 中性物种共存体系最需要的性质。必须纠正的错误引用流传中常把 eNRTL 说成 “Chen Song 2004”。事实是Chen Song (2004) 的论文标题为“Solubility Modeling with NonRandom Two-Liquid Segment Activity Coefficient Model”IECR 43, 8354–8362即NRTL-SAC 溶解度模型与 eNRTL 不是同一个模型 [C8]。本系列统一规范引用为Chen Evans (1986)可并列 1982 年单电解质原始文献。这条纠偏请写进你自己的文献笔记里因为一旦引错后面对参数回归第 13 篇的讨论会全线失焦。最新进展方面2024 年AIChE Journal有工作给出用于 MEA-H₂O-CO₂ 的modified eNRTL 含过量焓解析式[C9]可作为该方向仍在演进的旁证引用但它不是本系列的基础假设来源。从工程视角看eNRTL 的贵体现在两处。第一处是参数数量长程项需要离子尺寸与介电常数相关参数局部组成项需要一组相互作用参数这些参数不能凭空给出必须由 VLE、溶解度或量热数据回归得到。第二处是参数的可辨识性MEA-H₂O-CO₂ 体系的负载、温度、压力三个维度高度耦合如果实验数据覆盖的负载区间过窄回归出的参数一旦离开该区间就可能给出错误趋势。这正是第 13 篇要专门处理参数回归与残差诊断的原因。反过来看Kent-Eisenberg 之所以长期存在于工业现场正是因为它把这两笔成本压缩成一组表观常数用精度换取了可用性——这个取舍本身没有对错只有适用边界。1.3 Kent-Eisenberg半经验的代价Kent-Eisenberg及其 modified 版本源头是 Kent Eisenberg 1976 年Hydrocarbon Processing的 “Better data for amine treating”本系列经二次文献 [C6][C7] 转引属于半经验平衡模型用表观平衡常数去拟合实验负载计算快、参数少。代价也很清楚在高负载时精度相对下降[C6][C7]。它的优点是真实存在的——当你只需要快速扫大量工况、或者手里没有任何 VLE 数据做回归时K-E 类模型是能立刻跑起来的选择。但请不要把它当成严格热力学模型来写。1.4 Deshmukh-Mather长程项加相互作用参数Deshmukh-Mather源头是 Deshmukh Mather 1981 年Chemical Engineering Science的工作同样经二次文献 [C7] 转引走的是另一条路用Debye-Hückel 型长程项 相互作用参数描述电解质非理想性。文献给出的对照结论是在低 CO₂ 分压时Deshmukh-Mather 比 Kent-Eisenberg 更准[C7]而另一方面modified Kent-Eisenberg 在高负载下优于 Deshmukh-Mather[C6]。这两句话必须成对使用、并各自带工况低分压看 D-M高负载看 modified K-E。把两条结论各自截一半、拼成D-M 更准或K-E 更准的单一结论是本主题下最常见的写作错误。1.5 三者对照维度eNRTLKent-Eisenberg / modified K-EDeshmukh-Mather类型活度系数模型γ-φ严格热力学半经验平衡模型活度系数 / 电解质平衡模型核心构造同号离子互斥 局部电中性长程项 局部组成项表观平衡常数拟合负载Debye-Hückel 型 相互作用参数规范引用Chen Evans (1986)及 1982 原始文献Kent Eisenberg (1976)本系列经二次文献转引Deshmukh Mather (1981)本系列经二次文献转引优势工况需要全物种活度、需与速率基配套时快速扫描、参数少、高负载modified 版低 CO₂ 分压已知短板参数多、回归难无授权软件时难落地高负载精度下降高负载下不如 modified K-E软件支持Aspen Plus 内置 property method含 MEA 专用物性 insertK-E 类与 D-M 类可由用户/电解质模板实现 [C1][C2][S4]同上同上关于开源DWSIM 的自动化属性包注册表中没有 eNRTL 专用包有NRTLPropertyPackage、IdealElectrolytePropertyPackage、CAPEOPENPropertyPackage等[S16]因此禁止断言 DWSIM 支持 eNRTL若要在开源侧使用 eNRTL路径是自建或通过 CAPE-OPEN 接入外部属性包第 16 篇。1.6ln K(T)的形式与温度外推风险在 Aspen 中平衡常数随温度的变化采用ln K(T) a b/T c·lnT d·T形式 [C3]。这个形式本身很常见但它有明确的适用区间四个系数是在某个温度窗口内回归出来的把它外推到窗口之外例如从吸收温度外推到再生温度、或更高会得到没有物理意义的 K 值。因此吸收侧与再生侧的平衡常数参数最好分别在其各自温度窗口内确认涉及温度扫描的工况试验第 12 篇必须记录 K 参数的回归区间任何用一条 K(T) 曲线覆盖 40–120 ℃ 甚至更宽的写法都要标注为未经验证的外推。1.7 选型决策路径任务是否需要定量预测吸收/再生性能 ├─ 否只要量级、要快───────────────► Kent-Eisenberg半经验参数少 └─ 是 ├─ 是否有授权软件的内置电解质物性包 │ ├─ 是 ─► eNRTLChen Evans 1986若还要速率基塔必须配套传递物性 [C1] │ └─ 否 ─► 用 CAPE-OPEN 接入外部电解质包或自建第 16 篇 └─ 关注区间是哪一个 ├─ 低 CO₂ 分压 ─► Deshmukh-Mather 更准 [C7] └─ 高负载 ─► modified Kent-Eisenberg 优于 Deshmukh-Mather [C6]注意这条路径不是二选一之后互不相干速率基塔必须显式要求黏度、表面张力、导热系数、二元扩散系数等传递物性[C1]这意味着一旦你选择速率基第 05、08 篇eNRTL 类的完整物性包几乎是必选配置而不是可选项。二、完整代码与逐行剖析代码 1离子强度如何把活度系数传导到平衡分压# 03_code1_enrtl_activity_demo.py# 运行环境Python 3.10仅标准库 math# 用途演示 eNRTL 为什么必须由“长程项 局部组成项”两部分组成# 只用长程静电项活度系数会随离子强度单调下降高浓度区物理上不合理# 叠加局部组成项后出现“盐析”方向的反转并把 CO2 平衡分压整体抬高。# 【重要】本脚本用简化 Debye-Hückel 型长程项与线性局部项替代 eNRTL 的# Pitzer-Debye-Hückel 项与 NRTL 局部项物种分布、A_DH、B_ION、C_LOCAL# 全部是【示意参数】不能用于定量设计。真实计算须由 Aspen Plus 的# 电解质物性包含 eNRTL完成见第 04 篇。importmath# ---------------- 示意参数区 ----------------MEA_WT0.30# 底账基准 30 wt% MEA[C19][C14][C5]M_MEA61.08# MEA 摩尔质量, g/molM_H2O18.015# 水摩尔质量, g/molA_DH0.509# Debye-Hückel 常数25 ℃ 水溶液常见值示意B_ION1.5# 离子尺寸参数 B*a示意C_LOCAL0.35# 局部组成项示意系数对应 eNRTL 的 NRTL 相互作用参数LN10math.log(10.0)defmea_molality(wt:floatMEA_WT)-float:把 30 wt% MEA 换算成质量摩尔浓度 mol/kg 水示意换算。n_meawt/M_MEA# 每 1 g 溶液中 MEA 的摩尔数m_water_g1.0-wt# 每 1 g 溶液中水的质量, greturnn_mea/(m_water_g/1000.0)defionic_strength(loading:float,m_mea:float)-float:示意假设负载态 CO2 全部形成 MEACOO- 并伴随等量 MEAH 则 I 0.5*(m_cation m_anion) m_co2。 真实体系里 HCO3-/CO3^2-/OH- 都会参与必须由平衡求解给出。m_co2m_mea*loadingreturn0.5*(m_co2m_co2)deflog10_gamma_long_range(i_str:float)-float:长程静电项扩展 Debye-Hückellog10 γ± -A√I / (1 B a √I)。 单独使用时 γ± 随 I 单调下降——低浓度正确高浓度不成立。rootmath.sqrt(i_str)return-A_DH*root/(1.0B_ION*root)defln_gamma_local(i_str:float,c_local:floatC_LOCAL)-float:局部组成项NRTL 型的线性示意ln γ_local ≈ c * I。 真实 eNRTL 的局部项含相互作用参数由回归或电解质向导给出 这里的 c 只是示意值用来展示“叠加后趋势反转”这一结构性结论。returnc_local*i_strdefwater_activity(i_str:float,m_total:float)-float:用渗透系数近似水的活度ln a_w ≈ -φ * Σm_i * M_w / 1000。 为什么关心 a_w它进入水合路径 CO2 H2O ⇌ H HCO3- 的平衡 也决定气相水汽分压漏掉它会把再生塔的汽提蒸汽量算错。phi1.0-0.5*A_DH*math.sqrt(i_str)/(1.0math.sqrt(i_str))# 示意returnmath.exp(-phi*m_total*M_H2O/1000.0)defmain()-None:m_meamea_molality()print(f30 wt% MEA 的质量摩尔浓度示意{m_mea:.3f}mol/kg 水)print(f{负载:6s}{I(mol/kg):10s}{γ±(仅长程):12s}{γ±(长程局部):14s}f{a_w:8s}{p_CO2 相对倍数:14s})refNoneforloadingin(0.05,0.10,0.20,0.30,0.40,0.48):i_strionic_strength(loading,m_mea)g_lr10**log10_gamma_long_range(i_str)g_tot10**(log10_gamma_long_range(i_str)ln_gamma_local(i_str)/LN10)awwater_activity(i_str,i_str*2)relg_tot**2# 固定负载时 p_CO2 正比于 γ±^2 的示意关系ifrefisNone:refrelprint(f{loading:6.2f}{i_str:10.3f}{g_lr:12.3f}{g_tot:14.3f}f{aw:8.4f}{rel/ref:14.2f})print(\n判读只叠长程项时 γ± 随负载下降高浓度区不合理)print( 叠加局部组成项后 γ± 反转上升所需 CO2 平衡分压被整体抬高。)print( 忽略这一机制模型就会高估吸收能力、低估 SRD。)if__name____main__:main()预期输出与判读γ±(仅长程)随负载从约 0.69 降到约 0.56而γ±(长程局部)则从约 0.78 升到约 1.83p_CO2 相对倍数因此涨到 5 倍以上a_w从约 0.99 降到约 0.90。这段代码的核心结论有两条其一只叠长程项会给出物理上不合理的高浓度趋势这正是 eNRTL 必须叠加局部组成项的原因其二忽略活度修正的模型会系统性低估高负载下的平衡分压从而高估吸收能力、低估能耗。代码 2电解质反应集与物种表自洽性自检# 03_code2_reaction_species_check.py# 运行环境Python 3.10仅标准库 re# 用途对 MEA-H2O-CO2 体系的物种表与反应集做元素与电荷衡算自检。# 为什么这段代码重要不同热力学模型的差异首先来自“收了哪些物种、# 写了哪几个独立反应”。物种表都自洽不了后面的活度模型无从谈起。# 【说明】本代码不含任何拟合参数全部是可验证的化学计量衡算。importre# 物种表字符串为分子式省略 1第二个元素为电荷数。# 物种集合取自底账 [C5] 描述的平衡反应体系含 MEACOO- 与 HCO3- 生成。SPECIES{MEA:(C2H7NO,0),MEAH:(C2H8NO,1),MEACOO-:(C3H6NO3,-1),HCO3-:(CH1O3,-1),CO3^2-:(C1O3,-2),H:(H1,1),OH-:(H1O1,-1),H2O:(H2O1,0),CO2:(C1O2,0),}# 反应集负数反应物正数产物。系数为化学计量数。REACTIONS{R1 氨基甲酸盐生成:{MEA:-2,CO2:-1,MEACOO-:1,MEAH:1},R2 CO2 水合CO2 H2O:{CO2:-1,H2O:-1,H:1,HCO3-:1},R3 碳酸氢根二级解离:{HCO3-:-1,H:1,CO3^2-:1},R4 水的解离:{H2O:-1,H:1,OH-:1},R5 MEA 质子化:{MEA:-1,H:-1,MEAH:1},}_TOKENre.compile(r([A-Z][a-z]?)(\d*))defparse_formula(formula:str)-dict:把 C2H7NO 解析为 {C: 2, H: 7, N: 1, O: 1}。out{}forelement,countin_TOKEN.findall(formula):out[element]out.get(element,0)(int(count)ifcountelse1)returnoutdefcheck_reaction(name:str,stoich:dict)-bool:返回该反应是否同时满足元素衡算与电荷衡算。element_balance{}charge_balance0forsp,nuinstoich.items():formula,chargeSPECIES[sp]forelement,countinparse_formula(formula).items():element_balance[element]element_balance.get(element,0)nu*count charge_balancenu*charge bad_elements{k:vfork,vinelement_balance.items()ifv!0}ok(notbad_elements)andcharge_balance0detail通过ifokelsef不通过元素残差{bad_elements}, 电荷残差{charge_balance}print(f{name:26s}{detail})returnokdefmain()-None:print( 物种表 )forsp,(formula,charge)inSPECIES.items():print(f{sp:9s}{formula:8s}电荷{charge:d})print(\n 反应衡算自检 )all_okall(check_reaction(name,s)forname,sinREACTIONS.items())print(f\n全部反应自洽{all_ok})# 独立反应数 反应数本集合已手工去除线性相关组合。# 自由度提醒MEA-H2O-CO2 体系的平衡求解还需要每个独立反应的 K(T)# 而这些 K 在不同模型里来自“表观常数”或“活度模型 相平衡”两条完全不同的路径。print(f独立反应数本模型集{len(REACTIONS)})print(提示Kent-Eisenberg 用表观平衡常数合并表达上述反应)print( eNRTL / Deshmukh-Mather 则显式携带活度系数路径不同、结果自然不同。)if__name____main__:main()反直觉点很多人以为模型差异只出在活度系数公式上其实第一步差异是物种表与反应集。比如把水合路径写成CO2 OH⁻ ⇌ HCO₃⁻还是CO2 H₂O ⇌ H⁺ HCO₃⁻在显式活度模型里不影响平衡结果两者由水的解离关联但在半经验模型里会直接影响你拟合哪几个表观常数。代码 3热力学模型选型决策脚本# 03_code3_model_selection.py# 运行环境Python 3.10仅标准库# 用途把 1.7 节的决策路径落成可执行规则避免凭印象选模型。# 【说明】本脚本输出的是“选型建议与依据”不是计算结果# 所有依据均标注了出处未核实的地方一律回落到“以官方文档为准”。defchoose_model(quantitative:bool,has_electrolyte_pkg:bool,low_pressure_focus:bool,high_loading_focus:bool)-str:ifnotquantitative:return(Kent-Eisenberg半经验理由计算快、参数少注意其在【高负载】时精度相对下降[C6][C7])ifhas_electrolyte_pkg:base(eNRTL理由可计算全部离子与分子物种的活度系数规范引用 Chen Evans (1986)[C3]。若同时要做速率基塔还须配套黏度/表面张力/导热系数/扩散物性[C1])iflow_pressure_focusandhigh_loading_focus:base补充低分压可参考 Deshmukh-Mather[C7]高负载可参考 modified K-E[C6]eliflow_pressure_focus:base补充低 CO2 分压区间 Deshmukh-Mather 比 K-E 更准[C7]elifhigh_loading_focus:base补充高负载区间 modified K-E 优于 Deshmukh-Mather[C6]returnbasereturn(无内置电解质物性包可选路径用 CAPE-OPEN 接入外部电解质包或自建 eNRTL 并回归参数注意 DWSIM 的自动化属性包注册表中没有 eNRTL 专用包[S16]是否具备该能力请以官方文档为准)CASES[(快速量级扫描,dict(quantitativeFalse,has_electrolyte_pkgTrue,low_pressure_focusFalse,high_loading_focusTrue)),(有授权物性包的定量吸收塔设计,dict(quantitativeTrue,has_electrolyte_pkgTrue,low_pressure_focusFalse,high_loading_focusTrue)),(低分压烟气工况,dict(quantitativeTrue,has_electrolyte_pkgTrue,low_pressure_focusTrue,high_loading_focusFalse)),(仅开源无 eNRTL 环境,dict(quantitativeTrue,has_electrolyte_pkgFalse,low_pressure_focusFalse,high_loading_focusFalse)),]defmain()-None:fortitle,kwinCASES:print(f\n[{title}])print( 建议choose_model(**kw))if__name____main__:main()为什么把建议和结论分开选型脚本最容易被误用成自动决策器。它只负责把已知的适用区间整理成人能读的句子真正的判断仍取决于你有没有 VLE 数据、能不能拿到授权物性包、以及要不要上速率基。三、常见报错与排查现象 1文献综述里写eNRTL 由 Chen Song (2004) 提出被审稿人指出。根因误把 NRTL-SAC 当成 eNRTL。解法规范引用Chen Evans (1986)AIChE J. 32, 444–454及 1982 年单电解质原始文献 [C3][C4]Chen Song (2004) 是“Solubility Modeling with NRTL-SAC”IECR 43, 8354–8362[C8]。现象 2用 Kent-Eisenberg 算高负载工况结果明显偏离实测。根因K-E 类模型在高负载时精度相对下降 [C6][C7]。解法在高负载区改用 modified Kent-Eisenberg文献结论为在高负载下优于 Deshmukh-Mather [C6]或改用带活度系数的模型。现象 3把 Deshmukh-Mather 与 Kent-Eisenberg 的适用结论混用。根因两条结论方向不同——低 CO₂ 分压看 D-M [C7]高负载看 modified K-E [C6]。解法在文中成对引用并各自标注工况不要写成D-M 更准这种无条件的单句。现象 4温度扫描时ln K(T)外推到回归窗口之外K 值异常。根因ln K(T) a b/T c·lnT d·T的四个系数只在特定温度窗口内有效 [C3]。解法记录每个 K 参数的回归区间外推处显式标注未经验证的外推吸收侧与再生侧分别确认。现象 5在 DWSIM 里找 eNRTL 属性包找不到就以为软件装坏了。根因DWSIM 自动化属性包注册表中确实没有 eNRTL 专用包 [S16]。解法改用IdealElectrolytePropertyPackage做定性演示或通过CAPEOPENPropertyPackage接入外部包不要在文中断言 DWSIM 支持 eNRTL。四、动手练习练习 1长程项 vs 局部项运行代码 1先把C_LOCAL设为 0.0再设回 0.35观察γ±随负载的变化方向。成功判定C_LOCAL 0.0时γ±(长程局部)与γ±(仅长程)相同、随负载单调下降C_LOCAL 0.35时γ±(长程局部)随负载上升且p_CO2 相对倍数明显大于 1。请写下结论“eNRTL 必须同时含长程项与局部组成项缺一项趋势就会反。”练习 2反应集自检运行代码 2然后故意把 R1 的MEAH系数从 1 改成 2。成功判定程序应输出不通过元素残差 …, 电荷残差 …“。恢复后应全部显示通过”。这一步的目的是让你亲手确认反应集必须先通过元素与电荷衡算才谈得上拟合参数。练习 3选型演练运行代码 3为有授权物性包 低分压 高负载的工况手写一条新的CASES条目。成功判定输出应同时包含Deshmukh-Mather 更准与modified K-E 优于 D-M两句并各自带出处编号若只出现一句说明你把两个适用区间合并了。五、小结与下一篇预告本篇确定了热力学模型这道分水岭eNRTL 的假设与规范引用Chen Evans 1986、Kent-Eisenberg 与 Deshmukh-Mather 的适用区间低分压看 D-M、高负载看 modified K-E、ln K(T)形式的外推风险以及一条可执行的选型路径。三段代码分别演示了活度传导链、反应集自洽性自检与选型决策。下一篇04进入动手环节用 Aspen Plus 搭建本系列第一个平衡级 MEA 吸收塔——定义 MEA-H₂O-CO₂ 与离子物种、选择 eNRTL 路径、用 RadFrac 把填料高度抽象成理论级然后用Apwn.Document、Tree.FindNode和Node.Value把它脚本化并用贫液 49 ℃ 与 L/G 1.069 kg/kg 这两个已核实的量级做结果校核。本篇认知问题回显FAQQ1eNRTL 模型的假设是什么规范引用是哪一篇AeNRTL 基于同号离子互斥与局部电中性假设规范引用为 Chen Evans (1986)AIChE Journal 第 32 卷 444–454 页。Q2为什么 eNRTL 不能写成 Chen Song 2004A因为 Chen Song (2004) 是 NRTL-SAC 溶解度模型论文主题是溶解度建模与 eNRTL 不是同一个模型eNRTL 应引用 Chen Evans (1986)。Q3Kent-Eisenberg 与 Deshmukh-Mather 分别在什么区间更准A两组结论方向不同Deshmukh-Mather 在低 CO₂ 分压时更准modified Kent-Eisenberg 则在高负载下优于 Deshmukh-Mather两者必须成对引用。Q4Aspen 中平衡常数的温度关联式是什么形式AAspen 中平衡常数随温度采用 ln K(T) a b/T c·lnT d·T 形式四个系数仅在回归温度窗口内有效。Q5DWSIM 能否直接使用 eNRTL 计算 MEA 体系A不能直接使用。DWSIM 自动化属性包注册表中没有 eNRTL 专用包需自行开发或通过 CAPE-OPEN 接入外部电解质属性包。
燃烧后CO2捕集MEA吸收过程模拟教程03热力学模型的第一道分水岭——eNRTL、Kent-Eisenberg 与 Deshmukh-Mather版本声明块工具/软件Aspen Plus官方产品页未公示数字版本号以官方文档为准Python 3.10 及以上语言/环境Python仅使用标准库math、re、itertools无需第三方包本文目标掌握三类电解质热力学模型的假设、适用区间与选型判据并纠正一处流传甚广的错误引用一句话结论MEA-H₂O-CO₂ 体系的热力学模型分三道eNRTL基于同号离子互斥与局部电中性假设规范引用为 Chen Evans (1986), AIChE J. 32, 444–454、Kent-Eisenberg半经验表观平衡常数高负载精度下降、Deshmukh-MatherDebye-Hückel 型加相互作用参数低 CO₂ 分压时比 Kent-Eisenberg 更准选错哪一道都会让吸收塔的 CO₂ 平衡分压算错。〇、本篇要解决的认知问题eNRTL 的物理假设是什么它的规范引用到底是哪一篇文献为什么eNRTL Chen Song 2004这个说法是错的正确的对应关系是什么Kent-Eisenberg 与 Deshmukh-Mather 分别在哪一段工况更准能不能互相替代Aspen 中平衡常数的温度关联式ln K(T) a b/T c·lnT d·T意味着什么风险给定一个 MEA 吸收模拟任务怎么用一条决策路径选出热力学模型一、机制解析1.1 为什么热力学模型是第一道分水岭MEA 吸收塔的推动力来自气相 CO₂ 分压与液相 CO₂ 平衡分压之差。气相好算难的是液相溶液里同时存在 MEA、MEAH⁺、MEACOO⁻、HCO₃⁻、CO₃²⁻、H⁺、OH⁻、H₂O、CO₂ 等物种见第 01 篇。给定负载与温度液相 CO₂ 的平衡分压要由所有物种的活度共同决定。这就是分水岭同一个实验点换一套活度模型算出的平衡分压可以差出一大截进而把塔高、液气比、SRD 全部带偏。1.2 eNRTL假设、来源与规范引用eNRTLelectrolyte NRTL的核心假设是同号离子互斥与局部电中性它先由 Chen、Britt、Boston、Evans (1982) 针对单溶剂完全解离体系提出随后由Chen Evans (1986)推广到多组分水溶液体系 [C3][C4][C5]。它的能力边界是明确的可计算所有物种离子与分子、溶质与溶剂的活度系数——这正是 MEA 这种离子物种 中性物种共存体系最需要的性质。必须纠正的错误引用流传中常把 eNRTL 说成 “Chen Song 2004”。事实是Chen Song (2004) 的论文标题为“Solubility Modeling with NonRandom Two-Liquid Segment Activity Coefficient Model”IECR 43, 8354–8362即NRTL-SAC 溶解度模型与 eNRTL 不是同一个模型 [C8]。本系列统一规范引用为Chen Evans (1986)可并列 1982 年单电解质原始文献。这条纠偏请写进你自己的文献笔记里因为一旦引错后面对参数回归第 13 篇的讨论会全线失焦。最新进展方面2024 年AIChE Journal有工作给出用于 MEA-H₂O-CO₂ 的modified eNRTL 含过量焓解析式[C9]可作为该方向仍在演进的旁证引用但它不是本系列的基础假设来源。从工程视角看eNRTL 的贵体现在两处。第一处是参数数量长程项需要离子尺寸与介电常数相关参数局部组成项需要一组相互作用参数这些参数不能凭空给出必须由 VLE、溶解度或量热数据回归得到。第二处是参数的可辨识性MEA-H₂O-CO₂ 体系的负载、温度、压力三个维度高度耦合如果实验数据覆盖的负载区间过窄回归出的参数一旦离开该区间就可能给出错误趋势。这正是第 13 篇要专门处理参数回归与残差诊断的原因。反过来看Kent-Eisenberg 之所以长期存在于工业现场正是因为它把这两笔成本压缩成一组表观常数用精度换取了可用性——这个取舍本身没有对错只有适用边界。1.3 Kent-Eisenberg半经验的代价Kent-Eisenberg及其 modified 版本源头是 Kent Eisenberg 1976 年Hydrocarbon Processing的 “Better data for amine treating”本系列经二次文献 [C6][C7] 转引属于半经验平衡模型用表观平衡常数去拟合实验负载计算快、参数少。代价也很清楚在高负载时精度相对下降[C6][C7]。它的优点是真实存在的——当你只需要快速扫大量工况、或者手里没有任何 VLE 数据做回归时K-E 类模型是能立刻跑起来的选择。但请不要把它当成严格热力学模型来写。1.4 Deshmukh-Mather长程项加相互作用参数Deshmukh-Mather源头是 Deshmukh Mather 1981 年Chemical Engineering Science的工作同样经二次文献 [C7] 转引走的是另一条路用Debye-Hückel 型长程项 相互作用参数描述电解质非理想性。文献给出的对照结论是在低 CO₂ 分压时Deshmukh-Mather 比 Kent-Eisenberg 更准[C7]而另一方面modified Kent-Eisenberg 在高负载下优于 Deshmukh-Mather[C6]。这两句话必须成对使用、并各自带工况低分压看 D-M高负载看 modified K-E。把两条结论各自截一半、拼成D-M 更准或K-E 更准的单一结论是本主题下最常见的写作错误。1.5 三者对照维度eNRTLKent-Eisenberg / modified K-EDeshmukh-Mather类型活度系数模型γ-φ严格热力学半经验平衡模型活度系数 / 电解质平衡模型核心构造同号离子互斥 局部电中性长程项 局部组成项表观平衡常数拟合负载Debye-Hückel 型 相互作用参数规范引用Chen Evans (1986)及 1982 原始文献Kent Eisenberg (1976)本系列经二次文献转引Deshmukh Mather (1981)本系列经二次文献转引优势工况需要全物种活度、需与速率基配套时快速扫描、参数少、高负载modified 版低 CO₂ 分压已知短板参数多、回归难无授权软件时难落地高负载精度下降高负载下不如 modified K-E软件支持Aspen Plus 内置 property method含 MEA 专用物性 insertK-E 类与 D-M 类可由用户/电解质模板实现 [C1][C2][S4]同上同上关于开源DWSIM 的自动化属性包注册表中没有 eNRTL 专用包有NRTLPropertyPackage、IdealElectrolytePropertyPackage、CAPEOPENPropertyPackage等[S16]因此禁止断言 DWSIM 支持 eNRTL若要在开源侧使用 eNRTL路径是自建或通过 CAPE-OPEN 接入外部属性包第 16 篇。1.6ln K(T)的形式与温度外推风险在 Aspen 中平衡常数随温度的变化采用ln K(T) a b/T c·lnT d·T形式 [C3]。这个形式本身很常见但它有明确的适用区间四个系数是在某个温度窗口内回归出来的把它外推到窗口之外例如从吸收温度外推到再生温度、或更高会得到没有物理意义的 K 值。因此吸收侧与再生侧的平衡常数参数最好分别在其各自温度窗口内确认涉及温度扫描的工况试验第 12 篇必须记录 K 参数的回归区间任何用一条 K(T) 曲线覆盖 40–120 ℃ 甚至更宽的写法都要标注为未经验证的外推。1.7 选型决策路径任务是否需要定量预测吸收/再生性能 ├─ 否只要量级、要快───────────────► Kent-Eisenberg半经验参数少 └─ 是 ├─ 是否有授权软件的内置电解质物性包 │ ├─ 是 ─► eNRTLChen Evans 1986若还要速率基塔必须配套传递物性 [C1] │ └─ 否 ─► 用 CAPE-OPEN 接入外部电解质包或自建第 16 篇 └─ 关注区间是哪一个 ├─ 低 CO₂ 分压 ─► Deshmukh-Mather 更准 [C7] └─ 高负载 ─► modified Kent-Eisenberg 优于 Deshmukh-Mather [C6]注意这条路径不是二选一之后互不相干速率基塔必须显式要求黏度、表面张力、导热系数、二元扩散系数等传递物性[C1]这意味着一旦你选择速率基第 05、08 篇eNRTL 类的完整物性包几乎是必选配置而不是可选项。二、完整代码与逐行剖析代码 1离子强度如何把活度系数传导到平衡分压# 03_code1_enrtl_activity_demo.py# 运行环境Python 3.10仅标准库 math# 用途演示 eNRTL 为什么必须由“长程项 局部组成项”两部分组成# 只用长程静电项活度系数会随离子强度单调下降高浓度区物理上不合理# 叠加局部组成项后出现“盐析”方向的反转并把 CO2 平衡分压整体抬高。# 【重要】本脚本用简化 Debye-Hückel 型长程项与线性局部项替代 eNRTL 的# Pitzer-Debye-Hückel 项与 NRTL 局部项物种分布、A_DH、B_ION、C_LOCAL# 全部是【示意参数】不能用于定量设计。真实计算须由 Aspen Plus 的# 电解质物性包含 eNRTL完成见第 04 篇。importmath# ---------------- 示意参数区 ----------------MEA_WT0.30# 底账基准 30 wt% MEA[C19][C14][C5]M_MEA61.08# MEA 摩尔质量, g/molM_H2O18.015# 水摩尔质量, g/molA_DH0.509# Debye-Hückel 常数25 ℃ 水溶液常见值示意B_ION1.5# 离子尺寸参数 B*a示意C_LOCAL0.35# 局部组成项示意系数对应 eNRTL 的 NRTL 相互作用参数LN10math.log(10.0)defmea_molality(wt:floatMEA_WT)-float:把 30 wt% MEA 换算成质量摩尔浓度 mol/kg 水示意换算。n_meawt/M_MEA# 每 1 g 溶液中 MEA 的摩尔数m_water_g1.0-wt# 每 1 g 溶液中水的质量, greturnn_mea/(m_water_g/1000.0)defionic_strength(loading:float,m_mea:float)-float:示意假设负载态 CO2 全部形成 MEACOO- 并伴随等量 MEAH 则 I 0.5*(m_cation m_anion) m_co2。 真实体系里 HCO3-/CO3^2-/OH- 都会参与必须由平衡求解给出。m_co2m_mea*loadingreturn0.5*(m_co2m_co2)deflog10_gamma_long_range(i_str:float)-float:长程静电项扩展 Debye-Hückellog10 γ± -A√I / (1 B a √I)。 单独使用时 γ± 随 I 单调下降——低浓度正确高浓度不成立。rootmath.sqrt(i_str)return-A_DH*root/(1.0B_ION*root)defln_gamma_local(i_str:float,c_local:floatC_LOCAL)-float:局部组成项NRTL 型的线性示意ln γ_local ≈ c * I。 真实 eNRTL 的局部项含相互作用参数由回归或电解质向导给出 这里的 c 只是示意值用来展示“叠加后趋势反转”这一结构性结论。returnc_local*i_strdefwater_activity(i_str:float,m_total:float)-float:用渗透系数近似水的活度ln a_w ≈ -φ * Σm_i * M_w / 1000。 为什么关心 a_w它进入水合路径 CO2 H2O ⇌ H HCO3- 的平衡 也决定气相水汽分压漏掉它会把再生塔的汽提蒸汽量算错。phi1.0-0.5*A_DH*math.sqrt(i_str)/(1.0math.sqrt(i_str))# 示意returnmath.exp(-phi*m_total*M_H2O/1000.0)defmain()-None:m_meamea_molality()print(f30 wt% MEA 的质量摩尔浓度示意{m_mea:.3f}mol/kg 水)print(f{负载:6s}{I(mol/kg):10s}{γ±(仅长程):12s}{γ±(长程局部):14s}f{a_w:8s}{p_CO2 相对倍数:14s})refNoneforloadingin(0.05,0.10,0.20,0.30,0.40,0.48):i_strionic_strength(loading,m_mea)g_lr10**log10_gamma_long_range(i_str)g_tot10**(log10_gamma_long_range(i_str)ln_gamma_local(i_str)/LN10)awwater_activity(i_str,i_str*2)relg_tot**2# 固定负载时 p_CO2 正比于 γ±^2 的示意关系ifrefisNone:refrelprint(f{loading:6.2f}{i_str:10.3f}{g_lr:12.3f}{g_tot:14.3f}f{aw:8.4f}{rel/ref:14.2f})print(\n判读只叠长程项时 γ± 随负载下降高浓度区不合理)print( 叠加局部组成项后 γ± 反转上升所需 CO2 平衡分压被整体抬高。)print( 忽略这一机制模型就会高估吸收能力、低估 SRD。)if__name____main__:main()预期输出与判读γ±(仅长程)随负载从约 0.69 降到约 0.56而γ±(长程局部)则从约 0.78 升到约 1.83p_CO2 相对倍数因此涨到 5 倍以上a_w从约 0.99 降到约 0.90。这段代码的核心结论有两条其一只叠长程项会给出物理上不合理的高浓度趋势这正是 eNRTL 必须叠加局部组成项的原因其二忽略活度修正的模型会系统性低估高负载下的平衡分压从而高估吸收能力、低估能耗。代码 2电解质反应集与物种表自洽性自检# 03_code2_reaction_species_check.py# 运行环境Python 3.10仅标准库 re# 用途对 MEA-H2O-CO2 体系的物种表与反应集做元素与电荷衡算自检。# 为什么这段代码重要不同热力学模型的差异首先来自“收了哪些物种、# 写了哪几个独立反应”。物种表都自洽不了后面的活度模型无从谈起。# 【说明】本代码不含任何拟合参数全部是可验证的化学计量衡算。importre# 物种表字符串为分子式省略 1第二个元素为电荷数。# 物种集合取自底账 [C5] 描述的平衡反应体系含 MEACOO- 与 HCO3- 生成。SPECIES{MEA:(C2H7NO,0),MEAH:(C2H8NO,1),MEACOO-:(C3H6NO3,-1),HCO3-:(CH1O3,-1),CO3^2-:(C1O3,-2),H:(H1,1),OH-:(H1O1,-1),H2O:(H2O1,0),CO2:(C1O2,0),}# 反应集负数反应物正数产物。系数为化学计量数。REACTIONS{R1 氨基甲酸盐生成:{MEA:-2,CO2:-1,MEACOO-:1,MEAH:1},R2 CO2 水合CO2 H2O:{CO2:-1,H2O:-1,H:1,HCO3-:1},R3 碳酸氢根二级解离:{HCO3-:-1,H:1,CO3^2-:1},R4 水的解离:{H2O:-1,H:1,OH-:1},R5 MEA 质子化:{MEA:-1,H:-1,MEAH:1},}_TOKENre.compile(r([A-Z][a-z]?)(\d*))defparse_formula(formula:str)-dict:把 C2H7NO 解析为 {C: 2, H: 7, N: 1, O: 1}。out{}forelement,countin_TOKEN.findall(formula):out[element]out.get(element,0)(int(count)ifcountelse1)returnoutdefcheck_reaction(name:str,stoich:dict)-bool:返回该反应是否同时满足元素衡算与电荷衡算。element_balance{}charge_balance0forsp,nuinstoich.items():formula,chargeSPECIES[sp]forelement,countinparse_formula(formula).items():element_balance[element]element_balance.get(element,0)nu*count charge_balancenu*charge bad_elements{k:vfork,vinelement_balance.items()ifv!0}ok(notbad_elements)andcharge_balance0detail通过ifokelsef不通过元素残差{bad_elements}, 电荷残差{charge_balance}print(f{name:26s}{detail})returnokdefmain()-None:print( 物种表 )forsp,(formula,charge)inSPECIES.items():print(f{sp:9s}{formula:8s}电荷{charge:d})print(\n 反应衡算自检 )all_okall(check_reaction(name,s)forname,sinREACTIONS.items())print(f\n全部反应自洽{all_ok})# 独立反应数 反应数本集合已手工去除线性相关组合。# 自由度提醒MEA-H2O-CO2 体系的平衡求解还需要每个独立反应的 K(T)# 而这些 K 在不同模型里来自“表观常数”或“活度模型 相平衡”两条完全不同的路径。print(f独立反应数本模型集{len(REACTIONS)})print(提示Kent-Eisenberg 用表观平衡常数合并表达上述反应)print( eNRTL / Deshmukh-Mather 则显式携带活度系数路径不同、结果自然不同。)if__name____main__:main()反直觉点很多人以为模型差异只出在活度系数公式上其实第一步差异是物种表与反应集。比如把水合路径写成CO2 OH⁻ ⇌ HCO₃⁻还是CO2 H₂O ⇌ H⁺ HCO₃⁻在显式活度模型里不影响平衡结果两者由水的解离关联但在半经验模型里会直接影响你拟合哪几个表观常数。代码 3热力学模型选型决策脚本# 03_code3_model_selection.py# 运行环境Python 3.10仅标准库# 用途把 1.7 节的决策路径落成可执行规则避免凭印象选模型。# 【说明】本脚本输出的是“选型建议与依据”不是计算结果# 所有依据均标注了出处未核实的地方一律回落到“以官方文档为准”。defchoose_model(quantitative:bool,has_electrolyte_pkg:bool,low_pressure_focus:bool,high_loading_focus:bool)-str:ifnotquantitative:return(Kent-Eisenberg半经验理由计算快、参数少注意其在【高负载】时精度相对下降[C6][C7])ifhas_electrolyte_pkg:base(eNRTL理由可计算全部离子与分子物种的活度系数规范引用 Chen Evans (1986)[C3]。若同时要做速率基塔还须配套黏度/表面张力/导热系数/扩散物性[C1])iflow_pressure_focusandhigh_loading_focus:base补充低分压可参考 Deshmukh-Mather[C7]高负载可参考 modified K-E[C6]eliflow_pressure_focus:base补充低 CO2 分压区间 Deshmukh-Mather 比 K-E 更准[C7]elifhigh_loading_focus:base补充高负载区间 modified K-E 优于 Deshmukh-Mather[C6]returnbasereturn(无内置电解质物性包可选路径用 CAPE-OPEN 接入外部电解质包或自建 eNRTL 并回归参数注意 DWSIM 的自动化属性包注册表中没有 eNRTL 专用包[S16]是否具备该能力请以官方文档为准)CASES[(快速量级扫描,dict(quantitativeFalse,has_electrolyte_pkgTrue,low_pressure_focusFalse,high_loading_focusTrue)),(有授权物性包的定量吸收塔设计,dict(quantitativeTrue,has_electrolyte_pkgTrue,low_pressure_focusFalse,high_loading_focusTrue)),(低分压烟气工况,dict(quantitativeTrue,has_electrolyte_pkgTrue,low_pressure_focusTrue,high_loading_focusFalse)),(仅开源无 eNRTL 环境,dict(quantitativeTrue,has_electrolyte_pkgFalse,low_pressure_focusFalse,high_loading_focusFalse)),]defmain()-None:fortitle,kwinCASES:print(f\n[{title}])print( 建议choose_model(**kw))if__name____main__:main()为什么把建议和结论分开选型脚本最容易被误用成自动决策器。它只负责把已知的适用区间整理成人能读的句子真正的判断仍取决于你有没有 VLE 数据、能不能拿到授权物性包、以及要不要上速率基。三、常见报错与排查现象 1文献综述里写eNRTL 由 Chen Song (2004) 提出被审稿人指出。根因误把 NRTL-SAC 当成 eNRTL。解法规范引用Chen Evans (1986)AIChE J. 32, 444–454及 1982 年单电解质原始文献 [C3][C4]Chen Song (2004) 是“Solubility Modeling with NRTL-SAC”IECR 43, 8354–8362[C8]。现象 2用 Kent-Eisenberg 算高负载工况结果明显偏离实测。根因K-E 类模型在高负载时精度相对下降 [C6][C7]。解法在高负载区改用 modified Kent-Eisenberg文献结论为在高负载下优于 Deshmukh-Mather [C6]或改用带活度系数的模型。现象 3把 Deshmukh-Mather 与 Kent-Eisenberg 的适用结论混用。根因两条结论方向不同——低 CO₂ 分压看 D-M [C7]高负载看 modified K-E [C6]。解法在文中成对引用并各自标注工况不要写成D-M 更准这种无条件的单句。现象 4温度扫描时ln K(T)外推到回归窗口之外K 值异常。根因ln K(T) a b/T c·lnT d·T的四个系数只在特定温度窗口内有效 [C3]。解法记录每个 K 参数的回归区间外推处显式标注未经验证的外推吸收侧与再生侧分别确认。现象 5在 DWSIM 里找 eNRTL 属性包找不到就以为软件装坏了。根因DWSIM 自动化属性包注册表中确实没有 eNRTL 专用包 [S16]。解法改用IdealElectrolytePropertyPackage做定性演示或通过CAPEOPENPropertyPackage接入外部包不要在文中断言 DWSIM 支持 eNRTL。四、动手练习练习 1长程项 vs 局部项运行代码 1先把C_LOCAL设为 0.0再设回 0.35观察γ±随负载的变化方向。成功判定C_LOCAL 0.0时γ±(长程局部)与γ±(仅长程)相同、随负载单调下降C_LOCAL 0.35时γ±(长程局部)随负载上升且p_CO2 相对倍数明显大于 1。请写下结论“eNRTL 必须同时含长程项与局部组成项缺一项趋势就会反。”练习 2反应集自检运行代码 2然后故意把 R1 的MEAH系数从 1 改成 2。成功判定程序应输出不通过元素残差 …, 电荷残差 …“。恢复后应全部显示通过”。这一步的目的是让你亲手确认反应集必须先通过元素与电荷衡算才谈得上拟合参数。练习 3选型演练运行代码 3为有授权物性包 低分压 高负载的工况手写一条新的CASES条目。成功判定输出应同时包含Deshmukh-Mather 更准与modified K-E 优于 D-M两句并各自带出处编号若只出现一句说明你把两个适用区间合并了。五、小结与下一篇预告本篇确定了热力学模型这道分水岭eNRTL 的假设与规范引用Chen Evans 1986、Kent-Eisenberg 与 Deshmukh-Mather 的适用区间低分压看 D-M、高负载看 modified K-E、ln K(T)形式的外推风险以及一条可执行的选型路径。三段代码分别演示了活度传导链、反应集自洽性自检与选型决策。下一篇04进入动手环节用 Aspen Plus 搭建本系列第一个平衡级 MEA 吸收塔——定义 MEA-H₂O-CO₂ 与离子物种、选择 eNRTL 路径、用 RadFrac 把填料高度抽象成理论级然后用Apwn.Document、Tree.FindNode和Node.Value把它脚本化并用贫液 49 ℃ 与 L/G 1.069 kg/kg 这两个已核实的量级做结果校核。本篇认知问题回显FAQQ1eNRTL 模型的假设是什么规范引用是哪一篇AeNRTL 基于同号离子互斥与局部电中性假设规范引用为 Chen Evans (1986)AIChE Journal 第 32 卷 444–454 页。Q2为什么 eNRTL 不能写成 Chen Song 2004A因为 Chen Song (2004) 是 NRTL-SAC 溶解度模型论文主题是溶解度建模与 eNRTL 不是同一个模型eNRTL 应引用 Chen Evans (1986)。Q3Kent-Eisenberg 与 Deshmukh-Mather 分别在什么区间更准A两组结论方向不同Deshmukh-Mather 在低 CO₂ 分压时更准modified Kent-Eisenberg 则在高负载下优于 Deshmukh-Mather两者必须成对引用。Q4Aspen 中平衡常数的温度关联式是什么形式AAspen 中平衡常数随温度采用 ln K(T) a b/T c·lnT d·T 形式四个系数仅在回归温度窗口内有效。Q5DWSIM 能否直接使用 eNRTL 计算 MEA 体系A不能直接使用。DWSIM 自动化属性包注册表中没有 eNRTL 专用包需自行开发或通过 CAPE-OPEN 接入外部电解质属性包。
