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

    大尺度障礙物與泄爆面對(duì)天然氣內(nèi)爆炸的協(xié)同作用規(guī)律研究*

    2019-01-05 02:13:30呂鵬飛胡倩然
    關(guān)鍵詞:房間內(nèi)障礙物火焰

    楊 凱,呂鵬飛,胡倩然,龐 磊

    (1.北京石油化工學(xué)院 安全工程學(xué)院,北京 102617;2.北京市安全生產(chǎn)工程技術(shù)研究院,北京 102617)

    0 引言

    作為一種高效、清潔能源,天然氣在城市居民生活中得到廣泛應(yīng)用。近年來(lái),由于設(shè)備設(shè)施老化、違規(guī)使用等問(wèn)題誘發(fā)的天然氣爆炸事故已然成為城市居民生活中的典型事故類(lèi)型之一。一般而言,民用建筑內(nèi)天然氣爆炸多始發(fā)在廚房,而廚房的門(mén)窗在爆炸過(guò)程中往往會(huì)成為具有一定開(kāi)啟壓力的泄爆結(jié)構(gòu)。同時(shí),廚房?jī)?nèi)的櫥柜、吊柜等物品很可能加劇爆炸瞬態(tài)流場(chǎng)的復(fù)雜性。氣體爆炸屬于非點(diǎn)源爆炸形式,爆炸過(guò)程受約束條件的影響極為顯著。在典型的室內(nèi)天然氣爆炸事故中,具有一定開(kāi)啟壓力的泄爆面與內(nèi)置大尺度障礙物將共同作用于室內(nèi)的爆炸流場(chǎng)[1]。

    目前,關(guān)于泄爆結(jié)構(gòu)對(duì)氣體爆炸的影響規(guī)律方面已經(jīng)取得了較多的研究成果。Bao,Qi,F(xiàn)akandu,Guo等[2-5]分別借助不同形狀、不同尺度的實(shí)驗(yàn)裝置研究了開(kāi)啟壓力、泄爆面積、相對(duì)位置等泄爆面屬性對(duì)氣體爆炸及泄爆的影響,但由于事故過(guò)程、工藝特點(diǎn)等方面的差異,并沒(méi)有進(jìn)一步考查內(nèi)置障礙物的爆炸效應(yīng)。

    近年來(lái),部分學(xué)者采用實(shí)驗(yàn)手段或數(shù)值方法對(duì)內(nèi)置障礙物作用下氣體爆炸及泄爆效應(yīng)開(kāi)展了大量研究,但民用建筑內(nèi)的家具、電氣設(shè)備等障礙物都是以大尺度為主,有關(guān)大尺度障礙物與泄爆面對(duì)受限空間內(nèi)氣體爆炸協(xié)同的研究還比較缺乏。由于天然氣和液化石油氣是2種主要的居民生活用氣,因此很多文獻(xiàn)[6-15]中主要利用這2種氣體的主要成分甲烷或丙烷氣體作為爆源開(kāi)展研究。這些文獻(xiàn)中泄爆面和障礙物特征是主要考察的2個(gè)參數(shù)。由于研究中爆炸腔室體積尺寸的分布范圍較廣,如Hall等[9]和Tomlin等[14]分別采用體積為0.000 625 m3和182.25 m3的爆炸腔室,研究了液化石油氣和天然氣的約束爆炸過(guò)程,使得各文獻(xiàn)中泄爆面積和障礙物體積的差別較大,而約束泄爆過(guò)程具有尺度效應(yīng),小尺度泄爆實(shí)驗(yàn)結(jié)果不能直接用于指導(dǎo)民用建筑研究[1]?;跓o(wú)量綱分析發(fā)現(xiàn)文獻(xiàn)[6-15]中泄壓比和障礙物面積阻塞率的分布范圍較寬,分別為0.1~1和0%~70%。雖然這些文獻(xiàn)考察了各種泄壓比對(duì)泄爆過(guò)程的影響,但泄爆面的開(kāi)啟壓力相對(duì)較低,均在數(shù)kPa范圍內(nèi),而民用建筑房間玻璃的開(kāi)啟壓力通常在數(shù)kPa到幾十kPa之間。當(dāng)房間泄壓比在特定范圍時(shí),泄爆面開(kāi)啟壓力增加將加劇室內(nèi)爆炸超壓強(qiáng)度。盡管Bao等[2]借助大尺度實(shí)驗(yàn)裝置研究了泄爆面靜開(kāi)啟壓力對(duì)室內(nèi)爆炸超壓的影響,但并沒(méi)有考慮室內(nèi)障礙物的協(xié)同作用。障礙物的存在會(huì)增加室內(nèi)湍流源,從而促進(jìn)火焰速度和爆炸超壓。雖然文獻(xiàn) [6-15]考察了圓柱、立方體、長(zhǎng)方體等不同形狀障礙物對(duì)約束泄爆過(guò)程的影響,但這些障礙物通常以小尺度細(xì)長(zhǎng)形為主,并采用連續(xù)均勻方式布置,而民用建筑內(nèi)的家具、電氣設(shè)備等典型物品,特別是廚房?jī)?nèi)的障礙物的幾何尺度通常較大[14],障礙物布置較為緊湊。由此可見(jiàn),目前同時(shí)涉及障礙物與泄爆面的有關(guān)研究主要針對(duì)小尺度連續(xù)障礙物,且泄爆面開(kāi)啟壓力相對(duì)較低的情況,尚未關(guān)注房間內(nèi)的大尺度障礙物與約束泄爆面之間的協(xié)同作用,而針對(duì)小尺度障礙物的研究對(duì)民用建筑天然氣爆炸災(zāi)害研究的參考意義較為有限。

    基于此,以典型民用廚房?jī)?nèi)天然氣爆炸為例,研究了約束泄爆面與內(nèi)置大尺度障礙物對(duì)爆炸動(dòng)力學(xué)演化過(guò)程的協(xié)同作用機(jī)制,分析了泄爆面靜開(kāi)啟壓力和大尺度障礙物體積阻塞率對(duì)約束泄爆效應(yīng)的影響,力求為此類(lèi)事故的致因機(jī)理及事故防治提供科學(xué)依據(jù)。

    1 數(shù)值方法

    1.1 數(shù)值模型

    基于計(jì)算流體動(dòng)力學(xué)技術(shù)的氣體爆炸動(dòng)力學(xué)仿真軟件AutoReaGas適用于氣體爆炸及沖擊動(dòng)力學(xué)問(wèn)題,且該軟件的PRESSURFS邊界條件可以根據(jù)實(shí)際情況對(duì)特定面結(jié)構(gòu)設(shè)置特定的靜開(kāi)啟壓力,從而實(shí)現(xiàn)受限空間約束泄爆瞬態(tài)流場(chǎng)的求解,且該軟件在求解氣體爆炸泄壓流場(chǎng)的可靠性已經(jīng)在文獻(xiàn)[16-17]中得到了驗(yàn)證。因此,本文采用該軟件開(kāi)展受限空間內(nèi)預(yù)混氣體約束爆炸研究。軟件中假設(shè)可燃?xì)怏w爆炸為單步不可逆化學(xué)反應(yīng),并采用有限體積法數(shù)值對(duì)包括質(zhì)量守恒方程、動(dòng)量守恒方程和能量守恒方程等一系列方程組進(jìn)行求解來(lái)解決氣體爆炸過(guò)程中的氣體動(dòng)力學(xué)行為[18];通過(guò)k-ε模型來(lái)模擬氣體爆炸中的湍流行為。

    層流燃燒時(shí)的氣體燃燒率Rfu=C1ρRmin,其中,C1為層流燃燒率模型常數(shù),取系統(tǒng)默認(rèn)值;ρ為密度,kg/m3;Rmin為燃料質(zhì)量分?jǐn)?shù)、氧氣質(zhì)量分?jǐn)?shù)以及產(chǎn)物質(zhì)量分?jǐn)?shù)的最小值。

    層流燃燒速度Sb可表示為:

    Sb=Sl(1+FSRf)

    (1)

    式中:Sl為特定層流燃燒速度,m/s;Rf為球形火焰半徑,m;FS為層流火焰加速系數(shù),取值為0.15。

    湍流燃燒時(shí)的氣體燃燒率Rfu=Ctρ(St2/Γfu)Rmin;Ct為火焰速度常數(shù),取值為40;ρ為密度,kg/m3;Γfu為湍流耗散系數(shù);Rmin為燃料質(zhì)量分?jǐn)?shù)、氧氣質(zhì)量分?jǐn)?shù)以及產(chǎn)物質(zhì)量分?jǐn)?shù)的最小值。

    湍流燃燒速度St可以表示為:

    St=1.8ut0.412·Lt0.196·Sl0.784·ν-0.196

    (2)

    式中:ut為湍流強(qiáng)度;Lt為特征長(zhǎng)度尺寸,m;Sl為特定層流燃燒速度,m/s;v為未燃混合物的運(yùn)動(dòng)粘度,m2/s。

    1.2 網(wǎng)格獨(dú)立性驗(yàn)證

    AutoReaGas采用8節(jié)點(diǎn)結(jié)構(gòu)網(wǎng)格對(duì)空間離散,采用子網(wǎng)格技術(shù)描述小尺度實(shí)體對(duì)流場(chǎng)的作用[19]。為了考察網(wǎng)格尺寸對(duì)爆炸流場(chǎng)的精確捕捉的影響。選取尺寸為5 m×2 m×2.5 m房間作為研究對(duì)象。房間內(nèi)無(wú)障礙物。在房間較小的一面豎直墻體上設(shè)有1個(gè)尺寸為0.5 m×0.5 m的泄爆面,其開(kāi)啟壓力為60 kPa。在房間內(nèi)充滿濃度為9.5%且處于靜止?fàn)顟B(tài)的甲烷/空氣預(yù)混氣體。點(diǎn)火源位于房間后壁中心處。針對(duì)該物理模型,分別采用0.1 m×0.1 m×0.1 m(M1)和0.05 m×0.05 m×0.05 m(M2)的2套尺寸對(duì)其進(jìn)行網(wǎng)格劃分,并對(duì)比分析了各算例的峰值超壓。表1給出了2套網(wǎng)格的峰值超壓對(duì)比數(shù)據(jù)。從表1中可以看出峰值超壓的相對(duì)誤差在10%以?xún)?nèi),表明2種網(wǎng)格尺寸對(duì)本計(jì)算的影響較小,均可應(yīng)用本計(jì)算,為了節(jié)約計(jì)算時(shí)間,采用尺寸為0.1 m的網(wǎng)格開(kāi)展研究。

    表1 峰值超壓對(duì)比Table 1 Comparison of peak overpressure

    1.3 數(shù)值模型的有效性驗(yàn)證

    為了驗(yàn)證本文中數(shù)值模型的有效性,借助該數(shù)值模型對(duì)與本研究工況相類(lèi)似的Bao等[2]開(kāi)展的大尺度天然氣爆炸實(shí)驗(yàn)過(guò)程進(jìn)行數(shù)值模擬和對(duì)比分析。Bao等[2]采用2 m×2 m×3 m的大尺度實(shí)驗(yàn)裝置開(kāi)展了天然氣爆炸與泄爆實(shí)驗(yàn),并分析了6組不同靜開(kāi)啟壓力的泄爆面以及9組不同體積濃度的甲烷氣體的泄爆效應(yīng),泄爆面面積為0.64 m2,點(diǎn)火源位于爆室中心位置,超壓曲線取自爆室后壁幾何中心位置的測(cè)點(diǎn),爆室簡(jiǎn)化模型如圖1所示。

    圖1 對(duì)比實(shí)驗(yàn)中[2]爆炸室簡(jiǎn)化模型Fig.1 Simplified model of explosion chamber in contrastive experiment[2]

    在對(duì)應(yīng)的數(shù)值模擬中,參照實(shí)驗(yàn)參數(shù)設(shè)置泄爆結(jié)構(gòu)的開(kāi)啟壓力和泄爆面積,其他墻壁均設(shè)置為剛性固壁邊界。采用甲烷/空氣預(yù)混氣體作為爆源,其中甲烷體積濃度為9.5%。參照實(shí)驗(yàn)中的測(cè)試布置,將數(shù)值模擬中的測(cè)點(diǎn)布置在房間后壁幾何中心,距離后壁0.1 m。點(diǎn)火源位于房間中心位置,點(diǎn)火時(shí)環(huán)境初始?jí)毫统跏紲囟确謩e設(shè)定為1.013 25×105Pa和300 K?;趯?shí)驗(yàn)工況,對(duì)泄爆面開(kāi)啟壓力為1.5 ,12.1, 20.6, 31.4和62.7 kPa的5組泄爆工況分別進(jìn)行了數(shù)值模擬。不同開(kāi)啟壓力下數(shù)值模擬與實(shí)驗(yàn)所得的室內(nèi)峰值超壓如圖2所示。

    圖2 實(shí)驗(yàn)與數(shù)值模擬結(jié)果中爆炸超壓的對(duì)比Fig.2 Comparison of peak overpressure between numerical and experimental results

    對(duì)比發(fā)現(xiàn),實(shí)驗(yàn)和數(shù)值模擬的最大絕對(duì)誤差為3.4 kPa,最大相對(duì)誤差為4.2%。由于實(shí)驗(yàn)裝置自身、傳感器精度以及實(shí)驗(yàn)過(guò)程存在一定程度的不確定性,因此,其中的誤差可能由實(shí)驗(yàn)測(cè)試誤差與計(jì)算域離散等多方面原因綜合所致。實(shí)驗(yàn)與數(shù)值模擬的對(duì)比證實(shí),借助以上數(shù)值方法求解室內(nèi)天然氣爆炸與泄爆瞬態(tài)流場(chǎng)是可行的。

    2 研究方案

    由于民用建筑內(nèi)天然氣爆炸通常始發(fā)于廚房?jī)?nèi),因此基于中國(guó)民用建筑廚房典型幾何尺寸以及內(nèi)部裝修布置特征,選擇房間模型尺寸為5 m(長(zhǎng))×2 m(寬)×2.5 m(高)。模型四周及屋頂和地板均為剛性壁面,其中較小面積的一面墻體上設(shè)置1個(gè)尺寸為0.5 m×0.5 m的具有一定開(kāi)啟壓力的泄爆面。為準(zhǔn)確描述泄爆口外部二次爆炸對(duì)室內(nèi)爆炸災(zāi)害的影響,將計(jì)算域延泄爆口方向向外延伸1倍房間長(zhǎng)度,并將延伸的該部分計(jì)算域設(shè)置為自由出流邊界。采用尺寸為0.1 m的正方體結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算域進(jìn)行離散。點(diǎn)火源位于房間后壁墻體的幾何中心位置。由于天然氣的主要成分為甲烷,在房間內(nèi)均勻填充9.5%濃度且靜止?fàn)顟B(tài)的甲烷/空氣預(yù)混氣體,計(jì)算域內(nèi)的環(huán)境初始?jí)毫统跏紲囟确謩e設(shè)定為1.013 25×105Pa和300 K。

    根據(jù)民用建筑窗體玻璃類(lèi)型的不同,其靜開(kāi)啟壓力在7~60 kPa[2],為研究房間內(nèi)窗體不同靜開(kāi)啟壓力對(duì)室內(nèi)天然氣約束泄爆效應(yīng)的影響,在研究中改變泄爆面的靜開(kāi)啟壓力,分別設(shè)為20,40和60 kPa。同時(shí)為研究房間內(nèi)大尺度障礙物對(duì)天然氣約束泄爆的影響,在房間內(nèi)分別設(shè)置了體積阻塞率為0,0.15和0.30的障礙物。民用建筑廚房?jī)?nèi)的障礙物通常布置在房間四周,根據(jù)中國(guó)民用廚房結(jié)構(gòu)特征,可以分為“I型”、“Ⅱ型”、“U型”、“L型”4類(lèi),本文采用“Ⅱ型”的廚房布置形式,即沿著房間底部和頂部?jī)蓚?cè)長(zhǎng)度方向設(shè)置障礙物。房間內(nèi)障礙物尺寸參照中國(guó)典型廚房中的櫥柜、吊柜的尺寸和形式,位于下部的櫥柜和位于上部的吊柜橫截面尺寸分別為0.85 m(高)×0.5 m(寬)和0.65 m(高)×0.5 m(寬),櫥柜和吊柜的長(zhǎng)度均與房間長(zhǎng)度一致,此時(shí)房間體積阻塞率為0.3,而采用將櫥柜和吊柜的高度降低一半的方法來(lái)實(shí)現(xiàn)阻塞率為0.15的實(shí)驗(yàn)條件。根據(jù)房間長(zhǎng)度,在模型中設(shè)置了8個(gè)測(cè)點(diǎn),所有測(cè)點(diǎn)均位于房間中心軸線上,除了點(diǎn)火源附近和泄爆面兩側(cè)的測(cè)點(diǎn)外,其余測(cè)點(diǎn)均按照1 m等間距布置。廚房結(jié)構(gòu)以及測(cè)點(diǎn)布置如圖3所示。

    圖3 典型廚房及內(nèi)置大尺度障礙物的物理模型Fig.3 Physical model oftypical kitchen with built-in large scale obstacles

    3 結(jié)果與討論

    3.1 無(wú)障礙物時(shí)不同開(kāi)啟壓力對(duì)應(yīng)的室內(nèi)爆炸過(guò)程分析

    圖4給出了空房間內(nèi)不同開(kāi)啟壓力對(duì)應(yīng)的室內(nèi)軸線上的火焰速度。顯然,各組火焰速度整體呈增長(zhǎng)趨勢(shì),但火焰速度增長(zhǎng)相對(duì)較慢,且漲幅較小。同時(shí),最大火焰速度隨著泄爆面開(kāi)啟壓力的增加而不斷增大,表明開(kāi)啟壓力對(duì)火焰速度存在影響。這是由于開(kāi)啟壓力越大,泄爆時(shí)刻房間內(nèi)和房間外的壓差也越大,泄爆引發(fā)的氣體出流速度更快,進(jìn)而沿泄爆口的火焰速度增加最為顯著。

    圖4 無(wú)障礙物房間內(nèi)不同開(kāi)啟壓力對(duì)應(yīng)的火焰速度-距離Fig.4 Flame speed-distance in unobstructed rooms with different activation overpressures

    圖5給出了空房間內(nèi)不同開(kāi)啟壓力對(duì)應(yīng)的房間中心位置的爆炸超壓隨時(shí)間的分布,可見(jiàn),每條曲線均呈現(xiàn)2個(gè)明顯的超壓峰值,超壓峰值p1是由房間泄爆面開(kāi)啟造成的,而超壓峰值p2則是由室內(nèi)外壓力平衡以及外泄氣體二次爆炸等綜合作用所致。2個(gè)峰值均隨開(kāi)啟壓力的增加而顯著增大。隨著泄爆面開(kāi)啟壓力增加,室內(nèi)氣體的出流速度和湍流強(qiáng)度進(jìn)一步增加,火焰速度進(jìn)一步加快,進(jìn)而氣體燃燒速度和燃燒產(chǎn)物也隨之增加,泄爆面開(kāi)啟后房間內(nèi)外壓力平衡需求和外部二次爆炸對(duì)室內(nèi)流場(chǎng)的作用也更為強(qiáng)烈,最終導(dǎo)致室內(nèi)更大的超壓峰值p2。

    圖5 無(wú)障礙物房間內(nèi)不同開(kāi)啟壓力對(duì)應(yīng)的超壓-時(shí)間Fig.5 Overpressure-time at centric position in unobstructed rooms with different activation overpressures

    3.2 有障礙物時(shí)不同開(kāi)啟壓力對(duì)應(yīng)的室內(nèi)爆炸過(guò)程分析

    圖6(a)和圖6(b)分別為體積阻塞率為0.15和0.3對(duì)應(yīng)的室內(nèi)火焰速度分布。顯然,各組火焰速度總體上呈現(xiàn)增長(zhǎng)趨勢(shì),且沿傳播距離均可以分為3個(gè)類(lèi)似階段。第1階段是火焰緩慢增長(zhǎng)階段,出現(xiàn)在距房間后壁2 m范圍內(nèi),該范圍內(nèi)火焰速度小于10 m/s。第2階段是火焰快速增長(zhǎng)階段,出現(xiàn)在距房間后壁2 m至4 m范圍內(nèi),該范圍內(nèi)火焰速度介于10 m/s與70 m/s之間。第3階段是火焰波動(dòng)階段,該范圍內(nèi)火焰速度產(chǎn)生波動(dòng),但總體上處于較大的速度值。出現(xiàn)波動(dòng)是由于該區(qū)域靠近泄爆口,泄爆面開(kāi)啟作用于氣體流動(dòng),造成了火焰速度的不穩(wěn)定性。

    對(duì)比圖4和圖6發(fā)現(xiàn),相對(duì)于開(kāi)啟壓力,障礙物對(duì)促進(jìn)火焰速度增長(zhǎng)的影響更為顯著。由于障礙物是普通房間內(nèi)的主要湍流源,對(duì)爆炸災(zāi)害的傳播過(guò)程主要表現(xiàn)為剪切作用,誘發(fā)火焰鋒面前端的未燃?xì)怏w形成湍流,從而進(jìn)入火焰速度增加的正反饋機(jī)制。當(dāng)廚房?jī)?nèi)的大尺度障礙物采用“Ⅱ型”布置時(shí),由于障礙物沿房間長(zhǎng)度方向布置,并沒(méi)有改變火焰?zhèn)鞑サ目傞L(zhǎng)度。因此,障礙物體積阻塞率的增加雖然導(dǎo)致房間內(nèi)可燃?xì)怏w的體積減小,但同時(shí)也促使實(shí)際爆炸有效空間的長(zhǎng)徑比增加,造成氣體流速相對(duì)增加,最終導(dǎo)致火焰速度更早地進(jìn)入快速增長(zhǎng)階段。

    圖6 有障礙物時(shí)不同開(kāi)啟壓力對(duì)應(yīng)的火焰速度-距離Fig.6 Flame speed-distance in obstructed rooms with different activation overpressures

    圖7(a)和圖7(b)分別為體積阻塞率0.15和0.3時(shí)不同開(kāi)啟壓力對(duì)應(yīng)的超壓時(shí)間曲線。當(dāng)阻塞率為0.15時(shí),超壓峰值p1和p2均隨開(kāi)啟壓力的增加而增大,峰值到達(dá)時(shí)間也隨之增加。當(dāng)阻塞率為0.3時(shí),峰值p1越發(fā)不明顯。阻塞率的增加導(dǎo)致火焰速度增加更為迅速,促使氣體生成速率明顯大于氣體泄放速率,房間內(nèi)壓力在極短時(shí)間內(nèi)到達(dá)峰值,而泄爆面開(kāi)啟對(duì)室內(nèi)超壓的影響隨著測(cè)點(diǎn)遠(yuǎn)離泄爆口而發(fā)生延遲,進(jìn)而p1還未及時(shí)形成便進(jìn)入火焰快速增長(zhǎng)階段??傮w而言,隨著體積阻塞率增加,泄爆面開(kāi)啟壓力對(duì)峰值超壓p1的影響逐漸降低。

    圖7 有障礙物時(shí)不同開(kāi)啟壓力對(duì)應(yīng)的超壓-時(shí)間Fig.7 Overpressure-time at centric position in obstructed rooms with different activation overpressures

    3.3 泄爆面與大尺度障礙物的協(xié)同效應(yīng)分析

    圖8為房間內(nèi)不同阻塞率時(shí)最大火焰速度隨泄爆面開(kāi)啟壓力的變化。顯然,最大火焰速度出現(xiàn)于最大開(kāi)啟壓力和最大體積阻塞率同時(shí)存在時(shí),但相比于泄爆面開(kāi)啟壓力,障礙物體積阻塞率對(duì)火焰速度的影響更為顯著。從圖中可見(jiàn),當(dāng)室內(nèi)阻塞率為0.15時(shí)最大火焰速度增長(zhǎng)最快。可能是由于當(dāng)房間內(nèi)阻塞率為0時(shí),室內(nèi)湍流源較少,火焰速度相對(duì)較低,隨開(kāi)啟壓力的增長(zhǎng)也較緩慢;而當(dāng)阻塞率超過(guò)0.3時(shí),大體積障礙物的存在一方面會(huì)導(dǎo)致室內(nèi)可燃?xì)怏w量減少,另一方面也會(huì)對(duì)火焰加速傳播產(chǎn)生一定的阻礙作用。

    圖8 不同體積阻塞率對(duì)應(yīng)的最大火焰速度-開(kāi)啟壓力曲線Fig.8 Peak flame speed-activation overpressure with different volume blockages

    圖9 不同體積阻塞率對(duì)應(yīng)的峰值超壓-開(kāi)啟壓力曲線Fig.9 Peak overpressure-activation overpressure with different volume blockages

    圖9為房間內(nèi)不同阻塞率時(shí)峰值超壓隨泄爆面開(kāi)啟壓力的變化趨勢(shì)??梢?jiàn),峰值超壓出現(xiàn)于最大開(kāi)啟壓力和最大體積阻塞率同時(shí)存在時(shí),即阻塞率與開(kāi)啟壓力對(duì)室內(nèi)峰值超壓的增長(zhǎng)產(chǎn)生了協(xié)同作用,共同促進(jìn)其增長(zhǎng),但相比于泄爆面開(kāi)啟壓力,障礙物體積阻塞率對(duì)超壓發(fā)展的影響更為顯著。

    圖10為不同體積阻塞率時(shí)室內(nèi)峰值超壓到達(dá)時(shí)間隨開(kāi)啟壓力的變化曲線。

    圖10 不同體積阻塞率對(duì)應(yīng)的峰值超壓到達(dá)時(shí)間-開(kāi)啟壓力曲線Fig.10 Arrival time of peak overpressure-activation overpressure with different volume blockages

    由圖10可見(jiàn),峰值超壓到達(dá)時(shí)間現(xiàn)于最大開(kāi)啟壓力和最大體積阻塞率同時(shí)存在時(shí),這表明大尺度障礙物與開(kāi)啟壓力對(duì)峰值超壓到達(dá)時(shí)間也具有協(xié)同作用,共同促進(jìn)超壓快速發(fā)展,縮短峰值超壓到達(dá)時(shí)間,最短的峰值超壓到達(dá)時(shí)間出現(xiàn)于最大開(kāi)啟壓力和最大體積阻塞率同時(shí)存在時(shí)。同時(shí),相比開(kāi)啟壓力,障礙物體積阻塞率對(duì)峰值超壓到達(dá)時(shí)間的影響也更顯著。

    4 結(jié)論

    1)大尺度障礙物與泄爆面對(duì)室內(nèi)天然氣爆炸過(guò)程具有顯著的協(xié)同作用,共同促進(jìn)火焰速度與爆炸超壓的顯著增長(zhǎng),并縮短峰值超壓到達(dá)時(shí)間。在室內(nèi)天然氣爆炸事故調(diào)查分析過(guò)程中,應(yīng)對(duì)這一協(xié)同作用給予足夠重視,孤立分析單一因素容易低估爆炸強(qiáng)度。

    2)相比泄爆面開(kāi)啟壓力,大尺度障礙物體積阻塞率對(duì)泄爆效應(yīng)的影響更為明顯。盡管障礙物的存在顯著減少了參與爆炸的天然氣量,但從增加房間內(nèi)湍流源和房間相對(duì)長(zhǎng)徑比的角度來(lái)衡量,大尺度障礙物對(duì)爆炸過(guò)程具有加劇作用。而隨著開(kāi)啟壓力的增加,障礙物對(duì)泄爆效應(yīng)的影響更為強(qiáng)烈。在民用住宅天然氣爆炸事故調(diào)查時(shí),應(yīng)對(duì)室內(nèi)家具、電器等大尺度物品所形成的體積阻塞給予更多關(guān)注。

    3)火焰峰面是天然氣爆炸事故的主要致災(zāi)因素之一。大尺度障礙物與泄爆面協(xié)同作用下,室內(nèi)火焰速度沿泄爆方向呈現(xiàn)出緩慢增長(zhǎng)、快速增長(zhǎng)、波動(dòng)3個(gè)階段。這一特征加劇了室內(nèi)人員燒傷的復(fù)雜性和嚴(yán)重性,在災(zāi)害后果評(píng)估中應(yīng)予以考慮。

    猜你喜歡
    房間內(nèi)障礙物火焰
    《火焰》
    最亮的火焰
    越南·河內(nèi)
    給小狗的新毛衣
    高低翻越
    SelTrac?CBTC系統(tǒng)中非通信障礙物的設(shè)計(jì)和處理
    漂在水上的火焰
    穿過(guò)下午
    詩(shī)潮(2018年4期)2018-04-26 05:27:58
    智趣
    讀者(2017年10期)2017-04-26 18:33:14
    吹不滅的火焰
    學(xué)與玩(2017年6期)2017-02-16 07:07:22
    黑人巨大精品欧美一区二区mp4| 国产片内射在线| 午夜激情久久久久久久| 精品国产乱子伦一区二区三区 | 亚洲欧洲精品一区二区精品久久久| 国产成人av教育| 欧美精品人与动牲交sv欧美| av超薄肉色丝袜交足视频| 免费高清在线观看视频在线观看| 麻豆av在线久日| 亚洲精品自拍成人| 亚洲天堂av无毛| 老司机午夜十八禁免费视频| 精品国产乱码久久久久久小说| 欧美精品一区二区免费开放| 日韩有码中文字幕| 欧美黄色淫秽网站| 人妻久久中文字幕网| 高清在线国产一区| 日日爽夜夜爽网站| 一本大道久久a久久精品| 美女国产高潮福利片在线看| 欧美日韩精品网址| 欧美日韩黄片免| av线在线观看网站| 久久国产精品男人的天堂亚洲| 美女午夜性视频免费| 国产xxxxx性猛交| 欧美97在线视频| 午夜免费观看性视频| 人人妻人人爽人人添夜夜欢视频| 成人亚洲精品一区在线观看| 亚洲av日韩在线播放| 91国产中文字幕| 日韩制服骚丝袜av| 欧美中文综合在线视频| 国产又爽黄色视频| 一级a爱视频在线免费观看| 超碰成人久久| 欧美日韩精品网址| 成人av一区二区三区在线看 | 又大又爽又粗| 麻豆国产av国片精品| 久久久精品国产亚洲av高清涩受| 爱豆传媒免费全集在线观看| 欧美少妇被猛烈插入视频| 亚洲精品久久午夜乱码| 男男h啪啪无遮挡| 我的亚洲天堂| 法律面前人人平等表现在哪些方面 | 十八禁人妻一区二区| 美女大奶头黄色视频| 又大又爽又粗| 丝袜人妻中文字幕| 美女主播在线视频| 日本wwww免费看| 精品久久久久久电影网| 一本色道久久久久久精品综合| 一进一出抽搐动态| 国产97色在线日韩免费| 99久久精品国产亚洲精品| 精品久久久久久电影网| 最新在线观看一区二区三区| 亚洲第一欧美日韩一区二区三区 | 久久国产精品大桥未久av| 亚洲av日韩在线播放| 最黄视频免费看| 国产激情久久老熟女| 在线观看免费高清a一片| 777久久人妻少妇嫩草av网站| 午夜精品国产一区二区电影| 丁香六月天网| 国产成人系列免费观看| avwww免费| 热99久久久久精品小说推荐| 亚洲精品美女久久av网站| 欧美日韩亚洲综合一区二区三区_| 啪啪无遮挡十八禁网站| 操出白浆在线播放| 熟女少妇亚洲综合色aaa.| 高潮久久久久久久久久久不卡| 国产成+人综合+亚洲专区| 午夜两性在线视频| 乱人伦中国视频| 欧美变态另类bdsm刘玥| 国产一区二区三区在线臀色熟女 | 欧美一级毛片孕妇| 亚洲欧美日韩高清在线视频 | 国产精品久久久久成人av| 久久精品亚洲熟妇少妇任你| 777久久人妻少妇嫩草av网站| 性色av一级| 精品少妇一区二区三区视频日本电影| 亚洲精品一卡2卡三卡4卡5卡 | 日本欧美视频一区| 少妇精品久久久久久久| 国产淫语在线视频| 另类亚洲欧美激情| 日本精品一区二区三区蜜桃| 日韩三级视频一区二区三区| 免费高清在线观看日韩| 狠狠婷婷综合久久久久久88av| 久久精品亚洲av国产电影网| 日韩人妻精品一区2区三区| 男女之事视频高清在线观看| 久久免费观看电影| 国产黄频视频在线观看| 成人手机av| 国产熟女午夜一区二区三区| 亚洲国产av新网站| 色综合欧美亚洲国产小说| 国产精品免费大片| 久久午夜综合久久蜜桃| 亚洲av成人一区二区三| 成年人午夜在线观看视频| 欧美日韩av久久| 建设人人有责人人尽责人人享有的| 亚洲全国av大片| 一区福利在线观看| 99精品欧美一区二区三区四区| 免费观看人在逋| 99精品久久久久人妻精品| 国产成人啪精品午夜网站| 久久青草综合色| 免费在线观看视频国产中文字幕亚洲 | 国产精品久久久久成人av| 久热爱精品视频在线9| 亚洲全国av大片| 一区福利在线观看| 国产欧美日韩精品亚洲av| 美女脱内裤让男人舔精品视频| 一级,二级,三级黄色视频| 亚洲一区二区三区欧美精品| 最黄视频免费看| 91国产中文字幕| 一进一出抽搐动态| 老司机福利观看| 两个人看的免费小视频| 久久热在线av| 搡老熟女国产l中国老女人| 老熟妇乱子伦视频在线观看 | 精品一区二区三卡| 成人亚洲精品一区在线观看| 高清黄色对白视频在线免费看| 精品国产一区二区久久| 国产日韩一区二区三区精品不卡| 老司机深夜福利视频在线观看 | 亚洲中文av在线| 性色av乱码一区二区三区2| 91国产中文字幕| 日韩大片免费观看网站| 日本猛色少妇xxxxx猛交久久| 高潮久久久久久久久久久不卡| 搡老熟女国产l中国老女人| 久久国产精品影院| av在线app专区| 久久国产精品男人的天堂亚洲| 热99久久久久精品小说推荐| 久久精品人人爽人人爽视色| 高清黄色对白视频在线免费看| 后天国语完整版免费观看| e午夜精品久久久久久久| 成人影院久久| 人人妻人人添人人爽欧美一区卜| 欧美乱码精品一区二区三区| 亚洲成人免费电影在线观看| 无遮挡黄片免费观看| kizo精华| 亚洲九九香蕉| 青草久久国产| 老司机午夜福利在线观看视频 | 亚洲欧美一区二区三区黑人| 欧美成狂野欧美在线观看| 国产区一区二久久| 丰满迷人的少妇在线观看| 中文字幕精品免费在线观看视频| 国产日韩欧美视频二区| 欧美+亚洲+日韩+国产| 啦啦啦免费观看视频1| av有码第一页| 如日韩欧美国产精品一区二区三区| av一本久久久久| 久久久久网色| 国产野战对白在线观看| 国产一区二区激情短视频 | 亚洲欧美日韩高清在线视频 | 又黄又粗又硬又大视频| 久久久久国产精品人妻一区二区| 波多野结衣av一区二区av| 久久亚洲国产成人精品v| 十八禁网站网址无遮挡| 高清黄色对白视频在线免费看| 中文欧美无线码| 欧美一级毛片孕妇| 日日夜夜操网爽| 999精品在线视频| 久久精品国产a三级三级三级| 在线精品无人区一区二区三| 久久人妻熟女aⅴ| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美久久黑人一区二区| 老熟妇乱子伦视频在线观看 | 国产男女超爽视频在线观看| 99热国产这里只有精品6| 精品国产乱子伦一区二区三区 | 中文欧美无线码| 多毛熟女@视频| 手机成人av网站| 亚洲国产毛片av蜜桃av| 午夜福利免费观看在线| 国产成人精品久久二区二区91| 日韩有码中文字幕| 高潮久久久久久久久久久不卡| 性色av一级| 在线精品无人区一区二区三| 香蕉国产在线看| 久久热在线av| 国产精品国产三级国产专区5o| 久久毛片免费看一区二区三区| 99精国产麻豆久久婷婷| 三级毛片av免费| 真人做人爱边吃奶动态| 每晚都被弄得嗷嗷叫到高潮| 99re6热这里在线精品视频| 久久九九热精品免费| 啦啦啦中文免费视频观看日本| 每晚都被弄得嗷嗷叫到高潮| 高清视频免费观看一区二区| 一区二区三区精品91| 成人影院久久| 侵犯人妻中文字幕一二三四区| 成年动漫av网址| 在线看a的网站| 中文字幕色久视频| 超碰成人久久| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看日韩| 人人妻人人澡人人看| 日日夜夜操网爽| 午夜精品国产一区二区电影| 国产三级黄色录像| 亚洲人成电影免费在线| 日日夜夜操网爽| 亚洲精品国产av蜜桃| 成年人免费黄色播放视频| 欧美精品一区二区免费开放| 欧美少妇被猛烈插入视频| 久久国产精品影院| 男人舔女人的私密视频| 一本一本久久a久久精品综合妖精| 精品国产乱码久久久久久男人| 久久九九热精品免费| 欧美少妇被猛烈插入视频| 亚洲专区字幕在线| 国产精品久久久久久精品古装| 91麻豆精品激情在线观看国产 | 久久久久视频综合| 亚洲欧美清纯卡通| av在线老鸭窝| 日本a在线网址| 亚洲欧洲日产国产| 亚洲欧美精品综合一区二区三区| 欧美午夜高清在线| 伊人亚洲综合成人网| 丰满人妻熟妇乱又伦精品不卡| 曰老女人黄片| 国产精品久久久人人做人人爽| 嫁个100分男人电影在线观看| 久久综合国产亚洲精品| 性色av一级| 黄色视频不卡| 日韩制服骚丝袜av| 欧美人与性动交α欧美软件| 午夜福利乱码中文字幕| 18在线观看网站| 亚洲欧美日韩高清在线视频 | 97在线人人人人妻| 99国产综合亚洲精品| av片东京热男人的天堂| 亚洲精品国产av成人精品| 精品亚洲成a人片在线观看| 免费观看人在逋| 成人国产av品久久久| 久9热在线精品视频| 黑人欧美特级aaaaaa片| 动漫黄色视频在线观看| 亚洲精品国产精品久久久不卡| 日本wwww免费看| 一本大道久久a久久精品| 日韩电影二区| 国产精品成人在线| 91精品国产国语对白视频| 极品少妇高潮喷水抽搐| 性色av乱码一区二区三区2| 色婷婷av一区二区三区视频| 一本大道久久a久久精品| 国产淫语在线视频| 1024视频免费在线观看| 精品国内亚洲2022精品成人 | 99精品久久久久人妻精品| 女人爽到高潮嗷嗷叫在线视频| 成年动漫av网址| 香蕉国产在线看| 9191精品国产免费久久| 波多野结衣av一区二区av| 男女边摸边吃奶| 看免费av毛片| 欧美av亚洲av综合av国产av| 18在线观看网站| 看免费av毛片| 精品人妻熟女毛片av久久网站| 欧美精品av麻豆av| 久久天堂一区二区三区四区| 2018国产大陆天天弄谢| 日韩电影二区| 国产激情久久老熟女| 亚洲欧美色中文字幕在线| 久久精品亚洲av国产电影网| 天堂8中文在线网| 亚洲精品久久久久久婷婷小说| 搡老乐熟女国产| 91成人精品电影| 这个男人来自地球电影免费观看| 亚洲欧美清纯卡通| 亚洲精品美女久久av网站| 亚洲欧洲日产国产| 日韩 亚洲 欧美在线| 日韩电影二区| av超薄肉色丝袜交足视频| 777米奇影视久久| 两性午夜刺激爽爽歪歪视频在线观看 | 黄片播放在线免费| 丁香六月天网| 亚洲美女黄色视频免费看| 啦啦啦视频在线资源免费观看| 亚洲全国av大片| 91av网站免费观看| 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 亚洲七黄色美女视频| 国产免费视频播放在线视频| 亚洲视频免费观看视频| 亚洲免费av在线视频| 精品少妇一区二区三区视频日本电影| 欧美激情极品国产一区二区三区| 欧美亚洲 丝袜 人妻 在线| 又黄又粗又硬又大视频| 18禁国产床啪视频网站| 一区二区av电影网| 国产精品久久久久久精品电影小说| 亚洲av美国av| 狂野欧美激情性xxxx| 成人av一区二区三区在线看 | 99国产精品99久久久久| 午夜免费观看性视频| av福利片在线| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 一区二区三区四区激情视频| 国产真人三级小视频在线观看| 国产色视频综合| 岛国在线观看网站| 老汉色av国产亚洲站长工具| 亚洲综合色网址| 成人av一区二区三区在线看 | 午夜福利视频精品| 精品人妻在线不人妻| 欧美成狂野欧美在线观看| 午夜福利视频精品| 成年动漫av网址| 啪啪无遮挡十八禁网站| 我要看黄色一级片免费的| 免费在线观看完整版高清| 国产麻豆69| 一区二区三区激情视频| 成年美女黄网站色视频大全免费| 国产成人啪精品午夜网站| 亚洲成人免费av在线播放| 午夜福利乱码中文字幕| 操出白浆在线播放| 国产在视频线精品| 1024香蕉在线观看| 女警被强在线播放| 亚洲精品国产区一区二| 1024视频免费在线观看| 91精品伊人久久大香线蕉| 午夜福利,免费看| 久久ye,这里只有精品| 丝袜脚勾引网站| 日韩三级视频一区二区三区| 久久久国产一区二区| 中文字幕制服av| 亚洲激情五月婷婷啪啪| 在线 av 中文字幕| 99精国产麻豆久久婷婷| 啦啦啦中文免费视频观看日本| 日韩欧美一区二区三区在线观看 | 日韩欧美一区二区三区在线观看 | 欧美av亚洲av综合av国产av| 色综合欧美亚洲国产小说| 欧美av亚洲av综合av国产av| 法律面前人人平等表现在哪些方面 | 日韩大码丰满熟妇| 亚洲精品国产区一区二| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| 大码成人一级视频| 亚洲第一av免费看| 日本猛色少妇xxxxx猛交久久| 久热这里只有精品99| 国产高清国产精品国产三级| 亚洲国产日韩一区二区| 久久 成人 亚洲| 亚洲中文av在线| 午夜91福利影院| 成年美女黄网站色视频大全免费| 九色亚洲精品在线播放| 成年人黄色毛片网站| 中文字幕制服av| 国产一区二区在线观看av| 丝袜喷水一区| 丰满人妻熟妇乱又伦精品不卡| 亚洲,欧美精品.| 午夜免费鲁丝| 黑丝袜美女国产一区| 欧美性长视频在线观看| 丝袜人妻中文字幕| 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 久久精品国产亚洲av香蕉五月 | 国产精品免费大片| 精品乱码久久久久久99久播| 啦啦啦视频在线资源免费观看| 黄色视频,在线免费观看| 99久久99久久久精品蜜桃| 日本a在线网址| 久热爱精品视频在线9| 如日韩欧美国产精品一区二区三区| 亚洲精品国产av蜜桃| 欧美日韩成人在线一区二区| 老汉色∧v一级毛片| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 国产成人欧美在线观看 | 日本91视频免费播放| 国产亚洲欧美精品永久| 日本黄色日本黄色录像| 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯 | 不卡av一区二区三区| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久久精品古装| 国产成人一区二区三区免费视频网站| 久久国产精品男人的天堂亚洲| av一本久久久久| 国产一区二区三区综合在线观看| 一本久久精品| av超薄肉色丝袜交足视频| 啦啦啦啦在线视频资源| 色婷婷av一区二区三区视频| 在线观看免费视频网站a站| 91大片在线观看| 国产在线免费精品| 热99国产精品久久久久久7| 免费看十八禁软件| 91字幕亚洲| av又黄又爽大尺度在线免费看| 热99久久久久精品小说推荐| 国产免费视频播放在线视频| 亚洲av日韩在线播放| 亚洲成人免费电影在线观看| 久久这里只有精品19| 男人添女人高潮全过程视频| 精品一区二区三区av网在线观看 | 国产精品 国内视频| 亚洲精品中文字幕在线视频| 最新在线观看一区二区三区| 国产又色又爽无遮挡免| 精品高清国产在线一区| 欧美精品一区二区大全| 飞空精品影院首页| av又黄又爽大尺度在线免费看| 久久天堂一区二区三区四区| 亚洲中文日韩欧美视频| 亚洲黑人精品在线| 久久久久久久大尺度免费视频| 亚洲成国产人片在线观看| 久久毛片免费看一区二区三区| 国产99久久九九免费精品| 少妇裸体淫交视频免费看高清 | 亚洲黑人精品在线| 午夜视频精品福利| 午夜激情久久久久久久| 国产在线视频一区二区| 免费日韩欧美在线观看| 国产色视频综合| 亚洲少妇的诱惑av| 欧美国产精品一级二级三级| 大香蕉久久网| 国产精品久久久久久精品电影小说| 亚洲色图 男人天堂 中文字幕| 91成人精品电影| 国产人伦9x9x在线观看| 日韩视频一区二区在线观看| 侵犯人妻中文字幕一二三四区| 一进一出抽搐动态| 欧美黑人精品巨大| 免费一级毛片在线播放高清视频 | 老司机福利观看| 中文字幕人妻丝袜一区二区| 欧美人与性动交α欧美软件| 久久九九热精品免费| 男女午夜视频在线观看| 成人18禁高潮啪啪吃奶动态图| 日韩视频一区二区在线观看| 国产亚洲av高清不卡| 亚洲国产av新网站| 少妇被粗大的猛进出69影院| 丝袜美足系列| 久久精品成人免费网站| 丝瓜视频免费看黄片| 欧美国产精品一级二级三级| 国产免费现黄频在线看| 一二三四在线观看免费中文在| 视频区图区小说| 午夜免费观看性视频| 国产1区2区3区精品| 人成视频在线观看免费观看| 999久久久国产精品视频| 熟女少妇亚洲综合色aaa.| 黄色视频不卡| 精品福利观看| 国产精品秋霞免费鲁丝片| av视频免费观看在线观看| 国产成人a∨麻豆精品| a级毛片在线看网站| 日韩制服丝袜自拍偷拍| 正在播放国产对白刺激| 国产av精品麻豆| 日韩中文字幕欧美一区二区| 一本色道久久久久久精品综合| 精品乱码久久久久久99久播| 国产成+人综合+亚洲专区| 超色免费av| 国产亚洲av片在线观看秒播厂| 国产一区二区 视频在线| 肉色欧美久久久久久久蜜桃| 999精品在线视频| 午夜激情久久久久久久| 国产三级黄色录像| 一二三四在线观看免费中文在| 久久这里只有精品19| 国产不卡av网站在线观看| 欧美乱码精品一区二区三区| 成人手机av| 搡老岳熟女国产| 一区二区三区四区激情视频| 美女大奶头黄色视频| 日韩欧美国产一区二区入口| 亚洲欧洲精品一区二区精品久久久| 日韩视频一区二区在线观看| 性色av一级| 自拍欧美九色日韩亚洲蝌蚪91| www.av在线官网国产| 9热在线视频观看99| 一区二区av电影网| 久久中文看片网| 不卡av一区二区三区| av福利片在线| 国产精品久久久人人做人人爽| 天天躁日日躁夜夜躁夜夜| 成年人黄色毛片网站| 国产真人三级小视频在线观看| 日日摸夜夜添夜夜添小说| 1024视频免费在线观看| 亚洲成人免费av在线播放| 精品少妇久久久久久888优播| 纵有疾风起免费观看全集完整版| 自线自在国产av| 久久影院123| 国产免费福利视频在线观看| 午夜精品国产一区二区电影| 亚洲精品国产色婷婷电影| 91av网站免费观看| 午夜视频精品福利| 国产片内射在线| 中文字幕人妻丝袜制服| 美女视频免费永久观看网站| 夫妻午夜视频| 久久中文字幕一级| 999精品在线视频| 多毛熟女@视频| 国内毛片毛片毛片毛片毛片| 12—13女人毛片做爰片一| 一级毛片电影观看| 美女脱内裤让男人舔精品视频| 国产精品久久久久久精品古装| 免费久久久久久久精品成人欧美视频| 国产男女超爽视频在线观看| 丰满饥渴人妻一区二区三| 欧美日韩亚洲国产一区二区在线观看 | 黄色毛片三级朝国网站| 亚洲精品日韩在线中文字幕| 国产精品影院久久| 午夜免费鲁丝| 久久天躁狠狠躁夜夜2o2o| av在线老鸭窝| 少妇猛男粗大的猛烈进出视频| 亚洲精品久久久久久婷婷小说| 日本欧美视频一区| 大片免费播放器 马上看| 国产一区二区三区av在线| 日韩欧美一区视频在线观看|