方靜，張淑婷，李婷婷，李春利

?

離子液體用于燃油萃取脫硫的選擇與過程優(yōu)化

方靜，張淑婷，李婷婷，李春利

（河北工業(yè)大學(xué)化工學(xué)院，天津 300130）

針對(duì)基于COSMO-SAC模型分子設(shè)計(jì)方法的準(zhǔn)確性問題，采用離子液體脫硫機(jī)理分析和實(shí)驗(yàn)的方法對(duì)其進(jìn)行了驗(yàn)證，即對(duì)[HMIM][BF4]、[HMIM][PF6]、[BMIM][BF4]、[BMIM][PF6]、[EMIM][BF4]、[EMIM][PF6] 6種離子液體作萃取劑時(shí)的脫硫效果進(jìn)行了脫硫機(jī)理的分析和實(shí)驗(yàn)的驗(yàn)證，得到的脫硫性能排序與離子液體分子設(shè)計(jì)結(jié)果基本一致，且均認(rèn)為[HMIM][PF6]脫硫率較高。以[HMIM][PF6]為萃取劑，通過液相色譜法測(cè)定萃取后的液相組成，考察萃取時(shí)間、萃取溫度、劑油比3個(gè)因素對(duì)脫硫率、分配系數(shù)和選擇性系數(shù)的影響。通過正交實(shí)驗(yàn)設(shè)計(jì)確定了萃取時(shí)間40 min、萃取溫度20℃、劑油比2:1為較優(yōu)操作條件，單次脫硫率為72.74%，四級(jí)萃取可將模型油的含硫量由1200 μg·g?1降至6.98 μg·g?1，符合國(guó)Ⅴ標(biāo)準(zhǔn)。

離子液體；萃??；脫硫；模擬；正交實(shí)驗(yàn)

引 言

汽車尾氣中含硫氧化物是城市空氣的主要污染物之一。脫除燃油中的含硫化合物成為遏制汽車尾氣污染大氣的主要途徑[1-4]。

萃取脫硫（extraction desulfurization，EDS）是利用硫化物在萃取劑與燃料油中溶解度不同而進(jìn)行分離的脫硫工藝[5]。萃取脫硫的操作條件溫和（低溫、低壓或常溫、常壓），脫硫后油品的性質(zhì)不發(fā)生改變，不需要消耗氫氣、氧化劑等額外的物質(zhì)，設(shè)備簡(jiǎn)單，流程短，投資成本低。然而由于受到萃取劑性質(zhì)的限制，該方法的脫硫效率較低，而且傳統(tǒng)溶劑對(duì)環(huán)境的污染也較為嚴(yán)重[6-7]。因此，對(duì)于溶劑萃取脫硫工藝而言，開發(fā)一種高效、低污染的綠色溶劑是當(dāng)前的研究重點(diǎn)。

離子液體作為一種綠色溶劑，對(duì)傳統(tǒng)的加氫脫硫方法難以去除的噻吩類硫化物有較好的脫除效果，而且離子液體萃取脫硫可在室溫下進(jìn)行，操作簡(jiǎn)單，離子液體易于回收利用，克服了傳統(tǒng)脫硫方法工藝復(fù)雜、設(shè)備投資大等缺點(diǎn)[8-11]。因此，離子液體萃取脫硫具有良好的工業(yè)應(yīng)用前景[12]。

1995年Klamt[13]提出了COSMO-RS模型，即真實(shí)溶劑似導(dǎo)體屏蔽模型，該模型只需知道分子的結(jié)構(gòu)信息，便可利用分子表面屏蔽電荷密度分布計(jì)算得到體系的物理化學(xué)性質(zhì)，如飽和蒸氣壓和活度因子等[14-16]。其后Lin等[17-19]提出了與之類似的COSMO-SAC片段活度因子模型，通過引入Staverman-Guggenheim組合項(xiàng)克服了COSMO-RS在預(yù)測(cè)汽液平衡數(shù)據(jù)時(shí)不符合熱力學(xué)一致性原理的缺點(diǎn)[20]。但基于COSMO-SAC模型的離子液體設(shè)計(jì)方法是一種估算方法，在計(jì)算無(wú)限稀釋活度因子時(shí)存在一定誤差[21]。

本文通過離子液體脫硫機(jī)理和實(shí)驗(yàn)的方法對(duì)基于COSMO-SAC模型的分子設(shè)計(jì)方法的準(zhǔn)確性進(jìn)行驗(yàn)證，并就[EMIM][BF4]、[EMIM][PF6]、[BMIM][BF4]、[BMIM][PF6]、[HMIM][BF4]、[HMIM][PF6] 6種離子液體中模擬萃取效果最好的[HMIM][PF6]作為萃取脫硫溶劑，對(duì)不同操作條件下脫除燃油中的含硫化合物的效果進(jìn)行研究，旨在找出簡(jiǎn)單操作下能達(dá)到高脫硫率的最佳操作條件。

2 實(shí)驗(yàn)材料和方法

2.1 實(shí)驗(yàn)材料

苯并噻吩，純度98%，天津希恩思生化科技有限公司；正辛烷，分析純，天津福晨化學(xué)試劑廠；1-乙基-3-甲基咪唑四氟硼酸鹽[EMIM][BF4]，純度99%，上海成捷化學(xué)有限公司；1-乙基-3-甲基咪唑六氟磷酸鹽[EMIM][PF6]，純度99%，上海成捷化學(xué)有限公司；1-丁基-3-甲基咪唑四氟硼酸鹽[BMIM][BF4]，純度99%，上海成捷化學(xué)有限公司；1-丁基-3-甲基咪唑六氟磷酸鹽[BMIM][PF6]，純度99%，上海成捷化學(xué)有限公司；1-己基-3-甲基咪唑四氟硼酸鹽[HMIM][BF4]，純度99%，上海成捷化學(xué)有限公司；1-己基-3-甲基咪唑六氟磷酸鹽[HMIM][PF6]，純度99%，上海成捷化學(xué)有限公司。

選擇苯并噻吩為含硫化合物，正辛烷為溶劑，配制模型油。在精確度為0.0001 g的電子天平上準(zhǔn)確稱取1.0000 g苯并噻吩，在250 ml的錐形瓶中加入198 g正辛烷，與苯并噻吩配制成含硫量1200 μg·g?1的模型油。

2.2 儀器設(shè)備和方法

集熱式恒溫加熱磁力攪拌器，DF-101S型，上海翔雅儀器設(shè)備有限公司；高效液相色譜，SPD-M20A型，日本島津公司；旋轉(zhuǎn)蒸發(fā)器，RE-5299型，天津科諾儀器設(shè)備有限公司；電子天平，TC328A型，上海天平儀器廠。

實(shí)驗(yàn)中苯并噻吩和離子液體的測(cè)定采用液相色譜法。分別配置苯并噻吩含量為1、2、3、4、5 mg·L?1，[HMIM][PF6]含量為400、600、720、800、980 mg·L?1的標(biāo)準(zhǔn)溶液，在流動(dòng)相90%甲醇、10%水，流速1 ml·min?1下進(jìn)樣分析，以峰面積為縱坐標(biāo)，各物質(zhì)質(zhì)量濃度（mg·L?1）為橫坐標(biāo)，通過線性回歸，分別得到苯并噻吩和[HMIM][PF6]的線性回歸方程為

= 46088.7? 2705.9(= 0.99993)

= 12976.1+ 2.4(= 0.99912)

利用所得回歸方程，即可將檢測(cè)結(jié)果中的峰面積轉(zhuǎn)化為苯并噻吩和[HMIM][PF6]的濃度。

常壓下，在50 ml的錐形瓶中加入含硫量為1200 μg·g?1的模型油和離子液體，一定溫度下密閉磁力攪拌，萃取平衡后室溫下靜置，待模型油與離子液體萃取劑清晰分層。用注射器抽取0.1 ml上、下層溶液并稱重，再分別稀釋至50、100倍，即可用液相色譜儀檢測(cè)稀釋后苯并噻吩、離子液體的含量A和S，正辛烷含量B=?50A?100S。由此可得上下層溶液三相組成。

實(shí)驗(yàn)中離子液體在正辛烷中的溶解度同樣采用液相色譜法測(cè)定。配制25℃下離子液體-正辛烷飽和溶液，取0.1 ml并稱重，稀釋至100倍用液相色譜儀檢測(cè)稀釋后離子液體的含量S，正辛烷含量B=?100S，25℃下離子液體在正辛烷中的溶解度為P=100SB。

3 結(jié)果與討論

3.1 基于COSMO-SAC模型分子設(shè)計(jì)方法的準(zhǔn)確性驗(yàn)證

3.1.1 COSMO-SAC模型[20]及分子設(shè)計(jì) Lin等[17]在2002年提出COSMO-SAC活度系數(shù)模型，認(rèn)為混合物中組分的活度系數(shù)可以分解為兩個(gè)部分，即表面電荷片段對(duì)活度系數(shù)的貢獻(xiàn)和分子形狀大小差異對(duì)活度系數(shù)的貢獻(xiàn)。前者采用與COSMO-RS模型類似的片段活度系數(shù)計(jì)算方法表示，后者則采用Staverman-Guggenheim組合項(xiàng)來(lái)表示，即

式中，γ/S、分別為組分在混合物中的活度系數(shù)和Staverman-Guggenheim組合項(xiàng)活度系數(shù)；n為混合物中組分的物質(zhì)的量；p()為組分的表面電荷密度分布函數(shù)，即-Profile，該模型中純物質(zhì)的-Profile取自文獻(xiàn)[22]；s()、Г()分別為混合物和純物質(zhì)中表面電荷密度等于的片段對(duì)活度系數(shù)的貢獻(xiàn)值大小，均可按式(2)計(jì)算

(2)

式中，下標(biāo)S/表示混合物S或組分；Δ(σ，σ)為交換能，即從中性對(duì)中得到(σ，σ)荷電對(duì)所需的能量；為氣體常數(shù)；為體系溫度。

混合物中的表面電荷密度分布函數(shù)S()為所有純組分P()的加權(quán)和，即

式中，x為混合物中組分的摩爾分?jǐn)?shù)；A為組分中所有表面電荷的空穴表面積。

Staverman-Guggenheim組合項(xiàng)活度系數(shù)則由式(4)計(jì)算

式中，為配位數(shù)，通常取值為10；θ和φ分別為組分的表面積分?jǐn)?shù)和體積分?jǐn)?shù)

(5)

式中，r和分別為組分的量綱1分子體積和表面積，可以通過COSMO計(jì)算得到。l按式(6)計(jì)算

l=(/2)[(r-q)(r-1)] (6)

用COSMO-SAC模型對(duì)[EMIM][BF4]、[EMIM][PF6]、[BMIM][BF4]、[BMIM][PF6]、[HMIM][BF4]、[HMIM][PF6] 6種離子液體在30℃下進(jìn)行量子化學(xué)計(jì)算，得到其無(wú)限稀釋活度因子。選取分配系數(shù)A、選擇性系數(shù)、溶解度P作為萃取劑的性能指標(biāo)[23]，計(jì)算公式如式(7)～式(9)所示

(8)

(9)

式中，∞A,B為組分A在組分B中的無(wú)限稀釋活度因子；∞A,S為組分A在萃取劑S中的無(wú)限稀釋活度因子；∞B,S為組分B在萃取劑S中的無(wú)限稀釋活度因子。

分配系數(shù)A表示萃取達(dá)到平衡時(shí)，組分A在萃取劑S和組分B中的比值，比值越大，萃取能力越強(qiáng)；選擇性系數(shù)越大，萃取劑對(duì)兩組分的分離效果越好；溶解度P表示組分A在萃取劑S中的溶解能力，溶解度越大，分離效果也越好。

根據(jù)王秀麗等[24]提出的模糊綜合評(píng)判方法，對(duì)6種離子液體的脫硫性能進(jìn)行評(píng)價(jià)和排序。設(shè)萃取劑評(píng)判因素集= {A,,P}，對(duì)中所有因素均設(shè)評(píng)語(yǔ)集= {優(yōu)}。萃取劑的模糊綜合評(píng)判公式為== {1,2,…,}，其中為權(quán)重分配集，反映的是各因素的權(quán)重，這里考慮到選擇性系數(shù)、分配系數(shù)和溶解度的重要性基本相同，因此令= [1/3, 1/3, 1/3][21]；為綜合評(píng)判因子；為模糊轉(zhuǎn)換矩陣

是第個(gè)被評(píng)判對(duì)象中第個(gè)元素對(duì)的隸屬度。這里規(guī)定隸屬函數(shù)如式(11)～式(13)所示

(11)

(13)

基于COSMO-SAC模型的分子設(shè)計(jì)結(jié)果見表1，從6種離子液體的綜合評(píng)判因子大小可以看出，離子脫硫性能順序?yàn)椋篬HMIM]+> [BMIM]+> [EMIM]+，PF6?> BF4?，且離子液體中綜合性能較好的為[HMIM][PF6]。6種離子液體結(jié)構(gòu)如圖1所示。

表1 離子液體作為萃取脫硫劑的分子設(shè)計(jì)結(jié)果

3.1.2 離子液體脫硫機(jī)理驗(yàn)證 Zhang等[25]采用NMR光譜研究了噻吩在離子液體中的吸附作用，結(jié)果表明離子液體在溶解噻吩分子過程中進(jìn)行重新組裝，這種重新組裝的能力主要取決于噻吩和離子液體中陽(yáng)離子的相互作用，這主要?dú)w因于咪唑陽(yáng)離子與噻吩分子的芳香環(huán)間的電流作用[26]。

由Mochizuki等[27]、Nie等[28]的研究可知，在陰離子相同的情況下，隨著陽(yáng)離子上烷基鏈取代基中碳原子個(gè)數(shù)的增加，離子液體脫硫能力增強(qiáng)。這是因?yàn)椋环矫骐S著陽(yáng)離子上烷基鏈取代基的增長(zhǎng)，離子液體分子結(jié)構(gòu)增大，陰陽(yáng)離子之間的庫(kù)侖作用減弱，有利于噻吩分子插入離子液體的動(dòng)態(tài)分子結(jié)構(gòu)中，形成液相包合物，促進(jìn)萃取過程的進(jìn)行；另一方面，咪唑環(huán)的對(duì)稱性越差，陽(yáng)離子極化性就越強(qiáng)，與噻吩類硫化物的絡(luò)合作用就越強(qiáng)，并且烷基是給電子基，可以增強(qiáng)硫原子與咪唑環(huán)的π-π作用[29]。所以陽(yáng)離子脫硫效果順序?yàn)椋篬HMIM]+> [BMIM]+> [EMIM]+。

Mochizuki等[27]、胡松青等[30]也發(fā)現(xiàn)，在陽(yáng)離子相同的情況下，陰離子體積增大也會(huì)使脫硫率增大。對(duì)此的解釋是，硫化物在離子液體中的溶解度受空間效應(yīng)的影響，陰離子體積增大使得陽(yáng)離子不會(huì)靠得很近，離子液體與芳香環(huán)間的電流作用降低，“堆垛”結(jié)構(gòu)沒有那么緊密了，可插入的硫化物分子數(shù)增多，脫硫率提高[31]。所以，陰離子脫硫效果順序?yàn)镻F6?> BF4?。

由此可見，離子液體的結(jié)構(gòu)影響其萃取脫硫效果，在一定限度內(nèi)增大陰陽(yáng)離子的體積有助于提高脫硫率[31]。所以從脫硫機(jī)理來(lái)看，脫硫性能排序?yàn)閇HMIM]+> [BMIM]+> [EMIM]+，PF6?> BF4?。且6種離子液體中，[HMIM][PF6]的脫硫性能最好，這與基于COSMO-SAC模型的分子設(shè)計(jì)結(jié)果保持一致，證明基于COSMO-SAC模型的離子液體設(shè)計(jì)方法是有效的。但是需要注意的是，由于COSMO-SAC模型是一種估算方法，在計(jì)算無(wú)限稀釋活度因子時(shí)存在一定誤差，因此該萃取劑分子設(shè)計(jì)方法是一種初步選擇方法，萃取劑的綜合性能仍需要實(shí)驗(yàn)進(jìn)行驗(yàn)證[21]。

3.1.3 脫硫?qū)嶒?yàn)驗(yàn)證 在錐形瓶中分別加入等物質(zhì)的量的離子液體和一定量的模型油，在30℃下密閉磁力攪拌進(jìn)行脫硫?qū)嶒?yàn)，萃取達(dá)到平衡后室溫下靜置。用液相色譜檢測(cè)反應(yīng)前后模型油相中的硫含量，考察脫硫率（由于[EMIM][PF6]常溫下呈固態(tài)，其萃取溫度設(shè)為50℃）。6種離子液體的脫硫率數(shù)據(jù)見表2。

從表2中6種離子液體的脫硫率大小可以看出，其脫硫性能排序與基于COSMO-SAC模型的分子設(shè)計(jì)結(jié)果基本一致，均認(rèn)為[HMIM]+> [BMIM]+> [EMIM]+，PF6?> BF4?，且[HMIM][PF6]的脫硫效果最好。

表2 6種離子液體單次脫硫率

綜上所述，通過離子液體脫硫機(jī)理和實(shí)驗(yàn)兩種方法對(duì)基于COSMO-SAC模型分子設(shè)計(jì)方法的準(zhǔn)確性進(jìn)行驗(yàn)證，證明基于COSMO-SAC模型分子設(shè)計(jì)方法具有一定的準(zhǔn)確性。

6種離子液體在正辛烷中的溶解度數(shù)據(jù)見表3。從表中可以看出，6種離子液體在正辛烷中的溶解度都非常小，雖然[HMIM][PF6]在正辛烷中的溶解度較其他5種離子液體略大，但這并不影響分離過程中對(duì)萃取劑溶解度的要求。

表3 298.15 K離子液體在正辛烷中的溶解度

下面將采用分子設(shè)計(jì)方法得出的綜合評(píng)判因子最高的[HMIM][PF6]離子液體作為萃取劑，對(duì)脫除模型油中的苯并噻吩過程進(jìn)行優(yōu)化。

3.2 [HMIM][PF6]脫硫的過程優(yōu)化

3.2.1 萃取溫度對(duì)脫硫效果的影響 以質(zhì)量比1:1加入一定量的[HMIM][PF6]和模型油，在不同溫度下密封攪拌30 min后，室溫下靜置，測(cè)定上下層中的三相組成。分析萃取溫度對(duì)脫硫效果的影響，實(shí)驗(yàn)結(jié)果見表4。

表4 溫度對(duì)萃取效果的影響

從表4可以看出，從20℃到60℃，隨著溫度的升高，離子液體脫硫率逐漸降低，這是因?yàn)殡S著溫度的升高，兩相間互溶度增大，不利于萃取過程的進(jìn)行，且高溫條件會(huì)造成油品收率的降低。

3.2.2 萃取時(shí)間對(duì)脫硫效果的影響 以質(zhì)量比1:1分別加入一定量的模型油和離子液體，在30℃下密封磁力攪拌10、20、30、40、50、60 min后，室溫下靜置，測(cè)定上下層中的三相組成，分析隨萃取時(shí)間變化，脫硫率、分配系數(shù)及選擇性系數(shù)的變化。實(shí)驗(yàn)結(jié)果見表5。

表5 萃取時(shí)間對(duì)萃取效果的影響

從表5可以看出，從10 min到40 min，隨著時(shí)間的延長(zhǎng)，脫硫率逐漸增加，40 min時(shí)脫硫率達(dá)到最大值，萃取40 min后脫硫率不再變化。這是因?yàn)殡S著時(shí)間的延長(zhǎng)，含硫化合物不斷地從油品中萃取到離子液體中，導(dǎo)致離子液體中的含硫化合物逐漸趨于飽和，達(dá)到萃取平衡，從而使得脫硫率不再變化。因此，萃取時(shí)間選擇40 min為較優(yōu)萃取時(shí)間。

3.2.3 劑油比對(duì)脫硫效果的影響 分別以質(zhì)量比4:1、2:1、1:1、1:2加入一定量的模型油和離子液體，在30℃下密閉磁力攪拌30 min后室溫下靜置，測(cè)定上下層中的三相組成，分析隨劑油比的變化，脫硫率、分配系數(shù)及選擇性系數(shù)的變化。實(shí)驗(yàn)結(jié)果見表6。

表6 劑油比對(duì)脫硫率和萃取效果的影響

由表6可知，隨著離子液體與模型油質(zhì)量比的增大，離子液體脫硫率逐漸增大，這是因?yàn)殡S著萃取劑用量的增加，苯并噻吩的溶解的量增大，脫硫率升高，因此增加萃取劑用量有利于脫硫率的提高。

3.2.4 萃取級(jí)數(shù)對(duì)脫硫效果的影響 以質(zhì)量比1:2加入一定量的模型油和離子液體，在30℃下密閉磁力攪拌30 min，結(jié)束后室溫靜置，待清晰分層后，測(cè)定其中的含硫量，再對(duì)其重復(fù)多次萃取，分析隨萃取次數(shù)增加，脫硫效果的變化。模型油中含硫量和脫硫率隨萃取級(jí)數(shù)的變化如圖2所示。

由圖2可看出隨著萃取級(jí)數(shù)的增加，脫硫率顯著增大，一級(jí)萃取脫硫率為70.72%，含硫量降至351.36 μg·g?1，經(jīng)過四級(jí)萃取之后，脫硫率達(dá)到99.27%，模型油含硫量降至8.87 μg·g?1，符合國(guó)Ⅴ標(biāo)準(zhǔn)。因此，離子液體可用于燃油深度脫硫，并具有較好的效果。

3.2.5 離子液體的回收利用 將脫硫后的離子液體置于在旋轉(zhuǎn)蒸發(fā)儀中，90℃下抽真空（真空度0.1 MPa）蒸餾60 min，得到回收離子液體。用回收的離子液體進(jìn)行萃取實(shí)驗(yàn)，測(cè)定回收離子液體的脫硫效果，如圖3所示。

由圖3可以看出回收離子液體的脫硫效果并無(wú)很大的變化，重復(fù)使用4次脫硫率仍可達(dá)到51.89%。雖然離子液體價(jià)格略高，但其重復(fù)使用性提高了其綜合利用價(jià)值，這也使得離子液體優(yōu)于其他脫硫溶劑。值得注意的是，在本實(shí)驗(yàn)中，雖然多級(jí)萃取能夠?qū)⒋蟛糠至蛎摮?，但同時(shí)會(huì)造成一小部分模型油的損失，因此，如何在深度脫硫的同時(shí)保證燃油的收率值得繼續(xù)研究。

3.3 正交實(shí)驗(yàn)設(shè)計(jì)

通過正交試驗(yàn)來(lái)優(yōu)化萃取脫硫條件，設(shè)計(jì)三因素三水平的正交試驗(yàn)，以脫硫率為考察指標(biāo)，正交實(shí)驗(yàn)因素和水平見表7，實(shí)驗(yàn)及優(yōu)化結(jié)果見表8。

表7 正交實(shí)驗(yàn)影響因素和水平

由表8可以看出，各因素對(duì)脫硫率的影響順序?yàn)椋篊(劑油比) > A(萃取溫度) > B(攪拌時(shí)間)，較優(yōu)操作條件為A1B3C1，而該操作條件并未出現(xiàn)在表8的9次試驗(yàn)中，按照得出的較優(yōu)方案進(jìn)行實(shí)驗(yàn)驗(yàn)證，在劑油比為2:1，萃取溫度20℃，密閉攪拌40 min條件下，得到該條件下離子液體單次脫硫率為72.74%。經(jīng)過四級(jí)脫硫?qū)嶒?yàn)之后，模型油含硫量降至6.98 μg·g?1，符合國(guó)Ⅴ標(biāo)準(zhǔn)。

表8 正交實(shí)驗(yàn)結(jié)果及分析

4 結(jié) 論

（1）通過離子液體脫硫機(jī)理和實(shí)驗(yàn)兩種方法，對(duì)基于COSMO-SAC模型分子設(shè)計(jì)方法的準(zhǔn)確性進(jìn)行驗(yàn)證，即：對(duì)[HMIM][BF4]、[HMIM][PF6]、[BMIM][BF4]、[BMIM][PF6]、[EMIM][BF4]、[EMIM][PF6] 6種離子液體作萃取劑時(shí)的脫硫效果進(jìn)行了脫硫機(jī)理的分析和實(shí)驗(yàn)的驗(yàn)證，得到的離子液體脫硫效果排序都與基于COSMO-SAC模型的分子設(shè)計(jì)結(jié)果基本一致，證明該分子設(shè)計(jì)方法具有一定的準(zhǔn)確性。

（2）[HMIM][PF6]作為萃取脫硫劑，通過液相色譜法測(cè)定了萃取脫硫后的液相組成，同時(shí)考察了萃取溫度、萃取時(shí)間、劑油比3個(gè)因素對(duì)脫硫率、分配系數(shù)及選擇性系數(shù)的影響。

（3）通過正交試驗(yàn)設(shè)計(jì)確定了萃取溫度20℃、萃取時(shí)間40 min、劑油比2:1為較優(yōu)操作條件，得到該條件下離子液體單次脫硫率為72.74%，四級(jí)脫硫可將模型油中的含硫量由1200 μg·g?1降至6.98 μg·g?1，符合國(guó)Ⅴ標(biāo)準(zhǔn)。

符 號(hào) 說 明

A——分子表面積，nm2 b——模糊綜合評(píng)判因子 F——隨機(jī)變量 f——自由度 K——分配系數(shù) m——質(zhì)量，g n——物質(zhì)的量，mol P(σ)——表面電荷密度分布函數(shù) q——表面積參數(shù) R——?dú)怏w常數(shù)，J·mol-1·K-1 r——體積參數(shù) SP——溶解度，g·g?1 T——體系溫度，K ΔW——交換能，J x——萃余相中組分的質(zhì)量分?jǐn)?shù) xi——組分i的摩爾分?jǐn)?shù) y——萃取相中組分的質(zhì)量分?jǐn)?shù) z——配位數(shù) β ——選擇性系數(shù) Г(σ)——表面電荷密度為σ的片段活度系數(shù) γ——活度系數(shù) θ——表面積分?jǐn)?shù) σ——屏蔽電荷密度，e·nm?2 φ——體積分?jǐn)?shù) 上角標(biāo) SG——Staverman-Guggenheim組合項(xiàng) ∞——無(wú)限稀釋 下角標(biāo) A——苯并噻吩 B——正辛烷 i——組分i S——萃取劑

References

[1] 崔文廣, 周二鵬, 趙地順, 等. 催化裂化汽油脫硫技術(shù)研究進(jìn)展[J]. 河北工業(yè)科技, 2004, 21(2): 1-4+12. CUI W G, ZHOU E P, ZHAO D S,. Technology progress of desulfurization of FCC gasoline[J]. Hebei Journal of Industral Science & Technology, 2004, 21(2): 1-4+12.

[2] 龐新梅, 孫書紅, 高雄厚. 生產(chǎn)低硫汽油新型FCC催化劑研究進(jìn)展[J]. 石化技術(shù)與應(yīng)用, 2001, 19(6): 384-388+397-350. PANG X M, SUN S H, GAO X H. New FCC catalysts for gasoline desulfurization[J]. Petrochemical Technology & Application, 2001, 19(6): 384-388+397-350.

[3] 王軍民, 袁鐵. 超低硫清潔汽油的生產(chǎn)技術(shù)進(jìn)展[J]. 天然氣與石油, 2001, 19(4): 14-17+1. WANG J M, YUAN T. Research progress in technology for super low sulfur gasoline production[J]. Natural Gas and Oil, 2001, 19(4): 14-17+1.

[4] 王鵬, 傅軍, 何鳴元. 含噻吩烷烴在分子篩上裂化脫硫的研究[J]. 石油煉制與化工, 2000, 31(3): 58-62. WANG P, FU J, HE M Y. Cracking-desulfurization of thiophene contained paraffin over zeolites[J]. Petroleum Processing and Petrochemicals, 2000, 31(3): 58-62.

[5] 寇曉撒, 史麗娟. 汽油萃取脫硫技術(shù)進(jìn)展[J]. 廣州化工, 2013, 41(3): 30-31. KOU X S, SHI L J. The progress of gasoline extraction desulfurization technology[J]. Guangzhou Chemical Industry, 2013, 41(3): 30-31.

[6] 曾丹林, 胡義, 王可苗, 等. 燃料油萃取脫硫技術(shù)研究進(jìn)展[J]. 石油煉制與化工, 2012, 43(5): 98-102. ZENG D L, HU Y, WANG K M,. Research advance in desulfurization of transportation fuels by extraction process[J]. Petroleum Processing and Petrochemicals, 2012, 43(5): 98-102.

[7] 張偉偉, 侯俠. 汽油萃取脫硫技術(shù)的應(yīng)用研究與進(jìn)展[J]. 蘭州石化職業(yè)技術(shù)學(xué)院學(xué)報(bào), 2015, 15(4): 5-9. ZHANG W W, HOU X. Progress on research and application of gasoline extraction desulphurization technology[J]. Journal of Lanzhou Petrochemical College of Technology, 2015, 15(4): 5-9.

[8] SWATLOSKI R P, HOLBREY J D, ROGERS R D. Ionic liquids are not always green: hydrolysis of 1-butyl-3-methylimidazolium hexafluorophosphate[J]. Green Chemistry, 2003, 5(4): 361-363.

[9] WASSERSCHEID P, HAL R V, B?SMANN A. 1--Butyl-3 -methylimidazolium ([BMIM]) octylsulfate — an even ‘greener’ ionic liquid[J]. Green Chemistry, 2002, 4(4): 400-404.

[10] 周瀚成, 陳楠, 石峰, 等. 離子液體萃取脫硫新工藝研究[J]. 分子催化, 2005, 19(2): 94-97. ZHOU H C, CHEN N, SHI F,. Studies on desulfurization of gasoline by extracting with ionic liquids[J]. Journal of Molecular Catalysis (China), 2005, 19(2): 94-97.

[11] 馮婕, 李春喜, 孟洪, 等. 磷酸酯類離子液體在燃油深度脫硫中的應(yīng)用[J]. 石油化工, 2006, 35(3): 272-276. FENG J, LI C X, MENG H,. Application of phosphate ionic liquids in deep desulfurization of fuel[J]. Petrochemical Technology, 2006, 35(3): 272-276.

[12] 王坤, 劉大凡, 何愛珍, 等. 離子液體萃取脫硫的研究[J]. 石油化工, 2010, 36(6): 675-680. WANG K, LIU D F, HE A Z,. Extractive desulfurization of simulative oil with ionic liquids[J]. Petrochemical Technology, 2010, 36(6): 675-680.

[13] KLAMT A. Conductor-like screening model for real solvents: a new approach to the quantitative calculation of solvation phenomena[J]. Journal of Physical Chemistry, 1995, 99(7): 2224-2235.

[14] KLAMT A, VOLKER J, THORSTEN A B,. Refinement and parametrization of COSMO-RS[J]. Journal of Physical Chemistry A, 1998, 102(26): 5074-5085.

[15] KLAMT A, ECKERT F. COSMO-RS: a novel and efficient method for the a priori prediction of thermophysical data of liquids[J]. Fluid Phase Equilibria, 2000, 172(172): 43-72.

[16] KLAMT A, KROOSHOF G J P, TAYLOR R. COSMOSPACE: alternative to conventional activity-coefficient models[J]. AIChE Journal, 2002, 48(10): 2332-2349.

[17] LIN S T, SANDLER S I. A priori phase equilibrium prediction from a segment contribution solvation model[J]. Industrial & Engineering Chemistry Research, 2002, 41(5): 899-913.

[18] LIN S T, CHANG J, WANG S,. Prediction of vapor pressures and enthalpies of vaporization using a COSMO solvation mode[J]. Journal of Physical Chemistry A, 2004, 108(36): 7429-7439.

[19] HSIEH C M, SANDLER S I, LIN S T. Improvements of COSMO-SAC for vapor-liquid and liquid-liquid equilibrium predictions[J]. Fluid Phase Equilibria, 2010, 297(1): 90-97.

[20] 吳選軍, 余永富, 袁繼祖. COSMO結(jié)合狀態(tài)方程預(yù)測(cè)汽液平衡及可視化[J]. 武漢理工大學(xué)學(xué)報(bào)(信息與管理工程版), 2011, 33(3): 354-358. WU X J, YU Y F, YUAN J Z. Prediction and visualization of vapor -liquid equilibria based on COSMO model and state equation[J]. Journal of WUT (Information & Management Engineering), 2011, 33(3): 354-358.

[21] 李瑞, 崔現(xiàn)寶, 吳添, 等. 基于COSMO-SAC模型的離子液體萃取劑的選擇[J]. 化工學(xué)報(bào), 2013, 64(2): 452-469.LI R, CUI X B, WU T,. Selection of ionic liquid solvent for liquid-liquid extraction based on COSMO-SAC model[J]. CIESC Journal, 2013, 64(2): 452-469.

[22] MULLINS E, OLDLAND R, LIU Y A,. Sigma-profile database for using COSMO-based thermodynamic methods[J]. Industrial & Engineering Chemistry Research, 2006, 45(12): 4389-4415.

[23] PRETEL E J, LOPEZ P A, BOTTINI S B,. Computer-aided molecular design of solvents for separation process[J]. AIChE Journal, 1994, 48(8): 1349-1360.

[24] 王秀麗, 孫偉, 宋海華. 萃取精餾溶劑的模糊綜合評(píng)判[J]. 天津大學(xué)學(xué)報(bào), 2003, 36(4): 478-481.WANG X L, SUN W, SONG H H. Fuzzy comprehensive judgement on extractive distillation solvents[J]. Journal of Tianjin University, 2003, 36(4): 478-481.

[25] ZHANG S G, ZHANG Z C. Novel properties ofionic liquids in selective sulfur removal from fuels atroom temperature[J]. Green Chemistry, 2002, 4(4): 376-379.

[26] SU B M, ZHANG S G, ZHANG Z C. Structuralelucidation of thiophene interaction with ionic liquidsby multinuclear NMR spectroscopy[J]. Journal ofPhysical Chemistry B, 2004, 108(50): 19510-19517.

[27] MOCHIZUKI Y, SUGAWARA K. Removal of organic sulfur fromhydrocarbon resources using ionic liquids[J]. Energy & Fuels, 2008, 22(5): 3303-3307.

[28] NIE Y, LI C X, HONG M,.,-Dialkylimidazoliumdialkylphosphate ionic liquids: their extractive performance forthiophene series compounds from fuel oils versus the length of alkylgroup[J]. Full Processing Technology, 2008, 89(10): 978-983.

[29] 黃廷昌, 苑建. 離子液體應(yīng)用于燃油脫硫技術(shù)的研究進(jìn)展[J]. 河北化工, 2011, 34(5): 44-47+67.HUANG T C, YUAN J. Development on desulfurization of fuel oil by ionic liquids[J]. Hebei Huagong, 2011, 34(5): 44-47+67.

[30] 胡松青, 張軍, 劉冰, 等. 離子液體萃取脫硫的探索性研究[J]. 石油學(xué)報(bào)(石油加工), 2007, 23(1): 100-103. HU Q S, ZHANG J, LIU B,. Feasibility study of extraction desulfurization with ionic liquids[J]. ACTA Petroleisinica (Petroleum Processing Section), 2007, 23(1): 100-103.

[31] 林燕, 王芳, 張志慶, 等. 離子液體綠色脫硫機(jī)理及應(yīng)用進(jìn)展[J]. 化工進(jìn)展, 2013, 32(3): 549-557. LIN Y, WANG F, ZHANG Z Q,. Mechanism and application of ionic liquids in environmental friendly oil desulphurization[J]. Chemical Industry and Engineering Progress, 2013, 32(3): 549-557.

Selection and process optimization of ionic liquids for desulfurization

FANG Jing, ZHANG Shuting, LI Tingting, LI Chunli

(School of Chemical Engineering, Hebei University of Technology, Tianjin 300130, China)

To verify the reliability of ionic liquids’ molecular design method based on COSMO-SAC model, ionic liquids desulfurization mechanism analysis and experimental results are compared with COSMO simulation.By analyzing [HMIM][BF4], [HMIM][PF6], [BMIM][BF4], [BMIM][PF6], [EMIM][BF4] and [EMIM][PF6] on mechanism level, [HMIM][PF6] was considered as a better extractive desulfurization agent than others, and the order of desulfurization efficiency was consistent with the calculation results from COSMO-SAC model. The experimental result shows that the difference between the COSMO simulation and experiments is not significant. In this study, [HMIM][PF6] was used as an extractive desulfurization agent to remove benzothiophene from fuels, and the change of sulfur content in model oil was determined by liquid chromatography. Through considering the effect of extraction temperature, time and the ratio of solvent to oil, the optimal operating conditions, 20℃, 40 min and 2:1 respectively, were determined in orthogonal experiment with an efficiency over 72.74%. With [HMIM][PF6] used as the deep extraction desulfurization solvent, the sulfur content decreased from 1200 μg·g?1to 6.98 μg·g?1by 4-stages extraction, which met the Phase Ⅴ Emission Standards in China.

ionic liquids; extraction; desulfurization; simulation; orthogonal experiment

10.11949/j.issn.0438-1157.20170182

TQ 028.3

A

0438—1157（2017）09—3434—08

2017-02-24收到初稿，2017-05-31收到修改稿。

李春利。

方靜（1977—），女，博士，教授。

教育部博士點(diǎn)基金項(xiàng)目（20131317120014）。

2017-02-24.

Prof. LI Chunli, ctstfj@163.com

supported by the Ph.D Program Foundation of Ministry of Education of China (20131317120014).