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

蒙脫石對油頁巖干酪根熱解特性的影響

2016-03-29 07:10:28王擎張宏喜遲銘書崔達許祥成
化工進展 2016年3期

王擎,張宏喜,遲銘書,崔達,許祥成

(東北電力大學能源與動力工程學院,吉林 吉林 132012)

?

蒙脫石對油頁巖干酪根熱解特性的影響

王擎,張宏喜,遲銘書,崔達,許祥成

(東北電力大學能源與動力工程學院,吉林 吉林 132012)

摘要:選取中國樺甸、撫順、窯街3個地區油頁巖,酸洗得到各干酪根樣品。應用TG-FTIR技術對不同升溫速率下各干酪根與蒙脫石共熱解行為進行了研究,并利用Coats-Redfem積分法對升溫速率為10℃/min下的干酪根摻混樣品進行熱解反應動力學分析,獲得了蒙脫石對干酪根熱解產物析出規律的影響及熱解過程的活化能(E)和指前因子(A)。結果表明,隨著蒙脫石摻混比例增大,樺甸和撫順干酪根熱解失重率先升高后降低;在熱解過程中,440℃前干酪根與蒙脫石摻混樣品的熱解失重率較干酪根單獨熱解低,而440℃后其熱解失重率較干酪根熱解有所提高;蒙脫石摻混比例增大使得各干酪根熱解產物中CH3/CH2值有所增加;有蒙脫石摻混的干酪根熱解活化能相比干酪根熱解活化能有所降低。表明蒙脫石對油頁巖干酪根熱解具有物理吸附和裂解催化兩種作用,且隨著熱解過程的加深和蒙脫石配比的改變而不同。最后,對升溫速率為10℃/min的樣品進行了熱解反應動力學分析,進一步探討了蒙脫石對干酪根的熱解影響機理。

關鍵詞:蒙脫石;干酪根;熱解;油頁巖;動力學;催化裂解

油頁巖是一種具有無機礦物質骨架并含固體有機質的沉積巖。油頁巖經干餾生成頁巖油,對頁巖油加工可制取油品,因此油頁巖是石油的一種補充能源[1]。油頁巖有機質由干酪根和瀝青質組成,而干酪根是熱解生油的主要物質。油頁巖無機物主要有石英、高嶺土、黏土、云母、碳酸鹽以及硫鐵礦等礦物質以及銅、鎳、鈷、鉬、鈦、釩等化合物。對油頁巖熱解行為研究可以為工業化生產頁巖油提供理論指導,但油頁巖中礦物質含量高(60%~80%)、有機質與礦物質緊密交合、有機質與礦物質熱解同時發生等因素對油頁巖熱解機理的研究造成很大障礙,因此,針對礦物質對干酪根熱解影響的研究顯得尤為重要[2-3]。

諸多文獻表明礦物質中黏土礦物對有機質熱解影響較為顯著[4-6]。而蒙脫石作為黏土礦物之一對干酪根熱解影響作用最為顯著。劉曉艷[7]應用自制熱解模擬裝置研究了蒙脫石對褐煤熱解烴過程的影響,認為蒙脫石對有機質熱解存在重要作用,在低成熟度時,蒙脫石對有機質生油率顯著提高。JOHNS等[8]對脂肪酸與蒙脫石摻混的實驗樣品進行了熱解實驗,研究表明,蒙脫石對脂肪酸脫羧基形成長鏈烷烴以及長鏈烷烴進一步裂解生成短鏈烷烴具有催化作用。李術元等[9]對烴源巖催化生烴機制進行了綜述,認為在有水條件下,蒙脫石對有機質生烴催化作用最強,主要是通過改變反應的路徑降低反應的活化能而實現的。而ESPITALIE等[10]通過熱模擬實驗,對干酪根以及干酪根與礦物質(伊利石、蒙脫石、高嶺石、方解石)摻混樣品進行熱解實驗,研究表明礦物催化劑對干酪根熱解具有抑制作用。HETéNYI[11]利用熱解模擬實驗,研究了蒙脫石兩種類型干酪根熱解影響,發現蒙脫石對Ⅲ型干酪根熱解產物釋放具有抑制作用,但對Ⅰ型干酪根存在催化作用。由以上可知,蒙脫石催化機理基本都是從地球化學角度研究有機質生烴過程而來,即使是研究熱解反應過程中蒙脫石對干酪根產物的影響,也因為樣品單一及樣品處理方式的不同而存在正催化和反催化兩種結論,因此,深入蒙脫石對油頁巖干酪根熱解特性影響的研究具有實際意義。

本文應用TG-FTIR聯用技術對中國樺甸、撫順、窯街油頁巖干酪根進行不同升溫速率下蒙脫石摻混共熱解實驗研究,以期得到蒙脫石對油頁巖干酪根熱解行為的影響,為蒙脫石對油頁巖干餾產油影響提供基礎數據。

1 實驗部分

1.1 樣品采集

蒙脫石樣品采自吉林地區,并利用四分法進行取樣,經破碎研磨、篩選至0.2mm以下,置于

40℃恒溫烘干箱內干燥至恒重,裝瓶待用。

油頁巖樣品采自吉林樺甸(HD)、遼寧撫順(FS)及甘肅窯街(YJ)3個地區,經破碎研磨、篩選至0.2mm以下,置于40℃恒溫烘干箱內干燥至恒重,裝瓶待用。

1.2 干酪根的制備

依據沉積巖中干酪根分離方法GB/T 19144—2010的國家標準,制取油頁巖干酪根。取15g油頁巖樣品,首先對原頁巖樣進行索氏抽提去掉可溶瀝青,對去瀝青樣品加入150mL 20%的HCl溶液,置于恒溫70℃水浴振蕩器中勻速震蕩4h,將初次酸洗后的樣品(經)去離子水反復過濾,直至濾液用AgNO3滴定無白色沉淀。按照1∶1比例加入40%HF和20%HCl的混合酸溶液,重復上述試驗;然后加入20%的HNO3溶液,常溫酸洗2h,最終得到各地區油頁巖的干酪根樣品。各地區干酪根XRD譜圖如圖1所示,由圖1可知,各地區干酪根礦物質峰基本消失,表明干酪根符合實驗要求。

按照國標GB/T 476—2008方法對干酪根組成元素進行測定,C、H、N、S含量為3次平行樣的平均值,O含量用差減法獲得,測定結果見表1。

1.3 實驗儀器及方法

應用瑞士METTLER-TOLEDO公司生產的TGA/DSC熱重分析儀和美國ThermoFisher科技公司生產的NICOLET IS10傅里葉變換紅外光譜儀組成的TG-FTIR聯用分析儀進行試驗。選定熱解升溫速率分別為10℃/min、20℃/min、50℃/min,采用氮氣(99.999 %)作為載氣,流量為50mL/min,溫度區間50~850℃,并于常壓下運行。傅里葉紅外光譜儀的光譜測量范圍為4000~500cm?1,分辨率為4cm?1,掃描次數為16,掃描速度為1.8988cm/s。

圖1 各地區干酪跟XRD譜圖

表1 干酪根元素分析

1.4 數據處理

1.4.1 熱解失重率定量分析方法

由TG曲線可以得到w-t關系,將樣品在熱解前后的凈失重與樣品所含干酪根質量的比值作為熱解失重率?,如式(1)。

式中,?為干酪根熱解失重率,%;W1為樣品初始質量,g;Wτ為樣品在τ時刻質量,g;Ws為τ時刻蒙脫石單獨熱解失重量,g;α為干酪根所占樣品比例,%。

1.4.2 紅外定量分析方法

采用外標法測定熱解產物中甲基(CH3)、亞甲基(CH2)相對含量,所選標準物品為含100%脂肪鏈微晶石蠟[12-13]。為了確定標準曲線,重復進行5次實驗,所得標準曲線誤差為3.7%,符合實驗要求。

2 結果與討論

2.1 熱解產物析出特性

2.1.1 各樣品熱解失重率實驗

不同升溫速率下干酪根與蒙脫石共熱解產物失重率實驗結果見表2。由表可知,各干酪根單獨熱解時,升溫速率的改變并未引起熱解失重率顯著改變,但有蒙脫石摻混的HD干酪根熱解失重率隨著升溫速率提高略有提高,FS、YJ干酪根與蒙脫石摻混前后升溫速率對熱解失重率基本無影響;在同一升溫速率下,HD和FS干酪根熱解失重率隨著蒙脫石配比增加均呈現先增加后減少的趨勢,且YJ干酪根在蒙脫石摻混比為5時熱解失重率反而比干酪根單獨熱解有所下降;3個地區樣品熱解失重率達到最大值時蒙脫石摻混比也不同,HD和FS干酪根均在蒙脫石摻混比為1時轉化率達到最大,YJ干酪根在摻混比0.5時熱解失重率達到最大,該規律與張枝煥等[14]的研究結果一致。

表2 不同升溫速率下蒙脫石摻混的干酪根熱解失重率

對于有雙因素單獨觀測值的實驗,可以對所得數據進行方差分析得到各因素對目標值的顯著性檢驗[15]。各地區有蒙脫石摻混的干酪根熱解失重率方差分析結果見表3。

對各地區樣品方差分析可知,HD樣品FTemperature=5.08>F0.05(4,8)=3.84、FSmectite=10.07>F0.05(2,8)=4.46,FS樣品FTemperature=3.86>F0.05(4,8)=3.84、FSmectite=5.48>F0.05(2,8)=4.46,表明蒙脫石配比對HD、FS油頁巖干酪根熱解過程有顯著影響,而YJ干酪根FTemperature=1.43

表3 各地區有蒙脫石摻混的干酪根熱解失重率方差分析

2.1.2 不同熱解終溫的熱解失重率實驗

對升溫速率為10℃/min的干酪根樣品各升溫段熱解失重率進行分析,結果見表4。干酪根在320℃左右熱解產物開始析出,且在440~500℃之間析出速率迅速增大,有蒙脫石摻混的3個地區干酪根具有同樣的趨勢。在熱解過程中,440℃前干酪根與蒙脫石摻混樣品的熱解失重率較干酪根單獨熱解低,而440℃后其熱解失重率較干酪根熱解有所提高,表明在干酪根熱解早期蒙脫石層間吸附作用抑制了熱解產物的析出,但隨著熱解程度加深,在蒙脫石活化中心的催化作用下干酪根熱解產物析出速率增大,熱解失重率迅速提高。因此在干酪根的整個熱解過程中,蒙脫石對其熱解產物釋放具有先抑制后催化的作用[6]。

表4 10℃/min各地區樣品不同終溫熱解失重率

2.2 熱解產物紅外光譜

干酪根與有蒙脫石摻混的干酪根熱解產物紅外三維譜圖如圖2所示。由圖2可知,油頁巖干酪根單獨熱解與有蒙脫石摻混的干酪根熱解產物官能團種類基本相同,且主要以甲基(CH3)、亞甲基(CH2)為主。

圖3為HD干酪根480℃熱解產物紅外譜圖,由圖可知干酪根熱解產物中以脂肪烴為主,在2960cm?1、1460cm?1處分別為CH3非對稱伸縮振動及反對稱變形振動,2926cm?1、2860cm?1處依次為CH2對稱伸縮振動與反對稱伸縮振動,3015cm?1處為甲烷振動特征峰。

圖2 熱解產物紅外三維譜圖

圖3 HD干酪根480℃熱解產物紅外譜圖

CH3與CH2相對含量可以反映干酪根熱解產物中脂肪烴鏈長度與支化情況[16]。按照前述方法得到不同升溫速率下蒙脫石與干酪根共熱解產物中CH3與CH2相對含量,如表5所示。由表5可知,隨升溫速率的提高,各干酪根熱解產物中CH3/CH2值都有不同程度升高,表明隨著升溫速率的提高干酪根熱裂解更加劇烈,使得熱解產物中短鏈脂肪烴不斷增多。但不同干酪根其CH3/CH2值升高數值略有差異,因為不同干酪根熱解對升溫速率響應程度不同;此外,在相同升溫速率下各干酪根均表現隨著蒙脫石配比的增加CH3/CH2值也有所增加,所得熱解產物中脂肪烴鏈長度變短,表明蒙脫石在干酪根熱解過程中起到了催化裂解作用。張枝煥等[4]研究結果同樣驗證蒙脫石的添加使得干酪根熱解產物正構烷烴減少,異構烷烴增多。干酪根熱解主要是以自由基的機理進行,熱解產物中正構烷烴和烯烴較多,且有部分產生縮合和環化反應形成少量環烷烴和芳烴;而有蒙脫石摻混的干酪根熱解是一種催化降解過程,是以正碳離子反應機理進行的,蒙脫石含有Lewis酸和Bronsted酸,可以分別將烷烴和芳烴引發生成正碳離子,生成的正碳離子可以引發兩種反應,一是發生β位斷裂生成小分子烴,二是正碳離子可以引發分子內部氫遷移以及烷基和芳基的位置變動,從而引起各種碳骨架的異構化反應,而且比自由基反應機理異構化速度反應快得多[7,17-20]。因此,有蒙脫石摻混的干酪根樣品熱解產物中CH3/CH2值要比干酪根單獨熱解大,且這種作用隨蒙脫石配比加大而增強。

表5 熱解產物中CH3/CH2相對比值

2.3 熱解動力學分析

國內外諸多學者對油頁巖進行熱解動力學研究,得到多方面的結論[21-26]。根據熱重實驗數據探討了在10℃/min條件下了3個地區油頁巖干酪根的熱解動力學,運用Coats-Redfem積分法[27]求取各干酪根與蒙脫石摻混下的熱解反應動力學參數。其熱解速率函數如式(2)。

假設反應為一級反應動力學,可取f(α)=1?α,則利用用Coats-Redfem積分方法分離變量整理并取近似值可以得到式(3)。

圖4 各地區樣品干酪根與蒙脫石共熱解Arrhenius曲線

圖5為不同蒙脫石配比的干酪根熱解活化能與失重率對比。由圖5可知,HD和FS干酪根熱解活化能與失重率呈負相關,即當失重率升高時,在蒙脫石同一配比下的活化能則有所降低,當失重率降低時,在蒙脫石同一配比下的活化能則有所升高,這表明蒙脫石的催化使得熱解活化能降低,增大了熱解產物的析出量,進而熱解失重率升高;而YJ地區干酪根隨著蒙脫石配比增加其熱解活化能與失重率分別保持在同一水平,說明蒙脫石對YJ干酪根總的影響效果并不顯著。這是因為蒙脫石是2∶1型層狀硅鋁酸鹽,具有較強的吸附能力和離子交換能力,對干酪根熱解影響有兩方面,一是表面吸附和層間吸附作用,二是酸性活化中心的正催化作用。由表4實驗結果可知,蒙脫石對窯街樣品熱解產物具有吸附作用,并且從表5可以得到YJ干酪根熱解過程中CH3/CH2值隨著蒙脫石摻混比增加而升高,表明蒙脫石對YJ干酪根同樣存在催化裂解作用。因此,出現上述現象很可能與YJ地區油頁巖樣品有關,在酸洗過程中YJ酸洗液成黑色渾濁狀且在過濾紙上會殘留較多成黑色的有機質,表明YJ干酪根成熟度不夠高,有機大分子較易裂解成小分子,而生成的小分子較易被蒙脫石層間結構吸附,抵消了蒙脫石的催化作用,導致蒙脫石影響作用不夠顯著。因此,蒙脫石對油頁巖干酪根存在物理吸附和裂解催化兩種作用,并且兩種作用強度又與干酪根本身大分子結構的復雜性和差異性有關。

表6 各地區樣品干酪根與蒙脫石共熱解動力學參數

圖5 不同蒙脫石配比的干酪根熱解活化能與失重率

此外,隨著蒙脫石摻混比例增大,HD和FS干酪根表觀活化能呈先減少后增大趨勢且始終比干酪根單獨熱解小,說明當蒙脫石含量較少時由于蒙脫石酸性活化中心的作用,使得干酪根的催化裂解反應更易進行,反應活化能減少,熱解失重率增大;但隨著蒙脫石摻混比例增大,吸附焦化作用也在增強,使得反應活化能增大及釋放出的烴量減少,降低了干酪根熱解失重率。

總之,蒙脫石對油頁巖干酪根熱解同時具有催化和抑制兩種作用,這兩種作用隨著熱解過程的加深和蒙脫石摻混比的改變而不同。在熱解初期,蒙脫石對干酪根熱解產物的影響主要是吸附作用,而隨著熱解過程的加深,蒙脫石酸性活化中心對干酪根熱裂解催化作用較為顯著,熱解產物失重率變大;蒙脫石與干酪根低摻混比時蒙脫石催化裂解作用較強,熱解失重率較干酪根熱解大,但隨著蒙脫石摻混比例的增加,熱解產物的吸附作用增強,熱解失重率反而降低,且降低程度因干酪根樣品不同而有所差異。因此,蒙脫石配比對不同地區干酪根樣品熱解影響各異,但低摻混比對干酪根熱解均表現正催化。

3 結 論

(1)升溫速率對3個地區油頁巖干酪根樣品熱解失重率并無顯著影響,但是在同一升溫速率下3個地區干酪根隨著蒙脫石配比增加均呈現先增加后減少的趨勢,且HD和FS干酪根均在蒙脫石摻混比為1時熱解失重率最大,方差分析可知,升溫速率和蒙脫石添加對YJ地區干酪根樣品熱解并無顯著影響,這是由于YJ干酪根本身特征決定的。

(2)蒙脫石摻混的3個地區干酪根在熱解溫度達到440℃之前其熱解失重率較干酪根單獨熱解低,而440℃之后熱解失重率迅速提高并較干酪根單獨熱解有所提高。表明蒙脫石對其熱解產物釋放起到先抑制后催化的作用。

(3)每種干酪根隨升溫速率提高CH3/CH2比值都有不同程度升高,此外,在相同升溫速率下3個地區干酪根均表現隨著蒙脫石配比的增加CH3/CH2比值有所增加,表明蒙脫石對油頁巖干酪根熱解表現為催化異構。

(4)蒙脫石對油頁巖干酪根熱解產物存在物理吸附和裂解催化兩種影響,且隨著熱解過程的加深和蒙脫石配比的改變而不同,蒙脫石與干酪根摻混比例較小時對蒙脫石催化裂解作用較強,熱解失重率較干酪根單獨熱解大,但隨著蒙脫石摻混比例的增加熱解產物的吸附作用增強,使得熱解失重率降低。

參 考 文 獻

[1] 錢家麟,王劍秋,李術元. 世界油頁巖資源利用和發展趨勢[J]. 吉林大學學報(地球科學版),2006,36(6):878-887.

[2] MOHAMMAD A. Effect of demineralization and heating rate on the pyrolysis kinetics of Jordanian[J]. Fuel Processing Technology,2011,92(9):1805-1811.

[3] YAN J W,JIANG X M. A TG-FTIR investigation to the catalytic effect of mineral matrixin oil shale on the pyrolysis and combustion of kerogen[J]. Fuel,2013,104(3):307-317.

[4] 張枝煥,高先志. 黏土礦物對干酪根熱解產物的影響及其作用機理[J]. 石油大學學報,1995,19(5),11-17.

[5] GUO H F,JIA D L,YANG Y D. Effect of minerals on the self-heating retorting of oil shale:self-heating effect and shale-oil production[J]. Fuel,2014,118:186-193.

[6] 李忠. 試論油氣形成過程中黏土礦物的催化作用[J]. 石油實驗地質,1992,14(1):59-63.

[7] 劉曉艷. 黏土礦物對有機質演化的影響[J]. 天然氣地球科學,1995,6(3):23-26.

[8] JOHNS W D,SHIMOYAMA A. Clay minerals and petroleumforming reactions during burial and diagenesis[J]. AAPG Bulletin,1972,56(11):2160-2167.

[9] 劉洛夫,李術元. 烴源巖催化生烴機制研究進展[J]. 地質論評,2000,46(5):491-498.

[10] ESPITALIE J,MADEC M,TISSOT B. Role of mineral matrix in kerogen pyrolysis:influence on petroleum generation and migration[J]. AAPG Bulletin,1980,64(1):59-66.

[11] HETéNYI M. Simulated thermal maturation of type I and III kerogens in the presence,and absence,of calcite and montmorillonite [J]. Organic Geochemistry,1995,23(2):121-127.

[12] 翁詩甫. 傅里葉變換紅外光譜分析[M]. 2版. 北京:化學工業出版社,2010:263-264.

[13] 常建華,董綺功. 波譜原理及解析[M]. 2版. 北京:科學出版社,2005:107-109.

[14] 張枝煥,高先志. 黏土礦物對干酪根熱解生烴過程的影響[J]. 石油勘探與開發,1994,21(5):30-37.

[15] 張玲. 單因素及雙因素方差分析及檢驗的原理及統計應用[J]. 數學學習與研究(教研版),2010(7):92-94.

[16] 薛向欣,劉艷輝,李勇,等. 紅外光譜法研究油頁巖及干酪根的生油能力[J],東北大學學報(自然科學版),2010,31(9):1292-1295.

[17] EGLINTON T I,ROWLAND S J,CURTIS C D. Kerogen-mineral reaetions at raised temperatures in the presence of water[J]. Org. Geochem.,1986,10:1041-1052.

[18] 高先志,張萬選,張厚福. 礦物質對熱解影響的研究[J]. 石油實驗地質,1990,12(2):201-205.

[19] 房玄,雷懷彥. 黏土礦物對生物-熱催化過渡帶有機質熱解影響的研究[J]. 天然氣地球科學,1993,4(6):80-85.

[20] 吳平霄. 蒙脫石熱活化及其與微結構變化關系研究[D]. 廣州:華南理工大學,1998.

[21] TORRENTE M C,GALAN M A. Kinetics of the thermal decomposition of oil shale from Puertollano[J]. Fuel,2001,80(3):327-334.

[22] YUE S,LI C. Study of pyrolysis kinetics of oil shale[J]. Fuel,2003,82(3):337-342.

[23] KêK M V,PAMIR M R. Comparative pyrolysis and combustion kinetics of oil shales[J]. J. Anal. Appl. Pyrol.,2000,55(2):185-194.

[24] KHRAISHA Y H,SHABIB I M. Thermal analysis of shale oil using thermogravimetry and differential scanning calorimetry[J]. Energy Convers. Mgmt.,2002,43(2):229-239.

[25] JABER J O,PROBERT S D. Non-isothermal thermogravimetry and decomposition kinetics of two Jordanian oil shales under different processing conditions[J]. Fuel Process Tech.,2000,63(1):57-70.

[26] WANG Q,LIU H P,SUN B Z,et al. Study on pyrolysis characteristics of Huadian oil shale with isoconversional method[J]. Oil Shale,2009,26(2):148-162.

[27] COATS R. Kinetic parameters from thermo gravimetric data[J]. Nature,1964,201:68-69.

研究開發

Effect of smectite on the pyrolysis of kerogen isolated from oil shale

WANG Qing,ZHANG Hongxi,CHI Mingshu,CUI Da,XU Xiangcheng
(School of Energy and Power Engineering,Northeast Dianli University,Jilin 132012,Jilin,China)

Abstract:To examine the effect of smectite on oil generation of kerogen pyrolysis,we studied the kerogen isolated from oil shale of Huadian,Fushun and Yaojie by TG-FTIR and kinetic analysis. Results show that the weight loss of kerogen from Huadian and Fushun are reduced after the initial increase. Furthermore,the weight loss of kerogen which blended with smectite is lower than that of kerogen until the pyrolysis temperature reaches 440℃,but it then increases rapidly and exceeds that of kerogen after 440℃. With the increase of smectite,the ratio of CH3/CH2is also increased at same heating rate. The pyrolysis activation energy of kerogen with smectite is found lower than that of kerogen alone. Results show that the effect of smectite to kerogen pyrolysis is both physical adsorption and catalytic cracking,which will be different through the pyrolysis process and increasing blend of smectite. Finally,further discuss has been conducted for the mechanism of smectite by pyrolysis kinetics analysis of samples with a heating rate of 10℃/min。

Key words:smectite; kerogen; pyrolysis; oil shale; kinetics; catalytic cracking

基金項目:國家自然科學基金項目(51276034)。

收稿日期:2015-08-19;修改稿日期:2015-10-29。

DOI:10.16085/j.issn.1000-6613.2016.03.018

中圖分類號:TK 16

文獻標志碼:A

文章編號:1000–6613(2016)03–0766–07

第一作者及聯系人:王擎(1964—),男,博士,教授。E-mail rlx888@126.com。

主站蜘蛛池模板: 亚洲国产日韩一区| 久久亚洲精少妇毛片午夜无码| 无码粉嫩虎白一线天在线观看| 久久久精品久久久久三级| 国产免费一级精品视频 | 性视频久久| a毛片免费观看| 欧美、日韩、国产综合一区| 在线另类稀缺国产呦| 久青草国产高清在线视频| 曰AV在线无码| 91久久性奴调教国产免费| 最新日韩AV网址在线观看| 国产美女免费| 高潮爽到爆的喷水女主播视频 | 在线色国产| 国产精品偷伦在线观看| 97精品久久久大香线焦| 97精品伊人久久大香线蕉| 美女免费黄网站| 午夜在线不卡| 69av在线| 日韩中文精品亚洲第三区| 色综合成人| 国产精品黑色丝袜的老师| 欧美国产视频| 在线观看免费黄色网址| 日韩精品久久无码中文字幕色欲| 亚洲天堂久久久| 国产情侣一区| 久久婷婷国产综合尤物精品| 免费看黄片一区二区三区| 一级毛片基地| 九色综合视频网| 国产女人在线视频| 无码福利视频| 4虎影视国产在线观看精品| 国产午夜无码专区喷水| 99久久精彩视频| 91精品专区国产盗摄| 看国产一级毛片| 亚洲精品爱草草视频在线| 米奇精品一区二区三区| 欧美激情,国产精品| 97se亚洲| 99视频在线免费观看| 亚洲无码91视频| 亚亚洲乱码一二三四区| 中文字幕一区二区人妻电影| 免费在线国产一区二区三区精品| 91高清在线视频| 亚洲高清在线天堂精品| 在线精品自拍| 国产成人调教在线视频| 尤物国产在线| 国产精品福利在线观看无码卡| 另类综合视频| 亚洲综合一区国产精品| 亚洲欧美日韩中文字幕在线一区| 亚洲欧美另类久久久精品播放的| 色窝窝免费一区二区三区 | 久久免费看片| 九色国产在线| 国产18在线播放| 亚洲黄色片免费看| 中文字幕第1页在线播| 国产美女叼嘿视频免费看| 一级全免费视频播放| 日本尹人综合香蕉在线观看| 怡春院欧美一区二区三区免费| 国产精品专区第1页| 99视频精品全国免费品| 成人在线欧美| 国产 在线视频无码| 国产极品美女在线播放| 国产精品55夜色66夜色| swag国产精品| 国产高潮流白浆视频| 国产精品尤物铁牛tv| 成人国产免费| 88av在线| 亚国产欧美在线人成|