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

第4章 化工过程模拟的基本方法第4章 化工过程模拟的基本方法1主要内容过程单元过程单元过程单元的自由度过程单元的自由度过程单元的模型化与模拟过程单元的模型化与模拟过程系统过程系统过程系统的自由度过程系统的自由度系统结构的识别系统结构的识别系统的分解系统的分解序贯模块模拟法序贯模块模拟法原理原理寻找最佳断裂流股寻找最佳断裂流股断裂流股的收敛断裂流股的收敛面向方程模拟法及联立模块模拟法原理面向方程模拟法及联立模块模拟法原理过程系统的优化过程系统的优化常用商业化流程模拟软件简介常用商业化流程模拟软件简介主要内容过程单元2公司结构常减压催化焦化重整芳烃原油原油产品产品产品产品炼油厂换热器车间塔器车间钢材钢材产品产品产品产品机械厂石化公司不相关子系统不可分隔子系统回路系统公司结构常减压催化焦化重整芳烃原油产品产品炼油厂换热器塔器钢3过程系统模拟一般步骤：确定需模拟的系统识别不相关子系统（一级子系统，各子系统独立处理）一级子系统自由度分析找到分解成规模最小的二级子系统的决策变量各一级子系统内不可分隔子系统的识别（二级子系统，各二级子系统可依次求解）各二级子系统依次处理找回路、确定最佳断裂位置、确定计算顺序相关单元操作建模，确定迭代方法过程系统模拟一般步骤：确定需模拟的系统4第1节 流程的自由度分析1.一些基本概念过程(Process)：对原料进行某些物理或化学变换，使其性质发生预期的变化机械加工不能称为过程系统：由相互联系，相互作用的若干组成部分结合成的具有特定功能的总机体过程系统：由各种过程构成的系统第1节 流程的自由度分析1.一些基本概念5流程:描述化工生产的物料流向及能量流向及装置特点的过程REACTORFEEDRECYCLEREAC-OUTCOOLCOOL-OUTSEPPRODUCT流程:REACTORFEEDRECYCLEREAC-OUTC6模型：复杂的A简单的B来替代。研究B来预测A的行为 实物模型 /数学模型模拟：对某一描述实际过程的数学模型利用数学方法进行求解，并对结果作出解释模型：复杂的A简单的B来替代。研究B来预测A的行为72.物流自由度、单元自由度及系统自由度物流(stream)Dhem定理：对于一个已知每个组分初始质量的封闭体系，其平衡状态完全取决于两个独立变量，而不论该体系有多少个相，多少个组分或多少个化学反应 未知组成：C+2已知组成：2 2.物流自由度、单元自由度及系统自由度物流(stream)8过程单元自由度分析方法过程单元自由度:可改变单元操作状态的独立变量数目列方程法：列出现象方程、限制方程等 自由度=变量数-方程数l描述规则法：要完全描述一单元设备的操作，必须确定的自由度的数目必定等于能由设备结构确定或能用外部手段控制的变量的数目l公式分析法：过程单元自由度分析方法过程单元自由度:9常见过程单元自由度单元名称自由度数常规指定变量混合器0分流器S-1流量分配比闪蒸器2闪蒸温度、压力泵、节流阀1出口压力压缩机/膨胀机2绝热多变效率，出口压力换热器1某一物流出口温度反应器2+r反应程度，绝热，压降常规精馏塔5塔板数，进料位置，R,B,操作压力平衡级3侧线抽出比，压力，换热量常见过程单元自由度单元名称自由度数常规指定变量混合器0分流器10过程系统自由度过程系统自由度:确定系统状态的独立变量数目 系统自由度=单元自由度+进料自由度过程系统自由度过程系统自由度:11第节 过程系统结构的计算机识别主要内容 1.过程系统结构有向图2.过程系统结构的矩阵表示3.系统结构的识别第节 过程系统结构的计算机识别主要内容121.过程系统结构有向图节点边1.过程系统结构有向图节点边132.系统结构的矩阵表示节点-节点节点相邻矩阵节点-边关联矩阵边-边弧相邻矩阵2.系统结构的矩阵表示节点-节点143.系统结构的识别可及矩阵法不相关子系统不相关子系统过程系统：建模过程中可分别独立处理方程系统：写出方程组事件矩阵，Himmelblau算法识别不可分隔子系统不可分隔子系统可及矩阵法，索引矩阵法、图解法，Steward通路法3.系统结构的识别可及矩阵法不相关子系统15方程系统识别-Himmelblau算法：在mm事件矩阵M中，选出非零元素最多的列k。保留M中k列内每个零元素对应的行，k列中为1的元素所对应的行用布尔加法合并成一行排列在最后。得到的新的jm的布尔矩阵记做M(0)；重复，从而得到序列M，M(0)，M(N)；最终得到矩阵M(N)，其每一列只有一个非零元素，其每一行与原方程系统中的不相关子系统对应。方程系统识别-Himmelblau算法：在mm事件16例例17过程系统识别不相关子系统分析不可分隔子系统可及矩阵过程系统识别不相关子系统18公司结构常减压催化焦化重整芳烃原油原油产品产品产品产品炼油厂换热器车间塔器车间钢材钢材产品产品产品产品机械厂石化公司不相关子系统不可分隔子系统回路系统公司结构常减压催化焦化重整芳烃原油产品产品炼油厂换热器塔器钢19系统网络图 求解顺序：H-ABCDE-FG-I系统网络图 求解顺序：H-ABCDE-FG-I20第第3 3节节 序贯模块模拟法序贯模块模拟法主要内容：基本原理循环流股的断裂与迭代断裂变量的收敛第3节 序贯模块模拟法主要内容：21一、基本原理一、基本原理混合器Submix(F1,F2,P1)P1=F1+F2EndSub分割器SubSPLT(F1,P1,P2,ALFA)P1=F1*ALFAP2=F1*(1-ALFA)EndSub收敛单元S4S4单元子程序内容：基础数据 单元方程组 求解算法单元子程序功能：根据给定的进料 条件和设备参数，计算出单元输出 结果(预测型)一、基本原理混合器 收敛单元S4S4单元子程序内容：22设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,进行流程模拟计算23序贯模块法模块的特点 单向性单元操作输入输出设计规定流程计算收敛单元-循环流/设计规定积木式序贯模块法模块的特点 单向性单元操作输入输出设计规定流程计算24断裂位置的影响：如果如果S2S2与与S4S4的自由度不同，需迭代的变量数也将不同的自由度不同，需迭代的变量数也将不同断裂位置的影响：如果S2与S4的自由度不同，需迭代的变量数也25二、再循环流股的断裂二、再循环流股的断裂方法1方法2原则：原则：将所有闭合回路全部打开将所有闭合回路全部打开二、再循环流股的断裂方法1方法2原则：将所有闭合回路全部打开262.1 最优断裂准则(1)被切断的流股数最少；(2)被切断的流股变量数最少；(3)被切断的流股的权重因子之和最少；(4)回路切断的总次数最少。通常选择原则：满足（满足（4）的基础上，选（）的基础上，选（1）或（）或（2）2.1 最优断裂准则(1)被切断的流股数最少；27抉择依据：最少计算时间计算时间：计算时间：断裂方式；流程及变量灵敏度有效计算时间：流程计算；断裂流股迭代 迭代时间：迭代次数；收敛速度抉择依据：最少计算时间计算时间：282.2 回路矩阵 回路A:单元S2IIIS4回路B:单元IS1S2IIIS5I回路C:单元IS1S2IIIS3S6I回路D:单元S2IIIS3S72.2 回路矩阵 回路A:单元S2IIIS429回路矩阵 矩阵元素aij定义为回路矩阵 矩阵元素aij定义为302.3 UpadhyeGrens断裂法(II)1)有关术语 有效断裂组有效断裂组能够把全部简单回路至少切断一次的断裂流股的集合。S2,S1,S3,S4，S1,S2,S5 S1,S3 多余断裂组多余断裂组若从一个有效断裂组中至少可以除去一个流股，而且得到的断裂组仍为有效断裂组；或者存在着对一个回路的二次断裂。S1,S2,S5 S1,S3,S4 2.3 UpadhyeGrens断裂法(II)1)有关31非多余断裂组非多余断裂组除去多余断裂组以后的有效断裂组。断裂族断裂族具有相同计算顺序的有效断裂组的集合。非多余断裂组除去多余断裂组以后的有效断裂组。322)、替代规则令D1为一有效断裂组，Ai为全部输入流股均属于D1的单元，将Ai的所有输入流用Ai的全部输出流替代，构成新的断裂组D2，则D2也是有效断裂组对直接迭代，D2与D1具有相同的收敛性质 由于全部回路被打开，至少会有一个单元的全部输入条件由于全部回路被打开，至少会有一个单元的全部输入条件已知，可以作为计算的起点，计算出其输出流股。而该流股已知，可以作为计算的起点，计算出其输出流股。而该流股又将作为后续单元的输入流股，产生连锁效应！又将作为后续单元的输入流股，产生连锁效应！2)、替代规则 令D1为一有效断裂组，Ai为全部输333)、Westerberg算法从任何一有效断裂组开始，运用替代规则如果在任何一步中出现重复断裂组，则消去其中的重复流股，消去重复流股后形成的新断裂组作为新的起点重复、，直到没有重复断裂组出现，且每个树枝上的断裂组重复出现为止，从最后一个新的起点开始，其后出现的所有不重复的断裂组构成非多余断裂族非多余断裂族中总数最小的断裂组为最优断裂组3)、Westerberg算法从任何一有效断裂组开始，运用34例 例 35非多余断裂组非多余断裂组断裂变量数S29S1,S4,S72+3+2=7S3,S4,S52+3+3=8S4,S5,S6,S73+3+4+2=12准则1准则2全部满足准则4：每个回路都断裂，且只断裂每个回路都断裂，且只断裂1次次非多余断裂组非多余断裂组断裂变量数S29S1,S4,36三、断裂流股变量的收敛 三、断裂流股变量的收敛 37收敛单元功能 提供循环流x的初值x。根据初值x，及其它输入条件进行流程计算，求出循环流出口值y，比较x和y，若x与y之差满足精度要求则停止计算，否则转。根据x与y值以一定的方式产生x的新估计值x1，转 收敛单元功能 提供循环流x的初值x。38收敛单元的要求 对初值要求不高 易得，不易引起迭计算的发散提供的初值组数少例:初值的影响收敛单元的要求 对初值要求不高 易得，不易引起迭计算的发39数值稳定性好数值稳定性单调收敛数值稳定性好数值稳定性单调收敛40续数值稳定性振荡收敛续数值稳定性振荡收敛41续数值稳定性单调发散续数值稳定性单调发散42续数值稳定性振荡发散续数值稳定性振荡发散43收敛速度快 影响收敛的速度主要因素：迭代次数 流程计算次数(Pass)迭代中矩阵求逆次数收敛速度快 影响收敛的速度主要因素：44占内存小 占内存小 45四、常用迭代法 直接迭代法有界Wegstein法Broyden法主特征值法联立超松弛法 Newton-Raphon法Marquardt法四、常用迭代法 直接迭代法46阻尼直接迭代法 q取值的一般原则：稳定性较好流程，取小于零，加快收敛速度如合成氨过程取q=-0.75。收敛速度快稳定性差，可取01，以改善稳定性如反应器q=0.5。阻尼直接迭代法 q取值的一般原则：47几种收敛策略的比较 方法偏导数函数计算次数稳定性速度直接迭代法否1较好慢Newton-Raphon法是n较差快Broyden法否1较好适中Marquardt法是n好适中Wegstein法否1较好适中主特征值法否1较好适中几种收敛策略的比较 方法偏导数函数稳定性速度直接迭代法否1较48收敛判据 绝对误差 相对误差 收敛判据 绝对误差 49五、解决设计问题 估计反应单元的温度为估计反应单元的温度为T T估计再循环物流估计再循环物流S4S4依次计算混合单元、反应单元、分离单元，得到新的依次计算混合单元、反应单元、分离单元，得到新的S4S4的的比较比较S4S4与与S4S4，若两者相等则进行下一步，若不相等则返回，若两者相等则进行下一步，若不相等则返回 在在收收敛敛单单元元内内比比较较S5S5和和设设计计值值，若若两两者者不不相相等等则则返返回回，若若相等则计算结束相等则计算结束五、解决设计问题 估计反应单元的温度为T50同时收敛策略 控制模块的设置增加了迭代循环圈，导致计算量的增加D设计规定向量 H过程系统方程组p决策变量与系统参数向量同时收敛策略 控制模块的设置增加了迭代循环圈，导致计算量的增51第4节 常用流程模拟软件简介主要内容：流程模拟技术的进展流程模拟软件的结构第4节 常用流程模拟软件简介主要内容：52一、流股模拟技术的进展流股模拟技术序贯模块法联立方程法（面向方程法）联立模块法一、流股模拟技术的进展流股模拟技术531 联立方程法原理Sub flow(S1,ALFA,S2,S3,S4)S2=S1+S4 混合器模型 S3=S2*ALFA S4=S2*(1 ALFA)分割器模型End Sub 特点：特点：单元模型只提供方程，系统方程单元模型只提供方程，系统方程组一起求解组一起求解无循环物流及设计规定的特殊处无循环物流及设计规定的特殊处理问题理问题1 联立方程法原理Sub flow(S1,ALFA,S541.1模拟结构 循环物流方程设计规定方程仅提供方程仅提供关联式多!特殊处理1.1模拟结构 循环物流方程仅提供方程仅提供关联式多!特551.2 物性系统的处理占用7080%的计算量处理方法：1）不提供方程组，仅提供子程序2）提供 K,H 的方程组，其它物性由子程序调用1.2 物性系统的处理占用7080%的计算量56 1.3 序贯法与联立方程法的对比 1.3 序贯法与联立方程法的对比 571.3.1 序贯模块法优点：与实际过程的直观联系强，软件的建立、维护和扩充很方便；易于通用化；需要的计算机内存较小；易于诊断出错位置；缺点：循环物流及设计问题：慢1.3.1 序贯模块法优点：缺点：循环物流及设计问题581.3.2 面向方程法优点 便于实际用户的提出的质量性能要求。有设计规定的系统，仅增加几个方程 用空间换取时间。相同的过程单元出现多次，将占用更多的 存储空间，自由度的选择也可能各不相同 以空间换取求解的数值稳定性。有较多回流高度交互作用流程相对容易收敛1.3.2 面向方程法优点591.3.2 面向方程法缺点 要求提供较好的初值 难于自动产生流程的模型方程组，难以通用化 不容易诊断出现错误的地方。1.3.2 面向方程法缺点602.联立模块法联立模块法又称双层法，将过程系统的近似模型方程与单元模块交替求解联立模块法兼有序贯模块法和面向方程法的优点。既能使用序贯模块法积累的大量模块，又能将最费计算时间的流程收敛和设计约束收敛等迭代循环合并处理，通过联立求解达到同时收敛2.联立模块法联立模块法又称双层法，将过程系统的近似模型方61第4章-过程系统模拟的基本方法ppt课件62过程系统稳态模拟三种方法的比较过程系统稳态模拟三种方法的比较过程系统稳态模拟三种方法的比较63二、过程系统优化模型的建立过程模拟模型给定参数可调参数经济模拟模型经济参数最优化程序过程约束条件经济约束条件经济特性指标二、过程系统优化模型的建立过程模拟模型给定参数可调参数经济64三、过程系统模拟软件的结构流程模拟软件类型通用软件专用软件三、过程系统模拟软件的结构流程模拟软件类型651.流程模拟系统结构1.流程模拟系统结构662.模拟步骤2.模拟步骤67 输入流程拓朴结构信息 输入单元计算所需的必要信息，选择物性及热力学性质，传递性质计算方法 进行全流程的模拟计算 解释、考核计算结果单元名称进料混合器1反应器精馏塔混合器2分流器产品输入物流号1，7234，856，9输出物流号1，8234，596，7 输入流程拓朴结构信息单元名称进料混合器1反应器精馏塔混合68过程系统模拟一般步骤：确定需模拟的系统识别不相关子系统（一级子系统，各子系统独立处理）一级子系统自由度分析找到分解成规模最小的二级子系统的决策变量各一级子系统内不可分隔子系统的识别（二级子系统，各二级子系统可依次求解）各二级子系统依次处理找回路、确定最佳断裂位置、确定计算顺序相关单元操作建模，确定迭代方法过程系统模拟一般步骤：确定需模拟的系统69软件名称件名称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,0002000多种3000多组有特种工特种工艺系系统模型模型GLC气液接触模型、固定床反应器模型、溶液结晶模型、聚合过程模型、燃料电池模型多种PEP工艺模型（SRI的工艺包）与SPYRO软件接口炼油专用模块RefSYS,包括催化裂化，催化重整反应，异构化反应器，烷基化反应器，加氢裂化反应器，加氢精制/处理反应器炼油厂各种典型装置模型优化化计算算有有有有经济核算核算有有/有部分动态模模拟计算算均为动态(稳态时间方程为零)需AspenDynamics连接可直接转入动态模拟与Dynsim连接进中国市中国市场时间2002年80年代初90年代80年代末常见流程模拟软件对比软件名称gPROMSAspen PlusHysysPRO/I70第4章-过程系统模拟的基本方法ppt课件71第4章-过程系统模拟的基本方法ppt课件72