環(huán)己烷-正丙醇-乙二醇體系液液平衡數(shù)據(jù)的測定與關(guān)聯(lián)

摘 要：為了分離醫(yī)藥及化工領(lǐng)域排放的大量環(huán)己烷和正丙醇的混合廢液，采用Materials Studio的Dmol³模塊分析了環(huán)己烷、正丙醇和待選萃取劑的σ-profile圖和相互作用能，選定乙二醇為萃取劑.分別測定了303.15 K、313.15 K和323.15 K下環(huán)己烷-正丙醇-乙二醇的液液相平衡數(shù)據(jù)并繪制三相圖.采用NRTL熱力學模型對液液相平衡數(shù)據(jù)進行關(guān)聯(lián)，得到二元交互作用參數(shù).得到的均方根偏差（RMSD）小于1.0%，表明NRTL模型適用于三元體系液液相平衡數(shù)據(jù)的關(guān)聯(lián).用GUI-MATLAB軟件對回歸的二元交互作用參數(shù)進行了進一步的可靠性檢驗.證實了數(shù)據(jù)關(guān)聯(lián)回歸的有效性，可為萃取工藝模擬與優(yōu)化提供基礎(chǔ)數(shù)據(jù)支撐.

關(guān)鍵詞：相互作用能； 液液相平衡； 環(huán)己烷； 正丙醇； 乙二醇； NRTL模型

中圖分類號：O642.4+2

文獻標志碼： A

Determination and correlation of liquid-liquid balance data of

cyclohexane-n-propanol-ethylene glycol system

DAI Fang-fang^*， PENG Mei-yuan， CAO Jiang-ting， HAN Ling

（Key Laboratory of Auxiliary Chemistry amp; Technology for Chemical Industry， Ministry of Education， Shaanxi University of Science amp; Technology， Xi′an 710021， China）

Abstract：In order to separate a large amount of mixed waste liquid of cyclohexane and n-propanol discharged from the pharmaceutical and chemical fields，the σ-profile and interaction energy of cyclohexane，n-propanol and the selected extractants were analyzed using the Dmol³ model of Materials Studio，and ethylene glycol was selected as the extractants.The liquid-liquid equilibrium data of cyclohexane-n-propanol-ethylene glycol at 303.15 K，313.15 K and 323.15 K were determined and the three-phase diagram was drawn.The NRTL thermodynamic model was used to correlate the liquid-liquid equilibrium data，and the binary interaction parameters were obtained.The root mean square deviation （RMSD） obtained is less than 1.0%，indicating that the NRTL model is suitable for the correlation of liquid-liquid equilibrium data in ternary systems.The reliability test of the regression binary interaction parameters is carried out by using GUI-MATLAB software.The validity of data correlation regression is proved，which can provide basic data support for extraction process simulation and optimization.

Key words：interaction energy； liquid-liquid equilibrium； cyclohexane； n-propanol alcohol； ethylene glycol； NRTL model

0 引言

環(huán)己烷和正丙醇是重要的有機溶劑和有機合成助劑，廣泛應(yīng)用于化工領(lǐng)域.環(huán)己烷是一種無色液體，主要用于合成環(huán)己酮、環(huán)己醇、尼龍6等.少量用作丁基橡膠、油漆和清漆的溶劑，也用作農(nóng)藥制備過程的萃取劑^[1].正丙醇是合成酯的重要原料，同時也是一種優(yōu)良溶劑，被廣泛應(yīng)用于醫(yī)藥領(lǐng)域^[2]，還可用于制作食品添加劑和香料.醫(yī)藥及化工領(lǐng)域會排放大量環(huán)己烷和正丙醇的混合廢液^[3]，直接排放不僅會對自然生態(tài)系統(tǒng)造成破壞，還會威脅人類生命健康安全，同時造成資源浪費.因此，分離回收廢液中的環(huán)己烷和正丙醇十分重要.

正丙醇和環(huán)己烷的沸點相近，形成共沸物，無法用簡單精餾直接分離.邢慶浩等人^[3]選定[N1，8，8，8][Ac]為萃取劑，萃取精餾分離環(huán)己烷-正丙醇共沸體系.研究表明[N1，8，8，8][Ac]的加入會增大環(huán)己烷-正丙醇之間的相對揮發(fā)度.采用液液萃取的方法可以從廢液中對環(huán)己烷和正丙醇混合物進行分離和回收，降低了對環(huán)境的污染，還獲得了良好的經(jīng)濟效益.由于萃取劑眾多，需要進行初篩，才可以得到分離效果較好的萃取劑.使用類導體屏蔽模型的鏈段活度系數(shù)（COSMO-SAC）方法獲得的電荷密度分布函數(shù)（σ-profile）可確認兩個分子之間是否可形成氫鍵，從而定性地測定萃取劑和共沸物體系的氫鍵能力，判斷萃取效果^[4].COSMO-SAC模型是Lin等^[5-7]基于真實溶劑的類導體屏蔽模型（COSMO-RS）提出的.該模型借助量子力學計算分子的表面屏蔽電荷密度，進而表征分子之間的相互作用.其機理是將分子表面分割為多個片段，通過計算每一個片段的電荷密度及表面積，得到每一個片段的表面平均屏蔽電荷密度，統(tǒng)計這些片斷的電荷密度分布得到（σ-profile）.將量子化學計算和統(tǒng)計力學充分結(jié)合的COSMO模型被證明能夠快速篩選萃取劑.Zhang等^[8]基于COSMO-SAC模型和液液平衡法研究了離子液體分離異丙醚-丙酮體系，為共沸混合物萃取劑的選擇提供了新的理論指導.

本文采用COSMO-SAC模型繪制σ-profile圖，計算分子間相互作用能，篩選出乙二醇為萃取分離環(huán)己烷-正丙醇體系的潛在萃取劑.為了探究乙二醇對環(huán)己烷-正丙醇體系的分離特性，測定了環(huán)己烷-正丙醇-乙二醇體系在303.15 K、313.15 K和323.15 K下的LLE數(shù)據(jù)，且計算了分配系數(shù)（D）和分離因子（S）來分析乙二醇萃取分離環(huán)己烷-正丙醇體系的可行性，為分離回收廢液中的環(huán)己烷和正丙醇提供相平衡實驗數(shù)據(jù).采用NRTL熱力學模型關(guān)聯(lián)LLE數(shù)據(jù)，得到二元交互作用參數(shù).利用GUI-MATLAB對二元交互作用參數(shù)的可靠性進行了驗證.為萃取分離工藝和設(shè)計優(yōu)化提供實驗與理論支撐.

1 實驗部分

1.1 實驗試劑和儀器

實驗試劑和儀器如表1和表2所示，實驗中使用的所有試劑的純度通過氣相色譜法（GC）測定，未經(jīng)過進一步純化.

1.2 實驗過程及分析方法

用于獲得LLE數(shù)據(jù)的設(shè)備已在先前的研究中詳細描述^[9].實驗測定了101.325 kPa下，303.15 K、313.15 K和323.15 K時的環(huán)己烷-正丙醇-乙二醇體系的LLE數(shù)據(jù).于洗凈烘干的100 mL相平衡釜中加入體積比為1∶1的環(huán)己烷和乙二醇混合溶液，將溶液攪拌4小時至充分混合，然后靜置8小時至上下分層.達到平衡后，用氣相色譜法測定上、下兩相的組成.采用校準面積歸一法進行定量分析，無水乙醇作為校準物得到各物質(zhì)的相對校正因子，每一個樣品分析五次，取其平均值作為實驗測量值.通過分次向混合物中加入固定量的正丙醇和改變實驗溫度以獲得一系列液液相平衡數(shù)據(jù)^[10]，組成的標準不確定度按文獻[11]計算.色譜條件：Inert Cap WAX GC柱（30 × 0.25 mm × 0.25 μm）；氫火焰離子檢測器（FID）；進樣口溫度493.15 K；檢測器溫度523.15 K；載氣為高純氮氣（99.999%），空氣流量400 mL/min，氫氣發(fā)生器流量30 mL/min；程序升溫（初始323.15 K，以20 ℃/min的速率將柱溫升高至473.15 K并保持1 min）；進樣量為1 μL.

2 結(jié)果與討論

2.1 分子模擬篩選萃取劑

2.1.1" σ-profile分析

Yang等^[12]以1，2-丙二醇、1，4-丁二醇及二甲基亞砜為萃取劑從環(huán)己烷中萃取異丙醇，證明二醇類物質(zhì)和二甲基亞砜從環(huán)己烷中萃取正丙醇是可行的.本文使用Material Studio（MS）軟件中的Dmol³模塊對環(huán)己烷，正丙醇以及待選萃取劑（1，2-丙二醇，乙二醇和二甲基亞砜）進行量子化學計算.首先在Material Studio中繪制分子結(jié)構(gòu)，基于Dmol³模塊和密度泛函理論對初始分子進行結(jié)構(gòu)優(yōu)化，選擇密度泛函（DFT）中的VWN-BP泛函，基組選擇DNP v4.0.0，實際空間截斷值為5.5 ，物性中勾選Electron Density.結(jié)束結(jié)構(gòu)優(yōu)化后對分子結(jié)構(gòu)進行能量優(yōu)化，將分子表面劃分成若干個具有各自面積和屏蔽電荷密度的片段，生成COSMO文件.將COSMO文件中的三維屏蔽電荷密度分布投影到二維直方圖上，得到分子的屏蔽電荷分布函數(shù)（σ-profile），從而得到分子屏蔽表面上的理想電荷分布^[13].在COSMO-SAC模型（2002）^[14]中，純物質(zhì)的電荷密度分布ρ（σ）由式（1）計算：

ρi（σ）=ni（σ）ni=Ai（σ）Ai （1）

式（1）中：ni（σ）和Ai（σ）分別是電荷密度σ的片段數(shù)量和其片段的表面積，ni=∑σni（σ）是表面片段的總數(shù)，Ai=∑σAi（σ）則是空穴的總面積^[15].

圖1為計算的各組分σ-profile圖.譜圖被兩條直線（σ = ±0.008 5 e/²）分為三個區(qū)域，分別為氫鍵供體區(qū)（σlt;-0.008 5 e/²），氫鍵受體區(qū)（σgt;0.008 5 e/²）和非極性區(qū)（-0.008 5 e/²lt;σlt;0.008 5 e/²），若屏蔽電荷處于氫鍵供體區(qū)并且越遠離-0.008 5 e/²，表明物質(zhì)的氫鍵供體能力越強，相反，若屏蔽電荷處于氫鍵受體區(qū)并且越遠離0.008 5 e/² ，表明物質(zhì)的氫鍵受體能力越強^[16].

環(huán)己烷全部分布在非極性區(qū)，難以與萃取劑形成氫鍵.由于醇羥基的存在，正丙醇分布范圍為（-0.017， 0.016）e/²，在非極性區(qū)、氫鍵供體區(qū)和氫鍵受體區(qū)都有分布，與氫鍵供體區(qū)相比，正丙醇在氫鍵受體區(qū)有較大面積，表明有較強的接受氫鍵的能力.1，2-丙二醇的電荷密度分布范圍為（-0.017， 0.016）e/²，在非極性區(qū)和極性區(qū)都有分布，在氫鍵受體區(qū)峰值較大，但延伸較少.二甲基亞砜的電荷密度分布為（-0.01， 0.018）e/²，主要集中在非極性區(qū)域和氫鍵受體區(qū)，在氫鍵供體區(qū)幾乎沒有峰值.乙二醇的分布區(qū)域為（-0.018， 0.018）e/²，由于其具有雙醇羥基和極高的對稱結(jié)構(gòu)，其電荷密度分布在非極性區(qū)有較高的峰值，但整體面積遠小于環(huán)己烷，在氫鍵受體區(qū)和氫鍵供體區(qū)均有更多的延伸，可以與正丙醇形成較強的相互作用，是較為理想的萃取劑.

2.1.2 相互作用能

萃取分離的本質(zhì)是研究萃取劑與共沸體系之間相互作用的差異^[17].本文利用Material Studio軟件中的Dmol³模塊計算了環(huán)己烷，正丙醇與待選萃取的相互作用能.在對每個分子進行優(yōu)化后，將萃取劑分別與環(huán)己烷和正丙醇進行組合，并進一步優(yōu)化，得到各組合的鍵長.進而計算環(huán)己烷、正丙醇與萃取劑的相互作用能.為了使計算結(jié)果更加準確，使用平衡方法通過BSSE^[18]（基本集合疊加誤差）校正相互作用能，如等式（2）和（3）所示：

△E=EAB-EA-EB+EBSSE（2）

EBSSE=EA-E（A， AB）+EB-E（B， AB）（3）

式（2）、（3）中：EAB表示分子之間的相互作用能，EA和EB表示子系統(tǒng)A和B的能量；E（A， AB）表示A在整個系統(tǒng)AB中的能量；E（B， AB）表示B在整個系統(tǒng)AB中的能量，EBSSE是修正能量.優(yōu)化后能量最低的3D分子結(jié)構(gòu)如圖2所示.各組分間的相互作用能計算結(jié)果如表3所示，萃取劑與正丙醇的相互作用能均大于與環(huán)己烷的相互作用能.萃取劑與環(huán)己烷的相互作用能絕對值順序為1，2-丙二醇gt;二甲基亞砜gt;乙二醇，表明乙二醇與環(huán)己烷之間的相互作用能最小.萃取劑與正丙醇的相互作用能絕對值順序為：乙二醇gt;二甲基亞砜gt;1，2-丙二醇，表明乙二醇與正丙醇之間的相互作用能最強，這與σ-profile分析結(jié)果一致.最終選擇以乙二醇為萃取劑來分離環(huán)己烷-正丙醇體系.

2.2 實驗LLE數(shù)據(jù)

101.325 kPa下，測定環(huán)己烷-正丙醇-乙二醇體系在303.15 K、313.15 K和323.15 K下的LLE數(shù)據(jù)如表4所示.表中環(huán)己烷、正丙醇和乙二醇分別由組分1、2和3表示，表中所有物質(zhì)均以摩爾濃度表示.從表中可以看出，兩相中正丙醇的摩爾濃度隨著正丙醇的加入而逐漸增加，且在萃取相中的摩爾濃度高于萃余相中的摩爾濃度，說明正丙醇在乙二醇中有良好的溶解性.且乙二醇在萃余相中溶解度較低，在313.15 K下為最高（0.014 6），說明萃取時大部分萃取劑能夠回收.三元相圖如圖3～5所示，所有相圖均歸屬為第I類相圖，即萃取劑與溶劑部分互溶^[19].相圖左側(cè)的點表示萃取相的組成，右側(cè)的點表示萃余相的組成.環(huán)己烷-正丙醇-乙二醇體系在三個溫度下兩相區(qū)均較大，說明萃取劑的操作范圍較大，有利于環(huán)己烷與正丙醇的分離.隨著溫度的升高，兩相區(qū)減小，說明溫度升高使乙二醇的萃取能力下降.

分配系數(shù)（D）和分離因子（S）是評價溶劑萃取分離效果的重要指標^[20，21].D和S分別由等式（4）、（5）計算：

D=x2βx2α（4）

S=x2β/x1βx2α/x1α（5）

式（4）、（5）中：x2β、x1β分別為正丙醇和環(huán)己烷在萃取相中的摩爾分數(shù)；x2α，x1α分別為正丙醇和環(huán)己烷在萃余相中的摩爾分數(shù).

由表4可知，三個溫度下的S值均大于1，說明乙二醇表現(xiàn)出良好的萃取能力.各溫度下的S值隨著正丙醇摩爾濃度的增加而降低，萃取劑的萃取效果降低.這是由于正丙醇摩爾濃度的增加，兩相區(qū)面積減小導致的^[22].

圖6為S值與萃取相中正丙醇濃度的關(guān)系.由圖可以看出，乙二醇的萃取能力隨著溫度的升高而降低.D值反映了達到一定分離效果所需的萃取劑用量.分配系數(shù)越大，萃取分離過程所需的萃取劑的量越小^[23].

圖7為D值與萃取相中正丙醇濃度的關(guān)系.由圖可以看出，三個溫度下的D值均遠大于1，說明萃取過程需要的萃取劑用量較少.D值隨溫度的升高而降低，說明正丙醇在萃取劑中的溶解度隨溫度升高而降低，需要增大萃取劑用量才能達到分離目的.萃取劑乙二醇的分離因子和分配系數(shù)均較高，對正丙醇具有優(yōu)異的萃取效果，303.15 K為最優(yōu)萃取溫度.

2.3 二元交互作用參數(shù)的回歸

本文通過Aspen Plus 11軟件，使用NRTL熱力學模型對303.15 K、313.15 K和323.15 K時的環(huán)己烷-正丙醇-乙二醇三元體系LLE實驗數(shù)據(jù)進行回歸關(guān)聯(lián)，得到體系的二元交互作用參數(shù)，為分離過程和分離設(shè)備的設(shè)計與優(yōu)化提供了基礎(chǔ)數(shù)據(jù)^[24].NRTL方程的計算公式如下式^[25，26].

lnγi=∑jτjiGjiXj∑kGkiXk+∑jGijXj∑kGkiXk

（τij-∑kτkjGkjXk∑kGkjXk）

Gij=exp（-αijτij）（6）

Gji=exp（-αjiτji）

τij=aij+bij/T

αij=αji

式（6）中：i， j， k分別表示體系中每個組分的符號；aij， bij表示NRTL方程的每個二元交互作用參數(shù).

以式（7）計算均方根偏差（RMSD）來表示實驗數(shù)據(jù)和計算數(shù)據(jù)之間的差異，評估兩個結(jié)果之間的相關(guān)性，計算公式如下.

RMSD=∑^Mk=1∑²j=1∑³i=1（x^expijk-x^calijk）26M^1/2（7）

采用最小二乘法式（8）對實驗數(shù)據(jù)進行回歸擬合，目標函數(shù)（OF）最小時，得到模型的最優(yōu)二元交互作用參數(shù)^[27，28].

OF=∑^Mk=1∑²j=1∑³i=1（x^expijk-x^calijk）2（8）

式（8）中：M表示聯(lián)結(jié)線的總數(shù)，x^expijk表示實驗值，x^calijk表示NRTL模型的計算值.三元體系的NRTL模型計算數(shù)據(jù)繪制在上述相圖3～5中，與實驗得到的LLE數(shù)據(jù)吻合性較好.方程回歸得到的二元交互作用參數(shù)和RMSD值見表5所示，RMSD數(shù)值小于0.01，表明NRTL模型能很好地關(guān)聯(lián)環(huán)己烷-正丙醇-乙二醇三元體系的液液相平衡數(shù)據(jù).

2.4 數(shù)據(jù)一致性測試

為了進一步驗證回歸數(shù)據(jù)與實驗數(shù)據(jù)的一致性，采用MATLAB軟件編寫的圖形用戶界面方法來評估回歸數(shù)據(jù)^[29-31].圖8為101.325 kPa和303.15 K下通過使用NRTL模型對環(huán)己烷-正丙醇-乙二醇三元體系的拓撲分析.整個空間的G^M（L）/RT表面和實驗LLE聯(lián)絡(luò)線如圖8 （a）所示，其中所有實驗LLE數(shù)據(jù)和模擬數(shù)據(jù)聯(lián)絡(luò)線都在表面下方.

從圖8 （b）可以看出，子系統(tǒng)（環(huán)己烷-乙二醇）表現(xiàn)出部分混溶性，子系統(tǒng)（環(huán)己烷-正丙醇，正丙醇-乙二醇）表現(xiàn)出完全混溶性，證明環(huán)己烷-正丙醇-乙二醇體系的三元相圖屬于I型相圖^[19]，與三相圖結(jié)果一致.圖8 （c）～（f）顯示了實驗值和計算值的聯(lián)絡(luò)線（#1和#7）與G^M（L）/RT曲線部分的交點，聯(lián)絡(luò)線位于G^M（L）/RT平面上，且與平面相切，表明實驗體系具有最低的混合吉布斯自由能.這些結(jié)果證實了NRTL模型回歸的準確性和可靠性.

3 結(jié)論

本文使用COSMO-SAC模型繪制σ-profile圖，計算分子間相互作用能，選定乙二醇為分離共沸體系環(huán)己烷-正丙醇的萃取劑，測定了303.15 K、313.15 K和323.15 K下環(huán)己烷-正丙醇-乙二醇體系的液液相平衡數(shù)據(jù).三個溫度下的S值和D值都遠大于1，證明乙二醇可以萃取分離環(huán)己烷-正丙醇體系，且萃取劑可回收率較高.采用NRTL熱力學模型對實驗值和計算值進行擬合，得到的RMSD值小于1.0%，證明NRTL模型能成功關(guān)聯(lián)環(huán)己烷-正丙醇-乙二醇三元體系的LLE數(shù)據(jù).采用GUI-MATLAB方法進一步驗證了二元交互參數(shù)的可靠性.液液平衡相圖和二元交互作用參數(shù)可為萃取分離環(huán)己烷-正丙醇體系工藝設(shè)計和優(yōu)化提供重要依據(jù).

參考文獻

［1］ Ma S，Li J，Li L，et al.Liquid-liquid extraction of benzene and cyclohexane using sulfolane-based low transition temperature mixtures as solvents：Experiments and Simulation［J］.Energy amp; Fuels，2018，32（7）：8 006-8 015.

［2］Srinivasan C A，Saravabab C G.Emission reduction on ethanol-gasoline blend using fuel additives for an SI engine［J］.Energy Sources，Part A：Recovery，Utilization，and Environmental Effects，2013，35（12）：1 093-1 101.

［3］邢慶浩，李文秀，張 弢.［N1，8，8，8］［Ac］萃取精餾分離環(huán)己烷-正丙醇共沸體系研究［J］.山東化工，2023，52（12）：43-46，51.

［4］Klamt A.Conductor-like screening model for real solvents：A new approach to the quantitative calculation of solvation phenomena［J］.Physical Chemistry，1995，99（7）：2 224-2 235.

［5］Lin S T，Chang J，Wang S，et al.Prediction of vaporpressur-es and enthalpies of vaporization using a COSMO solvation model［J］.Physical Chemistry A，2004，108（36）：7 429-7 439.

［6］Lin S T，Wang L H，Chen W L，et al.Prediction ofmiscibility gaps in water/ether mixtures using COSMO-SAC model［J］.Fluid Phase Equilibria，2011，310（1）：9-24.

［7］Lin S T，Sandler S I.A priori phase equilibrium predictionfr-om a segment contribution solvation model［J］.Industrial amp; Engineering Chemistry Research，2002，41（5）：899-913.

［8］Zhang Y，Li H，Shan R，et al.Separation of isopropyl ether a-nd acetone using ionic liquids based on quantum chemistry calculation and liquid-liquid equilibrium［J］.Chemical Thermodynamics，2021，167：106 715.

［9］Dai F F，Xin K，Song Y H，et al.Liquid-liquid equilibria for the ternary system containing 1-Butanol + methoxy （methoxym-ethoxy） methane + water at temperatures of 303.15，323.15 and 343.15 K［J］.Fluid Phase Equilibria，2016，409：466-471.

［10］ 楊 沖，林旭楓，張金鋒，等.正己烷-異丙醇共沸體系液液相平衡數(shù)據(jù)測定及關(guān)聯(lián)［J］.化工學報，2020，71（7）：3 009-3 017.

［11］Wongsawa T，Hronec M，Soták T，et al.Ternary （liquid-liquid） equilibrium data of furfuryl alcohol with organic solvents at T=298.2 K：Experimental results and thermodynamic models［J］.Fluid Phase Equilibria，2014，369：88-96.

［12］Yang C，Cai Z，Chen H，et al.Measurement and correlation of liquid-liquid equilibria for extracting isopropanol from cyclohe-xane with different solvents［J］.Chemical amp; Engineering Data，2020，65（1）：140-145.

［13］Xu Y，Meng D，Li H，et al.Mechanism analysis for separation of cyclohexane and tert-butanol system via ionic liquids as e-xtractants and process optimization［J］.ACS Sustainable Chemistry amp; Engineering，2019，7（24）：19 984-19 992.

［14］Klamt A.COSMO-RS from quantum chemistry to fluid phase thermodynamics and drug design［J］.Chemical amp; Engineering Data，2005，51：1 480.

［15］楊 犁.基于COSMO-SAC模型和分子模擬預(yù)測復雜系統(tǒng)熱力學性質(zhì)［D］.上海：華東理工大學，2012.

［16］Wu Z，Liu C，Cheng H，et al.Tuned extraction andregenerati-on process for separation of hydrophobic compounds by aqueo-us ionic liquid［J］.Molecular Liquids，2020，308：113 032.

［17］曹江婷，陳布天，程曉亮，等.EG-EGDA分離萃取劑的篩選及其萃取工藝模擬［J］.陜西科技大學學報，2024，42（2）：103-109，115.

［18］Boys S F，Bernardi F.The calculation of small molecular inte-ractions by the differences of separate total energies.Some procedures with reduced errors［J］.Molecular Physics，1970，19（4）：553-566.

［19］Srensen J M，Magnussen T，Rasmussen P，et al.Liquid-liquid equilibrium data：Their retrieval，correlation and prediction Part II：Correlation［J］.Fluid Phase Equilibria，1979，3（1）：47-82.

［20］Yim J H，Park K W，Lim J S，et al.Liquid-liquid equilibrium measurements for the ternary system of water/2，3-butanediol/4-methyl-2-pentanol at various temperatures［J］.Chemical amp; Engineering Data，2019，64（9）：3 882-3 888.

［21］Zhu Z Y，Bai W，Xu Y，et al.Liquid-liquid extraction of met-hanol from its mixtures with hexane using three imidazolium-based ionic liquids［J］.Chemical Thermodynamics，2019，138：189-195.

［22］Alkhaldi K H A E，Al Tuwaim M S，F(xiàn)andary M S，et al.Separation of propylbenzene and n-alkanes from their mixtures using 4-methyl-N-butylpyridinium tetrafluoroborate as an ionic solvent at several temperatures［J］.Fluid Phase Equilibria，2011，309（1）：102-107.

［23］Dai F F，Cao J T，Liu N，et al.Determination and correlation of LLE data for n-hexane，ethyl acetate and different extractants［J］.Chemical Thermodynamics，2024，192：107 257.

［24］徐春明，殷夢凡，童 浩，等.離子液體萃取分離石腦油中芳烴的理論計算與試驗驗證［J］.中國石油大學學報（自然科學版），2023，47（5）：138-145.

［25］Renon H，Prausnitz J M.Local compositions in thermodynamic excess functions for liquid mixtures［J］.AIChE Journal，1968，14（1）：135-144．

［26］汪 勤，張冰劍，何 暢，等.環(huán)丁砜萃取精餾過程模擬分析及工藝參數(shù)優(yōu)化［J］.化工學報，2017，68（5）：1 969-1 976.

［27］莊志海，張建強，劉殿華.聚甲氧基二甲醚+水+正己烷三元體系的液液相平衡［J］.化工學報，2016，67（9）：3 545-3 551.

［28］ToledoI E P，F(xiàn)erreira Pinto L，Pedersen Voll F A，et al.Liquid-liquid equilibrium of the system {peanut biodiesel+glycerol+ethanol} at atmospheric pressure［J］.Chemical amp; Engineering Data，2019，64（5）：2 207-2 212．

［29］Díaz I，Rodríguez M，Gonz＇alez E J，et al.A simple and reli-able procedure to accurately estimate NRTL interaction parameters from liquid-liquid equilibrium data［J］.Chemical Engineering Science，2019，193：370-378.

［30］Marcilla A，Reyes Labarta J A，Olaya M M.Should we trust all the published LLE correlation parameters in phase equilibria necessity of their assessment prior to publication［J］.Fluid Phase Equilibria，2017，433：243-252.

［31］Reyes Labarta J A，Olaya M M，Velasco R，et al.Correlation of the liquid-liquid equilibrium data for specific ternary systems with one or two partially miscible binary subsystems［J］.Fluid Phase Equilibria，2009，278（1-2）：9-14.

【責任編輯：蔣亞儒】