張 圓，唐盛偉，張 濤

（四川大學(xué) 化學(xué)工程學(xué)院 四川省多相流傳質(zhì)與反應(yīng)工程重點(diǎn)實(shí)驗(yàn)室，四川 成都 610065）

預(yù)測(cè)乙烷/離子液體系統(tǒng)溶解平衡的UNIFAC模型

張 圓，唐盛偉，張 濤

（四川大學(xué) 化學(xué)工程學(xué)院 四川省多相流傳質(zhì)與反應(yīng)工程重點(diǎn)實(shí)驗(yàn)室，四川 成都 610065）

采用UNIFAC模型對(duì)乙烷在1-烷基-3-甲基咪唑雙三氟甲基磺酰亞胺鹽 （［RMIm］［Tf2N］，R=乙基（E）、丁基、己基、辛基（O）、癸基）離子液體中的113個(gè)溶解度數(shù)據(jù)進(jìn)行關(guān)聯(lián)，得到了乙烷與主基團(tuán)CH2—和［MIm］［Tf2N］間的交互參數(shù)，建立了預(yù)測(cè)乙烷在該類離子液體中溶解度的UNIFAC模型。利用該模型預(yù)測(cè)了乙烷在［EMIm］［Tf2N］和［OMIm］［Tf2N］兩種離子液體中的溶解度，預(yù)測(cè)值與文獻(xiàn)值吻合較好，平均誤差分別為4.01%和3.81%，最大誤差為9.86%。將該模型用于分析改變烷基鏈長對(duì)乙烷在［RMIm］［Tf2N］離子液體中溶解度的影響發(fā)現(xiàn)，在烷基鏈長較短時(shí)，增加其鏈長可有效提高乙烷的溶解度；但當(dāng)烷基鏈長較長時(shí)，其效果減弱。該模型可對(duì)乙烷在該類離子液體中的溶解度進(jìn)行有效預(yù)測(cè)，為其相關(guān)的傳質(zhì)分離過程提供相平衡數(shù)據(jù)。

乙烷；離子液體；溶解度；UNIFAC模型；1-烷基-3-甲基咪唑雙三氟甲基磺酰亞胺鹽

離子液體具有蒸汽壓極低、難揮發(fā)、良好的 化學(xué)穩(wěn)定性和熱穩(wěn)定性、不易燃及良好的導(dǎo)電性等優(yōu)良特性［1］，更重要的是其結(jié)構(gòu)可按實(shí)際要求進(jìn)行設(shè)計(jì)。離子液體已廣泛應(yīng)用于催化反應(yīng)、分離和電化學(xué)等眾多領(lǐng)域［2-8］。

離子液體系統(tǒng)的相平衡是對(duì)其在反應(yīng)、分離等涉及傳質(zhì)的過程中的應(yīng)用進(jìn)行研究、設(shè)計(jì)的基礎(chǔ)。如在離子液體催化的氣液反應(yīng)中，氣體在離子液體中的溶解度直接影響反應(yīng)的效率，進(jìn)一步會(huì)影響產(chǎn)物的收率；在氣體分離中，不同氣體組分在離子液體中的溶解度差異將決定其分離的難易程度。目前，已有許多學(xué)者測(cè)定了氣體在離子液體中的溶解度，但這遠(yuǎn)不能滿足實(shí)際應(yīng)用需求。單純從實(shí)驗(yàn)中獲得溶解度數(shù)據(jù)的工作量是巨大的，用相平衡模型預(yù)測(cè)氣體在離子液體中的溶解度對(duì)離子液體在化工過程中的應(yīng)用十分重要。UNIFAC基團(tuán)貢獻(xiàn)法［9］將分子劃分為基本基團(tuán)，只考慮基團(tuán)間的交互作用，特別適用于相平衡數(shù)據(jù)十分缺乏的離子液體系統(tǒng)。Lei等［10］將UNIFAC模型應(yīng)用于CO2/離子液體系統(tǒng)，回歸了CO2與22種離子液體基團(tuán)間的交互參數(shù)，并預(yù)測(cè)了CO2在混合離子液體中的溶解度，還回歸了CO和H2與離子基團(tuán)間的UNIFAC模型交互參數(shù)［11］。但在其他氣體/離子液體系統(tǒng)中，UNIFAC模型的應(yīng)用還較少。

本工作對(duì)乙烷/1-烷基-3-甲基咪唑雙三氟甲基磺酰亞胺鹽（［RMIm］［Tf2N］， R=乙基（E）、丁基（B）、己基（H）、辛基（O）、癸基（D））離子液體系統(tǒng)的UNIFAC模型基團(tuán)交互參數(shù)進(jìn)行了回歸，并預(yù)測(cè)了離子液體烷基鏈長對(duì)乙烷溶解度的影響。

1 UNIFAC模型

當(dāng)乙烷/離子液體系統(tǒng)處于相平衡時(shí)，有

式中，x1和y1分別為乙烷在液相和氣相中的摩爾分?jǐn)?shù)；φ1(T, p, y1)為乙烷在氣相中的逸度系數(shù)，由PR方程計(jì)算得到；p為系統(tǒng)壓力，kPa；為乙烷的飽和蒸氣壓，kPa；γ1為乙烷在離子液體中的活度系數(shù)，可由實(shí)驗(yàn)數(shù)據(jù)或UNIFAC模型求得。由于離子液體的蒸氣壓極低，氣相中幾乎不存在離子液體，所以氣相為純乙烷相，y1=1。

UNIFAC模型［9］是UNIQUAC模型與基團(tuán)分析法相結(jié)合而發(fā)展起來的活度系數(shù)模型，該模型將活度系數(shù)分為兩部分：組合項(xiàng)ln（是由分子大小和形狀引起的）和剩余項(xiàng)ln（是由分子間相互作用引起的）。其表達(dá)式為：組合項(xiàng)ln的計(jì)算式為：

式中，qi和ri分別為純組分分子的表面積參數(shù)和體積參數(shù)，其計(jì)算式為：

式中，xi為組分i的摩爾分?jǐn)?shù)；Xm為基團(tuán)m在溶液中的摩爾分?jǐn)?shù)；ψnm為基團(tuán)交互作用參數(shù)。

2 結(jié)果與討論

2.1 UNIFAC模型參數(shù)回歸

在應(yīng)用UNIFAC模型時(shí)，需將組分拆分為基本基團(tuán)。將氣態(tài)乙烷分子作為一個(gè)主基團(tuán)；對(duì)于離子液體，由于離子對(duì)之間有較強(qiáng)的靜電作用［12］，所以將陽離子甲基咪唑環(huán)和陰離子看作一個(gè)整體，將其作為主基團(tuán)。以［BMIm］［Tf2N］為例，將其拆分為1個(gè)［MIm］［Tf2N］、1個(gè)CH3—和3個(gè)CH2—共5個(gè)基團(tuán)（見圖1）。離子液體的基團(tuán)拆分見表1。

圖1 ［BMIm］［Tf2N］的基團(tuán)拆分示意圖Fig.1 Diagram of split groups of［BMIm］［Tf2N］.

表1 離子液體的基團(tuán)拆分Table 1 Groups of the ionic liquids for the UNIFAC model

UNIFAC模型計(jì)算中所需的基團(tuán)表面積參數(shù)（Q）和體積參數(shù)（R）取自文獻(xiàn)［12-13］，見表2。

表2 基團(tuán)體積參數(shù)和表面積參數(shù)Table 2 Volume parameters(R) and surface area parameters(Q) of the groups

UNIFAC模型交互參數(shù)回歸中，目標(biāo)函數(shù)OF定義為：

式中，xexp和xcal分別為乙烷在離子液體中的溶解度實(shí)驗(yàn)值和UNIFAC模型計(jì)算值；N為實(shí)驗(yàn)數(shù)據(jù)點(diǎn)數(shù)，N=113。

為提高交互參數(shù)回歸效率和減少交互參數(shù)回歸個(gè)數(shù)，CH2—（m）與［MIm］［Tf2N］（n）之間的交互參數(shù)直接使用文獻(xiàn)值［12］：amn=400.89，anm=145.8。C2H6與CH2—和［MIm］［Tf2N］之間的交互參數(shù)通過擬合乙烷在［EMIm］［Tf2N］［14］，［BMIm］［Tf2N］［15］，［HMIm］［Tf2N］［16-17］，［DMIm］［Tf2N］［15］離子液體中的溶解度實(shí)驗(yàn)數(shù)據(jù)得到。通過求解目標(biāo)函數(shù)式（10）得到的UNIFAC模型基團(tuán)交互參數(shù)回歸值見表3?；貧w時(shí)最小化目標(biāo)函數(shù)值為0.042 9。

乙烷在離子液體中的溶解度實(shí)驗(yàn)值與UNIFAC模型回歸值的比較見圖2。由圖2可見，絕大多數(shù)數(shù)據(jù)點(diǎn)的相對(duì)誤差在10%以內(nèi)，只有少數(shù)數(shù)據(jù)點(diǎn)的相對(duì)誤差大于10%，但未超過15%。

表3 UNIFAC模型基團(tuán)交互參數(shù)回歸值Table 3 Fitted group interaction parameters for the UNIFAC model

圖2 乙烷在離子液體中的溶解度實(shí)驗(yàn)值與UNIFAC模型回歸值的比較Fig.2 Comparison between the experimental data(xexp) and the values calculated by the UNIFAC model(xcal) for the solubility of ethane in the ionic liquids.

UNIFAC模型回歸的平均相對(duì)誤差見表4。由表4可知，在乙烷/［EMIm］［Tf2N］系統(tǒng)中平均相對(duì)誤差最小，為1.36%；在乙烷/［HMIm］［Tf2N］系統(tǒng)中平均相對(duì)誤差最大，為7.35%。考慮到溶解度數(shù)據(jù)來源于不同的文獻(xiàn)，存在實(shí)驗(yàn)誤差等因素，該預(yù)測(cè)結(jié)果是可以接受的，模型參數(shù)回歸效果較好。

表4 UNIFAC模型回歸的平均相對(duì)誤差Table 4 Average relative deviation(ARD) of regression with the UNIFAC model

2.2 利用溶解度數(shù)據(jù)驗(yàn)證UNIFAC模型參數(shù)

為驗(yàn)證回歸的基團(tuán)交互參數(shù)的有效性，采用乙烷在［EMIm］［Tf2N］［18］和［OMIm］［Tf2N］［15］中的溶解度數(shù)據(jù)與UNIFAC模型預(yù)測(cè)值進(jìn)行比較（見圖3），其相對(duì)誤差見表5。從表5可看出，UNIFAC模型預(yù)測(cè)值與實(shí)驗(yàn)值吻合得較好，最大相對(duì)誤差小于10%；乙烷/［EMIm］［Tf2N］系統(tǒng)的平均相對(duì)誤差為4.01%，乙烷/［OMIm］［Tf2N］系統(tǒng)的平均相對(duì)誤差為3.81%。UNIFAC模型中基團(tuán)交互參數(shù)的有效性與回歸所用實(shí)驗(yàn)數(shù)據(jù)的準(zhǔn)確性和數(shù)量有很大關(guān)系，實(shí)驗(yàn)數(shù)據(jù)越精確，數(shù)量越多，回歸得到的基團(tuán)交互參數(shù)越有效。目前乙烷在離子液體中的溶解度數(shù)據(jù)還較少，還需廣大科研工作者繼續(xù)完善。

表5 乙烷在離子液體中的溶解度實(shí)驗(yàn)值與UNIFAC模型預(yù)測(cè)值的相對(duì)偏差Table 5 Relative deviation of the values predicted by the UNIFAC model and the experimental data

圖3 乙烷在離子液體中的溶解度實(shí)驗(yàn)值與UNIFAC模型預(yù)測(cè)值的比較Fig.3 Comparison between the experimental data and the values predicted by the UNIFAC model for the solubility of ethane in the ionic liquids.

2.3 UNIFAC模型預(yù)測(cè)離子液體結(jié)構(gòu)對(duì)乙烷溶解度的影響

運(yùn)用回歸的基團(tuán)交互參數(shù)，可預(yù)測(cè)離子液體烷基鏈長與乙烷溶解性能的關(guān)系。以亨利常數(shù)表示乙烷在離子液體中的溶解度，其值越小，則表示乙烷的溶解度越大。亨利常數(shù)定義為：

式中，H（T）為溫度T下乙烷在離子液體中的亨利常數(shù)，MPa；為乙烷在離子液體中的無限稀釋活度系數(shù)，由UNIFAC模型計(jì)算；(T)為乙烷在溫度T下的飽和蒸氣壓，MPa；f1為乙烷的逸度，MPa；x1為乙烷在離子液體中的摩爾分?jǐn)?shù)。

預(yù)測(cè)了298.15 K時(shí)乙烷在不同烷基鏈長離子液體中的亨利常數(shù)，并與文獻(xiàn)值進(jìn)行比較（見圖4）。從圖4可看出，在考慮了實(shí)驗(yàn)方法不同以及實(shí)驗(yàn)裝置誤差的情況下，乙烷在［EMIm］［Tf2N］，［BMIm］［Tf2N］，［HMIm］［Tf2N］，［OMIm］·［Tf2N］，［DMIm］［Tf2N］等5種離子液體中的亨利常數(shù)的UNIFAC模型預(yù)測(cè)值與文獻(xiàn)值［14-19］吻合得很好，且其隨烷基鏈長的變化規(guī)律一致。這進(jìn)一步說明了UNIFAC模型能對(duì)乙烷在離子液體中的溶解度進(jìn)行比較準(zhǔn)確的預(yù)測(cè)。從圖4還可看出，離子液體烷基側(cè)鏈上的碳原子數(shù)目對(duì)乙烷在離子液體中的亨利常數(shù)影響很大。對(duì)于不同長度的烷基鏈，當(dāng)碳原子數(shù)目變化相同時(shí)，對(duì)亨利常數(shù)的影響程度差別很大。如當(dāng)烷基側(cè)鏈上碳原子數(shù)目由2增至4時(shí)，單位碳原子對(duì)亨利常數(shù)的改變量為2.12 MPa；而當(dāng)碳原子數(shù)目由18增至20時(shí)，單位碳原子對(duì)亨利常數(shù)的改變量僅為0.13 MPa?？傮w上，隨烷基鏈長的增加，乙烷在離子液體中的亨利常數(shù)呈降低趨勢(shì)，且降幅先快后慢，逐漸趨于平穩(wěn)。這表明隨烷基鏈長的增加，烷基鏈長變化對(duì)乙烷溶解度的影響逐漸減弱。所以在離子液體烷基鏈長較短時(shí)，增加其鏈長能有效提高乙烷在其中的溶解度；但當(dāng)烷基鏈長較長時(shí)，增加其鏈長對(duì)乙烷溶解度的影響減弱，此時(shí)需采取其他手段來提高乙烷的溶解度，如改變其陰離子結(jié)構(gòu)等。

圖4 298.15 K時(shí)UNIFAC模型預(yù)測(cè)乙烷在離子液體中的HFig.4 Henry constants(H) of ethane in ionic liquids at 298.15 K predicted by the UNIFAC model.

3 結(jié)論

1）對(duì)乙烷/［RMIm］［Tf2N］離子液體系統(tǒng)的UNIFAC模型基團(tuán)交互參數(shù)進(jìn)行了回歸，擬合效果較好，最小化目標(biāo)函數(shù)值為0.042 9。將回歸的基團(tuán)交互參數(shù)用于［EMIm］［Tf2N］和［OMIm］［Tf2N］離子液體中乙烷溶解度的模擬計(jì)算，計(jì)算結(jié)果與文獻(xiàn)值的平均相對(duì)誤差分別為4.01%和3.81%。

2）運(yùn)用擬合的UNIFAC模型，預(yù)測(cè)了298.15 K時(shí)乙烷在不同烷基鏈長離子液體中的亨利常數(shù)，預(yù)測(cè)值與實(shí)驗(yàn)值的吻合程度高。通過分析離子液體烷基鏈長對(duì)乙烷溶解度的影響發(fā)現(xiàn)，隨離子液體烷基鏈長的增加，改變烷基鏈長對(duì)乙烷溶解能力的影響逐漸減弱。

3）所建立的乙烷/［RMIm］［Tf2N］離子液體系統(tǒng)的UNIFAC模型可對(duì)乙烷在該類離子液體中的溶解度進(jìn)行有效預(yù)測(cè)，為其相關(guān)的傳質(zhì)分離過程提供相平衡數(shù)據(jù)。

符 號(hào) 說 明

a UNIFAC模型基團(tuán)交互參數(shù)

f1乙烷的逸度，MPa

H 亨利常數(shù)

N 實(shí)驗(yàn)數(shù)據(jù)點(diǎn)數(shù)

p 體系壓力，kPa

Q 基團(tuán)表面積參數(shù)

R 基團(tuán)體積參數(shù)

T 溫度，K

X 基團(tuán)在溶液中的摩爾分?jǐn)?shù)

x1乙烷在液相中的摩爾分?jǐn)?shù)

xexp乙烷在離子液體中的溶解度（摩爾分?jǐn)?shù)）

xcalUNIFAC模型計(jì)算得到的溶解度（摩爾分?jǐn)?shù)）

Γk基團(tuán)k的剩余活度系數(shù)

γ1乙烷在液相中的活度系數(shù)

θ 基團(tuán)表面積分?jǐn)?shù)

φ1乙烷在氣相中的逸度系數(shù)

ψ 基團(tuán)交互作用參數(shù)

下角標(biāo)

i，j 組分，i，j = 1，2，…

k，m，n 基團(tuán)，k，m，n=1，2，…

［1］ Rogers R D，Seddon K R. Ionic Liquids：Solvents of the Future?［J］. Science，2003，302（5646）：792 - 793.

［2］ Sheldon R. Catalytic Reactions in Ionic Liquids［J］. Chem Commun，2001（23）：2399 - 2407.

［3］ Berthod A，Ruiz-Angel M，Carda-Broch S. Ionic Liquids in Separation Techniques［J］. J Chromatogr A，2008，1184（1）：6 - 18.

［4］ Galiński M，Lewandowski A，St?pniak I. Ionic Liquids as Electrolytes［J］. Electrochim Acta，2006，51（26）：5567 - 5580.

［5］ 徐新，羅國華，王莉. 氯鋁酸離子液體催化苯與氯乙烷合成乙苯［J］. 石油化工，2012，41（1）：33 - 36.

［6］ 劉春萍，劉剛，胡玉才，等. 1-芐基-4-甲基吡啶鹽離子液體的合成及其電導(dǎo)性能的研究［J］. 石油化工，2011，40（7）：764 - 769.

［7］ 王吉林，王璐璐. 離子液體1-丁基-3-甲基咪唑磷酸二丁酯的制備及其催化酯化反應(yīng)的性能［J］. 石油化工，2011，40（6）：635 - 639.

［8］ 田征，崔現(xiàn)寶，彭艷枚，等. 離子液體萃取精餾分離乙酸甲酯-甲醇共沸物［J］. 石油化工，2013，42（5）：537 - 541.

［9］ Fredenslund A，Jones R L，Prausnitz J M. Group-Contribution Estimation of Activity Coeffcients in Nonideal Liquid Mixtures［J］. AIChE J，1975，21（6）：1086 - 1099.

［10］ Lei Zhigang，Dai Chengna，Wang Wei，et al. UNIFAC Model for Ionic Liquid-CO2Systems［J］. AIChE J，2014，60（2）：716 - 729.

［11］ Lei Zhigang，Dai Chengna，Yang Qian，et al. UNIFAC Model for Ionic Liquid-CO（H2） Systems：An Experimental and Modeling Study on Gas Solubility［J］. AIChE J，2014，60（12）：4222 - 4231.

［12］ Lei Zhigang，Zhang Jiguo，Li Qunsheng，et al. UNIFAC Model for Ionic Liquids［J］. Ind Eng Chem Res，2009，48（5）：2697 - 2704.

［13］ Banerjee T，Singh M K，Sahoo R K，et al. Volume，Surface and UNIQUAC Interaction Parameters for Imidazolium based Ionic Liquids via Polarizable Continuum Model［J］. Fluid Phase Equilib，2005，234（1）：64 - 76.

［14］ Hong Gao，Jacquemin J，Deetlefs M，et al. Solubility of Carbon Dioxide and Ethane in Three Ionic Liquids Based on the Bis{（Trifluoromethyl） Sulfonyl}Imide Anion［J］. Fluid Phase Equilib，2007，257（1）：27 - 34.

［15］ Gomes M F C，Pison L，Pensado A S，et al. Using Ethane and Butane as Probes to the Molecular Structure of 1-Alkyl-3-Methylimidazolium Bis［（Trifluoromethyl） Sulfonyl］Imide Ionic Liquids［J］. Faraday Discuss，2012，154：41 - 52.

［16］ Florusse L J，Raeissi S，Peters C J. High-Pressure Phase Behavior of Ethane with 1-Hexyl-3-Methylimidazolium Bis（Trifluoromethylsulfonyl）Imide［J］. J Chem Eng Data，2008，53（6）：1283 - 1285.

［17］ Costa G M. Low-Pressure Solubility and Thermodynamics of Solvation of Carbon Dioxide，Ethane，and Hydrogen in 1-Hexyl-3-Methylimidazolium Bis（Trifluoromethylsulfonyl）Amide Between Temperatures of 283 K and 343 K［J］. J Chem Eng Data，2007，52（2）：472 - 475.

［18］ Liu Xiangyang，Afzal W，Prausnitz J M. Solubilities of Small Hydrocarbons in Tetrabutylphosphonium Bis（2，4，4-Trimethylpentyl）Phosphinate and in 1-Ethyl-3-Methylimidazolium Bis（Trifuoromethylsulfonyl）Imide［J］. Ind Eng Chem Res，2013，52（42）：14975 - 14978.

［19］ Anthony J L，Anderson J L，Maginn E J，et al. Anion Effects on Gas Solubility in Ionic Liquids［J］. J Phys Chem B，2005，109（13）：6366 - 6374.

（編輯 安 靜）

UNIFAC Model for Prediction of Dissolution Equilibrium of Ethane/Ionic Liquid Systems

Zhang Yuan，Tang Shengwei，Zhang Tao
（Multi-Phases Transfer and Reaction Engineering Laboratory，College of Chemical Engineering，Sichuan University，Chengdu Sichuan 610065，China）

113 Solubility data of ethane in ionic liquids 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide(［RMIm］［Tf2N］，R=ethyl(E)，butyl，hexyl，octyl(O)，decyl) were fitted by the UNIFAC model and the interaction parameters between C2H6，and CH2— and［MIm］［Tf2N］ were obtained. The UNIFAC model was used to predict the solubilities of ethane in［EMIm］［Tf2N］ and ［OMIm］［Tf2N］. The results showed that the predicated data agreed with the experimental data well，the average relative deviations were 4.01% and 3.81% respectively，and the largest relative deviation was 9.86%. The effect of alkyl chain length on the ethane solubility in the ionic liquids was investigated by means of the UNIFAC model. It was found that increasing the alkyl chain length could enhance the dissolution effectively when the length was relatively short，but when the carbon number was larger than 10，the effects became weaker. The obtained UNIFAC model is useful for predicting the ethane solubility in［RMIm］［Tf2N］.

ethane；ionic liquids；solubility；UNIFAC model；1-alkyl-3-methylimidazolium bis(trifuoromethylsulfonyl)imide

1000 - 8144（2015）10 - 1212 - 06

TQ 013.1

A

2015 - 04 - 13；［修改稿日期］ 2015 - 06 - 27。

張圓（1989—），男，四川省武勝縣人，碩士，電郵 sccdzy2008@163.com。聯(lián)系人：唐盛偉，電話 028 - 85405201，電郵tangdynasty@scu.edu.cn。

國家自然科學(xué)基金項(xiàng)目（21276163）。