《PROII流程模拟的应用实例分解.doc》由会员分享，可在线阅读，更多相关《PROII流程模拟的应用实例分解.doc（19页珍藏版）》请在三一文库上搜索。
1、3.7 过程系统模拟的应用实例化工流程模拟软件的种类虽然较多，但针对化工过程模拟的基本思想和方法却相通。Aspen Plus应用非常广泛，下面我们以乙苯生产过程乙苯精馏塔生产实例为例介绍Aspen Plus稳态模拟的基本方法和基本思想。某干气制乙苯装置中，干气中的乙烯与苯催化反应生成乙苯，反应产物中同时包括反应副产品丙苯、二乙苯、多乙苯和没有反应完全的苯。反应物经过苯塔分离脱除苯后进入乙苯塔、丙苯塔获得乙苯、丙苯产品。已知进入乙苯塔物料含苯0.2%（摩尔分数，下同）、乙苯77%、丙苯10%、二乙苯12.8%，流量8679kg/h，压力1.45MPa（绝压，下同），温度281。乙苯塔为93块浮阀
2、塔盘，进料位置为第40块塔板，操作中，塔顶压力0.52MPa，冷凝器后压力0.5MPa，塔顶温度210，塔底温度254，回流量19000kg/h，塔顶产品中乙苯含量99.5%，釜液中乙苯含量1.5%，塔顶为全凝器，塔底为热虹吸式再沸器。生产中对乙苯产品指标有严格限制，要求乙苯产品纯度大于99.6%，同时为较小物耗要求乙苯塔釜液乙苯含量小于1%，试通过模拟分析提出改进方案。使用Aspen Plus进行模拟的基本步骤：（1）选择模板（2）选择运行类型（3）创建一个流程 （4）规定计算的全局信息（5）规定组分（6）选择物性方法（7）输入物流规定（8）输入模型规定（9）运行模拟（10）检查结果（11）
3、灵敏度分析（12）生成报告1. 建立一个新的运行当启动Aspen Plus并建立一个新的模拟时，可以从一个空白模拟着手或者从一个模板着手，见图3-27。模板设定了特定工业通常使用的缺省项包括 测量单位、 所要报告的物流组成信息和性质、 物流报告格式、 对游离水选项的缺省设置、 性质方法、 其它特定的应用缺省。Aspen Plus 内置以下列模板：l 空气分离l 化学工艺l 电解质l 气体加工l 一般工艺l 湿法冶金l 石油l 医药l 冶金l 固体l 特种化工对于每个模板，可以选择米制或英制作为缺省测量单位制，其它单位制也可用。 由于本分离工艺物系为芳烃分离，可以选择General with M
4、etric Units 模板作为模拟模板，见图3-27。图3-272. 选择运行类型当建立一个新运行时，必须在New对话框上的Run Type的列表框中选择运行类型。Flowsheet 运行类型为流程模拟计算的基本类型，见图3-27。表3-1 Aspen Plus运行类型表运行类型说明用途Flowsheet模拟计算的基本类型。对流程模拟包括灵敏度分析和优化，同时集成了性质常数估计、化验数据分析/虚拟组分产生、性质分析。化验数据分析一个独立的化验数据分析/虚拟组分产生运行。当你不想在同一个运行中进行流程模拟而分析化验数据时用数据回归一个独立的数据回归运行，能够包含性质常数估计和性质分析计算。将A
5、spen Plus所需的物性模型参数与所测量的纯组分、气液平衡、液液平衡和其它混合数据拟合，Aspen Plus在流程运行中不能执行数据回归。PROPERTIES PLUS一个PROPERTIES PLUS设置运行制备性质包，以便供Aspen Custom Modeler /ADVENT 使用性质分析一个独立的性质分析运行，能够包含性质常数估计和化验数据分析计算。当你不想在同一个运行中进行流程模拟时，通过生成物性值表格而进行性质分析。性质估计一个独立的性质常数估计运行。当你不想在同一个运行执行流程模拟时估计性质参数本算例的目的是对乙苯分离过程进行稳态模拟，因此选择Flowsheet 运行类型，
6、见图3-27。3创建一个流程Aspen Plus是基于流程图的通用流程模拟软件，因此要定义一个能够反应实际过程的流程图。定义一个流程用如下步骤：（1）在View 菜单下确认PFD 状态已经关闭，否则设置的模块和物流图形不能变成模拟模型的一部分，如图3-28；图3-28（2）选择单元操作模块并将它们放置到流程窗口。Aspen Plus 中过程单元操作模型包括8大类、57小类、349个单元操作模型。表3-2 单元操作模型类型模型说明混合器/分流器MixerFsplitSsplit物流混合物流分流子物流分流分离器Flash2Flash3DecanterSepSep2双出口闪蒸三出口闪蒸液-液倾析器多
7、出口组分分离器双出口组分分离器换热器HeaterHeatXMHeatXHetranAerotran加热器/冷却器双物流换热器多物流换热器与BJAC 管壳式换热器的接口程序与BJAC 空气冷却换热器的接口程序塔DSTWU Distl RadFrac Extract MultiFrac SCFrac PetroFrac Rate-Frac BatchFrac 简捷蒸馏设计简捷蒸馏核算严格蒸馏严格液-液萃取器复杂塔的严格蒸馏石油的简捷蒸馏石油的严格蒸馏连续蒸馏严格的间歇蒸馏反应器RStoic RYield REquil Rgibbs RCSTR RPlug RBatch 化学计量反应器收率反应器平衡
8、反应器平衡反应器连续搅拌罐式反应器活塞流反应器间歇反应器压力变送器Pump Compr Mcompr Pipeline Pipe Valve泵/液压透平压缩机/透平多级压缩机/透平多段管线压降单段管线压降严格阀压降手动操作器Mult Dupl ClChong 物流倍增器物流复制器物流类变送器固体Crystallizer Crusher Screen FabFl Cyclone Vscrub ESP HyCyc CFuge Filter SWash CCD 除去混合产品的结晶器固体粉碎器固体分离器滤布过滤器旋风分离器文丘里洗涤器电解质沉降器水力旋风分离器离心式过滤器旋转真空过滤器单级固体洗涤器固
9、体逆流倾析器用户模型User User2 用户提供的单元操作模型用户提供的单元操作模型本例中的任务是对乙苯塔进行模拟，因此选择严格蒸馏模型，如图3-29图3-29根据现场操作过程，选择FRACT2作为模拟模型，并添加乙苯塔位号T-303，见图3-30。图3-30（3）用物流线画出生产过程的物料流程，同时物流线将各单元模块连接起来。本案例中，进塔物流号为101，塔顶液相采出物流号为102，塔底液相采出物流号为103，如图3-31。图3-314规定计算的全局信息全局规定为整个运行设定缺省值，Aspen Plus 运行在输入详细的流程说明之前，要在Setup 窗体中规定全局信息，见图3-32。虽然在
10、选择模板时，Aspen Plus 已经规定了计算的默认值，但可以在任何时候返回这些表格和改变输入。图3-32表3-3 Setup 包含的内容内容说明用途Specification规定输入全局信息Simulation Options模拟选项规定计算闪蒸收敛系统选项时间和误差限制Stream Class物流类定义物流类和物流性质Substreams子物流定义子物流和属性Units-Sets单位集定义度量单位集Custom Units习惯单位确定用户习惯的单位Report Options报告选项规定报告选项5规定组分将模拟过程中所需的组分输入到组分列表中。本案例中所需要的组分为苯、乙苯、丙苯、二乙苯
11、，见图3-33、图3-34。图3-33图3-346选择物性方法选择接近的物性方法经常是决定模拟结果的精确度的关键步骤，同时也是模拟的核心工作。物性方法是一批方法和模型，Aspen Plus 用它们计算热力学物性和迁移性质。热力学物性包括逸度系数、焓、熵、吉布斯自由能、体积，迁移性质包括粘度、热导率、扩散系数、表面张力。Aspen Plus 包括大量内置的物性方法足以满足大部分的应用，但是同时用户可以建立新的适合模拟需要的物性方法。Aspen Plus为用户物性方法的选择提供了参考。表3-4 油和气产品应用推荐的物性方法储水系统PR-BM、RKS-BM板式分离PR-BM、 RKS-BM通过管线输
12、送油和气PR-BM、RKS-BM表3-5 炼油应用推荐的物性方法低压应用（最多几个大气压）、真空蒸馏塔、常压原油塔BK10、 CHAO-SEA 、GRAYSON中压应用（最多几十个大气压）Coker 主分馏器FCC 主分馏器CHAO-SEA 、GRAYSON、 PENG-ROB、 RK-SOAVE富氢的应用重整炉加氢器GRAYSON 、PENG-ROB 、RK-SOAVE润滑油单元、脱沥青单元 PENG-ROB、 RK-SOAVE表3-6 气体加工应用 推荐的物性方法烃分离脱甲烷塔C3 分离器PR-BM、 RKS-BM 、PENG-ROB、 RK-SOVAE深冷气体加工空气分离PR-BM 、R
13、KS-BM 、PENG-ROB、 RK-SOVAE表3-7 化工应用 推荐的物性方法乙烯装置初级分馏器轻烃串级分离器急冷塔CHAO-SEA、 GRAYSON、PENG-ROB 、RK-SOAVE芳香族环烃BTX 萃取WILSON、 NRTL 、UNIQUAC和它们的变化形式乙苯和苯乙烯装置PENG-ROB、 RK-SOAVE、WILSON、 NRTL 、UNIQUAC和它们的变化根据Aspen Plus为用户物性方法提供参考，考虑本案例为乙苯和苯乙烯装置，选用物性方法为PENG-ROB，见图3-35。图3-357输入物流规定对于一个模拟过程，根据自由度分析，必须规定过程系统进料物流的流率、组成
14、、热力学状态，还可以提供撕裂循环物流的初始估值。本案例中，乙苯塔物料含苯0.2%、乙苯77%、丙苯10%、二乙苯12.8%，流量8679kg/h，压力1.45MPa，温度281。将上述数据输入到物流101的物流信息表中，如图3-36。图3-368输入模型规定对于一个模拟过程，根据自由度分析，必须将模型的操作变量输入到模型规定中。不同的模型规定各不相同。本案例中，乙苯塔为93块浮阀塔盘，进料位置为第40块塔板，操作中，塔顶压力0.52MPa，冷凝器后压力0.5MPa，塔顶温度210，塔底温度254，回流量19000kg/h。塔顶为全凝器，塔底为热虹吸式再沸器。将上述信息输送到精馏塔的输入规定中，
15、如图3-37。注意，在模拟计算中，冷凝器和再沸器均按一块板来计，因此计算总板数为95块，进料位置为41块。图3-37值得注意的是塔板效率的选择。精馏塔模拟计算中，塔板效率常常是根据计算结果和实际分离结果比较选定的，塔板效率的初值通常根据塔板型式和分离物性选定。本案例中初选全塔效率60%，如图3-38和图39所示。 图3-38 图3-399运行模拟当所有相关信息添加完成，可以运行模拟10检查结果（1）检查运行的完成状态要检查有关运行收敛和完成状态的汇总信息，可使用Results Summary 页面。这个表指明了计算是否为正常完成。如果没能正常完成，说明模拟过程的输入需要调整，直至系统收敛。当系
16、统收敛后，可以检查计算结果，如图3-40。图3-40如图3-41，在计算结果表中，分别汇集了组分、物流、模块计算信息，根据实际装置运行情况检查模拟结果，如果模拟结果与实际运行结果偏差较大，不能反应实际过程的原貌，则需查找原因，采用调整物性方法等方法直至模拟结果接近实际运行结果。图3-41实际运行中塔顶产品中乙苯含量99.5%，釜液中乙苯含量1.5%。模拟结果顶产品中乙苯含量99.54%，釜液中乙苯含量1.455%，与实际运行结果基本吻合，同时塔顶、塔底温度与实际也基本吻合，因此可以认为该过程的模拟基本准确。11灵敏度分析灵敏度分析是检验一个过程如何对变化的关键操作变量和设计变量反应的一个工具。
17、可以用它改变一个或多个流程变量，并研究该变化对流程其它变量的影响。它是做工况研究的一个最有用的工具。被变化的流程变量必须是流程的输入参数，在模拟中计算出的变量不能被改变。可以用灵敏度分析来验证一个设计规定的解是否在操作变量的变化范围内，还可以用它做简单的过程优化，可以用灵敏度分析模块生成随进料物流模块输入参数或其它输入变量变化的模拟结果的表或图。灵敏度分析的步骤：（1） 建立一个灵敏度(Sensitivity)分析块；单击“new”，系统自动给出分析块的名字，如图3-42。 图3-42（2） 标识被采集的流程变量；对于本案例，乙苯精馏塔操作中，进料位置、回流比对分离结果有很大影响，因此有必要分
18、析进料位置的变化对塔顶乙苯浓度(EBTOP)、塔底乙苯浓度(EBBOT)的影响，回流比的变化对塔顶、塔底乙苯浓度的影响。单击“new”，输入EBTOP，系统弹出一定义变量对话框，如图3-43。 图3-43（3）标识要操纵的输入变量来生成一个表，如图3-44。图3-44（4）定义你要Aspen Plus 将什么制表，如图3-45；图3-45当所有相关信息添加完成，可以运行模拟。得出以下结果图3-46 选择进料板位置为横坐标，塔顶、塔底乙苯摩尔分数为纵坐标，用工具栏中“plot”做灵敏度分析图如图3-47：图3-47由于图3-37给的是回流量，因此做回流比灵敏度分析时，将回流量19000kg/h改
19、为回流比3，然后按上面同样步骤进行灵敏度分析，结果如图3-48所示。图3-48 通过灵敏度分析，如图3-47、图3-48可以发现适宜的进料位置6064，适宜的回流比3.2以上，因此本装置基本处于较优操作，但仍有进一步调优空间。生产中塔顶产品中乙苯含量99.5%，釜液中乙苯含量1.5%，均不能满足生产要求。通过图3-47可知，如果在保证现场回流比3不变的前提下，通过调整进料位置由40块至61块，可使塔顶产品中乙苯含量由99.5%提高至99.68%，釜液中乙苯含量1.5%降至0.9%，达到生产要求。通过图3-48可知，如果在保证进料位置40块不变的前提下，回流比由现场3增加至3.6（相当于回流量由
20、19000kg/h增加至23000kg/h），可使塔顶产品中乙苯含量由99.5%提高至99.67%，釜液中乙苯含量1.5%降至0.96%，达到生产要求。12生成报告生成的报告如表3-5所示。表3-5 乙苯塔现场工况模拟结果表Stage TemperaturePressurekPaLiquid flowkg/hrVapor flowkg/hr液相组成气相组成BENZEETHYLN-PROO-DIEBENZEETHYLN-PROO-DIE1207.60135001915500.0025970.9954420.0019611.93E-070.0069360.9917920.0012717.76E-0
21、82209.867452019510.13255400.0012380.9962110.002553.19E-070.0025970.9954420.0019611.93E-073209.9603520.219516.4925895.130.0006970.9961150.0031884.92E-070.0015730.9960210.0024052.88E-074210.0204520.419519.3525901.490.0004810.995640.0038797.31E-070.0011650.9959490.0028864.18E-075210.0682520.619520.8125
22、904.350.0003940.9949770.0046281.06E-060.0010020.9955910.0034065.99E-0739214.939527.919524.9325916.640.0003220.8484960.1033230.0478590.0008940.8973560.0747460.02700540215.6728528.119373.2125909.930.0003190.8294640.1069990.0632180.0008940.885430.0778460.0358341216.4958528.323019.2320732.350.0002120.80
23、98140.1083170.0816570.0004520.8842570.0782130.03707842216.5336528.623022.520725.230.0001090.8098580.1083750.0816580.0002320.8844260.0782750.03706857217.0547531.823060.8720763.984.97E-090.8031930.1147540.0820531.06E-080.8785660.0839830.03745158217.1203532.023067.3820766.872.55E-090.8015880.1162650.08
24、21475.43E-090.8771170.0853420.03754159217.1916532.223072.1620773.381.31E-090.7996220.1181110.0822672.79E-090.8753570.0869960.03764794249.6843539.725367.6223184.282.16E-210.0241180.5167070.4591758.33E-210.0344430.5737650.39179295252.9471540229423073.625.80E-220.0142930.4286370.557072.31E-210.0250830.5253490.449568