房 蕭 巨 峰 何 琪 付家亭

（1.中國(guó)礦業(yè)大學(xué)礦業(yè)工程學(xué)院，江蘇省徐州市，221116；2.煤炭資源與安全開采國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇省徐州市，221008）

千秋礦綜放面巨厚懸空礫巖層采動(dòng)應(yīng)力分布特征數(shù)值模擬研究＊

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

（1.中國(guó)礦業(yè)大學(xué)礦業(yè)工程學(xué)院，江蘇省徐州市，221116；2.煤炭資源與安全開采國(guó)家重點(diǎn)實(shí)驗(yàn)室，江蘇省徐州市，221008）

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

綜放工作面 巨厚懸空礫巖 采動(dòng)應(yīng)力 應(yīng)力分布 數(shù)值模擬

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

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

1 工作面概況

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

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

2 數(shù)值模型的建立

2.1 模擬參數(shù)確定

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

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

圍巖物理力學(xué)性質(zhì)參照該礦工作面實(shí)際巖體力學(xué)特性。節(jié)理特性考慮采動(dòng)影響，圍巖本構(gòu)關(guān)系采用莫爾－庫(kù)侖（Mohr－Coulomb）模型。

2.2 開挖過程的合理模擬

本數(shù)值模擬分析中的巖石工程是多次開挖完成的，由于巖石受力后，其內(nèi)部的應(yīng)力具有傳遞性，因此，前面的每次開挖都對(duì)后面的開挖產(chǎn)生影響，施工順序不同，開挖步驟不同，都有各自不同的最終力學(xué)效應(yīng)，也有不同的巖石工程穩(wěn)定性狀態(tài)。本次模擬的步驟為：

（1）建立走向模型，并進(jìn)行模型的原巖應(yīng)力平衡計(jì)算；

（2）走向模型進(jìn)行分步開挖，模型應(yīng)力平衡計(jì)算，進(jìn)行數(shù)據(jù)的提取與后處理；

（3）根據(jù)千秋礦21采區(qū)實(shí)際開采順序進(jìn)行開挖，先開挖21101和21121工作面，再開挖21181和21201工作面，最后開挖21141工作面。

2.3 測(cè)線與測(cè)點(diǎn)布置

在進(jìn)行初步的試運(yùn)算后，根據(jù)運(yùn)算情況在21141工作面的頂板中設(shè)置監(jiān)測(cè)線，以研究工作面回采過程中的應(yīng)力和位移變化規(guī)律。另外在巨厚礫巖層底部也布置監(jiān)測(cè)線來研究礫巖層的變形規(guī)律。

工作面開掘引起覆巖破斷與運(yùn)動(dòng)，數(shù)值分析過程的時(shí)間設(shè)置雖不能與實(shí)際開采影響的時(shí)間過程相對(duì)應(yīng)，但數(shù)值分析中不同時(shí)間段的應(yīng)力結(jié)果反映了實(shí)際開采過程中巖層應(yīng)力的演化過程。

3 模擬結(jié)果和分析

3.1 巨厚懸空礫巖層一次采動(dòng)應(yīng)力分布特征

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

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

3.2 巨厚懸空礫巖層二次采動(dòng)應(yīng)力分布特征

對(duì)21181和21201工作面進(jìn)行開挖，其對(duì)21141工作面覆巖運(yùn)動(dòng)變形影響如圖4所示。

圖4 二次采動(dòng)后21141工作面覆巖運(yùn)動(dòng)變形情況

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

圖5 二次采動(dòng)后監(jiān)測(cè)線應(yīng)力分布規(guī)律

3.3 不同開采階段應(yīng)力分布特征的比較

當(dāng)一次采動(dòng)后，采動(dòng)應(yīng)力達(dá)到24MPa，當(dāng)二次采動(dòng)后，隨著礫巖的進(jìn)一步彎曲，其下部受拉，應(yīng)力達(dá)到了34.98MPa，相比原巖應(yīng)力有較大的變化，巖層內(nèi)部由于彎曲產(chǎn)生的彈性能容易誘發(fā)沖擊礦壓。如圖6所示。

圖6 不同開采階段模型監(jiān)測(cè)線應(yīng)力變化曲線

當(dāng)煤層在巨厚懸空礫巖層下采動(dòng)時(shí)，煤層覆巖塑性區(qū)的分布規(guī)律（見圖7）可知，以目前21采區(qū)工作面的開采順序開采，21141工作面采空區(qū)上方的巨厚礫巖層產(chǎn)生了明顯的分層塑性變形，而兩側(cè)巨厚礫巖則發(fā)生整體性斷裂，形成相互擠壓的塊體結(jié)構(gòu)，這對(duì)21161工作面的開采是極為不利的。

圖7 21141工作面采后覆巖塑性區(qū)分布

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

圖8 21141工作面采后不同位置監(jiān)測(cè)線垂直應(yīng)力分布圖

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

4 結(jié)論

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

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

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

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

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

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

[3] 唐軍華，白海波，杜峰.采場(chǎng)支承壓力分區(qū)變異特征研究[J].采礦與安全工程學(xué)報(bào)，2011（2）

[4] 尹士獻(xiàn)，李德海，馬永慶.采動(dòng)影響下硐室群穩(wěn)定性預(yù)測(cè)研究[J].采礦與安全工程學(xué)報(bào)，2009（3）

[5] 楊官濤，李夕兵，王其勝等.地下采場(chǎng)失穩(wěn)的能量突變判斷準(zhǔn)則及其應(yīng)用[J].采礦與安全工程學(xué)報(bào)，2008（3）

[6] 徐學(xué)鋒，竇林名，曹安業(yè)等.覆巖結(jié)構(gòu)對(duì)沖擊礦壓的影響及其微震監(jiān)測(cè)[J].采礦與安全工程學(xué)報(bào)，2011（1）

[7] 劉義新，戴華陽，郭文兵.巨厚松散層下深部寬條帶開采地表移動(dòng)規(guī)律[J].采礦與安全工程學(xué)報(bào)，2009（3）

[8] 范鋼偉，張東升，馬立強(qiáng).神東礦區(qū)淺埋煤層開采覆巖移動(dòng)與裂隙分布特征[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào)，2011（2）

[9] 姚邦華，周海峰，陳龍.重復(fù)采動(dòng)下覆巖裂隙發(fā)育規(guī)律模擬研究[J].采礦與安全工程學(xué)報(bào)，2010（3）

[10] 姜福興，XUN Luo，楊淑華.采場(chǎng)覆巖空間破裂與采動(dòng)應(yīng)力場(chǎng)的微震探測(cè)研究[J].巖土工程學(xué)報(bào)，2003（1）

[11] 康紅普.煤礦井下應(yīng)力場(chǎng)類型及相互作用分析[J].煤炭學(xué)報(bào)，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，F(xiàn)u 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

國(guó)家自然科學(xué)基金面上項(xiàng)目（51074165）；國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目（50834004）

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

（責(zé)任編輯 張毅玲）