• <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)用
    亚洲一级一片aⅴ在线观看| 高清av免费在线| 国产视频首页在线观看| 我要看黄色一级片免费的| 另类精品久久| 精品人妻熟女毛片av久久网站| 中文字幕人妻丝袜制服| 日韩三级伦理在线观看| av卡一久久| 蜜桃国产av成人99| 国产在线免费精品| 五月伊人婷婷丁香| 久久97久久精品| 激情视频va一区二区三区| 日韩一区二区视频免费看| 亚洲第一区二区三区不卡| 亚洲美女视频黄频| av国产久精品久网站免费入址| 香蕉精品网在线| 久久亚洲国产成人精品v| 电影成人av| 熟女少妇亚洲综合色aaa.| 王馨瑶露胸无遮挡在线观看| 国产野战对白在线观看| 在线观看免费视频网站a站| 亚洲精品美女久久久久99蜜臀 | 亚洲av综合色区一区| 狠狠精品人妻久久久久久综合| 国产免费又黄又爽又色| 国产精品二区激情视频| 69精品国产乱码久久久| 少妇猛男粗大的猛烈进出视频| 一区二区三区激情视频| 久久鲁丝午夜福利片| 黄片无遮挡物在线观看| 久久97久久精品| 女人被躁到高潮嗷嗷叫费观| 亚洲成色77777| 如日韩欧美国产精品一区二区三区| 国产视频首页在线观看| 国产精品免费视频内射| 日韩电影二区| 男女免费视频国产| 日韩三级伦理在线观看| 2021少妇久久久久久久久久久| 高清在线视频一区二区三区| 成人亚洲精品一区在线观看| 又大又黄又爽视频免费| 久久99热这里只频精品6学生| 热re99久久精品国产66热6| 18禁国产床啪视频网站| 91午夜精品亚洲一区二区三区| 国产精品三级大全| 欧美日韩av久久| 精品视频人人做人人爽| 满18在线观看网站| 高清av免费在线| 下体分泌物呈黄色| 丝袜喷水一区| av.在线天堂| 精品一区二区免费观看| 精品少妇内射三级| 自线自在国产av| 国产在线免费精品| freevideosex欧美| av一本久久久久| 97精品久久久久久久久久精品| 久久婷婷青草| 中文天堂在线官网| 国产男人的电影天堂91| videossex国产| 国产野战对白在线观看| 免费女性裸体啪啪无遮挡网站| 日韩三级伦理在线观看| 在现免费观看毛片| 亚洲第一av免费看| 满18在线观看网站| 国产成人精品婷婷| 丝袜脚勾引网站| 秋霞伦理黄片| 一级片免费观看大全| 欧美中文综合在线视频| 少妇人妻精品综合一区二区| 久久久欧美国产精品| 美女大奶头黄色视频| 久久久精品94久久精品| 少妇人妻精品综合一区二区| kizo精华| 欧美精品一区二区免费开放| 日韩成人av中文字幕在线观看| 18禁动态无遮挡网站| 国产一区二区激情短视频 | 久久久久久久亚洲中文字幕| 黄频高清免费视频| 亚洲av福利一区| 曰老女人黄片| 久久精品aⅴ一区二区三区四区 | 色婷婷久久久亚洲欧美| 精品一区二区三卡| 成人毛片a级毛片在线播放| 日本猛色少妇xxxxx猛交久久| 欧美日韩国产mv在线观看视频| 最近的中文字幕免费完整| 人妻系列 视频| 妹子高潮喷水视频| 亚洲国产毛片av蜜桃av| 免费在线观看完整版高清| 99精国产麻豆久久婷婷| 一级片免费观看大全| 波多野结衣一区麻豆| 久久久久视频综合| 搡老乐熟女国产| 亚洲国产欧美网| 国产精品国产三级国产专区5o| 久久精品国产鲁丝片午夜精品| 一区二区日韩欧美中文字幕| 90打野战视频偷拍视频| 麻豆av在线久日| 日韩一卡2卡3卡4卡2021年| 9色porny在线观看| 亚洲国产精品成人久久小说| 亚洲精品一区蜜桃| 你懂的网址亚洲精品在线观看| 午夜激情久久久久久久| 97精品久久久久久久久久精品| 高清视频免费观看一区二区| 美女主播在线视频| 性色avwww在线观看| 成年女人在线观看亚洲视频| 亚洲av国产av综合av卡| 色视频在线一区二区三区| 男人舔女人的私密视频| 国产成人精品在线电影| 日韩三级伦理在线观看| 免费黄色在线免费观看| 深夜精品福利| 日韩伦理黄色片| 咕卡用的链子| 精品一区二区三卡| 日产精品乱码卡一卡2卡三| 最近中文字幕高清免费大全6| 亚洲图色成人| 男女免费视频国产| 国产高清国产精品国产三级| 超碰97精品在线观看| 亚洲成色77777| 又粗又硬又长又爽又黄的视频| 一级黄片播放器| 91精品三级在线观看| 亚洲伊人久久精品综合| 国产老妇伦熟女老妇高清| kizo精华| 国产极品粉嫩免费观看在线| 好男人视频免费观看在线| 青春草视频在线免费观看| 免费人妻精品一区二区三区视频| 亚洲精品av麻豆狂野| 精品人妻一区二区三区麻豆| 天天操日日干夜夜撸| 成人漫画全彩无遮挡| 欧美成人午夜免费资源| 久久99一区二区三区| 欧美精品av麻豆av| 麻豆av在线久日| 亚洲情色 制服丝袜| 天天躁夜夜躁狠狠躁躁| av.在线天堂| 欧美精品av麻豆av| 天天躁夜夜躁狠狠躁躁| 男女免费视频国产| 黄片播放在线免费| 日日啪夜夜爽| av网站免费在线观看视频| 久久久久久久久免费视频了| 韩国精品一区二区三区| 熟女电影av网| 丰满乱子伦码专区| 国产亚洲av片在线观看秒播厂| 麻豆精品久久久久久蜜桃| 蜜桃国产av成人99| 久久久久久久亚洲中文字幕| 亚洲,欧美,日韩| 午夜福利网站1000一区二区三区| 国产 一区精品| 看十八女毛片水多多多| 欧美精品人与动牲交sv欧美| 国产精品欧美亚洲77777| 久久久精品区二区三区| 九草在线视频观看| 成年女人毛片免费观看观看9 | 中文字幕另类日韩欧美亚洲嫩草| 熟女少妇亚洲综合色aaa.| 啦啦啦啦在线视频资源| 久久久国产精品麻豆| 777米奇影视久久| 少妇人妻精品综合一区二区| videosex国产| 亚洲精品国产一区二区精华液| 日韩中文字幕欧美一区二区 | 最黄视频免费看| 久久精品国产综合久久久| 亚洲国产看品久久| 国产精品无大码| av电影中文网址| 18+在线观看网站| 久久久久精品久久久久真实原创| 亚洲四区av| 久久亚洲国产成人精品v| 国产精品秋霞免费鲁丝片| 色婷婷久久久亚洲欧美| 99re6热这里在线精品视频| 极品人妻少妇av视频| 亚洲国产色片| 十八禁网站网址无遮挡| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品中文字幕在线视频| 精品国产露脸久久av麻豆| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 精品人妻偷拍中文字幕| 男男h啪啪无遮挡| 久久精品aⅴ一区二区三区四区 | 亚洲精品在线美女| 香蕉丝袜av| 青青草视频在线视频观看| tube8黄色片| 亚洲熟女精品中文字幕| 国产精品熟女久久久久浪| 精品少妇黑人巨大在线播放| 一本—道久久a久久精品蜜桃钙片| 亚洲三区欧美一区| 欧美激情高清一区二区三区 | 国产成人欧美| 寂寞人妻少妇视频99o| 看免费av毛片| 日韩精品免费视频一区二区三区| 久久精品国产自在天天线| 人妻少妇偷人精品九色| 久久久国产精品麻豆| 欧美精品av麻豆av| 久久久a久久爽久久v久久| 如日韩欧美国产精品一区二区三区| 国产精品嫩草影院av在线观看| 欧美激情极品国产一区二区三区| 午夜久久久在线观看| www.自偷自拍.com| 亚洲激情五月婷婷啪啪| 一级毛片我不卡| 深夜精品福利| 熟女少妇亚洲综合色aaa.| 丝瓜视频免费看黄片| 亚洲av中文av极速乱| 午夜老司机福利剧场| 亚洲在久久综合| a级毛片黄视频| 精品酒店卫生间| 久久ye,这里只有精品| 国产野战对白在线观看| 亚洲天堂av无毛| a 毛片基地| 亚洲av电影在线进入| 久久影院123| 久久久久久久亚洲中文字幕| 免费播放大片免费观看视频在线观看| 狂野欧美激情性bbbbbb| 大陆偷拍与自拍| 日本爱情动作片www.在线观看| 亚洲三级黄色毛片| 亚洲国产av新网站| 国产无遮挡羞羞视频在线观看| 人成视频在线观看免费观看| 丝袜美足系列| 99re6热这里在线精品视频| 一区二区av电影网| 最近手机中文字幕大全| 精品人妻熟女毛片av久久网站| 视频区图区小说| 亚洲综合精品二区| 大香蕉久久网| 国产日韩一区二区三区精品不卡| 天天影视国产精品| 国产一区二区在线观看av| 精品国产超薄肉色丝袜足j| 国产福利在线免费观看视频| 好男人视频免费观看在线| 亚洲在久久综合| 久久久精品94久久精品| 天天躁狠狠躁夜夜躁狠狠躁| 99热国产这里只有精品6| 国产日韩欧美视频二区| 久久精品国产亚洲av高清一级| 久久 成人 亚洲| 国产成人一区二区在线| 丰满乱子伦码专区| 久久国产精品男人的天堂亚洲| 视频区图区小说| 亚洲精品在线美女| 男女边摸边吃奶| 免费日韩欧美在线观看| 国产av精品麻豆| 亚洲,一卡二卡三卡| 蜜桃在线观看..| 免费久久久久久久精品成人欧美视频| 精品亚洲成国产av| 日韩伦理黄色片| 亚洲成人一二三区av| 男女边摸边吃奶| 女性被躁到高潮视频| 欧美日韩综合久久久久久| 一二三四中文在线观看免费高清| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 新久久久久国产一级毛片| 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| 久久国产精品男人的天堂亚洲| 欧美日韩精品成人综合77777| 欧美日韩av久久| 欧美激情极品国产一区二区三区| √禁漫天堂资源中文www| 欧美成人午夜免费资源| 亚洲av国产av综合av卡| 一本久久精品| 人人妻人人爽人人添夜夜欢视频| 各种免费的搞黄视频| 日韩中文字幕欧美一区二区 | 国产成人av激情在线播放| 天天影视国产精品| 婷婷色综合大香蕉| av不卡在线播放| 欧美日本中文国产一区发布| 一本色道久久久久久精品综合| 亚洲精品第二区| 精品酒店卫生间| 欧美精品国产亚洲| 国产无遮挡羞羞视频在线观看| 亚洲人成77777在线视频| 中文字幕色久视频| 亚洲三区欧美一区| 人体艺术视频欧美日本| 成人漫画全彩无遮挡| av在线app专区| 超色免费av| 久久午夜福利片| 18+在线观看网站| 一级毛片 在线播放| 国产无遮挡羞羞视频在线观看| 国产片特级美女逼逼视频| 老鸭窝网址在线观看| 亚洲综合色惰| 男人爽女人下面视频在线观看| 亚洲国产精品999| 国产一区二区在线观看av| 一二三四在线观看免费中文在| 国产又色又爽无遮挡免| 国产精品国产三级专区第一集| 高清不卡的av网站| 国产一区有黄有色的免费视频| 国产精品国产三级专区第一集| 久久韩国三级中文字幕| 亚洲国产看品久久| 美女主播在线视频| 亚洲一区中文字幕在线| 看免费成人av毛片| 午夜激情久久久久久久| 晚上一个人看的免费电影| 大话2 男鬼变身卡| 黄色怎么调成土黄色| 一级爰片在线观看| 亚洲天堂av无毛| 亚洲四区av| 大香蕉久久网| 亚洲精品国产色婷婷电影| 老女人水多毛片| 亚洲精品久久久久久婷婷小说| 九色亚洲精品在线播放| 久久免费观看电影| 精品少妇一区二区三区视频日本电影 | av线在线观看网站| 亚洲欧美成人综合另类久久久| 午夜精品国产一区二区电影| 最黄视频免费看| 国产爽快片一区二区三区| 久久精品久久久久久噜噜老黄| 午夜精品国产一区二区电影| 亚洲av成人精品一二三区| 国产深夜福利视频在线观看| 色婷婷av一区二区三区视频| 久久久国产欧美日韩av| 精品午夜福利在线看| 欧美在线黄色| 日本猛色少妇xxxxx猛交久久| 人妻 亚洲 视频| 欧美国产精品一级二级三级| 蜜桃在线观看..| 成人影院久久| 男人添女人高潮全过程视频| 亚洲天堂av无毛| 一级黄片播放器| 亚洲av国产av综合av卡| 在线观看www视频免费| 成人黄色视频免费在线看| 欧美 日韩 精品 国产| 国产欧美日韩综合在线一区二区| 久久精品国产a三级三级三级| 精品一区二区三区四区五区乱码 | 午夜福利影视在线免费观看| 制服丝袜香蕉在线| 水蜜桃什么品种好| 日产精品乱码卡一卡2卡三| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 亚洲av在线观看美女高潮| 七月丁香在线播放| 天天影视国产精品| 国产精品 欧美亚洲| 国产午夜精品一二区理论片| 久久人人爽av亚洲精品天堂| 看免费成人av毛片| 久久精品国产综合久久久| 久久人妻熟女aⅴ| 美女国产高潮福利片在线看| 这个男人来自地球电影免费观看 | 亚洲成av片中文字幕在线观看 | 国产成人午夜福利电影在线观看| 黄片播放在线免费| 成人免费观看视频高清| 一本久久精品| 天天操日日干夜夜撸| 精品少妇一区二区三区视频日本电影 | 在线天堂最新版资源| 亚洲欧美清纯卡通| 久久久久久久久久久久大奶| 欧美成人午夜精品| 五月伊人婷婷丁香| 伊人久久国产一区二区| 啦啦啦在线观看免费高清www| 少妇被粗大的猛进出69影院| 午夜免费鲁丝| 欧美最新免费一区二区三区| 国产精品熟女久久久久浪| 天天躁夜夜躁狠狠久久av| 韩国av在线不卡| 男人操女人黄网站| 亚洲图色成人| 亚洲av福利一区| 18+在线观看网站| 久久青草综合色| 午夜福利在线免费观看网站| 国产精品嫩草影院av在线观看| 日韩av在线免费看完整版不卡| 性色av一级| 精品国产乱码久久久久久小说| 亚洲图色成人| 国产亚洲欧美精品永久| 国产精品不卡视频一区二区| 一二三四中文在线观看免费高清| 日本免费在线观看一区| 久久精品夜色国产| 美女主播在线视频| 一个人免费看片子| 国产片内射在线| 欧美日本中文国产一区发布| 美国免费a级毛片| 男的添女的下面高潮视频| 在线观看www视频免费| 亚洲av成人精品一二三区| 一级爰片在线观看| 久久久久人妻精品一区果冻| 女人久久www免费人成看片| 高清欧美精品videossex| 欧美精品人与动牲交sv欧美| 国产精品久久久av美女十八| 多毛熟女@视频| 亚洲在久久综合| 丰满少妇做爰视频| 永久免费av网站大全| 日韩一区二区三区影片| 国产熟女午夜一区二区三区| tube8黄色片| 国产视频首页在线观看| 女人久久www免费人成看片| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产欧美网| 交换朋友夫妻互换小说| 亚洲国产欧美日韩在线播放| 99久久中文字幕三级久久日本| 成年女人在线观看亚洲视频| 日本爱情动作片www.在线观看| 夫妻性生交免费视频一级片| 亚洲欧美精品综合一区二区三区 | 99久久人妻综合| 在线观看三级黄色| 在线精品无人区一区二区三| 夜夜骑夜夜射夜夜干| 久久久亚洲精品成人影院| 蜜桃国产av成人99| 成年女人毛片免费观看观看9 | 国产日韩欧美视频二区| 久久久精品国产亚洲av高清涩受| 久久精品aⅴ一区二区三区四区 | 亚洲视频免费观看视频| 99久久综合免费| 久久久久久久久久久久大奶| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕免费大全7| 国产野战对白在线观看| 18禁裸乳无遮挡动漫免费视频| 国产亚洲精品第一综合不卡| av片东京热男人的天堂| 国产亚洲av片在线观看秒播厂| 少妇的丰满在线观看| 欧美97在线视频| 国产成人精品婷婷| 国产又爽黄色视频| 国产一区亚洲一区在线观看| tube8黄色片| 啦啦啦在线观看免费高清www| 亚洲第一av免费看| xxx大片免费视频| 美女午夜性视频免费| 亚洲国产欧美在线一区| 在线观看一区二区三区激情| 成人国语在线视频| 91成人精品电影| 一级毛片 在线播放| 在线观看人妻少妇| 亚洲成国产人片在线观看| 午夜福利视频在线观看免费| av又黄又爽大尺度在线免费看| tube8黄色片| 亚洲综合色网址| 久热这里只有精品99| 国产精品免费大片| 午夜福利乱码中文字幕| 男女高潮啪啪啪动态图| 18+在线观看网站| 黑人欧美特级aaaaaa片| 欧美黄色片欧美黄色片| 久久久久久久久免费视频了| 日本黄色日本黄色录像| av国产精品久久久久影院| 大香蕉久久网| 免费日韩欧美在线观看| 婷婷色综合大香蕉| 亚洲精品久久成人aⅴ小说| 日日爽夜夜爽网站| 久久久亚洲精品成人影院| 亚洲美女搞黄在线观看| 国产视频首页在线观看| 丝袜美腿诱惑在线| 妹子高潮喷水视频| 亚洲av福利一区| 黄色配什么色好看| 久久久久网色| 久久精品久久久久久噜噜老黄| 久久精品aⅴ一区二区三区四区 | 亚洲国产精品999| 久久国内精品自在自线图片| 免费观看a级毛片全部| 国产淫语在线视频| 国产精品久久久久久av不卡| 18在线观看网站| 精品亚洲成国产av| 校园人妻丝袜中文字幕| 天天操日日干夜夜撸| 免费黄网站久久成人精品| 91久久精品国产一区二区三区| 大陆偷拍与自拍| 精品国产超薄肉色丝袜足j| 久久久久国产一级毛片高清牌| 国产欧美亚洲国产| 一级a爱视频在线免费观看| 久久这里只有精品19| 欧美激情 高清一区二区三区| 春色校园在线视频观看| 午夜av观看不卡| 欧美97在线视频| 毛片一级片免费看久久久久| 午夜日本视频在线| 涩涩av久久男人的天堂| 日韩成人av中文字幕在线观看| 亚洲精品aⅴ在线观看| 久久精品国产亚洲av天美| 精品第一国产精品| 在线看a的网站| 男人舔女人的私密视频| 制服人妻中文乱码| 国产成人精品久久久久久| 精品99又大又爽又粗少妇毛片| 欧美精品av麻豆av| 久久久国产精品麻豆| 女性被躁到高潮视频| 亚洲精品,欧美精品| 韩国精品一区二区三区| 亚洲国产欧美在线一区| 1024香蕉在线观看| 美女福利国产在线| 中文字幕人妻丝袜制服| 黄网站色视频无遮挡免费观看| 日韩不卡一区二区三区视频在线| 国产色婷婷99| 黄片播放在线免费| 国产男女超爽视频在线观看| 在线 av 中文字幕| 久久99精品国语久久久| 亚洲精品国产av蜜桃| 高清av免费在线| 看免费成人av毛片| 在线 av 中文字幕| 亚洲av中文av极速乱| 在线亚洲精品国产二区图片欧美| 亚洲综合色惰| 国产精品国产av在线观看| av国产精品久久久久影院|