化工系统模拟与优化

化工系统模拟与优化丁二烯和二氧化硫合成丁二烯砜的案例研究徐国想指导教师：钱宇教授华南理工大学化工研究所一九九九年八月

目录序言0．输入信息间歇对连续流程图的输入－输出结构2．1．流程图的输入－输出结构2．2．设计变量和物料衡算的步骤流程图的循环结构3．1．过程循环结构的决策3．2．流程的物料衡算别离系统4．1．别离序列的综合4．2．别离系统的总体结构丁二烯制备丁二烯砜过程的全流程模拟5．1．丁二烯制备丁二烯砜过程的全流程图5．2．丁二烯制备丁二烯砜过程的全流程经济衡算参考文献附录PRO/ⅡKEYWORD输入文件

序言概念设计是依据开发性根底研究的结果、文献的数据、现有类似装置的操作数据和工程经验，按照所开发的新技术工业化规模而做出的预想设计。概念设计的目标就是寻找最正确工艺流程和估算最正确设计条件，用以指导过程研究及提出对开发性根底研究进一步的试验要求，所以，它是根底研究与过程研究的指南，是整个开发研究过程中十分关键的一个步骤。经验指出，如果要考虑所有的可能性，那么对任何过程的流程图一般都可以提出许多方案〔经常是104~109个〕，采用分层次决策的方法，就可以在很大程度上防止了漏掉某些重要方案选择的情况。概念设计的目标就是寻找最正确“方案〞，所以，需要尽快降低要予以考虑的方案数目，而基于数量级分析法的简捷算法常常精确得足以消去90％左右无效益的方案。采用简捷算法和分层次决策方法也可以尽快地把信息反响给正致力于开发该过程的化学家手中，将有助于化学家去收集那些有可能获取最大效益操作范围内的数据，并且终止掉那些在有效益操作范围之外的实验。分层次决策方法就是把设计问题简化为分层次的决策，其步骤如下：1.间歇对连续；2.流程图的输入输出结构；3.流程图的循环结构；4.别离系统的总体结构；5.热交换器网络。

输入信息反响信息a.可逆反响：丁二烯＋SO2丁二烯砜丁二烯＋SO2丁二烯砜HR=－48,000Btu/mol,Keq×10-11)exp[-36,940/T(0R)],k1×1015)exp[-52,200/T(0R)],k1=k-1Keq[mol/(ft3·h)],反响速率符合化学计量关系，并且采用CSTR反响器。b.反响条件：反响温度900F,反响器压力=150psia。c.假设反响的选择性为1：d.反响动力学：反响在产品沸点下有显著的逆反响速率，这时应选择适宜的反响温度、二氧化硫与丁二烯的摩尔比，不仅要求转化率高，还要求生产控制稳定、产品本钱低，一般选择二氧化硫与丁二烯的摩尔比r，丁二烯的转化率为40%e.液相，无催化剂2.丁二烯砜产量：P丁二烯砜=80mol/h3.丁二烯砜的产品纯度：xD4.物化数据：一般需要的信息是分子量、沸点、蒸汽压、热容、汽化热、反响热、液体密度和逸度系数〔或状态方程〕。5.费用数据：SO2＝0.064$/mol，丁二烯＝6.76$/mol，丁二烯砜＝8.50$/mol，反响器的年均建设费是3150[$/(ft3·a)]。

间歇对连续采用间歇操作颇为有利的因素产率：a．当产率低于10×106lb/a时，有时采用间歇；b．如果产率低于1×106lb/a时，一般采用间歇；c．多产品装置。市场的力量：a．季节性生产；b．产品的生存期短。放大问题：a．反响的时间非常长；b．要处理低流速的浆料；c．快速结垢性的物料。×106lb/a〔操作时数按8150h/a计〕，一般采用间歇，但根据反响的特点，决定选择一个连续的过程，操作费用和物流费用以年为基准，操作时数为8150h/a。

流程图的输入－输出结构2.1流程图的输入－输出结构1.净化进料物流情况：不需净化进料物流。2.可逆的副产品情况：无可逆副产物。3.气相循环。4.产品物流的数目。所有组分的去向在表2.1中给出，只有一股产品物流。表2.1反响工段各组分去向组分去向SO2循环C4H6循环C4H6SO2产品初步的流程图见图2.1：CC4H6，SO2C4C4H6C4H6SO2图2.1丁二烯制备丁二烯砜过程图2.1丁二烯制备丁二烯砜过程的输入输出结构2.2.设计变量和物料衡算的步骤可能出现在第2层次上的设计变量：对复杂反响有反响器转化率，反响物摩尔比，反响温度/或压力等；过量反响物的情况，过量反响物是指既不回收，又不是气体循环和放空的反响物。正常情况下，可由设计变量开发总物料平衡的公式，而根本不必考虑任何循环流量。只要总方程式定义缺乏，就必须寻找一个或多个设计变量来完善该问题的定义，而这些设计变量总是对应于“主要的过程优化〞问题。因此，最初的分析总是应该只集中于输入－输出的流量。

开发总物料平衡的步骤：从规定的产率入手由化学计量关系求出副产物流量和原料需用量〔以设计变量表示〕在反响物需要全部回收和循环的场合，计算进料物流内杂质的进口流率和出口流率在反响物既不回收也不循环的场合〔循环和放空或是空气或是水〕，以规定的过剩量表示，计算各反响物的出口流率计算随第四步中的反响物流带入的杂质的进口流率和出口流率。

流程图的循环结构过程循环结构的决策：需采用一套反响器系统只有一个主反响，只需要一个反响器系统。由于反响转化率不可能到达100%，故原料要循环回反响器。采用一股气相循环物流全流程有一个回到反响器去的循环物流。其流程图见图2.3。在反响器入口处无需过量反响物配置气体压缩机一台反响在高压下进行，需要一台进料气压缩机。反响器采用预热方式进料通过改变反响器的停留时间来调节平衡转化率。二氧化硫丁二烯二氧化硫丁二烯反响器二氧化硫丁二烯反响器二氧化硫丁二烯丁二烯砜别离系统图2.3丁二烯制备丁二烯砜过程循环结构3.2流程的物料衡算初期物料衡算是基于所有有价值物流全部回收的假设，衡算结果将为用PRO/II作流程模拟提供初值。严格的物料衡算待确定别离系统的细节后由PRO/II流程模拟结果给出。1.设计变量。将产物分布关联为有限反响物的单程转化率的函数，即设计变量为丁二烯的单程转化率。2.物料衡算。从图2.3看到，所有进入该过程的丁二烯都转化了〔没有丁二烯离开系统〕，这是基于丁二烯完全回收的假设，即忽略掉产品物流中有任何丁二烯的损失。当反响转化率为x、反响选择性为S、丁二烯砜产率为P丁二烯砜时，该过程的丁二烯新鲜进料量FFC4H6必须为(2-1)令进入反响器的丁二烯流率为FC4H6，当转化率是x时，离开反响器的丁二烯量是FC4H6(1-x)。对于完全回收的别离系统，离开反响器的流率将等于循环流率。在反响器前的混合器作一衡算，那么有FFC4H6+FC4H6(1-x)=FC4H6(2-2)于是反响器的进料量为FC4H6＝FFC4H6/x(2-3)

4别离系统别离序列的综合选用探试法选择简单塔塔序。探试法的搜索速度一般比数学规划法要高得多，如果探试法那么使用得当，不需要过多的计算便可迅速获得接近最优的别离序列。Nadgir和Liu提出的有序探试法，就是一种较好的别离序列综合的探试法。Nadgir和Liu把探试法那么分为四大类：别离方法探试〔M探试〕法那么；设计探试〔D探试〕法那么组分探试〔S探试〕法那么；组成探试〔C探试〕法那么。这四类探试法那么中包括如下探试法那么，这些探试法那么必须依次顺序采用，不可颠倒：规那么1〔M1法那么〕，尽量采用直接别离法〔只用能量别离剂的别离方法，如普通精馏〕，防止采用间接别离法〔需要用质量别离剂的别离方法〕。但当关键组分的相对挥发度〔或别离系数〕小于～时，不推荐用普通精馏，该采用质量别离剂，但别离剂应在下一个别离器中别离出来。规那么2〔M2法那么〕，防止温度和压力过于偏离环境条件。如果必须偏离，也宁可向高温或高压方向偏离，而尽量不向低温、低压方向偏离。如不个、得不采用真空蒸馏，可以考虑用液液萃取来代替。如果需要冷冻，可以考虑吸收等替代方案。规那么3〔D1法那么〕，倾向于产生的产品个数最少的别离序列，也就是防止别离那些在同一目标产物组中的组分。规那么4〔S1法那么〕，优先别离具有腐蚀性或毒性的成分。规那么5〔S2法那么〕，难于别离的组分最后别离。特别是关键组分的相对挥发度接近于时，应当在没有非关键组分存在的条件下进行别离。规那么6〔C1法那么〕，在相对挥发度允许的条件下，进料组成最大的组分优先别离。规那么7〔C2法那么〕，在进料组分的组成相差不大，且相对挥发度允许的条件下，倾向于将进料一分为二的别离，即将进料等摩尔地分成塔顶和塔底产品。假设难以判断哪一种别离最接近一分为二，那么可以把易别离系数(CES)值最大的别离点优先别离。易别离系数定义如下：(3-1)其中f为产品摩尔流率比，定义为(3-2)式中D—塔顶出料摩尔流率；B—塔底出料摩尔流率。式〔3-1〕中为(3-3)或(3-4)式〔3-3〕表示轻重关键组分间沸点差，式〔3-4〕中为轻重关键组分相对挥发度。易别离系数越大，表示轻、重关键组分越易被别离，所需别离费用一般也较低。4.2别离系统的总体结构从反响器的操作条件〔900F，150psia〕下出来的液相混合物，温度大约在900F，经过换热器进行间接升温至3400F进入别离系统，再通过分别设计蒸汽回收系统和液体别离系统来实现相的最终别离。丁二烯制备丁二烯砜过程的别离系统总体结构如图4.2.1所示：蒸汽别离蒸汽别离系统循环去反响循环去反响器液体蒸汽进料相别离器反响器系统液体蒸汽进料相别离器反响器系统过程别离系统的总体结构0F，再经过换热器降温至900F进反响器。丁二烯制备丁二烯砜过程别离系统的初步结构如图4.2.2所示：循环去循环去反响器3400F150psia0F15psia903400F150psia0F15psia900F150psia进料闪蒸器2减压阀反响器系统闪蒸器1减压阀反响器系统闪蒸器1产品液体产品液体别离系统图4.2.2丁二烯制备丁二烯砜过程别离系统的初步结构该别离系统可采用各种替代方案：升高混合物的温度或降低它的压力，利用闪蒸器来实现相的别离。升高混合物的温度，采用精馏塔来分割物流。升高混合物的温度，采用稳定塔来分割物流。以上各种替代方案按照费用递增的顺序排列，故希望采用排列靠前的方案。

5丁二烯制备丁二烯砜过程的全流程模拟过程的全流程图在以上分析的根底上，可进行全流程的模拟，利用PRO/II软件进行模拟计算，其流程如图5.1所示。811811M1M121新鲜进料C1F121新鲜进料C1F19910412763V1F2REACTOR10412763V1F2REACTOR135E3E2E1135E3E2E1过程的全流程图

过程的全流程经济衡算以年操作时数8150小时为基准，过程的经济潜力为：EP=产品的价值-原料费用-设备费用1．原料及产品费用SO2＝0.064$/mol，丁二烯＝6.76$/mol，丁二烯砜＝8.50$/mol2．反响器的费用模型采用CSTR反响器，反响器的年均建设费可按下式计算：3150[$/(ft3·a)]3．闪蒸器的费用模型一般闪蒸器的停留时间为15min左右〔0.25h〕，故闪蒸器的体积可按下式计算：式中V是闪蒸器的体积〔ft3〕，Qv是进入闪蒸器物流的体积流率〔ft3/h〕。取闪蒸器的L/D=1.5，那么闪蒸器的直径和长度为：；L=12D在引入投资归还因子1/3后，闪蒸器的年度建设费用可按下式计算：在本设计中M&S取792，Fc=3.67*1.0。换热器的费用模型换热器的热负荷均由PRO/II计算得到，换热所需传热面积由下式计算：式中Qc为换热器的热负荷〔Btu/hft20F〕，Uc为传热系数100Btu/hft20F，△Tm是换热器的对数温差，Ac为换热器所需换热面积〔ft2〕。在引入投资归还因子1/3后，换热器的年度建设费用可按下式计算：在本设计中M&S取792，当采用碳钢制造的浮头式换热器时Fc=1.0。另外，换热器的年度操作费用〔消耗的冷却水费用〕可按下式计算：式中冷却水的价格为0.06$/1000gal，△Tc是冷却水的进出口温差，年操作时数为8150小时。换热器年度费用=换热器年度建设费用+换热器年度操作费用压缩机的费用模型对于初期设计，假定压缩机效率为90%，压缩机理论功率hp由PRO/II流程模拟结果给出，那么制动马力为用Guthrie关联式计算压缩机的建设费用操作费用按压缩机电机效率为90%和电费为35$/kWh计算式中bhp单位为kW。一、费用数据二氧化硫费用$/kmol丁二烯费用$/kmol丁二烯砜价值8．50$/kmol电价35$/kWh冷却水×10-4$/kmol设备费用估算价格指数M&S=792,并引入投资归还因子1/3年，以便把建设费用置于年度的基准上。二、经济潜力估算以反响转化率x为例1．原料二氧化硫费用=(0.064$/mol)(80mol/h)(8150h/a)=104$/a原料丁二烯费用=(6.76$/mol)(80mol/h)(8150h/a)=106$/a产品丁二烯砜价值=(8.50$/mol)(80mol/h)(8150h/a)=106$/a2．反响器费用反响器进料体积流率v0=km3/h反响的速度常数k=1.045s-1(T=900F)×10-4h反响器容积V=0.065m3反响器年度费用＝3150[$/(ft3·a)]×35.315ft3/m3104$/a3．闪蒸器费用××35.315=2.709ft3式中V是闪蒸器的体积〔ft3〕，Qv是进入闪蒸器物流的体积流率〔ft3/h〕。取闪蒸器的L