蔣金泉 李小裕 丁 楠 徐傳偉

(1.山東科技大學(xué)礦山災(zāi)害預(yù)防控制省部共建國家重點(diǎn)實(shí)驗(yàn)室培育基地，山東省泰安市，271000；2.濟(jì)南城建集團(tuán)有限公司，山東省濟(jì)南市，250031)

斷層是煤礦開采中最為常見的地質(zhì)構(gòu)造之一，它的存在破壞了煤巖體的完整性，在其附近會(huì)產(chǎn)生地應(yīng)力的局部不連續(xù)及集中，尤其是在工作面開采推進(jìn)過程中遇到斷層時(shí)，工作面前方的煤柱受到斷層的應(yīng)力阻隔效應(yīng)會(huì)產(chǎn)生遠(yuǎn)遠(yuǎn)超出正常工作面超前支承應(yīng)力大小的應(yīng)力集中。這種應(yīng)力集中在淺埋深煤礦開采中體現(xiàn)得不太明顯，但是在深井煤礦開采中顯現(xiàn)得非常明顯，特別是在工作面接近斷層影響區(qū)時(shí)，開采活動(dòng)會(huì)對(duì)斷層煤柱的穩(wěn)定性造成擾動(dòng)，斷層煤柱存在失穩(wěn)的可能性，一旦斷層煤柱失穩(wěn)會(huì)瞬間釋放大量的能量，形成沖擊地壓，造成明顯的礦壓顯現(xiàn)。這些集聚的能量對(duì)工作面的支架等機(jī)械設(shè)備造成嚴(yán)重的破壞，如人員撤離不及時(shí)，會(huì)造成人員安全事故，產(chǎn)生巨大的經(jīng)濟(jì)損失，形成惡劣的社會(huì)影響。代進(jìn)等采用理論推導(dǎo)與現(xiàn)場實(shí)測(cè)的方法研究了上盤先開采和下盤先開采兩種情況下煤柱應(yīng)力演化的差異，得出斷層上盤開采時(shí)，載荷傳遞至下盤的能力強(qiáng)，下盤先開采時(shí)，載荷傳遞與斷層傾角有關(guān)。蔣金泉等采用三維數(shù)值模擬的方法研究了工作面向正斷層推進(jìn)、上盤工作面沿正斷層布置的采動(dòng)應(yīng)力演化特征。王濤等采用數(shù)值模擬和微震監(jiān)測(cè)結(jié)合的方式研究了正斷層條件下工作面推進(jìn)方式與斷層活動(dòng)的相關(guān)性。姜耀東等采用數(shù)值模擬的方法研究了工作面從斷層上盤和下盤向斷層逐漸回采過程中斷層接觸面的法向應(yīng)力和剪切應(yīng)力的時(shí)空演化規(guī)律，斷層上下盤的運(yùn)動(dòng)規(guī)律。上述研究均側(cè)重于對(duì)斷層這一地質(zhì)情況存在時(shí)工作面開采所表現(xiàn)出來的應(yīng)力應(yīng)變的奇異性，對(duì)于工作面埋深影響下斷層工作面開采時(shí)煤柱應(yīng)力演化規(guī)律的研究較少，因此有必要對(duì)不同埋深條件下工作面推進(jìn)至斷層影響區(qū)時(shí)煤柱的穩(wěn)定性進(jìn)行一些研究。本文通過三維數(shù)值模擬軟件FLAC3D模擬研究不同埋深下正斷層上盤工作面推進(jìn)至距斷層10 m、20 m、30 m、40 m、50 m不同寬度煤柱時(shí)煤柱垂直應(yīng)力演化特征及垂直應(yīng)力集中系數(shù)演化規(guī)律。

1 工程背景及數(shù)值模型

1.1 工程背景

楊柳礦10416工作面位于104采區(qū)東邊界，工作面中存在多條斷層，分別為F104-23、F104-43、F104-44、104DF52，工作面開采時(shí)先推過F104-23斷層。

1.2 數(shù)值模型

以楊柳礦10416工作面地質(zhì)情況為背景實(shí)例，采用FLAC3D數(shù)值模擬軟件建立相似三維數(shù)值模型，模型X方向長500 m，Y方向300 m，Z方向200 m，采用軟件內(nèi)置Interface功能建立接觸面模擬斷層，斷層傾角45°，為摩爾—庫倫力學(xué)模型。采用軟件內(nèi)置的Software應(yīng)變軟化力學(xué)本構(gòu)模型，模型水平方向四周及底面采用位移限定，模型頂面分別施加0 MPa、2.5 MPa、5 MPa、7.5 MPa、10 MPa、12.5 MPa、15 MPa、17.5 MPa的均布荷載，方向與Z坐標(biāo)軸負(fù)方向一致，模型四周施加側(cè)壓系數(shù)0.73的水平應(yīng)力，分別模擬工作面埋深150 m、250 m、350 m、450 m、550 m、650 m、750 m、850 m。建立三維模型如圖1所示，模型各巖體巖性物理力學(xué)參數(shù)見表1。

工作面位于斷層上盤，向斷層開挖推進(jìn)，單步開挖10 m，直至開挖至斷層，在工作面中間煤層中間隔5 m布置測(cè)點(diǎn)，監(jiān)測(cè)工作面前方煤柱垂直應(yīng)力隨工作面回采推進(jìn)的演化特征。

表1 巖體物理力學(xué)參數(shù)

1.3 計(jì)算方案

工作面布置于上盤，寬度100 m，以距離模型左邊界100 m作為上盤開采的起始位置，向斷層開采，直至開采到達(dá)斷層，表示實(shí)際開采中工作面由遠(yuǎn)到近逐漸向斷層靠近，考察不同采深條件下工作面在逐漸接近斷層影響區(qū)過程中工作面前方煤柱垂直應(yīng)力及應(yīng)力集中隨采深逐漸加深的演化規(guī)律。

圖1 巷道與斷層三維數(shù)值模型

2 斷層上盤開采垂直應(yīng)力演化規(guī)律

2.1 斷層煤柱垂直應(yīng)力演化特征

工作面不同推進(jìn)長度下前方不同寬度煤柱垂直應(yīng)力峰值隨采深變化演化曲線如圖2所示。由圖2可知，工作面前方的垂直應(yīng)力集中并不隨采深增加呈線性變化，而是呈現(xiàn)臺(tái)階形的階梯上升然后出現(xiàn)突然的下降，最后呈近似平直的直線小斜率緩慢上升。導(dǎo)致垂直應(yīng)力非線性變化的主要原因是：隨著埋深的增加，巖層集聚的彈性能與埋深的平方是一種正比關(guān)系增長；導(dǎo)致應(yīng)力突然下降是因?yàn)槁裆钜欢ǖ臈l件下，工作面向前推進(jìn)導(dǎo)致工作面前方煤柱的尺寸逐漸減小，煤柱不能承受工作面前方及斷層阻隔效應(yīng)帶來的交叉應(yīng)力疊加的作用，使得煤柱失穩(wěn)發(fā)生破壞，應(yīng)力因而突降。

圖2 不同煤柱垂直應(yīng)力峰值隨采深演化曲線

煤柱尺寸為50 m時(shí)，以此次模擬結(jié)果為例，當(dāng)埋深小于250 m時(shí)，曲線斜率較大，表明工作面前方煤柱超前垂直應(yīng)力增速隨著埋深增加增速很快，呈直線型線性增大；埋深250～750 m時(shí)，曲線斜率小于淺埋深斜率，近似于階梯型緩慢上升，說明隨著埋深的增加，垂直應(yīng)力逐階上升，但是增加速度要小于淺埋深。當(dāng)埋深大于750 m時(shí)，應(yīng)力出現(xiàn)斷崖式陡降，主要是因?yàn)椴缮畲?，垂直?yīng)力增加使得煤柱發(fā)生了塑性破壞，失去了承載能力，說明工作面埋深為750 m時(shí)，工作面開采影響范圍到達(dá)前方50 m，即斷層邊界。

煤柱尺寸為40 m和30 m時(shí)，垂直應(yīng)力演化曲線基本重疊，整體變化趨勢(shì)與煤柱為50 m時(shí)相同，淺埋深條件下為直線型上升，中等埋深條件下為階梯型上升，但是其峰值出現(xiàn)在埋深為650 m時(shí)，表明在當(dāng)前模擬條件下，當(dāng)工作面采深為650 m，煤柱保持承載能力的最窄寬度為30 m。在煤柱寬度為30 m條件下，埋深超過650 m時(shí)，煤柱發(fā)生塑性破壞，承載能力降低為煤體殘余支承應(yīng)力大小，埋深為550～650 m時(shí)，應(yīng)力增速加快，說明在煤柱在即將達(dá)到承載極限時(shí)對(duì)于埋深變化較為敏感。

煤柱尺寸為20 m，垂直應(yīng)力演化曲線整體變化趨勢(shì)與40 m和30 m煤柱相同，但是在中等埋深情況下呈近似直線型上升，表明煤柱對(duì)于工作面埋深的變化較為敏感，這是煤柱留設(shè)較窄、承載能力不足的表現(xiàn)。

煤柱尺寸為10 m，曲線變化趨勢(shì)與上述一致，但是煤柱在工作面埋深超過550 m時(shí)開始發(fā)生塑性失穩(wěn)破壞，承載能力逐漸降低，采深達(dá)到650 m及更大時(shí)降低為煤體殘余支承應(yīng)力。

2.2 斷層煤柱垂直應(yīng)力集中系數(shù)演化規(guī)律

工作面不同推進(jìn)長度下工作面前方煤柱垂直應(yīng)力集中系數(shù)隨采深變化演化曲線如圖3所示。

圖3 不同煤柱垂直應(yīng)力集中系數(shù)隨采深演化曲線

由圖3可知，曲線總體上呈先直線式增加，中期階梯式逐級(jí)下降，最后斷崖式下降并維持在一定數(shù)值附近。

以本模擬為例，埋深小于250 m時(shí)，工作面超前支承應(yīng)力集中系數(shù)線性直線上升，表明工作面前方煤柱垂直應(yīng)力集中增長速度要高于工作面埋深增加而引起的原巖應(yīng)力增加速度，此時(shí)應(yīng)力集中系數(shù)較大，但是因?yàn)槁裆钚?，原巖應(yīng)力較小，應(yīng)力集中數(shù)值相對(duì)偏小，工作面開采引起的煤體應(yīng)力集中影響要高于埋深引起的應(yīng)力增加，也就是說在淺埋深開采時(shí)，工作面發(fā)生沖擊地壓的主要因素并不是埋深。依據(jù)實(shí)際統(tǒng)計(jì)資料，工作面淺埋深時(shí)，引起沖擊地壓多是工作面開采布置的不合理及地質(zhì)構(gòu)造所致，因此淺埋深礦井開采應(yīng)該合理進(jìn)行工作面開采布置設(shè)計(jì)，在遇到特殊地質(zhì)構(gòu)造時(shí)加強(qiáng)針對(duì)性措施。

埋深大于250 m時(shí)，工作面超前支承應(yīng)力集中系數(shù)逐階降低，表明垂直應(yīng)力集中增長速度要低于工作面埋深增加引起的煤體應(yīng)力增加速度，此時(shí)應(yīng)力集中系數(shù)減小，但是埋深大，原巖應(yīng)力較高，應(yīng)力集中數(shù)值比較大，工作面開采引起的煤體應(yīng)力集中影響程度要低于埋深引起的應(yīng)力增加，即中等埋深及深井開采時(shí)，由于工作面上覆巖層厚度大的天然影響，工作面發(fā)生沖擊地壓的可能性要高于淺埋深礦井。

3 結(jié)論

通過數(shù)值計(jì)算模擬不同埋深條件下正斷層上盤工作面開采，得到了工作面與斷層之間留設(shè)不同寬度的煤柱時(shí)煤柱承載的工作面超前支承應(yīng)力及應(yīng)力集中系數(shù)隨埋深增加的演化特征規(guī)律。

(1)模擬研究數(shù)據(jù)表明，煤柱承載的支承應(yīng)力及應(yīng)力集中系數(shù)并不是隨埋深顯現(xiàn)線性增加或減小變化，煤柱支承應(yīng)力在淺埋深開采時(shí)呈直線型線性上升，中等及深埋深呈階梯型逐漸緩慢上升，淺埋深工作面開采時(shí)支承應(yīng)力增加速度快于中等及深埋深。

(2)煤柱支承應(yīng)力集中系數(shù)在淺埋深時(shí)線性快速上升，應(yīng)力集中系數(shù)較中等及深埋深大，但是應(yīng)力集中數(shù)值偏小，工作面開采影響是工作面沖擊地壓發(fā)生的主要影響因素；在中等及深埋時(shí)應(yīng)力集中系數(shù)呈現(xiàn)臺(tái)階型減小，但是應(yīng)力集中數(shù)值大，主要是由于上覆巖層較厚，表明埋深是工作面發(fā)生沖擊地壓的主要影響因素，且工作面易發(fā)生沖擊地壓。