晉牛礦1303綜放面采空區(qū)注氮方案研究及數(shù)值模擬*

賈寶山，汪　偉，祁　云，孫　勇，李守國1,

(1.遼寧工程技術(shù)大學(xué) 安全科學(xué)與工程學(xué)院，遼寧 阜新 123000；2.礦山熱動力災(zāi)害與防治教育部重點(diǎn)實(shí)驗室 遼寧工程技術(shù)大學(xué)，遼寧 葫蘆島 125105；3.煤科集團(tuán)沈陽研究院有限公司 煤礦安全技術(shù)國家重點(diǎn)實(shí)驗室，遼寧 沈陽 113112)

0　引言

煤礦火災(zāi)是影響礦井安全的一大災(zāi)害，而采空區(qū)遺煤自燃火災(zāi)在井下火災(zāi)總數(shù)的占比達(dá)到85%左右[1-2]。有資料表明，全國有約56%的礦井受到自然發(fā)火的威脅[3]。礦井火災(zāi)不僅使煤炭資源遭受破壞，損毀大量的井下設(shè)備，還會產(chǎn)生瓦斯、煤塵爆炸等繼發(fā)性災(zāi)害[4]。隨著綜采放頂煤技術(shù)在我國的廣泛應(yīng)用，礦井生產(chǎn)效率得到明顯提高，與此同時，也出現(xiàn)了工作面推進(jìn)速度降低，漏風(fēng)量增加，采空區(qū)遺煤較多等問題，使采空區(qū)自燃火災(zāi)發(fā)生幾率大大增加[5-7]。

目前，采空區(qū)防滅火技術(shù)主要有均壓通風(fēng)、灌漿、噴灑阻化劑、注氮等，亦或是這些技術(shù)的綜合使用。研究表明，如果封閉火區(qū)內(nèi)的氧氣濃度降到5%以下，火勢就會慢慢減弱最后熄滅[8]。注氮防滅火技術(shù)就是根據(jù)這一原理，將惰性氣體N2注入采空區(qū)，稀釋火區(qū)內(nèi)的氧氣濃度，滿足防火惰化要求。由于N2的注入可以在采空區(qū)產(chǎn)生正壓，使漏風(fēng)量減小，而且N2的密度要小于空氣，可以在采空區(qū)上浮、擴(kuò)散，能夠?qū)Σ煽諈^(qū)深部和高位火災(zāi)起到較好的治理效果，再加上注氮的經(jīng)濟(jì)成本比較低，因此注氮防滅火技術(shù)作為經(jīng)濟(jì)可靠的治理手段己經(jīng)被普遍采用[9-11]。而注氮流量、注氮口位置、注氮時間等參數(shù)制約著注氮效果的好壞，對此許多學(xué)者都進(jìn)行了相關(guān)研究[12-15]，但是尚未明確給出注氮量的模擬取值范圍的標(biāo)準(zhǔn)，沒有系統(tǒng)對氧化自燃帶寬度隨注氮量變化的關(guān)系進(jìn)行定量分析。為此本文結(jié)合前人的研究成果，針對晉牛煤礦1303綜放工作面開展試驗研究，通過在采空區(qū)進(jìn)、回風(fēng)側(cè)布置束管和溫度采樣測點(diǎn)，連續(xù)監(jiān)測采空區(qū)氣體濃度、溫度數(shù)據(jù)，劃分采空區(qū)自燃“三帶”分布區(qū)域，并基于CFD理論，利用流體力學(xué)COMSOL計算軟件，研究了不同注氮量、注氮位置參數(shù)下采空區(qū)自燃“三帶”分布規(guī)律，綜合現(xiàn)場實(shí)測、定量分析、數(shù)值模擬的方法系統(tǒng)地分析注氮參數(shù)對氧化自燃帶寬度變化的影響規(guī)律，結(jié)合工作面實(shí)際生產(chǎn)情況及模擬數(shù)據(jù)擬合式對注氮參數(shù)進(jìn)行優(yōu)化，從而為采空區(qū)遺煤自燃防治技術(shù)的進(jìn)一步提高提供理論指導(dǎo)。

1　采空區(qū)自燃“三帶”實(shí)測分析

1.1　 1303綜放工作面概況

晉牛煤礦1303綜放工作面位于1030水平11#煤層一采區(qū)，該工作面東起集中回風(fēng)大巷，西至井田邊界，南為1305工作面，北為1301工作面。工作面設(shè)計走向長度906.6 m，傾斜長度90 m，煤層厚度為5.24～7.30 m，平均厚度為6.17 m。采用綜采放頂煤開采技術(shù)，全部垮落法控制采空區(qū)頂板，工作面機(jī)采采高2.9 m，放頂煤高度3.3 m，采放比為1∶1.13。煤層節(jié)理發(fā)育，含有1～4層夾矸，夾矸以薄層泥巖為主，厚度變化不大，煤層傾角8～14°，平均為10°。開采煤層屬于Ⅰ類易自燃煤層，最短自然發(fā)火期為21 d，煤塵爆炸指數(shù)為45.79%，具有爆炸性。

1.2　采空區(qū)束管布置方案

針對1303綜放面的基本特征，在進(jìn)風(fēng)巷和回風(fēng)巷鋪設(shè)束管監(jiān)測系統(tǒng)，并利用晉牛煤礦提供的GC-2010型氣相色譜儀，對采空區(qū)內(nèi)各種氣體的含量隨著綜放面回采的變化規(guī)律進(jìn)行測定。采用進(jìn)、回風(fēng)巷同時布置測點(diǎn)的方式，測點(diǎn)相互間距為20 m，每側(cè)3個，具體布置方案如圖1所示，其中1～6#表示束管及溫度傳感器安設(shè)位置。

圖1　束管及溫度傳感器安設(shè)布置Fig.1　Beam tube and temperature sensor layout

測溫采用AD590型集成溫度傳感器，并使用自主研發(fā)的本安型數(shù)字測溫儀表。由于采空區(qū)頂板垮落，容易損壞測點(diǎn)的探頭，為防止頂板冒落及放頂煤時對束管產(chǎn)生巨大的沖擊，在現(xiàn)場測點(diǎn)布置時采用無縫鋼管對束管進(jìn)行保護(hù)。將束管管纜和測溫導(dǎo)線進(jìn)行固定后一起穿入套管內(nèi)，測溫導(dǎo)線預(yù)留一定的長度，防止測溫導(dǎo)線由于拉伸所導(dǎo)致斷裂，套管之間用快速接頭連接牢固。在各個測點(diǎn)處分別安裝上溫度傳感器及采樣頭，采樣頭和束管相連，溫度傳感器與測溫導(dǎo)線相連。束管鋪設(shè)前，將每個測點(diǎn)預(yù)鋪設(shè)的單芯束管與真空泵采氣處的三芯束管進(jìn)行對應(yīng)編號，以便識別監(jiān)測點(diǎn)。

1.3　采空區(qū)實(shí)測數(shù)據(jù)分析

利用兩個多月的時間對1303綜放面采空區(qū)的溫度及氣體進(jìn)行了實(shí)測，從監(jiān)測結(jié)果來看，氧氣濃度在上隅角監(jiān)測位置相對于其他監(jiān)測位置來說下降的更快些，這也正與氣體在采空區(qū)的實(shí)際運(yùn)移規(guī)律相符合。其中1#測點(diǎn)在被埋進(jìn)采空區(qū)48 m時，氧氣濃度開始降至18%，表明已進(jìn)入到氧化自燃帶；隨著推進(jìn)距離繼續(xù)增加，氧氣濃度也進(jìn)一步降低，到90 m的時候氧氣含量的監(jiān)測值已經(jīng)處于10%以下；煤體壓實(shí)程度加深，這一區(qū)域的漏風(fēng)量極小，推進(jìn)到106 m時，氧氣濃度低于8%進(jìn)入窒熄帶。2#，3#測點(diǎn)的記錄結(jié)果大體與1#測點(diǎn)的結(jié)果一致。5#測點(diǎn)的氧氣濃度在37 m時降至18%，埋進(jìn)95 m時降到8%，進(jìn)入窒息帶。4#，6#測點(diǎn)在推進(jìn)到39 m時，氧濃度開始低于18%，后期由于束管被堵不再監(jiān)測。整體上的監(jiān)測結(jié)果較為理想，于是依氧氣濃度劃分的原則，對1303綜放面采空區(qū)遺煤自燃“三帶”的分布區(qū)域進(jìn)行劃分，結(jié)果如表1所示。

表1　根據(jù)氧氣濃度劃分采空區(qū)自燃“三帶”Tab.1　“Three zones” of spontaneous combustion ingoaf according to oxygen concentration

2　1303工作面模型的建立

2.1　幾何模型的建立

根據(jù)晉牛煤礦1303綜放工作面資料及現(xiàn)場實(shí)測，利用COMSOL模擬軟件建立采空區(qū)幾何模型如圖2所示。其中工作面長度WG為90 m，寬度為7 m；采空區(qū)長度GL為260 m，寬度為90 m。邊界W1為進(jìn)風(fēng)巷的進(jìn)風(fēng)口，寬度為3 m；W2為回風(fēng)巷的回風(fēng)口，寬度為3 m；W3，W4，W7為工作面未采實(shí)體煤壁，W5為進(jìn)風(fēng)巷外側(cè)保護(hù)煤壁，W6回風(fēng)巷外側(cè)保護(hù)煤壁，G1為回風(fēng)巷采空區(qū)邊界，G2為進(jìn)風(fēng)巷采空區(qū)邊界，以頂板壓力穩(wěn)定區(qū)域邊界G3為計算區(qū)域的邊緣。

圖2　采空區(qū)幾何模型Fig.2　Geometric model of goaf

2.2　采空區(qū)氣體流動控制方程

由于流體運(yùn)動方程組的建立要滿足基本流體流動定律，當(dāng)氣體處在采空區(qū)深部低速流動，如果設(shè)定采空區(qū)氣體是不可壓縮的，可以將其看成二維層流流動的穩(wěn)態(tài)問題，此時不必考慮能量交換。所以采空區(qū)的風(fēng)流流動就只需要滿足動量守恒方程和連續(xù)性方程[16]：

(1)

(2)

(3)

式中：u,v分別為x,y方向上的風(fēng)流速度，m/s；ρ為礦井中空氣的密度，kg/m3；t為氣體流動時間，s；P為流體微元上的壓力，Pa；μ為采空區(qū)空氣的動力粘度，Pas；Su,Sv為自定義的源項。

2.3　邊界條件及計算參數(shù)設(shè)定

在選取主要計算條件和參數(shù)時，要遵循現(xiàn)場實(shí)際情況而定，將工作面進(jìn)風(fēng)巷設(shè)定為入口邊界；采空區(qū)出口邊界為自由邊界。將注氮口設(shè)定為速度入口，并取氮?dú)獾臐舛葹?7%。實(shí)測得到進(jìn)風(fēng)巷的風(fēng)流溫度是18.6 ℃，氧氣濃度是20.9%，工作面實(shí)際風(fēng)速是1.62 m/s；礦井平均空氣密度ρ=1.225 kg/m3，常溫下取空氣粘性系數(shù)η=1.789 4×10-5kg/ms，氣體的擴(kuò)散系數(shù)D=2.88×10-5m3/s，松散系數(shù)設(shè)置為1.5。采空區(qū)的孔隙率可以依據(jù)頂板冒落碎脹的實(shí)際情況由經(jīng)驗公式來得到：

(4)

采空區(qū)的滲透率由多孔介質(zhì)Carman公式進(jìn)行計算[17-18]：

(5)

式中：x為采空區(qū)距工作面距離，m；Dm為平均粒徑，m；n為孔隙率。

3　注氮量取值范圍的確定

首先要確定出一個合理的模擬方案，通過常用的采空區(qū)注氮設(shè)計方法，得到數(shù)值模擬所需注氮量的取值范圍。目前主要是按照采空區(qū)氧化升溫帶的氧氣濃度和工作面的煤炭日產(chǎn)量這2種參數(shù)計算方法來設(shè)計采空區(qū)注氮量。

3.1　按氧化升溫帶的氧氣濃度設(shè)計注氮量

按照氧化升溫帶的氧氣濃度來設(shè)計采空區(qū)所需注氮量的原理是，使氧化升溫帶內(nèi)原本的氧濃度值稀釋到滿足惰化要求時的氧濃度值，此時注氮量可用下式計算:

(6)

式中：QN為注氮量，m3/h；k為附加系數(shù)，取1.3；Q0為氧化升溫帶內(nèi)的漏風(fēng)量，取2.75 m3/h；C1為氧化升溫帶內(nèi)的平均氧氣濃度，按經(jīng)驗取15%；CN為注入氮?dú)獾臐舛?，?7%，C2為采空區(qū)達(dá)到防火要求時的氧氣濃度，這里為8%。代入到式(6)得出注氮量QN=300 m3/h。

3.2　按工作面煤炭日產(chǎn)量設(shè)計注氮量

按照工作面的煤炭日產(chǎn)量來設(shè)計采空區(qū)所需注氮量的原理是，通過注氮口注入氮?dú)猓瑢⒚刻煲驗榛夭啥斐傻牟煽諈^(qū)空間體積進(jìn)行惰化，稀釋氧化升溫帶內(nèi)氧氣使其濃度降到防火要求指標(biāo)以下，此時注氮量可用下式計算：

(7)

式中：A為煤炭的日產(chǎn)量，t；ρ為煤的密度，取1.28 t/m3；N1為管道輸?shù)?，?.9；N2為采空區(qū)注氮效率，取0.6；C3為空氣中的氧氣含量，取20.9%。代入式(7)得出注氮量QN=500 m3/h。通過以上計算，可以知道注氮量的取值范圍是300～500 m3/h。

4　注氮參數(shù)的數(shù)值模擬研究

4.1　注氮位置影響采空區(qū)流場仿真

因為不同注氮位置條件下，氮?dú)庠诓煽諈^(qū)中擴(kuò)散和運(yùn)移的范圍會產(chǎn)生變化，從而導(dǎo)致采空區(qū)“三帶”的分布區(qū)域也跟著發(fā)生改變，為了確定最佳注氮位置，當(dāng)注氮量為400 m3/h時，在注氮的位置與切頂線間距為10,20,30,50,70 m情況下，分別對采空區(qū)自燃“三帶”的分布情況進(jìn)行數(shù)值模擬。為了使模擬所得的“三帶”分布結(jié)果更貼合實(shí)際，考慮到現(xiàn)場的具體情況，這里采用的是綜合采空區(qū)流場風(fēng)速和氧氣濃度2種劃分指標(biāo)的方法，也就是分別把風(fēng)速0.004 m/s等值線和氧氣濃度8%等值線定為氧化帶的上界限和下界限。

由數(shù)值模擬得到不同注氮位置條件下采空區(qū)自燃“三帶”分布范圍如圖3所示。

(a)注氮口距離切頂線10 m

(b)注氮口距離切頂線20 m

(c)注氮口距離切頂線30 m

(d)注氮口距離切頂線50 m

(e)注氮口距離切頂線70 m圖3　不同注氮位置下采空區(qū)流場和自燃“三帶”分布Fig.3　Goaf flow field and spontaneous combustion three-zone distribution under different injection nitrogen position

可以看出，如果注氮位置發(fā)生改變，則采空區(qū)氧化自燃帶也會呈現(xiàn)不一樣的分布趨勢。隨著注氮口向采空區(qū)深部移動，影響的主要是氧氣濃度8%等值線，而漏風(fēng)速度0.004 m/s等值線的變化趨勢并不明顯。隨著注氮口位置不斷遠(yuǎn)離切頂線，采空區(qū)氧化自燃帶的寬度雖然開始縮小，但之后又逐漸增大。注氮口位置與切頂線間距為10，20，30，50，70 m的情況下，對應(yīng)氧化自燃帶寬度分別為44，37，28.5，26，31.6 m，氧化自燃帶寬度降低為未注氮條件下的75%，63%，49%，45%，54%，表明對采空區(qū)進(jìn)行注氮能夠使采空區(qū)氧化自燃帶寬度明顯縮小。

由于注氮位置參數(shù)會影響注氮效果及經(jīng)濟(jì)成本，因此要結(jié)合實(shí)際生產(chǎn)情況及數(shù)值模擬結(jié)果優(yōu)化出最佳的注氮位置。根據(jù)圖3結(jié)果運(yùn)用Origin軟件繪制了氧化自燃帶寬度與注氮位置變化關(guān)系曲線，如圖4所示。由圖4可知，隨著注氮位置越來越深入采空區(qū)，氧化自燃帶的寬度開始逐漸縮小，但當(dāng)注氮位置與切頂線的間距超過50 m以后氧化自燃帶寬度反而開始增長，這是因為注氮口埋入采空區(qū)過深時，與采空區(qū)窒息帶相距的比較近，對于采空區(qū)較淺位置的氧氣濃度不能起到很好的稀釋惰化作用，再加上注氮壓力的影響，注氮效果受到極大限制，所以結(jié)合現(xiàn)場實(shí)際與理論研究結(jié)果，采空區(qū)最佳注氮位置范圍應(yīng)該在距離工作面30～50 m之間，取30 m為最合適的注氮位置。

圖4　采空區(qū)氧化帶寬度隨注氮口位置變化Fig.4　Goaf oxidized zone width with injection nitrogen position change

4.2　注氮量影響采空區(qū)流場仿真

本次設(shè)計的注氮量模擬方案一共被分成5組，根據(jù)之前所得到的模擬值的取值范圍，當(dāng)注氮位置距離切頂線30 m時，分別在注氮量為300,350,400,450,500 m3/h條件下，對采空區(qū)內(nèi)的遺煤自燃“三帶”的分布情況進(jìn)行數(shù)值模擬。由數(shù)值模擬得到不同注氮量下采空區(qū)流場和自燃“三帶”分布范圍如圖5所示，可以看出變化較大的是氧氣濃度等值線，隨著注氮量加大，氧濃度8%等值線明顯向工作面移動，氧化自燃帶的范圍逐漸縮小，而漏風(fēng)速度等值線的變化趨勢不大。

通過圖5得到注氮量為300，350，400，450，500 m3/h的情況下采空區(qū)內(nèi)的氧化自燃帶的寬度，如表2所示。整體上氧化帶的上界限受注氮量的影響并不是十分顯著，只沿著回采方向移動了不到4 m，而氧化帶的下界限受注氮量的影響相對來說更為顯著，沿著回采方向移動了22 m左右。隨著注氮量由300 m3/h增大到500 m3/h，采空區(qū)內(nèi)的氧化自燃帶的范圍也跟著明顯的縮小，寬度也從41.5 m減小到了23.5 m；當(dāng)注氮量為500 m3/h時，氧化自燃帶的寬度較未對采空區(qū)進(jìn)行注氮時減小了約35 m。

(a)注氮量為300 m3/h

(b)注氮量為350 m3/h

(c)注氮量為400 m3/h

(d)注氮量為450 m3/h

(e)注氮量為500 m3/h圖5　不同注氮量下采空區(qū)流場和自燃“三帶”分布Fig.5　Goaf flow field and spontaneous combustion three-zone distribution under different injection nitrogen quantity

注氮量/(m3·h-1)氧化帶上限位置/m氧化帶下限位置/m寬度/m30046.58841.535045.880.334.54004573.528.54504471275004366.523.5

由于注氮量過大會使注氮成本增加，并且向工作面泄漏氮?dú)獾奈ｋU也逐漸增加，容易使工作面的氧氣濃度降到18.5%以下，從而對工人的身體健康造成危害，因此要結(jié)合實(shí)際生產(chǎn)情況及數(shù)值模擬結(jié)果優(yōu)化出最佳的注氮量。根據(jù)表2運(yùn)用Origin軟件對數(shù)據(jù)進(jìn)行擬合處理，得到采空區(qū)氧化自燃帶的寬度隨注氮量變化曲線如圖6所示。由圖可知，注氮量與氧化自燃帶的寬度基本符合指數(shù)關(guān)系，根據(jù)所得擬合式有：

vτ≥L=120.4exp(-0.004 7QN)+11.85

(8)

式中:v為工作面的回采速度；τ為自然發(fā)火期；晉牛煤礦1303綜放工作面的實(shí)際回采速度為1.6 m/d，為了考慮安全，取回采速度為v=1.5 m/d，自然發(fā)火期τ=21 d。通過計算可知最優(yōu)的氮?dú)庾⑷肓渴?86 m3/h，此時氧化自燃帶寬度為31.5 m。

5　結(jié)論

1)通過對各測點(diǎn)氧氣濃度隨工作面回采的變化趨勢的測定，劃分了晉牛煤礦1303綜放工作面采空區(qū)自燃“三帶”分布區(qū)域，得出采空區(qū)進(jìn)風(fēng)側(cè)距工作面48～106 m，回風(fēng)側(cè)距工作面37～95 m的范圍為氧化自燃帶。

2)利用COMSOL軟件模擬了注氮情況下采空區(qū)自燃“三帶”分布的變化規(guī)律，隨著注氮口位置向采空區(qū)深部移動，氧化自燃帶寬度開始縮小而后增大，最合適的注氮位置應(yīng)該在距離切頂線30 m左右。

3)通過數(shù)值模擬發(fā)現(xiàn)隨著注氮量的加大，氧化自燃帶的范圍逐漸縮短，結(jié)合工作面實(shí)際生產(chǎn)情況及模擬結(jié)果對注氮量參數(shù)進(jìn)行優(yōu)化，發(fā)現(xiàn)注氮量與氧化升溫帶寬度基本符合指數(shù)關(guān)系，由擬合式算出采空區(qū)最優(yōu)注氮量為386 m3/h，此時氧化帶的寬度相比未采取注氮措施時減小了27 m。

[1]王德明.礦井火災(zāi)學(xué)[M].徐州:中國礦業(yè)大學(xué)出版社,2011.

[2]段春生,雷炎云,李世雄,等.不同風(fēng)量下采空區(qū)自燃三帶判定及防滅火技術(shù)[J].煤炭科學(xué)技術(shù),2015,43(1):62-65.

DUAN Chunsheng,LEI Yanyun,LI Shixiong,et al. Spontaneous combustion three zones determination and fire control technology in goaf under different ventilation volume [J]. Coal Science and Technology,2015,43(1):62-65.

[3]馬東,秦波濤.綜放工作而采空區(qū)注氮量與氧化自燃帶分布關(guān)系[J].煤炭科學(xué)技術(shù),2016,44 (4):78-82.

MA Dong,QIN Botao. Relationship between oxidized spontaneous combustion zone and nitrogen injection quantity in goaf of fully-mechanized top coal caving mining face [J]. Coal Science and Technology,2016,44(4):78-82.

[4]周心權(quán).煤礦采空區(qū)自燃防治和應(yīng)急處置能力提高的建議[J].煤炭科學(xué)技術(shù),2013,41(9):151-153.

ZHOU Xinquan. Proposals on improvement of spontaneous combustion prevention and control in mining goaf and emergency handing capacity [J]. Coal Science and Technology,2013,41(9):151-153.

[5]周西華,白剛,李誠玉,等.天池礦102綜放孤島工作面控制風(fēng)量防滅火技術(shù)研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2015,11(3):105-111.

ZHOU Xihua,BAI Gang,LI Chengyu,et al. Study on fire prevention and control technology by controlling air volume on 102 fully-mechanized isolated island caving face of Tianchi Coal Mine [J]. Journal of Safety Science and Technology,2015,11(3):105-111.

[6]曹鏡清,鄔劍明,周春山,等.低位放頂煤采空區(qū)自燃區(qū)域劃分與注氮口位置確定[J].煤炭科學(xué)技術(shù),2017,45(2):89-94.

CAO Jingqing,WU Jianming,ZHOU Chunshan,et al. Division of spontaneous combustion dangerous regionand determination of nitrogen injection position in goaf of low level top coal caving mining face [J]. Coal Science and Technology,2017,45(2):89-94.

[7]張龍,袁樹杰,朱成濤.俯采工作面采空區(qū)“三帶”分布研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2013,9(5):67-71.

ZHANG Long,YUAN Shujie,ZHU Chengtao. Study on distribution of “Three Zones” in goaf of full mechanized underhand mining face [J]. Journal of Safety Science and Technology,2013,9(5):67-71.

[8]蔣仲安.綜采工作面采空區(qū)封閉與惰化防滅火技術(shù)的數(shù)值模擬[J].北京科技大學(xué)學(xué)報,2014,36(6):722-729.

JIANG Zhongan. Numerical simulation of goaf sealing and inerting fire-fighting technology for fully mechanized coal faces [J]. Journal of University of Science and Technology Beijing,2014,36(6):722-729.

[9]吳玉國,鄔劍明,張東坡,等.綜放工作面連續(xù)注氮下采空區(qū)氣體分布及“三帶”變化規(guī)律[J].煤炭學(xué)報,2011,36(6):964-967.

WU Yuguo,WU Jianming,ZHANG Dongpo,et al. Distribution law of gas and change rule of “three zones” in the goaf of fully mechanized top-coal caving working face under the continuous nitrogen injection [J]. Journal of China Coal Society,2011,36(6):964-967.

[10]羅新榮,唐冠楚,李亞偉,等. CFD模型下采空區(qū)瓦斯抽采與注氮防滅火設(shè)計[J].中國安全科學(xué)學(xué)報,2016,26(6):69-74.

LUO Xinrong,TANG Guanchu,LI Yawei,et al. CFD model based design of both gas drainage and nitrogen injection against fire in goaf [J]. China Safety Science Journal,2016,26(6):69-74.

[11]朱紅青.頂煤自燃危險性分析及注氮防火的理論研究[J].煤炭學(xué)報,2012,37(6):1015-1020.

ZHU Hongqing. Investigation into spontaneous combustion hazard and numerical simulation of nitrogen injection fire prevention effects in the roof coal [J]. Journal of China Coal Society,2012,37(6):1015-1020.

[12]劉星魁,王公忠,周愛桃.采空區(qū)注氮參數(shù)設(shè)計及三維流場分析[J].太原理工大學(xué)學(xué)報,2015,46(3):312-317.

LIU Xingkui,WANG Gongzhong,ZHOU Aitao. Nitrogen injection parameters design for goaf and its 3D flow field analysis [J]. Journal of Taiyuan University of Technology,2015,46(3):312-317.

[13]尹曉雷,戴廣龍,吳彬,等.綜采面動態(tài)注氮作用下采空區(qū)“三帶”分布及防滅火技術(shù)研究[J].中國安全生產(chǎn)科學(xué)技術(shù),2014,10(10):137-142.

YIN Xiaolei,DAI Guanglong,WU Bin,et al. Research on distribution law of “three zones” in goaf of full-mechanized mining face and fire prevention technology under the action of dynamic nitrogen injection [J]. Journal of Safety Science and Technology,2014,10(10):137-142.

[14]高科,劉劍,劉玉姣.急傾斜綜放面采空區(qū)注氮方案研究及數(shù)值模擬[J].安全與環(huán)境學(xué)報,2014,14(5):49-53.

GAO Ke,LIU Jian,LIU Yujiao. Innovated nitrogen-injection project and numerical simulation in gob at steep inclined fully-mechanized mining face [J]. Journal of Safety and Environment,2014,14(5):49-53.

[15]黃顯華,顧北方,靳曉華.注氮條件對采空區(qū)自燃“三帶”的影響[J].科技導(dǎo)報,2014,32(18):26-30.

HUANG Xianhua,GU Beifang,JIN Xiaohua. Effect of nitrogen injection condition on the spontaneous combustion “three zones” in goafs [J]. Science and Technology Review,2014,32(18):26-30.

[16]S J.The determination of a co-efficent of longwall gob per-meability [J]. Archives of Mining Sciences,2001,46(4):451-468.

[17]汪文革,袁奎.注氮條件下瓦斯抽采對采空區(qū)自燃“三帶”的影響[J].煤炭科學(xué)技術(shù),2014,42 (12):75-78,83.

WANG Wenge,YUAN Kui. Gas drainage affected to spontaneous combustion three zones of mining goaf under the condition of nitrogen injection [J]. Coal Science and Technology,2014,42(12):75-78,83.

[18]周西華,門金龍,李誠玉,等.綜放孤島工作面采空區(qū)自燃與爆炸危險區(qū)監(jiān)測及數(shù)值模擬[J].安全與環(huán)境學(xué)報,2016,16(1):24-28.

ZHOU Xihua,MEN Jinlong,LI Chengyu,et al. Numerical simulation and the monitoring activities for the spontaneous combustion in the goaf of the fully mechanized cave-mining isolated working face [J]. Journal of Safety and Environment,2016,16(1):24-28.