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

    正辛烷和1-辛烯在CeY分子篩上吸附的分子模擬研究

    2018-11-13 08:01:38任卓然秦玉才宋麗娟
    石油煉制與化工 2018年11期
    關鍵詞:辛烷辛烯等溫線

    任卓然,李 強,秦玉才,宋麗娟

    (遼寧石油化工大學遼寧省石油化工催化科學與技術重點實驗室,遼寧 撫順 113001)

    石油的重質(zhì)化趨勢使得其深度加工技術的重要性不言而喻[1]。催化裂化是目前常用的重油輕質(zhì)化技術[2-3]。在重質(zhì)油催化裂化過程中,C8以上的烷烴和烯烴都是十分重要的催化裂化原料[4]。其在催化劑上吸附擴散行為決定了催化裂化的效果,所以對大分子烴類在催化劑上吸附行為的研究顯得尤為重要。

    前人的研究較多關注重油催化裂化過程中重油組分的變化,以及烷烴、烯烴含量對辛烷值等性能指標的影響[5-6]。為了提高這些指標,不少的研究者對催化劑本身進行了大量的研究,包括新材料的合成和既有催化劑的改性,其中Ce離子改性的Y分子篩是很好的催化裂化催化劑。大量實驗研究表明Ce離子引入Y分子篩中,Ce離子是以Ce(OH)2+形式存在的。在稀土離子交換過程中,Ce(OH)2+首先通過物理吸附作用定位在分子篩的超籠中,繼而在焙燒過程中脫去水分子由超籠進入到方鈉石籠中,適度的離子交換度可提高分子篩的催化活性。且由于稀土離子的引入,在提高重油轉(zhuǎn)化和收率方面有顯著改善。在催化裂化過程中,這些指標實際上與原料在催化裂化催化劑中的內(nèi)擴散也密切相關。對于原料大分子在分子篩中的傳質(zhì)行為與重油轉(zhuǎn)化效率之間也存在一定的關系,目前這方面的研究尚未系統(tǒng)開展,催化劑中大分子傳質(zhì)行為的研究依然亟需深入。運用傳統(tǒng)的實驗方法去研究大尺寸、復雜分子在稀土改性Y分子篩上的吸附擴散性能,不僅存在成本高、操作復雜、消耗時間更多等缺點,而且又無法觀察和確定物質(zhì)在可以到達的分子篩孔道中的吸附和擴散行為[7-8]。與其相比,利用計算機模擬的方法不但成本更低、安全且不受各種反應過程的限制,準確性和重復性更好,又可以清楚地觀察到反應物在分子篩孔道內(nèi)的運動過程,并考察此過程的熱力學性質(zhì)和動力學性質(zhì)[9]。蒙特卡洛方法是模擬計算常用的方法,可應用于研究吸附分子在分子篩中的能量與位置以及其隨吸附量、溫度、壓力等的變化。并可以提供吸附等溫線,徑向分布函數(shù)等數(shù)據(jù)[10]。

    本課題以1-辛烯、正辛烷在催化裂化過程中Ce離子改性Y分子篩上的吸附為研究主體,利用巨正則蒙特卡洛(GCMC)的方法和分子動力學(MD)[11-12]方法探究正辛烷、1-辛烯在CeY分子篩上吸附擴散的行為,為實驗提供理論支持。

    1 計算模型和計算方法

    1.1 分子篩模型的建立

    利用Accelrys公司的Material Studio 5.5軟件數(shù)據(jù)庫建立Y型分子篩骨架結構,其結構參數(shù)為:空間群為Fd-3,晶胞參數(shù)a=b=c=2.502 8 nm,α=β=γ=90°。建立2×2×2的超胞模型(ultra cell,UC),利用perl程序?qū)⒛P椭械牟糠諷i原子替換為Al原子,Al原子的分布遵循Lowenstein規(guī)則[13],將分子篩骨架中的硅鋁比調(diào)整為2.62,此時,分子篩骨架中含Al原子424個。再用Sorption模塊向體系中加入相應的Na+和Ce(OH)2+以平衡由于Si被Al取代而產(chǎn)生的負電荷,其中Ce(OH)2+的交換度為45%,由下式得到:

    (1)

    式中:q為陽離子所帶電荷數(shù);m為陽離子個數(shù);x為Ce取代前的NaY中Na+的個數(shù);X為陽離子交換度,%。圖1為NaY分子篩及Ce(OH)2+質(zhì)量分數(shù)為45%的CeY 分子篩的結構模型。由式(1)計算可知,在CeY模型中含Ce(OH)2+共96個。

    圖1 NaY分子篩及Ce(OH)2+質(zhì)量分數(shù)為45%的CeY分子篩的結構模型—H; —O; —Si; —Na; —Al; —Ce

    1.2 計算方法

    在吸附模擬前,應用DMol3模塊,采用BLYP交換相關泛函數(shù)和DNP基組對吸附質(zhì)分子進行基于量子化學的結構優(yōu)化。在巨正則蒙特卡洛模擬吸附過程中,采用COMPASS力場[14],該力場基于量子力學從頭計算,將有機與無機力場相統(tǒng)一,更能夠準確模擬有機物分子與分子篩之間的作用。利用Sorption模塊模擬吸附質(zhì)分子在分子篩中的吸附擴散行為,隨機取樣方法為Metropolis取樣,即在模擬過程中,吸附質(zhì)分子的出現(xiàn)與消失具有相同的概率,分子篩與吸附質(zhì)分子以及吸附質(zhì)分子之間的靜電相互作用的處理采用Ewald加和方法,而范德華力的相互作用利用Atom based方法,截斷半徑為1.85 nm。計算過程中的模擬步數(shù)設為1.1×107步,為了使體系平衡,前1×106步設為平衡段,剩下的1×107步設為取樣段。

    1.3 吸附等溫線擬合模型

    為了說明Ce(OH)2+對NaY分子篩吸附的影響,對使用GCMC方法得到的吸附等溫線數(shù)據(jù),分別利用Langmuir和Langmuir-Freundlich兩種吸附模型對吸附等溫線進行擬合。

    Langmuir等溫方程在吸附擴散領域應用甚廣,其模型假設吸附劑表面均勻,吸附分子間無相互作用,且單分子層吸附,即一個活性位點只能占據(jù)一個分子。Langmuir方程擬合式為:

    (2)

    式中:Q為吸附量;Qm為飽和吸附量;b為吸附平衡常數(shù),即吸附系數(shù),它與溫度有關;p為組分的分壓。Langmuir-Freundlich等溫方程[15]是在結合Langmuir方程與Freundlich方程基礎上建立的方程,考慮了吸附位分布不均勻,以及吸附質(zhì)分子間存在相互作用時表面的吸附。Langmuir-Freundlich方程擬合式為:

    (3)

    式中,n代表均相指數(shù),n越趨近于1,吸附劑表面的吸附越均勻,吸附質(zhì)分子間作用力越小且分散越均勻,越符合理想的Langmuir等溫吸附方程。

    1.4 吸附能量和吸附熱

    吸附能在吸附過程中是一個重要的參數(shù),它是由體系中所有的吸附質(zhì)和吸附劑之間以及吸附質(zhì)與吸附質(zhì)之間的能量決定的。

    Ead=Eads-ads+Eads-zeo

    (4)

    式中:Ead為體系總能量;Eads-ads為吸附質(zhì)和吸附質(zhì)之間的相互作用能量;Eads-zeo為吸附質(zhì)和分子篩之間的相互作用能量。

    在吸附過程中分子放出熱量,從而達到相對穩(wěn)定狀態(tài),吸附熱計算式如下:

    (5)

    式中:Q1為吸附熱;R=8.314 J(molK),為氣體常數(shù);Ead為Eads-ads與Eads-Zeo之和;Eintra為吸附質(zhì)分子勢能;Nad為吸附質(zhì)個數(shù)[11]。

    1.5 吸附能量密度等高截面圖

    為了能夠直觀地觀察到吸附質(zhì)在CeY中的吸附情況,選取十二元環(huán)中心來進行截面[12],圖2a中藍線為截取位置,先對模型構建Connoly 表面,然后在對其進行截面。圖2b為十二元環(huán)結構圖,圖2c為圖2a左視圖,圖2d為吸附能量密度俯視圖,可以清晰地看到十二元環(huán)以及其它原子分布情況。

    圖2 CeY骨架結構及分子吸附能量密度分布—H; —O; —Si; —Na; —Al; —Ce

    2 結果與討論

    2.1 吸附等溫線

    為了探討正辛烷和1-辛烯在Ce離子改性Y型分子篩上吸附擴散行為,應用GCMC方法分別在303 K、423 K溫度下對其進行模擬研究。模擬壓力從1.0×10-14kPa到其飽和蒸氣壓,正辛烷和1-辛烯分子在不同溫度下的飽和蒸氣壓如表1所示。正辛烷和1-辛烯在CeY分子篩上的吸附等溫線如圖3所示,圖中散點為原始模擬數(shù)據(jù),線為對原始模擬數(shù)據(jù)進行的Langmuir-Freundlich模型擬合數(shù)據(jù)。由圖3可知:在溫度為303 K時,當壓力增加到1×10-7kPa左右,正辛烷和1-辛烯分子開始在CeY上吸附,隨后吸附量迅速上升,并在3.2×10-3kPa時吸附逐漸趨于平緩并最終達到吸附飽和;而在423 K時,正辛烷、1-辛烯分子在壓力達到2.49×10-3kPa時才開始進行吸附,隨后吸附量亦快速增加,并在4.48 kPa壓力下逐漸趨于吸附飽和。

    表1 正辛烷、1-辛烯在不同溫度下的飽和蒸氣壓 kPa

    圖3 不同溫度下正辛烷和1-辛烯在CeY分子篩上的吸附等溫線■—303 K-CeY-正辛烷; ▲—303 K-CeY-1-辛烯; ●—423 K-CeY-正辛烷; 辛烯

    在303 K溫度下,正辛烷和1-辛烯的飽和吸附量都大于423 K下的飽和吸附量。這是因為在低溫下,分子動能較低,CeY僅需一個較小的作用力就可以牢牢地將吸附質(zhì)吸附在活性位點上;而423 K下,分子熱運動加劇,分子動能較大,吸附質(zhì)與分子篩之間的范德華力被削弱,從而使得CeY分子篩無法捕捉到較多的氣體分子,導致飽和吸附量的下降;303 K、423 K溫度下,正辛烷的飽和吸附量都大于1-辛烯的飽和吸附量。

    為了進一步研究正辛烷和1-辛烯的吸附行為,分別用Langmuir以及Langmuir-Freundlich吸附模型對模擬出的吸附等溫線進行擬合,擬合數(shù)據(jù)如表2所示。

    根據(jù)表2中相關系數(shù)(R)的比較可以看出:兩種擬合方法的相關性都很高,接近于1;但Langmuir-Freundlich模型對吸附等溫線擬合的相關程度比利用Langmuir模型擬合的效果更好,說明正辛烷、1-辛烯的吸附行為,并非單分子層吸附,而且分子篩上吸附位分布并不均一,分子和分子之間存在相互作用力;303 K下,Langmuir-Freundlich模型擬合均相指數(shù)(n)較小,1-辛烯吸附的n值小于正辛烷的n值;這是由于Ce(OH)2+的引入造成分子篩吸附位分布不均一,且1-辛烯分子之間的相互作用力更強而導致的;Langmuir-Freundlich模型的飽和吸附量更接近原始數(shù)據(jù),相關性更好;隨著溫度升高,423 K下Langmuir模型R值和Langmuir-Freundlich模型R值都更接近1,Langmuir-Freundlich模型中n值變大,1-辛烯的n值略小于正辛烷的n值,但相比低溫二者n值更加接近。這說明溫度升高使得氣體分子動能增大,在體系中分散度更高,1-辛烯分子間相互作用也減弱。

    表2 不同溫度下CeY與正辛烷、1-辛烯的Langmuir和Langmuir-Freundlich模型擬合參數(shù)

    2.2 吸附質(zhì)與吸附劑之間的相互作用

    2.2.1吸附勢能曲線為了進一步探討正辛烷、1-辛烯與CeY之間的相互作用,從吸附能量角度來對其進行考察。圖4為不同溫度、不同吸附質(zhì)的吸附勢能變化曲線。圖5是303 K、423 K溫度下兩種吸附質(zhì)在吸附過程中吸附質(zhì)的分布密度。

    圖4 不同溫度下CeY分子篩與正辛烷和1-辛烯之間的吸附勢能曲線

    圖4a為303 K下不同壓力時正辛烷在CeY上吸附勢能曲線。由圖4a可以看出:在1.56×10-7kPa下,有少量正辛烷分子開始吸附;雖然曲線上看似只有一個峰,但是經(jīng)過擬合后發(fā)現(xiàn),在-87.94 kJmol處的主峰右側-82.5 kJmol的位置上有一個伴峰,在壓力小于3.25×10-3kPa的情況下都可以看出這2個峰。這說明在開始吸附到壓力較小時,正辛烷在CeY上存在兩種位置的吸附。由圖5a可以看出:在低壓吸附時,正辛烷在Ce(OH)2+附近存在吸附密度,當壓力逐漸增加后,正辛烷的吸附逐漸轉(zhuǎn)向Na+附近,Ce(OH)2+附近正辛烷出現(xiàn)的概率大大減弱;當吸附壓力高于3.25×10-3kPa時,吸附逐漸趨向飽和,勢能曲線的變化也趨向穩(wěn)定,最后在吸附壓力達到2.46 kPa時,正辛烷的吸附量達到飽和。此時,由圖5a可以看出,正辛烷只和分子篩骨架及Na+發(fā)生作用,吸附勢能曲線中只有一個勢能峰。但是當溫度升高至423 K時(如圖4c),無論壓力高低,正辛烷在CeY上吸附勢能均只有一種,同時,隨著吸附壓力的升高,峰位出現(xiàn)了明顯的位移,最可幾相互作用隨著吸附壓力的增大而變小,這與低溫時的變化情況相似。結合圖5b可以看出,高溫下吸附密度出現(xiàn)的區(qū)域都是在Na+附近,這說明正辛烷在Ce(OH)2+附近的吸附雖然比在Na+附近的吸附穩(wěn)定,但是受吸附溫度影響較大,溫度升高后,正辛烷在Ce(OH)2+附近的吸附能力變低。

    圖5 不同溫度下CeY分子篩與正辛烷和1-辛烯之間的吸附勢能密度截面分布紅藍色塊代表正辛烷或1-辛烯的分布概率,紅色為出現(xiàn)概率高,藍色為概率低

    對于1-辛烯分子的吸附,其情況與正辛烷稍有不同。結合圖4b和圖5c可以看出:在低溫下,吸附勢能曲線存在2個峰,即-87.78 kJmol、-96.14 kJmol,即兩種相互作用力,1-辛烯的吸附位同時出現(xiàn)在Na+和Ce(OH)2+附近;隨著吸附壓力的增加,-87.78 kJmol位置的勢能峰依舊存在,但是-96.14 kJmol位置的勢能峰則隨之迅速降低,直到4.89×10-6kPa飽和吸附壓力下幾乎消失。說明在壓力增大的過程中,隨著1-辛烯吸附量的上升,有一種相對較弱的相互作用力逐漸減小甚至消失;當溫度較高時,1-辛烯在CeY上的吸附也只有一種吸附勢能,且隨著吸附量的增加沒有較大的變化;與正辛烷吸附進行對比后發(fā)現(xiàn),1-辛烯的吸附勢能變化遠小于正辛烷。

    由以上分析可以看出,303 K下1-辛烯與正辛烷的吸附類似。均是在低的吸附量時才能穩(wěn)定地出現(xiàn)吸附質(zhì)分子與Ce(OH)2+的作用。隨著吸附量的增加,這種吸附作用逐漸減小。這也可以認為,吸附質(zhì)分子與CeY作用時,由于Ce(OH)2+的數(shù)量是固定的,無論吸附量增加多少,吸附在Ce(OH)2+的分子數(shù)量是不會發(fā)生變化的,所以大多數(shù)的分子都是在Na+附近和骨架旁邊吸附。從吸附勢能曲線的移動可以看出,正辛烷與CeY之間的相互作用隨著吸附壓力增加和吸附量的增加影響較大。反之,1-辛烯隨著吸附壓力增大與吸附量的增加對其相互作用影響較小。由此推斷,這可能是由于1-辛烯上的C=C鍵產(chǎn)生的π-π相互作用,導致1-辛烯分子與CeY之間存在比正辛烷分子與CeY間更強的相互作用而導致的。正辛烷、1-辛烯與CeY之間的主要相互作用在孔道中的Na+附近。

    根據(jù)飽和吸附壓力下的勢能分布曲線,整理出不同溫度下正辛烷、1-辛烯與CeY分子篩之間最可幾相互作用能量,如表3所示。

    表3 正辛烷、1-辛烯在CeY分子篩中飽和吸附時的最可幾相互作用能 kJmol

    表3 正辛烷、1-辛烯在CeY分子篩中飽和吸附時的最可幾相互作用能 kJmol

    溫度/KCeY-正辛烷CeY-1-辛烯303-88.67-86.87423-84.18-83.54

    由表3也可以看出,423 K下正辛烷、1-辛烯的最可幾相互作用能都大于303 K下的最可幾相互作用能。這是由于溫度的升高,分子動能增加,從而提高了體系的能量,使吸附態(tài)趨向較不穩(wěn)定的情況,從而降低了CeY與吸附質(zhì)之間的結合能,導致飽和吸附量的降低。且不論是高溫、低溫下,1-辛烯的最可幾相互作用能也都大于正辛烷,所以在不同溫度下,正辛烷的飽和吸附量要大于1-辛烯,也印證了吸附等溫線的結果。

    為了直觀地觀察正辛烷、1-辛烯在CeY上的擴散情況,探究了在不同溫度下飽和吸附情況下的單元晶胞上吸附質(zhì)的勢能密度分布情況,如圖6所示。

    圖6 不同溫度下正辛烷與1-辛烯在CeY分子篩上的勢能密度分布藍色越深說明分子勢能越高,紅色越深則說明分子勢能越低

    從圖6可以看出:在303 K下(圖6a、c)正辛烷和1-辛烯的分布范圍較為集中;在423 K下(圖6b、d)吸附質(zhì)分子的分布范圍更大。這是由于溫度的升高,分子動能增大,從而增大了分布范圍。從不同吸附質(zhì)來看,不論是低溫還是高溫,1-辛烯的飽和吸附量都明顯低于正辛烷,且1-辛烯的分布更為均勻。對比圖6a、b和圖6c、d可以看出,正辛烷的分布范圍隨溫度的變化相比較于1-辛烯的變化更大,分布更廣。結合表3的數(shù)據(jù)也可以看出,正辛烷的最可幾相互作用能隨溫度變化較大。

    2.2.2吸附熱圖7為模擬得到的兩種吸附質(zhì)分子在CeY中吸附的吸附熱隨吸附量的變化曲線。由圖7可以看出:正辛烷與1-辛烯在423 K下的吸附熱都高于303 K下的吸附熱,因為吸附過程是放熱過程,在高溫下分子動能更大,想要吸附在分子篩上則需要放出更多的能量來降低自身勢能從而達到一個相對穩(wěn)定的吸附態(tài);兩個溫度下,正辛烷的吸附熱都大于1-辛烯,這是由于1-辛烯內(nèi)部存在C=C鍵,分子本身的勢能高于正辛烷。

    圖7 CeY分子篩與正辛烷和1-辛烯之間的吸附熱■—303 K-CeY-正辛烷; 辛烯;●—423 K-CeY-正辛烷; ▲—423 K-CeY-1-辛烯

    從圖7還可以看出:正辛烷、1-辛烯的吸附熱,在吸附量較小的情況下變化較為緩慢,在吸附量增大時,吸附熱降低的速率也更快。這說明吸附量的大小對于吸附熱有一定的影響;隨著吸附量的增大,分子和分子之間相互作用也越來越明顯,分子與Ce(OH)2+之間的相互作用會減弱,新來的氣體分子只要釋放相對較少的能量就能吸附在CeY的活性位上;在303 K、423 K下的起始吸附階段,吸附量較小時,正是由于1-辛烯的吸附熱小于正辛烷,所以在起始吸附階段,1-辛烯先于正辛烷吸附在CeY上;而在吸附接近飽和時,正辛烷的吸附熱反而略低于1-辛烯,這也從吸附熱角度解釋了為什么正辛烷在不同溫度下的飽和吸附量會高于1-辛烯;且在吸附起始到飽和正辛烷的吸附熱變化幅度大于1-辛烯,也印證了圖3中正辛烷能量曲線變化幅度大于1-辛烯的結果。

    3 結 論

    (1)正辛烷、1-辛烯與CeY分子篩之間的吸附符合Langmuir-Freundlich吸附模型。Ce(OH)2+的引入使得分子篩上吸附活性位分布不均勻,且1-辛烯分子間由于C=C鍵的存在使其分子間相互作用力更強,因此,在303 K和423 K下,1-辛烯的飽和吸附量均低于正辛烷。

    (2)正辛烷、1-辛烯在CeY分子篩上的吸附量對吸附熱存在影響,隨著吸附量的增加,這種影響也增大。吸附作用中最主要的相互作用力為范德華力,吸附質(zhì)分子主要與孔道中Na+作用,與Ce(OH)2+之間的相互作用次之,且隨著吸附量的增加,分子與Ce(OH)2+之間的相互作用會減弱。

    猜你喜歡
    辛烷辛烯等溫線
    溶液法聚乙烯裝置1-辛烯的效能提升應用
    全氟辛烷磺酸對斑馬魚胚胎的急性毒性與致畸效應
    河蜆(Corbicula fluminea)對氧化石墨烯和全氟辛烷磺酸類物質(zhì)(PFOS)聯(lián)合脅迫的生理生化響應
    低溫吸附劑的低溫低壓吸附等溫線分布研究
    隨機環(huán)辛烷鏈的三類 Kirchhoff指數(shù)
    基于氟—氟相互作用的上轉(zhuǎn)換熒光法快速測定水中的全氟辛烷磺酸
    分析化學(2019年3期)2019-03-30 10:59:24
    辛烯共聚聚乙烯樹脂非等溫結晶動力學
    中國塑料(2015年6期)2015-11-13 03:02:48
    反式聚辛烯橡膠/丁腈橡膠并用膠的性能研究
    如何在新課改背景下突破等溫線判讀中的難點
    辛烯基琥珀酸魔芋葡甘聚糖酯的微波半干法制備工藝及性能
    食品科學(2013年14期)2013-03-11 18:25:01
    高清视频免费观看一区二区| 久热爱精品视频在线9| 免费日韩欧美在线观看| 午夜免费成人在线视频| 757午夜福利合集在线观看| 色综合婷婷激情| 精品福利观看| 男女高潮啪啪啪动态图| 久久人妻av系列| 亚洲精品中文字幕在线视频| 久久久久久久大尺度免费视频| 国产无遮挡羞羞视频在线观看| 久久精品国产a三级三级三级| 巨乳人妻的诱惑在线观看| 亚洲,欧美精品.| 精品国产乱码久久久久久男人| 欧美成狂野欧美在线观看| 菩萨蛮人人尽说江南好唐韦庄| 免费av中文字幕在线| 大片免费播放器 马上看| 中文亚洲av片在线观看爽 | 欧美国产精品va在线观看不卡| 亚洲熟女毛片儿| 嫩草影视91久久| av视频免费观看在线观看| 久久国产精品大桥未久av| 大香蕉久久成人网| 亚洲欧洲精品一区二区精品久久久| 色播在线永久视频| 亚洲精品在线美女| www.999成人在线观看| 久久免费观看电影| 国产精品久久电影中文字幕 | 亚洲中文日韩欧美视频| 中文字幕制服av| 国产一卡二卡三卡精品| 午夜福利在线观看吧| 久久久精品国产亚洲av高清涩受| 午夜免费鲁丝| 婷婷成人精品国产| 精品少妇黑人巨大在线播放| 一级片免费观看大全| 色婷婷av一区二区三区视频| 国精品久久久久久国模美| 精品一区二区三区av网在线观看 | 亚洲 欧美一区二区三区| 精品免费久久久久久久清纯 | 午夜视频精品福利| 蜜桃在线观看..| 久久人妻福利社区极品人妻图片| 菩萨蛮人人尽说江南好唐韦庄| 黄片播放在线免费| 国产精品一区二区精品视频观看| 欧美在线黄色| 久久精品熟女亚洲av麻豆精品| 动漫黄色视频在线观看| 国产精品久久久久久精品电影小说| 久9热在线精品视频| 久久免费观看电影| 黄色视频,在线免费观看| 老司机影院毛片| 99国产综合亚洲精品| 亚洲精品久久成人aⅴ小说| 日韩视频在线欧美| 性高湖久久久久久久久免费观看| 亚洲精品国产精品久久久不卡| 欧美日韩亚洲综合一区二区三区_| 十八禁网站网址无遮挡| 国产欧美亚洲国产| 两个人免费观看高清视频| 精品一区二区三区四区五区乱码| 精品欧美一区二区三区在线| 亚洲国产av影院在线观看| 在线av久久热| 亚洲人成电影观看| 国产一区二区三区在线臀色熟女 | 一区二区日韩欧美中文字幕| cao死你这个sao货| 一级片'在线观看视频| netflix在线观看网站| 黄色 视频免费看| 深夜精品福利| 日韩有码中文字幕| 成人国产av品久久久| 久久热在线av| 国产成人av激情在线播放| 99香蕉大伊视频| 老司机午夜福利在线观看视频 | 免费黄频网站在线观看国产| 满18在线观看网站| 久久久久精品国产欧美久久久| 黄色视频在线播放观看不卡| 欧美激情 高清一区二区三区| 超色免费av| 黄色丝袜av网址大全| 老司机亚洲免费影院| 国产成人精品久久二区二区免费| 色在线成人网| 新久久久久国产一级毛片| 午夜福利视频在线观看免费| 婷婷成人精品国产| 免费少妇av软件| 肉色欧美久久久久久久蜜桃| 亚洲欧洲日产国产| 亚洲熟女毛片儿| 人人妻人人爽人人添夜夜欢视频| 日韩有码中文字幕| 欧美国产精品va在线观看不卡| 美女午夜性视频免费| 国产福利在线免费观看视频| 久久人妻av系列| 精品一区二区三区视频在线观看免费 | 91成人精品电影| 久久午夜综合久久蜜桃| 久久午夜亚洲精品久久| 国产精品av久久久久免费| 黑人巨大精品欧美一区二区蜜桃| a级毛片在线看网站| 天天影视国产精品| 精品福利观看| 怎么达到女性高潮| 国产熟女午夜一区二区三区| 欧美av亚洲av综合av国产av| 日韩免费av在线播放| tube8黄色片| 国产伦人伦偷精品视频| 欧美精品av麻豆av| 狠狠精品人妻久久久久久综合| 日韩熟女老妇一区二区性免费视频| 欧美日韩视频精品一区| 一本色道久久久久久精品综合| 国产有黄有色有爽视频| 一区二区三区激情视频| 热re99久久国产66热| 少妇的丰满在线观看| 久久久久久久久久久久大奶| 国产在线视频一区二区| 日韩三级视频一区二区三区| 中亚洲国语对白在线视频| 十八禁人妻一区二区| 一夜夜www| 99精品在免费线老司机午夜| 18禁美女被吸乳视频| 日韩成人在线观看一区二区三区| 国产成人免费观看mmmm| av超薄肉色丝袜交足视频| videosex国产| 9热在线视频观看99| 18禁裸乳无遮挡动漫免费视频| 男女床上黄色一级片免费看| 亚洲av成人不卡在线观看播放网| 亚洲专区国产一区二区| 一区二区三区乱码不卡18| 亚洲熟女毛片儿| 十八禁高潮呻吟视频| 国产一卡二卡三卡精品| 亚洲精品美女久久久久99蜜臀| 久久av网站| 老汉色av国产亚洲站长工具| 国产视频一区二区在线看| e午夜精品久久久久久久| 亚洲国产精品一区二区三区在线| 国产精品九九99| 国产精品国产av在线观看| 岛国在线观看网站| 日韩制服丝袜自拍偷拍| 人人妻人人澡人人爽人人夜夜| 99精品久久久久人妻精品| 青青草视频在线视频观看| 亚洲精品一卡2卡三卡4卡5卡| 欧美精品高潮呻吟av久久| 女人久久www免费人成看片| 中文欧美无线码| 99香蕉大伊视频| bbb黄色大片| 这个男人来自地球电影免费观看| 丁香六月天网| 国产亚洲一区二区精品| 黄色视频不卡| 亚洲欧美精品综合一区二区三区| 99精品在免费线老司机午夜| 成人三级做爰电影| 久久午夜亚洲精品久久| 波多野结衣一区麻豆| 日韩人妻精品一区2区三区| 99国产精品一区二区蜜桃av | 久久久水蜜桃国产精品网| 国产一区二区三区在线臀色熟女 | 亚洲三区欧美一区| 十分钟在线观看高清视频www| 夜夜爽天天搞| 大陆偷拍与自拍| 日韩人妻精品一区2区三区| 精品久久蜜臀av无| 老司机在亚洲福利影院| 热99久久久久精品小说推荐| 久久久欧美国产精品| 超碰97精品在线观看| 午夜久久久在线观看| 欧美黄色淫秽网站| 午夜激情久久久久久久| 成人精品一区二区免费| av天堂在线播放| 91成年电影在线观看| 欧美日韩一级在线毛片| 午夜福利一区二区在线看| 午夜福利一区二区在线看| 成在线人永久免费视频| 人人澡人人妻人| 久久久久久人人人人人| www.999成人在线观看| 亚洲熟女精品中文字幕| 黄色怎么调成土黄色| 亚洲三区欧美一区| 欧美日韩中文字幕国产精品一区二区三区 | 欧美激情久久久久久爽电影 | 成年人黄色毛片网站| 色尼玛亚洲综合影院| 亚洲自偷自拍图片 自拍| 狠狠精品人妻久久久久久综合| 国产三级黄色录像| av一本久久久久| 大型黄色视频在线免费观看| 母亲3免费完整高清在线观看| 久久久国产一区二区| 老司机亚洲免费影院| 男女无遮挡免费网站观看| 老司机影院毛片| 操美女的视频在线观看| 国产在视频线精品| 不卡一级毛片| 久久精品亚洲精品国产色婷小说| 国产成人免费观看mmmm| 自拍欧美九色日韩亚洲蝌蚪91| 黑人猛操日本美女一级片| 亚洲欧洲精品一区二区精品久久久| 免费在线观看日本一区| 国产99久久九九免费精品| 美女扒开内裤让男人捅视频| 久久国产亚洲av麻豆专区| 99九九在线精品视频| 高清欧美精品videossex| 一边摸一边抽搐一进一出视频| 伊人久久大香线蕉亚洲五| 热99久久久久精品小说推荐| 人人妻人人爽人人添夜夜欢视频| 丝袜喷水一区| kizo精华| 亚洲精品国产精品久久久不卡| 午夜福利免费观看在线| 捣出白浆h1v1| 老汉色av国产亚洲站长工具| 国产伦人伦偷精品视频| 高清黄色对白视频在线免费看| 久久久久久久久久久久大奶| 精品国产亚洲在线| av网站免费在线观看视频| 亚洲色图 男人天堂 中文字幕| 操美女的视频在线观看| 成人国语在线视频| 久久精品国产综合久久久| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久成人av| 老汉色∧v一级毛片| 久久天堂一区二区三区四区| 免费高清在线观看日韩| 欧美日韩福利视频一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 精品久久久精品久久久| 飞空精品影院首页| 久久中文看片网| 母亲3免费完整高清在线观看| 在线观看免费午夜福利视频| 99香蕉大伊视频| 欧美精品av麻豆av| 69精品国产乱码久久久| 亚洲,欧美精品.| 国产精品一区二区在线观看99| 视频在线观看一区二区三区| 精品第一国产精品| 每晚都被弄得嗷嗷叫到高潮| 国产欧美亚洲国产| 国产成人免费观看mmmm| 777久久人妻少妇嫩草av网站| 宅男免费午夜| 女人精品久久久久毛片| 国产精品99久久99久久久不卡| 久久精品国产亚洲av香蕉五月 | 老汉色av国产亚洲站长工具| e午夜精品久久久久久久| 757午夜福利合集在线观看| 他把我摸到了高潮在线观看 | 国产精品久久久久久精品电影小说| 国产精品.久久久| 狂野欧美激情性xxxx| 成年版毛片免费区| 美女高潮喷水抽搐中文字幕| 国产精品美女特级片免费视频播放器 | 99香蕉大伊视频| 免费不卡黄色视频| 欧美黄色片欧美黄色片| 亚洲精品国产精品久久久不卡| 久久精品国产亚洲av高清一级| 日韩有码中文字幕| 国产精品一区二区免费欧美| 亚洲自偷自拍图片 自拍| 美女视频免费永久观看网站| 国产日韩欧美在线精品| 一二三四社区在线视频社区8| 久久中文字幕人妻熟女| 国产色视频综合| 亚洲第一欧美日韩一区二区三区 | 91成人精品电影| 亚洲,欧美精品.| 久久亚洲真实| 国产精品欧美亚洲77777| 中文亚洲av片在线观看爽 | 少妇粗大呻吟视频| 日韩欧美免费精品| 久久久久精品国产欧美久久久| 热re99久久精品国产66热6| 人人妻人人澡人人爽人人夜夜| 欧美日韩一级在线毛片| 欧美久久黑人一区二区| 成人18禁高潮啪啪吃奶动态图| 婷婷丁香在线五月| 国产av又大| 成人手机av| 日韩熟女老妇一区二区性免费视频| 老熟女久久久| 亚洲全国av大片| 在线观看www视频免费| 一区二区三区乱码不卡18| 丰满少妇做爰视频| av片东京热男人的天堂| 色精品久久人妻99蜜桃| 91av网站免费观看| 日韩熟女老妇一区二区性免费视频| 国产91精品成人一区二区三区 | 国产不卡av网站在线观看| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 狂野欧美激情性xxxx| 在线观看www视频免费| 日韩欧美免费精品| xxxhd国产人妻xxx| 一级a爱视频在线免费观看| 精品亚洲成a人片在线观看| 欧美变态另类bdsm刘玥| 亚洲人成77777在线视频| 别揉我奶头~嗯~啊~动态视频| 国产精品熟女久久久久浪| 天堂动漫精品| 国产精品 欧美亚洲| 久久人人97超碰香蕉20202| 久久久久久久久久久久大奶| 少妇猛男粗大的猛烈进出视频| 中文字幕精品免费在线观看视频| 国产野战对白在线观看| 俄罗斯特黄特色一大片| 91国产中文字幕| 高清在线国产一区| 大片免费播放器 马上看| 手机成人av网站| 无遮挡黄片免费观看| 美女午夜性视频免费| 国产av精品麻豆| 亚洲九九香蕉| 国产精品偷伦视频观看了| 99热网站在线观看| 国产亚洲精品一区二区www | 中文字幕人妻熟女乱码| 菩萨蛮人人尽说江南好唐韦庄| 50天的宝宝边吃奶边哭怎么回事| av国产精品久久久久影院| 久久精品亚洲熟妇少妇任你| 国产欧美日韩一区二区三区在线| a级片在线免费高清观看视频| 国产亚洲精品久久久久5区| 1024香蕉在线观看| 欧美精品一区二区大全| 男女之事视频高清在线观看| 国产在线视频一区二区| 久久中文字幕人妻熟女| 亚洲精华国产精华精| 欧美精品一区二区免费开放| 午夜成年电影在线免费观看| 国产一区二区 视频在线| 亚洲三区欧美一区| 国产精品影院久久| 国产淫语在线视频| 精品亚洲成国产av| 在线 av 中文字幕| 欧美 亚洲 国产 日韩一| 丝袜美足系列| 亚洲午夜精品一区,二区,三区| 女性被躁到高潮视频| 精品亚洲成国产av| 精品国产乱码久久久久久男人| 久久久久久亚洲精品国产蜜桃av| 亚洲精品中文字幕在线视频| 中国美女看黄片| 国产视频一区二区在线看| 日韩欧美免费精品| 亚洲色图综合在线观看| 国产又爽黄色视频| 亚洲av欧美aⅴ国产| 久久影院123| 妹子高潮喷水视频| 考比视频在线观看| 19禁男女啪啪无遮挡网站| 十八禁人妻一区二区| 操美女的视频在线观看| 中文字幕av电影在线播放| 亚洲精品中文字幕在线视频| 一区二区日韩欧美中文字幕| 亚洲av欧美aⅴ国产| 女人精品久久久久毛片| 丰满少妇做爰视频| www.999成人在线观看| 免费久久久久久久精品成人欧美视频| 国产野战对白在线观看| 欧美日韩福利视频一区二区| 亚洲熟女毛片儿| 大型av网站在线播放| 国产成人精品久久二区二区免费| 日韩三级视频一区二区三区| 国产97色在线日韩免费| 曰老女人黄片| 亚洲av第一区精品v没综合| av一本久久久久| 12—13女人毛片做爰片一| 天堂动漫精品| 大香蕉久久成人网| 美女扒开内裤让男人捅视频| 美女国产高潮福利片在线看| 国产男靠女视频免费网站| 精品少妇黑人巨大在线播放| 女警被强在线播放| 美女福利国产在线| 极品少妇高潮喷水抽搐| 黑丝袜美女国产一区| 麻豆av在线久日| 久久中文字幕一级| 精品少妇久久久久久888优播| av在线播放免费不卡| 亚洲国产av新网站| 狠狠精品人妻久久久久久综合| 国产色视频综合| 亚洲欧美精品综合一区二区三区| 大型av网站在线播放| 国产精品一区二区免费欧美| 1024香蕉在线观看| 欧美精品一区二区免费开放| av网站在线播放免费| av片东京热男人的天堂| 又黄又粗又硬又大视频| 国产真人三级小视频在线观看| 欧美在线一区亚洲| 精品第一国产精品| 欧美 日韩 精品 国产| 欧美日本中文国产一区发布| 99香蕉大伊视频| 午夜91福利影院| 欧美黑人精品巨大| 中国美女看黄片| 黄色怎么调成土黄色| 一个人免费看片子| 18禁美女被吸乳视频| 熟女少妇亚洲综合色aaa.| 国精品久久久久久国模美| 久久婷婷成人综合色麻豆| 午夜福利在线免费观看网站| 亚洲avbb在线观看| 麻豆乱淫一区二区| 他把我摸到了高潮在线观看 | 一个人免费在线观看的高清视频| 高清av免费在线| 男女无遮挡免费网站观看| 亚洲国产欧美在线一区| 免费少妇av软件| 色综合欧美亚洲国产小说| 国产日韩欧美亚洲二区| 首页视频小说图片口味搜索| 国产一区有黄有色的免费视频| 操美女的视频在线观看| 涩涩av久久男人的天堂| 亚洲伊人久久精品综合| 涩涩av久久男人的天堂| 国产精品秋霞免费鲁丝片| 大型黄色视频在线免费观看| 日本av免费视频播放| av免费在线观看网站| 国产在线视频一区二区| 久久九九热精品免费| 亚洲精品国产一区二区精华液| 亚洲欧美一区二区三区久久| 日韩制服丝袜自拍偷拍| 国产伦理片在线播放av一区| 国产精品 国内视频| 777米奇影视久久| 精品高清国产在线一区| 国产精品久久久久久精品古装| videos熟女内射| 久久久久久久久久久久大奶| 国产高清videossex| 伦理电影免费视频| 黄色丝袜av网址大全| 天天躁夜夜躁狠狠躁躁| 免费高清在线观看日韩| 少妇的丰满在线观看| 黑人巨大精品欧美一区二区mp4| avwww免费| a在线观看视频网站| 一级毛片精品| 女性生殖器流出的白浆| 精品国内亚洲2022精品成人 | 97在线人人人人妻| 成人国产av品久久久| 悠悠久久av| 国产黄频视频在线观看| 亚洲成人免费电影在线观看| 一区二区三区精品91| 国产精品98久久久久久宅男小说| 亚洲精品粉嫩美女一区| 怎么达到女性高潮| 久久狼人影院| 日韩大片免费观看网站| 国产真人三级小视频在线观看| 久久精品国产99精品国产亚洲性色 | 日韩欧美三级三区| 亚洲成av片中文字幕在线观看| 狠狠狠狠99中文字幕| 窝窝影院91人妻| 80岁老熟妇乱子伦牲交| 99国产精品一区二区蜜桃av | 一区二区日韩欧美中文字幕| 亚洲七黄色美女视频| 黄片播放在线免费| 国产一区有黄有色的免费视频| 如日韩欧美国产精品一区二区三区| 欧美激情久久久久久爽电影 | 欧美乱妇无乱码| 丁香六月欧美| 蜜桃在线观看..| 在线av久久热| 在线观看www视频免费| 深夜精品福利| 久久精品亚洲av国产电影网| 成人精品一区二区免费| 国产91精品成人一区二区三区 | 久久久精品94久久精品| 飞空精品影院首页| 久久人妻福利社区极品人妻图片| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区蜜桃| 国产欧美日韩综合在线一区二区| 999久久久国产精品视频| 三级毛片av免费| 久久精品国产a三级三级三级| e午夜精品久久久久久久| 日韩欧美国产一区二区入口| av福利片在线| 免费在线观看影片大全网站| 国产黄频视频在线观看| 国产精品98久久久久久宅男小说| 最新的欧美精品一区二区| 在线观看www视频免费| 国产黄色免费在线视频| 欧美日韩精品网址| 人人妻人人澡人人爽人人夜夜| 午夜91福利影院| 久久中文字幕一级| 成年人午夜在线观看视频| 色精品久久人妻99蜜桃| 在线观看66精品国产| 久久人妻熟女aⅴ| 国产亚洲精品一区二区www | 91麻豆精品激情在线观看国产 | 又大又爽又粗| 亚洲成a人片在线一区二区| 999精品在线视频| 精品国产乱码久久久久久男人| av国产精品久久久久影院| 天天影视国产精品| 成人手机av| 成年人午夜在线观看视频| 在线天堂中文资源库| 丰满人妻熟妇乱又伦精品不卡| 在线观看www视频免费| 久久久久久久久久久久大奶| 午夜91福利影院| 女性生殖器流出的白浆| 一本久久精品| 亚洲av日韩在线播放| 一级黄色大片毛片| 美女主播在线视频| 免费看a级黄色片| 大片电影免费在线观看免费| 蜜桃在线观看..| 99精品在免费线老司机午夜| 99精国产麻豆久久婷婷| 亚洲自偷自拍图片 自拍| 啪啪无遮挡十八禁网站| 久热爱精品视频在线9| 91国产中文字幕| 自线自在国产av| 亚洲成国产人片在线观看| 国产免费av片在线观看野外av| 搡老熟女国产l中国老女人| 免费女性裸体啪啪无遮挡网站| 美国免费a级毛片| 十八禁网站网址无遮挡| 成人三级做爰电影| 久久天堂一区二区三区四区|