千秋礦綜放面巨厚懸空礫巖層采動應力分布特征數值模擬研究＊

房 蕭 巨 峰 何 琪 付家亭

（1.中國礦業大學礦業工程學院，江蘇省徐州市，221116；2.煤炭資源與安全開采國家重點實驗室，江蘇省徐州市，221008）

千秋礦綜放面巨厚懸空礫巖層采動應力分布特征數值模擬研究＊

房 蕭1，2巨 峰1，2何 琪1，2付家亭1，2

（1.中國礦業大學礦業工程學院，江蘇省徐州市，221116；2.煤炭資源與安全開采國家重點實驗室，江蘇省徐州市，221008）

利用離散元軟件UDEC4.0對千秋礦21141綜放工作面在巨厚懸空礫巖層下的采動應力分布進行模擬研究。結果表明，隨推采面積的加大，礫巖在不斷彎曲變形，撓度增大，礫巖體對工作面前方煤體及切眼后方煤體的作用力逐漸增大，并隨煤體的開采而重新分布，因此工作面煤體內部支承應力將隨著煤體塑性范圍的擴大向深部轉移，并重新積聚產生新的沖擊危險。

綜放工作面 巨厚懸空礫巖 采動應力 應力分布 數值模擬

煤層頂板巖層特性是沖擊礦壓發生的重要影響因素。研究表明：上覆巖層的懸頂長度越大越容易發生沖擊，且隨著懸頂長度的增加，沖擊危險性迅速增高；堅硬、厚而完整的頂板突然斷裂，易誘發頂板型沖擊礦壓，且基本頂厚度越大，完整性越好，發生的沖擊強度也越大。但是，目前國內外對于巨厚懸空礫巖層采動應力分布研究并不明確，故有必要分別對巨厚懸空礫巖層在兩次采動影響下應力分布規律進行研究。

由于巖土材料結構的不連續性，離散元數值模擬軟件UDEC4.0已廣泛應用于工程巖體模擬研究。煤礦開采將使采場頂板巖體系統產生復雜的節理裂隙和變形塊體，并在頂板巖層變形、斷裂和下沉過程中產生應力擾動和釋放震動能量，本文以義馬煤業公司千秋煤礦21141綜采放頂煤工作面為例，采用UDEC數值模擬軟件對巨厚懸空礫巖影響下采場應力分布規律用離散單元法進行模擬研究。

1 工作面概況

21141綜采放頂煤工作面位于礦井西部二水平21采區下山西翼，北鄰21121采空區，南鄰未開采的21161工作面，西鄰千秋礦、耿村礦邊界煤柱，東鄰21采區下山煤柱。地面標高590.1～617.5m，工作面標高為－38.1～－84.7m，平均采深665.2m。對應地表位置為丘陵山地，地形復雜，沖溝發育，地表無村莊、無水體。工作面可采走向長度1298m，傾斜長度130m，傾角10～14°。該工作面回采二煤，平均煤厚約20m。煤層為黑色塊狀及粉末狀，具瀝青光澤，干燥、疏松破碎，極易自燃。煤層結構復雜，含2～4層夾矸，含矸巖性分別為細砂巖、粉砂巖、泥巖，局部煤體紊亂。一層夾矸為含礫細砂巖，礫石為細粒狀石英巖屑，平均厚度0.5m，基本穩定，特征明顯。二層夾矸為泥巖，較穩定。煤巖類型為半亮型—半暗型。

直接頂為深灰色泥巖，夾菱鐵質薄層，致密均一、斷口平坦，具隱蔽水平層理發育，厚度23.02～27.63m，分布較穩定，老頂為中侏羅系雜色礫巖，厚度較大，平均407m。底板自上而下為煤矸互層、碳質泥巖、含礫泥巖、細砂巖、粉砂巖礫巖，分布不穩定。

2 數值模型的建立

2.1 模擬參數確定

研究上覆巨厚礫巖層對采場回采時的影響，模型不可簡化，模型上邊界為地表，模型的左右邊界不小于一個完整工作面的開采及影響范圍。建立走向模型范圍為2000m×700m。圖1為21141工作面模型圖。

圖1 千秋礦21141工作面模型圖

圍巖物理力學性質參照該礦工作面實際巖體力學特性。節理特性考慮采動影響，圍巖本構關系采用莫爾－庫侖（Mohr－Coulomb）模型。

2.2 開挖過程的合理模擬

本數值模擬分析中的巖石工程是多次開挖完成的，由于巖石受力后，其內部的應力具有傳遞性，因此，前面的每次開挖都對后面的開挖產生影響，施工順序不同，開挖步驟不同，都有各自不同的最終力學效應，也有不同的巖石工程穩定性狀態。本次模擬的步驟為：

（1）建立走向模型，并進行模型的原巖應力平衡計算；

（2）走向模型進行分步開挖，模型應力平衡計算，進行數據的提取與后處理；

（3）根據千秋礦21采區實際開采順序進行開挖，先開挖21101和21121工作面，再開挖21181和21201工作面，最后開挖21141工作面。

2.3 測線與測點布置

在進行初步的試運算后，根據運算情況在21141工作面的頂板中設置監測線，以研究工作面回采過程中的應力和位移變化規律。另外在巨厚礫巖層底部也布置監測線來研究礫巖層的變形規律。

工作面開掘引起覆巖破斷與運動，數值分析過程的時間設置雖不能與實際開采影響的時間過程相對應，但數值分析中不同時間段的應力結果反映了實際開采過程中巖層應力的演化過程。

3 模擬結果和分析

3.1 巨厚懸空礫巖層一次采動應力分布特征

根據千秋煤礦21采區實際情況，首先對21101、21121工作面進行開挖。所以21141工作面在開采之前就深受巨厚礫巖的影響。將21101和21121工作面采后對21141工作面的影響稱為一次采動影響，將21181和21201工作面采后對21141工作面影響稱為二次采動影響。一次采動后，21141工作面覆巖運動變形影如圖2所示。

21101和21121工作面采完后應力分布如圖3所示。從左側煤壁到采空區，21101和21121工作面頂板中應力分布情況可劃分為4個區域：原巖應力區、支承應力影響區、應力降低區和應力恢復區。一次采動后巨厚礫巖層盡管由于其整體結構懸而不垮，但隨著下位巖層的垮落變形，巨厚礫巖層逐漸分層形成了塑性區。

3.2 巨厚懸空礫巖層二次采動應力分布特征

對21181和21201工作面進行開挖，其對21141工作面覆巖運動變形影響如圖4所示。

圖4 二次采動后21141工作面覆巖運動變形情況

巨厚礫巖層的懸而不垮使得已采工作面采空區覆巖不能充分壓實形成承載結構，只有采空區中部區域的覆巖觸矸穩定，工作面上下端頭區域均不承載。巖層承載結構因此發生了根本性的變化：二煤覆巖的自重不再是均載分布，而是由采區上端實體煤、采區下端實體煤和未開采的21141工作面承載。21141工作面二次采動后應力分布如圖5所示。

圖5 二次采動后監測線應力分布規律

3.3 不同開采階段應力分布特征的比較

當一次采動后，采動應力達到24MPa，當二次采動后，隨著礫巖的進一步彎曲，其下部受拉，應力達到了34.98MPa，相比原巖應力有較大的變化，巖層內部由于彎曲產生的彈性能容易誘發沖擊礦壓。如圖6所示。

圖6 不同開采階段模型監測線應力變化曲線

當煤層在巨厚懸空礫巖層下采動時，煤層覆巖塑性區的分布規律（見圖7）可知，以目前21采區工作面的開采順序開采，21141工作面采空區上方的巨厚礫巖層產生了明顯的分層塑性變形，而兩側巨厚礫巖則發生整體性斷裂，形成相互擠壓的塊體結構，這對21161工作面的開采是極為不利的。

圖7 21141工作面采后覆巖塑性區分布

工作面采后，覆巖逐層垮落變形，自下而上形成垮落帶、裂隙帶和彎曲下沉帶。煤層頂板中應力的分布直接受巨厚礫巖層的影響，煤層中最大支承應力與巨厚礫巖層中最大支承應力位于同一水平投影位置，如圖8所示。

圖8 21141工作面采后不同位置監測線垂直應力分布圖

綜上所述，21181和21201工作面采后巨厚礫巖層在采區上部出現了線狀的拉破斷區，使得巨厚礫巖層整體破斷。這種巨大的巖塊具有整體下滑、旋轉等失穩趨勢。而21101和21121采空區上方的礫巖層由于下位巖層的切頂作用，塑性區的分布呈現階段性的向上發展，具有分層而不是整體性的特點。巨厚礫巖層的這種局部塑性變形擴容、整體破斷垮落失穩的特點使得未開采的21141工作面的支承應力達到36.59MPa，為工作面原巖應力的2.01倍。

4 結論

上覆巨厚礫巖層的自重應力及下位巖層的沖擊載荷是21141工作面發生沖擊地壓的主要力源。

（1）隨工作面開采面積加大，礫巖不斷彎曲變形，撓度增大，礫巖體對工作面前方煤體及切眼后方煤體的作用力逐漸增大，因此工作面煤體內部支承應力將隨著煤體塑性范圍的擴大向深部轉移，并重新積聚產生新的沖擊危險。

（2）工作面直接頂初次垮落步距為40m左右，工作面推進60m后，老頂初次來壓。隨著工作面的推進，工作面周期來壓步距為50m左右。當工作面推進300m后，巨厚礫巖懸空距離達到最大，為14.06m。繼續向前推進10m后，巨厚礫巖底部開始分層垮落，說明310m是巨厚礫巖的極限懸空步距。此后，巨厚礫巖不再懸空開始觸底穩定并壓實下位巖層，采空區應力快速增大至原巖應力的1.6倍。

（3）由于巨厚礫巖層整體性好，分層厚度大，利于傳遞上覆巖層對下位煤層的沖擊載荷，對礫巖層形成的集中應力傳播速度較快，從而為沖擊地壓的產生提供了必要的條件。當支承應力傳遞到采場并向采場周圍空間釋放時，就會發生沖擊礦壓事故。

[1] 譚輔清，昝東峰，周楠，鄧雪杰.巨厚礫巖層下工作面過斷層覆巖運動規律研究及應用[J].中國煤炭，2011（9）

[2] 牟宗龍，竇林名，李慧民等.頂板巖層特性對煤體沖擊影響的數值模擬[J].采礦與安全工程學報，2009（1）

[3] 唐軍華，白海波，杜峰.采場支承壓力分區變異特征研究[J].采礦與安全工程學報，2011（2）

[4] 尹士獻，李德海，馬永慶.采動影響下硐室群穩定性預測研究[J].采礦與安全工程學報，2009（3）

[5] 楊官濤，李夕兵，王其勝等.地下采場失穩的能量突變判斷準則及其應用[J].采礦與安全工程學報，2008（3）

[6] 徐學鋒，竇林名，曹安業等.覆巖結構對沖擊礦壓的影響及其微震監測[J].采礦與安全工程學報，2011（1）

[7] 劉義新，戴華陽，郭文兵.巨厚松散層下深部寬條帶開采地表移動規律[J].采礦與安全工程學報，2009（3）

[8] 范鋼偉，張東升，馬立強.神東礦區淺埋煤層開采覆巖移動與裂隙分布特征[J].中國礦業大學學報，2011（2）

[9] 姚邦華，周海峰，陳龍.重復采動下覆巖裂隙發育規律模擬研究[J].采礦與安全工程學報，2010（3）

[10] 姜福興，XUN Luo，楊淑華.采場覆巖空間破裂與采動應力場的微震探測研究[J].巖土工程學報，2003（1）

[11] 康紅普.煤礦井下應力場類型及相互作用分析[J].煤炭學報，2008（12）

Numerical simulation about the stress distribution characteristics of the caving face under the thick impending mining conglomerate layer

Fang Xiao1，2，Ju Feng1，2，He Qi1，2，Fu Jiating1，2
（1.School of Mines，China University of Mining ＆Technology，Xuzhou，Jiangsu 221116，China；2.State Key Laboratory of Coal Resources ＆Mine Safety，Xuzhou，Jiangsu 221008，China）

The stress distribution characteristics of the 21141caving facing in Qianqiu coal mine under the thick impending mining conglomerate layer by discrete element software UDEC4.0have been studied.Result shows that as the increasing of the mining area，the conglomerate constant bending and deflection，the force of cutting the coal ahead and behind of the caving face are increasing and are redistributed.So the internal bearing stress of coal mining face will expand to deep layer as the plastic deformation of coal increasing，which formed new impact dangerous.

the caving face，the thick impending conglomerate，mining stress，stress distribution，numerical simulation

TD323

A

國家自然科學基金面上項目（51074165）；國家自然科學基金重點項目（50834004）

房蕭（1987－），男，遼寧凌源人，在讀碩士研究生，從事采動巖體控制及固體物充填開采等方面的研究。

（責任編輯 張毅玲）