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

    混合離子液體催化反應精餾合成乙酸正己酯

    2021-08-31 07:00:44宋振興崔現(xiàn)寶張纓張雪梅何杰馮天揚王紀孝
    化工學報 2021年8期
    關鍵詞:板數(shù)酯交換甲酯

    宋振興,崔現(xiàn)寶,張纓,張雪梅,何杰,馮天揚,王紀孝

    (天津大學化工學院,化學工程聯(lián)合國家重點實驗室,天津 300350)

    引 言

    乙酸甲酯(MeOAC)是生產(chǎn)PVA(聚乙烯醇)過程中的副產(chǎn)物,工業(yè)上乙酸甲酯的應用十分有限,因此常常通過酯交換反應將乙酸甲酯轉化為高附加值的酯[1-4]。乙酸正己酯(HeAC)是一種重要的溶劑[5],通過乙酸甲酯和正己醇進行酯交換反應,可得到乙酸正己酯產(chǎn)品,同時獲得的甲醇是生產(chǎn)PVA的原料,能夠進行循環(huán)利用。

    酯交換反應的傳統(tǒng)催化劑是無機酸,例如H2SO4、HF和H3PO4等,而無機酸具有較強的腐蝕性,難以回收,容易造成一系列的環(huán)境問題。為了解決這些問題,在酯交換反應中開始選用分子篩、離子交換樹脂等固體催化劑[6-11]。但是固體催化劑存在容易失活、裝填困難以及熱穩(wěn)定性差等問題[12]。離子液體作為一類同時具備優(yōu)良溶劑性能與催化活性的新興綠色催化劑,在各類反應中有著廣泛的應用[13]。Zhang等[14]報道了離子液體在催化酯交換生成生物柴油領域的應用,Deshmukh等[15]進行了離子液體催化碳酸二甲酯和苯酚酯交換反應的研究。混合離子液體,又稱“復鹽離子液體”[16-17],通過離子液體的混合來調(diào)控離子液體的密度、黏度、電導率以及溶解性,為離子液體的設計拓寬了渠道[18-22]。Larriba等[23]進行了混合離子液體萃取的研究,Taige等[24]報道了混合離子液體在電化學中的應用,Yang等[25]報道了單一和混合離子液體催化乙酸甲酯酯交換反應的動力學研究。將混合離子液體作為酯交換反應的催化劑,不僅起到調(diào)控催化劑物理性質的作用,還可以加快化學反應速率。離子液體由于具有極低的飽和蒸氣壓,僅通過簡單蒸餾過程便可回收重復利用,節(jié)省分離過程中的能耗,同時避免環(huán)境污染。

    離子液體催化反應精餾是一種新型的反應分離集成方法,它結合了反應精餾和離子液體各自的優(yōu)勢。反應精餾中反應過程和分離過程相互促進,且離子液體僅采用閃蒸分離便可回收,從而節(jié)省設備投資并降低能耗[26-27]。Peng等[26]報道了離子液體[Hpy][HSO4]催化反應精餾合成乙酸異丁酯。Yang等[28]研究了離子液體催化酯交換反應精餾合成異丙醇。至今未發(fā)現(xiàn)有關混合離子液體催化乙酸甲酯和正己醇(HeOH)的反應精餾研究。

    本文從5種具有不同陰、陽離子結構的質子酸離子液體中篩選出1-丙基磺酸-3-甲基咪唑三氟甲烷磺酸鹽([PSO3HMIm][OTf])作為乙酸甲酯和正己醇酯交換反應的催化劑,采用離子液體1-辛基-2,3-二甲基咪唑雙(三氟甲烷磺酰)亞胺鹽([OMMIm][Tf2N])與其混合獲得混合離子液體催化劑。測定了混合離子液體催化乙酸甲酯酯交換生成乙酸正己酯的反應動力學數(shù)據(jù),擬合實驗數(shù)據(jù)得到反應動力學方程。在反應動力學研究基礎上,建立乙酸甲酯酯交換反應精餾工藝流程,通過模擬計算探究了催化劑用量、回流比、持液量、反應段塔板數(shù)、進料位置及理論板數(shù)對混合離子液體反應精餾過程的影響,獲得了優(yōu)化的操作和設備參數(shù)。

    1 實驗部分

    1.1 實驗試劑

    實驗中所用的試劑如表1所示。所有試劑不需進一步處理,直接使用。所用離子液體的含水量均小于0.0002(經(jīng)卡爾-費休滴定法測定)。

    表1 實驗試劑Table 1 The reagents used in experiment

    1.2 實驗裝置和實驗步驟

    酯交換反應在一個帶有機械攪拌裝置以及裝有溫度計和冷凝器的四口圓底燒瓶中進行,通過恒溫水浴來控制反應體系的溫度,精度為±0.01 K。采用循環(huán)冷凍水進行冷卻,冷凝器的作用是確保揮發(fā)出來的反應組分冷凝并回流至反應體系中。實驗過程中從取樣口進行取樣分析。

    實驗步驟:將稱量好的正己醇和催化劑加入四口燒瓶中進行預熱,相應量的乙酸甲酯加入單口燒瓶中預熱,待預熱溫度均達到設定溫度時,將單口燒瓶中的乙酸甲酯迅速加入四口燒瓶中,同時開始計時。實驗過程中間隔一定時間便進行取樣分析,當四口燒瓶內(nèi)各個組分濃度趨于恒定值時該反應體系達到平衡,計時結束。

    1.3 分析測試方法

    實驗樣品采用氣相色譜儀(SP-1000)分析,采用修正面積歸一法定量。氣相色譜檢測條件為:FID檢測器,采用Wel-30毛細色譜柱,規(guī)格為50 m×0.32 mm×0.5μm。進樣器、柱箱、檢測器溫度分別為383.15、473.15、503.15 K。載氣氮氣、氫氣、空氣流量分別為30、30、300 ml/min。離子液體在氣相色譜中不出峰,因此氣相色譜得到的是扣除離子液體后各組分含量,根據(jù)這些數(shù)據(jù)計算乙酸甲酯轉化率。由于反應是液相反應,反應過程中體積變化很小,反應液的體積采用初始投料體積,離子液體的濃度采用式(1)計算:

    式中,c為離子液體濃度;n為離子液體物質的量;V為初始投料組分的體積。

    2 離子液體催化性能和反應動力學研究

    通過實驗評估,攪拌速率(100~300 r/min)對酯交換反應速率幾乎沒有影響,因此以下所有實驗的轉速值均設定為200 r/min。

    2.1 單一離子液體催化劑的篩選

    設定反應溫度為323.15 K,催化劑濃度為0.2 mol/L,反應物乙酸甲酯和正己醇初始摩爾比為1?;撬峁δ芑疊r?nsted酸性離子液體由于合成路線簡單得到了廣泛的應用[29-30],在酯交換反應中表現(xiàn)出良好的性能[25],因此本文選擇五種不同陰陽離子結構的磺酸功能化Br?nsted酸性離子液體,考察其對該酯交換反應的影響,結果見圖1。

    圖1 不同離子液體催化劑對酯交換反應的影響Fig.1 Effects of different ionic liquid catalysts on transesterification

    由圖1可以看出,不同結構的離子液體具有不同的催化活性,催化活性由大到小是:[PSO3HMIm][OTf]>[BSO3HMIm][OTf]>[BSO3HMIm][HSO4]>[BSO3HPy][OTf]>[BSO3HMIm][Cl]。一般來說,質子酸離子液體催化劑的活性與酸性強弱有關,以上五種離子液體的酸性由大到小依次是:[BSO3HMIm][HSO4]>[PSO3HMIm][OTf]>[BSO3HMIm][OTf]>[BSO3HPy][OTf]>[BSO3HMIm][Cl][25,31-33]。這五種離子液體催化劑的催化活性順序和酸性強弱順序并不完全相同,主要原因是這五種離子液體均不能完全溶解于反應體系當中,其催化活性不僅和催化劑的酸性強弱有關,也和催化劑與反應體系的溶解度相關[34]。離子液體[PSO3HMIm][OTf]酸性較強且與反應體系相溶性較好,呈現(xiàn)出最高的催化活性,因此選擇它為后續(xù)反應催化劑。

    2.2 混合離子液體催化的研究

    如2.1節(jié)所述,離子液體催化性能受到其在反應體系中溶解性的制約,[PSO3HMIm][OTf]不能夠完全溶解于反應體系中。因此,選用極性較弱和黏度較低的離子液體[OMMIm][Tf2N][35]與其混合來提高催化劑在反應體系中的溶解度。設定反應溫度為323.15 K,初始反應物摩爾比為1,催化劑的濃度為0.15 mol/L。[OMMIm][Tf2N]和[PSO3HMIm][OTf]的摩爾比由0.5∶1變化到2∶1,考察不同混合比對乙酸甲酯酯交換反應的影響,結果見圖2。

    由圖2可知,添加[OMMIm][Tf2N]對該反應的反應速率產(chǎn)生影響。當[OMMIm][Tf2N]和[PSO3HMIm][OTf]的摩爾比從0變化到1∶1時,乙酸甲酯酯交換的反應速率逐漸增大,而再增大摩爾比時,反應速率基本維持不變。離子液體[OMMIm][Tf2N]雖然不能夠解離出H+,不能起到催化作用,但是它可以促進催化劑[PSO3HMIm][OTf]在反應體系中的溶解,因此適當增加[OMMIm][Tf2N]的含量反應速率會增大。當[OMMIm][Tf2N]和[PSO3HMIm][OTf]的摩爾比達到1∶1時,催化劑[PSO3HMIm][OTf]能完全溶解于反應物系,因此再增加離子液體[OMMIm][Tf2N]的量時,反應速率不再增大。因此,選用摩爾比為1∶1的[OMMIm][Tf2N]和[PSO3HMIm][OTf]混合離子液體作為后續(xù)乙酸甲酯和正己醇酯交換反應的催化劑。另外需要說明的是,質子酸離子液體催化劑黏度一般很高,使用不便,混合離子液體不僅能夠提高催化活性,也起到降低催化劑體系黏度的作用。

    圖2 不同[OMMIm][Tf2N]∶[PSO3HMIm][OTf]摩爾比對酯交換反應的影響Fig.2 Effect of different molar ratio of[OMMIm][Tf2N]∶[PSO3HMIm][OTf]on transesterification

    2.3 混合離子液體催化合成乙酸正己酯反應動力學

    為方便起見,本節(jié)中設定混合離子液體催化劑濃度為0.145 mol/L,維持[OMMIm][Tf2N]和[PSO3HMIm][OTf]的摩爾比為1∶1,反應溫度為323.15 K,初始乙酸甲酯和正己醇摩爾比為1∶1作為參比反應條件。當研究其中一個參數(shù)的影響時,其他參數(shù)值的設定與參比反應條件保持一致。

    2.3.1 催化劑濃度的影響 催化劑濃度是影響反應速率的一個重要因素,因此實驗分別測定了催化劑濃度為0.06、0.10、0.145、0.20、0.25 mol/L條件下的反應轉化率變化曲線。如圖3所示,隨著催化劑濃度的升高,該反應體系達到平衡的時間縮短,由于催化劑不影響反應平衡,因此,不同濃度下的乙酸甲酯平衡轉化率維持不變。

    圖3 催化劑濃度對酯交換反應的影響Fig.3 Effects of catalyst concentration on transesterification

    本節(jié)研究了初始反應速率與催化劑濃度的關系,通過擬合初始反應速率和催化劑濃度關系方程,可以確定動力學方程中催化劑濃度的指數(shù)參數(shù)[36]。初始反應速率與催化劑濃度關系由式(2)表達,擬合得到該方程中的模型參數(shù)為E=0.3474 min-1,F(xiàn)=1。擬合結果如圖4所示。

    圖4 初始反應速率與催化劑濃度關系Fig.4 Initial reaction rate versus IL catalyst concentration

    式中,r0為初始反應速率,mol/(L·min);ci為組分i濃度,mol/L;vi為組分i化學計量系數(shù);ccat為催化劑[PSO3HMIm][OTf]的濃度,mol/L。

    2.3.2 反應溫度對酯交換的影響 選取313.15、318.15、323.15、327.65、333.15 K共5個溫度來考察反應溫度對乙酸甲酯反應速率和轉化率的影響,結果如圖5所示。

    由圖5可知,溫度對該反應的影響明顯,溫度越高,反應達到平衡的時間越短,反應速率越大。溫度對反應速率的影響可以由Arrhenius方程來解釋:溫度升高造成分子運動更加劇烈,因此發(fā)生了更多的有效碰撞使得反應速率加快。然而,由圖5可以看出,溫度對乙酸甲酯和正己醇酯交換反應的平衡轉化率影響較小,這主要是因為該反應的反應焓較低,因此受溫度影響較小。

    圖5 反應溫度對酯交換反應的影響Fig.5 Effects of reaction temperature on transesterification

    2.3.3 反應物摩爾比的影響 乙酸甲酯和正己醇酯交換是一個可逆反應,因此反應物初始摩爾比(正己醇∶乙酸甲酯)是影響反應的一個重要因素。本節(jié)考察了醇酯摩爾比0.25~4下的乙酸甲酯轉化率隨時間變化關系,結果如圖6所示。

    從圖6中可以看出,反應物初始摩爾比對平衡轉化率有著顯著的影響,乙酸甲酯的平衡轉化率隨著初始醇酯摩爾比R0的增大而增大,此現(xiàn)象產(chǎn)生主要是由于正己醇的過量會造成化學平衡向右移動,因此乙酸甲酯的轉化率會升高。

    圖6 正己醇和乙酸甲酯初始摩爾比對酯交換反應的影響Fig.6 Effect of initial molar ratio of n-hexanol to methyl acetate on the transesterification

    2.4 離子液體催化劑的回收

    對于離子液體催化酯交換反應來說,離子液體的可重復使用性能是至關重要的。因此,為了探究離子液體的復用性,設定反應條件為:反應溫度323.15 K,混合離子液體催化劑濃度0.25 mol/L,醇酯比1∶1?;旌想x子液體催化劑回收采取真空簡單蒸餾方法,由于離子液體的飽和蒸氣壓幾乎可以忽略,[OMMIm][Tf2N][35]和[PSO3HMIm][OTf][37]熱分解溫度較高,可以通過簡單蒸餾將反應體系的反應物和生成物去除以得到原來的催化劑。真空簡單蒸餾回收條件為:壓力500 Pa,溫度180℃。將催化劑重復回收5次進行實驗測定,結果如圖7所示,實驗證明催化劑在重復使用五次后活性無明顯降低,證明了離子液體催化劑具有良好的復用性能。

    圖7 混合離子液體催化劑的回收性能Fig.7 Performance of recovered mixed ionic liquid catalyst

    2.5 化學平衡的研究

    在動力學參比反應條件下分別測定了313.15、318.15、323.15、327.65、333.15 K 5個溫度下的化學平衡常數(shù)(Ka),乙酸甲酯和正己醇進行酯交換反應的化學平衡常數(shù)計算公式為:

    化學平衡常數(shù)由反應物和生成物的摩爾分數(shù)(x)和活度系數(shù)(γ)計算得到。其中,活度系數(shù)可由NRTL方程計算得到:

    式中,Bij為二元交互作用參數(shù);αij為非隨機參數(shù);T為溫度。乙酸甲酯與甲醇、甲醇與正己醇的模型參數(shù)來源于Aspen數(shù)據(jù)庫,[OMMIm][Tf2N]與乙酸甲酯和甲醇的模型參數(shù)來源于文獻值[38],COSMOSAC模型能夠較為準確地預測活度系數(shù)[39-40],并在離子液體精餾研究中被廣泛應用[27,41],因此其余缺失的NRTL模型參數(shù)通過COSMO-SAC模型預測得到。所用到的NRTL方程參數(shù)見表2。

    表2 NRTL模型參數(shù)Table 2 NRTL model parameters

    根據(jù)van’t Hoff等壓方程[式(5)],對Ka和溫度T進行擬合便可以得到ΔrH0,擬合結果如式(6)和圖8所示。

    圖8 反應溫度對化學平衡常數(shù)的影響Fig.8 Effect of reaction temperature on chemical equilibrium constant

    通過式(5)和式(6)的對比,可以得出ΔrH0的擬合值,標準反應焓ΔrH0也可以由反應體系四種物質的標準生成焓ΔfH0計算得到。四種物質的ΔfH0數(shù)值來源于NIST數(shù)據(jù)庫,計算和實驗兩種方法得到的值列于表3中。由表3可知,乙酸甲酯和正己醇進行的酯交換反應為吸熱反應,且反應焓的數(shù)值較小。

    表3 標準反應焓計算結果Table 3 Calculation results of standard reaction enthalpy

    2.6 反應動力學的擬合

    本課題組已有多篇文獻報道離子液體催化酯交換反應動力學[12,42-43],酯交換反應屬于二級反應,乙酸甲酯和正己醇酯交換反應動力學采用非理想擬均相模型(NIH)描述為:

    式中,ρ為反應體系的摩爾密度,ρ=n/V;α為對應物質的活度;k+、k-分別為正反應速率常數(shù)和逆反應速率常數(shù),可以通過計算Arrhenius方程得到:

    式中,A+、A-分別為正反應指前因子和逆反應指前因子;E+、E-分別為正反應活化能和逆反應活化能;R為氣體常數(shù),8.314 J/(mol·K);T為反應溫度,K。

    根據(jù)動力學NIH模型,將實驗測得的轉化率隨時間變化的數(shù)據(jù)點采用非線性最小二乘法回歸擬合得到反應動力學參數(shù),擬合中的目標函數(shù)為乙酸甲酯轉化率的平均相對偏差(MRD):

    式中,XMeOAC,cal,i為計算出的轉化率;XMeOAC,exp,i為實際轉化率;N為數(shù)據(jù)點數(shù)。最終得到混合離子液體催化乙酸甲酯和正己醇酯交換反應的動力學方程(MRD為0.0312):

    3 反應精餾過程模擬

    3.1 工藝流程設計及參比條件的選擇

    本文采用Aspen Plus中的RadFrac模塊和Flash模塊對整個工藝流程進行計算,RadFrac模塊采用平衡級模型,反應動力學方程見式(11)。由于Aspen Plus數(shù)據(jù)庫中沒有離子液體的物性參數(shù),需要采用Aspen的物性估算系統(tǒng)(PCES)根據(jù)離子液體的分子結構進行估算,其中離子液體的理想氣體熱容參數(shù)由Aspen Plus軟件采用Joback方法進行估算,氣相焓采用理想氣體熱容模型進行計算,離子液體的飽和蒸氣壓采用Extended Antoine飽和蒸氣壓模型進行計算,離子液體的臨界性質TC、PC、VC、ZC采用基團貢獻法進行預測得到。物性方法選用NRTL,涉及的模型參數(shù)見表2。

    由于反應物系中存在乙酸甲酯和甲醇的共沸物[44],因此采用雙塔流程。整個反應精餾流程如圖9所示(純度為摩爾分數(shù)),包括一個反應精餾塔(T2)和一個普通精餾塔(T1)。反應精餾塔是乙酸甲酯和正己醇進行化學反應和分離反應物系的場所,普通精餾塔的主要作用是分離出純度較高的產(chǎn)品甲醇,閃蒸罐主要是將催化劑和產(chǎn)品乙酸正己酯進行一次閃蒸分離,回收催化劑重復使用。

    整個流程的進料流股和初步模擬結果如圖9所示。接近共沸組成的乙酸甲酯和甲醇混合物從T1塔的第10塊板進料,混合離子液體和正己醇分別從T2塔的第5和第10塊板進料。T1塔分離得到的接近共沸組成的乙酸甲酯和甲醇混合物進入T2塔的第50塊板,T2塔頂分離出的混合物返回T1塔的第20塊板。各設備的操作條件分別是:T1常壓操作,回流比2,總理論板數(shù)30;T2常壓操作,總理論板數(shù)60,回流比3,持液量110 L;閃蒸罐溫度設定為180℃,壓力為500 Pa。通過模擬計算,分別在T1塔底和閃蒸罐得到純度為98.7%的甲醇和98.1%的乙酸正己酯。

    反應精餾塔T2是整個流程的核心,因此單獨研究了反應精餾塔的操作和設備參數(shù)對產(chǎn)品乙酸正己酯純度和塔底乙酸正己酯收率的影響。選取圖9中T2塔和閃蒸罐標示出的進料條件和初步模擬結果作為參比條件,當改變其中一個參數(shù)作為變量時,保持其他參數(shù)值不變。

    圖9 催化反應精餾流程圖Fig.9 Flow chart of catalytic reactive distillation

    3.2 催化劑用量的影響

    混合離子液體催化劑[OMMIm][Tf2N]和[PSO3HMIm][OTf](比例為1∶1)的用量對模擬結果的影響如圖10所示。從圖10中可以看出,乙酸正己酯的純度隨著催化劑用量的不斷增加呈現(xiàn)先增大后減小的趨勢。當催化劑用量小時,反應速率較低,反應物轉化不完全,因此正己醇轉化率較低,乙酸正己酯純度較低。當催化劑用量過大時,反應速率高,但是混合離子液體催化劑能夠增大乙酸甲酯相對于甲醇的相對揮發(fā)度(圖11),造成反應精餾塔塔頂餾出物中乙酸甲酯含量增加、甲醇含量減少,不利于反應平衡向產(chǎn)物方向移動,從而也造成反應不完全。反應精餾是反應過程和精餾過程的耦合過程,當催化劑用量小時化學反應為控制步驟,產(chǎn)品純度隨著催化劑用量變大而增大。但隨著催化劑用量變大,精餾分離成為控制步驟,因此該過程存在一個最佳的催化劑用量,其值為1 kmol/h。

    圖10 催化劑用量的影響Fig.10 The effect of catalyst flow rate

    圖11 混合離子液體催化劑不同濃度下乙酸甲酯-甲醇體系的y-x相圖Fig.11 The y-x phase diagram of methyl acetate methanol system with different concentrations of mixed ionic liquid catalyst

    3.3 理論板數(shù)的影響

    理論板數(shù)是精餾中的重要參數(shù),理論板數(shù)過小則達不到分離要求,過大會造成設備投資費用的增加,因此,選擇合適的理論板數(shù)是精餾塔設計中重要的一環(huán)。理論板數(shù)對反應精餾塔T2的影響如圖12所示。從圖12中可以看出:乙酸正己酯收率和產(chǎn)品純度隨著理論板數(shù)的增加而增加,當理論板數(shù)達到66后基本維持不變。因此,合適的理論板數(shù)為66。

    圖12 理論板數(shù)的影響Fig.12 The effect of theoretical plates

    3.4 反應段塔板數(shù)和進料位置的影響

    反應精餾塔可以分成三段:正己醇進料位置以上為精餾段,主要作用是進行甲醇濃縮;正己醇和乙酸甲酯進料位置之間為反應段,主要作用是進行酯交換反應;乙酸甲酯混合物進料位置之下的部分為提餾段,主要作用是濃縮乙酸正己酯。本節(jié)探究了不同進料位置和反應段塔板數(shù)(乙酸甲酯混合物與正己醇進料之間的塔板數(shù))對反應精餾的影響,分別考察了10種不同進料位置的情況,結果如表4所示。

    從表4中可以看出,反應段塔板數(shù)的增加會使得乙酸正己酯收率增大和產(chǎn)品純度增加,這主要是由于反應段塔板數(shù)增加,意味著反應體積增加,單位時間產(chǎn)物的生成量增加,反應更完全。因此選擇適宜的反應段塔板數(shù)為40?;旌想x子液體和正己醇無論是從不同塔板進料還是混合后從一塊板進料,對乙酸正己酯收率和純度幾乎沒有影響。若維持反應段塔板數(shù)為40不變,從表4中可以看出,當正己醇進料位置從第6塊理論板移動到第10塊理論板時,精餾段理論板數(shù)增加,反應段下移,提餾段理論板數(shù)減小,而乙酸正己酯收率和純度升高,但是當正己醇進料位置繼續(xù)下移時,乙酸正己酯收率和純度開始降低,因此反應物存在最佳進料位置。正己醇和乙酸甲酯混合物的最佳進料位置是第10和50塊板。

    表4 混合離子液體和反應物的進料位置及間隔板數(shù)的影響Table 4 The effect of feeding position and number of spacer plates

    3.5 回流比的影響

    對于普通精餾塔,回流比的增加能夠提高塔頂或塔底產(chǎn)品的純度。但是對于反應精餾來說,由于化學反應和精餾分離是同時進行的,回流比的變化對結果產(chǎn)生較為復雜的影響。本節(jié)研究了回流比(0.8~5.5)對反應精餾塔的影響,結果如圖13所示。從圖13中可以看出,該反應精餾過程存在一個最佳的回流比,這主要是因為在回流比較小時,塔頂正己醇和乙酸正己酯得不到徹底的分離,而當回流比太大時,催化劑被稀釋,反應速率減小,導致了產(chǎn)品收率降低。因此,最佳的回流比選3.75。

    圖13 回流比的影響Fig.13 The effect of reflux ratio

    3.6 持液量的影響

    塔板是化學反應的場所,因此塔板持液量的大小是影響反應精餾的一個重要因素。反應精餾塔持液量一般要求較大,為了保證反應精餾過程所需持液量,可以采用板式塔,例如篩板塔,塔板持液量可以在設計時調(diào)節(jié)溢流堰高度進行調(diào)整。持液量(60~126 L)對模擬結果的影響如圖14所示,從圖14中可以看出,乙酸正己酯收率和產(chǎn)品純度隨著持液量的增大而增大,并慢慢趨于平穩(wěn)。這主要是因為隨著持液量的不斷增加,反應精餾的控制步驟由反應過程向精餾過程過渡,因此持液量的影響最終趨于穩(wěn)定,該工藝條件下合適的持液量為120 L。

    圖14 塔板持液量的影響Fig.14 The effect of liquid holdup of tray

    3.7 優(yōu)化條件及模擬結果

    根據(jù)上述對各個操作和結構參數(shù)分析和優(yōu)化,得到了反應精餾塔的操作條件如下:混合離子液體[OMMIm][Tf2N]和[PSO3HMIm][OTf](比例為1∶1)用量為1 kmol/L,總理論塔板數(shù)為66,回流比為3.75,離子液體、正己醇和乙酸甲酯混合物分別從第5、10和50塊板進料(進料流量如圖9所示),塔板持液量為120 L。經(jīng)過模擬得到乙酸正己酯收率為0.9991,乙酸正己酯的純度為0.9993,塔內(nèi)液相各個組分隨塔板分布如圖15所示。

    圖15 反應精餾塔中各組分液相濃度分布Fig.15 Liquid concentration profiles in the reactive distillation column

    4 結 論

    通過乙酸甲酯和正己醇酯交換反應催化劑的篩選、反應動力學的研究以及反應精餾流程的設計和穩(wěn)態(tài)模擬,得到如下結論。

    (1)5種具有不同結構的離子液體[PSO3HMIm][OTf]、[BSO3HPy][OTf]、[BSO3HMIm][HSO4]、[BSO3HMIm][OTf]和[BSO3HMIm][Cl]均能催化乙酸甲酯與正己醇的酯交換反應,其中[PSO3HMIm][OTf]呈現(xiàn)最高的催化活性。[OMMIm][Tf2N]和[PSO3HMIm][OTf]形成的混合離子液體催化劑的催化活性比單一離子液體催化劑[PSO3HMIm][OTf]高?;旌想x子液體催化劑在循環(huán)使用五次后活性并無明顯降低。

    (2)測定了乙酸甲酯與正己醇酯交換反應的動力學數(shù)據(jù),建立了乙酸甲酯與正己醇的酯交換反應動力學模型,擬合得到反應動力學模型參數(shù)。

    (3)進行了反應精餾流程模擬計算并對反應精餾塔的影響因素進行了分析。乙酸正己酯收率和純度隨著理論板數(shù)、反應段板數(shù)、持液量的增大而增大,反應精餾過程存在最佳的回流比和催化劑用量。在優(yōu)化操作條件下,乙酸正己酯的收率為0.9991,乙酸正己酯的純度為0.9993。

    猜你喜歡
    板數(shù)酯交換甲酯
    二硫化碳-四氯化碳精餾工藝改進與模擬研究
    離子交換樹脂催化合成苯甲酸甲酯
    云南化工(2020年11期)2021-01-14 00:50:52
    Aspen Plus 軟件在低溫甲醇洗工藝中的應用
    化工管理(2020年19期)2020-07-28 02:59:02
    工業(yè)酶法和化學法酯交換在油脂改性應用中的比較
    中國油脂(2020年5期)2020-05-16 11:23:52
    用Excel進行精餾塔理論板數(shù)的計算
    山東化工(2018年19期)2018-10-29 08:12:20
    K/γ-Al2O3催化丙酸甲酯合成甲基丙烯酸甲酯
    化工進展(2015年3期)2015-11-11 09:07:41
    卡前列甲酯栓聯(lián)合鈣劑預防及治療產(chǎn)后出血的效果觀察
    超導可視聯(lián)合卡前列甲酯用于早早孕無痛人流術的效果觀察
    無溶劑體系下表面活性劑修飾的豬胰脂酶催化酯交換反應的研究
    碳基固體酸催化酯交換反應合成長碳鏈脂肪酸酯
    亚洲精品国产av蜜桃| 91久久精品国产一区二区成人| 免费观看a级毛片全部| 国产高潮美女av| 99视频精品全部免费 在线| 99精国产麻豆久久婷婷| 精品酒店卫生间| 在线天堂最新版资源| 日产精品乱码卡一卡2卡三| 少妇人妻久久综合中文| videossex国产| 国产精品久久久久久精品古装| av视频免费观看在线观看| 午夜福利在线观看免费完整高清在| 国产一区二区在线观看日韩| 国产免费视频播放在线视频| 如何舔出高潮| 大香蕉久久网| 国产伦精品一区二区三区视频9| 久久99热这里只频精品6学生| 免费黄色在线免费观看| 97超碰精品成人国产| 交换朋友夫妻互换小说| 丰满乱子伦码专区| 美女高潮的动态| 18+在线观看网站| 天天躁夜夜躁狠狠久久av| 自拍欧美九色日韩亚洲蝌蚪91 | 韩国高清视频一区二区三区| 高清视频免费观看一区二区| 日韩av在线免费看完整版不卡| 又粗又硬又长又爽又黄的视频| 一本久久精品| 男人添女人高潮全过程视频| freevideosex欧美| 免费大片18禁| 亚洲人成网站在线播| 久久鲁丝午夜福利片| 深夜a级毛片| 国产综合精华液| 久久99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 久久99蜜桃精品久久| 搡老乐熟女国产| 少妇人妻精品综合一区二区| 麻豆成人av视频| 久久99蜜桃精品久久| 青青草视频在线视频观看| 99热这里只有是精品在线观看| 欧美高清成人免费视频www| 亚洲精品视频女| 欧美另类一区| 超碰97精品在线观看| 国产伦在线观看视频一区| 国产精品国产三级专区第一集| 一级毛片电影观看| 亚洲av综合色区一区| 欧美人与善性xxx| 少妇人妻一区二区三区视频| 五月玫瑰六月丁香| 99九九线精品视频在线观看视频| 国产高清有码在线观看视频| 欧美zozozo另类| 18禁裸乳无遮挡免费网站照片| 男人和女人高潮做爰伦理| 黄色视频在线播放观看不卡| 热re99久久精品国产66热6| 亚洲欧美日韩另类电影网站 | 男女无遮挡免费网站观看| 亚洲欧美成人精品一区二区| 美女脱内裤让男人舔精品视频| 观看av在线不卡| 一边亲一边摸免费视频| 国产亚洲最大av| 又大又黄又爽视频免费| 久久久成人免费电影| 国产视频首页在线观看| 青青草视频在线视频观看| 18禁裸乳无遮挡免费网站照片| 99热这里只有是精品在线观看| 丝瓜视频免费看黄片| 人人妻人人看人人澡| 日产精品乱码卡一卡2卡三| 国产精品女同一区二区软件| 晚上一个人看的免费电影| 亚洲欧美日韩东京热| 精华霜和精华液先用哪个| 久久青草综合色| 欧美日韩亚洲高清精品| 亚洲电影在线观看av| 18禁裸乳无遮挡动漫免费视频| 熟妇人妻不卡中文字幕| 亚洲不卡免费看| 小蜜桃在线观看免费完整版高清| 久久女婷五月综合色啪小说| 精品人妻偷拍中文字幕| 色婷婷av一区二区三区视频| 青青草视频在线视频观看| 久久久久国产精品人妻一区二区| 丝袜喷水一区| 中文资源天堂在线| 三级国产精品欧美在线观看| 精品国产三级普通话版| 中国国产av一级| 日韩亚洲欧美综合| 又爽又黄a免费视频| 日韩国内少妇激情av| 国产爱豆传媒在线观看| 菩萨蛮人人尽说江南好唐韦庄| 18禁在线播放成人免费| 亚洲欧美清纯卡通| a级毛色黄片| 精品久久久噜噜| 国产成人精品福利久久| 亚洲欧美一区二区三区黑人 | 五月开心婷婷网| 日本免费在线观看一区| 亚洲美女黄色视频免费看| av又黄又爽大尺度在线免费看| 欧美成人a在线观看| 97精品久久久久久久久久精品| 亚洲av在线观看美女高潮| 熟女av电影| 精品熟女少妇av免费看| 国内揄拍国产精品人妻在线| 久久人人爽av亚洲精品天堂 | 国产av精品麻豆| 免费高清在线观看视频在线观看| 2018国产大陆天天弄谢| 啦啦啦中文免费视频观看日本| 欧美日韩国产mv在线观看视频 | 精品国产一区二区三区久久久樱花 | 欧美97在线视频| 免费看光身美女| 日日啪夜夜撸| 亚洲精品色激情综合| 亚洲四区av| 日韩欧美一区视频在线观看 | 综合色丁香网| 极品少妇高潮喷水抽搐| 男人舔奶头视频| 久久99热这里只频精品6学生| av线在线观看网站| 欧美精品一区二区大全| av国产精品久久久久影院| 国产欧美日韩一区二区三区在线 | 亚洲欧美日韩东京热| www.色视频.com| 中文字幕精品免费在线观看视频 | 国产一区二区在线观看日韩| 久久国产亚洲av麻豆专区| 国产人妻一区二区三区在| 小蜜桃在线观看免费完整版高清| 色婷婷av一区二区三区视频| 免费黄色在线免费观看| 老司机影院毛片| 亚洲国产欧美在线一区| 晚上一个人看的免费电影| 夜夜看夜夜爽夜夜摸| 免费久久久久久久精品成人欧美视频 | 一级毛片aaaaaa免费看小| 91午夜精品亚洲一区二区三区| 综合色丁香网| 一级毛片aaaaaa免费看小| 国产真实伦视频高清在线观看| 嫩草影院入口| 国产毛片在线视频| 日韩大片免费观看网站| 99热这里只有是精品在线观看| 女性生殖器流出的白浆| 啦啦啦在线观看免费高清www| 亚洲欧美中文字幕日韩二区| 国产精品人妻久久久久久| 精品国产乱码久久久久久小说| 大香蕉97超碰在线| 色5月婷婷丁香| 亚洲怡红院男人天堂| 男人和女人高潮做爰伦理| 久久av网站| 少妇人妻久久综合中文| 国产69精品久久久久777片| 国产高清不卡午夜福利| 一级毛片电影观看| 国产大屁股一区二区在线视频| 男女边摸边吃奶| 18禁在线播放成人免费| 在现免费观看毛片| 欧美成人午夜免费资源| 人人妻人人看人人澡| 久久久久久久久久久免费av| 中文字幕av成人在线电影| 3wmmmm亚洲av在线观看| av一本久久久久| 97热精品久久久久久| 国产精品久久久久久精品古装| 人人妻人人看人人澡| 中文字幕av成人在线电影| 亚洲自偷自拍三级| 18禁裸乳无遮挡免费网站照片| 高清在线视频一区二区三区| 久久久久国产网址| 国产午夜精品久久久久久一区二区三区| 久久久国产一区二区| 国产片特级美女逼逼视频| 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 亚洲国产av新网站| av卡一久久| 青青草视频在线视频观看| 亚洲国产色片| 高清毛片免费看| 看免费成人av毛片| 在线观看一区二区三区| 青春草亚洲视频在线观看| av女优亚洲男人天堂| 97在线人人人人妻| 欧美变态另类bdsm刘玥| 国内少妇人妻偷人精品xxx网站| 日韩不卡一区二区三区视频在线| 久久ye,这里只有精品| 久久精品夜色国产| 国产欧美日韩精品一区二区| 国产爱豆传媒在线观看| 国产成人精品婷婷| 中文在线观看免费www的网站| 久久99精品国语久久久| 亚洲国产色片| 大陆偷拍与自拍| 伦理电影免费视频| 成人美女网站在线观看视频| 国产精品秋霞免费鲁丝片| 韩国av在线不卡| 老司机影院毛片| 成人影院久久| 成人一区二区视频在线观看| 九九久久精品国产亚洲av麻豆| 狂野欧美激情性bbbbbb| 夜夜看夜夜爽夜夜摸| 天堂8中文在线网| 日本欧美视频一区| 99热国产这里只有精品6| 欧美xxxx性猛交bbbb| .国产精品久久| 欧美性感艳星| 日韩视频在线欧美| 日日撸夜夜添| 欧美极品一区二区三区四区| 亚洲欧美日韩无卡精品| 99热6这里只有精品| 日本免费在线观看一区| 精品久久久久久久末码| 亚洲欧美日韩东京热| 亚洲国产欧美人成| 免费少妇av软件| 婷婷色综合大香蕉| 麻豆乱淫一区二区| 欧美日韩综合久久久久久| 国产综合精华液| 成年女人在线观看亚洲视频| 99久久精品热视频| 女的被弄到高潮叫床怎么办| 亚洲美女黄色视频免费看| 国产精品久久久久久精品电影小说 | 少妇人妻 视频| 免费在线观看成人毛片| 在线观看av片永久免费下载| av不卡在线播放| 国产成人精品福利久久| 国产精品女同一区二区软件| 18禁动态无遮挡网站| 精品酒店卫生间| av国产免费在线观看| 一二三四中文在线观看免费高清| 尤物成人国产欧美一区二区三区| 九草在线视频观看| 国产一区二区三区av在线| 国产视频首页在线观看| 久久人人爽人人片av| 欧美激情国产日韩精品一区| 久久国产乱子免费精品| 在线观看免费高清a一片| 99热这里只有精品一区| 国产av码专区亚洲av| av又黄又爽大尺度在线免费看| 在线 av 中文字幕| 美女中出高潮动态图| 午夜福利在线在线| 综合色丁香网| 国产欧美日韩一区二区三区在线 | 亚洲精品成人av观看孕妇| 看十八女毛片水多多多| 大片电影免费在线观看免费| 欧美成人午夜免费资源| 高清黄色对白视频在线免费看 | 51国产日韩欧美| 亚洲国产精品999| 亚洲av男天堂| 久久精品国产鲁丝片午夜精品| 久久国产乱子免费精品| 国产精品久久久久久精品古装| 91午夜精品亚洲一区二区三区| 日本欧美国产在线视频| 在线天堂最新版资源| 国产成人a∨麻豆精品| 免费看光身美女| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品成人av观看孕妇| 永久网站在线| 婷婷色综合www| 久久女婷五月综合色啪小说| 精品久久久久久久末码| 国产精品一区二区三区四区免费观看| 91在线精品国自产拍蜜月| 日韩在线高清观看一区二区三区| 亚洲精品国产成人久久av| 亚洲第一av免费看| 久久久久视频综合| av在线老鸭窝| 精品午夜福利在线看| 日日啪夜夜爽| 亚洲国产精品专区欧美| 免费黄网站久久成人精品| 免费高清在线观看视频在线观看| 建设人人有责人人尽责人人享有的 | 亚洲内射少妇av| 国产精品麻豆人妻色哟哟久久| 日本wwww免费看| 亚洲精品国产av蜜桃| 九九在线视频观看精品| 多毛熟女@视频| 18禁动态无遮挡网站| 亚洲久久久国产精品| 99热国产这里只有精品6| 国产免费一区二区三区四区乱码| 高清在线视频一区二区三区| 日韩人妻高清精品专区| 人人妻人人添人人爽欧美一区卜 | 亚洲国产高清在线一区二区三| 免费av中文字幕在线| 啦啦啦啦在线视频资源| 老司机影院毛片| 啦啦啦在线观看免费高清www| 嘟嘟电影网在线观看| av.在线天堂| 精品一区在线观看国产| 精品99又大又爽又粗少妇毛片| 成人毛片60女人毛片免费| 观看美女的网站| 亚洲av成人精品一区久久| 看非洲黑人一级黄片| 亚洲婷婷狠狠爱综合网| 联通29元200g的流量卡| 亚洲熟女精品中文字幕| 日韩欧美精品免费久久| 新久久久久国产一级毛片| 男女啪啪激烈高潮av片| 久久影院123| 老熟女久久久| 成人影院久久| 妹子高潮喷水视频| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级 | 亚洲成人av在线免费| 欧美另类一区| 欧美高清性xxxxhd video| 在线看a的网站| 亚洲av电影在线观看一区二区三区| 欧美成人午夜免费资源| 自拍偷自拍亚洲精品老妇| 午夜老司机福利剧场| 免费少妇av软件| 日本色播在线视频| 在现免费观看毛片| 不卡视频在线观看欧美| 大码成人一级视频| 国产爽快片一区二区三区| 日韩av在线免费看完整版不卡| 午夜老司机福利剧场| av视频免费观看在线观看| 久久久久精品久久久久真实原创| 2021少妇久久久久久久久久久| 久久精品熟女亚洲av麻豆精品| 精品人妻视频免费看| 亚洲欧美中文字幕日韩二区| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级国产专区5o| 欧美日韩一区二区视频在线观看视频在线| 青青草视频在线视频观看| av专区在线播放| 少妇人妻久久综合中文| 麻豆精品久久久久久蜜桃| 身体一侧抽搐| 99久久精品国产国产毛片| 国产精品一区二区三区四区免费观看| 欧美精品国产亚洲| videos熟女内射| 国产精品秋霞免费鲁丝片| 欧美丝袜亚洲另类| 纵有疾风起免费观看全集完整版| 六月丁香七月| a级一级毛片免费在线观看| 少妇人妻久久综合中文| 久久精品夜色国产| 国产高清不卡午夜福利| 少妇裸体淫交视频免费看高清| 日韩欧美一区视频在线观看 | 日韩制服骚丝袜av| 乱系列少妇在线播放| 国产av码专区亚洲av| 久久 成人 亚洲| 赤兔流量卡办理| 亚洲精品色激情综合| 亚洲精品乱码久久久v下载方式| 国产黄频视频在线观看| 久久久久网色| 黄色日韩在线| 一级毛片久久久久久久久女| 国产av一区二区精品久久 | 97精品久久久久久久久久精品| 肉色欧美久久久久久久蜜桃| 91久久精品国产一区二区三区| 免费观看性生交大片5| 又黄又爽又刺激的免费视频.| 最黄视频免费看| 国产欧美日韩一区二区三区在线 | 亚洲成人手机| 看免费成人av毛片| 午夜免费鲁丝| 日韩成人伦理影院| 色哟哟·www| 99热这里只有是精品50| 亚洲美女视频黄频| 久久精品国产a三级三级三级| 国产极品天堂在线| 熟女人妻精品中文字幕| 午夜福利在线观看免费完整高清在| 大片电影免费在线观看免费| 一二三四中文在线观看免费高清| 97超碰精品成人国产| 18禁裸乳无遮挡动漫免费视频| 精品亚洲成a人片在线观看 | 大香蕉97超碰在线| 亚洲av二区三区四区| 国产国拍精品亚洲av在线观看| 新久久久久国产一级毛片| 国产真实伦视频高清在线观看| 丝袜脚勾引网站| 日本免费在线观看一区| 偷拍熟女少妇极品色| 直男gayav资源| 身体一侧抽搐| 插阴视频在线观看视频| 边亲边吃奶的免费视频| 成人二区视频| 色吧在线观看| 欧美变态另类bdsm刘玥| 亚洲aⅴ乱码一区二区在线播放| 国产有黄有色有爽视频| 99热6这里只有精品| 欧美日韩精品成人综合77777| 国产日韩欧美亚洲二区| 欧美少妇被猛烈插入视频| 一个人看视频在线观看www免费| 18禁在线无遮挡免费观看视频| 寂寞人妻少妇视频99o| 日韩一区二区三区影片| 国产伦精品一区二区三区视频9| 国产免费又黄又爽又色| 天堂8中文在线网| 国产欧美亚洲国产| 亚洲欧美日韩卡通动漫| 黄色一级大片看看| 国产av码专区亚洲av| 女人十人毛片免费观看3o分钟| 国产伦在线观看视频一区| 青春草亚洲视频在线观看| 精品一区二区三卡| 日韩中文字幕视频在线看片 | av.在线天堂| 免费观看无遮挡的男女| 中国美白少妇内射xxxbb| 有码 亚洲区| 欧美97在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 小蜜桃在线观看免费完整版高清| 高清av免费在线| 日韩视频在线欧美| 久久久精品94久久精品| 一个人看视频在线观看www免费| 国语对白做爰xxxⅹ性视频网站| 在线观看av片永久免费下载| 亚洲精品国产av蜜桃| 自拍欧美九色日韩亚洲蝌蚪91 | av免费观看日本| 免费黄频网站在线观看国产| 青春草亚洲视频在线观看| 欧美最新免费一区二区三区| 亚洲人与动物交配视频| 亚洲欧美成人综合另类久久久| 男女无遮挡免费网站观看| 视频中文字幕在线观看| 精品国产露脸久久av麻豆| 少妇 在线观看| 少妇丰满av| 青春草视频在线免费观看| 国产综合精华液| 插逼视频在线观看| a级毛色黄片| 99热网站在线观看| 欧美老熟妇乱子伦牲交| 精品久久久精品久久久| 免费人妻精品一区二区三区视频| 精品亚洲成a人片在线观看 | 久久久久久久久久人人人人人人| 午夜激情福利司机影院| 一区在线观看完整版| 国产成人a区在线观看| 搡老乐熟女国产| 久久久久久久久久人人人人人人| 99热6这里只有精品| 国语对白做爰xxxⅹ性视频网站| av天堂中文字幕网| 最后的刺客免费高清国语| 欧美区成人在线视频| 亚洲va在线va天堂va国产| 全区人妻精品视频| 国产精品熟女久久久久浪| 永久免费av网站大全| 国产 精品1| 国产精品国产三级专区第一集| 热re99久久精品国产66热6| 免费观看a级毛片全部| 内地一区二区视频在线| 七月丁香在线播放| 日韩电影二区| 亚洲精品国产色婷婷电影| 亚洲欧美成人精品一区二区| 美女福利国产在线 | 美女国产视频在线观看| 精品亚洲成a人片在线观看 | 黄色欧美视频在线观看| 国产男人的电影天堂91| 人妻 亚洲 视频| 欧美一区二区亚洲| 五月开心婷婷网| 欧美精品亚洲一区二区| 久久久成人免费电影| 国产精品不卡视频一区二区| 亚洲成人手机| 国产成人精品久久久久久| 最近的中文字幕免费完整| 这个男人来自地球电影免费观看 | 精品国产一区二区三区久久久樱花 | 精品久久久久久久末码| 熟女电影av网| 国产人妻一区二区三区在| 欧美亚洲 丝袜 人妻 在线| videos熟女内射| 在线观看美女被高潮喷水网站| 国产免费福利视频在线观看| 亚洲国产成人一精品久久久| 亚洲一级一片aⅴ在线观看| 精品视频人人做人人爽| 偷拍熟女少妇极品色| .国产精品久久| 在现免费观看毛片| 国产极品天堂在线| 亚洲精品456在线播放app| 观看av在线不卡| 久久久久久久久久久丰满| av又黄又爽大尺度在线免费看| 国产伦精品一区二区三区四那| 亚洲av电影在线观看一区二区三区| 亚洲最大成人中文| 亚洲国产最新在线播放| 1000部很黄的大片| 久热久热在线精品观看| 国产亚洲午夜精品一区二区久久| 黑丝袜美女国产一区| 国产无遮挡羞羞视频在线观看| 亚洲av福利一区| 国产成人精品福利久久| 少妇熟女欧美另类| 国产精品人妻久久久影院| 大香蕉久久网| 久久精品国产亚洲网站| 亚洲自偷自拍三级| 大香蕉久久网| 国产精品国产三级国产专区5o| 国产一区亚洲一区在线观看| 欧美激情极品国产一区二区三区 | 又大又黄又爽视频免费| 97超碰精品成人国产| 国内精品宾馆在线| 成人国产麻豆网| 久久av网站| av女优亚洲男人天堂| 少妇人妻久久综合中文| 伊人久久精品亚洲午夜| 最近的中文字幕免费完整| 亚洲精品乱久久久久久| 少妇高潮的动态图| a级毛色黄片| 少妇的逼好多水| 人体艺术视频欧美日本| 一级av片app| 亚洲欧洲日产国产| 日本免费在线观看一区| 久久久久久久国产电影| 亚洲av电影在线观看一区二区三区| 三级国产精品欧美在线观看| 久久久久久久久久久丰满| 校园人妻丝袜中文字幕| 久久精品人妻少妇|