第8章-过程系统模拟的基本方法ppt课件

《第8章-过程系统模拟的基本方法ppt课件》由会员分享，可在线阅读，更多相关《第8章-过程系统模拟的基本方法ppt课件（89页珍藏版）》请在装配图网上搜索。

1、第3章 化工过程模拟的基本方法第3章 化工过程模拟的基本方法1主要内容过程单元过程单元过程单元的自由度过程单元的自由度过程单元的模型化与模拟过程单元的模型化与模拟过程系统过程系统过程系统的自由度过程系统的自由度系统结构的识别系统结构的识别系统的分解系统的分解序贯模块模拟法序贯模块模拟法原理原理寻找最佳断裂流股寻找最佳断裂流股断裂流股的收敛断裂流股的收敛面向方程模拟法及联立模块模拟法原理面向方程模拟法及联立模块模拟法原理过程系统的优化过程系统的优化常用商业化流程模拟软件简介常用商业化流程模拟软件简介主要内容过程单元2公司结构常减压催化焦化重整芳烃原油原油产品产品产品产品炼油厂换热器车间塔器车间钢

2、材钢材产品产品产品产品机械厂石化公司不相关子系统不可分隔子系统回路系统公司结构常减压催化焦化重整芳烃原油产品产品炼油厂换热器塔器钢3过程系统模拟一般步骤：确定需模拟的系统识别不相关子系统（一级子系统，各子系统独立处理）一级子系统自由度分析找到分解成规模最小的二级子系统的决策变量各一级子系统内不可分隔子系统的识别（二级子系统，各二级子系统可依次求解）各二级子系统依次处理找回路、确定最佳断裂位置、确定计算顺序相关单元操作建模，确定迭代方法过程系统模拟一般步骤：确定需模拟的系统4第1节 流程的自由度分析1.一些基本概念过程(Process)：对原料进行某些物理或化学变换，使其性质发生预期的变化机械加

3、工不能称为过程系统：由相互联系，相互作用的若干组成部分结合成的具有特定功能的总机体过程系统：由各种过程构成的系统第1节 流程的自由度分析1.一些基本概念5流程:描述化工生产的物料流向及能量流向及装置特点的过程REACTORFEEDRECYCLEREAC-OUTCOOLCOOL-OUTSEPPRODUCT流程:REACTORFEEDRECYCLEREAC-OUTC6模型：复杂的A简单的B来替代。研究B来预测A的行为 实物模型 /数学模型模拟：对某一描述实际过程的数学模型利用数学方法进行求解，并对结果作出解释模型：复杂的A简单的B来替代。研究B来预测A的行为72.物流自由度、单元自由度及系统自由度

4、物流(stream)Dhem定理：对于一个已知每个组分初始质量的封闭体系，其平衡状态完全取决于两个独立变量，而不论该体系有多少个相，多少个组分或多少个化学反应 2.物流自由度、单元自由度及系统自由度物流(stream)8过程单元自由度分析方法过程单元自由度:可改变单元操作状态的独立变量数目列方程法：列出现象方程、限制方程等 自由度=变量数-方程数l描述规则法：要完全描述一单元设备的操作，必须确定的自由度的数目必定等于能由设备结构确定或能用外部手段控制的变量的数目l公式分析法：过程单元自由度分析方法过程单元自由度:9常见过程单元自由度单元名称自由度数常规指定变量混合器0分流器S-1流量分配比闪蒸

5、器2闪蒸温度、压力泵、节流阀1出口压力压缩机/膨胀机2绝热多变效率，出口压力换热器1某一物流出口温度反应器2+r反应程度，绝热，压降常规精馏塔5塔板数，进料位置，R,B,操作压力平衡级3侧线抽出比，压力，换热量常见过程单元自由度单元名称自由度数常规指定变量混合器0分流器10过程系统自由度过程系统自由度:确定系统状态的独立变量数目 系统自由度=单元自由度+进料自由度过程系统自由度过程系统自由度:11第节 过程系统结构的计算机识别主要内容 1.过程系统结构有向图2.过程系统结构的矩阵表示3.系统结构的识别第节 过程系统结构的计算机识别主要内容121.过程系统结构有向图节点边1.过程系统结构有向图节

6、点边132.系统结构的矩阵表示节点-节点节点相邻矩阵节点-边关联矩阵边-边弧相邻矩阵2.系统结构的矩阵表示节点-节点14节点相邻矩阵相邻矩阵定义为：式中，m为节点数目；sij为矩阵元素，定义为节点相邻矩阵相邻矩阵定义为：15163.系统结构的识别可及矩阵法不相关子系统不相关子系统过程系统：建模过程中可分别独立处理方程系统：写出方程组事件矩阵，Himmelblau算法识别不可分隔子系统不可分隔子系统可及矩阵法，索引矩阵法、图解法，Steward通路法3.系统结构的识别可及矩阵法不相关子系统17方程系统识别-Himmelblau算法：在mm事件矩阵M中，选出非零元素最多的列k。保留M中k列内每个零

7、元素对应的行，k列中为1的元素所对应的行用布尔加法合并成一行排列在最后。得到的新的jm的布尔矩阵记做M(0)；重复，从而得到序列M，M(0)，M(N)；最终得到矩阵M(N)，其每一列只有一个非零元素，其每一行与原方程系统中的不相关子系统对应。方程系统识别-Himmelblau算法：在mm事件18例例19过程系统识别不相关子系统分析不可分隔子系统可及矩阵过程系统识别不相关子系统20公司结构常减压催化焦化重整芳烃原油原油产品产品产品产品炼油厂换热器车间塔器车间钢材钢材产品产品产品产品机械厂石化公司不相关子系统不可分隔子系统回路系统公司结构常减压催化焦化重整芳烃原油产品产品炼油厂换热器塔器钢21可及

8、矩阵法识别不相关子系统的准则：可及矩阵法识别不相关子系统的准则：矩阵A的可及矩阵A*定义为：同属一个回路的节点当满足该条件的节点数量大于时，说明坐在公用一个节点的复合回路可及矩阵法识别不相关子系统的准则：矩阵A的可及矩阵A*定义为22可及矩阵法识别不可分隔子系统依据：可及矩阵法识别不可分隔子系统依据：Berge定理若用A表示某有向图的节点相邻矩阵，那么，矩阵H=A 中为1的元素hij表示从节点i沿弧的正方向经段弧可以到达节点j可及矩阵法识别不可分隔子系统依据：Berge定理23矩阵运算规则：矩阵运算规则：矩阵间：阵代数规则矩阵元素：布尔代数原则 布尔乘法布尔加法矩阵运算规则：矩阵间：阵代数规则

9、布尔乘法24例例25系统网络图 系统网络图 26例求解顺序：H-ABCDE-FG-I例求解顺序：H-ABCDE-FG-I27第第3 3节节 序贯模块模拟法序贯模块模拟法主要内容：基本原理循环流股的断裂与迭代断裂变量的收敛第3节 序贯模块模拟法主要内容：28一、基本原理一、基本原理混合器Submix(F1,F2,P1)P1=F1+F2EndSub分割器SubSPLT(F1,P1,P2,ALFA)P1=F1*ALFAP2=F1*(1-ALFA)EndSub收敛单元S4S4单元子程序内容：基础数据 单元方程组 求解算法单元子程序功能：根据给定的进料 条件和设备参数，计算出单元输出 结果(预测型)一、

10、基本原理混合器 收敛单元S4S4单元子程序内容：29设S1=10kmol/h，分割比=0.5,进行流程模拟计算解 设S4=0，进行MIX的模拟计算S2=S1+S4=10+0=0 进行分割器的模拟计算S4=S2=0.510=5 比较S4与S4 现假设S4=10，由MIX模块计算得到S2=S1+S4=10+10=20 进行分割器的模拟计算S4=S2=0.520=10 计算得S4与假设S4的数值相等，假设正确。由SPLT模块计算得S3=10。流程计算完成。设S1=10kmol/h，分割比=0.5,进行流程模拟计算30序贯模块法模块的特点 单向性单元操作输入输出设计规定流程计算收敛单元-循环流/设计规

11、定积木式序贯模块法模块的特点 单向性单元操作输入输出设计规定流程计算31断裂位置的影响：如果如果S2S2与与S4S4的自由度不同，需迭代的变量数也将不同的自由度不同，需迭代的变量数也将不同断裂位置的影响：如果S2与S4的自由度不同，需迭代的变量数也32二、再循环流股的断裂二、再循环流股的断裂方法1方法2原则：原则：将所有闭合回路全部打开将所有闭合回路全部打开二、再循环流股的断裂方法1方法2原则：将所有闭合回路全部打开332.1 最优断裂准则(1)被切断的流股数最少；(2)被切断的流股变量数最少；(3)被切断的流股的权重因子之和最少；(4)回路切断的总次数最少。通常选择原则：满足（满足（4）的基

12、础上，选（）的基础上，选（1）或（）或（2）2.1 最优断裂准则(1)被切断的流股数最少；34抉择依据：最少计算时间计算时间：计算时间：断裂方式；流程及变量灵敏度有效计算时间：流程计算；断裂流股迭代 迭代时间：迭代次数；收敛速度抉择依据：最少计算时间计算时间：352.2 回路矩阵 回路A:单元S2IIIS4回路B:单元IS1S2IIIS5I回路C:单元IS1S2IIIS3S6I回路D:单元S2IIIS3S72.2 回路矩阵 回路A:单元S2IIIS436回路矩阵 矩阵元素aij定义为回路矩阵 矩阵元素aij定义为372.3 UpadhyeGrens断裂法(II)1)有关术语 有效断裂组有效断裂

13、组能够把全部简单回路至少切断一次的断裂流股的集合。S2,S1,S3,S4，S1,S2,S5 S1,S3 多余断裂组多余断裂组若从一个有效断裂组中至少可以除去一个流股，而且得到的断裂组仍为有效断裂组；或者存在着对一个回路的二次断裂。S1,S2,S5 S1,S3,S4 2.3 UpadhyeGrens断裂法(II)1)有关38非多余断裂组非多余断裂组除去多余断裂组以后的有效断裂组。断裂族断裂族具有相同计算顺序的有效断裂组的集合。非多余断裂组除去多余断裂组以后的有效断裂组。392)、替代规则令D1为一有效断裂组，Ai为全部输入流股均属于D1的单元，将Ai的所有输入流用Ai的全部输出流替代，构成新的断

14、裂组D2，则D2也是有效断裂组对直接迭代，D2与D1具有相同的收敛性质 由于全部回路被打开，至少会有一个单元的全部输入条件由于全部回路被打开，至少会有一个单元的全部输入条件已知，可以作为计算的起点，计算出其输出流股。而该流股已知，可以作为计算的起点，计算出其输出流股。而该流股又将作为后续单元的输入流股，产生连锁效应！又将作为后续单元的输入流股，产生连锁效应！2)、替代规则 令D1为一有效断裂组，Ai为全部输403)、Westerberg算法从任何一有效断裂组开始，运用替代规则如果在任何一步中出现重复断裂组，则消去其中的重复流股，消去重复流股后形成的新断裂组作为新的起点重复、，直到没有重复断裂组

15、出现，且每个树枝上的断裂组重复出现为止，从最后一个新的起点开始，其后出现的所有不重复的断裂组构成非多余断裂族非多余断裂族中总数最小的断裂组为最优断裂组3)、Westerberg算法从任何一有效断裂组开始，运用41例 例 42非多余断裂组非多余断裂组断裂变量数S29S1,S4,S72+3+2=7S3,S4,S52+3+3=8S4,S5,S6,S73+3+4+2=12准则1准则2全部满足准则4：每个回路都断裂，且只断裂每个回路都断裂，且只断裂1次次非多余断裂组非多余断裂组断裂变量数S29S1,S4,43三、断裂流股变量的收敛 三、断裂流股变量的收敛 44收敛单元功能 提供循环流x的初值x。根据初值

16、x，及其它输入条件进行流程计算，求出循环流出口值y，比较x和y，若x与y之差满足精度要求则停止计算，否则转。根据x与y值以一定的方式产生x的新估计值x1，转 收敛单元功能 提供循环流x的初值x。45收敛单元的要求 对初值要求不高 易得，不易引起迭计算的发散提供的初值组数少例:初值的影响收敛单元的要求 对初值要求不高 易得，不易引起迭计算的发46数值稳定性好数值稳定性单调收敛数值稳定性好数值稳定性单调收敛47续数值稳定性振荡收敛续数值稳定性振荡收敛48续数值稳定性单调发散续数值稳定性单调发散49续数值稳定性振荡发散续数值稳定性振荡发散50收敛速度快 影响收敛的速度主要因素：迭代次数 流程计算次数

17、(Pass)迭代中矩阵求逆次数收敛速度快 影响收敛的速度主要因素：51占内存小 占内存小 52四、常用迭代法 直接迭代法有界Wegstein法Broyden法主特征值法联立超松弛法 Newton-Raphon法Marquardt法四、常用迭代法 直接迭代法53阻尼直接迭代法 q取值的一般原则：稳定性较好流程，取小于零，加快收敛速度如合成氨过程取q=-0.75。收敛速度快稳定性差，可取01，以改善稳定性如反应器q=0.5。阻尼直接迭代法 q取值的一般原则：54几种收敛策略的比较 方法偏导数函数计算次数稳定性速度直接迭代法否1较好慢Newton-Raphon法是n较差快Broyden法否1较好适中

18、Marquardt法是n好适中Wegstein法否1较好适中主特征值法否1较好适中几种收敛策略的比较 方法偏导数函数稳定性速度直接迭代法否1较55收敛判据 绝对误差 相对误差 收敛判据 绝对误差 56五、解决设计问题 估计反应单元的温度为估计反应单元的温度为T T估计再循环物流估计再循环物流S4S4依次计算混合单元、反应单元、分离单元，得到新的依次计算混合单元、反应单元、分离单元，得到新的S4S4的的比较比较S4S4与与S4S4，若两者相等则进行下一步，若不相等则返回，若两者相等则进行下一步，若不相等则返回 在在收收敛敛单单元元内内比比较较S5S5和和设设计计值值，若若两两者者不不相相等等则则

19、返返回回，若若相等则计算结束相等则计算结束五、解决设计问题 估计反应单元的温度为T57同时收敛策略 控制模块的设置增加了迭代循环圈，导致计算量的增加D设计规定向量 H过程系统方程组p决策变量与系统参数向量同时收敛策略 控制模块的设置增加了迭代循环圈，导致计算量的增58第4节 常用流程模拟软件简介主要内容：流程模拟技术的进展流程模拟软件的结构常用软件介绍第4节 常用流程模拟软件简介主要内容：59一、流股模拟技术的进展流股模拟技术序贯模块法联立方程法（面向方程法）联立模块法一、流股模拟技术的进展流股模拟技术601 联立方程法原理Sub flow(S1,ALFA,S2,S3,S4)S2=S1+S4

20、混合器模型 S3=S2*ALFA S4=S2*(1 ALFA)分割器模型End Sub 特点：特点：单元模型只提供方程，系统方程单元模型只提供方程，系统方程组一起求解组一起求解无循环物流及设计规定的特殊处无循环物流及设计规定的特殊处理问题理问题1 联立方程法原理Sub flow(S1,ALFA,S611.1模拟结构 循环物流方程设计规定方程仅提供方程仅提供关联式多!特殊处理1.1模拟结构 循环物流方程仅提供方程仅提供关联式多!特621.2 物性系统的处理占用7080%的计算量处理方法：1）不提供方程组，仅提供子程序2）提供 K,H 的方程组，其它物性由子程序调用1.2 物性系统的处理占用708

21、0%的计算量63 1.3 序贯法与联立方程法的对比 1.3 序贯法与联立方程法的对比 641.3.1 序贯模块法优点：与实际过程的直观联系强，软件的建立、维护和扩充很方便；易于通用化；需要的计算机内存较小；易于诊断出错位置；缺点：循环物流及设计问题：慢1.3.1 序贯模块法优点：缺点：循环物流及设计问题651.3.2 面向方程法优点 便于实际用户的提出的质量性能要求。有设计规定的系统，仅增加几个方程 用空间换取时间。相同的过程单元出现多次，将占用更多的 存储空间，自由度的选择也可能各不相同 以空间换取求解的数值稳定性。有较多回流高度交互作用流程相对容易收敛1.3.2 面向方程法优点661.3.

22、2 面向方程法缺点 要求提供较好的初值 难于自动产生流程的模型方程组，难以通用化 不容易诊断出现错误的地方。1.3.2 面向方程法缺点672.联立模块法联立模块法又称双层法，将过程系统的近似模型方程与单元模块交替求解联立模块法兼有序贯模块法和面向方程法的优点。既能使用序贯模块法积累的大量模块，又能将最费计算时间的流程收敛和设计约束收敛等迭代循环合并处理，通过联立求解达到同时收敛2.联立模块法联立模块法又称双层法，将过程系统的近似模型方68第8章-过程系统模拟的基本方法ppt课件69过程系统稳态模拟三种方法的比较过程系统稳态模拟三种方法的比较过程系统稳态模拟三种方法的比较70二、过程系统优化模型

23、的建立过程模拟模型给定参数可调参数经济模拟模型经济参数最优化程序过程约束条件经济约束条件经济特性指标二、过程系统优化模型的建立过程模拟模型给定参数可调参数经济71三、过程系统模拟软件的结构流程模拟软件类型通用软件专用软件三、过程系统模拟软件的结构流程模拟软件类型721.流程模拟系统结构1.流程模拟系统结构732.模拟步骤2.模拟步骤74 输入流程拓朴结构信息 输入单元计算所需的必要信息，选择物性及热力学性质，传递性质计算方法 进行全流程的模拟计算 解释、考核计算结果单元名称进料混合器1反应器精馏塔混合器2分流器产品输入物流号1，7234，856，9输出物流号1，8234，596，7 输入流程拓

24、朴结构信息单元名称进料混合器1反应器精馏塔混合75过程系统模拟一般步骤：确定需模拟的系统识别不相关子系统（一级子系统，各子系统独立处理）一级子系统自由度分析找到分解成规模最小的二级子系统的决策变量各一级子系统内不可分隔子系统的识别（二级子系统，各二级子系统可依次求解）各二级子系统依次处理找回路、确定最佳断裂位置、确定计算顺序相关单元操作建模，确定迭代方法过程系统模拟一般步骤：确定需模拟的系统76四、典型流程模拟软件1.著名通用流程模拟软件平台及适用范围20世纪50年代中期开始流程模拟软件的研制和开发1958年美国AI.W.Kellogg公司推出了世界上第一个化工模拟程序Flexible Flo

25、wsheeting在80年代，化工过程模拟进入成熟期，模拟软件的开发、研制走向了专业化、商业化到了21世纪初，流程模拟商业化软件更是到达了一个比较成熟的阶段，涌现出了大批的模拟软件四、典型流程模拟软件1.著名通用流程模拟软件平台及适用范围77常见的化工模拟软件ASPEN PLUS、HYSYS、PRO/IIgPROMS、CHEMCAD、VMGSimDesign II、ECSS 常见的化工模拟软件ASPEN PLUS、HYSYS、PRO/78ASPEN PLUS大型通用流程模拟系统，美国AspenTech公司出品全球各大化工、石化、炼油等过程工业制造企业及著名的工程公司都是Aspen Plus的用

26、户。在实际应用中，Aspen Plus 可以帮助工程师进行工艺过程的能量和质量平衡计算，用于装置标定；预测物流的流率、组成和性质；预测操作条件、设备尺寸；设计一个新的工艺过程、查找一个原油加工装置的故障或者优化一个乙烯全装置的操作等工程和操作的关键问题官方网站：http:/www.AspenT PLUS大型通用流程模拟系统，美国AspenTe79HYSYShysys是hyprotech公司的软件，原叫hysim，95年推出hysys，2002年Aspen收购了hyprotech公司全部股权，但2004年被美国联邦商务委员会判定为垄断经营。Honeywell和Aspen达成最终协议，由Honey

27、well收购hysys建模软件知识产权和操作培训仿真ots业务。条款包括：1、Aspen公司将保留hyprotech产品许可，以继续实施推广、销售和开发hyprotech的权利。2、Aspen公司在三年内不参与操作员培训服务竞争。3、两年支持协议，对转让的hyprotech产品，Aspen公司将继续为Honeywell公司提供新版产品的源代码。07年后Honeywell以Hysys源代码为基础开发UniSimHYSYShysys是hyprotech公司的软件，原叫hy80HYSYS可以利用HYSYS对工业装置进行安全分析和预测；进行操作规律和控制方案的研究及连锁控制调试；分析装置操作、生产过程

28、中的瓶颈问题；确定安全的开工方案；计算间歇生产过程的安全问题；计算特殊的非稳态过程；安全生产指导和调优；在线优化及先进控制；生产培训等。官方网站：http:/www.AspenT Process Modeling System）是英国PSE公司开发的通用工艺过程模拟系统特点是可建造任何反应过程、分离过程、多个过程的组合特别适用于任何新的工艺过程的研究开发，因为它可以做研发过程的实验设计，过程建模及参数的估值、数据的统计分析，如数据校验、一致性检验、残差分析等。由于它擅长于求解对时空变量的偏微分方程组，所以特别适用于动态过程的建模便于建立仿真培训系统，进行生产装置开工过程的培训、指导、控制，它可

29、以是在线的，也可以是离线的，并可以对动态和稳态过程进行优化控制已广泛用于化学工业、石油化工、石油和天然气加工、造纸、精细化工、食品工业、制药及生物制品技术等加工行业。官方网站：http:/ Process M83ChemCAD美国Chemstations公司开发的化工流程模拟软件。可以在计算机上建立与现场装置吻合的数据模型，并通过运算模拟装置的稳态或动态运行，为工艺开发、工程设计、优化操作和技术改造提供理论指导。官方网站：http:/ II美国WinSim Inc.公司开发的流程模拟软件经过近30年的开发和改进，Design II已经成为流程模拟变革的先驱；许多Design II的革新，如在线

30、Fortran和严格塔计算，均已确立了流程模拟的标准。官方网站：http:/www.WinS II美国WinSim Inc.公司开发的流程85VMGSim 加拿大Virtual Materials Group公司的产品VMGSim的核心成员是HYSIM/HYSYS的原始开发人员VMGSim 作为计算准确，功能强大以及高性价比的稳态流程模拟软件，可以详细预测工艺装置和工厂的性能，有效的帮助工程师达到提高操作效率，改善产品质量，实现节省投资，降低操作费用，提高效益，安全生产的目标。官方网站：http:/www.VMGS http:/ 加拿大Virtual Materials Gr86软件名称件名称

31、gPROMSAspen PlusHysysPRO/II历史背景史背景1997由英国帝国理工学院开发成功1981年在美国麻省理工学院MIT开发成功加拿大HyproTech从80年代开发SimulationScience公司开发(70年代)单元操作模型元操作模型库三相闪蒸多级分离塔模型反应器模型固体处理操作有8种8种可以有8种7种可以有可以自由组合5种/7种6种有物性物性计算系算系统纯化合物数目VLE双元交互系数电解质系统可用gPROMS、OLI、PC-SAFT、DIPPR、AspenPropertiesPlus等CAPE-OPEN物性数据库5941个40,000组有2000多种18,000200

32、0多种3000多组有特种工特种工艺系系统模型模型GLC气液接触模型、固定床反应器模型、溶液结晶模型、聚合过程模型、燃料电池模型多种PEP工艺模型（SRI的工艺包）与SPYRO软件接口炼油专用模块RefSYS,包括催化裂化，催化重整反应，异构化反应器，烷基化反应器，加氢裂化反应器，加氢精制/处理反应器炼油厂各种典型装置模型优化化计算算有有有有经济核算核算有有/有部分动态模模拟计算算均为动态(稳态时间方程为零)需AspenDynamics连接可直接转入动态模拟与Dynsim连接进中国市中国市场时间2002年80年代初90年代80年代末常见流程模拟软件对比软件名称gPROMSAspen PlusHysysPRO/I87第8章-过程系统模拟的基本方法ppt课件88第8章-过程系统模拟的基本方法ppt课件89