乙醇+丁醇+辛烷+壬烷+水+十四烷液液相平衡研究

袁慎峰, 包根来, 尹 红, 陈志荣

乙醇+丁醇+辛烷+壬烷+水+十四烷液液相平衡研究

袁慎峰, 包根来, 尹 红, 陈志荣

(浙江省化工高效制造技术重点实验室, 浙江大学 化学工程与生物工程学院, 浙江 杭州 310027)

针对乙醇水溶液萃取丁醇+辛烷+壬烷混合物后的萃取相中残留的辛烷和壬烷，提出了使用十四烷作为反萃剂对萃取相进行反萃，测定了常压、298.15K下乙醇+丁醇+辛烷或壬烷+水+十四烷五元体系和乙醇+丁醇+辛烷+壬烷+水+十四烷六元体系的液液相平衡数据。使用Othmer-Tobias方程检验了液液相平衡数据的准确性，线性相关系数2≥0.986 3。使用非随机两流体(NRTL)模型对五元体系的液液相平衡数据进行回归，得到了相关的二元交互作用参数；利用回归得到的二元交互作用参数对六元体系进行预测，结果表明，预测值与实验值吻合良好，这表明回归得到的二元交互作用参数具有较高的准确度。

丁醇+辛烷+壬烷混合物；萃取；液液相平衡；非随机两流体(NRTL)；Aspen Plus

1 前 言

费托合成反应[1-4](Fischer-Tropsch synthesis)以合成气为原料，在催化剂作用下生成烃类和含氧有机物等高附加值产物。费托合成产物低碳醇中的丁醇是一种重要的化工基础原料，广泛用于石油工业、医药、织物制造以及燃料[5-8]。费托产物中的丁醇与辛烷、壬烷易形成共沸物，普通精馏难以将丁醇分离出来，一般多采用共沸精馏、萃取精馏、加盐精馏、膜分离和液液萃取等方法[9-10]实现醇和烷烃的分离。

刘甲锋[11]选用质量分数为40%的乙醇水溶液作为萃取剂，从丁醇+辛烷+壬烷混合物中萃取丁醇，选用十四烷作为反萃剂，对萃取相中残留的少量辛烷和壬烷进行反萃，测定了常压、298.15K下乙醇+丁醇+辛烷+水和乙醇+丁醇+壬烷+水两组四元体系的液液相平衡数据，并使用非随机两流体 (non-random two liquids, 简称NRTL) 热力学模型对四元体系液液相平衡数据进行回归关联，得到了二元交互作用参数，利用Aspen Plus软件设计了一个从丁醇+辛烷+壬烷混合物中分离丁醇的工艺流程，并利用该回归得到的二元交互作用参数对全流程进行了模拟和优化，最终得到了质量分数为99.83%的丁醇产品，丁醇回收率达到99.45%。但是，前述工作未对与反萃剂十四烷相关的二元交互作用参数进行回归，也未见乙醇+丁醇+辛烷+壬烷+水+十四烷体系的液液相平衡数据的相关报道。

本研究拟测定常压、298.15K下乙醇+丁醇+辛烷+水+十四烷和乙醇+丁醇+壬烷+水+十四烷这两个五元体系以及乙醇+丁醇+辛烷+壬烷+水+十四烷六元体系的液液相平衡数据，通过Othmer-Tobias[12]方程对液液相平衡数据进行一致性检验，使用NRTL热力学模型对上述两个五元体系液液相平衡数据进行回归关联，得到与十四烷相关的二元交互作用参数，利用回归得到的二元交互作用参数对六元体系液液相平衡数据进行预测，并将实验值与预测值进行对比，以验证回归得到的二元交互作用参数的准确性。

2 实验部分

2.1 实验试剂

本文所用到的实验试剂的具体信息如表1所示。

表1 实验试剂

图1 液液平衡釜

2.2 实验装置及实验步骤

本研究使用平衡釜法测定液液相平衡数据，液液平衡釜如图1所示。实验仪器包括电子天平 (JINGHAI FA1004N，精度0.000 1g)、超级恒温水浴槽(CH-1015，精度0.01 K)、水银温度计(精度0.1 K)、气相色谱 (SHIMADZU GC-2014)、卡尔费休自动水分滴定仪 (ANTING ZSD-1)、磁力搅拌器 (78-1)、压力表 (DYM3-01)。

具体实验步骤如下：

(1) 使用电子天平准确称量不同配比所需的各个组分质量，并倒入平衡釜内；

(2) 将平衡釜放置在磁力搅拌器上，开启搅拌，接通超级恒温水浴槽，确保循环水充满平衡釜夹套。待超级恒温水浴槽的温度达到设定值且保持稳定时，持续搅拌2 h；

(3) 缓慢关闭搅拌，轻微震动平衡釜确保釜壁无液滴粘挂，静置8 h；

(4) 使用注射器从上、下取样口分别抽取澄清的烃相和醇相样品，密封保存样品等待分析；

(5) 通过气相色谱仪测定两液相样品的组分浓度，通过自动水分滴定仪测定两液相样品中水的含量。

2.3 分析方法

两相样品中有机物组分含量采用气相色谱、以仲丁醇为内标物进行测定，两相样品中水含量则利用卡尔费休容量法通过自动水分滴定仪进行测定。

气相色谱的分析条件如下：载气为氮气(体积分数99.999%)，色谱柱为HT-1毛细管色谱柱(30 m×0.32 mm×0.5 μm)，检测器为氢火焰离子化检测器(FID)，进样器温度为593.15 K，检测器温度为593.15 K，柱前压为0.1 MPa，初始柱温为333.15 K并保持3.5 min，再以升温速率20 K×min-1升至503.15 K后保持5 min。进样量为1 μL，体积分流比为50:1。

3 结果与讨论

3.1 五元体系液液相平衡数据

常压、298.15 K下乙醇(1)+丁醇(2)+辛烷(3)+水(5)+十四烷(6)和乙醇(1)+丁醇(2)+壬烷(4)+水(5)+十四烷(6) 两个五元体系的液液相平衡数据分别如表2、3所示。表中，1～6分别代表乙醇(1)、丁醇(2)、辛烷(3)、壬烷(4)、水(5)、十四烷(6)的质量分数。

表2 乙醇(1)+丁醇(2)+辛烷(3)+水(5)+十四烷(6)五元体系在常压和298.15K下的液液相平衡数据

aStandard uncertaintiesare() = 1kPa,() = 0.1K and() = 0.000 7.() is standard uncertainty of pressure,() is standard uncertainty of temperature,() is standard uncertainty of mass fraction.

表3 乙醇(1)+丁醇(2) +壬烷(4)+水(5)+十四烷(6)五元体系在常压和298.15K下的液液相平衡数据

aStandard uncertaintiesare() = 1kPa,() = 0.1K and() = 0.000 7.

3.2 六元体系液液相平衡数据

常压、298.15 K下乙醇(1)+丁醇(2)+辛烷(3)+壬烷(4)+水(5)+十四烷(6) 六元体系的液液相平衡数据如表4所示。

表4 乙醇(1)+丁醇(2)+辛烷(3)+壬烷(4)+水(5)+十四烷(6)六元体系在常压和298.15K下的液液相平衡数据

aStandard uncertaintiesare() = 1kPa,() = 0.1K and() = 0.000 7.

3.3 液液相平衡数据的一致性检验

五元体系和六元体系的液液相平衡数据的一致性可以通过下述Othmer-Tobias方程进行检验。

拟合结果如表5所示，线性关系如图2～4所示，表、图中的溶剂比均为质量比。由表5可见，拟合的线性相关系数2≥0.986 3，表明所测定的五元体系和六元体系的液液相平衡数据具有较好的一致性，用来回归二元交互作用参数具有较高的准确性。

表5 液液相平衡数据的Othmer-Tobias方程拟合结果

图2 乙醇(1)+丁醇(2)+辛烷(3)+水(5)+十四烷(6)五元体系的Othmer-Tobias方程拟合曲线

图3 乙醇(1)+丁醇(2)+壬烷(4)+水(5)+十四烷(6) 五元体系的Othmer-Tobias方程拟合曲线

3.4 二元交互作用参数回归

NRTL热力学模型的二元交互作用参数τ表达式见式(2)：

式中：A、B、E和F均为组分与的二元参数；为温度，K。

刘甲锋[11]已经回归得到了除反萃剂十四烷以外的其他组分间的二元交互作用参数(见表6)，本研究仅对与反萃剂十四烷相关的二元交互作用参数进行回归。为简化计算，A采用Aspen Plus自带的数据库中的默认参数，E和F的值设为0，仅对B进行回归，回归使用Britt-Luecke算法[13]，基于最大似然法则进行计算。与十四烷相关的默认参数和回归参数如表7所示。

图4 乙醇(1)+丁醇(2)+辛烷(3)+壬烷(4)+水(5)+十四烷(6)六元体系的Othmer-Tobias方程拟合曲线

表6 默认和前期工作回归得到的NRTL二元交互作用参数[11]

aDefault parameters.bRegressed parameters.

表7 默认和回归得到的NRTL二元交互作用参数

aDefault parameters.bRegressed parameters.

回归关联的准确性通过计算实验值与模拟值之间的均方根误差(RMSE)和平均绝对误差(AAD)来衡量，计算公式见式(3)、(4)。两个五元体系的RMSE和AAD值如表8所示。

图5 乙醇+丁醇+辛烷+水+十四烷五元体系的液液相平衡数据实验值与模拟值对比图

图6 乙醇+丁醇+壬烷+水+十四烷五元体系的液液相平衡数据实验值与模拟值对比图

表8 五元体系的RMSE和AAD

Table 8 RMSE and AAD in mass fraction for the quinary systems

式中：为数据点的个数，上标exp、cal分别代表实验值和模拟值。

从图5和6、表8可见，实验值与模拟值吻合良好，即两者具有很高的一致性。从表8可以看出，五元体系RMSE的最大值为0.0112，但考虑到此水相中乙醇含量较高，因此相对误差很小。

利用回归得到的二元交互作用参数τ对六元体系液液相平衡数据进行预测，实验值与预测值的对比如图7所示，图中上标est代表预测值；预测的准确性通过RMSE和AAD进行衡量，如表9所示。同时，使用Aspen Plus自带数据库中默认的二元交互作用参数对六元体系液液相平衡数据进行预测，并与实验值对比，如图8所示。

图7 乙醇+丁醇+辛烷+壬烷+水+十四烷六元体系的液液相平衡数据实验值与预测值对比图

图8 乙醇+丁醇+辛烷+壬烷+水+十四烷六元体系的液液相平衡数据实验值与Aspen Plus默认参数预测值对比图

表9 六元体系的RMSE和AAD

结合图7和表9可见，利用回归出来的二元交互作用参数得到的六元体系的液液平衡预测值与实验值吻合良好，具有很高的一致性。而从图8可见，利用Aspen Plus自带的数据库中默认的二元交互作用参数计算得到的液液平衡数据与实验值偏差较大。因此，本研究通过两个五元体系的液液相平衡实验数据回归得到的二元交互作用参数，可以更加准确地预测六元体系的液液相平衡，即本研究回归得到的二元交互作用参数可用于丁醇+辛烷+壬烷混合物分离工艺的模拟与优化。

4 结 论

本研究测定了常压、298.15K下乙醇+丁醇+辛烷或壬烷+水+十四烷体系的液液相平衡数据，使用NRTL热力学模型对该五元液液相平衡数据进行回归，得到了与反萃剂十四烷相关的二元交互作用参数，该参数可以更加准确地预测乙醇+丁醇+辛烷+壬烷+水+十四烷六元体系的液液相平衡，提高丁醇+辛烷+壬烷混合物分离工艺模拟结果的可靠性。

[1] ANDERSON R B. Fischer-Tropsch synthesis [M]. New York: Academic Press Inc, 1984.

[2] 孙启文, 吴建民, 张宗森. 费托合成技术及其研究进展 [J]. 煤炭加工与综合利用, 2020(2): 35-42, 4.

SUN Q W, WU J M, ZHANG Z S. Fischer-Tropsch synthesis technology and research progress [J]. Coal Processing & Comprehensive Utilization, 2020(2): 35-42, 4.

[3] 白天忠, 安良成, 张安贵. 费托合成油产品开发研究 [J]. 广东化工, 2020, 47(6): 141-142.

BAI T Z, AN L C, ZHANG A G. Research and development of Fischer-Tropsch synthetic oil [J]. Guangdong Chemical Industry, 2020, 47(6): 141-142.

[4] 薛丽娜, 刘岩. 费托合成工艺发展现状 [J]. 山西化工, 2014, 34(3): 36-38.

XUE L N, LIU Y. Development status of Fischer-Tropsch synthesis process [J]. Shanxi Chemical Industry, 2014, 34 (3): 36-38.

[5] 崔兆然. 正丁醇的合成进展简述 [J]. 中国石油和化工标准与质量, 2019, 39(7): 136-137.

CUI Z R. Brief introduction to the progress of synthesis of n-butanol [J]. China Petroleum and Chemical Standard and Quality, 2019, 39(7): 136-137.

[6] 刘博. 正丁醇添加于柴油-生物柴油对柴油机性能及颗粒物排放的影响研究 [J]. 河南科技, 2018(20): 129-132.

LIU B. Study on n-butanol addition to diesel-biodiesel effecting on engine performance and particulates emissions [J]. Henan Science and Technology, 2018(20): 129-132.

[7] 姜莉, 朱荣富, 石美玉, 等. 正丁醇作为生物燃料在生产和应用方面的进展 [J]. 中外能源, 2014, 19(7): 30-35.

JIANG L, ZHU R F, SHI M Y,. Progress in production and application of n-butanol as a biofuel [J]. Sino-Global Energy, 2014, 19(7): 30-35.

[8] 刘娅, 刘宏娟, 张建安, 等. 新型生物燃料—丁醇的研究进展 [J]. 现代化工, 2008, 28(6): 28-31, 33.

LIU Y, LIU H J, ZHANG J A,. Research progress in new biofuel butanol [J]. Modern Chemical Industry, 2008, 28(6): 28-31, 33.

[9] 张琪, 王涛, 张雪冰, 等. 费托合成产物中含氧组分资源化利用技术进展 [J]. 化工进展, 2020, 39(8): 1-15.

ZHANG Q, WANG T, ZHANG X B,. Advances in resource utilization of oxygenated compounds in Fischer-Tropsch synthesis products [J]. Chemical Industry and Engineering Progress, 2020, 39(8): 1-15.

[10] 马爱华, 云志. 费托合成水相副产物中具有共沸组成的低碳混合醇-水体系分离方法的研究进展 [J]. 石油学报(石油加工), 2013, 29(4): 738-743.

MA A H, YUN Z. Research progress of separation technologies for azeotropic mixture of lower alcohols-water system of the by-product in water of Fischer-Tropsch synthesis [J]. Acta Petrolei Sinica (Petroleum Processing Section), 2013, 29(4): 738-743.

[11] 刘甲锋. 从低碳醇烃混合物中分离丁醇的实验与模拟研究 [D]. 杭州: 浙江大学, 2017.

LIU J F. Study on the separation process of butanol from mixtures containing low-carbon primary alcohols and n-alkanes: Experiments and simulation [D]. Hangzhou: Zhejiang University, 2017.

[12] OTHMER D, TOBIAS P. Liquid-liquid extraction data – the line correlation [J]. Industrial & Engineering Chemistry Research, 1942, 34(6): 693-696.

[13] BRITT H I, LUECKE R H. The estimation of parameters in nonlinear, implicit models [J]. Technometrics, 1973, 15(2): 233-247.

Liquid-liquid equilibrium of ethanol + butanol + octane + nonane + water + tetradecane

YUAN Shen-feng, BAO Gen-lai, YIN Hong, CHEN Zhi-rong

(Zhejiang Provincial Key Laboratory of Advanced Chemical Engineering Manufacture Technology, College of Chemical and Biological Engineering, Zhejiang University, Hangzhou 310027, China)

Tetradecane was used as back-extractant to extract the residual octane and nonane in the extraction phase after the extraction of butanol + octane + nonane mixture with ethanol aqueous solution. The liquid-liquid equilibrium (LLE) data of the quinary systems ethanol+butanol+octane or nonane+water+tetradecane and the six-component system ethanol+butanol+octane+nonane+water+tetradecane were obtained under 298.15K and atmospheric pressure. The accuracy of the LLE data was examined by the Othmer-Tobias equation with the linear correlation2≥0.986 3. The LLE data of the quinary systems were regressed by the non-random two liquids (NRTL) model and related binary interaction parameters were obtained to predict the LLE data of the six-component system. The results show that the predicted values are in great agreement with experimental data, which indicates that the regressed parameters have high accuracy.

butanol+octane+nonane mixture; extraction; liquid-liquid equilibrium; non-random two liquids (NRTL); Aspen Plus

TQ413.2

A

10.3969/j.issn.1003-9015.2021.03.002

1003-9015(2021)03-0400-07

2020-04-27；

2020-06-19。

袁慎峰(1977-)，男，山东聊城人，浙江大学副教授，博士。

尹红，E-mail：yinh@zju.edu.cn