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

    天然氣脫蠟旋風(fēng)分離器分離效率的模擬

    2014-07-05 16:03:04袁惠新方毅付雙成
    化工進(jìn)展 2014年1期
    關(guān)鍵詞:切向速度段長度旋風(fēng)

    袁惠新,方毅,付雙成

    (常州大學(xué)機(jī)械工程學(xué)院,江蘇 常州 213016)

    天然氣脫蠟旋風(fēng)分離器分離效率的模擬

    袁惠新,方毅,付雙成

    (常州大學(xué)機(jī)械工程學(xué)院,江蘇 常州 213016)

    旋風(fēng)分離器分離效率高,不易堵塞,用于天然氣脫蠟效果顯著。通過CFD軟件Fluent模擬CYG-S型天然氣脫蠟旋風(fēng)分離器的兩相流場,得到了旋風(fēng)分離器內(nèi)的壓力、切向速度、軸向速度分布。對比了不同入口速度下的模擬與理論計算的分割粒徑x50,發(fā)現(xiàn)具有很好的吻合度,兩相模擬有一定的可靠性。結(jié)果表明:在旋風(fēng)分離器錐段底部靠近壁面處的石蠟液滴質(zhì)量濃度較高;隨著進(jìn)口流量的增加,旋風(fēng)分離器分離效率提高,當(dāng)進(jìn)口流量為1000 m3/h時,x50可以達(dá)到5.3 μm;大粒徑液滴的分離效果明顯,但在所研究的進(jìn)口流量范圍內(nèi),進(jìn)口流量的變化不能明顯地影響粒徑小于5 μm液滴的分離效率;柱段和錐段長度的增加使得旋風(fēng)分離器的整體長度增加,延長了液滴在旋風(fēng)分離器內(nèi)的停留時間,提高了旋風(fēng)分離器的分離效率。

    旋風(fēng)分離器;數(shù)值模擬;分離效率;進(jìn)口流量;柱段和錐段長度

    石蠟是由正構(gòu)烷烴、帶支鏈(或異構(gòu))烷烴和環(huán)烷烴的組成的復(fù)雜混合物,有很多品種的熔點(diǎn)范圍在常溫范圍(2~36 ℃)[1],以固態(tài)和半固態(tài)的的形式存在于天然氣中。由于石蠟具有黏結(jié)和流動性差的特性,極易在管道內(nèi)壁沉積,甚至堵塞設(shè)備,所以對于天然氣中石蠟的分離,傳統(tǒng)的過濾分離會非常容易堵塞過濾器,嚴(yán)重影響生產(chǎn)效率,存在安全隱患。

    袁惠新[2]發(fā)明的天然氣凈化裝置,對天然氣加熱至石蠟熔點(diǎn)以上,使石蠟以液滴狀態(tài)存在于天然氣中,再利用旋風(fēng)分離器進(jìn)行氣液分離,具有分離效率高、處理能力大、投資小、操作簡便可靠等優(yōu)點(diǎn),可有效防止堵塞的發(fā)生。目前該設(shè)備已在江漢油田投產(chǎn)應(yīng)用,效果顯著。

    為了進(jìn)一步優(yōu)化天然氣旋流脫蠟過程,考察旋流分離器內(nèi)兩相流動特性,推動該技術(shù)發(fā)展,本文作者選用CYG-S型天然氣脫蠟旋風(fēng)分離器,使用CFD軟件Fluent模擬旋風(fēng)分離器的內(nèi)部流場,并對不同入口流量、不同柱段和錐段長度下的分離效率進(jìn)行了模擬。

    1 模擬方法

    隨著計算機(jī)以及算法的不斷發(fā)展,數(shù)值模擬方法已成為重要的研究手段。目前常用的幾種湍流計算模型有標(biāo)準(zhǔn)k-ε模型[3]、RNG模型、雷諾應(yīng)力模型(RSM) 以及代數(shù)應(yīng)力模型(ASM)。RSM可以計算出雷諾應(yīng)力分量,因而更適用于強(qiáng)旋流旋風(fēng)分離器流場[4],已有不少學(xué)者采用了該模型并取得了較好的效果:王海剛等[5]通過對比3種模型的計算結(jié)果與實驗數(shù)據(jù),發(fā)現(xiàn)以RSM模型的預(yù)報結(jié)果最為合理。葛坡等[6]利用RSM模型,使用SIMPLEC算法,對對稱多入口型旋風(fēng)分離器進(jìn)行數(shù)值模擬,得到了不錯的模擬效果。

    當(dāng)石蠟液滴的體積分?jǐn)?shù)低于5%時,可采用離散相模型(DPM)模擬石蠟液滴的分離過程。將分散相的石蠟液滴看做離散顆粒,求解每一個離散顆粒的運(yùn)動方程;連續(xù)相運(yùn)動則通過N-S方程求解獲得;顆粒與連續(xù)相之間通過相互作用力模型耦合,具體步驟如下所述[7]。

    (1)計算連續(xù)相流場(此時計算域中不存在離散相)。

    (2)計算從每個噴射源開始的顆粒軌道,從而在計算域中引入離散相。

    (3)使用已經(jīng)得到的顆粒計算結(jié)果中的相間動量、質(zhì)量交換項重新計算連續(xù)相流場。

    (4)計算修正后的連續(xù)相流場中的顆粒軌跡。

    (5)重復(fù)上述兩個步驟,直到獲得收斂解。

    由此可知,DPM模型的計算過程是一個從單相的流場模擬(即第一步計算連續(xù)相流場)到兩相耦合計算的過程,而且單相模擬可得到旋風(fēng)分離器的內(nèi)部氣相的壓力、切向和軸向速度分布。將之與前人理論與實驗研究結(jié)果比較,可間接驗證模擬的可靠性,為接下來的兩相DPM模擬研究提供保障。

    旋風(fēng)分離器的分離效率計算方法見式(1)[8]。

    式中,Mo為該直徑顆粒在溢流口的質(zhì)量流率,kg/s;Mi為該直徑顆粒在進(jìn)料口的質(zhì)量流率,kg/s。通過Fluent模擬可監(jiān)測旋風(fēng)分離器入口和溢流口的液滴質(zhì)量流率,從而得出旋風(fēng)分離器的分離效率。

    1.1 幾何模型的建立及網(wǎng)格劃分

    圖1 模擬模型

    圖1(a)為天然氣脫蠟旋風(fēng)分離器的結(jié)構(gòu)簡圖,表1為旋風(fēng)分離器的主要結(jié)構(gòu)參數(shù),旋風(fēng)分離器直徑為300 mm,其入口形式為斜切式,提高了旋流器的流場穩(wěn)定性,減小了壓力損失。

    圖1(b)為gambit建立的模擬模型,整個三維模型的網(wǎng)格劃分采用六面體結(jié)構(gòu)化網(wǎng)格如圖1 (c),劃分了32698個網(wǎng)格。

    表1 旋風(fēng)分離器的主要結(jié)構(gòu)參數(shù)

    1.2 計算方法

    在流場的數(shù)值模擬中對于各個通用方程的離散,目前比較常用的差分格式主要有一階迎風(fēng)格式、二階迎風(fēng)格式、QUICK格式和中心差分格式等。當(dāng)流動具有強(qiáng)旋轉(zhuǎn)特性時,QUICK差分格式明顯具有較高的精度,能最準(zhǔn)確地預(yù)報旋流器內(nèi)的流場[9]。采用Fluent Pressure Based隱式求解器;計算模型選用RSM模型,控制方程中的壓力-速度耦合項選用協(xié)調(diào)一致的SIMPLEC算法。

    1.3 邊界條件及部分物性參數(shù)

    1.3.1 流體相的邊界條件

    (1)進(jìn)口邊界 旋風(fēng)分離器入口為流場的速度入口,速度方向垂直于入口界面指向入口內(nèi)側(cè),壓力為該處靜壓。

    (2)出口邊界 旋風(fēng)分離器的溢流口,為自由出口。底流口設(shè)為壁面,為液滴捕集面之一。

    (3)壁面 在旋風(fēng)分離器的近壁區(qū),湍流輸運(yùn)性質(zhì)會發(fā)生急劇的變化,應(yīng)用壁面函數(shù)法進(jìn)行近壁處理,解決計算精確度的問題。

    (4)基礎(chǔ)物性 旋風(fēng)分離器內(nèi)的流體為天然氣,分散相為石蠟液滴,其組成及相關(guān)的物性參數(shù)如表2所示。

    1.3.2 分散相邊界條件

    液滴粒徑分布通過對液滴尺寸進(jìn)行Rosin-Rammler分布擬合[10],所有的顆粒尺寸范圍被劃分為離散的尺寸組,每個尺寸組由該組射流源中的單個顆粒所代表,在計算運(yùn)動軌跡的時候所采用的就是這個平均直徑。圖2 為模擬采用的液滴的粒徑分布。

    溢流口設(shè)置為逃逸出口,當(dāng)液滴運(yùn)動到溢流口的時候會逃逸,此時停止對顆粒的追蹤,終止對該顆粒軌道的計算。

    壁面和底流口均設(shè)置為捕集面,液滴接觸壁面或到達(dá)底流口就認(rèn)為液滴得到分離,終止對該顆粒軌道的計算。

    表2 物料及相關(guān)物性參數(shù)

    圖2 Rosin-Rammler分布

    2 結(jié)果與分析

    2.1 壓力分布

    圖 3(a)為旋風(fēng)分離器在進(jìn)口流量Qi為500 m3/h的縱截面壓力分布云圖;圖3(b)為距離底流口900 mm處橫截面上的壓力分布。壓力呈對稱分布,由外至內(nèi)逐漸減小,壁面處壓力最高,軸心處壓力最低,有負(fù)壓產(chǎn)生。壓力分布比較穩(wěn)定,可以看出流場比較穩(wěn)定,有利于兩相分離。

    2.2 速度分布

    2.2.1 切向速度

    切向速度反映了流體在旋風(fēng)分離器內(nèi)旋轉(zhuǎn)的快慢,決定了離心力的大小,是使液滴沉降的驅(qū)動力。圖4(a)為旋風(fēng)分離器在進(jìn)口流量Qi為500 m3/h時的縱截面切向速度分布云圖;圖4(b)為距離底流口900 mm處橫截面上的切向速度分布。

    模擬發(fā)現(xiàn)切向速度滿足一般旋風(fēng)分離器的準(zhǔn)自由渦規(guī)律:從邊壁開始隨著半徑的減小,切向速度緩慢增大,在某一徑向位置達(dá)到最大值,然后急劇降低,在整個直徑范圍內(nèi),呈現(xiàn)“M”形駝峰分布。

    2.2.2 軸向速度

    軸向速度反映了流體向兩個軸向出口的流動情況。圖5(a)為旋風(fēng)分離器在進(jìn)口流量Qi為500 m3/h時的縱截面上的軸向速度分布云圖和零軸速包絡(luò)面;圖5(b)為距離底流口900 mm處橫截面上的軸向向速度分布??拷獗诘臑橄蛳碌耐庑?,速度逐漸減??;零軸速包絡(luò)面內(nèi)側(cè)則是向上的外旋流,軸向速度逐漸增大。在溢流管中心處軸向速度為負(fù),出現(xiàn)回流。

    圖3 壓力分布(Qi=500m3·h-1)

    圖4 切向速度分布(Qi=500 m3·h-1)

    圖5 軸向速度分布(Qi=500 m3·h-1)

    模擬所得的流場分布規(guī)律——壓力、切向和軸向速度分布規(guī)律,符合一般旋風(fēng)分離器的流場分布規(guī)律[11-12],說明了模擬的可靠性,為接下來的兩相模擬和分離效率的研究提供了保障。

    2.3 兩相分布

    當(dāng)石蠟液滴的體積濃度低于5%時,可采用離散相模型(DPM)模擬旋風(fēng)分離器內(nèi)石蠟液滴的分離過程。天然氣和石蠟液滴的混合流體由進(jìn)料口進(jìn)入旋風(fēng)分離器,在腔內(nèi)高速旋轉(zhuǎn)產(chǎn)生離心力場。在離心力的作用下,石蠟液滴發(fā)生離心沉降,遷移到四周,沿著壁面向下流動,最后經(jīng)底流口被分離出去。一些小粒徑的液滴這隨著氣相由溢流口排出。圖6為進(jìn)口流量Qi為500 m3/h時旋風(fēng)分離器內(nèi)部的石蠟液滴的質(zhì)量濃度分布,發(fā)現(xiàn)在其錐段底部靠近壁面處的液滴濃度較高,符合實際規(guī)律。

    2.4 分離效率

    分離效率是旋風(fēng)分離器最主要的性能參數(shù)。對于一個分離器,若分離某種物料的效率為95%,也許效率很高,但如果換一個細(xì)一點(diǎn)的進(jìn)料,效率就達(dá)不到95%;反之,若分級粗一些的物料效率就會高于95%。因此,一個分離器的分離性能應(yīng)當(dāng)用分離效率和粒度一起考慮,即得出粒級效率曲線。

    在粒級效率曲線上對應(yīng)于50%分離效率的粒度被定義為分割粒度,用x50表示。標(biāo)準(zhǔn)型旋風(fēng)分離器的分割粒度x50可用式(2)計算[13]。

    圖6 石蠟液滴質(zhì)量濃度分布(Qi=500m3·h-1)

    式中,D為旋風(fēng)分離器直徑,m;Vi為入口速度,m/s;μ為氣體黏度,Pa·s;ρp、ρg分別為液滴和氣體的密度。

    圖7是模擬的入口速度不同時的旋風(fēng)分離器的分割粒度和根據(jù)式(2)理論計算所得的分割粒度,可以發(fā)現(xiàn)兩者之間的誤差較小,且分布規(guī)律一致:隨著入口速度的增大,旋風(fēng)分離器的分割粒徑在減小。故兩相模擬具有一定的可靠性。

    2.4.1 進(jìn)口流量對分離效率的影響

    進(jìn)口流量Qi是旋風(fēng)分離器重要操作參數(shù),直接影響切向速度和壓力降。圖8是距離底流口900 mm截面上不同進(jìn)口流量Qi下的切向速度;圖9為不同進(jìn)口流量Qi下旋風(fēng)分離器的級效率曲線。

    圖7 不同入口速度下的理論和模擬的分割粒度

    圖8 z=900 mm截面處不同進(jìn)口流量下的切向速度分布

    圖9 不同進(jìn)口流量下的級效率曲線

    隨著進(jìn)口流量Qi的增加,級效率曲線向左偏移,分割粒徑x50減小,分級效率提高。這是因為隨著進(jìn)口流量Qi的增加,旋風(fēng)分離器內(nèi)流體的切向速度增大,產(chǎn)生更大的離心分離作用,使得液滴更快靠近邊壁而被分離出。

    在所研究的進(jìn)口流量范圍內(nèi),對于小于5 μm的液滴,流量的增加沒有顯著增加其分離效率。這是由于粒度小的液滴沉降速度小,來不及完成沉降,被直接夾帶出溢流口,未實現(xiàn)分離。

    當(dāng)進(jìn)口流量Qi大于800 m3/h后,級效率曲線變化不是很顯著。這是由于在此時旋流器內(nèi)的湍流強(qiáng)度增大,不利于液滴的分離,同時,湍流加速了液滴的破碎,使液滴粒徑變小,增加了分離難度。

    2.4.2 段長度L1和錐段長度L2對分離效率的影響

    柱段長度和錐段長度增加,增加了旋風(fēng)分離器的長度,延長了液滴在旋風(fēng)分離器內(nèi)的停留時間。流體在旋風(fēng)分離器內(nèi)的停留時間決定于流體在旋風(fēng)分離器內(nèi)運(yùn)動的速度和經(jīng)過的路徑長短,或流體的軸向速度W和路徑的軸向距離即旋風(fēng)分離器的長度L。旋風(fēng)分離器內(nèi)的物料的平均停留時間可以表達(dá)為式(3)。

    由式(3)可知,停留時間與進(jìn)口流量Qi、旋風(fēng)分離器直徑D和柱段長徑比L1/D以及錐段長徑比L2/D等有關(guān)。加大D、L1/D或L2/D,或者減小Qi都能延長停留時間,使更小的液滴完成沉降,從而提高分離效率。但是,加大D或減小Qi都會顯著地降低離心加速度,使得液滴沉降的驅(qū)動力不夠。因此,只有增加旋風(fēng)分離器的長度才是增加停留時間的有效手段。

    通過對不同柱段和錐段長度的旋風(fēng)分離器的DPM兩相模擬,增加旋風(fēng)分離器的長度可有效提高旋風(fēng)分離器的分離效率。

    保持錐段長度L2不變,圖10為不同柱段長徑比L1/D的旋風(fēng)分離器的級效率曲線。旋風(fēng)分離器的級效率曲線隨著柱段長徑比L1/D的增大向左偏移,x50減小,分離效率提高。

    保持柱段長度L1不變,隨著錐段長度L2的增加,如圖11,級效率曲線向左偏移,x50減小,分離效率增大。

    在實際工程中,更長的旋風(fēng)分離器的設(shè)計和制造成本更高,需要綜合考慮合適的長徑比。

    圖10 不同柱段長度L1的級效率曲線(Qi=500 m3·h-1)

    圖11 不同錐段長度L2的級效率曲線(Qi=500 m3·h-1)

    3 結(jié) 論

    (1)使用RSM湍流模型模擬天然氣脫蠟旋風(fēng)分離器,所得的流場分布規(guī)律符合一般旋風(fēng)分離器的流場分布規(guī)律:旋風(fēng)分離器內(nèi)部壓力呈對稱分布,由壁面向中心逐漸減小,在中心軸處壓力最低并有負(fù)壓產(chǎn)生;切向速度符合準(zhǔn)自由渦分布規(guī)律,呈現(xiàn)“M”形駝峰分布;軸向速度分布反映了向下的外旋流和向上的內(nèi)旋流的流動規(guī)律。

    (2)使用DPM模型模擬旋風(fēng)分離器內(nèi)的兩相流動,隨著進(jìn)口流量的增加,旋風(fēng)分離器的級效率曲線向左移動,x50減小,分離效率提高,與理論計算相符。當(dāng)入口流量為1000 m3/h時,x50可以達(dá)到5.3 μm。大粒徑液滴的分離效果明顯,但在所研究的處理量范圍內(nèi),處理量的變化不能明顯地影響小于5 μm液滴的分離效率。

    (3)柱段和錐段長度的增加增長了旋風(fēng)分離器的的總長度,延長了液滴在旋風(fēng)分離器內(nèi)的停留時間,旋風(fēng)分離器的級效率曲線向左偏移,x50減小,分離效率增大。

    [1] 馮異勇. 石蠟沉積及清除數(shù)值模擬研究[D]. 四川:西南石油學(xué)院,2002.

    [2] 袁惠新. 天然氣凈化裝置:中國,102786995 A[P]. 2012-08-31.

    [3] Hsieh K T,Rajamani R K. Mathematiacal model of the hydrocyclone based on physics of fluid flow[J]. AIChE Journal,1991,37(5):735-746.

    [4] 劉子紅,肖波,楊家寬. 旋風(fēng)分離器兩相流研究綜述[J]. 過濾與分離,2003,13(1):42-45.

    [5] 王海剛,劉石. 不同湍流模型在旋風(fēng)分離器三維數(shù)值模擬中的應(yīng)用和比較[J]. 熱能動力工程,2003,18(4):337-342.

    [6] 葛坡,袁惠新,付雙成. 對稱多入口型旋風(fēng)分離器的數(shù)值模擬[J].化工進(jìn)展,2012,31(2):296-299.

    [7] 牟春宇. 基于CFD的旋風(fēng)除塵器分離效率的數(shù)值模型研究[D].無錫:江南大學(xué),2008.

    [8] 許姸霞. 水力旋流分離過程數(shù)值模擬與分析[D]. 上海:華東理工大學(xué),2012.

    [9] 劉淑艷,張雅,王保國. 用RSM模擬旋風(fēng)分離器內(nèi)的三維湍流流場[J]. 北京理工大學(xué)學(xué)報,2005,25(5):377-379.

    [10] 張憲瑞,閆鵬,畢榮山. 氣液噴射器的不同液滴粒徑分布特性研究[J]. 科學(xué)技術(shù)與工程,2012,12(6):3933-3936.

    [11] 魏耀東,宋健斐,陳建義,等. 旋風(fēng)分離器內(nèi)氣相流場的相似?;治?Ⅰ) 流動參數(shù)[J]. 化工學(xué)報,2010,61(9):2265-2273.

    [12] 宋健斐,楊光福,陳建義,等. 旋風(fēng)分離器內(nèi)氣相流場的相似?;治?Ⅱ) 尺寸參數(shù)[J]. 化工學(xué)報,2010,61(9):2274-2279.

    [13] 袁惠新,馮骉. 分離工程[M]. 北京:中國石化出版社,2001.

    Numerical simulation of the separation efficiency of natural gas dewaxing cyclone separator

    YUAN Huixin,F(xiàn)ANG Yi,F(xiàn)U Shuangcheng
    (School of Mechanical Engineering,Changzhou University,Changzhou 213016,Jiangsu,China)

    The cyclone separator has high separation efficiency and is not easy to be blocked,it has good effect on the natural gas dewaxing. Numerical simulation of the two-phase flow of the CYG-S natural gas dewaxing cyclone separator was studied by using CFD software Fluent,and get the internal pressure,tangent and axial velocity distributions. The simulation and calculated cut size x50in different inlet velocity are contrasted,the results has a good degree of agreement,so the two-phase simulation are reliable. The results show that the mass concentration of paraffin droplet is higher at the bottom of the cone section. With the increasing of the inlet flow,the separation efficiency of cyclone separator was improved,under the 1000 m3/h inlet flow,the x50can reach 5.3 μm. However,in the studied range of inlet flow,the increase rate of the drop less than 5μm is not obvious. The lengthened of column and taper section increases the length of cyclone separator,extents the residence time of the droplets in the cyclone separator,and improves the separation efficiency of the cyclone separator.

    cyclone separator;numerical simulation;separation efficiency;inlet flow;column and taper section length

    TQ 051.8

    A

    1000-6613(2014)01-0043-07

    10.3969/j.issn.1000-6613.2014.01.008

    2013-07-16;修改稿日期:2013-08-20。

    及聯(lián)系人:袁惠新(1957—),男,教授,博士生導(dǎo)師,從事分離技術(shù)與設(shè)備研究。E-mail yuanhuixin2000@163.com。

    猜你喜歡
    切向速度段長度旋風(fēng)
    旋風(fēng)分離器內(nèi)氣相旋轉(zhuǎn)流不穩(wěn)定性的實驗研究
    雙旋流耦合式旋流反應(yīng)器內(nèi)切向速度分布研究
    過渡段長度對混合梁橋的受力影響
    腦力急旋風(fēng)
    腦力急旋風(fēng)
    腦力急旋風(fēng)
    腦力急旋風(fēng)
    高強(qiáng)鋼組合K型偏心支撐框架耗能梁段長度研究
    不同自由段長度預(yù)應(yīng)力錨索力學(xué)特性分析
    預(yù)應(yīng)力錨索不同錨固長度對錨索及巖體的影響
    河南科技(2013年7期)2013-08-14 05:28:02
    熟女电影av网| 听说在线观看完整版免费高清| 男女国产视频网站| 成人性生交大片免费视频hd| 精品国内亚洲2022精品成人| ponron亚洲| 国产私拍福利视频在线观看| 日本av手机在线免费观看| 麻豆乱淫一区二区| 高清视频免费观看一区二区 | 国产一级毛片七仙女欲春2| 久久久久精品久久久久真实原创| 最近中文字幕高清免费大全6| 秋霞在线观看毛片| 麻豆成人av视频| 亚洲精品自拍成人| 乱人视频在线观看| 女人被狂操c到高潮| 国产伦理片在线播放av一区| 亚洲欧美日韩无卡精品| 青春草亚洲视频在线观看| 人妻少妇偷人精品九色| 国产午夜精品一二区理论片| 日本五十路高清| 色综合站精品国产| 麻豆av噜噜一区二区三区| 赤兔流量卡办理| 成人毛片a级毛片在线播放| 97超视频在线观看视频| 国产老妇伦熟女老妇高清| 欧美高清性xxxxhd video| 国产精品99久久久久久久久| 汤姆久久久久久久影院中文字幕 | 久久综合国产亚洲精品| 18禁动态无遮挡网站| 日韩国内少妇激情av| av播播在线观看一区| or卡值多少钱| 免费看日本二区| 国产成人午夜福利电影在线观看| av播播在线观看一区| 一级二级三级毛片免费看| 国产中年淑女户外野战色| 看片在线看免费视频| 午夜福利高清视频| 18禁在线播放成人免费| 成人亚洲精品av一区二区| 久久鲁丝午夜福利片| 美女大奶头视频| 亚洲av一区综合| 建设人人有责人人尽责人人享有的 | 熟妇人妻久久中文字幕3abv| 国产亚洲5aaaaa淫片| 国产视频首页在线观看| 久久草成人影院| 伦精品一区二区三区| 国产三级在线视频| 国产免费视频播放在线视频 | 老司机影院毛片| 草草在线视频免费看| 国产成人a∨麻豆精品| 99久久人妻综合| 老司机影院毛片| 我的女老师完整版在线观看| 男插女下体视频免费在线播放| 国产91av在线免费观看| 看黄色毛片网站| 国产一级毛片七仙女欲春2| 久久久久精品久久久久真实原创| 久久久久久大精品| 久久久精品94久久精品| 国产免费视频播放在线视频 | 国产淫片久久久久久久久| 成年免费大片在线观看| 色视频www国产| 中文字幕人妻熟人妻熟丝袜美| 午夜亚洲福利在线播放| 一二三四中文在线观看免费高清| 亚洲国产最新在线播放| 成人性生交大片免费视频hd| 又爽又黄无遮挡网站| 啦啦啦韩国在线观看视频| 亚洲av中文字字幕乱码综合| 久久久久久久久久久免费av| 久久婷婷人人爽人人干人人爱| 午夜福利在线观看吧| 欧美极品一区二区三区四区| 成人亚洲欧美一区二区av| 七月丁香在线播放| 成人性生交大片免费视频hd| 亚洲国产日韩欧美精品在线观看| 久久久久精品久久久久真实原创| 久久午夜福利片| 精品人妻一区二区三区麻豆| 国产精品综合久久久久久久免费| 超碰av人人做人人爽久久| 欧美丝袜亚洲另类| 99久久九九国产精品国产免费| 亚洲欧美成人综合另类久久久 | 成人毛片a级毛片在线播放| 一级毛片电影观看 | 我的老师免费观看完整版| 高清av免费在线| 久久精品夜色国产| 看十八女毛片水多多多| 免费电影在线观看免费观看| 三级毛片av免费| 亚洲,欧美,日韩| 欧美日韩精品成人综合77777| 日韩av在线大香蕉| 日韩视频在线欧美| 91aial.com中文字幕在线观看| 两个人视频免费观看高清| 久久精品国产亚洲av天美| 99视频精品全部免费 在线| 黑人高潮一二区| 最近2019中文字幕mv第一页| 日韩,欧美,国产一区二区三区 | 久久这里有精品视频免费| 日产精品乱码卡一卡2卡三| 一夜夜www| 国产久久久一区二区三区| 亚洲中文字幕日韩| 最近视频中文字幕2019在线8| 久久精品国产99精品国产亚洲性色| 伦理电影大哥的女人| 免费电影在线观看免费观看| 在线免费观看不下载黄p国产| 内射极品少妇av片p| 午夜a级毛片| 午夜福利视频1000在线观看| 久久久久免费精品人妻一区二区| 国产一区二区在线观看日韩| 日韩高清综合在线| 国产亚洲精品av在线| 色噜噜av男人的天堂激情| 国模一区二区三区四区视频| 国产午夜精品一二区理论片| 久久久久久久久大av| 男人和女人高潮做爰伦理| 午夜免费男女啪啪视频观看| 日本wwww免费看| 日本欧美国产在线视频| 亚洲在久久综合| 国产伦精品一区二区三区视频9| 精品欧美国产一区二区三| 日韩欧美在线乱码| 亚洲伊人久久精品综合 | 一级黄色大片毛片| АⅤ资源中文在线天堂| 亚洲精品456在线播放app| 一级黄色大片毛片| 一级二级三级毛片免费看| 欧美成人午夜免费资源| 久久欧美精品欧美久久欧美| 丰满少妇做爰视频| 亚洲性久久影院| 热99在线观看视频| 亚洲成人久久爱视频| 麻豆成人av视频| 久久久成人免费电影| av女优亚洲男人天堂| 色5月婷婷丁香| 亚洲欧洲国产日韩| 99热6这里只有精品| 九草在线视频观看| 久久久久久大精品| 久久人妻av系列| 国产免费一级a男人的天堂| 国产视频首页在线观看| 丝袜喷水一区| 婷婷色综合大香蕉| 国产精品美女特级片免费视频播放器| .国产精品久久| 亚洲av成人av| 亚洲国产成人一精品久久久| 日韩人妻高清精品专区| 午夜福利在线观看吧| 看片在线看免费视频| 看非洲黑人一级黄片| 久久久色成人| 一个人看的www免费观看视频| 高清av免费在线| 国产v大片淫在线免费观看| 九九在线视频观看精品| 狂野欧美激情性xxxx在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产成人精品久久久久久| 最近视频中文字幕2019在线8| 亚洲不卡免费看| 国产免费一级a男人的天堂| 综合色丁香网| 国产精品电影一区二区三区| 精品国产三级普通话版| 在线免费十八禁| 亚洲最大成人中文| 日韩强制内射视频| 国产色婷婷99| 晚上一个人看的免费电影| 全区人妻精品视频| 国产黄片视频在线免费观看| 精品久久久噜噜| 日本色播在线视频| 国产成人aa在线观看| 大香蕉久久网| 久久综合国产亚洲精品| 国产又黄又爽又无遮挡在线| 在线观看一区二区三区| 国产精品一二三区在线看| 亚洲成人av在线免费| 91久久精品国产一区二区三区| 亚洲精品乱码久久久v下载方式| 日本黄大片高清| 免费观看人在逋| 国产色爽女视频免费观看| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 亚洲国产高清在线一区二区三| 超碰av人人做人人爽久久| 亚洲在线自拍视频| 搡女人真爽免费视频火全软件| 人妻少妇偷人精品九色| 精品熟女少妇av免费看| 亚洲av电影在线观看一区二区三区 | 观看美女的网站| 丝袜喷水一区| 3wmmmm亚洲av在线观看| 建设人人有责人人尽责人人享有的 | 美女黄网站色视频| 国产91av在线免费观看| 最近中文字幕2019免费版| 亚洲欧美日韩无卡精品| 久久精品熟女亚洲av麻豆精品 | 美女xxoo啪啪120秒动态图| 一级毛片我不卡| 欧美成人午夜免费资源| 欧美不卡视频在线免费观看| 免费看a级黄色片| 国产欧美另类精品又又久久亚洲欧美| 97热精品久久久久久| 亚洲图色成人| h日本视频在线播放| 日本色播在线视频| 中国国产av一级| 成人av在线播放网站| 国产精品麻豆人妻色哟哟久久 | 波多野结衣高清无吗| 中国国产av一级| 精品欧美国产一区二区三| 岛国毛片在线播放| 我要看日韩黄色一级片| 夜夜爽夜夜爽视频| 国产成人免费观看mmmm| 特大巨黑吊av在线直播| 人人妻人人澡欧美一区二区| 人妻少妇偷人精品九色| 麻豆久久精品国产亚洲av| 一区二区三区四区激情视频| 国产午夜精品论理片| 久久人人爽人人爽人人片va| 观看免费一级毛片| 久久久成人免费电影| 欧美丝袜亚洲另类| 身体一侧抽搐| 免费观看的影片在线观看| 亚洲最大成人中文| 观看美女的网站| 看片在线看免费视频| 国产精品国产三级国产av玫瑰| 亚洲精品456在线播放app| 国产乱人视频| 好男人在线观看高清免费视频| 色综合色国产| av在线天堂中文字幕| 欧美97在线视频| 国产精品av视频在线免费观看| 三级国产精品片| 亚洲av日韩在线播放| 性插视频无遮挡在线免费观看| 国产伦理片在线播放av一区| 我要看日韩黄色一级片| 国产精品人妻久久久久久| 国产成人精品婷婷| 日本一本二区三区精品| www.色视频.com| 蜜桃亚洲精品一区二区三区| 深夜a级毛片| 最近的中文字幕免费完整| 在线天堂最新版资源| 91久久精品国产一区二区成人| 久久久久免费精品人妻一区二区| 激情 狠狠 欧美| 啦啦啦韩国在线观看视频| 国产精品一及| 中文字幕人妻熟人妻熟丝袜美| 色噜噜av男人的天堂激情| 久久精品夜色国产| 男女啪啪激烈高潮av片| 亚洲欧美日韩高清专用| 日本-黄色视频高清免费观看| 亚洲国产精品成人综合色| av在线亚洲专区| 九色成人免费人妻av| 精品久久久久久久久久久久久| eeuss影院久久| 久久99热6这里只有精品| 夜夜看夜夜爽夜夜摸| 一二三四中文在线观看免费高清| 日韩国内少妇激情av| 乱人视频在线观看| 日韩欧美在线乱码| 激情 狠狠 欧美| a级一级毛片免费在线观看| 国产高清视频在线观看网站| av在线老鸭窝| 爱豆传媒免费全集在线观看| 又黄又爽又刺激的免费视频.| 久久这里只有精品中国| 高清毛片免费看| 婷婷色综合大香蕉| 欧美另类亚洲清纯唯美| 婷婷色综合大香蕉| 亚洲av电影不卡..在线观看| 欧美bdsm另类| 午夜福利成人在线免费观看| 国内揄拍国产精品人妻在线| 日本三级黄在线观看| 久久精品国产自在天天线| 韩国高清视频一区二区三区| 只有这里有精品99| 免费看av在线观看网站| 黄片无遮挡物在线观看| 女人被狂操c到高潮| 久久99精品国语久久久| 青春草国产在线视频| 天天躁夜夜躁狠狠久久av| 99久国产av精品国产电影| 中文亚洲av片在线观看爽| 卡戴珊不雅视频在线播放| 久久鲁丝午夜福利片| 午夜亚洲福利在线播放| 亚洲国产精品合色在线| 丝袜美腿在线中文| 不卡视频在线观看欧美| 一区二区三区四区激情视频| 三级经典国产精品| 久久韩国三级中文字幕| 亚洲成人av在线免费| 国产精品久久视频播放| 国产高清不卡午夜福利| 一区二区三区免费毛片| 在线a可以看的网站| 久久午夜福利片| 波野结衣二区三区在线| 亚洲精品国产av成人精品| 国产熟女欧美一区二区| 在现免费观看毛片| 在线观看av片永久免费下载| 久久精品国产亚洲网站| 亚洲自拍偷在线| 又爽又黄无遮挡网站| 亚洲自拍偷在线| 欧美3d第一页| 五月伊人婷婷丁香| 直男gayav资源| 草草在线视频免费看| 亚洲国产精品久久男人天堂| 搡女人真爽免费视频火全软件| 内地一区二区视频在线| 菩萨蛮人人尽说江南好唐韦庄 | 久久国产乱子免费精品| 精品99又大又爽又粗少妇毛片| 少妇丰满av| 精品久久久久久成人av| 国产成人免费观看mmmm| 成人综合一区亚洲| a级毛片免费高清观看在线播放| 久久这里有精品视频免费| 国产精品日韩av在线免费观看| 天堂av国产一区二区熟女人妻| 成年女人永久免费观看视频| 久久亚洲国产成人精品v| av在线天堂中文字幕| 亚洲,欧美,日韩| 久久99热这里只频精品6学生 | 亚洲精品国产成人久久av| 日韩欧美在线乱码| 日本爱情动作片www.在线观看| 国产成人freesex在线| 久久精品久久精品一区二区三区| 亚洲成色77777| 能在线免费看毛片的网站| 精品国产露脸久久av麻豆 | 三级男女做爰猛烈吃奶摸视频| 男女视频在线观看网站免费| 国产精品久久久久久久电影| 日韩亚洲欧美综合| 22中文网久久字幕| 麻豆av噜噜一区二区三区| 少妇人妻一区二区三区视频| 国产精品久久电影中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 欧美丝袜亚洲另类| 边亲边吃奶的免费视频| 激情 狠狠 欧美| 高清日韩中文字幕在线| 久久久久九九精品影院| 中文资源天堂在线| 黑人高潮一二区| 国产成人aa在线观看| 国产白丝娇喘喷水9色精品| 亚洲精品国产成人久久av| 老师上课跳d突然被开到最大视频| 边亲边吃奶的免费视频| 亚洲av熟女| 少妇的逼水好多| 我的女老师完整版在线观看| 男人狂女人下面高潮的视频| 天堂中文最新版在线下载 | 别揉我奶头 嗯啊视频| 国产av在哪里看| 日本欧美国产在线视频| 寂寞人妻少妇视频99o| av女优亚洲男人天堂| 午夜日本视频在线| 18禁动态无遮挡网站| 菩萨蛮人人尽说江南好唐韦庄 | 午夜a级毛片| 乱码一卡2卡4卡精品| 国产精品人妻久久久影院| 精品国内亚洲2022精品成人| 偷拍熟女少妇极品色| 九草在线视频观看| 最后的刺客免费高清国语| 亚洲电影在线观看av| 国产精品美女特级片免费视频播放器| 国产一区二区三区av在线| 免费黄网站久久成人精品| 久久久久久伊人网av| 啦啦啦观看免费观看视频高清| 淫秽高清视频在线观看| 岛国毛片在线播放| 国内精品宾馆在线| 国产亚洲5aaaaa淫片| 中文字幕亚洲精品专区| 国产成人免费观看mmmm| 精品久久久久久久人妻蜜臀av| 偷拍熟女少妇极品色| 免费无遮挡裸体视频| 最近手机中文字幕大全| 亚洲精品日韩在线中文字幕| 精品99又大又爽又粗少妇毛片| 国产精品久久久久久精品电影| 久久久久久久亚洲中文字幕| 91久久精品国产一区二区三区| 男女国产视频网站| 亚洲国产日韩欧美精品在线观看| 一区二区三区四区激情视频| 欧美三级亚洲精品| 亚洲精品国产av成人精品| 亚洲av二区三区四区| 校园人妻丝袜中文字幕| 三级经典国产精品| 亚洲欧美日韩卡通动漫| 国产成人一区二区在线| 听说在线观看完整版免费高清| 99久久九九国产精品国产免费| 国语对白做爰xxxⅹ性视频网站| 纵有疾风起免费观看全集完整版 | 韩国高清视频一区二区三区| 亚洲精品成人久久久久久| 日韩欧美国产在线观看| 国产亚洲午夜精品一区二区久久 | 国产单亲对白刺激| 亚洲av.av天堂| 国产淫语在线视频| 天美传媒精品一区二区| 中文亚洲av片在线观看爽| 国产黄色视频一区二区在线观看 | videos熟女内射| 亚洲内射少妇av| av在线亚洲专区| 亚洲高清免费不卡视频| 边亲边吃奶的免费视频| 极品教师在线视频| 久久久久久久久久黄片| av在线亚洲专区| videossex国产| 在线观看av片永久免费下载| 亚洲av电影不卡..在线观看| 长腿黑丝高跟| 亚洲最大成人中文| 精品欧美国产一区二区三| 亚洲欧美精品专区久久| 级片在线观看| 噜噜噜噜噜久久久久久91| 亚洲国产精品合色在线| 欧美色视频一区免费| 自拍偷自拍亚洲精品老妇| 高清在线视频一区二区三区 | 久久久a久久爽久久v久久| 国产精品爽爽va在线观看网站| 亚洲色图av天堂| 简卡轻食公司| 最新中文字幕久久久久| 一级毛片aaaaaa免费看小| 校园人妻丝袜中文字幕| 亚洲中文字幕一区二区三区有码在线看| 中文字幕精品亚洲无线码一区| 看非洲黑人一级黄片| 黄色欧美视频在线观看| 我要搜黄色片| 午夜福利在线观看免费完整高清在| 国产亚洲一区二区精品| 最近最新中文字幕免费大全7| 日韩在线高清观看一区二区三区| 成人无遮挡网站| 日韩精品青青久久久久久| 亚洲精品乱码久久久v下载方式| 亚洲欧洲日产国产| 三级经典国产精品| 午夜激情福利司机影院| 国产成人91sexporn| 国产精品人妻久久久久久| 久久午夜福利片| 午夜日本视频在线| 国产乱人偷精品视频| or卡值多少钱| 中文亚洲av片在线观看爽| 桃色一区二区三区在线观看| 在线免费观看的www视频| 亚洲欧美成人精品一区二区| 美女内射精品一级片tv| 久久精品人妻少妇| 亚洲不卡免费看| 淫秽高清视频在线观看| 美女脱内裤让男人舔精品视频| 99热全是精品| 美女黄网站色视频| 内射极品少妇av片p| 在现免费观看毛片| 3wmmmm亚洲av在线观看| 国产精品野战在线观看| 色播亚洲综合网| 日韩一本色道免费dvd| 蜜桃久久精品国产亚洲av| 国国产精品蜜臀av免费| 一级黄色大片毛片| 韩国高清视频一区二区三区| 亚洲自偷自拍三级| 国内精品一区二区在线观看| 熟女人妻精品中文字幕| 久久久成人免费电影| 亚洲国产日韩欧美精品在线观看| 美女cb高潮喷水在线观看| 国产精品福利在线免费观看| 亚洲在线观看片| 最近中文字幕高清免费大全6| 亚洲精品aⅴ在线观看| 成人亚洲精品av一区二区| 国国产精品蜜臀av免费| 成人亚洲欧美一区二区av| 韩国高清视频一区二区三区| 国产一区亚洲一区在线观看| 亚洲国产精品合色在线| 最近的中文字幕免费完整| 日日撸夜夜添| 一个人看视频在线观看www免费| 日韩亚洲欧美综合| 久久精品夜夜夜夜夜久久蜜豆| 亚洲人成网站在线观看播放| 性色avwww在线观看| 精品人妻视频免费看| 免费人成在线观看视频色| av卡一久久| 蜜臀久久99精品久久宅男| 亚洲av成人精品一区久久| 亚洲av熟女| 国产三级中文精品| 久久人妻av系列| av视频在线观看入口| 午夜福利高清视频| 日韩欧美国产在线观看| 搡老妇女老女人老熟妇| 中文字幕av成人在线电影| 久久久久久久久久久丰满| 免费大片18禁| 久久久久久九九精品二区国产| 尤物成人国产欧美一区二区三区| 欧美一区二区亚洲| 国产免费一级a男人的天堂| 亚洲图色成人| 综合色丁香网| 在线观看av片永久免费下载| 国产大屁股一区二区在线视频| 最近中文字幕高清免费大全6| 精品人妻偷拍中文字幕| 老女人水多毛片| 欧美一区二区亚洲| 成人毛片a级毛片在线播放| 精品人妻一区二区三区麻豆| 一夜夜www| 国产精品99久久久久久久久| 久久久国产成人免费| www.色视频.com| 国产成人freesex在线| 三级毛片av免费| 国产成年人精品一区二区| 亚洲色图av天堂| 免费看av在线观看网站| 亚洲欧洲日产国产| 日本黄色视频三级网站网址| 欧美丝袜亚洲另类| 天堂√8在线中文| 最近中文字幕2019免费版|