顧永強(qiáng)　田志超　劉業(yè)嬌

(內(nèi)蒙古科技大學(xué)建筑與土木工程學(xué)院)

?

邊角煤柱斷層沖擊地壓的數(shù)值模擬*

顧永強(qiáng)田志超劉業(yè)嬌

(內(nèi)蒙古科技大學(xué)建筑與土木工程學(xué)院)

摘要以石圪節(jié)煤礦3501孤島工作面為研究背景，通過數(shù)值模擬，發(fā)現(xiàn)工作面推進(jìn)不同距離時(shí)，斷層處等效應(yīng)力有很明顯的變化，沒有斷層處等效應(yīng)力基本保持穩(wěn)定，變化很小，造成邊角煤柱在回采時(shí)冒落嚴(yán)重，發(fā)生事故概率大，因此有必要采取預(yù)防措施，防止沖擊地壓發(fā)生。該數(shù)值模擬結(jié)果對于類似地質(zhì)條件的煤礦開采，具有一定的借鑒和指導(dǎo)意義。

關(guān)鍵詞邊角煤柱斷層力學(xué)模型沖擊地壓數(shù)值模擬

隨著煤礦資源開采規(guī)模的日益加大，回收利用率的要求逐漸提高,邊角煤是煤炭資源回收中的重要部分。然而在回收“孤島”型邊角煤柱時(shí)，礦難事故頻發(fā)，孤島煤柱發(fā)生沖擊地壓事故越來越嚴(yán)重，成為礦產(chǎn)資源開采過程中急需解決的難題，其中回收斷層中的預(yù)留煤柱時(shí)沖擊的安全防治顯得尤為迫切[1-3]。本文以石圪節(jié)煤礦3501孤島工作面為例，通過數(shù)值模擬，研究邊角煤柱的安全回收。

1　工程地質(zhì)概況

1.1采區(qū)位置及周邊關(guān)系

石圪節(jié)煤礦3501煤柱位于3501工作面停采線以西，TF4斷層以東，南鄰3309下分層工作面運(yùn)中巷，屬于三角形煤柱。煤柱傾斜寬20～80m，走向長約110m。3501工作面上、下分層分別在2003年7月、2007年7月回采完畢，3309上分層工作面在1996年4月回采完畢，TF4斷層以北的3503工作面上、下分層分別在2005年10月、2007年7月回采完畢。

1.2煤層賦存特征

開采煤層為3層煤，通過地質(zhì)資料分析，該煤柱范圍內(nèi)3層煤賦存穩(wěn)定，煤層厚度3.5～4.1m，平均3.8m。煤層傾角18°～22°，平均20°，以氣煤為主，內(nèi)生裂隙發(fā)育，普氏系數(shù)f=2.5。煤層頂?shù)装迩闆r見表1。

表1　煤層頂?shù)装迕枋?/p>

2　數(shù)值模擬過程分析

2.1數(shù)值計(jì)算模型

建立石圪節(jié)3#煤數(shù)值模型，工作面長度及兩側(cè)巷道和上下區(qū)段長度合計(jì)120m。模型四端固定水平運(yùn)動(dòng)，底端固定垂直運(yùn)動(dòng)，最上面墊層按等效載荷施加到模型頂端，各巖層物理參數(shù)如表2所示。

表2　巖層物理參數(shù)

工作面推進(jìn)時(shí)，距斷層不同距離時(shí)的計(jì)算簡圖見圖1。

圖1　工作面推進(jìn)時(shí)的簡化模型

2.2數(shù)值模擬計(jì)算結(jié)果分析

工作面推進(jìn)30m時(shí)距離斷層80m的模型，沒有材料屬性顯示的部分表示工作面推進(jìn)30m的采空區(qū)，采空區(qū)大小30m×75m，該模型總的大小為323m×78m，并定義各巖層屬性(見圖2)，然后進(jìn)行網(wǎng)格劃分，劃分成大小3m×2m微單元(見圖3)[4-5]。

圖2　距離斷層80 m的簡化模型

圖3　網(wǎng)格劃分后的模型

劃分好單元后，底端固定，邊界施加約束條件，左右兩邊不發(fā)生水平移動(dòng)，在上端施加2 500kN/m均布?jí)毫Α?/p>

2.2.1工作面推進(jìn)50 m

當(dāng)工作面推進(jìn)50m時(shí)，距離斷層60m的計(jì)算等效應(yīng)力云圖見圖4，煤層處等效應(yīng)力曲線見圖5。

圖4　等效應(yīng)力云圖

圖5　等效應(yīng)力曲線

當(dāng)工作面推進(jìn)50m時(shí)，等效應(yīng)力左邊為5.855MPa，右邊為6.915MPa，隨后兩旁等效應(yīng)力逐漸下降，但到斷層處的等效應(yīng)力上升，斷層處的等效應(yīng)力極大值左邊3.300MPa右邊3.200MPa，左邊比右邊大點(diǎn)。在斷層附近的等效應(yīng)力變化幅度相比開采30m時(shí)變化不大(變化幅值為0.9MPa，30m時(shí)為1.005MPa)。

2.2.2工作面推進(jìn)70 m

當(dāng)工作面推進(jìn)70m，距離斷層距離40m，等效應(yīng)力云圖見圖6、煤層處等效應(yīng)力曲線圖見圖7。

圖6　等效應(yīng)力云圖

圖7　等效應(yīng)力曲線

當(dāng)工作面推進(jìn)70m時(shí)，采空區(qū)兩旁的等效應(yīng)力更大，左邊小于右邊(左邊8.67MPa，右邊10.51MPa)。斷層處等效應(yīng)力極大值變大(左邊為4.625MPa、右邊為3.824MPa)，斷層處的等效應(yīng)力變化幅度比推進(jìn)30、50m時(shí)大(幅值為1.213MPa)。斷層結(jié)束后等效應(yīng)力緩慢下降，最后等效應(yīng)力趨向穩(wěn)定2.34MPa。

2.2.3工作面推進(jìn)90 m

當(dāng)工作面推進(jìn)90m，距離斷層距離20m，等效應(yīng)力云圖見圖8、煤層處等效應(yīng)力曲線見圖9。

圖8　等效應(yīng)力云圖

圖9　等效應(yīng)力曲線

當(dāng)工作面推進(jìn)90m時(shí)，采空區(qū)兩旁的等效應(yīng)力更大，采空區(qū)的左邊9.88MPa右邊13.24MPa，斷層處的等效應(yīng)力極大值變大(左邊6.98MPa，右邊5.220MPa)，變化幅度約0.9MPa，斷層后的等效應(yīng)力逐漸趨向穩(wěn)定2.42MPa。

2.2.4工作面推進(jìn)100 m

工作面推進(jìn)100m，距離斷層距離10m，等效應(yīng)力云圖見圖10、煤層處等效應(yīng)力曲線見圖11。

圖10　等效應(yīng)力云圖

當(dāng)工作面推進(jìn)100m時(shí)，采空區(qū)兩旁的等效應(yīng)力左邊13.320MPa、右邊15.0MPa，斷層處的等效應(yīng)力極大值變大，左邊12.22MPa、右邊6.75MPa，斷層處的等效應(yīng)力變化幅度約0.8MPa，斷層后的等效應(yīng)力趨向穩(wěn)定值2.450MPa。

圖11　等效應(yīng)力曲線

2.3比較分析

(1)工作面推進(jìn)不同的距離，即距離斷層不同的距離來看，4種情況下斷層處所的等效應(yīng)力是逐漸增大的，開始上升→下降→上升→下降→趨向穩(wěn)定，同時(shí)斷層處的變化幅度不大，基本保持在0.8～1.2MPa，工作面推進(jìn)的距離符合所研究的課題。

(2)工作面推進(jìn)不同距離時(shí)，等效應(yīng)力在斷層與無斷層時(shí)表現(xiàn)不同，斷層處等效應(yīng)力有很明顯的變化，等效應(yīng)力逐漸增大，左邊從3.32MPa增加到12.22MPa，右邊從3.16MPa增加到6.75MPa；沒有斷層處變化很小，等效應(yīng)力基本保持穩(wěn)定在2.42～2.45MPa。

3　結(jié)　論

回收邊角煤柱是提高煤炭資源利用率的重要途經(jīng)，通過數(shù)值模擬計(jì)算，發(fā)現(xiàn)工作面推進(jìn)不同距離時(shí)，等效應(yīng)力在有斷層與無斷層時(shí)表現(xiàn)不同，斷層處等效應(yīng)力逐漸增大，沒有斷層處變化很小，等效應(yīng)力基本保持穩(wěn)定。因此斷層在回采時(shí)很容易發(fā)生巖層冒落，而且隨著工作面推進(jìn)冒落地越來越嚴(yán)重，同時(shí)沖擊地壓事故在斷層處發(fā)生的概率相比其他大的多，有必要采取預(yù)防措施防止沖擊地壓發(fā)生。該數(shù)值模擬計(jì)算結(jié)果對于類似地質(zhì)條件的煤礦開采具有一定的借鑒和指導(dǎo)意義。

參考文獻(xiàn)

[1]田志超,劉業(yè)嬌,王榮超.非線性動(dòng)力學(xué)分岔理論在礦山開采應(yīng)用中的探討[J].現(xiàn)代礦業(yè),2012(11):54-56.

[2]吳教錕,趙旭生,李秋林.大灣煤礦保護(hù)層開采數(shù)值模擬研究[J],礦業(yè)安全與環(huán)保,2007,34(1):6-7.

[3]李樹清,龍祖根,羅衛(wèi)東.煤層群下保護(hù)層開采保護(hù)范圍的數(shù)值模擬[J].中國安全科學(xué)學(xué)報(bào),2012,22(6):34-39.

[4]田志超,劉業(yè)嬌,王傳朋.利用煤層中缺陷體防治沖擊地壓[J].煤炭工程，2010(1):72-74.

[5]閆鑫,田志超,劉業(yè)嬌.錨桿在深埋厚煤層巷道中的支護(hù)效果分析[J].現(xiàn)代礦業(yè)，2013(1):131-132.

(收稿日期2016-04-19)

NumericalSimulationofRockBurstinthePresenceofFaultsintheMarginalCoalPillar

GuYongqiangTianZhichaoLiuYejiao

(CollegeofArchitectureCivilEngineering,InnerMongoliaUniversityofScience&Technology)

AbstractTaking the 3501 isolated working face of Shiyijie coal mine as the engineering background,the numerical simulation results show that the working face advancing different distance,the fault equivalent stress is changed obviously,the equivalent stress is kept stable in the area without faults,which resulting in caving phenomenon of the corner of marginal coal pillar is more and more serious in the process of mining,the probability of the rock burst accident occurred on the faults is enlarged,therefore,it is necessary to take preventive measures to prevent the occurrence of rock burst.The numerical simulation results have some reference for the the mining of the coal mines under the similar geological conditions.

KeywordsMarginal coal pillar,Fault,Mechanical model,Rock burst,Numerical simulation

*內(nèi)蒙古自治區(qū)高等學(xué)?？茖W(xué)研究項(xiàng)目(編號(hào)：NJZY148)；內(nèi)蒙古科技大學(xué)創(chuàng)新基金資助項(xiàng)目(編號(hào):2001NCL033)。

顧永強(qiáng)(1970—)，男，副教授，014010 內(nèi)蒙古包頭市阿爾丁大街7號(hào)。