天津大学化工设计.ppt

《天津大学化工设计.ppt》由会员分享，可在线阅读，更多相关《天津大学化工设计.ppt（46页珍藏版）》请在装配图网上搜索。

1、3.5 化工流程MB和HB的联立计算（CB）,CB计算与仅做（单独做）MB计算时，有下面几点区别：,1、流程中所有流股的流股变量加上温度（T）；,2、流程中所有单元的单元变量加上描述传热的单元变量（dQ/dt或Q）；,3、流程中一个单元牵涉N种组份，就可以列出N个独立的MB方程（分流器除外），但一个单元无论其牵涉到多少种组份，都只能列出一个HB方程。,4、因为流股的焓值与流股的相态有关，做CB时还要注意流股的相态；,5、一个流股焓值的计算，大多数情况下可以考虑各组份的分摩尔焓按其在流股中的摩尔分率进行加和计算（即按理想状态，各分组份焓简单加和），少数情况下，要考虑采用偏摩尔焓加和；,6、dQ/

2、dt或Q，有两种表示方式： 、 dQ/dt或Q焓（进）焓（出） dQ/dt或Q0，则意味系统（或单元）向外散热，反之逆然。,、 dQ/dt或Q焓（出）焓（进） dQ/dt或Q0，则意味系统（或单元）需要向外界获取热量，反之逆然。,dQ/dt或Q0，表明系统（或单元）绝热。,7、焓是状态函数，流股焓的计算必须有基准态（参考态）：,焓值计算时参考态,标准态（25）,任意温度,流程中某个流股的温度,注意：在一个计算过程中，参考态一经选定，就不能再变更。,优点：数据来源方便,可能简化计算,8、MB和HB联立计算后，计算复杂程度加大，主要是因为焓值与温度的非线性关系，导致HB方程通常是非线性方程。,9、

3、在考虑CB问题时（不是仅做MB计算时），看一个流程设计条件是否给得正确，不能只看其过程的MB计算自由度，而必须看其过程的CB自由度。,例题： 有一个常压绝热猝冷过程（过程中无化学反应），所有组份含量均为摩尔分率，其流程示意图如下所示：,3,T2400 1000 mol/h C6H6 0.40 C6H5CH3 0.30 CH4 0.10 H2 0.20 （All in gas）,1,2,T120 C6H6 （L） 1.00,T3200 C6H6 C6H5CH3 CH4 H2 （All in gas）,绝热猝冷过程示意图,试作该流程的自由度分析，完成其物料衡算和热量衡算计算。,解：,流程自由度分析

4、表,1,0,0,4,4,0,9,MB,猝冷器,CB,0,自由度,0,已知其它关系式数,1,已知单元变量数,7(4+3),已知流股变量数,5(4+1),衡算方程数,1,单元变量数,12(9+3),流股变量数,流股变量加T,多了描述传热的 单元变量dQ/dt,多了一个HB方程,dQ/dt=0,流程自由度分析表,这个流程仅仅看MB计算自由度分析，不能就说它的设计条件给得不正确，而应该看器CB计算的自由度分析结果。,这种情况说明，该流程的MB不能单独求解，而必须与其HB方程联立起来，进行CB计算。,具体计算：,F31000F1 X3,C6H5CH3 F3300 X3,CH4 F3100 X3,H2 F

5、3200,4个MB方程,其中：5个未知变量，只有4个方程。,在列其HB方程之前，先考虑其参考态选什么温度？,如果热容数据好查！最简化计算的参考态温度应该选为3号流股的温度。,因为其未知的流量组成变量最多,流程（单元）的HB方程（选200为参考态温度）：,F2H2（400）H2（200）F1H1（200）H1（20）,列出详细HB方程如下：,苯在1atm，80.11沸点时的摩尔汽化热,注意Cp的量纲问题，及其适用温度范围。,例题：一个合成氨生产工艺流程图如下所示：已知 . 反应器1中N2的转化率为10； 图中组成均为摩尔分率； 5号流股为液相外，其余均为气相流股； . 除分离器外，其它单元均为绝

6、热操作； . 两反应器中均只有一个反应： N2 3H2 2NH3 要求： . 做过程的自由度分析； . 指出衡算计算的计算顺序。,解：,首先，根据流程特点，判断各流股的组份数。,NH3 N2 H2 -50,425,进料,1,50 N2 0.25 H2 0.75,反应器2,反应器1,混合器,换热器,分 离 器,2,8,3,7,425,6,N2 0.25 H2 0.75 50,4,dQ/dt,5,535,(3),(3),(3),(3),(3),(2),(2),(2),流程整体（OB）的特点：一个带一种反应的非绝热的虚拟反应器。,反应器1和反应器2自由度分析表,5,1,3,1,1,0,1,52,2,

7、3,1,11,2,0,1,6,1,3,0,0,0,4,8,2,3,0,1,2,1,3,混合器和换热器自由度分析表,4,4,4,自由度,0,0,0,已知附加式数,1,1,0,已知单元变量数,4,3,1,已知流股变量数,1,1,HB方程,3,3,MB方程数,1,1,0,单元变量,54,11,8,流股变量,CB,CB,MB,换热器,混合器,+2,2,5,4,自由度,0,0,0,0,已知附加式数,0,0,0,0,已知单元变量数,3,1,3,1,已知流股变量数,1,1,HB方程,3,3,3,3,MB方程数,2,1,1,0,单元变量,7,5,11,8,流股变量,CB,MB,CB,MB,整体OB,分离器,分

8、离器和整体（OB）自由度分析表,流程自由度分析结果一览表,0,0,0,0,0,0,0,0,0,0,0,0,已知附加式数,0,0,5,0,0,1,1,0,1,0,2,1,已知单元变量数,+2,+2,1,+5,+4,+4,+4,+4,+3,+4,+1,+1,自由度,3,1,8,3,1,4,4,1,2,0,2,1,已知流股变量数,1,5,1,1,1,1,1,HB方程数,3,3,12,3,3,3,3,3,3,3,3,MB方程数,2,1,7,1,0,1,1,0,2,1,2,1,单元变量数,7,5,24,11,8,9,11,8,8,6,7,5,流股变量数,CB,MB,CB,CB,MB,CB,CB,MB,C

9、B,MB,CB,MB,整体OB,过程,分离器,换热器,混合器,反应器2,反应器1,显然是一个弹性设计，必须挑选一个流股的流量作为计算基准。,最好选F71000mol/h为计算基准。,从自由度分析结果看，应该先解反应器1的MB，再解其CB（HB），而且完全求解。,解完反应器1后，应该接着解哪个单元？,完全求解反应器1的MB和HB意味着什么？,混合器和换热器的更新自由度分析表,4,+1,7,0,5,+3,分析结果：接着解混合器的CB，完全求解。,解完混合器后，应该接着解哪个单元？,反应器2和整体OB的更新自由度分析表,3,+1,5,0,2,4,+1,+1,分析结果：接着解反应器2的CB，完全求解。

10、,解完反应器2后，应该接着解哪个单元？,换热器的更新自由度分析表,8,0,分析结果：接着解换热器的CB（即HB），完全求解。,解完换热器后，应该接着解哪个单元？,分离器的更新自由度分析表,5,0,分析结果：接着解分离器的MB，最后解其HB，完全求解。,只剩分离器和整体OB，二者解其一即可。,8,0,综上分析结果，得到一个完全求解的衡算序列：,反应器1MB,反应器1HB,混合器CB,反应器2CB,换热器HB,分离器MB,分离器HB,具体求解：略,例题：一个甲烷生产流程初步设计的流程草图如下，试作流程的自由度分析，并指出物料及热量衡算计算顺序。 设计已知及要求为： 1、除了换热器1和分离器外，其它

11、单元均为绝热操作； 2、8号流股为纯水液相流股，其余流股均为气相流股； 3、整个流程处于100psia压力状态，压力对焓计算的影响可忽略； 4、产品流中甲烷的分摩尔流量为1000kmol/h, 图中组成均为摩尔分率； 5、反应器中只有：CO + 3H2 CH4 + H2O； 6、1号流股中只含3个组份，且CO和H2的摩尔比为1：2.9。,注：psia（英制压力单位）,混合器,反应器,换热器1,换热器2,分离器,分配器,dQ1/dt,dQ2/dt,1,2,进料 200 CO H2 H2O,3 1000,4 500,5 200,6 CO 0.12 H2 CH4 0.50 H2O 200,7 200

12、,H2O 100,8,H2O 0.01 100,9,10,(1),(4),(4),(4),(4),(4),(4),(4),(4),(3),解：,（一）混合器和反应器的自由度分析,+3,+3,+6,+5,自由度,0,0,1,1,已知附加式数,1,0,1,0,已知单元变量数,3,2,2,1,已知流股变量数,1,1,HB方程,4,4,4,4,MB方程数,2,1,1,0,单元变量数,10,8,14,11,流股变量数,CB,MB,CB,MB,反应器,混合器,（二）换热器1和换热器2的自由度分析,+2,0,0,4,1,1,6,CB,换热器1,+5,自由度,0,已知附加式数,1,已知单元变量数,6,已知流股

13、变量数,1,HB方程,MB方程数,1,单元变量数,12,流股变量数,CB,换热器2,（三）分流器和分离器的自由度分析,+2,+2,+1,+2,自由度,0,0,0,0,已知附加式数,0,0,1,已知单元变量数,6,3,6,3,已知流股变量数,1,1,HB方程,4,4,1,1,MB方程数,1,0,1,单元变量数,12,9,9,6,流股变量数,CB,MB,CB,MB,分离器,分流器,回自由度总表,（四）过程和整体（OB）的自由度分析,+1,+1,自由度,1,1,已知附加式数,0,0,已知单元变量数,6,3,4+8,已知流股变量数,1,HB方程,4,4,MB方程数,2,1,单元变量数,11,8,18+

14、10,流股变量数,CB,MB,CB,整体OB,过程,+1,+1,-1,+2,+2,+1,+2,+5,+2,+3,+3,+6,+5,自由度,1,1,1,0,0,0,0,0,0,0,0,1,1,已知附加式数,0,0,4,0,0,1,0,1,0,1,0,1,0,已知单元变量数,6,3,12,6,3,6,3,6,4,3,2,2,1,已知流股变量数,1,6,1,1,1,1,1,1,HB方程,4,4,13,4,4,1,1,4,4,4,4,MB方程数,2,1,7,1,0,1,0,1,1,2,1,1,0,单元变量数,11,8,28,12,9,9,6,12,6,10,8,14,11,流股变量数,CB,MB,CB

15、,CB,MB,CB,MB,CB,CB,CB,MB,CB,MB,整体OB,过程,分离器,分流器,换2,换1,反应器,混合器,流程自由度分析结果一览表,设计不对？问题在哪？,回自由度（三）,0,+2,0,5,问题：这个流程的衡算计算采用“变量置换法”有简化 计算的效果吗？,本章总结,化工流程衡算计算的大致步骤及每步相应的注意事项：,第一步：根据流程特点和已知设计条件确定流程中各流股含有的组份种类数。（如果是CB计算，还要注意各流股的相态。）,第二步：给流程中所有流股编号。（注意将分流器的有关流股编号相连。）,第三步：确定是否仅仅做MB计算（流股变量不考虑温度，单元变量不考虑dQ/dt或Q）；做CB

16、计算（流股变量考虑温度，单元变量考虑dQ/dt或Q）。注意多相流股的流股变量问题。,第四步：确定流程中所有反应器各自的独立反应个数。,第五步：如果做CB计算，应确认各单元的dQ/dt(或Q)是否已知。（绝热还是非绝热操作。）,第六步：根据流程特点判断整体（OB）的特征。注意： （1）、只要流程中有反应器（无论多少个反应器存在），则流程整体相当于一个大的虚拟反应器。进一步结合各反应器分别存在的反应个数，来判断OB的独立反应种类数。 （2）、只要流程中有单元不绝热，则OB是一个虚拟的非绝热单元。所有单元都绝热，OB才绝热。,第七步：做全流程的自由度分析。注意 （1）、间壁式换热器只有CB问题（其单

17、独MB无意义）。 （2）、分流器的自由度分析有特殊性（MB方程只有1个，还要注意其MB方程和HB方程的独立性）。 （3）、注意过程一栏中MB方程各单元（不包括OB）加和时的独立性。 （4）、注意过程一栏中已知附加式各单元（包括OB）加和时的独立性。,第八步：从全流程自由度分析结果及流程特点出发，判断是否是弹性设计，是否可用变量置换法简化衡算计算。（注意弹性设计时计算基准的挑选；注意变量置换后，将衡算计算结果换回原设计条件的流量要求。）,第九步：从全流程最后的自由度分析结果出发，逐步利用更新自由度分析表，判断衡算计算的步骤，确定衡算计算步骤。（各步骤单元部分求解序列居多；部分求解和完全求解混杂。

18、）,第十步：具体求解，并给出衡算结果一览表。（计算时注意数据的量纲。）,以上18步，即使是用计算机辅助衡算计算软件，也必须人工处理，而且还大有必要。计算机辅助衡算计算软件主要是帮助处理910步。,3.6 计算机辅助物料衡算与热量衡算,例题： 有一个常压绝热猝冷过程（过程中无化学反应），所有组份含量均为摩尔分率，其流程示意图如下所示：,我们来用PROCESS II 6.2版计算机辅助化工设计软件系统进行衡算计算。,第一步建立流程图 首先在PROCESS II软件空白流程界面中建立过程的流程示意图，输入所有单元设备及所有的流股，确定各流股的走向，并给所有流股编号。,第二步定义组份 给流程中的每个流

19、股分别说明其组份种类。同时，为了利用PROCESS II软件带有的物性数据库。第二步必须将流程中牵涉的所有组份种类输入PROCESS II软件。,第三步定义一个热力学方法 根据过程的压力、温度等情况在众多的热力学方程（如Pen-Robinson cubic equation of state、Ideal、NRTL、UNIQUAC等等）中，挑选一个适宜的热力学方法，软件本身能自行进行组份的物性数据计算。,第四步指定工艺装置和物流的已知数据 因为PROCESS II软件本身只能进行正向模拟计算，必须先将输入流股的已知数据输入软件中（这些已知输入流股参数在软件计算过程中始终保持不变），输入流股的一些

20、未知参数也必须输入初值（最好尽量合理一些，该初值可以由软件计算改变），输入各单元装置已知的设备参数。 在输入过程中可以改变参数的计量单位。,第五步输入设计约束条件 将输出流股的已知参数、其它设计约束条件以模拟计算控制模块（全部以已知方程方式）形式输入PROCESS II软件中。并设定模拟计算收敛精度和最大迭代计算次数。,第六步模拟计算 启动PROCESS II软件开始模拟计算，如果输入的数据满足了软件模拟计算的自由度要求，软件将开始模拟计算。随着PROCESS II软件模拟计算的进行，流程中各变量将不断变化，直到达到收敛精度或迭代次数，PROCESS II软件将给出计算结果。,第七步确定计算结

21、果的输出格式 软件本身有默认的计算结果输出格式，但可以按设计者爱好或设计计算要求，挑选确定计算结果的输出格式，确定数据的有效数据位数。,用PROCESS II软件计算以上绝热猝冷过程,用一台奔腾4 2.6（内存1M）的台式计算机，瞬间（CPU用时不到1秒）即可完成计算，计算的结果见下表。,有如下一个无反应的单精馏塔流程的设计方案:,按设计投产运行后，“查定”的该精馏塔实际运行结果：,精馏塔,1,2,3,100 kmol/h C3 0.20 i-C4 0.30 i-C5 0.20 C5 0.30,62.0087 kmol/h C3 0.3133 i-C4 0.40 i-C5 0.2742,37.

22、9913 kmol/h i-C4 0.1368 i-C5 0.0789 C5 0.7693,关于物料衡算不能平衡的几个情况,精馏塔,1,2,3,100 kmol/h C3 0.20 i-C4 0.30 i-C5 0.20 C5 0.30,62.0087 kmol/h C3 0.3133 i-C4 0.40 i-C5 0.2742,37.9913 kmol/h i-C4 0.1368 i-C5 0.0789 C5 0.7693,针对查定的实际运行结果，进行物料衡算：,C3的衡算： 2019.4273，误差为2.86,C5的衡算： 3029.2267，误差为2.58,问题：为什么衡算不起来，衡算误

23、差这么大？,其自由度104510,其物料衡算计算的所有方程：,原设计的物料衡算情况：,理论衡算结果如下：,精馏塔,1,2,3,100 kmol/h C3 0.20 i-C4 0.30 i-C5 0.20 C5 0.30,61.6667 kmol/h C3 0.3243 i-C4 0.40 i-C5 0.2757,38.3333 kmol/h i-C4 0.1391 i-C5 0.0783 C5 0.7826,精馏塔,1,2,3,100 kmol/h C3 0.20 i-C4 0.30 i-C5 0.20 C5 0.30,62.0087 kmol/h C3 0.3133 i-C4 0.40 i-

24、C5 0.2742,37.9913 kmol/h i-C4 0.1368 i-C5 0.0789 C5 0.7693,实际投产运行查定结果如下：,1、理论衡算结果和实际运行查定结果之间有差别； 2、实际运行查定结果衡算不起来。,原因何在？,精馏塔,1,2,3,100 kmol/h C3 0.20 i-C4 0.30 i-C5 0.20 C5 0.30,C3 i-C4 0.40 i-C5 C5 0.0125,C3 0.015 i-C4 i-C5 C5,设计要求： 1、进料中85的i-C5 从塔顶馏出； 2、组成均为摩尔分率；,其自由度12（102）47（52）10,其物料衡算计算的所有方程：,实

25、际运行过程的物料衡算情况：,精馏塔,1,2,3,100 kmol/h C3 0.20 i-C4 0.30 i-C5 0.20 C5 0.30,62.0087 kmol/h C3 0.3133 i-C4 0.40 i-C5 0.2742 C5 0.0125,37.9913 kmol/h C3 0.015 i-C4 0.1368 i-C5 0.0789 C5 0.7693,严格的物料衡算计算结果如下：,原因之一： 忽略某个或某些流股中含量较少的组份，导致衡算不起来或衡算误差。 但通过严格流股组份分析（包括含量极少的组份），衡算计算就没有问题了， 而且更接近实际工程情况。,原因之二：因为化工过程的复杂性，某些流股无法实现全组份分析，导致衡算计算的误差大。,原因之三：由于对化工过程研究不透彻，例如有意或无意忽略一些付反应过程，导致流股组份分析不全，也会导致衡算计算的误差大。,尤其是制药工艺，有些反应合成途径十分复杂，有时根本不知道付反应的反应式。这也是FDA监管一些药物时，要监管到药物的生产过程的原因。,