环丁砜芳烃抽提的流程模拟

环丁砜芳烃抽提的流程模拟*王强温晓明费维扬（清华大学化学工程系萃取实验室北京100084）摘要本文对环丁砜芳烃抽提的流程进行了全面深入的分析，建立了相应的单元操作模型和结构模型，并采用序贯模块法进行求解，开发了环丁砜芳烃抽提专用流程模拟软件。结合某厂环丁砜芳烃抽提装置的技术改造，利用本文所开发的软件对装置改造后的标定数据进行了核算，计算结果与操作数据吻合较好。关键词：芳烃抽提，化工流程模拟一、前言芳烃抽提是重要的石油化工过程，它采用萃取的方法分离加氢汽油、重整油等含芳原料中的芳烃和非芳烃。环丁砜作为芳烃抽提的溶剂具有溶解能力大、选择性高、稳定性好和易于回收等诸多优点，因此该工艺自六十年代初工业化以来得到了迅速的推广应用，我国也已先后引进了数套环丁砜芳烃抽提装置，取得了良好的经济效益[1]。但是由于存在着以下困难，环丁砜芳烃抽提的流程模拟一直未能很好地实现：1.环丁砜芳烃抽提体系复杂，非理想性严重。整个体系含有多达数十种的烷烃、环烷烃和芳烃，且沸点相差很大；又由于环丁砜和水的加入，使得体系具有严重的非理想性，计算比较困难。2.缺乏必要的基础数据，如环丁砜的物性、传递性质以及相关的热力学参数等。3.流程结构复杂。由于物料和能量的综合利用，使得流程中含有多条再循环回路，各设备之间相互联系、相互影响，增加了模拟的困难。正是由于以上原因，一些通用的化工流程模拟系统无法直接用于该流程的模拟计算，而在引进设备时，国外承包商也未能提供相应的数据和计算方法[1]。为了完成对引进设备的消化吸收，进行环丁砜芳烃抽提装置的优化操作和设计，有必要开发一套专用的环丁砜芳烃抽提流程模拟系统。二、流程概述*本文得到国家自然科学基金的资助和国家重点化学工程联合实验室的支持。环丁砜芳烃抽提装置主要由抽提塔、汽提塔、溶剂回收塔、水汽提塔、溶剂再生塔、抽余油水洗塔、抽余油分馏塔等分离设备，以及有关的热交换设备和流体输送设备组成。整套流程中含有近30个单元操作设备和近50条流线，包括烃循环、溶剂循环和水循环等三类循环回路，是一个复杂的分离过程，并且由于原料和分离要求的不同，各厂的实际流程还不尽相同。图1为某厂环丁砜芳烃抽提装置的流程示意图。图1环丁砜芳烃抽提流程示意图三、模型与算法清华大学萃取分离实验室在环丁砜芳烃抽提过程的热力学与动力学、抽提设备的性能和物性数据测定等方面进行了系统的实验研究[2-7]，并开展了有关分离塔设备计算机辅助设计的研究[8]。在此基础上，本文对环丁砜芳烃抽提流程进行了全面深入的分析，开展了环丁砜芳烃抽提专用流程模拟系统的研究开发工作。1.热力学模型考虑到环丁砜芳烃抽提装置在中低压下操作，选择两项VIRIAL方程描述气相的非理想性，并用Hayden和O’Connell的一般化方法计算第二VIRIAL系数[9]。液相非理想性的描述采用UNIFAC基团贡献法，由于流程中同时包含有精馏和萃取操作，在进行汽液平衡和液液平衡的计算时分别采用相应的UNIFAC基团交互作用参数。2.单元操作模型对环丁砜芳烃抽提流程中所包含的各类单元操作进行了详细的分析，分别建立了流股分割器、流股混合器、水冷（热）器、换热器、液液等温分相器、精馏塔和萃取塔等七类单元操作模型。经过分析，发现环丁砜芳烃抽提流程模拟结果的准确与否主要取决于其中塔设备的模拟计算情况，因此对流程中的各类精馏塔和萃取塔进行了详细的分析，建立了统一的多组份逆流分离的复杂塔模型，如图2所示。图2通用模型塔该模型塔描述如下：(1)塔具有N块板，编号从塔顶至塔底依次为1,2,…,N。对于有冷凝器和再沸器的情况，二者分别为第1和第N块板。塔处理的体系包含M个组份。(2)对第j块板，有一股进料Fj，其组成为Zfij；一股汽相（轻相）采出Gj；一股液相（重相）采出Sj。该板的加热量为Qj。规定G1=0，SN=0。(3)离开第j板的汽相（轻相）流量为Vj、组成为yij；离开第j板的液相（重相）流量为Lj、组成为xij。上面叙述中，j=1,2,…,N；i=1,2,…,M。12jj-1j+1N-1NjLj-1xij-1,Hlj-1Vjyij,HvjGjSjLjxij,HljVj+1yij+1,Hvj+1QjFjZfij,HfjV1F1S1FjGjSjFNGNLN根据不同的操作情况将塔设备的工艺计算分为四类：使用分凝器的精馏塔（A类）、使用全凝器的精馏塔（B类）、提馏塔（C类）和萃取塔（D类），如表1所示。对有冷凝器的精馏塔定义回流比R如下：RLVS=+111表1塔设备工艺计算的分类塔类型ABCD已进料情况各板进料流量Fj及进料组成Zfij(i=1,…,M;j=1,…,N)知侧线采出各板侧线采出流量Gj、Sj(j=1,…,N)的压力各板压力Pj(j=1,…,N)条加热情况Qj(j=2,…,N-1)Qj(j=1,…,N-1)Qj(j=1,…,N)件设计条件R,V1（或LN）R,V1=0V1（或LN）求流量分布Vj、Lj(j=1,…,N)解组成分布yij、xij(i=1,…,M;j=1,…,N)变温度分布Tj(j=1,…,N)量附加变量Q1、QNQN以平衡级模型为基础，建立了上述复杂塔的MESH方程组，并在模型中引入了组份的Murphree效率，使得程序对理论板和实际板均可进行模拟计算，同时可考虑不同组份在不同的板上具有不同的板效率，突破了平衡级模型的局限。在复杂塔模型的求解上采用联立方程同时求解的Newton-Raphson方法，避免了费时的泡、露点迭代计算，同时合理地给定初值，并结合一些计算技巧，既加快了塔计算的收敛速度，又保证了计算的稳定性。3.结构模型为了明确地表述系统的结构，采用信息流图和流线联接矩阵建立了系统的结构模型。软件采用序贯模块法进行环丁砜芳烃抽提流程的模拟计算。为了适应于不同的流程结构，采用可及向量法[10]进行系统分解，使得软件可以根据具体的流程结构自动确定切割流线集，同时排出单元模块的求解顺序。切割流线的迭代计算采用直接迭代法，并采用三个偏差指标来控制迭代的进行：各切割流线的各组份分率与上次迭代值的偏差；各切割流线的流量与上次迭代值的偏差；以及各切割流线的温度与上次迭代值的偏差。当上述三项指标同时达到收敛要求时，流程的迭代结束。为防止在迭代中出现振荡甚至发散，在更新切割流线的变量值时加入了阻尼因子，增加了收敛的稳定性。四、软件的结构与特点软件的设计以结构化程序设计思想为指导，采用主程序－子程序的结构。整个软件共分为9个模块：主体控制模块、热力学函数计算模块、数学计算模块、单元计算模块、单元接口模块、数据输入模块、系统分解模块、流程收敛控制模块和结果输出模块。主体控制模块用来安排程序运行时的数据结构，有序地调用软件的其它模块进行外界信息的输入、流程的分解、迭代计算和收敛判断，并输出模拟的结果，该模块的框图如图3所示。图3主体控制模块框图其它的模块均为主体控制模块所调用，分属于三个不同的层次：最底层为热力学函数计算模块和数学计算模块；单元水平上的为单元计算模块和单元接口模块；流程水平上的为数据输入模块、系统分解模块、流程收敛控制模块和结果输出模块。各模块分述如下：热力学函数计算模块用于计算各个组份的热力学函数（如逸度系数、活度系数、汽液平衡和液液平衡的分配系数等）、流股的焓值和泡、露点，以及有关参数对温度和组成的偏导数。偏导数的计算尽量使用解析式，对不易得到解析式的偏导数计算用差分代替。数学计算模块用于提供在流程模拟中所需用到的一些数学方法，包括矩阵相乘、高斯列主元消去法解线性方程组、三对角方程和块三对角矩阵方程求解以及流线中各物流变量的单位转换等。单元计算模块用于完成相应单元操作的物料衡算和热量衡算。单元接口模块用于完成数据在主体控制模块和各单元模块之间的传递。经过程序初始化，流程迭代次数置0输入物性数据和系统结构数据确定每类单元模块的个数，输入单元数据分析流程结构，确定单元的输入输出流线进行系统分解，确定单元计算顺序和切割流线集读入外界输入流股数据输入切割流线的变量初值迭代次数加1ForI=1To单元总数通过单元接口模块，调用相应单元计算模块求解第I单元调用流程收敛控制模块计算迭代偏差，更新切割流线的变量值Until迭代偏差≤偏差限调用结果输出模块进行结果输出结束仔细设计，使得各单元计算模块尽可能地独立于主体控制程序，有利于模块的修改、扩充、替换和增减。数据输入模块用于进行数据的输入。数据按其应用分为四个部分：纯组份及基团的性质数据；反映系统中各单元设备类型及确定流线联系矩阵的系统结构数据；单元设备的参数及设计变量数据；系统的外界输入流股数据和切割流股的变量初值。模块从数据文件读入数据，并将其转换为程序运行所需的形式。系统分解模块用于建立和求解系统的结构模型，确定切割流线集，并排出各单元模块的求解顺序。流程收敛控制模块用来控制切割流线的迭代计算。结果输出模块用于将流程模拟的结果转化成工程上所需要的形式，同时计算流股物性，以及塔设备和换热器内流体的物性，并以一定的格式输出。主要输出内容包括：流股的温度、压力、相态、泡露点、摩尔流量、质量流量、体积流量、摩尔分率、质量分率、平均分子量、密度、粘度、热容、摩尔焓、导热系数以及液体的表面张力、气体的压缩因子、绝热指数等；塔内各板的温度、压力、表面张力（界面张力），两相的摩尔流量、质量流量、体积流量、摩尔分率、质量分率、密度、粘度等；换热器的热负荷以及有相变情况下的流体焓与温度的关系等。整套软件具有清晰的结构，同时各模块又具有很强的通用性，使得软件便于维护、扩充和复用。在编码阶段采用标准C语言书写源程序，使得软件具有良好的可移植性，在微机和工作站上均可运行。在软件设计上对体系所含的组份数目和流程的大小不做限制（只受机器内存容量和基础数据的限制），在数据丰富后可用于其它分离过程的模拟计算。另外，对软件的输入信息进行了合理的安排，使得软件的使用简单、方便。五、软件的应用结合某厂环丁砜芳烃抽提装置的技术改造，利用该套软件对改造后的装置标定数据进行了核算。该套装置以加氢裂解汽油为原料，采用UOP技术的七塔流程（参见图1）。表2为原料加氢汽油的组成。表2加氢汽油组成组份苯甲苯C8芳烃C9以上芳烃非芳烃wt%38.821.7613.670.4125.36表3为标定时的主要操作参数。其中溶剂比为加入到T-101塔的贫溶剂与加氢汽油进料之比；第一反洗比为加入到T-101塔的第一反洗液（来自抽余油分馏塔塔顶）与进料之比；第二反洗比为加入到T-101塔的第二反洗液（来自汽提塔塔顶）与进料之比；T-105回流比为溶剂回收塔的回流量与抽提物（芳烃产品）之比。所有比值均为体积流量之比。表3环丁砜芳烃抽提装置主要操作参数操作条件溶剂比第一反洗比第二反洗比T-105回流比数据4.00.080.80.6在上述原料和操作条件下，对整套装置进行了流程模拟。由于不同的芳烃组份在抽提塔中具有不同的板效率，因此使用非平衡级模型对抽提塔中各组份在不同位置上的Murphree效率进行了计算[11]。流程模拟中对抽提塔使用了非平衡级模型的计算结果，即不同的组份具有不同的板效率。表4给出了检测点处流股的模拟计算结果与标定数据的比较（其中的抽余油是指抽提塔顶的出料流股，下同）。作为对比，对抽提塔内各组份的效率均与苯相同的情况进行了计算，计算结果列于表5中。表4模拟结果与标定数据的比较组分wt%流股抽余油抽提物反洗1反洗2计算标定计算标定计算标定计算标定苯.0763微52.256.80.1650.1232.531.7甲苯0.2950.2929.227.9.02090.0010.08.19C8以上芳烃1.491.4918.615.30.000.002.732.54C6非芳62.359.8327ppm204ppm98.492.551.050.6C7非芳31.934.01.407.383.756.90C8以上非芳3.503.670.000.00.0136微表5抽提塔内各组份效率相同的模拟结果流股苯甲苯C8芳烃C6非芳C7非芳C8非芳抽余油0.07630.1630.52763.032.33.54抽提物52