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

    石化管道丁烷氣體泄漏擴(kuò)散數(shù)值模擬

    2020-05-21 05:36:28王偉天黃維秋趙會(huì)軍
    安全與環(huán)境工程 2020年2期
    關(guān)鍵詞:丁烷射流組分

    周 寧,王偉天,陳 兵,李 雪,陳 力,黃維秋,趙會(huì)軍

    (1.常州大學(xué)石油工程學(xué)院,江蘇 常州 213164;2.中國(guó)安全生產(chǎn)科學(xué)研究院,北京 100012)

    丁烷作為燃料與原料大量用于工業(yè)生產(chǎn),為易燃、易爆、有毒的氣體,發(fā)生泄漏后形成氣云,易導(dǎo)致火災(zāi)、爆炸事故的發(fā)生,造成人員中毒傷亡。例如:2013年6月,云南硯山某公司的廠區(qū)管道破裂,發(fā)生丁烷氣體泄漏并引發(fā)火災(zāi),造成5人傷亡;2016年11月,位于美國(guó)路易斯安那州的Baton Rouge煉油廠發(fā)生了丁烷氣體泄漏,丁烷氣云遇明火源引發(fā)大面積火災(zāi),造成大量設(shè)備燒毀以及4人嚴(yán)重?zé)齻?。因此,研究丁烷氣體的泄漏擴(kuò)散規(guī)律對(duì)于安全生產(chǎn)意義重大。

    針對(duì)丁烷氣體泄漏擴(kuò)散行為的研究主要包括現(xiàn)場(chǎng)試驗(yàn)研究、物理模擬研究和數(shù)值模擬研究[1]。近年來(lái),為研究重氣泄漏擴(kuò)散特征,國(guó)內(nèi)外學(xué)者針對(duì)各類重質(zhì)氣體泄漏擴(kuò)散的一般規(guī)律進(jìn)行了研究[2-3]。如潘旭海等[4]從控制重氣云團(tuán)擴(kuò)散行為的微分方程入手,根據(jù)箱模型以及其他重氣擴(kuò)散模型采用向前插分和牛頓迭代的方法進(jìn)行了數(shù)值模擬,得到重氣云團(tuán)外形尺寸和空氣卷吸量隨時(shí)間的變化規(guī)律以及重氣云團(tuán)在下風(fēng)向固定點(diǎn)處地面的最大濃度值,較好地反映了重氣云團(tuán)特有的擴(kuò)散行為;吳玉劍等[5]利用計(jì)算流體力學(xué)的方法對(duì)Thorney Island Trial026實(shí)驗(yàn)進(jìn)行了數(shù)值模擬,驗(yàn)證了CFD軟件能較精確地模擬障礙物地形下的重氣擴(kuò)散過(guò)程;姜傳勝等[6]開展了重氣連續(xù)泄漏擴(kuò)散的風(fēng)洞模擬實(shí)驗(yàn),剖析了重氣連續(xù)泄漏擴(kuò)散的特點(diǎn),重點(diǎn)分析了環(huán)境風(fēng)速對(duì)重氣連續(xù)泄漏擴(kuò)散的影響;Chavez等[7]采用風(fēng)洞模擬試驗(yàn)和CFD軟件研究了六氟化硫在建筑物群的擴(kuò)散行為。目前,關(guān)于石化裝置及其相關(guān)工藝管道氣體泄漏擴(kuò)散規(guī)律方面的研究較多,而對(duì)化工園區(qū)周邊管廊管道氣體泄漏擴(kuò)散規(guī)律以及擴(kuò)散云團(tuán)對(duì)周圍環(huán)境的影響的研究較少,因此本文利用CFD軟件對(duì)某石化園區(qū)化工管廊管道丁烷氣體的泄漏擴(kuò)散規(guī)律進(jìn)行了數(shù)值模擬研究。

    1 數(shù)值模型的建立

    1. 1 幾何模型

    上海某化工園區(qū)的某段管廊全寬為9 m,總高為10.9 m,共三層管架,其中選作研究對(duì)象的編號(hào)55號(hào)管道位于第二層管架,距地面6.5 m,距管架左側(cè)邊緣2 m,管徑為0.078 m,壁厚為5.5 mm。根據(jù)石化園區(qū)化工管廊管道現(xiàn)場(chǎng)的實(shí)際情況并考慮泄漏氣體的充分流動(dòng),選取計(jì)算區(qū)域長(zhǎng)40 m、寬20 m,管道直徑為0.078 m,管道距離左側(cè)入口邊界為2 m,管道底部距離地面為6.5 m,管道位置可以保證流動(dòng)充分發(fā)展[8]。泄漏口位于管道頂部中間,泄漏口直徑為0.02 m。建立的幾何模型如圖1所示。

    圖1 幾何模型Fig.1 Geometric Model

    1. 2 數(shù)學(xué)模型

    丁烷氣體泄漏擴(kuò)散的過(guò)程遵循連續(xù)性方程、動(dòng)量守恒方程、能量守恒方程和湍流控制方程,考慮到丁烷氣體泄漏后在大氣中的擴(kuò)散行為,本次數(shù)值模擬計(jì)算將啟動(dòng)組分?jǐn)U散模型[9-10]。具體方程如下:

    (1) 連續(xù)性方程為

    (1)

    式中:ρ為流體的密度(kg/m3);t為時(shí)間(s);u為速度矢量(m/s)。

    (2) 動(dòng)量守恒方程為

    (2)

    (3) 能量守恒方程為

    (3)

    (4) 湍流控制方程Realizablek-ε模型比標(biāo)準(zhǔn)k-ε模型在濃度分布計(jì)算上有更好的精度[11],故對(duì)氣體擴(kuò)散湍流問(wèn)題采用Realizablek-ε模型,其湍流方程為

    (4)

    (5)

    式中:k為湍流動(dòng)能(J);ε為耗散率(%);xi、xj分別為各坐標(biāo)分量;C1ε、C3ε、C1為常數(shù);σk為湍流動(dòng)能k的湍流普朗特?cái)?shù),取值為1.0;σε為耗散率ε的湍流普朗特?cái)?shù),取值為1.2;C2為常數(shù),取值為1.9;Gk為因速度梯度產(chǎn)生的湍流動(dòng)能源項(xiàng)(J);Gb為因浮力產(chǎn)生的湍流動(dòng)能源項(xiàng)(J);YM為在可壓縮湍流中波動(dòng)擴(kuò)張引起的耗散項(xiàng)(%);?為運(yùn)動(dòng)黏度(m2/s)。

    (5) 組分?jǐn)U散模型為

    (6)

    1. 3 網(wǎng)格劃分與邊界條件設(shè)置

    管道泄漏口設(shè)置于管道頂部中間,泄漏口直徑為0.02 m,屬小孔泄漏。管內(nèi)溫度取管道運(yùn)行時(shí)的溫度即313.16 K,環(huán)境溫度為300 K,環(huán)境壓力為101 325 Pa。為了便于計(jì)算,管道氣體泄漏過(guò)程簡(jiǎn)化為連續(xù)泄漏且管道內(nèi)部壓力不隨泄漏過(guò)程改變。

    考慮到泄漏口附近壓力變化較大、流場(chǎng)復(fù)雜,為了保證計(jì)算的精度,計(jì)算模型在管道泄漏口及其附近進(jìn)行網(wǎng)格局部加密(見圖2),泄漏口處加密網(wǎng)格尺寸為0.008 mm×0.008 mm,其余網(wǎng)格尺寸為0.09 mm×0.09 mm,計(jì)算模型總網(wǎng)格數(shù)為314 420,90%的網(wǎng)格質(zhì)量不低于0.9,滿足本文的計(jì)算精度要求。

    泄漏口采用速度入口,左側(cè)入口邊界設(shè)置速度入口,計(jì)算區(qū)域頂部及右側(cè)邊界設(shè)置為壓力出口,管壁和底部地面為無(wú)滑移的恒溫壁面。計(jì)算采用Pressure-Base求解器求解,選擇瞬態(tài)計(jì)算,湍流模型選取Realizablek-ε模型,啟動(dòng)能量守恒方程,選擇組分輸運(yùn)模型。

    圖2 泄漏口處網(wǎng)格加密情況Fig.2 Grid encryption at the leak

    2 數(shù)值模擬結(jié)果與分析

    2. 1 無(wú)風(fēng)狀態(tài)下丁烷氣體的泄漏擴(kuò)散規(guī)律

    管道運(yùn)行壓力為1.8 MPa,與大氣壓力相差約16倍,當(dāng)丁烷氣體從泄漏口以射流形式進(jìn)入大氣后其與環(huán)境壓力迅速達(dá)到相對(duì)平衡,平衡區(qū)域較小且集中于泄漏口附近,故對(duì)整個(gè)計(jì)算區(qū)域造成的影響較小。無(wú)風(fēng)狀態(tài)下丁烷氣體泄漏穩(wěn)定后的速度和組分濃度分布云圖,見圖3。

    圖3 無(wú)風(fēng)狀態(tài)下丁烷氣體泄漏穩(wěn)定后的速度和組分 濃度分布云圖Fig.3 Cloud map of butane velocity and concentration distribution when the leakage is stable under the windless condition

    由圖3可見,當(dāng)丁烷氣體泄漏的自由擴(kuò)散穩(wěn)定后,速度和組分濃度分布均呈對(duì)稱關(guān)系。由于丁烷氣體泄漏的初始狀態(tài)為高壓力高濃度射流,在泄漏口處組分濃度和速度達(dá)到最大值,隨后逐漸擴(kuò)散,最終與大氣混合。

    根據(jù)丁烷氣體泄漏速度分布云圖[見圖3(a)]可以看出:由于丁烷氣體在泄漏初始狀態(tài)時(shí)速度最大,丁烷氣體以射流的形式從泄漏口噴射而出,并且在泄漏口上方形成錐形的射流流場(chǎng),當(dāng)丁烷氣體高速噴射到大氣后,迅速與大氣介質(zhì)均勻摻混,同時(shí)向周圍擴(kuò)散。

    此外,由圖3還可以看出,丁烷氣體泄漏為無(wú)固壁約束的自由湍流,由于存在壓力和速度差,丁烷氣體的湍性射流迅速與周圍邊界中活躍的湍流進(jìn)行混合,對(duì)周圍介質(zhì)產(chǎn)生吸卷,發(fā)生渦流現(xiàn)象。而在泄漏口的附近,丁烷氣體泄漏的射流作用使管壁兩側(cè)的壓力下降,形成局部低壓,產(chǎn)生“負(fù)壓區(qū)”(見圖4),最終導(dǎo)致丁烷泄漏氣體往射流邊界層發(fā)生橫向流動(dòng)(見圖5)。由于丁烷氣體進(jìn)入大氣后對(duì)空氣發(fā)生吸卷效應(yīng),使其射流的最大速度隨著距泄漏口距離的增大而逐漸減小,丁烷氣體泄漏射流的最大速度處于中心線上,其與距泄漏口的距離的平方根成反比。

    圖4 丁烷氣體泄漏近場(chǎng)負(fù)壓特征區(qū)Fig.4 Negative pressure feature area near field of the leakage of butane gas

    圖5 無(wú)風(fēng)狀態(tài)下丁烷氣體泄漏射流邊界層擴(kuò)大特征 云圖Fig.5 Profiling feature cloud map of jet boundary layer of butane gas leakage under the windless condition

    另外,由圖3(a)可以觀察到丁烷氣體泄漏后產(chǎn)生湍性射流的幾個(gè)特征:丁烷氣體與大氣介質(zhì)的均勻摻混呈現(xiàn)出自由射流的自模性[12],主要表現(xiàn)為在不同的截面上丁烷氣體與大氣介質(zhì)的混合長(zhǎng)度與射流寬度呈正比關(guān)系,氣流參數(shù)分布規(guī)律相似,而整體的氣體擴(kuò)散分布規(guī)律表現(xiàn)為線性漸進(jìn);氣體微團(tuán)由于受到湍性射流所產(chǎn)生的渦流的影響,發(fā)生橫向的質(zhì)量、熱量和動(dòng)量交換,伴隨著這些交換的發(fā)生而產(chǎn)生邊界層,邊界層的存在導(dǎo)致射流擴(kuò)展的同時(shí)更容易對(duì)周圍的大氣介質(zhì)產(chǎn)生吸卷[13];由于泄漏的丁烷氣體與環(huán)境大氣發(fā)生的吸卷效應(yīng),重力作用使射流出現(xiàn)了彎曲現(xiàn)象,這種現(xiàn)象在距泄漏口越遠(yuǎn)的地方越明顯,該射流屬于非等密度射流;氣體射流寬度與距泄漏口(即射流源)的距離呈正比關(guān)系[14]。

    丁烷氣體泄漏擴(kuò)散的組分濃度分布規(guī)律相似于速度分布規(guī)律,兩者均呈對(duì)稱分布,且存在邊界層。氣體射流中心區(qū)域?yàn)榧兌⊥闅怏w,射流邊界層以外為環(huán)境大氣介質(zhì),邊界層屬于丁烷氣體與大氣介質(zhì)的混合物。丁烷氣體的流動(dòng)性質(zhì)與純射流較為相近,表現(xiàn)在濃度分布中射流核心區(qū)域的射流氣體幾乎不受大氣介質(zhì)的影響,其濃度極高。由于射流寬度與距射流源的距離呈正比關(guān)系,伴隨射流高度的增加,射流氣體的濃度分布范圍也隨著射流寬度的增加而增大[15]。這是因?yàn)榇髿饨橘|(zhì)受到射流過(guò)程的卷吸而對(duì)其產(chǎn)生了阻滯作用,導(dǎo)致射流氣體速度下降,初始動(dòng)量減少,當(dāng)射流氣體受到浮力影響時(shí)其射流速度逐漸增大,擴(kuò)散范圍也逐漸擴(kuò)大,最終在重力的作用下發(fā)生沉降[16]。

    丁烷氣體爆炸上限為8.4%,爆炸下限為1.9%,取爆炸下限的二分之一為警戒濃度(0.95%),通過(guò)Tecplot軟件對(duì)丁烷氣體泄漏組分濃度分布云圖進(jìn)行處理,得到無(wú)風(fēng)狀態(tài)下丁烷氣體泄漏組分濃度等值線圖(見圖6),并對(duì)各區(qū)域面積進(jìn)行估算。

    圖6中,紅色區(qū)域?yàn)槎⊥闅怏w高濃度區(qū)域,丁烷氣體組分濃度高于其爆炸上限8.4%,區(qū)域面積約為80 m2;黃色區(qū)域?yàn)槎⊥闅怏w爆炸極限區(qū)域,丁烷氣體組分濃度區(qū)間為[1.9%,8.4%],介于其爆炸上、下限之間,區(qū)域面積約為35 m2;綠色區(qū)域?yàn)榫鋮^(qū)域,丁烷氣體組分濃度區(qū)間為[0.95%,1.9%],介于其警戒濃度與爆炸下限之間,區(qū)域面積約為16 m2。由于丁烷氣體泄漏的射流速度和壓力較大,其密度高于空氣密度,射流氣體的核心區(qū)域濃度變化緩慢,外圍區(qū)域氣體則迅速與大氣摻混,使其濃度降低。在初始動(dòng)能被大量消耗、氣體達(dá)到射流高度峰值時(shí),丁烷氣云仍然保持高濃度狀態(tài),產(chǎn)生渦卷現(xiàn)象的區(qū)域氣體濃度達(dá)到8.4%以上。

    2. 2 環(huán)境風(fēng)速對(duì)丁烷氣體泄漏擴(kuò)散過(guò)程的影響

    研究不同風(fēng)速對(duì)丁烷泄漏氣體濃度分布的影響對(duì)研究丁烷氣體泄漏擴(kuò)散規(guī)律有著重要的意義。通過(guò)數(shù)值模擬得到不同環(huán)境風(fēng)速(分別為3 m/s、5 m/s、8 m/s時(shí))、丁烷氣體泄漏初始速度為30 m/s狀態(tài)下其濃度分布云圖,見圖7。

    圖7 不同環(huán)境風(fēng)速下丁烷泄漏氣體組分濃度場(chǎng)分布云圖Fig.7 Clound map of concentration field distribution of butane gas leakage with different ambient wind speeds

    由圖7可見,丁烷氣體泄漏的濃度場(chǎng)往下風(fēng)向徑向偏斜率受環(huán)境風(fēng)速的影響,環(huán)境風(fēng)速增大使其射流偏離垂直方向的角度隨之增大;在泄漏口處丁烷氣體出流速度極高,遠(yuǎn)大于環(huán)境風(fēng)速,由于兩者速度差較大使泄漏口處的丁烷氣體射流動(dòng)壓力受風(fēng)流產(chǎn)生的動(dòng)壓力的影響微弱。

    不同環(huán)境風(fēng)速下管道泄漏孔處丁烷氣體泄漏射流無(wú)偏移特征區(qū)云圖,見圖8。

    圖8 不同環(huán)境風(fēng)速下管道泄漏孔處丁烷氣體泄漏 射流無(wú)偏移特征區(qū)云圖Fig.8 Cloud map of the no offset feature area of the butane gas leakage jet at the pipe leak hole with different ambient wind speeds

    由圖8可見,丁烷氣體位于泄漏口處的濃度場(chǎng)沒(méi)有發(fā)生明顯偏移,但隨著環(huán)境風(fēng)速的增大,無(wú)偏移的氣體射流縮短;而當(dāng)射流距離逐漸增大后,氣體射流速度下降,空氣阻力的作用逐漸超過(guò)氣體射流動(dòng)壓力的作用,同時(shí)其射流受到大氣湍流的影響也逐漸增加,氣體擴(kuò)散范圍擴(kuò)大,風(fēng)流的動(dòng)壓力逐步增強(qiáng)且影響作用更加明顯;在環(huán)境風(fēng)速的作用下丁烷氣體的濃度場(chǎng)往下風(fēng)向轉(zhuǎn)移,氣體擴(kuò)散受環(huán)境風(fēng)速的影響加強(qiáng),且丁烷氣云受到環(huán)境風(fēng)流的影響而連續(xù)性被破壞,高濃度區(qū)域的丁烷氣云被割裂成小型氣團(tuán),環(huán)境風(fēng)速越大,切割后形成的氣團(tuán)越小。

    采用Tecplot軟件對(duì)不同環(huán)境風(fēng)速下丁烷氣體泄漏濃度場(chǎng)分布云圖進(jìn)行處理,得到不同風(fēng)速下丁烷氣體泄漏組分濃度等值線圖(見圖9),并對(duì)各區(qū)域面積進(jìn)行估算。

    圖9 不同環(huán)境風(fēng)速下丁烷氣體泄漏組分濃度等值線圖Fig.9 Contours of concentration of butane gas leakage with different ambient wind speeds

    由圖9可見,當(dāng)環(huán)境風(fēng)速為3 m/s時(shí),丁烷氣體泄漏高濃度區(qū)域(丁烷氣體組分濃度>爆炸上限8.4%)面積約為66 m2,爆炸極限區(qū)域(丁烷氣體組分濃度介于爆炸上限8.4%、爆炸下限1.9%之間)面積約為128 m2,警戒區(qū)域(丁烷氣體組分濃度介于爆炸下限1.9%、警戒濃度0.95%之間)面積約為18 m2[見圖9(a)];當(dāng)環(huán)境風(fēng)速為5 m/s時(shí),丁烷氣體泄漏高濃度區(qū)域面積約為64 m2,爆炸極限區(qū)域面積約為234 m2,警戒區(qū)域面積約為30 m2[見圖9(b)];當(dāng)環(huán)境風(fēng)速為8 m/s時(shí),丁烷氣體泄漏高濃度區(qū)域面積約為60 m2,爆炸極限區(qū)域面積約為242 m2,警戒區(qū)域面積約為35 m2[見圖9(c)]。同時(shí),由圖9還可以發(fā)現(xiàn):隨著環(huán)境風(fēng)速的增大,丁烷云團(tuán)水平移動(dòng)距離增加,相同水平截面的云團(tuán)高度下降,說(shuō)明環(huán)境風(fēng)對(duì)氣云的擴(kuò)散有明顯的促進(jìn)作用,隨著環(huán)境風(fēng)速的增大,氣云整體擴(kuò)散面積增加,高濃度區(qū)域面積逐漸減小,而濃度相對(duì)較低的爆炸極限區(qū)域和警戒區(qū)域的面積不斷增大。

    2. 3 泄漏初始速度對(duì)丁烷氣體泄漏擴(kuò)散過(guò)程的影響

    對(duì)丁烷氣體泄漏初始速度分別為30 m/s、40 m/s、50 m/s和60 m/s 4工況下丁烷氣體的泄漏擴(kuò)散過(guò)程進(jìn)行了數(shù)值模擬,環(huán)境風(fēng)速設(shè)置為8 m/s,得到丁烷氣體泄漏濃度場(chǎng)分布云圖,見圖10。

    圖10 環(huán)境風(fēng)速為8 m/s時(shí)不同泄漏初始速度下丁烷 氣體泄漏組分濃度分布云圖Fig.10 Cloud map of concentration field distribution of butane gas leakage with different leaking initial velocities when the wind speed is 8 m/s

    由圖10可見,當(dāng)環(huán)境風(fēng)速相同時(shí),丁烷氣體泄漏初始速度越高,氣體射流動(dòng)壓力越大,氣體射流受風(fēng)流動(dòng)壓力的影響越微弱,氣體射流發(fā)生的軸向偏移量越小。這是因?yàn)楫?dāng)環(huán)境風(fēng)速不變時(shí),丁烷氣體泄漏初始速度較小,風(fēng)力對(duì)氣體射流產(chǎn)生束縛作用,氣體射流的動(dòng)壓力受到風(fēng)流動(dòng)壓力的作用更強(qiáng),使氣團(tuán)不易受到大氣湍流的影響,擴(kuò)散作用減弱,擴(kuò)散范圍縮小,徑向偏斜率增大,徑向擴(kuò)散距離增加;當(dāng)丁烷氣體泄漏初始速度增大,氣體射流將在較低的高度發(fā)生擴(kuò)散行為,徑向偏斜率和擴(kuò)散距離減小,此時(shí)丁烷氣體的湍性射流更容易與周圍大氣湍流混合,對(duì)靜止大氣介質(zhì)產(chǎn)生卷吸,發(fā)生渦流現(xiàn)象。

    采用Tecplot軟件對(duì)不同泄漏初始速度下丁烷組分濃度分布云圖進(jìn)行處理,得到不同泄漏初始速度下丁烷氣體泄漏組分濃度等值線圖(見圖11),并對(duì)各區(qū)域面積進(jìn)行估算。

    圖11 不同泄漏初始速度下丁烷氣體泄漏組分濃度 等值線圖Fig.11 Contours of concentration of butane gas leakage with different leaking initial velocities

    由圖11可見,當(dāng)丁烷氣體泄漏初始速度為30 m/s時(shí),丁烷氣體高濃度區(qū)域(丁烷氣體組分濃度高于爆炸上限8.4%)面積約為44 m2,爆炸極限區(qū)域(丁烷氣體組分濃度介于爆炸上限8.4%與爆炸下限1.9%之間)面積約為162 m2,警戒區(qū)域(丁烷氣體組分濃度介于爆炸下限1.9%與警戒濃度0.95%之間)面積約為30 m2[見圖11(a)];當(dāng)丁烷氣體泄漏初始速度為40 m/s時(shí),丁烷氣體高濃度區(qū)域面積約為60 m2,爆炸極限區(qū)域面積約為242 m2,警戒區(qū)域面積約為35 m2;當(dāng)丁烷氣體泄漏初始速度為50 m/s時(shí),丁烷氣體高濃度區(qū)域面積約為87 m2,爆炸極限區(qū)域面積約為284 m2,警戒區(qū)域面積約為40 m2;當(dāng)丁烷氣體泄漏初始速度為60 m/s時(shí),丁烷氣體高濃度區(qū)域面積約為128 m2,爆炸極限區(qū)域面積約為427 m2,警戒區(qū)域面積約為58 m2。同時(shí),由圖11還可以看出,隨著丁烷氣體泄漏初始速度的增大,氣云重量增加,氣云內(nèi)部分子壓強(qiáng)增大,風(fēng)力對(duì)云團(tuán)擴(kuò)散的影響不斷衰減,氣體自由擴(kuò)散作用隨之增強(qiáng)并逐漸成為云團(tuán)擴(kuò)散的主導(dǎo)因素,最終導(dǎo)致氣云的縱深增大,氣云同時(shí)向高空與地面蔓延,各濃度區(qū)域面積同時(shí)擴(kuò)大,其中爆炸極限區(qū)域面積的擴(kuò)大程度最強(qiáng)。由于高濃度區(qū)域空氣不足,丁烷氣體難以燃燒,低濃度區(qū)域丁烷氣體不足以引起燃燒,而爆炸極限區(qū)域的丁烷氣體易燃、易爆,與空氣混合后不穩(wěn)定,容易引發(fā)火災(zāi)、爆炸等事故,故爆炸極限區(qū)域?yàn)槎⊥闅怏w泄漏的主要危險(xiǎn)區(qū)域[14-15]。隨著丁烷氣體泄漏初始速度的增加,危險(xiǎn)區(qū)域面積增大并往地面貼近,因此丁烷氣體泄漏初始速度越大,氣云縱深越長(zhǎng),丁烷氣體泄漏危險(xiǎn)程度越高。該模擬結(jié)果符合氣體自由擴(kuò)散泄漏的表達(dá)式,即氣體泄漏質(zhì)量流量越大,擴(kuò)散面積越大。其表達(dá)式為

    (7)

    式中:Qm為泄漏氣體的質(zhì)量流量(kg/s);C0為泄漏氣體的流出系數(shù),無(wú)量綱;A為泄漏氣體的擴(kuò)散面積(m2);P0為容器內(nèi)壓(Pa);P為環(huán)境壓力(Pa);M為泄漏氣體的分子質(zhì)量,單位為1;Rg為理想氣體常數(shù)(J/mol·K);γ為泄漏氣體的絕熱系數(shù)(無(wú)量綱);T0為泄漏源的溫度(K)。

    3 結(jié) 論

    通過(guò)對(duì)某石化園區(qū)化工管廓管道丁烷氣體泄漏擴(kuò)散規(guī)律進(jìn)行數(shù)值模擬,得到以下結(jié)論:

    (1) 無(wú)風(fēng)狀態(tài)下,丁烷氣體從泄漏口以無(wú)固壁約束的自由射流形式進(jìn)入外部大氣環(huán)境,射流氣體的速度和濃度分布均呈對(duì)稱分布,分布規(guī)律較為相似。射流氣體離開泄漏口后受浮力作用,氣體流速表現(xiàn)出先增后減的規(guī)律,越靠近軸線該現(xiàn)象越明顯。

    (2) 環(huán)境風(fēng)速使丁烷氣體射流發(fā)生徑向偏斜,環(huán)境風(fēng)速越大,徑向偏斜率越大,射流氣體的偏斜點(diǎn)與泄漏口之間的距離越短;環(huán)境風(fēng)速對(duì)丁烷氣云擴(kuò)散有促進(jìn)作用,隨著環(huán)境風(fēng)速的增大,氣云水平擴(kuò)散距離增加,相同水平截面上氣云的擴(kuò)散高度減小;環(huán)境風(fēng)速對(duì)丁烷氣云產(chǎn)生割裂作用,環(huán)境風(fēng)速越大,氣云連續(xù)性受破壞的程度越大,割裂后形成的氣團(tuán)體積越小、數(shù)量越多;環(huán)境風(fēng)速與丁烷氣體泄漏擴(kuò)散引起的爆炸危險(xiǎn)區(qū)域面積呈正比關(guān)系,環(huán)境風(fēng)速增大,丁烷氣體泄漏擴(kuò)散引起的爆炸危險(xiǎn)區(qū)域的水平距離增加,其擴(kuò)散高度減小,且環(huán)境風(fēng)速越大,丁烷氣體泄漏擴(kuò)散引起的爆炸危險(xiǎn)區(qū)域主要分布位置往下風(fēng)向移動(dòng)距離越大,而其擴(kuò)散高度維持在5~15 m的高度上。

    (3) 丁烷氣體泄漏初始速度與丁烷氣體泄漏擴(kuò)散引起的爆炸危險(xiǎn)區(qū)域面積呈正比關(guān)系,丁烷氣體泄漏速率增大,丁烷氣體泄漏擴(kuò)散引起的爆炸危險(xiǎn)區(qū)域水平距離不受影響,但氣云縱深增加,其爆炸危險(xiǎn)區(qū)域面積擴(kuò)大。

    猜你喜歡
    丁烷射流組分
    深海逃逸艙射流注水均壓過(guò)程仿真分析
    低壓天然氣泄漏射流擴(kuò)散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    組分分發(fā)管理系統(tǒng)在天然氣計(jì)量的應(yīng)用
    一種難溶難熔未知組分板材的定性分析
    在線處理脫丁烷塔壓差高
    黑順片不同組分對(duì)正常小鼠的急性毒性
    中成藥(2018年8期)2018-08-29 01:28:26
    金雀花中黃酮苷類組分鑒定及2種成分測(cè)定
    中成藥(2018年2期)2018-05-09 07:20:09
    新型3-氨基氮雜環(huán)丁烷衍生物的合成
    正丁烷氧化制順酐催化劑的制備及其催化性能
    射流齒形噴嘴射流流場(chǎng)與氣動(dòng)聲學(xué)分析
    国产一区二区在线观看日韩 | 少妇裸体淫交视频免费看高清| 一进一出抽搐gif免费好疼| 久久人妻av系列| 欧美黑人巨大hd| 国产精品av视频在线免费观看| av福利片在线观看| svipshipincom国产片| 亚洲五月天丁香| av欧美777| 午夜两性在线视频| 黑人欧美特级aaaaaa片| 成人欧美大片| 啦啦啦韩国在线观看视频| xxx96com| 看黄色毛片网站| 国产毛片a区久久久久| 不卡一级毛片| 大型黄色视频在线免费观看| 亚洲av成人不卡在线观看播放网| 国产成人系列免费观看| 婷婷六月久久综合丁香| 精品乱码久久久久久99久播| 欧美日韩中文字幕国产精品一区二区三区| 制服丝袜大香蕉在线| 免费人成视频x8x8入口观看| 男插女下体视频免费在线播放| 久久这里只有精品19| 国产成人一区二区三区免费视频网站| 特大巨黑吊av在线直播| 两个人看的免费小视频| 99国产精品99久久久久| 久久久久国产精品人妻aⅴ院| 日韩中文字幕欧美一区二区| 日韩精品青青久久久久久| 亚洲国产看品久久| 久久中文看片网| 欧美+亚洲+日韩+国产| 一进一出抽搐gif免费好疼| 日韩成人在线观看一区二区三区| 一级a爱片免费观看的视频| 国产av不卡久久| 精品久久久久久久末码| 狂野欧美激情性xxxx| 九色国产91popny在线| 日韩国内少妇激情av| 三级国产精品欧美在线观看 | 午夜激情欧美在线| 国产激情久久老熟女| 国产欧美日韩精品一区二区| 国产av麻豆久久久久久久| 久久久久免费精品人妻一区二区| 在线观看免费视频日本深夜| 国产野战对白在线观看| 热99re8久久精品国产| 色av中文字幕| 中文字幕久久专区| 真人做人爱边吃奶动态| 欧美日韩国产亚洲二区| av天堂中文字幕网| 真人做人爱边吃奶动态| 757午夜福利合集在线观看| 全区人妻精品视频| 一区二区三区高清视频在线| 欧美日本亚洲视频在线播放| 午夜a级毛片| 国产91精品成人一区二区三区| 欧美黑人欧美精品刺激| 中文字幕熟女人妻在线| 亚洲欧美一区二区三区黑人| 99国产精品一区二区蜜桃av| 一本一本综合久久| 午夜激情欧美在线| bbb黄色大片| 国产av一区在线观看免费| 窝窝影院91人妻| 亚洲午夜理论影院| 国产精品野战在线观看| 免费一级毛片在线播放高清视频| 色老头精品视频在线观看| 悠悠久久av| 国产男靠女视频免费网站| 一级a爱片免费观看的视频| 国产真人三级小视频在线观看| 欧美又色又爽又黄视频| 免费看a级黄色片| 一区福利在线观看| 国产欧美日韩一区二区精品| 国产麻豆成人av免费视频| 少妇的丰满在线观看| 中文在线观看免费www的网站| 欧美高清成人免费视频www| 免费av不卡在线播放| 午夜福利高清视频| 丁香欧美五月| 午夜福利免费观看在线| 午夜精品在线福利| 9191精品国产免费久久| 国产精品久久久av美女十八| 亚洲黑人精品在线| 欧美另类亚洲清纯唯美| 亚洲电影在线观看av| 美女黄网站色视频| 这个男人来自地球电影免费观看| 久久久久久久久免费视频了| 成人av一区二区三区在线看| 91av网站免费观看| 久99久视频精品免费| 麻豆一二三区av精品| 好男人电影高清在线观看| 国产亚洲欧美在线一区二区| www国产在线视频色| 日韩国内少妇激情av| 国产一区在线观看成人免费| 丝袜人妻中文字幕| 欧美一级a爱片免费观看看| 国产极品精品免费视频能看的| 黄片小视频在线播放| 国产亚洲av嫩草精品影院| 午夜久久久久精精品| av福利片在线观看| АⅤ资源中文在线天堂| 亚洲成a人片在线一区二区| 亚洲午夜理论影院| 国产蜜桃级精品一区二区三区| 亚洲欧美一区二区三区黑人| 18禁裸乳无遮挡免费网站照片| 日韩成人在线观看一区二区三区| 色噜噜av男人的天堂激情| 天堂√8在线中文| 91久久精品国产一区二区成人 | 久久久久久大精品| 天天躁日日操中文字幕| 女人被狂操c到高潮| 久久精品国产综合久久久| 悠悠久久av| 欧美黑人欧美精品刺激| 久久久精品欧美日韩精品| 在线观看午夜福利视频| 精品无人区乱码1区二区| 亚洲精品456在线播放app | 亚洲国产色片| 午夜福利成人在线免费观看| 日韩大尺度精品在线看网址| 欧美绝顶高潮抽搐喷水| 日本黄大片高清| 最新中文字幕久久久久 | 国产成人欧美在线观看| 啦啦啦韩国在线观看视频| 婷婷精品国产亚洲av| av国产免费在线观看| 国产一区二区激情短视频| 这个男人来自地球电影免费观看| 亚洲av电影不卡..在线观看| 国产午夜精品论理片| 久久久久亚洲av毛片大全| 亚洲18禁久久av| 国产精品自产拍在线观看55亚洲| 黄片大片在线免费观看| 在线看三级毛片| 亚洲国产中文字幕在线视频| 欧美最黄视频在线播放免费| 久久久久国产精品人妻aⅴ院| 婷婷六月久久综合丁香| 免费在线观看成人毛片| 小说图片视频综合网站| 亚洲天堂国产精品一区在线| 又粗又爽又猛毛片免费看| а√天堂www在线а√下载| 国产aⅴ精品一区二区三区波| 91字幕亚洲| 免费观看人在逋| 久久久久精品国产欧美久久久| 一区二区三区高清视频在线| 国产精品 国内视频| 日韩欧美国产在线观看| 亚洲精品一区av在线观看| 午夜亚洲福利在线播放| 99热精品在线国产| 国产高清激情床上av| 午夜精品久久久久久毛片777| 黑人操中国人逼视频| 精品国产三级普通话版| 在线看三级毛片| 亚洲专区国产一区二区| 国产精品女同一区二区软件 | 亚洲熟妇熟女久久| 国产精品永久免费网站| 村上凉子中文字幕在线| 在线观看午夜福利视频| 超碰成人久久| 国产精品九九99| 美女高潮的动态| 好看av亚洲va欧美ⅴa在| 日韩免费av在线播放| 精品国产美女av久久久久小说| 欧美一区二区国产精品久久精品| 天天躁日日操中文字幕| 日韩有码中文字幕| 淫妇啪啪啪对白视频| 亚洲自偷自拍图片 自拍| 久久久久九九精品影院| 黄色女人牲交| 亚洲精华国产精华精| 亚洲国产精品久久男人天堂| 窝窝影院91人妻| 亚洲人成伊人成综合网2020| 欧美成人一区二区免费高清观看 | 美女大奶头视频| 久久天堂一区二区三区四区| 夜夜看夜夜爽夜夜摸| 亚洲自拍偷在线| 午夜a级毛片| 色综合站精品国产| 18禁美女被吸乳视频| 老鸭窝网址在线观看| 成人鲁丝片一二三区免费| 在线播放国产精品三级| 国产精品电影一区二区三区| 国产午夜福利久久久久久| 丁香欧美五月| 亚洲精品美女久久久久99蜜臀| 变态另类丝袜制服| 中文字幕熟女人妻在线| 亚洲七黄色美女视频| 嫁个100分男人电影在线观看| 桃色一区二区三区在线观看| 亚洲人成电影免费在线| 香蕉国产在线看| 国产精品亚洲av一区麻豆| 国产aⅴ精品一区二区三区波| 人妻丰满熟妇av一区二区三区| 天堂网av新在线| 成人高潮视频无遮挡免费网站| 欧美一区二区精品小视频在线| 搡老岳熟女国产| 日韩精品青青久久久久久| 久久中文看片网| 亚洲,欧美精品.| 亚洲精华国产精华精| 国产精品久久久久久人妻精品电影| 国产亚洲精品久久久com| 国产麻豆成人av免费视频| 成年女人毛片免费观看观看9| 日韩欧美国产在线观看| 欧美一级a爱片免费观看看| 手机成人av网站| 不卡一级毛片| 久久精品夜夜夜夜夜久久蜜豆| 日本a在线网址| 亚洲七黄色美女视频| 久久精品91无色码中文字幕| 免费观看的影片在线观看| 99久国产av精品| 欧美一级毛片孕妇| 三级毛片av免费| 亚洲一区高清亚洲精品| 色av中文字幕| 啦啦啦观看免费观看视频高清| 亚洲一区二区三区色噜噜| 精品国产美女av久久久久小说| 叶爱在线成人免费视频播放| 高清毛片免费观看视频网站| 亚洲国产精品久久男人天堂| 亚洲av成人一区二区三| 亚洲成av人片在线播放无| 中文资源天堂在线| 熟妇人妻久久中文字幕3abv| 欧美高清成人免费视频www| 麻豆久久精品国产亚洲av| 国产精品久久久久久亚洲av鲁大| 日本 av在线| 亚洲激情在线av| 免费人成视频x8x8入口观看| 亚洲欧美激情综合另类| 两个人看的免费小视频| 无限看片的www在线观看| 成年版毛片免费区| 一区二区三区高清视频在线| 国产精品免费一区二区三区在线| 极品教师在线免费播放| 1024手机看黄色片| 国产精品自产拍在线观看55亚洲| 给我免费播放毛片高清在线观看| 国产伦人伦偷精品视频| 免费观看的影片在线观看| 国产精品美女特级片免费视频播放器 | 亚洲午夜精品一区,二区,三区| 亚洲电影在线观看av| 精品一区二区三区视频在线观看免费| 男女午夜视频在线观看| 久久欧美精品欧美久久欧美| 国产高清激情床上av| 丰满的人妻完整版| 久久香蕉国产精品| 久久久久九九精品影院| 精品久久久久久成人av| 9191精品国产免费久久| 一区二区三区激情视频| 亚洲美女视频黄频| 91av网一区二区| 午夜福利免费观看在线| 天天躁日日操中文字幕| av视频在线观看入口| 90打野战视频偷拍视频| 国产精品女同一区二区软件 | 久久久久久人人人人人| 男女午夜视频在线观看| 亚洲精品在线观看二区| 黑人操中国人逼视频| 黑人巨大精品欧美一区二区mp4| 国产乱人伦免费视频| 精品国产超薄肉色丝袜足j| 美女免费视频网站| 国内揄拍国产精品人妻在线| 免费av毛片视频| 18禁国产床啪视频网站| 在线免费观看的www视频| 国产高清有码在线观看视频| 香蕉国产在线看| 成人18禁在线播放| 婷婷精品国产亚洲av| 中文字幕人成人乱码亚洲影| 国产成人精品无人区| 国内久久婷婷六月综合欲色啪| 中文字幕久久专区| 国产精品国产高清国产av| 欧美+亚洲+日韩+国产| 日韩欧美 国产精品| 在线永久观看黄色视频| 亚洲欧美一区二区三区黑人| 99热6这里只有精品| 一级毛片女人18水好多| 日本成人三级电影网站| 亚洲片人在线观看| 亚洲欧美一区二区三区黑人| 色视频www国产| 日韩欧美在线乱码| 黄色 视频免费看| 久久中文看片网| 中国美女看黄片| 午夜成年电影在线免费观看| 亚洲国产精品合色在线| 97超视频在线观看视频| 亚洲人成伊人成综合网2020| 12—13女人毛片做爰片一| 少妇裸体淫交视频免费看高清| 国产人伦9x9x在线观看| 久久久国产欧美日韩av| 91麻豆av在线| 亚洲专区中文字幕在线| cao死你这个sao货| 亚洲天堂国产精品一区在线| 大型黄色视频在线免费观看| av在线天堂中文字幕| 看片在线看免费视频| 少妇丰满av| 桃红色精品国产亚洲av| 制服人妻中文乱码| 久久这里只有精品中国| 老司机午夜十八禁免费视频| 亚洲成av人片在线播放无| 欧美午夜高清在线| 欧美色欧美亚洲另类二区| 久久精品人妻少妇| 在线观看午夜福利视频| 久久精品91蜜桃| 亚洲欧美一区二区三区黑人| 免费无遮挡裸体视频| 精品久久久久久久毛片微露脸| 色尼玛亚洲综合影院| 午夜视频精品福利| 欧美日韩福利视频一区二区| 成人特级黄色片久久久久久久| 黄色视频,在线免费观看| 在线观看舔阴道视频| 男女床上黄色一级片免费看| 国产在线精品亚洲第一网站| 一a级毛片在线观看| 99久国产av精品| 大型黄色视频在线免费观看| 国产成人av教育| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 亚洲色图av天堂| 亚洲性夜色夜夜综合| 国内久久婷婷六月综合欲色啪| 亚洲乱码一区二区免费版| 国产亚洲av嫩草精品影院| 搞女人的毛片| 亚洲精品美女久久久久99蜜臀| 亚洲精华国产精华精| 国产 一区 欧美 日韩| 无遮挡黄片免费观看| 1024手机看黄色片| 欧美又色又爽又黄视频| e午夜精品久久久久久久| 色哟哟哟哟哟哟| 成人av在线播放网站| 一个人看的www免费观看视频| 18禁美女被吸乳视频| 小蜜桃在线观看免费完整版高清| 午夜免费观看网址| 国产精品久久视频播放| 熟女电影av网| 国产高清有码在线观看视频| 少妇的逼水好多| www.999成人在线观看| av片东京热男人的天堂| 亚洲中文字幕日韩| 欧美成人一区二区免费高清观看 | 99国产综合亚洲精品| 一级毛片女人18水好多| 久久久久久久午夜电影| 日韩大尺度精品在线看网址| 久久精品国产亚洲av香蕉五月| e午夜精品久久久久久久| 亚洲电影在线观看av| 国产黄色小视频在线观看| 亚洲自偷自拍图片 自拍| 国产一区二区在线观看日韩 | 99久久精品国产亚洲精品| 99视频精品全部免费 在线 | 久久精品亚洲精品国产色婷小说| а√天堂www在线а√下载| 麻豆一二三区av精品| 精品国产美女av久久久久小说| 真实男女啪啪啪动态图| 他把我摸到了高潮在线观看| 亚洲第一欧美日韩一区二区三区| 国内揄拍国产精品人妻在线| 免费看日本二区| 国产极品精品免费视频能看的| 亚洲黑人精品在线| 欧美绝顶高潮抽搐喷水| 黄频高清免费视频| 国产精品 欧美亚洲| 亚洲精品在线美女| 色综合站精品国产| 又爽又黄无遮挡网站| 亚洲 国产 在线| 色精品久久人妻99蜜桃| 757午夜福利合集在线观看| 天天躁日日操中文字幕| 757午夜福利合集在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产高清激情床上av| 欧美中文日本在线观看视频| 国产高清激情床上av| 成年人黄色毛片网站| 亚洲专区中文字幕在线| 桃红色精品国产亚洲av| 亚洲 国产 在线| 亚洲国产色片| www日本在线高清视频| 日本 欧美在线| 欧美激情久久久久久爽电影| 叶爱在线成人免费视频播放| av福利片在线观看| 最近最新中文字幕大全免费视频| 在线a可以看的网站| 欧美zozozo另类| 不卡av一区二区三区| 久久久久久久久久黄片| 久久久国产欧美日韩av| 黄色视频,在线免费观看| 精品一区二区三区四区五区乱码| 9191精品国产免费久久| 精品久久久久久,| cao死你这个sao货| 欧美大码av| 男人舔女人的私密视频| 91av网站免费观看| 一本精品99久久精品77| 国产不卡一卡二| 99在线视频只有这里精品首页| 一本综合久久免费| 国产蜜桃级精品一区二区三区| 亚洲黑人精品在线| 免费观看精品视频网站| 久久久国产欧美日韩av| 免费在线观看成人毛片| 曰老女人黄片| 免费高清视频大片| 国产v大片淫在线免费观看| 我的老师免费观看完整版| 国产一区二区激情短视频| 国产视频内射| 欧美日韩精品网址| 亚洲18禁久久av| 狂野欧美白嫩少妇大欣赏| 国产精品野战在线观看| 国产精品1区2区在线观看.| 午夜免费观看网址| 午夜福利视频1000在线观看| 亚洲 欧美一区二区三区| 亚洲乱码一区二区免费版| 中亚洲国语对白在线视频| 波多野结衣巨乳人妻| 国产三级黄色录像| 美女扒开内裤让男人捅视频| 丁香欧美五月| 真人一进一出gif抽搐免费| 一区二区三区激情视频| 黄频高清免费视频| 无限看片的www在线观看| 美女高潮喷水抽搐中文字幕| 欧美丝袜亚洲另类 | 午夜免费激情av| 精品不卡国产一区二区三区| 热99re8久久精品国产| 成人国产一区最新在线观看| 99热这里只有精品一区 | 欧美性猛交╳xxx乱大交人| 网址你懂的国产日韩在线| 午夜精品在线福利| 我的老师免费观看完整版| 久久久国产成人精品二区| 成人三级黄色视频| 久久热在线av| 99热只有精品国产| 欧美日韩国产亚洲二区| 制服丝袜大香蕉在线| 国产亚洲av高清不卡| 国产成人一区二区三区免费视频网站| 一进一出抽搐gif免费好疼| 亚洲第一欧美日韩一区二区三区| 99re在线观看精品视频| 亚洲 欧美一区二区三区| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| 日日干狠狠操夜夜爽| 国产一区二区三区视频了| 精品一区二区三区视频在线 | 在线国产一区二区在线| 久久天躁狠狠躁夜夜2o2o| 网址你懂的国产日韩在线| 免费观看的影片在线观看| 国内精品美女久久久久久| 成年女人毛片免费观看观看9| 美女午夜性视频免费| 高潮久久久久久久久久久不卡| 天堂动漫精品| 成人高潮视频无遮挡免费网站| 免费人成视频x8x8入口观看| 亚洲 国产 在线| 欧美色欧美亚洲另类二区| 1024香蕉在线观看| netflix在线观看网站| 午夜亚洲福利在线播放| 曰老女人黄片| 99riav亚洲国产免费| 熟女电影av网| av女优亚洲男人天堂 | 两性夫妻黄色片| 精品久久久久久,| 亚洲av成人精品一区久久| 亚洲熟妇熟女久久| 网址你懂的国产日韩在线| 又大又爽又粗| 精品久久久久久,| 最新美女视频免费是黄的| 亚洲精品美女久久av网站| 免费看光身美女| 亚洲国产欧美一区二区综合| 久久欧美精品欧美久久欧美| 日韩欧美免费精品| 午夜两性在线视频| 又爽又黄无遮挡网站| 一区福利在线观看| 波多野结衣巨乳人妻| 哪里可以看免费的av片| 国内精品一区二区在线观看| 极品教师在线免费播放| 欧美日韩黄片免| 热99re8久久精品国产| 1024香蕉在线观看| 啦啦啦韩国在线观看视频| 天天躁日日操中文字幕| 久久精品91蜜桃| 婷婷丁香在线五月| 一进一出好大好爽视频| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 老司机在亚洲福利影院| 精品国产乱码久久久久久男人| 999精品在线视频| 非洲黑人性xxxx精品又粗又长| 免费电影在线观看免费观看| 啦啦啦观看免费观看视频高清| 夜夜爽天天搞| 一区二区三区激情视频| 中出人妻视频一区二区| avwww免费| 91在线观看av| 美女cb高潮喷水在线观看 | 亚洲人成伊人成综合网2020| 最好的美女福利视频网| 午夜精品一区二区三区免费看| 露出奶头的视频| 午夜两性在线视频| 亚洲欧美日韩高清专用| 在线观看一区二区三区| 久久热在线av| 亚洲av第一区精品v没综合| 精品电影一区二区在线| 性欧美人与动物交配| 国产精品九九99| 免费在线观看视频国产中文字幕亚洲| 99精品久久久久人妻精品| 999久久久精品免费观看国产| 99热这里只有精品一区 | 成熟少妇高潮喷水视频| 黄色女人牲交| 亚洲在线自拍视频|