2. Aspen Plus 流程模拟回顾

设计规定(Design Spec)和灵敏度分析(Sensitivity)之所以常被初学者用错,根源往往不在工具本身,而在前面的流程模拟没有"搭稳"。设计规定本质上是让 Aspen Plus 反算某个输入变量,使指定的输出达到目标值;灵敏度分析则是扫描一个或多个输入,观察输出如何变化。两者都建立在同一个前提上:物性方法选对、单元模型配置正确、收敛真实可靠、结果能被准确读取。本节把这条基础链梳理清楚,后续所有设计规定与灵敏度分析的操作才有意义。

物性方法:结果可信度的第一道关口

物性方法决定了气液平衡、焓、密度、黏度等性质的计算方式,选错方法,后面的收敛和优化都是在错误基础上做文章。选择的基本逻辑是:先判断体系是极性还是非极性、是否含电解质、是否在高压下、是否有强非理想性。常见对应关系如下:

体系特征推荐物性方法典型场景
烃类、轻气体,中低压SRK / PR(状态方程法)天然气处理、炼油分离
极性非电解质,低压NRTL / UNIQUAC / Wilson醇-水、酮-水精馏
强非理想、共沸体系NRTL + 二元交互参数共沸精馏、萃取精馏
电解质水溶液ELECNRTL酸碱中和、盐水浓缩
高压含氢体系PR / SRK(含 Henry 组分)加氢、合成气

实操中要特别注意两点:一是二元交互参数是否齐全,NRTL 等活度系数模型若缺少交互参数,软件会用估算值,结果可能明显偏离;二是Henry 组分设置,当体系中存在超临界轻气体(如 H₂、N₂、CO₂)时,必须在 Components 中将其设为 Henry 组分,否则气相溶解度计算会失真。物性方法一旦确定,后续设计规定中"以某塔顶产品纯度为目标反算回流比"这类任务才有物理意义。

单元操作模型:把工艺意图翻译成正确的模块

流程中的每个设备都要选对模型。精馏塔常用 RadFrac(严格平衡级),若只做简捷估算可用 DSTWU 或 Shortcut;换热器用 Heater 或 HeatX;反应器按动力学或平衡选择 RStoic、REquil、RGibbs 或 RPlug;泵、压缩机用 Pump、Compr。模型选错,自由度就错,设计规定也就无从谈起。

以 RadFrac 为例,其自由度由"理论板数、进料位置、回流比、馏出量、侧线采出"等共同决定。做设计规定时,通常固定理论板数和进料位置,让回流比或馏出量作为被调变量。这里有一个容易被忽略的细节:RadFrac 的收敛算法(如 Standard、Petroleum、Strongly Non-ideal)要与体系匹配,强非理想体系用 Standard 容易震荡,应改用 Newton 或 Non-ideal 算法。模型配置不当,灵敏度分析会出现"扫描到某一点突然不收敛"的假象。

收敛方法:让结果"真收敛"而不是"看起来收敛"

Aspen Plus 的收敛包括两层:单元模块内部收敛(如 RadFrac 的塔内迭代)和全流程撕裂流收敛(Tear Stream)。全流程默认用 Wegstein 加速,遇到循环回路多、强非理想或高度耦合的流程时,可切换为 Broyden 或 Newton。判断是否真收敛,不能只看"Convergence 显示 Converged",还要检查:

  • 残差是否降到容差以下(默认 1e-4 量级),而不是靠放宽容差"通过";
  • 关键变量(如塔顶温度、回流比、产品纯度)在最后几次迭代中是否稳定;
  • 有无"Convergence block not converged"警告被忽略;
  • 物料平衡与能量平衡是否闭合(可在 Report 中查看)。
做设计规定和灵敏度分析前,务必先让基础流程在不放宽容差的条件下干净收敛。否则设计规定会把"不收敛"误判为"无解",灵敏度曲线也会出现断裂或跳变。

报告与结果查看:读懂数据再动手优化

结果查看的核心是 Report。常用入口包括:Streams 的物流结果(组成、流量、温度、压力、焓)、Blocks 的模块结果(塔板温度/组成剖面、换热器热负荷、反应器转化率)、以及 Convergence 的收敛摘要。做设计规定前,先确认目标变量在哪个报告位置、单位是什么、是摩尔分率还是质量分率,避免目标值填错量纲。

下面这张图把"物性方法—单元模型—收敛—结果"这条基础链与设计规定、灵敏度分析的衔接关系画在一起,帮助建立整体认识:

流程模拟基础链 → 设计规定 / 灵敏度分析 ① 物性方法 SRK / PR / NRTL Henry 组分 · 交互参数 ② 单元操作模型 RadFrac / Heater / RStoic 自由度 · 理论板 · 进料位置 ③ 收敛方法 Wegstein / Broyden 残差 · 撕裂流 · 容差 ④ 报告与结果 Streams / Blocks 单位 · 量纲 · 平衡闭合 不收敛 → 回查模型与物性 基础流程干净收敛 设计规定 Design Spec 目标:塔顶纯度 = 0.995 被调变量:回流比 约束:馏出量固定 灵敏度分析 Sensitivity 扫描:回流比 1.2 → 3.0 观察:塔顶纯度、再沸负荷 步长:0.2,共 10 点 互为验证 设计规定给出"达到目标所需的操作条件",灵敏度分析给出"操作条件变化对目标的影响趋势" 两者都依赖:正确的物性方法 + 正确的模型 + 真实收敛 + 准确读取结果

一个可复现的示例:为设计规定准备基础流程

假设用 RadFrac 模拟一座脱丙烷塔,进料为丙烷/丁烷混合物,目标塔顶丙烷摩尔分率 0.995。基础配置可参考:

物性方法: PR(烃类中低压,含少量轻组分)
塔模型: RadFrac,理论板 20,进料板 10
操作: 冷凝器 Total,回流比初值 1.5,馏出量 D = 100 kmol/h
收敛: Standard 算法,容差 1e-4

先运行基础流程,确认收敛摘要无警告、塔顶温度与组成稳定、物料平衡闭合。此时若直接进入设计规定,把"塔顶丙烷摩尔分率"设为目标 0.995、被调变量设为回流比,软件会迭代反算。若基础流程本身不收敛或物性方法不当,反算可能报"Design spec not converged",这并不代表工艺无解,而是基础没打牢。

灵敏度分析同理:以回流比 1.2~3.0、步长 0.2 扫描,观察塔顶纯度与再沸器热负荷的变化曲线。若曲线在某点突然断裂,应回查该点附近是否出现塔内迭代不收敛,而不是简单认为"数据异常"。把设计规定的反算结果与灵敏度曲线上的对应点对照,两者应基本吻合,这就是对基础流程正确性的交叉验证。

建议养成习惯:每次做设计规定或灵敏度分析前,先保存一份"基础流程已收敛"的备份文件(.bkp),并在 Report 中记录关键物流的流量、组成、温度、压力。这样后续任何反算或扫描出现异常,都能快速判断是基础流程问题还是分析设置问题。

把物性方法、单元模型、收敛和结果查看这四件事做扎实,设计规定与灵敏度分析才是在"可信的模拟"上做"有意义的优化",而不是在不确定的结果上反复试错。