破頂斷面對(duì)圍巖穩(wěn)定性的影響數(shù)值模擬研究

方 勇，婁培杰，張必昇

(安徽理工大學(xué) 土木建筑學(xué)院，安徽 淮南 232001)

0 前 言

隨著當(dāng)前世界能源結(jié)構(gòu)形勢(shì)的不斷復(fù)雜變化，煤炭資源在我國(guó)的能源經(jīng)濟(jì)仍占據(jù)總能源50%以上，對(duì)我國(guó)的能源發(fā)展起著主導(dǎo)作用。然而國(guó)內(nèi)煤炭資源所處的工程地質(zhì)條件較為復(fù)雜，大傾角煤層和緩傾斜煤層占很大比例，其中更有上部存在力學(xué)性質(zhì)復(fù)雜的復(fù)合型頂板，導(dǎo)致開發(fā)難度大大增加。礦井在開采煤炭的過程中，各種問題日益突出，世界上主要采煤的國(guó)家對(duì)特殊地質(zhì)條件下的巷道支護(hù)技術(shù)研究較少，對(duì)大傾角煤層回采巷道的圍巖控制技術(shù)研究更少，通常采用理論分析、數(shù)值模擬以及相似模型等手段進(jìn)行近水平巷道失穩(wěn)機(jī)理研究[1-5]。雖然國(guó)內(nèi)外學(xué)者對(duì)大傾角煤層巷道失穩(wěn)機(jī)理與控制技術(shù)做了大量研究，但仍存在很多研究問題。

本文以桃園煤礦為研究對(duì)象，桃園煤礦Ⅱ1044工作面位于大傾角煤巖層。在機(jī)巷掘進(jìn)工作中為滿足一次成巷要求，降低施工難度及安全風(fēng)險(xiǎn)，礦上有兩種斷面方案，一種是梯形斷面，另一種是采取破“L”型內(nèi)嵌式頂板斷面。在兩種方案中，保持高低幫及跨度尺寸不變。本文以此為背景，運(yùn)用FLAC3D數(shù)值模擬方法，只考慮在巷道開挖后不支護(hù)情況下，不破頂梯形斷面與破頂后斷面對(duì)圍巖穩(wěn)定性(尤其高幫和頂?shù)装?影響的對(duì)比，分析巷道開挖后的圍巖應(yīng)力、塑性區(qū)和位移。

1 工程概況

桃園煤礦Ⅱ1044工作面10煤層屬二疊系山西組，煤層厚度為2.5～4.3 m，平均3.5 m，局部含有一層夾矸，厚0～0.4 m，平均0.2 m；煤層傾角28°～44°，平均35°。煤層直接頂為灰色泥巖，厚度為0～2.7 m，平均1 m，中厚層狀，平坦?fàn)?參差狀斷口，斷面富含植物化石碎屑；老頂為灰色至淺灰色細(xì)粒細(xì)砂巖，厚度為1.3～16.3 m，平均9.5 m，中厚層狀，成份以石英為主，硅泥質(zhì)膠結(jié)，頂部富聚菱鐵質(zhì)鮞粒，夾有粉砂巖薄層，中下部夾有粉砂質(zhì)包裹體，垂直裂隙，局部巖層破碎，斷面偶見植物化石碎屑；直接底為深灰色泥巖，厚度為0.3～2.1 m，平均1.2 m，中厚層狀，水平層理，層面見有少量白云母碎片及植物化石，性脆；老底為深灰色粉砂巖，厚度為8.3～11.7 m，平均10 m，中厚層狀，水平層理，層面見有少量白云母碎片及植物化石，性脆。埋深747～776 m，工程中按照最大埋深776 m計(jì)算。

Ⅱ1044工作面總體為走向南北、傾向東的單斜構(gòu)造。根據(jù)周邊巷道實(shí)際揭露以及地面三維地震勘探資料解釋，預(yù)計(jì)該工作面機(jī)巷在掘進(jìn)期間可能揭露的斷層共有2條，落差在0～6 m，其中DF12、BFN106為地面三維地震解釋斷層。

2 模型建立

根據(jù)桃園煤礦Ⅱ1044工作面機(jī)巷地質(zhì)采礦條件，確定模型計(jì)算所需的材料參數(shù)，建立FLAC3D數(shù)值分析模型。為簡(jiǎn)便應(yīng)力計(jì)算，模型尺寸建成104.7 m×24 m×132.8 m(長(zhǎng)×寬×高)，巷道尺寸為高幫5 m，低幫2.8 m，跨度4.6 m。模型四周和底部均采用位移限定邊界(即模型4個(gè)側(cè)面水平方向速度為0，底面豎直方向速度為0)，在模型頂面施加豎直向下的應(yīng)力邊界，水平應(yīng)力按靜水壓力設(shè)置。模型采用摩爾-庫倫(Mohr-Coulomb)屈服準(zhǔn)則。

根據(jù)桃園煤礦Ⅱ1044工作面機(jī)巷圍巖類別，最終將模型圍巖分為11層，模型以巷道中心為起點(diǎn)，巷道上方的巖層自下往上依次為泥巖、細(xì)砂巖、泥巖、細(xì)砂巖、粉砂巖，巷道下方的巖層自上往下依次為10煤、泥巖、粉砂巖、泥巖、粉砂巖、細(xì)砂巖。為更直觀地反映巷道附近位移、應(yīng)力和塑性區(qū)的變化，巷道附近設(shè)置為最小網(wǎng)格尺寸0.1 m×0.1 m×0.1 m(長(zhǎng)×寬×高)。所建立的數(shù)值模擬模型及巖層劃分情況如圖1所示。

模擬時(shí)所用的巖層力學(xué)參數(shù)根據(jù)實(shí)驗(yàn)室試驗(yàn)結(jié)果和礦上地質(zhì)資料加以確定。考慮到室內(nèi)試驗(yàn)條件與施工現(xiàn)場(chǎng)存在差異，對(duì)部分巖體力學(xué)參數(shù)做了一定調(diào)整。模擬所取的巖體力學(xué)參數(shù)見表1。

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

3 不同巷道斷面形式圍巖穩(wěn)定性

通過改變巷道的斷面形狀，可以有效控制巷道的變形。在桃園煤礦大傾角煤巖層中，為防破壞巖層層理結(jié)構(gòu)，且保證幫部和跨度尺寸滿足工程現(xiàn)場(chǎng)需求，斷面形狀及支護(hù)方式的選擇至關(guān)重要。就桃園煤礦Ⅱ1044機(jī)巷，探究破“L”型內(nèi)嵌式頂板斷面與不破頂梯形斷面下圍巖穩(wěn)定性的差異。根據(jù)礦上資料及其他大傾角煤礦工程實(shí)例顯示[6-10]，此類巷道主要失穩(wěn)位置在幫部和頂?shù)装?，以下只分析兩種斷面的高幫和頂?shù)装宸€(wěn)定性差異。兩種斷面高低幫、跨度尺寸均相同，破頂點(diǎn)與低幫等高，位于巷道中心線處。

圖2為兩種斷面形式下的垂直應(yīng)力分布情況，可以看出在不支護(hù)情況下開挖，由于上覆巖層荷載的作用，巷道周圍垂直應(yīng)力分布呈現(xiàn)明顯的非對(duì)稱特征。破頂后的巷道垂直應(yīng)力峰值小于未破頂梯形巷道，且出現(xiàn)應(yīng)力集中位置到幫側(cè)的距離不同，應(yīng)力集中范圍要大于梯形巷道。破頂巷道相比梯形巷道頂板面積更大，形成更大應(yīng)力釋放。

(a)梯形斷面

(b)破頂斷面

通過巷道高幫中點(diǎn)，在巷道左部布置50 m長(zhǎng)水平監(jiān)測(cè)點(diǎn)，得到兩種方案高幫指定水平線上的垂直應(yīng)力曲線如圖3所示。

圖3 兩種斷面形式高幫圍巖垂直應(yīng)力對(duì)比

從圖3可以看出破頂斷面在距高幫左側(cè)2.75 m處時(shí)垂直應(yīng)力達(dá)到峰值43.2 MPa，梯形斷面在3 m處達(dá)到峰值41.3 MPa，而后都逐漸趨向圍巖原始應(yīng)力。

為進(jìn)一步對(duì)比兩種斷面形式下的圍巖穩(wěn)定性，對(duì)兩方案圍巖塑性區(qū)進(jìn)行對(duì)比，見圖4。

(a)梯形斷面

(b)破頂斷面

由圖4可見，傾斜煤層梯形巷道塑性區(qū)出現(xiàn)在巷道四周，且沿煤層傾角方向分布，而破頂后巷道頂板上方呈現(xiàn)大范圍塑性區(qū)，巷道底板右下角塑性區(qū)減少，兩種斷面的塑性區(qū)破壞類型均為剪切破壞。

通過布置監(jiān)測(cè)點(diǎn)，可得兩種方案高幫及頂?shù)装逦灰品植记€，如圖5所示。

(a)

(b)

(c)

由圖5可以看出，與塑性區(qū)范圍變化相對(duì)應(yīng)，大傾角煤層破頂后高幫位移量小于未破頂情況下，頂板的位移量卻明顯大于未破頂情況下，梯形斷面高幫最大位移量135.7 mm，破頂斷面高幫最大位移量122 mm，兩種方案中的巷道高幫左側(cè)位移影響范圍在9.5 m左右。

梯形斷面頂板最大位移量47.7 mm，破頂斷面頂板最大位移量80.8 mm。兩種方案的底板位移量差異不大，都在110 mm左右。

4 結(jié) 論

通過對(duì)桃園煤礦Ⅱ1044工作面機(jī)巷破“L”型內(nèi)嵌式頂板斷面與不破頂梯形斷面數(shù)值模擬對(duì)比分析，得到以下結(jié)論。

1)相較于未破頂時(shí)的梯形斷面，破頂斷面低幫側(cè)垂直應(yīng)力峰值更小，但高幫側(cè)應(yīng)力集中峰值要略大，且出現(xiàn)應(yīng)力集中位置更靠近巷道。

2)兩種方案塑性區(qū)對(duì)比，主要差異出現(xiàn)在破頂斷面頂板上方出現(xiàn)大范圍塑性區(qū)。

3)相較于未破頂時(shí)的梯形斷面，破頂斷面高幫位移量減小，但頂板位移量顯著增大，底板位移量差異不明顯。

[ID：013809]