环氧乙烷合成模拟与优化(共28页).doc

《环氧乙烷合成模拟与优化(共28页).doc》由会员分享，可在线阅读，更多相关《环氧乙烷合成模拟与优化(共28页).doc（28页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上环氧乙烷合成模拟与优化 1 绪论1.1环氧乙烷的性质和用途。1.1.1环氧乙烷的性质环氧乙烷的分子式为C2H4O，分子量为44.05，简称EO。环氧乙烷常温下沸点10.40C,为无色气体。蒸汽易燃易爆，在空气中爆炸范围为3%100%。环氧乙烷可与水、醇、醚等有机溶剂混合。环氧乙烷的物理性质列于下表：环氧乙烷的重要物理性质性质数值性质数值分子量44.05闪点(开杯试验)/0C-18沸点（101.3KPa）/0C10.4凝固点/0C-112.5体膨胀系数（200C）/0C0.00161熔化热（250C）（KJ/mol）5.17临界压力/MPa7.19燃烧热（250C）（K

2、J/mol）1306.04临界温度/0C195.8常压下在水中溶解热（250C）（KJ/mol）6.3介电常数13.71电离势(实验值)/J(1.08-1.730*10-18(00C)偶极矩/（Cm）6.34*10-30折光指数1.3597在空气中爆炸极限（体积分数/%）上限100下限3同时，环氧乙烷的化学性质非常活泼，能与许多化合物发生开环加成反应。环氧乙烷能还原硝酸银，受热后易聚合，在有金属盐类或氧的存在下能分解。1.1.2 环氧乙烷的用途环氧乙烷是一种有毒的致癌物质，以前被用来制造。环氧乙烷易燃易爆，不易长途运输，因此有强烈的地域性。被广泛地应用于，等行业。在化工相关产业可作为的起始剂。

3、 环氧乙烷可杀灭（及其内孢子）、及，因此可用于消毒一些不能耐受高温消毒的物品。美国化学家Lloyd Hall在1938年取得以环氧乙烷消毒法保存香料的专利，该方法直到今天仍有人使用。环氧乙烷也被广泛用于消毒医疗用品诸如绷带、缝线及手术器具。 环氧乙烷有杀菌作用，对金属不腐蚀，无残留气味，因此可用材料的气体杀菌剂。 通常采用环氧乙烷-（两者之比为90：10）或环氧乙烷-的混合物，主要用于医院和精密仪器的消毒。环氧乙烷用熏蒸剂常用于粮食、食物的保藏。例如，干蛋粉的贮藏中常因受细菌的作用而分解，用环氧乙烷熏蒸处理，可防止变质，而蛋粉的化学成分，包括氨基酸等都不受影响。 环氧乙烷易与酸作用，因此可作为

4、抗酸剂添加于某些物质中，从而降你这些物质的酸度或者使用其长期不产生酸性。例如，在生产氯化丁基橡胶时，与异戊二烯共聚物的溶液在氯化前如果加入环氧乙烷，则成品即可完全不用碱洗和水洗。 由于环氧乙烷易燃及在空气中有广阔的爆炸浓度范围，它有时被用作燃料气化爆弹的燃料成份。 环氧乙烷自动分解时能产生巨大能量，可以作为火箭和喷气推进器的动力，一般是采用和环氧乙烷的混合物（60：40-95：5）。这种混合燃料燃烧性能好，凝固点低，性质比较稳定，不易引爆。总的来说，环氧乙烷的上述这灯直接用途消费量很少，环氧乙烷作为工业衍生物仅次于聚乙烯，为第二位的重要产品。其重要性主要是以其为原料生产的系列产业。由环氧乙烷衍

5、生的下游产品的种类远比各种乙烯衍生物多。环氧乙烷的毒性为的27倍，与氨的毒性相仿。在体内形成、乙二醇和，对中枢神经系统起麻醉作用，对粘膜有刺激作用，对细包原浆有毒害作用。 大部份的环氧乙烷被用于制造其它化学品，主要是乙二醇。乙二醇主要的最终用途是生产聚酯聚合物，也被用作汽车冷却剂及防冻剂。其次用于生产乙氧基化合物、乙醇胺、乙二醇醚、亚乙基胺、二甘醇、三甘醇、多甘醇、羟乙基纤维素、氯化胆碱、乙二醛、乙烯碳酸酯等下游产品。 环氧乙烷主要用于制造乙二醇（制涤纶纤维原料）、合成洗涤剂、非离子表面活性剂、抗冻剂、乳化剂以及缩乙二醇类产品，也用于生产增塑剂、润滑剂、橡胶和塑料等。广泛应用于洗染、电子、医药

6、、农药、纺织、造纸、汽车、石油开采与炼制等众多领域。1.2 环氧乙烷工艺发展历史及其特点环氧乙烷早期采用氯醇法工艺生产，20世纪20年代初，UCC公司进行了工业化生产，之后公司基于Lefort有关银催化剂的研究成果，使用银催化剂，推出空气法乙烯直接氧化生产环氧乙烷工艺。20世纪50年代末，Shell公司采用近乎纯氧代替空气作为生产环氧乙烷的氧原料，推出氧气法乙烯直接氧化生产环氧乙烷工艺，经过不断改进，目前较先进的生产方法是用银作催化剂，在列管式固定床反应器中，用纯氧与乙烯反应，采用乙烯直接氧化生产环氧乙烷。现就这几种方法进行分析比较。1.1 氯醇法环氧乙烷氯醇法生产分两步进行：首先氯气与水反应

7、生成次氯酸，再与乙烯反应生成氯乙醇；然后氯乙醇用石灰乳皂化生成环氧乙烷。这种方法存在的严重缺点大致有：1)消耗氧气，排放大量污水，造成严重污染；2)乙烯次氯酸化生产氯乙醇时，同时副产二氧化碳等副产物，在氯乙醇皂化时生产的环氧乙烷可异构化为乙醛，造成环氧乙烷损失，乙烯单耗高；3)氯醇法环氧乙烷，醛的质量分数很高，约为410-610-6最低也有210-6 。#氯醇法生产环氧乙烷，由于装置小、产量少、质量差、消耗高，因而成本也高，与大装置氧化法生产的高质量产品相比已失去了市场竞争能力。8 h4 Y8 F3 _. _1.2 直接氧化法乙烯直接氧化法，分为空气直接氧化法和氧气直接氧化法。空气直接氧化法。

8、空气直接氧化法用空气作氧化剂，因此生产中必须有空气净化装置，以防止空气中有害杂质带人反应器而影响催化剂的活性。空气法的特点是由两台或多台反应器串联，即主反应器和副反应器，为使主反应器催化剂的活性保持在较高水平(63%-75%)，通常以低转化率操作(20%-50%)。氧气直接氧化法此法是采用制备纯氧或有其他氧源作氧化剂。由于用纯氧作氧化剂，连续引入系统的惰性气体大为减少，未反应的乙烯基本上可完全循环使用。从吸收塔顶出来的循环气必须经过脱碳以除去二氧化碳，然后循环返回反应器，否则二氧化碳质量超过15%，将严重影响催化剂的活性。2 J/ F0 Q8 Z/ e3 * y: D n$ E! N%这两种氧

9、化方法均采用列管式固定床反应器。反应器是关键设备，与反应效果密切相关，其反应过程基本相同，包括反应、吸收、汽提和蒸馏精制等工序。但是氧气氧化法与空气氧化法相比前者具有明显的优越性，主要体现在以下几个方面：! 1 + k1 v) O! a/ h0 1)流程：空气法需要空气净化系统和二次反应器与吸收塔等，以及尾气催化转化器与热量回收系统；氧气法需要分离装置和二氧化碳脱除系统。氧气氧化法与空气氧化法相比，前者工艺流程稍短，设备较少，建厂投资少。 2)催化剂：生产中催化剂的用量直接影响到产品成本的高低。影响催化剂使用量多少主要有两个因素。一方面是催化剂的性能，如催化剂的选择性及转化率等。氧化法在这方面

10、比较有利，因为选择性高，催化剂需要量少。另一方面是催化剂中银的质量分数，德国Hals公司空气法催化剂银的质量分数约20%；而Shell公司氧气法催化剂银的质量分数仅为10%左右，同样体积的催化剂；含银量几乎相差50%。# Y* . k. j; Y R+ y2 c 3)反应器：在同样生产规模的前提下，氧气法需要较少的反应器，而且，反应器都是并联操作。空气法需要有副反应器，以及二次吸收和汽提塔等，增加了设备投资。收率和单耗：氧气法环氧乙烷收率高于空气法，而且乙烯单耗较低。7 j- % eR4 u. B. n7 i2 t8;综上所述，氧气氧化法无论是在生产工艺、生产设备、产品收率、反应条件上都具有明

11、显的优越性，因此目前世界上的EO/EG装置普遍采用氧气氧化法生产。但是由于氧气氧化法采用纯氧作原料，因此在氧气价格上涨时，对氧气法生产环氧乙烷的生产总费用会有一定的影响，而空气法就不存在氧气价格变动带来的总费用变动问题，而且氧气法对原料的纯度要求很高，如氧气纯度低，就会显著增加含烃放空气体的数量，造成乙烯单耗提高。尽管如此，通常氧气氧化法的生产成本要比空气氧化法低10%左右。1.3我国环氧乙烷生产与消费状况2008年，我国大陆共有 1 4家公司的 1 5套环氧乙烷装置开工，产能达到 2 3 0万 t ， 年，实际产量达到 210万t左右。从环氧乙烷原料路线上看，我国环氧乙烷以乙烯为原料，酒精-

12、乙烯路线在发展。从环氧乙烷商品能力和分布上看：2008年，我国环氧乙烷商品能力近80万t，实际商品量近60万t。从环氧乙烷商品制造商看：我国商品环氧乙烷能力约80出自于中国石化和中国石油两大公司。其中，上海金山石化公司53万t，南京扬子石化公司l1万t。2009年，环氧乙烷，乙二醇、PTA和聚酯产业前景不容乐观，全球环氧乙烷乙二醇增长看好。但是，中国作为全球最大的乙二醇消费国，已经受到中东2008年新增的280万t年的能力投运以及对中国市场销售的巨大挑战。因此，中国环氧乙烷的产品结构开始发生重大调整，中国的乙二醇装置将逐步转向环氧乙烷商品的生产，中国的环氧乙烷商品资源将更加充裕。 2环氧乙烷合

13、成工艺流程 环氧乙烷制备工艺流程是合成环氧乙烷中续原料气制备、净化之后的部分流程。从合成气开始包括环氧乙烷的合成、环氧乙烷的分离，而且带有物料循环回路。工艺流程图如下图所示。整个工艺流程操作稳定，可以认为是动态模拟系统，系统中各个参数均不随时间变化。2.1环氧乙烷合成工艺模拟环氧乙烷的制备分离装置环氧乙烷精制二氧化碳脱除环氧乙烷合成单元建立了一个合成塔的经验模型，通过入塔气体的流量与组成，计算得到出塔的环氧乙烷浓度；合成塔物料衡算单元是通过物料衡算在入塔气体流量、组成和出塔气体流量、组成，T、P已知的情况下，计算经过反应后出塔气体的流量及其余组分的浓度。环氧乙烷的分离单元包括环氧乙烷的汽提和再

14、吸收，二氧化碳的脱除，通过相平衡（气液平衡常数），物料平衡得到气体和液体的组成流量和循环比。我们只对合成过程进行分析.2.2工艺流程的分析 该模拟过程主要是对合成环氧乙烷反应器进行模拟。本模拟过程采用一维管式反应器模拟。在催化床轴向取高度为 d 的微元圆柱体, 通过物料衡算和热量衡算, 可以得到固定床反应器中 ET摩尔分、 二氧化碳摩尔分数和反应温度随床层轴向分布的微分方程组。使用龙格-库塔沿着轴向积分这三个方程。模拟出口物料的浓度随催化剂剂龄与入口温度的变化，从而找寻最佳的反应温度3.1环氧乙烷的合成流程模拟与主要参数3.1.1乙烯氧化合成环氧乙烷管式反应器的模拟目前，世界上ET 氧化合成E

15、O 反应器全部采用固定床列管式反应器，管内填装催化剂，管间流动的是处于沸点的冷却液。对管式反应器采用一维拟均相平推流模型进行模拟， 不考虑径向的含量差和温度差。反应器模型选用一维拟均相平推流模型，其依据如下：（l）固定床反应器中，雷诺准数Re 40，床层高度L / 催化剂颗粒当量直径ds大于l00 时，床层轴向返混可忽略；管内径dt / 颗粒直径dp大于8 时，可不计壁效应，可假定反应器中同一截面气流是均匀的7。本文模拟的反应器符合上述简化条件。（2）甘霖6采用二维非均相数学模型对装填YS- 6 型银催化剂的工业环氧乙烷合成反应器进行了模拟计算。从其结果看，气相径向浓度差、固相表面与气相浓度差

16、均非常小；气相径向温差小于0. 5 C；除反应器入口处的床层外，固相表面径向温差不高于6 C；除反应器入口处的床层外，同一处的固相表面与气相温度差不高于6 C，且大部分在4 C以下。（3）赖海明8、包力等2采用一维拟均相模型对工业环氧乙烷合成反应器进行了模拟计算，模拟结果与实际结果吻合很好。齐润红等9采用一维拟均相模型和二维拟均相模型分别对工业环氧乙烷合成反应器进行了模拟计算，模拟结果均与实际结果吻合很好。由于本文模拟的固定床反应器内的压降很小，可以忽略，从而假定反应器内的反应压力不随轴向而变化，为恒定值。4.1.2环氧乙烷合成的模拟计算 dy(ET)/dl=r2COR2+(1-0.5y(ET

17、)r1COR11-0.5y(ET)-0.25y(CO2） bA/1-0.5y0(ET)-0.25y0(CO2)NT,0，d y(CO2)/dl =2r 2COR2+0.5 y(CO2)r1COR1 1- 0.5 y(ET)-0.25y(CO2)b A/(1-0.5y0(ET)-0.25 y0(CO2) NT,0，dT/dl=-Hm,1bAr1COR1-Hm,2bAr 2COR2-Ktbmtdt（T-Tf)/(GCP,m)其中：r1=-dN1(ET)/dm(cat)= k1p(O2)p(ET)/1+K1p(O2)+K2p1/2(O2)p(CO2),r2=-d N2(ET)/dm(cat) =k2

18、p1/2(O2)p(ET)/1+K1p(O2)+K2p1/2(O2)p(CO2)k1=exp (12.39103-43.58572103/RgT)，k2=exp (18.59960-77.76320103/RgT)，K1=exp (-2.82805+18.32100103/RgT)，K2=exp (-1.00740+34.66055103/RgT)。A=1-0.5y(ET)-0.25y(CO2)/1-0.5y0(ET)-0.25y0(CO2) y(O2)=Ay0(O2 )+1.25y0(CO2)-0.5y0(ET)-1.25y(CO2)+0.5 y(ET)y(EO)=Ay0(EO)+y 0(E

19、T)+0.5y0(CO2)-0.5y(CO2)-y(ET)y(H2O)=Ay0(H2O)-y0(CO2)+y(CO2)y(CH4)=Ay0(CH4) 主副反应活性矫正系数催化剂剂龄的关系t/月COR1COR220.74270.411240.96150.534280.90090.4859100.85510.4801120.81040.4344140.77350.4187160.72100.3722180.69360.3503200.64560.3183将不同时期的催化剂活性校正系数回归成与催化剂龄的关系如下：COR1=0.7916+0.02220t-0.t2（R=0.8535）COR2=0.43

20、59+0.01365t-0.t2（R=0.6025）b=0.5801g/cm3=580.1kg/m3组分TC/KPC/MPaVc/(cm3/mol)Zc/DC2H404697.1940.2001400.2581.9CH2=CH2282.45.0360.0851290.2760CH4190.64.6000.008990.2880H20647.322.050.344560.2291.8CO2304.27.3760.225940.2740由于入口气体中甲烷和乙烯气体的比例占到了百分之八十五以上，故将入口气体简化为二组分混合。计算Kay规则适用判据： 均小于2，大于0.5。可以同Kay规则作为混合规则

21、。Tpc=0.294*282.4+0.706*190.6=217.6KPpc=0.294*5.036+0.706*4.600=4.728MPa=0.294*0.085+0.706*0.008=0.0306虚拟对比参数为Tpr=T0/Tpc=500/217.6=2.30Ppc=P0/Ppc=2.0/4.728=0.423用普遍化压缩因子法，查的Z0=0.9943 ， Z1=0.0354则：Z=Z0+Z1=0.9943+0.0306*0.0354=0.9954Vmol=0.9954*8.3145*500/2*106=2.069m3/kmolVR=m/b=31500/580.1=54.3m3Qv=S

22、v/VR=4500/54.3=82.87m3/hNt0=Qv/V=82.87/2.069=40.05kmol/hw=5.1(dt/dp)0.8(dt/l)0.1(dpG/m)0.46m/dtdt=3.1310-2mdp=8.310-3ml=8mmt=8821SR=VR/HR=54.3/7=7.76m2u=Qv/SR=82.87/7.76=10.68m/h=2.9710-3m/sii=i=0.089Vci1/3ij=ij(ij)1/2Mij=0.2942（0.8091291/3）3+0.2940.706（0.8091291/3）(0.809991/3)3/2+0.7062(0.809991/3)

23、3=44.45=0.2942（282.4/1.2593）（0.8091291/3）3+0.2940.706（282.4/1.2593）(190.6/1.2593)1/2（0.8091291/3）(0.809991/3)3/2 +0.7062(190.6/1.2593)(0.809991/3)3/44.45=7566.46/44.45=170.22Mm=0.2942（282.4/1.2593）（0.8091291/3）2281/2+0.2940.706（282.4/1.2593）(190.6/1.2593)1/2（0.8091291/3）(0.809991/3) 20.361/2+0.7062(

24、190.6/1.2593)(0.809991/3)2161/2 /（7566.46/44.451/3）2=17.28kg/kmolG=um=10.682.110317.28/(0.99548.314500)=93.66kg/(m2h)从以上计算中发现气体组分性质接近于CH4性质，且CH4含量大于50%，所以在下面的计算中，忽略其他组分，将气体看成纯甲烷系统。Tc=190.6K Pc=46.00bar Vc=99cm3/molM=16g/mol=(Tc/M3)1/6=(190.6/(163464)1/6=0.0884(P)-1V=Tr=T/190.6=5.246610-3TPr=2.1/4.6=

25、0.4565Q=(5.244610-3T)-3/2=12990.58776T-3/2=0.08664=7.538T-1Zk+1=1+-q=1+7.538T-1-97923.0T-5/2V=4.15710-6ZTm3/mol=4.157ZT cm3/molr=Vc/V=99/4.157ZT=23.82Z-1T-1(-0)+11/4=1.0230+0.23364r+0.58533Tc=190.6K Pc=46.00bar M=16 w=0.008=2.57/k=190.6（0.7915+0.16930.008）=150.12KT*=T/(/)=6.6610-3T Kv=m=10-63600=3.6

26、10-3 kg/(m.h)Re=dpG/mCp,m=A+BT+CT2+DT3=19.2380+0.T+0.T2-11.309410-9T3mer=0.38Re0.72(dt/dp)0.60m 对于CH4,在273,15K下，0=26.410-3kcal(mh)-1n=1.330=26.410-3=1.5210-5T1.33采用Stiel-Thodos法，可得到高压下的，计算方法如下：r=Vc/V=99/4.157ZT=23.82Z-1T-1当r0.5时，（-0）=1.2210-2exp(0.535r)-1当0.5r2.0时，（-0）=1.1410-2exp(0.67r)-1.069当2.0r2

27、.8时，（-0）=2.6010-3exp(1.155r)-2.016=210=210=156.96V=4.15710-6ZTm3/mol=4.157ZT cm3/molr=Vc/V=99/4.157ZT=23.82Z-1T-1Hm，l = -9.9900104-25.337 T +3.277210-2T2 -l.6936 10-5T3 +3.4l7 10-9T 4 Hm，2 = -l.327l03106 +l5.926 T -2.92310-3T2 -4.99l510-5T3 +l.036210-8T4 Tf=485.5K3.1.2环氧乙烷合成流程图NoYes输出y（EO），TYest20t=

28、t+2y（O2）0输入塔设备常量，气体入口条件，T0，t，以及运算步长利用龙格-库塔法计算轴长度为l时y（ET），y（EO），y（CO2），y（O2）TT=T+20T500The endYesNo3.2环氧乙烷合成的主要参数该工厂的EO 合成反应器的基本参数和操作条件如下：1)固定床中催化剂填充高度7 m；2)反应管内径31.3 mm，壁厚3.4 mm，每台反应器8 821 根列管；3)空速4 500 h-1，反应压力2.00 MPa，气包温度227.0 ；4)反应器入口各组分摩尔分数依次为：y(C2H4)=25.07%，y(O2)=78.0%，y(H2O)=0.2%，y(CH4)=60.21

29、%，y(CO2)=6.74%。Cp,m气体热容，kJ/(kmolK)；dp催化剂颗粒等外表面积当量直径，m；dt列管内径，m；G反应气质量流量，kgm- 2h- l；Hm,l反应（l）反应热，kJ/kmol；Hm，2反应（2）反应热，kJ/kmol ；kl反应（l）吸附平衡常数，MPa- l；k2反应（2）吸附平衡常数，MPa- l. 5；Kl反应（l）反应速率常数，kmolkg-lh-lMPa-2；K2反应（2）反应速率常数，kmolkg-lh-lMPa-l.5；l床层高，m；N物质的量，kmolh- l；Nt物质总量，kmolh- l；Nt0物质总量初值，kmolh- l；Pi组分分压，M

30、Pa；Rg气体常数，8. 3l4 kJkmol- lK- l；rl反应（l）反应速率，kmolkg- lh- l；r2反应（2）反应速率，kmolkg - lh - l；S列管内截面积，m2；T床层温度，上标0 表示初值，K 或；Ktb传热总系数，kJm- 2K - lh - l；u气体表观流速，mS- l；Wcat催化剂质量，kg；yi反应气体中组分摩尔分数；yi0反应气体中组分摩尔分数初值；m气体混合物导热系数，kJm-lK-lh-l；m气体混合物粘度，kgm- lh- l；b催化剂堆密度，kgm- 3；g气体密度，kgm- 3。4计算结果模拟催化剂的不同使用时间、入口的不同温度下所对应的

31、出口环氧乙烷浓度tT30032034036038040042044046048050020.0.0.04260.0.0.0.0.0.0.0.40.0.0.0.0.042710.0.0.0.0.0.60.0.0.042970.0.0.0.0.0.0.0.0252180.0.0.0.0.0.0.0.0.0.0.100.0.0.043120.0.0.0.0.0.0.0.120.0.0.0.0.043190.0430.0.0.0.0.140.0.0.0.0.0.0.0.0.0.0.160.0.0.0.0.0.0.0.0.037170.0.180.0.044080.0.0.0.0.0.0.0.0.200

32、.0.044750.0.0.0.0.0.0.0.0.在入口温度400K下，随着轴长L的变化轴面温度T及环氧乙烷轴面浓度Y（ET）变化图轴面温度T随轴长L的变化（横坐标为L，纵坐标T）轴面环氧乙烷浓度y（ET）随轴长L的变化（横坐标L 纵坐标y（ET）4.1.5计算结果分析1相同气体入口温度下.随着催化剂使用时间的增加，催化效率下降，但环氧乙烷的产率升高。2.相同催化剂投入使用时间下，随着气体入口温度T的增加，环氧乙烷产率下降。3.有效管式反应器长度为0.5米以内。4.随着轴向长度L的增加，温度升高，环氧乙烷浓度变大4.1.6实验结果分析 1.温度随轴向的变化 由于ET氧化合成EO反应是强放热反

33、应热点以前放热速率大于散热速率，而热点之后散热速率大于放热速率。另外大部分床层上的温度变化很小，在轴向距离进口大约2.5 m 处达到热点温度， 然后温度缓慢下降。模拟所得：温度随轴向长度l的增加而不断加大，并且在过拐点后直线上升。温度对产率的影响 通过对反应过程的化学热力学计算可知,各反应均可看作不可逆反应,因此怎样控制反应温度才能使反应在银催化剂的作用下,主要沿热力学趋势最小的主反应进行则显得相当重要。由化学动力学理论可知: 对于平行反应,如果主反应活化能大于副反应活化能,温度增高有利于选择性的增大;反之,降低温度有利于选择性的增大。因此要解决这个问题首先必须知道各主副反应的活化能。关于银催

34、化剂用不同的催化剂进行过测定,其数据见表5 。表5 在银催化剂上主、副反应的活化能kJ / mol研究人 主反应活化能 主要副反应活化能默里 50.21 62.76库里连柯 63.60 82.84波布柯夫 59. 83 89.54由表5 可知,主反应的活化能大于主要副反应的活化能,因此宜于低温操作。但是,温度低则反应速度低、转化率低、收率低、分离困难,不符合生产条件;而温度升高又存在副反应发生及产物的深度氧化等问题。另外,主、副反应均为强放热反应,因此怎样改进反应器的结构,如何更有效的撤走反应放出的热量,进而使所控制的温度对主反应更有利显然是非常重要的问题。模拟所得：温度升高确实反应速率加快而

35、选择性低，但低温条件下收率很大。3氧气在反应器内轴向很短的距离内完全反应。五结论模拟结果与实验结果存在着很大的偏差，我们有理由认为1. 反应气体混合性质较为复杂，对部分物理学常数无法简化为单组份或双组份混合体系。2. COR1，COR2的回归方程存在一定偏差，影响实验结果。3. 产率过低，实际合成过程需要循环气体操作虽然存在很大偏差，但我们建立模型的思想，简化问题的方式是符合工程设计的。做出的程序是可以准确模拟简化后的过程。而我们存在的一些暂时无法解决的问题会在今后的学习中不断完善。附录1MATLA辅助程序clearclc% t=20;% T0=400;h=0.0001;%XX=;YY=;ZZ

36、=;YEO=,;for t=2:2:20 for T0=300:20:500 l=0:h:8;pb=580.1;y0ET=0.6021;y0EO=0;y0CO2=0.0674;y0O2=0.078;y0H2O=0.002;NT0=40.05;M=17.28;Rg=8.314;mt=8821;dt=3.13*10(-2);dp=8.3*10(-3);Tf=485.5;G=93.66;p0=2;Z=0.9954;e=0.0884;o=2.57;F=156.96;nn=8/h;yCO2=;yET=;T=;yEO=;yEO(1)=y0EO;yCO2(1)=y0CO2;yET(1)=y0ET;T(1)=

37、T0;%for i=1:8/h ll=i*h; Cpm=19.2380+0.*T(i)+0.*T(i)2-11.3094*10(-9)*T(i)3;COR1=0.7916+0.02220*t-0.*t2;COR2=0.4359+0.01365*t-0.*t2;k1=exp(12.39103-43.58572*103/(Rg*T(i);k2=exp(18.59960-77.76320*103/(Rg*T(i);K1=exp(-2.82805+18.32100*103/(Rg*T(i);K2=exp(-1.00740+34.66055*103/(Rg*T(i);A=(1-0.5*yET(i)-0.

38、25*yCO2(i)/(1-0.5*y0ET-0.25*y0CO2);yO2=A*(y0O2+1.25*y0CO2-0.5*y0ET)-1.25*yCO2(i)+0.5*yET(i);yH2O=A*(y0H2O-y0CO2)+yCO2(i);pO2=yO2*p0;if pO20 break;endpCO2=yCO2(i)*p0;pET=yET(i)*p0;r1=k1*pO2*pET/(1+K1*pO2+K2*pO20.5*pCO2);r2=k2*pO20.5*pET/(1+K1*pO2+K2*pO20.5*pCO2);T_star=6.66*10(-3)*T(i);nv=1.155/(T_st

39、ar0.1462)+0.3945/(exp(0.6672*T_star)+2.05/(exp(T_star*2.168);n0=26.69*(M*T(i)0.5/(o2*nv);pr=23.82/Z/T(i);n=(1.0230+0.23364*pr+0.58533*pr2-0.40758*pr3+0.*pr4)-1)/e+n0;um=n*10(-3)*3.6;r0=1.52*10(-5)*T(i)1.33;if pr0.5|pr2.0 rm=1.14*10(-2)*(exp(0.67*pr)-1.069)/(F*Z5)+r0; else rm=2.60*10(-3)*(exp(1.155*p

40、r)-2.016)/(F*Z5)+r0; endendRe=dp*G/um;aw=5.1*(dt/dp)0.8*(dt/ll)0.1*(dp*G/um)0.46*rm/dt;rer=0.38*Re0.72*(dt/dp)0.60*rm;Ktb=1/(1/aw+dt/(8*rer);Hm1=-9.9900*104-25.337*T(i)+3.2772*10(-2)*T(i)2-1.6936*10(-5)*T(i)3+3.417*10(-9)*T(i)4;Hm2=-1.*106+15.926*T(i)-2.923*10(-3)*T(i)2-4.9915*10(-5)*T(i)3+1.0362*10

41、(-8)*T(i)4;f1=-(r2*COR2+(1-0.5*yET(i)*r1*COR1)*(1-0.5*yET(i)-0.25*yCO2(i)*pb*A/(1-0.5*y0ET-0.25*y0CO2)*NT0);f2=(2*r2*COR2+0.5*yCO2(i)*r1*COR1)*(1- 0.5*yET(i)-0.25*yCO2(i)*pb*A/(1-0.5*y0ET-0.25*y0CO2)*NT0);f3=(-Hm1*pb*A*r1*COR1-Hm2*pb*A*r2*COR2-3.*Ktb*mt*dt*(T(i)-Tf)/(G*Cpm);yET(i+1)=yET(i)+f1*h;yCO2

42、(i+1)=yCO2(i)+f2*h;T(i+1)=T(i)+f3*h;yEO(i+1)=A*(y0EO+y0ET+0.5*y0CO2)-0.5*yCO2(i+1)-yET(i+1);endii=t/2;jj=(T0-280)/20;YEO(ii,jj)=yEO(i)l=i*h;TT=T(i);X=0:h:l-h;plot3(X,yEO,T);hold on;grid; endendYEOxlswrite(YEO,YEO);附录2 引用文献化工工艺学 高等教育出版社 黄仲九 房鼎业主编乙烯及其部分衍生物工业基础 化工工业出版社 安钢主编乙烯氧化合成环氧乙烷管式反应器的模拟 2009 年第16 卷第6 期化工生产与技术Chemical Production and Technology 汪程刚 李涛* 应卫勇房鼎业流体与过程热力学第二版 郑丹星 编著流体热物性学-基本理论与计算 童景山 编著 环氧乙烷生产中工艺条件的改进方向 李红权万曰亮专心-专注-