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

三元水鹽體系冰點和共晶點的測定、相圖表達及計算

2017-09-16 06:04:35劉福云黃雪莉黃文婷王雪瑩潘毅
化工學報 2017年9期
關鍵詞:體系實驗

劉福云,黃雪莉,黃文婷,王雪瑩,潘毅

?

三元水鹽體系冰點和共晶點的測定、相圖表達及計算

劉福云,黃雪莉,黃文婷,王雪瑩,潘毅

(新疆煤炭潔凈轉化與化工過程重點實驗室,新疆大學化學化工學院,新疆烏魯木齊 830046)

針對Na+//Cl-,NO-3-H2O、Na+//Cl-,SO42--H2O和Na+,K+//Cl--H2O 3個三元水鹽體系,測定了不同組成時溶液的冰點、共晶點以及降溫過程的析鹽規律,繪制出三元體系溫度-組成立體相圖。結果表明:三元體系的冰點、共晶點、相關系與組成可用三棱柱表達,3個三元水鹽體系在共晶點溫度下均無復鹽存在,各含有1個兩鹽與冰的共晶點,3條單鹽與冰或鹽鹽的共晶線,1個冰和兩個單鹽結晶區;用三元體系溫度-組成立體相圖能更直觀地表達三元水鹽體系的結晶區、冰點、共晶點和溶液組成的關系,可以判斷溶液降溫過程中的鹽或冰的析出類型、溫度;同時建立了冰點經驗公式,應用經驗公式和過量吉布斯自由能模型預測三元體系的冰點和共晶點,計算值和實驗值吻合良好。

水鹽體系;共晶點;冰點;相圖;計算

引 言

鹽湖資源的開發利用主要是通過自然蒸發或強制蒸發、冷卻、結晶、過濾等傳統方法[1-2]加以分離,雖然利用了太陽能和風能,降低了能耗和成本,但也損失了大量的水,對于干旱缺水的鹽湖資源地區,這是一個嚴峻的問題,制約著鹽湖化工的發展。因此,開發節水新技術十分必要。

除了利用夏季高溫熱能蒸發鹵水,也可以考慮利用冬季冷能,在不損失或少損失水的條件下,實現鹵水的結晶和分離,但目前低溫下加工鹽湖資源[3-14]的技術很少。由于相平衡規律是鹽湖資源開發的理論基礎,因此,研究低溫條件下水鹽體系的冰點、共晶點以及鹽析規律,對利用冬季冷能開發低溫鹽湖化工工藝以及一些含鹽廢水的低溫處理技術[15-17]具有重要的理論指導意義。

對于大部分單鹽溶液,其冰點和共晶點數據都有報道[18-20],但對于三元及以上水鹽體系報道很少。Oakes等[21]研究了NaCl-CaCl2-H2O體系的冰點;Gibbard等[22]研究了NaCl-MgCl2-H2O體系的冰點;Desnoyers等[23]研究了NaCl-EtOH-H2O體系的冰點;王雪瑩等[24]研究了K+//Cl-,SO42-,NO-3-H2O和Na+,K+//SO42-,NO-3-H2O兩個四元水鹽體系及其相關的6個三元子體系的冰點、共晶點,提出了用簡單三棱柱圖形來表達三元水鹽體系冰點、共晶點與組成的關系,但未對其進行明確的相區劃分和應用研究。

一般來說,海水和絕大多數鹽湖體系均屬于四元以上水鹽體系,組分多,相關系復雜,如需對其研究,必須首先進行相關子體系的研究。Na+//Cl-,NO-3-H2O、Na+//Cl-,SO42--H2O和Na+,K+// Cl--H2O 3個三元體系為常見的基礎子體系,本文研究了這3個三元體系的冰點、共晶點和鹽析規律,為研究四元、五元或更復雜的鹽湖鹵水體系低溫相平衡打下基礎,也為實現利用冬季冷能處理廢鹽水、開發鹽湖資源提供理論支持。

1 實驗部分

1.1 實驗試劑和實驗儀器

試劑:NaCl、KCl、NaNO3、Na2SO4均為A.R. 級;實驗用水均為蒸餾水。

儀器:DHJF-4010型低溫(恒溫)攪拌反應浴(鄭州長城科工貿有限公司,使用范圍-40~100℃,控溫精度±0.2℃);JDW-3F型精密電子溫差測量儀(北京亞美福科技有限公司,使用范圍-50~100℃,精度±0.001℃);XTL-1型體視顯微鏡(南京江南光電股份有限公司);X 射線衍射儀(荷蘭Philip X公司)。

1.2 實驗方法

具體實驗裝置及操作方法詳見文獻[24-25],本文只作簡要說明。

1.2.1 三元體系共晶點溫度的測定 以三元體系Na+//Cl-,NO--3-H2O為例,首先,根據常溫下該體系的共飽點組成配制300 ml溶液,取約10 ml溶液于直徑2 cm的試管中,外加直徑為2.5 cm試管為套管,置于低溫恒溫裝置中降溫,外加電動增力攪拌器進行攪拌,記錄時間、溫度,作出時間-溫度關系圖,根據溫度波動點確定體系冰點、共晶點。

1.2.2 三元體系共晶點溶液組成的測定 將配制的剩余溶液放入平底燒瓶,置于低溫恒溫裝置,溫度設置為略低于已測出的共晶點或冰點溫度,攪拌降溫。當溶液到達共晶點溫度后,觀察到較多的冰析出時,停止攪拌靜止約10 min,冰浮于液面,鹽沉降于底部,取液相樣分析,確定共晶點溶液組成,同時分離固相進行鑒定。

1.3 化學分析方法

Cl-:硝酸銀容量法。K+:四苯硼鈉重量法。SO42-:比濁法。NO-3:重鉻酸鉀氧化法。Na+:差減法。固相鑒定采用體式顯微鏡或X射線晶體衍射綜合確定。

2 實驗結果與討論

2.1 三元體系研究結果

上述3個三元體系的共晶點的液相組成和溫度數據列于表1中。

由表1數據可知,3個三元體系在低溫下均為簡單體系,無復鹽存在;均含有1個共晶點,對應的共晶點溫度均低于所對應的兩種單鹽溶液的共晶點溫度。

表1 三元體系共晶點溫度及液相組成

Note: S10—Na2SO4·10H2O; Sy—KCl; NN—NaNO3; Hy—NaCl·2H2O.

2.2 三元體系冰點立體相圖

為了能直觀地表達三元水鹽體系的結晶區、冰點、共晶點和溶液組成的關系,采用三棱柱表達。底面等邊三角形為組成面,3個頂點為純組分;棱代表溫度;3個立面分別代表3個二元水鹽體系;三棱柱內部為三元水鹽體系。

以Na+//Cl-,NO-3-H2O三元體系為例進行具體說明,如圖1所示(為清晰起見,立體圖局部未嚴格按照比例繪制)。

由圖1可知:

(1)底面等邊三角形為NaCl、NaNO3、H2O的組成區,3個頂點分別代表NaCl、NaNO3、H2O 3個純組分。

(2)3個立面′′、′′和′′分別代表NaCl-H2O、NaNO3-H2O和NaCl-NaNO33個二元體系,其中NaCl-H2O、NaNO3-H2O是水鹽體系,NaCl-NaNO3是鹽鹽體系(′、′分別是NaCl、NaNO3熔點,分別為801℃、306.8℃[26]。NaCl-NaNO3鹽鹽體系的共熔點,一般低于對應的單鹽熔點,具體數據未見報道,本文因只研究水鹽體系,故未對其進行測定,設為點3,為示意點;由于3處的氯化鈉是以NaCl形式存在,而處的氯化鈉是以NaCl·2H2O形式存在,則3線上必有一個轉變點,是NaCl和NaCl·2H2O的共晶點,因測定困難,未進行測定,設為點5,為示意點)。

(3)3個空間曲面,代表了3個單固相的飽和溶液面。其中空間曲面′4153′是NaCl和NaCl·2H2O的飽和溶液面;空間曲面′253′是NaNO3的飽和溶液面;空間曲面′12′是冰的飽和溶液面。3個空間曲面上方的空間區域表示不飽和溶液區。

(4)3條空間曲線,代表3組雙固相的共晶溶液線。其中1是NaCl·2H2O和冰的共晶溶液線,2是NaNO3和冰的共晶溶液線,53是NaCl、NaCl·2H2O和NaNO3的共晶溶液線,1、2、3分別是二元體系中相應的共晶點。

(5)點是NaCl·2H2O、NaNO3、冰3個固相的共晶點,實驗測得該共晶點的冰點溫度為-26.91℃。

(6)由圖2可知:′3′、′4′、′354′、54、53、構成的多面體為NaCl的結晶區;由45、5415、41、、、5、1構成的多面體為NaCl·2H2O的結晶區;由′3′、′2′、′253′、35、2、構成的多面體為NaNO3的結晶區;由′2′、′1′、′12′、1、2、構成的多面體為冰的結晶區。每種鹽或冰的結晶區都有1個斷層面,落在底面的投影為線、、、′,這4條線左右的點雖然溶液組成接近,但鹽析出溫度差別比較大。

(7)由圖3可知:、、′、′、′、′、′′′、′′構成的多面體為NaCl和NaCl·2H2O的結晶區;、、′、′′、′構成的多面體為NaCl·2H2O和冰的結晶區;、、、、構成的多面體為NaNO3和冰的結晶區;、、、、構成的多面體為NaCl、NaCl·2H2O和NaNO3的結晶區。

(8)5個三角面、、、′、下方為全固相區域,液相完全消失。

2.3 三元體系變溫立體相圖的應用

根據繪制出的三元變溫立體相圖,可以確定任意組成的溶液在降溫過程中的冰鹽析出順序、冰點近似溫度、單鹽結晶溫度、共晶點溫度和液相組成。

以圖4中的點溶液降溫過程為例進行說明。

點位于不飽和溶液區,過點向下做垂線——降溫線,穿過NaNO3的結晶區,在底面中投影為4點。

當溶液從點降溫至1點時,NaNO3析出,1點即為NaNO3的結晶點,對應的縱坐標即為NaNO3的結晶溫度;由1點繼續降溫至2點時,冰也析出,對應的縱坐標即為NaNO3和冰的共晶點溫度;繼續降溫,NaNO3和冰不斷析出,液相點由點移向點,到達點后,NaCl·2H2O也飽和析出,總組分點由2降至3,液相全部消失,進入全固相區域。

1點的確定:過和4點畫直線1交2線于′點,過′點向上畫垂線2交三角面于″點,交2線于點,順著曲面′253′連接′為線3,與垂線4交于點1,1即是4降溫線與NaNO3的飽和溶液面′253′的交點。2點的確定:在曲面2上連接為線4,與垂線4交于點2,2即是4降溫線與NaNO3和冰的共晶面2的交點;3點的確定:連接″為線5,與垂線4交于點3,3即是4降溫線與三角面的交點。

利用該相圖還可以計算出溶液在降溫過程中每個階段的冰或鹽的析出量,確定溶液濃度的變化。因篇幅所限,本文不再詳述。

表2 三元體系不同組成溶液的冰點

表3 3個三元體系冰點實驗值與經驗公式計算值的比較

Note: Unit of concerntration is g·(100 g solution)-1.

2.4 三元體系的冰點及共晶點預測

2.4.1 經驗公式計算 三元體系變溫立體相圖的優勢在于能夠直觀地看出溶液在降溫過程中的相變化,缺點在于比較復雜,使用不便。由于海水、鹽湖鹵水以及一些含鹽廢水的冰點確定意義重大,本文通過補充測定一些位于圖1中′12′面上以及′1和′2(二元體系)邊上的溶液的冰點數據(見表2),結合上述測定的共晶點數據,擬合出了3個計算三元體系冰點的經驗公式,也即三元體系空間曲面′12′的方程,同時列出計算結果與實測冰點的比較,見表3。

總體預測結果較理想,個別值相對偏差較大,是因為這些點本身的冰點溫度(絕對值)較小,即使絕對偏差較小,相對偏差也較大。

2.4.2 熱力學模型計算 在工程應用上,經驗公式很實用,但缺乏理論意義。本文采用Alonso等[27]提出的以擴展UNIQUAC模型為基礎的過量吉布斯自由能模型來預測三元體系的冰點和共晶點。各離子的UNIQUAC模型參數由文獻[28-31]獲得。具體模型及計算過程見文獻[25, 32],在此不再詳述。模型計算的共晶點值和實驗值比較見表4。

表4 3個三元體系冰點實驗值與熱力學模型計算值的比較

總體來說,模型計算的共晶點值和實驗值吻合較好,計算值比實驗值略偏大,低濃度下計算結果更好。

3 結 論

(1)Na+//Cl-,NO-3-H2O、Na+//Cl-,SO42--H2O和Na+,K+//Cl--H2O 3個三元體系在冰點與共晶點溫度下無復鹽存在,各含有1個兩鹽與冰的共晶點,對應的共晶點溫度分別為-26.91℃、-21.31℃、-24.09℃,均低于所對應的兩種單鹽溶液的共晶點溫度,各含有兩個鹽結晶區(水合鹽與相應的單鹽合為1個),1個冰結晶區。

(2)采用三棱柱能直觀地表達三元水鹽體系的結晶區、冰點、共晶點和溶液組成的關系,可以確定任意組成的溶液在降溫過程中的冰鹽析出順序、冰點溫度、單鹽結晶點溫度、共晶點溫度和液相組成。

(3)應用本文所獲得的二次多項式經驗公式和文獻中過量吉布斯自由能模型來預測三元體系的冰點和共晶點,計算結果均較理想。

References

[1] SONG P S, LI W, SUN B,. Recent development on comprehensive utilization of salt lake resources[J]. Chinese Journal of Inorganic Chemistry, 2011, 27(5): 801-815.

[2] MA P. Sustainable exploitation and comprehensive utilization of salt lake resources in China[J]. Progress in Chemistry, 2009, 21(11): 2349-2357.

[3] SHEPARD M L, GOLDSTON C S, COCKS F H. The H2O-NaCl-glycerol phase diagram and its application in cryobiology[J]. Cryobiology, 1976, 13(1): 9-23.

[4] 朱巧麗, 黃雪莉.-15℃下Na+,K+,Mg2+//Cl-,SO42--H2O體系相平衡[J]. 化工學報, 2015, 66(4): 1252-1257. ZHU Q L, HUANG X L. Liquid-solid phase equilibrium of Na+,K+,Mg2+//Cl-,SO42--H2O system at-15℃[J].CIESC Journal, 2015, 66(4): 1252-1257.

[5] COCKS F H, BROWER W E. Phase diagram relationships in cryobiology[J]. Cryobiology, 1974, 11(4):340-358.

[6] PERALTA J M, RUBIOLO A C, ZORRILLA S E. Prediction of heat capacity, density and freezing point of liquid refrigerant solutions using an excess Gibbs energy model[J]. Journal of Food Engineering, 2007, 82(4): 548-558.

[7] 付振海, 張志宏, 馬艷芳, 等. 硫酸鹽型鹵水低溫處理及其液相蒸發析鹽規律的理論研究[J]. 無機鹽工業, 2013, 45(2): 29-32. FU Z H, ZHANG Z H, MA Y F,. Freeze sulfate-type brine at low temperature in stages and theoretical study on evaporation and crystallization behavior[J]. Inorganic Chemicals Industry, 2013, 45(2): 29-32.

[8] BEIER N, SEGO D, DONAHUE R,. Laboratory investigation on freeze separation of Saline mine waste water[J ]. Cold Regions Science and Technology, 2007, 48(3): 239-247.

[9] 張杰, 史學偉, 趙雙良, 等. 水鹽體系相平衡研究進展[J]. 化工學報, 2016, 67(2): 379-389. ZHANG J, SHI X W, ZHAO S L,.Progress in study on phase equilibria of salt-water systems[J]. CIESC Journal, 2016, 67(2): 379-389.

[10] GAYLE F W, COCKS F H, SHEPARD M L. The H2O-NaCl-sucrose phase diagram and applications in cryobiology[J]. Journal of Chemical Technology and Biotechnology, 1977, 27(5): 599-607.

[11] 伍倩, 鄭綿平, 乜貞, 等. 西藏當雄錯鹽湖鹵水冬季日曬蒸發實驗研究[J]. 地質學報, 2013, 87(3): 433-440. WU Q, ZHENG M P, NIE Z,. Experiment study of solar evaporation of brine form the Dangxiongcuo salt lake in winter[J]. Acta Geologica Sinica, 2013, 87(3): 433-440.

[12] GE X, WANG X. Estimation of freezing point depression, boiling point elevation, and vaporization enthalpies of electrolyte solutions[J]. Industrial & Engineering Chemistry Research, 2009, 48(4): 2229-2235.

[13] 李青海, 顧同欣, 于升松, 等. 南翼山油田鹵水低溫結晶過程研究[J]. 物理化學學報, 2011, 27(8): 1803-1808. LI Q H, GU T X, YU S S,. Study on the precipitation pathway of Nanyishan oilfield brine at subzero temperatures[J]. Acta Physico-Chimica Sinica, 2011, 27(8): 1803-1808.

[14] 張霞, 黃雪莉. 低溫下Na+//Cl-,SO42-,NO-3-H2O體系相平衡研究[J]. 化學通報, 2015, 78(4): 337-341. ZHANG X, HUANG X L. Study on the phase equilibrium of the quaternary system Na+//Cl-,SO42-,NO-3-H2O at low temperatures[J]. Chemistry Bulletin, 2015, 78(4): 337-341.

[15] SUN Y, ZHANG Y, QUAN X. Treatment of petroleum refinery wastewater by microwave-assisted catalytic wet air oxidation under low temperature and low pressure[J]. Separation & Purification Technology, 2008, 62(3): 565-570.

[16] GAO W, SMITH D W, SEGO D C. Ice nucleation in industrial wastewater[J].Cold Regions Science & Technology, 1999, 29(2): 121-133.

[17] LANGENHOFF A, STUCKEY D C. Treatment of dilute wastewater using an anaerobic baffled reactor: effect of low temperature[J]. Water Research, 2000, 34(15): 3867-3875.

[18] SILCOK H. Solubilities of Inorganic and Organic Compounds[M]. 3rd ed. Oxford: New York Pergamon Press, 1979: 1151.

[19] 牛自得, 程芳琴. 水鹽體系相圖及應用[M]. 天津: 天津大學出版社, 2002: 269. NIU Z D, CHENG F Q. The Phase Diagram of Salt-Water System and Its Application[M]. Tianjin: Tianjin University Press, 2002: 269.

[20] WEAST R C. CRC Handbook of Chemistry and Physics[M]. BocaRaton, USA: CRC Press, 1974: 2661.

[21] OAKES C S, BODNAR R J, SIMONSON J M. The system NaCl-CaCl2-H2O(Ⅰ): The ice liquidus at 1 atm total pressure[J]Geochimica Et Cosmochimica Acta, 1990, 54(3): 603-610.

[22] GIBBARD J R H F, GOSSMANN A F. Freezing points of electrolyte mixtures (Ⅰ): Mixtures of sodium chloride and magnesium chloride in water[J]. Journal of Solution Chemistry, 1974, 3(5): 385-393.

[23] DESNOYERS J E, BILLON M, LéGER S,. Salting out of alcohols by alkali halides at the freezing temperature[J]. Journal of Solution Chemistry, 1976, 5(10): 681-691.

[24] 王雪瑩, 黃文婷, 黃雪莉. 多元水鹽體系冰鹽共晶點的測定和圖形表達[J]. 化工學報, 2016, 67(5): 1687-1693. WANG X Y, HUANG W T, HUANG X L. Determination and graphics expression of ice-salt eutectic points of multicomponent salt-water systems[J]. CIESC Journal, 2016, 67(5): 1687-1693.

[25] 黃文婷. 水鹽體系冰點實驗研究與計算預測[D]. 烏魯木齊: 新疆大學, 2015. HUANG W T. The experimental determinatinon and prediction of freezing point of salt-water systems[D]. Urumqi: Xinjiang University, 2015.

[26] 張寒琦, 林英杰, 宋利珠, 等.實用化學手冊[M]. 北京: 科學出版社, 2001: 774. ZHANG H Q, LIN Y J, SONG L Z,. Handbook of Practical Chemistry[M]. Beijing: Science Press, 2001: 774.

[27] ALONSO H A T, PERALTA J M, RUBIOLO A C,. Prediction of the freezing point of multicomponent liquid refrigerant solutions[J]. Journal of Food Engineering, 2011, 104(1): 143-148.

[28] THOMSEN K, RASMUSSEN P, GANI R. Correlation and prediction of thermal properties and phase behaviour for a class of aqueous electrolyte systems[J]. Chemical Engineering Science, 1996, 51(14): 3675-3683.

[29] THOMSEN K. Aqueous electrolytes: model parameters and process simulation[D]. Denmark: Technical University of Denmark, 1997.

[30] GARCíA A V, THOMSEN K, STENBY E H. Prediction of mineral scale formation in geothermal and oilfield operations using the extended UNIQUAC model(Ⅰ): Sulfate scaling minerals[J]. Geothermics, 2005, 34(1): 61-97.

[31] GARCíA A V, THOMSEN K, STENBY E H. Prediction of mineral scale formation in geothermal and oilfield operations using the extended UNIQUAC model(Ⅱ): Carbonate-scaling minerals[J]. Geothermics, 2006, 35(3): 239-284.

[32] 黃文婷, 黃雪莉. 水鹽體系冰點的預測[J]. 鹽湖研究, 2015, 23(1): 38-43. HUANG W T, HUANG X L. The experimental determinatinon and prediction of freezing point of salt-water systems[J]. Journal of Salt Lake Research, 2015, 23(1): 38-43.

Determination, graphics expression and calculation of freezing point and eutectic point in ternary salt-water system

LIU Fuyun, HUANG Xueli, HUANG Wenting, WANG Xueying, PAN Yi

(Key Laboratory of Cleaner Transition of Coal & Chemicals Engineering of Xinjiang Uyghur Autonomous Region,College of Chemistry and Chemical Engineering, Xinjiang University, Urumqi 830046, Xinjiang, China)

In this paper, aiming at ternary salt-water systems of Na+//Cl-,NO-3-H2O, Na+//Cl-,SO42--H2O and Na+,K+//Cl--H2O, the freezing points, eutectic points and the crystallization regularities of salts were investigated; the three-dimensional phase diagrams of temperature-composition of the ternary systems were plotted according to experimental data. The results were as following: the relationships among the freezing point, eutectic point, phases and the compositions of solution for the ternary system can be expressed by triangular prism. In the three ternary systems, there were no double salts at the eutectic point temperature, and there was an eutectic point with two salts and ice, three eutectic lines with single salt and ice or two salts, a crystallization zone with ice and two single salt. The three-dimensional phase diagram of temperature-composition of solution can be used to describe the relationships among the crystalline regions, freezing points or eutectic points and the compositions of liquids for the ternary salt-water systems visually, and to determine the precipitation processes of salts or ice as cooling solution. Three empirical formulas about freezing points were established. The empirical formulas of freezing point and the excess Gibbs free energy model were used to predict the freezing point and eutectic point of the ternary systems, which turned out that the calculated values were in good agreement with the experimental values.

salt-water system; eutectic point; freezing point; phase diagram; calculation

10.11949/j.issn.0438-1157.20161809

O 642.5+4

A

0438—1157(2017)09—3336—07

2016-12-26收到初稿,2017-05-12收到修改稿。

黃雪莉。

劉福云(1990-),女,碩士研究生。

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

2016-12-26.

Prof. HUANG Xueli, xuelih@163.com

supported by the National Natural Science Foundation of China (21166022).

猜你喜歡
體系實驗
記一次有趣的實驗
微型實驗里看“燃燒”
構建體系,舉一反三
探索自由貿易賬戶體系創新應用
中國外匯(2019年17期)2019-11-16 09:31:14
做個怪怪長實驗
NO與NO2相互轉化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
如何建立長期有效的培訓體系
現代企業(2015年1期)2015-02-28 18:43:18
“曲線運動”知識體系和方法指導
《實驗流體力學》征稿簡則
主站蜘蛛池模板: 亚洲浓毛av| 国产在线观看91精品| 97视频精品全国在线观看| 亚洲综合色婷婷| www.日韩三级| 黄网站欧美内射| 欧美日韩另类国产| 国产情精品嫩草影院88av| 国产性精品| 国产成人高清精品免费软件| 99久久人妻精品免费二区| 综合五月天网| 老色鬼欧美精品| 伊人久久大香线蕉综合影视| 91视频区| 精品国产电影久久九九| 国产成人精品在线1区| 国产理论一区| 国产成人久久综合777777麻豆| 国产成人精品视频一区视频二区| 不卡无码h在线观看| 在线看AV天堂| 永久天堂网Av| 熟妇丰满人妻| 亚洲有无码中文网| 97视频免费在线观看| 热这里只有精品国产热门精品| 久热中文字幕在线| www.亚洲国产| 色偷偷一区| 午夜精品福利影院| 日本午夜三级| 亚洲第一视频网| 99九九成人免费视频精品| 亚洲av综合网| 欧美在线伊人| 免费人成视网站在线不卡| 精品99在线观看| 亚洲精品视频网| 日本道综合一本久久久88| 亚洲无码37.| 男女男免费视频网站国产| 日本91视频| 国产成人h在线观看网站站| 免费高清毛片| 国内精品手机在线观看视频| 国产精品密蕾丝视频| AV无码无在线观看免费| 欧美成人影院亚洲综合图| 日韩毛片免费| 久久青青草原亚洲av无码| 国产欧美视频一区二区三区| 色婷婷电影网| 四虎综合网| 91精品亚洲| 日本手机在线视频| 国产精品九九视频| 亚洲精品国偷自产在线91正片| 亚洲欧美日韩动漫| 麻豆精品在线播放| 国产一区二区三区免费| 激情无码视频在线看| 亚洲欧洲日本在线| 一级片一区| 精品人妻AV区| 国产福利免费视频| 日韩视频精品在线| 亚洲综合片| 一本二本三本不卡无码| 中文字幕色站| 99久久国产综合精品2020| 亚洲人成网站18禁动漫无码 | 男女猛烈无遮挡午夜视频| 91人妻日韩人妻无码专区精品| 欧美日韩国产精品综合| 婷婷亚洲最大| 日本在线国产| 国产乱人视频免费观看| 久996视频精品免费观看| 国产在线观看第二页| 亚洲日韩精品无码专区| 午夜久久影院|