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

渣油加氫反應過程中芳香性化合物轉化規律

2018-03-05 05:39:42張龍剛戴立順楊清河賈燕子
石油學報(石油加工) 2018年1期

張龍剛, 聶 紅, 戴立順, 楊清河, 賈燕子, 劉 濤

(中國石化 石油化工科學研究院, 北京 100083)

隨著輕質燃料需求的日益增加,渣油深度轉化成為煉油工業長期追求的目標。固定床渣油加氫和催化裂化組合工藝作為實現重油高效轉化的有效手段,在煉油廠得到了廣泛應用[1],其中固定床渣油加氫裝置的主要目的是脫除渣油中的硫、氮和金屬雜原子,降低殘炭值,從而為催化裂化裝置提供合格原料[2]。渣油中硫、氮原子多以芳香性雜環化合物形式存在,在后續加工過程中會造成催化劑失活。渣油中殘炭前身物主要存在于五環及五環以上的稠環芳烴類化合物中[3],其容易在催化裂化過程中生焦,并影響催化裂化裝置的產品分布。

囿于分析手段限制,原有對固定床渣油加氫轉化過程規律的研究[4-6]主要集中在渣油四組分和雜原子含量等性質上,而缺少從分子層面對渣油中芳香性化合物加氫轉化過程的有效認識。近年來傅里葉離子回旋共振質譜(FT ICR MS)在渣油烴類組成分析和應用方面取得了重大進展,采用電噴霧(ESI)電離源[7]可以選擇性電離微量的極性雜原子化合物,大氣壓光致電離(APPI)[8]可以有效電離非極性的烴類化合物和極性較弱的含硫類化合物。高分辨質譜可以獲得渣油中化合物類型、碳數和缺氫數分布等詳細信息[9-10],是認識渣油分子結構組成及加氫轉化規律的有效手段。

筆者在中型渣油加氫裝置上考察了兩種典型渣油的反應差異,采用傅里葉離子回旋共振質譜分析方法,比較沿反應器物流方向油品中各類型化合物的組成分布變化規律,旨在從分子水平逐步深入認識渣油中芳香性化合物的加氫轉化規律,從而為更高效的渣油加氫催化劑和工藝的開發提供科學參考。

1 實驗部分

1.1 原料油

選用兩種不同來源典型常壓渣油作為原料油,分別為儀長管輸原油常壓渣油(以下簡稱儀長渣油)和沙特阿拉伯輕質原油常壓渣油(以下簡稱沙輕渣油),其基本性質如表1所示。由表1可以看出,兩類渣油密度和殘炭值相近,其中儀長渣油屬于典型的高氮低硫類渣油,而沙輕渣油屬于高硫低氮類渣油。

1.2 試驗方法

渣油加氫試驗在中型加氫裝置上進行,其工藝流程簡圖如圖1所示。第一反應器R1中裝填加氫保護劑(RG-30B)和加氫脫金屬劑(RDM-35/RDM-32/RDM-33B),第二反應器R2中裝填加氫脫硫催化劑(RMS-30)和加氫脫殘炭催化劑(RCS-30/RCS-31B),所用催化劑均為中國石化石油化工科學研究院開發的已工業生產的第三代RHT系列催化劑。在裝置運行3個月催化劑活性進入穩定期后,分別通入儀長渣油和沙輕渣油,反應后采集第一反應器R1出口脫金屬后的油樣(記作R1樣品)和第二反應器R2出口加氫后油樣(記作R2樣品)。

圖1 渣油加氫工藝流程示意圖Fig.1 Flow diagram of the residue hydrotreating process

1.3 渣油分子組成分析方法

渣油分子組成分析采用apex-Qe型傅里葉變換離子回旋共振質譜儀(9.4 T),離子源為APPI大氣壓光致電離源。采用正離子模式,樣品溶解在甲苯(色譜純,美國迪馬公司生產)中,質量濃度為0.5 mg/mL,由注射泵以360 μL/h的速率送入儀器。噴霧氣溫度400℃,流速2.0 L/min。干燥氣溫度220℃,流速6.0 L/min。測試參數:相對分子質量范圍200~1200,射頻激發RF為250Vp-p,一級分離質量Q1 Mass為250,碰撞電壓為-4 V,采樣時掃描256次。

2 結果與討論

2.1 儀長和沙輕渣油加氫前后基本性質變化規律

在相同催化劑級配體系和工藝條件下,考察兩類渣油加氫前后基本性質的變化規律,結果如表1所示。從表1可以看出,沿物流方向不同反應器出口渣油密度、雜原子含量和殘炭值均呈現逐漸降低的趨勢,氫含量逐漸提高,說明隨加氫深度增加,渣油性質得到改善。

兩類渣油R1樣品中金屬鎳和釩原子大部分已經脫除,硫含量明顯降低,渣油殘炭值降低,儀長渣油脫殘炭率為21.2%,沙輕渣油脫殘炭率為22.5%,說明渣油殘炭前身物經過保護劑和加氫脫金屬劑反應后發生了部分轉化。R2樣品中雜原子含量進一步降低,儀長渣油和沙輕渣油加氫脫殘炭率分別提高到58.7%和67.5%,加氫脫硫率分別為86.6%和89.7%,加氫脫氮率分別為45.1%和40.0%。由此可見,R1樣品中儀長脫金屬渣油脫殘炭率略低于沙輕脫金屬渣油,R2樣品中兩類渣油脫殘炭率差距進一步加大,儀長加氫后渣油仍然呈現高氮低硫的特點,這一結果與文獻[11]報道一致。兩類渣油加氫脫硫率相差不大,但儀長渣油加氫脫氮率較高,加氫脫殘炭率相對較低,這可能是由于儀長渣油平均相對分子質量較大的原因,并且與渣油中高含量氮原子的存在形式也有關。

表1 沿反應器物流方向兩類渣油基本性質變化規律Table 1 Bulk properties of two kinds of residue along the reactor flow direction

Residue hydrotreating conditions:p=16.0 MPa;T=390℃; LHSV=0.20 h-1;V(H2)/V(Oil)=700

2.2 渣油加氫反應過程中芳香性化合物的組成分布變化規律

采用傅里葉離子回旋共振質譜對渣油加氫反應過程中芳香性化合物的組成和分布變化規律進行詳細分析,圖2為不同加氫深度下兩類渣油中芳香性化合物相對含量的變化規律。從圖2可以看出,儀長渣油中芳香性化合物以CH類化合物(無雜原子的烴類)、N1類化合物(含有1個氮原子的烴類,以下類推)和S1類化合物(含有1個硫原子的烴類)為主,分別占渣油總量的47.9%、33.3%和11.7%。沙輕渣油中以CH類化合物、S1類化合物和S2類化合物為主,分別占渣油總量的31.5%、47.1%和15.2%。在多雜原子化合物中儀長渣油以含有氮原子的N2類、NS類和NO類化合物居多,沙輕渣油中以含有硫原子的S2類和NS類化合物居多,與兩類渣油中硫、氮雜原子的含量特點相符。

圖2 不同反應器出口兩類渣油芳香性化合物分布 Fig.2 Aromatic compound distribution of two kinds of residue at different reactor outlets (a) Yichang AR; (b) Saudi-Arabian AR

R1樣品中儀長渣油CH類化合物占渣油中芳香性化合物總量的58.1%,沙輕渣油中CH類化合物占41.3%,較渣油原料中均有所提高。兩類渣油中S1類、S2類和NS類化合物含量有所下降,說明經過第1反應器R1后脫除了硫原子,含硫化合物發生了部分轉化。儀長渣油中N1類化合物含量有所下降,但下降幅度不大。沙輕渣油中N1類化合物相對含量反而有所增加,這可能與渣油中NS類和NO類等化合物脫除了1個硫原子和1個氧原子,轉變為N1類化合物有關。

R2樣品中儀長渣油CH類化合物含量為68.8%,沙輕渣油中為81.7%,兩類渣油中S1類化合物在第二反應器R2中大部分被轉化。隨著反應深度增加,兩類渣油殘炭值降低,對應CH類化合物含量升高,儀長加氫渣油中S1類和N1類化合物含量下降,而沙輕加氫渣油中S1類和S2類化合物含量下降。通過加氫反應,渣油中含雜原子的S1類和N1類等芳香性化合物加氫脫除硫、氮雜原子,進而轉化為CH類化合物,同時芳環縮合度降低,在此過程中渣油殘炭值逐步降低,進一步說明殘炭轉化過程與雜原子的脫除以及芳環加氫飽和密切相關。

對比圖2中兩類渣油的加氫轉化結果可以看出,R1樣品中兩類渣油中CH類化合物轉化程度差別較小,儀長渣油和沙輕渣油中CH類化合物含量分別提高了10.2和9.8百分點。R2樣品中不同類型化合物相對含量發生了顯著變化,儀長渣油中CH類化合物含量提高了20.9百分點,N1類化合物和S1類化合物含量分別下降了6.7和10.0百分點;而沙輕渣油中CH類化合物含量提高了50.2百分點,S1類和NS類化合物含量顯著下降,其中S1類化合物含量下降了40.8百分點。兩類渣油經過加氫后CH類化合物含量增加幅度不同,沙輕渣油中高含量的S1類化合物顯著降低,儀長渣油中高含量的N1類化合物降低幅度相對較小,這種差異影響了殘炭前身物的轉化,R1樣品中儀長渣油加氫脫殘炭率低于沙輕渣油,隨加氫深度增加,兩類渣油反應差異更加明顯。

2.3 兩類渣油加氫反應前后主要化合物碳數及DBE值變化規律

由前面的討論可知,兩類渣油加氫前后CH類、S1類和N1類化合物含量較高,是影響渣油雜原子脫除和殘炭前身物轉化的關鍵化合物,因此有必要對渣油中3類主要化合物的分子結構變化過程進行詳細分析。圖3、圖4和圖5分別為兩類常壓渣油原料、R1樣品和R2樣品中CH類、S1類和N1類化合物的碳數-DBE值(化合物分子等效雙鍵數目)分布氣泡圖,圖中氣泡的大小代表對應化合物含量多少。DBE值代表渣油分子的等效雙鍵數目,DBE值越大,說明化合物縮合度越高,DBE值高于15的化合物一般為五環以上的稠環芳烴類化合物,可以認為是主要殘炭前身物。

圖3 兩類常壓渣油原料中主要化合物中碳數-DBE分布Fig.3 Plots of DBE vs carbon number of different class compounds in two kinds of atmospheric residues (a), (b), (c) stands for CH, S1 and N1 compounds of Yichang AR and (d), (e), (f) stands for CH, S1 and N1 compounds of Saudi-Arabian AR

從圖3可以看出,兩類渣油原料中CH類化合物碳數分布在20~70之間,DBE值集中在4~16。S1類化合物碳數分布在17~74之間,DBE值集中在5~18。N1類化合物碳數分布在22~78之間,DBE值集中在8~23。圖3結果說明,渣油芳香性化合物中N1類化合物縮合度最高,其次是S1類化合物,CH類化合物縮合度相對最低。儀長渣油中CH類、S1類和N1類化合物的平均碳數和DBE值均高于沙輕渣油,說明儀長渣油分子更大,芳環縮合度更高。

從圖4和圖5可以看出,隨加氫深度增加,兩類渣油R1和R2樣品中CH類和S1類化合物逐漸向低碳數和低DBE值方向移動,N1類化合物碳數和DBE值的變化并不明顯。儀長渣油R1樣品中碳數高于60的S1類化合物含量較少,進一步說明一部分含硫化合物在經過加氫脫金屬劑后被脫除。沙輕渣油R2樣品中仍存在一些高碳數的S1類化合物,這部分化合物在經過第二反應器R2后大部分被轉化,這可能是沙輕渣油中S1類化合物本身含量很高,在第一反應器R1中未完全反應造成的。R1和R2樣品中高碳數N1類化合物含量降低幅度并不明顯,這是由于渣油中氮原子多存在于芳環中,而C=N雙鍵鍵能較高[12],無法通過芳環間橋鍵斷裂的方式由大分子轉化為小分子。可以看出,R2樣品中3類化合物中仍含有部分DBE值大于15的化合物,其中N1類化合物中比例最高,儀長渣油中占54.6%,沙輕渣油中占63.8%,說明兩類渣油N1類化合物中均含有較高比例的殘炭前身物。

圖4 兩類常壓渣油R1樣品中主要化合物中碳數-DBE分布Fig.4 Plots of DBE vs carbon number of different compounds in two residues at the outlet of 1st reactor (a), (b), (c) stands for CH, S1 and N1 compounds of Yichang AR and (d), (e), (f) stands for CH, S1 and N1 compounds of Saudi-Arabian AR

圖5 兩類常壓渣油R2樣品中主要化合物中碳數-DBE分布Fig.5 Plots of DBE vs carbon number of different compounds in two residues at the outlet of 2nd reactor (a), (b), (c) stands for CH, S1 and N1 compounds of Yichang AR and (d), (e), (f) stands for CH, S1 and N1 compounds of Saudi-Arabian AR

進一步對儀長渣油中高含量的N1類化合物和沙輕渣油中高含量的S1類化合物進行分析,詳細比較了兩類關鍵化合物在第二反應器R2中加氫脫硫劑和加氫脫殘炭劑加氫反應前后的分子結構變化,其碳數和DBE值變化分別如圖6和圖7所示。從圖6可以看出,經過第二反應器R2后,儀長渣油中N1類化合物中碳數大于60的化合物在加氫后比例有所上升,并且出現了一些R1樣品中沒有的碳數大于85的化合物,同時DBE值高于26的化合物加氫后未明顯轉化,一方面是由于高縮合度的含氮化合物本身較難加氫轉化,另一方面渣油中含氮化合物極性較強,容易吸附在催化劑表面生成積炭[13-14],不利于催化劑發揮加氫作用。

從圖7可以看出,沙輕渣油中S1類化合物經過第二反應器R2加氫處理后,碳數和DBE值發生了明顯變化。經過加氫脫硫劑和加氫脫殘炭劑后,沙輕渣油中碳數高于60的化合物絕大部分轉化,DBE值高于30的化合物基本全部消失,這是因為渣油中不同芳環結構間橋鍵硫原子鍵能低,硫原子容易加氫脫除[1],在此過程中渣油大分子轉化為小分子,DBE值和殘炭值隨之降低,這也是高硫低氮類沙輕渣油殘炭前身物容易轉化的一大原因。

圖7 沙輕渣油不同反應器出口S1類化合物碳數和DBE變化Fig.7 Distribution of carbon number and DBE of Saudi-Arabian AR S1 compounds at different reactor outlets (a) Carbon number; (b) DBE

由此可見,兩類渣油分子結構及其雜原子含量和存在形式的不同是導致其殘炭前身物轉化存在差別的主要原因。與沙輕渣油相比,儀長渣油平均碳數和平均縮合度均更高,加氫生成油中CH類、S1類和N1類化合物均含有部分DBE值高于15的化合物,其中N1類化合物中高碳數和高DBE值的化合物所占比例最高。這部分高碳數高縮合度N1類化合物的存在一方面增加了反應空間位阻,抑制了渣油分子與催化劑活性中心的接觸,導致含氮化合物本身反應速率較低;另一方面高縮合度的含氮化合物容易吸附在催化劑表面,無法及時脫附,從而造成積炭增多,引起催化劑失活加快。分析工業廢催化劑積炭量可以發現,國內煉油廠中加工儀長渣油和沙輕渣油的渣油加氫裝置上其廢催化劑上平均積炭沉積質量分數分別為23.5%和15.5%。

渣油中芳香性化合物轉化難易程度不同,其中N1類化合物在加氫過程中最難轉化。與氮原子脫除相比,渣油中芳環結構單元間以S—S鍵和C—S鍵形式存在的硫原子極易脫除,從而導致殘炭值降低,而芳環中存在的硫原子其脫除速率也要優于渣油芳環中的氮原子。因此,儀長渣油中氮含量高、硫含量低,N1類化合物中較高比例的DBE值大于15的含氮芳香性化合物,是其殘炭脫除率低于沙輕渣油的主要原因。

基于以上對兩種典型渣油中主要化合物分子結構特點和轉化規律的研究,認為針對高氮低硫類渣油,開發性能更優的渣油加氫脫殘炭脫硫催化劑不僅要側重于提高加氫性能,還要注重減小大分子殘炭前驅物的擴散限制,強化渣油分子和催化劑活性中心的可接近性,以及含氮芳香性化合物的加氫脫除性能。

3 結 論

(1)儀長渣油加氫脫殘炭率低于沙輕渣油,隨加氫深度增加,兩類渣油反應差異更加明顯,儀長渣油加氫后仍呈現高氮低硫的特點。

(2)儀長加氫生成油中主要是CH類和N1類化合物,沙輕加氫生成油中主要是CH類化合物。儀長加氫渣油中存在高含量的N1類化合物,其中DBE值大于15的化合物比例較高,是其殘炭脫除率低于沙輕渣油的主要原因。

(3)儀長渣油中CH類、S1類和N1類化合物的平均碳數和DBE值均略高于沙輕渣油,隨加氫深度增加,兩類渣油CH類和S1類化合物逐漸向低碳數和低DBE值方向移動,N1類化合物的變化并不明顯。經過加氫脫殘炭劑后,沙輕渣油中S1類化合物碳數高于60和DBE值大于30的化合物大部分轉化,儀長渣油中N1類化合物出現了少量高碳數和高DBE值的化合物。

(4)針對高氮低硫類渣油,開發性能更優的渣油加氫脫殘炭脫硫催化劑,不僅要提高催化劑加氫性能,還要強化渣油分子與催化劑活性中心的可接近性,以及含氮芳香性化合物的加氫脫除性能。

[1] 聶紅, 楊清河, 戴立順, 等. 重油高效轉化關鍵技術的開發及應用[J].石油煉制與化工, 2012, 43(1): 1-6.(NIE Hong, YANG Qinghe, DAI Lishun, et al. Development and commercial application of key technology for efficient conversion of heavy oil[J].Petroleum Processing and Petrochemicals, 2012, 43(1): 1-6.)

[2] 李大東, 聶紅, 孫麗麗. 加氫處理工藝與工程[M].北京: 中國石化出版社, 2015.

[3] 陳俊武, 曹漢昌. 殘炭前身化合物的結構及其在煉油過程中的作用[J].煉油設計, 1992, 22(1): 1-11.(CHEN Junwu, CAO Hanchang. The structure of carbon residue precursor and its role in petroleum refining[J].Petroleum Refinery Engineering, 1992, 22(1): 1-11.)

[4] 趙輝. 渣油加氫轉化規律研究[D].青島: 中國石油大學, 2009.

[5] 孫昱東, 楊朝合, 山紅紅, 等. 不同來源渣油加氫反應性能的對比[J].石油學報(石油加工), 2011, 27(1): 32-36. (SUN Yudong, YANG Chaohe, SHAN Honghong, et al. Hydrotreating reaction performance of different residua[J].Acta Petrolei Sinica (Petroleum Processing Section), 2011, 27(1): 32-36.)

[6] 張會成, 顏涌捷, 趙榮林, 等. 渣油加氫處理過程中金屬分布與脫除規律的研究[J].石油煉制與化工, 2006, 37(11): 7-10. (ZHANG Huicheng, YAN Yongjie, ZHAO Ronglin, et al. Study of the metals distribution and removal in residue hydrotreating[J].Petroleum Processing and Petrochemicals, 2006, 37(11): 7-10.)

[7] 史權, 趙鎖奇, 徐春明, 等. 傅立葉變化離子回旋共振質譜儀在石油組成分析中的應用[J].質譜學報, 2008, 29(6): 367-378. (SHI Quan, ZHAO Suoqi, XU Chunming, et al. Fourier transform ion cyclotron resonance mass spectrometry and its application in petroleum analysis[J].Journal of Chinese Mass Spectrometry Society, 2008, 29(6): 367-378.)

[8] PURCELL J M, JUYAL P, KIM D G, et al. Sulfur speciation in petroleum: Atmospheric pressure photoionization or chemical derivatization and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry[J].Energy Fuels, 2007, 21(5): 2869-2874.

[9] JEREMIAH M P, ISABELLE M, RYAN P R, et al. Stepwise structural characterization of asphaltenes during deep hydroconversion process determined by atmospheric pressure photoionization (APPI) Fourier Transform Ion Cyclotron Resonance (FT-ICR) mass spectrometry[J].Energy Fuels, 2010, 24(4): 2257-2265.

[10] 刁瑞. 渣油原料和加氫產品的分子水平認識[D].北京: 石油化工科學研究院, 2013.

[11] 董凱, 邵志才, 劉濤, 等. 儀長渣油加氫處理反應規律的研究Ⅰ. 儀長渣油性質特點及加氫反應特性[J].石油煉制與化工, 2015, 46(1): 1-5. (DONG Kai, SHAO Zhicai, LIU Tao, et al. Study on hydrotreating reactivity of Yichang residue Ⅰ Characterization and hydrotreating reactivity of Yichang residue[J].Petroleum Processing and Petrochemicals, 2015, 46(1): 1-5.)

[12] 李大東. 加氫脫氮的工藝及催化劑[J].石油煉制, 1983, 12: 14-23. (LI Dadong. Hydro-denitrification technology and catalyst[J].Petroleum Refining, 1983, 12: 14-23.)

[13] 邵志才, 高曉冬, 李皓光, 等. 氮化物對柴油深度和超深度加氫脫硫的影響: I. 氮化物含量的影響[J].石油學報(石油加工), 2006, 22(4): 12-17. (SHAO Zhicai, GAO Xiaodong, LI Haoguang, et al. Influence of nitrogen containing compounds on diesel ultra-deep hydrodesulfurization I Effect of nitrogen compounds content[J].Acta Petrolei Sinica(Petroleum Processing Section), 2006, 22(4): 12-17.)

[14] 崔瑞利, 趙愉生, 薛鵬, 等. 固定床渣油加氫處理裝置催化劑積炭分析及表征[J].石油煉制與化工, 2012, 43(1): 45-48. (CUI Ruili, ZHAO Yusheng, XUE Peng, et al. Characterization of coke on spent catalysts from fixed-bed hydrotreating unit processing residue feed[J].Petroleum Processing and Petrochemicals, 2012, 43(1): 45-48.)

主站蜘蛛池模板: …亚洲 欧洲 另类 春色| 凹凸国产熟女精品视频| 久久五月视频| 亚洲无码视频图片| 国产亚洲男人的天堂在线观看| 无码aⅴ精品一区二区三区| 老司国产精品视频| 国产亚洲精品无码专| 亚洲欧美一区二区三区图片| 国产男女XX00免费观看| 18禁影院亚洲专区| 久久男人资源站| 欧美色视频网站| 狠狠色丁香婷婷| 农村乱人伦一区二区| 在线无码私拍| 亚洲综合天堂网| 免费可以看的无遮挡av无码| 国产精品女人呻吟在线观看| 69综合网| 精品欧美一区二区三区久久久| 国产精品久久国产精麻豆99网站| 国产精品手机视频一区二区| 国产极品美女在线播放| 国产亚洲精久久久久久久91| 国产玖玖玖精品视频| 精品無碼一區在線觀看 | 免费A∨中文乱码专区| 精久久久久无码区中文字幕| 国产玖玖视频| 色亚洲激情综合精品无码视频| 日韩欧美91| 久久久久夜色精品波多野结衣| 亚洲国产中文在线二区三区免| 在线免费观看AV| 成人综合在线观看| 国产福利影院在线观看| 亚洲美女视频一区| 国产女人在线视频| 青青青国产视频| m男亚洲一区中文字幕| 成人国产精品一级毛片天堂| 色综合国产| 欧美成人综合在线| 干中文字幕| 久久人体视频| 亚洲伊人久久精品影院| 亚洲综合色婷婷中文字幕| 波多野结衣中文字幕一区二区| 国产 在线视频无码| 日韩在线欧美在线| 亚洲色偷偷偷鲁综合| 蜜臀av性久久久久蜜臀aⅴ麻豆| 高清无码手机在线观看| 国产网友愉拍精品| 内射人妻无套中出无码| 91最新精品视频发布页| 国产肉感大码AV无码| 亚洲视频欧美不卡| 色综合成人| 在线观看精品国产入口| 尤物国产在线| 真实国产乱子伦视频| 欧美精品1区2区| 国语少妇高潮| 日韩第九页| 国产乱论视频| 国产91丝袜在线播放动漫 | 国产国语一级毛片在线视频| 日韩第九页| 国产亚洲欧美在线人成aaaa| 日a本亚洲中文在线观看| 国产成人永久免费视频| 亚洲欧美日韩视频一区| 性网站在线观看| 欧美亚洲国产精品第一页| 国产成人久久综合一区| 成人一级黄色毛片| 国产精品久久久久久久久| 亚洲天堂网在线播放| 日韩福利在线视频| 国产精品毛片一区|