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

    管道內(nèi)障礙物對(duì)加氫甲烷爆炸特性的影響

    2016-12-14 09:02:13余明高袁晨樵鄭凱
    化工學(xué)報(bào) 2016年12期
    關(guān)鍵詞:傳播速度障礙物氫氣

    余明高,袁晨樵,鄭凱

    ?

    管道內(nèi)障礙物對(duì)加氫甲烷爆炸特性的影響

    余明高1,2,袁晨樵1,鄭凱2

    (1河南理工大學(xué)安全科學(xué)與工程學(xué)院,河南焦作 454003;2重慶大學(xué)煤礦災(zāi)害動(dòng)力學(xué)與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室,重慶 400044)

    在搭建小尺寸爆炸平臺(tái)的基礎(chǔ)上,研究了當(dāng)量比為1時(shí),氫氣體積分?jǐn)?shù)分別為0%、25%、50%、75%和100%時(shí),障礙物條件下添加氫氣對(duì)管道內(nèi)甲烷-空氣預(yù)混火焰?zhèn)鞑ヌ匦缘挠绊?。?shí)驗(yàn)結(jié)果表明:最大爆炸壓力、火焰?zhèn)鞑ニ俣榷紩?huì)隨著氫含量和障礙物數(shù)量的增加而增大,并且火焰顏色和火焰形狀也會(huì)隨之發(fā)生變化。但是,障礙物對(duì)火焰平均傳播速度的影響較為微弱,并且障礙物對(duì)甲烷和氫氣的火焰?zhèn)鞑ニ俣鹊挠绊懶Ч灿兴煌W畲蟊▔毫εc火焰顏色隨氫含量的變化趨勢(shì)隨著障礙物數(shù)量的增加越來(lái)越明顯。氫氣與甲烷相比,障礙物對(duì)前者爆炸特性的影響更為明顯。

    氫氣;障礙物;甲烷/空氣;爆炸特性

    引 言

    甲烷爆炸是煤礦重大惡性事故之一,給煤礦安全生產(chǎn)帶來(lái)了極大的威脅。為了有效防止煤礦甲烷事故的發(fā)生,許多學(xué)者都對(duì)甲烷的爆炸機(jī)理、爆炸特征及其影響因素進(jìn)行了研究,如鄭立剛等[1]研究了在管道口一端閉口一端開口條件下,點(diǎn)火源位置對(duì)甲烷-空氣預(yù)混氣爆燃超壓特征的影響。溫小萍等[2]基于火焰動(dòng)態(tài)傳播和超壓信號(hào)的高速同步采集,實(shí)驗(yàn)研究了不同湍流激勵(lì)條件下瓦斯爆炸火焰結(jié)構(gòu)與壓力波的耦合關(guān)系。然而,煤礦井下采空區(qū)、封閉火區(qū)等區(qū)域中的瓦斯除了甲烷氣體之外,還存在著一定量的氫氣[3]。同時(shí),煤體在一定條件下也會(huì)產(chǎn)生一定量的氫氣。例如,在常溫、常壓下,自然狀態(tài)的煤體可以通過脫附、解吸等作用下釋放出一定量的氫氣。根據(jù)前人研究結(jié)果[4-7]可知,添加氫氣可以顯著改變甲烷-空氣預(yù)混火焰?zhèn)鞑ヌ匦?,如提高火焰?zhèn)鞑ニ俣?、增加預(yù)混火焰溫度及爆炸反應(yīng)活性等。顯然,這些特性的改變?cè)黾恿送咚贡ㄎkU(xiǎn)性,因此有必要開展添加氫氣對(duì)瓦斯爆炸影響的研究。

    目前學(xué)者們對(duì)添加氫氣對(duì)瓦斯爆炸的影響進(jìn)行了一系列初步研究。如李增華等[3]研究了添加氫氣會(huì)降低甲烷的爆炸下限,增強(qiáng)其爆炸危害性。Ma等[8-9]分別通過球形容器和數(shù)值模擬研究了添加氫氣對(duì)甲烷/空氣爆炸特性的影響以及在開口和密閉條件下加氫對(duì)甲烷/空氣爆炸過程的影響,研究結(jié)果表明爆炸壓力會(huì)隨著氫含量的增加而增加,并且當(dāng)量比情況下壓力最大,以及在開口條件下加氫對(duì)甲烷/空氣的影響要強(qiáng)于密閉條件。Yu等[10]利用開口管道研究了加氫對(duì)瓦斯爆炸特征的影響,研究結(jié)果表明隨著氫含量的增加氣體爆炸壓力和火焰?zhèn)鞑ニ俣榷紩?huì)隨之提高,并且火焰結(jié)構(gòu)也會(huì)隨之發(fā)生變化。賈寶山等[11]通過數(shù)值模擬研究了受限空間內(nèi)瓦斯爆炸與氫氣促進(jìn)機(jī)理,結(jié)果表明,隨著混合氣中氫氣含量的增加,瓦斯引爆時(shí)間越來(lái)越短,其爆炸強(qiáng)度也隨之增大,且氫氣在一定程度上對(duì)有害氣體CO、CO2、NO、NO2的生成有很大影響。Porwski等[12]實(shí)驗(yàn)研究了在障礙物管道內(nèi)氫氣-甲烷-空氣混合氣體爆炸火焰的傳播、加速和爆轟,以及決定爆燃和爆轟的區(qū)域和火焰?zhèn)鞑サ乃俣?。Salzano等[13]研究了在密閉的柱形容器內(nèi)不同含量的氫氣/甲烷混合氣體和初始?jí)毫?duì)最大爆炸壓力、最大爆炸壓力上升速率和燃燒速率的影響。Woolley等[14]通過建立密閉和開放的數(shù)學(xué)模型,研究了氫氣/甲烷混合氣體的爆炸特性,結(jié)果表明加氫對(duì)氣體的爆炸壓力有顯著的影響。這些研究主要關(guān)注于添加氫氣對(duì)爆炸壓力、溫度等參數(shù)的影響,而對(duì)于預(yù)混火焰?zhèn)鞑ヌ卣鞯难芯枯^少,因此有必要開展添加氫氣對(duì)瓦斯爆炸預(yù)混火焰?zhèn)鞑ヌ卣鞯难芯俊?/p>

    煤礦井下不可避免地存在各種障礙物,研究表明,障礙物的存在能夠誘導(dǎo)火焰由層流向湍流轉(zhuǎn)變,從而提高爆炸反應(yīng)速度,增加爆炸危險(xiǎn)性。潘鵬飛等[15]通過密閉管道研究了條形障礙物對(duì)瓦斯爆炸的影響,研究表明隨著障礙物數(shù)量的增加,爆炸壓力、火焰?zhèn)鞑ニ俣取⒈▔毫ι仙俾室约氨ㄖ笖?shù)都會(huì)隨之增加。何學(xué)秋等[16]研究了障礙物對(duì)瓦斯爆炸火焰結(jié)構(gòu)及火焰?zhèn)鞑サ挠绊懀Y(jié)果表明障礙物會(huì)使火焰鋒面褶皺度增大,促進(jìn)了火焰加速。Alharbi等[17]研究了在障礙物管道內(nèi)氫氣的火焰?zhèn)鞑ィ瑢?shí)驗(yàn)表明障礙物阻塞率增大隨之氣體的峰值壓力以及變化率也會(huì)增大。除此之外,根據(jù)Ibrhaim等[18]的研究可知障礙物會(huì)對(duì)爆炸波的誘導(dǎo)、疊加、反射規(guī)律以及引起爆轟產(chǎn)生很大的影響。王海賓等[19]研究了水平管道內(nèi)障礙物對(duì)爆炸壓力的影響,結(jié)果表明隨著障礙物阻塞率的增加爆炸壓力會(huì)相應(yīng)增大,而改變障礙物的間距則對(duì)爆炸過程并無(wú)太大影響。余明高等[20]分別研究了交錯(cuò)障礙物對(duì)瓦斯爆炸特征的影響,認(rèn)為交錯(cuò)障礙物明顯增強(qiáng)了火焰的形變以及提高了火焰?zhèn)鞑ニ俣群捅▔毫?。林柏泉等[21]研究了障礙物對(duì)爆炸火焰以及爆炸波的影響,認(rèn)為障礙物增強(qiáng)了火焰?zhèn)鞑サ耐牧鳜F(xiàn)象以及爆炸波的變化幅度。尉存娟等[22]和Hall等[23]研究了管道內(nèi)障礙物數(shù)量對(duì)瓦斯爆炸過程的影響,前者研究表明爆炸壓力及其上升速率都隨障礙物數(shù)量增加呈先增大后減小的變化規(guī)律,而火焰?zhèn)鞑ニ俣葎t隨障礙物數(shù)量增加單調(diào)遞增,但其增幅較?。缓笳邉t認(rèn)為在一定范圍內(nèi),障礙物數(shù)量越多爆炸壓力越大。Gubba等[24]綜合研究了障礙物距離點(diǎn)火位置、障礙物數(shù)量以及障礙物間距對(duì)預(yù)混氣體爆炸的影響。Park等[25]在開放的燃燒室內(nèi)研究了不同長(zhǎng)徑比和阻塞率條件下,多組障礙物對(duì)氣體最大爆炸壓力和火焰?zhèn)鞑ニ俣鹊挠绊?。郭丹彤等[26]運(yùn)用流體力學(xué)軟件AutoRea Gas建立不同阻塞率和不同結(jié)構(gòu)的障礙物爆炸模型,模擬分析不同布置情況對(duì)氣體爆炸壓力場(chǎng)的影響程度和規(guī)律,研究表明同種障礙物結(jié)構(gòu)下,隨著阻塞率的增加,氣體爆炸壓力的增加程度在一定范圍內(nèi)呈現(xiàn)出先增大后減小的變化情況;相同阻塞率下,立體障礙物對(duì)氣體爆炸壓力場(chǎng)產(chǎn)生的影響明顯大于平面障礙物。

    然而,以上研究中還存在著一些缺陷,其研究對(duì)象主要還局限于甲烷/空氣等單一氣體,卻忽略了瓦斯中氫氣的存在。因此,本研究搭建了相關(guān)的實(shí)驗(yàn)平臺(tái),就障礙物管道內(nèi)添加氫氣對(duì)瓦斯爆炸特性的影響進(jìn)行了實(shí)驗(yàn)研究,研究對(duì)象主要為爆炸預(yù)混火焰?zhèn)鞑ヌ卣鳌㈩A(yù)混火焰?zhèn)鞑ニ俣燃氨▔毫Φ?,以期?duì)瓦斯爆炸的防治起到一定的理論幫助。

    1 實(shí)驗(yàn)裝置和實(shí)驗(yàn)方法

    1.1 實(shí)驗(yàn)裝置

    為了研究障礙物管道內(nèi)加氫對(duì)甲烷爆炸特性的影響,搭建了有關(guān)在障礙物管道內(nèi)測(cè)試當(dāng)量比為1時(shí)不同體積分?jǐn)?shù)混合氣體爆炸特性的實(shí)驗(yàn)系統(tǒng)。系統(tǒng)主要由水平管道、配氣系統(tǒng)、壓力測(cè)量系統(tǒng)、點(diǎn)火裝置、光電測(cè)量系統(tǒng)組成。其中,管道是實(shí)驗(yàn)的主體部分,所有實(shí)驗(yàn)都在此管道內(nèi)進(jìn)行。管道由截面為100 mm×100 mm,高為500 mm,管壁厚度為20 mm的有機(jī)玻璃制成,抗壓2 MPa以上。管道底部由不透鋼鋼板封閉,并且在鋼板上開了3個(gè)小孔,這是為了安置進(jìn)氣閥門、壓力傳感器和點(diǎn)火器。管道頂部由保鮮薄膜封閉,起到泄爆作用;配氣系統(tǒng)由空氣壓縮機(jī)、甲烷氣體瓶、氫氣氣體瓶、流量控制器以及三通管道構(gòu)成;壓力傳感器采用的是上海銘動(dòng)的高頻壓力傳感器,測(cè)試范圍為-0.1~0.1 MPa,然后通過數(shù)據(jù)采集卡采集,采集速率為15 kHz;點(diǎn)火裝置由電子點(diǎn)火器、穩(wěn)壓電源和點(diǎn)火控制器組成,電子點(diǎn)火器工作電壓為6 V;為了便于觀測(cè),光電傳感器安置在管道側(cè)壁的下方。除此之外,為了確保實(shí)驗(yàn)的安全性,本實(shí)驗(yàn)是在室外進(jìn)行的。同時(shí),為了確?;鹧鎴D像的準(zhǔn)確性和清晰度,實(shí)驗(yàn)一般在晚上進(jìn)行?;鹧?zhèn)鞑ミ^程由分辨率為1024×1024的德國(guó)LaVision高速攝像儀以5000 frame·s-1的速度進(jìn)行高速拍攝,以捕捉瓦斯爆炸過程中的火焰顏色和火焰形狀。

    為了對(duì)比分析不同數(shù)量的障礙物管道內(nèi)加氫對(duì)爆炸特性的影響,本實(shí)驗(yàn)分別在4個(gè)工況下進(jìn)行。其中,工況1為空管,工況2為1組障礙物管道,工況3為2組障礙物管道,工況4為3組障礙物管道,如圖1所示。

    圖1 實(shí)驗(yàn)工況

    1.2 實(shí)驗(yàn)方法

    分別對(duì)4個(gè)工況進(jìn)行實(shí)驗(yàn)對(duì)比,其中障礙物的阻塞率為0.3。實(shí)驗(yàn)過程中,先對(duì)氣體進(jìn)行配氣,其中燃?xì)庵袣錃怏w積分?jǐn)?shù)分別為0、25%、50%、75%和100%。通過流量控制器調(diào)配氣體的流量,經(jīng)三通閥將氣體預(yù)混通入管道內(nèi),在管道側(cè)壁上方靠近頂口處設(shè)有排氣口。依據(jù)文獻(xiàn)[11],調(diào)節(jié)好通氣流量,通氣總量大約為4倍的容器體積,通氣時(shí)間持續(xù)約為6~7 min,以保證排盡管道內(nèi)原有的氣體以及管道內(nèi)的氣體均勻混合。通氣完畢后,關(guān)閉流量控制器的電源以及管道進(jìn)出口閥門。待氣體靜止大約30 s后,點(diǎn)火引爆,對(duì)數(shù)據(jù)進(jìn)行采集。

    氫氣的體積分?jǐn)?shù)表示為

    燃?xì)猱?dāng)量比表示為

    式中,/為實(shí)際通入的燃?xì)馀c空氣體積比;(/)stoic表示為燃?xì)馀c空氣的化學(xué)計(jì)量比。當(dāng)量比小于1為貧燃,當(dāng)量比為1是化學(xué)當(dāng)量比,當(dāng)量比大于1為富燃。

    2 結(jié)果與分析

    2.1 對(duì)爆炸壓力的影響

    圖2表示在不同工況下不同氫含量與最大爆炸壓力的變化曲線。從圖中可以看出,無(wú)論是空管還是障礙物管道,加氫對(duì)氣體的最大爆炸壓力有著很明顯的影響,隨著氫含量的增加而增大,并且上升趨勢(shì)也越來(lái)越明顯。例如,在工況3下,氫含量為0時(shí),最大爆炸壓力為7.02 kPa。氫含量為25%時(shí),最大爆炸壓力為10.65 kPa,比純甲烷提高了51.7%。氫含量為50%時(shí),最大爆炸壓力為16.75 kPa,比純甲烷提高了137%。氫含量為75%時(shí),最大爆炸壓力為32.3 kPa,比純甲烷提高了360.1%。氫含量為100%時(shí),最大爆炸壓力為66.23 kPa,比純甲烷提高了843%。這是由于加氫可以提高反應(yīng)中自由基的濃度[11],加快了預(yù)混火焰的燃燒反應(yīng)速率并增加了反應(yīng)中產(chǎn)生的熱量,從而促進(jìn)了氣體爆炸。

    圖2 不同工況下氫含量與最大爆炸壓力的變化曲線

    從圖2中也可以看出,同一情況下,最大爆炸壓力隨著障礙物數(shù)量的增加而增加。例如,氫含量為0時(shí),工況1的最大爆炸壓力為5.01 kPa。然而,在障礙物的干擾作用下,火焰湍流強(qiáng)度增強(qiáng),爆炸壓力隨之增大[14]。工況2的最大爆炸壓力比工況1提高了22%,為6.11 kPa。工況3比工況1提高了40%,最大爆炸壓力為7.02 kPa。工況4比工況1提高了60%,最大爆炸壓力為8 kPa。同理,氫含量為25%、50%、75%和100%時(shí),其余3個(gè)工況的最大爆炸壓力相比于工況1都會(huì)有所提高。然而不同的是,氫含量為75%和100%時(shí),最大爆炸壓力隨障礙物數(shù)量的增長(zhǎng)幅度要明顯大于其余3種體積分?jǐn)?shù)。由此說明,氫含量小于50%時(shí),最大爆炸壓力隨障礙物數(shù)量的上升趨勢(shì)較為微弱,而氫含量大于50%時(shí),最大爆炸壓力隨障礙物數(shù)量的上升則較為明顯。換言之,障礙物對(duì)氫氣爆炸壓力的影響要強(qiáng)于它對(duì)甲烷的影響。

    分析上述原因:由于管道內(nèi)障礙物的存在,火焰波陣面會(huì)發(fā)生湍流現(xiàn)象,增大了火焰與氣體的接觸面積,前驅(qū)沖擊波對(duì)未燃混合物的加熱和壓縮的正反饋機(jī)理越顯著,使爆炸反應(yīng)速度和能量釋放速度隨之加快,氣體爆炸壓力和爆炸指數(shù)會(huì)隨之增大。障礙物數(shù)量的增加,湍流效應(yīng)增強(qiáng),隨之爆炸壓力和爆炸指數(shù)也會(huì)增大[14]。其中,根據(jù)文獻(xiàn)[27]可知,爆炸指數(shù)為工藝單元危險(xiǎn)系數(shù)和物質(zhì)系數(shù)的乘積。爆炸指數(shù)與危險(xiǎn)程度呈正比,爆炸指數(shù)越大,危險(xiǎn)程度越高。另外,隨著氫含量的增加火焰燃燒反應(yīng)速率增強(qiáng)[10]。同時(shí)在障礙物的擾動(dòng)作用下火焰鋒面褶皺現(xiàn)象加強(qiáng),從而極大促進(jìn)了火焰燃燒反應(yīng)速率,爆炸壓力也就隨之明顯增大。因此障礙物對(duì)氫氣爆炸壓力的影響更為突出。

    由圖2還可以看出,就整體而言,隨著障礙物數(shù)量的增加曲線斜率越來(lái)越大。由此表明,最大爆炸壓力隨氫含量的上升趨勢(shì)會(huì)隨著障礙物數(shù)量增加越來(lái)越明顯。在工況1下,氫含量從0逐漸變化到100%,最大爆炸壓力的平均增長(zhǎng)率為62.2%。由于障礙物的激勵(lì)作用[28],在工況2下,氫含量從0逐漸變化到100%,最大爆炸壓力的平均增長(zhǎng)率比工況1提高了6.23%,為68.43%。同理,在工況3和工況4下,氫含量從0逐漸變化到100%,最大爆炸壓力的平均增長(zhǎng)率相比于工況1也都有所提高,分別提高了13.8%和20.6%,為76%和82.8%。

    圖3為空管和兩組障礙物管道內(nèi)不同氫含量的時(shí)間-爆炸壓力關(guān)系。由圖中可以看出,無(wú)論是在空管還是障礙物管道內(nèi),氫含量小于50%時(shí)會(huì)出兩次壓力峰值,分別是泄爆壓力和最大爆炸壓力。而氫含量大于50%時(shí)則只出現(xiàn)了一次壓力峰值,也就是最大爆炸壓力。這是由于爆炸過程中薄膜破裂,部分未燃?xì)怏w從泄口排出,這就導(dǎo)致了壓力降低[9],從而出現(xiàn)了第1次壓力峰值,也就是泄爆壓力。隨著火焰的傳播,燃燒反應(yīng)速率漸漸超過了泄爆速率,壓力開始上升,直到燃燒反應(yīng)速率等于泄爆速率時(shí),第2次壓力峰值出現(xiàn),也就是最大爆炸壓力。然而,隨著氫含量的增加,層流火焰?zhèn)鞑ニ俣仍龃螅瑥亩鴮?dǎo)致燃燒反應(yīng)與泄爆之間的速率差增大[11]。從圖中也可以看出,隨著氫含量的增加,到達(dá)最大壓力峰值的時(shí)間縮短。但是,最大壓力峰值出現(xiàn)時(shí)刻并不完全與障礙物數(shù)量有關(guān)。圖4為氫含量為50%時(shí)不同工況下時(shí)間-爆炸壓力曲線關(guān)系。從圖中可以發(fā)現(xiàn),工況4時(shí)到達(dá)最大壓力峰值時(shí)刻最晚,1組障礙物管道時(shí)到達(dá)最大壓力峰值時(shí)刻最早。這就說明最大壓力峰值出現(xiàn)時(shí)刻也與火焰-湍流耦合效應(yīng)發(fā)生時(shí)刻有關(guān)。

    圖3 工況1與工況3下不同氫含量的時(shí)間-爆炸壓力關(guān)系

    圖4 氫含量為50%時(shí)不同障礙物管道內(nèi)時(shí)間-爆炸壓力關(guān)系

    2.2 對(duì)火焰顏色和火焰形狀的影響

    圖5為不同工況下不同氫含量的爆炸火焰圖像。從圖中可以看出,無(wú)論是空管還是障礙物管道,隨著氫含量的增加火焰變得越來(lái)越明亮且顏色也發(fā)生了變化。這主要是由于隨著氫含量的增加反應(yīng)氣體中自由基的濃度隨之增加,以及不同體積分?jǐn)?shù)的氫氣的光譜不同所導(dǎo)致的。除此之外,在空管和1組障礙物管道內(nèi),氫含量為50%時(shí)火焰顏色發(fā)生顯著變化變?yōu)榧t色。在2組和3組障礙物管道內(nèi),氫含量為25%時(shí)火焰顏色就變?yōu)榱思t色,氫含量為100%時(shí)火焰亮度變得十分耀眼。同時(shí)在障礙物管道內(nèi)火焰繞過障礙物后,火焰顏色會(huì)慢慢發(fā)生變化[29]。以上現(xiàn)象說明了由于障礙物的存在,火焰湍流強(qiáng)度增強(qiáng),火焰鋒面表面積增大,火焰內(nèi)部流場(chǎng)各反應(yīng)物質(zhì)熱運(yùn)動(dòng)能增大,從而引起燃燒反應(yīng)速率加快和熱釋放速率加快。因此隨著障礙物數(shù)量的增加火焰顏色也隨之變化,并且火焰顏色隨氫含量的變化趨勢(shì)也越明顯。同時(shí),由于甲烷和氫氣發(fā)生燃燒反應(yīng)時(shí)自由基濃度以及其光譜不同導(dǎo)致純甲烷和純氫在不同工況時(shí)火焰顏色的變化有明顯的區(qū)別,純甲烷時(shí)火焰顏色隨障礙物數(shù)量的增加由黃變紅,純氫氣時(shí)火焰顏色隨障礙物數(shù)量的增加由紅變得十分明亮耀眼。從圖中還可以發(fā)現(xiàn),燃?xì)鉃榧兗淄闀r(shí)隨著障礙物數(shù)量的增加,火焰顏色的變化較為微弱,但隨著氫含量的增加障礙物對(duì)火焰顏色的影響則較為明顯。這就表明了對(duì)于火焰顏色的影響,障礙物對(duì)氫氣的作用更為突出。這是由于隨著氫含量的增加火焰鋒面會(huì)出現(xiàn)褶皺現(xiàn)象,在障礙物作用下褶皺現(xiàn)象更為明顯,大大增加了火焰鋒面?zhèn)髋c氧氣的接觸面積,從而加快了火焰反應(yīng)速率和熱釋放速率,火焰顏色變化明顯。

    由圖5也可以看出,障礙物和氫氣都會(huì)對(duì)火焰的形狀產(chǎn)生很大的影響。隨著氫含量的增加火焰鋒面會(huì)出現(xiàn)褶皺現(xiàn)象。例如在空管內(nèi),氫含量小于50%時(shí)火焰鋒面較為平滑,當(dāng)氫含量大于50%時(shí)火焰鋒面開始有褶皺出現(xiàn)。在障礙物管道內(nèi),同一條件下,剛剛開始時(shí)火焰形狀與空管內(nèi)的火焰形狀一致。但當(dāng)火焰?zhèn)鞑サ秸系K物到附近時(shí),由于障礙物的存在,迫使火焰向中心發(fā)展,火焰開始發(fā)生形變。然后火焰鋒面從兩個(gè)障礙物之間的間隙噴出。當(dāng)火焰穿過障礙物后,火焰向兩側(cè)發(fā)展,火焰鋒面變寬,鋒面表面積增大[30]。當(dāng)火焰大約傳至管道中部時(shí),火焰鋒面開始出現(xiàn)褶皺現(xiàn)象。以上現(xiàn)象說明了隨著障礙物數(shù)量的增加火焰鋒面會(huì)出現(xiàn)明顯的褶皺現(xiàn)象,并且火焰形變也會(huì)越來(lái)越明顯。文獻(xiàn)[31]中表示火焰鋒面受到干擾后會(huì)改變其穩(wěn)定狀態(tài),從而使火焰鋒面產(chǎn)生褶皺。燃燒學(xué)認(rèn)為光滑的層流火焰表面表現(xiàn)出來(lái)任何變化都認(rèn)為是火焰湍流現(xiàn)象[26]。因此火焰鋒面因障礙物干擾或其他因素影響所導(dǎo)致的火焰鋒面褶皺現(xiàn)象就是火焰湍流化的表現(xiàn)。

    圖5 不同工況下不同氫含量的爆炸火焰圖像

    2.3 對(duì)火焰速度的影響

    已知管道長(zhǎng)度為500 mm,根據(jù)高速攝像機(jī)記錄到的火焰從點(diǎn)火到管口所用的時(shí)間,由此可以計(jì)算得到火焰平均傳播速度。火焰平均傳播速度可以展現(xiàn)出在不同工況下以及不同氫含量條件下火焰整體的速度變化。圖6所示的是不同工況下不同氫含量與火焰平均傳播速度的變化曲線。從圖中可以看出,在同一條件下,隨著障礙物數(shù)量的增加火焰平均傳播速度越來(lái)越快。例如,當(dāng)氫含量為0時(shí),工況1的火焰平均傳播速度為11.63 m·s-1。然而,在障礙物的干擾下,火焰鋒面出現(xiàn)褶皺,增大了火焰鋒面的表面積,火焰的傳播速度隨之增大。工況2的火焰平均傳播速度就比工況1提高了17.8%,為13.70 m·s-1。工況3比工況1提高了21.1%,為14.08 m·s-1。工況4比工況1高了36.5%,為15.87 m·s-1。同理,當(dāng)氫含量為25%、50%、75%和100%時(shí),相較于工況1,其余3個(gè)工況的火焰平均都會(huì)有所提高,并且隨著障礙物數(shù)量的增加火焰平均傳播速度越來(lái)越大。除此之外,在同一工況下,隨著氫含量的增加火焰平均傳播速度也越來(lái)越快。例如,在工況3下,氫含量為0時(shí),火焰平均傳播速度為14.08 m·s-1。然而,隨著氫含量的增加火焰燃燒反應(yīng)速率和熱釋放速率加快,火焰?zhèn)鞑ニ俣纫簿碗S之增大。氫含量為0.25時(shí),火焰平均傳播速度為17.54 m·s-1,比純甲烷提高了24.6%。氫含量為50%時(shí),火焰平均傳播速度為23.26 m·s-1,比純甲烷提高了65.2%。氫含量為75%時(shí),火焰平均傳播速度為41.67 m·s-1,比純甲烷提高了196%。氫含量為100%時(shí),火焰平均傳播速度為83.33 m·s-1,比純甲烷提高了492%。同理,其余3個(gè)工況下,火焰?zhèn)鞑ニ俣纫矔?huì)隨著氫含量的增加而增大。

    對(duì)比圖6和圖2還可以發(fā)現(xiàn),相較于障礙物對(duì)最大爆炸壓力的影響,火焰平均傳播速度隨障礙物數(shù)量的增長(zhǎng)趨勢(shì)較為微弱。這就表明障礙物對(duì)火焰平均傳播速度的影響較小。這是由于火焰在穿越障礙物時(shí)會(huì)有部分能量損失[22],從而導(dǎo)致了部分速度的損失。另一方面火焰在穿過障礙物之后,會(huì)向兩邊拉伸,使其縱向速度減小[32]。

    圖6 不同工況下氫含量與火焰平均傳播速度的變化曲線

    根據(jù)高速攝像機(jī)記錄得到的火焰鋒面位置隨時(shí)間的變化,以及測(cè)量出的從點(diǎn)火到火焰鋒面位置的距離,由此可以計(jì)算得出火焰?zhèn)鞑ニ俣??;鹧鎮(zhèn)鞑ニ俣瓤梢愿谜宫F(xiàn)出隨著火焰鋒面位置變化,火焰速度的變化情況。障礙物對(duì)含不同氫含量的瓦斯爆炸火焰?zhèn)鞑ニ俣鹊挠绊懶Ч兄黠@的不同。圖7為氫含量分別為0、50%以及100%時(shí),不同工況下火焰?zhèn)鞑ゾ嚯x與火焰?zhèn)鞑ニ俣戎g的曲線關(guān)系。從圖中可知,在相同條件下,火焰?zhèn)鞑コ跗?個(gè)工況的火焰速度基本一致。氫含量為0和50%時(shí)火焰?zhèn)鞑ニ俣入S著障礙物的添加出現(xiàn)了速度峰值,并且峰值出現(xiàn)的次數(shù)與障礙物數(shù)量是一致的。速度峰值依次出現(xiàn)在距離管道底部100~150、200~250以及300~350 mm。然而,氫含量為100%時(shí),4個(gè)工況的火焰?zhèn)鞑ニ俣扰c火焰?zhèn)鞑ゾ嚯x都基本呈線性關(guān)系,隨著火焰?zhèn)鞑ゾ嚯x的增加火焰速度越來(lái)越快。分析上述原因:火焰經(jīng)過障礙物時(shí),火焰的流通面積縮小,導(dǎo)致火焰?zhèn)鞑ニ俣仍龃?。?dāng)火焰繞過障礙物后被湍流卷吸而發(fā)生橫向拉伸,致使火焰縱向速度減小[32]。之后,在下一個(gè)障礙物的作用下,火焰再次獲得加速,如此下去,火焰速度不斷提高。但是,純氫的火焰?zhèn)鞑ニ俣确浅?欤梢詷O快地通過障礙物。

    圖7 不同工況下火焰?zhèn)鞑ゾ嚯x與火焰?zhèn)鞑ニ俣戎g的關(guān)系

    3 結(jié) 論

    (1)加氫和障礙物對(duì)最大爆炸壓力、火焰顏色和火焰?zhèn)鞑ニ俣榷紩?huì)產(chǎn)生一定的影響,隨著氫含量和障礙物數(shù)量的增加而增加,并且火焰形狀也會(huì)隨之發(fā)生變化。除此之外,純甲烷和純氫在不同工況時(shí)火焰顏色變化有明顯的區(qū)別。

    (2)當(dāng)氫含量小于50%,最大爆炸壓力隨障礙物數(shù)量的增長(zhǎng)趨勢(shì)較為微弱;當(dāng)氫含量大于50%時(shí),最大爆炸壓力隨障礙物數(shù)量的增長(zhǎng)趨勢(shì)較為明顯。同樣地,當(dāng)氫含量為0時(shí),火焰顏色隨障礙物數(shù)量的變化趨勢(shì)較為微弱,但隨著氫含量的增加障礙物對(duì)其影響則越來(lái)越明顯。以上表明障礙物對(duì)氫氣爆炸特性的影響更為明顯。除此之外,就整體而言,隨著障礙物數(shù)量的增加,最大爆炸壓力和火焰顏色隨氫含量的變化趨勢(shì)也越明顯。然而,火焰平均傳播速度隨障礙物數(shù)量的上升趨勢(shì)則較為微弱,這就表明障礙物對(duì)火焰平均傳播速度的影響較小。

    (3)障礙物對(duì)甲烷和氫氣的火焰?zhèn)鞑ニ俣扔绊懶Ч泻艽蟛煌?,甲烷在障礙物作用下會(huì)出現(xiàn)速度峰值。然而,氫氣在障礙物作用下與空管的速度變化趨勢(shì)基本一致,都隨著火焰?zhèn)鞑ゾ嚯x不斷增大。

    [1] 鄭立剛, 呂先舒, 鄭凱, 等. 點(diǎn)火源位置對(duì)甲烷/空氣爆燃超壓特征的影響 [J]. 化工學(xué)報(bào), 2015, 66 (7): 2749-2756. ZHENG L G, Lü X S, ZHENG K,. Influence of ignition position on overpressure of premixed methane-air deflagration [J]. CIESC Journal, 2015, 66 (7): 2749-2756.

    [2] 溫小萍, 武建軍, 謝茂昭. 瓦斯爆炸火焰結(jié)構(gòu)與壓力波的耦合規(guī)律 [J]. 化工學(xué)報(bào), 2013, 64 (10): 3871-3877. WEN X P, WU J J, XIE Z M. Coupled relationship between flame structure and pressure wave of gas explosion [J]. CIESC Journal, 2013, 64 (10): 3871-3877.

    [3] 李增華, 林柏泉, 張?zhí)m君, 等. 氫氣的生成及對(duì)瓦斯爆炸的影響 [J]. 中國(guó)礦業(yè)大學(xué)學(xué)報(bào), 2008, 37 (2): 147-151. LI Z H, LIN B Q, ZHANG L J,. Effects of hydrogen production on gas explosion [J]. Journal of China University of Ming &Technology, 2008, 37 (2): 147-151.

    [4] OKAFOR E C, HAYAKAWA A Y, KITAGAWA T. Effects of hydrogen concentration on premixed laminar flames of hydrogen-methane-air [J]. International Journal of Hydrogen Energy, 2014, 39 (24): 9-17.

    [5] CHEN Z, DAI P, CHEN S Y. A model for the laminar flame speed of binary fuel blends and its application to methane/hydrogen mixtures [J]. International Journal of Hydrogen Energy, 2012, 37 (103): 90-96.

    [6] DI SARLIA V, DI BENEDETTOB A. Laminar burning velocity of hydrogen-methane/air premixed flames [J]. International Journal of Hydrogen Energy, 2007, 32 (6): 37-46.

    [7] HUANG Z H, ZHANG Y, ZENG K,. Measurements of laminar burning velocities for natural gas-hydrogen-air mixtures [J]. Combustion Flame, 2006, 146 (3): 2-11.

    [8] MA Q J, ZHANG Q, CHEN J,. Effects of hydrogen on combustion characteristics of methane in air [J]. International Journal of Hydrogen Energy, 2014, 39: 11291-11298.

    [9] MA Q J, ZHANG Q, PANG L. Effects of hydrogen addition on the confined and vented explosion behavior of methane in air [J]. Journal of Loss Prevention in the Process Industries, 2014, 27: 65-73.

    [10] YU M G, ZHENG K, ZHENG L G. Effects of hydrogen addition on propagation characteristics of premixed methane/air flames [J]. Journal of Loss Prevention in the Process Industries, 2015, 34: 1-9.

    [11] 賈寶山, 溫海燕, 梁運(yùn)濤. 受限空間瓦斯爆炸與氫氣促進(jìn)機(jī)理研究 [J]. 中國(guó)安全科學(xué)學(xué)報(bào), 2012, 22 (2): 81-87. JIA B S, WEN H Y, LIANG Y T. Study on the methane explosion in an enclosed space and hydrogen promoting mechanism [J]. China Safety Science Journal, 2012, 22 (2): 81-87.

    [12] POROWSKI R, TEODORCZYK A. Experimental study on DDT for hydrogen-methane-air mixtures in tube with obstacles [J]. Journal of Loss Prevention in the Process Industries, 2013, 26: 374-375.

    [13] SALZANO E, CAMMAROTA F, DI BENEDETTO A,. Explosion behavior of hydrogen-methane/air mixtures [J]. Journal of Loss Prevention in the Process Industries, 2012, 25: 443-447.

    [14] WOOLLEY R M, FAIRWEATHER M, FALLE S A E G,. Prediction of confined, vented methane-hydrogen explosion using a computational fluid dynamic approach [J]. International Journal of Hydrogen Energy, 2013, 38: 6904-6914.

    [15] 潘鵬飛, 譚迎新, 王志青. 條形障礙物對(duì)瓦斯爆炸特性影響研究 [J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù), 2010, 6 (5): 71-76. PAN P F, TAN Y Z, WANG Z Q. Study on effects of strip-type obstacles on methane explosion [J]. Journal of Safety Science and Technology, 2010, 6 (5): 71-76.

    [16] 何學(xué)秋, 楊藝, 王恩元, 等. 障礙物對(duì)瓦斯爆炸火焰結(jié)構(gòu)及火焰?zhèn)鞑ビ绊懙难芯?[J]. 煤炭學(xué)報(bào), 2004, 2 (2): 186-189. HE X Q, YANG Y, WANG E Y,. Effects of obstacles on premixed flame microstructure and flame propagation in methane/air explosion [J]. Journal of China Coal Society, 2004, 2 (2): 186-189.

    [17] ALHARBI A, MASRI A R, IBRHAIM S. Turbulent premixed flames of CNG, LPG, and H2propagating past repeated obstacles [J]. Experimental Thermal and Fluid Science, 2014, 56: 2-8.

    [18] IBRHAIM S, MASRI A R. The effect of obstructions on overpressure resulting from premixed flame deflagration [J]. Journal of Loss Prevention in the Process Industries, 2001, 14: 213-221.

    [19] 王海賓, 尉存娟, 譚迎新. 水平管道內(nèi)障礙物對(duì)爆炸壓力影響的研究 [J].廣州化工, 2014, 39 (24): 43-46. WANG H B, WEI C J, TAN Y X. Effects of pressure of gas explosion in the pipeline of obstacles [J]. Guangzhou Chemical Industry, 2014, 39 (24): 43-46.

    [20] 余明高, 紀(jì)文濤, 溫小萍. 交錯(cuò)障礙物對(duì)瓦斯爆炸的實(shí)驗(yàn)研究 [J]. 中國(guó)礦業(yè)大學(xué)學(xué)報(bào), 2013, 42 (3): 349-354. YU M G, JI W T, WEN X P. Experimental study of the influence of staggered obstacles on gas explosion [J]. Journal of China University of MingTechnology, 2013, 42 (3): 349-354.

    [21] 林柏泉, 周世寧, 張仁貴. 障礙物對(duì)瓦斯爆炸過程中火焰和爆炸波的影響 [J]. 中國(guó)礦業(yè)大學(xué)學(xué)報(bào), 1999, 28 (2): 104-107. LIN B Q, ZHOU S N, ZHANG R G. Influence of barriers on flame transmission and explosion wave in gas explosion [J]. Journal of China University of Ming & Technology, 1999, 28 (2): 104-107.

    [22] 尉存娟, 譚迎新, 袁宏甦.水平管道內(nèi)障礙物數(shù)量對(duì)爆炸過程的影響研究 [J]. 中國(guó)安全科學(xué)學(xué)報(bào), 2012, 22 (6): 60-64. WEI C J, TAN Y X, YUAN H S. Study on the influence of obstacles quantity on gas explosion process in horizontal pipe device [J]. China Safety Science Journal, 2012, 22 (6): 60-64.

    [23] HALL S, MASRI A R, YAROSHCHYK P,. Effects of position and frequency of obstacles on turbulent premixed propagating flames [J]. Combustion and Flame, 2009, 156 (8): 439-446.

    [24] GUBBA S, IBRAHIM S, MALALASEKERA W,. Measures and LES calculations of turbulent premixed flame propagation past repeated obstacles [J]. Combustion and Flame, 2011, 158 (12): 2465-2481.

    [25] PARK D, LEE Y S, GREEN A R. Experiments on the effects of multiple obstacles invented explosion chambers [J]. Journal of Hazardous Materials, 2008, 153 (5): 340-350.

    [26] 郭丹彤, 呂淑然, 楊凱. 障礙物布置對(duì)氣體爆炸壓力場(chǎng)的影響效果研究 [J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù), 2015, 11 (9): 88-93. GUO D T, Lü S R, YANG K. Research on impact effect of obstacle arrangement to pressure filed of gas explosion [J]. Journal of Safety Science and Technology, 2015, 11 (9): 88-93.

    [27] 宋文化, 李小偉, 李冬梅. 道化學(xué)火災(zāi)爆炸指數(shù)評(píng)價(jià)法在合成氨轉(zhuǎn)化工序安全性評(píng)價(jià)中的應(yīng)用 [J]. 消防科學(xué)與技術(shù), 2008, 27 (5): 321-324. SONG W H, LI X W, LI D M. Application of Dow Chemical Fire and explosion index evaluation method in safety evaluation of synthetic ammonia conversion process [J]. Fire Science and Technology, 2008, 27 (5): 321-324.

    [28] 藺照東, 李如江, 劉恩良, 等. 障礙物對(duì)瓦斯爆炸沖擊波影響研究 [J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù), 2012, 8 (2): 28-32.LIN Z D, LI R J, LIU E L,. Study on effect of obstacles to shock wave of gas explosion [J]. Journal of Safety Science and Technology, 2012, 8 (2): 28-32.

    [29] WEN X P, XIE M Z, YU M G,. Porous media quenching behaviors of gas deflagration in the presence of obstacles [J]. Experimental Thermal and Fluid Science, 2013, 50 (6): 37-44.

    [30] 林柏泉, 葉青, 翟成, 等. 瓦斯爆炸在分叉管道中的傳播規(guī)律及分析 [J]. 煤炭學(xué)報(bào), 2008, 33 (2): 136-139. LIN B Q, YE Q, ZHAI C,. The propagation rule of methane explosion in bifurcation duct [J]. Journal of China Coal Society, 2008, 33 (2): 136-139.

    [31] 葉青, 賈真真, 林柏泉. 管內(nèi)瓦斯爆炸的火焰加速機(jī)理分析 [J]. 煤礦安全, 2008, (1): 78-80. YE Q, JIA Z Z, LIN B Q. Analysis on flame acceleration mechanism of gas explosion in tube [J]. Coal Mine Safety,2008, (1): 78-80.

    [32] 溫小萍. 瓦斯湍流爆燃火焰特性與多孔介質(zhì)淬熄抑爆機(jī)理研究[D]. 大連: 大連理工大學(xué), 2013. WEN X P. Mechanism study on flame characteristics and porous media quenching suppression of gas turbulent deflagration [D]. Dalian: Dalian University of Technology, 2013.

    Effects of hydrogen addition on explosion characteristics of gas under condition of obstacles

    YU Minggao1,2, YUAN Chenqiao1, ZHENG Kai2

    (1School of Safety Science Engineering, Henan Polytechnic University, Jiaozuo 454003, Henan, China;2State Key Laboratory of Coal Mine Disaster Dynamics Control, Chongqing University, Chongqing 400044, China)

    The effects of hydrogen addition on explosion characteristics of gas under the condition of obstacles were investigated. The hydrogen fraction in the methane-hydrogen mixture was varied from 0 to 1 at equivalence ratio of 1. The result indicated that the maximum pressure and propagation speed of flame both increased with increasing fraction of hydrogen and obstacle quantity. In addition, the hydrogen addition can also change flame color and flame front. However, the influence of the obstacle quantity on average propagation velocity of flame was slight, and the effects of the obstacle quantity on flame propagation speed of methane as well as hydrogen were different. Increasing trend of the maximum pressure and flame color with the fraction of hydrogen was more and more obvious as the obstacle quantity increased. The influence of the obstacle quantity on explosion characteristics of hydrogen was more obvious than that of methane.

    hydrogen; obstacle; methane/air; explosion characteristics

    date: 2016-05-11.

    Prof. YU Minggao, mg_yu@126.com

    10.11949/j.issn.0438-1157.20160645

    TD 712

    A

    0438—1157(2016)12—5311—09

    國(guó)家自然科學(xué)基金項(xiàng)目(U1361205,51574111);煤礦災(zāi)害動(dòng)力學(xué)與控制國(guó)家重點(diǎn)實(shí)驗(yàn)室自主研究課題重點(diǎn)項(xiàng)目(2011DA 105287*ZD201401)。

    supported by the Natural National Science Foundation of China (U1361205, 51574111).

    2016-05-11收到初稿,2016-09-05收到修改稿。

    聯(lián)系人及第一作者:余明高(1963—),男,教授。

    猜你喜歡
    傳播速度障礙物氫氣
    代謝綜合征患者臂踝脈搏波傳播速度與頸動(dòng)脈粥樣硬化的關(guān)系
    高低翻越
    SelTrac?CBTC系統(tǒng)中非通信障礙物的設(shè)計(jì)和處理
    新雷
    一類廣義canmassa—Holm方程的無(wú)限傳播速度與漸近行為
    一種用氫氣還原制備石墨烯的方法
    氫氣對(duì)缺血再灌注損傷保護(hù)的可能機(jī)制
    小天使·六年級(jí)語(yǔ)數(shù)英綜合(2015年2期)2015-01-22 04:50:59
    大型氫氣隔膜壓縮機(jī)的開發(fā)及應(yīng)用
    土釘墻在近障礙物的地下車行通道工程中的應(yīng)用
    亚洲熟女毛片儿| 99久久精品国产亚洲精品| 亚洲黑人精品在线| 一夜夜www| 俄罗斯特黄特色一大片| 国产日韩欧美在线精品| 国产精品影院久久| 成人国产一区最新在线观看| 中文字幕制服av| 不卡一级毛片| 999久久久精品免费观看国产| 亚洲成人免费电影在线观看| 飞空精品影院首页| 女同久久另类99精品国产91| 精品国产亚洲在线| 女性被躁到高潮视频| 香蕉久久夜色| 国产亚洲精品一区二区www | 久久性视频一级片| 色精品久久人妻99蜜桃| 亚洲七黄色美女视频| 男女之事视频高清在线观看| 成人亚洲精品一区在线观看| 日本精品一区二区三区蜜桃| 色老头精品视频在线观看| 色综合婷婷激情| 淫妇啪啪啪对白视频| 麻豆国产av国片精品| 国产成人精品在线电影| 国产精品国产av在线观看| 国产精品香港三级国产av潘金莲| 亚洲美女黄片视频| 侵犯人妻中文字幕一二三四区| 成年版毛片免费区| 欧美精品啪啪一区二区三区| 一区二区三区乱码不卡18| 午夜免费鲁丝| 青青草视频在线视频观看| av一本久久久久| 国产亚洲午夜精品一区二区久久| 成年女人毛片免费观看观看9 | 欧美日韩中文字幕国产精品一区二区三区 | 黄频高清免费视频| 十八禁网站网址无遮挡| 深夜精品福利| 国产精品二区激情视频| 午夜福利一区二区在线看| 国产av一区二区精品久久| 一区二区日韩欧美中文字幕| 国产男靠女视频免费网站| 精品国产超薄肉色丝袜足j| 丰满饥渴人妻一区二区三| 亚洲精品成人av观看孕妇| 国产精品免费视频内射| 精品福利观看| 国产精品久久久久久精品电影小说| 岛国在线观看网站| 一边摸一边做爽爽视频免费| 激情在线观看视频在线高清 | 国产老妇伦熟女老妇高清| 国产亚洲精品一区二区www | www日本在线高清视频| 热99国产精品久久久久久7| 老司机在亚洲福利影院| 亚洲五月色婷婷综合| 久久国产精品人妻蜜桃| 黄色视频,在线免费观看| 亚洲欧美激情在线| 国产成人av激情在线播放| 大香蕉久久成人网| 热99久久久久精品小说推荐| 亚洲欧美日韩另类电影网站| 欧美中文综合在线视频| 久久国产精品大桥未久av| 黄色 视频免费看| 久久人妻av系列| 最近最新免费中文字幕在线| 亚洲美女黄片视频| 777久久人妻少妇嫩草av网站| 亚洲欧洲精品一区二区精品久久久| av网站免费在线观看视频| 欧美在线黄色| 香蕉久久夜色| 国产精品电影一区二区三区 | 国产aⅴ精品一区二区三区波| 国产日韩一区二区三区精品不卡| 丰满人妻熟妇乱又伦精品不卡| 午夜成年电影在线免费观看| 国产成人免费观看mmmm| 国产精品成人在线| 久久九九热精品免费| 欧美精品一区二区免费开放| 一区在线观看完整版| 18禁美女被吸乳视频| 成在线人永久免费视频| 一个人免费在线观看的高清视频| 久久毛片免费看一区二区三区| 天天影视国产精品| 啦啦啦免费观看视频1| 啦啦啦免费观看视频1| 国产精品香港三级国产av潘金莲| 国产在线观看jvid| 欧美午夜高清在线| 亚洲精品一二三| 久久毛片免费看一区二区三区| h视频一区二区三区| 精品人妻在线不人妻| 亚洲成av片中文字幕在线观看| 久久av网站| 日韩熟女老妇一区二区性免费视频| 在线永久观看黄色视频| kizo精华| 深夜精品福利| 日韩熟女老妇一区二区性免费视频| 99精国产麻豆久久婷婷| 性色av乱码一区二区三区2| 啪啪无遮挡十八禁网站| 18在线观看网站| 两性夫妻黄色片| 18禁黄网站禁片午夜丰满| 视频区图区小说| 桃花免费在线播放| 老司机福利观看| 国产成人啪精品午夜网站| 欧美大码av| 黄色毛片三级朝国网站| 99久久99久久久精品蜜桃| 超碰成人久久| 九色亚洲精品在线播放| 999久久久精品免费观看国产| 日韩欧美国产一区二区入口| 在线观看人妻少妇| 亚洲国产看品久久| 国产一卡二卡三卡精品| 欧美激情 高清一区二区三区| 欧美激情极品国产一区二区三区| www.自偷自拍.com| 激情在线观看视频在线高清 | 99riav亚洲国产免费| 露出奶头的视频| 久久精品熟女亚洲av麻豆精品| 中文亚洲av片在线观看爽 | 18禁黄网站禁片午夜丰满| 国产伦人伦偷精品视频| 在线永久观看黄色视频| 国产高清视频在线播放一区| 99国产精品一区二区三区| 国产99久久九九免费精品| 亚洲国产av影院在线观看| 99国产精品一区二区蜜桃av | 成年人免费黄色播放视频| 亚洲av日韩在线播放| 麻豆国产av国片精品| 色综合欧美亚洲国产小说| 中文字幕色久视频| 两个人看的免费小视频| 国产主播在线观看一区二区| 精品一区二区三卡| 一本一本久久a久久精品综合妖精| 免费高清在线观看日韩| 天天操日日干夜夜撸| tocl精华| 亚洲国产欧美在线一区| 欧美国产精品一级二级三级| 一级毛片女人18水好多| 亚洲专区国产一区二区| 午夜精品国产一区二区电影| 少妇精品久久久久久久| 国产一卡二卡三卡精品| 91av网站免费观看| 免费人妻精品一区二区三区视频| 99热网站在线观看| www.熟女人妻精品国产| 王馨瑶露胸无遮挡在线观看| 亚洲一区中文字幕在线| 飞空精品影院首页| avwww免费| 18禁国产床啪视频网站| 国产在线免费精品| 亚洲精品自拍成人| 日本av免费视频播放| 久久久久视频综合| 国产aⅴ精品一区二区三区波| 脱女人内裤的视频| 美女福利国产在线| 中文字幕精品免费在线观看视频| 一边摸一边抽搐一进一出视频| 国产精品1区2区在线观看. | 欧美国产精品一级二级三级| 黄色怎么调成土黄色| 中文字幕精品免费在线观看视频| 久久久久久免费高清国产稀缺| 韩国精品一区二区三区| 两个人免费观看高清视频| 99在线人妻在线中文字幕 | 免费观看av网站的网址| 久久久久久久久久久久大奶| 久久午夜亚洲精品久久| 精品久久久久久电影网| 日韩有码中文字幕| 亚洲少妇的诱惑av| 搡老岳熟女国产| 亚洲国产毛片av蜜桃av| 黄色视频不卡| 一个人免费在线观看的高清视频| 黑丝袜美女国产一区| 国产精品久久久久久人妻精品电影 | 欧美国产精品一级二级三级| 高清黄色对白视频在线免费看| 欧美人与性动交α欧美精品济南到| 无人区码免费观看不卡 | 国产国语露脸激情在线看| 久久亚洲真实| 性少妇av在线| 一级片'在线观看视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成人免费av在线播放| 欧美在线黄色| 99re6热这里在线精品视频| 亚洲黑人精品在线| tocl精华| 国产一区二区在线观看av| 正在播放国产对白刺激| 黄色毛片三级朝国网站| av国产精品久久久久影院| 色播在线永久视频| 国产在视频线精品| 国产精品国产高清国产av | 激情在线观看视频在线高清 | www.熟女人妻精品国产| netflix在线观看网站| 国产精品久久久久成人av| 成人18禁在线播放| 精品一区二区三区av网在线观看 | 亚洲中文字幕日韩| 国产高清视频在线播放一区| 国产深夜福利视频在线观看| 宅男免费午夜| a级片在线免费高清观看视频| 成年动漫av网址| 亚洲人成77777在线视频| 性高湖久久久久久久久免费观看| 男女午夜视频在线观看| 欧美日韩亚洲高清精品| 丝瓜视频免费看黄片| 欧美黑人精品巨大| 久久久欧美国产精品| 欧美老熟妇乱子伦牲交| 国产一区二区三区在线臀色熟女 | 国产激情久久老熟女| 性少妇av在线| 免费观看a级毛片全部| 国产麻豆69| 淫妇啪啪啪对白视频| 又黄又粗又硬又大视频| 国产精品免费视频内射| 久久婷婷成人综合色麻豆| 国产精品香港三级国产av潘金莲| 一本一本久久a久久精品综合妖精| 国产一区二区激情短视频| 国产免费福利视频在线观看| 欧美日韩一级在线毛片| 18禁裸乳无遮挡动漫免费视频| 夫妻午夜视频| 国产精品电影一区二区三区 | √禁漫天堂资源中文www| 成人18禁在线播放| 精品国内亚洲2022精品成人 | 日韩三级视频一区二区三区| 亚洲黑人精品在线| 叶爱在线成人免费视频播放| 俄罗斯特黄特色一大片| 中亚洲国语对白在线视频| 一区在线观看完整版| 十分钟在线观看高清视频www| 亚洲第一欧美日韩一区二区三区 | 午夜福利影视在线免费观看| 欧美人与性动交α欧美精品济南到| 精品久久蜜臀av无| 国产一区二区 视频在线| 国产成人精品久久二区二区免费| 欧美日韩黄片免| 欧美+亚洲+日韩+国产| 久久精品熟女亚洲av麻豆精品| 久久中文字幕一级| 亚洲色图 男人天堂 中文字幕| 亚洲精品乱久久久久久| 午夜免费鲁丝| 中文亚洲av片在线观看爽 | 9色porny在线观看| 一本大道久久a久久精品| 国产亚洲精品第一综合不卡| 十八禁网站网址无遮挡| 亚洲黑人精品在线| 久久久精品国产亚洲av高清涩受| 精品国内亚洲2022精品成人 | 国产精品免费一区二区三区在线 | 亚洲av欧美aⅴ国产| 女警被强在线播放| 久久影院123| 亚洲av成人一区二区三| 亚洲熟女毛片儿| 国产精品国产av在线观看| videosex国产| 纵有疾风起免费观看全集完整版| 亚洲第一欧美日韩一区二区三区 | 啦啦啦视频在线资源免费观看| 一边摸一边抽搐一进一出视频| 丝袜喷水一区| 亚洲中文av在线| 桃红色精品国产亚洲av| 免费av中文字幕在线| 看免费av毛片| 一级黄色大片毛片| 国产在线观看jvid| 亚洲九九香蕉| 欧美日韩福利视频一区二区| 伦理电影免费视频| 另类精品久久| 国产在线观看jvid| 大香蕉久久网| 午夜福利在线免费观看网站| a在线观看视频网站| 亚洲第一欧美日韩一区二区三区 | 国产精品久久久久成人av| 欧美 日韩 精品 国产| 亚洲第一青青草原| 最近最新免费中文字幕在线| 成年女人毛片免费观看观看9 | 精品久久久久久久毛片微露脸| 国产欧美日韩一区二区精品| 性少妇av在线| netflix在线观看网站| 成人av一区二区三区在线看| avwww免费| 啪啪无遮挡十八禁网站| 91麻豆av在线| 高清欧美精品videossex| 亚洲美女黄片视频| 久久久国产精品麻豆| 99热国产这里只有精品6| 窝窝影院91人妻| 在线观看免费视频日本深夜| 精品一品国产午夜福利视频| 69av精品久久久久久 | 91字幕亚洲| 男女无遮挡免费网站观看| 中文字幕另类日韩欧美亚洲嫩草| 老司机午夜十八禁免费视频| 久久国产精品人妻蜜桃| 侵犯人妻中文字幕一二三四区| 亚洲av成人不卡在线观看播放网| 黑丝袜美女国产一区| 男女之事视频高清在线观看| 国产欧美日韩综合在线一区二区| 国产精品欧美亚洲77777| 我要看黄色一级片免费的| 一个人免费看片子| 国产一区二区激情短视频| 亚洲熟女精品中文字幕| 色综合婷婷激情| 黑人巨大精品欧美一区二区蜜桃| 日本精品一区二区三区蜜桃| 黄片大片在线免费观看| 久久久国产一区二区| netflix在线观看网站| 啦啦啦免费观看视频1| 91av网站免费观看| 人人妻人人澡人人看| 精品卡一卡二卡四卡免费| 国产男靠女视频免费网站| 成年动漫av网址| 国产xxxxx性猛交| 免费一级毛片在线播放高清视频 | 国产欧美日韩精品亚洲av| 国产野战对白在线观看| 两个人免费观看高清视频| 亚洲九九香蕉| 精品国产一区二区久久| 激情在线观看视频在线高清 | 搡老熟女国产l中国老女人| 老司机亚洲免费影院| 新久久久久国产一级毛片| 欧美精品人与动牲交sv欧美| 女警被强在线播放| 50天的宝宝边吃奶边哭怎么回事| 91字幕亚洲| 美女高潮到喷水免费观看| 午夜福利在线观看吧| 欧美精品av麻豆av| 一二三四在线观看免费中文在| 国产精品麻豆人妻色哟哟久久| 亚洲精品中文字幕一二三四区 | 精品国内亚洲2022精品成人 | av视频免费观看在线观看| 91麻豆精品激情在线观看国产 | √禁漫天堂资源中文www| 亚洲精品久久午夜乱码| 午夜激情av网站| 99在线人妻在线中文字幕 | 啦啦啦 在线观看视频| 欧美乱妇无乱码| 丰满少妇做爰视频| 免费少妇av软件| 狠狠狠狠99中文字幕| 91大片在线观看| 麻豆成人av在线观看| 啦啦啦免费观看视频1| 精品国产一区二区三区久久久樱花| 亚洲天堂av无毛| 国产精品一区二区免费欧美| 妹子高潮喷水视频| 欧美黄色片欧美黄色片| 国产精品偷伦视频观看了| 日本av手机在线免费观看| 日韩制服丝袜自拍偷拍| 五月开心婷婷网| 自拍欧美九色日韩亚洲蝌蚪91| 丰满少妇做爰视频| 乱人伦中国视频| 一级片免费观看大全| 国产精品 欧美亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 十八禁网站网址无遮挡| 久久精品国产亚洲av高清一级| 亚洲第一欧美日韩一区二区三区 | 精品高清国产在线一区| 免费黄频网站在线观看国产| 国产有黄有色有爽视频| 国产欧美日韩综合在线一区二区| av不卡在线播放| 久久久欧美国产精品| 高清毛片免费观看视频网站 | av一本久久久久| 18禁观看日本| 亚洲精品粉嫩美女一区| 亚洲国产av新网站| 国产精品1区2区在线观看. | 手机成人av网站| 成人永久免费在线观看视频 | 国产成人免费无遮挡视频| 国产成+人综合+亚洲专区| 国产日韩欧美亚洲二区| 亚洲 国产 在线| 久久国产精品男人的天堂亚洲| 精品国产一区二区三区久久久樱花| 电影成人av| 18禁黄网站禁片午夜丰满| 成人18禁高潮啪啪吃奶动态图| 91精品国产国语对白视频| 手机成人av网站| 欧美激情极品国产一区二区三区| 无人区码免费观看不卡 | 久久亚洲精品不卡| a在线观看视频网站| 在线观看舔阴道视频| 国产淫语在线视频| 最新的欧美精品一区二区| 久久影院123| 女人被躁到高潮嗷嗷叫费观| 亚洲国产欧美一区二区综合| 黄色怎么调成土黄色| 久久久久网色| 久久久国产精品麻豆| 午夜久久久在线观看| 欧美日韩精品网址| 性高湖久久久久久久久免费观看| 欧美成人免费av一区二区三区 | 一级片'在线观看视频| 国产深夜福利视频在线观看| 免费观看av网站的网址| 国产男女内射视频| 人人妻,人人澡人人爽秒播| 亚洲精品自拍成人| 中文字幕精品免费在线观看视频| 久久精品熟女亚洲av麻豆精品| 国产黄色免费在线视频| 日日夜夜操网爽| 国产野战对白在线观看| 51午夜福利影视在线观看| 久久ye,这里只有精品| 精品国产超薄肉色丝袜足j| 色94色欧美一区二区| 午夜激情久久久久久久| 亚洲伊人色综图| 久久久久久久久久久久大奶| 欧美乱码精品一区二区三区| 国产在线观看jvid| 色精品久久人妻99蜜桃| 国产精品久久久久成人av| 一二三四社区在线视频社区8| 多毛熟女@视频| 国产在线一区二区三区精| 中文字幕人妻丝袜制服| 大片免费播放器 马上看| 欧美变态另类bdsm刘玥| 大片电影免费在线观看免费| 美女视频免费永久观看网站| 国产91精品成人一区二区三区 | 亚洲国产欧美在线一区| 天堂8中文在线网| 亚洲专区中文字幕在线| 久久亚洲真实| 成人精品一区二区免费| 丝袜在线中文字幕| 丝袜人妻中文字幕| av一本久久久久| 男人舔女人的私密视频| 中文欧美无线码| 精品国产国语对白av| 国产欧美亚洲国产| 法律面前人人平等表现在哪些方面| 色在线成人网| 一区在线观看完整版| www日本在线高清视频| 露出奶头的视频| 淫妇啪啪啪对白视频| 国产伦理片在线播放av一区| 日本vs欧美在线观看视频| 丝袜在线中文字幕| 精品免费久久久久久久清纯 | 亚洲第一欧美日韩一区二区三区 | 一本久久精品| 丝袜喷水一区| 免费观看a级毛片全部| 国产老妇伦熟女老妇高清| 一夜夜www| 国产av国产精品国产| 高清视频免费观看一区二区| 丝袜美腿诱惑在线| 中文字幕最新亚洲高清| 精品亚洲成国产av| 女人精品久久久久毛片| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 色尼玛亚洲综合影院| 久久久久久久国产电影| 亚洲一卡2卡3卡4卡5卡精品中文| 美女国产高潮福利片在线看| 黄频高清免费视频| 精品高清国产在线一区| 国产在线免费精品| 亚洲情色 制服丝袜| 日本精品一区二区三区蜜桃| 日韩中文字幕视频在线看片| 精品人妻1区二区| 男人操女人黄网站| 一二三四在线观看免费中文在| 电影成人av| 国产亚洲精品一区二区www | 久久久久精品国产欧美久久久| 韩国精品一区二区三区| 国产精品美女特级片免费视频播放器 | 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av高清一级| 亚洲av片天天在线观看| 亚洲男人天堂网一区| 我的亚洲天堂| 国产91精品成人一区二区三区 | 我的亚洲天堂| 欧美日本中文国产一区发布| 免费不卡黄色视频| 97在线人人人人妻| 久久香蕉激情| 麻豆乱淫一区二区| 黄片大片在线免费观看| 五月天丁香电影| 国产成人精品久久二区二区91| 中文字幕精品免费在线观看视频| 久久人妻熟女aⅴ| 欧美日韩国产mv在线观看视频| 国产一区二区三区在线臀色熟女 | 美女高潮喷水抽搐中文字幕| 国产精品欧美亚洲77777| 亚洲,欧美精品.| 久久久久精品国产欧美久久久| kizo精华| 久久久精品94久久精品| 国产精品国产av在线观看| 国产成人影院久久av| 不卡av一区二区三区| 老鸭窝网址在线观看| 99国产精品一区二区三区| cao死你这个sao货| 国产亚洲精品一区二区www | 90打野战视频偷拍视频| 脱女人内裤的视频| videosex国产| 国精品久久久久久国模美| 成人免费观看视频高清| 精品国产乱子伦一区二区三区| 亚洲欧美色中文字幕在线| 两人在一起打扑克的视频| 国产成+人综合+亚洲专区| 国产亚洲精品久久久久5区| 精品国产国语对白av| 黄频高清免费视频| 超碰成人久久| av超薄肉色丝袜交足视频| 9热在线视频观看99| 欧美精品人与动牲交sv欧美| 不卡av一区二区三区| 国产精品 国内视频| av免费在线观看网站| 免费人妻精品一区二区三区视频| 亚洲精品在线美女| 黄色片一级片一级黄色片| 亚洲av成人不卡在线观看播放网| 日韩大片免费观看网站| 精品亚洲乱码少妇综合久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久狼人影院| 最近最新中文字幕大全电影3 | 黄频高清免费视频|