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

    拱頂罐內(nèi)油品蒸發(fā)損耗的數(shù)值模擬和實驗研究

    2020-08-14 07:13:18張佩宇黃維秋景海波孫憲航
    安全與環(huán)境工程 2020年4期
    關鍵詞:損耗量罐頂罐壁

    張佩宇,黃維秋,景海波,2,汪 城,孫憲航

    (1.常州大學石油工程學院,江蘇 常州 213164;2.合肥燃氣集團有限公司,安徽 合肥 230075)

    儲罐區(qū)揮發(fā)性有機物(VOCs)的排放會引起以光化學煙霧、城市霧霾為特征的復合型大氣污染問題。現(xiàn)階段有效控制儲罐區(qū)VOCs的排放及回收工作是我國大氣環(huán)境治理領域中的熱點問題[1-3]。為了降低儲罐區(qū)油品的蒸發(fā)損耗,國內(nèi)相關規(guī)范[4]規(guī)定儲存揮發(fā)性有機液體的儲罐應選用浮頂罐和集成油氣回收處理系統(tǒng)以及密閉排氣系統(tǒng)的常壓拱頂罐。拱頂罐與油氣回收集成工藝在油氣排放量的最終控制、整體效益上優(yōu)于單純的浮頂罐儲存工藝,是一種較為理想的輕質(zhì)油品儲存工藝[5]。通過對拱頂罐內(nèi)油氣和空氣的傳熱和傳質(zhì)機理進行分析并測算油罐排氣量,可為拱頂罐及其油氣回收集成系統(tǒng)的設計和管理以及降低拱頂罐小呼吸損耗提供理論依據(jù),以指導管理部門科學地制訂相關防治方案。

    國內(nèi)外的研究者對油品蒸發(fā)模型進行了一些探索。目前,數(shù)值模擬已成為油品蒸發(fā)損耗的主要研究方法[6-12]。如Wang等[13]和Li等[14]研究了大型浮頂罐在不同條件下儲罐油品溫度和溫降曲線的變化規(guī)律;Hata等[15]通過對汽車油箱進行呼吸損失測試,建立了一種汽車蒸發(fā)損耗量的估算模型;Saufi等[16]和Christian等[17]通過實驗分析對液滴蒸發(fā)過程進行了研究。

    上述研究均未從太陽輻射的角度對拱頂罐內(nèi)、罐壁溫度對儲罐油品蒸發(fā)損耗的影響進行研究,為了有效地降低拱頂罐內(nèi)油品的蒸發(fā)損耗,有必要對其蒸發(fā)損耗機理和評估方法展開研究。為此,本文建立了拱頂罐的非穩(wěn)態(tài)傳熱傳質(zhì)理論模型,自編Matlab程序并通過自行搭建實驗平臺驗證其可行性,對某液位為2 875 mm的1 000 m3拱頂罐在春分日6∶00~18∶00時間段儲罐內(nèi)油品的蒸發(fā)損耗過程進行了數(shù)值模擬,分析了儲罐內(nèi)溫度、傳熱系數(shù)、液相蒸發(fā)量的變化規(guī)律,并估計了儲罐的小呼吸損耗量,可為拱頂罐油品蒸發(fā)損耗的評估及其油氣收集回收系統(tǒng)的設計、管理提供依據(jù)。

    1 拱頂罐的非穩(wěn)態(tài)傳熱傳質(zhì)理論模型的建立

    影響儲罐內(nèi)油氣蒸發(fā)的因素很多,本文對拱頂罐非穩(wěn)態(tài)傳熱傳質(zhì)模型進行了必要的簡化處理,即忽略罐內(nèi)氣相、液相、罐壁和罐頂?shù)膶釤嶙韬蜌庀嗟膫髻|(zhì)阻力,分析了儲罐內(nèi)外的傳熱傳質(zhì)規(guī)律[18]。建立的拱頂罐非穩(wěn)態(tài)傳熱傳質(zhì)模型如圖1所示,外界向罐內(nèi)液相、氣相傳遞熱量,同時氣相向液相傳遞熱量;液相吸收的熱量全部用于自身溫升和蒸發(fā),液相蒸發(fā)即液相向氣相傳遞質(zhì)量,氣相吸收的熱量全部用于氣相溫升。

    圖1 拱頂罐的非穩(wěn)態(tài)傳熱傳質(zhì)理論模型Fig.1 Theoretical model of unsteady heat and mass transfer in the dome-roof tank

    1.1 外界與氣相和外界與液相間的傳熱

    外界向罐內(nèi)氣相傳遞的熱流量(QG)包括經(jīng)罐壁和罐頂傳入氣相的熱流量。外界向罐內(nèi)液相傳遞的熱流量(QL)包括經(jīng)罐壁和罐頂傳入液相的熱流量。由于經(jīng)罐底傳入罐內(nèi)液相的熱流量很小,所以可忽略不計。

    經(jīng)氣體空間罐壁傳入氣相的熱流量可按下式計算[18]:

    (1)

    (2)

    經(jīng)罐頂傳入氣相的熱流量可按下式計算:

    (3)

    (4)

    經(jīng)液體空間罐壁傳入液相的熱流量即外界向液相傳遞的熱流量(QL)可按下式計算:

    (5)

    (6)

    1.2 罐頂和罐壁吸收的太陽輻射量

    與太陽光照射方向垂直的單位面積罐頂或罐壁吸收的太陽輻射熱流密度(q0)可按下式計算:

    (7)

    式中:IC為太陽常數(shù),其值為1 372.57 W/m2;P為大氣透明系數(shù),其值為0.7~0.8;σ為與晝長有關的系數(shù);m1為與大氣質(zhì)量有關的系數(shù),m1=2/cosθ;θ為太陽正午時的天頂角(°);α為罐頂和罐壁的黑度系數(shù)。

    由蘭貝特定律可知,任意時刻罐壁和罐頂吸收的太陽輻射熱流密度可按下面公式計算:

    (8)

    (9)

    1.3 大氣環(huán)境溫度

    一天之內(nèi)大氣溫度變化用余弦函數(shù)表示。將起始時刻設定為日出時刻,任一時刻的大氣溫度可按下式計算:

    (10)

    1.4 氣液兩相間的傳熱

    氣液兩相間傳遞的熱流量(Qin)可按下式計算:

    Qin=Ainain(TL-TG)

    (11)

    式中:Qin為氣液兩相間傳遞的熱流量(W);ain為氣液兩相間的傳熱系數(shù)[W/(m2·K)];Ain為氣液兩相界面的面積(m2)。

    1.5 氣相溫度和液相蒸發(fā)量

    由氣相的熱平衡可計算得出氣相溫度(TG)隨時間的變化規(guī)律,其計算公式為:

    (12)

    式中:ρG為氣相密度(kg/m3);CPG為氣相比熱容[J/(kg·K)];VG為氣相體積(m3);τ、τ0分別為任意時刻和初始時刻(s);TG0為氣相初始時刻溫度(K)。

    由液相的熱平衡可計算得出液相蒸發(fā)量(m)隨時間的變化規(guī)律,其計算公式為:

    (13)

    式中:ρL為液相密度(kg/m3);CPG為液相比熱容[J/(kg·K)];VL為氣相體積(m3);h為液相汽化潛熱(J/kg);m0為初始時刻蒸發(fā)量(kg)。

    1.6 氣體膨脹量

    常壓罐呼出的混合氣體體積(Vh)由兩部分組成:第一部分為由于氣相溫升而引起的氣體體積膨脹;第二部分為逸入氣體空間的油氣所引起的氣體體積膨脹。氣體空間溫升而引起的氣體體積膨脹累積量Vh1可按下式計算:

    (14)

    式中:T1、T2分別為相鄰時刻的溫度(K);V為罐內(nèi)氣體空間體積(m3)。

    逸入氣體空間的油氣所引起的氣體體積膨脹累積量(Vh2)可按下式計算:

    (15)

    式中:m為某一時間段內(nèi)的液相蒸發(fā)質(zhì)量(kg);M為油氣摩爾質(zhì)量(kg/mol),其值為65 kg/mol;R為通用氣體常數(shù)[J/(mol·K)],其值為8.314 J/(mol·K);P為罐內(nèi)氣體壓強(Pa),其值為101 325 Pa。

    2 數(shù)值模擬和實驗驗證

    拱頂罐內(nèi)氣相溫度和液相蒸發(fā)量的數(shù)值模擬主要是基于變物性模型,通過建立拱頂罐的非穩(wěn)態(tài)傳熱傳質(zhì)理論模型,并對模型進行數(shù)值求解,具體數(shù)值模擬流程見圖2。

    圖2 拱頂罐內(nèi)氣相溫度和液相蒸發(fā)量的數(shù)值模擬流程Fig.2 Flow chart of numerical simulation of gas phase temperature and liquid evaporation in the dome-roof tank

    為了驗證建立的拱頂罐非穩(wěn)態(tài)傳熱傳質(zhì)理論模型的準確性,利用自制的儲罐實測罐內(nèi)氣相溫度變化并與數(shù)值模擬結(jié)果進行對比。自制儲罐的規(guī)格如下:罐徑D=500 mm,罐壁高度H1=500 mm,罐頂高度H2=100 mm,液位為100 mm,儲罐罐頂面積和氣液界面面積為0.196 3 m2,與氣體空間接觸的罐壁面積為0.628 m2,與液體空間接觸的罐壁面積為0.157 m2,氣相體積為0.085 m3。儲罐內(nèi)液體為92#汽油,密度為750 kg/m3,汽化潛熱為3.8×105J/kg,比熱容為1 821 J/(kg·K),油品氣相初始溫度、液相初始溫度和初始蒸發(fā)量分別為294 K、297.2 K和0 kg。

    圖3為儲罐的小呼吸損耗實驗原理圖。首先將熱電偶安裝在儲罐的罐內(nèi)、罐壁和罐頂測溫點上,再將熱電偶與溫度巡檢儀連接來讀取溫度值。當罐內(nèi)氣體空間溫度升高時,氣體膨脹導致罐內(nèi)氣體外溢;當罐內(nèi)氣體空間溫度降低時,氣體收縮導致外部空氣進入罐內(nèi),故可以通過氣體流量計來計量氣體流量。然后在罐頂開孔,利用直徑2 mm的聚乙烯管通過罐頂開孔伸入罐內(nèi)氣體取樣點進行油氣取樣,并利用氣相色譜來測定油氣樣本濃度??紤]到環(huán)境風速的影響,可以使用風速儀來測量風速,即在實驗開始后,每隔1個小時記錄1次溫度巡檢儀、氣體流量計和風速儀的讀數(shù),每隔2個小時對油氣采樣并測量其濃度。最后采用體積-濃度法按下式計算儲罐小呼吸損耗量:

    圖3 儲罐的小呼吸損耗實驗原理圖Fig.3 Experimental schematic diagram of breathing loss in the dome-roof tank

    (16)

    式中:ΔM為儲罐的小呼吸損耗量(kg);Vh,i為某時刻儲罐呼出的混合氣體積(m3);Cy,i為某時刻儲罐呼出的混合氣質(zhì)量濃度(g/m3)。

    在環(huán)境溫度和風速影響下,儲罐內(nèi)氣體溫度、氣體膨脹量和小呼吸損耗量實驗值與模擬值的對比,見圖4、圖5和表1。

    圖5 儲罐內(nèi)氣體膨脹量實驗值與模擬值的對比Fig.5 Comparison of volume expansion between experimental values and simulated results in the storage tank

    圖4 儲罐內(nèi)氣體溫度實驗值與模擬值的對比Fig.4 Comparison of gas temperature between experimental values and simulated values in the storage tank

    表1 儲罐的小呼吸損耗量實驗值與模擬值的對比Table 1 Comparison of breathing losses between experimental values and simulated values in the storage tank

    由圖4、圖5和表1可見,儲罐內(nèi)氣體溫度、氣體膨脹量和小呼吸損耗量實驗值與模擬值吻合得較好,從而驗證了本文所建立的拱頂罐非穩(wěn)態(tài)傳熱傳質(zhì)理論模型和計算方法的準確性。

    3 實例應用與分析

    本文利用上述模擬方法,以江蘇省常州地區(qū)某液位為2 875 mm的1 000 m3拱頂汽油罐為研究對象,對該拱頂罐在春分日6∶00~18∶00時間段儲罐內(nèi)油品的蒸發(fā)損耗過程進行了數(shù)值模擬,分析了儲罐內(nèi)氣相、罐頂、氣體空間罐壁和液體空間罐壁溫度、傳熱系數(shù)和液相蒸發(fā)量,并估算了儲罐的小呼吸損耗量。該拱頂儲罐數(shù)、汽油物性參數(shù)和外界環(huán)境條件見表2。其中,儲罐內(nèi)平均油品溫度一天之內(nèi)變化幅度不大,一般僅為1~3℃,因此假定罐內(nèi)油品溫度在春分日6∶00~18∶00時間段內(nèi)上升2℃[19]。

    表2 某拱頂儲罐、汽油物性參數(shù)和外界環(huán)境條件Table 2 Basic parameters of a dome-roof tank,gasoline physical properties and external environmental conditions

    3.1 儲罐溫度的變化規(guī)律

    儲罐溫度主要受太陽輻射和大氣環(huán)境溫度的影響而發(fā)生變化。在早上6∶00時大氣環(huán)境溫度近似等于晝夜間最低大氣溫度,最高大氣溫度一般出現(xiàn)在中午12∶00~14∶00。同時,在中午時,由于太陽高度角最大,太陽輻射經(jīng)過的大氣行程最短,受到的削弱也最小,到達地面的太陽輻射量也最多,儲罐吸收的太陽輻射量也最多,所以儲罐溫度也最高。圖6為儲罐內(nèi)氣相、罐頂、罐壁溫度和大氣環(huán)境溫度隨時間的變化趨勢圖。

    圖6 儲罐內(nèi)氣相、罐頂、罐壁溫度和大氣環(huán)境溫度隨 時間的變化趨勢圖Fig.6 Change curves of the temperature of gas phase, tank top,tank wall and the atmospheric environment with time 注:TG、Td、TGbi、TLbi、Tair分別為儲罐內(nèi)氣相、罐頂、氣體空間罐壁、液體空間罐壁和大氣環(huán)境的溫度。

    由圖6可見,儲罐內(nèi)氣相、罐頂、氣體空間罐壁和液體空間罐壁溫度的變化趨勢與大氣環(huán)境溫度的變化趨勢一致,呈先升高后降低,最高溫度出現(xiàn)在14∶00;由于氣相比熱容較小,其升溫最快,單位熱流量下的溫升最大;由于液相溫度較低且變化不大,故液體空間罐壁溫度小于氣體空間罐壁和罐頂?shù)臏囟惹易兓淮螅粴庀鄿囟染佑趦捱吔?罐壁和罐頂)的溫度和大氣溫度之間。

    3.2 儲罐傳熱系數(shù)的變化規(guī)律

    圖7為儲罐內(nèi)氣相、罐頂、罐壁自然對流換熱系數(shù)隨時間的變化趨勢圖。

    圖7 儲罐內(nèi)氣相、罐頂、罐壁自然對流換熱系數(shù) 隨時間的變化趨勢圖Fig.7 Change curves of the natural convection heat transfer coefficient of gas phase,tank top, tank wall with time 注:aG、aL、ad分別為儲罐內(nèi)氣相與罐壁間、液相與罐壁間、氣相與罐頂間的自然對流換熱系數(shù)。

    由圖7可見,儲罐內(nèi)氣相與罐壁間和氣相與罐頂間的自然對流換熱系數(shù)的變化波動很小,液相與罐壁間的自然對流換熱系數(shù)的變化波動較大。這是由于自然對流換熱系數(shù)與介質(zhì)的密度、導熱系數(shù)、運動黏度、比熱容和溫差有關,很顯然液相與氣相相比更易傳熱,因此液相與罐壁間的自然對流換熱系數(shù)大于氣相與罐壁、罐頂間的自然對流換熱系數(shù);又由于罐壁縱向自然對流強度大于罐頂橫向自然對流強度,所以氣相與罐壁間的自然對流換熱系數(shù)大于氣相與罐頂間的自然對流換熱系數(shù)。根據(jù)模擬結(jié)果可知,儲罐內(nèi)氣相與罐頂間自然對流換熱系數(shù)的變化范圍為1.40~1.61 W/(m2·K),氣相與罐壁間自然對流換熱系數(shù)的變化范圍為1.64~2.10 W/(m2·K),液相與罐壁間自然對流換熱系數(shù)的變化范圍為31.05~73.05 W/(m2·K)。

    圖8為儲罐內(nèi)氣相、罐頂、罐壁輻射換熱系數(shù)隨時間的變化趨勢圖。

    圖8 儲罐內(nèi)氣相、罐頂、罐壁輻射換熱系數(shù)隨時間的 變化趨勢圖Fig.8 Change curves of the radiation heat transfer coefficient of gas phase,tank top,tank wall with time 注:a2L、a2d、a2G分別為液體空間罐壁、罐頂、氣體空間罐壁的輻射換熱系數(shù)。

    由圖8可見,儲罐內(nèi)罐頂、氣體空間罐壁和液體空間罐壁的輻射換熱系數(shù)變化趨勢一致,變化幅度也都不大。這是由于儲罐邊界黑度恒定,輻射換熱系數(shù)的大小主要取決于儲罐邊界與大氣環(huán)境之間的溫差,由圖6可明顯看出,氣體空間罐壁與大氣環(huán)境之間的溫差最大,罐頂次之,液體空間罐壁最小,因此氣體空間罐壁的輻射換熱系數(shù)最大,罐頂次之,液體空間罐壁最小;同時中午時分儲罐邊界與大氣環(huán)境之間的溫差較大,所以中午時儲罐邊界的輻射換熱系數(shù)最大。根據(jù)模擬結(jié)果可知,儲罐內(nèi)罐頂?shù)妮椛鋼Q熱系數(shù)變化范圍為4.46~4.80 W/(m2·K),氣體空間罐壁的輻射換熱系數(shù)變化范圍為4.47~4.81 W/(m2·K),液體空間罐壁的輻射換熱系數(shù)變化范圍為4.06~4.18 W/(m2·K)。

    儲罐邊界的總傳熱系數(shù)由自然對流換熱系數(shù)、輻射換熱系數(shù)、強制對流換熱系數(shù)、儲罐邊界的厚度和導熱率決定。假設模擬的環(huán)境風速為定值(4 m/s),則儲罐內(nèi)罐頂、罐壁的強制對流換熱系數(shù)也為定值。儲罐內(nèi)罐壁的強制對流換熱系數(shù)近似于空氣橫掠單管時的換熱系數(shù),罐頂?shù)膹娭茖α鲹Q熱系數(shù)近似于空氣掠過平板時的換熱系數(shù),分別為2.30 W/(m2·K)和7.70 W/(m2·K)。儲罐邊界的厚度很小而導熱率又非常大,所以可忽略罐壁和罐頂?shù)臒嶙?。因此,根?jù)模擬可得到儲罐邊界的總傳熱系數(shù),具體見表3。

    表3 儲罐邊界的總傳熱系數(shù)Table 3 Total heat transfer coefficient of the storage tank boundary

    由表3可知,儲罐內(nèi)液體空間罐壁的總傳熱系數(shù)最大,其變化范圍為8.57~10.18 W/(m2·K),氣體空間罐壁的總傳熱系數(shù)次之,其變化范圍為1.45~1.80 W/(m2·K),罐頂?shù)目倐鳠嵯禂?shù)最小,變化范圍為1.16~1.31 W/(m2·K)。

    3.3 罐內(nèi)液相蒸發(fā)量和氣體膨脹量

    儲罐內(nèi)液相蒸發(fā)量主要由液相的氣化潛熱、外界與液相以及氣液兩相間的傳熱量決定。儲罐內(nèi)液相蒸發(fā)量的模擬結(jié)果見圖9。

    圖9 儲罐內(nèi)液相蒸發(fā)累積量與時間的變化趨勢圖Fig.9 Change curves of the liquid phase evaporation accumulation of gas phase with time

    由圖9可見,在春分日6∶00~18∶00時間段內(nèi)儲罐內(nèi)液相蒸發(fā)累積量逐漸增高,但其增加速率在早晨和下午時較大,中午時較小,這主要是由于早晨和下午時,太陽輻射主要集中于罐壁,液相吸收的太陽輻射熱流量較多,有利于液相的蒸發(fā),其蒸發(fā)量較大,而中午時,太陽輻射主要集中于罐頂,太陽輻射熱流量自罐頂向下傳遞,由于氣相熱阻較大,不利于傳熱,故液相的蒸發(fā)量較??;在初始時刻蒸發(fā)量為0的條件下,在春分日6∶00~18∶00時間段內(nèi),液位為2 875 mm的1 000 m3拱頂罐內(nèi)液相蒸發(fā)累積量為421.13 kg。已知罐內(nèi)氣體空間體積為850 m3,假設沒有油氣逸入外界環(huán)境,則罐內(nèi)平均油氣濃度為495.45 g/m3。另外,夜晚時,由于大氣環(huán)境溫度降低,且無太陽輻射,儲罐開始向外界大氣散熱,罐內(nèi)油氣冷凝,液相蒸發(fā)量會減少。

    根據(jù)式(14)和式(15)可計算得到氣相溫升引起的氣體累積體積膨脹量Vh1為56.62 m3,逸入氣體空間的油氣所引起的氣體體積膨脹累積量Vh2為157.44 m3,氣體體積總膨脹累積量Vh為214.06 m3。圖10為儲罐內(nèi)氣體體積膨脹量隨時間的變化趨勢圖。

    圖10 儲罐內(nèi)氣體體積膨脹累積量隨時間的變化 趨勢圖Fig.10 Change curves of the sum of gas volume expansion accumulation with time 注:Vh1、Vh2、Vh分別為氣相溫升引起的氣體體積膨脹量、油氣逸入氣體空間引起的氣體體積膨脹量、油罐呼出的混合氣體積。

    由圖10可見,Vh1主要受氣相溫升的影響,Vh1與氣相溫度變化趨勢一致。在早晨6∶00至7∶00時間段,由于氣相溫升最大,此時的Vh1增加速率也最大,隨后氣相溫度緩慢上升,Vh1也隨之緩慢增加,在中午14∶00時氣相溫度達到最高值,Vh1也達到最高值,14∶00之后,氣相溫度逐漸降低,Vh1也隨之逐漸降低。Vh2主要受液相蒸發(fā)量的影響,氣體體積膨脹累積量與液相蒸發(fā)累積量的變化趨勢一致。在早晨和下午時,液相蒸發(fā)速率較快,氣體體積膨脹累積量的增加速率也較快,在中午時,液相蒸發(fā)速率較慢,氣體體積膨脹累積量增加速率也較慢。

    3.4 儲罐小呼吸損耗量分析

    由于夜晚期間儲罐不再吸收太陽輻射,儲罐向外界傳熱,因此可認為夜晚儲罐的小呼吸損耗量為0,白天儲罐的小呼吸損耗量即為儲罐全天的小呼吸損耗量,為421.13 kg。另一方面,本模型很難得到罐內(nèi)呼出的混合氣體濃度,且其濃度必小于罐內(nèi)平均濃度,罐內(nèi)氣體空間體積為850 m3,假設沒有油氣逸入外界環(huán)境,則罐內(nèi)平均油氣濃度由c=Δm/V計算得到為495.45 g/m3;為此,可假設呼出的混合氣體油氣濃度分別為50 g/m3、100 g/m3、200 g/m3、250 g/m3、300 g/m3、400 g/m3,進而計算出儲罐相應的小呼吸損耗量,并與由公式(17)(美國API公式)計算得到的儲罐小呼吸損耗量進行比較[20]。美國API公式如下:

    (17)

    式中:ΔmAPI為固定頂罐一天的小呼吸損耗量(kg);PA為油品本體溫度下的蒸汽壓(kPa),取32 kPa;Pa為大氣壓(kPa),取100 kPa;ρL為液相密度(kg/m3);D為油罐直徑(m);H為氣體空間高度(m);ΔT為大氣環(huán)境溫度的平均日溫差(K);FP為涂漆系數(shù),取1;C為小罐修正系數(shù),取1;K1為單位換算系數(shù),取1;K2油品系數(shù),取1。

    基于API公式和體積-濃度法計算得到的儲罐小呼吸損耗量見表4。

    表4 基于API公式與體積-濃度法計算得到的儲罐小呼吸損耗量對比Table 4 Comparison of calculated breathing losses of the storage tank between API equation and volume-concentration method

    由表4可知,API公式計算得到的儲罐小呼吸損耗量為50.84 kg,通過對比發(fā)現(xiàn)假設該模型罐內(nèi)呼出混合氣體油氣濃度為250 g/m3時,基于體積-濃度法計算得到的儲罐小呼吸損耗量為53.50 kg,此時兩者是最為接近的。因此可以推斷出該模型的呼出混合氣體油氣濃度大約為250 g/m3時,儲罐小呼吸損耗量大約為50 kg,日損耗率約為0.23‰,降低儲罐的小呼吸損耗是必要的,可選用儲罐涂料等方法降低罐內(nèi)溫差,同時為響應國家“可持續(xù)性發(fā)展”的政策,各罐區(qū)應盡量進行油氣回收。

    4 結(jié) 論

    本文基于變物性模型,建立了拱頂罐的非穩(wěn)態(tài)傳熱傳質(zhì)理論模型,自編Matlab程序并通過自行搭建的實驗平臺驗證其可行性,并以江蘇省常州市地區(qū)某液位為2 875 mm的1 000 m3拱頂汽油罐為研究對象,研究在春分日6∶00~18∶00時間段內(nèi),儲罐的溫度、傳熱系數(shù)、液相蒸發(fā)量隨時間變化的規(guī)律,估算了儲罐的小呼吸損耗量,得到如下結(jié)論:

    (1) 儲罐內(nèi)氣相、罐頂和氣體空間罐壁溫度的變化趨勢與大氣環(huán)境溫度的變化趨勢一致,最大值均出現(xiàn)于14∶00左右;液體空間罐壁溫度最低且變化幅度很??;氣體空間罐壁和罐頂?shù)臏囟确浅=咏?,前者溫度高于后者溫度;氣相溫度居于儲罐邊?氣體空間罐壁和罐頂)溫度和大氣環(huán)境溫度之間。

    (2) 儲罐內(nèi)氣相與罐壁、罐頂間的自然對流換熱系數(shù)波動很小,而液相與罐壁間的自然對流換熱系數(shù)波動較大;罐壁與液相間的自然對流換熱系數(shù)最大,罐壁與氣相間次之,罐頂與氣相間的最小。輻射換熱系數(shù)中儲罐內(nèi)罐頂、氣體空間罐壁和液體空間罐壁的輻射換熱系數(shù)變化趨勢一致,變化幅度和數(shù)值大小相差不大;儲罐內(nèi)氣體空間罐壁的輻射換熱系數(shù)最大,罐頂次之,液體空間罐壁最小。總傳熱系數(shù)中與液相接觸的罐壁的總傳熱系數(shù)最大,與氣相接觸的罐壁的次之,罐頂?shù)淖钚 ?/p>

    (3) 儲罐內(nèi)油氣濃度為0的初始條件下,液相蒸發(fā)量為421.13 kg;氣相溫升引起的氣體體積膨脹累積量為56.62 m3,逸入氣體空間的油氣所引起的氣體體積膨脹累積量為157.44 m3,總氣體體積膨脹量為214.06 m3。假設呼出混合氣體油氣濃度分別為50 g/m3、100 g/m3、200 g/m3、250 g/m3、300 g/m3、400 g/m3,則通過體積、濃度法計算得到的儲罐小呼吸損耗量分別為10.70 kg、21.41 kg、42.81 kg、53.50 kg、64.23 kg、85.60 kg。將其與美國API公式進行比較,得出該模型呼出混合氣體油氣濃度大約為250 g/m3,小呼吸損耗量大約為50 kg。

    猜你喜歡
    損耗量罐頂罐壁
    噴射混凝土(濕噴工藝)在山嶺隧道損耗量控制
    科技資訊(2023年9期)2023-06-09 18:27:52
    大型LNG儲罐罐壁隔熱層保冷性能及其優(yōu)化對策探究
    基于AP-42的內(nèi)浮頂儲罐VOCs排放特性定量研究
    在役常壓儲罐罐壁完整性檢測與評價標準問題探討
    常壓儲罐罐頂結(jié)構力學性能分析
    GB50341和API650的罐壁設計比較
    大型LNG儲罐罐頂氣壓頂升作業(yè)動力系統(tǒng)的設計與控制
    TANKS軟件在儲罐大小呼吸計算中的應用
    下期要目
    化工設計(2017年1期)2017-03-03 06:39:56
    10萬m3LNG儲罐罐頂氣升技術措施及控制淺析
    上海煤氣(2016年3期)2016-11-09 06:06:44
    叶爱在线成人免费视频播放| 69av精品久久久久久 | 丝袜在线中文字幕| 老熟妇仑乱视频hdxx| 成人国产一区最新在线观看| 乱人伦中国视频| 日本五十路高清| 国产亚洲av片在线观看秒播厂| 国产淫语在线视频| 亚洲成人国产一区在线观看| 人人妻人人爽人人添夜夜欢视频| 19禁男女啪啪无遮挡网站| 五月开心婷婷网| 俄罗斯特黄特色一大片| 99国产精品一区二区蜜桃av | 国产亚洲精品第一综合不卡| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av香蕉五月 | 亚洲欧美一区二区三区久久| 色婷婷久久久亚洲欧美| 777久久人妻少妇嫩草av网站| 亚洲三区欧美一区| 男人舔女人的私密视频| 午夜福利视频在线观看免费| 午夜免费观看性视频| 深夜精品福利| 一个人免费在线观看的高清视频 | 久久久久国产一级毛片高清牌| 多毛熟女@视频| 亚洲欧美日韩高清在线视频 | 嫩草影视91久久| 欧美97在线视频| 在线精品无人区一区二区三| 亚洲人成77777在线视频| 满18在线观看网站| 亚洲五月婷婷丁香| 国产麻豆69| 熟女少妇亚洲综合色aaa.| 国产精品一区二区在线不卡| 免费一级毛片在线播放高清视频 | 国产区一区二久久| 老司机午夜福利在线观看视频 | 日韩有码中文字幕| 少妇的丰满在线观看| 亚洲精品国产一区二区精华液| 中文字幕最新亚洲高清| 中文字幕高清在线视频| 老司机亚洲免费影院| 一边摸一边做爽爽视频免费| 韩国精品一区二区三区| 丝袜喷水一区| 亚洲av美国av| 精品免费久久久久久久清纯 | 男人爽女人下面视频在线观看| 亚洲中文字幕日韩| 狠狠精品人妻久久久久久综合| 免费黄频网站在线观看国产| 国产精品久久久人人做人人爽| 国产精品二区激情视频| 99精国产麻豆久久婷婷| 精品福利永久在线观看| 欧美日韩视频精品一区| 亚洲伊人久久精品综合| 一级片'在线观看视频| 看免费av毛片| 欧美国产精品va在线观看不卡| 国产精品一二三区在线看| 国产日韩欧美在线精品| 在线永久观看黄色视频| 亚洲,欧美精品.| 美女高潮到喷水免费观看| 亚洲欧美精品自产自拍| 久久香蕉激情| av又黄又爽大尺度在线免费看| 亚洲av电影在线观看一区二区三区| 色播在线永久视频| 少妇 在线观看| av线在线观看网站| 国产又色又爽无遮挡免| 欧美激情极品国产一区二区三区| 在线看a的网站| 亚洲天堂av无毛| 大码成人一级视频| 看免费av毛片| 亚洲午夜精品一区,二区,三区| 亚洲精品中文字幕在线视频| 午夜影院在线不卡| 精品人妻熟女毛片av久久网站| 国产免费福利视频在线观看| 欧美日韩一级在线毛片| 2018国产大陆天天弄谢| 美女福利国产在线| 国产无遮挡羞羞视频在线观看| 十八禁高潮呻吟视频| 亚洲伊人色综图| 国产欧美日韩一区二区三 | 成人国语在线视频| 中国国产av一级| av不卡在线播放| 欧美日韩亚洲高清精品| 在线观看免费视频网站a站| 亚洲人成77777在线视频| 汤姆久久久久久久影院中文字幕| 性少妇av在线| 久久人人爽av亚洲精品天堂| 两性午夜刺激爽爽歪歪视频在线观看 | 高清视频免费观看一区二区| 19禁男女啪啪无遮挡网站| 男女免费视频国产| 麻豆av在线久日| 一进一出抽搐动态| √禁漫天堂资源中文www| 午夜福利视频精品| 99热全是精品| 97在线人人人人妻| 亚洲久久久国产精品| 国产免费福利视频在线观看| 国产色视频综合| 成年女人毛片免费观看观看9 | 国产伦理片在线播放av一区| 亚洲一码二码三码区别大吗| 最新在线观看一区二区三区| 淫妇啪啪啪对白视频 | 国产精品成人在线| 亚洲精品av麻豆狂野| 亚洲专区国产一区二区| 韩国精品一区二区三区| 91精品三级在线观看| 欧美日韩视频精品一区| 两个人看的免费小视频| 五月开心婷婷网| 人成视频在线观看免费观看| e午夜精品久久久久久久| 天堂8中文在线网| 9热在线视频观看99| 自线自在国产av| 亚洲,欧美精品.| 亚洲国产欧美一区二区综合| 精品久久久久久电影网| 我的亚洲天堂| 丝袜脚勾引网站| 三上悠亚av全集在线观看| 男人添女人高潮全过程视频| 日本黄色日本黄色录像| 久久精品成人免费网站| 18禁裸乳无遮挡动漫免费视频| 久久 成人 亚洲| 99热全是精品| 美女中出高潮动态图| 精品久久久久久久毛片微露脸 | 日韩有码中文字幕| 高清av免费在线| 国产成人免费无遮挡视频| 一边摸一边做爽爽视频免费| 一区福利在线观看| 黄色视频不卡| 黄频高清免费视频| 嫁个100分男人电影在线观看| 亚洲伊人久久精品综合| 精品一区二区三区av网在线观看 | 久久久久久久精品精品| 欧美日韩精品网址| 国产野战对白在线观看| 一二三四社区在线视频社区8| 最近最新中文字幕大全免费视频| 免费在线观看黄色视频的| √禁漫天堂资源中文www| 欧美性长视频在线观看| 乱人伦中国视频| 欧美激情 高清一区二区三区| 国产成人精品久久二区二区91| 亚洲国产av新网站| 亚洲精品av麻豆狂野| 欧美日韩黄片免| 国产欧美日韩一区二区三 | av天堂久久9| 老司机在亚洲福利影院| 天天影视国产精品| 久久中文字幕一级| 99热网站在线观看| av国产精品久久久久影院| 亚洲三区欧美一区| 他把我摸到了高潮在线观看 | 久久性视频一级片| 黑人猛操日本美女一级片| 欧美国产精品va在线观看不卡| 国产亚洲av片在线观看秒播厂| 日韩视频一区二区在线观看| 亚洲色图 男人天堂 中文字幕| www.精华液| 免费在线观看视频国产中文字幕亚洲 | 老汉色av国产亚洲站长工具| 亚洲精品久久久久久婷婷小说| 国产精品免费大片| 欧美国产精品va在线观看不卡| 无遮挡黄片免费观看| 成人亚洲精品一区在线观看| 久9热在线精品视频| 欧美久久黑人一区二区| 激情视频va一区二区三区| 亚洲第一青青草原| 黄片大片在线免费观看| 两人在一起打扑克的视频| 多毛熟女@视频| www日本在线高清视频| 精品人妻1区二区| 午夜福利免费观看在线| 日本一区二区免费在线视频| 久久99热这里只频精品6学生| 999久久久精品免费观看国产| 18禁国产床啪视频网站| 丝袜人妻中文字幕| 欧美大码av| 午夜两性在线视频| 91麻豆精品激情在线观看国产 | 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 巨乳人妻的诱惑在线观看| 久久中文字幕一级| 亚洲精品乱久久久久久| 国产精品久久久久成人av| 涩涩av久久男人的天堂| 亚洲九九香蕉| 一本色道久久久久久精品综合| 高清欧美精品videossex| 国产日韩欧美在线精品| 男女免费视频国产| 2018国产大陆天天弄谢| 国产精品一区二区精品视频观看| 成年动漫av网址| 美女大奶头黄色视频| 午夜成年电影在线免费观看| 国产精品久久久久久精品电影小说| 亚洲精品久久成人aⅴ小说| 成在线人永久免费视频| 亚洲天堂av无毛| 我的亚洲天堂| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看完整版高清| 91国产中文字幕| 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 欧美日本中文国产一区发布| 久久香蕉激情| 99国产精品99久久久久| 国产在线一区二区三区精| 精品国产一区二区三区久久久樱花| av片东京热男人的天堂| 天天躁日日躁夜夜躁夜夜| 国产在视频线精品| 啦啦啦免费观看视频1| 精品熟女少妇八av免费久了| 久久久久久免费高清国产稀缺| 男女免费视频国产| 久久99热这里只频精品6学生| 国产一级毛片在线| 日本精品一区二区三区蜜桃| 两性夫妻黄色片| 啦啦啦中文免费视频观看日本| 嫩草影视91久久| 大陆偷拍与自拍| 日韩大片免费观看网站| 97人妻天天添夜夜摸| 80岁老熟妇乱子伦牲交| 超色免费av| 午夜免费鲁丝| 亚洲专区中文字幕在线| 99精国产麻豆久久婷婷| 成人三级做爰电影| 侵犯人妻中文字幕一二三四区| 欧美黄色片欧美黄色片| 视频区图区小说| 欧美性长视频在线观看| 另类亚洲欧美激情| 99久久精品国产亚洲精品| 在线观看www视频免费| 亚洲七黄色美女视频| 国产成人a∨麻豆精品| 久久中文看片网| 一级片'在线观看视频| 亚洲精品一区蜜桃| 精品国产一区二区久久| 自线自在国产av| 如日韩欧美国产精品一区二区三区| 欧美乱码精品一区二区三区| 在线观看免费视频网站a站| 一级片免费观看大全| 久久影院123| 婷婷色av中文字幕| 91精品三级在线观看| 国产成人一区二区三区免费视频网站| 国产在视频线精品| 国产日韩欧美在线精品| 国产福利在线免费观看视频| 国产精品 欧美亚洲| 亚洲精品第二区| 久热这里只有精品99| 伊人亚洲综合成人网| 亚洲第一青青草原| 中文欧美无线码| 亚洲欧美日韩高清在线视频 | 黄网站色视频无遮挡免费观看| 久热这里只有精品99| 丝袜人妻中文字幕| 精品一区二区三区av网在线观看 | 欧美国产精品va在线观看不卡| 国产欧美日韩一区二区三 | 一本综合久久免费| 国产av一区二区精品久久| 美女主播在线视频| 美国免费a级毛片| 少妇被粗大的猛进出69影院| 最黄视频免费看| 久久精品亚洲熟妇少妇任你| 欧美老熟妇乱子伦牲交| 日本一区二区免费在线视频| 国产成+人综合+亚洲专区| 日本精品一区二区三区蜜桃| 黄色a级毛片大全视频| 国产精品国产av在线观看| 欧美一级毛片孕妇| 日本欧美视频一区| 高清av免费在线| 欧美av亚洲av综合av国产av| 亚洲国产av影院在线观看| 久久久久久人人人人人| 午夜精品久久久久久毛片777| 国产亚洲一区二区精品| 久久国产精品大桥未久av| 欧美一级毛片孕妇| 国产精品免费视频内射| 亚洲欧美精品自产自拍| 一本久久精品| 肉色欧美久久久久久久蜜桃| 我要看黄色一级片免费的| 亚洲av日韩精品久久久久久密| av在线播放精品| 亚洲精品美女久久久久99蜜臀| 2018国产大陆天天弄谢| 久热这里只有精品99| 亚洲avbb在线观看| 久久亚洲精品不卡| 高清在线国产一区| 1024香蕉在线观看| 高清av免费在线| 国产一级毛片在线| 国产精品av久久久久免费| 首页视频小说图片口味搜索| 欧美日韩精品网址| 成年人免费黄色播放视频| 国产精品一区二区免费欧美 | 男男h啪啪无遮挡| 亚洲中文av在线| 亚洲精品久久午夜乱码| 美女高潮到喷水免费观看| 窝窝影院91人妻| 日韩三级视频一区二区三区| 久久 成人 亚洲| 五月开心婷婷网| 免费一级毛片在线播放高清视频 | 日本vs欧美在线观看视频| 国产深夜福利视频在线观看| 免费人妻精品一区二区三区视频| 日日摸夜夜添夜夜添小说| 日韩欧美国产一区二区入口| 香蕉丝袜av| av一本久久久久| 精品国产超薄肉色丝袜足j| 亚洲精品美女久久av网站| 免费不卡黄色视频| 色婷婷久久久亚洲欧美| 国产不卡av网站在线观看| 一区二区三区精品91| 黄色毛片三级朝国网站| 成年动漫av网址| a在线观看视频网站| 女人爽到高潮嗷嗷叫在线视频| 午夜福利乱码中文字幕| 亚洲一区二区三区欧美精品| 日日摸夜夜添夜夜添小说| 亚洲国产毛片av蜜桃av| 欧美亚洲 丝袜 人妻 在线| 热99国产精品久久久久久7| 久久人妻熟女aⅴ| 中文字幕精品免费在线观看视频| cao死你这个sao货| 侵犯人妻中文字幕一二三四区| 久久久久国内视频| 热re99久久精品国产66热6| 如日韩欧美国产精品一区二区三区| 老汉色∧v一级毛片| av福利片在线| 亚洲欧美激情在线| 欧美日韩中文字幕国产精品一区二区三区 | 人妻人人澡人人爽人人| 国产一区二区三区综合在线观看| 国产男女内射视频| 国产免费av片在线观看野外av| 亚洲欧美精品自产自拍| 国产精品1区2区在线观看. | 欧美黄色片欧美黄色片| 亚洲精品av麻豆狂野| 女性生殖器流出的白浆| 亚洲专区字幕在线| 老鸭窝网址在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品高清国产在线一区| 美女高潮喷水抽搐中文字幕| 少妇猛男粗大的猛烈进出视频| 国产成人精品在线电影| 免费在线观看视频国产中文字幕亚洲 | 亚洲av电影在线观看一区二区三区| 人妻一区二区av| 一级毛片女人18水好多| 国产男女超爽视频在线观看| 久久精品亚洲熟妇少妇任你| 亚洲一区二区三区欧美精品| 婷婷成人精品国产| 少妇人妻久久综合中文| 日韩电影二区| 久久国产精品人妻蜜桃| 亚洲中文av在线| 91九色精品人成在线观看| 国产伦人伦偷精品视频| 久久国产精品影院| 别揉我奶头~嗯~啊~动态视频 | 婷婷丁香在线五月| 精品亚洲乱码少妇综合久久| 欧美精品一区二区免费开放| 在线十欧美十亚洲十日本专区| 亚洲精品国产精品久久久不卡| 久久午夜综合久久蜜桃| 精品一品国产午夜福利视频| 性少妇av在线| 免费在线观看影片大全网站| a在线观看视频网站| 一二三四在线观看免费中文在| 又大又爽又粗| 999久久久国产精品视频| 女人精品久久久久毛片| 国产主播在线观看一区二区| 色视频在线一区二区三区| 亚洲精品国产精品久久久不卡| 成人国产av品久久久| kizo精华| 色视频在线一区二区三区| 久久久国产精品麻豆| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 国产精品国产三级国产专区5o| 久久香蕉激情| 少妇精品久久久久久久| 国产人伦9x9x在线观看| 五月开心婷婷网| 国产精品久久久久成人av| 最近最新免费中文字幕在线| 12—13女人毛片做爰片一| 99久久精品国产亚洲精品| 欧美97在线视频| 青春草亚洲视频在线观看| 亚洲精品美女久久av网站| 咕卡用的链子| 午夜福利在线观看吧| 蜜桃国产av成人99| 亚洲av男天堂| 国产精品久久久人人做人人爽| 日本av手机在线免费观看| 精品久久久久久久毛片微露脸 | 俄罗斯特黄特色一大片| 国产精品一区二区精品视频观看| 欧美人与性动交α欧美精品济南到| 十八禁高潮呻吟视频| 国产男女内射视频| 捣出白浆h1v1| 啦啦啦视频在线资源免费观看| 日韩,欧美,国产一区二区三区| 9色porny在线观看| 国产av一区二区精品久久| 老司机福利观看| 极品少妇高潮喷水抽搐| 久9热在线精品视频| 亚洲专区国产一区二区| 国产成人精品久久二区二区91| 超碰成人久久| 欧美激情高清一区二区三区| av在线播放精品| 深夜精品福利| 成人18禁高潮啪啪吃奶动态图| 亚洲第一av免费看| 国产成人一区二区三区免费视频网站| 丰满饥渴人妻一区二区三| av网站免费在线观看视频| 一区二区三区乱码不卡18| 亚洲精品粉嫩美女一区| 成人影院久久| 精品国产一区二区三区四区第35| 亚洲免费av在线视频| 国产欧美日韩一区二区三区在线| 国产精品偷伦视频观看了| 久久国产精品影院| 亚洲精华国产精华精| 午夜福利影视在线免费观看| 一个人免费在线观看的高清视频 | 麻豆乱淫一区二区| 人成视频在线观看免费观看| 美国免费a级毛片| videosex国产| 精品国产乱子伦一区二区三区 | 国产男人的电影天堂91| 99久久人妻综合| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品国产三级国产专区5o| 亚洲黑人精品在线| 老司机影院成人| 国产伦人伦偷精品视频| 精品少妇内射三级| 国产在线免费精品| 中文字幕av电影在线播放| 在线精品无人区一区二区三| 中文字幕精品免费在线观看视频| 91国产中文字幕| 国产亚洲av片在线观看秒播厂| 欧美大码av| 大陆偷拍与自拍| 一边摸一边做爽爽视频免费| av片东京热男人的天堂| 欧美老熟妇乱子伦牲交| 一本综合久久免费| 两性夫妻黄色片| 国产精品麻豆人妻色哟哟久久| 国产精品一区二区免费欧美 | 午夜两性在线视频| 国产97色在线日韩免费| 亚洲男人天堂网一区| 午夜福利影视在线免费观看| 成人手机av| 欧美精品高潮呻吟av久久| 我的亚洲天堂| 亚洲一码二码三码区别大吗| 纯流量卡能插随身wifi吗| 国产激情久久老熟女| www.av在线官网国产| 亚洲一区中文字幕在线| 精品国产乱码久久久久久小说| avwww免费| 在线观看免费高清a一片| 一级片免费观看大全| 亚洲av成人不卡在线观看播放网 | 成人亚洲精品一区在线观看| 久久久久国产精品人妻一区二区| 三级毛片av免费| 热99久久久久精品小说推荐| 王馨瑶露胸无遮挡在线观看| 超色免费av| 欧美av亚洲av综合av国产av| 日本av手机在线免费观看| 高清黄色对白视频在线免费看| 在线亚洲精品国产二区图片欧美| 国产精品欧美亚洲77777| 一个人免费在线观看的高清视频 | 国产xxxxx性猛交| 国产成人欧美在线观看 | 老熟女久久久| 青春草亚洲视频在线观看| 肉色欧美久久久久久久蜜桃| 黄色 视频免费看| 日本猛色少妇xxxxx猛交久久| videosex国产| 亚洲中文字幕日韩| 成年女人毛片免费观看观看9 | 精品国产一区二区久久| 精品亚洲成国产av| 狂野欧美激情性xxxx| 两个人免费观看高清视频| 久久热在线av| 久久天堂一区二区三区四区| 性高湖久久久久久久久免费观看| 99久久99久久久精品蜜桃| 51午夜福利影视在线观看| 黄色视频在线播放观看不卡| 在线观看免费日韩欧美大片| 国产精品久久久久久精品古装| 精品国产乱码久久久久久小说| 国产xxxxx性猛交| 啦啦啦视频在线资源免费观看| 麻豆av在线久日| av在线老鸭窝| 啦啦啦在线免费观看视频4| 亚洲激情五月婷婷啪啪| 丝袜喷水一区| 黄色毛片三级朝国网站| 18在线观看网站| 久久国产精品男人的天堂亚洲| 真人做人爱边吃奶动态| 亚洲人成电影观看| 欧美97在线视频| 色播在线永久视频| 成年人黄色毛片网站| h视频一区二区三区| 国精品久久久久久国模美| 国产黄色免费在线视频| 国产日韩欧美在线精品| 午夜福利视频精品| 桃红色精品国产亚洲av| 国产精品久久久av美女十八| 91大片在线观看| 久久女婷五月综合色啪小说| 女性被躁到高潮视频| 午夜福利在线观看吧| 999精品在线视频| 99久久精品国产亚洲精品| 久久亚洲精品不卡| 国产精品一区二区在线不卡| 捣出白浆h1v1| 一本综合久久免费|