郭 銳

(汾西礦業(yè)集團(tuán)公司 安監(jiān)局,山西 介休 032000)

斷層是斷裂巖層有明顯相對(duì)位移的一種重要地質(zhì)構(gòu)造，在有斷層影響的地質(zhì)條件下，采掘工程破壞了斷層圍巖系統(tǒng)原有的應(yīng)力平衡狀態(tài)，從而使斷層發(fā)生滑移，直至穩(wěn)定，達(dá)到新的平衡狀態(tài)，即斷層活化。在斷層滑移過(guò)程中，圍巖系統(tǒng)釋放大量的能量，引發(fā)震級(jí)較高的礦震，這也是斷層沖擊地壓的特點(diǎn)。以河南躍進(jìn)煤礦為研究對(duì)象，采用UDEC數(shù)值模擬軟件，分析采掘工程對(duì)斷層滑移的影響，為煤礦的安全生產(chǎn)提供理論指導(dǎo)。

1 數(shù)值計(jì)算模型建立

由于斷層的滑移失穩(wěn)，斷層及圍巖系統(tǒng)釋放所積聚的彈性能，其中一部分能量在巖體裂縫動(dòng)力擴(kuò)展中消耗，另一部分則轉(zhuǎn)化為動(dòng)能，使得斷層兩盤(pán)突然卸載而發(fā)生劇烈的震動(dòng)，這是斷層影響沖擊地壓發(fā)生的一個(gè)方面。另一方面，斷層的存在影響了周?chē)鷳?yīng)力場(chǎng)的分布形式，導(dǎo)致工作面頂板局部應(yīng)力集中，而積聚的大量彈性能的釋放，則導(dǎo)致巷道頂板及圍巖的破壞。因此，在研究躍進(jìn)煤礦F16斷層對(duì)沖擊地壓發(fā)生的影響時(shí)，首先對(duì)斷層的滑移失穩(wěn)進(jìn)行研究，而采掘工程對(duì)于F16斷層滑移失穩(wěn)影響則是主要研究?jī)?nèi)容。

采用的數(shù)值模型軟件為Itasca公司的UDEC4.0軟件，是針對(duì)非連續(xù)介質(zhì)模型的二維離散元數(shù)值計(jì)算程序，主要模擬靜載或動(dòng)載條件下非連續(xù)介質(zhì)的力學(xué)行為特征。

1) 建立模型。

根據(jù)躍進(jìn)煤礦實(shí)際地質(zhì)情況建立模型。模型及邊界條件：將模型左右邊界、底部邊界固定，使邊界橫向位移與豎向位移都為“0”，模型上邊界施加17.5 MPa的垂直應(yīng)力，重力加速度設(shè)置為9.8 m/s2。模型范圍：300 m×100 m，斷層傾角45°。煤層頂板至模型上邊界58 m，至模型下邊界42 m，見(jiàn)圖1。

圖1 計(jì)算模型示意圖

2) 巖層屬性。

模型中各煤巖層的物理力學(xué)參數(shù)見(jiàn)表1。

2 開(kāi)采對(duì)斷層滑移數(shù)值模擬分析

2.1 開(kāi)采對(duì)斷層面應(yīng)力的影響分析

1) F16斷層下盤(pán)煤層開(kāi)采對(duì)斷層面應(yīng)力的影響。

表1 煤巖力學(xué)參數(shù)表

根據(jù)所建立模型，可以通過(guò)數(shù)值模擬方法計(jì)算得到，當(dāng)位于F16斷層下盤(pán)附近的工作面向斷層推進(jìn)時(shí)，工作面距離斷層不同位置所對(duì)應(yīng)的斷層面應(yīng)力狀態(tài)，見(jiàn)表2。利用計(jì)算所得數(shù)據(jù)，繪制出在回采工作面的擾動(dòng)下，F(xiàn)16斷層面上剪應(yīng)力與正應(yīng)力比值的時(shí)空關(guān)系曲線，見(jiàn)圖2。斷層結(jié)構(gòu)是否發(fā)生滑落失穩(wěn)，取決于斷層剪應(yīng)力和正應(yīng)力的比值。

表2 斷層面應(yīng)力狀態(tài)統(tǒng)計(jì)表

圖2 下盤(pán)開(kāi)采對(duì)斷層面應(yīng)力狀態(tài)影響的示意圖

從表2及圖2可以看出，當(dāng)工作面與F16斷層之間具有較遠(yuǎn)距離時(shí)，例如80 m、60 m時(shí)，F(xiàn)16斷層面表現(xiàn)出緩慢的應(yīng)力變化，正應(yīng)力的上升速度較為緩慢，應(yīng)力值約為30 MPa，而剪應(yīng)力則在0.1～0.7 MPa變化。當(dāng)工作面分別距離F16斷層40 m、30 m、20 m時(shí)，剪應(yīng)力上升變快，而正應(yīng)力上升依舊緩慢，導(dǎo)致兩者的比值急劇增加。因此，在F16斷層附近區(qū)域進(jìn)行采掘活動(dòng)時(shí)，會(huì)引起斷層面剪應(yīng)力上升速度大于斷層面正應(yīng)力的上升速度，容易導(dǎo)致斷層滑移失穩(wěn)現(xiàn)象的發(fā)生。

2) F16斷層上盤(pán)煤層開(kāi)采對(duì)斷層面應(yīng)力的影響。

利用相同模型，模擬工作面由F16斷層上盤(pán)向斷層推進(jìn)，得到工作面與F16斷層不同距離時(shí)斷層面的應(yīng)力狀態(tài)統(tǒng)計(jì)與斷層剪應(yīng)力與正應(yīng)力比值的時(shí)空關(guān)系變化曲線，見(jiàn)表3，圖3。

表3 斷層面應(yīng)力狀態(tài)統(tǒng)計(jì)表

圖3 上盤(pán)開(kāi)采對(duì)斷層面應(yīng)力狀態(tài)影響的示意圖

由表3及圖3分析可知，在F16斷層附近的工作面自斷層上盤(pán)向斷層推進(jìn)時(shí)，剪應(yīng)力變化趨勢(shì)為先減小后增大，而正應(yīng)力變化趨勢(shì)為逐漸增大。當(dāng)工作面距離斷層大于30 m時(shí)，斷層面剪切應(yīng)力與正應(yīng)力比值緩慢減小，說(shuō)明此時(shí)的斷層面處于一個(gè)穩(wěn)定的狀態(tài)。然而當(dāng)工作面自斷層上盤(pán)繼續(xù)向斷層推進(jìn)時(shí)，剪應(yīng)力上升趨勢(shì)明顯加快，正應(yīng)力變化微弱，斷層剪應(yīng)力與正應(yīng)力的比值隨剪應(yīng)力的變化而變大，此時(shí)，斷層容易發(fā)生滑移失穩(wěn)。從數(shù)據(jù)上可以看出，無(wú)論工作面是在斷層的下盤(pán)或者上盤(pán)，回采工作面向斷層面推進(jìn)都會(huì)對(duì)斷層面應(yīng)力狀態(tài)產(chǎn)生一定的影響。

2.2 開(kāi)采對(duì)斷層滑移的影響分析

1) F16斷層下盤(pán)煤層開(kāi)采對(duì)斷層滑移的影響。

煤層開(kāi)采是斷層沖擊地壓發(fā)生的直接原因。當(dāng)斷層圍巖系統(tǒng)處于穩(wěn)定平衡狀態(tài)時(shí)，遠(yuǎn)場(chǎng)微小擾動(dòng)就會(huì)打破這種平衡狀態(tài)，導(dǎo)致斷層兩盤(pán)發(fā)生相對(duì)移動(dòng)，此時(shí)斷層圍巖系統(tǒng)就會(huì)處于非穩(wěn)定狀態(tài)，繼而引發(fā)沖擊地壓。通過(guò)數(shù)值模擬，可以計(jì)算出由于采掘工程而引起的斷層滑移的量值：在斷層兩盤(pán)的工作面老頂位置布置兩個(gè)相鄰點(diǎn)，即A點(diǎn)與B點(diǎn)，在對(duì)不同位置時(shí)工作面模型進(jìn)行應(yīng)力平衡之前，將A、B兩點(diǎn)X、Y方向上的位移置零后，再對(duì)模型進(jìn)行計(jì)算，應(yīng)力平衡后提取兩點(diǎn)X、Y方向上的位移，再利用公式(1)對(duì)A、B兩點(diǎn)之間距離計(jì)算，計(jì)算出的值L即為斷層的滑移量：

(1)

工作面向斷層面方向推進(jìn)，斷層距離工作面不同位置時(shí)，監(jiān)測(cè)點(diǎn)X、Y方向上的位移及斷層滑移量，見(jiàn)表4，并獲得推進(jìn)后斷層模型X、Y方向上的位移云圖，可直接觀察到斷層面位移變化情況，見(jiàn)圖4。

表4 斷層滑移量統(tǒng)計(jì)表

圖4 下盤(pán)煤層開(kāi)采對(duì)斷層滑移影響的示意圖

根據(jù)斷層的滑移量就可得到開(kāi)采導(dǎo)致的斷層滑移曲線(圖4)。由圖4可以看出，當(dāng)位于F16斷層下盤(pán)的工作面向斷層推進(jìn)，并距離斷層80 m、60 m時(shí)，斷層滑移量變化不大。隨著工作面的推進(jìn)，距離斷層小于60 m時(shí)，斷層滑移量逐步增加，且增加幅度變大。工作面距離斷層20 m時(shí)的斷層滑移量是60 m時(shí)的10倍。斷層滑移量的變化表明，此時(shí)斷層的圍巖系統(tǒng)處于非穩(wěn)定狀態(tài)，很容易發(fā)生斷層沖擊地壓。

2) F16斷層上盤(pán)煤層開(kāi)采對(duì)斷層滑移的影響。

研究方法不變，分析工作面位于F16斷層上盤(pán)時(shí)，斷層滑移量在回采工作面影響下的變化情況。工作面推進(jìn)不同位置時(shí)的滑移量統(tǒng)計(jì)見(jiàn)表5，再根據(jù)計(jì)算出的斷層滑移量得到斷層滑移曲線，見(jiàn)圖5。

表5 斷層滑移量統(tǒng)計(jì)表

工作面向斷層面方向推進(jìn)，斷層距離工作面不同位置時(shí)，監(jiān)測(cè)點(diǎn)X、Y方向上的位移及斷層滑移量見(jiàn)表5，并獲得推進(jìn)后斷層模型X、Y方向上的位移云圖，可直接觀察到斷層面位移變化情況，見(jiàn)圖5。

圖5 上盤(pán)煤層開(kāi)采對(duì)斷層滑移量的影響

由圖5可知，距離斷層面相同，工作面在F16下盤(pán)開(kāi)采較上盤(pán)開(kāi)采對(duì)F16斷層滑移量的影響更大。無(wú)論工作面位于F16斷層上盤(pán)還是下盤(pán)，斷層滑移量均在工作面與斷層距離為20 m時(shí)達(dá)到最大值。當(dāng)工作面位于F16斷層上盤(pán)向斷層推進(jìn)時(shí)，斷層最大滑移量達(dá)到25.73 mm，僅為工作面位于F16斷層下盤(pán)時(shí)最大滑移量的52.5%。下盤(pán)開(kāi)采較上盤(pán)開(kāi)采更容易發(fā)生斷層沖擊地壓。

3 結(jié) 論

1) 當(dāng)躍進(jìn)礦25采區(qū)工作面在下盤(pán)開(kāi)采時(shí)，距離F16斷層60～80 m位置期間，斷層滑移量變化平穩(wěn)，維持在4.5～5 mm，而距離斷層20～60 m期間，斷層滑移量變化較大。當(dāng)工作面在上盤(pán)開(kāi)采時(shí)，距離F16斷層40～80 m位置期間，斷層滑移量變化平穩(wěn)，維持在4.6～4.8 mm，而距離F16斷層20～40 m期間，斷層滑移量變化較大。無(wú)論從下盤(pán)或上盤(pán)向斷層開(kāi)采，在距離F16斷層20 m時(shí)，斷層滑移量達(dá)到最大值。

2) 當(dāng)該礦25采區(qū)工作面在下盤(pán)開(kāi)采時(shí)，距F16斷層60 m位置，斷層面剪切應(yīng)力與主應(yīng)力比值上升幅度開(kāi)始增加，說(shuō)明工作面開(kāi)采到此位置時(shí)，采掘工程對(duì)斷層面的應(yīng)力環(huán)境的影響變大。而工作面在上盤(pán)開(kāi)采時(shí)，距F16斷層30 m位置，剪切應(yīng)力與正應(yīng)力的比值才發(fā)生如此變化，說(shuō)明下盤(pán)開(kāi)采更容易改變斷層面的應(yīng)力狀態(tài)。

3) 當(dāng)該礦25采區(qū)工作面在F16斷層上盤(pán)向斷層開(kāi)采時(shí)，與斷層面相同距離，下盤(pán)開(kāi)采較上盤(pán)開(kāi)采對(duì)斷層滑移量的影響更大。當(dāng)工作面位于F16斷層上盤(pán)開(kāi)采時(shí)，斷層最大滑移量達(dá)到25.73 mm，僅為下盤(pán)開(kāi)采時(shí)最大滑移量的52.5%，下盤(pán)開(kāi)采較上盤(pán)開(kāi)采更容易發(fā)生斷層沖擊地壓。

參 考 文 獻(xiàn)

[1]王經(jīng)明,董書(shū)寧,呂 玲，等.采礦對(duì)斷層的擾動(dòng)及水文地質(zhì)效應(yīng)[J].煤炭學(xué)報(bào),1997,17(4):27-31.

[2]潘一山,王來(lái)貴,章夢(mèng)濤，等.斷層沖擊地壓發(fā)生的理論與試驗(yàn)研究[J].巖石力學(xué)與工程學(xué)報(bào),1998,17(6):642-649.

[3]孟召平,彭蘇萍,黎 洪.正斷層附近煤的物理力學(xué)性質(zhì)變化及其對(duì)礦壓分布的影響[J].煤炭學(xué)報(bào),2001,26(6):561-566.

[4]勾攀峰,胡有光.斷層附近回采巷道頂板巖層運(yùn)動(dòng)特征研究[J].采礦與安全工程學(xué)報(bào),2006,23(3):285-288.