冀 凱 戴子棟
(河南龍宇能源股份有限公司陳四樓煤礦,河南 商丘 476600)
采空區(qū)遺煤易引起自然發(fā)火事故,采空區(qū)內(nèi)部火源點(diǎn)難以定位且隱蔽,是困擾礦井安全生產(chǎn)的主要災(zāi)害之一[1]。常規(guī)探測(cè)技術(shù)精確度不足,難以探測(cè)采空區(qū)深部。煤自燃“三帶”劃分是防止采空區(qū)自然發(fā)火最有效方法之一[2],精準(zhǔn)劃分采空區(qū)煤自燃“三帶”范圍對(duì)礦井各種防滅火措施的制定具有指導(dǎo)作用,極大促進(jìn)礦井安全開(kāi)采。
陳四樓煤礦21210 工作面煤層的厚度為0.6~4.8 m,平均厚度2.7 m。該工作面布置規(guī)整,適宜工作面的開(kāi)采。工作面中部過(guò)一處夾矸發(fā)育區(qū),結(jié)構(gòu)復(fù)雜,為較穩(wěn)定煤層。21210 工作面整體為一傾向NW 的單斜構(gòu)造,地層傾角為20°~26°,平均23°,采用綜采工藝具有較好的經(jīng)濟(jì)效益。正常情況下跟頂回采,不得任意留頂?shù)酌骸C簩雍穸刃∮?.7 m 時(shí),跟頂破底回采,采高不低于2.7 m。煤層厚度不小于2.7 m,不大于4 m 時(shí),跟頂跟底回采。當(dāng)煤層厚度大于4 m 時(shí),跟頂留底煤回采,采高不得超過(guò)4 m。
采煤工作面采空區(qū)自然發(fā)火“三帶”可劃分為散熱帶、氧化帶和窒息帶[3]。目前,劃分采空區(qū)自燃“三帶”的主要方法有:溫度法、漏風(fēng)強(qiáng)度法及氧氣濃度法[4]。溫度法劃分采空區(qū)自燃“三帶”有一定缺陷,采空區(qū)某處溫度較低時(shí)該處的遺留仍存在自然發(fā)火的可能性,不能全面衡量采空區(qū)的自燃危險(xiǎn)性;由于受測(cè)量?jī)x器的影響以及采空區(qū)風(fēng)流方向和大小的復(fù)雜性,漏風(fēng)強(qiáng)度法劃分采空區(qū)自燃“三帶”在井下現(xiàn)場(chǎng)應(yīng)用較為困難。因此,溫度法、漏風(fēng)強(qiáng)度法都不太適用。氧氣濃度對(duì)采空區(qū)煤炭自燃起決定性作用,目前普遍采用氧氣濃度法來(lái)劃分采空區(qū)自燃“三帶”,具有現(xiàn)場(chǎng)可操作性強(qiáng)、劃分指標(biāo)明確的優(yōu)點(diǎn)。因此,陳四樓煤礦21210 工作面選擇氧氣濃度法測(cè)試采空區(qū)自燃“三帶”分布。
21210 工作面采用“U”型通風(fēng)方式,進(jìn)風(fēng)巷道、回風(fēng)巷道各布置2 個(gè)采樣點(diǎn)。采樣點(diǎn)1 位于工作面進(jìn)風(fēng)隅角端頭,采樣點(diǎn)2 與采樣點(diǎn)1 間距約30 m。每個(gè)測(cè)點(diǎn)連接1 根束管,束管鋪設(shè)長(zhǎng)度各為150 m,束管布置于進(jìn)風(fēng)巷道緊貼巷幫一側(cè)。同樣的,回風(fēng)巷道采樣點(diǎn)布置與進(jìn)風(fēng)巷道采樣點(diǎn)布置距離、位置一致,具體如圖1。隨著工作面不斷推進(jìn),各采樣點(diǎn)逐漸進(jìn)入采空區(qū)內(nèi)部,在氣樣采集點(diǎn)抽取各束管氣體,即可監(jiān)測(cè)采空區(qū)不同位置的氣體成分。
隨著工作面的推進(jìn),1、3 號(hào)測(cè)點(diǎn)首先進(jìn)入采空區(qū),2、4 號(hào)測(cè)點(diǎn)后進(jìn)入采空區(qū)。對(duì)進(jìn)入采空區(qū)的測(cè)點(diǎn)現(xiàn)場(chǎng)取樣,周期為每天采樣一次,用實(shí)驗(yàn)室色譜分析儀對(duì)樣品進(jìn)行檢測(cè)分析。根據(jù)分析結(jié)果,將陳四樓煤礦采空區(qū)內(nèi)氧氣體積分?jǐn)?shù)距工作面的距離變化繪制成圖,采樣點(diǎn)1~采樣點(diǎn)4 的變化趨勢(shì)如圖2、圖3。
圖2 采樣點(diǎn)1、采樣點(diǎn)2 氧氣濃度變化曲線圖
圖3 采樣點(diǎn)3、采樣點(diǎn)4 氧氣濃度變化曲線圖
實(shí)際測(cè)定中,劃分散熱帶與氧化帶、氧化帶與窒息帶的氧氣濃度指標(biāo)分別是18%和5%,即氧氣濃度范圍大于18%為散熱帶,氧氣濃度范圍介于5%~18%之間為氧化自燃帶,氧氣濃度范圍低于5%的為窒息帶。依據(jù)O2濃度指標(biāo),21210 綜采工作面實(shí)測(cè)自燃“三帶”范圍見(jiàn)表1。
表1 21210 工作面采空區(qū)自燃“三帶”分布 m
由表1 的采樣點(diǎn)1、采樣點(diǎn)2 監(jiān)測(cè)結(jié)果可知,21210 工作面采空區(qū)進(jìn)風(fēng)側(cè)的散熱帶界線為采空區(qū)深部75~81 m,平均寬度為78 m;氧化自燃帶界線為采空區(qū)深部126~132 m,平均129 m;進(jìn)入采空區(qū)126~132 m 之后,即為窒息帶范圍。
由表1 的采樣點(diǎn)3、采樣點(diǎn)4 監(jiān)測(cè)結(jié)果可知,21210 工作面采空區(qū)回風(fēng)側(cè)的散熱帶界線為采空區(qū)深部41~44 m,平均寬度為42.5 m;氧化自燃帶界線為采空區(qū)深部90~95 m,平均92.5 m;進(jìn)入采空區(qū)90~95 m 之后,即為窒息帶范圍。
按照現(xiàn)場(chǎng)實(shí)際及數(shù)值模擬的要求,以陳四樓煤礦21210 工作面采空區(qū)實(shí)際尺寸及風(fēng)量為依據(jù),為了簡(jiǎn)化問(wèn)題、把握主要因素,將巷道、工作面以及采空區(qū)均視為長(zhǎng)方體,采空區(qū)長(zhǎng)×寬×高為450 m×193 m×2.9 m。在笛卡爾坐標(biāo)系下創(chuàng)建3D 物理模型,建立一個(gè)長(zhǎng)寬高為450 m×193 m×2.9 m 的模型,使用Ansys Mesh 組件對(duì)采空區(qū)三維模型進(jìn)行網(wǎng)格劃分,模型共有618.93 萬(wàn)個(gè)網(wǎng)格,網(wǎng)格步長(zhǎng)取0.5 m。采空區(qū)物理模型及網(wǎng)格劃分如圖4。
圖4 采空區(qū)三維物理模型網(wǎng)格圖
計(jì)算模型選擇穩(wěn)態(tài)模型(steady),壓力與速度之間的耦合選用基于交錯(cuò)網(wǎng)格的SIMPLE 算法[5]。模擬運(yùn)算中,進(jìn)風(fēng)巷道面邊界設(shè)置為速度入口,總進(jìn)風(fēng)為1176 m3/min,回風(fēng)巷道面邊界設(shè)置為壓力出口,視研究范圍內(nèi)的氣體為理想氣體,氧氣濃度為21%,其余均為氮?dú)?。巷道及采空區(qū)均為長(zhǎng)方體,內(nèi)部尺寸各有差異。為了便于計(jì)算,模擬期間各種物理量恒定,不受采空區(qū)推進(jìn)影響。模擬中考慮重力因素,大小為9.8 N/kg。
由圖5 可知,在“U”型通風(fēng)方式下,采空區(qū)高氧濃度區(qū)域近似呈三角形式分布。采空區(qū)進(jìn)風(fēng)側(cè)的高氧濃度范圍大于回風(fēng)側(cè),距工作面約88 m 的范圍內(nèi)氧氣濃度均高于18%。而在采空區(qū)回風(fēng)側(cè),采空區(qū)漏風(fēng)逐漸向工作面回風(fēng)隅角匯入,氧氣濃度分布范圍減小,逐漸向工作面方向收縮??傮w上,隨著進(jìn)入采空區(qū)深部,離工作面越遠(yuǎn),氧氣濃度越低,進(jìn)風(fēng)側(cè)氧氣濃度高于7%的范圍為0~122 m,回風(fēng)側(cè)氧氣濃度高于7%的范圍為0~85 m;沿采空區(qū)傾向,氧氣濃度在18%~7%的采空區(qū)范圍成長(zhǎng)條狀,且貫穿整個(gè)傾向。
圖5 “U”型通風(fēng)方式下采空區(qū)氧氣濃度分布平面圖
1)通過(guò)井下束管取樣,對(duì)21210 綜采工作面采空區(qū)氧氣濃度不間斷監(jiān)測(cè)分析,利用氧氣濃度指標(biāo)劃分21210 工作面采空區(qū)自燃“三帶”范圍。進(jìn)風(fēng)測(cè)范圍:散熱帶0~78 m,氧化升溫帶78~129 m,窒息帶大于129 m;回風(fēng)側(cè)范圍:散熱帶0~42.5 m,氧化升溫帶42.5~92.5 m,窒息帶大于92.5 m。數(shù)值模擬結(jié)果與現(xiàn)場(chǎng)測(cè)定的采空區(qū)自燃“三帶”范圍基本一致,差異在合理范圍內(nèi)。數(shù)值模擬劃分的采空區(qū)自燃“三帶”為,進(jìn)風(fēng)側(cè)范圍:散熱帶0~88 m,氧化升溫帶88~122 m,窒息帶大于122 m;回風(fēng)側(cè)范圍:散熱帶0~30 m,氧化升溫帶30~85 m,窒息帶大于85 m。通過(guò)對(duì)比,模擬結(jié)果比較真實(shí)地反映采空區(qū)的情況。
2)數(shù)值模擬分析采空區(qū)氧氣濃度變化情況,進(jìn)風(fēng)側(cè)的氧化自燃帶寬度較大,回風(fēng)側(cè)的氧化自燃帶寬度較小,模擬與現(xiàn)場(chǎng)實(shí)際吻合。