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

微波改性對活性炭及其甲醇吸附的影響

2014-11-28 17:41:57李立清梁鑫姚小龍李海龍馬衛武
湖南大學學報·自然科學版 2014年7期

李立清+梁鑫+姚小龍+李海龍+馬衛武+劉颯

摘要:分別在600,700和800 ℃下對活性炭進行微波輻照加熱改性.采用比表面積及孔徑分析儀、Boehm滴定、傅立葉轉換紅外光譜對活性炭的物化性質進行表征.并且在10 ℃下以甲醇為吸附質進行固定床吸附實驗.研究表明:微波改性后,活性炭的比表面積、總孔容小幅度減小,但微孔比表面積顯著增大;隨著溫度升高,活性炭表面酸性基團大量分解,堿性基團逐漸形成.Langmuir方程和Freundlich方程均能較好的描述甲醇在活性炭上的吸附.準二階動力學方程最適合描述甲醇的動態吸附過程,說明甲醇吸附是一個物理和化學復合的吸附過程,吸附受到活性炭表面官能團的影響.顆粒內擴散模型擬合結果分為3個線性階段:表面吸附階段、漸近吸附階段和吸附平衡階段.微波改性后活性炭對甲醇的吸附能均增大,吸附能與活性炭表面含氮官能團總量成正比.

關鍵詞:微波;改性活性炭;吸附;甲醇

中圖分類號:O613.71 文獻標識碼:A

甲醇是常見于農藥、醫藥、涂料、染料等有機化工產業產生的揮發性有機化合物(VOCs)之一,它對人體具有強烈的毒性,麻醉中樞神經系統引起急慢性中毒.由于其極易揮發,所以對空氣中存在的甲醇進行去除顯得非常必要.活性炭(AC)是一種具有巨大比表面積的多孔吸附材料,活性炭吸附也越來越廣泛地用于VOCs的去除\[1-2\].活性炭的表面化學性質及孔結構決定其吸附性能,通過各種途徑對活性炭進行改性,研究改性后活性炭的吸附性能一直是研究的熱點.

熱處理能改變活性炭的表面化學特性,同時也能有效調節活性炭的孔結構\[3\].微波加熱屬于內加熱,升溫迅速,微波場中無溫度梯度存在,受熱均勻,加熱效率高,能耗低,在活性炭改性研究中得到重視.有學者對比研究了微波加熱和傳統加熱對活性炭吸附行為的影響,發現微波改性活性炭具有更大的微孔孔容\[4\],基于微波加熱的優越性,本文對不同溫度微波改性活性炭吸附甲醇進行了實驗研究,多角度探討了微波改性對活性炭本身及其吸附甲醇的影響.

1實驗部分

1.1熱重實驗

利用熱重分析儀(LABSYS TG, SETARAM, France)在氮氣氛圍下測試原始未改性活性炭的TG特性曲線.

1.2活性炭改性

稱取50 g活性炭(RS5型,河南長葛利民活性炭有限公司)進行預處理:將活性炭置于燒杯中,采用電子萬用爐用水煮沸30 min(煮沸);煮沸后的活性炭用去離子水洗滌4~5次(漂洗);漂洗后的活性炭置真空干燥箱(DZF型,北京市永光明醫療儀器廠)中在110 °C下干燥12 h(烘干).這樣所得樣品為原始未改性活性炭,記為AC1.稱取10 g AC1置于微波高溫燒結爐(WZ3/2.45型,長沙隆泰科技有限公司)中,在流量為600 mL·min-1的高純N2氣氛下進行微波輻照,微波頻率為2.45 GHz,根據熱重實驗結果,分別選定在600,700和800 ℃ 3個溫度點下加熱30 min,冷卻后置于真空干燥箱中110 ℃下干燥24 h,所得樣品分別記為AC600,AC700和AC800.

1.3活性炭表征

1.3.1孔結構表征

活性炭的孔結構參數采用低溫氮氣吸附法進行測定.利用比表面積及孔徑分析儀SA3100 (BECKMAN COULTER, USA) 測定77 K下高純N2在活性炭上的吸附等溫線.BET比表面積SBET由標準BET法得到;微孔孔容Vmicro,微孔比表面積Smicro由tPlot方法計算得到;中孔孔容Vmeso,大孔孔容Vmacro基于BJH法計算得到,總孔容Vtotal由相對壓力為0.9814時的液氮吸附量換算成液氮體積得到;孔徑分布是基于BJH方法通過分析吸附或脫附等溫線分支得到.

1.3.2表面官能團表征

利用Boehm滴定定量測試活性炭表面總酸性基團、總堿性基團、羧基以及酚羥基的含量;采用傅立葉轉換紅外光譜儀NEXUS670 (Nicolet, USA) 表征活性炭表面某些特定結構的官能團.

1.4吸附實驗

固定床吸附實驗裝置如圖1所示,該裝置由配氣系統、恒溫系統、吸附床和測試系統組成.室內空氣經過硅膠干燥器干燥后,一部分直接進入混合器,另一部分穿過微型噴淋區和恒溫區得到飽和有機蒸汽,干燥空氣和飽和有機蒸汽在混合器中混合,混合氣體用真空泵吸取并泵入固定床進行吸附,尾氣經凈化后排放.吸附柱進氣濃度和出氣濃度由氣相色譜儀(SP6890型,山東魯南瑞虹化工儀器公司)測定,待出氣濃度與進氣濃度相等,并保持30 min左右,則認為吸附已達平衡\[5-6\].每次實驗中,將4 g活性炭樣置于內徑為1.1 cm石英管吸附柱,吸附溫度用恒溫水箱(DC1015型,上海天平儀器公司)保持在10 °C.分別在C1為6.11,C2為11.03, C3為13.69, C4為18.91和 C5為24.82 g·m-3 5種進氣濃度條件下,利用電子天平稱重法(JA1203N型,上海精密科學儀器有限公司)測得不同濃度下4種活性炭對甲醇(分析純 AR, 國藥集團化學試劑有限公司)的吸附量,從而得到甲醇的等溫吸附曲線.

① 干燥器;② 有機氣體發生器;③ 流量計;

④ 混合器;⑤ 真空泵;⑥ 固定吸附床;⑦ 恒溫水箱;

⑧ 氣相色譜儀;⑨ 尾氣處理

2結果與討論

2.1熱重分析

AC1的熱重曲線如圖2所示.曲線在50~100 °C范圍內呈現輕微的下降趨勢,這可能是活性炭中少量水分蒸發引起;在100~500 °C范圍時,隨溫度升高活性炭質量基本保持不變,這說明該階段活性炭內部的基團保持相對穩定;而在500~900 °C范圍時,曲線以約0.01 mg·°C-1的速度下滑,這說明此溫度范圍內,活性炭表面某些物質或者基團發生了熱分解.因此本文選定在600,700和800°C 3個溫度點下對活性炭進行微波加熱改性.

溫度/℃

2.2改性對孔結構的影響

表1為活性炭的孔結構參數結果.由表可知,AC1的BET比表面積、總孔容均為最大,分別為971.74 m2·g-1和0.517 5 mL·g-1.微波改性后,BET比表面積和總孔容隨著溫度的增加先減小后增大,其大小順序為:AC1>AC800>AC600>AC700,但是微孔比表面積的大小順序與之剛好相反,AC700具有最大的微孔比表面積,比AC1大156.27 m2·g-1.而且,與AC1相比,AC700微孔孔容占總孔容的比例Vmicro/Vtotal由72 %升至83 %,這些都說明,在600~700 ℃范圍內,形成了很多新的微孔\[7\].當溫度升至800 ℃時,微孔比表面積和微孔孔容又降低,可能是由于一部分較小微孔的堵塞,同時,中孔孔容和大孔孔容也有一定程度的降低,這可能是由于某些炭骨架的塌陷所引起.

2.3改性對表面官能團的影響

2.3.1Boehm滴定

表2為Boehm滴定結果.微波改性后,酚羥基、羧基的數量逐漸減小,溫度越高,酸性基團分解量越大,AC600和AC700酸性基團總量分別為AC1的40.8 %和17.7 %,溫度達800 ℃時,酸性基團完全分解.隨著溫度的升高,含氧官能團從活性炭表面去除之后,活性炭的堿性基團則逐漸生成\[8\],堿性基團總量的大小順序為:AC800>AC700>AC600>AC1,AC800的堿性基團總量約為AC1的2.6倍.

2.3.2傅立葉轉換紅外光譜(FTIR)

圖3為FTIR的測試結果.在某些相同波數段,改性活性炭和原始活性炭均出現明顯的吸收峰,同時,改性活性炭在不同波數段出現新的吸收峰,這表明,改性后既保留了原有的某些基團,又產生了新的基團.

由圖所示,2 845~2 975 cm-1處的吸收峰,表明存在C-H鍵的伸縮振動;AC1,AC600和AC700在1 620~1 670 cm-1處的吸收峰,則表明存在碳碳雙鍵的伸縮振動,而AC800無此峰,這說明800 ℃下碳碳雙鍵已經被破壞;AC1,AC600和AC700在3 200~3 670 cm-1處寬而尖銳的吸收峰,表明存在羥基的伸縮振動,而AC800在此波數段表現為雙峰,說明羥基在800 ℃下已完全分解,3 200和3 350 cm-1處的雙峰顯示伯酰胺氫鍵的締合,由N-H的伸縮振動引起.1 030~1 230 cm-1處的吸收峰表明存在脂肪胺C-N的伸縮振動.隨著溫度的升高,1 680~1 750 cm-1波數段的峰強明顯減小,表明活性炭表面羰基C=O,羧基COOH,酯鍵COO的數量逐漸減小.

2.4改性對等溫吸附線的影響

等溫吸附線的測量方法詳見1.4,實驗得到4種活性炭的等溫吸附線如圖4所示.根據國際純粹與應用化學聯合會(IUPAC) 分類法,其均屬于Ⅰ型等溫吸附線.微波改性后,吸附量由大到小順序為AC800>AC600>AC1>AC700.

相對壓力P/P0

由表3可知,Langmuir方程擬合的相關系數R2均大于0.991,表明甲醇在活性炭上的吸附等溫線能較好地用Langmuir模型來預測, Freundlich方程擬合的相關系數相對較小,說明吸附并不是發生在非均勻固體表面.據文獻\[9\],一般認為0.11.0時,吸附作用力變弱,不利于吸附.此處nf都位于0.22~0.31區間,并且微波改性炭的nf值均比AC1小.nf還是吸附線性偏差的一個度量,用來驗證吸附的類型,一般認為nf等于1.0時,吸附是線性的,如果nf小于1.0,則說明吸附是一個物理吸附,如果nf大于1.0,則說明吸附是一個化學吸附.這里最大的nf值為0.308 9,表明甲醇在4種活性炭上均以物理吸附為主,吸附時起主要作用的是范德華力而不是化學鍵力\[10\].

2.5改性對吸附動力學的影響

甲醇在4種活性炭上的動態吸附曲線如圖5所示.由圖可知,隨著吸附時間增大,活性炭的吸附速率減小,前20 min內,吸附速率約為6 mg·min-1,而后40 min內,吸附速率降至0.4 mg·min-1.

時間/min

圖5甲醇的動態吸附曲線

Fig.5The dynamic adsorption curves of methanol

為了全面研究4種活性炭對甲醇的吸附動力學特性,找到最適合描述此吸附過程的動力學模型,本文選用3種動力學模型來對圖5的數據進行擬合.準一階模型主要用于描述物理吸附過程,其動力學方程如下\[11\]:

ln(qe,exp-qt) = ln qe,cal-k1t .(3)

其中:qt和qe,exp分別是吸附時刻t和吸附平衡時刻的吸附量,mg·g-1;qe,cal為模型預測的平衡吸附量,mg·g-1;k1為準一階模型吸附速率常數,min-1.準二階模型主要用于描述物理和化學復合的吸附過程,其動力學方程如式(4)\[12\]:

t/qt=1/k2q2e,cal+t/qe,cal .(4)

其中:qt和qe,cal同上;k2為準二階模型吸附速率常數,g·mg-1·min-1.對圖5數據作ln(qe,exp-qt)~t以及t/qt~t的線性關系圖,通過計算直線的斜率及截距得到的擬合結果見表4.

由表4可知,4種活性炭吸附甲醇的過程用準一階模型和準二階模型擬合的線性相關系數均大于0.98,說明兩種模型均能很好地模擬甲醇的吸附動力學過程.但是,準二階動力學方程擬合的線性相關系數均比準一階動力學方程大,又說明4種活性炭吸附甲醇的過程是一個物理和化學復合的吸附過程,吸附受到活性炭表面官能團的影響,但是官能團的影響在吸附過程中不起主導作用.并且,對于AC600, AC700和AC800,兩種模型預測的平衡吸附量qe,cal的大小順序均與實驗得出的平衡吸附量qe,exp大小順序相吻合.

其中:qt為t時刻的吸附量,mg·g-1;KW為顆粒內擴散速率常數,mg·g-1·min-0.5;t為吸附時間,min;b為直線截距,其數值大小反應了邊界層效應.Weber和Morris認為,如果吸附量qt與t0.5所作直線通過原點,則表明吸附速率由內擴散控制;若不通過原點,截距b越大,膜擴散在速率控制中的影響也越大.利用圖5數據作qt~t0.5的線性關系圖,擬合曲線見圖6,擬合參數見表5.

t0.5/min0.5

圖6 顆粒內擴散模型擬合曲線

Fig.6Fitting curves of intraparticle diffusion model

由上圖可知,qt對t0.5呈現3個線性階段.第1階段為吸附質分子在吸附劑上的表面吸附階段,此階段在前2 min內迅速完成,吸附速率常數KW1較小,吸附速率主要受吸附質濃度的影響.第2階段為漸近吸附階段,此階段約在2~40 min左右完成,該階段吸附速率常數KW2的大小順序為:AC800>AC700>AC600>AC1,且截距b2不為0,說明吸附速率受顆粒內擴散和膜擴散的共同影響.第3階段為最終吸附平衡階段,40 min后活性炭的大部分的微孔被吸附質填充,接著開始在中孔和大孔表面形成多層吸附,吸附速率相比第2階段明顯減小,該階段b3值較大,說明膜擴散在此階段起到了重要作用.由表5還可知,3個階段的吸附速率常數由大到小為KW2,KW3,KW1,吸附質通過顆粒內擴散到吸附質內表面被吸附時的速率常數最大.3個階段線性擬合的相關系數均較大,擬合效果較好,甲醇在活性炭上的整個吸附過程受到顆粒內擴散和膜擴散的共同作用.

3結論

1) 微波輻照加熱改性后,活性炭的總比表面積、孔容較小幅度減小,但微孔比表面積顯著增大,形成許多新的微孔.隨著溫度升高,酸性基團大量分解,堿性基團逐漸形成.

2) 4種活性炭對甲醇的吸附量大小順序為AC800>AC600>AC1>AC700,說明最佳微波改性溫度為800 °C.

3) Langmuir方程和Freundlich方程均能較好的描述甲醇在活性炭上的吸附.

4) 準二階動力學方程最適合描述甲醇的動態吸附過程,說明甲醇吸附是一個物理和化學復合的吸附過程,吸附受到活性炭表面官能團的影響.顆粒內擴散模型擬合結果分為3個線性階段:表面吸附階段、漸近吸附階段和吸附平衡階段.甲醇在活性炭上的吸附受到顆粒內擴散和膜擴散的共同作用.

5) 微波改性后活性炭對甲醇的吸附能均增大,吸附能與活性炭表面含氮官能團總量成正比.

參考文獻

[1]CHIANG YuChun, CHIANG PenChi, HUANG ChinPao. Effects of pore structure and temperature on VOC adsorption on activated carbon\[J\]. Carbon, 2001, 39(4): 523-534.

\[2\]RAMOS M E, BONELLI P R, CUKIERMAN A L, et al. Adsorption of volatile organic compounds onto activated carbon cloths derived from a novel regenerated cellulosic precursor\[J\]. Journal of Hazardous Materials, 2010, 177(1/3): 175-182.

\[3\]MENENDEZ J A, MENENDEZ E M, IGLESIAS M J, et al. Modification of the surface chemistry of active carbons by means of microwaveinduced treatments\[J\]. Carbon, 1999, 37(7): 1115-1121.

\[4\]HUANG Lihui, SUN Yuanyuan, WANG Weiliang, et al. Comparative study on characterization of activated carbons prepared by microwave and conventional heating methods and application in removal of oxytetracycline (OTC)\[J\]. Chemical Engineering Journal, 2011, 171(3): 1446-1453.

\[5\]李立清,HansJürgen Rhm,WolfDieter Suβebach. 活性炭C40/4吸附甲苯和丙酮吸附性能研究\[J\]. 湖南大學學報:自然科學版,2003,30(5):47-50.

LI Liqing, HansJürgen Rhm, WolfDieter Suβebach. Study on adsorption isotherm of toluene and acetone onto activated carbon C40/4\[J\]. Journal of Hunan University: Natural Sciences, 2003, 30(5): 47-50. (In Chinese)

\[6\]李立清,石瑞,顧慶偉,等. 酸改性活性炭吸附甲苯的性能研究\[J\]. 湖南大學學報:自然科學版,2013,40(5):92-98.

LI Liqing, SHI Rui, GU Qingwei et al. Adsorption of toluene on activated carbons with acid modified\[J\]. Journal of Hunan University:Natural Sciences, 2013, 40(5): 92-98. (In Chinese)

\[7\]CHINGOMBE P, SAHA B, WAKEMAN R J. Surface modification and characterisation of a coalbased activated carbon\[J\]. Carbon, 2005, 43(15): 3132-3143.

\[8\]CARROTT P J M, NABAIS J M V, CARROTT M M L R, et al. Thermal treatments of activated carbon fibres using a microwave furnace\[J\]. Microporous and Mesoporous Materials, 2001, 47(2/3): 243-252.

\[9\]GERELI G, SEKI Y, KUSOGLU I M, et al. Equilibrium and kinetics for the sorption of promethazine hydrochloride onto K10 montmorillonite\[J\]. Journal of Colloid and Interface Science, 2006, 299(1): 155-162.

\[10\]VIMONSES V, LEI Shaomin, JIN B, et al. Kinetic study and equilibrium isotherm analysis of Congo Red adsorption by clay materials\[J\]. Chemical Engineering Journal, 2009, 148(2/3): 354-364.

\[11\]HAMEED B H, AHMAD A A, AZIZ N. Isotherms, kinetics and thermodynamics of acid dye adsorption on activated palm ash\[J\]. Chemical Engineering Journal, 2007,133(1/3):195-203.

\[12\]HO Y S, MCKAY G. Pseudosecondorder model for sorption process\[J\]. Process Biochemistry, 1999, 34(5): 451–465.

\[13\]WEBER W J, MORRIS J C. Kinetics of adsorption on carbon from solution\[J\]. J Sanit Eng Div Am Soc Civ Eng, 1963, 89(17): 31-60.

\[14\]CHEN S G, YANG R T. Theoretical basis for the potential theory adsorption isotherms. The DubininRadushkevich and DubininAstakhov equations\[J\]. Langmuir, 1994, 10(11): 4244-4249.

\[15\]WOOD G O. Affinity coefficients of the Polanyi/Dubinin adsorption isotherm equations: A review with compilations and correlations\[J\]. Carbon, 2001, 39(3): 343-356.

\[16\]MANGUN C L, DEBARR J A, ECONOMY J. Adsorption of sulfur dioxide on ammoniatreated activated carbon fibers\[J\]. Carbon, 2001, 39(11): 1689-1696.

\[7\]CHINGOMBE P, SAHA B, WAKEMAN R J. Surface modification and characterisation of a coalbased activated carbon\[J\]. Carbon, 2005, 43(15): 3132-3143.

\[8\]CARROTT P J M, NABAIS J M V, CARROTT M M L R, et al. Thermal treatments of activated carbon fibres using a microwave furnace\[J\]. Microporous and Mesoporous Materials, 2001, 47(2/3): 243-252.

\[9\]GERELI G, SEKI Y, KUSOGLU I M, et al. Equilibrium and kinetics for the sorption of promethazine hydrochloride onto K10 montmorillonite\[J\]. Journal of Colloid and Interface Science, 2006, 299(1): 155-162.

\[10\]VIMONSES V, LEI Shaomin, JIN B, et al. Kinetic study and equilibrium isotherm analysis of Congo Red adsorption by clay materials\[J\]. Chemical Engineering Journal, 2009, 148(2/3): 354-364.

\[11\]HAMEED B H, AHMAD A A, AZIZ N. Isotherms, kinetics and thermodynamics of acid dye adsorption on activated palm ash\[J\]. Chemical Engineering Journal, 2007,133(1/3):195-203.

\[12\]HO Y S, MCKAY G. Pseudosecondorder model for sorption process\[J\]. Process Biochemistry, 1999, 34(5): 451–465.

\[13\]WEBER W J, MORRIS J C. Kinetics of adsorption on carbon from solution\[J\]. J Sanit Eng Div Am Soc Civ Eng, 1963, 89(17): 31-60.

\[14\]CHEN S G, YANG R T. Theoretical basis for the potential theory adsorption isotherms. The DubininRadushkevich and DubininAstakhov equations\[J\]. Langmuir, 1994, 10(11): 4244-4249.

\[15\]WOOD G O. Affinity coefficients of the Polanyi/Dubinin adsorption isotherm equations: A review with compilations and correlations\[J\]. Carbon, 2001, 39(3): 343-356.

\[16\]MANGUN C L, DEBARR J A, ECONOMY J. Adsorption of sulfur dioxide on ammoniatreated activated carbon fibers\[J\]. Carbon, 2001, 39(11): 1689-1696.

\[7\]CHINGOMBE P, SAHA B, WAKEMAN R J. Surface modification and characterisation of a coalbased activated carbon\[J\]. Carbon, 2005, 43(15): 3132-3143.

\[8\]CARROTT P J M, NABAIS J M V, CARROTT M M L R, et al. Thermal treatments of activated carbon fibres using a microwave furnace\[J\]. Microporous and Mesoporous Materials, 2001, 47(2/3): 243-252.

\[9\]GERELI G, SEKI Y, KUSOGLU I M, et al. Equilibrium and kinetics for the sorption of promethazine hydrochloride onto K10 montmorillonite\[J\]. Journal of Colloid and Interface Science, 2006, 299(1): 155-162.

\[10\]VIMONSES V, LEI Shaomin, JIN B, et al. Kinetic study and equilibrium isotherm analysis of Congo Red adsorption by clay materials\[J\]. Chemical Engineering Journal, 2009, 148(2/3): 354-364.

\[11\]HAMEED B H, AHMAD A A, AZIZ N. Isotherms, kinetics and thermodynamics of acid dye adsorption on activated palm ash\[J\]. Chemical Engineering Journal, 2007,133(1/3):195-203.

\[12\]HO Y S, MCKAY G. Pseudosecondorder model for sorption process\[J\]. Process Biochemistry, 1999, 34(5): 451–465.

\[13\]WEBER W J, MORRIS J C. Kinetics of adsorption on carbon from solution\[J\]. J Sanit Eng Div Am Soc Civ Eng, 1963, 89(17): 31-60.

\[14\]CHEN S G, YANG R T. Theoretical basis for the potential theory adsorption isotherms. The DubininRadushkevich and DubininAstakhov equations\[J\]. Langmuir, 1994, 10(11): 4244-4249.

\[15\]WOOD G O. Affinity coefficients of the Polanyi/Dubinin adsorption isotherm equations: A review with compilations and correlations\[J\]. Carbon, 2001, 39(3): 343-356.

\[16\]MANGUN C L, DEBARR J A, ECONOMY J. Adsorption of sulfur dioxide on ammoniatreated activated carbon fibers\[J\]. Carbon, 2001, 39(11): 1689-1696.

主站蜘蛛池模板: 在线免费看黄的网站| 激情午夜婷婷| 国产精品久久久久久久久久98| 国产精品jizz在线观看软件| 久久久久亚洲Av片无码观看| 五月天久久婷婷| 亚洲一区第一页| 免费女人18毛片a级毛片视频| 91热爆在线| 亚洲国产日韩视频观看| 国产精品女熟高潮视频| 日韩a级片视频| 好久久免费视频高清| 香蕉在线视频网站| 天堂网国产| 午夜欧美在线| 亚洲视频一区| 国产精品妖精视频| 色精品视频| 国产素人在线| 日韩欧美中文字幕在线精品| 人人妻人人澡人人爽欧美一区| 亚洲精品视频免费观看| 97超碰精品成人国产| 综合天天色| 人妖无码第一页| 国产精品吹潮在线观看中文| 久久熟女AV| 国产情侣一区| 亚洲一区二区成人| www.精品国产| 国产网友愉拍精品视频| 日韩国产 在线| 婷婷成人综合| 国产三级a| 人妻丰满熟妇啪啪| 国产在线观看一区精品| 国产女人18水真多毛片18精品| 国产成人91精品免费网址在线 | 自慰网址在线观看| 亚洲人成人无码www| 亚洲二区视频| 欧美日韩亚洲国产| 午夜国产理论| 日韩黄色在线| аv天堂最新中文在线| 中文字幕亚洲另类天堂| 91精品国产一区自在线拍| 五月婷婷中文字幕| 国产精品手机视频| 成年人福利视频| 亚洲男人在线天堂| 精品国产欧美精品v| 免费播放毛片| 免费av一区二区三区在线| 成人中文字幕在线| 在线日韩一区二区| 日本高清在线看免费观看| 亚洲一区二区三区香蕉| 欧美国产视频| 欧美自拍另类欧美综合图区| 青青草国产一区二区三区| 91国内在线观看| 亚洲国产中文在线二区三区免| 欧美一级专区免费大片| a在线观看免费| 日本成人在线不卡视频| 91原创视频在线| 99久久精品免费视频| 国产日韩丝袜一二三区| 欧美日韩在线亚洲国产人| 中文字幕丝袜一区二区| 日韩在线播放中文字幕| av性天堂网| 亚洲欧洲自拍拍偷午夜色| www.亚洲一区| 影音先锋亚洲无码| 亚洲V日韩V无码一区二区| 色天天综合久久久久综合片| 亚洲美女一级毛片| 欲色天天综合网| 色成人亚洲|