汽油馏分蒸气压的计算

汽油馏分蒸气压的计算［摘要］推荐汽油组分蒸气压计算方法。技术人员一直追求更为简介准确的计算方法计算各种物性参数，譬如蒸气压，虽然已经有许多方法可以计算，但是有的简洁但结果误差相对较大，有的准确但计算复杂，本文推荐立方型状态方程计算气液相平衡，求取汽油馏分在任意指定温度下的饱和蒸气压，为汽油馏分的气液相平衡计算提供支持。方程求解通过编程实现，输入参数包括化合物的临界温度、临界压力和偏心因子，输出为对应温度下化合物的蒸气压。利用该方法可以计算汽油馏分中所包含的400多种烃化合物的蒸气压，比较了在378C下的蒸气压计算值和实测值，结果吻合较好，轻烃化合物蒸气压的平均相对误差低于1%，汽油全馏分范围内烃化合物蒸气压的平均相对误差低于4%。可以用于实际设计工作和生产核算工作中。蒸气压是化合物重要的热力学性质之一，指在一定温度下该物质的气液两相处于平衡状态时的气相压力，用于衡量流体的蒸发性能。化合物的蒸发趋势越大，蒸气压越大；蒸发性大的物质临界压力大，临界温度低，密度小，沸点低，称为轻组分，反之为重组分。在炼油和石油化工领域的安全、环保、优化设计、HAZOP分析，以及安全平稳操作等相关计算中都需要准确的蒸气压值。如：烃蒸发损失以及烃蒸气在空气中的爆炸极限等均与烃分子的蒸气压有关。挥发性高的戊烷掺入汽油中可以提高汽车发动机的点火性能，而挥发性低的化合物可以减少蒸发损失和降低气阻形成的倾向，汽油调和是轻、重组分含量的优化结果。纯烃化合物在不同温度下的蒸气压可以通过实验测定，在TRC热力学性质表中收集了部分实测数据，在API-TDB中也用图线方式表示出一些烃类和非烃类化合物蒸气压与温度的对应关系。文献中纯烃化合物蒸气压的计算方法大体分两类：一类是蒸气压与温度的直接经验关联式，另一类是普遍化关联式。如Antoine方程、Daubert方程这类直接经验关联式中一般包含3-5个关联常数，而且关联常数与烃分子一一对应，计算软件采用该类关联式计算化合物蒸气压；通用关联式是更实用的预测方法，即只需要化合物的基础性质如沸点等，简单形式如van-van关联式，复杂形式如Maxwell-Bonnel关联式，如果纯物质的临界性质如临界温度、临界压力和偏心因子3可以得到，也可以采用更详细的关联模型，如Lee-kesler关联式等。上述关联式都可以用来计算烃化合物的蒸气压，但是在实际应用中，一方面可得到的资源不一定能够满足对公式中关联常数的需求，另一方面多数关联式的适用温度范围受限，从而影响计算精度。推荐一种求解立方型状态方程的方法来计算烃化合物的蒸气压，减少对关联参数的依赖，从而满足更广泛的需求。汽油组分相关蒸气压计算方法汽油馏分分子组成模型是构建分子水平过程动力学模型的基础步骤，在过程模型中对原料和产物进行物性表征必不可少，因此为汽油馏分中的烃化合物提供蒸气压值是过程动力学模型的基础组成部分。依据石脑油、催化裂化汽油等样品的详细烃分析数据以及己经建立的石油烃类化合物分类方法，结合各种化合物的反应性能，为汽油馏分建立了分子组成模型。石脑油和重整油分子组成模型中包含200余种单体烃化合物，催化裂化汽油组成模型中包含300余种单体烃化合物，两类模型中烃化合物总数达400余种。如果利用前已述及的关联式计算这些烃化合物的蒸气压值，比如常用的Antoine方程：比卢=A_草ICinT\DT2[ (1)式中：A-E均为与烃分子自身 对应的关联常数；P、T分别为对应的温度和压力。其它关联式如Daubert方程：Inpv=( )(衣+方丁'3 '+4匚‘)(2)式中：a〜d为组分关联常数；pr、Tr分别为对比压力和比温度，T=1-Tr。上述关联常数在文献中仅能提供部分化合物的数据，因此数量上未能满足需要；如果采用如下的Lee-Kesler普通化关联式：In队=S92714—6罕扑一L288621n7；+

上rQI69347f"1&旺18-”一li472lln7；+Q43577V；) (3)式中，3为偏心因子，尽管使用方便，但是式中的关联系数或许会限制关联式的使用温度范围。以下将解析如何采用求解状态方程的方法预测汽油馏分中烃化合物的蒸气压。烃化合物蒸气压的计算模型a)立方型状态方程

作为一种描述流体PVT行为的模型，立方型状态方程既可以同时计算液相和气相的性质，又可以包含较宽的温度和压力范围。在工程界相平衡计算中，SRK方程和PR方程被广泛接受并经常使用。SRK状态方程基本形式及计算式如下：RT_ .P~V-h~V(V|b)a=a,a0.4274E砂史c= Pc九=Q48+L57U-Q176「 0.086(vlRTch= PR状态方程基本形式及计算式如下：'V~hV(V\b}-b(V~b)a= »Cl4572481^-T;E= 可 十「I•妇】「1/：„=Q37464+154226由一。26992; 0.0778RTcb= (4)(5)(6)(7)(8)(9)(10)(11)(12)O?(13)(U)两方程中都包含了参数a、b，(4)(5)(6)(7)(8)(9)(10)(11)(12)O?(13)(U)b)烃分子蒸气压的求解立方型状态方程是关于摩尔体积的三次方，故解方程可以直接求得气、液两相的摩尔体积。为方便求解，将V、a、b分别无因次化为Z、A、B，如下所示:(15)A——(15)当气、液两相处于平衡状态时，气相、液相的逸度f相等，由于两相压力也相等，因此又可以认为逸度系数6相等。利用上述无因次变量可以表示逸度系数，两方程中6的计算式分别如下:A7In皆=Z—1— —13—In (16)JlJ•ZL*jLJ1i y[ 1 /rji —x/2)me=Z一1—mrZ— In —(17)2-J2B7+B(l+也)式(16)用于SRK方程，式(17)用于PR方程。以SRK方程为例，说明由纯烃化合物的临界性质Tc、Pc、3计算在任意温度T下蒸气压p的计算过程。步骤如下：①输入已知量：Tc、Pc、3和T，利用式(5)〜式(9)计算状态方程参数a、b；②计算无因次参数A、B：据蒸气压p的假设值，利用式(15)计算A、B；③计算压缩因子：假定摩尔体积的初值，利用迭代法求解气、液相的摩尔体积V，并计算Z值；④计算逸度系数：利用式(16)分别计算气、液相的逸度系数6v和6l；⑤判断蒸气压的假设值是否合理：定义一个计算终止的误差判断式8=I61—6vI，当误差值小于一个设定的小值，则认为计算可以终止，此时的P值为最终的计算值；如果不满足该条件，则调整P值，重复步骤②〜⑤，继续迭代，直至满足要求为止。PR方程的求解过程与SRK方程类似。结果与讨论针对第1节的方程，通过查阅数据库和计算方法得到全部化合物的临界性质Tc、Pc、3，将其作为蒸气压计算程序的输入参数。求解立方型状态方程，计算400多种烃化合物在任意温度下的蒸气压值，状态方程可以在SRK方程和PR方程之间切换选择。以下通过对计算值与实测值的对比来考察各种计算方法的精度和可靠性。a)状态方程法与普遍化关联法的比较部分烃化合物在37.8°C下蒸气压的计算值与实验值的对比结果如图1〜图5所示。计算方法包括SRK方程、PR方程以及前述Lee-Kesler通用关联式。化合物的碳数范围是C3〜C10。化合物种类覆盖了汽油组分所含全部化合物类型，包括正、异构烷烃，正、异构烯烃，单取代基和多取代基的环烷烃和芳烃。蒸气压数值范围为0.5〜500kPa。由于在整个汽油馏程内所涉及的蒸气压范围较宽，所以如图1~5分段显示。

图1G〜G怪蒸气压计算值与实测值比较♦ ■—[>]<；•SRK.图2〜图5同2flIN1614IZ1085876S4328 10 12 14 16IB20蒸气压（实测值）/kP^图3G炷蒸气压计算值与实测值比较图4C*媛蒸气压计算值与实测值比较蒸气压（实测值）JkPa3.00 (L5 1.0 1.5 2,0 2.5 3.0蒸气压(实测值)图5G〜C”炷蒸气压计算值与实测值比较图1表示了C4〜C5烯烃组分的蒸气压，数值在100〜500kPa之间，3种方法的预测精度均较好，平均相对误差约为0.5%左右。图2表示了C5〜C6烃化合物的蒸气压，数值在20〜80kPa之间，3种方法的平均相对误差在1%左右，其中SRK方法比另外两者的预测精度更好，误差约为0.7%。图3表示了C7烃化合物的蒸气压，数值在7〜20kPa之间，3种方法的平均相对误差在2.5%左右，其中SRK方程的预测精度较好，误差约为2%，Lee-Kesler关联式的误差较大，约为3.7%。图4表示了C8烷烃化合物的蒸气压，数值在2〜6kPa之间，其中SRK方程的预测精度最好，误差约为2.5%，PR方程预测精度偏低，误差达6.7%。图5表示了C8〜C10烷烃和芳烃类化合物的蒸气压，数值较小，在0.3〜2kPa之间，相比之下PR方程关联式的误差偏大，约为7%左右。由上述分析结果可见，用3种方法预测低碳轻烃的蒸气压值均较可靠；在汽油全部馏程范围内，SRK方程预测效果最好，平均误差低于4%，Lee-Kesler关联式次之，误差低于6%，PR方程平均预测精度低于前两种方法，误差达7%。尽管PR方程预测液体摩尔体积的准确度较高，但是预测蒸气压未见优势。形式简单的Lee-Kesler式源自对三参数对应状态方程的关联，是计算气液相的p-V-T关系较准确的方法，也是API-TDB推荐的方法。b)状态方程法与关联常数法的比较分别采用SRK和PR两种状态方程以及前述Daibert方程、Antoine方程计算苯在127°C下的蒸气压。苯的关联常数可以从文献中查阅，式(3)的对比性质可以计算，将用4种方法得到的计算值与API的实测值进行比较，结果如表1所示。式(1)和Antoine方程2种关联常数法的计算误差分别为2.3%和2.1%，计算精度较好。关联常数通常是在指定的温度范围内回归得到，最理想的蒸气压关联式应该使该物质从三相点到临界点都有效，不过通常还没有关联式在这么宽的范围内都有效。通常在关联温度范围内关联式的计算精度较好，范围之外计算精度变差，而且不同来源的关联常数值可能也有差异，也是回归参数时的温度范围不同所致，通常温度越低，精度越高。这种方法简单方便，适合手工计算。表1中SRK和PR两种状态方程法预测值误差分别为2%和1.4%，都很可靠，预测精度也很接近。由于可查到的蒸气压的实验值数量有限，不便进行对各计算方法适用温度范围的考察。表1苯在127C蒸气压的几神预测值对比项目p.-kPa误基，切API实测值34&0方程预测值PR方程33Q01.3SRK方程3两02.0DauberL方程:瘀02.3Antoine方程35Z3匕14.结论应用SRK和PR两种状态方程计算气液相平衡，建立了纯烃化合物在任意指定温度下蒸气压的计算方法。输入数据为化合物的3个临界参数Pc、Tc和3,输出数据为对应温度下的蒸气压。求取了汽油馏分中所包含的近400种纯烃化合物的饱和蒸气压值。对37.8°C下蒸气压的计算值与实测值的比较结果表明：对于低碳烃化合物来说，这两种方程的蒸气压预测精度接近，平均相对误差均低于1%;