燃烧后CO2捕集MEA吸收过程模拟教程04第一个可跑的仿真——用 Aspen Plus 搭一个平衡级 MEA 吸收塔版本声明块工具/软件Aspen Plus官方产品页未公示数字版本号以官方文档为准RadFrac 为 Aspen Plus 最常用的平衡级多级分离模型速率基在其上扩展语言/环境Python 3.10代码 1 需 Windows pywin32 本机已安装 Aspen Plus 桌面版代码 2、3 纯标准库任何平台均可运行本文目标跑通本系列第一个平衡级 MEA 吸收塔并建立结果必须用已核实锚点校核的习惯一句话结论第一个可跑的 MEA 吸收塔应当由 RadFrac 平衡级模型给出初值脚本层用Apwn.Document打开 .apw、用Tree.FindNode(r\Data\Streams\FEED\Input\TOTFLOW)改流量、Reinit后运行并把贫液进塔温度 49 ℃ 与液气比 1.069 kg/kg 当作量级校核锚点而不是精确目标。〇、本篇要解决的认知问题MEA-H₂O-CO₂ 体系在 Aspen Plus 里要定义哪些组分与离子物种为什么物性方法要走 eNRTL 路径它和第 03 篇的分水岭结论怎么衔接怎么用 RadFrac 搭一个平衡级吸收塔把 15–24 m 填料高度抽象成理论级COM 自动化里Tree.FindNode、Node.Value、Reinit、Run2具体怎么用哪些写法不能被断言为唯一正确结果该拿什么校核为什么贫液 49 ℃ 与 L/G 1.069 kg/kg 只能用作量级参照一、机制解析1.1 建模路线总览组分定义 物性方法 单元模型 求解与校核 MEA / H2O / CO2 → eNRTL 电解质路径 → RadFrac 平衡级吸收塔 → 收敛 → 与 49 ℃ / 1.069 kg/kg 对照 MEACOO- / HCO3- 填料高度 → 理论级 ↓ 等离子物种 再切速率基第 08 篇这五步的顺序不能颠倒物种没定义全物性方法就无法建立物性方法没定RadFrac 的收敛结果没有意义没有校核锚点收敛了也不知道对不对。1.2 组分与电解质物种底账 [C5] 明确 MEA 水溶液吸收 CO₂ 时同时生成 MEACOO⁻ 与 HCO₃⁻并给出平衡反应式 2–6。因此组分定义至少要覆盖类别物种分子中性MEA、H₂O、CO₂离子MEACOO⁻、MEAH⁺、HCO₃⁻、CO₃²⁻、H⁺、OH⁻注意两点第一不要只定义 MEA / H₂O / CO₂ 三个分子就跑——没有离子物种MEA 与 CO₂ 反应生成氨基甲酸盐这件事在模型里就根本没有表达第二反应式必须先通过元素与电荷衡算第 03 篇代码 2 做的正是这件事否则后续的平衡计算会给出看似收敛、实则无意义的解。关于Aspen Plus 的电解质向导它的具体点击步骤与参数生成流程不在本系列可核实的资料范围内底账附录 A 已把电解质向导生成 eNRTL 参数的官方步骤文档列为待补方向本系列不做逐步描述请以官方文档为准。1.3 物性方法走 eNRTL 路径第 03 篇的结论在这里落地需要全物种活度系数时选 eNRTL规范引用为 Chen Evans (1986) [C3]。Aspen Plus 内置电解质物性包含 eNRTL 及 MEA 专用物性 insert [C1][C2][S4]。三点提示平衡常数在 Aspen 中采用ln K(T) a b/T c·lnT d·T形式 [C3]四个系数都有回归温度窗口吸收侧与再生侧不要共用一套外推如果后续要切速率基塔必须同时具备黏度、表面张力、导热系数、二元扩散系数等传递物性 [C1]这也是本系列建议一开始就选完整电解质物性包的原因CCSI 的MEA Steady State Model是 Aspen Plus 的附加模型含物性、水力学、界面面积、动力学四类子模型不是Aspen 内置的标准物性包 [C14]写文档时不要混淆。1.4 用 RadFrac 表达吸收塔填料高度 → 理论级RadFrac 是 Aspen Plus 最常用的平衡级多级分离模型速率基正是在其上扩展而来 [S3]。所以第一步先用平衡级搭塔是官方认可的路径AspenTech 官方课程明确包含把 RadFrac 平衡级 CO₂ 吸收塔转换为使用 Aspen Rate-Based Distillation 的塔[S15]。平衡级塔里塔高这一物理信息被压缩成理论级数。真实装置的填料高度是 15 mSherman FEED两塔并联与 24 mTCM三床 1266 m[C14][C19]。这里有一个非常有用的量级检验若把 15–24 m 填料折算到少数几个理论级上等价 HETP 会高达 5–8 m明显高于常规精馏填料塔的常见量级此为基于 [C14][C19] 数值的量级判断。这个反差本身就是平衡级抽象在 MEA 吸收塔上很粗糙的直接证据也是第 05、08 篇要走向速率基的动机。1.5 COM 对象模型真实名称与必须加的限定元素真实名称使用要点COM 对象 ProgIDApwn.Documentwin32com.client.gencache.EnsureDispatch(Apwn.Document)需在本机验证带版本后缀的 ProgID 未获官方确认打开文档InitFromArchive2(path)从 .apw / .bkp 启动节点树Application.TreeTree.FindNode(path)按路径取节点如\Data\Streams\FEED\Input\TOTFLOW遍历子元素Tree.Elements(Data).Elements(Blocks)/.Elements(name)用于探路避免硬编码输出路径读写数值Node.Value如\Data\Blocks\B1\Input\BASIS_RR运行Application.Engine.Run2()等价运行按钮Application.Run2()亦见于案例重置Application.Reinit()/Engine.Reinit(IAP_REINIT_SIMULATION)改动输入后使用官方 VB 示例给出go_Simulation.Engine.Reinit(IAP_REINIT_SIMULATION)另有IAP_REINIT_BLOCK、IAP_REINIT_CONVERGENCE停止Application.Stop()等价停止状态检查Application.COMPSTATUS/Node.COMPSTATUS输入完整性与可运行性检查弹窗抑制Application.SuppressDialogs True/False批量脚本必须设置可见性Application.Visible True/False关界面可提速导出Application.Export(code, filename)代码 1.bkp、2.rep、3.sum、4.inp、6运行消息、9/10动态报告、11DXF该组代码取自第三方 API 封装以官方 User Guide 第 38 章为准官方文档《Aspen Plus User Guide》第 38 章 Happ ActiveX interfaceAspenTech 支持站需账号关于Run2的写法真实案例中Application.Run2()与Application.Engine.Run2()都被使用过 [S17][S20]因此本系列不断言唯一正确写法写法上采用先试一种、异常时回退另一种并把最终结果交回本机验证。关于表格中方法的来源与时效Apwn.Document、InitFromArchive2、FindNode、Run2、导出代码这一组写法来自第三方可运行示例仓库随作者更新可能失效并有官方知识库示例佐证 [S17][S18][S19][S20]。官方完整的 ActiveX 方法清单在《Aspen Plus User Guide》第 38 章属于需要账号的支持文档 [S18][S21]因此本篇所有调用都标注为需在本机验证不得当作官方 SDK 的固定签名使用。关于工程坑COM 调用依赖本机安装的桌面版 Aspen Plusgencache.EnsureDispatch会生成类型库包装无桌面安装的服务器/容器环境不可用不设SuppressDialogs True时收敛失败弹窗会挂死脚本 [S17]。这些是第 11 篇批量自动化要系统解决的问题。1.6 校核锚点已核实、可作量级参照的锚点是吸收塔进液贫液温度49 ℃Sherman FEED35 wt% MEA[C14]液气比 L/G1.069 kg/kgSherman FEED[C14]富液负载 0.275–0.476、贫液负载 0.075–0.239 mol CO₂/mol MEA35 wt% MEA、CCGT 烟气[C14]。为什么只能当量级参照这三个值绑定的是 35 wt% MEA、特定烟气与特定装置配置Sherman FEED。如果你的模型用的是 30 wt% MEA 或不同烟气组成算出的 L/G 与 1.069 有差异是正常的差异必须回到工况层面解释。把锚点当必须命中的目标值去反向调参是这一篇最危险的操作。二、完整代码与逐行剖析代码 1用 Python COM 驱动 Aspen Plus需本机验证# 04_code1_aspen_com_absorber.py# 运行环境Windows Python 3.10 pywin32 本机已安装 Aspen Plus 桌面版# 前置条件准备一个已建好的 .apw其中含一台 RadFrac 吸收塔与一条名为 FEED 的进料流股# 【重要】Apwn.Document 与 Run2 / Reinit 的写法来自第三方可运行实现与官方示例佐证# 官方完整 ActiveX 清单在《Aspen Plus User Guide》第 38 章需账号。# 因此本脚本的每一处调用都需在本机验证不能断言唯一正确写法。importosimportwin32com.clientaswin32# 需 pip install pywin32APW_PATHrC:\path\to\mea_absorber.apw# 换成你自己的文件路径FEED_TOTAL_FLOW10000.0# 烟气总摩尔流率, kmol/h示意工况defget_node(aspen,path):按路径取节点取不到就抛出可读错误。 为什么先判 None树路径是硬编码的版本迁移后路径可能消失 直接 .Value 会抛出不带上下文的 COM 异常排查成本极高。nodeaspen.Tree.FindNode(path)ifnodeisNone:raiseRuntimeError(f节点不存在{path}升级后路径可能变化先做路径探测)returnnodedefmain():aspenwin32.gencache.EnsureDispatch(Apwn.Document)# 依赖本机注册的 COM 服务器try:aspen.InitFromArchive2(APW_PATH)# 从 .apw / .bkp 启动仿真aspen.VisibleFalse# 关界面批量跑时可提速aspen.SuppressDialogsTrue# 必须收敛失败弹窗会阻塞甚至挂死脚本# 输入完整性检查。不同接口封装下该属性的可读性不同读不到就跳过# 不让一个可选的诊断信息阻断整条流程。try:print(输入检查 COMPSTATUS ,aspen.COMPSTATUS)exceptExceptionasexc:print(COMPSTATUS 不可读以所用版本为准,exc)# 修改一个“有出处”的输入节点进料总流率flow_nodeget_node(aspen,r\Data\Streams\FEED\Input\TOTFLOW)flow_node.ValueFEED_TOTAL_FLOWprint(f已把 FEED 的 TOTFLOW 设为{flow_node.Value}kmol/h)aspen.Reinit()# 重置让输入改动生效# 运行Application.Run2() 与 Application.Engine.Run2() 两种写法都见诸案例# 这里先试 Engine.Run2()异常时回退到 Run2()最终以本机验证结果为准。try:aspen.Engine.Run2()exceptException:aspen.Run2()try:print(运行后 COMPSTATUS ,aspen.COMPSTATUS)exceptException:pass# 读结果不硬编码输出节点名而是遍历该流股 Output 下的子节点。# 为什么遍历Output 下的节点名随模型与版本变化硬编码是脚本升级后失效的头号原因。try:outaspen.Tree.FindNode(r\Data\Streams\RICH\Output)forelinout.Elements:print(f RICH.{el.Name}{el.Value})exceptExceptionasexc:print(遍历输出节点失败不同封装/版本遍历方式不同需在本机验证,exc)finally:# 无论成功失败都导出报告便于事后核对导出代码含义以官方 User Guide 第 38 章为准try:reportos.path.join(os.path.dirname(APW_PATH),mea_absorber.rep)aspen.Export(2,report)print(已导出报告,report)exceptExceptionasexc:print(导出报告失败,exc)if__name____main__:main()反直觉默认值SuppressDialogs的默认值是不抑制弹窗。这在交互界面上没问题但在脚本里一个是否继续的对话框就意味着无限等待。批量工况脚本必须把它显式设为True。代码 2不依赖 Aspen 的平衡级吸收塔教学模型# 04_code2_equilibrium_absorber.py# 运行环境Python 3.10仅标准库任何平台可运行# 用途用“平衡级 操作线”逐步追级法McCabe-Thiele 式估算 MEA 吸收塔的理论级数# 并与 15–24 m 真实填料高度对照理解“理论级”抽象在 MEA 体系里有多粗糙。# 【重要】平衡关系 y_eq C_ALPHA * x^2 / (X_MAX - x) 是【示意参数】形式# 不是 MEA-H2O-CO2 的真实平衡数据真实平衡分压须由第 03 篇的电解质模型给出。# 本脚本只用于理解“理论级 → 塔高”的换算逻辑。X_MAX1.0# 示意负载上限, mol CO2/mol MEAC_ALPHA0.041# 示意平衡关系系数R_MEA_GAS0.15# 示意MEA 摩尔流率 / 烟气总摩尔流率Y_IN0.035# 入口烟气 CO2 摩尔分数底账烟气量级 3.5-15% 的低端, [C14][C19]Y_OUT0.0035# 出口烟气 CO2 摩尔分数对应 90% 捕集底账经典基准 [C19][C14]X_LEAN0.23# 贫液负载, mol CO2/mol MEA底账实测 0.075-0.239 区间内, [C14]defy_eq(x:float)-float:气相 CO2 摩尔分数与液相负载的示意平衡关系单调递增x 趋近 X_MAX 时急剧上升。returnC_ALPHA*x**2/(X_MAX-x)defx_from_y(y:float,lo:float0.0,hi:floatX_MAX-1e-6,tol:float1e-10)-float:对 y_eq 做二分反解。 为什么要反解而不是解析解真实电解质模型的平衡关系没有解析式 工程代码一律用数值反解这里提前用最简单的函数把流程走通。for_inrange(200):mid0.5*(lohi)ify_eq(mid)y:lomidelse:himidifhi-lotol:breakreturn0.5*(lohi)defstages_needed(y_in:float,y_out:float,x_in:float)-tuple:从塔顶往下逐步追级直到进入塔底的气相浓度达到 y_in。 返回 (理论级数, 富液负载)。 操作线由全塔物料衡算给出y y_out R_MEA_GAS * (x - x_in) 其中斜率就是 MEA 摩尔流率与烟气摩尔流率之比——这一点常被写错成倒数。xx_from_y(y_out)# 第 1 级出口气与离开该级的液相平衡n1richxwhileTrue:y_nexty_outR_MEA_GAS*(x-x_in)# 进入下一级更靠下的气相组成ify_nexty_in:returnn,rich n1ifn100:# 防死循环超过 100 级说明工况不可行return-1,rich xx_from_y(y_next)richxdefmain()-float:n,richstages_needed(Y_IN,Y_OUT,X_LEAN)print( 平衡级吸收塔逐步追级示意平衡关系 )print(f入口 y_CO2 {Y_IN:.4f}出口 y_CO2 {Y_OUT:.4f}捕集率{(1-Y_OUT/Y_IN):.1%})print(f贫液负载 x_in {X_LEAN:.3f}mol CO2/mol MEA示意)print(f理论级数 N {n}富液负载 x_out {rich:.3f}mol CO2/mol MEA)# 用底账已核实的填料高度换算等价的 HETP作为“平衡级抽象有多粗糙”的量级检验forheightin(15.0,24.0):print(f若真实填料高度{height:.0f}m[C14][C19]由{n}个理论级承担f等价 HETP ≈{height/n:.1f}m/理论级)returnrichif__name____main__:main()预期输出与判读在示意参数下90% 捕集会给出3 个理论级与约0.49的富液负载。0.49 略高于底账实测的富液上限 0.476 mol CO₂/mol MEA [C14]——这不是程序算错而是提醒你示意平衡关系不是真实数据只能用来理解换算逻辑不能当设计结果。把 15 m 与 24 m 折算到 3 个理论级等价 HETP 约 5.0 与 8.0 m/理论级远高于常规精馏填料塔的量级。这不是模型错了而是**理论级这个抽象把反应吸收的传质阻力全部藏起来了**。代码 3结果校核脚本锚点对照不下结论# 04_code3_cross_check.py# 运行环境Python 3.10仅标准库# 用途把仿真或教学模型算出的关键量与底账已核实的锚点并排对照输出“在区间内/外”。# 【重要】锚点绑定的是 35 wt% MEA 与 Sherman FEED / CCGT 工况[C14]# 只作量级参照落在区间外不等于模型错必须回到工况差异解释禁止反向调参凑数。# 锚点全部来自底账 §4逐条带条件ANCHORS[(贫液进塔温度, ℃,49.0,49.0,Sherman FEED, 35 wt% MEA [C14]),(L/G, kg/kg,1.069,1.069,Sherman FEED [C14]),(富液负载, mol/mol,0.275,0.476,35 wt% MEA, CCGT 烟气 [C14]),(贫液负载, mol/mol,0.075,0.239,35 wt% MEA, CCGT 烟气 [C14]),(捕集率经典基准,0.90,0.90,90% 为经典设计基准 [C19][C14]),]defcheck(name:str,value:float)-None:rownext(aforainANCHORSifa[0]name)_,lo,hi,condrow insidelovaluehi flag在范围内ifinsideelse在范围外print(f{name:22s}计算值{value:8.4g}锚点 [{lo:g},{hi:g}]{flag}{cond})defmain()-None:# 下面这些“计算值”请替换为你自己模型的输出# 这里填的是代码 2 教学模型与代码 1 的示意工况结果。check(贫液进塔温度, ℃,49.0)check(L/G, kg/kg,1.023)check(富液负载, mol/mol,0.489)check(贫液负载, mol/mol,0.23)check(捕集率经典基准,0.90)print(\n提醒L/G 1.023 kg/kg 是按 30 wt% MEA 示意工况算出的值)print( 与 Sherman FEED 的 1.069 kg/kg 相差约 4%属工况差异不是模型错误)print( 富液负载 0.489 高于锚点上限 0.476原因是代码 2 的示意平衡关系不是真实数据)print( 真实值必须由 Aspen Plus 的电解质物性包算出。)if__name____main__:main()为什么单独写一个校核脚本工程上最容易出问题的不是算不出来而是算出来了但没人核对。把锚点固化成数据结构而不是散落在文档里每次模型输出都过一遍是最便宜的防错手段。三、常见报错与排查现象 1win32com.client.gencache.EnsureDispatch(Apwn.Document)报com_error。根因本机未安装 Aspen Plus 桌面版或 ProgID 未注册。Apwn.Document属第三方可运行实现 官方示例佐证的写法需在本机验证。解法确认桌面版已安装仍失败时改用后期绑定的win32com.client.Dispatch再试带版本后缀的 ProgID 未经官方确认不要直接写进代码。现象 2脚本改了TOTFLOW但结果完全没变。根因改输入后没有重置或重跑。解法按Node.Value xxx→Application.Reinit()→Application.Engine.Run2()的顺序执行Reinit也有Engine.Reinit(IAP_REINIT_SIMULATION)这类带参数的写法官方 VB 示例中出现可按需选择粒度。现象 3脚本无输出地卡死。根因收敛失败弹出对话框。解法进入脚本第一件事就是aspen.SuppressDialogs True同时用aspen.Visible False提速。现象 4升级 Aspen 后FindNode返回 None脚本全线报错。根因树路径硬编码版本迁移会改变节点结构AspenTech 官方把版本迁移当作需要官方介入的常规工程动作。解法把路径集中成常量表并给FindNode加 None 检查代码 1 的get_node升级后先跑一次路径探测再跑业务脚本。现象 5RadFrac 收敛了但出口 CO₂ 与手算或文献差很多。根因常见三类——物性方法选错例如用了不支持离子活度的模型、反应集未通过元素与电荷衡算、把填料的理论级数直接当成几米填料来理解。解法依次核对物性方法回到第 03 篇、反应集自洽性第 03 篇代码 2、以及理论级 → HETP → 塔高的换算关系本篇代码 2。四、动手练习练习 1锚点校核运行代码 3把 L/G 的计算值改成 1.069、富液负载改成 0.50再跑一次。成功判定L/G 应显示在范围内富液负载应显示在范围外0.50 0.476 上限。请写下结论“锚点只用于判断量级是否合理落在区间外时必须先解释工况差异。”练习 2理论级与塔高在代码 2 中把R_MEA_GAS从 0.15 改为 0.10相当于降低溶剂循环量重新运行。成功判定理论级数应从 3 增加到 5富液负载升到约 0.57等价的 HETP 从 5.0/8.0 m 降到约 3.0/4.8 m。请写下结论“减小溶剂循环量会让吸收塔需要更多理论级也就是更高或更多填料。”练习 3COM 路径探测在本机有 Aspen Plus 桌面版的前提下把代码 1 中的APW_PATH换成你自己的文件先用Tree.FindNode读取\Data\Streams\FEED\Input\TOTFLOW的当前值并打印。成功判定能打印出原值、能改写为新值、并在Reinit()Run2()后COMPSTATUS显示可收敛。若FindNode返回 None说明你的树路径与示例不同请回到界面上核对真实路径此步骤必须在本机完成。五、小结与下一篇预告本篇跑通了第一个平衡级 MEA 吸收塔组分覆盖 MEA、H₂O、CO₂ 与 MEACOO⁻、MEAH⁺、HCO₃⁻ 等离子物种物性走 eNRTL 路径用 RadFrac 把 15–24 m 填料抽象成理论级脚本层用Apwn.Document、Tree.FindNode、Node.Value、Reinit、Engine.Run2完成改参数—重置—运行—取数闭环并用贫液 49 ℃ 与 L/G 1.069 kg/kg 两个锚点做量级校核。入门四篇到此结束。从下一篇开始进入核心机制第 05 篇讲速率基与双膜理论解释为什么平衡级把传质阻力藏起来了第 08 篇把本篇搭好的 RadFrac 塔改造成速率基塔第 11 篇再把本篇的 COM 脚本扩展成带重试、超时与结果落盘的批量自动化框架。本篇认知问题回显FAQQ1用 Aspen Plus 模拟 MEA 吸收塔时组分里必须加哪些离子物种A除 MEA、H₂O、CO₂ 三个分子外还需加入 MEACOO⁻、MEAH⁺、HCO₃⁻、CO₃²⁻、H⁺、OH⁻ 等离子物种才能表达反应吸收。Q2MEA 吸收塔的物性方法为什么选 eNRTLA因为 eNRTL 可计算全部离子与分子物种的活度系数能反映高负载下平衡分压被抬高的效应其规范引用为 Chen Evans (1986)。Q3RadFrac 里的理论级与填料高度是什么关系ARadFrac 是平衡级模型用理论级数代替塔高把 15 m 或 24 m 填料折算到 3 个理论级等价 HETP 约为 5 m 与 8 m。Q4用 Python 驱动 Aspen Plus 时Run2 应该怎么写AApplication.Run2() 与 Application.Engine.Run2() 两种写法都见于真实案例需在本机验证不能断言唯一正确写法。Q5MEA 吸收塔仿真结果该拿什么做校核A可用贫液温度 49 ℃、液气比 1.069 kg/kg富液 0.275–0.476、贫液 0.075–0.239 mol CO₂/mol MEA 作量级参照绑定 35 wt% MEA 工况。
燃烧后CO2捕集MEA吸收过程模拟教程04第一个可跑的仿真——用 Aspen Plus 搭一个平衡级 MEA 吸收塔版本声明块工具/软件Aspen Plus官方产品页未公示数字版本号以官方文档为准RadFrac 为 Aspen Plus 最常用的平衡级多级分离模型速率基在其上扩展语言/环境Python 3.10代码 1 需 Windows pywin32 本机已安装 Aspen Plus 桌面版代码 2、3 纯标准库任何平台均可运行本文目标跑通本系列第一个平衡级 MEA 吸收塔并建立结果必须用已核实锚点校核的习惯一句话结论第一个可跑的 MEA 吸收塔应当由 RadFrac 平衡级模型给出初值脚本层用Apwn.Document打开 .apw、用Tree.FindNode(r\Data\Streams\FEED\Input\TOTFLOW)改流量、Reinit后运行并把贫液进塔温度 49 ℃ 与液气比 1.069 kg/kg 当作量级校核锚点而不是精确目标。〇、本篇要解决的认知问题MEA-H₂O-CO₂ 体系在 Aspen Plus 里要定义哪些组分与离子物种为什么物性方法要走 eNRTL 路径它和第 03 篇的分水岭结论怎么衔接怎么用 RadFrac 搭一个平衡级吸收塔把 15–24 m 填料高度抽象成理论级COM 自动化里Tree.FindNode、Node.Value、Reinit、Run2具体怎么用哪些写法不能被断言为唯一正确结果该拿什么校核为什么贫液 49 ℃ 与 L/G 1.069 kg/kg 只能用作量级参照一、机制解析1.1 建模路线总览组分定义 物性方法 单元模型 求解与校核 MEA / H2O / CO2 → eNRTL 电解质路径 → RadFrac 平衡级吸收塔 → 收敛 → 与 49 ℃ / 1.069 kg/kg 对照 MEACOO- / HCO3- 填料高度 → 理论级 ↓ 等离子物种 再切速率基第 08 篇这五步的顺序不能颠倒物种没定义全物性方法就无法建立物性方法没定RadFrac 的收敛结果没有意义没有校核锚点收敛了也不知道对不对。1.2 组分与电解质物种底账 [C5] 明确 MEA 水溶液吸收 CO₂ 时同时生成 MEACOO⁻ 与 HCO₃⁻并给出平衡反应式 2–6。因此组分定义至少要覆盖类别物种分子中性MEA、H₂O、CO₂离子MEACOO⁻、MEAH⁺、HCO₃⁻、CO₃²⁻、H⁺、OH⁻注意两点第一不要只定义 MEA / H₂O / CO₂ 三个分子就跑——没有离子物种MEA 与 CO₂ 反应生成氨基甲酸盐这件事在模型里就根本没有表达第二反应式必须先通过元素与电荷衡算第 03 篇代码 2 做的正是这件事否则后续的平衡计算会给出看似收敛、实则无意义的解。关于Aspen Plus 的电解质向导它的具体点击步骤与参数生成流程不在本系列可核实的资料范围内底账附录 A 已把电解质向导生成 eNRTL 参数的官方步骤文档列为待补方向本系列不做逐步描述请以官方文档为准。1.3 物性方法走 eNRTL 路径第 03 篇的结论在这里落地需要全物种活度系数时选 eNRTL规范引用为 Chen Evans (1986) [C3]。Aspen Plus 内置电解质物性包含 eNRTL 及 MEA 专用物性 insert [C1][C2][S4]。三点提示平衡常数在 Aspen 中采用ln K(T) a b/T c·lnT d·T形式 [C3]四个系数都有回归温度窗口吸收侧与再生侧不要共用一套外推如果后续要切速率基塔必须同时具备黏度、表面张力、导热系数、二元扩散系数等传递物性 [C1]这也是本系列建议一开始就选完整电解质物性包的原因CCSI 的MEA Steady State Model是 Aspen Plus 的附加模型含物性、水力学、界面面积、动力学四类子模型不是Aspen 内置的标准物性包 [C14]写文档时不要混淆。1.4 用 RadFrac 表达吸收塔填料高度 → 理论级RadFrac 是 Aspen Plus 最常用的平衡级多级分离模型速率基正是在其上扩展而来 [S3]。所以第一步先用平衡级搭塔是官方认可的路径AspenTech 官方课程明确包含把 RadFrac 平衡级 CO₂ 吸收塔转换为使用 Aspen Rate-Based Distillation 的塔[S15]。平衡级塔里塔高这一物理信息被压缩成理论级数。真实装置的填料高度是 15 mSherman FEED两塔并联与 24 mTCM三床 1266 m[C14][C19]。这里有一个非常有用的量级检验若把 15–24 m 填料折算到少数几个理论级上等价 HETP 会高达 5–8 m明显高于常规精馏填料塔的常见量级此为基于 [C14][C19] 数值的量级判断。这个反差本身就是平衡级抽象在 MEA 吸收塔上很粗糙的直接证据也是第 05、08 篇要走向速率基的动机。1.5 COM 对象模型真实名称与必须加的限定元素真实名称使用要点COM 对象 ProgIDApwn.Documentwin32com.client.gencache.EnsureDispatch(Apwn.Document)需在本机验证带版本后缀的 ProgID 未获官方确认打开文档InitFromArchive2(path)从 .apw / .bkp 启动节点树Application.TreeTree.FindNode(path)按路径取节点如\Data\Streams\FEED\Input\TOTFLOW遍历子元素Tree.Elements(Data).Elements(Blocks)/.Elements(name)用于探路避免硬编码输出路径读写数值Node.Value如\Data\Blocks\B1\Input\BASIS_RR运行Application.Engine.Run2()等价运行按钮Application.Run2()亦见于案例重置Application.Reinit()/Engine.Reinit(IAP_REINIT_SIMULATION)改动输入后使用官方 VB 示例给出go_Simulation.Engine.Reinit(IAP_REINIT_SIMULATION)另有IAP_REINIT_BLOCK、IAP_REINIT_CONVERGENCE停止Application.Stop()等价停止状态检查Application.COMPSTATUS/Node.COMPSTATUS输入完整性与可运行性检查弹窗抑制Application.SuppressDialogs True/False批量脚本必须设置可见性Application.Visible True/False关界面可提速导出Application.Export(code, filename)代码 1.bkp、2.rep、3.sum、4.inp、6运行消息、9/10动态报告、11DXF该组代码取自第三方 API 封装以官方 User Guide 第 38 章为准官方文档《Aspen Plus User Guide》第 38 章 Happ ActiveX interfaceAspenTech 支持站需账号关于Run2的写法真实案例中Application.Run2()与Application.Engine.Run2()都被使用过 [S17][S20]因此本系列不断言唯一正确写法写法上采用先试一种、异常时回退另一种并把最终结果交回本机验证。关于表格中方法的来源与时效Apwn.Document、InitFromArchive2、FindNode、Run2、导出代码这一组写法来自第三方可运行示例仓库随作者更新可能失效并有官方知识库示例佐证 [S17][S18][S19][S20]。官方完整的 ActiveX 方法清单在《Aspen Plus User Guide》第 38 章属于需要账号的支持文档 [S18][S21]因此本篇所有调用都标注为需在本机验证不得当作官方 SDK 的固定签名使用。关于工程坑COM 调用依赖本机安装的桌面版 Aspen Plusgencache.EnsureDispatch会生成类型库包装无桌面安装的服务器/容器环境不可用不设SuppressDialogs True时收敛失败弹窗会挂死脚本 [S17]。这些是第 11 篇批量自动化要系统解决的问题。1.6 校核锚点已核实、可作量级参照的锚点是吸收塔进液贫液温度49 ℃Sherman FEED35 wt% MEA[C14]液气比 L/G1.069 kg/kgSherman FEED[C14]富液负载 0.275–0.476、贫液负载 0.075–0.239 mol CO₂/mol MEA35 wt% MEA、CCGT 烟气[C14]。为什么只能当量级参照这三个值绑定的是 35 wt% MEA、特定烟气与特定装置配置Sherman FEED。如果你的模型用的是 30 wt% MEA 或不同烟气组成算出的 L/G 与 1.069 有差异是正常的差异必须回到工况层面解释。把锚点当必须命中的目标值去反向调参是这一篇最危险的操作。二、完整代码与逐行剖析代码 1用 Python COM 驱动 Aspen Plus需本机验证# 04_code1_aspen_com_absorber.py# 运行环境Windows Python 3.10 pywin32 本机已安装 Aspen Plus 桌面版# 前置条件准备一个已建好的 .apw其中含一台 RadFrac 吸收塔与一条名为 FEED 的进料流股# 【重要】Apwn.Document 与 Run2 / Reinit 的写法来自第三方可运行实现与官方示例佐证# 官方完整 ActiveX 清单在《Aspen Plus User Guide》第 38 章需账号。# 因此本脚本的每一处调用都需在本机验证不能断言唯一正确写法。importosimportwin32com.clientaswin32# 需 pip install pywin32APW_PATHrC:\path\to\mea_absorber.apw# 换成你自己的文件路径FEED_TOTAL_FLOW10000.0# 烟气总摩尔流率, kmol/h示意工况defget_node(aspen,path):按路径取节点取不到就抛出可读错误。 为什么先判 None树路径是硬编码的版本迁移后路径可能消失 直接 .Value 会抛出不带上下文的 COM 异常排查成本极高。nodeaspen.Tree.FindNode(path)ifnodeisNone:raiseRuntimeError(f节点不存在{path}升级后路径可能变化先做路径探测)returnnodedefmain():aspenwin32.gencache.EnsureDispatch(Apwn.Document)# 依赖本机注册的 COM 服务器try:aspen.InitFromArchive2(APW_PATH)# 从 .apw / .bkp 启动仿真aspen.VisibleFalse# 关界面批量跑时可提速aspen.SuppressDialogsTrue# 必须收敛失败弹窗会阻塞甚至挂死脚本# 输入完整性检查。不同接口封装下该属性的可读性不同读不到就跳过# 不让一个可选的诊断信息阻断整条流程。try:print(输入检查 COMPSTATUS ,aspen.COMPSTATUS)exceptExceptionasexc:print(COMPSTATUS 不可读以所用版本为准,exc)# 修改一个“有出处”的输入节点进料总流率flow_nodeget_node(aspen,r\Data\Streams\FEED\Input\TOTFLOW)flow_node.ValueFEED_TOTAL_FLOWprint(f已把 FEED 的 TOTFLOW 设为{flow_node.Value}kmol/h)aspen.Reinit()# 重置让输入改动生效# 运行Application.Run2() 与 Application.Engine.Run2() 两种写法都见诸案例# 这里先试 Engine.Run2()异常时回退到 Run2()最终以本机验证结果为准。try:aspen.Engine.Run2()exceptException:aspen.Run2()try:print(运行后 COMPSTATUS ,aspen.COMPSTATUS)exceptException:pass# 读结果不硬编码输出节点名而是遍历该流股 Output 下的子节点。# 为什么遍历Output 下的节点名随模型与版本变化硬编码是脚本升级后失效的头号原因。try:outaspen.Tree.FindNode(r\Data\Streams\RICH\Output)forelinout.Elements:print(f RICH.{el.Name}{el.Value})exceptExceptionasexc:print(遍历输出节点失败不同封装/版本遍历方式不同需在本机验证,exc)finally:# 无论成功失败都导出报告便于事后核对导出代码含义以官方 User Guide 第 38 章为准try:reportos.path.join(os.path.dirname(APW_PATH),mea_absorber.rep)aspen.Export(2,report)print(已导出报告,report)exceptExceptionasexc:print(导出报告失败,exc)if__name____main__:main()反直觉默认值SuppressDialogs的默认值是不抑制弹窗。这在交互界面上没问题但在脚本里一个是否继续的对话框就意味着无限等待。批量工况脚本必须把它显式设为True。代码 2不依赖 Aspen 的平衡级吸收塔教学模型# 04_code2_equilibrium_absorber.py# 运行环境Python 3.10仅标准库任何平台可运行# 用途用“平衡级 操作线”逐步追级法McCabe-Thiele 式估算 MEA 吸收塔的理论级数# 并与 15–24 m 真实填料高度对照理解“理论级”抽象在 MEA 体系里有多粗糙。# 【重要】平衡关系 y_eq C_ALPHA * x^2 / (X_MAX - x) 是【示意参数】形式# 不是 MEA-H2O-CO2 的真实平衡数据真实平衡分压须由第 03 篇的电解质模型给出。# 本脚本只用于理解“理论级 → 塔高”的换算逻辑。X_MAX1.0# 示意负载上限, mol CO2/mol MEAC_ALPHA0.041# 示意平衡关系系数R_MEA_GAS0.15# 示意MEA 摩尔流率 / 烟气总摩尔流率Y_IN0.035# 入口烟气 CO2 摩尔分数底账烟气量级 3.5-15% 的低端, [C14][C19]Y_OUT0.0035# 出口烟气 CO2 摩尔分数对应 90% 捕集底账经典基准 [C19][C14]X_LEAN0.23# 贫液负载, mol CO2/mol MEA底账实测 0.075-0.239 区间内, [C14]defy_eq(x:float)-float:气相 CO2 摩尔分数与液相负载的示意平衡关系单调递增x 趋近 X_MAX 时急剧上升。returnC_ALPHA*x**2/(X_MAX-x)defx_from_y(y:float,lo:float0.0,hi:floatX_MAX-1e-6,tol:float1e-10)-float:对 y_eq 做二分反解。 为什么要反解而不是解析解真实电解质模型的平衡关系没有解析式 工程代码一律用数值反解这里提前用最简单的函数把流程走通。for_inrange(200):mid0.5*(lohi)ify_eq(mid)y:lomidelse:himidifhi-lotol:breakreturn0.5*(lohi)defstages_needed(y_in:float,y_out:float,x_in:float)-tuple:从塔顶往下逐步追级直到进入塔底的气相浓度达到 y_in。 返回 (理论级数, 富液负载)。 操作线由全塔物料衡算给出y y_out R_MEA_GAS * (x - x_in) 其中斜率就是 MEA 摩尔流率与烟气摩尔流率之比——这一点常被写错成倒数。xx_from_y(y_out)# 第 1 级出口气与离开该级的液相平衡n1richxwhileTrue:y_nexty_outR_MEA_GAS*(x-x_in)# 进入下一级更靠下的气相组成ify_nexty_in:returnn,rich n1ifn100:# 防死循环超过 100 级说明工况不可行return-1,rich xx_from_y(y_next)richxdefmain()-float:n,richstages_needed(Y_IN,Y_OUT,X_LEAN)print( 平衡级吸收塔逐步追级示意平衡关系 )print(f入口 y_CO2 {Y_IN:.4f}出口 y_CO2 {Y_OUT:.4f}捕集率{(1-Y_OUT/Y_IN):.1%})print(f贫液负载 x_in {X_LEAN:.3f}mol CO2/mol MEA示意)print(f理论级数 N {n}富液负载 x_out {rich:.3f}mol CO2/mol MEA)# 用底账已核实的填料高度换算等价的 HETP作为“平衡级抽象有多粗糙”的量级检验forheightin(15.0,24.0):print(f若真实填料高度{height:.0f}m[C14][C19]由{n}个理论级承担f等价 HETP ≈{height/n:.1f}m/理论级)returnrichif__name____main__:main()预期输出与判读在示意参数下90% 捕集会给出3 个理论级与约0.49的富液负载。0.49 略高于底账实测的富液上限 0.476 mol CO₂/mol MEA [C14]——这不是程序算错而是提醒你示意平衡关系不是真实数据只能用来理解换算逻辑不能当设计结果。把 15 m 与 24 m 折算到 3 个理论级等价 HETP 约 5.0 与 8.0 m/理论级远高于常规精馏填料塔的量级。这不是模型错了而是**理论级这个抽象把反应吸收的传质阻力全部藏起来了**。代码 3结果校核脚本锚点对照不下结论# 04_code3_cross_check.py# 运行环境Python 3.10仅标准库# 用途把仿真或教学模型算出的关键量与底账已核实的锚点并排对照输出“在区间内/外”。# 【重要】锚点绑定的是 35 wt% MEA 与 Sherman FEED / CCGT 工况[C14]# 只作量级参照落在区间外不等于模型错必须回到工况差异解释禁止反向调参凑数。# 锚点全部来自底账 §4逐条带条件ANCHORS[(贫液进塔温度, ℃,49.0,49.0,Sherman FEED, 35 wt% MEA [C14]),(L/G, kg/kg,1.069,1.069,Sherman FEED [C14]),(富液负载, mol/mol,0.275,0.476,35 wt% MEA, CCGT 烟气 [C14]),(贫液负载, mol/mol,0.075,0.239,35 wt% MEA, CCGT 烟气 [C14]),(捕集率经典基准,0.90,0.90,90% 为经典设计基准 [C19][C14]),]defcheck(name:str,value:float)-None:rownext(aforainANCHORSifa[0]name)_,lo,hi,condrow insidelovaluehi flag在范围内ifinsideelse在范围外print(f{name:22s}计算值{value:8.4g}锚点 [{lo:g},{hi:g}]{flag}{cond})defmain()-None:# 下面这些“计算值”请替换为你自己模型的输出# 这里填的是代码 2 教学模型与代码 1 的示意工况结果。check(贫液进塔温度, ℃,49.0)check(L/G, kg/kg,1.023)check(富液负载, mol/mol,0.489)check(贫液负载, mol/mol,0.23)check(捕集率经典基准,0.90)print(\n提醒L/G 1.023 kg/kg 是按 30 wt% MEA 示意工况算出的值)print( 与 Sherman FEED 的 1.069 kg/kg 相差约 4%属工况差异不是模型错误)print( 富液负载 0.489 高于锚点上限 0.476原因是代码 2 的示意平衡关系不是真实数据)print( 真实值必须由 Aspen Plus 的电解质物性包算出。)if__name____main__:main()为什么单独写一个校核脚本工程上最容易出问题的不是算不出来而是算出来了但没人核对。把锚点固化成数据结构而不是散落在文档里每次模型输出都过一遍是最便宜的防错手段。三、常见报错与排查现象 1win32com.client.gencache.EnsureDispatch(Apwn.Document)报com_error。根因本机未安装 Aspen Plus 桌面版或 ProgID 未注册。Apwn.Document属第三方可运行实现 官方示例佐证的写法需在本机验证。解法确认桌面版已安装仍失败时改用后期绑定的win32com.client.Dispatch再试带版本后缀的 ProgID 未经官方确认不要直接写进代码。现象 2脚本改了TOTFLOW但结果完全没变。根因改输入后没有重置或重跑。解法按Node.Value xxx→Application.Reinit()→Application.Engine.Run2()的顺序执行Reinit也有Engine.Reinit(IAP_REINIT_SIMULATION)这类带参数的写法官方 VB 示例中出现可按需选择粒度。现象 3脚本无输出地卡死。根因收敛失败弹出对话框。解法进入脚本第一件事就是aspen.SuppressDialogs True同时用aspen.Visible False提速。现象 4升级 Aspen 后FindNode返回 None脚本全线报错。根因树路径硬编码版本迁移会改变节点结构AspenTech 官方把版本迁移当作需要官方介入的常规工程动作。解法把路径集中成常量表并给FindNode加 None 检查代码 1 的get_node升级后先跑一次路径探测再跑业务脚本。现象 5RadFrac 收敛了但出口 CO₂ 与手算或文献差很多。根因常见三类——物性方法选错例如用了不支持离子活度的模型、反应集未通过元素与电荷衡算、把填料的理论级数直接当成几米填料来理解。解法依次核对物性方法回到第 03 篇、反应集自洽性第 03 篇代码 2、以及理论级 → HETP → 塔高的换算关系本篇代码 2。四、动手练习练习 1锚点校核运行代码 3把 L/G 的计算值改成 1.069、富液负载改成 0.50再跑一次。成功判定L/G 应显示在范围内富液负载应显示在范围外0.50 0.476 上限。请写下结论“锚点只用于判断量级是否合理落在区间外时必须先解释工况差异。”练习 2理论级与塔高在代码 2 中把R_MEA_GAS从 0.15 改为 0.10相当于降低溶剂循环量重新运行。成功判定理论级数应从 3 增加到 5富液负载升到约 0.57等价的 HETP 从 5.0/8.0 m 降到约 3.0/4.8 m。请写下结论“减小溶剂循环量会让吸收塔需要更多理论级也就是更高或更多填料。”练习 3COM 路径探测在本机有 Aspen Plus 桌面版的前提下把代码 1 中的APW_PATH换成你自己的文件先用Tree.FindNode读取\Data\Streams\FEED\Input\TOTFLOW的当前值并打印。成功判定能打印出原值、能改写为新值、并在Reinit()Run2()后COMPSTATUS显示可收敛。若FindNode返回 None说明你的树路径与示例不同请回到界面上核对真实路径此步骤必须在本机完成。五、小结与下一篇预告本篇跑通了第一个平衡级 MEA 吸收塔组分覆盖 MEA、H₂O、CO₂ 与 MEACOO⁻、MEAH⁺、HCO₃⁻ 等离子物种物性走 eNRTL 路径用 RadFrac 把 15–24 m 填料抽象成理论级脚本层用Apwn.Document、Tree.FindNode、Node.Value、Reinit、Engine.Run2完成改参数—重置—运行—取数闭环并用贫液 49 ℃ 与 L/G 1.069 kg/kg 两个锚点做量级校核。入门四篇到此结束。从下一篇开始进入核心机制第 05 篇讲速率基与双膜理论解释为什么平衡级把传质阻力藏起来了第 08 篇把本篇搭好的 RadFrac 塔改造成速率基塔第 11 篇再把本篇的 COM 脚本扩展成带重试、超时与结果落盘的批量自动化框架。本篇认知问题回显FAQQ1用 Aspen Plus 模拟 MEA 吸收塔时组分里必须加哪些离子物种A除 MEA、H₂O、CO₂ 三个分子外还需加入 MEACOO⁻、MEAH⁺、HCO₃⁻、CO₃²⁻、H⁺、OH⁻ 等离子物种才能表达反应吸收。Q2MEA 吸收塔的物性方法为什么选 eNRTLA因为 eNRTL 可计算全部离子与分子物种的活度系数能反映高负载下平衡分压被抬高的效应其规范引用为 Chen Evans (1986)。Q3RadFrac 里的理论级与填料高度是什么关系ARadFrac 是平衡级模型用理论级数代替塔高把 15 m 或 24 m 填料折算到 3 个理论级等价 HETP 约为 5 m 与 8 m。Q4用 Python 驱动 Aspen Plus 时Run2 应该怎么写AApplication.Run2() 与 Application.Engine.Run2() 两种写法都见于真实案例需在本机验证不能断言唯一正确写法。Q5MEA 吸收塔仿真结果该拿什么做校核A可用贫液温度 49 ℃、液气比 1.069 kg/kg富液 0.275–0.476、贫液 0.075–0.239 mol CO₂/mol MEA 作量级参照绑定 35 wt% MEA 工况。
