• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    玻璃鋼與含CO2流體界面模型與MS模擬

    2016-08-01 02:18:25張東興王冠輝王世剛陳雨時(shí)
    材料科學(xué)與工藝 2016年2期
    關(guān)鍵詞:玻璃鋼環(huán)氧樹脂

    張東興,張 兵,2,王冠輝,王世剛,賈 近,陳雨時(shí)

    ?

    玻璃鋼與含CO2流體界面模型與MS模擬

    張東興1,張兵1,2,王冠輝1,王世剛3,賈近1,陳雨時(shí)1

    (1.哈爾濱工業(yè)大學(xué)材料科學(xué)與工程學(xué)院,哈爾濱150001;2.斯達(dá)玻璃鋼有限公司,哈爾濱150001;3.中國石油吉林油田公司勘察設(shè)計(jì)院,吉林松原138000)

    摘要:在CO2驅(qū)油項(xiàng)目中,由于CO2是以超臨界狀態(tài)與水交替注入井下,局部形成了具有腐蝕性的CO2和水的混合物.為研究玻璃鋼(GFRP)與這種混合物的界面間相互作用,本文采用Material Studio(MS)模擬軟件建立了環(huán)氧樹脂與鏈烴、水、CO2、H2CO34種純凈物流體之間的界面模型,以及樹脂與鏈烴/水/CO2、鏈烴/H2CO32種混合流體間的界面模型,并通過分子動(dòng)力學(xué)計(jì)算得到了6種界面模型的最低能量結(jié)構(gòu)和界面能.結(jié)果表明,鏈烴與樹脂界面幾乎不吸附,CO2、水與樹脂界面吸附較差,H2CO3與樹脂界面能較高,能與樹脂中的—NH2、—OH基團(tuán)形成氫鍵具有較好的吸附作用,即H2CO3是混合體系中與環(huán)氧樹脂作用的主要物質(zhì).

    關(guān)鍵詞:玻璃鋼;環(huán)氧樹脂;分子模擬;界面模型;CO2

    CO2驅(qū)油是一種高效環(huán)保的驅(qū)采油技術(shù)[1],不但可以大幅降低稠油粘度提高石油采收率,還可以填埋大量工業(yè)廢氣減少溫室氣體的排放.玻璃鋼(GFRP)由于具有良好的力學(xué)性能和耐腐蝕性能,在石油化工領(lǐng)域廣泛應(yīng)用.在CO2驅(qū)油項(xiàng)目中,由于CO2是以超臨界狀態(tài)與水交替注入井下的[2-3],局部形成了油、水和CO2共同存在的混合物,而這類混合物具有一定腐蝕性,因而研究GFRP與混合物的界面間相互作用具有重要意義.

    超臨界CO2具有優(yōu)良的溶解性能[4],可將稠油溶解并與水形成低粘度混合液,增加了原油的流動(dòng)性,從而提高低滲油田的采收率,具有非常廣闊的應(yīng)用前景.在采出液中,同時(shí)存在原油、水和CO2,而CO2水溶液對(duì)金屬管材具有一定的腐蝕性,使得輸油管道壽命大大低于設(shè)計(jì)壽命[5-6].采用纖維纏繞成型工藝制造的玻璃鋼管道具有良好的力學(xué)性能和水利特性,因具有比金屬管道更好的耐腐蝕性能而在石油化工行業(yè)得以快速推廣和應(yīng)用,但是玻璃鋼不是絕對(duì)惰性材料,在使用中也存在著腐蝕現(xiàn)象,但國內(nèi)外對(duì)此進(jìn)行的研究很少,其作用機(jī)理還不明確.

    由于玻璃鋼樹脂基體結(jié)構(gòu)以及采出液成分的復(fù)雜性,實(shí)驗(yàn)研究其腐蝕機(jī)理具有一定的局限性,而利用分子模擬技術(shù)則可以較快獲得基體與不同組分的相互作用特征.分子模擬技術(shù)不僅可以模擬分子的靜態(tài)結(jié)構(gòu),也可以模擬分子體系的動(dòng)態(tài)行為(如氫鍵的締合與解締、吸附、擴(kuò)散等)[7-8].

    本文以分子模擬作為手段,通過Accelrys公司的Material Studio(MS)建立玻璃鋼的表面樹脂及鏈烴、水、CO2、H2CO3、鏈烴/水/CO2、鏈烴/ H2CO36種流體結(jié)構(gòu)模型,并對(duì)其進(jìn)行分子力學(xué)優(yōu)化,從而建立樹脂與幾種流體的界面模型并進(jìn)行分子動(dòng)力學(xué)計(jì)算,進(jìn)而得到各種界面模型的最低能量結(jié)構(gòu),探討含CO2流體不同組分與玻璃鋼基體的相互作用機(jī)理.

    1 環(huán)氧樹脂界面模型建立與優(yōu)化

    由于玻璃鋼管材表面的環(huán)氧樹脂對(duì)內(nèi)部結(jié)構(gòu)具有良好的包覆性,而分子模擬的研究以分子尺度相互作用為主,故忽略環(huán)氧樹脂缺陷及界面的影響,將整體界面簡化為石油采出液與環(huán)氧樹脂的界面模型.其中,與采出液接觸的管材內(nèi)表面成分為胺類固化雙酚A型環(huán)氧樹脂,固化劑為間苯二胺.采用Accelrys公司的Material Studio(MS)軟件中的“Amorphous cell”模塊,建立三維網(wǎng)狀結(jié)構(gòu)的胺類固化環(huán)氧樹脂的模型,如圖1所示.模型的晶胞尺寸為:a=3.5 nm,b=3.6 nm,c=5.6 nm,α=β=γ=90°.

    圖1 三維網(wǎng)狀結(jié)構(gòu)的環(huán)氧樹脂晶胞模型

    初始模型的結(jié)構(gòu)優(yōu)化:首先采用“smart minimizer”方法對(duì)模型進(jìn)行能量最小化優(yōu)化,匯聚水平設(shè)為costomized,最大迭代次數(shù)設(shè)為10 000次.其中“牛頓方法”的迭代原子數(shù)設(shè)為200個(gè),模擬時(shí)的能量最大偏移為默認(rèn)值5 000 kcal/mol,熱力學(xué)狀態(tài)描述選擇Berendsen方法.然后,在正則系綜(NVT)下對(duì)初步優(yōu)化后的結(jié)構(gòu)進(jìn)行分子動(dòng)力學(xué)計(jì)算,模擬溫度定為598 K,模擬時(shí)間定為50萬步,計(jì)算的時(shí)間步長為1 fs,總共的模擬時(shí)間為500 ps,保存軌跡選擇全部的運(yùn)動(dòng)及坐標(biāo)參數(shù),同時(shí)設(shè)定每5 000步輸出一個(gè)模擬結(jié)果.模擬計(jì)算中的長程相互作用力的截?cái)嗑嚯x為1.0 nm,在分子力學(xué)及分子動(dòng)力學(xué)模擬中采用的力場是Compass力場.將動(dòng)力學(xué)計(jì)算得到的環(huán)氧樹脂穩(wěn)定結(jié)構(gòu)作為樹脂/流體的界面模型的起始表面模型.

    2 環(huán)氧/流體界面模型與分子模擬

    石油采出液是一種由水、氣和各種烴類有機(jī)物組成的復(fù)雜混合物,其中烴類有機(jī)物的主要成分為烷烴.本文為方便計(jì)算將含CO2的采出液成分簡化為液態(tài)鏈烴、水、CO2及H2CO34種物質(zhì),而液態(tài)鏈烴則選取碳原子數(shù)最多、且具有多支鏈結(jié)構(gòu)的2-甲基-3-乙基戊烷進(jìn)行計(jì)算.采用MS軟件分別建立鏈烴、水、CO2、H2CO3、鏈烴/水/CO2混合流體、鏈烴/H2CO3混合流體6種流體的周期性邊界模型.

    首先對(duì)初始界面模型進(jìn)行分子力學(xué)計(jì)算,然后對(duì)界面模型進(jìn)行模擬退火計(jì)算.模擬退火的目的是使分子能夠跨越較高的勢壘,找到分子的最低能量構(gòu)型[9-10].模擬退火法的具體過程為:初始溫度10 K,終止溫度1 010 K,每2 000步的分子動(dòng)力學(xué)模擬,溫度相應(yīng)變化40 K,循環(huán)次數(shù)為20.每種模型均采用多個(gè)不同的初始結(jié)構(gòu),分別進(jìn)行退火模擬.將模擬退火法得到的數(shù)百個(gè)較低能量的構(gòu)象中,根據(jù)界面作用能式(1),找出其中能量最低的結(jié)構(gòu).然后,利用分子力學(xué)方法將此結(jié)構(gòu)的模型體系進(jìn)行優(yōu)化.經(jīng)過上述處理后的模型體系基本已處于平衡狀態(tài),最后在298 K下采用NVT系綜分子動(dòng)力學(xué)方法對(duì)此體系進(jìn)行計(jì)算.界面相互作用能公式為[11]

    式中:EInteraction為界面作用能;Etotal為玻璃鋼表面與流體界面體系的總能量;Esurface為玻璃鋼表面的能量(不含流體);Efluent為流體的能量(不含玻璃鋼表面).

    3 結(jié)果與討論

    3.1 環(huán)氧與4種純凈物的界面模型結(jié)構(gòu)分析

    固體表面與內(nèi)部的結(jié)構(gòu)和化學(xué)組成是不完全相同的.環(huán)氧樹脂基體內(nèi)部原子、分子間相互作用并達(dá)到一個(gè)平衡,而表面上的分子與原子則有自發(fā)降低表面能的傾向[12].由于固化后的環(huán)氧樹脂無法通過收縮來降低表面能,只能通過吸附其他物質(zhì)分子降低界面張力,從而降低表面能.因此,固體表面對(duì)其他物質(zhì)的吸附強(qiáng)弱在一定程度上也反映出兩者間發(fā)生反應(yīng)的難易程度.

    采用MS模擬得到了鏈烴、水、CO2和H2CO34種物質(zhì)與樹脂界面吸附模型的最低能量結(jié)構(gòu),如圖2所示.鏈烴分子在界面處分布無序而且距離樹脂界面的較遠(yuǎn),樹脂/鏈烴鏈界面模型最低能量結(jié)構(gòu)的界面能幾乎為0,屬于一種非常弱的相互作用界面,因此兩者之間發(fā)生反應(yīng)的幾率很低,與很多文獻(xiàn)和資料中飽和烷烴鏈和玻璃鋼不發(fā)生反應(yīng)的結(jié)論是一致的[13-15].樹脂/水、樹脂/CO2界面模型的最低能量結(jié)構(gòu)比較類似,2種分子都均勻聚集在距離樹脂界面不遠(yuǎn)的位置,分子間分布比較均勻,這說明小分子物質(zhì)與樹脂表面相互吸附,卻又沒有產(chǎn)生較強(qiáng)的作用而相互接近或穿插,屬于一種比較弱的相互作用.

    H2CO3是CO2與水生成的一種不穩(wěn)定弱酸,存在2個(gè)很活潑的—OH,容易與環(huán)氧樹脂中的有機(jī)基團(tuán)相互作用.圖2(d)為樹脂/H2CO3界面模型的最低能量結(jié)構(gòu)圖,可以看出,以分子形式存在的H2CO3與環(huán)氧樹脂界面出現(xiàn)了相互的穿插,H2CO3分子深入到樹脂內(nèi)部,H2CO3結(jié)構(gòu)中不穩(wěn)定的—OH與樹脂表面的—NH2及—OH基團(tuán)產(chǎn)生氫鍵(如圖3中淺藍(lán)色虛線所示),使樹脂表面對(duì)H2CO3吸附能力增強(qiáng),從而形成了樹脂/H2CO3界面間的較強(qiáng)相互作用.氫鍵的產(chǎn)生,不但可以提高混合物的酸性,促進(jìn)酯類的水解,而且還影響苯環(huán)上的親核取代反應(yīng),甚至可以在一些化學(xué)合成中充當(dāng)催化劑[16],因此對(duì)樹脂的化學(xué)穩(wěn)定性有一定的影響.

    圖2 界面模型最低能量結(jié)構(gòu)圖

    圖3 樹脂/H2CO3界面間形成的氫鍵

    3.2 環(huán)氧與2種混合物的界面模型結(jié)構(gòu)分析

    由于玻璃鋼管材輸送介質(zhì)是以鏈烴、水、CO2為主要成分的混合物,因此建立樹脂/(鏈烴/水/ CO2)和樹脂/(鏈烴/H2CO3)2種界面模型更具有實(shí)際意義,如圖4所示.

    圖4 樹脂/(鏈烴/水/CO2)、樹脂/(鏈烴/H2CO3)界面模型最低能量結(jié)構(gòu)圖

    由圖4可以看出,樹脂/(鏈烴/水/CO2)界面中的鏈烴分子依然無序地存在于遠(yuǎn)離界面的區(qū)域,而CO2和H2O分子距離環(huán)氧樹脂距離則與其各自單組分界面模型相似,之間沒有很強(qiáng)的吸附.另外,由于鏈烴分子處于固-液界面的外側(cè),因此也沒有表現(xiàn)出對(duì)界面吸附的明顯阻礙作用.樹脂/(鏈烴/H2CO3)界面最低能量結(jié)構(gòu)與圖2中的樹脂/H2CO3界面最低能量結(jié)構(gòu)類似,H2CO3與環(huán)氧樹脂界面出現(xiàn)了相互的穿插,H2CO3分子深入到樹脂內(nèi)部,鏈烴分子依然處于固-液界面的外側(cè).這說明在混合體系下,環(huán)氧樹脂與混合物的強(qiáng)相互作用主要是由H2CO3起作用,其他的幾種組分則作用不明顯,這也是實(shí)際應(yīng)用中H2CO3環(huán)境下易造成腐蝕的主要原因.

    3.3 6種界面結(jié)構(gòu)的界面能對(duì)比

    樹脂與鏈烴、水、CO2、H2CO3、鏈烴/水/CO2、鏈烴/H2CO36種流體按照1~6的順序表示,其界面結(jié)構(gòu)最低能量模型的界面能如圖5所示.

    由圖5可以看出,6種流體與樹脂的界面能相差很大.樹脂與鏈烴的界面能非常小,基本可以忽略它們之間的相互作用;樹脂/水和樹脂/CO2界面能相差不大,接近-300 kcal/mol,這也表明水、CO2對(duì)樹脂表面的吸附性較弱.相比鏈烴、水和CO23種物質(zhì),H2CO3與樹脂界面的界面能最高,接近-2 500 kcal/mol,這表明H2CO3極易吸附在樹脂表面,與樹脂表面官能團(tuán)相作用,并可能產(chǎn)生進(jìn)一步的化學(xué)反應(yīng).樹脂/(鏈烴/水/CO2)模型的界面能介于2種單組分界面模型的界面能之間,說明這2種組分之間對(duì)界面吸附?jīng)]有明顯的相互促進(jìn)或者抑制作用.

    樹脂/(鏈烴/H2CO3)界面能遠(yuǎn)遠(yuǎn)高于樹脂/ H2CO3的界面能,這是因?yàn)殒湡N阻礙了樹脂表面對(duì)H2CO3的吸附行為有關(guān).另外,相比樹脂/(鏈烴/水/CO2)界面,樹脂/(鏈烴/H2CO3)結(jié)構(gòu)界面能更高,證明了H2CO3具有比H2O與CO2混合物更好的界面活性,對(duì)環(huán)氧樹脂具有更大的腐蝕性.

    綜上所述,環(huán)氧樹脂與4種流體的相互作用具有較大的差異性,H2CO3易與—OH、—NH形成界面間的相互作用,從而進(jìn)一步與樹脂發(fā)生反應(yīng)形成腐蝕.

    4 結(jié) 論

    1)鏈烴與樹脂界面幾乎不吸附,CO2、水單組份與樹脂界面吸附較差,CO2和水的混合物對(duì)樹脂界面的吸附介于2種單一組分之間,沒有明顯的相互促進(jìn)和抑制作用.

    2)H2CO3與樹脂表面可以形成較好的吸附結(jié)構(gòu),并能與樹脂中的—NH2、—OH基團(tuán)形成氫鍵,這也表明H2CO3是混合體系中與環(huán)氧樹脂作用的主要物質(zhì).

    3)鏈烴的存在使樹脂與混合介質(zhì)的界面能降低,對(duì)其他組分有一定的阻礙作用,即鏈烴阻礙了樹脂表面對(duì)H2CO3的吸附行為.

    參考文獻(xiàn):

    [1] 喬紅軍,任曉娟,閆鳳平,等.低滲透儲(chǔ)層水氣交替注入方式室內(nèi)試驗(yàn)研究[J].石油天然氣學(xué)報(bào),2013,35(7):114-117.QIAO Hongjun,REN Xiaojuan,YAN Fengping,et al. Experimental study of low permeability reservoir with water alternate gas injection method[J].Journal of Oil and Gas Technology,2013,35(7):114-117.

    [2] 高云叢,趙密福,王建波,等.特低滲油藏CO2非混相驅(qū)生產(chǎn)特征與氣竄規(guī)律[J].石油勘探與開發(fā),2014,41(1):79-85. GAO Yuncong,ZHAO Mifu,WANG Jianbo,et al. Performance and gas breakthrough during CO2immiscible flooding in ultra-low permeability reservoirs[J]. Petroleum Exploration and Development,2014,41 (1):79-85.

    [3] QAMAR M M,ISLAM M R.CO2Injection in the weyburn field of canada:optimization of enhanced oil recovery and greenhouse gas storage with horizontal wells [C]//SPE/DOE Improved Oil Recovery Symposium. Oklahoma:[s.n.],2000:1-17.

    [4] LIU B,WONG-FOY A G,MATZGER A J.Rapid and enhanced activation of microporous coordination polymers by flowing supercritical CO2[J].Chem Commun,2013,49(14):1419-1421.

    [5] 李天成.榆濟(jì)天然氣管道內(nèi)腐蝕原因[J].腐蝕與防護(hù),2014,35(3):282-285. LI Tiancheng.Internal corrosion causes of Yu-Ji natural gas transmission pipeline[J].Corrosion and Protection,2014,35(3):282-285.

    [6] 王世杰.原油含水率對(duì)油氣管材超臨界CO2腐蝕行為的影響[J].腐蝕科學(xué)與防護(hù)技術(shù),2015,1(1): 73-77. WANG Shijie.Effect of water content of crude oil on corrosion behavior of oil and gas tubular goods steel in supercritical CO2system[J].Corrosion Science and Protection Technology,2015,1(1):73-77.

    [7] NADLER R,SREBNIK S.Molecularsimulation of polyamide synthesis by interfacial polymerization[J]. Journal of Membrane Science,2008,315(1/2): 100-105.

    [8] YANG Q,WEN X P,SUN C M,et al.Molecular simulation of adsorption mechanism of silicone rubber interfacial systems[J].Advanced Research on Energy Chemistry Materials&Information III,2014,1023: 15-18.

    [9] 崔守鑫,胡海泉,肖效光,等.分子動(dòng)力學(xué)模擬基本原理和主要技術(shù)[J].聊城大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,18(1):30-34. CUI Shouxin,HU Haiquan,XIAO Xiaoguang,et al. The basic principles and methods of molecular dynamics simulation[J].Journal of Liaocheng University (Nat Sci),2005,18(1):30-34.

    [10]樊康旗,賈建援.經(jīng)典分子動(dòng)力學(xué)模擬的主要技術(shù)[J].微納電子技術(shù),2005(3):133-138. FAN Kangqi,JIA Jianyuan.An overview on classical molecular dynamics simulation[J].MEMS Device and Technology,2005(3):133-138.

    [11]江浩,岳紅,趙建勇,等.復(fù)合材料界面相互作用的分子模擬研究與進(jìn)展[J].中國膠粘劑,2012 (7):50-54. JIANG Hao,YUE Hong,ZHAO Jianyong,et al.Research and progress on molecular simulation of interfacial interaction in composite[J].China Adhesives,2012(7):50-54.

    [12]周靜,易勇,胡麗,等.氫及其同位素在聚苯乙烯與Cu界面的擴(kuò)散[J].強(qiáng)激光與粒子束,2013,25 (5):1179-1183. ZHOU Jing,YI Yong,HU Li,et al.Diffusion of hydrogen and its isotopes in polystyrene-Cu interface[J]. High Power Laser and Particle Beams,2013,25(5): 1179-1183.

    [13]HUANG G U.Tensile behaviour of some high performance filaments after sulphuric acid treatment[J]. Materials and Design,2008,29(8):1671-1674.

    [14]HIROYUKI K,RIVASTAVA V K.The effect of an acidic stress environment on the stress-intensity factor for GRP laminates[J].Composites Science and Technology,2001,61:1109-1114.

    [15]FARSHAD M,NECOLA A.Strain corrosion of glass fiber-reinforced plastic pipe[J].Polymer Testing,2004(2):517-521.

    (編輯 程利冬)

    中圖分類號(hào):TB332

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1005-0299(2016)02-0058-05

    doi:10.11951/j.issn.1005-0299.20160207

    收稿日期:2015-05-22.

    基金項(xiàng)目:國家科技重大專項(xiàng)(2011ZX05016-004).

    作者簡介:張東興(1961—),男,教授,博士生導(dǎo)師.

    通信作者:張東興,E-mail:dongxingzhang@163.com.

    The interfacial structure and MS simulation of GFRP in transfering fluid containing CO2

    ZHANG Dongxing1,ZHANG Bing1,2,WANG Guanhui1,WANG Shigang3,JIA Jin1,CHEN Yushi1

    (1.School of Materials Science and Engineering,Harbin Institute of Technology,Harbin 150001,China;2.Star Fiberglass Harbin Co.LTD,Harbin 150088,China;3.Survey&Design Institute,China National Petroleum Corporation Jilin Oilfield,Songyuan 138000,China)

    Abstract:It is significant to focus on the corrosion of GFRP in CO2flooding project because the corrosive mixture of CO2and water are formed when CO2is injected underground alternatively with water at supercritical state.The interfacial models between GFRP and hydrocarbon,H2O,CO2,H2CO3,hydrocarbons/H2O/CO2and hydrocarbons/H2CO3were established by Material Studio to investigate the effects of interface.The interfacial energy and minimum energy structure for six kinds of interface models are obtained through the molecular dynamics simulation.The results show that the hydrocarbon is merely adsorbed by epoxy resin,and CO2and H2O are absorbed at very limited level,respectively.However the adsorption between the resin and H2CO3is remarkable compared with other fluent,and hydrogen bond was formed between—NH2and—OH.It can be deduced that the H2CO3is the main substance act with amine cured GFRP in the mixture.

    Keywords:GFRP;epoxy;molecular simulation;interfacial model;CO2

    猜你喜歡
    玻璃鋼環(huán)氧樹脂
    環(huán)氧樹脂基復(fù)合材料加筋板結(jié)構(gòu)吸濕行為研究
    YLZ3500玻璃鋼人防隔斷門數(shù)控液壓機(jī)的研制
    裝配式玻璃鋼閥門井的研發(fā)及應(yīng)用
    煤氣與熱力(2021年3期)2021-06-09 06:16:18
    大型纏繞成型玻璃鋼容器的優(yōu)化設(shè)計(jì)
    河北中復(fù)玻璃鋼有限公司
    中國釀造(2019年9期)2019-10-08 05:44:04
    環(huán)氧樹脂灌漿技術(shù)在EPP項(xiàng)目中的應(yīng)用
    粉末涂料用環(huán)氧樹脂的合成研究
    上海建材(2017年2期)2017-07-21 14:02:10
    碳納米管陣列/環(huán)氧樹脂的導(dǎo)熱導(dǎo)電性能
    玻璃鋼快艇隨動(dòng)液壓舵“跑舵”問題淺析
    可膨脹石墨對(duì)環(huán)氧樹脂的阻燃改性
    午夜福利视频1000在线观看| 欧美日韩综合久久久久久 | 国产精品久久久久久人妻精品电影| 麻豆成人午夜福利视频| 久久久久久久久免费视频了| 一二三四在线观看免费中文在| 无遮挡黄片免费观看| 久久久久久国产a免费观看| 男人的好看免费观看在线视频| 亚洲男人的天堂狠狠| 欧美乱色亚洲激情| 丁香欧美五月| 岛国在线免费视频观看| 成年免费大片在线观看| 欧美日韩国产亚洲二区| 精品人妻1区二区| 亚洲精品粉嫩美女一区| 女警被强在线播放| 国产伦精品一区二区三区四那| 精品国产亚洲在线| 精品乱码久久久久久99久播| 午夜福利视频1000在线观看| 淫妇啪啪啪对白视频| 欧美不卡视频在线免费观看| 婷婷六月久久综合丁香| 性色avwww在线观看| 欧美中文综合在线视频| 女同久久另类99精品国产91| 亚洲精品色激情综合| 一边摸一边抽搐一进一小说| 久久人妻av系列| 黄色日韩在线| 久久中文看片网| 欧美成人一区二区免费高清观看 | 国产高清视频在线观看网站| 欧美成狂野欧美在线观看| 天堂影院成人在线观看| 欧美日韩黄片免| 国产私拍福利视频在线观看| 一本久久中文字幕| 高潮久久久久久久久久久不卡| 夜夜看夜夜爽夜夜摸| 1024香蕉在线观看| 国产欧美日韩一区二区精品| 啦啦啦韩国在线观看视频| 欧美三级亚洲精品| 深夜精品福利| 毛片女人毛片| 国产1区2区3区精品| 一进一出抽搐gif免费好疼| www.999成人在线观看| 禁无遮挡网站| 高潮久久久久久久久久久不卡| 不卡一级毛片| 深夜精品福利| 国产美女午夜福利| 亚洲成av人片在线播放无| 国产精品影院久久| 亚洲av第一区精品v没综合| 欧美黄色淫秽网站| 国产精品1区2区在线观看.| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片| www.熟女人妻精品国产| 在线观看免费视频日本深夜| АⅤ资源中文在线天堂| 国内精品久久久久久久电影| 99国产精品一区二区三区| 伊人久久大香线蕉亚洲五| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 亚洲午夜理论影院| 搡老岳熟女国产| 亚洲欧美一区二区三区黑人| 黄色丝袜av网址大全| 又紧又爽又黄一区二区| 国产精品免费一区二区三区在线| 在线观看免费午夜福利视频| 日日摸夜夜添夜夜添小说| 18禁观看日本| 观看美女的网站| 国产v大片淫在线免费观看| 亚洲欧美日韩东京热| 国产又色又爽无遮挡免费看| 男插女下体视频免费在线播放| 淫妇啪啪啪对白视频| 国产高清视频在线播放一区| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人免费电影在线观看| 99在线人妻在线中文字幕| 在线观看免费午夜福利视频| 香蕉久久夜色| 久久婷婷人人爽人人干人人爱| 欧美xxxx黑人xx丫x性爽| 成人18禁在线播放| 999久久久精品免费观看国产| 亚洲国产欧洲综合997久久,| 色视频www国产| 免费在线观看视频国产中文字幕亚洲| avwww免费| www.www免费av| 欧美一级a爱片免费观看看| 不卡一级毛片| 午夜亚洲福利在线播放| 日本三级黄在线观看| 日本成人三级电影网站| 琪琪午夜伦伦电影理论片6080| 免费在线观看成人毛片| 亚洲精品一卡2卡三卡4卡5卡| 精品国产美女av久久久久小说| 毛片女人毛片| 国产伦一二天堂av在线观看| 黄片小视频在线播放| 后天国语完整版免费观看| 国产1区2区3区精品| 午夜影院日韩av| 欧美zozozo另类| 国产精品永久免费网站| 日本免费一区二区三区高清不卡| 99国产综合亚洲精品| 成年女人毛片免费观看观看9| 色老头精品视频在线观看| 国产免费av片在线观看野外av| 国产成人影院久久av| 亚洲在线观看片| 大型黄色视频在线免费观看| 69av精品久久久久久| 男人的好看免费观看在线视频| 午夜久久久久精精品| 人妻夜夜爽99麻豆av| 欧美性猛交黑人性爽| 亚洲一区高清亚洲精品| 亚洲专区国产一区二区| 国产精品爽爽va在线观看网站| 国产精品亚洲av一区麻豆| 国产毛片a区久久久久| 色综合亚洲欧美另类图片| 欧美乱妇无乱码| 麻豆一二三区av精品| 午夜福利18| 国产精品久久久久久人妻精品电影| 日本在线视频免费播放| 大型黄色视频在线免费观看| 精品久久久久久久末码| 人妻丰满熟妇av一区二区三区| 18禁观看日本| 亚洲国产精品成人综合色| 久久人妻av系列| 又粗又爽又猛毛片免费看| 午夜亚洲福利在线播放| 午夜精品一区二区三区免费看| 久久精品人妻少妇| 国产精品久久视频播放| 久久精品国产综合久久久| 搡老岳熟女国产| 俺也久久电影网| 国产1区2区3区精品| 美女免费视频网站| 99在线视频只有这里精品首页| 黑人操中国人逼视频| 午夜福利欧美成人| 亚洲欧美日韩高清专用| 日韩精品中文字幕看吧| 国产午夜精品久久久久久| 亚洲一区高清亚洲精品| 麻豆久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 欧美日韩亚洲国产一区二区在线观看| 不卡一级毛片| 欧美激情在线99| 久久人妻av系列| 国产精品影院久久| 国产成人影院久久av| 成人精品一区二区免费| av欧美777| 久久久国产成人免费| 欧美绝顶高潮抽搐喷水| 九色成人免费人妻av| 国产97色在线日韩免费| 成年女人看的毛片在线观看| 制服丝袜大香蕉在线| 99久久成人亚洲精品观看| ponron亚洲| 男女之事视频高清在线观看| 天堂网av新在线| 午夜影院日韩av| 日韩高清综合在线| 亚洲专区字幕在线| 亚洲精品美女久久av网站| 精品久久蜜臀av无| 国产高潮美女av| 午夜亚洲福利在线播放| 国产激情偷乱视频一区二区| 久久精品国产清高在天天线| 欧美色欧美亚洲另类二区| 网址你懂的国产日韩在线| 床上黄色一级片| 久久久国产成人精品二区| 黑人巨大精品欧美一区二区mp4| 亚洲精品粉嫩美女一区| 欧美中文综合在线视频| 久久人人精品亚洲av| 在线观看午夜福利视频| 欧美日本视频| 波多野结衣高清作品| 国产亚洲欧美98| 日韩中文字幕欧美一区二区| 日日干狠狠操夜夜爽| 精品久久久久久久久久免费视频| 日韩 欧美 亚洲 中文字幕| 后天国语完整版免费观看| 超碰成人久久| 综合色av麻豆| www.熟女人妻精品国产| 一级a爱片免费观看的视频| 久久午夜亚洲精品久久| 手机成人av网站| 久久久水蜜桃国产精品网| 特大巨黑吊av在线直播| 99国产极品粉嫩在线观看| 天天躁日日操中文字幕| 在线播放国产精品三级| av天堂中文字幕网| 少妇人妻一区二区三区视频| 偷拍熟女少妇极品色| 99久久精品国产亚洲精品| 亚洲av成人一区二区三| 亚洲av成人不卡在线观看播放网| 欧美日韩中文字幕国产精品一区二区三区| 操出白浆在线播放| 亚洲av中文字字幕乱码综合| 日本黄大片高清| 国产蜜桃级精品一区二区三区| 黄频高清免费视频| 老司机福利观看| 国产精品98久久久久久宅男小说| 丰满的人妻完整版| 美女cb高潮喷水在线观看 | 精品国产三级普通话版| 亚洲黑人精品在线| 99久久成人亚洲精品观看| 午夜两性在线视频| 中文字幕人妻丝袜一区二区| 麻豆成人午夜福利视频| 小蜜桃在线观看免费完整版高清| 国产男靠女视频免费网站| 亚洲九九香蕉| 成熟少妇高潮喷水视频| 免费av不卡在线播放| 在线永久观看黄色视频| 午夜精品在线福利| 精品不卡国产一区二区三区| av视频在线观看入口| 99国产精品99久久久久| 男女床上黄色一级片免费看| 级片在线观看| 午夜a级毛片| 国产精品美女特级片免费视频播放器 | 一区福利在线观看| 亚洲第一欧美日韩一区二区三区| 欧美一区二区精品小视频在线| 最近最新中文字幕大全免费视频| 国产亚洲精品av在线| 丰满人妻一区二区三区视频av | 可以在线观看的亚洲视频| 老熟妇乱子伦视频在线观看| 精品免费久久久久久久清纯| 日韩欧美在线乱码| 国产精品日韩av在线免费观看| 级片在线观看| 日本 av在线| 非洲黑人性xxxx精品又粗又长| 国产视频一区二区在线看| 91av网站免费观看| 99久久99久久久精品蜜桃| 色吧在线观看| 欧美极品一区二区三区四区| 黄色视频,在线免费观看| 国产精品一区二区免费欧美| 亚洲激情在线av| 99久国产av精品| 久久久久久九九精品二区国产| 男女视频在线观看网站免费| 欧美乱码精品一区二区三区| 淫秽高清视频在线观看| 看免费av毛片| avwww免费| 好男人在线观看高清免费视频| 最近最新中文字幕大全电影3| 视频区欧美日本亚洲| 色视频www国产| 黄色日韩在线| 久久久久久久久久黄片| 好男人在线观看高清免费视频| 在线十欧美十亚洲十日本专区| 最新中文字幕久久久久 | 很黄的视频免费| 成人永久免费在线观看视频| 成人av一区二区三区在线看| 在线免费观看的www视频| 黄片大片在线免费观看| 国产av麻豆久久久久久久| 午夜福利欧美成人| 国模一区二区三区四区视频 | 久久久久久大精品| 国产97色在线日韩免费| 精品不卡国产一区二区三区| 宅男免费午夜| 69av精品久久久久久| 久久久久久久久免费视频了| 久久久国产欧美日韩av| 亚洲欧美一区二区三区黑人| 99久久精品国产亚洲精品| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 这个男人来自地球电影免费观看| 伦理电影免费视频| 在线观看免费视频日本深夜| 欧美午夜高清在线| x7x7x7水蜜桃| 熟妇人妻久久中文字幕3abv| 国产精品影院久久| 天堂动漫精品| 熟女少妇亚洲综合色aaa.| 岛国在线观看网站| 美女免费视频网站| 天天躁日日操中文字幕| 色av中文字幕| 欧美一区二区精品小视频在线| 久久久精品欧美日韩精品| 日本黄色片子视频| 日韩国内少妇激情av| av天堂在线播放| 最近最新免费中文字幕在线| 最近最新中文字幕大全电影3| 日本黄色视频三级网站网址| 老熟妇乱子伦视频在线观看| 男人舔女人的私密视频| 听说在线观看完整版免费高清| 欧美一级毛片孕妇| 天天一区二区日本电影三级| 国产真人三级小视频在线观看| 嫁个100分男人电影在线观看| 午夜福利在线观看吧| 男人和女人高潮做爰伦理| 亚洲在线观看片| 三级国产精品欧美在线观看 | 一进一出抽搐动态| 国产成人av教育| 99在线视频只有这里精品首页| 久久久久久久久中文| 亚洲av成人不卡在线观看播放网| 淫秽高清视频在线观看| 成人亚洲精品av一区二区| 婷婷精品国产亚洲av| 18禁国产床啪视频网站| 一级作爱视频免费观看| 国产av一区在线观看免费| 午夜精品久久久久久毛片777| 99热这里只有精品一区 | 欧美3d第一页| 亚洲人成网站高清观看| 欧美日韩综合久久久久久 | 女人高潮潮喷娇喘18禁视频| 一个人免费在线观看电影 | 日本 av在线| 国产欧美日韩一区二区精品| 国产成人精品久久二区二区免费| 日韩欧美精品v在线| 亚洲精品久久国产高清桃花| 国产欧美日韩精品一区二区| 村上凉子中文字幕在线| 日本 欧美在线| 亚洲欧美日韩高清在线视频| 午夜福利18| 桃红色精品国产亚洲av| 国产午夜精品论理片| 九九久久精品国产亚洲av麻豆 | 波多野结衣巨乳人妻| 一二三四社区在线视频社区8| 人妻丰满熟妇av一区二区三区| 午夜免费观看网址| 日韩欧美国产在线观看| 日韩欧美一区二区三区在线观看| 黑人欧美特级aaaaaa片| 99热只有精品国产| 亚洲人成网站高清观看| 哪里可以看免费的av片| 国产精品野战在线观看| 琪琪午夜伦伦电影理论片6080| 免费观看精品视频网站| 久久久久久国产a免费观看| 少妇的逼水好多| 女生性感内裤真人,穿戴方法视频| 99热这里只有是精品50| 一二三四社区在线视频社区8| 亚洲在线观看片| 在线免费观看的www视频| 九色国产91popny在线| 中文字幕人妻丝袜一区二区| 中文字幕久久专区| 床上黄色一级片| 看片在线看免费视频| 露出奶头的视频| 色综合婷婷激情| 午夜福利免费观看在线| 国产成+人综合+亚洲专区| 欧美性猛交╳xxx乱大交人| 日韩欧美在线乱码| 成在线人永久免费视频| 成人国产综合亚洲| 久久天躁狠狠躁夜夜2o2o| 欧美日韩精品网址| 一本久久中文字幕| 久久久久久九九精品二区国产| av女优亚洲男人天堂 | 中亚洲国语对白在线视频| 香蕉久久夜色| 51午夜福利影视在线观看| 国产精品乱码一区二三区的特点| 麻豆av在线久日| 国产精品1区2区在线观看.| 国产精品 国内视频| 久久国产精品影院| 一本久久中文字幕| 日本成人三级电影网站| 18禁裸乳无遮挡免费网站照片| 身体一侧抽搐| 一进一出抽搐动态| 精品人妻1区二区| 又爽又黄无遮挡网站| 少妇丰满av| 色av中文字幕| 日韩欧美一区二区三区在线观看| 一级毛片女人18水好多| 亚洲熟妇熟女久久| 色综合婷婷激情| 性色av乱码一区二区三区2| 成人国产综合亚洲| 18美女黄网站色大片免费观看| 一区二区三区激情视频| 女警被强在线播放| 免费看a级黄色片| 成人亚洲精品av一区二区| 久久精品国产亚洲av香蕉五月| 黄色丝袜av网址大全| 久久久久久久久久黄片| 色综合婷婷激情| 香蕉丝袜av| 国产毛片a区久久久久| 日韩欧美 国产精品| 性色av乱码一区二区三区2| 国产精品乱码一区二三区的特点| 久久天躁狠狠躁夜夜2o2o| 亚洲 欧美 日韩 在线 免费| 在线观看日韩欧美| 2021天堂中文幕一二区在线观| 欧美极品一区二区三区四区| 丰满的人妻完整版| 国产成人av激情在线播放| 综合色av麻豆| 国产免费男女视频| 老鸭窝网址在线观看| 午夜福利在线观看免费完整高清在 | 97人妻精品一区二区三区麻豆| 亚洲精品456在线播放app | 国产高清三级在线| 长腿黑丝高跟| 老司机深夜福利视频在线观看| 夜夜爽天天搞| 日韩中文字幕欧美一区二区| 国内精品久久久久久久电影| 国产真人三级小视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品免费一区二区三区在线| 国产激情偷乱视频一区二区| 久久中文字幕一级| 亚洲专区字幕在线| 露出奶头的视频| 国内久久婷婷六月综合欲色啪| 成人亚洲精品av一区二区| 此物有八面人人有两片| 亚洲精品久久国产高清桃花| 国产亚洲精品av在线| 国产精品免费一区二区三区在线| 成年女人毛片免费观看观看9| 日韩精品青青久久久久久| 黄片大片在线免费观看| 白带黄色成豆腐渣| 日本a在线网址| 亚洲av片天天在线观看| 老汉色av国产亚洲站长工具| 国产久久久一区二区三区| 最好的美女福利视频网| 给我免费播放毛片高清在线观看| 黄色日韩在线| 国产黄色小视频在线观看| 国产伦精品一区二区三区四那| 18禁国产床啪视频网站| 精品久久久久久成人av| 亚洲自拍偷在线| 美女 人体艺术 gogo| 人妻夜夜爽99麻豆av| 天天添夜夜摸| 成人欧美大片| 88av欧美| 夜夜躁狠狠躁天天躁| 久久中文看片网| 亚洲欧洲精品一区二区精品久久久| 国产免费av片在线观看野外av| 色尼玛亚洲综合影院| 在线观看一区二区三区| 色av中文字幕| 久久这里只有精品中国| 亚洲avbb在线观看| 99re在线观看精品视频| 熟女少妇亚洲综合色aaa.| 五月伊人婷婷丁香| 久久久久国产一级毛片高清牌| 免费看a级黄色片| 欧美黄色淫秽网站| 久久久国产欧美日韩av| 久久久国产成人精品二区| 无遮挡黄片免费观看| 成熟少妇高潮喷水视频| 在线视频色国产色| 欧美黄色淫秽网站| 97人妻精品一区二区三区麻豆| 成人鲁丝片一二三区免费| 日本与韩国留学比较| 精品无人区乱码1区二区| av天堂在线播放| aaaaa片日本免费| 国产男靠女视频免费网站| 精品一区二区三区四区五区乱码| 女人高潮潮喷娇喘18禁视频| 97超级碰碰碰精品色视频在线观看| 久久久久亚洲av毛片大全| 最新在线观看一区二区三区| av天堂在线播放| 日本黄色片子视频| 亚洲美女黄片视频| 特大巨黑吊av在线直播| 国产精品香港三级国产av潘金莲| www.精华液| 一进一出好大好爽视频| 高清在线国产一区| 久久精品国产亚洲av香蕉五月| 久久久久精品国产欧美久久久| 婷婷六月久久综合丁香| 高清在线国产一区| 99精品欧美一区二区三区四区| 淫妇啪啪啪对白视频| 亚洲熟妇中文字幕五十中出| 国产激情欧美一区二区| 久久久国产精品麻豆| 特级一级黄色大片| 色播亚洲综合网| 精品不卡国产一区二区三区| 观看美女的网站| 99精品久久久久人妻精品| 在线观看66精品国产| 欧美xxxx黑人xx丫x性爽| 国产99白浆流出| 非洲黑人性xxxx精品又粗又长| 97人妻精品一区二区三区麻豆| 美女扒开内裤让男人捅视频| 欧美黄色片欧美黄色片| 天堂av国产一区二区熟女人妻| 制服丝袜大香蕉在线| 国产乱人伦免费视频| 午夜福利在线观看吧| 日本成人三级电影网站| 精品欧美国产一区二区三| 波多野结衣高清无吗| 国产v大片淫在线免费观看| 久久99热这里只有精品18| 亚洲 欧美一区二区三区| av在线天堂中文字幕| 岛国在线观看网站| 久久九九热精品免费| 亚洲 国产 在线| 成人亚洲精品av一区二区| 国产高清激情床上av| 人妻丰满熟妇av一区二区三区| 国产精品女同一区二区软件 | 搞女人的毛片| 欧美成人性av电影在线观看| 亚洲午夜精品一区,二区,三区| 亚洲精华国产精华精| 国产精品久久久久久久电影 | 欧美3d第一页| 国产私拍福利视频在线观看| 日韩大尺度精品在线看网址| 视频区欧美日本亚洲| 一个人看视频在线观看www免费 | 国产一区二区在线观看日韩 | 少妇丰满av| 国产精品亚洲av一区麻豆| 日韩精品青青久久久久久| 国产精华一区二区三区| 久久香蕉精品热| 麻豆av在线久日| 一级黄色大片毛片| 日本 欧美在线| 欧美黑人巨大hd| 熟妇人妻久久中文字幕3abv| 噜噜噜噜噜久久久久久91| 中文字幕av在线有码专区| 国产一区二区激情短视频| 真人做人爱边吃奶动态| 久久精品综合一区二区三区| 看免费av毛片| 亚洲精品粉嫩美女一区|