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

納米流體中CO2水合物生成特性實驗研究

2015-12-22 08:18:24劉妮張亞楠柳秀婷由龍濤
制冷學報 2015年2期
關鍵詞:實驗

劉妮 張亞楠 柳秀婷 由龍濤

(上海理工大學能源與動力工程學院 上海 200093)

納米流體中CO2水合物生成特性實驗研究

劉妮 張亞楠 柳秀婷 由龍濤

(上海理工大學能源與動力工程學院 上海 200093)

在自行設計的小型氣體水合物反應裝置上進行了納米流體中CO2水合物生成特性的實驗研究,探討了納米粒子的種類、粒徑和質量分數對CO2水合物生成特性的影響。研究發現,與純水相比,納米粒子CuO和SiO2增加了CO2耗氣量,但延長了氣體水合物生成的誘導時間。金屬納米粒子Cu和金屬氧化物納米粒子Al2O3對CO2水合物生成的誘導時間和耗氣量有明顯改善。對不同粒徑的Al2O3納米粒子對氣體水合物生成特性的研究發現,30 nm的Al2O3納米流體對水合物生成特性影響最大。與純水相比,0.1%-30 nm-Al2O3納米流體中水合物生成的誘導時間縮短了76.9%,耗氣量增加了23.2%。CO2水合物耗氣量隨著Cu粒子質量分數的增加先增加后減少。最后進行了納米粒子對CO2水合物生成特性的影響的理論分析。

蓄冷;CO2水合物;納米流體;誘導時間;耗氣量

隨著氣體水合物這種新型蓄冷工質的研究發展,CO2水合物因其環保與高儲能密度等特點引起研究者的廣泛關注[1]。氣體水合物是一種特殊的包絡化合物,主體分子即水分子間以氫鍵相互結合形成的籠型孔穴將客體分子包絡在其中所形成的非化學計量的化合物,客體分子和主體分子之間通過范德華力連接,形成了具有熱力學穩定性的氣體水合物。

在蓄冷空調的使用中,CO2氣體水合物的生成緩慢是目前面臨的一個重要問題[2],因此對CO2氣體水合物生成過程的熱力學平衡與生成動力學及其微觀反應機理進行深入的研究非常重要。由于氣體水合物的生成過程涉及復雜的氣液固多相傳熱傳質過程,氣相與液相之間難于充分混合,相界面處擴散速度很慢,導致水合反應的誘導時間長、過冷度大、水合物生成速度緩慢。因此,如何提高水合物的生成速率和儲氣密度、縮短誘導時間是蓄冷空調實現高效節能運行,并對其開展技術應用的關鍵,而其核心問題則是要解決氣液固多相之間的傳熱傳質阻力。研究學者常用的方法有攪拌[3]、氣泡[4]、噴霧[5]和加入添加劑[6-7]等,雖然這些方法的作用突出,但存在能耗大、不利于實際應用等特點。因此,尋找新的提高水合物生成特性的方法是學者的研究方向。

1995年,Choi SU S[8]首次提出了“納米流體”的概念,即將1~100 nm的金屬或者非金屬粒子懸浮在基液中形成的穩定懸浮液。近十幾年來,隨著納米技術的發展,納米流體已成為一種最有吸引力的強化傳熱傳質介質。一些研究者提出了納米流體改善水合物生成過程中的誘導時間、生成量和生成速率的設想,并進行了相關實驗研究。Sung-Seek Park等[9]研究了甲烷水合物在多壁碳納米管溶液中的生成特性,實驗發現在質量分數為0.004%的多壁碳納米管溶液中的甲烷水合物生成量最多,此時甲烷的耗氣量是在純水中的4.5倍。李金平[10]研究了制冷劑水合物在TiO2、Cu和Ag等納米流體中的靜態生成過程。實驗結果表明納米粒子的加入能夠加速氣體水合物的結晶和生長。彭浩等[11]研究了THF水合物在不同納米碳管質量分數(0.1%~10%)、不同溫度(258.15~270.15 K)時的熱導率。實驗結果顯示,含納米碳管的THF水合物的熱導率高于純THF水合物的熱導率,隨著納米碳管質量分數和溫度的增大,THF水合物的熱導率大幅升高。Krishnamurthy S等[12]的實驗研究工作表明懸浮的納米粒子明顯地增強了兩元混合液體中的擴散過程,表明納米流體可以強化傳質過程。Abolfazl Mohammadi等[13]實驗研究結果表明,SDS和銀納米顆粒并沒有顯著降低誘導時間和增加CO2水合物的存儲量,而SDS和銀納米顆粒的混合物顯著增加了CO2的存儲量。Mostafa Keshavarz Moraveji等[14]的研究表明納米粒子積聚在甲烷和水的接觸界面可以降低表面和界面張力。因此,與純水相比,添加納米粒子的溶液增強了甲烷的溶解度。

縱觀各類文獻,目前關于納米流體對二氧化碳水合物生成的熱力學及其對二氧化碳水合物生成過程的影響和作用機理的研究國內外鮮有報道。本文在小型水合物反應裝置上,以CO2氣體水合物為研究對象,探討了納米流體中CO2氣體水合物生成過程中的誘導時間、耗氣量和生成速率的變化,并對其作用機理進行了分析。

1 實驗部分

1·1 實驗材料

實驗所用CO2氣體的純度為99.99%,由上海偉創氣體有限公司提供。分散劑采用十二烷基苯硫酸鈉,即 SDBS,分子式為 C18H29NaO3S,其純度大于92%,由上海凌峰化學試劑有限公司提供。實驗選用四種納米粒子Cu、CuO、Al2O3和SiO2的純度均為99.99%,來自上海超威納米科技有限公司。實驗用水為自制的一次蒸餾水。

1·2 實驗裝置

實驗裝置如圖1所示。核心部件高壓反應釜,由316 L不銹鋼材料加工而成,有效容積為680 mL,耐壓20 MPa。為了增加氣體和溶液的接觸面積,反應釜頂部安裝了磁力攪拌器,其轉速在0~1000 r/min內可調。反應釜所需要的溫度由兩臺獨立的精密低溫恒溫槽提供,精度為0.01℃,所用冷媒為濃度40%~60%的乙二醇溶液。反應釜內溫度和壓力通過Pt100鉑電阻溫度傳感器和NS-I/10 MPa擴散硅式壓力傳感器來監測。實驗采用Agilent 34970A數據采集儀記錄數據。

圖1 實驗裝置系統示意圖Fig·1 SchematiCdiagramof experimental system

1·3 實驗過程

整個實驗在恒容條件下進行。實驗采用兩步法制備納米流體。先在蒸餾水中添加0.8 mg/mL 的SDBS作為分散劑,再加入相應納米粒子配成納米流體。對配置好的溶液進行30 min以上的超聲振蕩。然后,用粒徑測試儀和靜置法對納米流體的分散穩定性進行表征。實驗過程中先用蒸餾水清洗反應釜2~3次并擦拭干凈,再加入配置好的200 mL的納米流體并密封。將反應釜放入20℃的恒溫水浴中,用真空泵對系統抽真空,并檢查系統的氣密性。然后充入CO2氣體,開啟磁力攪拌器(轉速200 r/min)。當釜內溫度恒定且壓力穩定在預定壓力之后關閉供氣閥,將反應釜快速置入預先設定反應溫度(低于相平衡溫度)的恒溫槽中,降溫冷卻進行CO2水合物的生成反應,并記錄數據。反應結束后,保存實驗數據,排盡釜內剩余氣體。實驗過程中的反應溫度設為2℃,初始壓力為3 MPa左右。

2 結果與分析

2·1 不同種類的納米粒子對CO2水合物生成特性的影響

圖2是分別添加納米粒子 Cu、CuO、Al2O3和SiO2的溶液中典型的CO2水合物生成曲線,粒子濃度均為0.1%。圖2(a)中,隨著時間的變化,溫度突升,表示水合物開始大量生成。由于添加的納米粒子的導熱系數不同,曲線斜率存在差別,表示溶液降溫速率不同。水合物的平均生成速率用一段時間內反應器中CO2氣體的壓降來表示[15]。從圖2(b)中看出,一開始不同粒子的納米流體中CO2氣體的溶解速率基本一致。隨著CO2水合物的大量生成,不同的納米流體體系表現出異樣的壓降特性。Al2O3納米流體中CO2水合物的平均生成速率最大。

圖2 納米流體中CO2水合物生成曲線Fig·2 CO2hydrate formation in different nanofluids

在相同的實驗條件下,每種納米流體體系中CO2水合物的生成實驗重復三次。表1是不同納米流體體系中CO2水合物生成的誘導時間和耗氣量的實驗結果。本文中誘導時間定義為實驗中水合物開始出現明顯結晶的時間和水合物相平衡曲線與系統冷卻曲線相交點的時間之差[16]。與純水中CO2水合物的生成特性相比,Cu和Al2O3納米流體體系均縮短了CO2水合物的誘導時間,增加了CO2耗氣量。Al2O3納米流體中CO2水合物生成的誘導時間最短,與純水相比,縮短了76.9%。Cu納米流體中CO2耗氣量最大,與純水相比,提高了33.6%。分析認為,一方面納米添加劑降低了接觸界面的表面和界面張力,提高了CO2在水中的溶解速率,同時納米粒子的存在提供了大量的成核點,增加了氣液有效接觸面積,增強了溶液中的傳質過程。另一方面添加表面活性劑SDBS的Cu和Al2O3納米流體體系均有較高的導熱率[17-18],能夠增強換熱,提供CO2水合物生成所需要的有利條件。于是,分別進行了不同粒徑下的Al2O3和不同質量分數下的Cu納米流體體系中CO2水合物生成特性的實驗研究。

表1 CO2水合物生成實驗結果Tab·1 Experimental results of CO2hydrate formation

圖3為實驗反應結束后所拍攝的水合物照片。由圖可知,水合物圍繞金屬壁面生成。實驗中,反應釜從20℃的恒溫水浴進入2℃的恒溫水浴,金屬壁面和溶液之間形成了一個較大的溫度差,滿足了水合物成核需要的過冷度,使得CO2水合物快速生成。同時添加的納米粒子具有較大的表面自由能,為CO2水合物生成提供了大量的成核點。李金平[10]研究了添加鐵絲后納米流體中制冷劑氣體水合物生成現象,

也得出了相同的水合物生成位置及納米粒子有利于水合物生成的結論。

圖3 反應釜內生成的CO2水合物照片Fig·3 Photos of CO2hydrate in cell

2·2 不同粒徑Al2O3納米粒子對CO2水合物生成特性的影響

實驗中使用質量分數0.1%和1.0%,粒徑為10 nm、30 nm和50 nm的Al2O3納米流體進行研究。圖4是質量分數為0.1%情況下,不同粒徑Al2O3納米流體中CO2水合物生成過程中典型的p-t曲線。由圖可知,30 nm的Al2O3納米流體中CO2水合物生成所需誘導時間最短。表2是0.1%和1.0%質量分數下,不同粒徑的Al2O3納米流體中CO2水合物生成特性的實驗結果。從表中看出,粒徑為30 nm的Al2O3納米流體中CO2耗氣量最大,粒徑的變化對水合物生成的誘導時間和耗氣量影響顯著。一方面由于納米粒子粒徑越小,微運動強度越大,在溶液中的移動越頻繁,使得溶液內部能量快速傳遞,同時增加CO2氣體在水中的溶解速率。另一方面是溶液中熱傳遞過程發生液體和粒子界面,且小粒徑的納米粒子具有較大的比表面積,因此小粒徑的納米粒子更有利于水合物的生成。但10 nm的納米流體中水合物生成所需的誘導時間和耗氣量均比30 nm的納米流體低,分析認為粒徑越小,單位體積內粒子數目越多,粒子間距越小,相互作用力越大,使得黏度增加,溶液傳熱傳質速率降低。所以,存在一個最佳粒徑能有效改善納米流體中水合物的生成特性。

2·3 不同質量分數Cu納米粒子對CO2水合物的生成特性的影響

實驗使用的Cu納米粒子粒徑為30 nm,質量分數分別為0.1%、0.5%、1.0%和1.5%。表3是含不同質量分數的Cu納米流體對CO2水合物生成特性的影響。從表中看出,質量分數為0.5%的Cu納米流體中CO2水合物誘導時間最短,與純水相比,縮短了54.7%。CO2水合物耗氣量隨著Cu納米質量分數的增加先增加后降低,如圖5所示。1%質量分數的Cu納米流體中水合物生成的耗氣量最多,與純水相比,增加了37.6%。溶液中添加一定量的納米粒子,可以增加納米流體的導熱率,能夠及時降低溶液溫度,增加水合物成核驅動力。但隨著納米粒子質量分數的增加,溶液中納米粒子容易發生團聚沉淀,降低流體內部的能量傳遞,即降低溶液導熱率,從而不利于CO2水合物的生成。另一方面由于納米流體黏度發生了變化,導致隨著粒子質量分數的增加水合物生成量減少,如Liu Yang等[19]的研究。

圖4 Al2O3納米粒徑對誘導時間的影響Fig·4 Effects of Al2O3nanoparticle size on induction time

表2 不同的Al2O3納米粒徑下水合物生成實驗結果Tab·2 Experimental results of hydrate formation in differentAl2O3nanoparticle size

3 結論

本文對納米流體體系中CO2水合物的生成特性進行了實驗研究,探討了納米粒子的類別、粒徑和質量分數對CO2水合物生成特性的影響。研究發現,與純水相比,Cu和Al2O3納米流體中CO2水合物生成的誘導時間和耗氣量均有明顯改善。與純水相比,CuO和SiO2納米流體中CO2水合物生成的耗氣量增加,但誘導時間沒有縮短。納米粒子粒徑對CO2水合物生成特性影響顯著。存在最佳的粒徑值使得CO2水合物生成特性最佳。在30 nm的Al2O3納米流體中,CO2水合物生成所需的誘導時間最短,耗氣量最多。CO2水合物生成的耗氣量隨著Cu納米粒子質量分數的增加先增加后降低。Cu納米粒子質量分數為1%時,納米流體中水合物耗氣量最多,與純水相比,增加了37.6%。不同質量分數下的Cu納米流體均縮短了CO2水合物生成的誘導時間。綜合納米粒子對CO2水合物生成特性的影響,分析認為一是納米粒子可以增強溶液的傳熱傳質過程,二是納米粒子會增加溶液的黏度,當兩者處于耦合最佳狀態時,納米粒子對CO2水合物的生成特性的影響最佳。

表3 不同質量分數的Cu納米流體中CO2水合物生成實驗結果Tab·3 Experimental results of CO2hydrate formation in differentmass fraction of Cu nanofluid

圖5 Cu質量分數對CO2水合物耗氣量的影響Fig·5 Effects ofmass fraction of Cu nanoparticle on CO2hydrate consumption

[1] 魏晶晶,謝應明,劉道平.水合物在蓄冷及制冷(熱泵)領域的應用[J].制冷學報,2009,30(6):36-43.(Wei Jingjing,Xie Yingming,Liu Daoping.Application of Hydrate in Thermal Storage and Refrigeration(Heat Pump) System[J].Journal of Refrigeration,2009,30(6):36-43.)

[2] 謝振興,謝應明,舒歡,等.蓄冷用CO2水合物漿的強化制備研究進展[J].制冷學報,2013,34(5):82-89.(Xie Zhenxing,Xie Yingming,Shu Huan,et al.Research Progress of CO2Hydrate Slurry Formation Enhancement in Cool Storage Application[J].Journal of Refrigeration,2013,34(5):82-89.)

[3] 劉妮,劉道平,謝應明.水合物法高效儲存二氧化碳氣體的實驗研究[J].中國電機工程學報,2009,29(14):36-40.(Liu Ni,Liu Daoping,Xie Yingming.Experimental Study on CO2Storage by Hydrate Crystallization[J].Proceedings of the CSEE,2009,29(14):36-40.)

[4] Brij B Maini,PR Bishnoi.Experimental Investigation of Hydrate Formation Behavior ofANatural Gas Bubble inASimulated DeePSea Environment[J].Chemical Engineering Science,1981,36(1):183-189.

[5] 李剛,劉道平,肖楊,等.CO2水合物噴霧合成的生長特性實驗研究[J].上海理工大學學報,2009,31(3):213-217.(LiGang,Liu Daoping,Xiao Yang,etal.Experimental Research on Growth Characteristics of CO2Hydrate Composed by Water Spray[J].J.University of Shanghai for Science and Technology,2009,31(3):213-217.)

[6] Asheesh Kumar,Tushar Sakpal,Praveen Linga,et al.Influence of ContactMediumand Surfactants on Carbon Dioxide Clathrate Hydrate Kinetics[J].Fuel,2013,105:664-671.

[7] Mingjun Yang,Yongchen Song,Weiguo Liu,etal.Effects of Additive Mixtures(THF/SDS)on Carbon Dioxide Hydrate Formation and Dissociation in Porous Media[J].Chemical Engineering Science,2013,90:69-76.

[8] Choi SU S.Enhancing Thermal Conductivity of Fuids with Nanoparticles in Developments and Applications of non-Newtonian[J].ASME FED,1995,231(66):99-103.

[9] Sung-Seek Park,Eoung-Jin An,Sang-Baek Lee,et al.Characteristics of Methane Hydrate Formation in Carbon Nanofluids[J].Journal of Industrial and Engineering Chemistry,2012,18(1):443-448.

[10]李金平.納米流體制冷劑氣體水合物生成動力學和相平衡特性研究[D].上海:上海交通大學,2005.

[11]彭浩,樊栓獅,孫始財.納米碳管強化水合物導熱的研究[J].制冷學報,2006,28(2):39-43.(Peng Hao,Fan Shuanshi,Sun Shicai.Research on Hydrate Thermal Conductivity Enhanced by Carbon Nanotube[J].Journal of Refrigeration,2006,28(2):39-43.)

[12]Krishnamurthy S,Bhattacharya P,Phelan PE.Enhanced Mass Transport in Nanofluids[J].Nano Letters,2006,63: 419-423.

[13]Abolfazl Mohammadi,Mehrdad Manteghian,Ali Haghtalab,et al.KinetiCStudy of Carbon Dioxide Hydrate Formation in Presence of Silver Nanoparticles and SDS[J].Chemical Engineering Journal,2014,237:387-395.

[14]Mostafa Keshavarz Moraveji,Milad Golkaram,Reza Davarnejad,et al.Biosynthesized Silver Nanofluid Effect on Methane Dissolution in Water[J].Journal ofMolecular Liquids,2013,184:1-3.

[15]劉妮,張國昌,R E羅杰斯.二氧化碳氣體水合物生成特性的實驗研究[J].上海理工大學學報,2007,29(4): 405-408.(Liu Ni,Zhang Guochang,R E Rogers.Experimental Study of CO2Gas Hydrates Formation[J].J.University of Shanghai for Science and Technology,2007,29 (4):405-408.)

[16]Liu Ni,Xuan Xiaobo,Liu Daoping,et al.CO2Reduction by Clatchrate Hydrate Crystallization[J].Environmental Engineering and Management Journal,2011,10(6):837-841.

[17]Dongsheng Zhu,Xinfang Li,Nan Wang,et al.Dispersion Behavior and Thermal Conductivity Characteristics of Al2O3-H2O Nanofluids[J].Current Applied Physics,2009, 9(1):131-139.

[18]Xinfang Li,Dongsheng Zhu,Xianju Wang,et al.Thermal Conductivity Enhancement Dependent pH and Chemical Surfactant for Cu-H2O Nanofluids[J].Thermochimica Acta,2008,469(1/2):98-103.

[19]Liu Yang,Kai Du,Xiaofeng Niu,et al.Experimental Study on Enhancement of Ammonia-water Falling FilmAbsorption by Adding Nanoparticles[J].International Journal of Refrigeration,2011,34(3):640-647.

About the corresponding author

Liu Ni,female,Ph.D./associate professor,School of Energy and Power Engineering,University of Shanghai for Science and Technology,+86 21-55271619,E-mail:liu_ni@163.com.Research fields:control of CO2emission and cool storage technology by gas hydrate.

Experimental Study on Characteristics of CO2Hydrate Formation in Nanofluids

Liu Ni Zhang Yanan Liu Xiuting You Longtao

(School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai,200093,China)

The characteristics of CO2hydrate formation in nanofluids were investigated experimentally inAsmall scale reaction system. The effects of type,particle size and mass fraction of nanoparticles on CO2hydrate formation characteristics were stuided.Compared with pure water,CuO and SiO2nanoparticles increase gas consumption,but extend induction time of gas hydrate formation.Results show that metal nanoparticle Cu andmetaloxide nanoparticle Al2O3both haveAsignificant impacton the induction time and gas consumption of CO2hydrate formation.Al2O3showsAremarkable effcet on characteristics of gas hydrate formation withAparticle size of 30 nm.Compared with purewater,the induction time during hydrate formation in 0.1%-30 nm-Al2O3nanofluid is reduced by 76.9%,and gas consumption is increased by 23.2%.With the increase of Cu nanoparticle mass fraction,CO2gas consumption is first increased and then decreased.Themechanismof the effect of nanoparticles on characteristics of CO2hydrate formation was discussed.

cool storage;carbon dioxide hydrate;nanofluids;induction time;gas consumption

TB64;TB383

A

0253-4339(2015)02-0041-06

10.3969/j.issn.0253-4339.2015.02.041

簡介

劉妮,女,博士,副教授,上海理工大學能源與動力工程學院,(021)55271619,E-mail:liu_ni@163.com。研究方向:二氧化碳氣體排放控制及水合物蓄能技術應用研究。

國家自然科學基金(50706028)資助項目。(The projectwas supported by the National Natural Science Foundation of China(No. 50706028).)

2014年8月8日

猜你喜歡
實驗
我做了一項小實驗
記住“三個字”,寫好小實驗
我做了一項小實驗
我做了一項小實驗
記一次有趣的實驗
有趣的實驗
小主人報(2022年4期)2022-08-09 08:52:06
微型實驗里看“燃燒”
做個怪怪長實驗
NO與NO2相互轉化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
主站蜘蛛池模板: www.国产福利| 99性视频| 亚洲欧美日韩色图| 久久99热66这里只有精品一| 国产女人喷水视频| 欧美在线精品怡红院| 亚洲一级无毛片无码在线免费视频| 免费在线一区| 国产视频一二三区| 日韩欧美亚洲国产成人综合| 国产乱人视频免费观看| 欧美va亚洲va香蕉在线| 亚洲毛片网站| 日韩高清一区 | 黄色网址免费在线| 久久这里只有精品23| 欧美区在线播放| 无码高清专区| 国产免费好大好硬视频| 国产成人精品亚洲日本对白优播| 亚洲中文字幕无码mv| 91原创视频在线| 色噜噜中文网| 亚洲精品自拍区在线观看| 毛片久久网站小视频| 国产主播在线观看| 一级片一区| 天堂网亚洲系列亚洲系列| 国产一区二区影院| 国产喷水视频| 国产精品成人久久| 亚洲最新在线| 99热这里只有成人精品国产| 久久窝窝国产精品午夜看片| 依依成人精品无v国产| 2021精品国产自在现线看| 国产本道久久一区二区三区| 国产视频一二三区| 国产成人精品高清在线| 色婷婷国产精品视频| 久久黄色小视频| 國產尤物AV尤物在線觀看| 亚洲国产精品无码AV| 五月天福利视频| 国产中文一区二区苍井空| 激情国产精品一区| 亚洲精品卡2卡3卡4卡5卡区| 国产成人h在线观看网站站| 国产中文一区二区苍井空| 欧美日韩中文国产| 亚洲黄色视频在线观看一区| 国产精品微拍| 中文无码毛片又爽又刺激| 国产自产视频一区二区三区| 美女国内精品自产拍在线播放 | 国产欧美日韩精品第二区| 精品久久人人爽人人玩人人妻| 538国产视频| 日韩无码黄色网站| 久久综合AV免费观看| 中文字幕在线免费看| 色偷偷一区二区三区| 亚洲中文字幕97久久精品少妇| 久久无码高潮喷水| 精品国产免费人成在线观看| 精品欧美一区二区三区在线| 亚洲A∨无码精品午夜在线观看| 亚洲人成成无码网WWW| 欧美日韩国产精品va| 国产SUV精品一区二区6| 狠狠干综合| 亚洲第一天堂无码专区| 欧美精品成人一区二区在线观看| 国产杨幂丝袜av在线播放| 999精品视频在线| vvvv98国产成人综合青青| 国产在线观看一区二区三区| a级毛片毛片免费观看久潮| 亚洲伊人天堂| 日本道中文字幕久久一区| 色窝窝免费一区二区三区| 免费国产高清精品一区在线|