劉 寅,馮 飛,屈 潔
(中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司,西安 710065)
非洲某水電站母線(xiàn)豎井圍巖穩(wěn)定分析
劉 寅,馮 飛,屈 潔
(中國(guó)電建集團(tuán)西北勘測(cè)設(shè)計(jì)研究院有限公司,西安 710065)
結(jié)合非洲某工程實(shí)際,利用PHASE2軟件,通過(guò)平面計(jì)算和軸對(duì)稱(chēng)計(jì)算,對(duì)此項(xiàng)目長(zhǎng)達(dá)142.3 m的母線(xiàn)豎井進(jìn)行圍巖穩(wěn)定分析。根據(jù)計(jì)算結(jié)果,確定不同類(lèi)別圍巖的支護(hù)形式,為豎井圍巖穩(wěn)定性評(píng)價(jià)和支護(hù)設(shè)計(jì)方案的選擇提供了一種新的方法。關(guān)鍵詞:母線(xiàn)豎井;圍巖穩(wěn)定;支護(hù)形式;有限元計(jì)算;PHASE2
隨著水利水電工程的發(fā)展。作為一種典型的廠(chǎng)房結(jié)構(gòu),地下廠(chǎng)房已經(jīng)日益普遍地出現(xiàn)在各類(lèi)水電站之中。其中圍巖穩(wěn)定問(wèn)題是地下洞室中最為突出的問(wèn)題之一。關(guān)于地下廠(chǎng)房主變的布置,常規(guī)有2種形式,一種將主變布置在平行于主廠(chǎng)房的主變洞內(nèi),另一種則通過(guò)母線(xiàn)豎井將母線(xiàn)引出至地面,接入位于地面的主變。本文以非洲某水電站長(zhǎng)達(dá)142.3 m的豎井為模型,通過(guò)有限元圍巖穩(wěn)定分析,確定其支護(hù)形式。
PHASE2是由加拿大Rocscience公司推出的巖土專(zhuān)用計(jì)算軟件,該軟件主要用于計(jì)算邊坡及地下洞室開(kāi)挖后巖土應(yīng)力、應(yīng)變和位移狀況。軟件特性:可分析平面應(yīng)變問(wèn)題和空間軸對(duì)稱(chēng)問(wèn)題;材料的本構(gòu)模型包括了彈性模型、Mohr-Coulomb模型[2]和Hoek-Brown模型[3-4]等備用的模型;可模擬多種具有不同力學(xué)特性的材料;可以進(jìn)行不同開(kāi)挖階段的模擬[5];可考慮均勻場(chǎng)力單獨(dú)作用[6]、重力單獨(dú)作用和二者同時(shí)作用下的地應(yīng)力分布[7];可以對(duì)開(kāi)挖面噴混凝土、開(kāi)挖面施加錨桿、錨索等支護(hù)加固措施的效果進(jìn)行分析計(jì)算[8]。
以往常用的平面計(jì)算只能在二維層面考慮問(wèn)題,這種情況下,在空間上的問(wèn)題得不到解決。例如,對(duì)于一個(gè)豎井,平面計(jì)算只能考慮豎井在某一高程、一個(gè)特定斷面的結(jié)果,但是無(wú)法考慮此斷面高程以上或以下范圍內(nèi)的影響,往往會(huì)造成計(jì)算結(jié)果偏小。
因此對(duì)于圓形的豎井,且四周?chē)鷰r類(lèi)別基本成水平分布,可以考慮用軸對(duì)稱(chēng)計(jì)算來(lái)模擬此類(lèi)情況。軸對(duì)稱(chēng)計(jì)算是用二維模型來(lái)模擬一個(gè)關(guān)于某個(gè)軸圓周對(duì)稱(chēng)的三維模型??梢杂幂^少的計(jì)算資源(二維單元),模擬三維問(wèn)題,考慮空間上整個(gè)結(jié)構(gòu)的穩(wěn)定性。此計(jì)算方法雖然在適用條件上不如平面模型廣泛,但只要滿(mǎn)足其邊界條件,計(jì)算結(jié)果較平面模型更加趨于實(shí)際情況。
本工程母線(xiàn)豎井總長(zhǎng)142.3 m,外徑7 m,下部連接廠(chǎng)房?jī)?nèi)部母線(xiàn)室,上部連通主變平臺(tái),母線(xiàn)豎井部位地表堆積殘坡積塊石碎石土,結(jié)構(gòu)松散,局部架空,厚度小于15 m,需清除;母線(xiàn)豎井自上而下穿過(guò)強(qiáng)風(fēng)化石英黑云母片麻巖和石英巖、中風(fēng)化石英巖以及微風(fēng)化~新鮮的石英黑云母片麻巖,片麻理發(fā)育,片麻理產(chǎn)狀因巖體揉皺而變化較大,主要有NE70°~83°SE∠20°~30°。母線(xiàn)豎井縱向需要穿越強(qiáng)風(fēng)化、中風(fēng)化和弱風(fēng)化3類(lèi)圍巖。根據(jù)地質(zhì)鉆孔資料顯示,在41.5 m以上為強(qiáng)風(fēng)化;強(qiáng)風(fēng)化與中風(fēng)化之間,有厚度約為1 m的軟弱夾層帶;中風(fēng)化厚度約為33 m,以下均為弱風(fēng)化圍巖。各類(lèi)圍巖的參數(shù)見(jiàn)表1。對(duì)于計(jì)算中各種支護(hù)的參數(shù)取值,見(jiàn)表2。
表1 圍巖力學(xué)參數(shù)表
表2 支護(hù)參數(shù)表
初步設(shè)計(jì)對(duì)于強(qiáng)風(fēng)化圍巖(Ⅳ類(lèi)圍巖)采取布置長(zhǎng)3 m、間距1.5 m、直徑25 mm的錨桿支護(hù),并設(shè)有100 mm厚掛網(wǎng)噴混凝土;對(duì)于中風(fēng)化圍巖(Ⅲ類(lèi)圍巖)采取布置長(zhǎng)3 m、間距2 m、直徑25 mm的錨桿支護(hù),并設(shè)有100 mm厚掛網(wǎng)噴混凝土;對(duì)于弱風(fēng)化圍巖(Ⅱ類(lèi)圍巖)僅進(jìn)行100 mm厚掛網(wǎng)噴混凝土。
計(jì)算使用PHASE2軟件進(jìn)行平面計(jì)算,采用摩爾庫(kù)倫非線(xiàn)性模型進(jìn)行巖體參數(shù)的定義,計(jì)算模型四周采用法向約束。根據(jù)以往的計(jì)算經(jīng)驗(yàn),當(dāng)計(jì)算模型四周的厚度達(dá)到一定程度,就不會(huì)對(duì)計(jì)算結(jié)果產(chǎn)生影響。因此模型在高度方向和長(zhǎng)度方向上均取為5倍的豎井開(kāi)挖直徑,即為40 m。整個(gè)模型采用三角形單元?jiǎng)澐?,共有?jié)點(diǎn)2 097個(gè),單元4 354個(gè),計(jì)算模型見(jiàn)圖1。
為了計(jì)算開(kāi)挖后圍巖穩(wěn)定性,需要考慮初始應(yīng)力。在平面模型中,計(jì)算的初始應(yīng)力由豎向地應(yīng)力得到。根據(jù)圍巖深度的不同,不同類(lèi)別圍巖計(jì)算時(shí)的初始應(yīng)力場(chǎng)有所不同。
圖1 母線(xiàn)豎井縱剖面圖
根據(jù)圍巖深度,圍巖的豎向應(yīng)力大小為:
(1)
式中:σv為豎向應(yīng)力;γ為巖石容重;H為巖體深度。
水平向圍巖壓力為:
(2)
式中:σh為豎向應(yīng)力;μ為巖石泊松比。
計(jì)算模型輸入各類(lèi)荷載和支護(hù)條件后,其計(jì)算結(jié)果分析如下:
(1) 強(qiáng)風(fēng)化圍巖在有支護(hù)的情況下,出現(xiàn)的最大壓應(yīng)力為0.76 MPa,遠(yuǎn)小于圍巖的強(qiáng)度。豎井的最大位移值為0.45 mm,此數(shù)值很小。因此認(rèn)為對(duì)于強(qiáng)風(fēng)化圍巖,其支護(hù)形式的選取能夠滿(mǎn)足要求,圍巖在支護(hù)后穩(wěn)定。
(2) 中風(fēng)化圍巖在有支護(hù)的情況下,出現(xiàn)的最大壓應(yīng)力為1.17 MPa,遠(yuǎn)小于圍巖的強(qiáng)度。豎井的最大位移值為0.28 mm,此數(shù)值很小。因此認(rèn)為對(duì)于中風(fēng)化圍巖,其支護(hù)形式的選取能夠滿(mǎn)足要求,圍巖在支護(hù)后穩(wěn)定。
(3) 弱風(fēng)化圍巖在有支護(hù)的情況下,出現(xiàn)的最大壓應(yīng)力為1.8 MPa,遠(yuǎn)小于圍巖的強(qiáng)度。豎井的最大位移值為0.3 mm,此數(shù)值很小。因此認(rèn)為對(duì)于弱風(fēng)化圍巖,其支護(hù)形式的選取能夠滿(mǎn)足要求,圍巖在支護(hù)后穩(wěn)定。
根據(jù)以上計(jì)算結(jié)果分析,圍巖開(kāi)挖后,無(wú)論是第1主應(yīng)力還是第2主應(yīng)力的數(shù)值都不大,且圍巖的變形值很低,并且塑性區(qū)幾乎不可見(jiàn)。因此認(rèn)為支護(hù)設(shè)計(jì)基本合理,圍巖在開(kāi)挖后不存在失穩(wěn)的危險(xiǎn)。
軸對(duì)稱(chēng)計(jì)算模型,可以用二維模型來(lái)模擬一個(gè)關(guān)于某個(gè)軸圓周對(duì)稱(chēng)的三維模型。這種計(jì)算比二維平面計(jì)算要更加準(zhǔn)確,特別是能夠反映在不同高度下圍巖的穩(wěn)定特性。
根據(jù)實(shí)際設(shè)計(jì)的母線(xiàn)豎井,建立軸對(duì)稱(chēng)模型,對(duì)稱(chēng)軸為母線(xiàn)豎井中心線(xiàn)。
計(jì)算的材料依然使用摩爾庫(kù)倫模型,模型底部和兩側(cè)都使用法向約束。
通過(guò)計(jì)算結(jié)果可以看出,最大壓應(yīng)力為5 MPa,遠(yuǎn)小于圍巖的強(qiáng)度。豎井的最大位移值為3 mm,此數(shù)值很小。因此認(rèn)為對(duì)于強(qiáng)風(fēng)化圍巖,其支護(hù)形式的選取是能夠滿(mǎn)足要求的,圍巖在支護(hù)后是穩(wěn)定的。
對(duì)比平面計(jì)算結(jié)果,軸對(duì)稱(chēng)模型計(jì)算的結(jié)果相對(duì)較大,但總體數(shù)值比較接近,所得結(jié)論也是一致的。此類(lèi)計(jì)算結(jié)果同圍巖參數(shù)、豎井的深度和直徑關(guān)系密切,可以看出,如果僅考慮平面計(jì)算結(jié)果,在一些地質(zhì)條件較差、地應(yīng)力較大的區(qū)域,也許會(huì)得到與實(shí)際情況相違背的結(jié)果。但通過(guò)軸對(duì)稱(chēng)計(jì)算,可以充分地反映出三維空間上圍巖的穩(wěn)定性,其計(jì)算結(jié)果要更加真實(shí)。
根據(jù)平面和軸對(duì)稱(chēng)模型的計(jì)算結(jié)果,對(duì)于圍巖進(jìn)行的支護(hù)設(shè)計(jì)是合適的,能夠滿(mǎn)足安全要求。此工程中,雖然2種計(jì)算所得結(jié)果相同,但可以看出軸對(duì)稱(chēng)模型計(jì)算結(jié)果較平面模型計(jì)算結(jié)果大。這種大是考慮了三維空間影響得到的結(jié)論,是合理的,更加趨向于真實(shí)的情況。
平面計(jì)算模型可以在任何條件下使用,其適應(yīng)性較強(qiáng)。而對(duì)于軸對(duì)稱(chēng)模型,首先必須是圓形隧洞,并且在周邊是對(duì)稱(chēng)的,一旦在計(jì)算范圍內(nèi)有其他洞室存在,就無(wú)法考慮這些因素,會(huì)使得計(jì)算的邊界條件不準(zhǔn)確,無(wú)法保證計(jì)算結(jié)果的真實(shí)性。
因此,在后續(xù)工程中,在滿(mǎn)足軸對(duì)稱(chēng)模型計(jì)算邊界條件要求的前提下,使用此種計(jì)算能夠更加真實(shí)地反映實(shí)際情況。如果邊界條件不能夠滿(mǎn)足,可以用平面模型作為基本計(jì)算,軸對(duì)稱(chēng)計(jì)算作為復(fù)核計(jì)算,在相互驗(yàn)證之后,找到最為經(jīng)濟(jì)和合理的支護(hù)設(shè)計(jì),從而收到良好的經(jīng)濟(jì)效益。
[1] 華東水利學(xué)院.巖石力學(xué)[M].北京:水利電力出版社,1986.
[2] 韓瑞庚.地下工程新奧義[M].北京:科學(xué)出版社,1987.
[3] HOEK E,BROWN E T. Practical estimates of rack mass strength[J].International Journal of Rock Mechanics and Mining Science,1998,38(08):1165-1186.
[4] R.K.Goel,Bhawani Singh.Rock mass classification [M].ELsevier Science Ltd,1999.
[5] 李杰,郭海燕,劉凌云,牟桂杰.丈八口隧道圍巖穩(wěn)定分析[J].巖土力學(xué),2004(S2): 537-541.[6] 王立忠,胡亞元,王百林,陳云敏.崩塌松散圍巖隧道施工穩(wěn)定性分析與監(jiān)控[J].巖石力學(xué)與工程學(xué)報(bào),2003,22(04):79-85.[7] 何欣,曹懷園,劉永智,劉潔玉.Phase2軟件在隧洞開(kāi)挖圍巖支護(hù)時(shí)機(jī)中的應(yīng)用[J].西北水電,2015(03):57-61.
[8] 溫森,楊圣奇.基于Hoek-Brown準(zhǔn)側(cè)的隧洞圍巖變形研究[J].巖土力學(xué),2011,32(01):70-76.
[9] 高凱,王林維,郭建宏.云南某水電站地下洞室圍巖穩(wěn)定性分析[J].西北水電,2015(05):33-37.
[10] 石高峰.基于Phase2的地下廠(chǎng)房及尾調(diào)室開(kāi)挖與支護(hù)設(shè)計(jì)[J].水力發(fā)電,2015(12):63-67.
Analysis on Stability of Surrounding Rock of Basbar Shaft
LIU Yin, FENG Fei, QU Jie
(Northwest Engineering Corporation Limited, Xi'an 710065,China)
In combination with the practice of one African project and by application of PHASE 2 software, stability of the surrounding rock of the 142.3 m long busbar shaft of the project is analyzed through plan calculation and axis symmetry calculation. According to the calculation results, supports modes for different surrounding rocks are finalized, this provides stability assessment and support design scheme of the shaft surrounding rock with a new method.Key words:basbar shaft; stability of surrounding rock; support mode; calculation by finite element method; PHASE2
1006—2610(2016)06—0094—03
2016-11-18
劉寅(1987- ),男,河北省秦皇島市人,助理工程師,主要從事水利水電工程設(shè)計(jì)工作.
TV554;TU457
A
10.3969/j.issn.1006-2610.2016.06.024