盧義玉　夏彬偉　郭世杰（1.重慶大學(xué)煤礦災(zāi)害動力學(xué)與控制國家重點實驗室，重慶400030；2.重慶大學(xué)復(fù)雜煤氣層瓦斯抽采國家地方聯(lián)合工程實驗室，重慶400030）

?

堅硬頂板運(yùn)動對特厚煤層工作面瓦斯涌出的影響研究

盧義玉［1，2］夏彬偉［1，2］郭世杰［1，2］
（1.重慶大學(xué)煤礦災(zāi)害動力學(xué)與控制國家重點實驗室，重慶400030；2.重慶大學(xué)復(fù)雜煤氣層瓦斯抽采國家地方聯(lián)合工程實驗室，重慶400030）

摘要針對頂板運(yùn)動瓦斯異常涌出的問題，基于空氣動力學(xué)原理及牛頓第二定律，分析瓦斯涌出與頂板載荷、采空區(qū)高度、頂板面積等影響因素的關(guān)系，建立了頂板垮落—瓦斯涌出的耦合模型，借助COMSOL多物理場耦合軟件，以大同塔山煤礦二盤區(qū)8214工作面為例，根據(jù)其頂板垮落特征分析了上覆頂板分層垮斷下的采空區(qū)瓦斯涌出規(guī)律，并通過現(xiàn)場監(jiān)測數(shù)據(jù)進(jìn)行了對比分析。研究表明：頂板垮斷時，采空區(qū)瓦斯大部分將涌入工作面上隅角和高抽巷，其瓦斯?jié)舛瘸尸F(xiàn)短時間內(nèi)急劇增大后減小的變化，且涌出量的大小與頂板載荷、采空區(qū)高度、頂板面積等有著直接關(guān)系；頂板垮落時上層比下層造成更大量的瓦斯涌出，且上隅角瓦斯?jié)舛茸兓葹槊黠@。

關(guān)鍵詞特厚煤層；堅硬頂板；頂板垮落；工作面來壓；瓦斯涌出

“十三五”期間，我國煤炭行業(yè)將繼續(xù)向整合兼并重組方向調(diào)整，特大型煤炭生產(chǎn)基地將占主導(dǎo)地位，特別是特厚煤層礦區(qū)。以大同塔山礦為代表的特厚煤層綜放開采過程中，采空區(qū)空間大，頂板堅硬在回采過程中容易出現(xiàn)頂板整體垮落致使采空區(qū)瓦斯瞬間涌出造成瓦斯災(zāi)害，嚴(yán)重威脅著礦井的安全生產(chǎn)。因此探明頂板垮落下的采空區(qū)瓦斯涌出運(yùn)移規(guī)律對礦井的安全生產(chǎn)有著重要意義。

目前國內(nèi)外學(xué)者針對采動圍巖活動下的瓦斯涌出進(jìn)行了一些研究，學(xué)者認(rèn)為堅硬頂板控制下的特厚煤層在回采過程中，容易造成頂板的整體性垮落，使采空區(qū)氣體在沖擊作用下涌入工作面，造成工作面瓦斯災(zāi)害［1-3］。翁明月［4］對屯留煤礦S2205綜放工作面進(jìn)行了煤巖破壞的微震監(jiān)測、工作面礦壓顯現(xiàn)監(jiān)測和瓦斯涌出的實測研究，得到了綜放工作面煤巖破壞、礦壓顯現(xiàn)和瓦斯涌出三者之間的關(guān)系；李化敏等［5］以平煤13礦12071工作面為例，對支架支護(hù)阻力、超前支承壓力以及工作面瓦斯體積分?jǐn)?shù)進(jìn)行了實時監(jiān)測，分析了采動過程與瓦斯涌出之間的關(guān)系；顧鐵風(fēng)［6］建立了封閉采空區(qū)條件下颶風(fēng)災(zāi)害的理論模型，給出了差分解法；息金波［7］推導(dǎo)了長壁采煤法懸頂垮落沖擊力的計算公式；然而目前針對工作面瓦斯多以常態(tài)環(huán)境下的工作面靜態(tài)流場研究為主，對瓦斯異常涌出的瓦斯來源分析較少，且頂板垮落多以氣體沖擊災(zāi)害為研究對象，沒有對頂板垮落下的動態(tài)瓦斯涌出運(yùn)移規(guī)律進(jìn)行研究。

本文針對特厚煤層堅硬頂板整體垮落造成瓦斯異常涌出的問題，結(jié)合牛頓第二定律及空氣動力學(xué)原理，分析瓦斯涌出與頂板載荷、采空區(qū)高度、頂板面積等影響因素的關(guān)系，建立了頂板垮落—瓦斯涌出的耦合模型，借助comsol多物理場耦合軟件，以大同塔山礦生產(chǎn)工作面為例，對上覆堅硬巖層分層垮斷下的采空區(qū)瓦斯涌出規(guī)律進(jìn)行解算分析，結(jié)合現(xiàn)場監(jiān)測數(shù)據(jù)進(jìn)行驗證，為防治礦井瓦斯災(zāi)害提供依據(jù)。

1　頂板垮斷采空區(qū)瓦斯涌出模型研究

采空區(qū)頂板垮落沖擊是導(dǎo)致采空區(qū)（氣體）瓦斯涌出的主要原因［4］，在壓縮過程中，使采空區(qū)氣體壓力逐漸上升，造成采空區(qū)瓦斯瞬間涌向工作面，為了研究采空區(qū)頂板垮落下的氣體涌出，建立如2頁圖1頂板垮落氣體涌出模型并做出以下假設(shè)：（1）頂板堅硬完整性好，在回采過程中會出現(xiàn)頂板的整體垮落，造成采空區(qū)瓦斯在短時間內(nèi)涌向工作面及回風(fēng)巷；（2）將采空區(qū)視為規(guī)則長方體容器，采空區(qū)深部已被壓實，采空區(qū)內(nèi)空氣為理想氣體，頂板垮落氣體變化為等溫過程；（3）頂板垮落為頂板載荷及采空區(qū)氣體控制的運(yùn)動；（4）采空區(qū)氣體匯聚進(jìn)風(fēng)巷氣體流入回風(fēng)巷。

圖1　頂板垮落氣體涌出模型

頂板大面積垮落前，容器內(nèi)的氣體壓力為P0，采空區(qū)頂板面積A0，空氣的流速為v0，設(shè)定采空區(qū)空氣密度為ρ0。則：頂板重力為G=ρdabhg采空區(qū)氣體作用在頂板向上的力為F=abPZ，其中，ρd為上浮巖層塊體密度kg/m3，a，b，h分別為塊體的長寬高，m，PZ為塊體下降過程中采空區(qū)氣體的壓強(qiáng)，Pa，z為塊體下降高度，m，則頂板下降過程中，有G-F=maZ，則頂板下落加速度：

頂板在下降距離為z時，有氣體壓縮方程可知，P0V0=PZVZ，其中P0、V0為頂板垮落前采空區(qū)氣體壓強(qiáng)及體積，V0=Hab，下降z時VZ=（H-Z）ab，則有

頂板垮落過程中，在t時刻的壓縮體積為：

假設(shè)頂板壓縮氣體全部有回風(fēng)巷進(jìn)行排出，采空區(qū)瓦斯分布均勻，濃度為c1，進(jìn)風(fēng)巷為無瓦斯風(fēng)流，斷面面積A1，風(fēng)流速度v1，回風(fēng)巷巷道斷面面積A2，風(fēng)流速度風(fēng)速為v2，采空區(qū)壓縮氣體與進(jìn)風(fēng)巷風(fēng)流混合后流向回風(fēng)巷，假設(shè)采空區(qū)及回風(fēng)巷氣體密度相同，則根據(jù)可壓縮氣體的連續(xù)性（質(zhì)量守恒）方程在時間t時刻，流經(jīng)回風(fēng)巷氣體體積為：

則涌出瓦斯總量為VZc1，回風(fēng)巷瓦斯?jié)舛葹閏2，則：

將式（4）帶入（6）可得（7）：

公式中，頂板垮落過程是一個變加速過程，其速度是與頂板面積、采空區(qū)高度、頂板載荷、頂板下落位置以及氣體壓強(qiáng)5個因素相關(guān)的函數(shù)，想要直接求出vZ再求出瓦斯涌出情況是十分困難的，可通過數(shù)值模擬的方法進(jìn)行計算分析。

2　堅硬頂板垮落下的采空區(qū)瓦斯運(yùn)移模擬

本文以大同塔山煤礦堅硬頂板綜放工作面的布置情況為例，建立模型分析頂板垮落沖擊作用下的采空區(qū)瓦斯運(yùn)移。

2.1煤層及瓦斯涌出概況

塔山煤礦設(shè)計生產(chǎn)能力1 500萬t/a，目前開采的3-5#煤層平均煤厚12 m，屬特厚煤層，煤層瓦斯含量1.78 m3/t。礦井采用綜放低位放頂煤一次采全高的方法，采煤厚度3.5 m，放煤厚度8.5 m，通風(fēng)采用一進(jìn)兩回的通風(fēng)方式，進(jìn)風(fēng)風(fēng)量為3 200 m3/min，回風(fēng)風(fēng)量為2 200 m3/min，高抽巷風(fēng)量為1 000 m3/min。當(dāng)前礦井工作面正常回采期間絕對瓦斯涌出量保持在30 m3/ min～50 m3/min之間，瓦斯涌出異常期間，工作面絕對瓦斯涌出量可達(dá)到50 m3/min～75 m3/min之間，甚至更高。3-5#煤層頂板中存在近距離且厚度相差較大的堅硬頂板，下層堅硬頂板厚度4 m，上層堅硬頂板厚度20 m，兩層堅硬頂板之間夾4 m厚的2#煤層。堅硬巖層的普氏硬度為Ⅳ甲類，屬相當(dāng)堅硬巖石，頂?shù)装鍘r層柱狀圖如圖2所示。

圖2　巖層柱狀

2.2堅硬頂板垮落規(guī)律

通過平板理論［8-9］計算可知，塔山下層堅硬頂板斷裂步距為37.8 m，上層堅硬頂板斷裂步距為69.4 m，由現(xiàn)場檢測可知，工作面推進(jìn)到35 m～40 m時，下層堅硬頂板垮落，推進(jìn)到60 m～70 m左右時，兩層堅硬頂板可能呈現(xiàn)整體回轉(zhuǎn)垮落現(xiàn)象，并有明顯的大小來壓現(xiàn)象。頂板斷裂步距與大小來壓現(xiàn)象與現(xiàn)場觀測均較為一致。（見圖3）

圖3　支架工作阻力曲線

2.3堅硬頂板分層破斷頂板載荷計算

由煤層柱狀圖2，根據(jù)由關(guān)鍵層載荷［10］計算公式：

式中E1——第i層的彈性模量；hi——第i層的厚度；γi——第i層的容重；n——關(guān)鍵層所控巖層數(shù)；即q為全部（n層）控制巖層的自重載荷。則頂板巖層的力學(xué)模型如圖4。

圖4　頂板巖層的力學(xué)模型

其中q1為下層堅硬頂板所受單位面積載荷，q2為上層堅硬頂板所受單位面積載荷。根據(jù)巖層柱狀圖及關(guān)鍵層理論知上堅硬頂板控制高度達(dá)54.1 m，至砂質(zhì)泥巖層，下堅硬頂板控制高度即為2#煤層厚度，為4 m。

由煤巖層力學(xué)參數(shù)通過公式（8）計算出q1=0.158 MPa，q2=1.85 MPa。

2.4計算模型及參數(shù)

根據(jù)煤層頂板巖性以及垮斷步距運(yùn)用comsol軟件計算塔山煤礦8214綜放工作面采空區(qū)三維狀態(tài)下頂板垮落后采空區(qū)多孔介質(zhì)孔隙瓦斯壓力變化及瓦斯?jié)舛茸兓哼\(yùn)濤等［10］建立了垮落帶介質(zhì)孔隙率和滲透率的非均勻連續(xù)分布模型，分析可知，采空區(qū)初始孔隙率設(shè)置為0.4，由1.2章分析可知，設(shè)下層堅硬頂板垮斷步距40 m，上層堅硬頂板垮斷步距70 m，根據(jù)現(xiàn)場實際建立模型，模型只考慮采空區(qū)瓦斯涌出情況，忽略工作面煤壁瓦斯涌出狀況：下層堅硬頂板破斷（圖5）：頂板垮落后作用于采空區(qū)上部載荷為0.158 MPa，采空區(qū)高度12 m，工作面進(jìn)風(fēng)3 200 m3/min；上層堅硬頂板破斷（圖6）：頂板垮斷后作用于采空區(qū)上部載荷為1.85 MPa，采空區(qū)高度20 m，工作面進(jìn)風(fēng)3 200 m3/min。

圖5　下層堅硬頂板垮落作用采空區(qū)計算模型

圖6　上層堅硬頂板垮落作用采空區(qū)計算模型

2.5計算結(jié)果及分析

通過分析可知，頂板垮落的沖擊力主要是在一定時間內(nèi)改變采空區(qū)的瓦斯流場狀態(tài)，而采空區(qū)瓦斯流場狀態(tài)決定著回風(fēng)巷以及上隅角瓦斯?jié)舛鹊淖兓?，所以對采空區(qū)在頂板垮落沖擊力作用下的瓦斯流場進(jìn)行研究有著重要意義。根據(jù)以上假設(shè)條件以及控制方程，通過數(shù)值計算可以得到采空區(qū)在頂板垮落影響下的不同時間的瓦斯?jié)舛确植紙D，如圖7所示。

2.5.1下層堅硬頂板垮斷采空區(qū)瓦斯涌出計算結(jié)果

圖7　計算模型在不同時刻的瓦斯?jié)舛确植?/p>

圖8　頂板垮落沖擊作用后高抽巷、上隅角及回風(fēng)巷瓦斯?jié)舛茸兓€

2.5.2上層堅硬頂板垮斷采空區(qū)瓦斯涌出計算結(jié)果

圖9　計算模型在不同時刻的瓦斯?jié)舛确植?/p>

圖10　頂板垮落沖擊作用后高抽巷、上隅角及回風(fēng)巷瓦斯?jié)舛茸兓€

由圖7、圖9可以看出頂板垮落瞬間對采空區(qū)造成擠壓，使采空區(qū)瓦斯在風(fēng)流作用下向回風(fēng)巷以及高抽巷一側(cè)運(yùn)移，由于上隅角處于風(fēng)流盲區(qū)，在無采取防超措施下瓦斯?jié)舛茸钕仍龃?，結(jié)合圖8、圖10瓦斯?jié)舛茸兓€，可以看出，頂板垮落造成瓦斯?jié)舛仍诙虝r間內(nèi)急劇增大然后減小，與頂板的垮落呈現(xiàn)一定的相關(guān)性，在頂板垮落8 s～12 s時間內(nèi)瓦斯?jié)舛冗_(dá)到最高，隨即迅速降低，維持時間較短；回風(fēng)巷瓦斯?jié)舛茸兓雍笊嫌缃牵?0 s范圍內(nèi)達(dá)到最高，由于風(fēng)流稀釋作用，濃度低于上隅角；而高抽巷瓦斯?jié)舛戎饾u增大，在頂板垮落47 s～53 s時間內(nèi)瓦斯?jié)舛冗_(dá)到最高，高抽巷瓦斯?jié)舛绕毡楸３州^高狀態(tài)，且瓦斯?jié)舛茸兓^緩慢，分析可知，采空區(qū)瓦斯主要通過高抽巷進(jìn)行排放。

對比圖8、圖10兩層頂板垮落瓦斯?jié)舛茸兓€，下層堅硬頂板斷裂后，監(jiān)測點最大濃度上隅角為4.8%回風(fēng)巷為3.4%，高抽巷為5.6%；上層堅硬頂板破斷后，監(jiān)測點最大濃度上隅角濃度為9.3%，回風(fēng)巷為6.1%，高抽巷為8.3%；結(jié)合公式（7）可以看出，回風(fēng)側(cè)瓦斯?jié)舛扰c頂板面積、頂板載荷呈現(xiàn)正相關(guān)性，上層寬厚頂板破斷更容易造成瓦斯災(zāi)害，其中，上隅角瓦斯?jié)舛茸兓顬槊黠@。

由于現(xiàn)場已對上隅角及回風(fēng)巷進(jìn)行了防瓦斯超限處理，瓦斯?jié)舛茸兓幻黠@，高抽巷由于位置的原因受人為因素影響較小，能夠很好的反應(yīng)采空區(qū)瓦斯涌出變化，同時支架阻力的變化能夠較好的反映頂板圍巖活動，高抽巷瓦斯?jié)舛扰c支架加權(quán)平均阻力對應(yīng)關(guān)系如圖11所示。根據(jù)現(xiàn)場高抽巷檢測數(shù)據(jù)與模擬計算結(jié)果進(jìn)行對比，具有很好的統(tǒng)一性。

圖11　高抽巷瓦斯?jié)舛扰c支架加權(quán)平均阻力對應(yīng)關(guān)系

由圖11可知，支架加權(quán)平均阻力呈現(xiàn)大小周期變化，大周期為11天～13天（按日推進(jìn)度5.6 m計算，即對應(yīng)步距61.6 m～72.8 m），小周期6天～7天（對應(yīng)步距33.6 m～39.2 m），與計算頂板垮落步距相吻合，同時對應(yīng)高抽巷瓦斯?jié)舛茸兓?，亦呈現(xiàn)大小周期變化，其變化周期與支架阻力變化相一致，很好的反映上下層頂板垮斷下的瓦斯涌出關(guān)系，其中大周期最大濃度可達(dá)7.9%，根據(jù)高抽巷風(fēng)量1 000 m3/min計算，采空區(qū)瓦斯涌出將大于79 m3/min，小周期最大濃度可達(dá)6.1%，采空區(qū)瓦斯涌出將大于61 m3/min，且與模擬計算結(jié)果最大瓦斯?jié)舛认鄬?yīng)，對頂板垮落下的采空區(qū)瓦斯運(yùn)移模型進(jìn)行了驗證。

3　結(jié)論

1、結(jié)合空氣動力學(xué)理論及牛頓第二定律，建立了頂板垮落—瓦斯涌出的耦合模型，分析了瓦斯涌出與頂板載荷、采空區(qū)高度、頂板面積等影響因素的關(guān)系；

2、運(yùn)用多物理場耦合計算揭示了頂板垮落沖擊作用下的采空區(qū)瓦斯動態(tài)運(yùn)移規(guī)律，對比兩種工況下的監(jiān)測點瓦斯?jié)舛茸兓€，頂板垮落時采空區(qū)瓦斯涌入工作面上隅角和高抽巷，其中在上層頂板垮落時，瓦斯?jié)舛让黠@增大，上隅角瓦斯?jié)舛韧辉龇扔葹槊黠@，因此對采空區(qū)及上隅角瓦斯應(yīng)重點治理，結(jié)合理論公式以及塔山頂板狀況，分析可知采空區(qū)體積、頂板面積及載荷與瓦斯?jié)舛茸兓尸F(xiàn)一定的正相關(guān)性；

3、現(xiàn)場實測支架阻力及高抽巷瓦斯?jié)舛茸兓?，瓦斯?jié)舛茸兓c支架阻力相吻合，均呈現(xiàn)大小周期變化，且其變化步距與頂板垮落步距一致，對頂板垮落下的采空區(qū)瓦斯運(yùn)移模型進(jìn)行了驗證，同時，可通過高抽巷瓦斯?jié)舛茸兓从稠敯宓倪\(yùn)動狀況。

參考文獻(xiàn)

［1］孟憲銳，張文超，賀永強(qiáng).高瓦斯綜放面瓦斯涌出特征研究［J］.采礦與安全工程學(xué)報，2006，23（4）：419-422.

［2］張瑞林，李東印，魏軍，等.瓦斯涌出影響因素及其變化特征研究［J］.煤炭科學(xué)技術(shù)，2005，33（10）：20-22.

［3］李樹剛.綜放開采圍巖活動及瓦斯運(yùn)移［M］.徐州：中國礦業(yè)大學(xué)出版社，2000：60-72.

［4］翁明月，徐金海，李沖.綜放工作面煤巖破壞及礦壓顯現(xiàn)與瓦斯涌出關(guān)系的實測研究［J］.煤炭學(xué)報，2011，36（10）：1709-1714.

［5］李化敏，王文，熊祖強(qiáng).采動圍巖活動與工作面瓦斯涌出關(guān)系［J］.采礦與安全工程學(xué)報，2008，25（1）：11-16.

［6］顧鐵鳳，宋選民.封閉采空區(qū)頂板垮落-空氣沖擊耦合模型與差分解法［J］.煤炭學(xué)報，2008，11:1211-1215.

［7］息金波.采場頂板大面積垮落的颶風(fēng)災(zāi)害的理論研究［D］.太原理工大學(xué)，2006.

［8］楊耀乾，平板理論［M］.中國鐵道出版社，1980:211-212.

［9］錢鳴高，廖協(xié)興，許家林，等.巖層控制的關(guān)鍵詞理論［M］.中國礦業(yè)大學(xué)出版社，2003：11.

［10］梁運(yùn)濤，張騰飛，王樹剛，等.采空區(qū)孔隙率非均質(zhì)模型及其流場分布模擬［J］.煤炭學(xué)報，2009，09：1203-1207.

Study on the Influence on Gas Emission in Extreme Thick coal Seam Working Face of the Hard Roof Movement

Lu Yiyu1，2Xia Binwei1，2Guo Shijie1，2
（1.State Key Laboratory for Coal Mine Disaster Dynamics and Control，Chongqing University，Chongqing 400044；2.National and Local Joint Engineering Laboratory of Gas Drainage in Complex Coal Seam，Chongqing University，Chongqing 400044）

Abstract:Aiming at the problem of gas abnormal emission from roof movement，based on Newton secondary law and air dynamics principle，analyzed the relationship between gas emission and roof load，height of goaf，roof area and other factors，the model of gas emission with roof caving was built. In this paper，the 8214 working face of Erpan in Tashan was selected to analyze the gas emission law in the goaf under the overlying strate with the aid of COMSOL multi-physical field coupling software；The result shows that when the layered roof rock caving，the most of mining goaf gas will be poured into the working face corner and the high pumping lane，the gas concentration in the short time increases sharply and then decreases，and the amount of emission is directly related to the roof load，the height of the mined out area and the area of the roof；Hard roof caving upper caused more gas emission than the lower，and the upper corner gas concentration change is obvious.

Key words：extra thick coal seam；hard top slab；top slab caving；the pressure of working face；gushing gas

中圖分類號TD712

文獻(xiàn)標(biāo)識碼A

文章編號1000-4866（2016）03-0001-05

基金項目：長江學(xué)者和創(chuàng)新團(tuán)隊發(fā)展計劃資助（IRT13043）；國家自然科學(xué)基金項目（51374258）；國家自然科學(xué)基金項目（51104191）

作者簡介

盧義玉，男，教授，博士生導(dǎo)師，長江學(xué)者特聘教授，煤礦災(zāi)害動力學(xué)與控制國家重點實驗室主任，重慶大學(xué)資源及環(huán)境科學(xué)學(xué)院副院長。

收稿日期：2016-06-01