(能源化工行业)化工过程分析与合成作业

（能源化工行业）化工过程分析与合成作业

化工过程分析和合成

（三级闪蒸的模拟）

组员：0815010115尹小莉

0815010114杨秋静

0815010113张志莉

0815010119罗浩

0815010130李家翔

0815010131杜秦

壹、序贯模块法简介：

序贯模块法的基本原理

序贯模块法的基础是单元模块（子程序），通常单元模块和过程单元是壹壹对应的。过程单元的输入物流变量即为单元模块的输入，单元模块的输出即为过程单元的输出物流变量，单元模块是依据相应过程单元的数学模型和求解算法编制而成的子程序。单元模块具有单向性特点，给定其输入物流变量及参数可计算出相应的输出物流变量，但不能由输出变量计算输入变量，也不能由输入、输出变量计算模块参数。

序贯模块法的基本思想

从系统入口物流开始，经过接受该物流变量的单元模块的计算得到输出物流变量，这个输出物流变量就是下壹个相邻单元的输入物流变量。依次逐个的计算过程系统中的各个单元，最终计算出系统的输出物流。

3.序贯模块法的优点：

和实际过程的直观联系强；模拟系统软件的建立、维护和扩充都很方便，易于通用化；计算出错时易于诊断出错位置。其主要缺点是计算效率较低，尤其是解决设计和优化问题时计算效率更低。计算效率低是由于序贯模块法本身的特点所决定的。对于单元模块来说，信息的流动方向是固定的，只能根据模块的输入物流信息计算输出物流信息，而且在进行系统模拟过程中，对物料、单元模块计算、断裂物流收敛计算等，将进行三重嵌套迭代。虽然如此，序贯模块法仍不失为壹种优秀的方法。但在处理过程中设计和优化问题时，由于其循环迭代嵌套甚至可高达五层，以至其求解效率就太低了。

二、三级闪蒸模拟

1、三级闪蒸过程如下图：

图1

图中入Array料流量

为

453.6mol

/h，入料

组成为

丁烷

30%，戊

烷40%,

己烷

30%。入

料温度

为

121.1℃，压力位1723.7kPa。三个闪蒸器的压力均为709.5kPa。当闪蒸器温度分别为以下值时

闪蒸器1:106.9℃

闪蒸器2:98.9℃

闪蒸器3:114℃

闪蒸器1:107.2℃

闪蒸器2:94.0℃

闪蒸器3:119.2℃

分别用直接迭代法和阻尼直接迭代法计算气相和液相产品的流量和组成。阻尼因子分别取值为0.5,0.3，-0.2，-0.3，-0.7，-0.9。

依据闪蒸条件设该闪蒸过程为理想体系，三个闪蒸器均为等温闪蒸过程，建成相应的单元模块且将图1绘成图2的三级闪蒸过程模拟模块流程。

按图2所示的模拟流程，对俩组闪蒸条件进行模拟计算

三级闪蒸过程

图2

b.

从上例可见，阻尼因子q值的选取具有较大的任意性和经验性。1958年Wegstien提出了壹种简便的方法，能够弥补这种阻尼因子取值困难的弱点。

用联立模块法对三级闪蒸过程进行稳态模拟

建立简化模型

严格单元模块的输入流股变量向量x和输出流股变量y之间有严格模型：

y=G(x)

上式的壹阶台劳展开式为：

即） 令，，，便可得到严格模型的线性增量简化模型： 利用上式分别对每个过程单元写出其简化模型： 混合器 闪蒸器1：

闪蒸器2： 闪蒸器3：

由于混合器的严格模型为线性模型，且系统入料流量为给定值，所以有：

把上述线性简化模型写成矩阵形式的迭代格式，得到壹个稀疏线性方程组：

显而易见，上式是壹个维数不高的稀疏线性方程组。 从严格模块计算简化模型的系数

上式的系数矩阵可通过对严格模块的扰动计算得到。前面我们假定:,也就得到了A 。而A 是从壹阶台劳展开式得到的。偏离点后便会产生偏差，因此要不断进行修正。

分别对每个单元建立简化模型，然后把单元简化模型、联结方程、设计规定方程集合到壹起组成过程系统的简化模型，由于切断了全部联结物流，描述整个过程的系统的简化模型方程数为：

如果不考虑设计规定方程，用线性增量模型作为单元的简化模型，则过程系统的简化模型如下：

单元模型: 联结方程：

把上式写成矩阵形式：

对于壹个包括100各联结流股，每个流股有8各组分，十个设计规定系统，其系统简化模型数为： 由此可见对于较大的系统简化方程数很大，消去物流联结方程，从而使简化模型的维数大大减少： 此时简化模型方程数为： 三、以不可分割子系统为基本单位建立简化模型，回路切割方式 （1）将带有循环回路的虚拟单元作为简化模型的基本单位。 （2）虚拟单元所包含的各单元间的联结流股变量不出当下简化模型中，大大降低了简化模

型的维数。 （3）简化方程组数为： （4）壹些简单的方程组求解技术就能够处理这样的流程模型。 用回路切割模式建立三级闪蒸过程的线性增量简化模型。 最佳切割物流是S2

简化线性模型： 用矩阵表示为：

闪蒸模拟的面向方程法

面向方程法简介

把描述过程的所有数学模型汇集在壹起，形成壹个非线性方程进行求解，即： F(X,W)=0

式中：X------状态变量向量

2

332S A S ∆=∆4

664S A S ∆=∆I

A A A ===2126256

52S I S I S ∆+∆=∆0876243214321ϖ=⎥⎥⎥⎥⎥

⎥⎥⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢

⎢⎣⎡∆∆∆∆∆∆∆∆⎥⎥⎥⎥

⎥⎥⎥⎥⎥⎥⎦⎤

⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡--------k

k x x x x x x x x I I I I I I I I I A I

A I A I A ()u N i i eq N C M st ++=∑

=1

322

882S A S ∆=∆*

S S 22∆=∆