联立模块法在流程模拟中的应用
李玉刚周传光郑世清韩方煜
( 青岛化工学院计算机与化工研究所,  青岛266042)
摘要:  探讨了联立模块法在流程模拟中的应用, 提出了简化模型的建立和求解方法。求解用P 矩阵表达的各单元模块的简化模型时, 只需联立求解几
个基本混合点的物流, 其它的物流即可串联地解出。
关键词:  联立模块法; 流程模拟; 计算机
中图法分类号:  TQ  018
0 前言
随着计算机工业的发展和流程模拟软件的广泛应用, 流程模拟技术的研究也越来越深
入。目前开发应用的模拟方法有序贯模块法, 联立方程法和联立模块法  1 。其中序贯模块法是
把流程的各个单元表示成模块的形式, 并把模块的模型编成通用的子程序。在程序运行时, 一个模块接一个模块地调用。因此, 这种方法所需存贮量小, 查错方便, 可以根据需要添加程序块, 并且能够继承和有效地利用过去开发的软件包。这种方法在国内外开发的过程模拟软件中被广泛采用。但是当流程中存在循环流或者遇到设计型及优化型问题时, 应用这种方法需要反复迭代, 计算效率不高。联立方程法是将描述流程各单元的方程组合在一起, 联立求解这个大的非线性方程组。这样以来, 不论是模拟型问题还是设计型问题, 不论流程中有没有回路, 都无需不断地分级迭代, 一举可求得方程的全部解。计算效率高。然而目前求解大的非线性方程组还没有十分有效的方法, 致使联立方程法仍停留在原型试验阶段。联立模块法是介于以上两者之间的一种模拟方法。与序贯模块法相比, 在求解流程中含有
回路或者是设计型和优化型问题时, 具有更大的灵活性和较高的计算效率; 与联立方程法相比, 它只处理流程水平的简化方程组, 因而能处理更大规模的问题。尽管联立模块法的研究和应用发展很快, 国内外许多学者在从事这方面的研究工作, 但目前还未取得实质性的成果。该方法在大系统中的应用还很少报道, 其通用性和实用性需要进一步地探讨  2  。作者在前人工作
1 模型的建立与求解
联立模块法的思想最早是由R o sen 3    提出的, 其基本思想是: 流程中的每一单元都建立两
个模型, 一个是反映单元真实情况的严格模型, 另一个是根据严格模型提供的信息而形成的简化模型。
简化模型一般是线性的或者是非线性近似的工程模型。这样, 整个流程的模拟计算就可在两个不同的水平上进行, 即模块水平和流程水平。在模块水平上, 应用严格模型求出简化
收稿日期: 1996202202
18 青岛化工学院学报18 卷
模型的各参数。在流程水平上, 联立简化模型的线性方程和约束方程, 求出各单元上的输入、输出数据, 通过这两种模型的交替调用, 最后达到整个流程的收敛。整个收敛过程如图1 所示。
图1    联立模块法求解过程示意图
从图1 可以看出, 联立模块法同时需要建立严格模型和简化模型, 关于严格模型的建立和求解同序贯模块等方法相同, 可参考文献4, 5  。文中重点讨论简化模型的建立与求解。
2. 1 简化模型的建立
对流程中每一单元模块来说, 其进口流在一定的波动范围内, 它的出口流与进口流具有近似的线性关系。可表示为:  Y i  =  P i X 。
式中, X  为进料中各组份的流率向量, 对具有多股进料的单元模块, 可先经过一个虚拟的混合器合并为一股进料。Y i 为第i 股出料的各组份流率向量, P i 为第i 股出料的P 矩阵, 其中元素P i [ j ] 表示第j 个组份在该股出物流中的流量占进料量中其含量的分数。对一个单元模块有几个出口流就有几个P 矩阵, 在已知进口流的情况下, 各出口流就可由各自的P  矩阵求出。因此, 这些P 矩阵就构成了单元模块的简化模型。不同的单元模块具有不同的P  矩阵, 下面分别介绍。在以下公式中, n 表示组份数。
1. 1. 1 反应单元模块
假设反应单元只有一股进料和一股出料, 对于多股出料的情况可假设再经过一个虚拟的分离器进行分离。设反应单元中具有以下反应:
n
∑v i A i  =  0
i=  1
v i  为化学反应的计量系数, 对产物v i  >  0, 对反应物v i  <  0。其它组份v i  =
假设A K  为关键组份,  其转化率为ΓK , 则该反应单元的P 矩阵元素为:
0。
1
1 -  ΝK v iΝK ƒv K 0 i =  j ≠k
i =  j =  k
i ≠k,i =  1, 2  N
其它
(1)
P ij  =
1. 1. 2 分离单元模块
对一股进料, m  股出料的分离单元, 第i 股输出流的P 矩阵为:
Γi1
Γi2
P i j  = i =  1, 2, m(2)
Γin
19
第  1  期联立模块法在流程模拟中的应用
式中,  Γi j  为分配系数, 它们满足:
m
∑Γij  =  1j =  1, 2. . . n
i=  1
本模型适用于各种分离单元, 如分离塔、闪蒸器等。
1. 1. 3 分流器单元
分流器是分离器的特殊情况, 它是把输入流分成若干股输出流, 而各股输出流的组成与输入流的相同, 即
Γi1  =  Γi2  =
m
= Γin i =  1, 2, , m
∑Γij  =  1 J  =  1, 2 (3)
n
i=  1
1. 1. 4 混合器单元模块
liyugang混合器的作用是将多股流混合成一股流, 即
Y = IX 1  +  IX 2  +  IX 3  + . . . (4)其中I 是单位矩阵。
由上面讨论可知, 对反应单元只要知道关键组份的转化率即可代入( 1) 式求出该单元的P 矩阵。各反应单元模块的关键组份转化率可通过各自的严格模型求出。
对其它单元模块, 只要知道各进口物流的信息, 即可由各严格模型求出各出口流, 再由( 2~4) 式求出P 矩阵。
这样, 知道每个模块的P 矩阵, 即可建立起整个流程的简化模型。
1. 2 简化模型的求解
各个单元模块的P 矩阵求出以后, 在流程水平上就可以求解由P 矩阵表达的线性模型来代替严格模型。对这些线性模型的求解方法有两种, 一是联立求解这些线性模型和约束方程。这要求计算机贮存一个庞大的线性方程组, 还要采用稀疏矩阵的技术。另一种方法是保持原来的流程结构, 采用序贯模块法思想, 一个模块一个模块地进行计算, 这样虽然节省了内存, 但是需要劈分迭代技术。
为能既利用方程的线性本质, 又能利用流程网络结构的特性, 通过分析探索, 提出了一种新的求解简化模型方法。该方法的特点在于只需联立几个特殊点, 其它的点就可串联解出。如果将每股物流看作元素, 各个单元模块的P 矩阵看作是一个元素到另一个元素的弧( 映射) ,这样就可得到整个流程的信息拓朴图。在信息拓朴图中, 由连续的弧组成的弧列称为路。多于一个输入弧的点叫混合点。有些混合点, 它至少有路通向两个以上的混合点, 这类点称为基本混合点, 也就是需要联立求解的特殊点。求基本混合点的方法是, 首先去掉所有的进出物流(指整个系统的进料和出料)。然后将所有的简单路(没有分叉的路) 变成弧, 如果几条弧的起点和终点相同, 就合并为一条弧。这样剩下的元素就是基本混合点。这样只需联立求解这些基本混合点和某些约束条件, 其他点可串联解出, 从而使问题简化  2  。
2 应用实例
本文提出的方法已成功地应用于专用模拟软件“尿素流程模拟软件U R P S”6    和通用模拟软件“数据驱动的流程模拟软件DO P S”7    之中。在专用模拟软件U R P S 中, 由于P 矩阵法
20 青岛化工学院学报18 卷
循环回路和约束条件的收敛是采用简化模型实现的, 大大减少了复杂严格模型的调用次数, 运行速度显著提高。文中通过选用不同的初值和收敛判据对这两种方法进行了比较, 达到同样的精度要求和计算结果, 前者比后者~4 倍, 比较结果见表1。
表1 U R P S 中两种方法的比较结果
方法调用严格模型次数耗时(A ST ƒ386 微机)
序贯模块法联立模块法30~50
10~20
40~60 m in
10~20 m in
在通用模拟软件DO P S 中, 由于本方法的应用, 大大缩短了模拟和优化时间。尤其是在进行优化计算时, 需要几百次甚至上千次地调用各单元模块, 采用P 矩阵表达的简化模型时, 其效果更加突出, 一般能缩短计算时间10 倍以上。
实践证明, 本文采用的以P 矩阵表达简化模型的方法, 具有较好的实用性, 只要各单元模块的P 矩阵对进口流的变化不十分敏感, 就会取得非常好的效果。另外, 我们还发现, 在大多数情况下, 在计算开始, 只要用严格模型循环3~5 次, 各单元模块的P 矩阵就基本稳定, 可以启用简化模型进行计算。另外本文对简化模型求解的新方法, 充分利用了问题的线性本质和流程网络结构的特点, 速度快占内存少, 适于在微机上实现。
参考文献
1 2 3 4 5 6 7 朱建业, 李有润, 彭秉璞. 联立模块法研究进展. 化工进展, 1989, ( 1)
李玉刚. 青岛化工学院硕士论文. 1992
R o sen E M . C h en E ng P r o g , 1962, 58: 10
韩方煜编. 过程工程的模型化和模拟计算. 青岛化工学院讲义. 1985
杨友麒编. 化工系统工程. 北京: 化学工业出版社, 1989
青岛化工学院. 尿素生产流程模拟软件U R P S. 技术手册. 青岛化工学院. 1993 青岛化工学院. 数据驱动的流程模拟软件DO P S. 技术手册. 1995
A pp l ica t i on of S i m ul tan e ous-M odula r A pproa c h  in F l ow s hee t i n g
L i Y u g a n g Z h o u C h u a n g u a n g Z h e n g S h iq in g H a n  F a n gy u
(R e sea rch  I n st i tu te o f Com p u te r and C h e m ica l E ng inee r i ng, Q ingdao  I n st itu te o f C h e m ica l T ech no lo gy, Q ingdao 266042) )
A bstra c t:  T h e app li ca t i o n o f s i m u ltan e o u s2m o d u l a r app ro a ch i n f l ow sh e e t i n g is s tu d i ed.  It a r ise s how  to  s i m p li fy  th e m o de l b u ild i n g an d  th e so lu t i o n .  A  p rac t i ca l m e tho d o f  so l v i n g li n ea r equ a t i o n s is p re sen ted, an d  th is m e tho d can  save m u ch com p u te r m em o ry an d ca l cu a2 t i o n t i m e.
Key W ord s:  f l ow sh e e t i n g; s i m u ltan e o u s2m o d u l a r app ro a ch