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

    Fe,Co,Ni改性SAPO-34分子篩上MTO反應(yīng)的理論研究

    2018-10-15 09:49:42李對春劉紅艷王寶俊李瑞豐
    潔凈煤技術(shù) 2018年5期
    關(guān)鍵詞:丙烯烯烴催化活性

    李對春,邢 斌,劉紅艷,劉 平,王寶俊,李瑞豐

    (1.太原理工大學(xué) 化學(xué)化工學(xué)院,山西 太原 030024;2.大同大學(xué) 化學(xué)與環(huán)境工程學(xué)院,山西 大同 037009;3.中國科學(xué)院 山西煤炭化學(xué)研究所,山西 太原 030001;4.太原理工大學(xué) 煤科學(xué)與技術(shù)教育部和山西省重點實驗室,山西 太原 030024;5.太原理工大學(xué) 能源化工催化研究中心,山西 太原 030024)

    0 引 言

    低碳烯烴中乙烯、丙烯是現(xiàn)代化學(xué)工業(yè)的基本有機(jī)原料,隨著化工產(chǎn)業(yè)的飛速發(fā)展,其需求量將越來越大。制取乙烯、丙烯的傳統(tǒng)路線是通過石腦油裂解生產(chǎn),但由于石油是不可再生資源,儲量有限,價格起伏較大,尋求替代傳統(tǒng)的非石油路線生產(chǎn)低碳烯烴勢在必行[1-4]。甲醇制烯烴(methanol-to-olefins,MTO)是替代傳統(tǒng)路線的新工藝,是環(huán)境催化領(lǐng)域的研究熱點[2]。目前,甲醇的制備路線已成熟完善,不僅可通過生物質(zhì)發(fā)酵制得,還可以通過煤、天然氣經(jīng)由合成氣(CO+H2)制得。而MTO是該替代路線的關(guān)鍵環(huán)節(jié)。

    目前,共提出20多種甲醇制烯烴反應(yīng)機(jī)理,但至今沒有得出與反應(yīng)完全符合的機(jī)理[2],其中,最具代表性的為直接反應(yīng)機(jī)理和間接反應(yīng)機(jī)理。根據(jù)反應(yīng)過程中產(chǎn)生的特殊中間體,直接反應(yīng)機(jī)理又被分為卡賓機(jī)理[1,5]、氧鎓葉立德機(jī)理[6-7]、碳正離子機(jī)理[8]、自由基機(jī)理[9]等。Dahl等[10-12]提出了間接反應(yīng)機(jī)理(烴池機(jī)理),即限域在分子篩孔道內(nèi)的烴池物種與分子篩作為共催化劑,通過與甲醇發(fā)生連續(xù)的甲基化,使碳鏈增長,經(jīng)低碳烯烴消除反應(yīng)生成產(chǎn)物烯烴(乙烯和丙烯),同時,烴池物種和分子篩的B酸性位點再生,完成一個催化循環(huán),烴池機(jī)理得到了研究者們的普遍認(rèn)可。以烯烴為基礎(chǔ)的循環(huán)起源于30年前由Dessau[13-14]提出的烯烴甲基化和裂化路線。王傳明等[15]通過密度泛函理論證明了SAPO-34分子篩上發(fā)生的MTO反應(yīng)除了聚甲基苯外,烯烴本身也是主要的烴池物種,即在MTO催化反應(yīng)中同時存在芳烴路線和烯烴路線。

    在MTO反應(yīng)中,SAPO-34分子篩表現(xiàn)出優(yōu)異的催化性能。SAPO-34分子篩具有8元環(huán)構(gòu)成的橢球形籠和三維孔道結(jié)構(gòu),這種分子篩的小孔結(jié)構(gòu)、中等酸性強(qiáng)度及良好的水熱穩(wěn)定性使其具有良好的MTO催化性能[16-17]。不同的分子篩具有不同的骨架拓?fù)浣Y(jié)構(gòu)及酸性質(zhì),包括酸性位點的分布、酸性位點的密度及強(qiáng)度等,這些因素決定了其催化活性和產(chǎn)物的選擇性。其中,分子篩催化劑酸性位點的強(qiáng)度是影響MTO反應(yīng)催化性能的重要因素[18-20],但通過試驗手段無法明確單個因素對MTO反應(yīng)催化性能的影響。因此,本文采用色散力矯正的密度泛函理論方法,構(gòu)建SAPO-34分子篩模型,以2,3-二甲基-2-丁烯(iso-C6)為烴池物種[21],選擇不同金屬雜原子Fe、Co和Ni同晶取代Al原子對SAPO-34分子篩催化劑進(jìn)行改性,考察了不同金屬原子Fe、Co和Ni改性的SAPO-34 分子篩中B酸位點強(qiáng)度的變化,并分析改性SAPO-34分子篩對MTO反應(yīng)的催化活性和產(chǎn)物選擇性的影響。

    1 計算方法與模型

    1.1 計算方法

    采用VASP軟件包的密度泛函理論方法(DFT)進(jìn)行結(jié)構(gòu)優(yōu)化和總能量計算[22-24],電子的交換關(guān)聯(lián)勢采用廣義梯度近似(GGA)下的PBE泛函來描述[25],用PAW方法描述離子實,價電子用平面波基組擴(kuò)展[26-27],同時考慮了庫倫相互作用的色散效應(yīng)[28]。設(shè)定截斷能為400 eV,通過CI-NEB方法尋找過渡態(tài)[29],并通過頻率分析確認(rèn)過渡態(tài),在計算過程中,還考慮了零點能(ZPE)校正[30-31]。對于SAPO-34分子篩及Fe、Co和Ni改性的SAPO-34分子篩,布里淵帶使用1×1×1的Monkhost-Pack篩狀的k點空間取樣[32]。計算過程中基于自旋極化方法考慮了Fe、Co和Ni原子的磁性,結(jié)構(gòu)優(yōu)化時晶胞參數(shù)保持不變,直至每個原子上的力小于0.02 eV/Angst。

    1.2 計算模型

    SAPO-34分子篩具有菱沸石(CHA)骨架結(jié)構(gòu)單元[33],屬于小孔分子篩[34]。CHA晶胞單元中含有36個T原子(1個Si原子,18個Al原子,17個P原子)和72個O原子。由于Si、P原子價態(tài)不同,需要用額外的一個質(zhì)子(H+)補(bǔ)償晶胞中產(chǎn)生的一個負(fù)電荷來保持整體結(jié)構(gòu)的電中性,該位點(Si-O(H)-Al)即為Br?nsted酸性中心[35]。SAPO-34的周期性結(jié)構(gòu)如圖1所示,優(yōu)化的晶胞參數(shù)為a=b=13.85×10-10m,c=15.08×10-10m,與國際分子篩協(xié)會數(shù)據(jù)庫中的晶胞參數(shù)(a=b=13.68×10-10m,c=14.77×10-10m)基本相符[36]。

    金屬雜原子Fe、Co和Ni通過同晶取代Al原子的方式進(jìn)入分子篩骨架,得到MeAPO-34(Me=Fe、Co、Ni)催化劑。由于Fe原子與Al原子的價態(tài)相同,所以Fe(+3)原子替換Al(+3)原子時不會產(chǎn)生電荷差,而Co(+2)原子和Ni(+2)原子與Al(+3)原子的價態(tài)不同,所以替換Al原子時,同樣用質(zhì)子(H+)補(bǔ)償晶胞中產(chǎn)生一個負(fù)電荷。FeAPO-34,CoAPO-34和NiAPO-34分子篩優(yōu)化后的穩(wěn)定結(jié)構(gòu)如圖2所示,計算模型中所有的原子保持弛豫。

    2 結(jié)果與討論

    2.1 金屬雜原子引入SAPO-34分子篩后對酸性強(qiáng)度的影響

    分子篩催化劑的B酸位點的強(qiáng)度是影響MTO反應(yīng)催化性能的重要因素[19-20,37]。B酸定義為任何化學(xué)物種(分子或離子)能夠失去或貢獻(xiàn)出氫離子(質(zhì)子)。去質(zhì)子化能是評價酸性強(qiáng)度的標(biāo)準(zhǔn)之一,根據(jù)B酸的定義,本文采用去質(zhì)子化能EDPE(從Si-O(H)-Al位點上失去質(zhì)子(H+)所需要的能)描述B酸性位點的酸性強(qiáng)度[38-40]。計算公式[39-40]為

    圖2 3種金屬原子Fe、Co、Ni改性的SAPO-34結(jié)構(gòu)示意Fig.2 Configuration of FeAPO-34,CoAPO-34 and NiAPO-34

    EDPE=EMeAPO-34-+EH+-EMeAPO-34

    其中,EMeAPO-34為MeAPO-34分子篩的結(jié)構(gòu)能量;EMeAPO-34-為MeAPO-34結(jié)構(gòu)中失去一個質(zhì)子H+后的結(jié)構(gòu)能量;EH+為一個質(zhì)子H+的能量。EDPE越小,對應(yīng)的B酸位點越容易貢獻(xiàn)出H+質(zhì)子,B酸位點的酸性越強(qiáng)。EDPE的計算結(jié)果見表1,其與Jones.等[41]的計算結(jié)果一致。B酸位點周圍有4個O位點,為了比較不同MeAPO-34結(jié)構(gòu)中B酸的酸性強(qiáng)度,分別計算4個O位點的EDPE值和EDPE平均值(Eav-DPE)(表1)。

    Eav-DPE=(EDPE-1+EDPE-2+EDPE-3+EDPE-4)/4

    從表1可知,當(dāng)3種金屬雜原子Fe、Co和Ni引入SAPO-34結(jié)構(gòu),NiAPO-34結(jié)構(gòu)的Eav-DPE值比未改性SAPO-34結(jié)構(gòu)的Eav-DPE值大,說明NiAPO-34結(jié)構(gòu)中的B酸位點的強(qiáng)度較弱。但FeAPO-34和CoAPO-34結(jié)構(gòu)的Eav-DPE值均比未改性SAPO-34結(jié)構(gòu)的Eav-DPE值小,說明FeAPO-34和CoAPO-34結(jié)構(gòu)中的B酸位點的強(qiáng)度較強(qiáng)。因此,改性前后SAPO-34結(jié)構(gòu)的B酸位點的酸性強(qiáng)度順序為FeAPO-34>CoAPO-34>SAPO-34>NiAPO-34,表明金屬雜原子引入SAPO-34結(jié)構(gòu)中能改變分子篩的B酸位點的酸性強(qiáng)度。

    表1SAPO-34結(jié)構(gòu)及MeAPO-34結(jié)構(gòu)中的EDPE和Eav-DPE
    Table1CalculatedEDPEandEav-DPEofSAPO-34andMeAPO-34structure

    結(jié)構(gòu)位點nEDPE/(kJ·mol-1)Eav-DPE/(kJ·mol-1)11 173.5SAPO-34-Si-OHn-Al21 166.71 168.131 158.541 173.711 037.4FeAPO-34-Si-OHn-Al21 009.01 017.431 006.341 016.911 092.4CoAPO-34-Si-OHn-Al21 082.41 073.431 039.141 079.811 191.9NiAPO-34-Si-OHn-Al21 159.81 174.031 174.841 169.6

    2.2 金屬雜原子引入SAPO-34分子篩后對MTO反應(yīng)性能的影響

    在MeAPO-34穩(wěn)定結(jié)構(gòu)的基礎(chǔ)上,基于有機(jī)活性中心為iso-C6分子烯烴循環(huán)的烴池機(jī)理,即iso-C6分子與分子篩作為共催化劑,與甲醇發(fā)生連續(xù)的甲基化,使碳鏈增長,然后發(fā)生低碳烯烴消除反應(yīng),生成產(chǎn)物乙烯和丙烯,同時,烴池物種iso-C6分子和分子篩的B酸位點再生,形成一個完整的催化循環(huán)。反應(yīng)路線如圖3所示。

    反應(yīng)步驟可分為:

    1)甲基化反應(yīng)

    (M1)

    (CH3)3C+—C(CH3)(CH2CH3)+H2O+Z-

    (M2)

    (M3)

    圖3 MTO反應(yīng)的烯烴循環(huán)路徑Fig.3 Reaction pathway of olefin-based for MTO conversion

    2)過渡反應(yīng)

    (CH3)3C—C(CH3)(CH2)+H2O+HZ

    (D1)

    (CH3)3C—C(CH3)(CHCH3)+H2O+HZ

    (D2)

    H2O+Z-

    (C1)

    C+(CH3)2(C3H7)+H2O+Z-

    (C2)

    3)烯烴脫除反應(yīng)

    (E1)

    (E2)

    各個反應(yīng)步驟的活化能(Ea)和反應(yīng)熱(ΔH)見表2。本文選擇甲基化反應(yīng)步驟(M1、M2、M3)和低碳烯烴(乙烯和丙烯)消除步驟(E1、E2)來評價分子篩的酸性強(qiáng)度對MTO反應(yīng)催化性能的影響。

    表2MeAPO-34結(jié)構(gòu)中不同反應(yīng)步驟的活化能
    (Ea)和反應(yīng)熱(ΔH)
    Table2CalculatedEaandΔHofdifferentreactionstepsoverMeAPO-34structurekJ/mol

    MTO反應(yīng)的第1步是甲基化反應(yīng)(M1),其反應(yīng)物、產(chǎn)物和過渡態(tài)(S1、S2、TS1)結(jié)構(gòu)如圖4所示。在FeAPO-34,CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生M1反應(yīng)的反應(yīng)熱分別為26.1、-13.5和-28.9 kJ/mol,需要克服的能壘分別為104.2、99.4和112.9 kJ/mol。在NiAPO-34結(jié)構(gòu)中,M1反應(yīng)放熱最多,需要克服的能壘最高,說明Ni的存在不利于CH3OH的活化;但在FeAPO-34和CoAPO-34結(jié)構(gòu)中,M1反應(yīng)需要克服的能壘較接近,說明FeAPO-34和CoAPO-34對CH3OH活化的催化活性相近。

    第2步甲基化反應(yīng)(M2)的反應(yīng)物、產(chǎn)物和過渡態(tài)(S4、S5、TS3)結(jié)構(gòu)與M1反應(yīng)的S1、S2、TS1結(jié)構(gòu)相似,其鍵長見表3。在FeAPO-34、CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生M2反應(yīng)的反應(yīng)熱分別為-25.1、34.7和-39.6 kJ/mol,需要克服的能壘分別為98.4、100.3和148.6 kJ/mol。NiAPO-34結(jié)構(gòu)仍不利于CH3OH的活化,F(xiàn)eAPO-34和CoAPO-34結(jié)構(gòu)對CH3OH活化的催化活性相近。

    乙烯消除反應(yīng)步驟(E1)的反應(yīng)物、產(chǎn)物和過渡態(tài)(S6,S7,TS5)結(jié)構(gòu)如圖5所示。在FeAPO-34、CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生乙烯消除反應(yīng)的反應(yīng)熱分別為-12.5、5.8和-54.0 kJ/mol,需要克服的能壘分別為86.8、130.3和83.0 kJ/mol。結(jié)果表明:FeAPO-34和NiAPO-34表現(xiàn)出較好的催化活性,而CoAPO-34的催化活性較差,與Kang等[42]通過快速水熱晶化法合成的Fe、Co、Ni三種金屬改性的SAPO-34分子篩上對乙烯選擇性的順序一致。

    FeAPO-34、CoAPO-34和NiAPO-34的生成產(chǎn)物丙烯的反應(yīng)路徑中,甲基化反應(yīng)(M3)的反應(yīng)熱分別為32.8、-6.8和-56.0 kJ/mol,需要克服的能壘分別為104.2、93.6和138.0 kJ/mol,與M1、M2反應(yīng)的計算結(jié)果相似,即NiAPO-34結(jié)構(gòu)對甲基化反應(yīng)的催化活性最差,F(xiàn)eAPO-34和CoAPO-34結(jié)構(gòu)對甲基化反應(yīng)的催化活性相近。M3甲基化反應(yīng)步驟中的反應(yīng)物、產(chǎn)物和過渡態(tài)(S9,S10,TS7)結(jié)構(gòu)的鍵長見表3。

    丙烯消除反應(yīng)(E2)的反應(yīng)物、產(chǎn)物和過渡態(tài)(S11,S12,TS9)結(jié)構(gòu)與E1反應(yīng)中的S6,S7,TS5結(jié)構(gòu)相似,其鍵長見表4。在FeAPO-34,CoAPO-34和NiAPO-34結(jié)構(gòu)上發(fā)生丙烯消除反應(yīng)的反應(yīng)熱分別為-9.6、-34.7和46.3 kJ/mol,需要克服的能壘分別為38.6、41.5和131.2 kJ/mol,可知,對于丙烯消除反應(yīng)E2,F(xiàn)eAPO-34和CoAPO-34結(jié)構(gòu)催化活性明顯比NiAPO-34結(jié)構(gòu)的催化活性好。

    圖4 不同MeAPO-34結(jié)構(gòu)中M1的反應(yīng)物(S1)、產(chǎn)物(S2)和過渡態(tài)(TS1)結(jié)構(gòu)示意Fig.4 Configurations of reactant (S1),product (S2) and transition state (TS1) of M1 reaction step over different MeAPO-34 structures

    反應(yīng)步驟MeAPO-34結(jié)構(gòu)結(jié)構(gòu)鍵長/10-10 mC1—C2C2—O1H1—O1H1—O2FeAPO-343.361.451.431.06CoAPO-34反應(yīng)物(S4)3.401.461.271.17NiAPO-343.621.451.601.03FeAPO-342.242.041.021.40M2CoAPO-34過渡態(tài)(TS3)2.251.991.011.76NiAPO-342.251.991.011.76FeAPO-341.463.101.041.62CoAPO-34產(chǎn)物(S5)1.583.970.981.97NiAPO-341.583.970.981.97FeAPO-343.531.461.421.08CoAPO-34反應(yīng)物(S9)3.781.451.341.12NiAPO-343.261.462.241.04FeAPO-342.322.041.021.61M3CoAPO-34過渡態(tài)(TS7)2.351.991.011.71NiAPO-342.192.191.011.45FeAPO-341.604.010.983.00CoAPO-34產(chǎn)物(S10)1.623.200.991.88NiAPO-341.594.371.002.09

    注:C1為有機(jī)分子的C原子;C2為CH3OH分子中的C原子;O1為CH3OH分子中的O原子;H1為分子篩骨架上B酸位點的H原子;O2為分子篩骨架上B酸位點的O原子。

    表4不同MeAPO-34結(jié)構(gòu)上丙烯消除反應(yīng)E2的反應(yīng)物,產(chǎn)物和過渡態(tài)結(jié)構(gòu)中的鍵長
    Table4Correspondingbondlengthofreactant,productandtransitionstateofE2reactionstepsoverdifferentMeAPO-34structures

    MeAPO-34結(jié)構(gòu)結(jié)構(gòu)鍵長/10-10 mC1—C2'C3—H1'H1'—O1'H2—O1'H2—O2FeAPO-341.581.102.770.991.87CoAPO-34反應(yīng)物(S11)1.571.102.510.982.04NiAPO-341.561.102.520.992.25FeAPO-342.891.291.381.051.52CoAPO-34過渡態(tài)(TS9)3.211.341.341.011.73NiAPO-342.901.261.431.012.06FeAPO-343.233.020.981.551.03CoAPO-34產(chǎn)物(S12)3.251.951.011.261.18NiAPO-343.591.951.011.261.18

    注:C2′為丙烯分子中的C原子;C3為丙烯分子中的C原子;H1′為丙烯分子中的H原子;O1′為H2O分子中的O原子;H2為H2O分子中的H原子。

    圖5 不同MeAPO-34結(jié)構(gòu)中乙烯消除反應(yīng)E1的反應(yīng)物(S6)、產(chǎn)物(S7)和過渡態(tài)(TS5)結(jié)構(gòu)示意Fig.5 Configurations of reactant (S6),product (S7) and transition state (TS5) of E1 reaction step over different MeAPO-34 structures

    由表4可知,在3種不同的MeAPO-34結(jié)構(gòu)中,F(xiàn)eAPO-34和CoAPO-34上丙烯消除反應(yīng)E2需要克服的反應(yīng)能壘比乙烯消除反應(yīng)E1要低,而NiAPO-34上丙烯消除反應(yīng)E2需要克服的反應(yīng)能壘比乙烯消除反應(yīng)E1高,說明在本文的烯烴循環(huán)反應(yīng)路徑中,F(xiàn)eAPO-34和CoAPO-34分子篩對丙烯的選擇性好,NiAPO-34分子篩對乙烯的選擇性好。

    3 結(jié) 論

    1)通過比較不同MeAPO-34結(jié)構(gòu)中Si-O(H)-Al酸性位點的去質(zhì)子化能可知,MeAPO-34分子篩的B酸位的酸性強(qiáng)度順序為FeAPO-34>CoAPO-34>SAPO-34>NiAPO-34,不同金屬雜原子的引入對SAPO-34分子篩B酸位點的酸性強(qiáng)度起調(diào)節(jié)作用。

    2)以iso-C6分子為烴池物種,研究烯烴循環(huán)反應(yīng)路徑對MTO反應(yīng)中的甲基化反應(yīng)步驟(M1、M2、M3)和烯烴消除反應(yīng)(E1、E2)的催化活性,結(jié)果表明,F(xiàn)e、Co和Ni改性的SAPO-34分子篩、 M1、M2、M3反應(yīng)的催化活性順序為FeAPO-34≈CoAPO-34>NiAPO-34;對乙烯生成反應(yīng)(E1)的催化活性順序為FeAPO-34≈NiAPO-34>CoAPO-34;對丙烯生成反應(yīng)(E2)的催化活性順序為FeAPO-34≈CoAPO-34>NiAPO-34。

    3)NiAPO-34結(jié)構(gòu)中B酸位的酸性強(qiáng)度比FeAPO-34和CoAPO-34結(jié)構(gòu)中的弱,表明強(qiáng)酸性位點對甲基化反應(yīng)(M1、M2、M3)和丙烯消除反應(yīng)E2的催化活性好,而弱酸性位點對乙烯消除反應(yīng)E1的催化活性好。

    4)對于烯烴循環(huán)反應(yīng)路徑而言,F(xiàn)eAPO-34和CoAPO-34分子篩對丙烯的選擇性好,NiAPO-34分子篩對乙烯的選擇性好,即強(qiáng)酸位點有利于丙烯產(chǎn)物的生成,而弱酸位點有利于乙烯產(chǎn)物的生成。計算結(jié)果為MTO催化反應(yīng)中的分子篩催化劑改性提供理論依據(jù),從而制備出產(chǎn)物選擇性和催化活性更好的分子篩催化劑。

    猜你喜歡
    丙烯烯烴催化活性
    錢愛康
    《共生》主題系列作品
    大眾文藝(2022年16期)2022-09-07 03:08:04
    煤制烯烴副產(chǎn)混合碳四的綜合利用
    云南化工(2021年8期)2021-12-21 06:37:38
    烯烴不稀罕
    苯丙烯菌酮
    液化氣中的丙烯含有的雜質(zhì)對丙烯聚合反應(yīng)的影響
    MTO烯烴分離回收技術(shù)與烯烴轉(zhuǎn)化技術(shù)
    稀土La摻雜的Ti/nanoTiO2膜電極的制備及電催化活性
    環(huán)化聚丙烯腈/TiO2納米復(fù)合材料的制備及可見光催化活性
    過硼酸鈉對高碳烯烴的環(huán)氧化
    欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品一区av在线观看| 国产在线观看jvid| 久久久精品欧美日韩精品| 国产精品久久视频播放| 国产成人欧美在线观看| 中文字幕高清在线视频| 一边摸一边做爽爽视频免费| 国产精品久久久av美女十八| 日韩欧美免费精品| 久久久久国产一级毛片高清牌| 高潮久久久久久久久久久不卡| 亚洲最大成人中文| 天天躁狠狠躁夜夜躁狠狠躁| 美女高潮到喷水免费观看| 一区二区三区精品91| 亚洲欧美精品综合久久99| 亚洲中文日韩欧美视频| 身体一侧抽搐| 欧美激情 高清一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 夜夜夜夜夜久久久久| 欧美最黄视频在线播放免费| 视频区欧美日本亚洲| 国产一区在线观看成人免费| 国产亚洲精品久久久久5区| 嫁个100分男人电影在线观看| 老司机午夜福利在线观看视频| 精品卡一卡二卡四卡免费| 欧美日韩黄片免| 他把我摸到了高潮在线观看| 一区二区三区精品91| 中文字幕最新亚洲高清| 亚洲欧美精品综合一区二区三区| 午夜福利一区二区在线看| 久久中文字幕一级| 狂野欧美激情性xxxx| 午夜福利影视在线免费观看| 巨乳人妻的诱惑在线观看| 亚洲九九香蕉| 中出人妻视频一区二区| 中文字幕精品免费在线观看视频| 国产真人三级小视频在线观看| 香蕉丝袜av| 操美女的视频在线观看| 女性被躁到高潮视频| 精品久久久久久,| 国产精品精品国产色婷婷| 国产成人av教育| 欧美一级a爱片免费观看看 | 久久 成人 亚洲| 又紧又爽又黄一区二区| 老汉色∧v一级毛片| 一区二区三区国产精品乱码| 老熟妇仑乱视频hdxx| 深夜精品福利| a在线观看视频网站| 十分钟在线观看高清视频www| 亚洲欧美日韩高清在线视频| 亚洲三区欧美一区| 亚洲成人国产一区在线观看| 男女下面插进去视频免费观看| 国产成人一区二区三区免费视频网站| 多毛熟女@视频| 成人精品一区二区免费| 久久婷婷成人综合色麻豆| 精品人妻在线不人妻| 欧美日韩乱码在线| 欧美色欧美亚洲另类二区 | 国产在线观看jvid| 国产三级黄色录像| 国产精华一区二区三区| 日本三级黄在线观看| 超碰成人久久| 12—13女人毛片做爰片一| 一级a爱片免费观看的视频| 天堂√8在线中文| 久久亚洲真实| 老司机午夜福利在线观看视频| 精品卡一卡二卡四卡免费| 色哟哟哟哟哟哟| avwww免费| 国产精品精品国产色婷婷| 欧美日韩福利视频一区二区| 我的亚洲天堂| 精品不卡国产一区二区三区| 久久国产精品男人的天堂亚洲| 亚洲精品粉嫩美女一区| 亚洲人成电影观看| 久久影院123| 午夜两性在线视频| 国产精品免费一区二区三区在线| 嫁个100分男人电影在线观看| 亚洲情色 制服丝袜| 亚洲国产高清在线一区二区三 | 日韩国内少妇激情av| 精品久久久精品久久久| 在线观看www视频免费| 久久狼人影院| 两个人视频免费观看高清| 一级,二级,三级黄色视频| aaaaa片日本免费| 免费在线观看日本一区| 女人被躁到高潮嗷嗷叫费观| 亚洲激情在线av| 99国产精品一区二区蜜桃av| 日韩一卡2卡3卡4卡2021年| 九色国产91popny在线| 欧美日韩精品网址| 波多野结衣av一区二区av| 黄色 视频免费看| 久久精品aⅴ一区二区三区四区| 美女高潮到喷水免费观看| 午夜福利视频1000在线观看 | 亚洲最大成人中文| 欧美激情久久久久久爽电影 | 免费高清视频大片| 日本撒尿小便嘘嘘汇集6| 久久午夜综合久久蜜桃| 麻豆成人av在线观看| 黄色女人牲交| 香蕉丝袜av| 高清毛片免费观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 又黄又爽又免费观看的视频| 亚洲国产精品999在线| 啦啦啦 在线观看视频| 午夜免费激情av| 国产欧美日韩综合在线一区二区| 国产1区2区3区精品| 亚洲成人久久性| 欧美av亚洲av综合av国产av| 欧美日韩乱码在线| 亚洲成国产人片在线观看| 女人爽到高潮嗷嗷叫在线视频| 悠悠久久av| 岛国视频午夜一区免费看| 桃红色精品国产亚洲av| 制服人妻中文乱码| 亚洲电影在线观看av| 两个人免费观看高清视频| 最近最新中文字幕大全免费视频| 亚洲精品国产一区二区精华液| 国产精品亚洲美女久久久| 无遮挡黄片免费观看| 国产精品自产拍在线观看55亚洲| 真人一进一出gif抽搐免费| 久久热在线av| 美女国产高潮福利片在线看| 日韩精品免费视频一区二区三区| 又大又爽又粗| 在线观看日韩欧美| 91九色精品人成在线观看| 韩国精品一区二区三区| 91成年电影在线观看| 99久久国产精品久久久| 成人三级做爰电影| 欧美成人午夜精品| 宅男免费午夜| 久9热在线精品视频| 人妻丰满熟妇av一区二区三区| 夜夜看夜夜爽夜夜摸| 精品人妻在线不人妻| 国产成人精品无人区| 亚洲精华国产精华精| 成人三级做爰电影| 一级毛片高清免费大全| 日韩欧美国产一区二区入口| 亚洲九九香蕉| 色综合欧美亚洲国产小说| 一级毛片高清免费大全| 999久久久国产精品视频| av有码第一页| 91国产中文字幕| 性欧美人与动物交配| 丰满人妻熟妇乱又伦精品不卡| 国产精品国产高清国产av| 亚洲国产看品久久| 又大又爽又粗| 真人做人爱边吃奶动态| 久久人人精品亚洲av| 国产黄a三级三级三级人| 久久这里只有精品19| 亚洲av成人一区二区三| 丁香欧美五月| 国产成人影院久久av| 久久影院123| 欧美不卡视频在线免费观看 | 一卡2卡三卡四卡精品乱码亚洲| 好看av亚洲va欧美ⅴa在| 亚洲欧美精品综合一区二区三区| 亚洲精品美女久久av网站| 视频在线观看一区二区三区| 两个人看的免费小视频| 老司机靠b影院| 国产成年人精品一区二区| 一个人观看的视频www高清免费观看 | 午夜激情av网站| 一区福利在线观看| 69av精品久久久久久| 国产激情久久老熟女| 国产精品免费视频内射| 欧美黑人欧美精品刺激| 制服诱惑二区| 国产成人啪精品午夜网站| 亚洲精品av麻豆狂野| 亚洲av五月六月丁香网| 午夜福利,免费看| 女警被强在线播放| 国产精品一区二区在线不卡| 一进一出抽搐gif免费好疼| 亚洲中文av在线| 国产区一区二久久| 欧美乱码精品一区二区三区| 精品福利观看| 一二三四社区在线视频社区8| 视频在线观看一区二区三区| 99精品在免费线老司机午夜| 欧美+亚洲+日韩+国产| 老汉色av国产亚洲站长工具| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品合色在线| 岛国视频午夜一区免费看| 久久这里只有精品19| 九色亚洲精品在线播放| 久久久水蜜桃国产精品网| 久久九九热精品免费| 免费观看精品视频网站| 搡老岳熟女国产| 亚洲精品中文字幕一二三四区| 久久精品成人免费网站| 高潮久久久久久久久久久不卡| 成人三级黄色视频| 国产成人影院久久av| 久久 成人 亚洲| 亚洲精品一区av在线观看| 免费看十八禁软件| 免费在线观看日本一区| 涩涩av久久男人的天堂| 免费人成视频x8x8入口观看| 黑人欧美特级aaaaaa片| 涩涩av久久男人的天堂| 亚洲av熟女| 一级毛片女人18水好多| 亚洲视频免费观看视频| 18禁美女被吸乳视频| 国产人伦9x9x在线观看| 伊人久久大香线蕉亚洲五| av视频在线观看入口| 欧美色欧美亚洲另类二区 | 久久久国产欧美日韩av| 午夜福利视频1000在线观看 | 国产真人三级小视频在线观看| 国产精品 国内视频| 精品久久蜜臀av无| 麻豆久久精品国产亚洲av| 黄色视频不卡| 极品人妻少妇av视频| 老司机深夜福利视频在线观看| 亚洲国产精品999在线| 男人舔女人下体高潮全视频| 国产精品美女特级片免费视频播放器 | 亚洲狠狠婷婷综合久久图片| 国产片内射在线| 亚洲av成人不卡在线观看播放网| 亚洲熟妇熟女久久| 亚洲一区高清亚洲精品| 成年版毛片免费区| 亚洲精品粉嫩美女一区| 欧美中文日本在线观看视频| 久久青草综合色| 国产三级在线视频| 丝袜人妻中文字幕| 国产成人免费无遮挡视频| 一a级毛片在线观看| 成人手机av| 色综合婷婷激情| 在线观看一区二区三区| 99在线视频只有这里精品首页| 国产日韩一区二区三区精品不卡| 国产精品影院久久| 精品国产乱码久久久久久男人| 久久久久精品国产欧美久久久| 久久久久久人人人人人| 亚洲九九香蕉| 亚洲欧洲精品一区二区精品久久久| 国产一级毛片七仙女欲春2 | 国产精品,欧美在线| 在线观看免费日韩欧美大片| 又黄又爽又免费观看的视频| 亚洲男人天堂网一区| 午夜老司机福利片| 91在线观看av| 女生性感内裤真人,穿戴方法视频| 欧美黄色淫秽网站| 亚洲五月色婷婷综合| 亚洲熟妇熟女久久| www.熟女人妻精品国产| 亚洲av成人不卡在线观看播放网| 午夜福利视频1000在线观看 | 黄片小视频在线播放| 欧美色欧美亚洲另类二区 | 国产精品久久视频播放| 国产1区2区3区精品| 久久精品aⅴ一区二区三区四区| 欧美人与性动交α欧美精品济南到| 亚洲专区国产一区二区| 欧美日韩亚洲综合一区二区三区_| 制服人妻中文乱码| 国产不卡一卡二| aaaaa片日本免费| 一级黄色大片毛片| 大陆偷拍与自拍| 亚洲av电影不卡..在线观看| 变态另类成人亚洲欧美熟女 | 精品一区二区三区四区五区乱码| 97人妻天天添夜夜摸| 日韩三级视频一区二区三区| 两性夫妻黄色片| 久热爱精品视频在线9| 国产三级黄色录像| 久久久久精品国产欧美久久久| 亚洲少妇的诱惑av| 国产精品久久视频播放| 欧美绝顶高潮抽搐喷水| 国产精品二区激情视频| bbb黄色大片| 日日爽夜夜爽网站| 亚洲五月色婷婷综合| 国产精品日韩av在线免费观看 | 精品不卡国产一区二区三区| 他把我摸到了高潮在线观看| 久久久久久亚洲精品国产蜜桃av| 久久精品91蜜桃| 日日干狠狠操夜夜爽| 久久久久精品国产欧美久久久| www.熟女人妻精品国产| 999久久久国产精品视频| 99国产精品一区二区蜜桃av| 黄色丝袜av网址大全| 女人被躁到高潮嗷嗷叫费观| 国产一区二区在线av高清观看| 欧美中文综合在线视频| 波多野结衣高清无吗| 亚洲精品久久国产高清桃花| 亚洲自偷自拍图片 自拍| 丰满的人妻完整版| 自拍欧美九色日韩亚洲蝌蚪91| 日本黄色视频三级网站网址| 欧美激情高清一区二区三区| 亚洲天堂国产精品一区在线| 国产高清视频在线播放一区| 露出奶头的视频| 午夜成年电影在线免费观看| 国产成人免费无遮挡视频| 久久国产乱子伦精品免费另类| 免费在线观看完整版高清| 最近最新免费中文字幕在线| 50天的宝宝边吃奶边哭怎么回事| 侵犯人妻中文字幕一二三四区| 亚洲熟妇中文字幕五十中出| 制服丝袜大香蕉在线| 婷婷丁香在线五月| 亚洲欧美日韩高清在线视频| 日韩成人在线观看一区二区三区| 黄频高清免费视频| 制服人妻中文乱码| 国产精品一区二区精品视频观看| 美女午夜性视频免费| 两个人免费观看高清视频| 国产av精品麻豆| 中文字幕人妻熟女乱码| 亚洲欧美日韩高清在线视频| www.自偷自拍.com| www.www免费av| 夜夜躁狠狠躁天天躁| 97人妻精品一区二区三区麻豆 | 99国产极品粉嫩在线观看| av在线播放免费不卡| 国产麻豆成人av免费视频| 在线观看www视频免费| 色播亚洲综合网| 天天一区二区日本电影三级 | 久久久久九九精品影院| 变态另类成人亚洲欧美熟女 | 国产成人精品久久二区二区免费| 国产成人系列免费观看| 国产单亲对白刺激| 波多野结衣av一区二区av| 操出白浆在线播放| 激情视频va一区二区三区| 露出奶头的视频| 制服人妻中文乱码| 涩涩av久久男人的天堂| 午夜精品在线福利| 成人国产综合亚洲| 国产一区二区三区综合在线观看| 亚洲人成伊人成综合网2020| 女性被躁到高潮视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久精品国产欧美久久久| 国产一区二区三区视频了| 午夜日韩欧美国产| 精品久久久精品久久久| 天堂影院成人在线观看| 欧美一级毛片孕妇| 精品卡一卡二卡四卡免费| 美女国产高潮福利片在线看| 久久天堂一区二区三区四区| 亚洲男人的天堂狠狠| 日韩欧美国产一区二区入口| 一边摸一边抽搐一进一小说| 好看av亚洲va欧美ⅴa在| 精品国产国语对白av| 老汉色av国产亚洲站长工具| 精品无人区乱码1区二区| 黄色 视频免费看| 黄频高清免费视频| 久久午夜综合久久蜜桃| 久久精品国产清高在天天线| 桃色一区二区三区在线观看| 男人舔女人的私密视频| 一本久久中文字幕| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 美女午夜性视频免费| 一个人观看的视频www高清免费观看 | 91麻豆av在线| 久久欧美精品欧美久久欧美| 老鸭窝网址在线观看| 人妻丰满熟妇av一区二区三区| 免费看a级黄色片| 国产精品久久久久久人妻精品电影| 日韩欧美三级三区| 黄色片一级片一级黄色片| 成人特级黄色片久久久久久久| 宅男免费午夜| 欧美黄色片欧美黄色片| 国产精品精品国产色婷婷| 国产伦人伦偷精品视频| 精品国产一区二区久久| www日本在线高清视频| 高清黄色对白视频在线免费看| 丝袜美足系列| 19禁男女啪啪无遮挡网站| av天堂久久9| 中文字幕最新亚洲高清| 国产又爽黄色视频| x7x7x7水蜜桃| 在线视频色国产色| 极品人妻少妇av视频| 久久久久久免费高清国产稀缺| 亚洲国产精品成人综合色| 69精品国产乱码久久久| 男男h啪啪无遮挡| 两个人免费观看高清视频| 亚洲国产精品久久男人天堂| 美女大奶头视频| 99riav亚洲国产免费| 97人妻天天添夜夜摸| 免费在线观看影片大全网站| 欧美人与性动交α欧美精品济南到| 变态另类成人亚洲欧美熟女 | 美女高潮喷水抽搐中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲专区中文字幕在线| 制服诱惑二区| 女人爽到高潮嗷嗷叫在线视频| 精品人妻1区二区| 一区在线观看完整版| 久久精品亚洲熟妇少妇任你| 日韩三级视频一区二区三区| 亚洲天堂国产精品一区在线| 一区福利在线观看| 亚洲av电影不卡..在线观看| 看免费av毛片| 亚洲熟妇中文字幕五十中出| 国产av一区二区精品久久| 久久久久精品国产欧美久久久| 人人妻,人人澡人人爽秒播| 99精品久久久久人妻精品| 女性被躁到高潮视频| 黄网站色视频无遮挡免费观看| 国产精华一区二区三区| 国产精品99久久99久久久不卡| a在线观看视频网站| 手机成人av网站| 亚洲av第一区精品v没综合| 午夜福利成人在线免费观看| 久久狼人影院| 色av中文字幕| av在线播放免费不卡| 欧美成人一区二区免费高清观看 | 国产精品秋霞免费鲁丝片| 日韩中文字幕欧美一区二区| 亚洲九九香蕉| 精品国产国语对白av| 国产成+人综合+亚洲专区| 国产精品二区激情视频| 久久精品国产清高在天天线| 999久久久精品免费观看国产| 久久人妻福利社区极品人妻图片| 黑人操中国人逼视频| 欧美激情 高清一区二区三区| 夜夜躁狠狠躁天天躁| 日韩三级视频一区二区三区| 亚洲av电影不卡..在线观看| 99热只有精品国产| 久久伊人香网站| 精品乱码久久久久久99久播| 久久人人爽av亚洲精品天堂| 一进一出好大好爽视频| aaaaa片日本免费| 少妇 在线观看| 18美女黄网站色大片免费观看| 国内久久婷婷六月综合欲色啪| av视频免费观看在线观看| 久久久久久久精品吃奶| 精品福利观看| 一级,二级,三级黄色视频| 777久久人妻少妇嫩草av网站| 久久中文字幕一级| 国产伦一二天堂av在线观看| 国产亚洲精品久久久久5区| av电影中文网址| 亚洲av熟女| 久久精品国产综合久久久| АⅤ资源中文在线天堂| 窝窝影院91人妻| 宅男免费午夜| 成人欧美大片| 日韩成人在线观看一区二区三区| 欧美国产日韩亚洲一区| 日本欧美视频一区| 黄色丝袜av网址大全| 波多野结衣av一区二区av| 自拍欧美九色日韩亚洲蝌蚪91| 巨乳人妻的诱惑在线观看| 久久久久国产一级毛片高清牌| av在线播放免费不卡| 午夜久久久久精精品| 亚洲狠狠婷婷综合久久图片| 少妇熟女aⅴ在线视频| 午夜精品在线福利| 999精品在线视频| 波多野结衣一区麻豆| 麻豆久久精品国产亚洲av| 在线观看免费午夜福利视频| 日本免费一区二区三区高清不卡 | 夜夜夜夜夜久久久久| 亚洲少妇的诱惑av| 美女 人体艺术 gogo| 波多野结衣av一区二区av| 亚洲第一欧美日韩一区二区三区| 亚洲av五月六月丁香网| 18禁美女被吸乳视频| 欧美日本视频| 免费人成视频x8x8入口观看| 亚洲精品一卡2卡三卡4卡5卡| 变态另类丝袜制服| 国产精品电影一区二区三区| 国产精品一区二区三区四区久久 | 亚洲成人免费电影在线观看| 禁无遮挡网站| 日本黄色视频三级网站网址| 一进一出抽搐动态| 黑人巨大精品欧美一区二区蜜桃| 亚洲五月婷婷丁香| 亚洲 欧美 日韩 在线 免费| 国产欧美日韩综合在线一区二区| 亚洲欧美日韩无卡精品| 国产高清videossex| 国产精品免费一区二区三区在线| 女人爽到高潮嗷嗷叫在线视频| 91av网站免费观看| 又紧又爽又黄一区二区| 国产熟女午夜一区二区三区| 日本五十路高清| 校园春色视频在线观看| 99热只有精品国产| 免费在线观看日本一区| 色播亚洲综合网| 亚洲精品一卡2卡三卡4卡5卡| 真人一进一出gif抽搐免费| 最好的美女福利视频网| 国产在线精品亚洲第一网站| 欧美激情极品国产一区二区三区| 国产精品电影一区二区三区| 国产午夜精品久久久久久| 午夜福利影视在线免费观看| 亚洲色图 男人天堂 中文字幕| 一进一出抽搐gif免费好疼| 国产精品久久电影中文字幕| www.999成人在线观看| 美女国产高潮福利片在线看| 桃红色精品国产亚洲av| 日本欧美视频一区| 久久久久久亚洲精品国产蜜桃av| 黄片小视频在线播放| 精品卡一卡二卡四卡免费| 大香蕉久久成人网| 亚洲成人久久性| 国产麻豆成人av免费视频| 国产精品香港三级国产av潘金莲| 亚洲男人的天堂狠狠| 在线观看66精品国产| 亚洲七黄色美女视频| 久久香蕉精品热| 精品一区二区三区四区五区乱码| 亚洲少妇的诱惑av| 夜夜躁狠狠躁天天躁| 欧美激情高清一区二区三区| 在线观看免费视频日本深夜| 狠狠狠狠99中文字幕| 成人国产一区最新在线观看|