基于不同正斷層角度的上盤(pán)工作面突水?dāng)?shù)值模擬分析

劉晉予

(山西煤炭進(jìn)出口集團(tuán) 蒲縣豹子溝煤業(yè)有限公司,山西 臨汾 041204)

礦井水害嚴(yán)重影響著煤礦安全生產(chǎn),給國(guó)家和人民帶來(lái)巨大的損失。對(duì)于含有斷層的工作面在開(kāi)采過(guò)程中更易發(fā)生突水事故,研究表明,在諸多斷層種類(lèi)中正斷層更易導(dǎo)水,而在正斷層中不同的斷層傾角對(duì)工作面突水的影響又有所不同,本文針對(duì)某礦41112工作面位于正斷層上盤(pán)這一情況,模擬不同斷層傾角條件下的工作面推進(jìn)情況。研究不同斷層傾角對(duì)工作面突水的影響。

1 工作面概況

某礦41112工作面,煤層結(jié)構(gòu)復(fù)雜,煤層厚度變化較小,煤層中間有夾石,夾石層小于2層,煤層厚度平均1.8 m,視密度1.33 t/m,煤層埋深約740 m,底板與四灰距離35 m,四灰的厚度約8 m。工作面偏北側(cè)為F21斷層,北高南低,落差8 m～12 m,平均傾角59.6°。

2 數(shù)值模擬

2.1 FLAC3D原理

FLAC3D的原理是拉格朗日差分法,這種方法擁有和有限元法一樣的優(yōu)點(diǎn),可以運(yùn)用它來(lái)分析含有孔隙水壓力的流固耦合的問(wèn)題,它擁有很強(qiáng)的滲流運(yùn)算能力,可以模擬多孔介質(zhì)中的流體運(yùn)動(dòng),可以單獨(dú)運(yùn)算流體模型,只考慮滲流作用,同時(shí)也可以進(jìn)行流固耦合計(jì)算[1-3]。

2.2 模型建立

根據(jù)工作面的相關(guān)水文地質(zhì)條件建立斷層傾角分別為30°、45°、60°,長(zhǎng)、寬、高分別為160 m、120 m、100 m的三種模型，見(jiàn)圖1。力學(xué)邊界條件為下部邊界固定,設(shè)定埋深740 m,載荷均布在模型上部,上覆巖層的設(shè)定密度為2 600 kg/m3,水壓力邊界固定,底部含水層水壓設(shè)為7.5 MPa。含水層巖體飽和度設(shè)為1,工作面底板內(nèi)初始水壓力是依據(jù)梯度水壓的變化進(jìn)行定義,工作面回采后采空區(qū)成為滲流邊界,所以采空區(qū)邊界取固定的水壓力值為0。

工作面推進(jìn)方向與斷層走向平行,在斷層的邊界處留設(shè)20 m的煤柱,模型的邊界留設(shè)20 m以防邊界效應(yīng),工作面寬60 m,一次采全高。從y=20 m處開(kāi)始推進(jìn),分5步開(kāi)挖,每次10 m,最后推進(jìn)到y(tǒng)=80 m處結(jié)束。分別選取y=40 m、y=50 m、y=60 m、y=70 m、y=80 m時(shí)的塑性區(qū)、孔隙壓力的分布圖研究它們?cè)谕七M(jìn)過(guò)程中的變化規(guī)律。

圖1 不同斷層傾角時(shí)的模型Fig.1 Models at different fault dips

2.2.1 斷層傾角為30°的模型

1)塑性區(qū)分布規(guī)律。分析不同推進(jìn)距離下的塑性區(qū)分布圖(見(jiàn)圖2)可知:在剛開(kāi)始推進(jìn)時(shí),斷層帶已經(jīng)出現(xiàn)了塑性破壞,這正印證了正斷層帶容易被破壞的特征;隨著推進(jìn)繼續(xù)進(jìn)行,在煤層頂板的塑性破壞區(qū)越來(lái)越大,留設(shè)的斷層煤巖柱塑性破壞區(qū)域也逐漸增大直至完全被破壞,底板巖層的塑性破壞變化速度不如頂板,但底板破壞區(qū)域和深度也有所增大,靠近斷層附近底板塑性區(qū)域有向著斷層破碎帶繼續(xù)延伸的趨勢(shì),并且從底板塑性區(qū)整體來(lái)看,靠近斷層側(cè)塑性區(qū)范圍大于遠(yuǎn)離斷層的另一側(cè);當(dāng)工作面推進(jìn)結(jié)束時(shí),煤層頂板、底板的塑性最大破壞深度分別達(dá)24 m、12 m。

圖2 不同推進(jìn)距離下塑性區(qū)分布圖Fig.2 Distribution diagram of plastic zones at different advancing distances

2)孔隙壓力分布規(guī)律。分析不同推進(jìn)距離下的孔隙壓力分布圖(見(jiàn)圖3)可知:在工作面剛開(kāi)始推進(jìn)時(shí),處于斷層下盤(pán)的承壓水已經(jīng)開(kāi)始有沿著斷層裂隙向上導(dǎo)升的趨勢(shì),斷層的上下盤(pán)導(dǎo)水帶尚未溝通;當(dāng)工作面推進(jìn)到y(tǒng)=50 m時(shí),位于斷層上、下盤(pán)處的導(dǎo)水帶開(kāi)始導(dǎo)通,繼續(xù)推進(jìn)后下部的承壓水開(kāi)始出現(xiàn)出現(xiàn)較為明顯的導(dǎo)升;當(dāng)工作面推進(jìn)到y(tǒng)=

70 m時(shí),承壓水的導(dǎo)升趨勢(shì)開(kāi)始減緩,此時(shí)承壓水區(qū)域及斷層含水區(qū)域有向工作面下方底板區(qū)域擴(kuò)散擴(kuò)大的趨勢(shì);當(dāng)工作面推進(jìn)到y(tǒng)=80 m時(shí),承壓水繼續(xù)導(dǎo)升且向煤層底板擴(kuò)散范圍擴(kuò)大。與圖2中的塑性區(qū)分布特征進(jìn)行對(duì)比,分析原因,主要是工作面推進(jìn)中,靠近斷層附近底板塑性區(qū)范圍及破壞深度增加,承壓水沿裂隙滲透。

圖3 不同推進(jìn)距離下孔隙壓力分布圖Fig.3 Distribution diagram of pore pressure at different advancing distances

2.2.2 斷層傾角為45°的模型

1)塑性區(qū)分布規(guī)律。分析不同推進(jìn)距離下的塑性區(qū)分布圖(見(jiàn)圖4)可知:工作面推進(jìn)時(shí),斷層帶發(fā)生塑形破壞;隨著工作面的推進(jìn),斷層帶和煤層頂板的塑性區(qū)范圍持續(xù)變大,而底板在工作面推進(jìn)到y(tǒng)=50 m時(shí)塑形區(qū)范圍趨于穩(wěn)定,破壞范圍不在發(fā)生顯著變化,而靠近斷層側(cè)底板的塑性區(qū)范圍比遠(yuǎn)離斷層側(cè)明顯偏大。將煤層工作面推進(jìn)結(jié)束時(shí),煤層頂板、底板的塑性最大破壞深度達(dá)到24 m、10 m;隨推進(jìn)距離增大,煤巖柱塑形破壞加大,且有與斷層接通的趨勢(shì)。

圖4 不同推進(jìn)距離下塑性區(qū)分布圖Fig.4 Distribution diagram of plastic zones at different advancing distances

2)孔隙壓力分布規(guī)律。分析不同推進(jìn)距離下的孔隙壓力分布圖(見(jiàn)圖5)可知:隨著工作面推進(jìn),斷層上、下盤(pán)的含水帶接通,并且接通范圍隨推進(jìn)明顯增大;含水層水沿?cái)鄬訋蛏蠈?dǎo)升,并且導(dǎo)升趨勢(shì)及高度明顯,當(dāng)工作面推進(jìn)到y(tǒng)=60 m時(shí),導(dǎo)升已經(jīng)接近煤層;工作面推進(jìn)結(jié)束后,承壓水到達(dá)煤柱附近且有向煤柱滲透的趨勢(shì),同時(shí)導(dǎo)升繼續(xù)升高的趨勢(shì)。如若此時(shí)煤柱因塑性破壞嚴(yán)重,極有可能使得斷層內(nèi)水進(jìn)入工作面引起突水。

圖5 不同推進(jìn)距離下孔隙壓力分布圖Fig.5 Distribution diagram of pore pressure at different advancing distances

2.2.3 斷層傾角為60°的模型

1)塑性區(qū)分布規(guī)律。分析不同推進(jìn)距離下的塑性區(qū)分布圖(見(jiàn)圖6)可知:從圖中可以看出,當(dāng)工作面推進(jìn)到y(tǒng)=40 m時(shí),斷層帶整體發(fā)生塑形破壞;頂板的塑性區(qū)的范圍隨著繼續(xù)推進(jìn)而增大,底板處的塑性破壞區(qū)域基本不變,但靠近斷層的范圍比遠(yuǎn)離斷層的范圍大,隨著推進(jìn)的進(jìn)行,煤巖柱的塑形破壞程度變大,有和斷層相溝通的趨勢(shì)。工作面推進(jìn)結(jié)束時(shí)時(shí),煤層頂、底板的最大塑性破壞深度達(dá)到24 m、10 m。

圖6 不同推進(jìn)距離下塑性區(qū)分布圖Fig.6 Distribution diagram of plastic zones at different advancing distances

2)孔隙壓力分布規(guī)律。分析不同推進(jìn)距離下的孔隙壓力分布圖(見(jiàn)圖7)可知:隨工作面推進(jìn),斷層上、下盤(pán)的含水帶快速貫通,承壓水沿?cái)鄬訋?dǎo)升速度十分明顯;工作面推進(jìn)到y(tǒng)=50 m時(shí),承壓水到達(dá)工作面附近且向煤層滲透,范圍比較明顯,且在下盤(pán)煤層出現(xiàn)滲透;推進(jìn)結(jié)束時(shí),承壓水繼續(xù)向上導(dǎo)升,且斷層上、下盤(pán)煤層區(qū)域滲透范圍有變大的趨勢(shì)。

圖7 不同推進(jìn)距離下孔隙壓力分布圖Fig.7 Distribution diagram of pore pressure at different advancing distances

綜合不同斷層傾角的模擬結(jié)果分析可知,底板塑性破壞范圍及深度隨斷層傾角的增大有減小趨勢(shì),其中靠近斷層側(cè)塑性區(qū)變化最為明顯,破壞深度由傾角30°時(shí)的12 m減小到60°時(shí)的10 m;斷層傾角為30°時(shí),斷層邊界煤巖柱完全發(fā)生塑形破壞,而45°、60°時(shí),煤巖柱塑形區(qū)變大,卻未完全破壞;隨斷層傾角增大,孔隙壓力從2 MPa增加到4 MPa,承壓水的導(dǎo)升高越明顯,從傾角30°時(shí)的5 m增大到60°時(shí)的12 m,增加率達(dá)140%,但30°時(shí),承壓水在底板滲透范圍較大[4-5]。

3 結(jié)論

1)斷層傾角越小,底板及煤巖柱塑性破壞越嚴(yán)重,底板斷層及含水層越容易與底板塑性區(qū)導(dǎo)通引發(fā)突水。

2)斷層傾角越大,斷層越容易活化成為承壓水導(dǎo)升通路,且導(dǎo)升幅度不斷增加,此時(shí)斷層內(nèi)水容易通過(guò)煤柱塑性區(qū)與工作面導(dǎo)通引起突水。

[1] 卜萬(wàn)奎,茅獻(xiàn)彪.斷層傾角對(duì)斷層活化及底板突水的影響研究[J].巖石力學(xué)與工程學(xué)報(bào),2009,28(2):386-394. BU Wankui,MAO Xianbiao.Research on Effect of Fault Dip on Fault Activation and Water Inursh of Coal Floor[J]. Chinese Journal of Rock Mechanics and Engineering,2009,28(2):386-394.

[2] 苗旺,張培森,王浩,等.斷層傾角對(duì)逆斷層活化規(guī)律影響的模擬分析[J].煤礦安全,2016,47(9):187-189,193. MIAO Wang,ZHANG Peisen,WANG Hao,etal.Simulation Analysis for Effect of Fault Dip on Reverse Fault Activation [J].Safety in Coal Mines,2016,47(9):187-189,193.

[3] 吳佐漢,金愛(ài)兵,吳順川,等.斷層傾角對(duì)煤層底板突水影響的機(jī)理研究[J].煤炭技術(shù),2014,33(7):57-60. WU Zuohan,JIN Aibing,WU Shunchuan,etal.Study of Mechanism About Fault Dip’s Impact on Water Inrush of Coal Floor [J].Coal Technology,2014,33(7):57-60.

[4] 張寧.采動(dòng)條件下斷層組合突水?dāng)?shù)值模擬及防治研究[D].青島:山東科技大學(xué),2015.

[5] 陳文國(guó).工作面走向斷層突水模擬及防治研究[D].青島:山東科技大學(xué),2016.