饒運(yùn)章,吳國興,羅 歡

(江西理工大學(xué)資源與環(huán)境工程學(xué)院,江西贛州 341000)

礦山巷道地壓形成機(jī)理及分布規(guī)律的研究

饒運(yùn)章,吳國興,羅 歡

(江西理工大學(xué)資源與環(huán)境工程學(xué)院,江西贛州 341000)

通過對國內(nèi)外學(xué)者研究巷道地壓問題論述的分析,總結(jié)出礦山巷道地壓形成的一般機(jī)理,并運(yùn)用數(shù)值分析軟件和有限差分的基本原理,對半圓拱形巷道在一定的側(cè)壓系數(shù)條件下開挖前后圍巖應(yīng)力狀態(tài)進(jìn)行數(shù)值模擬與分析,總結(jié)出半圓拱形巷道開挖后圍巖應(yīng)力的分布規(guī)律。

巷道地壓;數(shù)值模擬;有限差分;側(cè)壓系數(shù);圍巖應(yīng)力

1 前言

為了開采礦石,常常需要掘進(jìn)各種斷面形狀的巷道,例如:礦山開拓中的階段運(yùn)輸巷道,采準(zhǔn)切割中的的切割巷道,礦石開采中的天井和鑿巖巷道。巷道一旦開挖,原巖應(yīng)力將受到開挖繞動,導(dǎo)致應(yīng)力的重新分布,形成巷道地壓。

這種圍巖應(yīng)力與原巖應(yīng)力相比,表現(xiàn)出很大的應(yīng)力集中現(xiàn)象,例如:巷道周邊環(huán)向應(yīng)力將大幅升高,徑向應(yīng)力明顯降低。對于不支護(hù)巷道,巷道周邊徑向應(yīng)力為0。巷道圍巖中這種重新分布的應(yīng)力組合是否超過巖體的強(qiáng)度極限或某種強(qiáng)度條件,將直接關(guān)系到礦山地下工程的穩(wěn)定性。例如經(jīng)常出現(xiàn)的“地壓顯現(xiàn)”可以表現(xiàn)為巷道底鼓、片幫、冒頂、斷面收縮、支架破壞,甚至垮塌。所以對礦山巷道開挖后地壓的研究顯得相當(dāng)重要。

2 巷道地壓的形成機(jī)理

從根本上講,地應(yīng)力是礦山巷道地壓顯現(xiàn)的根本來源。地應(yīng)力是存在于地層中的天然應(yīng)力,也稱原巖應(yīng)力。在沒有開挖工程擾動的情況下,巖體處于原始平衡狀態(tài)。礦山巷道的開挖,打破了原始平衡狀態(tài),導(dǎo)致地應(yīng)力的釋放,從而引起巖體的變形和向自由面的位移,引起圍巖應(yīng)力的重新分布。圍巖的過量位移和應(yīng)力集中將導(dǎo)致圍巖局部的或整體的失穩(wěn)和破壞。這是地壓發(fā)生的過程和機(jī)理。

因此從本質(zhì)上來定義,地壓就是巖體因受開挖擾動而產(chǎn)生的力學(xué)效應(yīng)。它與巖體的受力狀態(tài)、巖體結(jié)構(gòu)和質(zhì)量、巖體物理力學(xué)性質(zhì)、工程地質(zhì)條件以及時間等因素有關(guān)。

按表現(xiàn)形式,地壓可分為四類:散體地壓、變形地壓、沖擊地壓和膨脹地壓。下面以圓形巷道為例說明開挖引起的巷道圍巖地壓分布規(guī)律。

地壓巖體和金屬材料不同,它不是線彈性體,而是復(fù)雜的非線性體,多數(shù)可簡化為彈塑性體。若以彈塑體來考慮,并假設(shè)原巖應(yīng)力為均勻應(yīng)力場,其值等于P0,那么巷道開挖后一般將在圍巖中形成三個區(qū),即塑性區(qū)(A/B)、彈性區(qū)(C)和原巖應(yīng)力區(qū)(D),其中塑性區(qū)通常又分為兩個區(qū),即松動區(qū)(A)和塑性應(yīng)力升高區(qū)(B),如圖1所示。

圖1 巷道周邊圍巖

地下巷道開挖前,巖體處于原始平衡狀態(tài)。巷道開挖后出現(xiàn)圍巖應(yīng)力重新分布,圍巖的一部分彈性變形產(chǎn)生的壓力為彈性變形地壓。當(dāng)二次應(yīng)力場使彈性圍巖局部地區(qū),也是在巷道周邊地區(qū)的應(yīng)力狀態(tài)達(dá)到屈服條件時,這部分圍巖就是進(jìn)入塑性狀態(tài),產(chǎn)生塑性區(qū)。應(yīng)力再次重新分布,塑性區(qū)的應(yīng)力一部分向彈性區(qū)轉(zhuǎn)移。另一部分因塑性變形而釋放,同時產(chǎn)生塑性變形地壓,塑性變形的進(jìn)一步發(fā)展就會引起巖體的破壞,同時巷道周圍區(qū)域可能出現(xiàn)的拉應(yīng)力也會導(dǎo)致巖體的破裂,從而在最靠近巷道周邊的區(qū)域形成松動區(qū)。上述就是A,B,C三個區(qū)的形成過程。各區(qū)的應(yīng)力、圍移狀態(tài)分析如下。

塑性區(qū):

假設(shè)巖體滿足摩爾-庫侖屈服準(zhǔn)則,根據(jù)屈服條件和平衡方程,有(1),(2)式即可求出塑性區(qū)內(nèi)的應(yīng)力大小。

由公式(1)～(3)可以看出,塑性區(qū)的應(yīng)力狀態(tài)與原巖應(yīng)力P0無關(guān),而只與巖體的強(qiáng)度指標(biāo)(c,φ值)有關(guān)。塑性區(qū)通常也是極限平衡區(qū)。

彈性區(qū):

彈性區(qū)內(nèi)的應(yīng)力大小可由彈性力學(xué)的后壁筒公式,(4)、(5)式求得。

由(4)式、(5)式可以看出來,彈性區(qū)內(nèi)的應(yīng)力大小不僅與巖體強(qiáng)度指標(biāo)(c,φ值)有關(guān),而且與原巖應(yīng)力有關(guān)。

塑性區(qū)的半徑可根據(jù)彈、塑性交界面上的應(yīng)力、位移連續(xù)的條件求得:

根據(jù)廣義胡克定律和伊留申彈塑性小變形理論,求得巷道圍巖中的徑向位移u。因?yàn)樵搯栴}滿足軸對稱條件,故周向位移為零。

3 半圓拱形巷道數(shù)值模擬試驗(yàn)描述

3.1 試驗(yàn)條件

對于圓形和橢圓形等斷面形狀的巷道,在假定巖石(巖體)均質(zhì)、連續(xù)、各向同性的前提下,利用彈性力學(xué)理論,可以求出其圍巖應(yīng)力分布的精確分析解。但是無論是礦山巷道,還是公路和鐵路隧道,其斷面形狀絕大多數(shù)都是直墻拱形,而這種直墻拱形斷面巷道圍巖應(yīng)力分布,到不前為止,還沒有辦法利用彈性或彈塑性理論求得解析解。而隨著計算機(jī)技術(shù)的發(fā)展,數(shù)值模擬已經(jīng)成為一種近似求解的工具,本試驗(yàn)是用數(shù)值模擬軟件flac3D來對半圓形巷道進(jìn)行模擬分析。FLAC3D是一個三維有限差分程序,采用了一個“顯式解”方案,沒有必要存儲剛度矩陣,計算更快捷。對模擬塑性破壞和塑性流動采用的是“混合離散法”。這種方法比有限元法中通常采用的“離散集成法”更為準(zhǔn)確、合理。

3.2 試驗(yàn)?zāi)Ｐ?/p>

假設(shè)研究的巷道所處深度相對于開挖高度大很多,在巷道所處位置原巖應(yīng)力為確定值,即視原巖垂直應(yīng)力和水平應(yīng)力在巷道及其圍巖所處空間范圍內(nèi)為常數(shù),并忽略圍巖自重的影響。深度H處的巖石所受應(yīng)力,根據(jù)海姆(Heim)的結(jié)論有:

式中:γ為巖石的容重;H為上覆蓋巖的深度;μ為泊松比;σν、σh為原巖的垂直和水平應(yīng)力。在巖體中開挖巷道,應(yīng)力重分布的范圍是有限的。彈性理論分析和現(xiàn)場測試表明,對于巷道開挖的應(yīng)力應(yīng)變,僅在其周圍距開挖空間中心3～5倍巷道開挖寬度(或高度)的范圍內(nèi)存在實(shí)際影響。在3倍寬度處的應(yīng)力變化一般在10%以下,在5倍寬度處的應(yīng)力變化一般在3%以下。所以設(shè)計的模型邊界寬度是開挖巷道寬度的5倍。建立20m×10m×20m的模型,巷道寬4m,墻高2m,半圓半徑為2m,建立的模型網(wǎng)格劃分見圖2。

圖2 巷道模型

4 模擬實(shí)驗(yàn)結(jié)果

(1)開挖前后巷道周邊的位移變化見圖3～8。

圖3 開挖前拱頂Z方向位移變化

(3)巷道圍巖開挖前后各部分監(jiān)測數(shù)據(jù)表。

如圖3、4、5、6、7和8可知,監(jiān)測點(diǎn)的垂直位移測試數(shù)據(jù)如表1。

表1 不同監(jiān)測點(diǎn)巷道圍巖的垂直位移

如圖9、10、11、12、13和14可知,監(jiān)測點(diǎn)的徑向應(yīng)力測試數(shù)據(jù)如表2。

表2 不同監(jiān)測點(diǎn)巷道圍巖的徑向應(yīng)力/pa

監(jiān)測點(diǎn)的周向應(yīng)力測試數(shù)據(jù)如表3。

表3 不同監(jiān)測點(diǎn)巷道圍巖的周向應(yīng)力/pa

5 數(shù)值模擬分析

(1)由圖3、4可知,巷道開挖前在自重應(yīng)力下達(dá)到平衡狀態(tài)過程中,拱頂垂直位移隨著時間的增長,不斷增大,最后曲線趨向于水平達(dá)到穩(wěn)定平衡。巷道開挖后,忽略巷道平衡產(chǎn)生的位移,考慮開挖的影響,位移開始走負(fù)方向發(fā)展。圖3、4同樣說明了巷道底部中間類似的情況。從這些可知,巷道受開挖的影響變形非常明顯。

(2)從圖7、8可知,巷道圍巖各點(diǎn)在開挖前后的位移變化。在開挖后,巷道圍巖的應(yīng)力集中相對明顯,位移加速變化較大。表1列舉了拱頂中間,兩邊墻中間,底部中間三點(diǎn)的位移,可以看出在清零后的位移變化還是比較大的。拱頂中間相對其他兩點(diǎn)的位移變化較大。

(3)從圖11、12、13和14可知,巷道圍巖在開挖前都是以壓應(yīng)力為主,但在開挖后許多部位應(yīng)力集中比較明顯,產(chǎn)生了拉應(yīng)力,因巖石的抗拉強(qiáng)度較小,容易破壞。

(4)巷道兩邊墻中間處的周向應(yīng)力、拱頂中間的徑向應(yīng)力和底部中間的徑向應(yīng)力基本上始終為0。這與巷道圍巖理論中的開挖空間周邊處徑向應(yīng)力始終為0的結(jié)論是相符的。

(5)巷道周邊產(chǎn)生塑性變形,隨著時間的推移,塑性變形越來越大,最終形成流變。特別是,拱頂中間、兩邊墻中間、底部中間等應(yīng)力集中比較明顯的地方。后處理工程中,應(yīng)該加強(qiáng)巷道支護(hù),防止塑性流變的發(fā)生。

6 結(jié)語

運(yùn)用數(shù)值模擬或數(shù)值實(shí)驗(yàn)方法,分析計算半圓拱形巷道在一定原巖應(yīng)力條件下開挖前后圍巖應(yīng)力場和位移場的變化和分布情況,并準(zhǔn)確確定巷道周邊某些點(diǎn)周向和徑向應(yīng)力的大小,無論從技術(shù)角度還是經(jīng)濟(jì)方面,都是十分可行的。

從理論上說,彈性介質(zhì)受力后其內(nèi)部的應(yīng)力分布是與介質(zhì)本身的力學(xué)性能無關(guān)。flac3D軟件都是建立在有限差原理基礎(chǔ)之上的,而有限差分計算是先通過動態(tài)顯示計算出結(jié)點(diǎn)位移的。因而對于具體的工程進(jìn)行數(shù)值模擬或數(shù)值實(shí)驗(yàn)之前,物理方法測定巖石介質(zhì)的彈性模量、泊松比、單軸抗壓強(qiáng)度等力學(xué)參數(shù)是必要的。

由于巖石是脆性材料,在其內(nèi)部各種宏觀或微觀節(jié)理裂隙作用下,其抗拉強(qiáng)度幾乎為0。對于確定的巷道工程,斷面形狀與尺寸初定后,根據(jù)原巖應(yīng)力側(cè)壓系數(shù)和巖石的力學(xué)性能參數(shù)等進(jìn)行數(shù)值模擬試驗(yàn),并修改設(shè)計,確保圍巖不產(chǎn)生周向拉應(yīng)力,在技術(shù)、安全和經(jīng)濟(jì)方面都是具有重要意義的。

REFERENCES

[1] 蔡美峰.金屬礦山采礦設(shè)計優(yōu)化與地壓控制——理論與實(shí)踐.北京:科學(xué)出版社,2001,4.

[2] 彭文斌.FLac 3D實(shí)用教程.北京:機(jī)械工業(yè)出版社,2007,8.

[3] 劉加冬,陸文,路洪斌,何偉.半圓拱形巷道圍巖應(yīng)力分布規(guī)律的研究[J].礦業(yè)快報,2008,3:22-25.

[4] 劉君,孔憲京.節(jié)理巖體中隧道開挖與支護(hù)的數(shù)值模擬[J].巖土力學(xué),2007,28(2):321-326.

Ground Pressure Formation M echan ism ofM in ing Tunnels and Research on StressD istribution

RAO Yun-zhang,WU Guo-xing,LUO Huan
(JiangxiUniversity of Science&Technology,Resources&Environment Engineering Institute ,Ganzhou 341000,Jiangxi,China)

Through analysis on dissertation of tunnel ground pressure issues by scholars at home and abroad,the article summarizes the for mingmechanis m ofmining tunnel ground pressure.Numerical simulation and analysis are made on surrounding rock stress status of the semicircle arch tunnel before and after excavation under the condition of the some lateralpressure coefficients by numerical analysis software and basic principle of the finite difference.The distribution rules of surrounding rock stress are illustrated after excavation of semicircle arch tunnels.

roadways pressure;numerical simulation;the finite difference;lateral pressure coefficient;surrounding rock stress

TD3

:A

:1009-3842(2010)01-0016-05

2009-11-24

饒運(yùn)章(1963-),男,江西會昌人,博士,教授,從事采礦工程與環(huán)境巖土工程方面的教學(xué)、研究工作。