《化工过程数学模型与计算机模拟》课程案例研究.docx

1、化工过程数学模型与计算机模拟课程案例研究化工过程数学模型与计算机模拟课程案例研究 之一甲醇 二甲醚+ 水前言概念设计又称为“预设计”，在根据开发基础研究成果、文献的数据、现有类似的操作数据和工作经验,按照所开发的新技术工业化规模而作出的预设计,用以指导过程研究及提出对开发性的基础研究进一步的要求，所以它是实验研究和过程研究的指南,是开发研究过程中十分关键的一个步骤。概念设计不同于工程设计,因而不能作为施工的依据,但是成功的概念设计不但可以节省大量的人力和物力，而且又可以加快新技术的开发速度，提高开发的水平和实用价值。即使一个很普通的单一产品的生产过程,也可能有1109个方案可供选择。如何从技术

2、、经济的角度把最有希望的方案设计出来,是作为强化研究开发工作的方向，这是一种系统化的分级决策过程,也正是概念设计的真谛。概念设计是设计者综合开发初期收集的技术经济信息,通过分析研究之后。对开发项目作出一种设想的方案,其主要内容包括：原料和成品的规格,生产规模的估计，工艺流程图机简要说明，物料衡算和热量衡算,主要设备的规模，型号和材质的要求,检测方法，主要技术和经济指标，投资和成本的估算，投资回收预测,三废治理的初步方案以及对中试研究的建议。随着计算技术和计算机技术的发展,化工流程过程模拟软件也越来越成熟,计算机辅助设计也日趋广泛。在进行概念设计时，采用流程系统模拟物料衡算和热量衡算，投资和成本

3、估算等问题以及采用流程模拟软件进行整体优化业越来越普遍。本文采用国际上最成功和最流行的过程模拟软件之一的ASPLEN PS作为辅助设计的主要工具。与过程有关的物料和能量的衡算基本上有该软件给出,并从设计流程计算的收敛与否来检验该流程是否可行。本文通过概念设计,其目标是寻找最佳工艺流程（即:选择过程单元以及这些单元之间的相互连接)和估算最佳设计条件。采用分层次决策的方法和简捷设计能消去大量无效益的方案。本文按照以下基本步骤进行设计计算: 1.间歇对连续; 2. 流程图的输入-输出结构; 3.流程图的循环结构； . 分离系统的总体结构; a.蒸气回收系统； b. 液体回收系统。5. 热交换器网络。

4、1 输入信息11市场信息二甲醚工业生产的兴起是同氟氯烷的限制和禁止使用紧密相连的。0年代初国际上气雾剂制品得到了迅速发展,气雾剂生产中，气雾抛射剂主要采用氟氯烷。近年来,发现氟氯烷对地球大气臭氧层有严重的破坏作用,要限制和禁止使用。鉴于二甲醚的饱和蒸汽压等物理性质和二氟二氯甲烷相近,以及其优良的环保性能和无毒,使之成为氟氯烷的理想替代品。自80年代以来,二甲醚作为一种安全的气雾剂得到突飞猛进地发展，目前,气雾剂制品已成为二甲醚最重要的应用市场。二甲醚不仅可以做致冷剂和气雾剂，而且可以做液体燃料。低压下的二甲醚变为液体，与石油液化有相似之处。二甲醚也可以做醇醚燃料,与甲醇按一定比例混合后，可克服

5、单一液态燃料的缺点,从而改善燃料性能，具有清洁、使用方便等优点。据市场调查,二甲醚市场应用前景广阔，国内需求量远远超过供给量。现市场上甲醇价格为1502000元吨，二甲醚价格为00元/吨。以甲醇为原料，经催化脱水得到二甲醚，是一条新兴的工艺。本设计按照概念设计的思路，寻找从甲醇催化脱水生产二甲醚的最佳工艺流程和估算最佳设计条件。1.2反应信息a.反应方程式：CH3O(CH3）O + H2Ob.反应热：HR（2）=-11770/kmc.反应条件：温度T= 25C370C，反应压力11brd.选择性:该反应为催化脱氢，催化剂为12%硅酸处理的无定性氧化铝。在400C以下时，该反应过程为单一、不可逆

6、、无副产品的反应,选择性S=1。e.反应为气相反应。f.甲醇的转化率在以上。g.二甲醚产率:130kml/hh.二甲醚产品纯度:9.5w%i.原料：常温下工业级甲醇. 间歇对连续 选择一个连续的过程，操作费用和物流费用以年为基准,操作时数为810 a。.流程图的输入输出结构和循环结构3净化进料物流原料是工业级的甲醇,内含少量水,和极少量杂质，而水是反应的产物,所以不需要对原料进行净化处理。3.副产物在所选催化剂和反应条件下,主反应的选择性几乎为，副产物极少，在整个设计过程中可以忽略副产物的存在。.3循环和放空由于反应的转化率只有8%，所以需要一股循环物流将未反应的甲醇循环。反应产生的水，含有机

7、杂质很少,可以不经处理而排放,对环境无污染。3.4物料平衡和物流的费用a.物料平衡二甲醚产量PME130kml/,转化率=0.8，选择性S=1,则根据反应方程式有：排放水130kml/,原料甲醇F20ol设循环甲醇量为R km有（260+R）(10.8)=R循环甲醇量=6mo。b.物流费用甲醇价格200元吨，二甲醚价格700元/吨,作为废液水的价格为c.经济潜力EP2=(70001346/100+0-20026030/1000）150=20元/年3. 第二层次的替代方案本方案所生产的二甲醚纯度高达9.5w%，主要用做气雾剂。但二甲醚用做液体燃料和致冷剂时，纯度不需这样,但价格也便宜。根据市场的

8、要求，可以开发同时几种质量要求的二甲醚产品，当然可以在分离上少一些设备。所以这两种方案需要进行评价和比较。.流程的循环结构4.1设计的决策a.只有一个主反应,故只需一台反应器。反应为催化反应，需要催化。反应有大量的反应热放出。采用绝热固定床式反应器。b.有一股循环物流。反应的转化率只有0,还有大量的甲醇没有反应，分离出来的未反应的甲醇需要循环回反应器继续反应。4.2循环的物料平衡在前面，根据转化率和甲醚的产量，已经算出甲醇的循环量为R=65mol/h。4.3反应器的热效应为了作出关于反应器的热效用的决策,首先要估算出反应器的热负荷和绝热的温度变化。这些计算可能提供一些解决反应器热效用的疑难指南

9、。同样,我们也要注意设计问题所隐含的各种温度的限制。根据反应的限制条件,反应催化剂不能在400C的温度之上工作,如果温度超过4000C,主反应的选择性大为降低。而当温度在60370C之间时,反应不仅有接进1的选择性，而且反应速率也较高。此反应为一放热反应,反应热为R11770k/kmo,如果采用绝热反应器,采用出口温度Tot=50,在前面对于转化率x=0.8和相应的物流条件下,由Aspn 模拟反应器，得出反应器进口温度Tin=50。由此可见,可以采用绝热反应器。.4反应器的设计和费用甲醇催化脱水反应器一般采用绝热固定床反应器。在工程上要确定反应器的几何尺寸,首先得确定出一定生产能力下所需的催化

10、剂容积,再根据高径比确定反应器几何尺寸。本文采用实际生产中常见的空速来设计反应器，取7500 h-1 , 于是催化剂容积为：式中 V 催化剂容积,米3 标准状态下甲醇的容积流率,米/小时 U 空速，标米3甲醇/米3催化剂h-知道反应器尺寸和操作条件，由tri的费用关系式计算反应器的费用,如下式:式中 D反应器直径，f H反应器高度，ft FcFm*Fp 取 &S=79进入反应器的物流量为325kol/，则标准状态下甲醇的容积流率为：m/h所需催化剂体积为:V=7280/00.64m绝热固定床反应器有径向的传热传质影响，设计反应器时，为了尽可能避免径向的影响，取反应器的长径比L/D为10。根据一

11、般反应器尺寸,取D.m262f，则催化剂床层高度L=8,此为反应器中催化剂的实际填充高度。反应器的时间高度为：反应器高度L=催化剂高度反应器空闲高度L= + 2 =1 m=32.8ft则反应器费用为:$03106元按投资偿还因子为1/3年计算，则反应器的年度建设费0.391061/30.31元年5分离系统从反应器中出来的气体含有二甲醚、未反应的甲醇、水等物质，它们都是以气体。在进入分离塔之前,要将气体冷却成液体或者气液两相共存。三组分的混合体系，采用两个精馏塔,即一个二甲醚塔和一个甲醇塔来将三种物质分离。51塔序为了清晰的分割了混合物,可通过先回收最轻的组分,也可以先回收最重的组分。当组分数增

12、多时,替代方案数量急剧上升。因而在排定蒸馏塔的塔序时，人们得到了两组推理法则。 表 排定塔序的通用推理法则1.尽快脱出腐蚀性组分2尽快脱出反应性组分或单体3.以溜出物移出产品4.以溜出物移出循环物流，如果它们 是循环送回填料床反应器尤要这样 表2 排定塔序的推理法则1.流量最大的优先.最轻的优先3高收率的分离最后4.分离困难的最后.等摩尔的分割优先6.下一个分离应该是最便宜的在根据上述推理法则，三组分中二甲醚的流量最大，而且也最轻,所以本设计中塔的分离顺序如下图所示:.2二甲醚分离塔二甲醚在常压下的沸点相当低(240C）,如果在常压下进行精馏,则塔顶冷凝器处大量低温冷却介质，必须要增加制冷设备

13、，显然不可取。只有在加压条件下进行精馏,在一定的压力下，可以用常温下的水把精馏塔馏出物冷凝下来。经Asenplu模拟,在压力为0.a下加压精馏，有较好的效果。a.塔的主体设计根据Aspelus的模拟结果（详见附录）,塔采用严格计算有:回馏比R=6358塔板数=进料塔板位置 1冷凝器热负荷Q冷凝器=10424wat再沸器热负荷Q再沸器=274watt实际设计的过程中,取塔板效率E=0，则有实际的塔板数为N实际22/05=4,板间距取0.5，两端共加上2m，则塔的高度为:此塔的横截面积可由下式计算：按塔顶处计算,有：V=（+）D=（0.38+1)130=2261kmolh4Tb70=7.F代入公式

14、可以算出：根据可以计算出塔径：D3f0945由capco可以算出费用为:307071b.冷凝器的设计根据ASPPUS的模拟结果,塔顶冷凝器的热负荷为Q10420wt由此处用室温下的水冷却，可取Tm=0C=68FU可取为100u/h.f2.F则由cacos计算费用为:193$c.再沸器的设计塔底的热负荷为2257wat由取Tm=7F，U=00tu/h.t2.F于是其面积为：由capcot计算费用为：66$.3甲醇分离塔根据Apenpl的模拟结果（详见附录2),塔采用严格计算有:回馏比R=1.32塔板数N=26进料塔板位置1冷凝器热负荷Q冷凝器=16629w再沸器热负荷Q再沸器=15650wat实

15、际设计的过程中,取塔板效率E=05，则有实际的塔板数为N实际6/.=52,板间距取.5m，两端共加上m，则塔的高度为:a.塔主体的设计按塔顶处设计，塔顶处122.7C=52.F,M=30V=（R+1）D=(1.8421）65=1ol/h由计算塔的横截面积为：根据可以计算出塔径：D=2.6f=0.84m由cacost计算费用为：3225$b.冷凝器的设计根据ASPENPS的模拟结果，塔顶冷凝器的热负荷为Q=166290watt由此处用室温下的水冷却,可取T=70oC=58FUC可取为Btu/.ft2.F则由capcost计算费用为:2485$c.再沸器的设计塔底的热负荷为19650at由取Tm=

16、F，00t/h.2.F于是其面积为：由cpcos计算费用为：3$6. 热交换器网络在过程设计中节能总是重要的。所以，普遍采用在反应器和蒸馏塔的周围安装进料和出料的换热器。反应器中的原料供给系统在高于环境温度下操作，反应炉进料需要加热,出料要冷却到塔的进料温度。通过给出了需要加热或冷却的流股，取最小允许温差,将热流股的初、终温度分别减去最小允许温差,与冷流股的初终温度一起排序,这样把原问题划分为多个温度区间。对每个温区进行流股焓平衡计算，以确定净热需求量: 式中 输入到第i个温区的热量； 从第i个温区输出的热量; 温区端点温度; 热容流率。根据温度区间之间热传递的特性,并假定各温区均与外界不发生

17、热量交换,则有： 通过求狭点之上狭点匹配温区热流数和冷流数,以满足H NC若采用两两匹配,则需对热流股进行分割，这里选择多流股换热器来换热。温区净热需求量为负值，在狭点之上使用外部冷却器会使总公用工程消耗增大，为避免使用外部冷却器，将一些流股分出一个冷流股和热流股进行换热。狭点之下温区也采用多流股换热器，不分割热流股。根据温区内流股热量平衡的原则对其它温区子网络进行设计,把所有子网络合并便得到换热网络的初始方案。从子网络结构中挑选那些能构成相同的相邻匹配的子网络组合成换热网络，然后合并相同的匹配，从而减少换热单元数。加热器可从低温部位向高温部位迁移，与处于高温部位的加热器合并。冷却器可从高温部

18、位向低温部位迁移,与低温部位的冷却器合并。在组合过程中需对某些流股的最小传热温差进行松弛。参考文献 美J.道格拉斯著,蒋楚生等译，化工过程的概念设计，化学工业出版社,1994年2 吴指南等编著,基本有机化工工艺学，化学工业出版社,990年3 杨冀宏，麻德贤编著,过程系统工程导论，烃加工出版社，9年 陈甘棠等编著,化学反应工程,化学工业出版社,1990年92年5 APLE LUS L.附录1 甲醇催化脱水生产二甲醚工艺流程图1-增压泵；1-混合器;H1-加热汽化器， 2换热器,R1反应器；C1-冷却器；T1-二甲醚分离塔；V1-降压阀;2甲醇塔；C2冷却器；V2-降压阀 附录2 EPLUS模拟结

19、果FLOWSEET SECIONFLOWSHT CONNECITY Y TEAMTRA SOUCE ET STEA SOURC DST 1 - P1 2 1 M1 3 M1 1 4 H1 H2 7 H2 C1 5 H2 R1 6 R1 H2 8 C1 T1 1 - 0 T1 V 12 T2 M1 13 T2 C2 4 2 V2 V2 - 11 T2 FLWSHEETCONNTIVTY B OKS CK INLET OUTES P1 M1 12 3 H1 3 4 2 64 75 R 6 C1 7 8 T1 9 10 T 11 13 2 13 14 V2 14 15 V1 10 1 CONVEEN

20、CE STTUS MRYTE TRE SUMMAYSREA MAXIM MXIMU VARIALE C ID ERO TOERACE ERR/ ID STA BLCK 12 0504E-16 0.608-09 0.42962E-7 WAE MOLELOW# OLVER01 6 018405 .363505 0.185 DIET-01MOLFOW # $LVER01 = ONVERGED =OT CONVRGEDCONVERGNE BLO： $OVER01 Ta Stam ： 12 Toleranused: .10D-03 0.03 Traeolefac: 0100D-0 0.10D-0 MAX

21、IT30 WI TEATIONS EFORE ACELERATING QMX =00E+0 IN=. ETHOD： WEGTN STATUS：ONVERGED T UMER OF ITERATIONS: 13 FLOWSHEET SECIN CONVERGEE BOCK: $OLVR01 （COTINED） * FAL LUES* ARBLE VALUE PREV VALUE ERR/TL TOTALMOEFLOW KO/HR 0.650496D02 0.54961D+02 0.68029D-11 AL OEFLO KMOL/R 0.3250496D+03 .320330D+03 0.5090

22、D00 ETHA01MOLFLOW KMOL/H .6490+02 0.6490653D+02 0.0000000D+0 DIE-01MOLEFOW KMOH 0.30061+0 0.300861D+00 00000D+00 WTER MLEFLW MOLR 0.1299177-01 0.12991770 .29619D-07 PRESSURE AR 0.73000D+01 0.700000+ 0.00000D+00 MA THALPY JKG -.15655 -.7562D-6394-08 ETHA-1MOEFLOW KML/HR 6498131D+02 .649780+02 0.806D0

23、 DIME0MOFLW MOL 0.300927D+0 0.10861D+03 .5105100 WATER MOEFLWKMOL/H 0.12975D03 0.960+0 0.5090064D+00PESR AR .110000D+0 0.1100000+02 0.000D+00 MAS ENTALPY J/G -.592485D+07 .524885D+07 0.4057939D-03 ITERAON ISTRY * TEARSTREAM: ITEATIO AER/TL STREAM ID VARABLE 1 .000E+05 6 TAL MOLEO 2 0.100E+05 2 ASS E

24、NTHALY 3 208. 6 DMET-01MEFL 4 2008. DIMET-0MOLEFLO 5 331 6 DIT-01MLEFLO 6 334. 12 IMT-01MEFOW 7 6.58 DIE-1MOLOW 8 6.58 IM01LEFO 9 1.2 6 IET-01MOLEFLW 10 122 12 DIMET01MOEFLOW .58 DMET01MOLOW 12 .558 DIE-01LEFL 13 0.510 6 DIE-1MOLFLOWCOMPUTTINA SEQUCE SEUNEUED WAS:P1 $OVE01 M1 C1 T V12 R1 (ETUN $OVR01） CV2 OVERA FLWSHEET ALANE * MAS ADERY BALANCE * CONVENTIONAL IN OUT GENERATO REAVE DIFF.MPONETS （KML/HR ） META-1 26.00 7169E-01-259925 0.758