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

    汽油裝罐油氣擴散排放的實驗測定及數(shù)值模擬

    2016-10-13 18:52:02黃維秋王兆利紀虹趙晨露呂愛華徐先陽王翊紅
    化工學報 2016年12期
    關鍵詞:油口損耗率油罐

    黃維秋,王兆利,紀虹,趙晨露,呂愛華,徐先陽,王翊紅

    ?

    汽油裝罐油氣擴散排放的實驗測定及數(shù)值模擬

    黃維秋1,王兆利1,紀虹1,趙晨露2,呂愛華1,徐先陽1,王翊紅1

    (1江蘇省油氣儲運技術重點實驗室(常州大學),江蘇常州 213016;2上海河圖工程股份有限公司,上海 200000)

    正確預測儲油罐收油作業(yè)時罐內(nèi)油氣擴散排放規(guī)律對研究油品蒸發(fā)損耗及污染控制具有重要意義,因而兩個關鍵參數(shù)(油罐排放氣液比及損耗率)被重點考慮?;赩OF模型、擴散傳質(zhì)模型和RNG湍流模型,對汽油噴濺式裝油損耗進行數(shù)值模擬研究。分析比較了不同裝油口高度、不同裝油速度、不同油罐初始油氣濃度條件下罐內(nèi)的油氣擴散規(guī)律,并建立汽油裝罐蒸發(fā)損耗實驗平臺驗證以上數(shù)值模擬,模擬值與實驗值吻合良好。研究結果表明:裝油口位置越高,汽油損耗率越大,隨裝油速度增大,高裝油口損耗率最大約為0.34%,低裝油口約為0.025%。增大油罐初始油氣濃度,高裝油口損耗率最大約為0.44%,中裝油口約為0.21%,低裝油口約為0.043%。最后建議固定頂油罐大呼吸API損耗評估計算公式考慮裝油速度及油罐初始油氣濃度對蒸發(fā)損耗的影響,并盡量采用低裝油口及清洗油罐裝油。

    汽油;蒸發(fā);擴散;裝油速度;初始油氣濃度;數(shù)值模擬

    引 言

    石油開采、煉制、儲運及銷售過程中,都會造成一定的蒸發(fā)損耗,不僅會造成資源浪費,對環(huán)境也會產(chǎn)生極大的污染。其中儲運環(huán)節(jié)中裝油過程是主要的損耗階段,裝油口位置對油品損耗及油氣排放濃度(油氣排放因子)的影響甚大。盡管目前各種油罐、罐車、油船基本采用頂部浸沒式或底部裝油,裝油過程中油品的湍流運動能夠得到有效控制,從而產(chǎn)生較少的損耗量[1]。然而油品通過油罐底部裝油口裝油時,在油品浸沒裝油口之前仍然會產(chǎn)生強烈的對流擴散傳質(zhì)。目前,汽油等輕質(zhì)油品以內(nèi)浮頂罐儲存為主。帶有油氣回收處理系統(tǒng)的常壓固定頂罐已經(jīng)成為一種新的儲油方式[2-4]。為此,拱頂罐收油作業(yè)時罐內(nèi)油氣擴散排放的機理分析并正確預測其排放規(guī)律,對研究油品蒸發(fā)損耗及污染控制具有重要意義,也可為目前國內(nèi)正在實施的VOCs排污監(jiān)管及收費提供理論依據(jù)。同時,該研究成果還可為罐車、油船、加油站地下儲罐的油氣排放評估做參考。

    研究油罐噴濺式裝油需要模擬空氣和汽油的兩相流動。主要使用VOF模型捕捉氣液界面從而計算傳質(zhì),數(shù)值模擬在該現(xiàn)象的研究上有廣泛應用[5-7]。Waheed等[8]利用有限元法對處在自然對流或強迫對流中的液滴的傳質(zhì)規(guī)律進行了解算。Haelssig等[9]采用VOF模型追蹤氣液相流的動態(tài)界面,研究了動態(tài)界面的傳質(zhì)規(guī)律。朱玲等[10]針對93#汽油的蒸發(fā)損耗過程進行了實驗研究,分別模擬了上裝浸沒式進油、噴濺式進油、油罐靜儲存和地面敞口靜置4種不同的油品蒸發(fā)狀態(tài)。Hassanvand等[11]利用VOF模型研究了臥式油罐噴射裝油,考察了溫度、裝油速度等操作條件對汽油蒸發(fā)質(zhì)量和蒸發(fā)速率的影響。

    上述研究大都針對自由液面的動態(tài)傳質(zhì)規(guī)律及收油過程中油品的蒸發(fā)率,并未對收油作業(yè)時油罐排放氣液比、損耗率等變化規(guī)律進行研究,以及未對汽油等油品蒸發(fā)損耗進一步核實對比。為此本文通過對油罐噴濺式收油作業(yè)現(xiàn)象進行分析,依據(jù)計算流體力學理論基礎,采用有限體積法求解流場的偏微分方程,選取VOF模型、RNG湍流模型及傳質(zhì)模型,并添加氣液相界面質(zhì)量源項,建立油罐噴濺式裝油模型并采用PISO算法求解。通過比較裝油口高度、裝油速度、初始油氣濃度條件下罐內(nèi)油品蒸發(fā)擴散規(guī)律,為油罐收油損耗評估提供重要參考。

    1 模型建立

    1.1 數(shù)學模型

    油罐噴濺式裝油的流動過程可以用VOF多相流模型計算,在VOF模型中不同的流體組分共用一套動量方程,計算時在整個計算區(qū)域的每個計算單元內(nèi),都記錄下各相組分所占有的體積率或體積分數(shù)。所以,VOF模型通過計算體積分數(shù)的連續(xù)性方程捕捉氣液相界面。該方程為

    式中,v分別對應(,,)3個方向的運動距離,m;v為(,,)3個方向的速度分量,m·s-1;為時間,s;代表體積分數(shù)。

    當計算網(wǎng)格內(nèi)全部充滿氣相時=1,相反則=0。氣液相交界面可以在0<<1時被捕捉到。

    連續(xù)性方程為

    式中,v是所求相的密度,kg·m-3;為所求相的速度矢量,m·s-1。

    動量守恒方程為

    式中,x、x分別表示在軸和軸方向的運動距離,m;,分別表示在軸和軸方向的速度矢量,m·s-1;是動力黏度,Pa·s;是重力加速度,m·s-2;是絕對壓力,Pa;為混合密度,kg·m-3;σ是體積力,N·(m3·s)-1。

    湍流方程:油罐噴濺式裝油屬于湍流流動,其中RNG-模型特別適用于VOF多相流模型[12],因此,本模擬過程選擇RNG-模型。該模型對流動分離和二次流有很好的表現(xiàn),與標準-模型相比,有兩個主要不同點:RNG-模型為湍流黏性增加了一個公式;為耗散率增加了新的傳輸方程,更準確地反映了主流的時均應變率[13]。

    RNG-模型

    式中,是湍動能,m2·s-2;是耗散率,m2·s-3;常數(shù)σ=1.0,σ=1.2,=0.09,1ε=1.44,2ε=1.9;G是由于層流速度梯度產(chǎn)生的湍流動能,b是由于浮力而產(chǎn)生的湍流動能;eff是有效動力黏度,Pa·s;是運動黏度,m2·s-1。

    汽油油氣是多組分混合物,傳質(zhì)過程較為復雜,為簡化物理模型,將油氣看成單一組分進行解算,氣相組分輸運標量方程為[14-15]

    eff是有效組分擴散系數(shù),m2·s-1,它等于分子擴散系數(shù)和湍流擴散系數(shù)之和

    式中,AB是通過富勒定律計算的二元擴散系數(shù),取8.9×10-6 m2·s-1;eff是有效動力黏度,Pa·s;t是紊流Schmidt數(shù)。根據(jù)單相傳質(zhì)規(guī)律,傳質(zhì)通量與有效擴散系數(shù)和濃度梯度有關[16]

    式中,是質(zhì)量通量,kg·(m2·s)-1;為油氣濃度,kg·m-3;c是飽和油氣質(zhì)量分數(shù);是氣液相界面面積,m2。

    1.2 物理模型

    油罐模型如圖1所示,油罐具體規(guī)格參數(shù)如下:罐直徑60 cm,壁高89.5 cm,總高為120 cm,高裝油口高度為117 cm,中裝油口高度為56.5 cm低裝油口高度為10.5 cm,鶴管直徑1.5 cm,罐頂通氣孔直徑5 cm。由于油罐屬于三維軸對稱模型,為便于計算,簡化為二維軸對稱模型來計算。為了使計算結果更精確,采用結構化網(wǎng)格進行劃分,并對油罐速度進口和壓力出口以及油罐底部汽油噴濺區(qū)域進行加密,各裝油口油罐網(wǎng)格數(shù)為280000左右,圖2為低裝油口附近網(wǎng)格分布情況。邊界條件分別為速度入口和壓力出口以及無滑移壁面邊界條件。

    由于本文模擬及實驗油罐尺寸較小,為考察驗證裝油速度、初始油氣濃度對油罐收油作業(yè)時油氣擴散排放規(guī)律影響的放大效應,進一步對大型拱頂油罐收油過程進行數(shù)值模擬研究,油罐具體規(guī)格參數(shù)如下:油罐內(nèi)直徑2370 cm,罐壁高1270 cm,油罐總高為1530 cm,裝油口直徑為20 cm,裝油口距罐底的距離為40 cm,罐頂通氣孔直徑為25 cm。

    2 數(shù)值模擬

    2.1 條件假設

    (1)由于裝油時間較短,忽略了溫度變化對油氣擴散排放的影響,將溫度設為常溫300 K。

    (2)由于汽油是一種復雜的混合物,包含許多不同的組分,如果計算汽油中的所有組分,將會產(chǎn)生巨大的計算量(事實上,目前還難以實現(xiàn)),所以為了減少計算量并側重研究汽油裝罐時罐內(nèi)流體的擾動及油氣濃度分布狀態(tài),假設裝油過程中氣液界面蒸發(fā)擴散的組分保持不變,并把油氣和空氣作為兩種組分。實際油罐收油作業(yè)時罐內(nèi)油品會隨著時間的變化,損耗量不斷增加,即蒸發(fā)出的油氣組分會隨著損耗量、溫度等其他影響因素的變化而變化。

    2.2 求解方法

    控制方程[式(1)~式(5)]用有限體積法進行離散,其中動量方程中對流項利用高階差分格式離散[17],擴散項通過上一時間步的計算值進行計算。為了計算氣液接觸面的曲率和表面張力,Brackbill等[18]曾使用B樣條插值函數(shù)對容積率進行光順處理,本文使用同樣的處理方式。連續(xù)性方程式(1)中容積率的值通過高階離散格式進行計算。

    輸運方程式(7)中的對流項采用二階迎風格式,非穩(wěn)態(tài)項采用一階顯示時間積分方式。輸運方程式(7)只在氣相進行計算,并且汽油的蒸發(fā)損耗率取決于氣液相界面的源項,由于源項是交界面組分梯度的函數(shù),而液相并不計算輸運方程,所以會導致交界面組分梯度的不連續(xù)性。為解決以上問題,需要定義一個中間變量用來計算交界面處的組分梯度[11]。中間變量的值在氣相中等于組分質(zhì)量分數(shù),在液相中等于平均質(zhì)量分數(shù)。具體通過以下步驟來計算汽油的蒸發(fā)損耗[7]:

    (1)通過VOF方程確定所有單元內(nèi)液體的體積分數(shù);

    (2)計算氣相傳輸特性;

    (3)在含有部分或全部液體的單元內(nèi)對中間變量指定飽和氣體質(zhì)量分數(shù);

    (4)在含有部分或全部氣體的單元內(nèi)對中間變量指定氣體質(zhì)量分數(shù)[式(7)];

    (5)確定每一單元體積內(nèi)的蒸發(fā)率[式(9)];

    (6)求解控制方程[式(1)~式(5)]。

    3 實驗驗證

    建立小型汽油裝罐蒸發(fā)損耗實驗平臺如圖3所示。汽油從工作油罐上部裝入,工作罐排放的氣體由氣體流量計計量。裝油口也同樣設置有高裝油口(高度為117 cm),中裝油口(高度為56.5 cm)、低裝油口(高度為10.5 cm),在工作罐內(nèi)部=40 cm、方向上,設置5個不同高度的油氣采樣點(T1~T5),如圖4。采樣點T1~T5(圖4)及總排放管的采樣點T6(圖3)采集的油氣-空氣混合氣采用帶有FID的島津氣相色譜儀2010-Plus進行油氣濃度分析。

    4 結果與分析

    4.1 罐內(nèi)油氣濃度及流場分布

    油罐收油作業(yè)時,裝油口高度對罐內(nèi)油氣濃度分布影響較大。圖5為裝油速度為0.89 m·s-1時不同高度裝油口在540 s時罐內(nèi)油氣濃度分布實驗值與模擬值對比(其中橫坐標為罐內(nèi)氣體空間點與罐底的距離,縱坐標為罐內(nèi)氣體空間點油氣濃度質(zhì)量分數(shù));圖6為不同高度裝油口在0.89m·s-1裝油速度下不同時刻的罐內(nèi)油氣濃度分布云圖,油罐初始油氣濃度為0。裝油口的裝油速度是根據(jù)油罐裝油流量及裝油口直徑計算出的油流速度。由圖6(a)可知,高裝油口裝油時,罐內(nèi)氣體空間油氣濃度較高,當=540 s時基本達到飽和狀態(tài),結合圖5可知,此時罐頂區(qū)域油氣質(zhì)量分數(shù)約為0.35。主要是由于高裝油口裝油時,汽油從裝油口噴濺流出,并快速分散成多種液滴,形成較大的蒸發(fā)面積,同時由于油流沖擊油面引起汽油的攪動作用,以及液面上升引起的罐內(nèi)氣體空間氣流流動,使汽油蒸發(fā)和對流擴散加劇,罐內(nèi)氣體空間油氣濃度很快趨于飽和。此時,氣體空間油氣-空氣之間的傳質(zhì)機理是由強制對流而實現(xiàn)的。由圖6(b)可知,中裝油口裝油時,進油口以下氣體空間油氣濃度值較高,主要是因為油品還沒有沒過裝油口前,罐內(nèi)氣體空間油氣-空氣傳質(zhì)以對流擴散為主,同時,由于汽油對罐內(nèi)液相的沖擊與攪拌作用以及罐內(nèi)氣體空間氣流的流動,使蒸發(fā)和對流擴散逐漸加劇。裝油口以上氣體空間受對流影響較小,所以濃度值偏低,由圖5可知,540 s時罐頂區(qū)域油氣質(zhì)量分數(shù)為0.15左右。此時,氣體空間油氣-空氣之間的傳質(zhì)機理是先由強制對流而實現(xiàn)的,后段(裝油口浸沒油品后)對流傳質(zhì)的強度逐漸減弱。圖6(c)為低裝油口在0.89 m·s-1裝油速度下罐內(nèi)不同時刻油氣濃度分布云圖,低裝油口裝油時,油品很快沒過進油口,罐內(nèi)氣體空間油氣-空氣之間傳質(zhì)是對流傳質(zhì)和分子擴散共存,因而油氣空間的濃度增加速度相對緩慢,總體上罐內(nèi)油氣濃度值不高。通過圖5數(shù)據(jù)對比可以看出模擬曲線變化趨勢與實驗值吻合良好,說明了數(shù)值模擬的準確性,并很好地揭示了裝油過程罐內(nèi)氣體空間油氣-空氣之間的傳質(zhì)機理。

    油罐收油作業(yè)時,隨著裝油速度增加,汽油蒸發(fā)和對流現(xiàn)象逐漸加劇,罐內(nèi)氣體空間油氣濃度越快到達飽和狀態(tài)。圖7為不同高度裝油口在1.2 m·s-1裝油速度下罐內(nèi)油氣濃度分布云圖,油罐初始油氣濃度為0。通過對比圖6可知,隨著裝油速度增大,相同液位時中、高裝油口罐內(nèi)油氣濃度值明顯增大,由于低裝油口裝油時受對流影響較小,由圖7(c)可知,隨著裝油速度增大,罐內(nèi)油氣濃度值變化較小。

    圖8為不同高度裝油口在裝油速度為1.2 m·s-1時不同時刻罐內(nèi)流場分布。由圖8可以看出,汽油從裝油口噴出后,通過擴散卷吸周圍空氣,使空氣沿軸中心位置向下方兩側擾動,到達罐壁后,沿罐壁向上方運移。當汽油到達罐底,與罐底表面發(fā)生撞擊,并在罐底軸中心兩側形成渦流,加劇了汽油蒸發(fā)速率,結合圖7可以看出,=30 s時罐底兩側位置濃度較高。隨著罐內(nèi)汽油液面的上升,液面以下位置由于不斷受到后面流體的擠壓及攪動而出現(xiàn)渦流,流體做湍流流動。液面以上部分由于液面上升引起的氣流變化及汽油噴射引起的卷吸作用,從而形成不同程度大小的渦流,結合圖7可以看出,近液面處渦流位置油氣濃度值偏高。從而,圖8也進一步印證了上述的裝油過程罐內(nèi)氣體空間油氣-空氣之間傳質(zhì)機理的分析。

    油罐裝油過程中,油品蒸發(fā)擴散排放規(guī)律會受罐內(nèi)初始油氣濃度影響。設定常溫下油罐初始油氣質(zhì)量分數(shù)為0、0.1、0.3,進油速度為1.5 m·s-1。圖9為中裝油口收油作業(yè)時罐內(nèi)油氣濃度隨油罐初始油氣濃度變化規(guī)律,可以看出隨初始油氣濃度增加,罐內(nèi)油氣濃度梯度越小,達到飽和狀態(tài)的時間越短。

    4.2 油罐排放口油氣濃度分布

    圖10為裝油速度在0.89 m·s-1時不同裝油口高度、不同初始油氣濃度條件下油罐排放口油氣濃度分布模擬和實驗對比。初始油氣濃度為0時,由圖中可以看出高裝油口油罐排放口油氣濃度迅速上升,約在濃度達到1.3 g·L-1之后,增速平緩,形成一個向上的拐點;中裝油口油罐排放口在0~5 min時間內(nèi)增速很小,濃度在0.08 g·L-1左右,隨后排放口油氣濃度逐漸增大,主要是油品浸沒裝油口之前以對流擴散傳質(zhì)為主,油品沒過裝油口之后排放口油氣濃度值約為1.2 g·L-1,之后罐內(nèi)油氣擴散傳質(zhì)以分子擴散為主,濃度曲線較為平緩。低裝油口裝油時,油氣-空氣傳質(zhì)主要以分子擴散為主,所以0~15 min內(nèi)排放口油氣濃度值較低,隨著罐內(nèi)液面的不斷上升,導致氣體空間變小,混合氣體受到壓縮,排放口的油氣濃度隨著裝油過程的進行而逐漸增加。油罐初始油氣濃度對油罐排放口油氣濃度分布影響較大。結合圖9可知,油罐在收油作業(yè)時,初始油氣濃度越大,罐內(nèi)油氣濃度梯度越小,達到飽和狀態(tài)時間越短。由圖10可以看出,初始油氣濃度設為0.3 g·L-1時,中、高裝油口油罐排放口油氣濃度比初始油氣濃度為0時較快達到飽和狀態(tài),低裝油口油罐排放口油氣濃度達到飽和狀態(tài)時間相差不大。

    4.3 裝油過程中氣液比變化規(guī)律

    圖11為油罐排放氣液比隨裝油速度、裝油口高度變化規(guī)律的實驗和模擬數(shù)據(jù)對比。從圖中可以看出,低裝油口裝油時,油罐排放氣液比隨裝油速度增加而略微減小,對比圖6、圖7可知,當裝油速度從0.89 m·s-1增大到1.2 m·s-1時罐內(nèi)油氣濃度增幅較小,并且由于油罐裝滿時間變短,排放口排氣量減少,所以氣液比變小[19-20]。中、高裝油口裝油時,對流傳質(zhì)逐漸占主導地位,從而對氣液比產(chǎn)生影響。由圖11可以看出高裝油口裝油時,氣液比隨裝油速度增大,先快速增大后緩慢減小,因為裝油速度越大,罐內(nèi)對流擴散及油品的附加蒸發(fā)逐漸加劇;裝油速度小于1.2 m·s-1時氣液比隨裝油速度增大而增大,隨著裝油速度繼續(xù)增大,雖然罐內(nèi)對流擴散及油品附加蒸發(fā)仍然很劇烈,但油罐裝滿時間變短,罐內(nèi)液面位置高濃度層油氣來不及擴散至排放口,所以裝油速度大于1.2 m·s-1時,油罐排放氣液比略有減小。

    圖12為油罐排放氣液比隨裝油口高度、初始油氣濃度變化規(guī)律模擬實驗對比??梢钥闯?,隨罐內(nèi)初始油氣濃度的增加值減小。由于初始油氣濃度的增加會減小組分質(zhì)量分數(shù)梯度從而抑制汽油的蒸發(fā)速率,導致汽油附加蒸發(fā)量減小,因此排放出的混合氣也相應減小。如果初始油氣濃度越接近飽和濃度,油罐排放氣液比越接近1。所以,當罐內(nèi)初始油氣濃度從0增加到0.9 g·L-1時,值逐漸減小。

    4.4 裝油過程中損耗率變化規(guī)律

    圖13、圖14分別為裝油損耗率隨裝油速度、初始油氣濃度變化規(guī)律實驗模擬數(shù)據(jù)對比。由圖13可以看出,裝油速度相同時,低裝油口裝油損耗率最低,約為0.025%。高裝油口裝油損耗率較大,最大值約為0.34%。中裝油口裝油損耗率居于高、低裝油口之間,裝油速度為1.5 m·s-1時約為0.1%,并且總體來看損耗率較低。由于高裝油口裝油時,液面一直在裝油口以下,罐內(nèi)油品蒸發(fā)擴散主要是油流對液面的攪動及罐內(nèi)氣體空間氣液兩相對流引起的,因此高裝油口裝油時損耗率較大。中裝油口裝油時,先以噴淋裝油為主,罐內(nèi)油品蒸發(fā)以對流擴散傳質(zhì)為主,油品沒過裝油口之后,蒸發(fā)損耗主要是由分子擴散及液面上升引起的。低裝油口裝油時,油品很快沒過裝油口,罐內(nèi)油品蒸發(fā)損耗以分子擴散為主,所以排放口油氣濃度較低,損耗量較小。由圖14可以看出,當罐內(nèi)初始油氣濃度從0增加到0.3 kg·m-3時,值也相應增加。1.8 m·s-1裝油速度下,高裝油口損耗率約為0.44%,中裝油口約為0.21%,低裝油口約為0.043%。雖然初始油氣濃度的增加會抑制汽油的蒸發(fā)速率,但油罐排放的油氣-空氣混合氣總累積量要大于清洗罐排放氣的總累積量,所以值會隨著罐內(nèi)初始油氣濃度的增加而增加。

    4.5 大型油罐模擬結果及分析

    為驗證小型油罐數(shù)值模擬計算結果的實際意義及放大效應,進一步對大型拱頂油罐收油過程進行數(shù)值模擬研究,但由于篇幅所限,具體研究成果擬另文發(fā)表,本文只簡要給出部分研究內(nèi)容與小型油罐數(shù)值模擬結果相比對。圖15為大型拱頂油罐在1.5 m·s-1裝油速度下3.8 h時罐壁油氣濃度分布圖及流場分布,由圖15(a)可以看出,=3.8 h時油氣剛剛擴散至罐頂氣體空間,但分布并不均勻,根據(jù)計算結果,此時排放口油氣濃度質(zhì)量分數(shù)約為0.003。圖15(b)為=3.8 h時罐內(nèi)流場分布,此時油品已經(jīng)沒過裝油口,罐內(nèi)氣體流動主要是由于油流攪動及液面上升產(chǎn)生的影響。結合圖15(a)可以看出,渦流處油氣濃度值相對較高。

    圖16為大型拱頂油罐損耗率模擬值與理論值對比,由圖中可以看出,模擬曲線變化趨勢與理論值吻合良好,說明了以上數(shù)值模擬方法的正確性并可應用于大型油罐對收油作業(yè)時油品擴散排放規(guī)律的研究,為實際生產(chǎn)應用中油品損耗評估提供參考。

    4.6 對API損耗評估公式的討論

    美國環(huán)保局2006年最新修正的固定頂油罐API收油損耗計算公式為[21]

    式中,w為年工作損失,kg·a-1;為油罐年周轉(zhuǎn)系數(shù);LX為最大液位高度,m;為油罐直徑,m;N為工作損失周轉(zhuǎn)(飽和)因子;P為油品系數(shù);V為油氣密度,kg·m-3;B為呼吸閥設定校正因子,呼吸閥開啟壓力設置在±207 Pa范圍內(nèi)時,B=1。

    由式(10)可以看出,固定頂油罐收油損耗與油罐直徑及罐內(nèi)液位高度相關,但式中并未給出裝油速度和初始油氣濃度對收油損耗的影響,而根據(jù)模擬結果顯示,裝油速度以及初始油氣濃度對油品損耗率有較大影響,尤其是中、高裝油口收油過程。因此,建議API固定頂油罐大呼吸損耗計算公式考慮裝油速度和油罐初始油氣濃度對蒸發(fā)損耗的影響。

    5 結 論

    通過CFD軟件綜合計算連續(xù)性方程、動量方程、質(zhì)量傳遞方程,并結合實驗實測數(shù)據(jù),根據(jù)不同裝油口高度、不同裝油速度、不同初始油氣濃度來研究汽油從噴濺式裝油產(chǎn)生蒸發(fā)損耗的過程。

    (1)增加裝油速度,中、低裝油口油罐排放氣液比均減小,高裝油口油罐排放氣液比先增大后略減??;增大油罐初始油氣濃度,高、中、低裝油口油罐排放氣液比均減小,因為初始油氣濃度的增加會減小組分質(zhì)量分數(shù)梯度從而抑制汽油的蒸發(fā)速率,導致汽油附加蒸發(fā)量減小,因此排放出的混合氣也相應減小。

    (2)裝油口位置越高,汽油損耗率越大,隨裝油速度增大,高裝油口損耗率最大約為0.34%,低裝油口約為0.025%。增大油罐初始油氣濃度,高裝油口損耗率最大約為0.44%,中裝油口約為0.21%,低裝油口約為0.043%。

    (3)通過對大型拱頂油罐數(shù)值模擬研究發(fā)現(xiàn),油罐收油損耗率與理論計算值吻合良好,說明本數(shù)值模擬方法同樣適用于大型拱頂油罐,并可為油罐收油損耗評估提供一定參考。

    (4)裝油速度、油罐初始油氣濃度對排放口混合氣濃度、損耗率、氣液比等影響較大,建議固定頂油罐大呼吸API計算公式考慮裝油速度及油罐初始油氣濃度的影響,裝油過程中根據(jù)裝油口高度適當調(diào)整裝油速度大小并盡量使用清洗罐裝油。

    [1] 黃維秋, 趙書華, 呂愛華, 等. 油氣回收基礎理論及其應用[M]. 北京: 中國石化出版社, 2011: 89. HUANG W Q, ZHAO S H, Lü A H,. Fundamental Theory of Oil Vapor Recovery and Its Application[M]. Beijing: China Petrochemical Press, 2011: 89.

    [2] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局, 環(huán)境保護部. 石油煉制工業(yè)污染物排放標準: GB 31570—2015[S]. 北京: 中國環(huán)境科學出版社, 2015. General Administration of Quality Supervision, Environment Protection Department. Emission standard of pollutants for petroleum refining industry : GB 31570—2015[S]. Beijing: China Environment Science Press, 2015.

    [3] 中華人民共和國住房和城鄉(xiāng)建設部. 石油庫設計規(guī)范: GB 50074—2014[S]. 北京: 中國計劃出版社, 2015. Ministry of Housing and Urban-Rural Construction of the People’s Republic of China. Code for design of oil depot: GB 50074—2014[S]. Beijing: China Planning Press, 2015.

    [4] 黃維秋, 饒原剛, 呂艷麗. 固定頂油罐與油氣回收集成工藝研究[J]. 油氣儲運, 2009, 28(11): 20-23, 81, 4. HUANG W Q, RAO Y G, Lü Y L. Research on integrated process of fixed roof tanks with vapor recovery system[J]. Oil & Gas Storage and Transportation, 2009, 28(11): 20-23, 81, 4.

    [5] YAMADA H. Contribution of evaporative emissions from gasoline vehicles toward total VOC emissions in Japan[J]. Science of the Total Environment, 2013, 449: 143-149.

    [6] KOYNOV A, KHINAST J G. Micromixing in reactive, deformable bubble and droplet swarms[J]. Chemical Engineering & Technology, 2006, 29(1): 13-23.

    [7] DAVIDSON M R, RUDMAN M. Volume of fluid calculations of heat or mass transfer across deforming interfaces in two-fluid flow[J]. Numer. Heat Transfer B, 2002, 41(3/4): 291-308.

    [8] WAHEED M A, HENSCHKE M, ANDREAS P. Mass transfer by free and forced convection from single spherical liquid drops[J]. International Journal of Heat and Mass Transfer, 2002, 45(22): 4507-4514.

    [9] HAELSSIG J B, ANDRE Y T, THIBAULT J,. Direct numerical simulation of interphase heat and mass transfer in multicomponent vapor-liquid flows[J]. International Journal of Heat and Mass Transfer, 2010, 53(19/20): 3947-3960.

    [10] 朱玲, 陳家慶, 柳巖, 等. 汽油揮發(fā)過程的試驗分析[J]. 車用發(fā)動機, 2011, (2): 68-72. ZHU L, CHEN J Q, LIU Y,. Experimental analysis of gasoline evaporation process[J]. Vehicle Engine, 2011, (2): 68-72.

    [11] HASSANVAND A, HASHEMABADI S H, BAYAT M. Evaluation of gasoline evaporation during the tank splash loading by CFD techniques[J]. International Communications in Heat and Mass Transfer, 2010, 37(7): 907-913.

    [12] BANERJEE R, ISAAC K M. Evaluation of turbulence closure schemes for stratified two-phase flow[C]//2003 ASME International Mechanical Engineering Congress and Exposition. Washington DC, 2003.

    [13] LATEB M, MASSON C, STATHOPOULOS T,. Comparison of various types of-models for pollutant emissions around a two-building configuration[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2013, (115): 9-21.

    [14] 趙晨露, 黃維秋, 鐘璟, 等. 外浮頂罐油氣泄漏的數(shù)值模擬[J]. 化工學報, 2014, 65(10): 4203-4209. ZHAO C L, HUANG W Q, ZHONG J,. Numerical simulation of oil vapor leakage from external floating-roof tank[J]. CIESC Jorunal, 2014, 65(10): 4203-4209.

    [15] 趙晨露, 黃維秋, 石莉, 等. 內(nèi)浮頂罐中油氣擴散運移的數(shù)值模擬[J]. 安全與環(huán)境學報, 2015, 15(3): 72-77. ZHAO C L, HUANG W Q, SHI L,. Numerical simulation of oil vapor diffusion and migration from internal floating-roof tank[J]. Journal of Safety and Environment, 2015, 15(3): 72-77.

    [16] KAYS W, CRAWFORD M, WEIGAND B. Convective Heat and Mass Transfer [M]. New York: McGraw-Hill, 2005: 121.

    [17] SAMIMI ABIANEH O, CHEN C P, MAHALINGAM S. Modeling of multi-component droplets collision in dense evaporating and non-evaporating sprays[J]. International Journal of Automotive Technology, 2014, 15(7): 1091-1100.

    [18] BRACKBILL J U, KOTHE D B, ZEMACH C. A continuum method for modeling surface tension[J]. Journal of Computational Physics, 1992, 100(2): 335–354.

    [19] 黃維秋, 鐘秦. 非穩(wěn)態(tài)油品蒸發(fā)的數(shù)值分析(Ⅰ)數(shù)學模型的建立[J]. 石油學報(石油加工), 2003, 19(6): 58-63. HUANG W Q, ZHONG Q. Numerical analysis of petroleum products evaporation under unsteady conditions(Ⅰ): Establishment of mathematic model[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2003, 19(6): 58-63.

    [20] 黃維秋, 鐘秦. 非穩(wěn)態(tài)油品蒸發(fā)的數(shù)值分析(Ⅱ)數(shù)值分析及應用實例[J]. 石油學報(石油加工), 2004, 20(1): 52-57. HUANG W Q, ZHONG Q. Numerical analysis of petroleum products evaporation under the unsteady conditions(Ⅱ):Numerical analysis and application examples[J]. Acta Petrolei Sinica (Petroleum Processing Section), 2004, 20(1): 52-57.

    [21] U. S. Environmental Protection Agency Office. AP 42 Fifth Edition, Volume I, Chapter: Liquid Storage Tanks[EB/OL]. https://www3.epa.gov/ttn/chief/ap42/ch07/final/c07s01.pdf.

    Experimental determination and numerical simulation of vapor diffusion and emission in loading gasoline into tank

    HUANG Weiqiu1, WANG Zhaoli1, JI Hong1, ZHAO Chenlu2, Lü Aihua1, XU Xianyang1, WANG Yihong1

    (1Jiangsu Key Laboratory of Oil & Gas Storage and Transportation Technology (Changzhou University), Changzhou 213016, Jiangsu, China; 2Shanghai Hoto Engineering Inc., Shanghai 200000, China)

    It was of great significance correctly to predict the vapor diffusion and emission inside the tank for the researches of gasoline evaporation loss and vapor pollution control during the operation of loading gasoline into a tank. Then, two key parameters of the volumetric ratioof the displacement mixture gas of the vapor-air to the loaded gasoline and the qualitative ratioof the evaporation loss to the loaded gasoline were mainly considered. Based on the models of volume of fluid (VOF), mass transfer and RNG-turbulence, the evaporation loss was numerically simulated and experimentally investigated during the splash loading operation, and the oil vapor-air diffusion was analyzed and compared for the different loading exit heights, different loading velocity and the different initial vapor mass fraction. Meanwhile, an experimental system of evaporation loss in loading into a tank was built up to verify the numerical simulation, and the results of the numerical simulation were agreed well with the experimental data. The simulation results furthermore showed that the higher of the loading exit, the greater the qualitative ratioof the evaporation loss to the loaded gasoline. The qualitative ratio of high exit was at around 0.34% and the qualitative ratio of low exit at around 0.025% by the increase of the loading velocity. The qualitative ratioof high exit was at around of 0.44%, the qualitative ratio of mid exit at around 0.21% and the qualitative ratio of low exit at around 0.043% by increasing the initial vapor mass fraction. It was recommended that the effect of loading velocity and the initial vapor mass fraction should be considered in API loss formula by using a clean tank and low exit when loading and reducing loading speed appropriately before the loading pipe exit was submerged.

    gasoline; evaporation; diffusion; loading velocity; initial vapor mass fraction; numerical simulation

    date: 2016-04-26.

    Prof. HUANG Weiqiu, hwq213@cczu.edu.cn

    10.11949/j.issn.0438-1157.20160550

    TE 85

    A

    0438—1157(2016)12—4994—12

    國家自然科學基金項目(51574044);江蘇省高?!扒嗨{工程”資助項目(SCZ1409700002);江蘇省油氣儲運重點實驗室項目(SCZ1211200004/004)。

    supported by the National Natural Science Foundation of China (51574044), the Qing Lan Project of the Jiangsu Higher Education Institutions of China (SCZ1409700002) and the Major Research Plan of the Oil and Gas Storage and Transportation Laboratory of Jiangsu Province (SCZ1211200004/004).

    2016-04-26收到初稿,2016-09-19收到修改稿。

    聯(lián)系人及第一作者:黃維秋(1965—),男,教授。

    猜你喜歡
    油口損耗率油罐
    高溫重油罐區(qū)超壓分析總結
    云南化工(2021年6期)2021-12-21 07:31:30
    一種全液壓助力制動叉車液壓系統(tǒng)
    安徽科技(2021年3期)2021-04-06 08:50:16
    把握主動權,提高油罐火災撲救成功率
    水上消防(2020年4期)2021-01-04 01:51:20
    損耗率高達30%,保命就是保收益!這條70萬噸的魚要如何破存活率困局?
    淺談鋼煤斗制作原材料損耗率的研究
    Clinical outcomes of endoscopic management of pancreatic fluid collections in cirrhotics vs non-cirrhotics: Α
    脫水撬脫水效果分析與評價
    KD424:小型挖掘機用泵
    某分餾廠油罐區(qū)設計
    帶有轉(zhuǎn)閥換向機構的手動油泵
    科技資訊(2016年23期)2016-05-30 18:00:34
    熟女电影av网| 色综合婷婷激情| 禁无遮挡网站| 久久精品91无色码中文字幕| 性色av乱码一区二区三区2| 美女高潮喷水抽搐中文字幕| 91av网一区二区| 色在线成人网| 亚洲精品一区av在线观看| 久久久久久久久免费视频了| 亚洲国产精品久久男人天堂| 真人一进一出gif抽搐免费| 欧美xxxx黑人xx丫x性爽| 国产精品免费一区二区三区在线| 亚洲精品一区av在线观看| 亚洲片人在线观看| 精华霜和精华液先用哪个| 欧美日本亚洲视频在线播放| 人人妻,人人澡人人爽秒播| 国产精品98久久久久久宅男小说| 超碰成人久久| www日本在线高清视频| 熟妇人妻久久中文字幕3abv| 91在线观看av| 午夜日韩欧美国产| 国产熟女xx| 国产精品亚洲美女久久久| 无限看片的www在线观看| 搡老熟女国产l中国老女人| 中文资源天堂在线| 亚洲国产日韩欧美精品在线观看 | 亚洲国产中文字幕在线视频| 一卡2卡三卡四卡精品乱码亚洲| 黄色丝袜av网址大全| 久久精品国产99精品国产亚洲性色| 国产成人欧美在线观看| 久久精品国产99精品国产亚洲性色| 欧美黄色淫秽网站| 欧美黄色淫秽网站| 欧美中文日本在线观看视频| 可以在线观看的亚洲视频| 欧美丝袜亚洲另类 | 香蕉av资源在线| av国产免费在线观看| 久久精品国产亚洲av香蕉五月| 他把我摸到了高潮在线观看| 91麻豆精品激情在线观看国产| 精品乱码久久久久久99久播| 热99在线观看视频| 国产91精品成人一区二区三区| 国产精品一及| 亚洲av成人av| 亚洲精品国产精品久久久不卡| 亚洲精品粉嫩美女一区| 99久久精品热视频| 美女大奶头视频| 亚洲美女视频黄频| 色综合亚洲欧美另类图片| 精品久久久久久,| www.自偷自拍.com| av国产免费在线观看| 男女下面进入的视频免费午夜| 国产三级在线视频| 日本与韩国留学比较| 999久久久国产精品视频| 不卡一级毛片| 大型黄色视频在线免费观看| 日韩成人在线观看一区二区三区| 免费看a级黄色片| bbb黄色大片| 精品无人区乱码1区二区| 狂野欧美激情性xxxx| 成熟少妇高潮喷水视频| 人妻久久中文字幕网| 99久久成人亚洲精品观看| 男女那种视频在线观看| 中文字幕av在线有码专区| 成年版毛片免费区| 中文亚洲av片在线观看爽| 1000部很黄的大片| 国产三级在线视频| 精品国产美女av久久久久小说| АⅤ资源中文在线天堂| 中文字幕人妻丝袜一区二区| 精品久久久久久久毛片微露脸| 一个人免费在线观看电影 | 丁香六月欧美| 久久久久久久久中文| 久久精品国产综合久久久| 欧美绝顶高潮抽搐喷水| 九九在线视频观看精品| 91字幕亚洲| 久久国产乱子伦精品免费另类| 亚洲国产精品sss在线观看| 在线永久观看黄色视频| 18禁国产床啪视频网站| 国产成人福利小说| 搡老妇女老女人老熟妇| 丝袜人妻中文字幕| 最近在线观看免费完整版| 国产成人欧美在线观看| 观看免费一级毛片| 看片在线看免费视频| 日本熟妇午夜| 国产精品亚洲一级av第二区| 国产乱人视频| av国产免费在线观看| 1024香蕉在线观看| 女生性感内裤真人,穿戴方法视频| 国产精品女同一区二区软件 | 久久久久国内视频| 国产成人一区二区三区免费视频网站| 99久久综合精品五月天人人| 日本一本二区三区精品| 亚洲自偷自拍图片 自拍| 精品不卡国产一区二区三区| 97碰自拍视频| 男女做爰动态图高潮gif福利片| 欧美黄色片欧美黄色片| 国产精品女同一区二区软件 | aaaaa片日本免费| xxxwww97欧美| 制服人妻中文乱码| 国产精品野战在线观看| 欧美乱色亚洲激情| 变态另类成人亚洲欧美熟女| 18禁观看日本| 在线视频色国产色| 亚洲av免费在线观看| 久久久国产精品麻豆| 51午夜福利影视在线观看| 一边摸一边抽搐一进一小说| 亚洲成人久久爱视频| 成人三级做爰电影| 国产av麻豆久久久久久久| 在线观看舔阴道视频| 熟女电影av网| 欧美av亚洲av综合av国产av| 国产欧美日韩精品亚洲av| 69av精品久久久久久| 精品欧美国产一区二区三| 熟女少妇亚洲综合色aaa.| 亚洲人成网站高清观看| 国产探花在线观看一区二区| 欧美乱码精品一区二区三区| 手机成人av网站| 不卡一级毛片| 最新中文字幕久久久久 | 亚洲成av人片免费观看| 女人被狂操c到高潮| 黑人欧美特级aaaaaa片| 成人永久免费在线观看视频| 欧美日韩综合久久久久久 | 亚洲一区二区三区不卡视频| 国产精品亚洲一级av第二区| 国产精品久久久久久亚洲av鲁大| 国产一区在线观看成人免费| 精品乱码久久久久久99久播| 综合色av麻豆| 国产成人精品久久二区二区91| 黄色片一级片一级黄色片| 伦理电影免费视频| 小蜜桃在线观看免费完整版高清| 色老头精品视频在线观看| 动漫黄色视频在线观看| 啦啦啦免费观看视频1| 床上黄色一级片| 午夜免费激情av| 欧美另类亚洲清纯唯美| 99久久99久久久精品蜜桃| 国产极品精品免费视频能看的| 在线视频色国产色| 宅男免费午夜| 欧美一区二区国产精品久久精品| 亚洲av第一区精品v没综合| 精品欧美国产一区二区三| 美女扒开内裤让男人捅视频| 99国产精品99久久久久| 亚洲熟女毛片儿| www日本黄色视频网| 九九热线精品视视频播放| 99久久无色码亚洲精品果冻| 成人av在线播放网站| 国产精品永久免费网站| 免费观看精品视频网站| xxx96com| 搡老熟女国产l中国老女人| 超碰成人久久| 亚洲精品一卡2卡三卡4卡5卡| 97碰自拍视频| 亚洲av成人av| 久久久久国内视频| 91麻豆av在线| 国产激情欧美一区二区| 免费av不卡在线播放| 91av网一区二区| 免费电影在线观看免费观看| 亚洲欧美一区二区三区黑人| 两性夫妻黄色片| 色尼玛亚洲综合影院| 亚洲最大成人中文| 最好的美女福利视频网| 国产欧美日韩精品亚洲av| 日日干狠狠操夜夜爽| 在线免费观看不下载黄p国产 | 中文字幕av在线有码专区| 一区二区三区激情视频| 亚洲国产看品久久| 久久中文字幕人妻熟女| 国产成人aa在线观看| 国产高清videossex| 老司机深夜福利视频在线观看| 999久久久国产精品视频| 日本精品一区二区三区蜜桃| 一个人免费在线观看电影 | 国产高清视频在线播放一区| 97超视频在线观看视频| 国产精品女同一区二区软件 | bbb黄色大片| 中文字幕最新亚洲高清| 无限看片的www在线观看| 日本免费一区二区三区高清不卡| 日韩成人在线观看一区二区三区| 亚洲国产欧美人成| 一进一出抽搐动态| 欧美一级a爱片免费观看看| 国语自产精品视频在线第100页| 成人国产一区最新在线观看| 亚洲av美国av| 这个男人来自地球电影免费观看| 亚洲一区高清亚洲精品| www.熟女人妻精品国产| 真人做人爱边吃奶动态| 国产99白浆流出| 草草在线视频免费看| 精品免费久久久久久久清纯| 国产成人欧美在线观看| 欧美日本视频| 1000部很黄的大片| 久久久久久久久久黄片| 国产精品久久久久久人妻精品电影| 久久久国产成人精品二区| 少妇的逼水好多| 一区二区三区激情视频| 亚洲美女视频黄频| 亚洲黑人精品在线| 三级毛片av免费| 中文字幕av在线有码专区| www.自偷自拍.com| 久久久久久国产a免费观看| 日本一二三区视频观看| 国产高清三级在线| 国产精品亚洲一级av第二区| 国产男靠女视频免费网站| 波多野结衣高清作品| 亚洲五月婷婷丁香| 亚洲avbb在线观看| 亚洲欧美日韩东京热| 国产三级在线视频| 精品国产乱码久久久久久男人| aaaaa片日本免费| 在线国产一区二区在线| 熟女电影av网| 九九久久精品国产亚洲av麻豆 | 97超级碰碰碰精品色视频在线观看| 两个人看的免费小视频| 熟女电影av网| 日韩欧美三级三区| a在线观看视频网站| 精品久久久久久久久久免费视频| 国产午夜福利久久久久久| 欧美黑人欧美精品刺激| 国产精品 国内视频| 天堂网av新在线| 国产精品影院久久| 亚洲精华国产精华精| 国产探花在线观看一区二区| 18禁黄网站禁片免费观看直播| 久久天躁狠狠躁夜夜2o2o| 日日摸夜夜添夜夜添小说| 国产野战对白在线观看| 国内毛片毛片毛片毛片毛片| 后天国语完整版免费观看| 欧美丝袜亚洲另类 | 久久香蕉精品热| 天堂动漫精品| 偷拍熟女少妇极品色| 国产成人aa在线观看| 狂野欧美白嫩少妇大欣赏| 免费一级毛片在线播放高清视频| 国产极品精品免费视频能看的| 国产欧美日韩精品亚洲av| 亚洲国产色片| 午夜精品一区二区三区免费看| 亚洲 欧美 日韩 在线 免费| 欧美黄色片欧美黄色片| 好男人在线观看高清免费视频| 亚洲 欧美一区二区三区| 久久香蕉精品热| 亚洲18禁久久av| 一本一本综合久久| 看免费av毛片| 久久久久久人人人人人| 国产精品自产拍在线观看55亚洲| 国内精品久久久久久久电影| 韩国av一区二区三区四区| 亚洲乱码一区二区免费版| 一个人免费在线观看电影 | www日本黄色视频网| 69av精品久久久久久| 丰满人妻熟妇乱又伦精品不卡| 国产99白浆流出| 一个人免费在线观看的高清视频| 桃色一区二区三区在线观看| 国产亚洲精品av在线| 巨乳人妻的诱惑在线观看| www.www免费av| xxx96com| 精品人妻1区二区| av视频在线观看入口| 国产精品 欧美亚洲| 九九在线视频观看精品| 亚洲精品色激情综合| 欧美日韩一级在线毛片| 欧美乱码精品一区二区三区| 国产精品久久久av美女十八| 男女之事视频高清在线观看| 色尼玛亚洲综合影院| 久久久久久人人人人人| 极品教师在线免费播放| 精品午夜福利视频在线观看一区| 免费无遮挡裸体视频| 久久久久久久久久黄片| 精品欧美国产一区二区三| 中文字幕av在线有码专区| 国产成人系列免费观看| 一个人观看的视频www高清免费观看 | 欧美日韩综合久久久久久 | ponron亚洲| 叶爱在线成人免费视频播放| 蜜桃久久精品国产亚洲av| 一个人看的www免费观看视频| 1024香蕉在线观看| 久久精品91无色码中文字幕| 真实男女啪啪啪动态图| 成人国产一区最新在线观看| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清| 国产精品久久久久久亚洲av鲁大| 久久久久久久精品吃奶| 最近最新中文字幕大全电影3| 亚洲中文字幕一区二区三区有码在线看 | 久久精品aⅴ一区二区三区四区| 久久久久久国产a免费观看| 成年女人毛片免费观看观看9| 亚洲专区国产一区二区| 岛国视频午夜一区免费看| bbb黄色大片| 白带黄色成豆腐渣| 久久99热这里只有精品18| 91久久精品国产一区二区成人 | 色播亚洲综合网| 亚洲 欧美一区二区三区| av国产免费在线观看| 久久午夜亚洲精品久久| 国产伦精品一区二区三区四那| 日韩成人在线观看一区二区三区| 这个男人来自地球电影免费观看| 天堂av国产一区二区熟女人妻| 两个人的视频大全免费| av国产免费在线观看| 中文字幕精品亚洲无线码一区| 色老头精品视频在线观看| 久久久色成人| 夜夜爽天天搞| 麻豆av在线久日| 男人舔女人下体高潮全视频| 亚洲国产精品sss在线观看| 亚洲片人在线观看| 中国美女看黄片| 99国产精品99久久久久| 18禁黄网站禁片免费观看直播| 国产亚洲av嫩草精品影院| 久久久久久久久中文| 男女下面进入的视频免费午夜| 亚洲五月天丁香| 成人欧美大片| 亚洲国产精品sss在线观看| 久久性视频一级片| netflix在线观看网站| 国产伦精品一区二区三区视频9 | 宅男免费午夜| 桃色一区二区三区在线观看| 三级男女做爰猛烈吃奶摸视频| a级毛片a级免费在线| 午夜成年电影在线免费观看| 日韩免费av在线播放| 少妇裸体淫交视频免费看高清| 国产又黄又爽又无遮挡在线| 999精品在线视频| 久久久久久久精品吃奶| 看免费av毛片| 一级a爱片免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品av在线| 国产蜜桃级精品一区二区三区| 国产淫片久久久久久久久 | 日韩免费av在线播放| 69av精品久久久久久| 看黄色毛片网站| 97超视频在线观看视频| 国产成人欧美在线观看| 亚洲国产精品sss在线观看| 欧美在线一区亚洲| 成人欧美大片| 男女做爰动态图高潮gif福利片| 中文字幕精品亚洲无线码一区| 精品国产超薄肉色丝袜足j| 欧美在线一区亚洲| 久久精品国产综合久久久| 18禁国产床啪视频网站| 久久精品国产99精品国产亚洲性色| 99久久久亚洲精品蜜臀av| 色精品久久人妻99蜜桃| 啦啦啦韩国在线观看视频| 婷婷六月久久综合丁香| 国产成人aa在线观看| 国产精品久久久久久久电影 | 久久久国产精品麻豆| 两个人看的免费小视频| 日本与韩国留学比较| 久久久国产成人免费| 亚洲欧美精品综合一区二区三区| 99在线人妻在线中文字幕| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 一边摸一边抽搐一进一小说| 国产伦在线观看视频一区| 一进一出好大好爽视频| 99久久精品国产亚洲精品| 亚洲精品在线美女| 亚洲黑人精品在线| 91九色精品人成在线观看| 深夜精品福利| 人妻久久中文字幕网| 我的老师免费观看完整版| 母亲3免费完整高清在线观看| 国产1区2区3区精品| 免费电影在线观看免费观看| 中出人妻视频一区二区| 一进一出抽搐gif免费好疼| 久久久久九九精品影院| 久久精品国产清高在天天线| 亚洲一区二区三区色噜噜| 一夜夜www| 国产aⅴ精品一区二区三区波| 一a级毛片在线观看| 久久性视频一级片| 婷婷亚洲欧美| 国产伦人伦偷精品视频| 精品久久久久久久人妻蜜臀av| 少妇丰满av| 黄色片一级片一级黄色片| 精品国产亚洲在线| 99国产精品一区二区蜜桃av| 一本综合久久免费| 视频区欧美日本亚洲| 精品一区二区三区av网在线观看| 亚洲片人在线观看| 欧美一级毛片孕妇| 欧美日本亚洲视频在线播放| 狂野欧美白嫩少妇大欣赏| 精品国产超薄肉色丝袜足j| 亚洲五月天丁香| 国内精品一区二区在线观看| 嫩草影院精品99| 成人欧美大片| 国产成年人精品一区二区| 狂野欧美激情性xxxx| 久久精品人妻少妇| 给我免费播放毛片高清在线观看| 国产精品一区二区三区四区免费观看 | 色尼玛亚洲综合影院| 天天躁狠狠躁夜夜躁狠狠躁| 美女cb高潮喷水在线观看 | 亚洲欧美一区二区三区黑人| 亚洲av成人一区二区三| 久久久国产精品麻豆| 丁香六月欧美| 变态另类丝袜制服| 精品一区二区三区av网在线观看| 97超视频在线观看视频| 黄色丝袜av网址大全| 日本 av在线| 日本撒尿小便嘘嘘汇集6| 国产av麻豆久久久久久久| 国产成年人精品一区二区| 级片在线观看| 久久99热这里只有精品18| www国产在线视频色| 亚洲美女黄片视频| xxxwww97欧美| 男人舔女人的私密视频| 一本一本综合久久| 成人特级黄色片久久久久久久| 嫁个100分男人电影在线观看| 亚洲欧美精品综合一区二区三区| 免费搜索国产男女视频| 欧美激情久久久久久爽电影| 岛国在线免费视频观看| netflix在线观看网站| 久久久久久大精品| 看黄色毛片网站| 成人特级av手机在线观看| 脱女人内裤的视频| 丰满人妻一区二区三区视频av | 欧美成狂野欧美在线观看| 一区二区三区激情视频| 嫁个100分男人电影在线观看| 又紧又爽又黄一区二区| 亚洲精华国产精华精| 老熟妇乱子伦视频在线观看| 性色av乱码一区二区三区2| 99久久99久久久精品蜜桃| 亚洲av电影在线进入| 99国产精品99久久久久| 国产精品日韩av在线免费观看| 中文字幕最新亚洲高清| 国产精品一及| 亚洲第一电影网av| 久久精品国产清高在天天线| 国产久久久一区二区三区| av欧美777| 99久久成人亚洲精品观看| 精华霜和精华液先用哪个| 成人三级做爰电影| 国产不卡一卡二| 欧美日韩亚洲国产一区二区在线观看| 欧美又色又爽又黄视频| x7x7x7水蜜桃| 午夜福利18| xxxwww97欧美| 午夜福利18| 欧美一级毛片孕妇| 美女午夜性视频免费| 成人三级做爰电影| 久久精品国产亚洲av香蕉五月| 怎么达到女性高潮| 岛国视频午夜一区免费看| 成人一区二区视频在线观看| 久久精品综合一区二区三区| 午夜精品一区二区三区免费看| 深夜精品福利| 亚洲在线自拍视频| 床上黄色一级片| 久久久精品大字幕| 亚洲av成人av| 免费一级毛片在线播放高清视频| 女生性感内裤真人,穿戴方法视频| 可以在线观看毛片的网站| 综合色av麻豆| 精品久久久久久久人妻蜜臀av| 久久久精品欧美日韩精品| 9191精品国产免费久久| 我要搜黄色片| 老熟妇乱子伦视频在线观看| 日本 欧美在线| 亚洲专区国产一区二区| 免费观看的影片在线观看| av在线蜜桃| 国产成人一区二区三区免费视频网站| ponron亚洲| 亚洲一区高清亚洲精品| 免费看a级黄色片| 亚洲色图 男人天堂 中文字幕| 我的老师免费观看完整版| 欧美乱妇无乱码| 国产毛片a区久久久久| 舔av片在线| 91在线观看av| 偷拍熟女少妇极品色| 欧美不卡视频在线免费观看| 欧美激情在线99| 婷婷精品国产亚洲av| 国产精华一区二区三区| ponron亚洲| 很黄的视频免费| 色在线成人网| www.精华液| 久久久久久久久久黄片| 天堂av国产一区二区熟女人妻| 亚洲九九香蕉| 亚洲av五月六月丁香网| 高清毛片免费观看视频网站| 我的老师免费观看完整版| 日本撒尿小便嘘嘘汇集6| 巨乳人妻的诱惑在线观看| 亚洲精品乱码久久久v下载方式 | 国产三级黄色录像| 亚洲av片天天在线观看| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| 老鸭窝网址在线观看| 日本成人三级电影网站| av女优亚洲男人天堂 | 首页视频小说图片口味搜索| 在线观看一区二区三区| 国产综合懂色| 综合色av麻豆| 欧美精品啪啪一区二区三区| 十八禁人妻一区二区| av福利片在线观看| av女优亚洲男人天堂 | 少妇的丰满在线观看| 久久亚洲精品不卡|