閆永波　陸　文　張喜濤　王俊峰

(1.西南科技大學(xué)環(huán)境與資源學(xué)院;2.西昌學(xué)院工程技術(shù)學(xué)院)

?

最大水平主應(yīng)力對巷道圍巖穩(wěn)定性的影響研究

閆永波1陸文1張喜濤1王俊峰2

(1.西南科技大學(xué)環(huán)境與資源學(xué)院;2.西昌學(xué)院工程技術(shù)學(xué)院)

摘要基于金川三礦實測地應(yīng)力和某分段巷道,采用理論分析，得出了σHv型應(yīng)力場下巷道的最優(yōu)布置夾角,針對最大水平主應(yīng)力與巷道軸線呈不同夾角時的巷道布置，采用FLAC3D數(shù)值模擬軟件進行數(shù)值模擬,對其應(yīng)力分布、塑性區(qū)分布和圍巖位移進行對比研究。結(jié)果表明，在σHv型應(yīng)力場下，巷道的最優(yōu)布置方式是與最大水平主應(yīng)力成某一角度;當夾角達到一定程度時,巷道頂?shù)装鍟霈F(xiàn)大小相當?shù)沫h(huán)向應(yīng)力;引起巷道底板變形破壞的重要原因是水平應(yīng)力。

關(guān)鍵詞最大水平主應(yīng)力數(shù)值模擬圍巖位移

地下工程的穩(wěn)定性取決于圍巖性質(zhì)、原巖應(yīng)力場的分布、工程斷面形狀、以及工程支護方式等因素。地應(yīng)力是引起地下工程圍巖變形和破壞的主要因素,其大小和方向直接關(guān)系著地下工程圍巖的穩(wěn)定與否[1-4]。相關(guān)資料[5-6]顯示,巖層中的水平應(yīng)力在多數(shù)情況下大于垂直應(yīng)力,且具有明顯的方向性，巷道底鼓、頂板沉降和兩幫內(nèi)擠與其密切相關(guān)。鄭書兵[7]、張延新等人[8]的研究也都說明了巷道的正確布置有利于其穩(wěn)定。

在實際施工過程中，巷道布置常常難以與原巖應(yīng)力場相適應(yīng)，一是巷道的布置決定于礦床開拓方式與開采方法,難以保證巷道的走向與最大主應(yīng)力方向始終一致或垂直；二是原巖在漫長的地質(zhì)年代中經(jīng)歷多次地質(zhì)構(gòu)造運動,原巖應(yīng)力是自重應(yīng)力與歷次構(gòu)造運動產(chǎn)生的殘余應(yīng)力疊加的結(jié)果,最大主應(yīng)力的大小和方向是變化的。因此,原巖最大主應(yīng)力的方向與巷道走向(軸線方向)的夾角對于巷道圍巖穩(wěn)定的影響,便成了重要的研究課題。本文基于金川集團地下礦山實測原巖應(yīng)力,通過理論分析和對某分段巷道進行數(shù)值模擬,探索巷道軸線與最大水平主應(yīng)力的夾角對巷道穩(wěn)定性的影響。

1巷道的最優(yōu)布置理論分析

原巖應(yīng)力實測結(jié)果表明：金川集團三礦區(qū)地下巖層的水平應(yīng)力一般大于垂直應(yīng)力,且水平應(yīng)力具有明顯的方向性;最大水平主應(yīng)力與最小水平主應(yīng)力的比值為1.5～2.5。

巷道頂?shù)装宓姆€(wěn)定性與水平應(yīng)力之間有以下特點：①當巷道軸線與最大水平主應(yīng)力平行時,頂?shù)装宸€(wěn)定性最好;②巷道軸線與最大水平主應(yīng)力相交時,其頂板變形破壞偏向巷道某一幫;③巷道軸線與最大水平主應(yīng)力垂直時,頂?shù)装宸€(wěn)定性最差。根據(jù)彈性力學(xué)的孔口理論,巷道截面尺寸要與原巖的正應(yīng)力大小相適應(yīng),最好實現(xiàn)等應(yīng)力軸比。金川公司地下礦山采場分段巷道為直墻拱形,寬度與高度均為5.0 m左右,因此巷道最優(yōu)布置應(yīng)該使原巖水平正應(yīng)力與垂直主應(yīng)力的比值等于1,即σn/σv=1。

按照最大水平主應(yīng)力σH、最小水平主應(yīng)力σh和垂直主應(yīng)力σv三者數(shù)值的大小關(guān)系,將巖層中的地應(yīng)力場分成3種形式,即σH型,σH>σh>σv;σHv型,σH>σv>σh;σv型,σv>σH>σh。張重遠[9]等人在金川三礦區(qū)的地應(yīng)力測量成果表明,三礦區(qū)最大水平主應(yīng)力約是最小水平主應(yīng)力的2倍,最大水平主應(yīng)力σH=14.7 MPa,最小水平主應(yīng)力σh=7.2 MPa,平均垂直主應(yīng)力σv=11.9 MPa,σH>σv>σh，屬于σHv型，可以滿足σn/σv=1。任意平面上的法向應(yīng)力與最小和最大水平主應(yīng)力之間的表達式為

(1)

式中,α為巷道軸線與最大水平主應(yīng)力方向的夾角，(°)。

將σn/σv=1代入式(1),可解得巷道布置的最優(yōu)夾角：

(2)

將測量的地應(yīng)力結(jié)果代入式(2),解得金川集團三礦區(qū)巷道的最優(yōu)布置夾角為53°。

為進一步尋找巷道布置的最優(yōu)夾角與原巖地應(yīng)力場之間的關(guān)系,現(xiàn)設(shè)垂直主應(yīng)力為12 MPa,通過式(2)，可計算出在不同的應(yīng)力場狀態(tài)下時，巷道在巖層中布置的最優(yōu)夾角,結(jié)果見表1。

表1　不同應(yīng)力狀態(tài)時巷道布置最優(yōu)夾角(σv =12 MPa)

當最大水平主應(yīng)力固定不變時,巷道布置的最優(yōu)夾角隨著最小水平應(yīng)力的減小而增大,即巷道布置的最優(yōu)夾角偏離最大水平主應(yīng)力方向。

2數(shù)值模擬及結(jié)果分析

2.1數(shù)值模型建立與模擬方案

以金川集團三礦某分段巷道為模擬對象,運用美國Itasca公司開發(fā)的三維有限差分FLAC3D軟件進行數(shù)值模擬。巷道斷面形狀為直墻拱形,巷道尺寸(寬×高)為5.0 m×4.6 m。根據(jù)彈性力學(xué)孔口理論,地下巷道開挖引起巖體中應(yīng)力重分布的范圍為4倍的開挖跨度和高度。因此綜合考慮巷道穩(wěn)定性的各種影響因素[10],以巷道拱圓心為原點,確定所建立的計算模型尺寸(寬×長×高)為50 m×30 m×50 m,共劃分59 520個單元,62 031個節(jié)點。模型采用Mohr-coulomb破壞準則,考慮到巷道埋深遠遠大于20倍巷道半徑,故采用應(yīng)力邊界控制,忽略所建模型內(nèi)巖體自重。

為了模擬最大水平主應(yīng)力與巷道軸線不同夾角(0°～90°)下巷道圍巖的應(yīng)力分布和變形破壞規(guī)律,每隔10°建立一組模型,共10組模型。

2.2數(shù)值模擬參數(shù)選擇

根據(jù)張重遠等人的地應(yīng)力測量成果,金川三礦數(shù)值模擬應(yīng)力場參數(shù)選取為最大水平主應(yīng)力為14.7 MPa,最小水平主應(yīng)力為7.2 MPa,垂直主應(yīng)力為11.9 MPa。根據(jù)實驗室所做試驗和前人的模擬經(jīng)驗選擇巖石力學(xué)參數(shù),見表2。

表2　巖石力學(xué)參數(shù)

2.3模擬結(jié)果分析

圓形巷道軸對稱問題的線彈性解析表明,巷道周邊將產(chǎn)生環(huán)向應(yīng)力的最大值和徑向應(yīng)力的最小值。直墻拱形巷道與圓形巷道在形狀上存在一定的差距,但在彈塑性條件下，直墻拱形巷道圍巖的水平應(yīng)力和垂直應(yīng)力的數(shù)值解與圓形巷道的解析解有一定類似之處,因此在分析應(yīng)力時可將頂?shù)装鍑鷰r中水平應(yīng)力以及幫部圍巖中垂直應(yīng)力近似看作環(huán)向應(yīng)力;頂?shù)讎鷰r中垂直應(yīng)力和幫部圍巖中水平應(yīng)力可近似看作徑向應(yīng)力。

2.3.1圍巖應(yīng)力分析

在不同角度下巷道布置的數(shù)值模擬結(jié)果見圖1～圖3。

圖1　不同夾角下圍巖水平應(yīng)力的分布

圖2　不同夾角下圍巖垂直主應(yīng)力的分布

圖3　不同夾角下巷道圍巖應(yīng)力集中系數(shù)

(1)由圖1可得，隨著最大水平主應(yīng)力與巷道軸線夾角的增大,0°～60°底板產(chǎn)生的環(huán)向拉應(yīng)力數(shù)值不斷減小,受拉區(qū)域也持續(xù)減小,由連續(xù)受拉變?yōu)椴糠质芾?并且在底板呈類似于反拱形，由于拉伸作用使得巷道圍巖整體性能降低,在巷道底板產(chǎn)生徑向裂隙,造成底鼓破壞；70°～90°底板所受的環(huán)向應(yīng)力逐漸變大,在90°時環(huán)向拉應(yīng)力具有最大值(0.25 MPa)；幫部只在最大水平主應(yīng)力與巷道軸線平行和垂直時產(chǎn)生徑向拉應(yīng)力,數(shù)值分別為0.026 MPa和0.02 MPa；在整個數(shù)值模擬過程中頂板沒有產(chǎn)生環(huán)向拉應(yīng)力,一直處于受壓狀態(tài)。

(2)分析圖2可看出,在整個數(shù)值模擬過程中兩幫沒有產(chǎn)生環(huán)向拉應(yīng)力,頂板也沒有產(chǎn)生徑向拉應(yīng)力;底板一直受到徑向拉應(yīng)力作用,這是因為頂板為拱形,有利于降低應(yīng)力集中,底板為直長邊,容易產(chǎn)生較大的應(yīng)力集中。隨著最大水平主應(yīng)力與巷道軸線夾角的逐漸增大,拉應(yīng)力先減小后增大,在50°時最小,為0.011 MPa,0°和90°時分別為 0.043和0.042 MPa。

(3)由圖3可以看出，頂?shù)装瀛h(huán)向應(yīng)力集中系數(shù)隨著最大水平主應(yīng)力與巷道軸線夾角的增大而持續(xù)減小。當夾角為50°時,應(yīng)力集中系數(shù)為1.06,隨夾角的繼續(xù)增大到70°時，應(yīng)力集中系數(shù)急劇上升到1.59,在此之后較為穩(wěn)定,沒有出現(xiàn)大幅度增加；兩幫環(huán)向應(yīng)力集中系數(shù)總體上呈先減小后小幅上升的趨勢,在60°時達到最小值為1.34,但垂直時的應(yīng)力集中系數(shù)值小于平行時的應(yīng)力集中系數(shù)值。

綜合來看,在σHv型應(yīng)力場下,巷道圍巖的穩(wěn)定性隨著夾角(0°～90°)的改變呈現(xiàn)出復(fù)雜的變化；0°～50°時,最小水平主應(yīng)力對頂?shù)装宓沫h(huán)向應(yīng)力和兩幫的徑向應(yīng)力起決定性作用,這時頂?shù)装逅艿沫h(huán)向集中應(yīng)力隨著夾角的增大,其影響作用逐漸減小,巷道圍巖的穩(wěn)定性較小夾角時得到了提高；50°～90°時,最大水平主應(yīng)力在巷道圍巖的穩(wěn)定性中開始逐漸占據(jù)主導(dǎo)地位,此時頂?shù)装迨艿降沫h(huán)向集中應(yīng)力隨著最大水平主應(yīng)力與巷道軸線夾角的增大開始增大,頂?shù)装宓乃苄詤^(qū)也隨著增大。

2.3.2塑性區(qū)

當最大水平主應(yīng)力與巷道軸線夾角為0°～90°時,塑性區(qū)會出現(xiàn)以下幾種特征：0°～20°時，幫部的塑性區(qū)大于頂?shù)撞克苄詤^(qū),但總體巷道周邊塑性區(qū)沒有出現(xiàn)大的變化;30°～60°時，隨著夾角的增加，幫部塑性區(qū)呈增大趨勢,頂?shù)撞克苄詤^(qū)逐漸變小;70°～90°時，幫部以及頂?shù)撞克苄詤^(qū)變化情況恰好與30°～60°相反,但底板塑性區(qū)要大于頂板塑性區(qū)。

2.3.3巷道圍巖位移

在巷道頂、底板和兩幫中點布置監(jiān)測點,以確定巷道圍巖在最大水平主應(yīng)力與巷道軸線不同夾角下的收斂量變化情況,根據(jù)監(jiān)測結(jié)果繪出圍巖收斂量與不同夾角的關(guān)系曲線,見圖4。

圖4　圍巖收斂量與不同夾角的關(guān)系

整體分析來看, 幫部和頂?shù)装宓膰鷰r收斂量隨著夾角的增加呈先減小后增大的趨勢。當0°～10°時,圍巖收斂量變化并不明顯;當10°～50°時,巷道周邊圍巖收斂量持續(xù)減小,但幫部收斂量變化速率比頂?shù)装宓氖諗苛孔兓俾蚀?這主要是因為原巖水平主應(yīng)力值逐漸趨近于垂直主應(yīng)力值,巷道圍巖穩(wěn)定性比其他夾角下的穩(wěn)定性好，在夾角為50°時巷道圍巖收斂量最小,兩幫收斂量為6.07 cm,頂?shù)装迨諗苛繛?8.55 cm；隨著夾角的繼續(xù)增大,巷道周邊圍巖收斂量急劇增大，當達到70°～90°時,巷道收斂量增幅減緩,到90°時兩幫和頂?shù)装宓氖諗苛糠謩e為43.82，52.86 cm。

從以上分析可知,金川三礦巷道最優(yōu)布置方式應(yīng)該是最大水平主應(yīng)力與巷道軸線成50°附近,這和理論分析的53°最優(yōu)布置夾角一致。頂?shù)装宓氖諗孔冃屋^兩幫嚴重,在后期建設(shè)過程中應(yīng)加強頂?shù)装逯ёo。

3結(jié)論

(1)在σHv型應(yīng)力場中,巷道的最優(yōu)布置夾角并不是巷道軸線與最大水平主應(yīng)力成0°,而是由最大水平主應(yīng)力、最小水平主應(yīng)力和垂直主應(yīng)力三者的大小決定。在本次計算模型中,當巷道軸線與最大水平主應(yīng)力成53°時,圍巖的穩(wěn)定性最好。

(2)在σHv型應(yīng)力場中,當夾角小于巷道布置最優(yōu)夾角時,最小水平主應(yīng)力對頂?shù)装宓沫h(huán)向應(yīng)力分布起決定性作用；當夾角大于巷道布置最優(yōu)夾角時,對頂?shù)装瀛h(huán)向應(yīng)力分布起決定性作用的則為最大水平主應(yīng)力，且頂?shù)装鍟霈F(xiàn)大小相當?shù)沫h(huán)向應(yīng)力集中,造成巷道底鼓和頂板沉降。

(3)金川公司三礦巷道的變形與破壞主要表現(xiàn)為底鼓與頂板沉降,因此應(yīng)注重加強巷道頂?shù)装逯ёo。

參考文獻

[1]宋志敏,程增慶,張生華.構(gòu)造應(yīng)力區(qū)軟巖巷道圍巖變形與控制[J].礦山壓力與頂板管理,2005,22(4):48-50.

[2]Kushwaha A, Murali M, Singh S,etal. Effect of in situ horizontal stresses on roadway stability[J].Journal of Mines, Metals and Fuels, 2003,51(3):134-142.

[3]王永才,康紅普.金川礦山深井高應(yīng)力開采潛在的問題[J].中國礦業(yè),2010,9(12):52-55.

[4]魯巖,鄒喜正,劉長友,等.構(gòu)造應(yīng)力場中的巷道布置[J].采礦與安全工程學(xué)報,2008,25(2):144-149.

[5]倪興華.地應(yīng)力研究與應(yīng)用[M].北京:煤炭工業(yè)出版社,2007.

[6]康紅普.煤礦井下應(yīng)力場類型及相互作用分析[J].煤炭學(xué)報,2008,33(12):1329-1335.

[7]鄭書兵.寺河煤礦三維地應(yīng)力場分布和巷道布置優(yōu)化[J].煤炭學(xué)報,2010(5):717-722.

[8]張延新,蔡美峰,歐陽振華.地應(yīng)力與巷道布置關(guān)系的理論研究[J].巖土工程技術(shù),2005,19(2):93-97.

[9]張重遠,吳滿路,廖椿庭.金川三礦地應(yīng)力測量及應(yīng)力狀態(tài)特征研究[J].巖土力學(xué),2013,34(11):3254-3260.

[10]謝文兵,陳曉祥,鄭百生.采礦工程問題數(shù)值模擬研究與分析[M].徐州:中國礦業(yè)大學(xué)出版社,2005.

(收稿日期2015-11-01)

Influence of the Maximum Horizontal Principal Stress to the Stability of the Roadway Surrounding Rock

Yan Yongbo1Lu Wen1Zhang Xitao1Wang Junfeng2

(1.School of Environment and Resource, Southwest University of Science & Technology;2.School of Engineering Science,Xichang College)

AbstractBased on the actual measured ground stress of Jinchuan third Mine and a segment roadway,by theoretical analysis method, the optimal placement angle of roadway under the σHv-type stress filed is obtained. Based on FLAC3Dsimulation software,the roadway arrangement is simulated under the condition of the difference angle of maximum horizontal principal stress and roadway axis, besides that,the stress distribution, plastic zone distribution and surrounding rock displacement are conducted contrast analysis. The result show that in the σHvtype stress field, the optimal arrangement of the roadway is forming a certain angle with the maximum horizontal stress; when the angle reaches a certain level,the roof and floor of the roadway will appear the hoop stress, the horizontal stress is the main factor that cause the deformation and destruction of the roadway floor.

KeywordsMaximum horizontal principal stress, Numerical simulation, Surrounding rock displacement

閆永波(1990—),男,碩士研究生,621010 四川省綿陽市涪城區(qū)青龍大道中段59號。