999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

醋酸戊酯差壓熱耦合反應(yīng)精餾過(guò)程的節(jié)能優(yōu)化

2017-04-27 03:14:43張青瑞李永磊
化學(xué)與生物工程 2017年4期

郭 通,張青瑞,李永磊,禹 超

(青島科技大學(xué)化工學(xué)院,山東 青島 266042)

?

醋酸戊酯差壓熱耦合反應(yīng)精餾過(guò)程的節(jié)能優(yōu)化

郭 通,張青瑞*,李永磊,禹 超

(青島科技大學(xué)化工學(xué)院,山東 青島 266042)

為了實(shí)現(xiàn)醋酸戊酯反應(yīng)精餾系統(tǒng)反應(yīng)-分離-能量集成的優(yōu)化設(shè)計(jì),建立了差壓熱耦合反應(yīng)精餾(DPT-RD)新工藝。采用Aspen軟件對(duì)DPT-RD流程進(jìn)行模擬,考察了塔板持液量、降壓分餾塔壓力、HAc進(jìn)料位置、兩塔耦合方式等對(duì)DPT-RD過(guò)程的影響。結(jié)果表明,在塔板持液量為0.25 m3、降壓分餾塔壓力為0.16 bar、HAc第6塊板進(jìn)料、常規(guī)分餾塔塔板數(shù)為24的最優(yōu)條件下,DPT-RD比常規(guī)反應(yīng)精餾節(jié)能約78.22%,節(jié)能效果顯著。

醋酸戊酯;差壓熱耦合反應(yīng)精餾;節(jié)能優(yōu)化

醋酸戊酯(C7H14O2)是一種重要的化工原料,作為溶劑、稀釋劑用于制造香精、化妝品、人造革、膠卷、火藥等[1]。醋酸戊酯主要是由醋酸和戊醇在催化劑的作用下酯化生成,在該反應(yīng)體系中反應(yīng)物與產(chǎn)物之間存在多種共沸物,共沸現(xiàn)象嚴(yán)重,后續(xù)分離復(fù)雜且能耗大。由于精餾塔能同時(shí)進(jìn)行反應(yīng)和分離操作,及時(shí)分離產(chǎn)物,避免共沸物的形成,大大降低能耗,因此醋酸戊酯反應(yīng)精餾技術(shù)得到了快速發(fā)展。

Lee等[2]對(duì)以Dowex50Wx8-100為催化劑的醋酸和戊醇合成醋酸戊酯的酯化反應(yīng)動(dòng)力學(xué)進(jìn)行了研究。Chiang等[3]研究了醋酸戊酯反應(yīng)器與塔耦合、反應(yīng)精餾2種設(shè)計(jì)結(jié)構(gòu),并對(duì)其進(jìn)行了動(dòng)態(tài)模擬。Luyben等[4]對(duì)醋酸戊酯反應(yīng)精餾流程進(jìn)行了模擬和討論。潘晶[5]討論了醋酸戊酯反應(yīng)精餾體系內(nèi)部物質(zhì)耦合,并研究了內(nèi)部物質(zhì)耦合對(duì)其動(dòng)態(tài)特性的影響。Wu等[6-7]對(duì)正戊醇、正己醇混合物酯化反應(yīng)的隔壁塔反應(yīng)精餾及熱耦合反應(yīng)精餾過(guò)程的設(shè)計(jì)與控制進(jìn)行了深入研究。李洪等[8]對(duì)醋酸戊酯反應(yīng)精餾的動(dòng)態(tài)控制過(guò)程進(jìn)行了討論,并對(duì)不同擾動(dòng)下溫度控制和組成控制的效果進(jìn)行了比較。

目前,在反應(yīng)精餾的基礎(chǔ)上已開(kāi)發(fā)出了多種節(jié)能方案[9-11],如隔壁塔反應(yīng)精餾、側(cè)反應(yīng)精餾、熱耦合等,但醋酸戊酯差壓熱耦合反應(yīng)精餾(DPT-RD)的研究還未見(jiàn)報(bào)道。作者在此分析了差壓熱耦合(DPT)技術(shù)應(yīng)用于醋酸戊酯反應(yīng)精餾(RD)過(guò)程的可行性,重點(diǎn)研究了DPT-RD流程中各個(gè)變量對(duì)產(chǎn)品質(zhì)量及節(jié)能效果的影響,并對(duì)其進(jìn)行了優(yōu)化。

1 節(jié)能方案的選定

1.1 醋酸戊酯反應(yīng)精餾系統(tǒng)

醋酸戊酯是由醋酸和戊醇在反應(yīng)精餾塔中酯化生成,反應(yīng)式為:

(1)

該過(guò)程是一個(gè)可逆反應(yīng),體系中存在醋酸(HAc)、戊醇(AmOH)、醋酸戊酯(AmAc)和水(H2O) 4 種組分,這些組分構(gòu)成了強(qiáng)液相非理想體系,用NRTL方程可以很好地描述其熱力學(xué)特性,二元交互參數(shù)如表1所示[3]。

表1 醋酸戊酯體系的二元交互參數(shù)

Tab.1 Binary parameters of n-amyl acetate system

注:bij、bji、cij為NRTL模型中與溫度相關(guān)的二元參數(shù),可通過(guò)VLE或LLE數(shù)據(jù)回歸確定,用于活度的求取。

本實(shí)驗(yàn)動(dòng)力學(xué)模型采用以酸性離子交換樹(shù)脂Amberlyst 15為催化劑的擬均相可逆反應(yīng)的動(dòng)力學(xué),如式(2)~(4)所示:

r=mcat(kFcHAccAmOH-kBcAmAccH2O)

(2)

kF=31.1667exp(-51740/RT)

(3)

kB=2.2533exp(-45280/RT)

(4)

式中:r為反應(yīng)速率,kmol·s-1;kF、kB分別為正、逆向反應(yīng)速率常數(shù),m6·kmol-1·kg-1·s-1;mcat為催化劑的質(zhì)量,kg;cHAc、cAmOH、cAmAc、cH2O分別為醋酸、戊醇、醋酸戊酯和水的濃度,kmol·m-3。

1.2 常規(guī)反應(yīng)精餾過(guò)程

在醋酸戊酯反應(yīng)精餾過(guò)程中,醋酸和戊醇以物質(zhì)的量比1∶1進(jìn)料,戊醇作為較重組分在反應(yīng)段上部進(jìn)料,醋酸在反應(yīng)段中上部進(jìn)料。由于存在共沸物,一部分未反應(yīng)的反應(yīng)物與水一同進(jìn)入塔頂作為氣相采出后進(jìn)入全凝器進(jìn)行冷凝,然后進(jìn)入傾析器進(jìn)行分相;水相全部采出,油相全部回流至反應(yīng)精餾塔內(nèi)進(jìn)行進(jìn)一步反應(yīng)和精餾,產(chǎn)品醋酸戊酯作為重組分在塔底采出。

根據(jù)以上流程,采用Aspen中Radfrac模塊對(duì)該過(guò)程的穩(wěn)態(tài)流程進(jìn)行模擬,基本操作參數(shù)如圖1所示。

圖1 醋酸戊酯常規(guī)反應(yīng)精餾流程圖Fig.1 The flow chart of conventional reactive distillation of n-amyl acetate

1.3 DPT-RD可行性分析

DPT-RD技術(shù)是將反應(yīng)精餾塔分為兩個(gè)塔,通過(guò)壓縮機(jī)和節(jié)流閥連接。兩塔間壓力不同,常規(guī)分餾塔的壓力高于降壓分餾塔的壓力,當(dāng)常規(guī)分餾塔塔頂氣相物料溫度高于降壓分餾塔塔底液相物料并存在足夠的傳熱動(dòng)力時(shí),兩者將進(jìn)行換熱,從而降低常規(guī)分餾塔冷凝器和再沸器的熱負(fù)荷,達(dá)到節(jié)能效果[12-14]。DPT-RD流程如圖2所示。

以常規(guī)反應(yīng)精餾過(guò)程的基本操作參數(shù)為基礎(chǔ),得到較為合適的參比操作條件,并根據(jù)建立的DPT-RD流程模型在Aspen中進(jìn)行模擬,結(jié)果如表2所示。

由表2可知,常規(guī)分餾塔塔頂溫度高于降壓分餾塔塔底溫度,足夠的傳熱溫差使兩股物流有大的傳熱動(dòng)力;塔頂冷凝器熱負(fù)荷大于塔底再沸器熱負(fù)荷,說(shuō)明塔頂氣相物料有足夠的熱量與塔底液相物料進(jìn)行換熱,多余的氣相物料熱量依靠輔助冷凝器實(shí)現(xiàn)進(jìn)一步降溫;換熱后物料出口溫度仍較高,可用常溫水冷凝,降低冷凝費(fèi)用。這表明差壓熱耦合技術(shù)的確可以應(yīng)用于醋酸戊酯反應(yīng)精餾過(guò)程,且較常規(guī)反應(yīng)精餾過(guò)程有一定的節(jié)能效果。

1.常規(guī)分餾塔 2.降壓分餾塔 3.節(jié)流閥 4.壓縮機(jī) 5.主換熱器 6.輔助冷凝器 7.輔助再沸器 8.傾析器圖2 DPT-RD流程圖Fig.2 The flow chart of DPT-RD

表2 DPT-RD過(guò)程參比操作條件和模擬計(jì)算結(jié)果

Tab.2 Reference operating conditions and simulation results of DPT-RD process

2 過(guò)程參數(shù)優(yōu)化及結(jié)果分析

與普通精餾過(guò)程不同,常規(guī)反應(yīng)精餾過(guò)程中各參數(shù)對(duì)流程有很大的影響,差壓熱耦合技術(shù)與反應(yīng)精餾過(guò)程的結(jié)合更增加了流程的復(fù)雜性。以產(chǎn)品純度及流程總能耗為目標(biāo)函數(shù),在固定其它條件不變,分別對(duì)塔板持液量、降壓分餾塔壓力、HAc進(jìn)料位置、兩塔耦合方式等幾個(gè)變量進(jìn)行分析,考察它們對(duì)目標(biāo)函數(shù)的影響并對(duì)其進(jìn)行優(yōu)化。

2.1 塔板持液量、降壓分餾塔壓力的影響

降壓分餾塔壓力對(duì)出口溫度的影響如圖3所示。

圖3 降壓分餾塔壓力對(duì)出口溫度的影響Fig.3 Effect of pressure of pressure-reducing distillation tower on outlet temperature

由圖3可以看出,隨著降壓分餾塔壓力的增大,再沸器的出口溫度不斷升高,在0.19 bar時(shí)與冷凝器出口溫度相同;若繼續(xù)增大壓力,再沸器出口溫度將高于冷凝器出口溫度,此時(shí)將無(wú)法保證差壓熱耦合技術(shù)正常運(yùn)行。為了保證交換器的正常換熱,再沸器與冷凝器間需有一定的溫差(最小溫差為10 ℃),因此,降壓分餾塔壓力不宜超過(guò)0.16 bar。

塔板持液量的大小一定程度上決定了反應(yīng)的進(jìn)行程度。為了考察塔板持液量和降壓分餾塔壓力對(duì)DPT-RD過(guò)程的影響,在其它條件不變的情況下,考察塔板持液量、降壓分餾塔壓力對(duì)產(chǎn)品純度的影響,結(jié)果如圖4所示。

圖4 塔板持液量和降壓分餾塔壓力對(duì)產(chǎn)品純度的影響Fig.4 Effect of tray liquid holdup and pressure of pressure-reducing distillation tower on product purity

由圖4可以看出,三維圖像呈現(xiàn)出一個(gè)光滑的曲面,在塔板持液量一定的情況下降壓分餾塔壓力對(duì)產(chǎn)品純度的影響很小,而持液量的改變卻對(duì)產(chǎn)品純度影響很大。在持液量低于0.20 m3時(shí),產(chǎn)品純度隨持液量的增加迅速上升;持液量高于0.20 m3時(shí)升幅趨緩,在持液量達(dá)到0.22 m3時(shí)產(chǎn)品純度基本可以達(dá)到0.99。考慮到持液量的增加會(huì)增加催化劑用量,導(dǎo)致成本增加,因此,持液量以0.25 m3為宜。

圖5是塔板持液量和降壓分餾塔壓力對(duì)總能耗的影響,由于常規(guī)分餾塔塔頂冷凝器提供的熱量足夠保證降壓分餾塔塔底再沸器的熱量交換,因此總能耗僅由壓縮機(jī)產(chǎn)生。

圖5 塔板持液量和降壓分餾塔壓力對(duì)總能耗的影響Fig.5 Effect of tray liquid holdup and pressure of pressure-reducing distillation tower on total energy consumption

由圖5可以看出,隨著持液量的增加及壓力的增大總能耗均減小,在持液量超過(guò)0.25 m3時(shí)總能耗降幅減緩,且在0.16 bar下總能耗始終處于一個(gè)較低的狀態(tài),因此,降壓分餾塔壓力以0.16 bar為宜。

2.2 HAc進(jìn)料位置、兩塔耦合方式的影響

圖6是HAc進(jìn)料位置和兩塔耦合方式對(duì)產(chǎn)品純度的影響。

圖6 HAc進(jìn)料位置和兩塔耦合方式對(duì)產(chǎn)品純度的影響Fig.6 Effect of HAc feed location and two towers coupling mode on product purity

由圖6可以看出,HAc進(jìn)料位置對(duì)產(chǎn)品純度的影響較大,在高于第5塊板處進(jìn)料時(shí)產(chǎn)品純度迅速下降,這是由于進(jìn)料位置過(guò)高,離塔頂較近,導(dǎo)致反應(yīng)不夠完全;同時(shí)在第5塊板以下的位置進(jìn)料時(shí)產(chǎn)品純度變化不大,在第6、7塊板進(jìn)料時(shí)產(chǎn)品純度最大。相對(duì)于進(jìn)料位置,兩塔耦合方式對(duì)產(chǎn)品純度的影響較小,但隨著常規(guī)分餾塔塔板數(shù)的增加,即反應(yīng)段分配在常規(guī)分餾塔塔板數(shù)的增加,產(chǎn)品純度逐漸升高;當(dāng)常規(guī)分餾塔塔板數(shù)超過(guò)23時(shí),即反應(yīng)段完全分配在常規(guī)分餾塔內(nèi),增加塔板數(shù)產(chǎn)品純度升幅不明顯。在HAc進(jìn)料位置在5~10塊板處、常規(guī)分餾塔塔板數(shù)超過(guò)23時(shí)產(chǎn)品純度基本可達(dá)0.99。

圖7是HAc進(jìn)料位置和兩塔耦合方式對(duì)總能耗的影響。

圖7 HAc進(jìn)料位置和兩塔耦合方式 對(duì)總能耗的影響Fig.7 Effect of HAc feed location and two towers coupling mode on total energy consumption

由圖7可以看出,HAc進(jìn)料位置升高總能耗降低,在高于7塊板時(shí)總能耗變化不大;同時(shí)隨著常規(guī)分餾塔塔板數(shù)的增加,總能耗不斷降低,塔板數(shù)超過(guò)23時(shí)總能耗變化很小基本維持恒定。在保證產(chǎn)品純度達(dá)到一定要求的情況下,選擇較低能耗的操作方式,即進(jìn)料位置為第6塊板處,常規(guī)分餾塔塔板數(shù)為24。

2.3 最優(yōu)條件下的節(jié)能效果分析

除了降壓分餾塔壓力對(duì)產(chǎn)品純度影響很小外,其它變量對(duì)目標(biāo)函數(shù)均影響較大,但由于降壓分餾塔壓力的改變導(dǎo)致再沸器出口溫度發(fā)生變化,從而影響到DPT-RD過(guò)程的進(jìn)行。所以在實(shí)際操作中應(yīng)該綜合考慮各個(gè)變量,以保證流程平穩(wěn)運(yùn)行,實(shí)現(xiàn)高效高產(chǎn)。

綜合過(guò)程參數(shù)對(duì)產(chǎn)品純度和總能耗的影響,得到其優(yōu)化工藝參數(shù)為:持液量0.25 m3、降壓分餾塔壓力0.16 bar、HAc第6塊板進(jìn)料、常規(guī)分餾塔塔板數(shù)24。在該工藝條件下,對(duì)醋酸戊酯常規(guī)反應(yīng)精餾流程及DPT-RD流程進(jìn)行總能耗的對(duì)比,結(jié)果見(jiàn)表3。

表3 兩種反應(yīng)精餾流程的總能耗對(duì)比/kW

Tab.3 Comparison of total energy consumption for two kinds of reactive distillation processes/kW

由表3可知,相對(duì)于常規(guī)反應(yīng)精餾,DPT-RD由于再沸器能耗完全由冷凝器供給,總能耗即壓縮機(jī)的能耗,總能耗減少了762.5 kW,節(jié)能約78.22%。所以,DPT-RD技術(shù)具有明顯的節(jié)能優(yōu)勢(shì)。

3 結(jié)論

通過(guò)對(duì)DPT-RD的可行性分析,發(fā)現(xiàn)差壓熱耦合技術(shù)在醋酸戊酯常規(guī)反應(yīng)精餾過(guò)程的應(yīng)用是可行的。研究了塔板持液量、降壓分餾塔壓力、HAc進(jìn)料位置、兩塔耦合方式等對(duì)DPT-RD過(guò)程的影響,得到最優(yōu)工藝條件為:塔板持液量0.25 m3、降壓分餾塔壓力0.16 bar、HAc第6塊板進(jìn)料、常規(guī)分餾塔塔板數(shù)24。對(duì)比醋酸戊酯常規(guī)反應(yīng)精餾流程及DPT-RD流程的總能耗,發(fā)現(xiàn)差壓熱耦合技術(shù)的應(yīng)用對(duì)降低能耗有很大幫助,節(jié)能約78.22%,節(jié)能效果顯著。

[1] ZHU F H,HUANG K J,WANG S F,et al.Towards further internal heat integration in design of reactive distillation columns.Part Ⅳ:application to a high-purity ethylene glycol reactive distillation column[J].Chemical Engineering Science,2009,64(15):3498-3509.

[2] LEE M J,WU H T,LIN H M.Kinetics of catalytic esterification of acetic acid and amyl alcohol over Dowex[J].Industrial & Engineering Chemistry Research,2000,39(11):4094-4099.

[3] CHIANG S F,KUO C L,YU C C,et al.Design alternatives for the amyl acetate process:coupled reactor/column and reactive distillation[J].Industrial & Engineering Chemistry Research,2002,41(13):3233-3246.

[4] LUYBEN W L,YU C C.Reactive Distillation Design and Control[M].John Wiley & Sons,2008:147-177.

[5] 潘晶.乙酸戊酯反應(yīng)精餾塔的綜合設(shè)計(jì)與控制-強(qiáng)化內(nèi)部物質(zhì)耦合的影響[D].北京:北京化工大學(xué),2011.

[6] WU Y C,LEE H Y,LEE C H,et al.Design and control of reactive divided wall column for esterification with mixedn-amyl alcohol andn-hexanol feed[J].Computer Aided Chemical Engineering,2012,31:1582-1586.

[7] WU Y C,LEE H Y,LEE C H,et al.Design and control of thermally-coupled reactive distillation system for esterification of an alcohol mixture containingn-amyl alcohol andn-hexanol[J].Industrial & Engineering Chemistry Research,2013,52(48):17184-17197.

[8] 李洪,孟瑩,李鑫鋼,等.乙酸戊酯酯化反應(yīng)精餾過(guò)程系統(tǒng)控制模擬及分析[J].化工進(jìn)展,2015,34(12):4165-4171.

[9] LONG N V D,LE Q M,LE C N,et al.A novel self-heat recuperative dividing wall column to maximize energy efficiency and column throughput in retrofitting and debottlenecking of a side stream column[J].Applied Energy,2015,159:28-38.

[10] 安維中,林子昕,江月,等.考慮內(nèi)部熱集成的乙二醇反應(yīng)精餾系統(tǒng)設(shè)計(jì)與優(yōu)化[J].化工學(xué)報(bào),2013,64(12):4634-4640.

[11] 黃玉鑫,湯吉海,陳獻(xiàn),等.不同溫度反應(yīng)與精餾集成生產(chǎn)醋酸叔丁酯的過(guò)程模擬[J].化工學(xué)報(bào),2015,66(10):4039-4046.

[12] IWAKABE K,NAKAIWA M,HUANG K J,et al.Energy saving in multicomponent separation using an internally heat-integrated distillation column (HIDiC)[J].Applied Thermal Engineering,2006,26(13):1362-1368.

[13] 李洪,李鑫鋼,羅銘芳.差壓熱耦合蒸餾節(jié)能技術(shù)[J].化工進(jìn)展,2008,27(7):1125-1128.

[14] 孫蘭義,寧亞南,劉雪暖,等.反應(yīng)精餾應(yīng)用于差壓熱耦合精餾技術(shù)的模擬[J].化工進(jìn)展,2011,30(S2):149-153.

Energy-Saving Optimization on Differential Pressure Thermally Coupled Reactive Distillation Process ofn-Amyl Acetate

GUO Tong,ZHANG Qing-rui*,LI Yong-lei,YU Chao

(CollegeofChemicalEngineering,QingdaoUniversityofScienceandTechnology,Qingdao266042,China)

Inordertorealizetheoptimaldesignofreaction-separation-energyintegrationforareactivedistillationsystemofn-amylacetate,anovelprocessofdifferentialpressurethermallycoupledreactivedistillation(DPT-RD)wasestablished.TheDPT-RDprocesswassimulatedwithAspensoftware,andtheeffectsoftheparametersonDPT-RDprocesswereinvestigated,includingthetrayliquidholdup,thepressureofpressure-reducingdistillationtower,theHAcfeedlocation,andthetwotowerscouplingmode.Resultsshowedthat,thetotalenergyconsumptionofDPT-RDprocessreducedby78.22%comparedwiththeconventionalreactivedistillationundertheoptimalconditionsasfollows:thetrayliquidholdupwas0.25m3,thepressureofpressure-reducingdistillationtowerwas0.16bar,theHAcfeedlocationwasthesixthplate,andtheplatenumberofatmosphericdistillationtowerwas24.

n-amylacetate;differentialpressurethermallycoupledreactivedistillation;energy-savingoptimization

山東省自然科學(xué)基金資助項(xiàng)目(ZR2013BM001)

2016-10-12

郭通(1993-),男,河南新鄉(xiāng)人,碩士研究生,研究方向:過(guò)程系統(tǒng)工程,E-mail:2207190830@qq.com;通訊作者:張青瑞,副教授,博士,E-mail:zhangqingrui@qust.edu.cn。

10.3969/j.issn.1672-5425.2017.04.014

TQ028.31 TQ051.81

A

1672-5425(2017)04-0058-05

郭通,張青瑞,李永磊,等.醋酸戊酯差壓熱耦合反應(yīng)精餾過(guò)程的節(jié)能優(yōu)化[J].化學(xué)與生物工程,2017,34(4):58-62.

主站蜘蛛池模板: 亚洲精品国产乱码不卡| 国产菊爆视频在线观看| 色色中文字幕| 日本午夜影院| 四虎永久在线| 人妻中文字幕无码久久一区| 成人国产精品网站在线看| 精品91视频| 岛国精品一区免费视频在线观看 | 一级做a爰片久久免费| 最新加勒比隔壁人妻| 久久99国产视频| 色欲色欲久久综合网| 一级高清毛片免费a级高清毛片| 国产精品开放后亚洲| 伊人丁香五月天久久综合| 国产一级精品毛片基地| 干中文字幕| 五月天久久综合国产一区二区| 国产成人av一区二区三区| 国产成年无码AⅤ片在线| www.国产福利| 亚洲欧美成人| 亚洲av无码久久无遮挡| 沈阳少妇高潮在线| 亚洲国产看片基地久久1024 | 亚洲综合色婷婷| 午夜无码一区二区三区| 国内自拍久第一页| 老司机aⅴ在线精品导航| 欧美激情第一欧美在线| 国产精品爽爽va在线无码观看| 久久黄色毛片| 国产高清在线丝袜精品一区 | 欧美有码在线观看| 欧美日韩91| 亚洲欧美在线综合一区二区三区| 在线国产91| 四虎永久免费网站| 国产欧美性爱网| 亚洲三级a| 亚洲成年人片| 色婷婷成人网| 国产真实乱子伦视频播放| 青青青伊人色综合久久| 99这里只有精品免费视频| 玩两个丰满老熟女久久网| 91精品日韩人妻无码久久| 内射人妻无码色AV天堂| 欧美在线中文字幕| 亚亚洲乱码一二三四区| 狠狠色狠狠综合久久| 亚洲乱伦视频| 国产乱人伦偷精品视频AAA| 色欲色欲久久综合网| 亚洲国产精品久久久久秋霞影院| 免费无码AV片在线观看国产| www.99在线观看| 亚洲天堂网2014| 国产精品人人做人人爽人人添| 九色在线观看视频| 久久一色本道亚洲| 亚洲无码在线午夜电影| 欧美色丁香| 秋霞一区二区三区| 视频一区亚洲| 日韩高清中文字幕| 欧洲欧美人成免费全部视频| 欧美日韩国产精品综合| 亚洲免费人成影院| 91久久青青草原精品国产| 色综合中文| 日本www在线视频| 秘书高跟黑色丝袜国产91在线| 凹凸国产熟女精品视频| 色妞www精品视频一级下载| 狠狠亚洲婷婷综合色香| 国产欧美日韩18| 国产一区二区精品福利| 亚洲最大看欧美片网站地址| 国产在线精品美女观看| 国产欧美日韩专区发布|