亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        不同粒徑Pd/Al2O3催化乙炔加氫反應微觀動力學分析

        2016-10-14 00:06:16龐洪強隋志軍朱貽安周興貴
        化工學報 2016年9期
        關鍵詞:基元乙炔微觀

        龐洪強,隋志軍,朱貽安,周興貴

        ?

        不同粒徑Pd/Al2O3催化乙炔加氫反應微觀動力學分析

        龐洪強,隋志軍,朱貽安,周興貴

        (華東理工大學化學工程聯(lián)合國家重點實驗室,上海 200237)

        采用等量浸漬法制備了α-Al2O3負載的系列Pd催化劑,運用BET、XRD、ICP-AES、CO化學吸附、TEM等手段對催化劑進行了表征;根據(jù)部分析因實驗設計方案進行動力學實驗,采用微觀反應動力學方法模擬和分析了所獲穩(wěn)定期本征動力學實驗結果。結果發(fā)現(xiàn),制備所得催化劑Pd顆粒的平均粒徑分別為1.6、3.4、5.5 nm,CO化學吸附所測定達到活性穩(wěn)定期后的催化劑表面Pd原子數(shù)與Hardeveld模型計算的Pd(111) 表面原子數(shù)一致;模擬結果表明該微觀動力學模型可以很好地模擬不同粒徑催化劑上的動力學結果,在所研究范圍內表面最豐物種為C2H4*和C2H3*,通過微觀與宏觀動力學的特征判斷3種催化劑上乙炔加氫的速率控制步驟為乙烯基加氫生成乙烯。

        多相反應;乙炔加氫;微觀動力學模型;粒徑效應;催化劑

        引 言

        蒸汽裂解烴類原料生產(chǎn)聚合級乙烯是重要的石油化工過程之一,是乙烯生產(chǎn)最重要的方法[1]。通過這種方法生產(chǎn)的乙烯中含有0.5%~2%的乙炔,會導致后續(xù)聚合過程的催化劑中毒,因此必須將乙烯產(chǎn)品中的乙炔含量降低至10-6量級[2-4]。目前,工業(yè)上去除乙炔應用最廣泛的方法是在負載Pd系催化劑上進行乙炔選擇性加氫[5-7]。

        由于乙烯工業(yè)的重要性,乙炔選擇性加氫反應已得到廣泛研究[8-9]。近年來量化計算的應用加深了人們對于這一反應的認識。Norskov等[10]針對乙炔加氫反應成功發(fā)展了基于DFT計算和微觀反應動力學分析的催化劑理性設計方法;Hu等[11]系統(tǒng)模擬研究了Pd催化劑表面C和H等物種對加氫、表面聚合等反應的影響。

        以上量化計算研究積累了大量乙炔加氫基元步驟動力學數(shù)據(jù),有利于進行理想催化劑表面上這一反應的微觀動力學分析,但是這些數(shù)據(jù)是否適用于實際催化劑表面仍需進一步檢驗。對于乙炔加氫反應而言,通過微觀反應動力學分析建立動力學模型具有顯著的優(yōu)勢。由于在反應過程中乙炔含量的變化會跨越4個數(shù)量級(由10-2降低到10-6量級),對于傳統(tǒng)反應動力學模型的建立方法而言這意味著極大的實驗工作量和模型參數(shù)回歸,而微觀反應動力學模擬時無須使用速率控制步驟等假設,因此模型具有更為廣泛的適用性。理論上微觀反應動力學的模擬結果不僅可以在實驗范圍以外外推使用,更可以外推至其他催化體系使用[12]。

        在王新宇等[13]的研究中,通過微觀反應動力學模擬發(fā)現(xiàn),Norskov等[10]的計算結果經(jīng)過熱力學一致性校驗和不超過0.1 eV的調整可以與實驗值符合良好。在本研究中,將進一步檢驗所建微觀反應動力學模型描述不同粒徑Pd催化乙炔加氫反應的有效性和適用性。

        1 實驗部分

        1.1 催化劑的制備

        (1)載體的制備 稱取22 g擬薄水鋁石和0.22 g硝酸鑭,加入80 ml去離子水,磁力攪拌條件下逐滴加入濃硝酸直至懸浮液變?yōu)槿榘咨该魅苣z。將制得的溶膠在80℃下干燥12 h后研磨過篩,將樣品置于馬弗爐中1200℃焙燒4.5 h,制得α-Al2O3。

        (2)Pd/α-Al2O3催化劑的制備 分別以H2PdCl4和α-Al2O3為前體和載體,利用等量浸漬法制備了3種不同粒徑的Pd催化劑,按照負載量由小到大順序分別記作Cat.1、Cat.2、Cat.3。制備方法如下:稱取適量已制備的α-Al2O3載體,攪拌條件下逐滴滴加一定濃度的H2PdCl4前體溶液至飽和吸水量,在室溫下老化12 h后于110℃下烘干12 h,然后在500℃下焙燒2 h,再將焙燒后的樣品在體積比1:1的H2/N2混合氣中150℃下還原2 h即得最終催化劑。

        1.2 催化劑的表征

        采用N2物理吸附(ASAP 2020,Micrometrics)測定載體的BET比表面積;采用D/max2550V型X射線衍射(XRD)測定載體的晶相結構;采用等離子體原子發(fā)射光譜法(ICP-AES)測定Pd的負載量;采用CO化學吸附(AutoChem 2920,Micrometrics)分別測定反應前的粒徑大小和進入穩(wěn)定期催化劑上Pd的分散度,化學吸附條件在此前的研究中給出[13];采用JEOL JEM 2100型透射電鏡(TEM)表征Pd的粒徑。

        1.3 動力學實驗

        反應裝置為Altamira公司的Microactivity全自動控制微型固定床反應器,不銹鋼反應管的長度和內徑分別為30 cm和0.8 cm。催化劑裝填量為0.03 g,與0.4 g石英砂催化劑混合均勻后加入反應管中,催化劑在總壓2.1 MPa、H2含量1.6%、C2H2含量0.9%、C2H4含量75%(其余氣體為N2)、溫度80℃條件下進行了空白實驗及內外擴散排除實驗。結果表明:催化劑反應管壁和石英砂對加氫反應的影響可以排除;當催化劑顆粒粒徑小于0.2 mm、/小于2.7×10-4 g·min-1·ml-1時,內外擴散可以排除。本研究所用催化劑顆粒粒徑小于0.015 mm,氣體流速為130 ml·min-1,滿足以上條件,因此所獲數(shù)據(jù)為本征動力學數(shù)據(jù)。采用正交實驗法設計動力學實驗,在每組催化劑進行動力學實驗前,催化劑均在上述條件下反應10 h以保證催化劑活性達到穩(wěn)定[8]。反應前后的氣體濃度采用在線四通道Inficon 3000微型氣相色譜測定,每個樣品點取樣5次,取平均值作為實驗數(shù)據(jù)。

        1.4 乙炔加氫微觀反應動力學模型

        乙炔加氫的微觀動力學模型根據(jù)順序加氫機理建立,這一機理主要包含八步基元反應步驟(ST1~ST8)。各基元步驟活化能的數(shù)據(jù)來源于Norskov等[10]基于Pd(111)面的DFT模擬結果,并根據(jù)Cat.2催化劑上的動力學實驗結果對吸脫附活化能進行優(yōu)化,使模擬結果與實驗值一致。指前因子數(shù)據(jù)來源中,ST1、ST5、ST8使用Sellers等[14]的數(shù)據(jù),ST2使用Salmeron等[15]和Sellers[16]的數(shù)據(jù)計算得出,其余的數(shù)據(jù)來源于Dumesic[17]的數(shù)據(jù)。各基元反應步驟和對應動力學參數(shù)如表1所示,其中a和分別代表活化能和指前因子,下角標+和-表示正向和逆向反應。各基元步驟正逆反應速率常數(shù)分別為kk,使用Arrhenius公式計算得到,如式(1)所示

        速率方程表達式r如表1所示,其中y代表氣體的分壓與標準壓力的比值;θ代表物種的表面覆蓋度,θ代表空白位點的表面覆蓋度,滿足

        穩(wěn)態(tài)實驗中反應器滿足平推流反應器假設,以單位質量催化劑為物料衡算單元,可以列出反應器模型如下。

        氣相物種的物料衡算式為

        催化劑表面物種的物料衡算式為

        式中,為氣體摩爾流率,mol·s-1;為催化劑質量,g;表示量綱1時間;site是催化劑活性位濃度,mol·(g cat)-1。

        表1 各基元反應步驟及參數(shù) Table 1 Elementary reaction steps and parameters

        2 實驗結果與討論

        2.1 催化劑的表征結果

        N2物理吸附結果表明所制備Al2O3的比表面積為14.9 m2·g-1,孔容為0.086 cm3·g-1,平均孔徑為20.0 nm。XRD表征結果顯示Al2O3衍射峰與PDF#81-1667卡片號吻合,說明本研究所用Al2O3為α-Al2O3。

        采用等離子體原子發(fā)射光譜法(ICP-AES)對3種催化劑的負載量進行了測試,通過透射電鏡(TEM)和CO化學吸附對催化劑Pd粒徑的大小進行了表征,結果如表2所示。

        表2 催化劑Pd粒徑表征結果 Table 2 Particle size analysis of catalysts

        ① TEM. ② CO chemisorption.

        圖1為不同粒徑的催化劑形貌及Pd粒徑分布。從TEM圖可以看出,通過調節(jié)負載量可以制備出不同粒徑的Pd催化劑,并且TEM所統(tǒng)計的粒徑大小與CO化學吸附表征結果接近。

        2.2 不同粒徑催化劑上的微觀動力學模擬

        2.2.1 活性位數(shù)的確定

        在進行微觀反應動力學模擬前需要確定參與反應的活性位總數(shù)site。本文研究的是Pd催化劑上乙炔加氫的穩(wěn)態(tài)動力學,催化劑在進行動力學實驗前經(jīng)過10 h反應使其活性達到穩(wěn)定。與初始時潔凈的Pd催化劑相比,此時Pd表面可能累積一定的含碳物種以及在次表面生成含碳、含氫物種[8,18],因此表面暴露的Pd原子數(shù)量將減少,同時也可以形成具有不同催化性質的活性位。由于微觀反應動力學基于平均場理論的假定,因此在本文中不考慮這種不同活性位變化的影響。本文通過CO化學吸附測定達到穩(wěn)定期后各催化劑表面Pd原子總數(shù)at,s,如表3所示。

        由于本文各基元步驟動力學參數(shù)來自基于Pd(111)面上的DFT計算結果,因此有必要檢驗at,s與Pd(111)面原子數(shù)的關系。根據(jù)Hardeveld等[19]的報道,催化劑表面金屬顆粒各晶面原子數(shù)比例與金屬顆粒的直徑均存在一定的函數(shù)關系,假定Pd顆粒為截角八面體結構,則可以根據(jù)測得的催化劑粒徑大小直接計算出一定質量催化劑上對應Pd(111)面的Pd原子個數(shù)at,(111),計算結果同樣列入表3。由表中數(shù)據(jù)可知,通過CO化學吸附測得的表面Pd原子數(shù)與Hardeveld模型預測的表面Pd(111)面原子數(shù)很接近。雖然這并不是Pd(111)作為穩(wěn)定期Pd催化劑活性中心的直接證據(jù),但也旁證了采用Pd(111)面上的DFT計算結果的合理性。本文采用CO化學吸附測量的at,s作為動力學模型中參與反應的活性位總數(shù)。

        2.2.2 動力學模擬結果

        根據(jù)Cat.2催化劑的動力學實驗結果,對原文獻的參數(shù)值進行了微調,結果列入表1。采用這些動力學參數(shù)分別模擬3種催化劑在對應反應條件下的結果,模擬結果如圖2和圖3所示。從圖中可以看出,微觀動力學模型不僅可以很好地模擬Cat.2催化劑在不同H2和C2H2分壓下的乙炔出口濃度值,也同樣可以模擬Cat.1和Cat.3催化劑上的動力學實驗結果。這說明了微觀反應動力學模型的廣泛適用性。

        2.2.3 微觀動力學模擬分析

        利用表1的動力學參數(shù)進行模擬,可以得到不同條件下各基元步驟的正逆反應速率及催化劑表面物種覆蓋度等信息。以Cat.1和Cat.3為例,不同H2和C2H2分壓下微觀動力學模擬的催化劑表面覆蓋度如表4所示。從表中可以看出,不同粒徑催化劑在相同反應條件下各表面物種覆蓋度差別很小,并且C2H2*和C2H3*隨著C2H2分壓的增加而增加,而隨H2分壓的增加而減小。

        表4 微觀動力學表面覆蓋度模擬結果 Table 4 Microkinetics analysis results of surface coverage

        綜合所有計算結果,在本文反應條件下3種催化劑表面最豐物種均為C2H4*?;襟EST4的反應速率最小,以Cat.3為例,在C2H20.9%、H22.1%、C2H475%、80℃反應條件下,其值為1.08 s-1;這一數(shù)值與ST1和ST3步驟的正反應速率接近,對應條件下其值均為2.05 s-1左右,而這些基元步驟的速率均遠小于其他步驟速率(約為20 s-1)。

        2.3 微觀與宏觀動力學的聯(lián)系

        由于ST1、ST3和ST4基元步驟速率接近,進一步通過微觀與宏觀動力學的關系以檢驗該反應的速率控制步驟。

        2.3.1 動力學實驗結果

        以考評條件(總壓2.1 MPa、H2含量1.6%、C2H2含量0.9%、C2H4含量75%、溫度80℃)為中心,分別改變C2H2和H2分壓進行動力學實驗,考察不同粒徑催化劑上乙炔消失速率分別與乙炔濃度和氫氣濃度的關系。以Cat.1和Cat.3為例,結果如圖4和圖5所示。從圖中可以看出,乙炔消失速率對乙炔的反應級數(shù)接近0.6,而氫氣對乙炔消失反應的級數(shù)均在0.7左右。

        2.3.2 反應級數(shù)的微觀解釋

        假設1:乙炔加氫生成乙烯反應過程中的決速步驟為ST4,因此總包反應速率可近似由ST4基元步驟的正反應速率代替

        由表4數(shù)據(jù)可知(C2H2*+H*)<< (C2H3*+C2H4*),因此式(2)可以寫成

        通過擬穩(wěn)態(tài)假設可以推導出各物種表面覆蓋度,即

        將式(14)~式(17)代入式(13),可得

        假設2:ST1為乙炔加氫生成乙烯反應過程中的決速步驟,推導過程類似,可以得到

        假設3:ST3為乙炔加氫生成乙烯反應過程中的決速步驟,則

        根據(jù)表4的計算結果,由假設1可以計算得到對應的C2H2和H2的反應級數(shù)分別為0.6和0.8左右,這與實驗結果基本一致;而對假設2和假設3,由式(22)和式(24)計算可知H2的反應級數(shù)分別為0級和0.3級左右,與實驗值不符。因此ST4更有可能是乙烯加氫生成乙烯過程的決速步驟。

        3 結 論

        (1)采用CO化學吸附測定活性達到穩(wěn)定后Pd催化劑的表面原子數(shù),結果與理論模型預測的Pd(111)表面原子數(shù)接近,這一數(shù)值作為微觀反應動力學模擬中催化活性位的總數(shù)是合理的。

        (2)利用優(yōu)化后的DFT計算結果作為微觀動力學參數(shù)(表1),在本文反應條件下,模型適用于Pd粒徑在1~5 nm的不同催化劑催化乙炔加氫反應,表現(xiàn)出很好的適用性。

        (3)在本文研究范圍內,微觀動力學模擬結果顯示3種粒徑催化劑上的乙炔加氫反應過程的微觀信息基本一致:表面最豐物種均為C2H4*,次豐物種為C2H3*,決速步驟為乙烯基加氫生成乙烯,但該步驟與乙炔吸附和亞乙烯基加氫速率接近。

        References

        [1] HUANG W, MCCORMICK J R, LOBO R F,. Selective hydrogenation of acetylene in the presence of ethylene on zeolite- supported bimetallic catalysts [J]. Journal of Catalysis, 2007, 246(1): 40-51.

        [2] 田亮, 蔣達, 錢鋒. 乙炔加氫反應系統(tǒng)操作優(yōu)化策略 [J]. 化工學報, 2015, 66(1):373-377.
        TIAN L, JIANG D, QIAN F. Reactor system switch strategy for acetylene hydrogenation process [J]. CIESC Journal, 2015, 66(1): 373-377.

        [3] 張齊, 戴偉, 穆瑋, 等. 原位紅外光譜法研究Pd-Ag/Al2O3和Pd/Al2O3乙炔加氫催化劑 [J]. 化工學報, 2011, 62(1): 71-77.
        ZHANG Q, DAI W, MU W,. Pd-Ag/Al2O3and Pd/Al2O3catalysts for selective hydrogenation of ethyne with-DRIFTS [J]. CIESC Journal, 2011, 62(1): 71-77.

        [4] KIM W J, KANG J H, AHN I Y,. Effect of potassium addition on the properties of a TiO2-modified Pd catalyst for the selective hydrogenation of acetylene [J]. Applied Catalysis A: General, 2004, 268(1): 77-82.

        [5] HUANG W, LI A, LOBO R F,. Effects of zeolite structures, exchanged cations, and bimetallic formulations on the selective hydrogenation of acetylene over zeolite-supported catalysts [J]. Catalysis Letters, 2009, 130(3/4): 380-385.

        [6] MOLERO H, BARTLETT B F, TYSOE W T. The hydrogenation of acetylene catalyzed by palladium: hydrogen pressure dependence [J]. Journal of Catalysis, 1999, 181(1): 49-56.

        [7] SARKANY A, BECK A, HORVáTH A,. Acetylene hydrogenation on sol-derived Pd/SiO2[J].Applied Catalysis A: General, 2003, 253(1): 283-292.

        [8] BORODZINSKI A, BOND G C. Selective hydrogenation of ethyne in ethane-rich streams on palladium catalysts(Ⅰ): Effect of changes to the catalyst during reaction [J]. Catalysis Reviews, 2006, 48(2): 91-144.

        [9] BORODZINSKI A, BOND G C. Selective hydrogenation of ethyne in ethane-rich streams on palladium catalysts(Ⅱ): Steady-state kinetics and effects of palladium particle size, carbon monoxide, and promoters [J]. Catalysis Reviews, 2008, 50(3): 379-469.

        [10] STUDT F, ABILD-PEDERSEN F, BLIGAARD T,. Identification of non-precious metal alloy catalysts for selective hydrogenation of acetylene [J]. Science, 2008, 320(5881): 1320-1322.

        [11] YANG B, BURCH R, HARDACRE C,. Influence of surface structures, subsurface carbon and hydrogen, and surface alloying on the activity and selectivity of acetylene hydrogenation on Pd surfaces: a density functional theory study [J]. Journal of Catalysis, 2013, 305(9): 264-276.

        [12] DUMESIC J A, TREVINO A A, MILLIGAN B A,. A kinetic modeling approach to the design of catalysts: formulation of a catalyst design advisory program [J]. Industrial & Engineering Chemistry Research, 1987, 26(7): 1399-1407.

        [13] 王新宇, 隋志軍, 朱貽安, 等. Pd/Al2O3催化乙炔選擇性加氫微觀反應動力學 [J]. 化學反應工程與工藝, 2015, 31(4): 322-329.
        WANG X Y, SUI Z J, ZHU Y A,. Microkinetics of acetylene selective hydrogenation over Pd/α-Al2O3catalyst [J]. Chemical Reaction Engineering and Technology, 2015, 31(4): 322-329.

        [14] SELLERS H, GISLASON J. Adsorption and desorption rate constants for small molecules on metal surfaces: an example of Trouton's rule [J]. Surface Science, 1999, 426(2): 147-153.

        [15] MITSUI T, ROSE M K, FOMIN E,. Hydrogen adsorption and diffusion on Pd(111) [J]. Surface Science, 2003, 540(1): 5-11.

        [16] SELLERS H. Entropies of desorption from temperature programmed desorption data: trends and applications to rate constant determinations [J]. The Journal of Physical Chemistry B, 2003, 107(37): 10206-10208.

        [17] DUMESIC J A. 多相催化微觀動力學 [M]. 沈儉一, 譯. 北京: 國防工業(yè)出版社, 1998: 160.
        DUMESIC J A. The Microkinetics of Heterogeneous Catalysis [M]. SHEN J Y, trans. Beijing: National Defense Industry Press, 1998: 160.

        [18] 盧根民, 陳誦英, 彭少逸. 動態(tài)分析解析解方法研究可逆吸附及不可逆吸附在多相催化反應中的作用 [J]. 催化學報, 1991, 12(4): 301-313.
        LU G M, CHEN S Y, PENG S Y. Study on the effect of reversible and irreversible adsorption on catalytic reaction by analytic solution model of dynamic method [J]. Journal of Catalysis, 1991, 12(4): 301-313.

        [19] HARDEVELD R V, HARTOG F. The statistics of surface atoms and surface sites on metal crystals [J]. Surface Science, 1969, 15(2): 189-230.

        [20] LYNGGAARD H, ANDREASEN A, STEGELMANN C,. Analysis of simple kinetic models in heterogeneous catalysis [J]. Progress in Surface Science, 2004, 77(3): 71-137.

        Microkinetics analysis of acetylene hydrogenation over Pd/Al2O3catalyst with different particle sizes

        PANG Hongqiang, SUI Zhijun, ZHU Yi’an, ZHOU Xinggui

        (State Key Laboratory of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China)

        Pd/α-Al2O3with different particle sizes were prepared by the incipient wetness method and characterized by BET, XRD, ICP-AES, CO chemisorption and TEM. The kinetic characteristics over these catalysts were explored following a fractional factorial designs of experiment and the experiment results were analyzed by microkinetic method. The size of three catalysts obtained was 1.6 nm, 3.4 nm and 5.5 nm, respectively, as verified by TEM and CO chemisorption results. After stable performance of the catalyst, the exposed surface Pd number of catalysts measured from CO chemisorption agreed well with that of Pd(111) surface calculated from Hardeveld’s model. The microkinetic analysis results showed that the microkinetic model can fit the kinetic experiment results quite well on all catalysts. C2H4* and C2H3* were the most abundant surface species independent of the reaction conditions. The rate determining step was the hydrogenation of vinyl group to form ethylene according to the relationship between microkinetic information and macrokinetic characteristics.

        multiphase reaction; acetylene hydrogenation; microkinetic model; particle size effects; catalyst

        Foundation item: supported by the National Basic Research Program of China (2012CB720500) and the National Natural Science Foundation of China (21376076).

        date: 2016-03-10.

        SUI Zhijun, zhjsui@ecust.edu.cn

        TQ 032.4

        A

        0438—1157(2016)09—3692—07

        10.11949/j.issn.0438-1157.20160283

        國家重點基礎研究發(fā)展計劃項目(2012CB720500);國家自然科學基金項目(21376076)。

        2016-03-10收到初稿,2016-06-16收到修改稿。

        聯(lián)系人:隋志軍。第一作者:龐洪強(1989—),男,碩士研究生。

        猜你喜歡
        基元乙炔微觀
        關注基元反應的考查
        烷美無炔
        人體細胞內存在全新DNA結構
        科學導報(2018年30期)2018-05-14 12:06:01
        一種新的結合面微觀接觸模型
        超導電乙炔炭黑DENKA BLACK(電池用途)
        電源技術(2016年9期)2016-02-27 09:05:25
        微觀的山水
        詩選刊(2015年6期)2015-10-26 09:47:10
        Numerical Modeling and Analysis of Gas Entrainment for the Ventilated Cavity in Vertical Pipe*
        微觀中國
        浙江人大(2014年8期)2014-03-20 16:21:15
        微觀中國
        浙江人大(2014年6期)2014-03-20 16:20:46
        在線激光乙炔分析儀在醋酸乙烯合成中的應用
        自動化博覽(2014年8期)2014-02-28 22:33:06
        一区五码在线| 亚洲av综合永久无码精品天堂| 久久精品国产亚洲av大全| 国产一区二区三区啪| 人妻尤物娇呻雪白丰挺| 亚洲精品有码日本久久久| 一二三四日本中文在线| 亚洲 欧美 综合 另类 中字| 中文无码制服丝袜人妻AV| 国产精品国产三级国产剧情| 国产欧美日韩精品丝袜高跟鞋| 亚洲中文字幕无码久久2020| A亚洲VA欧美VA国产综合| 青青草手机视频免费在线播放| 粗大的内捧猛烈进出少妇| 午夜无码片在线观看影院| 精品高清国产乱子伦| 东风日产车是不是国产的 | 日产精品一区二区在线| 久久久免费看少妇高潮| 欧洲熟妇色xxxxx欧美老妇伦| 久久99精品久久久久久齐齐百度| 国产精品日本中文在线| 国产精品亚洲精品日韩已方 | 国产情侣真实露脸在线| 亚洲av久久无码精品九九| 2022精品久久久久久中文字幕| 按摩师玩弄少妇到高潮av| 国产精品人妻一码二码| 一本大道东京热无码中字 | 国产精品人伦一区二区三| 亚洲国产精品毛片av不卡在线| 国产午夜亚洲精品不卡福利| 亚洲综合av一区在线| 无遮掩无码h成人av动漫| 亚洲欧洲中文日韩久久av乱码| 精品黄色av一区二区三区| 亚洲一区二区三区中文字幕网| 影音先锋女人av鲁色资源网久久| 欧美性一区| 亚洲本色精品一区二区久久|