第 06 篇|Aspen Plus 进阶与工程案例:精馏塔严格建模与 RadFrac 配置
版本声明块
- 软件:aspenONE V15 版 Aspen Plus(Windows 64 位桌面版;截至 2026-10 公开渠道最新大版本为 V15)
- 单位体系:示例以 METCBAR 为基线并逐处显式标注单位
- 本文目标:把「RadFrac 到底该填哪些数、为什么」讲透——用平衡级方程解释「两条操作规定」的来历,用可解性判据解释「为什么有的组合不许用」,用算法矩阵解释「该选哪种收敛方案」
- 事实来源:AspenTech 官方培训课程页(Aspen Plus: Distillation Modeling;Modeling and Optimizing Chemical Processes Using Aspen Plus)、官方产品页(Distillation Modeling in Aspen Plus、Aspen Plus)、公开镜像的高校精馏课件;界面页签命名在不同版本间有差异,凡未在本机核实的都以「以本机界面为准」标注
一句话结论:RadFrac 把每一层塔板都当作一个平衡级(equilibrium stage),对一台有冷凝器和再沸器的简单塔,在理论级数、进料(位置/流量/组成/状态/焓)与压力剖面都给定了之后,整塔恰好还留下两个自由度——冷凝器与再沸器各自的规格已锁掉这两个设备的自由度,剩余的两条必须由两条线性独立的操作规定补足(典型如「回流比 + 塔顶采出量」);这两条规定应分别管住「分离强度」与「产品去向」,而任何一对线性相关的规定(如进料固定时的「塔顶采出 + 塔底采出」)都会被判为不可解。
〇、本篇要解决的认知问题
第 05 篇把「不收敛」拆成了可定位的诊断,你也学会了在控制面板里读 Max Err/Tol、在 Run Status 里找报错。但有一个更靠前的问题一直没回答:你填进 RadFrac 的那些数,为什么是那些数?很多人建塔的方式是「照抄别人的模板」或「随便试两个数,能收敛就行」——而一旦塔不收敛、或分离指标对不上,就完全不知道从哪改起。本篇围绕五个问题展开:
- RadFrac 的每一级到底有几个自由度?为什么给定了理论级数、进料与压力之后,整塔「恰好剩下两个自由度」需要两条操作规定补足?多给或少给一个会怎样?
- 冷凝器类型(Total / Partial-Vapor / Partial-Vapor-Liquid / None)与再沸器类型(Kettle / Thermosyphon / None)的组合,分别对应什么真实设备?为什么它们各自会「吃掉」模型里的一个自由度,从而决定后续要补几条规定?
- 两条操作规定怎么选才可解?为什么「回流比 + 塔顶采出量」是最常见的一对,而「塔顶采出量 + 塔底采出量」在进料固定时是重复的?官方那句「Some combinations are not allowed」背后的数学是什么?
- 压力剖面、进料 convention(Above-Stage / On-Stage)、侧线采出(side draw)、中间回流(pumparound)、三相(Vapor-Liquid-Liquid)、塔板效率,各自占了模型里的哪一部分、会改变自由度计数吗?
- 收敛方案Standard / Strongly Non-ideal Liquid / Petroleum-Wide-Boiling / Azeotropic / Cryogenic各自的底层算法与初始化方法是什么,什么场合该用哪一个?
这五个问题对应三类真实故障:计数类(第 1、2 问,自由度给错了,塔要么无解要么过约束)、选型类(第 3、4 问,规定组合或附加设备设错,物理上就不可解)、数值类(第 5 问,算法与初始化选错,明明可解却收敛不了)。它们的共同特征是:界面不报语法错,塔能给出一组「结果」,但那股结果要么对不上工程期望,要么根本是另一个操作点。
一、机制解析
1.1 平衡级与 MESH 方程:RadFrac 的基本假设
RadFrac 是严格多级分离模型,它的基本假设是:每一层理论板都是一个平衡级——离开该级的汽相与液相达到相平衡。官方产品页称其为「Rigorous two-three phase fractionation(most popular)」,即最常用的严格两相/三相分馏模型。
对单个平衡级iii,描述它的是那一组经典方程,业界取其首字母合称MESH:
- Material balance(物料衡算):每个组分进出该级的物料守恒;
- Equilibrium(相平衡):yi,j=Ki,j xi,jy_{i,j}=K_{i,j}\,x_{i,j}yi,j=Ki,jxi,j,其中KKK为相平衡常数(取决于物性方法,见第 02 篇);
- Summation(归一):∑jxi,j=1\sum_j x_{i,j}=1∑jxi,j=1、∑jyi,j=1\sum_j y_{i,j}=1∑jyi,j=1;
- Heat balance(能量衡算):进出该级的焓守恒。
一台塔把这 N 层板的 MESH 方程联立起来,未知量包括每级的温度、汽液相组成、以及级间的汽液相流率。未知量数与方程数之差,就是这台塔的自由度——而 RadFrac 要你补的输入,正是用来「吃掉」这些剩余自由度。
塔顶冷凝器 (Condenser) ← 一个带规格的平衡级 │ 回流 L v ═══════════ 第 2 级 ═══════════ 第 3 级 ... ⋮ ──► 进料 F ──► 第 f 级(进料位置,给定) ... ⋮ ═══════════ 第 N-1 级 ═══════════ 第 N 级 │ 塔底再沸器 (Reboiler) ← 一个带规格的平衡级 │ 塔底采出 B 给定的:级数 N、进料(F, z, T, P, 位置 f)、压力剖面 待补的:两条操作规定 → 整塔自由度归零,可解1.2 为什么恰好是「两个自由度」
这个「恰好为 2」不是拍脑袋定的,它来自「哪些量已经给定、哪些设备已经自带规格」的账。官方与高校课件给出的 RadFrac 基本输入清单是:级数(Number of stages)+ 冷凝器类型 + 再沸器类型 + 两条操作规定。所以账要这么算:
- 已经给定的(不占自由度):理论级数 N、进料流量与组成与温度压力、进料位置、压力剖面;
- 设备规格已锁定的:冷凝器是一个平衡级,它的类型(全凝/分凝)本身就规定了塔顶蒸气如何冷却、如何在冷凝器温度下达到相平衡,等于把冷凝器这一级的自由度锁死;再沸器同理(Kettle 釜式、Thermosyphon 热虹吸各有一种固定的汽化-返流方式),也自带规格。冷凝器与再沸器各自吃掉一个自由度。
- 剩下待补的:把上面这些全部代入 N 层 MESH 方程后,残余的自由度正好是 2,需要用户再用两条操作规定补足。
这里有一个必须点破的物理含义:这两条规定不是「随便挑两个数」,而是把两个独立的「操作维度」钉死。
| 维度 | 物理含义 | 典型规定 |
|---|---|---|
| 分离强度 | 塔内汽液负荷有多大、回流/沸腾有多强,决定能达到多高的分离程度 | 回流比(Reflux ratio)、沸腾比(Boil-up ratio)、冷凝器负荷、再沸器负荷 |
| 产品去向 | 进料在塔顶与塔底之间如何分配 | 塔顶采出量(Distillate rate)、塔底采出量、D/F、B/F |
工程直觉:一台塔要同时回答「分得多干净」和「产品各拿多少」两件事,这两件事恰好对应两个独立自由度。只给强度不给去向,塔顶塔底流量无法确定;只给去向不给强度,分离程度无法确定。
给多与给少是两类症状完全不同的硬错,先把它俩认出来,就能在输入层面排除一半问题:
- 多给一个(过约束):既给了回流比、又给了再沸器负荷、还给了塔顶采出量。此时 MESH 方程组互相矛盾,RadFrac 通常直接报错「规定冲突」,或迭代时两三个变量轮流成为「本轮最大误差」、Max Err/Tol 长期驻留在某一量级不下降。
- 少给一个(欠约束):只给了一条操作规定。此时方程组欠定、解不唯一,不同初值会收敛到不同的操作点。这类症状更阴险:塔「能收敛」,可你两次运行、或换个初值,分离度就变了——「能收敛但结果漂移」的背后常常就是它。
由此得到一条判定自由度的口诀:输入清单里必须有且仅有两条操作规定,且二者线性独立。规定数错了,后面所有收敛努力都建立在错误的地基上。
1.3 冷凝器 / 再沸器组合矩阵
冷凝器与再沸器的类型决定「这台塔的边界在哪、边界上又锁掉了什么」。这张矩阵是建模第一步:
| 冷凝器类型 | 物理设备 | 额外要给的规格 | 说明 |
|---|---|---|---|
| Total | 全凝器 | 无(默认饱和液体馏出物与饱和回流) | 最常用;塔顶全为液相产品 |
| Partial-Vapor | 分凝器 | 冷凝器温度或馏出物汽化分率 | 塔顶有汽相产品(如含不凝气时) |
| Partial-Vapor-Liquid | 分凝器 + 分液 | 冷凝器温度或馏出物汽化分率 | 塔顶同时出汽相与液相(三相 VLL 计算时的典型配置) |
| None | 无冷凝器 | — | 塔顶直接出汽相(如汽提塔、无塔顶冷凝的吸收塔) |
| 再沸器类型 | 物理设备 | 说明 |
|---|---|---|
| Kettle | 釜式再沸器 | 常用;相当于一个理论级 |
| Thermosyphon | 热虹吸再沸器 | 常用;循环式 |
| None | 无再沸器 | 塔底不出产品汽化回流(如无再沸吸收塔) |
官方课件对冷凝器规格的补充是:若冷凝器是 Partial-Vapor-Liquid,需要给冷凝器温度或馏出物汽化分率;而对液相馏出物与回流,可以给过冷温度或过冷度,默认是饱和液体馏出物与饱和回流。三条来自官方课件的最佳实践值得记牢:
- 回流比与沸腾比绝不能设为 0。要零回流就用
Condenser = None(那是另一台塔),不要在回流比栏里填 0。 - 含不凝气时用分凝器,并把馏出物汽化分率设到足以把不凝气全部带进汽相馏出物。
- 蒸汽进料、塔顶出汽相产品的塔,正确做法是塔顶建汽相产品股 + 冷凝器选 None,而不是给一个全凝器。
1.4 两条操作规定的可解性判据
RadFrac 允许从十个量里挑两个:塔顶采出量 D、塔底采出量 B、回流量 R、沸腾量 Vb、回流比 RR、沸腾比 BR、D/F、B/F、冷凝器负荷 Qc、再沸器负荷 Qr。官方原文强调一句:「Some combinations are not allowed.」为什么?因为挑出的两个量必须在数学上线性独立,否则方程组要么矛盾、要么仍欠定。判据可以归纳成两条经验法则:
- 一「去向」+ 一「强度」通常可解:如「D + RR」「D/F + RR」「B + Qr」「D + BR」等。这是最常见也最稳的组合。
- 两个「同类」量通常不可解或冗余:
- D + B:在进料固定时二者之和必须等于进料中的产品分配,稳态下不独立(差一个物料衡算关系)——所以应改用其中一个 + 一个强度量;
- D + D/F:D/F 就是 D 与进料 F 的比值,F 已给定,二者本质是同一个量,重复;
- Qc + Qr:在能量平衡闭合的场合,二者通过全塔焓平衡耦合,也常常不独立。
官方课件给出的取向偏好:「Flow rate specifications are in general preferred over duty specifications, especially for wide-boiling systems.」——流率规格通常优于热负荷规格,尤其在宽沸程体系。原因是宽沸程下热负荷对温度极不线性,给热负荷容易把塔推进不收敛区;给流率则稳定得多。而D/F、B/F的独特价值在于:当进料流量本身不确定时它们是便利规格,且便于给初值和上下界。
为什么「线性相关」会让求解直接失败,而不是「解不唯一、随便给个解」?因为 RadFrac 内部是把整塔 MESH 方程联立后迭代求解的(Standard 用 inside-out、Newton 方案用牛顿法),每加一条操作规定,就等于往方程组里再添一个方程。若两条规定在数学上线性相关(如进料固定时的 D 与 B),对应方程的导数矩阵就会出现秩亏(奇异),迭代法无法在该方向上稳定更新变量——现场表现就是「规定冲突」报错,或两三个变量轮流成为本轮最大误差、Max Err/Tol 长期驻留不降。反过来说,「两条规定必须线性独立」的数学含义,正是它们拼出的方程组雅可比必须满秩;秩亏时无论把迭代次数加多大都救不回来(呼应第 05 篇「结构性问题调数值无效」)。这就是官方提示「Some combinations are not allowed」的底层原因。
1.5 附加自由度:压力、进料 convention、侧线、pumparound、三相、效率
上面算的是「一台最简单的塔」。真实装置还会引入若干「附加自由度」,它们要么改变计数、要么改变方程的边界条件,必须逐个认账:
(1)压力剖面。Setup / Pressure页有三种给法:Top/Bottom Pressure(给塔顶/塔底压力,中间按柱压降处理;不填压降时第 1 级/冷凝器压力即塔压)、Pressure profile(逐级给压力,未填的级按已填值线性插值或外推)、Section pressure drop(按段给压降,塔顶压力必填)。压力是「给定」的,不占两条操作规定的名额——因为压力由设备结构(回流罐、塔板压降)决定,是设计输入而非操作自由度。
(2)进料 convention(进料引入方式)。Setup / Streams页的Convention有四档:
| convention | 行为 |
|---|---|
| Above-Stage(默认) | RadFrac 把进料流引入相邻两级之间:液相部分流到指定级,汽相部分流到指定级上一级 |
| On-Stage | 汽相与液相部分都进入指定级 |
| On-Stage-Liquid / On-Stage-Vapor | 同 On-Stage,但不做闪蒸,进料被当作完全处于指定相 |
官方提示:塔板数很多时两种 convention 结果相近;但当已知进料为单相时优选 On-Stage——它能省掉闪蒸计算,并避免超临界体系的闪蒸问题。
(3)侧线采出(side draw)。Setup / Streams页可给「产品级位置 + 采出流量」。侧线是从模型里往外拿走一股,它会改变该级的物料衡算——若采出流量给定,则它又额外消耗一个自由度;因此有侧线的塔常常需要更多规格(用 Vary/Design Spec 处理)。实务上的常见做法是:把侧线流量设为操纵变量(Vary),用一个设计规定(如侧线产品纯度或某关键组分的回收率)去反算它——这正是「规格数 = 操纵变量数」在侧线上的一次典型应用(详见 1.7 节)。
(4)中间回流 / 中段回流(pumparound)。官方描述:RadFrac 可处理从任一级到任一级的 pumparound;一次 pumparound 可以是某级液相、第一液相、第二液相或汽相的部分或全部抽出。若为部分抽出,必须再给两个规格(从流量、温度等中选)。pumparound 是炼油常压/减压塔的核心特征(如原油塔的中段回流),它会显著改变塔内热量分布,是「附加自由度」里最容易漏账的一项。
pumparound 为什么最容易漏账?它把某一级的液相抽走、换热后再送回(可送回同一级,也可送回另一级)。抽走的液相量与换热后的返流温度,都是新引入的独立变量;若抽走的是「部分液相」却只给了一个规格,模型仍然欠定——所以官方要求再给两个规格。另一个高频踩坑是级号约定:官方原油塔示例特别提醒,指定抽出级与返回级时要加 1——因为冷凝器算作第 1 级。这与「Nstages = Ntrays + 2」是同一套计数约定;弄错级号会让中段回流挂到错误的塔板上,全塔热分布整段跑偏。
(5)三相(Vapor-Liquid-Liquid)。Setup / Configuration的Valid phases有Vapor-Liquid(默认)与Vapor-Liquid-Liquid。选三相时,模型不假定两个液相的性质,且倾析器(Decanter)可挂在任意一级上。三相计算会引入第二液相的组成变量——自由度账要重算,这也是共沸精馏(如乙醇脱水用苯夹带剂形成两液相)必须用它的原因。
(6)塔板效率(Efficiencies)。RadFrac 默认每级完全平衡;对非理想分离可用汽化效率(Vaporization)与 Murphree 效率修正,二者可在Convergence或Setup下按组分、按级或按段给。官方建议优先用 Murphree 效率,因为汽化效率会改变温度剖面、从而改变相对挥发度。注意一个实用的对应关系:对板式塔,Nstages = Ntrays + 2(多出的两级是冷凝器与再沸器)。效率不改变自由度计数,但会改变分离结果。
1.6 收敛算法矩阵:六种方案各自的底层算法与初始化
RadFrac 的收敛方案由两块组成:底层算法与初始化方法(给温度剖面与汽液组成剖面的初值)。方案可在Setup / Configuration页的下拉框中选,也可在Convergence文件夹里细化(该文件夹在 V15 与旧版间页签命名有差异,如 Basic / Advanced,以本机界面为准)。
| 方案 | 底层算法 | 初始化方法 | 适用场合 |
|---|---|---|---|
| Standard(默认) | Standard 算法(Aspen Technology 的 inside-out 算法) | Standard | 应始终首先尝试 |
| Strongly Non-ideal Liquid | Nonideal 算法 | Standard | 强非理想两相塔且 Standard 收敛慢时 |
| Petroleum / Wide-boiling | Sum-Rates 算法 | Standard | 石油/石化、宽沸程、多组分、多设计规定 |
| Azeotropic | Newton 算法 | Azeotropic | 两相共沸精馏(如乙醇脱水用苯作夹带剂) |
| Cryogenic | Standard 算法 | Cryogenic | 深冷应用(如空气分离) |
| Custom | 用户混搭 | 用户选 | 高级用户按需组合 |
配合理解Standard 初始化策略:把所有进料合并成一股复合流;对该复合流做闪蒸得到汽/液组成,用这组 V/L 组成作为所有塔板的初值;再对复合流做泡点与露点计算,取T_top = T_bubble、T_bottom = T_dew,中间各级温度在这之间线性插值。
初值的三种量(T、x-y、V/L)通过Convergence / Estimates文件夹给。官方指南的边界很明确:温度和组成估值一般不需要给;吸收/汽提塔应给温度估值以覆盖默认值;吸收/汽提塔可能还需要液/汽流量估值;强非理想或极宽沸程体系可能需要组成估值。
1.7 Design Spec / Vary:把「规定」变成「操纵」
一台塔除了两条基本操作规定,工程上还常有「产品纯度必须达到 X」这类性能规定。RadFrac 内置了设计模式(等价于外部设计规定),用Design Specs + Vary两个页签实现:每条设计规定配一个操纵变量(manipulated variable),且操纵变量必须是 RadFrac 输入里已经被指定过的量(典型:回流比、塔顶采出量、进料量)。官方原文的关键句:「For each design target, there should be a corresponding manipulated variable」,以及「The number of Specs should, in most cases, be equal to the number of varies」(规格数通常等于操纵变量数)。
这不是小事:规格数与操纵变量数不相等,就是过约束或欠约束,塔要么不收敛、要么解不唯一。而操纵变量的初始填值会被当作初值使用(如例子里把回流比 1.50 当初始猜测,最终解出 1.895)。这与第 05 篇讲收敛块嵌套直接相关:设计规定的操纵变量如果由外层回路决定,就要把规定嵌到外层(Sequencing = Outside),否则规定在内层反复调也白费。
1.8 收敛 ≠ 分离可行(铁律 4 落点)
再次钉死铁律 4:收敛 ≠ 算对。RadFrac 收敛只保证 MESH 方程残差达标、塔内剖面自洽,它不保证这台塔在工程上分离得出来。一台「收敛」的塔完全可能意味着:进料位置错了(本该在中部进料却从塔顶进)、回流比设到了不可行的极端、或者规定落在精馏边界之外(分离根本不可能)。官方课件给出的 RadFrac 收敛注意项第一条就是「Ensure that column operating conditions are feasible」——先确认操作条件可行,再谈数值收敛。收尾必须做三件事:质量平衡校验(D + B + 侧线 = 进料)、能量平衡校验(Qc 与 Qr 与进出料焓差闭合)、以及分离可行性判断(关键组分回收率、产品纯度是否落在相对挥发度允许的范围内)。
二、完整操作与逐项剖析
2.1 从零建一台 RadFrac:界面操作顺序
界面操作步骤:
- 在
Blocks页拖入Columns / RadFrac,连好进料、塔顶产品、塔底产品与(可选)侧线、pumparound 回路。 - 进
Setup / Configuration:Calculation type选Equilibrium(平衡级)或Rate-Based(速率基,需传质系数);填Number of stages(理论级数);选Condenser与Reboiler类型;选Valid phases(Vapor-Liquid 或 Vapor-Liquid-Liquid);在下方Operating specifications里填两条操作规定。 - 进
Setup / Streams:给进料的级位置与Convention(默认 Above-Stage);给塔顶/塔底产品与侧线的级位置与流量。 - 进
Setup / Pressure:给塔顶压力(及压降)——压力是设备输入,不占操作规定名额。 - 若为分凝器,进
Setup / Condenser给冷凝器温度或馏出物汽化分率、以及过冷度(默认饱和)。 - 进
Convergence文件夹:选收敛方案(先 Standard),必要时给Estimates(吸收塔给温度、非理想体系给组成)。 - 运行;按第 05 篇的方法读控制面板与 Run Status,收敛后再做 2.3 的平衡校验。
关键设定表:
| 页签 | 必填项 | 缺省/常见值 |
|---|---|---|
Configuration | Calculation type、级数、冷凝器、再沸器、Valid phases、两条操作规定 | Equilibrium / V-L |
Streams | 进料级位置与 Convention、产品级位置、(侧线流量) | Above-Stage |
Pressure | 塔顶/塔底压力或压力剖面 | 塔压、压降 |
Condenser | 分凝器的冷凝器温度或馏出物汽化分率、过冷度 | 饱和液体馏出物/回流 |
Convergence | 收敛方案、迭代上限、Estimates | Standard |
期望输出:塔结果页Profiles给出逐级的温度、压力、汽液流量与组成剖面(单调或符合精馏常识);Summary给出塔顶/塔底温度、冷凝器与再沸器负荷、回流比;控制面板显示塔块收敛(Max Err/Tol < 1.0,呼应第 05 篇)。
2.2 三条不能违反的配置纪律
- 回流比与沸腾比不为 0:零回流用
Condenser = None表达,不要在数值栏填 0。 - 含不凝气用分凝器:并把馏出物汽化分率设到足以把不凝气全部带进汽相馏出物。
- 流率规格优于热负荷规格:宽沸程体系尤其如此;给热负荷容易把塔推入不收敛区。
2.3 收敛后的三项校验(铁律 4 落地)
| 校验项 | 做法 | 通过标准 |
|---|---|---|
| 质量平衡 | D + B + Σ(侧线) 与进料逐组分比对 | 各组分相对残差远小于收敛容差 |
| 能量平衡 | 冷凝器负荷 + 再沸器负荷 + 进出料焓差闭合 | 残差与容差同量级 |
| 分离可行性 | 关键组分回收率、产品纯度是否合理 | 落在相对挥发度允许范围内,无「一级塔分馏」式谬误 |
三、常见报错与排查
下面五条覆盖 RadFrac 建模高频现场,每条按「现象 → 根因 → 解法」给。
1. 塔一直不收敛,Max Err/Tol上下振荡或缓慢爬升。
- 现象:控制面板里塔块迭代到上限仍不收敛,
Tear History无下降趋势。 - 根因:最常见是操作条件不可行——两条规定把塔推到一个物理上不存在的操作点,或进料位置/级数不足以完成该分离。官方收敛注意项第一条即「确认操作条件可行」。
- 解法:按官方 RadFrac 收敛注意项的顺序排查——① 确认操作条件可行;② 检查物性是否在整个浓度-温度区间都算得对(铁律 1、5、10:物性先定案、不外推、可溯源);③ 若 Err/Tol 在稳定下降只是慢,就加大最大迭代数(在 RadFrac 的
Convergence设定里);④先不加 Design Spec 与 Vary,让塔裸收敛,收敛后再逐步加规格与操纵变量;⑤ 换收敛方案(见下条)。
2. 强非理想体系(如醇-水)用 Standard 收敛极慢。
- 现象:塔能收敛但迭代数很高、耗时很长,温度剖面初值像「不听话」。
- 根因:Standard 的 inside-out 算法对强非理想汽液平衡的适应不如专用算法;且组成的线性插值初值离真实剖面太远。
- 解法:把方案换成Strongly Non-ideal Liquid(底层 Nonideal 算法、Standard 初始化),并按官方指南在
Convergence / Estimates给组成估值(针对非理想体系)。若体系接近共沸,则改选Azeotropic(Newton 算法 + Azeotropic 初始化)。
3. 吸收塔 / 汽提塔不收敛,或收敛到明显不合理的温度剖面。
- 现象:吸收塔(无再沸器或极端操作)默认初始化给出的温度剖面不符合物理。
- 根因:默认 Standard 初始化是从进料泡露点线性插值的「精馏型」初值,对吸收/汽提这类传质主导、温差小的塔不适用。
- 解法:按官方指南在
Convergence / Estimates / Temperature页给关键级的温度估值(吸收/汽提塔应当给温度估值以覆盖默认值);必要时再给液/汽流量估值。这一步常常是吸收塔收敛的胜负手。
4. 报错提示两条操作规定冲突或某规定无法满足。
- 现象:界面提示规定组合非法,或塔收敛到明显矛盾的解。
- 根因:要么选了线性相关的两个规定(如进料固定下的「D + B」),要么在多规格时规格数与操纵变量数不等(官方要求二者基本相等)。
- 解法:把两条规定改成「一去向 + 一强度」的独立组合(经验法则);有设计规定时保证Specs 数 = Vary 数,且每个操纵变量必须是 RadFrac 输入里已指定的量(如回流比、塔顶采出量、进料量)。别用两个同类量硬凑。
5. 改了进料位置或级数后,塔「能收敛但结果跳变」。
- 现象:只动了进料位置或级数,分离指标突然大幅变化,甚至不再收敛。
- 根因:进料位置与级数是结构参数,改动后塔板剖面的初值与旧的收敛解不再对应;若还叠加了 pumparound 或侧线,旧规定可能彼此冲突(参考第 05 篇:结构参数变更后相关规定必须同步修改并重置)。
- 解法:改结构参数后重置(Reinitialize)再跑;给进料一个基于「泡点/露点位置」的合理级号(进料应放在其组成对应泡点温度与塔内温度剖面相交附近);有 pumparound 时同步核对抽出级与返回级(注意级号以冷凝器为第 1 级,实际数级要加 1——官方原油塔示例的做法)。另可先退回简单模型(DSTWU)出一次解作初值,再切回 RadFrac(官方策略)。
四、动手练习
练习 1:做一次 RadFrac 自由度账目核对。
做法:取一台有全凝器(Total)与釜式再沸器(Kettle)的二元塔,记录你给的全部输入:级数、进料流量/组成/温度/压力/位置、压力、两条操作规定。然后逐项判定每条输入属于「给定(不占自由度)」「设备规格(占掉自由度)」还是「操作规定(补足自由度)」。
判定标准:产出一张账目表,明确标出「共给出 2 条操作规定」,且这 2 条必须分属「分离强度」与「产品去向」两个维度(如回流比 + 塔顶采出量)。若你选的两条同属一个维度(如 D + B),必须判为不可解并指出正确替代。账目表还须逐条标注每条输入的归属(给定 / 设备规格 / 操作规定)并给出计数:给定项应为级数、进料参数与压力;设备规格项应为冷凝器与再沸器类型;操作规定项应恰为两条。若计数结果不是「恰好两条」,则无论塔是否收敛,都判为自由度设定有误,必须改到恰好两条为止。
练习 2:用两种收敛方案比较同一台强非理想塔的收敛行为。
做法:取一台强非理想(如醇-水或醇-酯)两相塔,先在缺省的Standard方案下跑,记录迭代数与是否收敛;再换成Strongly Non-ideal Liquid(必要时在Estimates给组成估值),在相同输入下重跑。
判定标准:产出一张对照表,含「方案、底层算法、初始化、最大迭代数、是否收敛、最终 Max Err/Tol、耗时」。若体系含共沸,再补一列Azeotropic方案。结论必须写清:该体系下哪个方案是首选,并说明判断依据(Err/Tol 是否稳定下降、迭代数差异)。
练习 3:验证「流率规格优于热负荷规格」,并完成平衡校验。
做法:对同一台塔,方案 A 用「回流比 + 塔顶采出量」,方案 B 用「回流比 + 再沸器负荷」;比较二者的收敛难易与结果一致性。随后对方案 A 的收敛结果做质量平衡(D + B = 进料,逐组分)与能量平衡(Qc + Qr 与焓差闭合)校验。
判定标准:产出一张两方案对照表(含「是否收敛、迭代数、塔顶纯度、回流比、Qc、Qr」);方案 A 必须收敛更稳或至少不差。平衡校验表三行残差必须给「绝对值、相对值、是否在容差内」。若能量平衡相对残差显著大于收敛容差,判定为「收敛但不可用」并指出最可能的输入错误。
五、小结与下一篇预告
这一篇把 RadFrac 从「填数就行」升级成了「知道每个数为什么这么填」。核心结论有四条:
第一,RadFrac 每级是一个平衡级,整塔自由度由 MESH 方程决定;在级数、进料与压力给定、冷凝器与再沸器自带规格之后,整塔恰好留下两个自由度,必须由两条操作规定补足——这两条应分别管「分离强度」与「产品去向」。
第二,冷凝器/再沸器组合决定塔的边界与规格:Total / Partial-Vapor / Partial-Vapor-Liquid / None 与 Kettle / Thermosyphon / None 各有对应的物理设备;回流比与沸腾比不得填 0,含不凝气用分凝器,单相进料优选 On-Stage convention。
第三,两条规定必须线性独立:流率规格通常优于热负荷规格;「D + B」「D + D/F」「Qc + Qr」这类同类组合往往不可解;有设计规定时保证Specs 数 = Vary 数,操纵变量必须是已指定的输入量。
第四,收敛方案分三类用途:Standard 首试,强非理想用 Strongly Non-ideal Liquid,石油宽沸程用 Petroleum,共沸用 Azeotropic,深冷用 Cryogenic;初始化方法要与体系匹配,吸收/汽提塔必须给温度估值。而无论用哪种方案,收敛之后都要做质量/能量平衡与分离可行性校验(铁律 4)。
再把两张「最容易记混」的清单并列收尾,方便日后自查。第一张是「给定 vs 待补」:级数、进料(位置/流量/组成/状态/焓)、压力剖面属于给定,不占名额;冷凝器与再沸器类型属于设备规格,各自锁死一个自由度;剩下恰好两条操作规定由你补。第二张是「占自由度 vs 不占自由度」:侧线(给定流量时)、pumparound(部分抽出时)、设计规定(每条配一个操纵变量)都会消耗自由度,需要相应地增加规格数;而进料 convention、三相 Valid phases、塔板效率则属于改方程或改结果但不改变基本自由度账的设定。把这两张清单钉在草稿纸上,RadFrac 的输入就不会再「凭感觉乱填」。
留下的问题是那些「非平衡级」的单元:本篇整篇都在平衡级的框架里打转,而真实装置里换热器的温差与面积、反应器里转化率与动力学、泵与压缩机的效率,都不是「一个平衡级」能描述的——尤其当反应与分离耦合(反应精馏)时,RadFrac 的自由度账还要再叠一层。
下一篇:《第 07 篇|Aspen Plus 进阶与工程案例:反应器建模从计量平衡到动力学》——RStoic(计量)、RYield(收率)、REquil(平衡)、RGibbs(吉布斯最小化)、RCSTR/RPlug(动力学)五类反应器的适用边界,转化率/平衡/动力学三种输入方式的取舍,以及铁律 3 与铁律 10 的落点:没有动力学数据就不要用动力学反应器编参数。
本篇认知问题回显(FAQ)
Q1:Aspen Plus 的 RadFrac 为什么只要给两条操作规定塔就能解,这两条的物理含义是什么?
A:因为在给定理论级数、进料与压力、且冷凝器与再沸器各自自带规格之后,整塔剩余的待定自由度恰好为 2;这两条应分别管住「分离强度」(回流/沸腾)与「产品去向」(采出量/采出比)两个独立维度。
Q2:Aspen Plus 的 RadFrac 里冷凝器与再沸器类型分别有哪些,各自会带来什么必填规格?
A:冷凝器有 Total、Partial-Vapor、Partial-Vapor-Liquid、None,分凝器需给冷凝器温度或馏出物汽化分率;再沸器有 Kettle、Thermosyphon、None;含不凝气用分凝器,无冷凝器或无再沸器时选 None 而不是把回流比填 0。
Q3:Aspen Plus 的 RadFrac 两条操作规定该怎么选,为什么有的组合不被允许?
A:两条规定必须线性独立,常见可解组合是「一回流比 + 一采出量」或「一采出比 + 一热负荷」;「D + B」「D + D/F」「Qc + Qr」这类同类或在物料/能量衡算上相关的组合会被判不可解,官方因此提示 Some combinations are not allowed。
Q4:Aspen Plus 的 RadFrac 里 pumparound、side draw 和三相分别该怎么配,会改变自由度吗?
A:pumparound 可从任一级抽出到任一级,部分抽出时需再给两个规格;side draw 在给定采出流量时会额外消耗一个自由度;三相把 Valid phases 设为 Vapor-Liquid-Liquid,引入第二液相变量并允许倾析器挂在任意级,自由度账需重算。
Q5:Aspen Plus 的 RadFrac 六种收敛方案该怎么选,各自的底层算法与初始化是什么?
A:Standard(inside-out 算法 + Standard 初始化)应首试,Strongly Non-ideal Liquid(Nonideal + Standard)用于强非理想两相,Petroleum/Wide-boiling(Sum-Rates + Standard)用于石油宽沸程,Azeotropic(Newton + Azeotropic)用于共沸精馏,Cryogenic(Standard + Cryogenic)用于深冷,Custom 供高级用户混搭算法与初始化。