甲醇合成工段工艺设计(仅供参考).doc

一、 甲醇的国内外发展简况及甲醇的用途甲醇作为极其重要的有机化工原料，是碳一化学工业的基础产品，在国民经济中占有重要地位。长期以来，甲醇都是被作为农药，医药，染料等行业的工业原料，但随着科技的进步与发展，甲醇将被应用于越来越多的领域。1. 生产的发展1) 世界甲醇工业的发展总体上说，世界甲醇工业从90年代开始经历了19911998的供需平衡，19981999的供大于求，从2000年初至今的供求基本平衡三个基本阶段。1据Nexant Chem Systems公司的最新统计，全球2004年甲醇生产能力为42265万ta2以下是近几年的甲醇需求统计。表1 全球主要地区甲醇消费构成2001年2002年2003年2004年按用途分甲醛940(31)970(32)1010(32)1050(33)MTBE830(28)810(26)780(25)760(22)（其中美国）470(16)430(14)340(11)270(8)醋酸270(9)290(9)300(10)310(10)MMA90(3)90(3)100(3)100(3)其它880(29)900(29)930(30)970(30)需求合计3020(100)3060(100)3100(100)3180(100)按地区分亚洲920(30)940(31)990(32)1040(33)北美1000(33)1000(33)980(31)970(30)西欧630(21)640(21)650(21)670(21)其它470(16)480(16)490(16)500(16)需求合计3020(100)3060(100)3110(100)3180(100) 从上表可以看出，到2004年为止，甲醇仍主要用于制造甲醛和MTBE。用于制造甲醛的甲醇用量随年份成增长趋势，而MTBE的需求量则逐年降低。亚洲需求量增长比较迅速，与此相反，北美地区需求则在减少。2) 我国甲醇工业发展我国的甲醇工业经过十几年的发展，生产能力得到了很大提高。1991年，我国的生产能力仅为70万吨，截止2004年底，我国甲醇产能已达740万吨 ，ll7家生产企业共生产甲醇44065万 生产甲醇，2005年19月甲醇量达39322万t，因此增长222，进口量991万t，因此下降31。22. 生产技术的发展1) 装置大型化 与上世纪末相比，现在新建甲醇规模超过百万吨的已不在少数。在20042008年新建的l4套甲醇装置中平均规模为134万ta，其中卡塔尔二期工程项目高达230万ta。最小规模的是智利甲醇项目，产能也达84万t/a ，一些上世纪末还称得上经济规模的6O万t/a 装置因失去竞争力而纷纷关闭。2) 二次转化和自转化工艺 合成气发生占甲醇装置总投资的5060，所以许多工程公司将其视为技术改进重点。已经形成的新工艺主要是Synetix(前ICI)的先进天然气加热炉转化工艺(AGHR)、Lurgi的组合转化工艺(CR)和Tops e的自热转化工艺(ATR)3) 新甲醇反应器和合成技术大型甲醇生产装置必须具备与其规模相适应的甲醇反应器和反应技术。传统甲醇合成反应器有ICI的冷激型反应器、Lungi的管壳式反应器、Topsdpe的径向流动反应器等 近期出现的新合成甲醇反应器有日本东洋工程(TEC)的MRFZ反应器等，而反应技术方面则出现了Lurgi推出的水冷一气冷相结合的新流程。4) 引入膜分离技术的反应技术通常的甲醇合成工艺中，未反应气体需循环返回反应器，而KPT则提出将未反应气体送往膜分离器，并将气体分为富含氢气的气体，前者作燃料用，后者返回反应器。5) 液相合成工艺 传统甲醇合成采用气相工艺，不足之处是原料单程转化率低、合成气净化成本高、能耗高。相比之下，液相合成由于使用了比热容高、导热系数大的长链烷烃化合物作反应介质，可使甲醇合成在等温条件下进行。33. 甲醇的用途甲醇是重要的化工原料，甲醇主要用于生产甲醛，其消耗量约占甲醇总量的3040，其次作为甲基化剂，生产甲胺、甲烷氯化物、丙烯酸甲酯、甲基丙烯酸甲酯、对苯二甲酸二甲酯等，甲醇羰基化可生产醋酸、醋酐、甲酸甲酯、碳酸二甲酯等。从甲醇低压羰基化生产醋酸，近年来发展很快。随着碳一化工的发展，由甲醇出发合成乙二醇、乙醛、乙醇等工艺正在日益受到重视。甲醇作为重要原料在敌百虫、甲基对硫磷、多菌灵等农药生产中，在医药、染料、塑料、合成纤维等工业中有着重要的地位。甲醇还可经生物发酵生成甲醇蛋白，用作饲料添加剂。 甲醇不仅是重要的化工原料，而且还是性能优良的能源和车用燃料。它可直接用作汽车燃料，也可与汽油渗合使用，它可直接用于发电站或柴油机，或经zsM5分于筛催化剂转化为汽油，它可与异丁烯反应生成甲基叔丁基醚，用作汽油添加剂。4二、 甲醇合成方法 当今甲醇生产技术主要采用中压法和低压法两种工艺，并且以低压法为主，这两种方法生产的甲醇约占世界甲醇产量的80以上。 高压法：（19.629.4MPa）是最初生产甲醇的方法，采用锌铬催化剂，反应温度360400，压力19.629.4MPa。高压法由于原料和动力消耗大，反应温度高，生成粗甲醇中有机杂质含量高，而且投资大，其发展长期以来处于停顿状态。 低压法（5.08.0MPa）是20世纪60年代后期发展起来的甲醇合成技术。低压法基于高活性的铜系催化剂，其活性明显高于锌铬催化剂，反应温度低（240-270）。在较低压力下可获得较高的甲醇收率，且选择性好，减少了副反应，改善了甲醇质量，降低了原料的消耗。此外，由于压力低，动力消耗降低很多，工艺设备制造容易。中压法（9.812.0MPa）：随着甲醇工业的大型化，如采用低压法势必导致工艺管道和设备较大，因此在低压法的基础上适当提高合成压力，即发展成为中压法。中压法仍采用高活性的铜系催化剂，反应温度与低压法相同，但由于提高了压力，相应的动力消耗略有增加。比较以上三者的优缺点，以投资成本，生产成本，产品收率为依据，选择中压法为生产甲醇的工艺路线。三、 甲醇合成的生产原理及影响因素1. 合成系统中的反应 CO2 + H2 CH3OH +102.37kJ/mol CO + H2 (CH3)2O+H2O +200.39kJ/mol CO + 3H2 CH4+H2O +115.69kJ/mol CO2 + H2 CO+H2 -42.92kJ/mol CO + H2 C4H9OH+3H2O +49.62kJ/mol2. 甲醇合成生产的影响因素1) 温度用来调节甲醇合成反应过程的工艺参数中，温度对于反应混合物的平衡和速率，都有很大影响。由一氧化碳加氢生成甲醇的反应和由二氧化碳加氢生成甲醇的反应，均为可逆的放热反应。i) 一氧化碳加氢合成甲醇的最佳温度对于可逆放热反应而言，存在最佳反应温度曲线。当甲醇含量较低时，由于平衡的影响相对的很小，最佳温度高，随着反应的进行，甲醇含量升高，平衡影响增大，最佳温度就低。反应器操作时，沿着最佳温度线进行，则反应速率最高而可得最大产量，但温度不能超过催化剂耐热允许温度，对于铜基催化剂一般不超过300。最佳温度值受操作压力，反应气体混合物的初始组成以及所用催化剂的颗粒大小等因素影响。ii) 有二氧化碳参加反应时的最佳温度 由于两个反应的速率常数和平衡常数不同，同一组成下，两个反应的最佳温度不同。故使两个反应速率之和最大时的温度，才是该反应系统的最佳温度。实现最佳温度的方法，工业上采用连续换热式或多段冷激式催化剂筐，使催化床的温度尽可能沿最佳温度分布。2) 压力 压力也是甲醇合成反应过程的重要工艺条件之一。甲醇合成反应时分子数减少，因此增加压力对平衡有利，由于压力高，组分的分压提高，因而催化剂的生产强度也提高。操作压力的选用与催化剂的活性温度范围有关。对于锌铬催化剂，其起始活性温度在320，由于反应平衡的限制，只能选用2530mpa。再高的压力则会导致副反应加速，造成温度猛升而损坏催化剂。而铜系催化剂的活性温度范围在200300，操作压力可降至5MPa，若能解决移热问题，则可设计高压下铜基催化剂合成甲醇的系统。对于现有合成塔的操作，催化剂使用初期，活性好，操作压力可较低;催化剂使用后期，活性降低，往往采用较高压力，以保持一定的生产强度。3) 原料配比 实验证明，最大反应速率时的组成并不是化学计量组成，由动力学方程确定的最佳氢碳比表明，在远离平衡的情况下，ZnCr催化剂最大反应速率的组成H2：CO4，据报道，氢含量还可以减少副反应以及降低催化剂的中毒程度。以上所指组成的氢碳比是合成反应器内的组成，近似于循环气组成，至于合成系统的进口新鲜气组成，则应根据整个合成系统的物料平衡而定，以维持系统的稳定生产。CO2的存在对于甲醇合成是有益的，一般维持在26即可，新鲜气中CO2含量还可略高些。新鲜气的组成，一般受上游流程的制约，但也要尽量满足本系统的要求。按化学计量要求（H2-CO2）与(CO+CO2)的摩尔比一般为2.052.15。惰性气体的含量也能影响反应速率，含量太高，降低反应速率，生产单位产量的动力消耗也大；维持低惰性气含量，则放空量加大，多损失有效气体。一般来说，惰性气体含量要根据具体情况来定，而且这也是调节工况的手段之一。4) 催化剂颗粒尺寸 由动力学研究可知，催化剂颗粒大小对甲醇合成的宏观速率有显著的影响。催化剂颗粒小，内表面利用率大，从而宏观反应速率大，可减少催化剂用量；粒度减小，会使床层压力降增大，从而增加动力消耗。因此催化剂的最佳颗粒尺寸尚需视气流和床层的特性即有关的具体情况而定。较合理的情况是，反应器上部装小颗粒，下部则装大颗粒催化剂。5) 空速 对于甲醇合成过程，若采用较低空速，则反应速率变化较大，反应过程中气体混合物的组成与平衡组成较接近，反应速率较低，催化剂的生产强度较低，但单位甲醇产品所需循环气量小，气体循环的动力消耗小，预热未反应气体到催化剂进口温度所需换热面积较小，且离开反应器气体的温度较高，热能利用价值高。若采用较高空速，则反应速率变化较小，催化剂的生产强度提高，但增大了预热所需传热面积，热能利用率降低，增大了循环气体通过设备的压力降及动力消耗，并且由于气体中反应产物的浓度降低，增加了分离反应产物的费用。还应注意，空速增大到一定程度后，催化剂床层温度不能维持。因此，必须综合多方面因素来考虑最佳空速的问题。6) 反应器结构 甲醇合成塔内件的型式繁多，内件的核心为催化剂筐，他的设计好坏直接影响合成塔的产量和消耗定额。一个好的催化剂筐设计应满足以下要求：i) 能保证催化剂在升温、还原过程中操作正常、还原充分，尽可能地提高催化剂的活性，达到最大的生产强度。ii) 能有效地移去反应热，合理地控制催化剂层的温度分布，使其逼近最佳操作温度线，提高甲醇净值和催化剂的使用寿命。 iii) 能保证气体均匀地通过催化剂层，阻力小，气体处理量大，甲醇产量高。iv) 充分利用高压空间，尽可能多装催化剂，提高容积利用系数。v) 操作稳定、调节方便，能适应各种操作条件的变化。vi) 结构简单，运转可靠，装卸催化剂方便，制造、安装和维修容易等。vii) 妥善处理各个内件的连接与保温，避免产生热应力，使内件在塔内能自由胀缩。如上所述，合成塔内件型式很多，有各自的特点和适用场合，因此不能用简单的方法肯定一种塔型或全部否定一种塔型。传统的高压法甲醇合成多用连续换热的三套管并流式和单管并流式，中低压法流程则多用多层冷激式合成塔和管式合成塔以及两者的改进型合成塔。无论在何种压力下操作，为减少阻力面应用径向合成塔或轴、径向复合式合成塔。四、 设计依据1. 设计裕量：102. 生产规模：3万吨/年3. 年工作时间：8000小时4. 粗甲醇分离收率：87.55. 甲醇衡算数据：见下表表2 甲醇衡算数据表组分组成合成塔出口气中甲醇含量3.95新鲜气（mol）驰放气（mol）粗甲醇kg/100kg溶解气NM3/tCO20.11640.06800.85224.9280CH40.01960.1560-0.1008N20.0112-0.5824H20.85280.66440.00013.1808CO0.10500.00680.5712CH3OH-0.005679.11400.6944(CH3)2O-0.00030.54800.0560C4H9OH-0.1990-H2O-0.000719.24090.1008五、 工艺流程说明工艺流程可简略的表示为下图： 合成塔循环压缩冷却冷凝甲醇分离闪蒸1新鲜气23驰放气5溶解（闪蒸）气 86粗甲醇 7图1 甲醇合成工艺流程示意图 本工艺采用中压法，新鲜气与由压缩机产生的循环压缩气混合，形成入塔气，其温度为40，压力为10MPa。在进入合成塔前，入塔气与出塔气进行热交换，温度升到150，压力仍为10.0MPa,经铜系催化剂催化反应后，温度变为220，压力变为9.6MPa。出塔气与入塔气热交换后，进入冷却冷凝器，此过程大部分甲醇冷凝，温度降至25。少量未冷凝的甲醇与其它气体一起进入循环压缩机，然后再与新鲜气混合。在进入压缩机前，有一排气口，将少部分冷凝后的气体排放，以减少惰性气体在体系中的累积。冷凝后的甲醇进入闪蒸罐，得到粗甲醇和闪蒸气。闪蒸气排空，粗甲醇进入精馏系统精制。六、 工艺流程设计1. 物料衡算以1吨粗甲醇为基础，对整个系统作总的物料衡算，物料衡算图如图2所示。 第5流股驰放气第1流股新鲜气第8流股溶解（闪蒸）汽第7流股粗甲醇 图2 总物料衡算图CO2CH4N2H2COCH3OH(CH3)2OC4H9OHH2Osum质量/kg/h8.522000.0010.068791.145.481.99192.409999.61摩尔流量/Kmol/h0.194000.00050.00224.7230.1100.02710.68935.746摩尔分率5.42E-03001.40E-056.79E-056.92E-013.07E-037.52E-042.99E-011由甲醇衡算数据表，可计算出第7流股中各组分质量，摩尔流量，摩尔分率表3 第7流股各组分构成情况以CO2为例：m CO2=0.8522108.522kg/hn CO2=8.522/44=0.194kmol/hx CO2=0.194/xini=5.42e-03 表4 第8流股各组成构成情况CO2CH4N2H2COCH3OH(CH3)2OC4H9OHH2Osum体积/NM34.9280.1010.5823.1810.5710.6940.05600.101摩尔流量/Kmol/h0.2200.0050.0260.1420.0260.0310.00300.0050.456摩尔分率0.4820.0100.0570.3110.0560.0680.00500.0101.0001) 新鲜气、驰放气摩尔流量设新鲜气摩尔流量为F，CO的摩尔含量为x。驰放气摩尔流量为Q，其中甲烷含量为y。根据原子平衡，可得以下方程：C原子：F(0.136+x)=25.532+Q(0.1792+y) 1O原子：F(0.2328+x)=36.442+Q0.2476 2N原子：0.0112F=0.026+Q(0.156-y) 3H原子：F(1.784-2x)=121.649+Q(1.354+4y) 4由方程3得Qy=0.026+0.156Q-0.0112F 5将5代入1，2，4，得：0.1472+Fx=25.558+0.3352Q 60.2328+Fx=36.442+0.2476Q 7 0.9144F-Fx=86.4345+0.989Q 8由6，7，8，解得 1.0616F=86.4345+1.3242Q 0.0856F=10.884-0.0876Q 故 F=106.5436 kmol/h Q=20.1421 kmol/h由方程2，3得 x =0.156y=0.0982) 循环比的计算设循环比为R,则循环流量为(R106.5436)kmol/h,得到下列方程：(106.5436R+20.1421)0.0056+36.20160.6838 =0.0395 106.5436R+20.1421+36.2016故R=6.2688再根据上述方法，对每个系统作物料衡算，可得到各流股的摩尔流量，摩尔分率。结果见下表：表5 物料衡算表流股项目CO2CH4N2H2COCH3OH(CH3)2OC4H9OHH2Osum1摩尔流量12.4022.088 1.193 74.235 16.626 0000106.544 摩尔分率0.116 0.020 0.011 0.697 0.156 00001.000 2摩尔流量57.819 67.574 39.900 517.988 86.755 3.740 0.200 0.000 0.468 774.444 摩尔分率0.075 0.087 0.052 0.669 0.112 0.005 0.000 0.000 0.001 1.000 3摩尔流量47.201 67.465 39.900 457.278 72.272 28.607 0.319 0.027 11.176 724.244 摩尔分率0.065 0.093 0.055 0.631 0.100 0.039 0.000 0.000 0.015 1.000 4摩尔流量0.414 0.005 0.026 0.143 0.028 24.754 0.112 0.027 10.694 36.202 摩尔分率0.011 0.000 0.001 0.004 0.001 0.684 0.003 0.001 0.295 1.000 5摩尔流量1.370 1.975 1.167 13.382 2.115 0.113 0.006 0.000 0.014 20.142 摩尔分率0.068 0.098 0.058 0.664 0.105 0.006 0.000 0.000 0.001 1.000 6摩尔流量45.417 65.486 38.707 443.753 70.130 3.740 0.200 0.000 0.468 667.901 摩尔分率0.068 0.098 0.058 0.664 0.105 0.006 0.000 0.000 0.001 1.000 7摩尔流量0.194 0.000 0.000 0.001 0.002 24.723 0.110 0.027 10.689 35.746 摩尔分率0.005 0.000 0.000 0.000 0.000 0.692 0.003 0.001 0.299 1.000 8摩尔流量0.220 0.005 0.026 0.142 0.026 0.031 0.003 0.000 0.005 0.456 摩尔分率0.482 0.010 0.057 0.311 0.056 0.068 0.005 0.000 0.010 1.000 备注：摩尔流量的单位为kmol/h。3) 转化率的计算 单程转化率 X单=(57.819+86.755)-(47.201+72.272)/(57.819+86.755)=0.174 全程转化率 X全=(57.819+86.755)-(47.201+72.272)/(12.402+16.626)=0.8652. 热量衡算1) 入塔气与出塔气的热交换过程热量衡算，下图为热交换过程示意图。 入塔气10.0MPa,40出塔气9.6MPa,22010MPa,1509.6MPa,T图3 入塔气与出塔气的热交换示意图i) 入塔气焓变计算r.g.10.0Mpa,40i.g.10.0Mpa,40r.g.10.0Mpa,150i.g.10.0Mpa,150H1H2H3H查手册5，得到入塔气各组分的物性数据，如下表：表6 入塔气各组分的物性数据组成CO2CH4N2H2COCH3OH(CH3)2OC4H9OHH2O摩尔分率0.075 0.087 0.052 0.669 0.112 0.005 0.000 0.000 0.001 Pc/MPa7.3824.5993.3911.3133.4998.8043.6444.41422.055Tc/K304.19190.56126.133.18132.92512.5466.7563647.13利用公式a= 0.42748R2Tc2.5/Pc b=0.08664RTc/Pc可计算出各组分的a,b,再利用混合规则amix1/2=yiai1/2,bmix=biyi ，可得到混合物的a,b.如下表所示：表7 混合气的a,b的计算组分abyiai1/2biyiCO26.45992.9682E-050.18982.216E-06CH43.22072.9847E-050.15662.604E-06N21.55602.6786E-050.06431.38E-06H20.14271.8203E-050.25271.218E-05CO1.72022.7364E-050.14693.065E-06CH3OH21.73424.5666E-050.02252.205E-07(CH3)2O38.15509.2255E-050.00162.387E-08C4H9OH50.34719.1877E-050.00000H2O14.27272.1136E-050.00231.276E-086.45992.9682E-050.83662.17E-05所以： amix=0.6999, bmix=2.1710-5由RK方程p=RT/(V-b)-a/T0.5V(V+b),求T=313.15K,P=10MPa时的摩尔体积V。解非线性方程的方法用Newton法。即： x(k+1)=x(k)-f(x)/f(x)求解过程如下： f(x)=RT/(Vm-b)-a/T0.5Vm(Vm+b)-pf(x)=RT-1/(Vm-b)2-a/T0.5-1/Vm(Vm+b)(2Vm+b)取初值Vm=2.0010-4m3；迭代计算表如下：Vm(k)Vm(k+1)f(x)f(x)2.00E-043.71E+06-7.34E+102.51E-042.51E-047.97E+05-4.53E+102.68E-042.68E-045.53E+04-3.92E+102.70E-042.70E-043.05E+02-3.88E+102.70E-042.70E-049.39E-03-3.88E+102.70E-04故入塔气的Vm=2.7010-4m3 .由此，可得到：H1=RT-pVm+3a/(2bT0.5)ln(1+b/Vm) 8.314313.15-101062.7010-4+30.6999/（22.1710-5313.150.5）ln（1+2.1710-5/2.7010-4） 119.534J/(mol*h)同理，可求得H3=-1.45J/(mol*h), T=423.15K,P=10MPa下的Vm计算过程如下表：Vm(k)Vm(k+1)f(x)f(x)2.00E-04-8.96E+06-1.03E+112.87E-042.87E-042.89E+06-4.75E+103.48E-043.48E-045.31E+05-3.17E+103.64E-043.64E-042.56E+04-2.87E+103.65E-043.65E-046.53E+01-2.85E+10-3.65E-04-由H2=T1T2CpdT可知，求解H2的关键为找到Cp与T的关系。查表5，得到入塔气各组分的Cp与T的关系： Cp/J/(mol*k)=a0+a1T+a2T2+a3T3+a4T4组分a0a1a2a3a4CO23.2591.36E-031.50E-05-2.37E-081.06E-11CH44.568-8.98E-033.63E-05-3.41E-081.09E-11N23.5385-2.61E-047.40E-081.57E-09-9.89E-13H22.88333.68E-03-7.72E-066.92E-09-2.13E-12CO3.912-3.91E-031.18E-05-1.30E-085.15E-12CH3OH4.714-6.99E-034.21E-05-4.44E-081.54E-11(CH3)2O4.3616.07E-032.90E-05-3.58E-081.28E-11C4H9OH4.4671.64E-026.69E-05-9.69E-083.86E-11H2O4.395-4.19E-031.41E-05-1.56E-086.32E-12利用上表，可求得入塔气的Cp与T的关系。组分摩尔分率a0a1a2a3a4CO20.0747 0.2433 1.012E-041.121E-06-1.772E-097.884E-13CH40.0873 0.3986 -7.831E-043.168E-06-2.973E-099.519E-13N20.0515 0.1823 -1.345E-053.813E-098.109E-11-5.094E-14H20.6689 1.9285 2.462E-03-5.16E-064.625E-09-1.421E-12CO0.1120 0.4382 -4.383E-041.324E-06-1.459E-095.769E-13CH3OH0.0048 0.0228 -3.374E-052.034E-07-2.146E-107.413E-14(CH3)2O0.0003 0.0011 1.570E-067.501E-09-9.265E-123.317E-15C4H9OH000000H2O0.0006 0.0027 -2.527E-068.482E-09-9.442E-123.815E-151.0000 3.2175 1.293E-036.733E-07-1.731E-099.263E-13故入塔气的Cp=3.2175+1.29310-3T+6.73310-7T2+(-1.73110-9)T3+9.26310-13T4H2=3.2175(423.15-313.15)+ 1.29310-3/2(423.152-313.152)+6.73310-7/3(423.153-313.153)+(-1.73110-9)/4(423.154-313.154)+9.26310-13/5(423.155-313.155) 408.661 J/(mol*h)H(H1+H2+H3)F2 = (119.534+408.661-1.45)774.4441000 4.08108J/hii) 出塔气经热交换后的温度 出塔气的H1的求解方法与入塔气相同，利用相同的公式，把流股3的各组分的摩尔分率代入，即可求得 amix=1.0330, bmix=2.25910-5r.g.9.6Mpa,220i.g.9.6Mpa,220r.g.10.0Mpa,Ti.g.9.6Mpa,TH1H2H3H出塔气换热前的摩尔体积Vm(k)Vm(k+1)f(x)f(x)2.00E-041.25E+07-1.20E+113.04E-043.04E-044.52E+06-4.89E+103.96E-043.96E-041.10E+06-2.80E+104.35E-044.35E-041.03E+05-2.30E+104.40E-044.40E-041.08E+03-2.26E+104.40E-04 Vm=4.4010-4m3H1=RT-pVm+3a/(2bT0.5)ln(1+b/Vm) 8.314493.15-9.61064.4010-4+31.0330/（22.25910-5493.150.5）ln（1+2.25910-5/4.4010-4） =33.3 J/(mol*h) 同理，求得出塔气的Cp与T的关系为： Cp=3.2998+8.32410-4T+2.56610-6T2+(-3.58310-9)T3+1.53310-12T4 假设出塔气经还热后温度T=100,求得该温度下的摩尔体积，即可求得H2,H3,即可求得H，比较H与H之间的差值，再利用二分法调整T。当达到所需精度要求后，此时的T即为要求的T。流程图如下所示：设T=100利用RK方程求得Vm求得H2,H3求得HH+H得到T是否，利用二分法重新取TT=100时，Vm(k)Vm(k+1)f(x)f(x)2.00E-04-6.69E+06-8.72E+102.77E-042.77E-041962858.9-4.36E+103.22E-043.22E-042.87E+05-3.18E+103.31E-043.31E-048.15E+03-3.00E+103.31E-043.31E-046.95E+00-2.99E+10-3.31E-04- Vm=3.3110-4m3H2=3.2998(373.15-493.15)+8.32410-4/2(373.152-493.152)+2.56610-6/3(373.153-493.153)+(-3.58310-9)/4(373.154-493.154)+1.53310-12/5(373.155-493.155) =-468.587 J/(mol*h)H3=RT-pVm+3a/(2bT0.5)ln(1+b/Vm)8.314373.15-9.61063.3110-4+31.0330/（22.25910-5373.150.5）ln（1+2.25910-5/3.3110-4） =-158.414 J/(mol*h)H =(H1+H2+H3)F3 =(33.3-468.587-158.414)724.24421000 -4.30108J/h =H+H=-4.30108+4.08108=-2.20107 J/h迭代过程如下：次数T/KVmH1/J/mol*hH2/J/mol*hH3/J/mol*hH/J/hH/J/h/J/h1373.153.311E-0433.3-468.587-158.414-4.30E+084.08E+08-2.20E+072383.153.403E-0433.3-430.389-144.179-3.92E+084.08E+081.59E+073378.153.357E-0433.3-449.508-151.177-4.11E+084.08E+08-2.97E+064380.653.380E-0433.3-439.954-147.649-4.01E+084.08E+086.51E+065379.153.366E-0433.3-445.687-149.759-4.07E+084.08E+088.26E+056378.653.361E-0433.3-447.598-150.467-4.09E+084.08E+08-1.07E+067379.053.365E-0433.3-446.069-149.900-4.07E+084.08E+084.47E+058378.853.363E-0433.3-446.834-150.184-4.08E+084.08E+08-3.11E+05故T=105.7.2) 合成器冷却器的热量衡算从合成塔出来的气体经两段冷却后，只考虑甲醇冷凝，针对此过程作热量衡算，求得各段换热器所需的冷却介质的流量。i) 上段换热器的冷却介质为软水，初温与终温如下图所示：流股3105.7，9.6MPa流股370，9.6MPa软水30软水50HH 由上述计算可知流股3： Cp=3.2998+8.32410-4T+2.56610-6T2+(-3.58310-9)T3+1.53310-12T4 故H3.2998(343.15-378.15)+ 8.32410-4/2(343.152-378.152)+2.56610-6/3(343.153-378.153)+(-3.58310-9)/4(343.154-378.154)+1.53310-12/5(343.155-378.155) -135.403 J/(mol*h) 查表5，得水的Cp与T的关系： Cp=92.053+(-0.039953)T+(-2.1110-4)T2+5.3510-7T3 H=92.053(323.15-303.15)+(-0.039953)/2(323.152-303.152)+ (-2.1110-4)/3(323.153-303.153)+5.3510-7/4(323.154-303.154) =1506 J/(mol*h)因此所需软水的摩尔流量 Q1-(-135.403)724.244210004.714/1506 = 3.07105mol/hii) 下段换热器的冷却介质为循环水，初温与终温如下图所示：流股370，9.6MPa流股325，9.6MPa循环水19循环水30HH 查表得5，甲醇的沸点tb64.7。通过前面的物料衡算，知冷凝后的流股4的甲醇摩尔流量为24.754kmol/h,只考虑甲醇冷凝，可得经下段换热器冷却后的气体组成及其摩尔流量，摩尔分率，从而求得该气体流股的Cp。 流股3经下段换热器冷却后，甲醇气体的摩尔流量变为(28.607-24.754)=3.853kmol/h,其它组分流量不变，计算结果如下：组分摩尔流量/kmol/h摩尔分率a0a1a2a3a4CO247.2010.06750.21999.15E-051.01E-06-1.60E-097.13E-13CH467.4650.09640.4406-8.66E-043.50E-06-3.29E-091.05E-12N239.9000.05700.2018-1.49E-054.22E-098.98E-11-5.64E-14H2457.2780.65371.88492.41E-03-5.05E-064.52E-09-1.39E-12CO72.2720.10330.4042-4.04E-041.22E-06-1.35E-095.32E-13CH3OH3.8530.00550.0260-3.85E-052.32E-07-2.45E-108.46E-14(CH3)2O0.3190.00050.00202.76E-061.32E-08-1.63E-115.84E-15C4H9OH0.0273.84E-051.72E-046.30E-072.57E-09-3.73E-121.49E-15H2O11.1760.01600.0702-6.69E-052.24E-07-2.50E-101.01E-13699.4901.0003.24981.11E-031.17E-06-2.14E-091.04E-12故Cp=3.2498+1.1110-3T+1.1710-6T2+(-2.1410-9)T3+1.0410-12T4 气体流股冷却的摩尔焓变： Ha=3.24