許盛運(yùn) 李國(guó)君

(遼寧鐵法能源有限責(zé)任公司, 遼寧 112700)

近年來, 計(jì)算流體力學(xué) (CFD) 技術(shù)已應(yīng)用于煤礦安全的許多領(lǐng)域。澳大利亞、英國(guó)、南非等國(guó)家利用計(jì)算流體力學(xué) (CFD) 數(shù)值模擬, 結(jié)合現(xiàn)場(chǎng)情況, 已經(jīng)開展了大量的工作面瓦斯、粉塵等方面的研究工作, 為煤礦災(zāi)害防治提供了理論依據(jù)。

1 采空區(qū)瓦斯流動(dòng)分布規(guī)律的CFD 模擬研究

采空區(qū)瓦斯流動(dòng)特征是復(fù)雜的, 它受多種因素的影響, 包括通風(fēng)、瓦斯密度、浮力, 以及采空區(qū)滲透等等。針對(duì)鐵法礦區(qū)的具體情況, 本文研究選擇了礦區(qū)兩個(gè)具有代表性的長(zhǎng)壁工作面, 即大興礦N1405、曉明礦N2410 工作面, 建立了三維CFD 模型, 對(duì)工作面采空區(qū)瓦斯流動(dòng)和分布規(guī)律進(jìn)行了模擬, 在此基礎(chǔ)上研究了瓦斯流動(dòng)運(yùn)移對(duì)地面垂直鉆井布井位置和抽采方式的影響。

1．1 大興礦N1405 工作面CFD 模擬研究

大興礦N1405 工作面的有效走向長(zhǎng)度為570m,傾斜長(zhǎng)度為150m。煤層傾角南西向, 最大傾角15°, 最小2°, 平均8．5°。煤層厚度在2．1～2．7m之間, 一般為2．55m, 瓦斯涌出量預(yù)計(jì), 上煤組瓦斯含量為11～27m3/t , 瓦斯壓力為98～401．8 ×104Pa。根據(jù)N1403 工作面瓦斯涌出量實(shí)際測(cè)定,實(shí)際為11．6m3/t。全負(fù)壓U 型通風(fēng), 配風(fēng)量600～900m/min。現(xiàn)場(chǎng)資料顯示, 正?；夭善陂g采空區(qū)瓦斯涌出量約為30～50m3/min (500～1000l/s) , 抽采瓦斯采用井下埋管、頂板鉆孔及地面采空區(qū)井相結(jié)合的方式, N1405 工作面布置了三口地面垂直鉆井。

根據(jù)現(xiàn)場(chǎng)N1405 工作面的基本情況, 建立了三維CFD 模型, 并進(jìn)行了解算。表1 給出工作面CFD模型基本參數(shù)及邊界條件, 圖1 和2 分別給出工作面CFD 模型幾何特征。模型包括了三個(gè)地面垂直鉆井, 其布置方式按實(shí)際位置確定。

1．2 N1405 工作面CFD 數(shù)值模擬結(jié)果

利用以上CFD 數(shù)值模型, 對(duì)N1405 工作面推進(jìn)到距開切眼600m 位置時(shí)采空區(qū)氧氣及瓦斯的分布情況進(jìn)行了模擬研究, 模擬結(jié)果表明：

(1) 采空區(qū)進(jìn)風(fēng)順槽側(cè)氧氣濃度由工作面至開切眼呈遞減趨勢(shì), 尤其是在采空區(qū)瓦斯涌出量較低(300l/s) 的情況下, 在工作面后進(jìn)風(fēng)順槽側(cè)100～200m 范圍內(nèi)基本維持在15%左右的高濃度氧氣;

表1 大興礦N1 采區(qū)405 工作面模型CFD 模型基本參數(shù)及邊界條件

圖1 工作面CFD 模型幾何特征- 平面圖

圖2 工作面CFD 模型幾何特征- 立體圖

(2) 工作面開切眼處氧氣濃度基本維持在10%左右, 表明工作面進(jìn)風(fēng)順槽側(cè)開切眼處存在漏風(fēng)情況;

(3) 回風(fēng)側(cè)氧氣濃度沿開切眼至上隅角氧氣濃度呈遞減趨勢(shì), 直至工作面回風(fēng)流;

(4) 采空區(qū)瓦斯沿回風(fēng)側(cè)匯集, 并沿上隅角方向漸呈上升趨勢(shì), 沿回風(fēng)側(cè)形成瓦斯富集帶, 濃度達(dá)80%以上;

(5) 3 口地面垂直鉆井的孔底位置基本落在瓦斯富集帶內(nèi), 有利于抽采高濃度瓦斯, 整體上, 1號(hào)地面井瓦斯抽采濃度將低于2 號(hào)和3 號(hào)地面井;

(6) 采空區(qū)瓦斯沿開采水平至頂板冒落裂隙帶呈上升趨勢(shì), 沿工作面采空區(qū)深部, 形成瓦斯富集帶;

(7) 隨著工作面的推進(jìn)和采空區(qū)面積的增大,高濃瓦斯富集帶會(huì)隨之?dāng)U大, 在采空區(qū)回風(fēng)側(cè)后部進(jìn)行瓦斯抽采是有利的;

(8) N1405 工作面屬于俯采, 采用地面垂直鉆井對(duì)抽放采空區(qū)瓦斯控制非常有利; 孔位布置非常重要, 除了鉆井間距和回風(fēng)側(cè)距離, 孔底距回采煤層的垂直距離也是決定能否抽采高濃度瓦斯的重要參數(shù), 孔底距回采煤層太近, 將降低抽采瓦斯?jié)舛? 增加采空區(qū)漏風(fēng)強(qiáng)度, 引發(fā)采空區(qū)遺煤自燃;

(9) 在低瓦斯涌出情況下, 工作面采空區(qū)自惰情況差, 采空區(qū)始終處于高濃度氧氣狀態(tài), 尤其是采空區(qū)兩側(cè)和工作面開切眼處, 具備自燃發(fā)火的供氧條件;

(10) 靠近工作面上方的頂板冒落裂隙帶, 瓦斯?jié)舛群凸┭蹩赡芫邆渫咚谷紵龡l件, 應(yīng)采取足夠的安全措施。

圖3 分別給出工作面后50m, 100m, 200m,300m 和400m 位置沿工作面傾向 (橫跨采空區(qū)) 頂板瓦斯?jié)舛鹊哪M變化曲線。結(jié)果顯示, 工作面后100m 內(nèi)采空區(qū)瓦斯?jié)舛茸兓^大, 50m 內(nèi)瓦斯?jié)舛炔怀^10%; 尤其是采空區(qū)進(jìn)風(fēng)側(cè), 該區(qū)域瓦斯運(yùn)移主要受工作面通風(fēng)影響; 工作面后100m 采空區(qū)中部瓦斯?jié)舛茸兓饾u趨于穩(wěn)定并升高, 受風(fēng)壓影響減小, 濃度逐漸達(dá)到60%, 可以看到, 高濃度瓦斯沿回風(fēng)巷側(cè)逐漸形成, 有利于地面井(孔) 的抽采。

圖3 采空區(qū)不同位置沿工作面傾向(橫跨采空區(qū)) 頂板瓦斯?jié)舛茸兓€

圖4 采空區(qū)中間沿工作面走向 (工作面至切眼)頂板不同高度瓦斯?jié)舛茸兓€

圖4 分別給出采空區(qū)中間沿工作面走向 (工作面至切眼) 頂板不同高度 (0m, 10m, 20m 和30m)瓦斯?jié)舛茸兓M曲線。結(jié)果顯示, 在工作面后100m, 頂板30m 范圍內(nèi)瓦斯?jié)舛茸兓€是較大的,瓦斯運(yùn)移主要受工作面通風(fēng)影響, 濃度較低, 不利于頂板孔抽采; 150m 后頂板瓦斯?jié)舛融呌诜€(wěn)定,濃度可達(dá)70%以上, 有利于地面井 (孔) 的抽采。分析模擬結(jié)果顯示, 利用地面井抽采采空區(qū)卸壓瓦斯, 鉆井宜布置在靠近回風(fēng)側(cè)30～80m, 距切眼(第一口井) 50m 范圍內(nèi), 有利抽采高濃度瓦斯。

2 曉明礦N2 采區(qū)410 工作CFD 模擬研究

2．1 曉明礦N2 采區(qū)410 工作面地質(zhì)條件及模型建立

曉明礦N2410 工作面的有效走向長(zhǎng)度為1386m, 長(zhǎng)度為194m。煤層厚度一般為2．54m。全負(fù)壓U 型通風(fēng), 配風(fēng)量1100～1500m3/min?，F(xiàn)場(chǎng)資料顯示, 正常回采期間采空區(qū)瓦斯涌出量約為30m3/min (500l/s) 以上。抽采瓦斯采用井下埋管,頂板及穿層鉆孔及地面采空區(qū)鉆井相結(jié)合的方式,值得注意的是, N2410 工作面回采第一階段為俯采, 而第二階段為仰采。

表2 曉明礦N2 采區(qū)410 工作面模型CFD 模型基本參數(shù)及邊界條件

根據(jù)現(xiàn)場(chǎng)情況, 建立了N2410 工作面三維CFD模型, 模擬工作面推進(jìn)至距切眼500m 位置第一階段俯采以及第二階段時(shí)工作面及采空區(qū)氧氣和瓦斯分布情況。表2 給出工作面CFD 模型基本參數(shù)及邊界條件。該模擬的主要目的是分析對(duì)比工作面開采走向在俯采或仰采時(shí)采空區(qū)氧氣和瓦斯的運(yùn)移規(guī)律。

圖5 給出了N24102 工作面CFD 模型幾例特征及計(jì)算網(wǎng)格方式。本模型采用了約30 萬個(gè)計(jì)算單元。模型包括了1 個(gè)地面垂直轉(zhuǎn)孔, 其布置方式按實(shí)際位置確定。

圖5 N2410 工作面CFD 模型幾何特征及計(jì)算網(wǎng)格劃分

2．2 N2410 工作面CFD 數(shù)值模擬結(jié)果

2．2．1 第一階段： 工作面走向俯采

在工作面走向俯采, N2410 工作面CFD 數(shù)值模擬結(jié)果表明：

(1) 采空區(qū)瓦斯沿回風(fēng)側(cè)匯集, 沿采空區(qū)高位區(qū)域形成瓦斯富集帶, 濃度可達(dá)90%以上;

(2) 地面垂直鉆井的孔底位置落在瓦斯富集帶內(nèi), 優(yōu)于405 地面井布置方式, 更有利于抽采高濃度瓦斯 (可達(dá)90%以上) ;

(3) 410 工作面屬于俯采, 回風(fēng)順槽高于進(jìn)風(fēng)順槽20m 以上; 有利于地面垂直鉆井抽采采空區(qū)瓦斯;

(4) 采空區(qū)瓦斯沿回風(fēng)側(cè)采空區(qū)深部高位點(diǎn)匯集, 對(duì)工作面幾上隅角瓦斯管理更為有利;

(5) 工作面采空區(qū)自惰情況較好, 有利于抑制采空區(qū)自發(fā)火。

2．2．2 第二階段： 工作面走向仰采

在工作面走向仰采時(shí), 模擬結(jié)果表明：

(1) 采空區(qū)瓦斯雖然沿回風(fēng)側(cè)匯集, 沿采空區(qū)深部中間區(qū)域形成較小的高濃度瓦斯富集帶;

(2) 地面垂直鉆井的孔底位置落在瓦斯富集帶邊緣, 不利于地面垂直鉆孔抽采采空區(qū)瓦斯;

(3) 受通風(fēng)、空氣和瓦斯密度的影響, 與走向俯采時(shí)相比, 采空區(qū)瓦斯沿回風(fēng)側(cè)向工作面和上隅角匯集, 極易造成瓦斯超限, 對(duì)瓦斯管理極為不利;

(4) 工作面上方的頂板冒落裂隙帶, 聚積的瓦斯可能具備燃燒條件, 且更接近工作面;

(5) 采空區(qū)一定深度區(qū)域始終處于高濃度氧氣狀態(tài), 自惰情況差, 尤其是采空區(qū)兩側(cè), 具備遺煤自燃的供氧條件。

2．3 地面垂直鉆孔瓦斯抽采的CFD 模擬

利用以上CFD 模型, 研究了地面垂直鉆井布孔位置和抽采方式對(duì)采空區(qū)瓦斯抽采影響。2．3．1 N1采區(qū)405 工作面地面垂直鉆井瓦斯

給出在工作面采空區(qū)瓦斯涌出量為500l/s 時(shí)采用單一鉆孔或多個(gè)鉆孔在不同抽采強(qiáng)度 (混合流量) 下采空區(qū)瓦斯的分布情況, 模擬結(jié)果表明：

(1) 采空區(qū)深部 (靠近開切眼附近) 地面井1瓦斯抽采濃度較低, 在抽采混合流量100l/s (6m3/min) 時(shí), 可以達(dá)到較好的抽采效果;

(2) 抽采流量在200l/s 以上時(shí), 瓦斯抽采濃度會(huì)較大幅度降低, 并增加采空區(qū)富氧帶和空氣涌入深度;

(3) 工作面采空區(qū)中部地面井2 瓦斯抽采效果較好, 在同樣抽采強(qiáng)度 (混合流量200l/s 即12m3/min) 條件下, 瓦斯抽采濃度較高, 可達(dá)到70%以上;

(4) 距工作面一定距離 (150～200m) 的淺部地面井3 瓦斯抽采效果較好, 更利于上隅角瓦斯管理, 但應(yīng)嚴(yán)格控制抽采強(qiáng)度, 不宜超過200l/s 以上, 以防增加采空區(qū)富氧帶和空氣涌入深度;

(5) 與無地面垂直鉆孔抽采瓦斯相比較, 地面鉆孔的利用改變了采空區(qū)局部瓦斯分布, 將上隅角高濃度瓦斯向采空區(qū)深部移動(dòng), 有利于降低上隅角瓦斯?jié)舛? 即便是采空區(qū)深部鉆孔抽采, 也會(huì)影響上隅角瓦斯積聚。

(6) 采用多個(gè)鉆孔中等抽采流量同時(shí)抽采, 尤其是深淺部鉆孔聯(lián)合抽采, 可在采空區(qū)內(nèi)形成較寬區(qū)域的低壓帶, 優(yōu)于采用單一鉆孔高強(qiáng)度集中抽采。

除抽采流量 (混合流量) 變化外, 模擬結(jié)果表明類似以上結(jié)論：

(1) 嚴(yán)格控制抽采流量, 以防增加采空區(qū)富氧帶和空氣涌入深度; 聯(lián)合抽采單孔強(qiáng)度高于12m3/min 時(shí), 采空區(qū)富氧帶和空氣涌入深度會(huì)增加;

(2) 應(yīng)盡量采用多個(gè)鉆孔中等抽采流量同時(shí)抽采, 尤其是深淺部鉆孔聯(lián)合抽采;

N1405 工作面地面垂直鉆井瓦斯現(xiàn)場(chǎng)考察結(jié)果表明, 鉆井混合流量6．11m3/min, 純量4．21m3/min, 平均濃度瓦斯70%, 與以上CFD 模擬結(jié)果較吻合。

2．3．2 N2采區(qū)410 工作面地面垂直鉆井抽采瓦斯

N2采區(qū)410 工作面布置了2 口地面垂直鉆井,第一口鉆井距切眼約110m, 距回風(fēng)順槽70m, CFD模擬分析表明, 由于工作面采用走向俯采, 地面垂直鉆井的孔底位置落在瓦斯富集帶內(nèi), 有利于抽采高濃度瓦斯 (可達(dá)70～90%) ; 第二口鉆井距切眼約800m, 由于煤層傾角的變化, 工作面變?yōu)檠霾?采空區(qū)瓦斯的分布規(guī)律不同于初期俯采; 利用CFD模型, 對(duì)工作面第二階段開采走向仰采時(shí)地面垂直鉆孔進(jìn)行了模擬。模擬結(jié)果表明：

(1) 模型中地面垂直鉆井的孔底位置落在瓦斯富集帶邊緣, 不利于地面垂直鉆井抽放采空區(qū)瓦斯; 抽采瓦斯?jié)舛葘⒌陀诘谝豢阢@井; 第二口鉆井實(shí)際布置位置靠近采空區(qū)中線, 較模型中孔位為宜;

(2) 與無地面垂直鉆井抽采時(shí)相比, 地面垂直鉆井抽放采空區(qū)瓦斯可以將上隅角高瓦斯向采空區(qū)深部移動(dòng), 有利于降低上隅角瓦斯超限, 但影響幅度有限;

(3) 地面垂直鉆井抽采大大增加采空區(qū)富氧帶和空氣涌入深度。

3 結(jié)論

計(jì)算流體力學(xué) (CFD) 模擬已經(jīng)成為一種必要的技術(shù)手段, 用于研究采空區(qū)瓦斯、自然發(fā)火及長(zhǎng)壁工作面粉塵等現(xiàn)象及其控制措施, 其應(yīng)用日益廣泛并在實(shí)踐中得到很好的驗(yàn)證。

采空區(qū)瓦斯涌出是工作面瓦斯涌出的主要來源, 采空區(qū)深部是瓦斯的富集區(qū), 掌握采空區(qū)瓦斯瓦斯運(yùn)移規(guī)律, 特別是高瓦斯工作面, 對(duì)制定采取合理的瓦斯抽排管理措施, 對(duì)提高工作面的生產(chǎn)安全, 提高瓦斯資源利用, 具有特別重要的意義。

針對(duì)鐵法礦區(qū)的具體情況, 本分析研究選擇了本礦區(qū)兩個(gè)具有代表性的長(zhǎng)壁工作面, 即大興礦N1405 工作面、曉明礦N2410 工作面, 建立了三維CFD 模型, 對(duì)工作面采空區(qū)瓦斯流動(dòng)和分布規(guī)律進(jìn)行了的模擬, 在此基礎(chǔ)上, 研究了瓦斯流動(dòng)和分布規(guī)律對(duì)地面垂直鉆井布孔位置和抽采方式的影響。

CFD 模擬揭示了長(zhǎng)壁工作面采空區(qū)瓦斯運(yùn)移多方面的本質(zhì)運(yùn)移規(guī)律, 模擬結(jié)果對(duì)本分析研究的現(xiàn)場(chǎng)實(shí)踐和具體瓦斯抽采方法, 尤其是采空區(qū)地面鉆井抽采技術(shù)相關(guān)參數(shù)優(yōu)化設(shè)計(jì), 具有重要的指導(dǎo)意義。