在化工过程模拟领域Aspen Plus 和 Aspen HYSYS 是工程师进行流程设计和优化的核心工具。实际工业流程中物料循环和能量回收是提升经济性的关键手段但带有循环的流程模拟也是新手最容易遇到收敛问题的场景。本文将以一个典型的带有循环的工艺流程为例详细讲解在 Aspen Plus 中建立模型、设置撕裂流股、选择收敛方法以及分析结果的完整步骤。1. 理解流程循环与收敛的基本原理1.1 为什么化工流程中会存在循环化工生产追求原料的充分利用和能量的高效回收。例如在反应系统中未反应的原料需要从产物中分离并返回反应器入口在精馏系统中塔顶冷凝器的一部分液体作为回流返回塔顶塔釜的一部分液体也可能作为再沸器的进料。这些都将流程的输出端与输入端连接起来形成了循环。从数学角度看一个带有循环的流程实际上是在求解一个非线性方程组。循环流股通常称为“撕裂流股”或“循环流股”的变量温度、压力、流量、组成既是方程组的输入也是输出。模拟软件需要不断迭代直到输入的猜测值与计算出的输出值在允许的误差范围内相等。1.2 Aspen Plus 求解循环流程的核心机制Aspen Plus 采用序贯模块法进行模拟。这意味着计算必须按照模块之间的输入输出关系有顺序地进行。当流程中存在循环时这个顺序就被打破了因为没有一个模块可以作为计算的起点。为了解决这个问题Aspen Plus 引入了“撕裂流股”的概念。软件会选择一个或多个流股作为撕裂流股并为其变量提供初始猜测值。然后模拟沿着一个预设的计算顺序进行一轮计算得到撕裂流股的新值。将新值与初始猜测值进行比较如果差异大于收敛容差则用新值或某种插值方法得到的新猜测值作为下一轮迭代的初始值重复此过程直到收敛。整个过程的核心是“收敛模块”它负责管理迭代过程。用户的关键任务在于合理选择撕裂流股。为撕裂流股提供尽可能接近真实解的初始值。选择合适的收敛算法。2. 案例流程描述与模拟目标本案例模拟一个简化的带有循环的化学反应精馏系统。流程描述如下一股新鲜的原料流FEED包含反应物 A 和少量惰性组分 I。FEED与一股循环流RECYCLE混合形成反应器进料流REAC-FEED。REAC-FEED进入一个平衡反应器REACTOR发生反应 A - B。反应未完全转化率为 70%。反应器出口流股REAC-OUT进入一个精馏塔COLUMN。精馏塔塔顶采出主要含未反应物 A 和惰性组分 I 的流股DIST-OUT该流股一部分作为产品排出PRODUCT另一部分作为循环流RECYCLE返回与新鲜原料混合。精馏塔塔釜采出为主要产物 B 的流股BOTTOM。模拟目标建立并收敛该带有物料循环的流程。确定在给定新鲜进料和塔顶采出比的情况下循环流股RECYCLE的流量和组成。分析系统的物料平衡。3. 在 Aspen Plus V12 中构建流程模型3.1 创建新模拟与物性方法选择启动 Aspen Plus选择“模板”新建一个模拟。在“化学”类别下选择“通用 with Metric Units”模板。这将设置一套常用的公制单位。进入模拟环境后首先需要定义组分和物性方法。点击导航栏的ComponentsSpecificationsSelection。在组分列表中输入本案例所需的组分A例如 Benzene、B例如 Toluene和 I例如 n-Hexane。您可以使用真实组分或自定义的虚拟组分。点击MethodsSpecificationsGlobal选择物性方法。对于本例中的烃类混合物选择“PENG-ROB”Peng-Robinson 状态方程是合适的。完成物性分析MethodsReviewResults查看是否有警告或错误。3.2 搭建流程图切换到Simulation环境。从模型库中拖放以下模块到流程模拟窗口Mixer混合器命名为MIXER。RStoic化学计量反应器命名为REACTOR。RadFrac严格精馏塔命名为COLUMN。Fsplit流股分割器命名为SPLITTER。按照以下顺序连接流股新建一股进料流股连接到MIXER的进口命名为FEED。从MIXER的出口连接到REACTOR的进口流股命名为REAC-FEED。从REACTOR的出口连接到COLUMN的进料板例如第 5 板流股命名为REAC-OUT。从COLUMN的塔顶出口Vapor Draw连接到SPLITTER的进口流股命名为DIST-OUT。从SPLITTER的第一个出口连接到MIXER的第二个进口流股命名为RECYCLE。这就是形成循环的关键流股。从SPLITTER的第二个出口作为产品采出流股命名为PRODUCT。从COLUMN的塔釜出口作为产品采出流股命名为BOTTOM。至此带有循环的流程图搭建完成。Aspen Plus 会自动识别出循环并将RECYCLE流股标记为撕裂流股。3.3 定义流股与模块参数进料流股FEED设定双击FEED流股在Specifications页面输入温度90 °C压力300 kPa总流量100 kmol/hr组成A 0.95摩尔分数I 0.05。反应器REACTOR设定双击REACTOR模块。在Specifications页面设定操作条件如温度、压力与进料相同。在Reactions页面点击New...创建一个新的反应集。定义反应A - B设定反应类型为“Power law”并指定转化率为 0.7070%。精馏塔COLUMN设定双击COLUMN模块使用“Setup Wizard”进行快速设置。塔板数10包括冷凝器和再沸器。进料板位置第 5 板。冷凝器类型全凝器Total。再沸器类型釜式Kettle。规定回流比2塔顶采出流量DIST-OUT30 kmol/hr。塔釜采出流量BOTTOM由物料平衡计算得出通常不直接规定确保塔内平衡即可。这里可以先输入一个估计值如 50 kmol/hr后续迭代会修正。分流器SPLITTER设定双击SPLITTER模块。在Specifications页面设定分流比。例如设定RECYCLE与PRODUCT的流量比为 21。这意味着进入SPLITTER的DIST-OUT流股有 2/3 进入RECYCLE1/3 进入PRODUCT。4. 收敛策略撕裂流股与算法选择4.1 查看与修改收敛设定首次运行模拟按 F5很可能会因为循环不收敛而报错。我们需要检查并优化收敛设置。点击顶部菜单Flowsheeting OptionsConvergenceConvergence。在Default收敛方法下可以看到 Aspen Plus 已经自动将RECYCLE流股添加为撕裂流股。点击Tear Streams标签页可以查看和修改撕裂流股的初始估计值。为撕裂流股提供良好的初始值至关重要。4.2 为撕裂流股提供初始估计双击RECYCLE流股在Estimate页面提供合理的初始值。即使不精确也比完全空白的默认值好。温度80 °C压力280 kPa总流量150 kmol/hr根据物料平衡粗略估计组成A 0.90I 0.10。4.3 选择合适的收敛算法在ConvergenceMethods页面有三种主要的收敛算法Wegstein默认且最常用的方法适用于大多数单变量循环问题收敛速度快。对于本案例通常首选此法。Direct直接迭代法非常稳定但可能很慢。当 Wegstein 法振荡或不收敛时可以尝试此法。Broyden适用于大型、强非线性系统比 Wegstein 法更强大但设置更复杂。对于本案例保持使用Wegstein算法并可以适当调整参数如将加速参数的上限Qmax从默认的 0 调整为 5 或 10有时能改善收敛性。5. 运行模拟与结果分析完成以上设置后再次运行模拟F5。如果一切设置正确状态窗口会显示“Results Available”表示流程已收敛。查看循环流股结果查看RECYCLE流股的结果报告ResultsStream ResultsRECYCLE。你会看到收敛后的流量、组成、温度、压力等详细信息。这表明软件成功求解了循环方程组。查看全局物料平衡通过比较总进料FEED和总出料PRODUCTBOTTOM的流量验证模拟的准确性。理论上在稳态下输入应等于输出。6. 常见收敛问题与排查路径循环流程不收敛是常见问题排查应遵循以下路径问题现象可能原因检查与解决步骤“Tear Stream ... has not converged”1. 撕裂流股初始估计值太差。2. 模块规定过于严格或相互冲突。3. 物性方法选择不当导致计算不稳定。1.提供更好的初始估计根据工艺知识给出更接近真实值的温度、压力和流量。2.放松收敛容差在ConvergenceOptions中暂时增大撕裂流股的收敛容差如从 0.0001 改为 0.001收敛后再改回。3.检查模块规定确保精馏塔的规定回流比、采出量是合理的不会导致塔内干板或液泛。计算在几个值之间振荡Wegstein 法加速过度。1.限制 Wegstein 参数在收敛方法设置中将 Qmax 设为 0相当于直接迭代法或设为一个较小的正数如 5。2.使用 Direct 方法切换至 Direct 迭代法虽然慢但稳定。物性计算错误在循环计算的某个中间状态温度、压力或组成进入了物性模型的非适用区间。1.检查物流结果查看迭代过程中的中间结果找到物性计算报错的位置。2.调整物性方法考虑更换更稳健的物性方法或使用物性分析工具检查混合物的气液平衡行为。流程设计本身不合理例如循环量远大于新鲜进料量导致系统对初始值极其敏感。1.简化流程先断开循环分别模拟反应器和精馏塔获得各流股的合理估计值再连接循环。2.逐步增加循环使用“计算器”模块或设计规定逐步将循环流股引入系统。7. 循环流程模拟的最佳实践先简化后复杂对于复杂循环流程先断开循环流股将其作为进料或产品使流程变为简单的开式流程。分别收敛各个单元模块获得所有流股尤其是计划作为撕裂流股的流股的合理估计值。然后连接循环并利用已获得的估计值进行初始化。良好的初始值是成功的一半花费时间提供准确的撕裂流股初始值远比盲目调整收敛算法有效。理解工艺是根本对工艺流程的深刻理解有助于判断初始值是否合理以及当模拟不收敛时问题是出在模型设置上还是工艺设计本身就有缺陷。善用诊断工具Aspen Plus 的ConvergenceHistory可以查看迭代历史帮助判断是单调发散还是振荡从而有针对性地调整算法参数。考虑使用“设计规定”本案例中塔顶采出比是固定的。在实际设计中你可能希望通过调节采出比来控制产品纯度。这时就需要使用“设计规定”功能它会引入另一个控制循环增加收敛难度需要更谨慎地设置。掌握带有循环的流程模拟是有效使用 Aspen Plus 进行复杂化工过程设计的关键一步。从理解原理出发通过规范的建模步骤、合理的初始估计和系统的排查方法可以解决大多数收敛挑战从而获得可靠的计算结果用于工程分析。