化工过程分析与合成-3

化工过程分析与合成化工过程分析与合成（三）（三）重庆大学化学化工学院重庆大学化学化工学院化学工程与工艺专业化学工程与工艺专业TSHY2/35第三章第三章 过程系统分解过程系统分解本章内容本章内容3.1 过程系统的分隔过程系统的分隔 3.2 不可分隔回路的分解不可分隔回路的分解 3.3 过程系统分解举例过程系统分解举例3.4 流程描述与分解应用（讨论）流程描述与分解应用（讨论）TSHY3/353.1 过程系统的分隔过程系统的分隔1. 1. 过程系统分隔方法过程系统分隔方法2. 2. 单元搜索法单元搜索法3. 3. 邻接矩阵法邻接矩阵法4. 4. 关联矩阵法关联矩阵法TSHY4/353.1.1 过程系统分隔方法过程系统分隔方法（1）循环回路（环路）循环回路（环路） （2）分隔的基本要求）分隔的基本要求 1）用信息流图描述过程系统； 2）分隔出每个部分可独立求解。（3）分隔的原则方法）分隔的原则方法 从一个节点出发，沿物流方向搜索，当流线回到出发的节点，就找到一个环路。一个环路原则上可以独立求解。ABCDEF10123456789TSHY5/353.1.2 单元搜索法单元搜索法 从一个单元（节点）开始，沿流线方向搜索，当流线回到起点单元时，所经过的单元就构成一个环路。如：ABCDEF（无环路） ACDEF（无环路） BCB（B、C节点经3、8流线构成环路） DEFD（D、E、F节点经5、6、9流线构成环路）ABCDEF10123456789TSHY6/353.1.3 邻接矩阵法邻接矩阵法1.1.可及原理可及原理 从某个单元出发，沿流线经一次前移可到达下一个单元，经若干次前移可到达流线上的某单元，这就是信息流图上的可及移动法。 在数学上，将邻接矩阵做布尔乘方运算，每乘方一次，代表各单元沿流线前移一次。经若干次乘方后，从某个单元也可到达另外的单元。当对角线元素为1时，表示前移到自己上，出现了环路。TSHY7/353.1.3 邻接矩阵法邻接矩阵法2.2.可及算法可及算法（1）构成初始邻接矩阵B 按邻接为1，不邻接为0，构成初始邻接矩阵BABCDEF101234567890 00 01 10 00 00 01 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 00 00 00 00 01 11 10 0B BA B C D E FABCDEFTSHY8/353.1.3 邻接矩阵法邻接矩阵法2.2.可及算法可及算法（2）一次可及B2 做布尔运算： B2=BBABCDEF101234567890 01 10 00 00 00 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 01 10 00 00 01 11 11 10 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 00 00 00 00 01 11 10 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 00 00 00 00 01 11 10 0B B2 2TSHY9/353.1.3 邻接矩阵法邻接矩阵法2.2.可及算法可及算法（3）二次可及B3 做布尔运算： B3=B2B 原则上，最长线路上有多少个单元，就需要做多少次可及。ABCDEF101234567891 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 00 00 01 10 01 10 01 10 00 01 10 01 10 00 00 01 11 11 11 10 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 00 00 00 00 01 11 10 00 01 10 00 00 00 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 01 10 00 00 01 11 11 10 0B B3 3TSHY10/353.1.3 邻接矩阵法邻接矩阵法3.3.求循环回路求循环回路（1）对各阶可及矩阵做并运算：R=BB2B3 并运算：1+1=1，1+0=1，0+1=1，0+0=01 11 11 10 00 00 01 11 11 10 00 00 01 11 11 10 00 00 01 11 11 11 11 10 00 01 11 11 11 10 00 01 11 11 11 10 01 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 00 00 01 10 01 10 01 10 00 01 10 01 10 00 00 01 11 11 11 10 00 01 10 00 00 00 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 01 10 00 00 01 11 11 10 00 00 01 10 00 00 01 10 00 00 00 00 00 01 10 00 00 00 00 00 01 10 01 10 00 00 00 01 10 00 00 00 00 01 11 10 0R RABCDEF10123456789TSHY11/353.1.3 邻接矩阵法邻接矩阵法3.3.求循环回路求循环回路（2）将并矩阵R与其转置矩阵RT做交运算：W=RRT 交运算：1+1=1，1+0=0，0+1=0，0+0=0 交运算矩阵中为1的方阵就构成循环回路。1 11 11 10 00 00 01 11 11 10 00 00 01 11 11 10 00 00 00 00 00 01 11 10 00 00 00 01 11 10 00 00 00 00 00 00 01 11 11 11 10 00 01 11 11 11 11 11 11 11 11 11 11 11 10 00 00 01 11 11 10 00 00 01 11 11 10 00 00 00 00 00 01 11 11 10 00 00 01 11 11 10 00 00 01 11 11 10 00 00 01 11 11 11 11 10 00 01 11 11 11 10 00 01 11 11 11 10 0WWABCDEF10123456789TSHY12/353.1.4 关联矩阵法关联矩阵法1. 建立关联矩阵建立关联矩阵 省略出入系统流线。有两个输出的划分为两个单元。ABCDEF102345689S2S3S4S5S6S8S9S10A-1A-1B1-11C1-1C-11D1-11E1-1F1-1TSHY13/353.1.4 关联矩阵法关联矩阵法2. 标记和调整关联矩阵标记和调整关联矩阵 在每一行标记一个元素，使每列只有一个元素被标记。 调整列，使标记元素都处于主对角线上。删除无标记列。S2S3S4S5S6S8S9S10A-1A-1B1-11C1-1C-11D1-11E1-1F1-1S2S10S3S8S4S5S6S9A-1-1B1-11C1-1C-1D1-11E1-1F1-1TSHY14/353.1.4 关联矩阵法关联矩阵法3. 通路搜索找环路通路搜索找环路（1）消去仅有对角线元素的行，以及对应的列，标记消去行的顺序。（2）若同一行上有两个或多个符号不同的元素，则在相应列中寻找是否有相反元素。如果有，表示存在环路，找出环路的所有行和列消去，并作为并集标记。（3）重复查找。 本例得到三个独立求解块： 单元 A 单元B、C 单元D、E、FS2S10S3S8S4S5S6S9标记标记A-1A-1B1-11(BC)C1-1C1-1D1-11(DEF)E1-1F1-1TSHY15/353.2 不可分隔回路的分解不可分隔回路的分解1. 1. 环路切割原则环路切割原则2. 2. 环路矩阵法环路矩阵法3. 3. 无多余切割法无多余切割法4. 4. 切割流线变量最少法切割流线变量最少法TSHY16/353.2.1 环路切割原则环路切割原则 寻找出来的环路通常由两个或两个以上的单元构成，起整体本身无法直接计算，需要将环路打开后进行迭代计算。对单个环路打开哪条流线都无关紧要，但对复杂系统，特别是有嵌套的环路系统，选择打开环路的流线就十分重要。因此，在选择切割流线时，需遵循以下基本原则： (1) 切割流线总数最少切割流线总数最少 可减少迭代量可减少迭代量 (2) 切割流线所含变量数最少切割流线所含变量数最少 可减少计算量可减少计算量 (3) 切割流线的总灵敏度最小切割流线的总灵敏度最小 迭代容易收敛迭代容易收敛 (4) 环路切割的总次数最少环路切割的总次数最少 迭代容易收敛迭代容易收敛TSHY17/353.2.2 环路矩阵法环路矩阵法以如下信息流图为例：以如下信息流图为例：其中共包含了四个环路：其中共包含了四个环路：ABCDE81234567BCDE2346AC85ABC125DE47环路环路d环路环路a环路环路c环路环路bTSHY18/353.2.2 环路矩阵法环路矩阵法1. 构成环路矩阵构成环路矩阵 行代表环路行代表环路 列代表流线列代表流线 1表示环路有该流线表示环路有该流线 0表示环路无该流线表示环路无该流线 为流频，流线经过的环为流频，流线经过的环路数路数 R为秩，环路中的流线数为秩，环路中的流线数BCDE2346AC85ABC125DE47环路环路d环路环路a环路环路c环路环路bS1S2S3S4S5S6S7S8Ra000010012b011101004c110010003d00010010212122111TSHY19/353.2.2 环路矩阵法环路矩阵法2. 消去不独立流线消去不独立流线 同一行中，若有流线的流频为：同一行中，若有流线的流频为： k j，则则k流线相对于流线相对于j流线不流线不独立，表示为：独立，表示为：Sk Sj。 在本环路矩阵中，在本环路矩阵中， S1S3 S2， S6S7 S4， S8 S5 消去不独立流线，剩下的流线则为侯选切割流线。消去不独立流线，剩下的流线则为侯选切割流线。S1S2S3S4S5S6S7S8Ra000010012b011101004c110010003d00010010212122111S2S4S5Ra0012b1104c1013d0102222TSHY20/353.2.2 环路矩阵法环路矩阵法3. 确定切割流线确定切割流线 在消去不独立流线后的环路矩阵中，秩为在消去不独立流线后的环路矩阵中，秩为1的环路上的流线即为的环路上的流线即为切割流线切割流线。 本例中的切割流线确定为本例中的切割流线确定为S4和和S5。 切割这两条流线后，所有环路都将打开。切割点就为迭代点。切割这两条流线后，所有环路都将打开。切割点就为迭代点。S2S4S5Ra0011b1102c1012d0101222ABCDE8123456754流程的计算顺序为：流程的计算顺序为： AEBCD（迭代流线迭代流线4-4，5-5）TSHY21/353.2.3 无多余切割法无多余切割法 为避免迭代计算的繁琐或难以收敛，对一个环路应为避免迭代计算的繁琐或难以收敛，对一个环路应避免进行多次切割。无多余切割法就是根据这一原则选避免进行多次切割。无多余切割法就是根据这一原则选择切割流线的。择切割流线的。 以如下信息流图为例：以如下信息流图为例：ABCEF81234567D9TSHY22/353.2.3 无多余切割法无多余切割法1. 建立环路矩阵建立环路矩阵 在矩阵最下方在矩阵最下方增加一行判别标增加一行判别标记记。ABCEF81234567D9S1S2S3S4S5S6S7S8S9a1111b11111c1111d11121331211TSHY23/353.2.3 无多余切割法无多余切割法2. 寻找切割流线寻找切割流线（1）选择）选择最大的流线为切最大的流线为切割流线，标记其的割流线，标记其的为为1。（2）与该流线处于同一环路）与该流线处于同一环路上的流线上的流线标记为标记为0。 切割流线切割流线4，打开，打开a、b、c三个环路。三个环路。ABCEF81234567D9S1S2S3S4S5S6S7S8S9a1111b11111c1111d11121331211000100010（3）重复以上步骤，切）重复以上步骤，切割流线割流线8，打开环路，打开环路d。TSHY24/353.2.3 无多余切割法无多余切割法3. 确定计算顺序确定计算顺序 根据已知输入计算输出的原则，起计算顺序为：根据已知输入计算输出的原则，起计算顺序为：FEDABC（迭代迭代4-4，8-8流线）流线）ABCEF81234567D948 本例中如果选择本例中如果选择切割流线切割流线5和流线和流线9也是一样的。也是一样的。TSHY25/353.2.4 切割流线变量最少法切割流线变量最少法 如果已知上例中各条流线的变量数，则可以在无多余如果已知上例中各条流线的变量数，则可以在无多余切割的基础上在考虑选择流线变量上的流线进行切割，切割的基础上在考虑选择流线变量上的流线进行切割，以减少迭代计算量。以减少迭代计算量。 例如已知各流线的变量数为：例如已知各流线的变量数为：ABCEF81234567D9S1S2S3S4S5S6S7S8S9v842635112TSHY26/351. 建立环路矩阵建立环路矩阵 在矩阵中在矩阵中的的下方增加一行变下方增加一行变量数量数v。ABCEF81234567D9S1S2S3S4S5S6S7S8S9a1111b11111c1111d11121331211v8426351123.2.4 切割流线变量最少法切割流线变量最少法TSHY27/352. 选择切割流线选择切割流线（1）首先用无多余切）首先用无多余切割法选择切割流线。割法选择切割流线。（2）在）在相同的情况相同的情况下，再选择下，再选择v小的流线。小的流线。S1S2S3S4S5S6S7S8S9a1111b11111c1111d11121331211v8426351120000100103.2.4 切割流线变量最少法切割流线变量最少法 因此，本例中选择切割流线因此，本例中选择切割流线5和流线和流线8。计算顺序为：。计算顺序为：DABCFEABCEF81234567D958TSHY28/353.3 过程系统分解举例过程系统分解举例环氧乙烷提纯系统的分解环氧乙烷提纯系统的分解1. 环氧乙烷提纯流程图环氧乙烷提纯流程图ABCDEFGHIJKLMNOTSHY29/352. 用信息流图描述流程用信息流图描述流程ABCDEJIHGFKLMNO22123456789101112131415161718192021混合器混合器分流器分流器塔设备塔设备换热器换热器省略了输入和输出系统的流线省略了输入和输出系统的流线TSHY30/353. 系统分隔系统分隔ABCDEJIHGFKLMNO22123456789101112131415161718192021寻找出所有循环回路：寻找出所有循环回路： 环路环路a：BCHB 环路环路b：BCDEB 环路环路c：BCDEFJ KNOAB 环路环路d：BCDEJIGB 环路环路e：DEFJ KMD 环路环路f： DEFJKNOLMD 环路环路g：JKJ 环路环路h：KNOLK TSHY31/354. 环路切割环路切割建立环路矩阵，选择切割流线：建立环路矩阵，选择切割流线：2、5、10、20 S12345678910111213141516171819202122a111b1111c1111111111d11111111e111111f111111111g11h111141135144312111233121110100100001000000000100TSHY32/355. 计算顺序计算顺序ABCDEJIHGFKLMNO22123456789101112131415161718192021EFJIGKNOABCLMD TSHY33/353.4 流程描述与分解实例流程描述与分解实例讨论任务讨论任务1. 1. 合成氨变换流程合成氨变换流程(8)(8)2. 2. SOSO2 2四段两次转化流程四段两次转化流程(9)(9)讨论要求讨论要求1. 1. 流程图的信息流图描述流程图的信息流图描述 2. 2. 系统分隔（建议邻接矩阵法）系统分隔（建议邻接矩阵法）3. 3. 切割流线及运算顺序的确定切割流线及运算顺序的确定TSHY34/353.4.1 中型合成氨变换流程中型合成氨变换流程 TSHY35/353.4.2 SO2四段两次转化流程四段两次转化流程