• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    萃取精餾復(fù)合溶劑的復(fù)配選擇

    2017-11-22 10:20:23方靜軒碧涵李曉春李春利
    化工學(xué)報(bào) 2017年11期
    關(guān)鍵詞:環(huán)己烷苯胺乙酸乙酯

    方靜,軒碧涵,李曉春,李春利

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

    萃取精餾復(fù)合溶劑的復(fù)配選擇

    方靜,軒碧涵,李曉春,李春利

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

    提出按溶劑之間復(fù)配形成氫鍵的情況和修正的UNIFAC模型相結(jié)合的方法選擇分離共沸體系的萃取精餾復(fù)合溶劑。準(zhǔn)確測量了共沸體系與萃取精餾單溶劑和復(fù)合溶劑的汽液平衡數(shù)據(jù),結(jié)果表明采用的萃取精餾單溶劑和復(fù)合溶劑的選擇方法具有較好的準(zhǔn)確度,模擬值與實(shí)驗(yàn)值比較,相對偏差均小于 9%。通過綜合比較得到分離乙酸乙酯-乙醇、乙醇-水和環(huán)己烷-苯3個(gè)共沸體系的最佳復(fù)合溶劑及配比。在最佳復(fù)合溶劑組成下,原共沸體系的相對揮發(fā)度達(dá)到最大值,該值優(yōu)于組成復(fù)合溶劑的單溶劑所達(dá)到的分離效果。在最佳復(fù)合溶劑組成下,改變原料體系組成,復(fù)合溶劑組成對原料體系相對揮發(fā)度的影響趨勢基本相同,這表明復(fù)合溶劑在萃取精餾塔內(nèi)不同原料組成下均能較好地增加體系的相對揮發(fā)度,為工業(yè)上在萃取精餾過程中使用復(fù)合溶劑提供了理論依據(jù)。

    萃取精餾;復(fù)合溶劑;共沸(混合)物;相平衡;分離

    引 言

    決定萃取精餾過程分離效率高低的關(guān)鍵是溶劑的選擇合理與否。用于萃取精餾的溶劑從組成上分為單一溶劑和復(fù)合溶劑兩大類,在實(shí)踐應(yīng)用過程中,有些單一溶劑往往存在自身選擇性與溶解性的矛盾,或是在萃取精餾塔內(nèi)溶劑比過大造成生產(chǎn)過程中的能耗大,同時(shí)溶劑回收時(shí)的能耗也相應(yīng)增加。

    單一溶劑的這些弱點(diǎn)限制了萃取精餾技術(shù)的應(yīng)用及發(fā)展,于是人們開始了使用混合溶劑進(jìn)行萃取精餾的研究,并取得了良好的效果?;旌先軇┚褪窃趩我蝗軇┗A(chǔ)上加入一種或一種以上的組分構(gòu)成新的溶劑,加入的組分稱為副溶劑[1],其主要作用就是增大組分間揮發(fā)度的差異,從而降低分離的難度,一方面可以改善單一溶劑的選擇性,另一方面可以避免當(dāng)溶劑比過大時(shí)溶劑與待分離體系形成不完全互溶的兩個(gè)液相,既能保證溶劑在萃取精餾塔內(nèi)保持較高的濃度,同時(shí)又使整個(gè)萃取精餾混合物液相處于均相狀態(tài),有利于傳質(zhì)和傳熱。目前,關(guān)于萃取精餾溶劑國內(nèi)外的研究進(jìn)展主要針對離子液體、加鹽溶劑和復(fù)合溶劑3類。

    Orchillés 等[1-2]測定了離子液體 1-乙基-3-甲基咪唑二氰胺([EMIM][DCA])和1-乙基-3-甲基咪唑三氟甲磺酸鹽([EMIM][CF3SO3])存在時(shí)乙醇-水物系的常壓汽液平衡數(shù)據(jù),發(fā)現(xiàn)了共沸現(xiàn)象消失。劉艷升等[3]以咪唑?yàn)殛栯x子,以氟磷酸或鹵素為陰離子,合成了1種離子液體或者2種以上離子液體的復(fù)合溶劑用于分離苯和環(huán)己烷。Domańska等[4-8]測定了烷烴、環(huán)烷烴、芳香烴等多種溶質(zhì)在離子液體 4-甲基-N-丁基吡啶([bmpy][Tf2N])、1-丁基 3-甲基吡啶硫氰酸鹽([BMPY][SCN])、1-丙基-1-甲基吡啶雙三氟甲磺酰亞胺鹽([PMPIP][Tf2N])、1-丁基-1-甲基吡啶硫氰酸鹽([BMPIP][SCN])、1-丁基-1-甲基吡咯烷酮硫氰酸鹽([BMPYR][SCN])、1-丁基-1-甲基吡咯烷酮三氟甲磺酸鹽([BMPYR][CF3SO3])中的無限稀釋活度系數(shù),計(jì)算出離子液體對苯-正庚烷和苯-環(huán)己烷物系的選擇性系數(shù),以及對苯的容量,并與傳統(tǒng)的N-甲基吡咯烷酮(NMP)和環(huán)丁砜溶劑進(jìn)行了比較。發(fā)現(xiàn)離子液體對苯-正庚烷和苯-環(huán)己烷的選擇性系數(shù)高于NMP和環(huán)丁砜的選擇性系數(shù),離子液體對苯的容量也高于環(huán)丁砜溶劑,表明離子液體可替代傳統(tǒng)溶劑用于這兩種物系的分離。Jongmans等[9-10]對離子液體分離乙苯-苯乙烯的過程進(jìn)行了研究,考察了離子液體結(jié)構(gòu)對其性能的影響。發(fā)現(xiàn)由烷基鏈較短的陽離子和不存在電子離域的陰離子構(gòu)成的離子液體的選擇性較高;由烷基鏈較長的陽離子和存在電子離域的陰離子構(gòu)成的離子液體的容量較大。國內(nèi)外在離子液體萃取精餾方面做了大量工作,也取得了一定的進(jìn)展。但離子液體存在價(jià)格昂貴、制備過程復(fù)雜、動(dòng)力學(xué)和熱力學(xué)數(shù)據(jù)缺乏等問題,其回收再生的報(bào)道較少,目前僅限于實(shí)驗(yàn)室研究[11]。

    針對N,N-二甲基甲酰胺(DMF)溶劑萃取精餾分離 C4或 C5時(shí)溶劑用量大、能耗高的現(xiàn)狀,Lei等[12]測定了DMF-C4物系和加鹽DMF-C4物系的汽液平衡數(shù)據(jù)。實(shí)驗(yàn)結(jié)果顯示,鹽的加入使物系的相對揮發(fā)度提高了1.4倍。李東風(fēng)等[13]對加鹽萃取精餾分離C5的過程進(jìn)行了研究,發(fā)現(xiàn)NaSCN的存在能提高物系的相對揮發(fā)度至 1.889,減少溶劑用量65%。為了脫除粗苯中的正庚烷等雜質(zhì),Dong等[14]采用加鹽DMF溶劑,對苯-正庚烷的分離進(jìn)行了研究,發(fā)現(xiàn)KSCN的加入使物系的相對揮發(fā)度的提高最大,達(dá)到3以上??傮w來看,加鹽溶劑可以改善物系的相對揮發(fā)度,減少溶劑用量,降低生產(chǎn)能耗,但在含鹽的非電解質(zhì)混合溶液中,鹽效應(yīng)情況十分復(fù)雜,理論研究還不夠深入,仍需做進(jìn)一步的工作[11]。

    一些研究者采用復(fù)合溶劑作為萃取劑在一些萃取精餾過程中取得了良好的效果。Borwn等[15-17]通過對 2,3-二甲基戊烷和環(huán)己烷、環(huán)戊烷和 2,2-二甲基丁烷、苯乙烯和鄰二甲苯分離體系的系列實(shí)驗(yàn)研究。研究表明混合溶劑比單一溶劑在相同條件下有更高的選擇性,混合溶劑的選擇性不僅與主溶劑密切相關(guān)而且還與輔助溶劑有著直接的關(guān)系?;旌先軇┧瑔我惠腿┎煌瑫r(shí),其選擇性差異明顯,當(dāng)混合溶劑含有相同的單一萃取劑但其配比不同時(shí),其選擇性也有很大的差異。新型混合溶劑提高了分離過程的選擇性、降低了理論板數(shù),在實(shí)際的工業(yè)領(lǐng)域也有較好的應(yīng)用。

    祝石華等[18]采用 N-甲酰嗎啉(NFM)添加一種助溶劑(COS)構(gòu)成復(fù)合溶劑作萃取精餾溶劑,從裂解加氫汽油苯餾分中回收高純苯。研究了助溶劑含量、溶劑比等因素對萃取精餾效果的影響,適宜的助溶劑質(zhì)量分?jǐn)?shù)為5%~20%,溶劑比為4.5~6.0。采用復(fù)合溶劑,苯產(chǎn)品收率達(dá)到 99.8%,溶劑回收溫度降至190℃以下,與NFM單溶劑萃取精餾工藝相比,具有產(chǎn)品質(zhì)量好、苯收率高、操作苛刻度低的優(yōu)點(diǎn)。

    曹新波等[19]利用N-甲酰嗎琳(NFM)和甲乙酮(MEK)作為混合溶劑萃取精餾分離正丁烯與丁烷,并用模擬計(jì)算考察了溶劑比、回流比、溶劑組成等參數(shù)對NFM-MEK混合溶劑萃取精餾提濃正丁烯過程的影響。結(jié)果表明合理溶劑比、溶劑組成的選擇應(yīng)基于經(jīng)濟(jì)權(quán)衡,而回流比則存在一個(gè)最佳值,所得出的結(jié)論對工業(yè)生產(chǎn)具有一定的指導(dǎo)意義。

    華超等[20]利用改進(jìn)的 UNIFAC模型來選擇萃取精餾苯-甲基環(huán)戊烷過程中的混合萃取劑,對混合溶劑N-甲酰嗎啉(NFM)-乙二醇(EG)對萃取精餾苯-甲基環(huán)戊烷的過程進(jìn)行了詳細(xì)的理論分析和實(shí)驗(yàn)研究,重點(diǎn)考察了混合溶劑配比和溶劑比對分離過程可行性的影響。研究結(jié)果表明,混合溶劑NFM-EG具有良好的選擇性和溶解性,為粗苯精制過程中苯-甲基環(huán)戊烷的分離提供了必要的基礎(chǔ)數(shù)據(jù)和指導(dǎo)。Yin等[21]開展了類似的研究工作,發(fā)現(xiàn)助溶劑EG的質(zhì)量分?jǐn)?shù)為20%時(shí),苯-環(huán)己烷物系的分離效果最佳。張志剛等[22]采用N,N-二甲基甲酰胺和二甲亞砜作為復(fù)合溶劑萃取精餾分離苯-環(huán)己烷體系,考察了不同萃取劑在不同回流比及萃取溶劑加入速率情況下對分離苯-環(huán)己烷共沸體系的影響。結(jié)果表明,二元混合溶劑能夠解決單一溶劑的選擇性與溶解性相矛盾的問題,且在同等條件下,綜合性能優(yōu)于單一溶劑。

    與離子液體和加鹽溶劑相比,復(fù)合溶劑具有產(chǎn)品純度和收率高、使用已有的工藝設(shè)備、操作條件簡單等優(yōu)點(diǎn),因此目前復(fù)合溶劑更具實(shí)用價(jià)值,具有極好的開發(fā)應(yīng)用前景。

    為了提高萃取精餾溶劑的分離性能,減少精餾塔內(nèi)的溶劑比,降低精餾能耗,提出了按溶劑間形成氫鍵的情況和修正UNIFAC模型相結(jié)合的萃取精餾復(fù)合溶劑的復(fù)配方法,選擇出分離性能高于單溶劑的復(fù)合溶劑萃取精餾分離共沸體系,對復(fù)合溶劑的復(fù)配理論、選擇方法、配比原則、復(fù)合溶劑組成對原料組分相對揮發(fā)度的影響等進(jìn)行探索性研究,并以乙酸乙酯-乙醇、乙醇-水和環(huán)己烷-苯3種工業(yè)上典型的物系為研究對象,建立汽液平衡測量裝置,準(zhǔn)確測量汽液平衡數(shù)據(jù),將實(shí)驗(yàn)值與模擬值進(jìn)行比較并進(jìn)行結(jié)果討論,驗(yàn)證這種方法的可靠性與實(shí)用性,為工業(yè)上在萃取精餾過程中使用復(fù)合溶劑提供了理論依據(jù)。

    1 萃取精餾溶劑的復(fù)配

    1.1 復(fù)配理論研究

    在萃取精餾過程中,溶劑與待分離體系組分之間的作用力主要是較弱的范德華力和較強(qiáng)的氫鍵力。范德華力包括感生力、取向力及色散力,但它們相對于氫鍵力而言都是很小的。從生成機(jī)理上說,氫鍵是由一個(gè)氫原子和兩個(gè)電負(fù)性原子結(jié)合構(gòu)成的,如C—H…D,此處C和D是電負(fù)性原子,“…”表示氫鍵,氫鍵是一種方向性很強(qiáng)的作用力。溶劑之間形成氫鍵的必要條件是:第一,要有電子接受體即質(zhì)子給體(proton donor);第二,要有電子給予體即質(zhì)子受體[23](proton acceptor)。按形成氫鍵的情況,可將溶劑分為以下4種類型。

    (1)N型溶劑:它們不能形成氫鍵,分子電負(fù)性非常小。如苯、環(huán)己烷、四氯化碳等。

    (2)A型溶劑:即受電子溶劑。如鹵代烴等,可形成氫鍵。

    (3)B型溶劑:即給電子溶劑。如酮、醛、醚、酯等,可形成氫鍵。

    (4)A BⅠ型溶劑:給受電子溶劑。如水、醇、酸、胺等,可形成氫鍵。

    (5)A BⅡ型溶劑:生成內(nèi)氫鍵的給受電子溶劑。如鄰硝基苯酚等,已形成內(nèi)氫鍵,故性質(zhì)同N型或B型溶劑相似。

    根據(jù)上述分類將常用的溶劑分為質(zhì)子給體和質(zhì)子受體,當(dāng)兩種溶劑復(fù)配時(shí)形成的氫鍵比混合以前更強(qiáng)或是數(shù)目更多的時(shí)候則利于互溶,對被分離的共沸體系的分離效果也將更好。

    1.2 混合物熱力學(xué)性質(zhì)的估算

    基于熱力學(xué)模型,建立定性定量選擇萃取劑的的理論方法,已報(bào)道了Prausnitz & Anderson理論[24]、正規(guī)溶液理論Scatchard-Hildebrnad法[25]、無限稀釋活度系數(shù)法[26]、官能團(tuán)UNIFAC法[27]、官能團(tuán)ASOG法[28]、溶解度參數(shù)MOSCED法[29]等多種利用各類熱力學(xué)模型計(jì)算的方法。由于萃取精餾過程中存在的物系大都為非理想性極強(qiáng)的體系,所以大多數(shù)不以官能團(tuán)作為基礎(chǔ)模型參數(shù)的熱力學(xué)方法在選擇萃取溶劑的過程中,局限性很強(qiáng)。因此,這些方法發(fā)展速度較慢或不再進(jìn)行。對于有機(jī)溶劑而言,它們雖然種類繁多,但基本是由幾十種基團(tuán)所構(gòu)成。目前,通過應(yīng)用官能團(tuán)相互作用參數(shù)值預(yù)測組分和混合物比較流行的方法是UNIFAC法。

    本文采用Gmehling等[30-32]修正的UNIFAC模型來計(jì)算組分的活度系數(shù)及無限稀釋活度系數(shù),求組分i在組分j中的無限稀釋活度系數(shù)時(shí),模型中的公式可以簡化為如下形式

    利用這些簡化后的公式,可以比較便捷地求得體系中組分的無限稀釋活度系數(shù),為以后分離過程溶劑性能的計(jì)算提供了便利。

    1.3 萃取精餾溶劑復(fù)配的流程

    萃取精餾溶劑的復(fù)配流程如圖1所示。

    圖1 萃取精餾復(fù)合溶劑復(fù)配流程Fig.1 Combination of mixed solvent in extractive distillation

    2 復(fù)合溶劑的應(yīng)用實(shí)例

    2.1 乙酸乙酯-乙醇體系

    從乙酸乙酯和乙醇分子極性來看,這兩種分子均具有極性,乙酸乙酯分子為質(zhì)子受體,易與質(zhì)子給體形成氫鍵,乙醇分子既是質(zhì)子給體又是質(zhì)子受體,易與質(zhì)子給體和受體物質(zhì)形成氫鍵。乙酸乙酯和乙醇的極性及第一電離能見表1。分子間由極性引起的誘導(dǎo)力、靜電力及色散力見式(13)~式(15)。由于乙酸乙酯和乙醇的偶極矩、第一電離能基本較為接近,誘導(dǎo)作用貢獻(xiàn)一般都很小,從表1及式(13)~式(15)來看,分子之間產(chǎn)生靜電力及色散力基本相等,故靜電力及色散力影響比較小。

    表1 乙酸乙酯、乙醇的偶極矩及第一電離能[33]Table 1 Dipole moment and first ionization energy of ethyl acetate and ethanol

    式中,εp1、εp2、εp3分別為分子間的誘導(dǎo)力、靜電力、色散力;αi、αj為分子的極化率;μi、μj為分子的偶極矩;r為分子間距離;K為Boltzmann因子;Ii、Ij為第一電離能。

    分離乙酸乙酯-乙醇體系的常用溶劑列于表2中。

    2.1.1 不同復(fù)合溶劑對乙酸乙酯-乙醇體系相對揮發(fā)度的影響 由表2可以看出只有二甲亞砜(DMSO)與乙苯為質(zhì)子受體,其他溶劑都是質(zhì)子給體與受體同體的溶劑,其中多元醇類溶劑,由于黏度較大,易于聚合,在實(shí)際生產(chǎn)中應(yīng)用困難,因此只將乙二醇(EG)作為醇類溶劑用作復(fù)合溶劑的復(fù)配。同時(shí)由于乙二醇和乙苯能形成共沸物不能復(fù)配組合,N-甲基吡咯烷酮(NMP)和二甲基亞砜(DMSO)缺少交互作用參數(shù),無法用UNIFAC方程預(yù)測復(fù)配結(jié)果。將余下的溶劑用修正的UNIFAC方程經(jīng)Matlab軟件編程計(jì)算,有復(fù)合溶劑效應(yīng)的有以下復(fù)配組合:

    A.DMSO(a)+ DMF(b)

    B.DMSO(a)+ EG(b)

    C.NMP(a)+乙苯(b)

    D.乙苯(a)+ DMF(b)

    并將 A、B、C、D 4種復(fù)合溶劑與共沸組成(xEA=0.69)下的乙酸乙酯-乙醇體系在溶劑質(zhì)量比為1:1的條件下,經(jīng)修正的UNIFAC方程計(jì)算得到的隨復(fù)合溶劑中 a組分的增加乙酸乙酯-乙醇相對揮發(fā)度的變化曲線繪制成圖2。

    表2 萃取精餾分離乙酸乙酯-乙醇體系單溶劑備選Table 2 Single-solvent of separating ethyl acetate-ethanol system by extractive distillation

    圖2 溶劑復(fù)配對乙酸乙酯-乙醇的相對揮發(fā)度的影響Fig.2 Influence of mixed solvent content for ethyl acetate-ethanol on relative volatility

    從圖中可以看出這幾組溶劑復(fù)配時(shí)對乙酸乙酯-乙醇共沸體系的相對揮發(fā)度均有影響,并且都存在一個(gè)最佳的溶劑配比。各組復(fù)合溶劑對體系的影響不相同,從圖中可以看出DMSO和DMF復(fù)配時(shí)對乙酸乙酯-乙醇體系的影響最大。溶劑復(fù)配的最佳值見表3。

    表3 萃取精餾分離乙酸乙酯-乙醇體系復(fù)合溶劑備選Table 3 Mixed-solvent of separating ethyl acetate-ethanol system by extractive distillation

    從表中可以看出60%DMSO和40%DMF復(fù)配時(shí)得到的乙酸乙酯和乙醇的相對揮發(fā)度最大為3.0266,大于單溶劑分離乙酸乙酯-乙醇共沸體系的最大相對揮發(fā)度。

    2.1.2 不同組成下的復(fù)合溶劑(DMSO+DMF)對乙酸乙酯-乙醇體系相對揮發(fā)度的影響 在溶劑比(質(zhì)量比)為1:1時(shí),改變復(fù)合溶劑DMSO和DMF的組成,并改變乙酸乙酯-乙醇體系中乙酸乙酯的質(zhì)量分?jǐn)?shù),模擬計(jì)算乙酸乙酯-乙醇的相對揮發(fā)度(圖3)。

    圖3 溶劑組成(DMSO+DMF)對乙酸乙酯-乙醇相對揮發(fā)度的影響Fig.3 Influence of DMSO+DMF ratio in binary mixed solvent for ethanol-ethyl acetate on relative volatility

    從圖3可以看出對于乙酸乙酯-乙醇體系,改變體系中乙酸乙酯的質(zhì)量分?jǐn)?shù)時(shí)復(fù)合溶劑組成對體系相對揮發(fā)度的影響趨勢基本相同,即在不同的原料組成下(xEA=0.3,xEA=0.5,xEA=0.69)復(fù)合溶劑的最佳組成都在xDMSO=0.6(質(zhì)量分?jǐn)?shù)),這說明如果在萃取精餾操作的過程中,即使在萃取精餾段中乙酸乙酯的組成隨塔高變化,使用適宜配比的復(fù)合溶劑,在一定的溶劑比下,都可以保持溶劑對乙酸乙酯-乙醇體系的高分離效果。

    表4 萃取精餾分離乙醇-水體系單溶劑備選Table 4 Single-solvent of separating ethanol-water system by extractive distillation

    2.1.3 不同溶劑比對乙酸乙酯-乙醇體系相對揮發(fā)度的影響 在單溶劑(DMSO、DMF)、復(fù)合溶劑(DMSO+DMF)的最佳組合下(xDMSO=0.6)改變?nèi)軇┍龋玫揭宜嵋阴?乙醇體系中(xEA=0.69)相對揮發(fā)度隨溶劑比的變化趨勢,結(jié)果如圖4所示。

    圖4 不同溶劑比對乙酸乙酯-乙醇體系相對揮發(fā)度的影響Fig.4 Influence of solvent ratio for ethanol-ethyl acetate on relative volatility

    從圖4中可以看出隨著溶劑比的增加,不管是單一溶劑還是復(fù)合溶劑,體系的相對揮發(fā)度都逐漸增加,在溶劑比較小時(shí)(小于 4)體系的相對揮發(fā)度增加很快,隨著溶劑比的逐漸增加,體系的相對揮發(fā)度增加速度變化緩慢。這主要是因?yàn)楫?dāng)溶劑比較小時(shí),隨著溶劑比的增加,溶劑改變乙酸乙酯-乙醇相對揮發(fā)度的作用增強(qiáng)效果較為明顯。但是,當(dāng)溶劑比較大時(shí),增加溶劑比,也會(huì)使溶劑對原料組分的稀釋作用增強(qiáng),因此相對揮發(fā)度的增加幅度減小。

    2.2 乙醇-水體系

    乙醇分子和水分子都是質(zhì)子給體和質(zhì)子受體同體的分子,分離乙醇-水體系的常用溶劑列于表4中。

    2.2.1 不同復(fù)合溶劑對乙醇-水體系相對揮發(fā)度的影響 由表4可以看出只有二甲亞砜(DMSO)為質(zhì)子受體,其他溶劑都是質(zhì)子給體與受體同體的溶劑,其中多元醇類溶劑,由于黏度較大,易于聚合,在實(shí)際生產(chǎn)中應(yīng)用困難,因此只將乙二醇作為醇類溶劑用作復(fù)合溶劑的復(fù)配。經(jīng)修正的UNIFAC方程的計(jì)算,得到以下復(fù)配結(jié)果:

    A.DMSO(a)+DMF(b)

    B.DMSO(a)+乙二醇(b)

    C.DMSO(a)+N,N-二甲基乙酰胺(DMAC)(b)

    并將 A、B、C 3種復(fù)合溶劑與共沸組成(xEOH=0.95)下的乙醇-水體系在溶劑質(zhì)量比為1:1的條件下,經(jīng)修正的UNIFAC方程計(jì)算得到的隨復(fù)合溶劑中 a組分的增加乙醇-水相對揮發(fā)度的變化曲線繪制成圖5。

    圖5 溶劑復(fù)配對乙醇-水體系的相對揮發(fā)度的影響Fig.5 Influence of mixed solvent content for ethanol-water on relative volatility

    從圖5中可以看出這幾組溶劑復(fù)配時(shí)對乙醇-水共沸體系的相對揮發(fā)度的影響均有一組最佳的溶劑配比,但是各個(gè)組合對體系的影響不同,其中復(fù)合溶劑 A(DMSO+DMF)與復(fù)合溶劑 C(DMSO+DMAC)對乙醇-水共沸體系相對揮發(fā)度的影響變化關(guān)系相似,而復(fù)合溶劑A在最佳復(fù)配組合下得到較高的相對揮發(fā)度1.65,高于單溶劑的最高相對揮發(fā)度。當(dāng)二甲基亞砜與乙二醇復(fù)配時(shí),沒有出現(xiàn)明顯的復(fù)配效果,隨著DMSO在復(fù)合溶劑中質(zhì)量分?jǐn)?shù)的增加,乙醇-水體系的相對揮發(fā)度大體上呈遞減的變化趨勢。對于A、B、C 3種復(fù)配溶劑的最佳配比見表5。

    表5 萃取精餾分離乙醇-水體系復(fù)合溶劑備選Table 5 Mixed-solvent of separating ethanol-water system by extractive distillation

    從表5中可以看出溶劑70%DMSO+30%DMF復(fù)配時(shí)乙醇-水的相對揮發(fā)度為 1.6504為最佳的溶劑復(fù)配結(jié)果。

    2.2.2 不同組成下的復(fù)合溶劑(DMSO+DMF)對乙醇-水體系相對揮發(fā)度的影響 在溶劑質(zhì)量比為1:1時(shí),改變復(fù)合溶劑DMSO和DMF的組成并改變乙醇-水體系中乙醇的質(zhì)量分?jǐn)?shù),所得的乙醇-水的相對揮發(fā)度結(jié)果如圖6所示。

    圖6 溶劑組成(DMSO+DMF)對乙醇-水相對揮發(fā)度的影響Fig.6 Influence of DMSO+DMF ratio in binary mixed solvent for ethanol-water on relative volatility

    從圖6中可以看出對于乙醇-水體系改變乙醇的質(zhì)量分?jǐn)?shù)時(shí),復(fù)合溶劑(DMSO+DMF)對體系相對揮發(fā)度的改變趨勢相似,均存一最佳配比,即在DMSO在復(fù)合溶劑中質(zhì)量分?jǐn)?shù)為0.7左右時(shí),體系的相對揮發(fā)度達(dá)到最大值。這表明如果在萃取精餾操作的過程中,即使在萃取精餾段中乙醇的組成隨塔高變化,使用適宜配比的復(fù)合溶劑,在一定的溶劑比下,都可以保持溶劑對乙醇-水體系的高分離效果。

    2.2.3 不同溶劑比對乙醇-水體系相對揮發(fā)度的影響 在單溶劑(DMSO、DMF)、復(fù)合溶劑(DMSO+DMF)的最佳組合下(xDMSO=0.7)改變?nèi)軇┍?,得到乙?水共沸體系中相對揮發(fā)度隨溶劑比的變化趨勢,結(jié)果如圖7所示。

    圖7 不同溶劑比對乙醇-水體系相對揮發(fā)度的影響Fig.7 Influence of solvent ratio on relative volatility of ethanol-ethyl acetate

    從圖7中可以看出隨著溶劑比的增加,不管是單一溶劑還是復(fù)合溶劑,體系的相對揮發(fā)度都逐漸增加,在溶劑比較小時(shí)(<3)體系的相對揮發(fā)度增加很快,隨著溶劑比的逐漸增加體系的相對揮發(fā)度增加速度變化緩慢。

    2.3 環(huán)己烷-苯

    環(huán)己烷分子是飽和烴,為非極性化合物,既不是質(zhì)子給體又不是質(zhì)子受體。苯分子是具有π電子系統(tǒng)的質(zhì)子受體。分離環(huán)己烷-苯的常用溶劑列于表6中。

    2.3.1 不同復(fù)合溶劑對環(huán)己烷-苯體系相對揮發(fā)度的影響 由表6可以看出吡啶、苯乙酮、DMSO為質(zhì)子受體,苯胺、DMF是質(zhì)子給體與受體同體的溶劑,應(yīng)將這兩種溶劑進(jìn)行組合得到復(fù)配溶劑。其中DMF作為溶劑,雖然其選擇性相對較高,但是它與環(huán)己烷的互溶度較低,溶劑出現(xiàn)分層現(xiàn)象,影響萃取精餾的效果。2003年張志剛[34]在對苯-環(huán)己烷萃取精餾復(fù)配溶劑對研究時(shí)也發(fā)現(xiàn),當(dāng)原料中環(huán)己烷含量稍大時(shí),由于混合液中存在 DMF出現(xiàn)分層現(xiàn)象,萃取精餾效果受到影響,因此復(fù)配時(shí)不考慮DMF。同時(shí)由于苯胺-二甲基亞砜均缺少交互作用參數(shù),將余下的溶劑進(jìn)行復(fù)配,經(jīng)過修正的UNIFAC方程的計(jì)算,得到以下復(fù)配結(jié)果:

    表6 萃取精餾分離環(huán)己烷-苯體系單溶劑備選Table 6 Single-solvent of separating ethyl cyclohexane-benzene system by extractive distillation

    A.苯胺(a)+ 苯乙酮(b)

    B.苯胺(a)+吡啶(b)

    將A、B兩種復(fù)合溶劑與共沸組成(xbex=0.55)下的環(huán)己烷-苯體系在溶劑質(zhì)量比為1:1的條件下,經(jīng)修正的UNIFAC方程計(jì)算得到的隨復(fù)合溶劑中a組分的增加環(huán)己烷-苯相對揮發(fā)度的變化曲線繪制成圖8。

    圖8 溶劑復(fù)配對環(huán)己烷-苯體系的相對揮發(fā)度的影響Fig.8 Influence of mixed solvent content for cyclohexane-benzene on relative volatility

    從圖8中可以看出這幾組溶劑復(fù)配時(shí)對環(huán)己烷-苯共沸體系的相對揮發(fā)度均有影響,并且都存在一個(gè)最佳的溶劑配比。各組復(fù)合溶劑對體系的影響不相同,從圖中可以看出苯乙酮和苯胺復(fù)配時(shí)對體系環(huán)己烷-苯的影響最大。溶劑復(fù)配的最佳值見表7。

    表7 萃取精餾分離環(huán)己烷-苯體系復(fù)合溶劑備選Table 7 Mixed-solvent of separating cyclohexane-benzene system by extractive distillation

    從表7中可以看出60%苯胺(aniline)+40%苯乙酮(hypnone)復(fù)配時(shí)得到的乙酸乙酯和乙醇的相對揮發(fā)度最大為 2.9404,大于單溶劑分離環(huán)己烷-苯共沸體系的最大相對揮發(fā)度。

    2.3.2 不同組成下的復(fù)合溶劑(苯胺+苯乙酮)對環(huán)己烷-苯體系相對揮發(fā)度的影響 在溶劑質(zhì)量比為1:1時(shí),改變復(fù)合溶劑苯胺和苯乙酮的組成并改變環(huán)己烷-苯體系中環(huán)己烷的質(zhì)量分?jǐn)?shù),所得的環(huán)己烷-苯的相對揮發(fā)度結(jié)果如圖9所示。

    圖9 溶劑組成(苯胺+苯乙酮)對環(huán)己烷-苯相對揮發(fā)度的影響Fig.9 Influence of mass ratio in binary mixed solvent on relative volatility of cyclohexane-benzene

    從圖9中可以看出對于環(huán)己烷-苯體系改變環(huán)己烷的質(zhì)量分?jǐn)?shù)時(shí),復(fù)合溶劑(苯胺+苯乙酮)對體系相對揮發(fā)度的改變趨勢相似,均存一最佳配比,即在苯胺在復(fù)合溶劑中質(zhì)量分?jǐn)?shù)為60%左右時(shí),體系的相對揮發(fā)度達(dá)到最大值。這表明如果在萃取精餾操作的過程中,即使在萃取精餾段中環(huán)己烷的組成隨塔高變化,使用適宜配比的復(fù)合溶劑,在一定的溶劑比下,都可以保持溶劑對環(huán)己烷-苯體系的高分離效果。

    2.3.3 不同溶劑比對環(huán)己烷-苯體系相對揮發(fā)度的影響 在單溶劑(苯乙酮、苯胺)、復(fù)合溶劑(苯乙酮+苯胺)的最佳組合下(xaniline=0.6)改變?nèi)軇┍?,得到環(huán)己烷-苯共沸體系中相對揮發(fā)度隨溶劑比的變化趨勢,結(jié)果如圖10所示。

    圖10 不同溶劑比對環(huán)己烷-苯體系相對揮發(fā)度的影響Fig.10 Influence of solvent ratio on relative volatility of cyclohexane-benzene

    從圖10中可以看出隨著溶劑比的增加,不管是單一溶劑還是復(fù)合溶劑,體系的相對揮發(fā)度都逐漸增加,在溶劑比較小時(shí)(<3)體系的相對揮發(fā)度增加很快,隨著溶劑比的逐漸增加體系的相對揮發(fā)度增加速度變化緩慢。

    3 溶劑-共沸物系的汽液平衡實(shí)驗(yàn)驗(yàn)證

    為了驗(yàn)證選擇萃取精餾復(fù)合溶劑理論的可靠性,本節(jié)采用平衡釜法結(jié)合氣相色譜分析直接測定常壓下復(fù)合溶劑與被分離體系(乙酸乙酯-乙醇,乙醇-水,環(huán)己烷-苯)的汽液平衡數(shù)據(jù),計(jì)算體系的相對揮發(fā)度,并將實(shí)驗(yàn)值與模擬值進(jìn)行比較并進(jìn)行結(jié)果討論。

    3.1 實(shí)驗(yàn)裝置

    實(shí)驗(yàn)中所用汽液平衡釜采用改進(jìn)的 Othmer循環(huán)釜,其結(jié)構(gòu)如圖11所示。

    圖11 汽液平衡實(shí)驗(yàn)測定裝置Fig.11 Sketch of VLE determination system

    3.2 汽液平衡實(shí)驗(yàn)結(jié)果與分析

    3.2.1 乙酸乙酯-乙醇體系 在101.325 kPa下,采用上述實(shí)驗(yàn)裝置和實(shí)驗(yàn)方法對單溶劑(DMF、DMSO)對乙酸乙酯-乙醇共沸體系相對揮發(fā)度的影響進(jìn)行實(shí)驗(yàn)測定,以驗(yàn)證提出的分子設(shè)計(jì)方法及所用到的修正 UNIFAC方程,結(jié)果分別如圖12、圖13所示。

    圖12 DMF 對乙酸乙酯-乙醇共沸體系相對揮發(fā)度的影響Fig.12 Influence of DMF on relative volatility of ethyl acetate-ethanol azeotropic system

    圖13 DMSO 對乙酸乙酯-乙醇共沸體系相對揮發(fā)度的影響Fig.13 Influence of DMSO on relative volatility of acetate-ethanol azeotropic system

    從圖12、圖13中可以看出,對于乙酸乙酯-乙醇體系的分離,隨著單一溶劑DMSO和DMF溶劑比的增加,乙酸乙酯-乙醇之間的相對揮發(fā)度逐漸增大,DMSO對該體系的相對揮發(fā)度的影響要強(qiáng)于DMF。實(shí)驗(yàn)測定的汽液平衡數(shù)據(jù)及實(shí)驗(yàn)值和模擬值的結(jié)果比較如表8所示,其中分別是液相與汽相達(dá)平衡時(shí)乙酸乙酯脫除溶劑后的質(zhì)量分?jǐn)?shù)。相對偏差δ用式(16)表示

    表8 實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較(共沸體系:乙酸乙酯-乙醇,溶劑:DMF,DMSO)Table 8 Comparison of calculated value and experimental data(azeotropic system: ethyl acetate-ethanol, solvent: DMF,DMSO)

    由表8可以看出單溶劑對乙酸乙酯-乙醇體系相對揮發(fā)度的影響,其實(shí)驗(yàn)值與計(jì)算值較為接近,最大偏差為0.0945,最小偏差為0.0082,說明模擬值與實(shí)驗(yàn)值吻合良好。

    常壓下,溶劑比(質(zhì)量比)為1的條件下,改變復(fù)合溶劑(DMSO+DMF)組成,通過實(shí)驗(yàn)測定并計(jì)算出乙酸乙酯-乙醇共沸體系相對揮發(fā)度的改變與1.2小節(jié)中修正的UNIFAC模型得到的計(jì)算值進(jìn)行比較,結(jié)果如圖14所示。

    圖14 混合溶劑DMSO+DMF組成對乙酸乙酯-乙醇共沸體系相對揮發(fā)度之間的影響Fig.14 Influence of mixed solvent (DMSO+DMF) on relative volatility of ethyl acetate-ethanol azeotropic system

    從圖14中可以看出,乙酸乙酯和乙醇的相對揮發(fā)度隨著溶劑 DMSO含量的增加呈現(xiàn)先增加后減小的變化,存在一最佳復(fù)配組合即當(dāng)DMSO在混合溶劑中的質(zhì)量含量為0.6時(shí),復(fù)合溶劑的分離效果最好,這一點(diǎn)模擬值和實(shí)驗(yàn)值基本吻合。實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較見表9。

    表9 實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較(共沸體系:乙酸乙酯-乙醇,溶劑:DMF+DMSO)Table 9 Comparison of calculated value and experimental data (azeotropic system: ethyl acetate-ethanol, solvent:DMF+DMSO)

    從表8和表9中可以看出,兩種單一溶劑DMF和 DMSO對體系的相對揮發(fā)度的實(shí)驗(yàn)值和用修正的UNIFAC模型的計(jì)算值比較,兩者的相對誤差均在 9%以內(nèi),兩種溶劑復(fù)配時(shí)實(shí)驗(yàn)值與計(jì)算值比較相對誤差也在9%以內(nèi),所以修正的UNIFAC模型可以較好地應(yīng)用于萃取精餾分離乙酸乙酯-乙醇體系單溶劑和復(fù)合溶劑的選擇。

    3.2.2 乙醇-水體系 在常壓下,對乙醇-水共沸體系,通過實(shí)驗(yàn)測定的單一溶劑DMF和DMSO對乙醇-水的相對揮發(fā)度的影響和用修正的 UNIFAC模型得到的計(jì)算值進(jìn)行比較,結(jié)果分別如圖15、圖16所示。

    從圖15、16中可以看出,對于乙醇-水體系無論是DMF還是DMSO,隨著溶劑比的增加輕重關(guān)鍵組分之間的相對揮發(fā)度都在增加,使體系更易于分離。在相同的溶劑比條件下,乙醇-水體系在DMSO中的相對揮發(fā)度均要高于 DMF。在溶劑比較小時(shí)體系的相對揮發(fā)度隨溶劑比增加較快。當(dāng)溶劑比較大時(shí),相對揮發(fā)度的變化相對緩慢。實(shí)驗(yàn)結(jié)果和計(jì)算結(jié)果增加的趨勢基本相同。實(shí)驗(yàn)測定的汽液平衡數(shù)據(jù)及實(shí)驗(yàn)值和模擬值的結(jié)果比較見表10。其中分別是液相與汽相達(dá)平衡時(shí)乙醇脫除溶劑后的質(zhì)量分?jǐn)?shù)。

    圖15 DMF對乙醇-水體共沸系相對揮發(fā)度的影響Fig.15 Influence of DMF on relative volatility of ethanol-water azeotropic system

    圖16 DMSO對乙醇-水共沸體系相對揮發(fā)度的影響Fig.16 Influence of DMSO on relative volatility of ethanol-water azeotropic system

    表10 實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較(共沸體系:乙醇-水,溶劑:DMF、DMSO)Table 10 Comparison of calculated value and experimental data(azeotropic system: ethanol-water, solvent:DMF, DMSO)

    在101.325 kPa下,在溶劑比為1(質(zhì)量比)的條件下,改變復(fù)合溶劑(DMSO+DMF)的組成得到乙醇-水共沸體系的相對揮發(fā)度的實(shí)驗(yàn)數(shù)據(jù)和用修正的 UNIFAC模型的計(jì)算值比較,結(jié)果如圖17所示。

    圖17 復(fù)合溶劑(DMSO+DMF)組成與乙醇-水共沸體系相對揮發(fā)度的影響Fig.17 Influence of mixed solvent (DMSO+DMF) on relative volatility of water-ethanol azeotropic system

    從圖17中可以看出兩種溶劑均能改變乙醇和水的相對揮發(fā)度,在兩種溶劑DMF和DMSO單獨(dú)作用下乙醇和水的相對揮發(fā)度分別為 1.2050和1.4315,兩者復(fù)合后隨著溶劑中DMSO含量的增加乙醇和水的相對揮發(fā)度逐漸增加,當(dāng)復(fù)合溶劑中DMSO質(zhì)量含量為70%時(shí),乙醇與水的相對揮發(fā)度的實(shí)驗(yàn)值達(dá)到最大為1.6809,之后隨著DMSO的含量增加乙醇-水的相對揮發(fā)度降低直至逐漸與純DMSO接近。計(jì)算值與實(shí)驗(yàn)值的變化趨勢基本相同。實(shí)驗(yàn)測定的汽液平衡數(shù)據(jù)及實(shí)驗(yàn)值和模擬值的結(jié)果比較見表11。

    表11 實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較(共沸體系:乙醇-水,溶劑:DMF+DMSO)Table 11 Comparison of calculated value and experimental data (azeotropic system: ethanol-water,solvent: DMF+DMSO)

    從表10、表11中可以看出,兩種單一溶劑DMF和DMSO的實(shí)驗(yàn)值和用修正的UNIFAC模型得到的計(jì)算值比較兩者的相對誤差均在 8%以內(nèi),兩種溶劑復(fù)配時(shí)實(shí)驗(yàn)值與計(jì)算值的相對誤差均在 6%以內(nèi),所以修正的 UNIFAC模型可以較好地應(yīng)用于萃取精餾分離乙醇-水體系單溶劑和復(fù)合溶劑的選擇。

    以上兩種共沸體系均為極性很強(qiáng)的體系,UNIFAC模型在建模過程中未考慮分子特性,如分子締合等,所以對于強(qiáng)極性混合物系的計(jì)算難免會(huì)出現(xiàn)誤差。但對于兩種體系所有的實(shí)驗(yàn)和計(jì)算誤差均在9%以內(nèi),驗(yàn)證了本文所選擇的模型的可行性,也說明了本文提出的復(fù)合溶劑復(fù)配方法的可行性。

    3.2.3 環(huán)己烷-苯體系 在101.325 kPa下,對環(huán)己烷-苯共沸體系,通過實(shí)驗(yàn)測定的單一溶劑苯乙酮和苯胺對環(huán)己烷-苯的相對揮發(fā)度的影響和用修正的UNIFAC模型得到的計(jì)算值進(jìn)行比較,結(jié)果分別如圖18、圖19所示。

    圖18 苯乙酮對環(huán)己烷-苯共沸體系相對揮發(fā)度的影響Fig.18 Influence of hypnone on relative volatility of cyclohexane-benzene azeotropic system

    圖19 苯胺對環(huán)己烷-苯共沸體系相對揮發(fā)度的影響Fig.19 Influence of aniline on relative volatility of cyclohexane-benzene azeotropic system

    從圖18、圖19中可以看出,實(shí)驗(yàn)結(jié)果和計(jì)算結(jié)果增加的趨勢基本相同。實(shí)驗(yàn)測定的汽液平衡數(shù)據(jù)及實(shí)驗(yàn)值和模擬值的結(jié)果比較見表12。其中xshex、yshex分別是液相與汽相達(dá)平衡時(shí)環(huán)己烷脫除溶劑后的質(zhì)量分?jǐn)?shù)。

    表12 實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較Table 12 Comparison of calculated value and experimental data (azeotropic system: cyclohexane-benzene,solvent: aniline, hypnone)

    在101.325 kPa下,在溶劑比為1(質(zhì)量比)的條件下,改變復(fù)合溶劑(苯胺+苯乙酮)的組成得到環(huán)己烷-苯共沸體系的相對揮發(fā)度的實(shí)驗(yàn)數(shù)據(jù)和用修正的UNIFAC模型的計(jì)算值,結(jié)果如圖20所示。

    圖20 復(fù)合溶劑(苯胺+苯乙酮)組成對乙醇-水共沸體系相對揮發(fā)度的影響Fig.20 Influence of mixed solvent (aniline +hypnone)on relative volatility of cyclohexane-benzene azeotropic system

    從圖20中可以看出,環(huán)己烷和苯的相對揮發(fā)度隨著溶劑中苯胺含量的增加呈現(xiàn)先增加后減小的變化,存在一最佳復(fù)配組合即當(dāng)苯胺在混合溶劑中的質(zhì)量含量為 0.6時(shí),復(fù)合溶劑的分離效果最好,這一點(diǎn)實(shí)驗(yàn)值和模擬值基本吻合。實(shí)驗(yàn)測定的汽液平衡數(shù)據(jù)及實(shí)驗(yàn)值和模擬值的結(jié)果比較見表13。

    表13 實(shí)驗(yàn)值和計(jì)算值的結(jié)果比較Table 13 Comparison of calculated value and experimental data(azeotropic system: cyclohexane-benzene,solvent: aniline+hypnone)

    從表12、表13中可以看出,兩種單一溶劑苯胺和苯乙胺的實(shí)驗(yàn)值和用修正的UNIFAC模型得到的計(jì)算值比較兩者的相對誤差均在 7%以內(nèi),兩種溶劑復(fù)配時(shí)試驗(yàn)值與計(jì)算值的相對誤差均在 7%以內(nèi),所以修正的UNIFAC模型可以較好地應(yīng)用于萃取精餾分離極性較弱的環(huán)己烷-苯體系單溶劑和復(fù)合溶劑的選擇。

    4 結(jié) 論

    提出按溶劑之間復(fù)配形成氫鍵的情況和修正的UNIFAC模型相結(jié)合的方法,以乙酸乙酯-乙醇、乙醇-水、環(huán)己烷-苯共沸體系為例,預(yù)測出萃取精餾分離這些共沸體系的復(fù)配溶劑及復(fù)合溶劑的最佳配比如下:

    (1)分離乙酸乙酯-乙醇共沸體系的最佳復(fù)合溶劑及配比(質(zhì)量分?jǐn)?shù))為60%DMSO+40%DMF;

    (2)分離乙醇-水共沸體系的最佳復(fù)合溶劑及配比(質(zhì)量分?jǐn)?shù))為70%DMSO+30%DMF;

    (3)分離環(huán)己烷-苯共沸體系的最佳復(fù)合溶劑及配比(質(zhì)量分?jǐn)?shù))為60%苯胺+40%苯乙酮。

    利用汽液平衡釜測定了在101.325 kPa下,在單溶劑和復(fù)合溶劑作用下乙酸乙酯-乙醇、乙醇-水、環(huán)己烷-苯3個(gè)體系的汽液平衡數(shù)據(jù),將實(shí)驗(yàn)得到的相對揮發(fā)度與采用修正UNIFAC方程估算得到的相對揮發(fā)度值進(jìn)行了比較,結(jié)果表明計(jì)算值與實(shí)驗(yàn)值吻合較好,3個(gè)體系的相對揮發(fā)度估算值最大相對偏差均小于9%。

    在最佳復(fù)合溶劑組成下,原共沸體系的相對揮發(fā)度達(dá)到最大值,該值優(yōu)于組成復(fù)合溶劑的單溶劑所達(dá)到的分離效果。隨著復(fù)合溶劑的溶劑比增加,原料體系的相對揮發(fā)度的增加,增加幅度均呈現(xiàn)先大后小的趨勢。在最佳復(fù)合溶劑組成下,改變原料體系組成,復(fù)合溶劑組成對原料體系相對揮發(fā)度的影響趨勢基本相同,這表明復(fù)合溶劑在萃取精餾塔內(nèi)不同原料組成下均能較好增加體系相對揮發(fā)度,為工業(yè)上在萃取精餾中使用復(fù)合溶劑提供了理論依據(jù)。

    符 號(hào) 說 明

    a,b,c,ψ——基團(tuán)相互作用參數(shù)

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

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

    Sp——溶解性

    T——溫度

    X——混合物中的摩爾分?jǐn)?shù)

    x——液相中的質(zhì)量分?jǐn)?shù)

    γS,∞——在溶劑S中無限稀釋活度系數(shù)

    α——組分的相對揮發(fā)度

    Γ——過?;疃认禂?shù)

    Γ(i)——在純組分i中的過?;疃认禂?shù)

    δ——相對偏差

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

    v(i)k——分子i中基團(tuán)k的個(gè)數(shù)

    上角標(biāo)

    s——平衡

    下角標(biāo)

    cal——模擬值

    exp——實(shí)驗(yàn)值

    i,j——組分i、組分j

    k,m,n——基團(tuán)k、基團(tuán)m、基團(tuán)n

    [1] ORCHILLéS A V, MIGUEl P J, LLOPIS F J,et al. Isobaric vapor-liquidequilibria for the extractive distillation of ethanol+water mixturesusing 1-ethyl-3-methylimidazolium dicyanamide[J]. J. Chem.Eng. Data, 2011, 56(12): 4875-4880.

    [2] ORCHILLéS A V, MIGUEL P J, VERCHER E,et al. Using 1-ethyl-3-methylimidazolium trifluoromethanesulfonate as an entrainer for the extractive distillation of ethanol+water mixtures[J]. J.Chem. Eng. Data, 2010, 55(4): 1669-1674.

    [3] 劉艷升, 胡玉峰, 曹睿, 等. 一種離子液體作為萃取精餾分離苯和環(huán)己烷的溶劑的應(yīng)用: CN200710064448.X[P]. 2008-09-17.LIU Y S, HU Y F, CAO R,et al. Application of ionic liquid as solvent for extractive distillation separation of benzene and cyclohexane:CN200710064448.X[P]. 2008-09-17.

    [4] DOMA?SKA U, REDHI G G, MARCINIAK A,et al. Activity coefficients at infinite dilution measurements for organic solutes and water in the 1-hexyloxymethyl-3-methyl-imidazolium and 1,3-dihexyloxymethyl-imidazolium bis(trifluoromethylsulfonyl)-imide ionic liquids—the cation influence[J]. Fluid Phase Equilibria, 2009,286(2): 154-161.

    [5] DOMA?SKA U, LASKOWSKA M. Measurements of activity coefficients at infinite dilution of aliphatic and aromatic hydrocarbons,alcohols, thiophene, tetrahydrofuran, MTBE, and water in ionic liquid[BMIM][SCN] using GLC[J]. Journal of Chemical Thermodynamics,2009, 41(5): 645-650.

    [6] DOMA?SKA U, PADUSZY?SKI K. Measurements of activity coefficients at infinite dilution of organic solutes and water in 1-propyl-1-methylpiperidinium bis{(trifluoromethyl)sulfonyl}imide ionic liquid using GLC[J]. Journal of Chemical Thermodynamics,2010, 42(11): 1361-1366.

    [7] DOMA?SKA U, KRóLIKOWSKA M. Measurements of activity coefficients at infinite dilution in solvent mixtures with thiocyanate-based ionic liquids using GLC technique[J]. J. Phys.Chem. B, 2010, 114(25): 8460-8466.

    [8] DOMA?SKA U, KRóLIKOWSKA M. Measurements of activity coefficients at infinite dilution for organic solutes and water in the ionic liquid 1-butyl-1-methylpiperidinium thiocyanate[J]. J. Chem.Eng. Data, 2011, 56(1): 124-129.

    [9] JONGMANS M T G, SCHUUR B, HAAN A B D. Ionic liquid screening for ethylbenzene/styrene separation by extractive distillation[J]. Ind. Eng. Chem. Res., 2011, 50(18): 10800-10810.

    [10] JONGMANS M T G, SCHUUR B, HAAN A B D. Binary and ternary LLE data of the system (ethylbenzene+ styrene+ 1-ethyl-3-methylimidazolium thiocyanate) and binary VLE data of the system(styrene+ 1-ethyl-3-methylimidazolium thiocyanate)[J]. Journal of Chemical Thermodynamics, 2012, 47: 234-240.

    [11] 唐曉東, 袁嬌陽, 李晶晶, 等. 萃取精餾溶劑的研究應(yīng)用進(jìn)展[J].現(xiàn)代化工, 2013, 33(12): 31-35.TANG X D, YUAN J Y, LI J J,et al. Research and application progresses in solvents for extractive distillation[J]. Modern Chemical Industry, 2013, 33(12): 31-35.

    [12] LEI Z G, WANG H Y, ZHOU R Q,et al. Influence of salt added to solvent on extractive distillation[J]. Chemical Engineering Journal,2002, 87(2): 149-156.

    [13] 李東風(fēng), 劉智信, 程建民, 等. 一種碳五分離用混合溶劑:CN201110178859. 8[P]. 2013-01-02.LI D F, LIU Z X, CHENG J M,et al. Mixed solvent for separating C5:CN201110178859.8[P]. 2013-01-02.

    [14] ZHANG J, YUE G J, DONG H X,et al. Liquid-liquid equilibria determination and correlation for systen of benzene + cyclohexane+N,N-dimethylformamide + potassium thiocyanate[J]. Journal of Chemical Engineering of Chinese Universities, 2013, 27(2):330-333.

    [15] LEE F M. Extractivedistillation: close-boiling-point, adding co-solvents can make tough separations technically and economically feasible[J]. Chemical Engineering, 1998, 105 (12): 112-120.

    [16] LEE F M, BROWN R E. Extractive distillation of hydrocarbons employing solvent mixture: US4921581[P]. 1990-05-01.

    [17] BROWN R E, LEE F M. Extractive distillation of hydrocarbon feed employing mixed solvent: US 4954224[P]. 1990-09-04.

    [18] 祝石華, 唐文成, 田龍勝. NFM-COS復(fù)合溶劑萃取精餾分離苯[J].石油化工, 2002, 31(1): 24-27.ZHU S H, TANG W C, TIAN L S. Extractive distillation process for benzene recovery with NFM-COS solvent[J]. Petrochemical Technology, 2002, 31(1): 24-27.

    [19] 曹新波, 湯志剛. NFM-MEK混合溶劑萃取精餾分離正丁烯/丁烷的流程模擬與參數(shù)分析[J]. 計(jì)算機(jī)與應(yīng)用化學(xué), 2005, 22(5):366-369.CAO X B, TANG Z G. Process simulation and parameter analysis for separatingn-butene and butane by extractive distillation with the mixture ofN-formylmorpholine (NFM) and methyl ethyl ketone(MEK) as solvent[J]. Computers and Applied Chemistry, 2005, 22(5):366-369.

    [20] 華超, 李鑫鋼, 徐世民, 等.N-甲酰嗎啉-乙二醇混合溶劑間歇萃取精餾苯-甲基環(huán)戊烷的研究[J]. 天然氣化工(C1化學(xué)與化工), 2006,31(6): 22-27.HUA C, LI X G, XU S M,et al. Study on the batch extractive distillation of azerotropic mixture benzene-methyl cyclopentane by mixed solventN-formymorpholine-ethylene glycol[J]. Natural Gas Chemistry, 2006, 31(6): 22-27.

    [21] YIN W, DING S H, XIA S Q,et al. Cosolvent selection for benzene-cyclohexane separation in extractive distillation[J]. Journal of Chemical & Engineering Data, 2010, 55(9): 3274-3277.

    [22] 張志剛, 徐世民, 李鑫鋼. 常規(guī)間歇萃取精餾分離苯-環(huán)己烷的研究[J]. 化學(xué)工程, 2006, 34(4): 5-8.ZHANG Z G, XU S M, LI X G. Study on separating benzene-cyclohexane by regular batch extractive distillation[J].Chemical Engineering, 2006, 34(4): 5-8.

    [23] GMEHLING J, LI J D, SCHILLER M. A Modified UNI FAC Model(Ⅱ): Present parameter matrix and results for different thermodynamic properties[J]. Ind. Eng. Chem. Res., 1993, 32(1):178-193.

    [24] PRAUSNITZ J M, ANDERSON R. Thermodynamics of solvent selectivity in extractive distillation of hydrocarbons[J]. AIChE Journal,1961, 7(1): 96-101.

    [25] KYLE B G, LENG D E. Solvent selection for extractive distillation[J].Ind. Eng. Chem., 1965, 57(2): 92-107.

    [26] TOMAS E R, ECHER C A. Prediction of limiting activity coefficients by a modified separation of cohensive energy density model and UNIFAC model[J]. Ind. Eng. Chem. Prog. Des., 1984,23(2): 194-209.

    [27] SKJOLDJORGENSEN S, KOLBE B, GMEHLING J,et al.Vapor-liquid equilibria by UNIFAC group contribution. revision and extension[J]. Industrial & Engineering Chemistry Process Design &Development, 1979, 18(4): 2352-2355.

    [28] TOCHIGI K, TIEGS D, GMEHLING J,et al. Determination of new ASOG parameters[J]. Chem. Eng. Japan, 1990, 23(4): 453-463.

    [29] ABRAMS D S, PRAUSNITZ J M. Statistical thermodynamics of liquid mixtures: a new expression for the excess Gibbs energy of partly or completely miscible systems[J]. AIChE Journal, 1975, 21(1):116-128.

    [30] WEIDLICH U, GMEHLING J. A modified UNIFAC model (Ⅰ):Prediction of VLE,hE, andγ∞[J]. Industrial & Engineering Chemistry Research, 1987, 26(7): 1372-1381.

    [31] GMEHLING J, LOHMANN J, JAKOB A. A modified UNIFAC model (Ⅲ): Revision and extension[J]. Ind. Eng. Chem. Res., 2008,37(12): 4876-4882.

    [32] GMEHLING J, WITTIG R, LOHMANN J,et al. A modified UNIFAC (Dortmund) model (Ⅳ): Revision and extension[J].Industrial & Engineering Chemistry Research, 2002, 41(6):1678-1688.

    [33] 盧煥章. 石油化工基礎(chǔ)數(shù)據(jù)手冊[M]. 北京: 化學(xué)工業(yè)出版社,1982: 985-1081.LU H Z. Handbook of Basic Data in Petrochemical Engineering[M].Beijing: Chemical Industry Press, 1982: 985-1081

    [34] 張志剛. 混合溶劑間歇萃取精餾的研究[D]. 天津: 天津大學(xué),2003.ZHANG Z G. Study on batch extractive distillation by using mixed solvent[D]. Tianjin: Tianjin University, 2003.

    date:2017-02-21.

    Prof. LI Chunli, lichunli_hebut@126.com

    supported by the Research Foundation for the Doctoral Program of China (20131317120014).

    Combination of mixed solvent in extractive distillation

    FANG Jing, XUAN Bihan, LI Xiaochun, LI Chunli
    (School of Chemical Engineering,Hebei University of Technology,Tianjin300130,China)

    Extractive distillation is an effective separation process to deal with the separation of close-boiling compounds and azeotropic systems. The selection of appropriate solvent is the key factor for feasibility and economy in extractive distillation. A new method of screening mixed solvent of extractive distillation was proposed, which was according to the situation of forming the hydrogen bond and modified UNFAC model. To take the polar systems (ethyl acetate-ethanol, ethanol-water) and non-polar system (cyclohexane-benzene) as examples, the appropriate mixed-solvents and the optimal proportion were sought out. The vapor-liquid equilibrium data of solvents (single-solvents and mixed-solvents) and azeotropic system were determined by modified Othmer kettle. The results showed that the calculated values were in good agreement with the experimental data and the relative error was less than 9%. By taking into the synthetical evaluation to the solvents,the optimal mixed solvents to separate ethyl acetate-ethanol, ethanol-water, and cyclohexane-benzene were obtained respectively. Under the action of optimal mixed-solvent, the relative volatility of feed components reached maximum which was higher than the value by using any single solvent. At different feed components the trend of volatility of feed components changed by a mixed solvent was almost the same, which indicated that the mixed solvent could improve the volatility of feed components in any tray of extractive distillation column efficiently.

    extractive distillation; compound solvents; azeotrope; phase equilibria; separation

    TQ 028

    A

    0438—1157(2017)11—4186—15

    10.11949/j.issn.0438-1157.20170159

    2017-02-21收到初稿,2017-05-15收到修改稿。

    聯(lián)系人:李春利。

    方靜(1977—),女,博士,教授。

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

    猜你喜歡
    環(huán)己烷苯胺乙酸乙酯
    一種有效回收苯胺廢水中苯胺的裝置
    能源化工(2021年6期)2021-12-30 15:41:26
    密蒙花乙酸乙酯萃取層化學(xué)成分的分離與鑒定
    環(huán)己烷催化氧化生產(chǎn)環(huán)己酮催化劑的專利技術(shù)綜述
    化工管理(2021年7期)2021-05-13 00:45:22
    吸收法處理有機(jī)尾氣技術(shù)研究
    環(huán)己烷冷液傾析器腐蝕與選材
    廣西莪術(shù)乙酸乙酯部位的抗血栓作用
    中成藥(2018年6期)2018-07-11 03:01:04
    澤漆乙酸乙酯提取物對SGC7901/DDP多藥耐藥性的逆轉(zhuǎn)及機(jī)制
    中成藥(2017年8期)2017-11-22 03:19:25
    鎖陽乙酸乙酯提取物的雌激素樣作用研究
    抗氧劑壬基二苯胺的合成及其熱穩(wěn)定性
    基于環(huán)己烷甲酸根和2,2′-聯(lián)吡啶配體的雙核錳(Ⅱ)配合物的合成與表征
    久久人妻熟女aⅴ| 一区二区日韩欧美中文字幕| 国产亚洲av片在线观看秒播厂| 国产日韩欧美视频二区| 最近的中文字幕免费完整| videos熟女内射| 岛国毛片在线播放| 最黄视频免费看| 国产精品欧美亚洲77777| 老司机深夜福利视频在线观看 | 欧美中文综合在线视频| 女性被躁到高潮视频| 九九爱精品视频在线观看| 中文字幕人妻丝袜一区二区 | 亚洲国产欧美日韩在线播放| a级片在线免费高清观看视频| 一级片免费观看大全| 各种免费的搞黄视频| 国产99久久九九免费精品| 水蜜桃什么品种好| 热99久久久久精品小说推荐| 日韩一本色道免费dvd| 啦啦啦在线观看免费高清www| 国产免费一区二区三区四区乱码| 在线亚洲精品国产二区图片欧美| kizo精华| 超色免费av| 男女高潮啪啪啪动态图| 国产精品一二三区在线看| h视频一区二区三区| 国产成人精品在线电影| 一级毛片 在线播放| 色播在线永久视频| 免费黄网站久久成人精品| 蜜桃在线观看..| 欧美在线一区亚洲| 男女国产视频网站| 丝袜脚勾引网站| 亚洲国产看品久久| 视频在线观看一区二区三区| 国产视频首页在线观看| 韩国高清视频一区二区三区| 1024香蕉在线观看| 国产精品成人在线| 日日爽夜夜爽网站| 人人妻,人人澡人人爽秒播 | 老鸭窝网址在线观看| 狂野欧美激情性xxxx| 人妻人人澡人人爽人人| 精品亚洲成a人片在线观看| 国产免费现黄频在线看| 可以免费在线观看a视频的电影网站 | 2021少妇久久久久久久久久久| 国产在线一区二区三区精| 亚洲av日韩精品久久久久久密 | 18禁裸乳无遮挡动漫免费视频| 国产免费一区二区三区四区乱码| 涩涩av久久男人的天堂| 男女免费视频国产| 亚洲欧洲日产国产| 天堂俺去俺来也www色官网| 久久天堂一区二区三区四区| 亚洲人成电影观看| 一区二区三区四区激情视频| 午夜精品国产一区二区电影| 国产在线免费精品| 爱豆传媒免费全集在线观看| 国产精品嫩草影院av在线观看| 九九爱精品视频在线观看| 国产免费现黄频在线看| xxx大片免费视频| 无限看片的www在线观看| 极品少妇高潮喷水抽搐| 国产男女超爽视频在线观看| 七月丁香在线播放| 校园人妻丝袜中文字幕| 在线免费观看不下载黄p国产| 亚洲伊人色综图| 精品视频人人做人人爽| 在线精品无人区一区二区三| 大香蕉久久网| 久久久亚洲精品成人影院| 色婷婷av一区二区三区视频| 久久久精品94久久精品| 无限看片的www在线观看| 黄色毛片三级朝国网站| 亚洲中文av在线| 一级片'在线观看视频| 亚洲欧美一区二区三区久久| 老熟女久久久| 国产成人精品久久久久久| 国产成人欧美在线观看 | 亚洲四区av| 亚洲国产看品久久| 如日韩欧美国产精品一区二区三区| 欧美久久黑人一区二区| 丝瓜视频免费看黄片| 久久国产精品大桥未久av| 日韩一卡2卡3卡4卡2021年| 高清欧美精品videossex| 狠狠精品人妻久久久久久综合| 国产淫语在线视频| 成人免费观看视频高清| 欧美日韩视频精品一区| 久久午夜综合久久蜜桃| 欧美日韩精品网址| 极品少妇高潮喷水抽搐| videosex国产| 午夜av观看不卡| 黑人猛操日本美女一级片| 欧美日韩av久久| 日本欧美视频一区| 天堂中文最新版在线下载| 99热网站在线观看| 久久久久久久精品精品| 又黄又粗又硬又大视频| 国产精品一二三区在线看| 精品亚洲乱码少妇综合久久| 99国产精品免费福利视频| 香蕉丝袜av| 中文欧美无线码| 久久久久视频综合| 亚洲国产欧美网| 建设人人有责人人尽责人人享有的| 人妻一区二区av| 别揉我奶头~嗯~啊~动态视频 | 欧美精品av麻豆av| 久久久久久久国产电影| 老司机影院毛片| 中文字幕色久视频| 日韩中文字幕欧美一区二区 | 看非洲黑人一级黄片| 好男人视频免费观看在线| 欧美av亚洲av综合av国产av | 中文字幕制服av| 亚洲成人av在线免费| 各种免费的搞黄视频| 午夜激情av网站| 国产av一区二区精品久久| 午夜福利网站1000一区二区三区| 免费在线观看黄色视频的| 国产高清国产精品国产三级| 夜夜骑夜夜射夜夜干| 中文字幕制服av| 丝袜在线中文字幕| 男人舔女人的私密视频| 老汉色av国产亚洲站长工具| 色婷婷av一区二区三区视频| 在线观看免费高清a一片| 亚洲五月色婷婷综合| 国产精品久久久人人做人人爽| 成人免费观看视频高清| 国产亚洲欧美精品永久| 麻豆精品久久久久久蜜桃| 久久久久精品人妻al黑| 在线观看三级黄色| 另类精品久久| 免费观看人在逋| 欧美中文综合在线视频| 亚洲国产成人一精品久久久| 亚洲欧美精品自产自拍| 两个人看的免费小视频| 国产乱人偷精品视频| 满18在线观看网站| 久久97久久精品| 一个人免费看片子| 91国产中文字幕| 国产视频首页在线观看| 久久久精品国产亚洲av高清涩受| 久久精品亚洲av国产电影网| 一区在线观看完整版| xxx大片免费视频| 丝袜美足系列| 女人爽到高潮嗷嗷叫在线视频| 一边摸一边做爽爽视频免费| 亚洲av在线观看美女高潮| 熟女av电影| 婷婷色麻豆天堂久久| 最新在线观看一区二区三区 | 成人手机av| 男男h啪啪无遮挡| 免费黄频网站在线观看国产| 十八禁人妻一区二区| 中文字幕高清在线视频| 一区二区日韩欧美中文字幕| 国产欧美日韩综合在线一区二区| www日本在线高清视频| 一区福利在线观看| 日本一区二区免费在线视频| 欧美日韩av久久| 亚洲国产日韩一区二区| 亚洲综合色网址| 精品少妇一区二区三区视频日本电影 | 亚洲国产精品一区二区三区在线| 丰满迷人的少妇在线观看| 黑人巨大精品欧美一区二区蜜桃| 午夜免费鲁丝| 国产免费又黄又爽又色| 又大又黄又爽视频免费| 老司机在亚洲福利影院| 国产高清不卡午夜福利| 久久热在线av| 高清在线视频一区二区三区| 亚洲av中文av极速乱| 在线观看免费高清a一片| 精品视频人人做人人爽| 少妇人妻久久综合中文| 韩国高清视频一区二区三区| 肉色欧美久久久久久久蜜桃| 91老司机精品| 午夜免费观看性视频| 欧美人与性动交α欧美精品济南到| av电影中文网址| 18禁国产床啪视频网站| 亚洲国产欧美在线一区| 熟妇人妻不卡中文字幕| 天天躁夜夜躁狠狠久久av| 中文字幕av电影在线播放| 丝袜喷水一区| 国产一区二区三区av在线| 又黄又粗又硬又大视频| 女性生殖器流出的白浆| 亚洲av男天堂| 婷婷色综合www| 三上悠亚av全集在线观看| 亚洲中文av在线| 国产成人91sexporn| 黄色怎么调成土黄色| 蜜桃国产av成人99| 老汉色av国产亚洲站长工具| 欧美激情高清一区二区三区 | 黄网站色视频无遮挡免费观看| 日韩电影二区| a级毛片在线看网站| 一二三四中文在线观看免费高清| 秋霞伦理黄片| 七月丁香在线播放| 婷婷成人精品国产| 男女边摸边吃奶| 色播在线永久视频| 十八禁网站网址无遮挡| 亚洲欧洲精品一区二区精品久久久 | 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www| 麻豆精品久久久久久蜜桃| 免费观看人在逋| 亚洲精品久久午夜乱码| 国产精品国产av在线观看| 午夜老司机福利片| 欧美国产精品va在线观看不卡| 精品一区二区三卡| 亚洲人成77777在线视频| 国产精品无大码| 成人影院久久| 2021少妇久久久久久久久久久| 欧美精品高潮呻吟av久久| 国产伦理片在线播放av一区| 亚洲av国产av综合av卡| 可以免费在线观看a视频的电影网站 | 国产无遮挡羞羞视频在线观看| 国产不卡av网站在线观看| 亚洲精品,欧美精品| 2021少妇久久久久久久久久久| 欧美日本中文国产一区发布| 亚洲欧洲精品一区二区精品久久久 | 国产免费现黄频在线看| 黄片播放在线免费| 青春草视频在线免费观看| 午夜老司机福利片| 久久久久久久大尺度免费视频| 美女大奶头黄色视频| 久久 成人 亚洲| 80岁老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 曰老女人黄片| 国产在线免费精品| 免费观看性生交大片5| 两个人免费观看高清视频| 99久久精品国产亚洲精品| 国产精品久久久久久精品古装| 亚洲精品aⅴ在线观看| 91精品伊人久久大香线蕉| 久久亚洲国产成人精品v| 色婷婷av一区二区三区视频| 狠狠精品人妻久久久久久综合| 最新的欧美精品一区二区| 亚洲四区av| 亚洲av电影在线进入| 丰满迷人的少妇在线观看| 久久狼人影院| 777米奇影视久久| 日日啪夜夜爽| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品中文字幕在线视频| 久久99精品国语久久久| 少妇 在线观看| 亚洲欧美精品自产自拍| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 欧美日韩视频高清一区二区三区二| avwww免费| 18禁裸乳无遮挡动漫免费视频| 男人舔女人的私密视频| 99热国产这里只有精品6| 久久久精品区二区三区| 在线天堂最新版资源| www.精华液| 亚洲av综合色区一区| 亚洲成av片中文字幕在线观看| 亚洲av成人精品一二三区| netflix在线观看网站| 国产精品香港三级国产av潘金莲 | 久久久久人妻精品一区果冻| 考比视频在线观看| 又粗又硬又长又爽又黄的视频| 极品少妇高潮喷水抽搐| 桃花免费在线播放| 伊人久久国产一区二区| 精品一区二区免费观看| 久久久久久久大尺度免费视频| 国产精品偷伦视频观看了| 一级a爱视频在线免费观看| 久久99一区二区三区| 日韩欧美精品免费久久| 美女脱内裤让男人舔精品视频| 在线观看免费高清a一片| 国产精品二区激情视频| 人人妻人人添人人爽欧美一区卜| 午夜福利,免费看| 亚洲精品国产色婷婷电影| 亚洲精品久久成人aⅴ小说| 成人国语在线视频| 欧美成人精品欧美一级黄| 最近最新中文字幕大全免费视频 | 亚洲欧美精品综合一区二区三区| 男女床上黄色一级片免费看| 可以免费在线观看a视频的电影网站 | 久久天堂一区二区三区四区| 国产97色在线日韩免费| 欧美成人精品欧美一级黄| 男女之事视频高清在线观看 | 在线观看三级黄色| 韩国精品一区二区三区| 亚洲精品自拍成人| 又黄又粗又硬又大视频| videosex国产| 丰满迷人的少妇在线观看| 亚洲av福利一区| 少妇人妻精品综合一区二区| 狠狠婷婷综合久久久久久88av| 精品少妇久久久久久888优播| 一区福利在线观看| 操美女的视频在线观看| 老司机靠b影院| 午夜日本视频在线| 国产精品久久久久成人av| 啦啦啦 在线观看视频| www.熟女人妻精品国产| 亚洲国产成人一精品久久久| 午夜日本视频在线| 老司机影院成人| 老司机影院毛片| 啦啦啦中文免费视频观看日本| 日韩视频在线欧美| 国产精品秋霞免费鲁丝片| 韩国高清视频一区二区三区| 国产精品久久久久久人妻精品电影 | av又黄又爽大尺度在线免费看| 亚洲精品国产av蜜桃| 这个男人来自地球电影免费观看 | 国产福利在线免费观看视频| 91精品国产国语对白视频| 久久亚洲国产成人精品v| 色综合欧美亚洲国产小说| www.自偷自拍.com| 日本av免费视频播放| 国产精品一区二区精品视频观看| 成人亚洲精品一区在线观看| 搡老乐熟女国产| 日本91视频免费播放| 国产精品无大码| 国产黄色视频一区二区在线观看| 亚洲av电影在线进入| 久久精品久久精品一区二区三区| 成年动漫av网址| 日本wwww免费看| 蜜桃国产av成人99| 久久精品国产亚洲av高清一级| 久久99热这里只频精品6学生| 巨乳人妻的诱惑在线观看| 香蕉国产在线看| 女性被躁到高潮视频| 黄色一级大片看看| √禁漫天堂资源中文www| 精品人妻一区二区三区麻豆| 大香蕉久久成人网| 国产免费一区二区三区四区乱码| 啦啦啦 在线观看视频| 国产1区2区3区精品| 精品酒店卫生间| 中文字幕高清在线视频| 肉色欧美久久久久久久蜜桃| 9热在线视频观看99| 一级,二级,三级黄色视频| 久久人人爽人人片av| 久久精品亚洲av国产电影网| 久久久久久人人人人人| 性少妇av在线| 亚洲欧美色中文字幕在线| 久久久国产一区二区| 午夜福利乱码中文字幕| 亚洲情色 制服丝袜| 国产毛片在线视频| 国产免费视频播放在线视频| 亚洲成人av在线免费| 国产精品久久久久成人av| 午夜福利网站1000一区二区三区| 国产精品久久久久成人av| 国产男女超爽视频在线观看| 亚洲精品av麻豆狂野| 亚洲精品久久成人aⅴ小说| 菩萨蛮人人尽说江南好唐韦庄| 免费看av在线观看网站| 麻豆av在线久日| 久久 成人 亚洲| 亚洲成人国产一区在线观看 | 操美女的视频在线观看| 丝袜在线中文字幕| 视频区图区小说| 午夜精品国产一区二区电影| 亚洲美女视频黄频| 免费少妇av软件| 免费久久久久久久精品成人欧美视频| a级片在线免费高清观看视频| 久久久久久久精品精品| 亚洲成人av在线免费| 日本欧美国产在线视频| 18在线观看网站| 国产成人欧美| 亚洲免费av在线视频| 男女无遮挡免费网站观看| 国产人伦9x9x在线观看| 岛国毛片在线播放| 婷婷色综合www| 老鸭窝网址在线观看| 国产精品一二三区在线看| 婷婷色综合大香蕉| 久久久国产欧美日韩av| www.精华液| 国产毛片在线视频| 精品卡一卡二卡四卡免费| 亚洲欧美一区二区三区黑人| 在线免费观看不下载黄p国产| 色94色欧美一区二区| 好男人视频免费观看在线| 久久女婷五月综合色啪小说| 男人爽女人下面视频在线观看| 一级,二级,三级黄色视频| 亚洲成人一二三区av| 亚洲美女视频黄频| 美女中出高潮动态图| 国产精品一区二区在线不卡| 国产精品.久久久| 最近手机中文字幕大全| 亚洲精品日韩在线中文字幕| 女的被弄到高潮叫床怎么办| 国产97色在线日韩免费| 少妇被粗大的猛进出69影院| 高清不卡的av网站| 国产成人精品久久久久久| 亚洲精品第二区| 少妇被粗大猛烈的视频| 亚洲国产毛片av蜜桃av| 亚洲av电影在线观看一区二区三区| 国产在线视频一区二区| 国产在视频线精品| 国产黄色免费在线视频| 18禁观看日本| 亚洲一级一片aⅴ在线观看| 80岁老熟妇乱子伦牲交| 日日摸夜夜添夜夜爱| 人体艺术视频欧美日本| 久久久久国产精品人妻一区二区| 丝袜人妻中文字幕| 伊人久久大香线蕉亚洲五| 各种免费的搞黄视频| 国产又色又爽无遮挡免| 午夜福利视频精品| 交换朋友夫妻互换小说| 亚洲精品中文字幕在线视频| 一本—道久久a久久精品蜜桃钙片| 久热这里只有精品99| 成人漫画全彩无遮挡| 波野结衣二区三区在线| 九色亚洲精品在线播放| 最黄视频免费看| 搡老岳熟女国产| 精品第一国产精品| 最近中文字幕2019免费版| 日韩 亚洲 欧美在线| 日本午夜av视频| 国产成人系列免费观看| 两个人看的免费小视频| 国产精品av久久久久免费| 热99国产精品久久久久久7| 超色免费av| 亚洲精品一区蜜桃| 欧美国产精品一级二级三级| 大香蕉久久网| 下体分泌物呈黄色| 看非洲黑人一级黄片| 久久国产精品男人的天堂亚洲| 两个人看的免费小视频| 国产熟女欧美一区二区| 免费黄色在线免费观看| 又大又黄又爽视频免费| 欧美黑人精品巨大| 午夜福利视频精品| 亚洲国产欧美一区二区综合| 亚洲国产成人一精品久久久| 久久久国产一区二区| 十八禁网站网址无遮挡| 免费日韩欧美在线观看| 亚洲天堂av无毛| 欧美av亚洲av综合av国产av | 午夜福利免费观看在线| 丝袜脚勾引网站| 在线观看免费午夜福利视频| 日韩精品免费视频一区二区三区| 啦啦啦在线免费观看视频4| 中文字幕人妻丝袜制服| 青春草国产在线视频| 少妇人妻久久综合中文| 18禁国产床啪视频网站| 人妻人人澡人人爽人人| 久久性视频一级片| 亚洲av成人不卡在线观看播放网 | 国产97色在线日韩免费| 国产伦理片在线播放av一区| 最近手机中文字幕大全| 青草久久国产| 国产日韩欧美在线精品| 啦啦啦中文免费视频观看日本| 满18在线观看网站| 亚洲国产成人一精品久久久| 日本猛色少妇xxxxx猛交久久| 精品酒店卫生间| 国产精品99久久99久久久不卡 | 国产精品香港三级国产av潘金莲 | 久久人人97超碰香蕉20202| 男女边吃奶边做爰视频| 亚洲av电影在线进入| 秋霞在线观看毛片| 国产97色在线日韩免费| 欧美人与善性xxx| 老司机亚洲免费影院| 别揉我奶头~嗯~啊~动态视频 | 久久久久久久国产电影| svipshipincom国产片| 国产精品一区二区在线观看99| 国产精品熟女久久久久浪| 观看av在线不卡| 精品一品国产午夜福利视频| 丝袜美足系列| 欧美亚洲 丝袜 人妻 在线| 精品亚洲成a人片在线观看| 9色porny在线观看| 中文字幕另类日韩欧美亚洲嫩草| 婷婷色av中文字幕| 成人三级做爰电影| 成人免费观看视频高清| 精品少妇一区二区三区视频日本电影 | 久久国产亚洲av麻豆专区| 欧美精品一区二区大全| 亚洲av在线观看美女高潮| 久久午夜综合久久蜜桃| 国产伦人伦偷精品视频| 18禁动态无遮挡网站| 欧美国产精品一级二级三级| 99精国产麻豆久久婷婷| 亚洲成人一二三区av| 大片电影免费在线观看免费| 十八禁网站网址无遮挡| 水蜜桃什么品种好| 国产精品无大码| 国产野战对白在线观看| 一级毛片电影观看| 欧美激情极品国产一区二区三区| 亚洲av中文av极速乱| 亚洲欧美激情在线| 一区在线观看完整版| 国产精品久久久av美女十八| 成人18禁高潮啪啪吃奶动态图| 亚洲第一青青草原| 丝袜美足系列| 电影成人av| 伊人久久大香线蕉亚洲五| 大片电影免费在线观看免费| 亚洲久久久国产精品| 中文字幕精品免费在线观看视频| 亚洲av中文av极速乱| 亚洲专区中文字幕在线 | 纯流量卡能插随身wifi吗| 国产精品人妻久久久影院| 久久亚洲国产成人精品v| 丝袜美足系列| 美女大奶头黄色视频| 亚洲精品国产av蜜桃| 啦啦啦在线观看免费高清www| 中文字幕人妻丝袜一区二区 | 中文字幕人妻丝袜制服| 少妇 在线观看| 久久这里只有精品19| 午夜福利乱码中文字幕| 看免费av毛片|