上鄰近保護(hù)層開采Y型通風(fēng)采空區(qū)瓦斯分布規(guī)律及控制
——平煤六礦戊8-32010工作面案例研究

薛俊華，馬 騫，李延河，李洪彪，袁占棟

(1.西安科技大學(xué) 安全科學(xué)與工程學(xué)院，陜西 西安710054；2.平頂山天安煤業(yè)股份有限公司 六礦，河南 平頂山 467000)

0 引 言

高瓦斯煤層在井工開采過程中，工作面一直面臨著采空區(qū)和上隅角的瓦斯積聚問題[1-2]，尤其是近距離上保護(hù)層開采工作面，由于開采擾動(dòng)后下鄰近層卸壓瓦斯竄入采空區(qū)，導(dǎo)致采空區(qū)和上隅角瓦斯?jié)舛染痈卟幌拢瑖?yán)重妨礙著礦井安全生產(chǎn)。Y型通風(fēng)系統(tǒng)可以消除上隅角瓦斯積聚，解決上隅角瓦斯易于超限的問題，是一條實(shí)現(xiàn)安全、高效、經(jīng)濟(jì)一體化生產(chǎn)的有效途徑[3]。

雖然Y型通風(fēng)系統(tǒng)解決了上隅角瓦斯超限問題，但同時(shí)改變了采空區(qū)瓦斯?jié)舛葓鯷4-5]，所以在施工過程中不能套用以前U型通風(fēng)的瓦斯?jié)舛葓龇植家?guī)律以及瓦斯治理技術(shù)[6-7]。只有摸清楚瓦斯?jié)舛鹊姆植家?guī)律，找出適合該工作面的系統(tǒng)風(fēng)量及2條進(jìn)風(fēng)巷的風(fēng)量分配，確定出高濃度瓦斯區(qū)域，在留巷內(nèi)向此區(qū)域施工抽采鉆孔才是解決采空區(qū)瓦斯?jié)舛雀叩目茖W(xué)方法[8-9]。

為此，眾多學(xué)者從理論研究、數(shù)值模擬、實(shí)驗(yàn)室試驗(yàn)和現(xiàn)場試驗(yàn)等方面進(jìn)行了大量研究，獲得了Y型通風(fēng)采空區(qū)瓦斯?jié)舛葓龇植家?guī)律[10]，偏Y型通風(fēng)采空區(qū)瓦斯涌出規(guī)律[11]，開發(fā)了Y型通風(fēng)風(fēng)控瓦斯技術(shù)、留巷鉆孔法抽采瓦斯技術(shù)等[12]。但是由于理論研究將現(xiàn)場復(fù)雜的條件做了理想簡化及假設(shè)，一般不能很好的指導(dǎo)現(xiàn)場試驗(yàn)[13-14]。數(shù)值模擬試驗(yàn)在理論研究的基礎(chǔ)上也做了一些條件假設(shè)，會(huì)與現(xiàn)場的試驗(yàn)有些差異，但是其適用性較強(qiáng)[15-17]。只有結(jié)合理論研究和數(shù)值模擬試驗(yàn)，才可以較好的指導(dǎo)現(xiàn)場工業(yè)性實(shí)驗(yàn)。

為了研究平頂山礦區(qū)煤層群特征下的上鄰近保護(hù)層開采Y型通風(fēng)采空區(qū)瓦斯?jié)舛葓龇植家?guī)律及控制技術(shù)，文中基于理論研究，建立了采空區(qū)瓦斯涌出模型，并導(dǎo)入到FLUENT數(shù)值模型之中，研究了上鄰近保護(hù)層Y型通風(fēng)采空區(qū)瓦斯分布規(guī)律，并以進(jìn)風(fēng)巷風(fēng)量分配為自變量，研究了其對(duì)瓦斯分布規(guī)律的影響，結(jié)合現(xiàn)場試驗(yàn)，確定最合適的進(jìn)風(fēng)巷風(fēng)量分配，并在現(xiàn)場進(jìn)行了驗(yàn)證，對(duì)比分析了導(dǎo)致數(shù)值模擬試驗(yàn)和現(xiàn)場試驗(yàn)結(jié)果之間差異的原因，最后利用留巷鉆孔法對(duì)抽采效果進(jìn)行考察，以期能為同類型煤層條件下的采空區(qū)瓦斯治理提供一定的參考依據(jù)。

1 礦井及工作面概況

平頂山天安煤業(yè)股份有限公司六礦主采煤層從上至下分別為丁5-6煤層、戊8煤層和戊9-10煤層，丁5-6煤層和戊8煤層平均層間距71 m，戊8煤層和戊9-10煤層平均層間距5 m，經(jīng)煤層突出危險(xiǎn)性鑒定后，決定將戊8煤層作為首采層。

戊8-32010工作面位于戊8煤層，平均煤厚2.2 m，設(shè)計(jì)可采長度2 300 m，傾向長度220 m，標(biāo)高-633～-768 m，埋深873～938 m，通風(fēng)系統(tǒng)結(jié)合高位瓦斯抽采巷及機(jī)巷切頂留巷布置，采用兩進(jìn)一回Y型通風(fēng)方式，機(jī)巷、風(fēng)巷進(jìn)風(fēng)，機(jī)巷留巷和機(jī)巷高抽巷回風(fēng)，煤層為近水平煤層，遺煤平均厚度為0.2 m，可解吸瓦斯含量為4.3 m3/t，煤壁瓦斯初始涌出速率為0.003 56 m3/(m2·min)。工作面布置示意圖如圖1所示。

圖1 煤層層位關(guān)系及工作面布置示意Fig.1 Relation of coal layers and working face

2 采空區(qū)瓦斯?jié)舛确植家?guī)律及抽采數(shù)值模擬

2.1 條件假設(shè)

1)采空區(qū)為冒落矸石和巖塊構(gòu)成的帶有空隙的多孔介質(zhì)，并將多孔介質(zhì)視為各向同性。

2)采空區(qū)氣體不可壓縮，不產(chǎn)生化學(xué)反應(yīng)，其流動(dòng)近似為穩(wěn)態(tài)流動(dòng)、等溫過程。

3)忽略瓦斯的浮升效應(yīng)，瓦斯充滿整個(gè)采空區(qū)。

4)采空區(qū)瓦斯涌出源來自鄰近層和遺煤。

5)遺煤和鄰近層產(chǎn)生的瓦斯都是瞬間釋放的，忽略解吸時(shí)間。

2.2 模型建立及網(wǎng)格劃分

模擬對(duì)象為戊8-32010工作面，工作面長度220 m，采空區(qū)深度300 m，據(jù)現(xiàn)場考察，豎向裂隙帶高度穩(wěn)定在57 m，所以采空區(qū)高度設(shè)置為57 m，工作面寬度10 m，工作面高度4 m，風(fēng)巷、機(jī)巷長度20 m，寬度4 m，高度3.5 m，建成的模型如圖2所示。

圖2 戊8-32010工作面數(shù)值模型Fig.2 Numerical simulation of Fifth8-32010 working face

風(fēng)巷、機(jī)巷和留巷網(wǎng)格間距為1 m×1 m，采空區(qū)模塊網(wǎng)格尺寸為2 m×2 m，網(wǎng)格劃分結(jié)果如圖3所示。

圖3 網(wǎng)格劃分Fig.3 Net divsion

2.3 采空區(qū)瓦斯涌出源模型

2.3.1 采空區(qū)空隙率分布模型

采空區(qū)空隙率與滲透率分布規(guī)律符合“O”形圈覆巖垮落形態(tài)，采空區(qū)空隙率與滲透率分布的數(shù)學(xué)表達(dá)如公式(1)和(2)[18]。

n(x,y,z)=

(1)

(2)

式中，n(x，y，z)為空隙率；L為工作面長度，m；k(x，y，z)為滲透率，m2；Dp為采空區(qū)冒落巖塊平均粒徑，取250 mm[19].

2.3.2 采空區(qū)遺煤瓦斯放散模型

采空區(qū)不同位置處的瓦斯涌出速率qc(y)可用公式組(3)求得[20]。

(3)

式中a為瓦斯涌出初始強(qiáng)度，m3/min；b為瓦斯涌出衰減系數(shù)，min-1；y為遺煤位置到工作面的距離，m；Vx為工作面推進(jìn)速度，m/min；Qt為單位體積遺煤可解吸瓦斯放散體積，m3/m3；T為解吸時(shí)間，取288 000 min(即200 d).

通過上述方程組結(jié)合工作面遺煤放散條件，即可求出遺煤瓦斯涌出速率。

2.3.3 下鄰近層瓦斯涌出模型

戊8-32010工作面開采時(shí)，戊9-10卸壓瓦斯涌入工作面采空區(qū)，下鄰近層瓦斯涌出量可由公式(4)求得，其中工作面煤壁絕對(duì)瓦斯涌出量參數(shù)，可由公式(5)求得，瓦斯排放率可由公式(6)求得，最終得出下鄰近層的瓦斯涌出量[11]。

(4)

(5)

(6)

式中Qs為戊9-10煤層量瓦斯涌出量，m3/min；mi為戊9-10煤層厚度，3.1 m；Qb為戊8-32010工作面絕對(duì)瓦斯涌出量，m3/min；hi為鄰近層層間距，m；其余薄煤層瓦斯涌出量忽略不計(jì)；ηs為戊9-10煤層的瓦斯排放率；L為戊8-32010工作面長度，m；h為開采厚度，2.2 m；S為開采層與鄰近層層間距，13 m；q0為工作面煤壁剛暴露時(shí)單位面積煤壁的瓦斯涌出強(qiáng)度，m3/(m2·min)；V為工作面采煤機(jī)割煤時(shí)期平均牽引速度，m/min，實(shí)測值為10.08 m/min；Tc為工作面正常生產(chǎn)班期間，采煤機(jī)完成一次割煤周期所用時(shí)間，min，實(shí)測值為44 min。

2.4 采空區(qū)瓦斯運(yùn)移方程組選擇

選擇的氣體組分為甲烷-空氣混合氣體，選擇湍流方程、連續(xù)性方程和動(dòng)量守恒方程作為混合氣體在采空區(qū)內(nèi)的流動(dòng)方程。假設(shè)混氣擴(kuò)散符合Fick定律，其中混合氣體在采空區(qū)中的流動(dòng)方程的統(tǒng)一表達(dá)式如式(7)所示，F(xiàn)ick擴(kuò)散定律的表達(dá)式如式(8)所示[21]。

(7)

式中ρ為氣體密度，kg/m3；φ為通用變量；t為時(shí)間，s；Γ為廣義擴(kuò)散系數(shù)；S為廣義源項(xiàng)。

(8)

式中r為極坐標(biāo)半徑，m；D為擴(kuò)散系數(shù)，m2/s；c為瓦斯?jié)舛龋琺ol/m3。

2.5 邊界條件及相關(guān)參數(shù)設(shè)置

為研究系統(tǒng)風(fēng)量為2 830 m3/min時(shí)，兩進(jìn)風(fēng)巷進(jìn)風(fēng)配比對(duì)工作面及留巷的瓦斯?jié)舛鹊挠绊懀苍O(shè)置了2個(gè)進(jìn)風(fēng)入口和一個(gè)出口。2個(gè)進(jìn)風(fēng)入口，分別表示戊8工作面機(jī)巷和風(fēng)巷，工作面機(jī)巷風(fēng)量與機(jī)巷風(fēng)量之比分別為：1∶3，1∶2和1∶1.5，一個(gè)出口設(shè)置為自然出流。

2.6 抽采條件設(shè)置

在留巷向里距離工作面切頂線10 m處增加抽放匯對(duì)采空區(qū)瓦斯進(jìn)行抽采，抽采管直徑為500 mm，鉆孔高度為巷道底部2 m處。

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

3.1 瓦斯?jié)舛确植家?guī)律實(shí)驗(yàn)結(jié)果及分析

在總風(fēng)量2 830 m3/min的基礎(chǔ)上，對(duì)工作面機(jī)、風(fēng)巷風(fēng)量分配進(jìn)行調(diào)整試驗(yàn)，分別對(duì)機(jī)巷、風(fēng)巷風(fēng)量比例1∶3，1∶2和1∶1.53種狀態(tài)對(duì)工作面瓦斯質(zhì)量分?jǐn)?shù)的影響進(jìn)行研究，結(jié)果如圖4所示。

圖4 采空區(qū)瓦斯?jié)舛葓瞿M結(jié)果Fig.4 Simulation results of gas concentration field in goaf

由圖4可得，配風(fēng)比例為1∶3時(shí)，通過工作面的風(fēng)量較大，上下隅角壓差增大，風(fēng)巷端口與留巷壓差也隨之增大，導(dǎo)致工作面漏風(fēng)嚴(yán)重，造成采空區(qū)防滅火難度加大；配風(fēng)量為1∶1.5時(shí)，采空區(qū)和留巷處瓦斯?jié)舛葏^(qū)域面積明顯增大，基本可以看出增大風(fēng)巷風(fēng)量比例，工作面漏風(fēng)加劇，減少風(fēng)巷風(fēng)量，采空區(qū)和留巷內(nèi)瓦斯?jié)舛葧?huì)增大，而當(dāng)采用配風(fēng)比例1∶2時(shí)，一方面可以令工作面漏風(fēng)情況得以控制，另一方面也可以緩解采空區(qū)和留巷內(nèi)瓦斯積聚，所以配風(fēng)比例為1∶2時(shí)更適合于戊8-32010工作面。

在Y型通風(fēng)條件下，瓦斯向采空區(qū)深部不斷積聚，且采空區(qū)尾部風(fēng)巷側(cè)瓦斯體積分?jǐn)?shù)高于尾部留巷側(cè)，由卸壓鉆孔法抽采卸壓瓦斯機(jī)理可知[12]，采空區(qū)和覆巖裂隙區(qū)相連通，在留巷后部向采空區(qū)布置瓦斯抽采鉆孔，可以持續(xù)高濃度抽采瓦斯資源。

3.2 瓦斯抽采實(shí)驗(yàn)結(jié)果及分析

采用配風(fēng)比例為1∶2的通風(fēng)模式建立戊8-32010工作面采空區(qū)抽采模型并布置尾抽鉆孔，抽采模型如圖5所示，抽采結(jié)果如圖6所示。

圖5 戊8-32010工作面抽采模型Fig.5 Extraction model of Fifth8-32010 working face

圖6 戊8-32010工作面瓦斯抽采模擬結(jié)果Fig.6 Simulation results of gas extraction at Fifth8-32010 working face

由圖6可以看出，利用尾抽鉆孔抽采瓦斯后，采空區(qū)瓦斯體積分?jǐn)?shù)明顯降低，說明在留巷深部向采空區(qū)布置鉆孔可以高效抽采瓦斯。

4 討 論

(9)

式中Ma為空氣摩爾質(zhì)量，kg/mol；Mg為瓦斯摩爾質(zhì)量，kg/mol。

在總風(fēng)量2 830 m3/min的基礎(chǔ)上，對(duì)工作面機(jī)、風(fēng)巷風(fēng)量分配進(jìn)行調(diào)整試驗(yàn)，分別對(duì)機(jī)巷、風(fēng)巷風(fēng)量比例1∶3，1∶2和1∶1.5這3種狀態(tài)時(shí)，對(duì)工作面3個(gè)測點(diǎn)和留巷內(nèi)2個(gè)測點(diǎn)共5個(gè)測點(diǎn)的瓦斯?jié)舛冗M(jìn)行對(duì)比分析，測點(diǎn)布置如圖7所示，現(xiàn)場試驗(yàn)測點(diǎn)與對(duì)應(yīng)測試位置見表1。

圖7 瓦斯?jié)舛葴y點(diǎn)布置Fig.7 Layout of gas concentration measuring dots

表1 試驗(yàn)測點(diǎn)與對(duì)應(yīng)測試位置

先將現(xiàn)場試驗(yàn)和數(shù)值模擬實(shí)驗(yàn)5個(gè)測點(diǎn)位置的瓦斯體積分?jǐn)?shù)統(tǒng)計(jì)、并繪圖，瓦斯體積濃度對(duì)比結(jié)果如圖8所示。

由圖8(a)可得，機(jī)、風(fēng)巷風(fēng)量配比1∶3時(shí)工作面瓦斯?jié)舛仍?.16%～0.3%之間，留巷段瓦斯?jié)舛仍?.32%～0.69%之間；配比為1∶1.5時(shí)工作面瓦斯?jié)舛仍?.24%～0.5%之間，留巷段瓦斯?jié)舛仍?.51%～0.64%之間。由圖8(b)可得，機(jī)、風(fēng)巷風(fēng)量配比1∶3時(shí)工作面瓦斯?jié)舛仍?.02%～0.05%之間，留巷段瓦斯?jié)舛仍?.05%～0.3%之間；配比為1∶1.5時(shí)工作面瓦斯?jié)舛仍?.05%～0.07%之間，留巷段瓦斯?jié)舛仍?.04%～0.28%之間。

綜上可得，風(fēng)巷風(fēng)量比例越小，工作面瓦斯?jié)舛仍礁撸粝锒瓮咚節(jié)舛仍礁撸伙L(fēng)巷風(fēng)量比例越大，工作面瓦斯?jié)舛仍降停F(xiàn)場試驗(yàn)和數(shù)值模擬試驗(yàn)采空區(qū)瓦斯?jié)舛鹊淖兓厔莼疽恢隆５ぷ髅嫣幒土粝飪?nèi)瓦斯?jié)舛戎挡町愝^大，造成這一差異的主要原因是數(shù)值模擬中無法控制工作面向采空區(qū)內(nèi)漏風(fēng)和采空區(qū)向留巷內(nèi)漏風(fēng)的過程，后續(xù)筆者會(huì)就這一問題做進(jìn)一步研究。

最后在留巷段向里距離工作面切頂線10 m處施工了瓦斯尾抽鉆孔，抽采了采空區(qū)瓦斯，在此期間，工作面和采空區(qū)瓦斯都處于安全值以內(nèi)[23-24]，相比于上個(gè)相鄰U型工作面推進(jìn)速度，本工作面推進(jìn)速度明顯提高，如圖9為日產(chǎn)量與抽采濃度之間的關(guān)系圖。

圖8 瓦斯體積濃度對(duì)比Fig.8 Comparison of gas volume and itss concentration

圖9 工作面日產(chǎn)量與瓦斯?jié)舛鹊年P(guān)系Fig.9 Relationship between day-products and gas concentration in working face

從圖9可以看出，利用尾抽抽采采空區(qū)尾部瓦斯，瓦斯?jié)舛榷荚?0%以上，說明在尾部布置瓦斯抽采效果較好。在抽采效果考察后期，日產(chǎn)量有短時(shí)間的相對(duì)減小，這是由于此時(shí)煤層上覆巖層破斷垮落，裂隙帶初步形成，將尾抽為主的抽采方式轉(zhuǎn)化為高抽巷裂隙帶瓦斯為主抽采方式，一部分瓦斯通過裂隙帶進(jìn)行抽采，暫時(shí)降低了日產(chǎn)量。

5 結(jié) 論

1)在采空區(qū)瓦斯涌出和下鄰近層瓦斯涌出強(qiáng)度一定時(shí)，近距離上保護(hù)層Y型通風(fēng)的瓦斯易于積聚在采空區(qū)尾部，越靠近尾部，瓦斯體積分?jǐn)?shù)增加速率越大，采空區(qū)尾部留巷側(cè)瓦斯體積分?jǐn)?shù)低于尾部風(fēng)巷側(cè)，在留巷尾部一定位置處布置瓦斯抽采鉆孔能高效抽采瓦斯，并且此濃度場與單煤層Y型通風(fēng)瓦斯?jié)舛葓龇植家?guī)律基本一致。

2)增大風(fēng)巷風(fēng)量比例，會(huì)導(dǎo)致采空區(qū)瓦斯大量涌入留巷；減小風(fēng)巷風(fēng)量比例，會(huì)導(dǎo)致工作面瓦斯增加，因此改變2條進(jìn)風(fēng)巷配風(fēng)量，工作面和留巷的瓦斯?jié)舛葧?huì)隨之變化。