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

GSS362水溶液動態表面張力及吸附動力學研究

2017-07-19 13:09:06喻紅梅朱國華劉駿峰孔德存
石油化工 2017年1期
關鍵詞:擴散系數

喻紅梅,華 平,朱國華,劉駿峰,孔德存,劉 倩

(1. 南通大學 化學化工學院,江蘇 南通 226019;2. 南通大學 紡織服裝學院,江蘇 南通 226019)

GSS362水溶液動態表面張力及吸附動力學研究

喻紅梅1,華 平1,朱國華1,劉駿峰2,孔德存1,劉 倩1

(1. 南通大學 化學化工學院,江蘇 南通 226019;2. 南通大學 紡織服裝學院,江蘇 南通 226019)

用最大泡壓法測定了磺酸鹽型雙子表面活性劑1,3-丙二醇雙子琥珀酸二異辛酯磺酸鈉(GSS362)在不同溫度和濃度下的動態表面張力;采用Rosen模型對動態表面張力進行了表征,同時用Ward-Tordai方程對GSS362的吸附機理進行了探討。實驗結果表明,GSS362的濃度越大,動態表面張力降低得越快,且介平衡值越接近平衡值;溫度升高有利于GSS362在新鮮表面的吸附;GSS362在新鮮表面的吸附前期為擴散控制吸附,吸附后期為混合動力控制吸附。

琥珀酸酯磺酸鹽;雙子表面活性劑;動態表面張力;吸附機理

表面活性劑降低水溶液的表面張力是由于表面活性劑分子不斷取代表面的水分子,隨著表面活性劑分子不斷取代水分子,表面張力逐漸降低,最終在表面達到飽和,此時的張力為平衡表面張力,而在表面活性劑分子達到飽和前的表面張力為動態表面張力[1]。在某些場合,對溶液動態吸附性能的研究比平衡狀態下吸附性能的研究更重要,也更能與表面活性劑的實際應 用相聯系[2-3]。三次采油中,驅油體系在油層的毛細孔道中運行,與原油瞬時接觸,了解表面活性劑水溶液與原油接觸瞬間的表面張力,可更好地研究如何提高采油率[4]。這些應用均與表面活性劑的動態吸附性能有關。最常見的動態吸附性能測定方法是動態表面張力的測定。

雙磺基琥珀酸酯磺酸鈉(GSS)通過柔性連接基團將兩個傳統的琥珀酸酯鹽連接,連接基團使得GSS分子中兩個離子頭基間的靜電斥力大大削弱,同時加強了碳氫鏈間的疏水結合力,使GSS具有比單鏈琥珀酸酯磺酸鹽更低的臨界表面張力[5-6]。鮮有文獻報道GSS的動態表面張力,因此有必要對其動態表面性能進行研究,為其實際應用提供理論基礎。

本工 作測定了1,3-丙二醇雙子琥珀酸二異辛酯磺酸鈉(GSS362)的靜態表面張力,同時用最大泡壓法研究了GSS362的動態表面張力,并對濃度和溫度的影響進行了探討;采用Rosen模型對動態表面張力進行表征,同時采用Ward-Tordai方程對動態吸附機理進行了探討。

1 實驗部分

1.1 試劑和儀器

GSS362:參照文獻[7]的方法合成,并參照文獻[8]的方法進行提純,經測定表面張力曲線無最低點,說明不含高活性雜質[9];實驗用水均為二次去離子水經亞沸重蒸處理。

JK99C型全自動表面張力儀(測量精度0.05 mN/m):上海中晨數字技術設備有限公司;SHA-C型水浴恒溫振蕩器(控溫精度±0.5 ℃):金壇市恒豐儀器廠;BP100型動態表面張力儀:德國Kruss公司。

1.2 表面性能的測定

1.2.1 系列濃度表面活性劑水溶液的配制

待測GSS362溶液的濃度范圍設定為1×10-6~1×10-2mol/L。首先配制一定量濃度為1×10-2mol/L的GSS362水溶液作為母液備用,然后移取適量母液逐級稀釋,得到一系列不同濃度的溶液,在設定溫度的振蕩水浴鍋中恒溫30 min。

1.2.2 靜態表面張力的測定

采用JK99C型表面張力儀由低濃度至高濃度測定恒溫后系列溶液的表面張力(γ),做γ-lgc(c為濃度)曲線。測試均在帶有恒溫夾套的玻璃杯中進行,外接數控超級恒溫槽,溫度控制在預設溫度,溫度選擇為20,30,40 ℃。

1.2.3 動態表面張力的測定

用BP100型動態表面張力儀由低濃度至高濃度測定恒溫后系列溶液的動態表面張力。測試均在帶有恒溫夾套的玻璃杯中進行,外接數控超級恒溫槽,溫度控制在預設溫度,溫度選擇20,30,40 ℃。

2 結果與討論

2.1 GSS362水溶液表面張力的測定

不同溫度下GSS362水溶液的靜態表面張力見圖1,介平衡表面張力與平衡表面張力見表1。

圖1 不同溫度下GSS362水溶液的靜態表面張力Fig.1 Surface tension(γ) vs. lgc of 1,3-propanediol bis-isooctyl sulfo succinate(GSS362) aqueous solutions at different temperature.c:concentration of GSS362 in aqueous solution. Temperature/℃:■ 20;● 30;▲ 40

表1 不同溫度下的介平衡表面張力與平衡表面張力Table 1 Medium equilibrium surface tension(γm) and equilibrium surface tension(γeq) of the GSS362 aqueous solutions at different temperature

不同溫度下GSS362的動態表面張力見圖2。

由圖2可知,不同濃度下,GSS362溶液動態表面張力下降的趨勢不一樣,濃度較低時(lgc= -4.52),圖1的靜態表面張力遠低于溶劑表面張力,而測得的溶液動態表面張力等于溶劑表面張力,因為,在濃度較低時,溶液相中表面活性劑分子濃度較低,分子擴散至表面的速率較慢,因而動態表面張力降低需要較長時間。在lgc=-3.00時,可看出較完整的動態表面張力降低曲線,開始時表面張力接近溶劑的表面張力,隨時間的延長表面張力先快速下降,降至一定數值后,降幅趨緩,在一段時間內基本保持不變。因為氣/液界面新生成時,GSS362尚未在表面吸附,界面張力為溶劑的表面張力,隨著時間的延長,表面活性劑分子開始在表面吸附,氣/液界面張力降低,隨著吸附的進行,表面逐步趨于飽和,表面活性劑分子在表面吸附的速率下降,因而動態表面張力下降速率變慢,趨于平衡。而表面活性劑分子在表面達到飽和吸附量還需要進一步的吸附與定向排列過程,因而觀察到動態表面張力曲線中的張力值較對應濃度下的靜態表面張力值高。而水溶液中GSS362濃度較高時(c>cmc,cmc為臨界膠束濃度),表面活性劑分子很快吸附到表面,開始測定時,新生成的氣/液界面吸附較多的表面活性劑分子,因而動態表面張力低于溶劑的表面張力,且隨著濃度的增加,初始動態表面張力降低。當c>cmc時,為表面活性劑膠束溶液,在溶液表面上吸附的表面活性劑分子達到飽和,達到平衡張力并保持不變。膠束似乎不具備表面活性,事實上膠束存在時表面活性劑溶液表面的吸附情況比無膠束時復雜得多。從分子角度看,溶液中存在膠束形成和分離的過程,從而導致了溶液中表面活性劑單體在體相和表面間的傳遞過程,即在表面活性劑的膠束溶液中,單體與膠束都參與了動態吸附。

圖2 不同溫度下GSS362的動態表面張力Fig.2 Dynamic surface tension(γ(t)) of GSS362 aqueous solutions at different temperature.lgc:■ -2.00;● -2.30;▲ -2.52;▼ -2.70;◆ -3.00;? -3.52;? -4.00;★ -4.52

2.2 動態表面張力的表征

Rosen等[10]研究動態表面張力隨時間變化時提出,許多表面活性劑溶液的動態表面張力隨時間的變化曲線可分為4個階段:誘導區、表面張力快速下降區、介平衡區和平衡區。前3個區域的動態表面張力可用以下經驗方程表示:

式中,γ0為溶劑的表面張力,mN/m;γ(t)為表面活性劑溶液t時刻的動態表面張力,mN/m;γm為介平衡時的表面張力,mN/m(部分濃度下的介平衡表面張力見表1);t為時間,s;n與t*均為常數,t*的單位與t的單位相同,n為無因次量。對式(1)兩邊取對數可得:

測得的動態表面張力數據經式(2)處理后,令lg[(γ0- γ(t))/(γ(t) - γm)]= lgK,以lgK對lgt做圖,見圖3。由圖3可看出,隨t的增大,γ(t)減小,lgK增大,這些均反映了能量的變化。從直線的斜率和截距可得到動態表面張力的特性參數n與t*,并可進一步得到誘導區結束時間ti、介平衡開始時間tm及t1/2(體系的表面壓達到介平衡時表面壓二分之一的時間)時表面張力的下降速率R1/2[11]:

不同溫度和濃度下GSS362溶液的動態表面張力參數見表2。從表1和表2可看出,溫度一定時,表面活性劑溶液濃度越大,γm越小,越接近于平衡張力γeq,t*,ti,tm越小,R1/2越大,n變化不大。濃度一定時,溫度越高,γm越小,t*,ti,tm越小,R1/2越大,n越小。

GSS362濃度越大,擴散速率越大,分子在氣/液界面吸附越快,則GSS362的動態表面張力下降越快,誘導區結束越早(ti越?。?,GSS362分子在氣/液界面越易達到介平衡(tm越短);溫度越高,分子擴散越快,則GSS362分子擴散至新鮮表面的速率越快,動態表面張力下降越快,誘導區結束越早(ti越小),GSS362分子在氣/液界面越易達到介平衡(tm越短)。這與圖3的曲線趨勢和計算得到的tm值一致。

n值反映了吸附初期(極短的時間內,t→0)表面活性劑分子從本體溶液擴散到面下層的過程,對于離子表面活性劑,n越小,擴散勢壘越小,表面活性劑擴散越快。GSS362濃度范圍為lgc= -2.00~-2.70時,表面活性劑溶液為膠束溶液,溶液中GSS362單體的濃度保持不變,表面活性劑擴散速率較接近。溫度升高,分子運動加快,GSS362擴散速率增大,而n隨溫度的升高而減小,即擴散勢壘減小。

圖3 GSS362 在不同溫度下的lgK~lgt曲線Fig.3 lgK vs. lgt at different temperature.lgc:■ -2.00;● -2.30;▲ -2.52;▼ -2.70;◆ -3.00

t*反映了吸附后期(t→∞)表面活性劑從面下層擴散到溶液表面的過程,t*越小,吸附勢壘越大。當lgc=-3.00~-2.00時,GSS362濃度越大,相同時間內擴散至表面的GSS362分子越多,GSS362繼續在表面吸附所需克服的能壘越大,這與濃度越高t*越小相吻合。

溫度對動態表面張力的影響主要反映在:溫度越高,表面活性劑分子運動越劇烈,分子擴散至新鮮界面的速率越快,因而動態表面張力的降低速度越快。

R1/2越大表明動態表面活性越大[12]。濃度越大,溫度越高,R1/2越大。

表2 不同溫度和濃度下GSS362溶液的動態表面張力參數Table 2 Dynamic surface tension parameters of the GSS362 aqueous solutions with different concentration at different temperature

2.3 動態吸附模型

表面活性劑分子在新鮮表面上的吸附分兩步:1)分子從本體溶液到次表面上的遷移;2)分子在次表面和表面間的吸附平衡。在上述動態吸附模型的基礎上提出了兩種主要的動態吸附機理:擴散控制吸附、混合動力控制吸附[13]。

對于c<cmc的表面活性劑溶液,Ward和Tordai根據擴散控制機理提出了描述表面吸附的動力學模型[14]:

式中,Γ(t)為t時刻新鮮表面上表面活性劑的吸附量,mmol/m2;c0為表面活性劑溶液本體溶液濃度,mol/L;D為表面活性劑分子擴散系數,m2/s;cs為次表面層中表面活性劑溶液濃度,mol/L,為時間的函數;τ為虛擬變量。其中表示分子從體相到面下層的遷移表示隨著面下層濃度的增加分子由面下層擴散回體相。由于反擴散項積分無法計算,對于離子型表面活性劑,Fainerman等用漸進的方法對Ward-Tordai方程進行了修正[15]。

式中,Гeq表示平衡時的表面吸附量,mmol/m2。對于雙子表面活性劑,n=3,Гeq由平衡時的γ~c曲線經計算機處理得到。

對于混合動力控制吸附模型,Liggieri和Raveva引進的校正擴散系數D*與根據式(6)得到的D的關系[16]為:

當εa→0時,D*→D,混合動力控制過程趨近于擴散控制過程?;旌蟿恿^程可看成存在一定能壘的擴散過程,修正后的方程為:

式中,Da=D*2/D=Dexp(-2εa/RT),當吸附后期為擴散控制吸附時,Da= D。

由式(5)和(6)可看出,溶液的γ(t)(t→0)與t1/2、γ(t)(t→∞)與t-1/2呈線性關系。

對20 ℃下的動態表面張力曲線的吸附機理進行討論。以γ(t)分別對t1/2和t-1/2做圖,見圖4~5,直線的截距為對應濃度下的平衡張力,分子擴散系數可由直線部分斜率求出,計算結果見表3。

從表3可知,對于c<cmc的GSS362溶液(lgc = -4.00~-3.00),在吸附初期,隨表面活性劑濃度的增大,分子擴散系數增大。這是因為溶液濃度增大,溶液中游離的表面活性劑分子數增加,分子擴散加快,因而在吸附初期表面活性劑分子的擴散速率快。在吸附后期,隨表面活性劑濃度的增加,溶液的平衡表面張力降低,分子擴散系數降低,且較吸附初期的擴散系數低并有數量級的差異。這與單一擴散控制吸附機理不符,說明有吸附能壘的存在。對于離子型表面活性劑,在吸附后期,溶液表面被相當數量的表面活性劑分子占據,形成吸附膜,帶電的極性頭基排列緊密,電荷密度較大,阻礙了新的表面活性劑分子從面下層到表面的吸附;溶液表面的表面活性劑分子越多,則剩余空間越小,對擴散產生影響;此外,表面的分子處于吸附和脫附的動態平衡中,吸附層分子數增多使得脫附速率增大[11],即表面活性劑在表面的脫附作用使得表面層與次表面層間的分子交換存在能壘。這與Rosen方程的處理結果,即隨著濃度的增加,表面活性劑分子的吸附勢壘增加相一致。所以可以推斷,GSS362分子在吸附初期為擴散控制吸附,在吸附后期為混合動力控制吸附。

圖4 不同濃度的GSS362溶液的γ(t)~t1/2曲線Fig.4 Dynamic surface tensions of the GSS362 solutions vs. t1/2.lgc:■ -2.00;● -2.30;▲ -2.52;▼ -2.70

圖5 不同濃度的GSS362溶液的γ(t)~t-1/2曲線Fig.5 Dynamic surface tension of the GSS362 solutions vs. t-1/2.lgc:■ -2.00;● -2.30;▲ -2.52;▼ -2.70

表3 不同GSS362濃度下體系的擴散系數Table 3 Diffusion coefficients of the GSS362 solutions with different concentration

3 結論

1)GSS362濃度越大,誘導區結束越早(ti越小),GSS362分子在氣/液界面越易達到介平衡(tm越?。?。

2)溫度越高,分子擴散越快,GSS362分子擴散至新鮮表面的速率加快,動態表面張力下降越快,誘導區結束越早(ti越?。珿SS362分子在氣/液界面越易達到介平衡(tm越?。?。

3)在吸附初期,濃度越大,表面活性劑擴散速率越快;在吸附后期,表面活性劑的擴散存在吸附能壘,濃度越大,能壘越大。GSS362在吸附初期為擴散控制吸附,在吸附后期為混合動力控制吸附。

符 號 說 明

c 濃度,mol/L

c0表面活性劑溶液本體溶液濃度,mol/L

cs次表面層中表面活性劑溶液濃度,mol/L

D 分子擴散系數,m2/s

Da表觀擴散系數,m2/s

D*重整化擴散系數,m2/s

n 常數,無因次量

R 氣體常數,J/(mol·K)

R1/2t1/2時表面張力的下降速率,mN/(m·s)

t 吸附時間,s

t*常數,s

ti誘導區結束時間,s

tm介平衡開始時間,s

t1/2體系的表面壓達到介平衡時表面壓二分之一的時間,s

Γ(t) t時刻新鮮表面上表面活性劑的吸附量,mmol/m2

Γeq平衡時的表面吸附量,mmol/m2

γ(t) 表面活性劑溶液在t時刻的動態表面張力,mN/m

γ0溶劑的表面張力,mN/m

γeq平衡表面張力,mN/m

γm介平衡時的表面張力,mN/m

εa吸附能壘,J/mol

τ 虛擬變量

[1] 趙國璽,朱?瑤. 表面活性劑作用原理[M].北京:中國輕工業出版社,2003:138.

[2] 姚志剛,胡艾希,黎小武,等. ELENMOL JS-2水溶液的動態表面張力及吸附動力學[J].華南理工大學學報:自然科學版,2007,35(7):82-87.

[3] 李干佐,牟建海,隋華. 表面活性劑溶液的動態表面張力與吸附動力學研究[J].日用化學工業,1999(5):28-33.

[4] Tang Xiaolan,Dong Jinfeng,Li Xuefeng. A comparison of spreading behaviors of Silwet L-77 on dry and wet lotus leaves[J].J Colloid Interface Sci,2008,325(1):223-227.

[5] 喻紅梅,李建華,華平,等. 1,4-丁二醇雙子琥珀酸聚醚(3)異辛基混合雙酯磺酸鈉的合成[J].精細化工,2013,30(11):1217-1222.

[6] Rosen M J. Geminis:A new generation of surfactants:These materials have better properties than conventional ionic surfactants as well as positive synergistic effects with non-ionics[J]. Chemtech,1993,11(11):887-892.

[7] 戴寶江,華平,黃艷芳. 1,3-丙二醇雙磺基琥珀酸異辛酯鈉的合成及性能研究[J].精細化工中間體,2015,45(3):67-69.

[8] 于艷春,肖鶴鳴. 琥珀酸二油脂磺酸鈉的合成、結構及水合作用[J].物理化學學報,2009,25(1):30-34.

[9] 凌錦龍,曹楓,徐敏虹,等. C16mimBr/Triton X-100混合體系的表面性質和膠團化行為[J].化學學報,2011,69(22):2658-2664.

[10] Rosen M J,Song L D. Dynamic surface tension of aqueous surfactant solutions:8. Effect of spacer on dynamic properties of Gemini surfactant solutions[J].J Colloid Interface Sci,1996,179(1):261-268.

[11] 陳文君,李干佐,柴金嶺,等. 十八烷基二甲基氯化銨水溶液動態表面張力及吸附動力學研究[J].化學學報,2002,60(4):669-673.

[12] 王曉春,張磊,宮清濤,等. 不同結構烷基苯磺酸鈉水溶液的泡沫性能及動態表面張力[J].高等學?;瘜W學報,2007,28(11):2118-2123.

[13] Eastoe J,Dalton J S. Dynamic surface tension and adsorption mechanisms of surfactants at the air-water interface[J].Adv Colloid Interface Sci,2000,85(2/3):103-144.

[14] Ward A F H,Tordai L. Time-dependence of boundary tensions of solutions:Ⅰ. The role of diffusion in time-effects[J].J Chem Phys,1946,14(7):453-461.

[15] Fainerman V B,Makievski A V,Miller R. The analysis of dynamic surface tension of sodium alkyl sulphate solutions,based on asymptotic equations of adsorption kinetic theory[J]. Colloids Surf,A,1994,87(1):61-75.

[16] Liggieri L,Ravera F,Passerone A. A diffusion-based approach to mixed adsorption kinetics[J].Colloids Surf,A,1996,114(16):351-359.

(編輯 王 萍)

Dynamic surface tension and adsorption kinetics of GSS362 aqueous solutions

Yu Hongmei1,Hua Ping1,Zhu Guohua1,Liu Junfeng2,Kong Decun1,Liu Qian1
(1. School of Chemistry and Chemical Engineering,Nantong University,Nantong Jiangsu 226019,China;2. School of Textile and Clothing,Nantong University,Nantong Jiangsu 226019,China)

The dynamic surface tensions of 1,3-propanediol bis-isooctyl sulfo succinate(GSS362) aqueous solutions with different concentration were investigated at different temperature using the maximum bubble pressure method,and characterized by means of the Rosen model. The adsorption mechanism of GSS362 was studied by means of the Ward-Tordai equation. The results showed that,the higher the concentration was,the faster the dynamic surface tension of GSS362 decreased. When the concentration was higher than the critical micelle concentration,the medium equilibrium surface tension was nearly equal to the equilibrium surface tension. Temperature rise was beneficial to the adsorption of GSS362 on the fresh surface. The adsorption was controlled by diffusion in the beginning of the adsorption and controlled by mixed diffusion-kinetics in the later stage.

sulfosuccinate;Gemini surfactant;dynamic surface tension;adsorption mechanism

1000-8144(2017)01-0068-07

TQ 423.1

A

10.3969/j.issn.1000-8144.2017. 01.010

2016-06-25;[修改稿日期]2016-10-26。

喻紅梅(1977—),女,江蘇省如東縣人,碩士,副教授。聯系人:華平,電話 15950805899,電郵 hua.p@ntu.edu.cn。

江蘇省科技計劃項目(BY2015047-12)。

猜你喜歡
擴散系數
石墨烯厚膜熱擴散系數與微觀結構的關系
表觀擴散系數值與肝細胞癌分級的相關性以及相關性與腫瘤大小關系的分析
傳染病信息(2022年6期)2023-01-12 08:57:54
磁共振表觀擴散系數對肝轉移瘤化療療效評估應用
甘肅科技(2020年20期)2020-04-13 00:30:42
一類具有變擴散系數的非局部反應-擴散方程解的爆破分析
表觀擴散系數與乳腺浸潤性導管癌預后因素的相關性分析
磁共振成像(2015年9期)2015-12-26 07:20:32
非腫塊型強化的乳腺癌磁共振成像表觀擴散系數值與HER-2表達的相關性分析
磁共振成像(2015年7期)2015-12-23 08:53:09
非腫塊型乳腺癌的MR表觀擴散系數及腫瘤大小與Ki-67表達的相關性研究
腫瘤影像學(2015年3期)2015-12-09 02:38:51
基于Sauer-Freise 方法的Co- Mn 體系fcc 相互擴散系數的研究
上海金屬(2015年5期)2015-11-29 01:13:59
FCC Ni-Cu 及Ni-Mn 合金互擴散系數測定
上海金屬(2015年6期)2015-11-29 01:09:09
非時齊擴散模型中擴散系數的局部估計
主站蜘蛛池模板: 国产精品视频白浆免费视频| 亚洲色中色| 国产美女免费网站| 特级毛片8级毛片免费观看| 国产免费黄| 人人看人人鲁狠狠高清| 国产一区免费在线观看| 国产精品美女免费视频大全 | 美女扒开下面流白浆在线试听| 国产精品嫩草影院视频| 午夜福利网址| 天天综合网色中文字幕| 婷婷开心中文字幕| 99久久99这里只有免费的精品| 午夜性爽视频男人的天堂| 东京热高清无码精品| 亚洲精品自拍区在线观看| 亚洲欧美成aⅴ人在线观看| 国产成人高清亚洲一区久久| 六月婷婷激情综合| 亚洲AV无码乱码在线观看裸奔 | 人妻丰满熟妇αv无码| 好久久免费视频高清| 亚洲最猛黑人xxxx黑人猛交 | julia中文字幕久久亚洲| 国产经典免费播放视频| 激情五月婷婷综合网| 久久天天躁狠狠躁夜夜2020一| 欧美日本在线| 午夜日b视频| 亚洲精品视频网| 国产va在线观看免费| 国产无码精品在线播放| 国产成人综合久久精品下载| 99热这里只有免费国产精品 | 成年人国产网站| 幺女国产一级毛片| 国产高清色视频免费看的网址| 亚洲天堂2014| 成年人午夜免费视频| 精品国产成人a在线观看| 国产精品私拍99pans大尺度 | 亚洲精品欧美重口| 国产亚洲欧美在线视频| a在线亚洲男人的天堂试看| 久久青草免费91观看| 国产亚洲现在一区二区中文| 欧美综合区自拍亚洲综合绿色 | 免费看a毛片| 精品91自产拍在线| 日本人妻丰满熟妇区| 午夜一区二区三区| 在线精品自拍| 国产乱人伦偷精品视频AAA| 国产又爽又黄无遮挡免费观看 | 中文字幕无码电影| 亚洲AV人人澡人人双人| 日韩在线视频网站| 五月综合色婷婷| 国产av无码日韩av无码网站| 久久久久久久97| 婷婷亚洲视频| 中文字幕免费播放| 免费高清a毛片| 波多野结衣中文字幕久久| 免费观看成人久久网免费观看| 欧美亚洲欧美| 波多野一区| 在线观看无码a∨| 国产欧美在线观看一区| 欧美日韩第三页| 国产精品网址在线观看你懂的| 免费毛片全部不收费的| 亚洲精品日产精品乱码不卡| 亚洲综合二区| 亚洲国产中文欧美在线人成大黄瓜 | 亚洲制服丝袜第一页| 国产丝袜精品| 国产一区二区三区精品久久呦| 四虎永久免费地址| 91成人在线免费视频| 中文字幕无码中文字幕有码在线|