如何使用ASPEN软件模拟完成精馏的设计和控制马后炮.docx
《如何使用ASPEN软件模拟完成精馏的设计和控制马后炮.docx》由会员分享，可在线阅读，更多相关《如何使用ASPEN软件模拟完成精馏的设计和控制马后炮.docx（25页珍藏版）》请在冰点文库上搜索。
如何使用ASPEN软件模拟完成精馏的设计和控制马后炮
ModifiedbyJEEPonDecember26th,2020.
如何使用ASPEN软件模拟完成精馏的设计和控制马后炮
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
第6章：
使用稳态计算选择控制结构
Steadt-stateCalculationsforControlStructureSelection
在我们转入将稳态模拟转化为动态模拟细节讨论之前，要先讨论一些重要的稳态模拟计算方法。
因为经常被用于精馏
设计中帮助为其选择一个实用且高效的控制结构，。
故此类讨论可能是一定意义的。
绝大部分精馏塔的设计是为了将两种关键组分分离获得指定的分离效果。
通常是两个设计自由度指定为馏出物中重关
键组分的浓度和塔底产品中轻关键组分的浓度。
因此，在精馏塔的操作和控制中，“理想的”控制结构需测定两股产品的组成
并操控两输入变量（如，回流流量和再沸器的输入热量），从而能够达到两股产品中关键组分的纯度要求。
然而，由于一些现实的原因，很少有精馏塔使用这种理想的控制结构。
组分检测仪通常购价昂贵且维修成本高，其可
靠性对连续在线控制而言，有时略显不足。
如果使用色层法，还会在控制回路中引入死时间。
此外，不使用直接测量组分
法，通常也有可能取得非常高效的控制效果。
温度测量被广泛应用于组分的推理控制。
温度传感器廉价而又可靠，在控制回路上只有很小的测量滞后。
对恒压二元
体系，温度与组成是一一对应相关的。
这在多组分体系中不适用，但精馏塔中合适位置的温度通常能够相当准确地提供关
于关键组分浓度的信息。
在单端控制结构中，只需控制某块塔板的温度；选择剩下的“控制自由度”时应使产品质量可变性最小。
例如，确定一定
的回流比RR或者固定回流与进料流量的比值R/F。
有时候，需要控制两个温度（双温控制系统）。
我们将在本章中讨论这
些被选方案。
如果选择使用塔板温度控制，那么问题便是选择最佳一块或数块塔板，该处的温度保持恒定。
在精馏文献中，这个问
题已讨论了半个世纪以上，且提出了一些可选择的方法。
我们将一一审视这些方法，并举例说明其在各个系统中的有效性。
需要重点关注的是，所有这些方法都仅使用稳态信息，因此，如AspenPlus之类的稳态过程模拟器可便捷地用于计算。
这些方法均要求恒定某些变量的同时将另一些变量变化。
例如，两股产品的组成或是某块塔板温度及回流流量恒定不变，
而进料组成变化。
在AspenPlus中，“DesignSpec/Vary”功能可以用来使期望的自变量恒定不变，计算所有其余应变量的
值。
在一些方法中，变化的变量是进料组成。
但对于任何一种方法，均不考虑进料流量。
这是因为进料流量的扰动可以直
接通过固定受控变量的流量与进料量的比值来处理。
当然，这需要假设整个塔的塔板效率固定不变。
同时，还需要假设每
个塔板的压力均不变。
这很少见，因为当气液流率变化时，塔板压降及塔板持液高度也会发生变化。
但是，这些影响均小
到不足以对控制系统造成很大的不利影响。
方法概要
斜率判据
满足斜率判据，关键在于选择相邻塔板之间温差最大的那块塔板。
绘制出在设计条件下的温度剖面图，研究剖面图的斜率，寻找斜率最大的那块塔板。
相邻塔板之间温度变化大，说明
该区域内重要成分的组成发生了变化。
控制此位置的塔板温度不变，则应该可以维持此精馏塔的组成剖面，防止轻组分流
向塔底、重组分窜入塔顶。
灵敏度判据
满足灵敏度判据的重点在于寻找由于一个受控变量的变化引起最大温度变化的那块塔板。
改变某一个受控变量（比如，回流流量），使其发生很小的变化（设计值的%）。
研究产生的塔板温度变化，观察哪
块塔板的温度变化最大。
对于其他受控变量（如再沸器热量输入），重复这一过程。
塔板温度的变化值除以受控变量的变化
值，就是这个塔板温度与此受控变量之间的开环稳态增益。
温度变化最大的塔板即是最“灵敏”的，故选择控制它。
增益较大，
说明此塔板的温度可以由相应的受控变量有效地控制。
增益较小说明阀门饱和态易于发生，且操作区域受到限制。
奇异值分解判据
Moore曾详尽地研究了稳态增益矩阵中奇异值分解（SingularValueDecomposition）问题。
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第2页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
在上节，我们计算了所有塔板温度和两个受控变量之间的稳态增益，从而形成了一个增益矩阵K。
它有NT行（塔板个
数）和2列（受控变量――回流比R、再沸器热量输入Q――的个数）。
利用标准奇异值分解程序（如Matlab中的svd函
数,svd（x,0）），将此矩阵分解为三个矩阵：
K=UσVT（分解的结果:
U为NT×2矩阵,σ为2×2矩阵,VT为2×2矩阵）。
以塔盘数
为座标轴,绘制两个矢量U1值和U2值的曲线。
对应着U量值最大的一块或数块塔板则指出了塔中最有效的控制位置。
σ
是2×2的对角矩阵，（其对角线元素）是K矩阵的奇异值。
（σ的对角元素中的）较大值和较小值的比值即为条件数，可以用来
评估二元温度控制方案的可行性。
条件数较大（或者最小奇异的值较小）时表明该系统难以控制。
调节器就是装置增益矩
阵的逆矩阵，假如奇异值为”0”则说明这是一个退化矩阵，不可转秩。
1.增益矩阵的概念:
按照控制理论,在多变量耦合控制系统中,选择其中第i个受控变量,当只有uj作用时,即只改变uj,使其他各受控变量uk
（k=1,2,3….n,k≠j）保持不变,当uj变化Δuj时,所得的控制变量yi的变化量与uj的变化量之比,称为uj到yi通道的第一放大系数,
Δy
也有称为开路增益的,K=|（u,k=1,2,3...,n,k≠j）此处稳态开路增益的概念与自控常用的相对增益矩阵的概念有区别.
i
ijΔu
j
2.奇异值分解的实际例子:
上述的描述过于抽象,现在按照Luyben先生1997年的”EssentialsofProcessControl”第458页的讲述,
试举例如下:
某9块塔盘的精馏塔用于分离水和异丙醇,对于其两个受控变量回流比R及热输入Q,各塔盘温度有如下的稳态增益:
塔盘编号ΔTn/ΔRΔTn/ΔQ
9
8
7
6
5
4
3
2
1
上表内各项即为该塔的稳态增益矩阵K的各元素,K矩阵有9行2列
ΔT
ΔT
=
.
.
ΔT
R
K
Q
对K做奇异值分解,得到以下3个矩阵:
K=UΣVT
-Σ=9VT=
-
U=
?
.34520
--
0
--
-
--
--
-
-
-
?
?
?
?
?
则矩阵U的第一列中最大的元素是,对应的塔盘数为6,因此,SVD法的结论是受控变量为回流量时,可以控制第
6块塔盘的温度.
必须指出,一般的对于一个m×n阶矩阵A的奇异值分解,得到的是一个m×m阶正交矩阵U,m×n阶对角阵Σ和n×n阶正交
矩阵V,这时所用的MATLAB命令是[U,S,V]=svd（X）,此处所用的MATLAB命令是[U,S,V]=svd（X,0）,即求”经济
尺寸”的奇异值分解,与常规的SVD得到的U,Σ,V的矩阵阶数并不相同.
详情请参见MATLAB主页:
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第3页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
3.条件数的概念:
在数值计算领域,条件数是与特定问题使用数字计算的可控制性（适定性）相关的一个度量参数.问题的
条件数较小,则称其为良态的,假如条件数较大则称为病态的.对于矩阵的条件数而言,并没有关于良态和病态的界定数
值.
在我们的讨论所涉及范围内,由条件数的定义推导出:
矩阵的条件数κ（A）,有:
σ（A）
κ（A）=
σ（A）
其中σmax（A）和σmin（A）分别是矩阵A奇异值的极大和极小值.
扩展阅读资料:
1.有关奇异值分解在增益矩阵上的应用，参见：
C.F.Moore,Selectionofcontrolledandmanipulatedvariables,in
PracticalDistillationControl,VanNostrand-Reinhold,1992,Chapter8.或Luyben的EssentialsofProcessControl,
（1997）第13章;或Luyben的PlantwideProcessControl（1998）第6章.
2.有关奇异值分解本身的线性代数定义、意义及其在数值计算中的作用及MATLAB中的基本使用，建议阅读：
《数值
分析与科学计算》，第2章等相关内容，JefferyJ.Leader着，清华大学出版社2008年5月第1版。
3.有关开路稳态增益、增益矩阵、相对增益矩阵及其解析解法的实际应用，增益值的对比和对实际控制结构选择上的
指导意义，参见《自动化仪表与过程控制，第四版》第8章及以下相关内容，作者：
施仁、刘文江等，电子工业出版
社，2009年2月第1版。
恒定温度判据判据
保持塔的馏出物和塔底产品的纯度不变，我们在预期的范围里改变进料组成值。
选择其温度不随进料组成变化而变化
那块塔板。
这个方法的困难性在于，对于进料组成的所有变化，可能并不存在保持恒定温度的塔板。
对于多组分体系，尤其如此，
非关键组分可能的改变会严重影响塔板温度，尤其在靠近塔两端处。
产品变化最小判据
满足产品变化最小判据的要点在于选择面对进料组成波动，塔盘的温度保持不变，引起的产品纯度变化最小的塔板。
选择几个候选塔板位置。
固定某一块塔板的温度，及第二控制自由度，如回流比或回流流量。
然后，在预期范围内改
变进料组成，计算产品的组成。
对于其他塔板位置，重复上述过程。
选择面对进料组成波动，其温度保持不变，引起产品
纯度变化最小的那块塔板。
小结
上述各节中，我们描述了五个最常用的判据。
有时候，这些判据推荐的控制塔板位置相同。
而在其他情况下，它们推
荐不同的控制塔板位置。
下一部分，我们将这些判据应用到几种典型的工业精馏体系中，评估其相对有效性。
丙烷/异丁烷二元体系
我们研究的第一个分离体系是丙烷和异丁烷的二元混合物。
进料流率为1Kmol/s，设计进料组成为丙烷摩尔分数40%。
按照习惯的符号标注，进料组成记为z，馏出物组成为xD，塔底产品组成为x（B均以丙烷的摩尔分率计）。
塔压设定在，
冷凝水可以用于冷却冷凝器（回流罐的温度为315K）。
这个塔使用Aspen的编号习惯给塔板编号，从回流罐直到塔釜有37
块理论板。
进料由第16块理论板引入。
馏出物纯度设定为摩尔分率为98%的丙烷，塔底产品的杂质含量设定为丙烷的摩尔分率为2%。
实现上述纯度所需的回
流为。
斜率判据
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第4页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
我们研究的第一个分离体系是丙烷和异丁烷的二元混合物。
进料流率为1Kmol/s，设计进料组成为丙烷摩尔分数40%。
按照习惯的符号标注，进料组成记为z，馏出物组成为xD，塔底产品组成为x（B均以丙烷的摩尔分率计）。
塔压设定在，
冷凝水可以用于冷却冷凝器（回流罐的温度为315K）。
这个塔使用Aspen的编号习惯给塔板编号，从回流罐直到塔釜有37
块理论板。
进料由第16块理论板引入。
馏出物纯度设定为摩尔分率为98%的丙烷，塔底产品的杂质含量设定为丙烷的摩尔分率为2%。
实现上述纯度所需的回
流为。
灵敏度判据
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第5页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
如图所示，上图给出了塔板温度和两个受控变量--回流流量R和再沸器热量输入QR--之间的开环回路增益。
实线
对应回流流量的变化，虚线对应再沸器热量输入的变化。
这两个输入信号在稳态值的基础上有很小（+%）的增加。
正如
所预料的，塔板温度针对回流之间的增益是负的，而塔板温度对再沸器热量输入之间的增益是正的。
两曲线说明了第8块理论板对回流流量的变化敏感，而第8块理论板和第29块理论板均对再沸器热量输入敏感。
因此，
第8块理论板的温度既可由回流流量控制也可由再沸器热量输入控制；而第29块理论板温度只能由再沸器热量输入控制。
应该记住的是，这些都是稳态结果，并没有反映任何动态信息。
所有塔板的温度由热量输入的变化而快速变化，因此将
热量输入和任何塔板的温度组对在动态上均是可行的。
然而，回流流量的变化需要相当长的时间去影响靠近塔釜塔板的温度，这是因为液相的流体滞后（每块塔板约3-6s）。
因此，将回流流量与离塔顶较远的塔板温度组对，可以预期只会出现糟糕的控制效果。
奇异值分解判据
如图下图所示，给出了奇异值分解（SVD）分析中U1和U2向量的值。
前者用实线表示，与回流有关。
后者用虚线表
示，与热量输入有关。
奇异值分解的结果类似于灵敏度的结果。
它们说明第8块理论板可由回流控制，而第29块理论板由热量输入控制。
稳
态增益矩阵的奇异值σ1=，σ2=，故条件数CN=σ1/σ2=。
这表明这两个温度彼此相当独立，因此二元温度控
制模式是可行的，至少从稳态角度如此。
温度恒定判据
Figure丙烷/异丁烷;保持产品纯度不变，由进料组成改变引起的温度剖面图的变化
图给出了在稳态值（丙烷40mol%）的正负两侧两个进料组成下温度剖面的变化。
实线代表了组成为35mol%丙烷
的进料，虚线则代表了组成为45mol%丙烷的进料。
对于两种进料组成，馏出物和塔底产品的组成分别设为98mol%和
2mol%。
正如所预期的，在二元恒压体系中恒定组成就确定了温度。
因此，塔顶和塔釜的温度均不发生变化。
理论上讲，控制
这些端部温度可以取得恒定的产品纯度。
然而，实际上，少量其他组分或者压力变化将导致使用端部温度控制效果不佳。
在讨论多组分体系时，这个结论将被证实。
产品变化最小判据
图给出了当某一特定塔板温度不变，进料组成变化时，产品纯度如何变化。
在这个图中，恒定的第二控制自由度是
回流流量。
表给出了选择恒定的回流流量而不是恒定回流比的理由。
这里用到在上一部分获得的信息。
产品纯度一定，进料组
成变化。
不用观察温度剖面图，我们注意到回流流量和回流比的所需要的变化。
正如表所清楚表明的，回流流量所需变
化比回流比所需变化要小的多。
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第6页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
Figure（a）丙烷/异丁烷：
保持回流和塔盘温度不变，由进料组成改变引起的馏出物纯度变化;（b）保持回流和塔盘温度
不变，由进料组成改变引起的塔底产品杂质含量变化.
表恒定产品纯度条件下,进料组分对回流量和回流比的影响
体系
产品纯度或杂质含量
（mol%）
z（mol%）
回流量
（kmol/s）
回流比
D/P98/235C3
40C3
45C3
苯/甲苯/二甲苯25B
30B
35B
多组分98/235C3/34IC4
40C3/29IC4
45C3/24IC4
醋酸甲酯25MeAc
30MeAc
33MeAc
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第7页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
因此，在图中，回流流量固定在s，固定一个塔板上的温度（理论板2号，8号，12号，20号，29
号，和37号）。
表中横坐标表示进料中丙烷的摩尔分率；纵坐标为馏出物的纯度xD，塔底产品的杂质含量xB。
图中显示了一些不合常理的结果。
控制靠近塔顶的第8块理论板的温度比控制靠近塔釜的第29块理论板温度能够
获得更好的塔底产品纯度。
塔底产品的杂质含量很接近、甚至低于预期的2mol%丙烷。
另一方面，控制靠近塔釜的第29
块理论板的温度能够获得更好的馏出物纯度，等于甚至高于预期的98mol%丙烷。
而传统知识认为，控制靠近产品物流的
塔板使得其纯度更恒定。
这些结果说明在这二元系统中，运用单回路温度控制方法，控制第8块和第29块理论板均能够获得相当好的产品纯
度。
如果使用二元温度控制方法，即控制精馏塔两端的温度，在压力恒定时，可以使产品纯度完全维持在稳态条件下预期
的值。
苯/甲苯/二甲苯三元体系
下一个的研究的分离体系是苯/甲苯/邻二甲苯三元混合物。
进料量1Kmol/s，进料组成（摩尔分数）设定为，苯30%，
甲苯40%，邻二甲苯40%。
操作目标是从重关键组分甲苯中分离出轻关键组分苯。
当然，比重关键组分甲苯更重的二甲苯
会随着甲苯进入塔底产品。
塔压为常压。
塔的理论板数为32块，原料由第16块理论板引入。
馏出物杂质含量为%（摩
尔分数）的甲苯，塔底产品的杂质含量为%（摩尔分数）的苯。
实现上述纯度所需的回流比为。
斜率判据
如图所示，上图给出了在设计条件下温度的剖面图，下图显示了邻近塔板之间温度的差异。
在进料板处有一个很大
的变化。
由于在该处引入进料，所以此处不适于温度控制。
在靠近塔底位置同样有一个大的温度变化，这是因为比重关键
组分更重的二甲苯在该处的聚集。
由于我们竭力探究苯和甲苯的组成，因此此处也不是温度控制的好位置。
斜率分析建议
第21块理论板作为温度控制的位置。
灵敏度判据
如图所示，上表给出了塔板温度和两个被控制变量之间的开环回路稳态增益。
这些曲线表明第21块理论板对热量
输入变化敏感，第22块理论板对回流流量变化敏感。
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第8页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
奇异值分解判据
如图所示，下表给出了奇异值分解分析中U1和U2的值。
前者由实线表示，与回流量相关。
后者由虚线表示，与
热量输入相关。
奇异值分解的结果与灵敏度的结果类似。
它们表明第21块理论板可由回流流量控制，而第23块理论板可由热量输入
控制。
稳态增益矩阵的奇异值σ1=，σ2=，故条件数CN=σ1/σ2=。
这表明这两个温度并不像丙烷/异丁烷体系
中的两个温度那样独立无关，因此二元温度控制体系可能并不有效。
这很容易理解，因为第21块理论板和第23块理论板
靠的太近以至于不能相互独立。
恒定温度判据
Figure苯/甲苯/二甲苯体系:
保持产品纯度不变时,由苯和二甲苯进料组成的变化引起的温度分布曲线变化
译者免责声明：
译者已经竭尽所能地确保译文正确完整地传达原作的意旨。
然而文中所论及的方法在工程中的具体使用，其使用责任完全在于使用人员。
本文仅为
学习了解所用，一切版权归于JohnWiley&Sons,Inc.请于下载后的24小时之内将此删除，译者不承担由此引起的一切法律责任。
第9页共18页
如何使用ASPENTM软件模拟完成精馏的设计和控制威廉·L·鲁平博士
图给出了这个三元体系中3个不同进料组成下的温度剖面图。
进料组成（摩尔分数）苯/甲苯/二甲苯（BTX）设定
为30%/30%/40%。
塔底产品和馏出物的杂质含量分别设定为苯%（摩尔分数），甲苯%（摩尔分数）。
实线表示BTX
进料组成（摩尔分数）为25%/35%/40%；虚线表示BTX进料组成（摩尔分数）为35%/25%/40%；点线虚线表示BTX
进料组成（摩尔分数）为25%/25%/50%。
对于进料中苯和甲苯比例的变化，得到的结果表明在产品纯度恒定时，第27块理论板的温度保持恒定。
因此，如果预
期的进料组成波动是这种类型，控制第27块理论板应该有效。
然而，对于进料中二甲苯浓度的变化，第27块理论板的温度几乎变化了3K。
Figure苯/甲苯/二甲苯体系:
保持回流和塔盘温度不变，由苯和二甲苯进料组成改变引起的馏出物纯度变化.
最小产品波动判据
图表明了当一块特定塔板的温度恒定，及进料组成变化时产品的纯度如何变化。
在本图中固定的第二控制自由
升级会员 