亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        水合物聚集體受力分析及臨界聚集體粒徑的計(jì)算

        2014-06-07 05:57:08劉海紅李玉星王武昌
        石油化工 2014年1期
        關(guān)鍵詞:模型

        劉海紅,李玉星,王武昌

        (中國(guó)石油大學(xué)(華東) 儲(chǔ)運(yùn)與建筑工程學(xué)院,山東 青島 266580)

        水合物聚集體受力分析及臨界聚集體粒徑的計(jì)算

        劉海紅,李玉星,王武昌

        (中國(guó)石油大學(xué)(華東) 儲(chǔ)運(yùn)與建筑工程學(xué)院,山東 青島 266580)

        分析了油包水體系水合物聚集體的受力特性,作用在水合物聚集體上的力有范德華力、液橋力、固橋力、靜電力、碰撞力、剪切力、重力等,通過(guò)受力比較確定了水合物聚集體在流動(dòng)體系中主要的聚集力和分離力分別是液橋力和剪切力。用液橋力代替聚集力,對(duì)水合物臨界聚集體粒徑的計(jì)算模型進(jìn)行了改進(jìn)。驗(yàn)證結(jié)果表明,在合理選取或準(zhǔn)確測(cè)量水合物聚集體表面接觸角的前提下,用液橋力代替聚集力是合理的,改進(jìn)的水合物臨界聚集體粒徑的計(jì)算模型能準(zhǔn)確計(jì)算水合物聚集體的臨界粒徑,該模型能為水合物漿的安全流動(dòng)提供技術(shù)指導(dǎo),促進(jìn)水合物漿液技術(shù)的工程應(yīng)用。

        油包水體系;水合物;聚集力;臨界聚集體粒徑

        近年來(lái)隨著深海油氣田的開(kāi)發(fā),水合物問(wèn)題日益突出,從20世紀(jì)90年代,動(dòng)力學(xué)抑制劑和防聚劑的開(kāi)發(fā),水合物冷流技術(shù)成為水合物防治的研究重點(diǎn)[1],實(shí)現(xiàn)在不抑制水合物生成的前提下保證水合物顆粒分散在油相中安全輸送。水合物顆粒聚集是影響水合物漿液穩(wěn)定輸送的關(guān)鍵因素,但目前對(duì)水合物聚集過(guò)程的受力特性和聚集機(jī)理尚不明確,亟需開(kāi)展針對(duì)水合物顆粒受力特性和聚集過(guò)程的研究[1-2]。國(guó)內(nèi)外學(xué)者通?;谑芰Ψ治鼋⒕蹐F(tuán)受力平衡模型來(lái)預(yù)測(cè)臨界聚集體粒徑[3-6],作為工業(yè)生產(chǎn)判斷顆粒聚集的重要指標(biāo), Camargo等[7]提出了適用于水合物聚集體的受力平衡模型,該模型是基于水合物聚集體所受聚集力和分離力的平衡,并結(jié)合水合物特有的疏松多孔特性建立的水合物臨界聚集體粒徑計(jì)算模型,在模型中提到了水合物聚集體間的聚集力,但沒(méi)有給出聚集力的明確計(jì)算式或數(shù)值。目前國(guó)內(nèi)外學(xué)者對(duì)聚集力的確定大都是基于微機(jī)械測(cè)力設(shè)備直接測(cè)定水合物聚集體間的黏附力[8-11],而從理論計(jì)算角度分析水合物顆粒受力的研究工作較少。

        本工作通過(guò)分析水合物聚集體的受力情況,比較各種力的大小確定水合物臨界聚集體粒徑計(jì)算模型中的主要聚集力,以使該模型能更方便的計(jì)算水合物臨界聚集體的粒徑,為油氣水管道的安全流動(dòng)提供指導(dǎo)。

        1 油包水體系水合物聚集體上的主要作用力

        油包水體系中水合物聚集體受重力、浮力、范德華力、靜電力、液橋力、固橋力、碰撞力以及流體剪切力等力的作用,根據(jù)各作用力在水合物聚集過(guò)程中發(fā)揮的效果不同可以將作用力大致分為3類[12]:促使水合物聚集的力為范德華力、液橋力、固橋力、靜電力;促使水合物分離的力為碰撞力、流體剪切力;促使水合物聚集體旋轉(zhuǎn)的力為滑動(dòng)摩擦力和滾動(dòng)摩擦力等。

        只考慮水合物聚集體對(duì)心碰撞的情況,不考慮聚集體碰撞后的旋轉(zhuǎn),水合物聚集體受力的示意圖見(jiàn)圖1。

        圖1 水合物聚集體受力的示意圖Fig.1 Schematic diagram of forces between hydrate agglomerates.Fvw:Van der Waals force;Fcap:liquid bridge force;Fs:solid bridgeforce;Fe:electrostatic force;Fv:buoyant force;Fshear:shear force;Fc:collision force;G:gravity.

        1.1 范德華力

        水合物聚集體Ⅰ和聚集體Ⅱ之間的范德華力可由式(1)計(jì)算。

        1.2 液橋力

        水合物表面具有強(qiáng)親水性,當(dāng)兩個(gè)聚集體相互靠近時(shí)在兩者之間會(huì)形成液橋(液橋的輪廓見(jiàn)圖2),從而產(chǎn)生液橋力改變水合物的受力特性,水合物聚集體間液橋力由式(2)計(jì)算。

        圖2 水合物聚集體間液橋的輪廓Fig.2 Outline of liquid bridge between hydrate agglomerates.

        1.3 固橋力

        基于水合物顆粒接觸誘導(dǎo)聚集機(jī)理[13],水滴和水合物顆粒或聚集體接觸,隨著水合反應(yīng)的進(jìn)行,水滴就會(huì)固化而轉(zhuǎn)化為水合物顆粒,假如水滴存在于兩個(gè)或多個(gè)水合物顆?;蚓奂w之間,水滴變形吸附在顆?;蚓奂w表面上并固化形成連接多個(gè)水合物顆?;蚓奂w的固橋,也就產(chǎn)生了固橋力,固橋力沒(méi)有固定的計(jì)算公式,一般通過(guò)實(shí)驗(yàn)得到。通常固橋形成后在沒(méi)有強(qiáng)外力(強(qiáng)攪拌等)的作用下,固橋是不會(huì)破裂的,因此在這里并不把固橋力作為導(dǎo)致顆粒發(fā)生聚集的一種作用力,而是把固橋視為水合物聚集體的一部分,一旦固橋形成,就認(rèn)為水合物顆粒發(fā)生聚集或聚集體體積增大形成更大塊的聚集體,在其他力的計(jì)算時(shí)用新聚集體的粒徑代替形成固橋前的粒徑。

        1.4 靜電力

        水合物聚集體在運(yùn)動(dòng)中相互間摩擦、碰撞以及和壁面的摩擦、碰撞都會(huì)導(dǎo)致聚集體表面產(chǎn)生靜電力,但由于體系中存在導(dǎo)電性能良好的水,因此水合物聚集體表面的靜電力可以迅速消除[14],所以水合物聚集過(guò)程中不考慮靜電力的作用。

        1.5 碰撞力

        普通固體聚集體發(fā)生碰撞時(shí),碰撞力是判斷兩聚集體碰撞后聚集與否的關(guān)鍵,但水合物聚集體碰撞過(guò)程卻表現(xiàn)出了不同的特性,這主要是源于水合物表面的強(qiáng)親水性,聚集體間的液橋阻礙了兩者的碰撞變形,傳統(tǒng)意義上的物理碰撞變形作用力在水合物聚集體碰撞中不存在。

        1.6 流體剪切力

        剪切力正比于粒子與周圍流體間的局部速率差,其計(jì)算方法見(jiàn)式(3)。

        1.7 凈重力

        水合物聚集體的凈重力是聚集體的重力與浮力之差,由式(4)計(jì)算。

        2 主要聚集力和分離力的確定

        2.1 基本參數(shù)的選取

        基于水合物聚集體間各作用力的定義,參照文獻(xiàn)[15]選擇以下基本參數(shù)取值,作為油包水體系水合物聚集力和分離力的計(jì)算基礎(chǔ)。

        水合物聚集體在油相中的Hamaker常數(shù)A取值為5×10-21J;油水界面張力(0.02~0.04 N/m[16])取值為0.035 N/m;油相黏度μ0取值為0.135 Pa·s。參考甲烷水合物晶體的密度(910 kg/m3)[17],對(duì)水合物聚集體的密度ρa(bǔ)取值為800 kg/m3。

        2.2 油包水體系水合物聚集力和分離力的計(jì)算分析

        水合物聚集體上4種作用力的計(jì)算結(jié)果見(jiàn)圖3。從圖3可看出,水合物聚集體粒徑在幾到幾百微米的范圍內(nèi),液橋力遠(yuǎn)大于范德華力和凈重力,Turner等[18]在油包水體系利用聚焦光束反射測(cè)量?jī)x觀測(cè)發(fā)現(xiàn)分散相水滴直接轉(zhuǎn)化成初始水合物顆粒,轉(zhuǎn)化后水合物顆粒粒徑和分散相水滴粒徑一致,而在油包水體系中形成的初始分散水滴的直徑一般都在幾微米的數(shù)量級(jí),即使是水滴形成的水合物顆粒發(fā)生聚集,聚集體粒徑也不過(guò)幾百微米,由于在該粒徑范圍內(nèi)液橋力比其他聚集力大得多,所以水合物聚集體間的主要聚集力是液橋力;分離力則比較容易判斷,剪切力是主要考慮的分離力。

        圖3 水合物聚集體上4種作用力的計(jì)算結(jié)果Fig.3 Calculation results of four forces between hydrate agglomerates.■ Liquid brige force;● Van der Waals force;▲ Net gravity;▲Shear force

        3 水合物臨界聚集體粒徑計(jì)算模型的改進(jìn)

        Camargo等[7]基于水合物聚集體受力平衡,提出水合物聚集體粒徑計(jì)算方法,見(jiàn)式(5)。

        考慮到水合物聚集體間的聚集力中液橋力遠(yuǎn)大于范德華力,因此,將式(5)中的水合物聚集體間的聚集力用液橋力代替,得到式(6)。

        式(6)即為改進(jìn)的水合物臨界聚集體粒徑的計(jì)算模型。由式(6)計(jì)算得到的dA實(shí)際上就是水合物最大聚集體粒徑dAmax。初始形成的水合物聚集體較小,此時(shí)dA<dAmax,水合物聚集體受到的聚集力大于分離力,聚集體會(huì)繼續(xù)增大,直到水合物聚集體dA接近或等于dAmax,水合物的聚集過(guò)程才會(huì)停止。在水合物聚集過(guò)程中,若水合物漿液在管道內(nèi)能保持正常流動(dòng),將不會(huì)發(fā)生水合物堵塞事故。

        根據(jù)式(6),只要已知水合物體積分?jǐn)?shù)(φ)、分形維數(shù)(f)、初始水合物顆粒直徑(dp)、油相黏度(μ0)、剪切速率(γ)、油水界面張力(γll)、接觸角(θ),就可以計(jì)算得到相應(yīng)的最大聚集粒徑dAmax。但如果計(jì)算得到的dAmax<dp,只需取dAmax=dp。

        4 水合物臨界聚集體粒徑計(jì)算模型的驗(yàn)證

        為了驗(yàn)證改進(jìn)的水合物臨界聚集體粒徑的計(jì)算模型(見(jiàn)式6),參照文獻(xiàn)[11]對(duì)相關(guān)參數(shù)取值:f=2.5,dp=1.5 μm,μ0=0.06 Pa·s,φ=0.3,其他參數(shù)見(jiàn)表1。

        表1 水合物懸浮液環(huán)道實(shí)驗(yàn)數(shù)據(jù)[11]Table 1 Experimental data for hydrate suspension in a loop[11]

        水合物懸浮液黏度與水合物聚集體粒徑的關(guān)系見(jiàn)式(7)。

        由式(7)計(jì)算得到水合物聚集體的粒徑,進(jìn)而根據(jù)式(6)計(jì)算得到聚集體間的聚集力,計(jì)算結(jié)果見(jiàn)表2。

        表2 改進(jìn)模型中水合物聚集體間聚集力的計(jì)算結(jié)果Table 2 Agglomeration force between hydrate agglomerates calculated by the improved model

        根據(jù)式(1)和式(2)分別計(jì)算得到的范德華力和接觸角分別為0°和54°時(shí)的液橋力見(jiàn)表3。

        對(duì)比表2和表3可看出,由式(2)計(jì)算得到的液橋力數(shù)值更接近式(6)計(jì)算得到的聚集力數(shù)值,兩力的數(shù)量級(jí)均約為10-7N,而式(1)計(jì)算得到的范德華力數(shù)量級(jí)在10-9N左右,其值遠(yuǎn)小于聚集力,大約是聚集力的1/20。對(duì)該計(jì)算有兩點(diǎn)說(shuō)明:第一,計(jì)算范德華力時(shí)分離距離取值1 nm,但實(shí)際上由于水合物表面的粗糙性及水合物表面易吸附油中瀝青質(zhì)等特點(diǎn)使得分離距離大于1 nm,根據(jù)式(1)可知水合物聚集體間的范德華力數(shù)量級(jí)必定小于10-9N,進(jìn)一步說(shuō)明水合物聚集體間的范德華力不是導(dǎo)致水合物聚集的主要聚集力;第二,表3中接觸角取值0°時(shí)計(jì)算得到的液橋力是式(6)計(jì)算得到的聚集力的兩倍左右,實(shí)際上環(huán)道實(shí)驗(yàn)的油中含有芘乙酸,它的存在會(huì)改變水合物表面的潤(rùn)濕特性,使得聚集體表面的接觸角增大,這里無(wú)法測(cè)量接觸角的實(shí)際值,但若將接觸角取值54°時(shí)計(jì)算得到的水合物聚集體間的液橋力數(shù)值與模型計(jì)算的聚集力數(shù)值一致,因此通過(guò)上面的分析可以說(shuō)明改進(jìn)的水合物臨界聚集體粒徑的計(jì)算模型是正確的。

        水合物臨界聚集體粒徑是保證管道內(nèi)流體穩(wěn)定流動(dòng)的最大水合物聚集體極限粒徑,改進(jìn)的水合物臨界聚集體粒徑的計(jì)算模型可為工業(yè)生產(chǎn)提供指導(dǎo)。

        表3 水合物聚集體間范德華力和液橋力的計(jì)算結(jié)果Table 3 Calculated Van der Waals force and liquid bridge force between hydrate agglomerates

        5 結(jié)論

        1)在流動(dòng)體系中,水合物聚集體間聚集力主要有液橋力、范德華力,分離力主要是剪切力。水合物聚集過(guò)程起主導(dǎo)作用的聚集力是液橋力,分離力是剪切力。

        2)用液橋力代替聚集力對(duì)水合物臨界聚集體粒徑的計(jì)算模型進(jìn)行改進(jìn)。改進(jìn)的水合物臨界聚集體粒徑的計(jì)算模型可以計(jì)算相應(yīng)工況下水合物臨界聚集體粒徑。

        3)水合物表面接觸角的正確選擇和準(zhǔn)確測(cè)量對(duì)驗(yàn)證模型是非常重要的,也是應(yīng)用改進(jìn)模型準(zhǔn)確計(jì)算油包水體系水合物臨界聚集體粒徑的關(guān)鍵。

        符 號(hào) 說(shuō) 明

        A Hamaker常數(shù),J

        dA水合物聚集體粒徑,m

        dAmax水合物聚集體最大臨界粒徑,m

        dp初始水合物顆粒的直徑,m

        f 分形維數(shù),無(wú)因次

        Fa水合物聚集體間聚集力,N

        Fc水合物聚集體間碰撞力,N

        Fcap水合物聚集體間液橋力,N

        Fe水合物聚集體間靜電力,N

        Fs水合物聚集體間固橋力,N

        Fshear水合物聚集體間剪切力,N

        Fv水合物聚集體間浮力,N

        Fvw水合物聚集體間范德華力,N

        G 水合物聚集體的凈重力,N

        g 重力加速度,m/s2

        Ra水合物聚集體半徑,m

        S 水合物聚集體表面間距,m

        γ 剪切速率,1/s

        γll油水界面張力,N/m

        θ 接觸角,°

        μ 水合物懸浮液黏度,Pa·s

        μ0分散相流體黏度,Pa·s

        ρa(bǔ),ρf水合物聚集體及連續(xù)相的密度 ,kg/m3

        φ 水合物體積分?jǐn)?shù),%

        φmax水合物聚集體最大的填充率,%

        [1] Lund A. Comments to Some Preliminary Results from the Exxon Hydrate Flow Loop[J]. Ann N Y Acad Sci,1994,715:447 -449.

        [2] Zhou Tao,Li Hongzhong,Shinohara K. Agglomerating Fluidization of Group C Particles:Major Factors of Coalescence and Breakup of Agglomerates[J]. Adv Powder Technol,2006,17(2):159 - 166.

        [3] Zhou Tao,Li Hongzhong. Force Balance Modelling for Agglomerating Fluidization of Cohesive Particles[J]. Powder Technol,2000,111(1/2):60 - 65.

        [4] Iwadate Y,Horio M. Prediction of Agglomerate Sizes in Bubbling Fluidized Beds of Group C Powders[J]. Powder Technol,1998,100(2/3):223 - 236.

        [5] 蘇吉,朱慶山. 聲場(chǎng)流化床中超細(xì)顆粒聚團(tuán)受力與尺寸[J].過(guò)程工程學(xué)報(bào),2010,10(3):431 - 437.

        [6] 周濤,李洪鐘. 粘性顆粒流化床中聚團(tuán)大小的計(jì)算模型[J].化學(xué)反應(yīng)工程與工藝,1999,15(1):44 - 50.

        [7] Camargo R,Palermo T. Rheological Properties of Hydrate Suspensions in an Asphaltenic Crude Oil[C]//Proceedings of the Fourth International Conference on Gas Hydrates,Yokohama,2002:880 - 885.

        [8] Yang S O,Kleehammer D M,Huo Z,et al. Temperature Dependence of Particle-Particle Adherence Forces in Ice and Clathrate Hydrates[J]. J Colloid Interface Sci,2004,277(2):335 - 341.

        [9] Taylor C J,Dieker L E,Miller K T,et al. Micromechanical Adhesion Force Measurements Between Tetrahydrofuran Hydrate Particles[J]. J Colloid Interface Sci,2007,306(2):255 - 261.

        [10] Dieker L E,Aman Z M,George N C,et al. Micromechanical Adhesion Force Measurements Between Hydrate Particles in Hydrocarbon Oils and Their Modifcations[J]. Energy Fuels,2009,23(12):5966 - 5971.

        [11] Aman Z M,Aspenes G,Sloan E D,et al. The Effect of Chemistry and System Conditions on Hydrate Interparticle Adhesion Forces Toward Aggregation and Hydrate Plug Formation[C]//SPE International Symposium on Oilfeld Chemistry,Texas:Society of Petroleum Engineers,2011.

        [12] 張文斌,祁海鷹,由長(zhǎng)福,等. 碰撞誘發(fā)顆粒團(tuán)聚及破碎的力學(xué)分析[J]. 清華大學(xué)學(xué)報(bào):自然科學(xué)版,2002,42(12):1639 - 1643.

        [13] Fidel-Dufour A,Gruy F,Herri J M. Rheology of Methane Hydrate Slurries During Their Crystallization in a Water in Dodecane Emulsion Under Flowing[J]. Chem Eng Sci,2006,61(2):505 - 515.

        [14] 鄭建祥. 粘附性顆粒動(dòng)理學(xué)及氣固兩相流動(dòng)力特性的研究[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2008.

        [15] Anklam M R,York D,Helmerich L,et al. Effects of Antiagglomerants on the Interactions Between Hydrate Particles[J]. AIChE J,2007,54(2):565 - 574.

        [16] 楊會(huì)麗,王業(yè)飛,任熵,等. 原油/水界面張力的影響因素[J]. 承德石油高等??茖W(xué)校學(xué)報(bào),2006,8(1):1 - 3.

        [17] 陳光進(jìn),孫長(zhǎng)宇,馬慶蘭. 氣體水合物科學(xué)與技術(shù)[M]. 北京:化學(xué)工業(yè)出版社,2007:37 - 38.

        [18] Turner D J,Miller K T,Sloan E D. Direct Conversion of Water Droplets to Methane Hydrate in Crude Oil[J]. Chem Eng Sci,2009,64(23):5066 - 5072.

        (編輯 李治泉)

        ·技術(shù)動(dòng)態(tài)·

        石科院裂解汽油苯乙烯抽提成套新工藝技術(shù)通過(guò)鑒定

        石油化工科學(xué)研究院等單位承擔(dān)的“裂解汽油抽提蒸餾回收苯乙烯成套技術(shù)(STED)開(kāi)發(fā)”項(xiàng)目通過(guò)鑒定。采用STED成套工藝技術(shù),在廣東新華粵石化股份有限公司建成了30 kt/a苯乙烯抽提裝置,實(shí)現(xiàn)了該技術(shù)的首次工業(yè)應(yīng)用,并實(shí)現(xiàn)了長(zhǎng)周期穩(wěn)定運(yùn)行。運(yùn)行結(jié)果表明,各項(xiàng)技術(shù)指標(biāo)均達(dá)到了設(shè)計(jì)要求,產(chǎn)品質(zhì)量好、回收率高,能耗、物耗低。標(biāo)定考核結(jié)果表明,苯乙烯純度達(dá)到99.87%(w),鉑-鈷色號(hào)小于10,抽提回收率98%,總回收率達(dá)94.3%,主要技術(shù)指標(biāo)顯著優(yōu)于國(guó)外同類技術(shù)。STED工藝技術(shù)已申請(qǐng)發(fā)明專利10件,獲授權(quán)3件,具有創(chuàng)新性和自主知識(shí)產(chǎn)權(quán)。

        Analysis of Forces Between Hydrate Agglomerates and Calculation of Critical Agglomerate Size

        Liu Haihong,Li Yuxing,Wang Wuchang
        (College of Pipeline and Civil Engineering,China University of Petroleum,Qingdao Shandong 266580,China)

        Forces between hydrate agglomerates in a water-in-oil system,including Van der Waals force,liquid bridge force,solid bridge force,electrostatic force,buoyant force,shear force,collision force and gravity,were analyzed. The results showed that the dominant agglomeration force and separating force between the hydrate agglomerates in the fow system were the liquid bridge force and shear force,respectively. Based on the force analysis,the calculating model for the critical hydrate agglomerate size was improved by substituting the liquid bridge force for the agglomeration force. The results indicated that the substitution was reasonable under the condition of correctly choosing or accurately measuring the surface contact angle of the hydrate agglomerate. It was showed that using the improved model could calculate the critical diameter of the hydrate agglomerate accurately.

        water-in-oil system;hydrate;agglomeration force;critical agglomerate size

        1000 - 8144(2014)01 - 0046 - 05

        TE 88

        A

        2013 - 06 - 30;[修改稿日期] 2013 - 09 - 28。

        劉海紅(1989—),女,山東省東營(yíng)市人,碩士生,電話 13465867764,電郵 upcliuhaihong@163.com。聯(lián)系人:李玉星,電郵 lyxstar@126.com。

        國(guó)家自然科學(xué)基金項(xiàng)目(51006120);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20110133110004);中國(guó)石油大學(xué)(華東)研究生創(chuàng)新工程資助項(xiàng)目(CX-1241);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助(11CX04048A)。

        猜你喜歡
        模型
        一半模型
        一種去中心化的域名服務(wù)本地化模型
        適用于BDS-3 PPP的隨機(jī)模型
        提煉模型 突破難點(diǎn)
        函數(shù)模型及應(yīng)用
        p150Glued在帕金森病模型中的表達(dá)及分布
        函數(shù)模型及應(yīng)用
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
        3D打印中的模型分割與打包
        日本一区二区三区免费播放| 亚洲人av毛片一区二区| 久久综合九色综合久久久| 亚洲精品人成中文毛片| 国产特级毛片aaaaaa| 成人区视频| 手机在线观看亚洲av| 亚洲AV色欲色欲WWW| 国产丝袜美腿在线视频| 奶头又大又白喷奶水av| 国产人在线成免费视频| 91精品91久久久久久| 黄色潮片三级三级三级免费| 老色鬼在线精品视频| 亚洲精品无码久久久久sm| 亚洲日本国产乱码va在线观看| 深夜黄色刺激影片在线免费观看| 国产国产人免费人成免费视频| 好男人视频在线视频| 亚洲国产视频精品一区二区| 国产在线观看视频一区二区三区| 亚洲av永久无码精品网站在线观看| 毛片网站视频| 欧美精品久久久久久三级| 暴露的熟女好爽好爽好爽| 欧洲成人一区二区三区| 国产精品福利视频一区| 欧亚精品无码永久免费视频| 五月激情四射开心久久久| 人妻中文字幕无码系列| 亚洲自拍另类制服在线| 一区二区三区在线视频爽| 久草青青91在线播放| 人妻无码一区二区三区免费| 日本丶国产丶欧美色综合| 亚洲国产精品国自产拍av在线| 东北女人一级内射黄片| 蜜臀av 国内精品久久久| 青青青爽国产在线视频| 精品女同一区二区三区在线播放器 | 国产免费a∨片在线软件|