翟小偉，蔣上榮，王 博

(1.西安科技大學(xué) 安全科學(xué)與工程學(xué)院，陜西 西安 710054；2.西部礦井開采及災(zāi)害防治教育部重點實驗室，陜西 西安 710054；3.陜西省煤火災(zāi)害防治重點實驗室，陜西 西安 710054)

煤炭作為我國的主體能源，在我國國民經(jīng)濟(jì)發(fā)展中占有重要的戰(zhàn)略地位。同時，我國有70%的礦井有煤自然發(fā)火隱患，礦井火災(zāi)作為礦井五大災(zāi)害之一，嚴(yán)重影響著礦井的安全生產(chǎn)和工人的生命安全。采空區(qū)煤自燃的“三帶”(即散熱帶、氧化升溫帶和窒息帶)是判斷煤層采空區(qū)自燃危險性的重要指標(biāo)[1]。針對煤層采空區(qū)危險區(qū)域劃分方法及判定依據(jù)，國內(nèi)外學(xué)者做了大量的研究，現(xiàn)階段關(guān)于煤自燃“三帶”的劃分方法主要集中在：采空區(qū)遺煤升溫速率、漏風(fēng)風(fēng)速和氧氣濃度三個方面[2-4]。以采空區(qū)煤體升溫速率為“三帶”劃分指標(biāo)的方法在實際中很少運(yùn)用，一是由于煤體是熱的不良導(dǎo)體，對于采空區(qū)溫度的監(jiān)測難度很大；二是采空區(qū)條件極其復(fù)雜，關(guān)于遺煤升溫規(guī)律以及采空區(qū)溫度場分布特征的研究還不明確。此外，由于條件限制，在采空區(qū)內(nèi)進(jìn)行風(fēng)速測定十分困難，采用漏風(fēng)風(fēng)速劃分自燃“三帶”依然處于理論分析階段，在實際運(yùn)用中存在較大的局限性。

目前最常用的采空區(qū)“三帶”劃分方法是，通過預(yù)埋管路或打鉆等方式在采空區(qū)布置測點，抽取氣體進(jìn)行分析，依據(jù)氧氣體積分?jǐn)?shù)來進(jìn)行劃分。研究表明，當(dāng)采空區(qū)氧濃度高于18%時，采空區(qū)漏風(fēng)嚴(yán)重會帶走大量熱能，熱量無法積聚，煤體很難發(fā)生自燃；當(dāng)采空區(qū)氧濃度低于下限氧濃度Cmin[5，6]時，煤體由于缺氧氧化產(chǎn)生的熱量無法達(dá)到自燃，該值可以通過絕熱氧化試驗測算得到?，F(xiàn)階段，關(guān)于采空區(qū)煤自燃“三帶”劃分的研究主要是基于U型通風(fēng)條件下，采用多種手段和方式對不同開采方法、開采條件進(jìn)行的[7-13]，對于較復(fù)雜通風(fēng)條件下煤層采空區(qū)的危險區(qū)域判定研究較少，對于U+L型通風(fēng)方式，因其多出一條進(jìn)風(fēng)巷或者回風(fēng)巷，造成采空區(qū)的漏風(fēng)比較嚴(yán)重，影響了采空區(qū)煤自燃“三帶”的劃分。本文以陜西杭來灣煤礦U+L型通風(fēng)工作面為例，通過絕熱氧化試驗確定了該煤層氧化自燃的下限氧濃度，進(jìn)而通過現(xiàn)場觀測與數(shù)值模擬相結(jié)合的方法研究采空區(qū)“三帶”區(qū)域分布特征，確定該采空區(qū)煤自燃“三帶”。研究內(nèi)容對類似工作面建立完備的自燃火災(zāi)預(yù)測預(yù)報與防控體系、防止采空區(qū)遺煤自燃具有一定的指導(dǎo)意義。

1 采空區(qū)煤自燃“三帶”劃分基礎(chǔ)

1.1 采空區(qū)煤自燃危險性分析

該礦30105工作面主采3#煤層為易自燃煤層，采用U+L通風(fēng)方式。工作面煤層厚度為8.8～9.6m，平均厚度9.2m，設(shè)計采高4.5m，工作面回采率不低于93%。由于采用大采高分層法開采，且回采過程中頂煤留設(shè)不穩(wěn)定，造成采空區(qū)浮煤較多，這為采空區(qū)煤炭自燃留下隱患。另外，在巷道設(shè)計過程中采用輔運(yùn)巷以減輕主運(yùn)巷的運(yùn)輸壓力，輔運(yùn)巷緊鄰采空區(qū)，采空區(qū)受輔運(yùn)巷的影響漏風(fēng)嚴(yán)重，使得遺煤氧化自燃最終形成火災(zāi)成為可能。故綜上所述，該工作面采空區(qū)具有一定量的遺煤，且因為其通風(fēng)方式的原因，有較大的漏風(fēng)，有一定的遺煤自燃危險性。

1.2 采空區(qū)煤自燃“三帶”劃分依據(jù)

通過絕熱氧化試驗進(jìn)行研究，得到了極限浮煤厚度、下限氧濃度、極限漏風(fēng)強(qiáng)度等煤自燃極限特征參數(shù)，為采空區(qū)煤自燃“三帶”劃分提供了依據(jù)。根據(jù)遺煤自燃特性，本文選取氧濃度作為采空區(qū)煤自燃“三帶”的劃分指標(biāo)，其中散熱帶氧濃度大于18%，窒息帶氧濃度小于絕熱氧化試驗測算出的下限氧濃度Cmin，氧化升溫帶氧濃度介于下限氧濃度Cmin和18%之間。

1.2.1 試驗裝置

運(yùn)用煤自然發(fā)火實驗臺進(jìn)行絕熱氧化試驗測試。該實驗臺由爐體、氣路及控制檢測三部分組成。爐體呈圓柱形，總裝煤量1.5～2t，外部安裝有爐體保溫層，使?fàn)t內(nèi)煤體處于良好的蓄熱環(huán)境，爐體頂、底部均有氣流緩沖層，使氣流由下向上均勻通過實驗煤體，空氣預(yù)熱之后從爐體底部送入，爐內(nèi)布置了測溫探頭和氣體采樣點。實驗裝置的結(jié)構(gòu)如圖1所示。

圖1 絕熱氧化實驗臺結(jié)構(gòu)圖

1.2.2 試驗過程及數(shù)據(jù)

采集工作面新鮮煤樣2t，密封包裝后運(yùn)至自然發(fā)火實驗室，運(yùn)用鄂式破碎機(jī)破碎并裝入實驗爐體，實驗煤樣使用混合粒度，經(jīng)稱重，共裝煤1561.51kg。預(yù)先向試驗爐體內(nèi)通入干空氣，排出爐內(nèi)氣體后開始實驗，實驗過程抽取不同溫度點的氧化氣體，由SP-2120色譜進(jìn)行分析。具體實驗條件見表1，實驗所得的原始數(shù)據(jù)曲線如圖2所示，并用于下一步計算。

表1 實驗條件

圖2 絕熱氧化實驗原始數(shù)據(jù)曲線

由絕熱氧化升溫試驗，得到了煤自然發(fā)火過程中的煤體溫度、氧氣濃度隨著反應(yīng)時間的變化曲線，并通過計算得到了煤體放熱強(qiáng)度隨著反應(yīng)時間的變化曲線，為后續(xù)的下限氧濃度計算提供了基礎(chǔ)參數(shù)。

1.2.3 工作面浮煤厚度測算

由于采空區(qū)浮煤厚度分布的復(fù)雜性和采空區(qū)的封閉性，目前還沒有準(zhǔn)確測量浮煤厚度的方法，因此，浮煤厚度是通過煤層厚度、回采率以及孔隙率等經(jīng)過經(jīng)驗計算的方法得到的。

30105工作面平均采高4.5m，工作面中部回采率在93%左右，采空區(qū)內(nèi)部空隙率根據(jù)經(jīng)驗考慮為30%。則工作面采空區(qū)浮煤厚度可推算如下：

浮煤厚度=采高×(1-回采率)/(1-孔隙率)

(1)

將30105工作面的各參數(shù)帶入上式，求得采空區(qū)浮煤厚度為0.45m。

1.2.4 下限氧濃度計算結(jié)果

下限氧濃度與煤的氧化放熱性、松散煤體堆積厚度、周圍散熱條件以及煤體溫度有關(guān)。近似計算式為[14]：

該工作面推進(jìn)速度以17.7m/d來計算，當(dāng)采空區(qū)長度為400m時，需推進(jìn)約23d，由煤自然發(fā)火實驗數(shù)據(jù)可知，此時煤樣溫度升至60℃，根據(jù)實驗結(jié)果測算出煤體在60℃時，極限浮煤厚度與下限氧濃度的對應(yīng)關(guān)系見表2。

表2 不同浮煤厚度時的下限氧濃度

根據(jù)計算所得的采空區(qū)浮煤厚度為0.45m，選取該浮煤厚度時下限氧濃度為12.90%，則該工作面采空區(qū)煤自燃“三帶”劃分標(biāo)準(zhǔn)見表3。

表3 “三帶”劃分標(biāo)準(zhǔn)

2 采空區(qū)煤自燃“三帶”現(xiàn)場觀測研究

2.1 采空區(qū)氣體觀測方法

采空區(qū)的氣體觀測采用兩工作面巷道預(yù)埋束管的方法，采用負(fù)壓抽氣泵取樣進(jìn)行色譜分析。設(shè)計采空區(qū)測定范圍為325m，每75m布置一個測點，兩工作面巷道各3個測點，編號1#—6#，同步觀測。具體分布如圖3所示。

圖3 工作面采空區(qū)“三帶”觀測測點布置圖

2.2 采空區(qū)O2濃度分布

根據(jù)現(xiàn)場氣體觀測數(shù)據(jù)，以氧氣濃度作為劃分依據(jù)，分別對6個測點監(jiān)測的O2濃度數(shù)據(jù)進(jìn)行曲線擬合，得出進(jìn)、回風(fēng)側(cè)氧氣濃度隨埋深變化趨勢。同時采用超前預(yù)測的方法，按照三個測點中氧化升溫帶范圍最大時來進(jìn)行劃分，從而對煤自燃危險進(jìn)行提前預(yù)防，標(biāo)注如圖4、圖5所示。

圖4 進(jìn)風(fēng)側(cè)O2濃度隨埋深變化規(guī)律

圖5 回風(fēng)側(cè)O2濃度隨埋深變化規(guī)律

由圖4可以看出，三條曲線變化規(guī)律趨于一致，進(jìn)風(fēng)側(cè)氧氣濃度隨埋深增加而降低，且降低速度隨埋深增加而加快。1#、2#、3#測點分別在埋深達(dá)到260m、230m、190m左右時氧氣濃度才降低至18%，進(jìn)入氧化升溫帶；埋深分別達(dá)到330m、335m、310m時，氧濃度降低至下限氧濃度12.9%，進(jìn)入窒息帶。說明該采空區(qū)進(jìn)風(fēng)側(cè)存在較為嚴(yán)重的漏風(fēng)現(xiàn)象，主要是因為受30105工作面緊鄰的輔助運(yùn)輸巷及聯(lián)絡(luò)巷的影響，工作面推進(jìn)后，由于煤柱的存在使得采空區(qū)進(jìn)風(fēng)側(cè)頂板在很長一段距離內(nèi)無法垮落，臨近巷道側(cè)的采空區(qū)頂板垮落不嚴(yán)實，導(dǎo)致采空區(qū)孔隙較大；另外采空區(qū)覆蓋的聯(lián)絡(luò)巷密閉性不良，二者的綜合作用是導(dǎo)致采空區(qū)漏風(fēng)嚴(yán)重的主要因素。

由圖5可以看出，4#、5#、6#三個測點氧濃度降低速度隨著埋深的增加而加快，與采空區(qū)進(jìn)風(fēng)側(cè)具有相同的趨勢。由4#、5#、6#三個測點可以得到，在埋深分別為34m、68m、92m時，氧濃度降低至18%，進(jìn)入氧化升溫帶；埋深分別為120m、135m、141m時，氧濃度降低至12.9%，進(jìn)入窒息帶；由此可知，該工作面回風(fēng)側(cè)氧濃度也較高，但隨著埋深的增加，回風(fēng)側(cè)氧濃度降低速度快于進(jìn)風(fēng)側(cè)，可見回風(fēng)側(cè)漏風(fēng)量要低于進(jìn)風(fēng)側(cè)，氧化升溫帶寬度明顯大于進(jìn)風(fēng)側(cè)采空區(qū)。

3 采空區(qū)煤自燃“三帶”數(shù)值模擬研究

3.1 幾何模型及邊界條件

以工作面實際為原型，利用ANSYS Fluent模擬軟件建立采空區(qū)幾何模型并進(jìn)行網(wǎng)格劃分，如圖6、圖7所示。采空區(qū)長度為400m，寬299.1m，浮煤厚度取0.5m，浮煤上為4.5m厚的巖石，采空區(qū)旁每隔14.3m有長60m，寬19.5m的煤柱，煤柱間留有與輔運(yùn)巷相連的19.5m寬，14.3m長的聯(lián)絡(luò)巷。計算區(qū)域劃分網(wǎng)格為非結(jié)構(gòu)化網(wǎng)格，浮煤中網(wǎng)格在X、Y、Z三個方向上步長為1m，巖石中步長為0.7m，共劃分網(wǎng)格3058605個。

圖6 幾何模型

圖7 網(wǎng)格劃分示意圖

邊界條件設(shè)置應(yīng)根據(jù)工作面實際工況而定，風(fēng)流入口邊界類型均設(shè)置為速度入口(velocity- inlet)，根據(jù)現(xiàn)場測試得到主運(yùn)、輔運(yùn)風(fēng)速分別為1.39m/s、1.72m/s，氣流溫度300K；遺煤孔隙率為0.3；巖體孔隙率設(shè)置為0.4；數(shù)值模擬中，煤巖體密度、比熱等材料參數(shù)可以由《杭來灣煤礦地質(zhì)勘探報告》的力學(xué)性質(zhì)試驗結(jié)果以及相關(guān)公式計算[15]得到，分別為：遺煤密度1600kg/m3，比熱1200J/(kg·K)；巖體密度1800kg/m3，比熱920J/(kg·K)；風(fēng)流出口邊界類型設(shè)置為自由出流(outflow)。

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

通過ANSYS Fluent模擬軟件進(jìn)行數(shù)值求解，得到采空區(qū)氧濃度分布規(guī)律，采空區(qū)氧氣濃度分布等值線與氧化升溫帶分布云圖分別如圖8、圖9所示，進(jìn)、回風(fēng)側(cè)采空區(qū)氧氣濃度隨埋深變化規(guī)律圖分別如圖10、圖11所示。

圖8 采空區(qū)氧氣濃度分布等值線圖

圖9 氧化升溫帶分布圖

圖10 進(jìn)風(fēng)側(cè)O2濃度隨埋深變化規(guī)律

圖11回風(fēng)側(cè)O2濃度隨埋深變化規(guī)律

從圖10、11中可以看出，在采空區(qū)回風(fēng)側(cè)，氧氣濃度最高值約為18.5%，最低值約為2.5%，而在采空區(qū)進(jìn)風(fēng)側(cè)，氧氣濃度從21%逐漸降低至3%，采空區(qū)進(jìn)風(fēng)側(cè)氧氣濃度下降速度要慢于回風(fēng)側(cè)，同一深度氧濃度差別較大。在采空區(qū)進(jìn)風(fēng)側(cè)埋深較淺處，氧濃度存在先減小后增加的波動區(qū)域，這是由于工作面進(jìn)風(fēng)口向采空區(qū)的漏風(fēng)較大，且輔運(yùn)巷通過聯(lián)絡(luò)巷向采空區(qū)漏風(fēng)也較為嚴(yán)重，使得距離工作面較近的區(qū)域及靠近聯(lián)絡(luò)巷的進(jìn)風(fēng)側(cè)氧氣濃度偏高，這為遺煤自燃提供了條件，且難以判定其危險性，對礦井的安全生產(chǎn)會產(chǎn)生較為嚴(yán)重的影響。

3.3 采空區(qū)煤自燃“三帶”劃分結(jié)果

通過絕熱氧化升溫試驗得到采空區(qū)煤自燃下限氧濃度，以18%、12.9%為上、下劃分界限對現(xiàn)場監(jiān)測結(jié)果與數(shù)值模擬結(jié)果進(jìn)行處理，以超前預(yù)測為原則得到該U+L型工作面三帶危險區(qū)域劃分結(jié)果見表4。

表4 “三帶”劃分結(jié)果 m

通過現(xiàn)場觀測與數(shù)值模擬結(jié)合的方法對采空區(qū)煤自燃“三帶”危險區(qū)域進(jìn)行劃分，得到采空區(qū)自燃氧化的危險區(qū)域，結(jié)果表明由于輔運(yùn)巷漏風(fēng)作用的影響，采空區(qū)氧化升溫帶區(qū)域較大，由其對進(jìn)風(fēng)側(cè)采空區(qū)影響嚴(yán)重，出現(xiàn)了煤自燃“三帶”區(qū)域的波動，這部分位置也有發(fā)生煤層自燃的可能，現(xiàn)場觀測過程中會由于束管監(jiān)測密度較小導(dǎo)致出現(xiàn)危險區(qū)域判定的準(zhǔn)確性降低，忽略小范圍氧化升溫帶預(yù)測，自燃危險相對較大。

4 結(jié) 論

通過現(xiàn)場觀測與數(shù)值模擬結(jié)合的方法，對大采高、U+L型工作面采空區(qū)危險區(qū)域進(jìn)行研究，主要得到以下結(jié)論：

1)通過絕熱氧化升溫試驗，確定該工作面采空區(qū)在不同浮煤厚度下，煤體氧化自燃的下限氧濃度，并根據(jù)浮煤厚度測算結(jié)果，得到該采空區(qū)自燃的下限氧濃度值為12.9%。

2)根據(jù)現(xiàn)場觀測與數(shù)值模擬結(jié)果，確定了采空區(qū)O2濃度隨埋深的變化關(guān)系，掌握了采空區(qū)的煤自燃規(guī)律及特點，進(jìn)而對采空區(qū)進(jìn)行了“三帶”危險區(qū)域劃分。

3研究結(jié)果表明U+L型通風(fēng)方式在進(jìn)風(fēng)側(cè)由于輔運(yùn)巷漏風(fēng)作用，會存在O2濃度的波動，煤自燃“三帶”區(qū)域出現(xiàn)非穩(wěn)態(tài)現(xiàn)象，增加了自燃危險性。應(yīng)增加測點密度，加強(qiáng)觀測，及時采取防滅火措施對采空區(qū)進(jìn)行治理。

4)利用現(xiàn)場觀測和數(shù)值模擬相結(jié)合的研究方法進(jìn)行采空區(qū)自燃危險區(qū)域的劃分，能夠更加準(zhǔn)確地得到采空區(qū)煤自燃“三帶”的實際分布規(guī)律，為制定有針對性的綜合防火措施提供科學(xué)依據(jù)和理論基礎(chǔ)。