單 威, 李建東

(內(nèi)蒙古自治區(qū) 水利水電勘測設(shè)計院,內(nèi)蒙古 呼和浩特 010020)

隆盛水庫溢洪道高邊坡FLAC3D數(shù)值模擬穩(wěn)定性分析

單 威, 李建東

(內(nèi)蒙古自治區(qū) 水利水電勘測設(shè)計院,內(nèi)蒙古 呼和浩特 010020)

隆盛水庫是內(nèi)蒙古大黑河一級支流拐角鋪河上在建的一座中型控制性水利樞紐工程,是調(diào)洪蓄水和防止泥沙進入大黑河的重要水利工程。因順坡向緩傾角節(jié)理裂隙及混合巖化作用發(fā)育,水庫溢洪道右岸高邊坡在開挖后修坡過程中,發(fā)生了兩次大規(guī)?；?在收集該地區(qū)大量的野外調(diào)查資料的基礎(chǔ)上,運用FLAC3D軟件模擬分析隆盛水庫溢洪道高邊坡的穩(wěn)定性,對工程區(qū)邊坡開挖和修坡條件下的土體變形及位移情況進行模擬計算,獲得了邊坡在開挖和修坡條件下的最大、最小主應(yīng)力圖和一系列觀測點位移圖,為工程設(shè)計提供參考。

高邊坡;FLAC3D模擬計算;穩(wěn)定性分析;隆盛水庫

1 研究背景

隆盛水庫是內(nèi)蒙古大黑河一級支流拐角鋪河上在建的一座中型控制性水利樞紐工程。工程的任務(wù)是灌溉、城鎮(zhèn)及工業(yè)供水、攔沙、防洪等綜合利用,根據(jù)調(diào)節(jié)計算每年灌溉可供水量為658萬m3,城鎮(zhèn)及工業(yè)供水可供水量為380萬m3。工程等別為Ⅲ等,主要建筑物包括大壩、溢洪道、泄洪排沙洞及供水洞等,主要建筑級別為3級建筑物,次要建筑物為4級,施工臨時建筑物級別為5級。

因順坡向緩傾角節(jié)理裂隙及混合巖化作用發(fā)育,水庫溢洪道右岸高邊坡在開挖后修坡過程中,發(fā)生了兩次大規(guī)?；?造成一、二級馬道嚴重損毀。本文以分析確定該高邊坡潛在不穩(wěn)定滑動面及應(yīng)力集中區(qū)域分布范圍為主要研究目的,運用FLAC3D軟件建模分析,對工程區(qū)邊坡開挖和修坡條件下的土體變形及位移情況進行模擬計算,獲得了邊坡在開挖和修坡條件下的最大、最小主應(yīng)力圖和一系列觀測點位移圖,為工程設(shè)計邊坡加固處理措施提供參考。

2 研究區(qū)概況

2.1 工程概況

隆盛水庫溢洪道設(shè)計方式為開敞式溢洪道,由進水渠段、控制段、泄槽段、挑流消能段組成。目前,溢洪道開挖進水渠段已至設(shè)計高程1 461.50 m,由控制段至泄槽段現(xiàn)狀底板高程為1 458.78～1 448.68 m,右部邊坡分為五級馬道,一—四級馬道邊坡高度15 m左右,四—五級邊坡高度約30 m,五級馬道以上邊坡高度10～15 m,工程開挖坡高在100 m左右,屬工程開挖高邊坡。每級邊坡坡比為1∶0.75,二級以上邊坡已完成噴錨處理。

溢洪道右岸一級邊坡開挖修坡過程中發(fā)生較大的邊坡滑塌兩次,一次是發(fā)生在2015年8月15日,塌坡范圍樁號0+080～0+120,高度為一級馬道邊坡坡腳(高程1471.0 m左右)—二級馬道頂部(高程1 486.0 m)。另一次較大的滑塌發(fā)生在2015年8月31日。在2015年8月30日發(fā)現(xiàn)一級馬道樁號0+020～0+040發(fā)生局部塌落,樁號0+060～0+100一級馬道發(fā)生縱向鋸齒狀拉裂縫,裂縫寬度10～15 cm,地質(zhì)設(shè)代向監(jiān)理反映了相關(guān)情況,并建議停止開挖。次日,樁號0+060～0+100段發(fā)生較大邊坡滑塌,塌坡高度為一級馬道頂部(高程1 471.0 m左右)—一級馬道邊坡坡腳(高程1 458.0 m左右)。

2.2 右岸高邊坡巖性特征

一級邊坡底部巖性為下太古界烏拉山群大理巖組(Ar1wl2)黑云母角閃斜長片麻巖,分布厚度8～10 m,片麻理非常發(fā)育,節(jié)理發(fā)育,巖體呈碎裂結(jié)構(gòu)。上部巖性0-020～0+005段為下太古界烏拉山群大理巖組(Ar1wl2)白云石大理巖,分布厚度為3～5 m,節(jié)理較發(fā)育,巖體呈次塊狀結(jié)構(gòu);0+005～0+028段為下太古界烏拉山群大理巖組(Ar1wl2)綠泥石英巖、石英片巖及混合巖化巖,分布厚度5～7 m,節(jié)理較發(fā)育,巖體呈碎裂結(jié)構(gòu)。

二級邊坡根據(jù)當(dāng)時地質(zhì)編錄巖性為下太古界烏拉山群大理巖組(Ar1wl2)黑云母角閃斜長片麻巖,片麻理非常發(fā)育,節(jié)理發(fā)育,巖體呈碎裂結(jié)構(gòu),現(xiàn)已噴錨支護。

2.3 右岸高邊坡卸荷裂隙拉裂縫的分布和性質(zhì)特征

右壁發(fā)育有6條拉裂縫L2-1～L2-6。其中L2-1、L2-2位于樁號0+005～0+030之間一級、二級馬道坡頂,L2-1發(fā)育長度約25.0 m,裂縫寬度5～7 mm,L2-2發(fā)育長度約20.0 m,裂縫寬度5～7 mm,兩條拉裂縫延伸較長,均沿溢洪道軸向水平發(fā)育;L2-3、L2-4位于樁號0+065～0+090、0+105～0+110的二級馬道坡頂處,L2-3發(fā)育長度約25.0 m,裂縫寬度5～7 mm,L2-4發(fā)育長度約5.0 m,裂縫寬度5～7 mm,均沿溢洪道軸向水平發(fā)育;L2-5、L2-6位于樁號0+130～0+140間的二級馬道坡頂—三級馬道坡頂處,發(fā)育長度7.5～15.0 m,裂縫寬度3～8 mm,垂直溢洪道軸向順坡面縱向發(fā)育。

分析認為,L2-1為卸荷裂隙結(jié)構(gòu)面,結(jié)構(gòu)面與邊坡呈外傾關(guān)系,即結(jié)構(gòu)面傾向與邊坡坡向相同,結(jié)構(gòu)面傾角38°～43 °,小于邊坡坡度角(1∶0.75,即53°)。結(jié)構(gòu)面平整,傾角上陡下緩,結(jié)構(gòu)面延伸遠,控制范圍在樁號0+000～0+130。結(jié)構(gòu)面以上巖體已發(fā)生滑塌,滑塌范圍在樁號0+000～0+130?；^嚴重段,滑塌高度由一級馬道坡角至二級馬道坡頂,高程在1 464.00～1 486.50 m。L2-1卸荷裂隙結(jié)構(gòu)面,是造成右邊坡失穩(wěn)滑塌控制性的主結(jié)構(gòu)面。6條裂縫推測均為卸荷裂隙,其下部巖體為潛在不穩(wěn)定巖體,已建議進行邊坡加固處理。

3 FLAC3D數(shù)值模擬高邊坡穩(wěn)定性分析

3.1 FLAC3D概述

FLAC3D是在三維連續(xù)介質(zhì)中應(yīng)用的快速拉格朗日差分分析方法。它在分析巖土工程結(jié)構(gòu)的彈塑性力學(xué)行為、模擬施工過程等方面有其獨到的優(yōu)點。尤其在發(fā)生塑性流動或失穩(wěn)的情況下,此法可以很方便地用于模擬結(jié)構(gòu)從彈性到塑性屈服、失穩(wěn)破壞直至大變形的全過程,這是其它一些模擬計算方法無法比擬的[1]。

3.1.1 空間導(dǎo)數(shù)的有限差分近似

在FLAC3D中采用了混合離散法,區(qū)域被劃分為常應(yīng)變六面體單元的集合體,而在計算過程中,程序內(nèi)部又將每個六面體分為以六面體角點為角點的常應(yīng)變的集合體,變量均在四面體上進行計算,六面體單元的應(yīng)力、應(yīng)變?nèi)≈禐槠鋬?nèi)四面體的體積加權(quán)平均值[2]。

3.1.2 計算循環(huán)

FLAC3D的計算循環(huán)是FLAC3D程序?qū)嬎銌卧g的不平衡力,將此不平衡力重新添加到各節(jié)點上,再進行下一步的迭代運算,直到不平衡力足夠小或者各個節(jié)點的位移趨于平衡為止。

3.2 邊坡模型的建立

考慮到建立三維模型的復(fù)雜性,在建立三維模型時,對一些因素進行了簡化。同時,本文主要模擬邊坡開挖后周邊的應(yīng)力變化及位移對邊坡的影響,所以主要對邊坡區(qū)進行模擬計算。模型計算的基本參數(shù)根據(jù)研究區(qū)域的巖層柱狀圖和現(xiàn)場測得的巖層物理力學(xué)的實驗結(jié)果來確定。結(jié)合工程具體情況,建模時將工況分為:開挖前、開挖后、修坡三個階段予以考慮[3]。

此三維模型的建立采用分區(qū)組合的方法進行構(gòu)造,首先按層面和模型外邊界進行分區(qū),分別劃分每個區(qū)的網(wǎng)格,然后再將各分區(qū)組合在一起,形成計算模型。模型外圍規(guī)格1 600 m×400 m×630 m,共劃分出15 420個單元,16 544個節(jié)點(未開挖前的網(wǎng)格示意圖如圖1)。具體建模過程如下:按溢洪道右岸高邊坡工程地質(zhì)剖面圖將圖形剖分,以海拔高程1 400 m作為模型坐標中心,向水平方向和垂直方向建模。為減小計算誤差,滿足FLAC3D計算精度,假設(shè)模型沿邊坡走向長400 m,形成一個三維模型。其中X,Y,Z正方向依照右手定律,X以右方向為正,Z以垂直向上為正。

形成的初始狀態(tài)未開挖前的網(wǎng)格圖見圖1。

3.3 驗證模型合理性

模型建立完畢后,首先通過生成模型的最大不平衡力圖形確定邊坡天然條件下的穩(wěn)定性,驗算邊坡模型是否最終達到平衡。經(jīng)模擬計算,高邊坡最大不平衡力隨時間推移不斷減小,說明在天然條件下,模型經(jīng)過一段時間的應(yīng)力平衡過程后,各向應(yīng)力隨時間延續(xù)逐漸平衡,說明模型逐漸趨于穩(wěn)定,所建模型合理(見圖2)[4]。

3.4 開挖后邊坡穩(wěn)定性分析

高邊坡開挖后(見圖3),應(yīng)力分布發(fā)生改變,模型重新進行應(yīng)力計算并趨于平衡,最大不平衡力隨時間推移最終趨于定值(見圖4),說明開挖后模型達到平衡,應(yīng)力重新分布結(jié)束[5-6]。

圖1 邊坡開挖前的網(wǎng)格圖Fig.1 Grid pattern before slope excavation

圖2 開挖前最大不平衡力圖Fig.2 Maximum unbalanced graph before excavation

圖3 開挖后邊坡模型Fig.3 Slope model after excavation

開挖后,應(yīng)力分布改變,通過平衡后的模型最大主應(yīng)力和最小主應(yīng)力分布圖對比開挖前、后的不平衡力分布狀態(tài)。其中,由開挖前、后的最大、最小主應(yīng)力分布對比圖(見圖5、圖6),在高邊坡坡腳位置出現(xiàn)應(yīng)力集中,該位置可能出現(xiàn)較大規(guī)模位移,因此對該部分單元選取考察點觀測其位移變化。

圖4 開挖后最大不平衡力圖Fig.4 Maximum unbalanced graph after excavation

本次模擬計算中,設(shè)計在應(yīng)力集中區(qū)域內(nèi),布設(shè)位移考察點9個(見圖7),其中PK5、PK4、PK3、PK2、PK1點為邊坡開挖后考察點,P2、P3、P10、P15點為修坡后考察點。

通過各點開挖后應(yīng)力重分布狀態(tài)下的最大位移量計算,生成各考察點X、Y、Z軸向位移量分析圖,5個考察點PK5、PK4、PK3、PK2、PK1中,以PK2點位移最大,且隨時間推移該點位移最后趨于穩(wěn)定,即卸荷裂隙停止發(fā)展,說明開挖后邊坡處于穩(wěn)定狀態(tài),最大位移量為6.94 cm。詳見圖8開挖后點PK2(230,0,330)的X向位移。

3.5 修坡后邊坡穩(wěn)定性分析

修坡后邊坡穩(wěn)定性分析過程與開挖后邊坡穩(wěn)定性分析過程相同。首先根據(jù)軟件生成的最大不平衡力圖判斷模型是否穩(wěn)定,再通過分析軟件所生成的最大、最小主應(yīng)力圖分析邊坡是否出現(xiàn)應(yīng)力集中,最后通過各監(jiān)測點的位移圖判斷模型是否發(fā)生變形破壞。

圖5 最大主應(yīng)力分布對比圖Fig.5 Comparison chart of maximum principal stress

圖6 最小主應(yīng)力分布對比圖Fig.6 Comparison chart of minimum principal stress

圖7 考察點位置示意圖Fig.7 Sketch map of position of observation point

圖8 開挖后點PK2(230,0,330)的X向位移Fig.8 X directiomal displacement of PK2(230,0,330) after excavation

3.6 模擬計算成果分析

(1) 開挖前,模型中的最大不平衡力隨時間延續(xù)逐漸減小,最后趨近于零,這與天然條件下邊坡的穩(wěn)定情況相同,說明所建模型合理。

(2) 開挖后,應(yīng)力發(fā)生重分布,對比開挖前和開挖后模型最大主應(yīng)力和最小主應(yīng)力分布圖,在邊坡區(qū)部分位置出現(xiàn)應(yīng)力集中,說明該位置可能出現(xiàn)較大位移。所有觀察點在邊坡開挖后首先隨時間延續(xù)位移逐漸增大,但最終達到穩(wěn)定狀態(tài)。在所有觀察點中以PK2點的X方向位移最大,達到6.94 cm,但最終位移趨于一定值,不再出現(xiàn)更大位移,其他觀測點位移均較小。由此推之,邊坡整體在開挖后,只在PK1、PK2點附件發(fā)生較大位移,位移大小7 cm左右,其他位置位移很小,所以邊坡不會發(fā)生破壞,邊坡穩(wěn)定,無須加以支護措施[7]。

(3) 修坡后,分析其最大主應(yīng)力和最小主應(yīng)力分布圖可知,在邊坡位置除P2、P3、P10、P15點附近出現(xiàn)應(yīng)力集中外,其他位置應(yīng)力分布基本不發(fā)生改變。而P2、P3、P10、P15點周圍的應(yīng)力集中可能導(dǎo)致單元體發(fā)生大規(guī)模位移,使巖體破壞,邊坡失穩(wěn)。在所有觀察點中以P3點的X方向位移最大,達到10.84 cm,且所有觀察點均會發(fā)生較大規(guī)模位移,位移大小8～10 cm左右,且所有觀察點在邊坡開挖修坡后隨時間延續(xù)位移逐漸增大。由此說明,修坡后在邊坡修坡段會出現(xiàn)大規(guī)模位移,可能導(dǎo)致邊坡失穩(wěn),需加支護。加錨桿后,邊坡最大位移最大處為6.36 cm,比未加錨桿邊坡最大位移減小4 cm。

4 結(jié)論及高邊坡加固處理措施建議

(1) 優(yōu)化設(shè)計方案,建議在滿足水工泄洪順暢條件下,不再繼續(xù)向下開挖,不破壞邊坡坡腳,不產(chǎn)生臨空面,以保障現(xiàn)開挖邊坡相對穩(wěn)定條件。

(2) 對潛在不穩(wěn)定邊坡采用不同方式錨固處理。據(jù)FLAC3D模擬分析及水平孔揭露,建議于高邊坡坡腳考察點P2-P3范圍的應(yīng)力集中區(qū)及卸荷區(qū)域進行水平錨固,貫穿深度應(yīng)>22.0 m。

(3) 做好排水措施,于考察點P2水平向布設(shè)排水孔,減少節(jié)理裂隙內(nèi)的積水并降低再次滑塌的可能。

(4) 對考察點P2、P3、P10、P15等模擬位移量較大區(qū)域,做好邊坡巖體的長期變形觀測工作。

[1] 楊立強,張中杰,林舸,等.FLAC基本原理及其在地學(xué)中的應(yīng)用[J].地學(xué)前緣,2003,10(1):1-3.

[2] 龔紀文,崔建軍,席先武,等.FLAC數(shù)值模擬軟件及其在地學(xué)中的應(yīng)用[J].大地構(gòu)造與成礦學(xué),2002,21(3):321-325.

[3] 康紅普.回采巷道錨桿支護影響因素的FLAC分析[J].巖土力學(xué)與工程學(xué)報,1999,18(5):497-502.

[4] 胡斌,張倬元,黃潤秋,等.FLAC3D前處理程序的開發(fā)及仿真效果檢驗[J].巖石力學(xué)與工程學(xué)報,2002,21(9):1387-1391.

[5] Yang H Huang.土坡穩(wěn)定分析[M].包承綱,王清友,等譯.北京:清華大學(xué)出版社,1988:117-126.

[6] 殷宗澤,郭志平,徐鴻江,等.條分法土坡穩(wěn)定分析計算程序(SLP)[M]//姜弘逍.水工結(jié)構(gòu)工程與巖土工程的現(xiàn)代計算方法及程序.南京:河海大學(xué)出版社,1992:343-351.

[7] 張雄,陸明萬.巖土結(jié)構(gòu)穩(wěn)定性分析的塊體—夾層模型[J].計算力學(xué)學(xué)報,1997(S):397-401.

[8] 馮樹仁,豐定詳,葛修潤,等.邊坡穩(wěn)定性的三維極限平衡分析方法及應(yīng)用[J].巖土工程學(xué)報,1999,21(6):657-661.

(責(zé)任編輯：陳姣霞)

FLAC3D Numerical Stability Analysis of High Slope in Longsheng Reservoir

SHAN Wei, LI Jiandong

(WaterResourcesandHydropowerSurveyandDesignInstituteofInnerMongoliaAutonomousRegion,Hohhot,InnerMongolia010020)

Longsheng reservoir is medium controlling key water-control project,also is the flood water and preventing the sediment into the large key water conservancy project in Heihe.Because of llow-angle joint cracks,the development of migmatization,two large-scale collapse occur in the process of the high slope after excavation on the right bank slope of reservoir spillway. The paper has collected a large number of field survey data in the region,analyzed reservoir spillway slope stability by FLAC3D. Maximum and minimum principal stress diagram and a series of observation point displacement diagram by simulating calculation under hydro-geologic conditions of slope excavation and slope under the condition of the soil deformation and displacement,provide reference for the engineering design.

Longsheng reservoir; high slope; FLAC3D numerical simulation; stability analysis

2016-04-22;改回日期:2016-05-06

內(nèi)蒙古自治區(qū)烏蘭察布市卓資縣隆盛水庫工程地質(zhì)勘察(S090044105010)。

單威(1983-),男,工程師,水文與水資源工程專業(yè),從事水文地質(zhì)、水文水資源工作。E-mail:42392534@qq.com

TU457

A

1671-1211(2016)03-0315-05

10.16536/j.cnki.issn.1671-1211.2016.03.016

數(shù)字出版網(wǎng)址:http://www.cnki.net/kcms/detail/42.1736.X.20160505.1531.032.html 數(shù)字出版日期:2016-05-05 15:31