• <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
    国产男靠女视频免费网站| 亚洲熟妇中文字幕五十中出| 色精品久久人妻99蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 一级毛片高清免费大全| 12—13女人毛片做爰片一| 久久精品国产自在天天线| 午夜福利免费观看在线| 精品熟女少妇八av免费久了| 成人一区二区视频在线观看| 亚洲专区中文字幕在线| 国产三级黄色录像| 久久久国产成人免费| 观看免费一级毛片| 夜夜看夜夜爽夜夜摸| 丰满的人妻完整版| 午夜免费激情av| 欧洲精品卡2卡3卡4卡5卡区| 在线a可以看的网站| 亚洲精品成人久久久久久| 欧美性感艳星| 亚洲精品456在线播放app | 1024手机看黄色片| 成人一区二区视频在线观看| 97碰自拍视频| 国内毛片毛片毛片毛片毛片| 老司机午夜十八禁免费视频| 亚洲成av人片在线播放无| 国内少妇人妻偷人精品xxx网站| 99热只有精品国产| av女优亚洲男人天堂| 国产在线精品亚洲第一网站| 国产精品综合久久久久久久免费| 国产一区二区激情短视频| 久久久色成人| 99精品在免费线老司机午夜| 欧美+亚洲+日韩+国产| 亚洲精品国产精品久久久不卡| 免费观看的影片在线观看| 国内久久婷婷六月综合欲色啪| 欧美绝顶高潮抽搐喷水| 国产不卡一卡二| 狠狠狠狠99中文字幕| 很黄的视频免费| 欧洲精品卡2卡3卡4卡5卡区| 女同久久另类99精品国产91| 亚洲专区中文字幕在线| avwww免费| 琪琪午夜伦伦电影理论片6080| 亚洲狠狠婷婷综合久久图片| 美女高潮喷水抽搐中文字幕| 综合色av麻豆| 欧美一区二区国产精品久久精品| 丝袜美腿在线中文| 亚洲成av人片免费观看| 在线免费观看不下载黄p国产 | 日韩欧美一区二区三区在线观看| 亚洲国产欧洲综合997久久,| 日韩高清综合在线| 特大巨黑吊av在线直播| 精品久久久久久久久久免费视频| 美女被艹到高潮喷水动态| 欧美乱码精品一区二区三区| 久久久久久久亚洲中文字幕 | 精品国内亚洲2022精品成人| netflix在线观看网站| 国产伦一二天堂av在线观看| 老司机在亚洲福利影院| 日韩免费av在线播放| 国产男靠女视频免费网站| 一区福利在线观看| 免费高清视频大片| 老汉色av国产亚洲站长工具| 有码 亚洲区| 亚洲黑人精品在线| 熟女人妻精品中文字幕| 欧美成人免费av一区二区三区| 国产激情欧美一区二区| 国产三级在线视频| 欧美一区二区亚洲| 美女高潮的动态| 欧美绝顶高潮抽搐喷水| 91麻豆精品激情在线观看国产| a在线观看视频网站| 真实男女啪啪啪动态图| 神马国产精品三级电影在线观看| av视频在线观看入口| 亚洲国产精品久久男人天堂| 久久香蕉精品热| 亚洲av免费在线观看| 好男人在线观看高清免费视频| 啦啦啦韩国在线观看视频| 日韩人妻高清精品专区| 99国产精品一区二区蜜桃av| 精品久久久久久成人av| 麻豆久久精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| aaaaa片日本免费| 国产亚洲精品一区二区www| 国产真实伦视频高清在线观看 | 久久久久久久午夜电影| 99久久综合精品五月天人人| 国产精品久久久久久久电影 | 国产午夜精品论理片| 97人妻精品一区二区三区麻豆| 亚洲av熟女| 国产久久久一区二区三区| 国产视频一区二区在线看| 色吧在线观看| x7x7x7水蜜桃| 麻豆成人av在线观看| 国产精品免费一区二区三区在线| 亚洲成人久久爱视频| 午夜视频国产福利| 精品国产三级普通话版| 国产真实乱freesex| 日韩欧美在线乱码| 日韩欧美一区二区三区在线观看| 亚洲人与动物交配视频| 国产精品国产高清国产av| 色综合亚洲欧美另类图片| 一本一本综合久久| 人人妻人人看人人澡| 亚洲精品日韩av片在线观看 | 特级一级黄色大片| 国产精品,欧美在线| 精品一区二区三区视频在线 | 亚洲精品影视一区二区三区av| 在线播放国产精品三级| 久久欧美精品欧美久久欧美| 亚洲欧美一区二区三区黑人| 99热这里只有是精品50| 亚洲成人精品中文字幕电影| 国产精华一区二区三区| 亚洲一区二区三区不卡视频| 精品乱码久久久久久99久播| 熟女人妻精品中文字幕| 熟妇人妻久久中文字幕3abv| 日韩欧美三级三区| 性欧美人与动物交配| 国语自产精品视频在线第100页| 国产精品免费一区二区三区在线| 欧美3d第一页| 亚洲中文字幕一区二区三区有码在线看| 亚洲av日韩精品久久久久久密| 欧美中文日本在线观看视频| 色哟哟哟哟哟哟| 国产精品精品国产色婷婷| 精品99又大又爽又粗少妇毛片 | 国语自产精品视频在线第100页| 男女床上黄色一级片免费看| av片东京热男人的天堂| 亚洲 欧美 日韩 在线 免费| 观看美女的网站| 中文字幕人妻熟人妻熟丝袜美 | 婷婷丁香在线五月| 老熟妇乱子伦视频在线观看| 国产成人系列免费观看| 亚洲内射少妇av| 高潮久久久久久久久久久不卡| 国产极品精品免费视频能看的| 亚洲成人中文字幕在线播放| 极品教师在线免费播放| 99久久精品国产亚洲精品| 亚洲午夜理论影院| 日韩欧美国产一区二区入口| 国产高清视频在线观看网站| 亚洲av一区综合| 午夜精品久久久久久毛片777| 亚洲国产精品久久男人天堂| 欧美日本视频| 每晚都被弄得嗷嗷叫到高潮| 午夜a级毛片| 最近最新中文字幕大全免费视频| or卡值多少钱| 一个人观看的视频www高清免费观看| 免费在线观看亚洲国产| 两个人看的免费小视频| 又紧又爽又黄一区二区| 男女下面进入的视频免费午夜| 我要搜黄色片| 免费看日本二区| 一级毛片高清免费大全| 久久久久久九九精品二区国产| 欧美性感艳星| 99热精品在线国产| 久久这里只有精品中国| 三级毛片av免费| 成年版毛片免费区| 叶爱在线成人免费视频播放| 脱女人内裤的视频| 丰满的人妻完整版| 99久久综合精品五月天人人| 亚洲精品亚洲一区二区| 国产亚洲av嫩草精品影院| 91九色精品人成在线观看| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩高清专用| 欧美激情久久久久久爽电影| 精品久久久久久久毛片微露脸| 校园春色视频在线观看| 国产精品影院久久| 18美女黄网站色大片免费观看| 少妇的逼水好多| 搡女人真爽免费视频火全软件 | 国产伦精品一区二区三区四那| 91av网一区二区| 久久久色成人| 少妇高潮的动态图| 极品教师在线免费播放| 99热精品在线国产| 国产私拍福利视频在线观看| 老司机午夜福利在线观看视频| 日本黄色视频三级网站网址| 亚洲av熟女| 成人欧美大片| 国产精品久久久久久精品电影| x7x7x7水蜜桃| 12—13女人毛片做爰片一| 免费观看人在逋| 国产成人a区在线观看| 国产成人系列免费观看| 国产中年淑女户外野战色| 国语自产精品视频在线第100页| 亚洲av日韩精品久久久久久密| 成人18禁在线播放| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| www日本在线高清视频| 亚洲国产欧洲综合997久久,| 99久久成人亚洲精品观看| 午夜免费成人在线视频| 别揉我奶头~嗯~啊~动态视频| 成人国产综合亚洲| 噜噜噜噜噜久久久久久91| 国产aⅴ精品一区二区三区波| 午夜福利在线在线| 欧美成狂野欧美在线观看| 99国产精品一区二区三区| 国产国拍精品亚洲av在线观看 | 国产高清有码在线观看视频| 看黄色毛片网站| 久久中文看片网| 少妇的丰满在线观看| 一区二区三区高清视频在线| 十八禁人妻一区二区| 欧美+日韩+精品| 欧美黑人欧美精品刺激| 尤物成人国产欧美一区二区三区| 男女下面进入的视频免费午夜| 非洲黑人性xxxx精品又粗又长| 麻豆成人av在线观看| 欧美日韩乱码在线| 久久久久久久久中文| 精品乱码久久久久久99久播| 色尼玛亚洲综合影院| 国产野战对白在线观看| 非洲黑人性xxxx精品又粗又长| 在线天堂最新版资源| 日本a在线网址| 国产激情偷乱视频一区二区| 亚洲,欧美精品.| 特大巨黑吊av在线直播| 国产精品久久久久久亚洲av鲁大| 亚洲国产色片| 欧美日韩综合久久久久久 | 国产一区在线观看成人免费| 两人在一起打扑克的视频| 国产一级毛片七仙女欲春2| 国产精品久久久久久人妻精品电影| 天堂av国产一区二区熟女人妻| 国产免费av片在线观看野外av| 欧美丝袜亚洲另类 | 欧美不卡视频在线免费观看| 午夜亚洲福利在线播放| 搡老熟女国产l中国老女人| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 无遮挡黄片免费观看| 成人性生交大片免费视频hd| 亚洲色图av天堂| 一二三四社区在线视频社区8| 性色avwww在线观看| 国产精品久久视频播放| www国产在线视频色| 午夜福利在线观看免费完整高清在 | 最后的刺客免费高清国语| 日韩欧美国产在线观看| 久久久久久人人人人人| 香蕉丝袜av| 日本在线视频免费播放| 免费在线观看影片大全网站| 国产一区二区三区在线臀色熟女| 亚洲av不卡在线观看| 国产精品久久久久久久电影 | 丁香六月欧美| 美女大奶头视频| 午夜免费激情av| 免费一级毛片在线播放高清视频| 在线视频色国产色| 亚洲五月天丁香| 精品国内亚洲2022精品成人| 色在线成人网| 日韩高清综合在线| 看黄色毛片网站| 国内揄拍国产精品人妻在线| 亚洲五月婷婷丁香| 小蜜桃在线观看免费完整版高清| 日韩成人在线观看一区二区三区| 亚洲精品影视一区二区三区av| 免费在线观看日本一区| netflix在线观看网站| 婷婷精品国产亚洲av在线| 亚洲国产精品成人综合色| 中文字幕久久专区| 亚洲美女黄片视频| 日韩高清综合在线| 午夜日韩欧美国产| 国产高清激情床上av| 亚洲欧美日韩高清专用| 国产精品99久久久久久久久| www.色视频.com| 欧美日韩福利视频一区二区| 日本熟妇午夜| www.999成人在线观看| 亚洲成a人片在线一区二区| 欧美最新免费一区二区三区 | 国产精品久久久人人做人人爽| 亚洲人成电影免费在线| 日韩精品中文字幕看吧| 一二三四社区在线视频社区8| 免费看日本二区| 国产高潮美女av| 亚洲国产精品成人综合色| 午夜免费男女啪啪视频观看 | 国产精品嫩草影院av在线观看 | АⅤ资源中文在线天堂| tocl精华| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 亚洲aⅴ乱码一区二区在线播放| 国产av一区在线观看免费| 亚洲国产日韩欧美精品在线观看 | 国产黄色小视频在线观看| 精品熟女少妇八av免费久了| 国产精品野战在线观看| 色老头精品视频在线观看| 丰满的人妻完整版| 在线观看日韩欧美| 日韩国内少妇激情av| 99热只有精品国产| 欧美丝袜亚洲另类 | 又爽又黄无遮挡网站| 国产精品亚洲一级av第二区| 九九热线精品视视频播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人久久性| 嫩草影院入口| 免费av观看视频| 香蕉av资源在线| 国产精品永久免费网站| 国产精品,欧美在线| 国产aⅴ精品一区二区三区波| 久久久久久久久大av| 级片在线观看| 亚洲,欧美精品.| 最近最新中文字幕大全电影3| 国产私拍福利视频在线观看| 亚洲在线观看片| 中文字幕人妻丝袜一区二区| 日韩大尺度精品在线看网址| 国产成人福利小说| 中文字幕av在线有码专区| 亚洲无线在线观看| 国产伦一二天堂av在线观看| 亚洲av熟女| 老熟妇乱子伦视频在线观看| 欧美成人a在线观看| 一个人看的www免费观看视频| 男女视频在线观看网站免费| 久久香蕉精品热| 禁无遮挡网站| 日韩人妻高清精品专区| 97超级碰碰碰精品色视频在线观看| 国产精品影院久久| 国产精品久久久久久久久免 | 亚洲18禁久久av| 国产成人a区在线观看| www.999成人在线观看| 日韩大尺度精品在线看网址| 日韩欧美一区二区三区在线观看| 狂野欧美激情性xxxx| 夜夜看夜夜爽夜夜摸| 国内毛片毛片毛片毛片毛片| 男女那种视频在线观看| 精华霜和精华液先用哪个| 熟女人妻精品中文字幕| 国产精品1区2区在线观看.| 国产欧美日韩精品亚洲av| 国产97色在线日韩免费| 我的老师免费观看完整版| 国产精品久久久久久人妻精品电影| 黄色丝袜av网址大全| 国产精品影院久久| 午夜免费激情av| 久久精品国产亚洲av涩爱 | 久久欧美精品欧美久久欧美| 男插女下体视频免费在线播放| 在线播放无遮挡| 男女床上黄色一级片免费看| 又爽又黄无遮挡网站| 一区二区三区高清视频在线| 国产精品久久视频播放| 香蕉久久夜色| 亚洲国产欧洲综合997久久,| 法律面前人人平等表现在哪些方面| 人人妻人人澡欧美一区二区| 亚洲成av人片免费观看| 99久久九九国产精品国产免费| or卡值多少钱| 午夜精品一区二区三区免费看| 色噜噜av男人的天堂激情| 在线看三级毛片| 久久久久久人人人人人| 亚洲成a人片在线一区二区| 97人妻精品一区二区三区麻豆| 精品久久久久久成人av| 在线免费观看的www视频| 精品国内亚洲2022精品成人| av在线天堂中文字幕| 99国产精品一区二区三区| 夜夜爽天天搞| 国产一区二区三区视频了| 亚洲成人久久爱视频| 人人妻人人澡欧美一区二区| 亚洲久久久久久中文字幕| 色哟哟哟哟哟哟| 麻豆一二三区av精品| 丰满的人妻完整版| 99久国产av精品| 欧美+亚洲+日韩+国产| 桃红色精品国产亚洲av| 特级一级黄色大片| 亚洲av美国av| 亚洲人成网站在线播| 色在线成人网| 一进一出抽搐动态| 亚洲国产高清在线一区二区三| 国产精品久久久久久人妻精品电影| 好看av亚洲va欧美ⅴa在| 成年女人永久免费观看视频| 高清日韩中文字幕在线| 一个人免费在线观看的高清视频| 精品午夜福利视频在线观看一区| 亚洲 国产 在线| 免费看十八禁软件| h日本视频在线播放| 亚洲国产精品成人综合色| 人妻丰满熟妇av一区二区三区| 亚洲精品色激情综合| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 午夜激情福利司机影院| 国产爱豆传媒在线观看| 99精品久久久久人妻精品| 亚洲成人久久性| 青草久久国产| 在线a可以看的网站| 精品免费久久久久久久清纯| 看黄色毛片网站| 国产成人福利小说| 欧美国产日韩亚洲一区| 久久精品国产亚洲av香蕉五月| 色尼玛亚洲综合影院| 亚洲激情在线av| 午夜久久久久精精品| 久久中文看片网| 久久人妻av系列| 亚洲激情在线av| 国产成人av激情在线播放| 岛国视频午夜一区免费看| 一本精品99久久精品77| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 婷婷精品国产亚洲av| а√天堂www在线а√下载| 91麻豆av在线| 老熟妇仑乱视频hdxx| 舔av片在线| 欧美一级a爱片免费观看看| 亚洲av免费在线观看| 男女那种视频在线观看| 国产欧美日韩一区二区三| xxx96com| 波多野结衣高清作品| 国产高清videossex| 高清日韩中文字幕在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久久久久久久久| 国产精品一及| 国产高清videossex| 一进一出好大好爽视频| 欧美丝袜亚洲另类 | av黄色大香蕉| 一进一出好大好爽视频| 国产aⅴ精品一区二区三区波| 久久精品国产清高在天天线| 在线播放国产精品三级| 丰满的人妻完整版| 中文字幕人妻熟人妻熟丝袜美 | 真实男女啪啪啪动态图| 久久中文看片网| 亚洲天堂国产精品一区在线| 国产高清激情床上av| 校园春色视频在线观看| av欧美777| 校园春色视频在线观看| 精品国产三级普通话版| 中文字幕人妻熟人妻熟丝袜美 | 很黄的视频免费| 一区二区三区高清视频在线| tocl精华| 色哟哟哟哟哟哟| 久久午夜亚洲精品久久| av天堂中文字幕网| 真人一进一出gif抽搐免费| 在线十欧美十亚洲十日本专区| 日日干狠狠操夜夜爽| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| www日本在线高清视频| 一个人观看的视频www高清免费观看| 欧美一级a爱片免费观看看| 欧美日韩国产亚洲二区| 亚洲18禁久久av| 精品久久久久久成人av| 久久精品夜夜夜夜夜久久蜜豆| 老司机午夜福利在线观看视频| 亚洲成人中文字幕在线播放| 亚洲久久久久久中文字幕| 精品国产美女av久久久久小说| 亚洲人成伊人成综合网2020| 欧美+亚洲+日韩+国产| 很黄的视频免费| 99国产极品粉嫩在线观看| www国产在线视频色| 日本 av在线| 久久精品国产99精品国产亚洲性色| 亚洲精品成人久久久久久| 天天躁日日操中文字幕| 麻豆久久精品国产亚洲av| 在线看三级毛片| 在线观看一区二区三区| 婷婷丁香在线五月| 黄片小视频在线播放| 国产欧美日韩精品一区二区| 亚洲av一区综合| 欧美区成人在线视频| 在线观看日韩欧美| 欧美+亚洲+日韩+国产| 欧美黄色淫秽网站| 成人三级黄色视频| 国产成人欧美在线观看| 一进一出好大好爽视频| 小蜜桃在线观看免费完整版高清| 成年女人看的毛片在线观看| tocl精华| 日韩欧美精品v在线| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费一区二区三区在线| 久久久精品大字幕| 极品教师在线免费播放| 伊人久久精品亚洲午夜| 桃红色精品国产亚洲av| 嫩草影视91久久| 婷婷丁香在线五月| 亚洲av免费在线观看| a级一级毛片免费在线观看| 琪琪午夜伦伦电影理论片6080| 日本熟妇午夜| 给我免费播放毛片高清在线观看| 午夜两性在线视频| x7x7x7水蜜桃| www.色视频.com| 久久久久久久久中文| 欧美黑人巨大hd| 老汉色∧v一级毛片| 狂野欧美白嫩少妇大欣赏| 一a级毛片在线观看| 国产久久久一区二区三区| 国产精品三级大全| 成人高潮视频无遮挡免费网站| 黄片大片在线免费观看| 日韩精品青青久久久久久| 国产97色在线日韩免费| 日韩欧美一区二区三区在线观看| 中文在线观看免费www的网站| 九九热线精品视视频播放| 黄色日韩在线| 欧美成人性av电影在线观看| 一区二区三区激情视频| 久久久久性生活片| 亚洲激情在线av| 成人精品一区二区免费| а√天堂www在线а√下载| 国产一级毛片七仙女欲春2| 中文字幕精品亚洲无线码一区| 午夜a级毛片| 免费在线观看日本一区| 久久久久久久午夜电影| 久久久久免费精品人妻一区二区| 深爱激情五月婷婷| 网址你懂的国产日韩在线| 亚洲欧美日韩高清专用| 亚洲一区二区三区不卡视频|