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

吡啶硝酸鹽與醇二元混合體系物化性能研究

2012-12-26 06:58:32王建英梁麗亞昝昊搏胡永琪
河北科技大學學報 2012年5期
關鍵詞:體系

王建英,梁麗亞,昝昊搏,胡永琪

(1.河北科技大學化學與制藥工程學院,河北石家莊 050018;2.唐山職業技術學院臨床醫學部,河北唐山 063000)

吡啶硝酸鹽與醇二元混合體系物化性能研究

王建英1,梁麗亞1,昝昊搏2,胡永琪1

(1.河北科技大學化學與制藥工程學院,河北石家莊 050018;2.唐山職業技術學院臨床醫學部,河北唐山 063000)

采用MDY-2電子密度儀和德國DCAT21全自動表面張力儀,在298.15 K下測定了離子液體N-丁基吡啶硝酸鹽(N-butylpyridinium nitrate,[BuPy]NO3)與有機溶劑甲醇、乙醇、正丁醇在全濃度范圍內(離子液體在有機溶劑中的摩爾分數為0~1)的密度和表面張力,計算了二元體系{[BuPy]NO3+CH3OH},{[BuPy]NO3+C2H5OH}和{[Bu Py]NO3+C4H9OH}的超額摩爾體積V E和表面張力偏差δγ,分別利用Redlich-Kister方程對二元體系的超額摩爾體積V E和表面張力偏差δγ進行了擬合。結果表明:二元體系{[BuPy]NO3+CH3OH}和{[BuPy]NO3+C2H5OH}的超額摩爾體積V E在全摩爾組成范圍內均為負值,而{[BuPy]NO3+C4H9OH}二元體系在低離子液體摩爾分數時為正值,隨著離子液體摩爾分數的增大由正值變為負值。在全摩爾組成范圍內,{[BuPy]NO3+CH3OH}的表面張力偏差δγ為正值,{[BuPy]NO3+C2H5OH}和{[BuPy]NO3+C4H9OH}的δγ為負值,隨著醇碳鏈的增長,表面張力偏差出現由正到負的變化趨勢。

離子液體;密度;表面張力;超額摩爾體積

離子液體的物理性質數據是相關工業設計和開發的重要基礎[1]。離子液體的一些獨特的性質與它們的結構密切相關。一方面,通過改變陰、陽離子的不同組合來調節它們的物理性質,如密度、表面性能、極性、溶解性能、熔點等;另一方面,有機溶劑的加入也能顯著改變離子液體的物理性質,從而影響其性能。離子液體與有機溶劑二元體系在很多方面顯示出不同于分子溶劑的特性,因而在有機合成、催化反應等方面表現出了相當的優勢[2]。例如:離子液體1-丁基-3-甲基咪唑四氟硼酸鹽([BMIM][BF4])和1-丁基-3-甲基咪唑六氟磷酸鹽([BMIM][PF6])與醇類混合溶劑體系中進行親電氟化反應,其反應速率及選擇性比傳統溶劑有很大的提高[3];離子液體[BMIM][BF4]和二氯甲烷的混合溶液是硝酸鈰銨催化1,3-二羰基化合物和取代加成反應的優秀溶劑,反應速率比在乙腈中有大幅提高[4]。因此,開展離子液體與分子溶劑二元體系的物理性能研究對于拓寬與離子液體有關的工業應用具有非常重要的意義。

目前,有大量文獻報道了陰離子為[PF6],[BF4],[MeSO4],[SCN],[EtSO4],[MetSO4],[triflate],[NTf2],[dca],[EtOSO3],[CH3SO4],[OcSO4][5-8]的咪唑類離子液體和有機溶劑如乙醇、丁醇、戊醇、己醇、苯甲醛和N,N二甲基甲酰胺等有機溶劑二元混合體系的表面張力、密度、黏度、折光率、比熱容及其熱力學性質等。但是這些文獻大都是基于咪唑類離子液體和有機溶劑二元體系的研究,關于吡啶類離子液體和有機溶劑混合物的研究卻少見報道。

筆者系統測定了298.15 K下離子液體[BuPy]NO3和甲醇、乙醇、正丁醇混合的二元體系在全濃度范圍內的密度和表面張力,計算了3個二元體系的超額摩爾體積V E和表面張力偏差δγ,并對3個二元混合物的體積性質和表面性質進行了討論。

1 實驗材料與方法

1.1 儀器及試劑

MDY-2電子密度儀(上海方瑞儀器有限公司提供);DCAT21全自動表面張力儀(德國Dataphysics公司提供);FC204型電子天平(上海精科天平儀器廠提供);DC-3005A型低溫恒溫槽(寧波天恒儀器廠提供)。

N-丁基吡啶硝酸鹽離子液體([BuPy]NO3,純度為99%,購自上海成捷化學有限公司);甲醇、乙醇、正丁醇,均為色譜醇(購自天津市永大化學試劑開發中心)。

1.2 離子液體與醇二元體系密度及表面張力測定

在溫度為298.15 K、高純氮氣保護下,采用MDY-2電子密度儀和德國DCAT21全自動表面張力儀,測定了二次蒸餾水的密度和表面張力,測定值與文獻值一致。在高純氮氣保護下,用低溫恒溫槽進行循環水控制溫度,測量了離子液體[BuPy]NO3和甲醇、乙醇、正丁醇以不同摩爾比混和的二元體系的密度和表面張力。

2 結果與討論

2.1 離子液體[BuPy]NO3和甲醇、乙醇、正丁醇二元混合體系的密度及體積性質

從圖1可以看出,[BuPy]NO3和醇二元混合體系的密度值隨著離子液體醇溶液中離子液體摩爾分數的增大而呈現緩慢增加的趨勢。因而可以通過添加一定比例的醇來改變離子液體的物性,進而滿足其在不同領域的應用。

超額摩爾體積的研究有助于了解[BuPy]NO3和醇二元混合體系中的分子間相互作用信息[9]。由實驗測得[BuPy]NO3和醇二元混合體系的密度值,通過式(1)計算得到超額摩爾體積:

式中:xIL為混合物中[BuPy]NO3的摩爾分數;ρ,ρ1,ρ2分別是二元混合體系的密度、純離子液體[BuPy]NO3的密度和溶劑醇的密度;M1和M2分別是[BuPy]NO3和溶劑醇的分子質量。超額摩爾體積由 Redlich-Kister[10-12]方程(簡稱 R-K 方程)進行擬合,即

式中:Z代表超額性質;xIL是離子液體[BuPy]NO3的摩爾分數;B p是擬合參數;M是用Marquardt多項式優化的多項式擴張階數。

由式(2)得到的實驗值和理論值的標準偏差σ為

式中:Zexp為實驗值;Zcal為計算值;ndat為實驗組數。

298.15 K時,[BuPy]NO3和醇二元體系超額摩爾體積V E隨離子液體摩爾分數變化曲線見圖2。

圖1 T=298.15 K,[BuPy]NO3和醇二元體系密度ρ隨離子液體摩爾分數變化曲線Fig.1 Curve of densityρof binary mixture of[BuPy]NO3 and alcohol against the molar fraction at T=298.15 K

從圖2可以看出,{[BuPy]NO3+ 甲醇}和{[BuPy]NO3+乙醇}二元混合體系的超額摩爾體積V E在全摩爾組成范圍內均為負值。這種現象與二元體系{[BuPy]BF4+甲醇},{[BuPy]BF4+ 乙醇}的V E變化趨勢相似[4]。對二元體系{[BuPy]NO3+甲醇}來說,在xIL<0.43時,超額摩爾體積隨著離子液體摩爾組成的增大而減小,在xIL=0.43時達到最小值;在xIL>0.43時,超額摩爾體積呈現出隨離子液體摩爾分數增大而增大的趨勢。而對于{[BuPy]NO3+乙醇}體系,在xIL<0.59時,超額摩爾體積隨著離子液體摩爾分數的增大而減小,在xIL=0.59處出現最小值,當xIL>0.59時,超額摩爾體積隨著摩爾分數的增加而增大。由圖2可知二元體系{[BuPy]NO3+乙醇}超額摩爾體積隨摩爾分數變化曲線是幾乎對稱的,而[BuPy]BF4和醇混合的二元體系的V E曲線呈現出不同的形狀[4]。

與{[BuPy]NO3+ 甲醇}和{[BuPy]NO3+ 乙醇}2個體系相比,{[BuPy]NO3+正丁醇}二元體系超額摩爾體積在全摩爾組成范圍內呈現出不同的變化趨勢。由圖2可以看出,V E隨摩爾分數的變化曲線在全摩爾組成范圍內呈倒S型,在xIL為0.2和0.74時分別出現最大值和最小值,在xIL=0.44時,{[BuPy]NO3+正丁醇}的超額摩爾體積由正變負,曲線的變化趨勢和吡啶離子液體[Pyrr][C7CO2]與正丁醇二元混合體系的變化趨勢相同[13]。

圖2 [BuPy]NO3和醇二元體系超額摩爾體積V E隨離子液體摩爾分數變化曲線Fig.2 Excess molar volumes V E of binary mixture of[BuPy]NO3 and alcohol against the molar fraction

{[BuPy]NO3+甲醇}和{[BuPy]NO3+乙醇}體系的超額摩爾體積在全摩爾組成范圍內呈現負值,說明不同分子間的交叉締合強于同分子間的自締合,組分間分子大小的不對稱產生填充[18],故超額摩爾體積為負。超額摩爾體積為正值即二元混合體系與理想溶液形成正偏差,這表明在混合過程中,有機溶劑醇降低甚至打破了純離子液體的陰離子和陽離子之間強烈的靜電吸引作用。對于3個二元混合體系來說,{[BuPy]NO3+甲醇}混合物的超額摩爾體積最小,{[BuPy]NO3+正丁醇}的超額摩爾體積最大,說明離子液體[BuPy]NO3與甲醇分子間的作用力最大,與乙醇的次之,與正丁醇的分子間作用力最小。二元體系超額摩爾體積出現最小值的原因是在這個摩爾分數時,二元體系分子間相互作用力最大,從而導致體積最小。

2.2 離子液體[BuPy]NO3和甲醇、乙醇、正丁醇二元混合體系的表面張力

圖3為298.15 K下,[BuPy]NO3和甲醇、乙醇、正丁醇3個二元混合體系的表面張力(γ)隨離子液體摩爾分數變化的結果,表面張力偏差隨離子液體摩爾分數變化的曲線見圖4。

圖3 [BuPy]NO3和醇二元混合體系表面張力隨離子液體摩爾分數變化曲線Fig.3 Curve of surface tension of binary mixture of[BuPy]NO3 and alcohol against the molar fraction

圖4 [BuPy]NO3和醇二元混合體系表面張力偏差隨離子液體摩爾分數變化曲線Fig.4 Surface tension deviations of binary mixture of[BuPy]NO3 and alcohol against the molar fraction

從圖3可以看出,[BuPy]NO3與甲醇、乙醇、正丁醇二元混合體系在全摩爾組成范圍內的表面張力值介于純離子液體與醇之間,且隨離子液體摩爾分數的增加而增大。對于{[BuPy]NO3+甲醇},{[Bu Py]NO3+乙醇}和{[BuPy]NO3+正丁醇}這3個二元混合體系來說,隨著脂肪醇碳鏈的增長,表面張力曲線隨離子液體摩爾分數的增加趨勢變緩,這說明分子間作用力對表面張力起著決定性的作用。

[BuPy]NO3與乙醇混合的表面張力隨離子液體摩爾分數的增加幾乎呈線性變化,這表明在全摩爾組成范圍內,乙醇分子和離子液體在表面相為均勻混合。而[BuPy]NO3和甲醇、正丁醇二元混合體系的表面張力隨離子液體摩爾分數變化的趨勢與{[BuPy]NO3+乙醇}二元混合體系不同,{[BuPy]NO3+甲醇}二元混合體系的表面張力在xIL<0.3時增加較快,但在xIL>0.3后增長趨勢變緩。而{[BuPy]NO3和正丁醇}二元混合體系的表面張力在xIL<0.6時增加較緩,而在xIL>0.6后迅速增長。即{[BuPy]NO3+甲醇}和{[BuPy]NO3+正丁醇}這2個二元混合體系的表面張力隨著離子液體摩爾組成變化曲線均出現拐點。這種表面張力隨離子液體摩爾組成變化曲線出現拐點的現象與[EMIM][CH3SO4]和甲醇、乙醇及正丁醇混合的二元混合體系是相似的[6]。

表面張力偏差(δγ)通過式(4)計算得出:

式中:γ和γi分別代表混合體系的表面張力和組分i的表面張力;x i代表組分i的摩爾分數。

T=298.15 K,[BuPy]NO3和醇二元混合體系表面張力偏差隨離子液體摩爾分數變化曲線見圖4。圖4是用R-K方程擬合的表面張力偏差,擬合方程如下:

由圖4可以看出,二元混合體系{[BuPy]NO3+甲醇}的表面張力偏差在全摩爾組成范圍內均為正值,且在xIL=0.3處出現最大值;二元混合體系{[BuPy]NO3+乙醇}和{[BuPy]NO3+正丁醇}的表面張力偏差在全摩爾組成范圍內均為負值,且在xIL=0.79處出現最小值。由此可見:隨著脂肪醇碳鏈的增長,離子液體[Bu Py]NO3和脂肪醇混合的二元體系的表面張力偏差呈現由正到負的變化趨勢,這種現象與二元混合體系[EMIM][CH3SO4]和甲醇、乙醇、正丁醇的表面張力偏差隨摩爾組成的變化趨勢相似[6]。

表面張力偏差的擬合參數和標準偏差列于表1。

表1 T=298.15 K,[BuPy]NO3和醇二元混合體系的擬合參數和標準偏差Tab.1 Fitting parameters and root mean square deviations for binary mixture of[BuPy]NO3 and alcohol at T =298.15 K

3 結 語

二元體系[BuPy]NO3和甲醇、乙醇、正丁醇混合物的密度和表面張力隨著摩爾組成的增加而增大。由密度值計算得到3個二元混合體系的超額摩爾體積V E,并用Redlich-Kister方程進行了擬合。二元混合體系{[BuPy]NO3+甲醇}和{[BuPy]NO3+乙醇}的超額摩爾體積V E在全摩爾組成范圍內均為負值,{[BuPy]NO3+正丁醇}在低摩爾分數時為正值,且隨著離子液體摩爾分數的增大而由正值變為負值。

二元混合體系{[BuPy]NO3+乙醇}的表面張力隨著摩爾組成幾乎呈線性變化。{[BuPy]NO3+甲醇}和{[BuPy]NO3+正丁醇}二元混合體系表面張力隨離子液體摩爾分數的變化出現拐點。由二元混合體系的表面張力值計算得到表面張力偏差δγ,并用Redlich-Kister方程擬合。在全摩爾組成范圍內,{[BuPy]NO3+甲醇}的表面張力偏差δγ為正值,{[BuPy]NO3+乙醇}和{[BuPy]NO3+正丁醇}的表面張力偏差δγ為負值。隨著醇碳鏈的增長,表面張力偏差出現由正到負的變化趨勢。

[1]ANOUTI M,VIGEANT A,JACQUEMIN J,et al.Volumetric properties,viscosity and refractive index of the protic ionic liquid,pyrrolidinium octanoate,in molecular solvents[J].J Chem Thermodynamics,2010,42(7):834-845.

[2]GARCIA-MARDONES M,PEREZ-GREGORIO V,GUERRERO H,et al.Thermodynamic study of binary mixtures containing 1-butylpyridinium tetrafluoroborate and methanol,or ethanol[J].J Chem Thermodynamics,2010,42(12):1 500-1 505.

[3]BAUDOUX J,SALIT A F,CAHARD D,et al.Ionic liquids as solvents of choice for electrophilic fluorination:fluorination of Indoles by F-TEDA-BF4[J].Tetra Lett,2002,43(37):6 573-6 574.

[4]BAR G,PARSONS A F,THOMAS C B.Manganese(Ⅲ)acetate mediated radical reactions in the presence of an ionic liquid[J].Chem Commun,2001,15:1 350-1 351.

[5]DOMANSKA U,KROLIKOWSKA M.Effect of temperature and composition on the surface tension and thermodynamic properties of binary mixtures of 1-butyl-3-methylimidazolium thiocyanate with alcohols[J].Journal of Colloid and Interface Science,2010,348(2):661-667.

[6]WANG Jian-ying,ZHAO Feng-yun,LIU Yu-min,et al.Thermophysical properties of pure 1-ethyl-3-methylimida zolium methylsulphate and its binary mixtures with alcohols[J].Fluid Phase Equilibria,2011,305(2):114-120.

[7]de CASTRO C A N,LANGA E,MORAIS A L,et al.Studies on the density,heat capacity,surface tension and infinite dilution diffusion with the ionic liquids[C4mim][NTf2],[C4mim][dca],[C2mim][EtOSO3]and[Aliquat][dca][J].Fluid Phase Equilibria,2010,294(1/2):157-179.

[8]DOMANSKA U,POBUDKOWSKA A,ROGALSKI M.Surface tension of binary mixtures of imidazolium and ammonium based ionic liquids with alcohols,or water:Cation,anion effect[J].Journal of Colloid and Interface Science,2008,322(1):342-350.

[9]ZHANG Suo-jiang,LI Xian,CHEN Hui-ping,et al.Determination of physical properties for the binary system of 1-ethyl-3-methylimidazolium tetrafluoroborate+H2O[J].J Chem Eng Data,2004,49(4):760-764.

[10]GOMEZ E,GONZALEZ B,CALVAR N,et al.Physical properties of pure 1-ethyl-3-methylimidazolium ethylsulfate and its binary mixtures with ethanol and water at several temperatures[J].J Chem Eng Data,2006,51(6):2 096-2 102.

[11]PEREIRO A B,RODRíGUEZ A.Thermodynamic properties of ionic liquids in organic solvents from(293.15 to 303.15)K[J].J Chem Eng Data,2007,52(2):600-608.

[12]GONZALEZ E J,GONZALEZ B,CALVAR N,et al.Physical properties of binary mixtures of the ionic liquid 1-ethyl-3-methylimidazolium ethyl sulfate with several alcohols atT=(298.15,313.15,and 328.15)K and atmospheric pressure[J].J Chem Eng Data,2007,52(5):1 641-1 648.

[13]ANOUTI M,VIGEANT A,JACQUEMIN J,et al.Volumetric properties,viscosity and refractive index of the protic ionic liquid,pyrrolidinium octanoate,in molecular solvents[J].Journal of Chemical Thermodynamics,2010,42(7):834-845.

[14]ZAFARANI-MOATTAR M T,SHEKAARI H.Application of prigogine-flory-patterson theory to excess molar volume and speed of sound of 1-n-butyl-3-methylimidazolium hexafluorophosphate or 1-n-butyl-3-methylimidazolium tetrafluoroborate in methanol andacetonitrile[J].J Chem Thermodynamics,2006,38(11):1 377-1 384.

[15]魏珺芳,王延吉,勵 娟,等.單分散亞微HZSM-5分子篩的制備及其環已烯水合性能[J].河北工業大學學報(Journal of Hebei University of Technology),2011,40(4):20-25.

[16]郭 斌,高 亮,段二紅,等.[bmim][BF4]、[bmim][PF6]離子液體及其與DMF二元體系在293.15~353.15 K的密度和體積性質[J].河北科技大學學報(Journal of Hebei University of Science and Technology),2012,33(2):135-141.

[17]劉寶友,王園園,魏福祥.離子液體走向工業化應用的兩個重要理論問題研究[J]河北工業科技(Hebei Journal of Industrial Science and Technology),2009,26(2):124-126.

Physical and chemical properties of binary mixture ofN-butylpyridinium nitrate and alcohol

WANG Jian-ying1,LIANG Li-ya1,ZAN Hao-bo2,HU Yong-qi1
(1.College of Chemical and Pharmaceutical Engineering,Hebei University of Science and Technology,Shijiazhuang Hebei 050018,China;2.Clinical Medicine Department,Tangshan Vocational and Technical College,Tangshan Hebei 063000,China)

ionic liquid;density;surface tension;excess molar volume

O645.13

A

In this paper,MDY-2 electronic densitometer and Germany DCAT21 automatic surface tension meter are used to measure the density and surface tension over the whole concentration range(the molar fraction of ionic liquids in oganic solvents is 0~1)for the binary mixtures of ionic liquidsN-butylpyridinium nitrate([BuPy]NO3)and organic solvents such as methanol,ethanol and butanol at the temperature of 298.15 K.The excess molar volumesV Eand the surface tension deviationsδγof binary mixtures{[BuPy]NO3+CH3OH},{[BuPy]NO3+C2H5OH}and{[BuPy]NO3+C4H9OH}are determined.VEandδγare fitted by using the Redlich-Kister equation.The results show that the excess molar volumes of binary mixtures of{[BuPy]NO3+CH3OH}and{[BuPy]NO3+C2H5OH}are negative over the whole composition range,while theV Eof binary mixture{[BuPy]NO3+C4H9OH}has positive values at low molar fraction of ionic liquid,passing through a maximum and then decreases and becomes negative,showing minimum at higher ionic liquid molar fraction.The surface tension deviationsδγof{[BuPy]NO3+CH3OH}systems are positive,but those of{[BuPy]NO3+C2H5OH}and{[BuPy]NO3+C4H9OH}systemsare negative over the entire molar fraction range.It seems that surface tension deviationsδγvarys from positive to negative with the increase of alkyl chain length of alcohols.

1008-1542(2012)05-0400-06

2012-05-04;

2012-09-05;責任編輯:張士瑩

國家自然科學基金資助項目(20776037);河北省自然科學基金資助項目(B2012208084)

王建英(1973-),女,河北靈壽人,副教授,博士,主要從事離子液體方面的研究。

胡永琪教授。E-mail:huyq@hebust.edu.cn

猜你喜歡
體系
TODGA-TBP-OK體系對Sr、Ba、Eu的萃取/反萃行為研究
“三個體系”助力交通安全百日攻堅戰
杭州(2020年23期)2021-01-11 00:54:42
構建體系,舉一反三
探索自由貿易賬戶體系創新應用
中國外匯(2019年17期)2019-11-16 09:31:14
常熟:構建新型分級診療體系
中國衛生(2015年12期)2015-11-10 05:13:40
如何建立長期有效的培訓體系
現代企業(2015年1期)2015-02-28 18:43:18
E-MA-GMA改善PC/PBT共混體系相容性的研究
汽車零部件(2014年5期)2014-11-11 12:24:28
“曲線運動”知識體系和方法指導
加強立法工作 完善治理體系
浙江人大(2014年1期)2014-03-20 16:19:53
日本終身學習體系構建的保障及其啟示
主站蜘蛛池模板: 91精品免费久久久| 欧美一道本| 在线精品自拍| 国产极品美女在线| 欧美午夜在线视频| 九九精品在线观看| 中国一级特黄大片在线观看| AV网站中文| 久久久久亚洲Av片无码观看| 欧美午夜视频在线| 在线看片国产| 超碰aⅴ人人做人人爽欧美 | 亚洲久悠悠色悠在线播放| 一本一道波多野结衣一区二区 | 久久国产精品77777| 亚洲高清在线天堂精品| 中文字幕在线永久在线视频2020| 美女高潮全身流白浆福利区| 国产亚洲精品资源在线26u| 国产在线视频福利资源站| 国产aaaaa一级毛片| 伊人久久青草青青综合| 国产精品jizz在线观看软件| 99中文字幕亚洲一区二区| 性色一区| 99在线视频免费| 亚洲欧美成人在线视频| 精品视频在线观看你懂的一区| 青青青国产精品国产精品美女| 在线欧美国产| 制服丝袜国产精品| 呦女精品网站| 成AV人片一区二区三区久久| 在线a视频免费观看| 免费人成网站在线观看欧美| 毛片网站观看| 亚洲国产高清精品线久久| 亚洲欧美不卡中文字幕| 欧美精品v| 国产精品尤物铁牛tv| 国产成人亚洲毛片| 国产成人精品视频一区二区电影 | 国产成人综合在线观看| 国产99欧美精品久久精品久久| 久久青青草原亚洲av无码| 青青草原国产精品啪啪视频| 精品国产免费观看| 久久成人免费| 国产麻豆aⅴ精品无码| 91蜜芽尤物福利在线观看| 啊嗯不日本网站| 国产粉嫩粉嫩的18在线播放91| 一本一道波多野结衣一区二区| 国产无码制服丝袜| 日韩欧美在线观看| 熟妇无码人妻| 欧美区一区| 99色亚洲国产精品11p| 2021精品国产自在现线看| 凹凸国产分类在线观看| 亚洲一级毛片| 香蕉视频国产精品人| 久久无码高潮喷水| 激情六月丁香婷婷四房播| 综合色亚洲| 人妻无码AⅤ中文字| 婷婷色在线视频| a天堂视频在线| 久久亚洲日本不卡一区二区| 大香伊人久久| 91探花国产综合在线精品| 亚洲欧美在线综合图区| 亚洲人成色在线观看| 色综合综合网| 综1合AV在线播放| 久久精品亚洲中文字幕乱码| 国产美女91呻吟求| 国产美女视频黄a视频全免费网站| 黑色丝袜高跟国产在线91| 亚洲中久无码永久在线观看软件| 国产亚洲精品91| 国产精品3p视频|