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

    硼和氟摻雜ZSM-5分子篩的制備及其甲醇制丙烯反應(yīng)性能

    2020-09-10 07:41:00胡曉燕閆新龍武建軍
    關(guān)鍵詞:中國(guó)礦業(yè)大學(xué)化工學(xué)院建軍

    馮 銳 胡曉燕 閆新龍 武建軍

    (中國(guó)礦業(yè)大學(xué)化工學(xué)院,徐州 221116)

    Propylene,as one of the most important intermediates in petrochemical industry,is traditionally produced by the steam cracking of naphtha or ethane,as well as the fluid catalytic cracking(FCC)process.To meet the ever-increasing demand of propylene and reduce its dependence on crude oil resource,new strat-egies for propylene production from cheaper and/or renewable feedstocks have attracted an increasing interest in recent decades.Among the newly-developed technologies,methanol to olefins(MTO)and methanol to propylene(MTP)technologies have been paid much attention due to the easy availability of methanol and high selectivity to light olefins(C2=~C4=).Specifically,the MTP process can optimize the propylene yield over ZSM-5 catalyst due to its unique microporous structure and surface acidity.However,the propylene selectivity and catalytic lifetime need to be enhanced.Many studies have shown that the physiochemical properties of ZSM-5 zeolite such as surface acidity,pore structure,and crystal sizes,play important roles in improving propylene selectivity and catalyst lifetime.For example,the acid sites in ZSM-5 zeolites are the active sites not only for methanol conversion and propylene formation,but also the undesired secondary reactions such as hydrogen transfer,cyclization and aromatization[1].Generally,high-density acid sites induce the rapid deactivation of catalysts due to coke deposition;and the intrinsic micropores of ZSM-5 zeolites increase the diffusion resistance of propylene,the secondary reactions of propylene decrease its selectivity[2-3].Therefore,in the last decade,efforts have been devoted to the ZSM-5 zeolites:(a)synthesizing the hierarchical or nanosized ZSM-5 to increase the diffusion of propylene[3-4];(b)modifying the surface acid strength and acid density to decrease the adsorption and hydrogen transfer of propylene on acid sites[5,7].For example,Liu et al.[8]modified high-silicon ZSM-5 with 0.1% phosphorus to reduce the acid strength of Br?nsted acid sites,which increased the propylene selectivity by 10% in the MTP reaction.In our previous study,the Lewis acid sites were found to be important sites for coke formation in hydrocarbon cracking reactions[9-11].However,the effect of acid sites type of ZSM-5 zeolite in the methanol conversion process has not been fully considered yet.In this study,to elucidate the role of Lewis and Br?nsted acid sites on the activity and selectivity in the MTP reaction,ZSM-5 zeolites with varied acid properties were synthesized with boron and/or fluorine doping and tested in the MTP reaction.

    1 Experimental

    1.1 Raw materials

    Ammonium fluoroborate(NH4BF4,97.00%(w/w))was purchased from Shanghai SSS Reagent Co.,Ltd.Ammonium fluoride(NH4F),boric acid(H3BO3),tetraethyl orthosilicate(TEOS),aluminium isopropoxide,and sodium hydroxide(NaOH)with analytical purity werepurchasedfrom Sinopharm GroupChemical Reagent Co.,Ltd.Tetrapropyl ammonium hydroxide(TPAOH,25%(w/w)aqueous solution)was purchased from Sinopharm Group Chemical Reagent Co.,Ltd.

    1.2 Sample preparation

    By varying the synthesis temperature during the crystallization,the nucleation and growth process should be controlled to obtain ZSM-5 zeolites with small and uniform crystals[12].Here,boron and/or fluorine doped ZSM-5 zeolites were synthesized using a two-stage crystallization method.Appropriate amounts of TEOS,TPAOH,NaOH,and aluminum isopropoxide were mixed with deionized water and stirred at room temperature for 6 h.Then,NH4BF4,NH4F,or H3BO3was added into the above mixture as modifiers and stirred for another 12 h.The final gel compositions of synthesizing systemnadditive)were 50∶1.0∶2:8∶3 000∶x.Here,three samples were marked as Z5-BF1,Z5-BF2,Z5-BF3 when NH4BF4was used as additive withx=3,6,and 9,in turn;two samples were marked as Z5-F1 and Z5-F2 when NH4F was used as additive withx=12 and 36,respectively;two samples were marked as Z5-B1 and Z5-B2 when H3BO3was used as additive withx=3 and 9,respectively.In comparison,the sample was marked as Z5 whenx=0.The above gel was then transferred into an autoclave for pre-crystallization at 110℃for 3 h and further crystallization at 170℃for 48 h.Thereafter,the as-synthesized solid products were filtered,washed with deionized water,dried overnight,and calcined at 550℃for 4 h in air to remove the organic templates.The H-form ZSM-5 zeolites were finally obtained after twice ion-exchanges with 1.0 mol·L-1of NH4Cl solution and once calcination at 550℃for 4 h in air.

    1.3 Characterization

    The crystalline phases were measured by X-ray diffraction(XRD)on a Bruker D8 Advance Diffractometer,using a CuKαradiation(λ=0.154 06 nm)operating at 40 kV and 30 mA,and the scanning angle of 5°~50°.The surface area and pore volume of zeolites were derived from the nitrogen sorption curves on a Micrometrics Tristar 3000 analyzer.Solid-state magic angle rotation nuclear magnetic resonance(MAS NMR)spectra were obtained from a Bruker AdvanceⅢHD 600 MHz instrument.The29Si NMR spectra were used to analyze the framework silica-to-alumina ratios(nSiO2/nAl2O3)of zeolites.Acid sites of all samples were measured on a Thermo Nicolet iS5 Fourier transform infrared(FT-IR)spectrometer,using pyridine as a probe molecule.Acid strength distributions were measured by NH3temperature-programmed desorption(NH3-TPD)on a Quanta chrome ChemStarTMinstrument.About 100 mg of samples were activated in flowing Ar at 350℃for 1 h,then cooled to room temperature and exposed to flow of 6.5%(V/V)of NH3/Ar gas.Physically adsorbed NH3was removed in Ar flow at 100℃for 2 h before collecting data at temperature ramped up to 650℃ at a rate of 10°C·min-1in helium flow of 30 mL·min-1.Field emission scanning electron microscopy(SEM)images were collected by a Quanta 250 instrument with an acceleration voltage of 15 kV.

    1.4 Evaluation of MTP performance

    The methanol to propylene reaction tests were performed in a fixed bed of a quartz tubular reactor with inner diameter of 6 mm.ZSM-5 catalysts with particle sizes of 40~60 mesh were prepared by pressing,crushing and sieving process.100 mg of catalyst was placed in the constant temperature zone of the reactor and pretreated at 550℃for 1 h under Ar gas before the reaction.Methanol was injected into the reactor by a constant flux pump using a 50 mL·min-1of Ar as carrier gas.The reaction temperature was set to 450℃and the weight hourly space velocity(WHSV)of methanol was 4.0 h-1.The reaction products were analyzed by an on-line gas chromatograph(GC-2014C,Shimadzu GC),which is equipped with a TG-BONG Q column,a thermal conductivity detector(TCD)and a flame ionization detector(FID).The maximum absolute error(<4%)was measured by triplicate experiments.The methanol conversion was defined as usual:XMeOH=(mMeOH,in-mMeOH,out)/mMeOH,in,while selectivity toiproduct was defined as follows:Si=mi/∑mi,wheremiis the mass of compoundi,and∑miis the total mass of all products.

    2 Results and discussion

    The X-ray diffraction(XRD)pattern in Fig.1 shows the characteristic peaks of MFI structure of ZSM-5 zeolite,corresponding to the(101),(020),(501),(151)and(303)crystal faces,respectively(PDF No.01-085-1208)[13-14].It shows that all samples exhibited a typical MFI phase and no peaks of crystallized impurity appeared.Fig.1b shows the shifts of two diffraction peaks to higher 2θangles with increasing the NH4BF4usages,indicative of the increase of framework SiO2/Al2O3ratios(nSiO2/nAl2O3).The diffraction peak shift of H3BO3modified samples such as Z5-B2 was more obvious than NH4F modified ones,indicating that H3BO3does more to increase the framework SiO2/Al2O3ratio.It might be attributed to the substitution of Al atoms by B atoms,since an increased B content was detected with increasing the boron usages[15-16].The relative crystallinities(RC)were calculated by comparing the peak areas of the characteristic peaks at 22°~25°and listed in Table 1.It shows that the crystallinity of ZSM-5 zeolites changed a little with heteroatom doping,and the relative crystallinity values of ZSM-5 zeolites decreased with increasing the NH4BF4usages.

    The29Si NMR spectra in Fig.2 exhibits a strong resonance at about-113(chemical shift,the same below),which corresponds to the Si(OSi)4species in the zeolites framework.A weak shoulder resonance at about-106 was ascribed to the Si(OSi)3(AlO)species.A weak resonance at about-102 was ascribed to the(SiO)3SiOH species.No resonance was observed below-100,indicating that there were no other silicate species[17-18].The framework SiO2/Al2O3ratios were calculated from the deconvoluted profileusing a Gaussian-Lorentzian mixed function and listed in Table 1.It shows that the variation of framework SiO2/Al2O3ratios was in consistent with that derived from XRD pattern.The doping of B and/or F increased the framework SiO2/Al2O3ratios compared with Z5 sample.The27Al NMR spectra of four selected ZSM-5 zeolites in Fig.3a exhibits two distinct resonances at 52 and-3,corresponding to a tetra-coordinated framework and a hexa-coordinated extra-framework of the Al species,respectively.No shoulderresonance ofa pentacoordinated Al species was observed at approximately 35.The results show that combining F and B for Z5-BF2 remarkably decreases the percentage of the extraframework Al species,but the percentages of extraframework Al species with B and F alone changed small.

    Fig.1 XRD patterns of as-synthesized ZSM-5 zeolites

    Fig.229Si NMR spectra of as-synthesized ZSM-5 zeolites with revolved resonances

    Table 1 Physiochemical properties of as-synthesized ZSM-5 zeolites

    The11B NMR spectra in Fig.3b displays two main peaks at-3.9 and 13.8.The former resonance is attributed to the tetrahedral coordinated framework boron,and the latter broad band is assigned to the trigonal coordinated extraframework boron species[15,19].The small shoulder between these two bands is ascribed to the second-order quadrupole broadening of the trigonal coordinated framework boron,which is easily derived from the dehydration of tetrahedral coordinated framework boron species[16].This result confirms that the majority of B species is incorporated into the framework of zeolite.The19F NMR spectra in Fig.3c shows that there are five peaks for zeolites Z5-BF2 and Z5-F1.The peaks from-66 to-87 are attributed to F-anions located in the cages of silica zeolites[20-21].The peak at around-123 is assigned to the presence of the F-ion in zeolite channels as counter ion of balanced cations such as H+,Na+or NH4+[22-23].The peaks at-141 and-142 are associated with the extraframework Al species such as AlF63-.The peaks at-158 and-161 are related to Si-F groups which are caused by the replacement of hydroxyl group in the silanols by fluo-rine atom,or a terminal fluorine atom at the surface,indicative of the fluorination in the form of(SiO3)Si-F groups[9,23].The peaks at-177 and-180 are assigned to the partially hydrated AlF3species[24-25].The NMR data indicate that the boron and fluorine exist in the form of framework incorporation and chemical bonding with extraframework Al species.

    Fig.3 27Al NMR spectra(a),11B NMR spectra(b),and19F NMR spectra(c)of selected as-synthesized ZSM-5 zeolites

    Fig.4 shows the N2sorption isotherms and pore size distributions of the ZSM-5 zeolites.The isotherms are identified as typeⅣ,which is characteristic of micropore and mesopore mixed structure of ZSM-5 materials.The Brunauer-Emmet-Teller(BET)specific surface areas of the samples were calculated from N2isotherms and listed in Table 2.It shows that with increasing the NH4BF4usages,microporous surface areas of ZSM-5 zeolites increased,but their mesoporous surface areas decreased obviously.As a result,the total surface areas gradually decreased,from 322 m2·g-1for Z5 to 272 m2·g-1for Z5-BF3.The B and F doping contributes to the formation of microspore but impaired the mesopore,being in consistent with that of H3BO3and NH4F modified samples.The same tendency has also been observed for the pore volumes,i.e.,the micropore volumes of ZSM-5 zeolites containing B and F were higher than that of Z5.The pore size distribution obtained from the desorption branch indicates the reduction of mesoporous structure for B and F doped samples,particularly for the F doped samples.

    Fig.4 N2sorption isotherms(a)and pore size distributions(b)of as-synthesized ZSM-5 zeolites

    Table 2 Textural properties of as-synthesized ZSM-5 zeolites

    FT-IR analysis was carried out to characterize the acid amounts of ZSM-5 zeolites.As highlighted in Fig.5a,it is generally recognized that the bands at 1 447 and 1 546 cm-1are ascribed to the pyridine molecules coordinated to Lewis acid sites and pyridinium ions formed by protonation of pyridine on Br?nsted acid sites,respectively[10,26-27].The band at 1 490 cm-1is assigned to pyridine associated with both the two kinds of acid sites.The band at 1 597 cm-1is ascribed to the hydrogen bonded pyridine.The acid amount was calculated using a semi-quantitative method proposed by Emeis et al[28].The results in Table 3 show that with the increase of NH4BF4usages,the amounts of Br?nsted acid sites of Z5-BFxobviously increased;in contrast,the amounts of Lewis acid sites remarkably decreased.It is caused by the doping of boron and fluorine atoms,which could effectively increase the amount of Br?nsted acid sites but decrease the amount of Lewis acid sites,in view of the results of H3BO3and NH4F used alone.

    The FTIR spectra of OH stretching region in Fig.5b is used to characterize the effect of B and F doping on OH groups.The band at 3 742 cm-1is ascribed to the OH stretching vibration(nOH)of free silanols which is thought to be located at the external surface of zeolite,namely SiOHs[26,29].It shows that combing B and F increases the amount of external surface SiOHs which is indicative of the external surface dealumination,being in consistent with the results of27Al NMR analysis in Fig.3.The bands at 3 691 and 3 604 cm-1are assigned to the vibration of extraframework Al-OH and bridging Si-OH-Al groups,respectively[30-31].Obviously,B and/or F doping could decrease the amounts of extraframework Al-OH and bridging Si-OH-Al species.Those species play important roles in constituting the acid sites,especially Br?nsted acid sites,which is consistent with the variation tendency of Br?nsted acid amounts with increasing the heteroatom doping(Table3).

    Fig.5 FT-IR spectra of pyridine absorbed by as-synthesized ZSM-5 zeolites

    Table 3 Surface acidity of as-synthesized ZSM-5 zeolites

    Fig.6 NH3-TPD profiles of as-synthesized ZSM-5 zeolites

    NH3-TPD tests were also performed to characterize the acid strength distributions of ZSM-5 zeolites.As shown in Fig.6,all samples exhibited two resolved desorption peaks:the low-temperature around 250℃and high-temperature peak higher than 300℃,corresponding to the weak and strong acid sites,respectively.It is obvious that the peak of strong acid sites on B and/or F doped samples shifted to higher temperatures compared to Z5.The relative amounts of acid sites derived from the resolved overlapped curves are listed in Table 3.It can be seen that the weak acid sites of B and/or F doped samples decreased compared with Z5 sample,which is consistent with the change of Lewis acid amount;however,the change of strong acid amount is in keeping with the Br?nsted acid amount(Table 3).The decrease of Lewis acid sites is probably attributed to the isomorphous substitution of Al atom by B atom or the bond combination of B and Al atoms;in contrast,the strong electronegativity of F atom reduces the electron cloud on Si-OH-Al group,enhancing the strength of Br?nsted acid sites[9-10,32].

    Fig.7 SEM images of as-synthesized ZSM-5 zeolites

    SEM images in Fig.7 show that all samples possessed relatively uniform particles.The pre-crystallization at low temperature contributes to the formation of crystal nucleus and further crystallization at high temperature are beneficial to the synthesis of small sized ZSM-5 zeolite with good uniformity in a short time.Sample Z5 had well dispersed spherical shape parti-clesca.2 μm.With increasing the NH4BF4usages,the particle sizes of Z5-BFxdecreased gradually.Z5-BF1 had cubic and round cornered particles with diameter ofca.500 nm,while Z5-BF3 had the smallest particles with diameter of approximately 200 nm.Compared to NH4BF4,H3BO3and NH4F used alone had the reverse effect on the particle size of ZSM-5 zeolites.At low usage of modifiers,the particle sizes of Z5-B1 and Z5-F1 were smaller than sample Z5,while they were much larger than Z5 at higher usage.Notably,high concentration of fluorine alone could not only accelerate the crystallization rate of zeolite but also etch the crystal surface(Fig.7g),as previously reported[33-34].In contrast,boron alone made the ZSM-5 particle surface coarser than Z5(Figs.7h and 7i).

    TEM images of typical Z5 and Z5-BF2 samples is shown in Fig.8.Both samples possessed relatively uniform crystals with smooth surface,being in consistent with that measured by SEM images.Some cavities are observed on the border and inside the crystals,indicating that the mesoporous structures detected in the N2sorption isotherms are mainly intra-crystal mesopores.

    The products distributions of all ZSM-5 catalysts in MTP reaction are listed in Table 4.All ZSM-5 catalysts exhibited nearly full methanol conversion during the initial reaction period,indicative of their high initial activity.By defining 85% of methanol conversion as the point of catalyst deactivation,the boron and/or fluorine doped catalysts showed longer catalytic lifetime than Z5.With the increases of NH4BF4usages,the lifetimes of Z5-BFxincreased first and then decreased.Z5-BF2 possessed the longest lifetime of 23 h.By contrast,catalysts with boron used alone had better performance in improving the catalyst lifetime compared with that used fluorine alone.As for product selectivity,the Z5-BF2 had the highest propylene selectivity up to 41.5%.In contrast,the propylene selectivity on catalysts with fluorine alone decreased greatly compared with Z5 catalyst.

    Methanol conversion and selectivity to light olefins(C2=,C3=and C4=)as a function of time on stream(TOS)over three selected catalysts are displayed in Fig.9.It shows that the selectivity of light olefins(C2=~C4=)on ZSM-5 catalysts was high at the initial period,and then decreased to a stable stage with slight fluctuation.Of all catalysts,Z5-BF2 possessed the highest average selectivity to propylene and light olefins,up to 41.5% and 63.6%,respectively.In the methanol conversion on zeolites,aromatic-based cycle and olefinbased cycle simultaneously occur,and ethylene mainly comes from the former while propylene and other heavi-er olefins mainly come from the latter[35].Taking Z5 and Z5-BF2 for examples,decreasing acid site density and the amount of Lewis acid sites could reduce the rates of aromatization and hydrogen transfer,resulting in lower ethylene selectivity but higher propylene selectivity for Z5-BF2.In contrast,the higher ethylene selectivity of Z5-BF3 catalyst is attributed to its improved amount of Br?nsted acid sites,compared to Z5 catalyst.Similarly,high acid site density and high amount of Lewis acid site could increase the selectivity of C5+products(mainly aromatic species),due to the aromatization reactions[18,36].

    Fig.8 TEM images of Z5(a,b)and Z5-BF2(c,d)

    Table 4 MTP reaction results on H-form ZSM-5 catalysts with boron and/or fluorine dopinga

    Fig.9 Conversion and product selectivity for the MTP reaction vs time-on-stream(TOS)over catalysts:(a)Z5;(b)Z5-BF2;(c)Z5-B1

    From the foregoing,boron alone could improve the lifetime of ZSM-5 catalyst while fluorine alone had limited effect on prolonging the catalyst lifetime.The combination of boron and fluorine on Z5-BFxcatalysts improved both catalyst lifetime and selectivity of light olefins,propylene in particular.Generally,the deactivation of MTP catalysts is caused by the accumulation of coke species on the surface of ZSM-5 zeolites,which cover the surface acid sites or plug pores and suppress the access of reactants to acid sites[7,37-38].However,it seems that mesopore did not do much in reducing the coke deactivation,considering that the decrease of mesoporous volumes of Z5-BFxdid not affect their lifetimes compared with Z5 catalyst.

    The crystal size is important in improving propylene selectivity and catalyst lifetime[39-42].Generally,nanosized ZSM-5 zeolites have been reported to show superior catalytic stability in the MTP reaction due to containing short diffusion paths.Despite prolonging the catalytic lifetime for nanosized ZSM-5 zeolites in most cases,their exposed external surface acid sites are more likely to absorb light olefins and coke precursors and contribute little to improving propylene selectivity because of their weakened shape selectivity.In this study,with increasing the NH4BF4usages,the particle sizes of Z5-BFxdecreased gradually.However,propylene selectivity and catalyst lifetime did not monotonically increase or decrease with the variation of crystal sizes.The superior catalytic performance of Z5-BF2 more reasonably results from the optimization of acid sites in view of the significant changes of surface acidity compared to Z5 in Table 3.

    By contrast,acid sites,especially the strong acid sites,are thought to be the predominant reason for coke deposition[38,43-46].The strong adsorption of light olefins on strong acid sites(the Br?nsted acid sites in this work)and further polymerization with other light olefins should be responsible for the catalyst deactivation[2,39].In view of the acid properties of catalysts in Table 3,the decrease of Lewis acid amount of Z5-BF1 and Z5-BF2 obviously improved the catalyst lifetimes compared with Z5 catalyst;while the sharp increase of Br?nsted acid amount of Z5-BF3 reduced the catalyst lifetime compared with Z5-BF2.Therefore,Lewis acid sites and strong acid sites,especially strong Br?nsted acid sites,are unfavorable due to the fast coking deactivation,as previously reported[11,47-48].The amount of strong acid sites on ZSM-5 catalyst should be controlled to a relatively low level.

    Some references related with boron or fluorine modified ZSM-5 and their catalytic performances in the MTP reactions are listed in Table 5.In the previous studies,most of the B or F modified ZSM-5 catalysts possessed higher propylene selectivity and longer catalyst lifetimes than unmodified ZSM-5 catalysts.As shown in Table 5,all catalysts had relatively high propylene selectivity more than 40%,however,their catalyst lifetimes differed widely.The big differences of cat-alyst lifetimes are caused by many factors such as the acidity,pore structure,and crystal sizes.It is noted that the reaction conditions such as catalyst loading amount,water addition into methanol feed,and even inner diameter of the reactor are also important factor influencing the results.In our two studies[36,49],the commercial MTP catalysts which was tested under two different reaction conditions exhibited totally different catalytic lifetimes(25 h vs>400 h).In this study,we used the prepared ZSM-5 catalyst in the same reactor with that used in our previous studies[3,36].23 h of catalyst lifetime was comparable to that of the commercial MTP catalyst.More importantly,we believe that it is more meaningful to compare the catalytic performances of catalysts in the same reactor and same reaction conditions.

    Table 5 Comparison of MTP performance on B or F modified ZSM-5 catalysts

    3 Conclusions

    A series of ZSM-5 catalysts were successfully synthesized in the presence of boron and/or fluorine containing reagents.The structural and morphological characterization results reveal that the appropriate amount of NH4BF4increases the micropores but reduces the Lewis acid sites without obvious increase of Br?nsted acid sites.Compared to the conventional Z5 catalyst,Z5-BF2 exhibits a superior catalytic performance in the methanol to propylene reaction with higher propylene selectivity and longer lifetime.Moreover,reducing Lewis acid sites and controlling the acid strength of Br?nsted acid sites is expected to be suitable for improving the performance of MTP reaction.

    Acknowledgements:This work was financially supported by the Fundamental Research Funds for the Central Universities(Grant No.2018QNB04).

    猜你喜歡
    中國(guó)礦業(yè)大學(xué)化工學(xué)院建軍
    使固態(tài)化學(xué)反應(yīng)100%完成的方法
    慶祝建軍95周年
    國(guó)家開放大學(xué)石油和化工學(xué)院學(xué)習(xí)中心列表
    【鏈接】國(guó)家開放大學(xué)石油和化工學(xué)院學(xué)習(xí)中心(第四批)名單
    Spatio-temporal evolution characteristics and pattern formation of a gas–liquid interfacial AC current argon discharge plasma with a deionized water electrode
    無論等多久
    《化工學(xué)報(bào)》贊助單位
    高校學(xué)生評(píng)教的問題與對(duì)策——以中國(guó)礦業(yè)大學(xué)為例
    中國(guó)礦業(yè)大學(xué)教育培訓(xùn)工作簡(jiǎn)介
    慶祝建軍八十三周年
    久久人人精品亚洲av| 精品一区二区三区av网在线观看| 日韩欧美精品免费久久 | 亚洲国产高清在线一区二区三| 少妇高潮的动态图| 亚洲中文日韩欧美视频| 精品一区二区三区av网在线观看| 免费观看精品视频网站| 中文在线观看免费www的网站| 免费看日本二区| 一进一出抽搐动态| 欧美日韩一级在线毛片| 精品熟女少妇八av免费久了| 麻豆国产97在线/欧美| 久久久久久国产a免费观看| 亚洲18禁久久av| 亚洲18禁久久av| 美女黄网站色视频| 亚洲精品美女久久久久99蜜臀| 99热只有精品国产| 18禁国产床啪视频网站| 国产精品久久久人人做人人爽| 村上凉子中文字幕在线| 老汉色∧v一级毛片| 757午夜福利合集在线观看| 亚洲无线观看免费| 成年版毛片免费区| 亚洲精品一卡2卡三卡4卡5卡| 在线国产一区二区在线| 麻豆一二三区av精品| 淫妇啪啪啪对白视频| 午夜精品一区二区三区免费看| 国产精品 欧美亚洲| 日本一本二区三区精品| 首页视频小说图片口味搜索| 一边摸一边抽搐一进一小说| 制服丝袜大香蕉在线| 国产成人福利小说| 国产成年人精品一区二区| 午夜a级毛片| 日韩欧美精品v在线| 老汉色∧v一级毛片| 老汉色∧v一级毛片| 99国产极品粉嫩在线观看| 我的老师免费观看完整版| 欧美日韩综合久久久久久 | 有码 亚洲区| 十八禁网站免费在线| 亚洲精品日韩av片在线观看 | 国产真人三级小视频在线观看| 欧美+亚洲+日韩+国产| 国产三级中文精品| 国内精品美女久久久久久| 亚洲成av人片在线播放无| 亚洲天堂国产精品一区在线| 一级黄片播放器| 天天躁日日操中文字幕| 黄色视频,在线免费观看| 国产一区二区三区在线臀色熟女| 两性午夜刺激爽爽歪歪视频在线观看| 小蜜桃在线观看免费完整版高清| 91在线观看av| 成人永久免费在线观看视频| 在线观看舔阴道视频| 可以在线观看毛片的网站| 最近最新中文字幕大全电影3| 国产精品一区二区免费欧美| 精品一区二区三区av网在线观看| 久久久久精品国产欧美久久久| 高清在线国产一区| www.熟女人妻精品国产| 怎么达到女性高潮| 两个人看的免费小视频| 国产色婷婷99| 国产午夜精品论理片| 黄色日韩在线| 亚洲av二区三区四区| 国产精品国产高清国产av| 亚洲av成人精品一区久久| svipshipincom国产片| 午夜精品久久久久久毛片777| 精品免费久久久久久久清纯| 久久久久久久久久黄片| 日韩 欧美 亚洲 中文字幕| 亚洲精品日韩av片在线观看 | 五月伊人婷婷丁香| 免费看a级黄色片| 午夜福利18| 亚洲人成网站在线播放欧美日韩| 麻豆成人av在线观看| 国产色爽女视频免费观看| 欧美色欧美亚洲另类二区| 禁无遮挡网站| 最近最新中文字幕大全免费视频| 国产亚洲欧美在线一区二区| 香蕉av资源在线| 日本与韩国留学比较| 免费在线观看影片大全网站| 国产黄片美女视频| 人人妻人人看人人澡| 小蜜桃在线观看免费完整版高清| 99久久99久久久精品蜜桃| 人妻丰满熟妇av一区二区三区| 久久国产精品人妻蜜桃| 国产精品一区二区免费欧美| 国产成人av激情在线播放| 精品99又大又爽又粗少妇毛片 | 中文亚洲av片在线观看爽| 国产毛片a区久久久久| 天天一区二区日本电影三级| 欧美最黄视频在线播放免费| 法律面前人人平等表现在哪些方面| 特大巨黑吊av在线直播| e午夜精品久久久久久久| 国产精品嫩草影院av在线观看 | 国产v大片淫在线免费观看| 女人被狂操c到高潮| x7x7x7水蜜桃| 成人av在线播放网站| 国产精品美女特级片免费视频播放器| 桃色一区二区三区在线观看| 亚洲人成电影免费在线| 久久性视频一级片| 成年免费大片在线观看| 黄色丝袜av网址大全| 看免费av毛片| 一本综合久久免费| 日韩精品中文字幕看吧| 亚洲精品在线美女| 成人国产一区最新在线观看| 久久久久亚洲av毛片大全| 婷婷六月久久综合丁香| 亚洲av成人av| 亚洲欧美一区二区三区黑人| 乱人视频在线观看| 国产主播在线观看一区二区| 亚洲va日本ⅴa欧美va伊人久久| 小蜜桃在线观看免费完整版高清| 亚洲人成网站高清观看| 亚洲男人的天堂狠狠| 热99re8久久精品国产| av女优亚洲男人天堂| 欧美成狂野欧美在线观看| 18美女黄网站色大片免费观看| 夜夜爽天天搞| 亚洲五月天丁香| 亚洲av日韩精品久久久久久密| 97人妻精品一区二区三区麻豆| 亚洲精品成人久久久久久| 男女做爰动态图高潮gif福利片| 在线视频色国产色| 国产高清视频在线播放一区| 伊人久久大香线蕉亚洲五| 国产高清激情床上av| 高清日韩中文字幕在线| 天堂√8在线中文| 日韩欧美免费精品| 91在线精品国自产拍蜜月 | 久久久精品大字幕| 99国产精品一区二区三区| 欧美一区二区精品小视频在线| 九九久久精品国产亚洲av麻豆| 18禁在线播放成人免费| 男人和女人高潮做爰伦理| 国内精品久久久久久久电影| 亚洲avbb在线观看| 亚洲成av人片在线播放无| 色哟哟哟哟哟哟| av国产免费在线观看| 无遮挡黄片免费观看| 黄片小视频在线播放| 在线观看美女被高潮喷水网站 | www日本黄色视频网| 88av欧美| 观看美女的网站| 日本一二三区视频观看| 国产单亲对白刺激| 欧美最黄视频在线播放免费| 国产精品1区2区在线观看.| 精品一区二区三区视频在线观看免费| 欧美日韩黄片免| 精品国产三级普通话版| 床上黄色一级片| 狠狠狠狠99中文字幕| 99久久无色码亚洲精品果冻| 亚洲精品久久国产高清桃花| 男人和女人高潮做爰伦理| 欧美中文日本在线观看视频| 午夜久久久久精精品| 人人妻人人看人人澡| 久久久久久人人人人人| 日韩 欧美 亚洲 中文字幕| 男女床上黄色一级片免费看| 一夜夜www| 午夜免费激情av| 精品日产1卡2卡| 日韩大尺度精品在线看网址| 手机成人av网站| 久久精品人妻少妇| 亚洲欧美日韩无卡精品| 国产精品日韩av在线免费观看| а√天堂www在线а√下载| 成人性生交大片免费视频hd| 啦啦啦观看免费观看视频高清| 午夜免费男女啪啪视频观看 | 熟女电影av网| 在线观看舔阴道视频| 国产精品99久久99久久久不卡| 美女大奶头视频| 亚洲狠狠婷婷综合久久图片| 男插女下体视频免费在线播放| 少妇的丰满在线观看| 精品国内亚洲2022精品成人| 国产探花极品一区二区| 熟女电影av网| 国产三级中文精品| 精品乱码久久久久久99久播| 国产精品三级大全| 在线观看午夜福利视频| 色视频www国产| 国产高清三级在线| 国产亚洲精品一区二区www| 久久久久国产精品人妻aⅴ院| 黄色成人免费大全| 日韩欧美国产一区二区入口| www.色视频.com| 午夜a级毛片| 十八禁网站免费在线| 99热只有精品国产| 国产一区二区三区在线臀色熟女| 岛国在线免费视频观看| 美女 人体艺术 gogo| 制服人妻中文乱码| 好看av亚洲va欧美ⅴa在| 日韩欧美国产一区二区入口| 久久久久久久久久黄片| 最后的刺客免费高清国语| 18禁黄网站禁片午夜丰满| 精品人妻一区二区三区麻豆 | 高清毛片免费观看视频网站| 午夜激情福利司机影院| 成人鲁丝片一二三区免费| av黄色大香蕉| 国产亚洲精品综合一区在线观看| 99久久99久久久精品蜜桃| 在线十欧美十亚洲十日本专区| 成人av一区二区三区在线看| 搡老岳熟女国产| 国产精品日韩av在线免费观看| 亚洲精品久久国产高清桃花| 人妻久久中文字幕网| 成人精品一区二区免费| 亚洲精品在线美女| 在线视频色国产色| 国产一区二区激情短视频| 青草久久国产| 久久久久久久亚洲中文字幕 | 国产欧美日韩一区二区三| 精品久久久久久久毛片微露脸| 久久天躁狠狠躁夜夜2o2o| 色视频www国产| 国产激情偷乱视频一区二区| 一级毛片高清免费大全| 成年版毛片免费区| 久久久久免费精品人妻一区二区| 国产亚洲精品一区二区www| 亚洲精品一卡2卡三卡4卡5卡| 日韩欧美国产一区二区入口| 国产免费av片在线观看野外av| 深爱激情五月婷婷| 日韩亚洲欧美综合| 久久久久久久亚洲中文字幕 | 小蜜桃在线观看免费完整版高清| 午夜两性在线视频| 久久亚洲精品不卡| 一个人免费在线观看电影| 熟女人妻精品中文字幕| 午夜精品一区二区三区免费看| 一区二区三区激情视频| 欧美午夜高清在线| 最近最新免费中文字幕在线| 亚洲第一欧美日韩一区二区三区| 可以在线观看毛片的网站| 免费无遮挡裸体视频| 在线观看一区二区三区| 日本a在线网址| 国产成人av教育| 两个人看的免费小视频| 国语自产精品视频在线第100页| 久久精品人妻少妇| 国内揄拍国产精品人妻在线| 国内毛片毛片毛片毛片毛片| 97超级碰碰碰精品色视频在线观看| 国产精品一区二区免费欧美| 国产乱人伦免费视频| 日本与韩国留学比较| 亚洲国产欧美人成| 欧美日韩一级在线毛片| av黄色大香蕉| 亚洲欧美日韩东京热| 亚洲精品亚洲一区二区| 九色成人免费人妻av| 日本 欧美在线| 日本撒尿小便嘘嘘汇集6| 亚洲欧美日韩高清在线视频| 国产老妇女一区| 在线观看一区二区三区| 日本免费a在线| 十八禁人妻一区二区| 亚洲av不卡在线观看| 日韩国内少妇激情av| 亚洲精品一区av在线观看| 成年女人永久免费观看视频| 欧美日韩亚洲国产一区二区在线观看| 欧美不卡视频在线免费观看| 日韩欧美三级三区| 日本黄色片子视频| 成人特级黄色片久久久久久久| 在线观看一区二区三区| 精品熟女少妇八av免费久了| 欧美精品啪啪一区二区三区| 一区福利在线观看| 熟女人妻精品中文字幕| 看黄色毛片网站| 3wmmmm亚洲av在线观看| 一夜夜www| 黄色丝袜av网址大全| 日本黄大片高清| 九九热线精品视视频播放| 久99久视频精品免费| 国产精品99久久久久久久久| 国产伦精品一区二区三区视频9 | 首页视频小说图片口味搜索| 欧美不卡视频在线免费观看| 深夜精品福利| 国产美女午夜福利| 手机成人av网站| 亚洲18禁久久av| 欧美黄色淫秽网站| 国产精品亚洲一级av第二区| 床上黄色一级片| 欧美成人免费av一区二区三区| 久久香蕉精品热| 欧美一区二区亚洲| 欧美丝袜亚洲另类 | 欧美一区二区亚洲| 亚洲国产中文字幕在线视频| 亚洲中文字幕一区二区三区有码在线看| 哪里可以看免费的av片| 免费看日本二区| 亚洲中文字幕日韩| 中文字幕熟女人妻在线| 国内久久婷婷六月综合欲色啪| 日韩免费av在线播放| 偷拍熟女少妇极品色| 成人18禁在线播放| ponron亚洲| 国产精品精品国产色婷婷| 欧洲精品卡2卡3卡4卡5卡区| 我的老师免费观看完整版| 欧美性猛交╳xxx乱大交人| 国产精品亚洲av一区麻豆| 一个人看视频在线观看www免费 | 免费看日本二区| 波野结衣二区三区在线 | 国产精品美女特级片免费视频播放器| 日韩国内少妇激情av| 免费电影在线观看免费观看| 欧美色视频一区免费| 国产伦在线观看视频一区| 国产免费av片在线观看野外av| av女优亚洲男人天堂| 91九色精品人成在线观看| av国产免费在线观看| 国模一区二区三区四区视频| 久久久久久久亚洲中文字幕 | 亚洲国产精品999在线| 亚洲成av人片免费观看| 欧美绝顶高潮抽搐喷水| 88av欧美| 亚洲欧美激情综合另类| 精品一区二区三区视频在线观看免费| 90打野战视频偷拍视频| 国产精品女同一区二区软件 | av国产免费在线观看| 老鸭窝网址在线观看| 日本在线视频免费播放| 一级黄色大片毛片| 成年女人永久免费观看视频| 国产在线精品亚洲第一网站| 久久国产乱子伦精品免费另类| 性色avwww在线观看| 国产高清三级在线| 国产精品久久久久久亚洲av鲁大| 日本免费a在线| 人妻丰满熟妇av一区二区三区| 亚洲欧美日韩无卡精品| 亚洲欧美日韩高清专用| 99国产精品一区二区三区| 日本精品一区二区三区蜜桃| 欧美成人免费av一区二区三区| 老汉色∧v一级毛片| 欧美成人a在线观看| 久久精品亚洲精品国产色婷小说| 久久香蕉国产精品| 亚洲人成网站在线播| 亚洲成人久久性| 国产伦人伦偷精品视频| 制服人妻中文乱码| 啦啦啦免费观看视频1| 女警被强在线播放| 欧美绝顶高潮抽搐喷水| 99国产极品粉嫩在线观看| 免费观看的影片在线观看| 757午夜福利合集在线观看| 亚洲av日韩精品久久久久久密| 久久久久久久亚洲中文字幕 | 在线免费观看的www视频| 国产日本99.免费观看| 99视频精品全部免费 在线| 久久久国产精品麻豆| 久久国产精品人妻蜜桃| 97人妻精品一区二区三区麻豆| 老鸭窝网址在线观看| 99久久综合精品五月天人人| 三级男女做爰猛烈吃奶摸视频| 国产精品女同一区二区软件 | 久久香蕉国产精品| 欧洲精品卡2卡3卡4卡5卡区| 97碰自拍视频| 成年人黄色毛片网站| 免费无遮挡裸体视频| 舔av片在线| 欧美一级a爱片免费观看看| 久久精品国产亚洲av香蕉五月| 亚洲美女黄片视频| 一进一出好大好爽视频| 欧美日韩福利视频一区二区| 国语自产精品视频在线第100页| 听说在线观看完整版免费高清| 国产在线精品亚洲第一网站| 久久国产乱子伦精品免费另类| 国产精品嫩草影院av在线观看 | 国产成人福利小说| 日本三级黄在线观看| 2021天堂中文幕一二区在线观| 亚洲最大成人手机在线| 老熟妇乱子伦视频在线观看| 一个人免费在线观看的高清视频| 高清在线国产一区| 国产男靠女视频免费网站| 亚洲黑人精品在线| 国产精品野战在线观看| 婷婷六月久久综合丁香| 欧美黄色淫秽网站| 日本黄大片高清| 老鸭窝网址在线观看| 波多野结衣巨乳人妻| 久久久成人免费电影| 狂野欧美激情性xxxx| 夜夜看夜夜爽夜夜摸| 久久久久久久精品吃奶| 国产精品99久久久久久久久| 亚洲国产欧洲综合997久久,| 免费看日本二区| 国产高潮美女av| 性色avwww在线观看| 内射极品少妇av片p| АⅤ资源中文在线天堂| 日本在线视频免费播放| 亚洲精品在线观看二区| 日本一二三区视频观看| 国产视频一区二区在线看| 一个人免费在线观看电影| 午夜日韩欧美国产| 国产国拍精品亚洲av在线观看 | 看黄色毛片网站| 国产主播在线观看一区二区| 1024手机看黄色片| 国产亚洲精品av在线| 国产精品免费一区二区三区在线| 三级毛片av免费| 激情在线观看视频在线高清| 亚洲狠狠婷婷综合久久图片| 国产高清videossex| 婷婷精品国产亚洲av在线| 在线观看午夜福利视频| 欧美乱色亚洲激情| 日本精品一区二区三区蜜桃| 日本撒尿小便嘘嘘汇集6| 99视频精品全部免费 在线| 国产精品爽爽va在线观看网站| 一本精品99久久精品77| 国产高清激情床上av| 又爽又黄无遮挡网站| 国产野战对白在线观看| 国产中年淑女户外野战色| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久久久人妻精品电影| 搡老熟女国产l中国老女人| 国产一区二区亚洲精品在线观看| 久久国产精品人妻蜜桃| 免费av毛片视频| 成人亚洲精品av一区二区| 欧美+亚洲+日韩+国产| 成年人黄色毛片网站| 精品一区二区三区视频在线观看免费| 亚洲精华国产精华精| 婷婷亚洲欧美| 69av精品久久久久久| 手机成人av网站| 亚洲av熟女| 国产亚洲精品一区二区www| 舔av片在线| 国产亚洲精品av在线| 一个人免费在线观看电影| 国内精品一区二区在线观看| 色播亚洲综合网| 国产精品1区2区在线观看.| 国产精品久久视频播放| 国产黄片美女视频| 亚洲精品久久国产高清桃花| 国产精品1区2区在线观看.| 午夜福利在线观看吧| 怎么达到女性高潮| 欧美黑人欧美精品刺激| 两个人看的免费小视频| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 国产主播在线观看一区二区| 成熟少妇高潮喷水视频| 精品人妻偷拍中文字幕| 无人区码免费观看不卡| 精品电影一区二区在线| 欧美日韩福利视频一区二区| 国内少妇人妻偷人精品xxx网站| 国产真实乱freesex| 999久久久精品免费观看国产| 国产v大片淫在线免费观看| 成人三级黄色视频| 精品日产1卡2卡| 国产蜜桃级精品一区二区三区| 国产综合懂色| 天堂动漫精品| 国产精品野战在线观看| 亚洲av免费高清在线观看| 动漫黄色视频在线观看| 亚洲黑人精品在线| 欧美日韩精品网址| 午夜福利在线观看吧| 夜夜看夜夜爽夜夜摸| 国产淫片久久久久久久久 | 日本熟妇午夜| 亚洲aⅴ乱码一区二区在线播放| 老熟妇乱子伦视频在线观看| 久久亚洲真实| 在线观看免费视频日本深夜| 国产精品嫩草影院av在线观看 | 男女之事视频高清在线观看| 亚洲最大成人中文| 一个人观看的视频www高清免费观看| 国产黄片美女视频| 男人的好看免费观看在线视频| 国产高潮美女av| 最新在线观看一区二区三区| 亚洲片人在线观看| 男女那种视频在线观看| 女生性感内裤真人,穿戴方法视频| a级一级毛片免费在线观看| 乱人视频在线观看| 老司机在亚洲福利影院| 欧美乱妇无乱码| 国内精品美女久久久久久| 男人和女人高潮做爰伦理| 亚洲精品亚洲一区二区| 亚洲一区高清亚洲精品| 国产三级在线视频| 99久国产av精品| 免费av毛片视频| 国产成人a区在线观看| 国产欧美日韩精品亚洲av| 亚洲av中文字字幕乱码综合| or卡值多少钱| 少妇熟女aⅴ在线视频| 亚洲五月天丁香| 日本黄大片高清| 少妇人妻精品综合一区二区 | 精品日产1卡2卡| www日本黄色视频网| 国产成人啪精品午夜网站| 香蕉av资源在线| 高潮久久久久久久久久久不卡| 成人特级av手机在线观看| 亚洲一区二区三区不卡视频| 成人午夜高清在线视频| 一级黄色大片毛片| 亚洲人与动物交配视频| 国产野战对白在线观看| 亚洲av免费高清在线观看| 人妻夜夜爽99麻豆av| 日本黄色片子视频| 午夜a级毛片| 国产在视频线在精品| 亚洲av成人av| 欧美日韩国产亚洲二区| 男女视频在线观看网站免费| 99久久久亚洲精品蜜臀av| 老司机午夜福利在线观看视频| 91麻豆精品激情在线观看国产| 亚洲av电影不卡..在线观看| 一区福利在线观看| 国内精品久久久久精免费| 舔av片在线| 欧美不卡视频在线免费观看| www.www免费av|