1、1Chemcad 在化工课程设计中的应用摘要：介绍了化工流程模拟软件 Chemcad 的特点、作用，以精馏过程的工艺设计为实例说明其应用。关键词：Chemcad 软件，设计，应用化工课程设计是学习化工设计基础知识、培养学生化工设计能力的重要教学环节。计算机辅助化工过程设计以成为化工设计的基本手段，有效地利用化工模拟设计软件进行化工设计工作可以极大地提高工作效率。将模拟软件用于过程设计或过程模拟，对于当代的化工过程工程师已是一件很普通的工作。基于此，本院购进美国 Chemstations 公司化工流程模拟软件 Chemcad。本文主要介绍该软件特点、模块功能作用，并以精馏过程的工艺设计实例较详细
2、地对该软件的应用进行说明。1.CHEMCAD 软件概述：CHEMCAD 系列软件是美国 Chemstations 公司开发的化工流程模拟软件。使用它，可以在计算机上建立与现场装置吻合的数据模型，并通过运算模拟装置的稳态或动态运行，为工艺开发、工程设计、优化操作和技术改造提供理论指导。1.1 使用 CHEMCAD 可以做的工作主要有以下几项：A.设计更有效的新工艺和设备使效益最大化 B.通过优化/脱瓶颈改造减少费用和资金消耗 C.评估新建/旧装置对环境的影响 D.通过维护物性和实验室数据的中心数据库支持公司信息系统1.2CHEMCAD 中的单元操作：CHEMCAD 提供了大量的操作单元供用户选择
3、，使用这些操作单元，基本能够满足一般化工厂的需要。 对反应器和分离塔，提供了多种计算方法。ChemCAD 可以模拟以下单元操作：蒸馏、汽提、吸收、萃取、共沸、三相共沸、共沸蒸馏、三相蒸馏、电解质蒸馏、反应蒸馏、反应器、热交换器、压缩机、泵、加热炉、控制器、透平、膨胀机等 50 多个单元操作。1.3 热力学物性计算方法：CHEMCAD 提供了大量的最新的热平衡和相平衡的计算方法，包含 39 种 K 值计算方法，和 13 种焓计算方法。K 值方法主要分为活度系数法和状态方程法等四类，其中活度系数法包含有 UNIFAC 、UPLM (UNIFAC for Polymers)、 Wilson 、T.
4、K. Wilson 、HRNM Modified Wilson 、Van Laar 、Non-Random Two Liquid (NRTL) 、Margules 、GMAC (Chien-Null) 、Scatchard-Hildebrand (Regular Solution)等。焓计算方法包括 Redlich-Kwong 、Soave-Redlich-Kwong 、Peng-Robinson 、API Soave-Redlich-Kwong 、Lee-Kesler 、Benedict-Webb-Rubin-Starling 、Latent Heat 、Electrolyte 、Heat
5、of Mixing by Gamma 等。1.4 使用 CHEMCAD 主要有以下步骤：A.画出您的流程图 B.选择组分 2C.选择热力学模型 D.详细指定进料物流E.详细指定各单元操作 F.运行 G.计算设备规格H.研究费用评估方案 I.评定环境影响 J.分析结果/按需优化 K.生成物料流程图/报告1.5CHENCAD 的应用领域A、蒸馏/萃取(间歇 蒸馏塔,吸收塔,再沸吸收塔,和分馏塔.副产品和副热/冷器也可以出现在 SCDS 的严格模型中,默佛里塔板效率也常在SDCS 中模拟与输出，SCDS 操作塔根据无限层和 5 个输入流,4 个副产品,这里没有限制SCDS 在流程图里的数目.SCDS
6、 强调规格的多样性,比如,总摩尔流比率,热负荷,回流比,蒸出速率,温度,分馏摩尔数,分馏回收,成份比率及产品中部分物流的比率.这些模型能够模拟严格的两相或三相蒸馏系统.如果计算三相系统,用户可以选择将一相液体加入冷凝器,然后回流其它相.SCDS 的主要构成是模拟化学系统中非理想 K 值.它用于 Newton-Raphon 数学收敛和计算并引出严格的方程式,包括 DK/DX(引出主要混合物的 K 值).方面,即有效模拟化工系统.SDCS 的运行时间通常大于其它特殊 TOWR 模组,当组分超过 10 种且较为复杂时.32.2TOWR 简介：TOWR 是一种精确的多级气液平衡模块，它可以模拟任何单塔
7、计算，其中包括蒸馏塔、吸收塔、再沸吸收塔和汽提塔。TOWR 还可以精确的模拟测线产品和侧面供热和降温。TOWR 最多可以处理 5 种进料流和 4 种侧线产品。在工艺流程图中，对 TOWR 的数量没有限制。TOWR 提供了各种技术规范便于用户使用。用户可以指定冷凝器，再沸器或塔板条件。各种技术规范，诸如整体摩尔流动速度，热负荷，回流比，沸腾比，温度，摩尔分率，回收率，相对流动速度，产品的重量和分子量等在 TOWR 中都有明确的说明。通常，TOWR 的收敛速度比主板计算模块 SCDS 要快。2.3 TOWR.PLUS 简介：当所模拟的装置或物流数量较多，就可以使用 TOWR.PLUS。它可模拟带有
8、侧线汽提塔，泵回路，侧线供热、急冷和侧线产品的塔。它是为模拟塔而设计的，但它也可以模拟任何单塔计算，其中包括蒸馏塔，吸收塔，再沸吸收塔和汽提塔。侧线汽提塔，泵回路都看作是 TOWR.PLUS 模块的一部分并且与主塔计算同时进行，不需要再循环计算。冷凝器中自由水的分离可以由用户设定从任意塔板上流出。一个工艺流程图中最多允许 50 个塔单元操作。TOWR.PLUS 提供了一系列技术规范以方便用户使用。用户可以指定冷凝器，再沸器或任何塔板条件。各种技术规范，诸如整体摩尔流动速度，热负荷，回流比，再沸比，温度，体积流速，摩尔流速，质量流速，任何体积分数中的 TBP/D86 温度，溢流，气液比，产品分数
9、，分子质量，气体压力，沸点，闪电和回收率都可以在 TOWR.PLUS 模块中指定。2.4 SHORTCUT 简介：SHORTCUT 蒸馏模块是用 Fenske-Underwood-Gilliland 方法模拟具有一股进料和两股出料的简单蒸馏塔。进料位置可以通过 Fenske 或 Kirkbride 平衡计算得之。但是SHORTCUT 蒸馏模块计算方法不适合塔设计，在这种情况下多用 TOWR 和 TOWR.PLUS模块3 应用举例（一）3.1 设计条件处理量(万吨/年)（8000h）16 （16 /8000 = 2 kg/h）710410进料组成(质量分率) 52% (苯) 48%(甲苯)操作压
10、力 常压进料状态 露点分离要求 塔顶苯含量95% 塔釜苯含量2%完成日期3.2 设计要求根据化工原理课程设计的要求，板式精馏装置的设计应包括以下主要内容：A 设计方案的说明：对所给或选定的整个精馏装置的流程、操作条件和主要设备的型式等进行简要的论述。B 精馏塔的设计计算：确定精馏塔所需的塔板数以及塔的主要尺寸。C 装置的辅助设备，如再沸器、冷凝器等的选择或计算。D 描绘精馏装置的工艺流程图和精馏塔的的设备工艺条件图，编写板式塔精馏装置设计的说明书。3.3 应用软件模拟设计43.3.1 步骤（1） 保存文件*.ccx 点击程序左上角的工具栏“file” ，在出现的菜单中点击“New”，即可出现如
11、下的对话框在对话框中，填入文件名 “chemical”，点击“保存” ，即可出现以下绘画画面（2）画流程图 先进行简捷计算。单击第七行第七列的塔（shortcut 模块） ，然后在左面空白出单击左键，即可画得一精馏塔。同理，画出进料管线“红色箭头”与出料管线“蓝色管线” ，用标有“stream ”的折线连接，如下图 （3）设置单位将“Simulation/Graphic”锁定在 Simulation，下拉菜单“ Format” “Engineering Units” 可出现如下画面：5如选择“Si”制，将“Presure ”设置为“kpa” ，单击“ok ”。（4）选择组分点击工具栏“Ther
12、mophysical” “component list” ，出现以下画面：在“search for”中填入“” ，点击“add ”即可在左面的空白区看见 “ benzene”，表示已经选中第一组分苯。同理选中第二组分甲苯 41，点击“ok” 。 （5）选择热力学模型 点击“thermophysical” “k-values” ，出现画面：如选 SRK 方程，点击 “ok”，即可。同样,点击“thermophysical” “enthalpy” ，选择 SRK 方程，点击“ok”即可。（6）详细定义进料物流 6双击红色箭头后方框内的流股 1，填写数据如下所示： 其中温度、压力和气象组成三者可任选
13、两个填入，注意：不可三者都填。将苯和甲苯的进料流量组成填入，设置完后单击即可。（7）详细定义各单元操作 双击塔 1（设备 1） ，可得如下画面，其中在“Light key split”中填入 ，在xD“Heavy key split”中填入 ，在 R/Rmin 中填入一个倍数，如 1.5。单击 OK。xD（8）运行模拟计算，并显示计算结果按住圆内的，点右键，选择运行：再双击圆内的，出现如下画面，显示计算结果7可知运算的的理论板数为块，进料板是第块，单击 Overall Mass Balance kmol/h kg/hInput Output Input OutputBenzene 133.33
14、8 133.338 10415.602 10415.602Toluene 104.342 104.342 9614.132 9614.132Total 237.680 237.680 20029.734 20029.734Shortcut Distillaton SummaryEquip. No. 1Name Mode 2 Light key component 1.0000Light key split 0.9573Heavy key component 2.0000Heavy key split 0.0427Reflux ratio 1.5000Number of stages 161.0
15、055Min. No. of stages 8.0463Feed stage 81.0028Condenser duty MJ/h -10267.8535Reboiler duty MJ/h 2547.0806Reflux ratio, minimum 1.8834Calc. Reflux ratio 1.5000（9）严格计算简捷计算所得结果可作为严格计算的初值选择第七行第四列的模块CDS column #1，画出流程图8以下操作步骤完全与上述相同，只是第（7）步画面不同于前者，如下：点击 Specifications，例如选择冷凝器 6，填入组分 Specification（ ） ；选择再
16、沸xD器 6，填入组分 Specification（ ） 。xW点击 General ，填入 No.of stages（理论板数）14 ，Feed stages （进料板）7。OK（10）运行 列出物料衡算表Scds Rigorous Distillation SummaryEquip. No. 2Name No. of stages 14 1st feed stage 7Condenser mode 6Condenser spec. 0.9573Cond. comp i 1Reboiler mode 6Reboiler spec. 0.0235Rebl. comp i 1Cond duty
17、MJ/h -16170.7314Reblr duty MJ/h 8458.3379Reflux mole kmol/h 382.3582Reflux ratio 2.79489Reflux mass kg/h 30096.5410Column diameter m 1.9812Tray space m 0.6096No of sections 1No of passes (S1) 1 Weir side width m 0.2159Weir height m 0.0508System factor 1.0000（11）设备尺寸计算 点击 Sizing Trays，在 Select unitop
18、s 中填入 2，OK 。选择塔板类型，Valve Tray（浮阀塔板） ，OK。出现如下画面：如选板间距 0.6m，点击 ok，出现如下数据（只列出第二块板的设计结果）CHEMCAD 5.2.0 Page 1Job Name: chemical2 Date: 04/18/2003 Time: 15:04:35Vapor load is defined as the vapor from the tray below.Liquid load is defined as the liquid on the tray.Equip. 2 Tray No. 2Tray Loadings Vapor Li
19、quid40835.939 kg/h 30067.238 kg/h 14673.854 m3/h 37.018 m3/h Density 2.783 kg/m3 812.229 kg/m3 System factor 1.000Valve type : V-1Valve material : S.S.Valve thickness, gage 12.00010Deck thickness, gage 14.000Tower internal diameter, m 1.981Tray spacing, m 0.610No. of tray liquid passes 1Downcomer di
20、mension, Width m Length m Area m2Side 0.216 1.235 0.182Avg. weir length m 1.235Weir height, m 0.051Flow path length m 1.549Flow path width m 1.755Tray area, m2 3.083Tray active area m2 2.719% flood 69.004Hole area m2 0.517Approx # of valves 436Tray press loss, m 0.104Tray press loss, atm 0.008Dry pr
21、ess drop, m 0.061Downcomer clearance m 0.044Downcomer head loss m 0.006Downcomer backup m 0.184Downcomer residence time, sec 3.265Liquid holdup m3 0.151Liquid holdup kg 122.406Design pressure, atm 1.000Joint efficiency 0.850Allowable stress atm 932.230Corrosion allowance m 0.001Column thickness m 0.
22、000Bottom thickness m 0.0003.3.2 TPXY图 XY图 点击Plot TPXY，出现如下画面：填入组分如上，第一组分苯，第二组分甲苯，大气压 1atm，No. of point 10 点OK ，即可出现如下：11温度组成图苯甲苯组成图（X Y 相图）4 应用举例（二）4.1 设计任务及条件 某厂以苯和乙烯为原料，通过液相烷基化反应生成含乙苯和苯的混合物。经水解、水洗等工序获得烃化液。烃化液经过精馏分离出的苯返回循环使用，而从脱除苯的烃化液中分离出乙苯作为生产苯乙烯的原料。现要求设计一采用常规精馏的方法，从烃化液中分离出苯的精馏装置。要求从塔顶馏出的苯液中，乙苯的含
23、量低于.（摩尔分数，下同） 。釜液中苯含量要求小于.。塔顶馏出液及釜液要求降至 ，其它条件如下所述。C04.1.1 烃化液的流量及组成烃化液参数如下表所示。12烃化液流量及组成表 组分名称 苯 甲苯 乙苯 二甲苯 二乙苯 三乙苯 焦油 水组成摩尔分数 0.5982 0.00602 0.3091 0.00906 0.06278 0.003110.00556 0.00608质量流量（kg/h） 14082.9流量摩尔流量/（kmol/h） 154.63温度/ C0 60 C04.1.2 公用工程条件A 加热蒸汽等级：.Mpa(绝)。 B 循环冷却水： 。0C 电容量可满足需要。4.2 分离工艺流程
24、通过考察实际生产装置，结合设计条件和要求拟定原则，工艺流程方案如图-所示。由- 泵将烃化液从中间罐 -中引出， 送至进料预热器-，在 - 中与釜液换热回收热量被预热至 ，进入 -塔中。-塔所C0需热量由再沸器-加入驱动精馏过程后，其热量由冷凝器- 从塔顶移出，使塔顶蒸汽全部冷凝。凝液经回流泵-一部分送至 -塔顶作为回流，余下部分作为产品送至冷却器-，冷却至 ，然后进入产品罐- 中。 0- 塔排出的釜液经 -与进料换热后，由 - 再次冷却至 ，然C0后进入产品罐-中。4.3 过程模拟在所画流程图中，3 号塔（苯乙苯塔）塔顶分离出苯，塔底分离出乙苯；5 号塔（乙苯二乙苯）顶分离出乙苯，塔底分离出二
25、乙苯流程图13乙苯二乙苯塔模拟结果塔顶馏液 釜液物料名称 进料生产值 Aspen plus 计算值Chemcad 计算值生产值 Aspen plus 计算值Chemcad计算值苯 1.195142E-7 0.00028 1.5E-7 1.153661E-7 4.49E-15 1.809147E-20甲苯 0.001301378 0.00184 5.3E-5 0.001673199 2.32E-9 2.226303E-10乙苯 0.7920002 0.99652 0.99994 0.9960001 0.11548 0.1109 0.07799914二甲苯0.02325865 0.00129 2.
26、51E-6 0.002326426 0.09606 0.1004 0.09652166二乙苯0.1611802 0.0007 5.2E-8 1.078943E-12 0.69276 0.69332 0.7253126三乙苯0.007984558 1.5E-20 1.47715E-21 0.03428 0.034346 0.0359306焦油 0.01427464 8.9E-30 2.636403E-20 0.06142 0.061403 0.06423603组成（摩尔分率）水 1.999823E-7 1.3E-7 2.571199E-7 2.17E-13 1.398E-184.3.2 3 塔结
27、果尺寸设计结果如下（只列出塔顶第一块板的结果） Equip. 3 Tray No. 2Tray Loadings Vapor Liquid10803.227 kg/h 3471.463 kg/h 3343.239 m3/h 4.295 m3/h Density 3.231 kg/m3 808.327 kg/m3 System factor 1.000Valve type : V-1Valve material : S.S.Valve thickness, gage 12.000Deck thickness, gage 14.000Tower internal diameter, m 1.067
28、Tray spacing, m 0.610No. of tray liquid passes 1Downcomer dimension, Width m Length m Area m2Side 0.083 0.570 0.032Avg. weir length m 0.570Weir height, m 0.051Flow path length m 0.902Flow path width m 0.921Tray area, m2 0.894Tray active area m2 0.830% flood 59.010Hole area m2 0.15814Approx # of valv
29、es 133Tray press loss, m 0.073Tray press loss, MPa 0.001Dry press drop, m 0.044Downcomer clearance m 0.044Downcomer head loss m 0.000Downcomer backup m 0.134Downcomer residence time, sec 3.571Liquid holdup m3 0.029Liquid holdup kg 23.157Design pressure, MPa 0.120Joint efficiency 0.850Allowable stres
30、s MPa 94.458Corrosion allowance m 0.001Column thickness m 0.002Bottom thickness m 0.0094.3.3 5 号塔顶板 Equip. 5 Tray No. 2Tray Loadings Vapor Liquid55563.566 kg/h 50590.883 kg/h 17039.765 m3/h 66.503 m3/h Density 3.261 kg/m3 760.725 kg/m3 System factor 1.000Valve type : V-1Valve material : S.S.Valve th
31、ickness, gage 12.000Deck thickness, gage 14.000Tower internal diameter, m 2.286Tray spacing, m 0.600No. of tray liquid passes 1Downcomer dimension, Width m Length m Area m2Side 0.311 1.568 0.335Avg. weir length m 1.568Weir height, m 0.051Flow path length m 1.664Flow path width m 2.064Tray area, m2 4
32、.104Tray active area m2 3.434% flood 72.248Hole area m2 0.652Approx # of valves 551Tray press loss, m 0.11415Tray press loss, MPa 0.001Dry press drop, m 0.065Downcomer clearance m 0.044Downcomer head loss m 0.012Downcomer backup m 0.206Downcomer residence time, sec 3.744Liquid holdup m3 0.237Liquid holdup kg 180.619Design pressure, MPa 0.100Joint efficiency 0.850Allowable stress MPa 94.458Corrosion allowance m 0.001Column thickness m 0.000Bottom thickness m 0.0005 结束语 CHEMCAD 工程模拟软件为辅助教学、化工原理课程设计、毕业设计等提供了可靠的帮助。6 参考文献匡国柱，史启才编化工单元过程及设备课程设计北京：北京化学工业出版社，2002
