賈泉升,司榮軍,李潤之,王 磊
(1.瓦斯災(zāi)害監(jiān)控與應(yīng)急技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,重慶 400037;2.中煤科工集團(tuán)重慶研究院有限公司,重慶 400037)
瓦斯爆炸過程中,能量的瞬間大量釋放會產(chǎn)生高溫環(huán)境。煤礦井下瓦斯爆炸位置一般發(fā)生在掘進(jìn)巷道、回采工作面[1],爆炸能量不易消散,在巷道中產(chǎn)生持續(xù)性高溫,給井下人員、設(shè)備造成一定損害。受初始瓦斯?jié)舛扔绊?,爆炸產(chǎn)生的高溫不同。因此,研究不同初始瓦斯?jié)舛葘Ρㄟ^程中溫度特性的影響,可降低煤礦瓦斯爆炸事故給人員、設(shè)備帶來的危害,為煤礦瓦斯爆炸事故溫度研究提供重要依據(jù)。
目前,國內(nèi)外學(xué)者針對瓦斯爆炸過程中溫度變化展開研究:Olsen[2]通過建立數(shù)學(xué)物理方程計算瓦斯爆炸溫度;文獻(xiàn)[3-9]通過利用FLACS、FLUENT、AutoReaGas、CHEMKIN等數(shù)值模擬軟件,研究爆炸在定容或管道傳播條件下的溫度變化規(guī)律,認(rèn)為瓦斯爆炸過程中的溫度最高達(dá)2 400 K以上,爆炸過程中的溫度變化具有振蕩特征,且初始瓦斯?jié)舛仍?.5%時溫度上升速率最高;部分學(xué)者通過模擬瓦斯室內(nèi)爆炸溫度變化規(guī)律,認(rèn)為爆炸過程中同一時刻室內(nèi)超壓分布基本均勻,空間差異較小,最高溫度可達(dá)2 226 K[10]。數(shù)值模擬研究一般在等溫、絕熱條件下進(jìn)行,與實(shí)際實(shí)驗(yàn)數(shù)據(jù)存在一定誤差,不能準(zhǔn)確模擬實(shí)際瓦斯爆炸過程中溫度變化。因此,部分學(xué)者在實(shí)驗(yàn)條件下對爆炸過程中的溫度變化進(jìn)行研究:王從銀等[11]通過電壓信號表示瓦斯爆炸在管道中傳播時的溫度變化發(fā)現(xiàn),瓦斯爆炸過程中管道中上部釋放的熱量大于下部;文獻(xiàn)[12-13]通過研究爆炸在管道中傳播時不同位置處的溫度變化規(guī)律發(fā)現(xiàn),管道上部火焰溫度高于管道下部火焰溫度;崔建霞等[14]利用R型微細(xì)熱電偶測試小型管道瓦斯爆炸溫度最高為1 920 K;文獻(xiàn)[15-16]利用C2-7-K、C2-1-K型熱電偶研究爆炸傳播過程中的溫度變化規(guī)律,測試溫度最高達(dá)1 292.27 K;劉貞堂等[17]通過分析管道爆炸傳播火焰與溫度關(guān)系發(fā)現(xiàn),溫度升高能夠促進(jìn)火焰的傳播;Nie等[18]利用二維溫度場的輻射測溫法研究爆炸火焰周圍溫度變化規(guī)律。
目前,針對瓦斯爆炸溫度的研究方法主要以數(shù)值模擬為主,實(shí)驗(yàn)條件下的研究也主要以不同測點(diǎn)處的溫度變化為主,針對不同初始瓦斯?jié)舛葪l件下爆炸溫度特性及溫度、壓力之間的相互作用關(guān)系的研究相對較少。因此,本文基于爆炸實(shí)驗(yàn),對不同初始瓦斯?jié)舛葪l件下爆炸過程中的溫度變化以及溫度、壓力之間的相互作用關(guān)系進(jìn)行研究,研究結(jié)果可為完善瓦斯爆炸溫度變化機(jī)理及災(zāi)害防控提供依據(jù)。
整個實(shí)驗(yàn)系統(tǒng)示意如圖1所示。實(shí)驗(yàn)配氣系統(tǒng)主要由電磁閥、配氣壓力傳感器及激光CH4濃度傳感器組成。實(shí)驗(yàn)過程中借助激光CH4濃度傳感器對罐體中的CH4濃度進(jìn)行準(zhǔn)確測試,通過微調(diào)使罐體中的瓦斯?jié)舛冗_(dá)到設(shè)定值,精度為0.1%。
圖1 實(shí)驗(yàn)系統(tǒng)示意Fig.1 Schematic diagram of experimental system
數(shù)據(jù)采集系統(tǒng)由測試傳感器和瞬態(tài)信號測試系統(tǒng)組成(采樣頻率為100 kHz),測試傳感器包括溫度、壓力和火焰?zhèn)鞲衅?如圖1中6所示,從上至下依次是壓力、溫度、火焰?zhèn)鞲衅?;溫度傳感器采用C型熱電偶,探測端距離罐體壁103 mm,傳感器精度0.1 K,量程273~2 373 K,響應(yīng)時間為ms級,可準(zhǔn)確測試爆炸過程中的溫度;壓力傳感器采用CYG400系列的高頻動態(tài)壓力傳感器(固有頻率200 kHz),可以對爆炸過程中的壓力信號進(jìn)行準(zhǔn)確測定,能夠測試罐體壁面的壓力;火焰?zhèn)鞲衅鳛镃KG100型,響應(yīng)時間≤10 μs,可對爆炸罐體中的火焰進(jìn)行采集,并以火焰信號采集通道為內(nèi)觸發(fā)通道,傳感器接收到火焰信號時,經(jīng)瞬態(tài)信號測試系統(tǒng)處理后,同時觸發(fā)溫度、壓力傳感器對爆炸罐體中的溫度、壓力數(shù)據(jù)進(jìn)行采集、分析和存儲。
實(shí)驗(yàn)首先抽真空,并借助壓力配比法將爆炸罐體中的瓦斯-空氣混合氣配制到0.15 MPa,靜置15 min后將罐體的排氣孔通過氣體流量計和激光CH4濃度傳感器連接,調(diào)節(jié)氣體的流量使其穩(wěn)定在200~300 mL/min,2 min后讀取CH4濃度。通過向罐體中加入空氣進(jìn)行微調(diào),直至罐體中瓦斯氣體濃度達(dá)到實(shí)驗(yàn)所需值,然后將罐體中多余的混合氣體排出,最后借助20 L球形實(shí)驗(yàn)系統(tǒng)進(jìn)行點(diǎn)火,并采集爆炸罐中的溫度、壓力數(shù)據(jù),實(shí)驗(yàn)結(jié)束后對罐體中的氣體置換2~3次,待罐體中的溫度恢復(fù)常溫后進(jìn)行下一次實(shí)驗(yàn)。實(shí)驗(yàn)步長為0.5%,點(diǎn)火能量10 J,環(huán)境溫度293 K。
瓦斯爆炸過程中產(chǎn)生的高溫傷害性極大,當(dāng)初始瓦斯?jié)舛葹?.5%時,爆炸過程中溫度隨時間的變化規(guī)律如圖2所示。
圖2 初始瓦斯?jié)舛葹?.5%時爆炸過程中溫度變化Fig.2 Curve of temperature change during 9.5% gas explosion
溫度延遲時間指火焰出現(xiàn)至溫度開始上升的時間。由圖2可知,爆炸過程出現(xiàn)二次升溫現(xiàn)象,即火焰出現(xiàn)后經(jīng)歷41 ms溫度延遲后急劇上升,在爆炸時間t1即331 ms時,出現(xiàn)第1個溫度峰值T1為909.22 K,在爆炸時間t2即450 ms時降至894.35 K,然后溫度逐漸升高,在爆炸時間t2為在570 ms時,達(dá)到第2個溫度峰值T2為937.66 K,隨后溫度逐漸降低,主要原因是爆炸溫度測試點(diǎn)位于距離罐體中心103 mm處,受爆炸罐體約束作用及爆炸壓力等因素影響,火焰在20 L罐體中出現(xiàn)一定振蕩、疊加,火焰燃燒波的疊加在一定程度上促進(jìn)溫度的二次升高。
爆炸最高溫指在某初始瓦斯?jié)舛戎禃r,所測爆炸過程中溫度的最大值。爆炸極限范圍內(nèi)爆炸最高溫隨初始瓦斯?jié)舛鹊淖兓?guī)律如圖3所示。
圖3 爆炸最高溫隨初始瓦斯?jié)舛茸兓疐ig.3 Change of maximum explosion temperature with gas concentration
由圖3可知,爆炸最高溫與初始瓦斯?jié)舛瘸收袷幾兓?guī)律。隨初始瓦斯?jié)舛壬?,爆炸最高溫呈升?降低-升高-降低-升高-降低的現(xiàn)象。初始瓦斯?jié)舛仍?.5%、9.5%、12%時,爆炸最高溫出現(xiàn)極大值EV1、EV2、EV3分別為995,932,1 153 K,其中當(dāng)初始濃度為12%時,爆炸最高溫相對最大。原因是瓦斯爆炸過程中的能量主要以氣體膨脹和熱能等形式釋放[19],實(shí)驗(yàn)中能量主要存在于高壓膨脹氣體、高溫環(huán)境、反應(yīng)產(chǎn)物以及熱傳導(dǎo)中,初始瓦斯?jié)舛炔煌?,對?yīng)爆炸前反應(yīng)能量不同。爆炸過程中產(chǎn)生的高溫環(huán)境是能量釋放形式之一,爆炸最大壓力隨瓦斯?jié)舛鹊淖兓?guī)律如圖4所示。由圖4可知,當(dāng)瓦斯?jié)舛葹?.5%時,產(chǎn)生的爆炸壓力相對最高,爆炸以高壓膨脹氣體形成釋放的能量大于瓦斯?jié)舛葹?.5%、12%時釋放的能量,導(dǎo)致瓦斯?jié)舛葹?.5%時溫度相對最低;初始瓦斯?jié)舛葹?2%時,由于參與反應(yīng)的CH4較多,爆炸釋放的總能量高于初始瓦斯?jié)鉃?.5%釋放的總能量,導(dǎo)致初始瓦斯?jié)舛葹?2%時的爆炸最高溫高于初始瓦斯?jié)舛葹?.5%時的爆炸最高溫。不同初始瓦斯?jié)舛缺ㄟ^程中的火焰持續(xù)時間不同,使爆炸最高溫隨初始瓦斯?jié)舛瘸收袷幾兓?/p>
圖4 爆炸最大壓力隨瓦斯?jié)舛茸兓€Fig.4 Change of maximum explosion pressure with gas concentration
點(diǎn)火源位于球形爆炸罐中心,受溫度傳感器安裝位置影響,當(dāng)點(diǎn)火源引爆周圍瓦斯氣體后,溫度傳感器沒有即刻探測到溫度的變化,出現(xiàn)時間延遲。初始瓦斯?jié)舛炔煌?,從火焰出現(xiàn)至溫度上升時間長短不同。溫度延遲時間隨初始瓦斯?jié)舛鹊淖兓?guī)律如圖5所示。升溫時間即升溫開始至達(dá)到最高溫所用時間,一定程度可反映爆炸溫度變化快慢。
圖5 溫度延遲時間隨初始瓦斯?jié)舛茸兓€Fig.5 Change of temperature delay time with initial gas concentration
由圖5可知,溫度延遲時間與初始瓦斯?jié)舛茸兓€近似呈U型。初始瓦斯?jié)舛仍诒ㄉ?、下限附近范圍?nèi)溫度延遲時間較長;當(dāng)初始瓦斯?jié)舛葹?.5%時,溫度延遲時間相對較短,為46 ms;當(dāng)初始瓦斯?jié)舛葹?.5%~6.0%時,溫度延遲時間出現(xiàn)斷崖式縮短;當(dāng)初始瓦斯?jié)舛葹?1.5%~14%時,溫度延遲時間快速增長,最大值達(dá)到550 ms。從鏈?zhǔn)椒磻?yīng)(著火的電理論)角度分析,瓦斯?jié)舛仍诒ㄉ?、下限附近變化時,能夠較大程度改變爆炸反應(yīng)的化學(xué)活性,從而進(jìn)一步影響爆炸強(qiáng)度和溫度延遲時間;初始瓦斯?jié)舛仍?.5%附近時,爆炸強(qiáng)度較大,爆炸熱傳遞速度相對最快,爆炸過程中的溫度延遲時間受瓦斯?jié)舛扔绊懴鄬^小。
升溫時間隨初始瓦斯?jié)舛茸兓鐖D6所示。由圖6可知,初始瓦斯?jié)舛仍诒ㄉ?、下限附近時,升溫時間較長,為1 200 ms;當(dāng)初始瓦斯?jié)舛仍?%~12%時,升溫時間介于475~625 ms,主要原因是初始瓦斯?jié)舛容^低,CH4分子間的距離相對較大,當(dāng)點(diǎn)火發(fā)生時,參與反應(yīng)的CH4分子較少,致使爆炸反應(yīng)的鏈引發(fā)、傳遞較為緩慢,溫度開始上升時間較晚;當(dāng)初始瓦斯?jié)舛戎饾u升高時,點(diǎn)火源周圍CH4分子量增加,爆炸反應(yīng)敏感性迅速增強(qiáng),極易形成爆炸反應(yīng)鏈,溫度快速升高,導(dǎo)致升溫時間在爆炸下限出現(xiàn)斷崖式縮短。當(dāng)瓦斯?jié)舛仍?%~12%時,點(diǎn)火源周圍的CH4分子運(yùn)動相對活躍,受瓦斯?jié)舛扔绊?,升溫時間變化不明顯;當(dāng)初始瓦斯?jié)舛葹楸ㄉ舷迺r,隨瓦斯?jié)舛壬撸鯘舛冉档?,不足以提供CH4反應(yīng)需要的氧化物,成為限制爆炸發(fā)生的主要因素之一,且對爆炸影響較大,因此升溫時間在爆炸上限也存在斷崖式變化。
圖6 升溫時間隨初始瓦斯?jié)舛茸兓€Fig.6 Change of heating time with initial gas concentration
溫升速率是描述瓦斯爆炸過程中溫度變化快慢的物理量。對爆炸過程中溫度數(shù)據(jù)進(jìn)行分析得出,低濃度(6.5%)、當(dāng)量濃度(9.5%)和高濃度(12%)瓦斯爆炸時溫升速率變化規(guī)律,如圖7所示。由圖7可知,初始瓦斯?jié)舛葹?.5%時,瓦斯爆炸過程中溫升速率最高為2.7 K/ms;初始瓦斯?jié)舛葹?.5%時,在97 ms時溫升速率達(dá)到最大值5.41 K/ms,后逐漸降為0 K/ms,然后出現(xiàn)小幅度升高,再緩慢降低至-0.5 K/ms,并最終穩(wěn)定在0 K/ms附近。爆炸溫升速率初始階段急劇上升,說明初始瓦斯?jié)舛葹?.5%附近爆炸劇烈,溫度上升較快;當(dāng)初始瓦斯?jié)舛葹?2%時,溫升速率最高達(dá)2.8 K/ms,小于瓦斯?jié)舛葹?.5%時的溫升速率。
圖7 不同初始瓦斯?jié)舛葪l件下溫升速率與溫度的關(guān)系Fig.7 Relationship between temperature rise rate and temperature under different gas concentration conditions
不同瓦斯?jié)舛葪l件下溫升速率對比如圖7(d)所示。由圖7(d)可知,初始瓦斯?jié)舛葹?.5%時爆炸最先開始升溫,說明爆炸過程中溫度上升延遲時間較短,且爆炸過程中溫升速率相對最高;低濃度(6.5%)和高濃度(12%)瓦斯爆炸過程中溫升速率和溫度延遲時間小于當(dāng)量瓦斯?jié)舛?,說明初始瓦斯?jié)舛仍诋?dāng)量濃度時,分子運(yùn)動最為活躍,且各反應(yīng)物相對充足,爆炸極易發(fā)生且爆炸強(qiáng)度相對最大。
以初始瓦斯?jié)舛葹?.5%時的溫度和壓力數(shù)據(jù)為例,研究溫度與壓力間的關(guān)系,如圖8所示。由圖8可知,壓力在115 ms達(dá)到峰值0.899 MPa;溫度在331 ms時第1次達(dá)到峰值為909.22 K,與壓力峰值時間間隔t1為216 ms;溫度在570 ms第2次達(dá)到峰值936.34 K,與壓力峰值時間間隔t2為455 ms。說明瓦斯在爆炸過程中,壓力波的傳播比火焰?zhèn)鞑サ目?,因此溫度達(dá)到第1個峰值的時間較壓力波達(dá)到峰值的時間晚,爆炸壓力波傳至爆炸罐體壁后,受罐體壁反向力作用,使爆炸壓力波反向傳播,同時與相對的火焰波匯合,促使火焰波反向傳播,正向火焰波與反向火焰波疊加,使火焰溫度進(jìn)一步升高。因此,在爆炸強(qiáng)度較大時,瓦斯爆炸溫度會出現(xiàn)二次升溫現(xiàn)象。
圖8 初始瓦斯?jié)舛葹?.5%時爆炸壓力與溫度關(guān)系Fig.8 Relationship between gas explosion pressure and temperature under 9.5% concentration
以初始瓦斯?jié)舛葹?.5%,12%為例,爆炸過程中溫度隨時間變化規(guī)律如圖9所示。由圖9可知,在爆炸過程中沒有出現(xiàn)二次升溫現(xiàn)象,主要原因是受初始瓦斯?jié)舛扔绊?,在爆炸上、下限附近時爆炸強(qiáng)度降低,爆炸過程中壓力上升速率減慢,壓力達(dá)到峰值的時間較長,壓力波到達(dá)爆炸罐體后反向傳播能力減弱。壓力、溫度峰值的時間間隔t3、t4較t2明顯縮短,爆炸過程中溫度峰值和壓力峰值的時間間隔較小,不能促進(jìn)火焰波的反向傳播,導(dǎo)致溫度上升過程中的二次升溫現(xiàn)象減弱,甚至無二次升溫現(xiàn)象。
圖9 初始瓦斯?jié)舛仍?.5%和12%條件下爆炸壓力與溫度關(guān)系Fig.9 Relationship between explosion pressure and temperature under different gas concentration conditions
1)爆炸溫度隨初始瓦斯?jié)舛壬叱收袷幮妥兓?,?2%、6.5%時出現(xiàn)極大值,分別為1 153,995 K,初始濃度為9.5%時,爆炸最高溫為932 K。
2)不同初始瓦斯?jié)舛缺ㄟ^程中,溫度延遲時間、升溫時間隨初始瓦斯?jié)舛瘸蔝型變化,在濃度為6.5%時最高分別達(dá)到550,900 ms,然后呈斷崖式降低;在當(dāng)量濃度附近分別縮短至46,475 ms,在濃度為12.5%時又快速升高。
3)當(dāng)量(9.5%)濃度瓦斯爆炸過程中,壓力波傳播速度大于火焰?zhèn)鞑ニ俣龋瑝毫Σ〞龠M(jìn)火焰燃燒波的反向傳播,出現(xiàn)2個溫度峰值;在爆炸上、下限附近時,爆炸強(qiáng)度降低,壓力波促進(jìn)火焰燃燒波的反向傳播能力減弱,爆炸過程中僅出現(xiàn)1個峰值。