趙淵，王亮清，周鵬

(1.武漢市勘察設(shè)計(jì)有限公司，湖北 武漢 430022； 2.中國(guó)地質(zhì)大學(xué)(武漢)，湖北 武漢 430074；3.中交第二航務(wù)工程局有限公司技術(shù)中心，湖北 武漢 430040)

基于立體投影和3DEC的某水電站裂隙巖體邊坡穩(wěn)定性研究

趙淵1*，王亮清2，周鵬3

(1.武漢市勘察設(shè)計(jì)有限公司，湖北 武漢 430022； 2.中國(guó)地質(zhì)大學(xué)(武漢)，湖北 武漢 430074；3.中交第二航務(wù)工程局有限公司技術(shù)中心，湖北 武漢 430040)

裂隙巖體邊坡穩(wěn)定性研究是國(guó)內(nèi)外巖土工程領(lǐng)域研究的熱點(diǎn)問題，傳統(tǒng)方法在評(píng)價(jià)裂隙巖體邊坡穩(wěn)定性時(shí)具有諸多限制。本文以某水電站溢洪道開挖邊坡為例，提出了考慮確定性結(jié)構(gòu)面與隨機(jī)性結(jié)構(gòu)面，基于立體投影和3DEC的裂隙巖體邊坡穩(wěn)定性分析方法。在水電站工程地質(zhì)調(diào)研的基礎(chǔ)上，確定了優(yōu)勢(shì)結(jié)構(gòu)面，開展了結(jié)構(gòu)面力學(xué)參數(shù)的現(xiàn)場(chǎng)試驗(yàn)，采用立體投影方法進(jìn)行了設(shè)計(jì)邊坡角下邊坡可能滑動(dòng)模式分析；依據(jù)可能的滑動(dòng)模式，建立了三維非連續(xù)地質(zhì)模型與數(shù)學(xué)力學(xué)模型，進(jìn)行了自然、暴雨、地震三種工況下的位移與穩(wěn)定性分析。該成果為裂隙巖體邊坡穩(wěn)定性研究提供了一種新的高效的方法，對(duì)相似工程具有借鑒作用。

節(jié)理巖體；破壞模式；邊坡穩(wěn)定性；立體投影；非連續(xù)數(shù)值模擬

1 前 言

隨著經(jīng)濟(jì)的發(fā)展，大批水電工程、公路、鐵路工程在工程地質(zhì)條件的裂隙巖體山區(qū)修建，而工程建設(shè)中將會(huì)破壞原有自然斜坡的穩(wěn)定性，威脅著生命財(cái)產(chǎn)安全。因此如何高效準(zhǔn)確的評(píng)價(jià)由結(jié)構(gòu)面控制的裂隙巖體邊坡的穩(wěn)定性有著重要的意義[1，2]。結(jié)構(gòu)面在邊坡穩(wěn)定性研究中包括數(shù)量?jī)?yōu)勢(shì)或規(guī)模優(yōu)勢(shì)的確定性結(jié)構(gòu)面與規(guī)模較小或數(shù)量較少的隨機(jī)性結(jié)構(gòu)面，前者主要影響邊界條件或幾何模型，后者主要影響計(jì)算參數(shù)。

對(duì)于隨機(jī)裂隙巖體而言，國(guó)內(nèi)外常用的邊坡穩(wěn)定性研究方法有工程地質(zhì)類比法[3，4]、極限平衡法[5，6]、模型試驗(yàn)法[7，8]和數(shù)值模擬法[9，10]。工程地質(zhì)類比法雖能綜合考慮各影響邊坡穩(wěn)定的因素，但類比條件因地而異、經(jīng)驗(yàn)性強(qiáng)，一般有經(jīng)驗(yàn)的工程師才能給出較好的判斷。極限平衡法無法反映邊坡變形破壞的過程，需要首先確定假想的滑動(dòng)面，而假想的滑動(dòng)面局限性、隨機(jī)性很大[1]。模型試驗(yàn)法花費(fèi)較大，僅用于大型工程，很難準(zhǔn)確反映研究區(qū)所處的環(huán)境條件。在數(shù)值模擬方法方面，常用連續(xù)介質(zhì)和非連續(xù)介質(zhì)的方法分析。連續(xù)介質(zhì)分析大多考慮不同風(fēng)化程度的巖體邊坡的變形破壞過程，當(dāng)考慮結(jié)構(gòu)面系統(tǒng)對(duì)巖體邊坡變形破壞的影響時(shí)存在一些局限。而非連續(xù)介質(zhì)法，特別是離散元法，可反映復(fù)雜的巖體結(jié)構(gòu)，具有宏觀上的不連續(xù)性，不僅可以評(píng)價(jià)節(jié)理巖體邊坡穩(wěn)定性，還可較真實(shí)地刻畫出節(jié)理巖體邊坡的變形變化特征及破壞過程[11～14]。

離散元法由Peter Cundall(1971)首先提出[15，16]。目前常用的離散元軟件為ITASCA 開發(fā)的UDEC 及3DEC 軟件[17]。離散單元法是當(dāng)前研究邊坡穩(wěn)定性比較流行的非連續(xù)變形方法之一，國(guó)內(nèi)外一些學(xué)者從不同的角度采用離散單元法進(jìn)行了邊坡穩(wěn)定性研究。李世海等(2002)用3DEC模擬了三峽永久船閘高邊坡開挖后沿節(jié)理面的滑移過程。張沖等(2006)用3DEC對(duì)拱壩―壩肩進(jìn)行了三維可變形離散元整體穩(wěn)定分析。寧宇等(2007)利用3DEC軟件結(jié)合強(qiáng)度折減法對(duì)某水電站高邊坡進(jìn)行了穩(wěn)定性分析。孟國(guó)濤等(2007) 用3DEC對(duì)一典型的反傾向?qū)訝罱Y(jié)構(gòu)巖質(zhì)高邊坡進(jìn)行了模擬，分析了進(jìn)水口邊坡在開挖、加固作用下的變形狀況和穩(wěn)定性。鄭文棠等(2007)用3DEC進(jìn)行了復(fù)雜邊坡三維地質(zhì)可視化和數(shù)值模型構(gòu)建，探討加錨方法、三維合位移等值線云圖和塑性區(qū)可視化實(shí)現(xiàn)。巨能攀等(2009)基于塊體理論，用3DEC研究了邊坡的穩(wěn)定性。朱永生等(2011)基于3DEC分析了白鶴灘水電站復(fù)雜塊體穩(wěn)定性，并給出了加固建議[18～23]。M.Brideau等(2011)用3DEC對(duì)加拿大阿爾伯塔省亞省鴉巢峽南峰的邊坡的穩(wěn)定性進(jìn)行了分析，考慮了個(gè)別的不連續(xù)面、尺寸效應(yīng)和內(nèi)摩擦角的影響。G.Firpo等(2011) 通過3DEC法用地面數(shù)字測(cè)量法分析了巖石邊坡的穩(wěn)定性。M.Brideau等(2012) 用3DEC模型對(duì)平面平移下的邊坡破壞機(jī)制進(jìn)行了運(yùn)動(dòng)學(xué)的控制評(píng)價(jià)。K.S.Kalenchuk等(2013)考慮原位數(shù)據(jù)和三維數(shù)值模擬模型(3DEC)闡述了唐尼邊坡的地質(zhì)力學(xué)解釋。M.Curtaz等(2014) 用地面攝影測(cè)量法和3DEC數(shù)值模擬法對(duì)高山地區(qū)的巖石邊坡進(jìn)行了穩(wěn)定性分析[24～28]。

綜合以上研究成果：采用3DEC進(jìn)行邊坡穩(wěn)定性研究仍存在以下不足：①基于具體工程和工程地質(zhì)條件的差別，進(jìn)行分區(qū)分段，采用3DEC開展研究的成果較少；②在3DEC建模分析時(shí)，對(duì)確定性結(jié)構(gòu)面與隨機(jī)結(jié)構(gòu)面的分析不夠深入；③基于工程地質(zhì)分區(qū)分段，考慮確定性結(jié)構(gòu)面與隨機(jī)性結(jié)構(gòu)面，采用立體投影和3DEC進(jìn)行裂隙巖體邊坡穩(wěn)定性分析的成果目前尚未有報(bào)道。

基于此，本文以某水電站溢洪道開挖邊坡為例，提出了在工程地質(zhì)分區(qū)分段的基礎(chǔ)上，考慮確定性結(jié)構(gòu)面與隨機(jī)性結(jié)構(gòu)面，基于立體投影和3DEC法的裂隙巖體邊坡穩(wěn)定性研究方法。即對(duì)邊坡區(qū)優(yōu)勢(shì)結(jié)構(gòu)面及力學(xué)參數(shù)進(jìn)行研究，采用立體投影法進(jìn)行了邊坡可能滑動(dòng)模式的確定，建立了真實(shí)的三維地質(zhì)模型與數(shù)學(xué)力學(xué)模型，采用3DEC法進(jìn)行了自然、暴雨、地震工況下考慮不同滑動(dòng)模式邊坡的變形與穩(wěn)定性研究。研究成果為類似邊坡穩(wěn)定性研究提供了一種新的高效的方法。

2 地質(zhì)背景

2.1地質(zhì)條件

某水電站下壩址溢洪道開挖邊坡由弱風(fēng)化為主的英安巖巖體組成。根據(jù)設(shè)計(jì)資料可知，左坡坡長(zhǎng) 48 m，坡高 40 m，坡角53°，設(shè)計(jì)開挖兩級(jí)，開挖坡面呈階梯狀，坡面傾向?yàn)?8°。右坡較左坡大，坡長(zhǎng) 101.5 m，坡高 80 m，坡角55°，設(shè)計(jì)四級(jí)開挖，坡面傾向?yàn)?48°。溢洪道開挖邊坡工程地質(zhì)剖面圖如圖1所示。

圖1 溢洪道開挖邊坡工程地質(zhì)剖面圖

2.2工程分區(qū)及優(yōu)勢(shì)結(jié)構(gòu)面的確定

通過平硐編錄資料及現(xiàn)場(chǎng)精細(xì)測(cè)量的數(shù)據(jù)進(jìn)行結(jié)構(gòu)面統(tǒng)計(jì)分析，根據(jù)結(jié)構(gòu)面走向節(jié)理玫瑰花圖及節(jié)理極點(diǎn)等值線圖對(duì)結(jié)構(gòu)面進(jìn)行優(yōu)勢(shì)分組(見圖2)，確定了溢洪道開挖邊坡區(qū)優(yōu)勢(shì)結(jié)構(gòu)面的產(chǎn)狀，如表1所示。

溢洪道巖體優(yōu)勢(shì)結(jié)構(gòu)面及長(zhǎng)大裂隙產(chǎn)狀 表1

圖2 溢洪道巖體結(jié)構(gòu)面走向節(jié)理玫瑰花圖及裂隙極點(diǎn)等值線圖

在確定優(yōu)勢(shì)結(jié)構(gòu)面產(chǎn)狀后，區(qū)分結(jié)構(gòu)面屬于確定性結(jié)構(gòu)面還是隨機(jī)性結(jié)構(gòu)面。確定性結(jié)構(gòu)面指工程中的長(zhǎng)大裂隙、產(chǎn)狀類似數(shù)量眾多的優(yōu)勢(shì)結(jié)構(gòu)面，隨機(jī)性結(jié)構(gòu)面指工程中分布較少且不均的小型非貫通裂隙。對(duì)于工程巖體而言，邊坡區(qū)常見的結(jié)構(gòu)面分布情況為：①全為隨機(jī)性結(jié)構(gòu)面；②有一組為確定性結(jié)構(gòu)面，其余為隨機(jī)性結(jié)構(gòu)面；③有兩組或以上的結(jié)構(gòu)面為確定性結(jié)構(gòu)面，其余為隨機(jī)性結(jié)構(gòu)面。確定性與隨機(jī)性結(jié)構(gòu)面的不同組合情況如圖3所示，其中②可能沿確定性結(jié)構(gòu)面發(fā)生平面滑動(dòng)，③可能沿兩組確定性結(jié)構(gòu)面發(fā)生楔形體滑動(dòng)。

圖3不同結(jié)構(gòu)面的組合情況示意圖

結(jié)合現(xiàn)場(chǎng)地質(zhì)情況，綜合確定本區(qū)優(yōu)勢(shì)結(jié)構(gòu)面在之后穩(wěn)定性研究中按照確定性結(jié)構(gòu)面考慮。

根據(jù)工程地質(zhì)條件與工程特點(diǎn)進(jìn)行分區(qū)分段。鑒于實(shí)際工程的復(fù)雜性及規(guī)模性，要想系統(tǒng)明確的對(duì)其分析，有必要從地層巖性、工程方向、風(fēng)化程度、優(yōu)勢(shì)結(jié)構(gòu)面產(chǎn)狀等幾方面對(duì)工程區(qū)進(jìn)行分區(qū)分段討論。而對(duì)本工程而言，根據(jù)溢洪道邊坡區(qū)的巖性、溢洪道走向、風(fēng)化程度和結(jié)構(gòu)面產(chǎn)狀特點(diǎn)，可將溢洪道開挖邊坡分成一段進(jìn)行穩(wěn)定性評(píng)價(jià)，典型地質(zhì)剖面見圖1。

2.3結(jié)構(gòu)面與巖體力學(xué)性質(zhì)參數(shù)的確定

結(jié)構(gòu)面力學(xué)參數(shù)的確定，包括確定性結(jié)構(gòu)面、隨機(jī)性結(jié)構(gòu)面及假想結(jié)構(gòu)面。確定性結(jié)構(gòu)面力學(xué)參數(shù)需要根據(jù)室內(nèi)外試驗(yàn)資料進(jìn)行擬合，根據(jù)分布特點(diǎn)確定力學(xué)參數(shù)；隨機(jī)性結(jié)構(gòu)面力學(xué)參數(shù)可以根據(jù)現(xiàn)場(chǎng)試驗(yàn)、經(jīng)驗(yàn)估算或連續(xù)(有限單元法) 與非連續(xù)數(shù)值試驗(yàn)(離散單元法與顆粒流方法)確定表征單元體體積與相應(yīng)的力學(xué)參數(shù)。在利用3DEC進(jìn)行非貫通結(jié)構(gòu)面建模的過程中需要引入輔助面——假想結(jié)構(gòu)面，“假想結(jié)構(gòu)面”由Kulatilake等(1992)[29]，Wang和Kulatilake(1993)[30]提出，是指力學(xué)性質(zhì)與完整巖塊一致的結(jié)構(gòu)面。假想結(jié)構(gòu)面的參數(shù)可用實(shí)驗(yàn)法及估算法確定，Kulatilake等(1992)[29]給出了假想結(jié)構(gòu)面的建議取值方法，吳瓊[31]改進(jìn)了假想結(jié)構(gòu)面變形參數(shù)的估算公式。

由上，依據(jù)經(jīng)驗(yàn)類比與參數(shù)估算方法綜合確定了水電站裂隙巖體邊坡巖體的力學(xué)參數(shù)(如表2所示)。根據(jù)現(xiàn)場(chǎng)原位試驗(yàn)、室內(nèi)試驗(yàn)與經(jīng)驗(yàn)類比綜合確定了確定性結(jié)構(gòu)面的力學(xué)參數(shù)(如表3所示)。假想結(jié)構(gòu)面力學(xué)參數(shù)的選取賦值，主要依據(jù)Kulatilake等(1992)[29]提出的理論與經(jīng)驗(yàn)估算方法[31]進(jìn)行綜合取值(如表3所示)。

巖體力學(xué)參數(shù)取值表 表2

結(jié)構(gòu)面力學(xué)參數(shù)取值表 表3

3 基于立體投影的邊坡破壞模式確定

根據(jù)立體投影方法及工程設(shè)計(jì)的幾何參數(shù)確定發(fā)生平面滑動(dòng)與楔形體滑動(dòng)的可能性。對(duì)于以確定性結(jié)構(gòu)面為主的邊坡，發(fā)生平面滑動(dòng)和楔形體滑動(dòng)的可能性較大；對(duì)于以隨機(jī)性結(jié)構(gòu)面為主的邊坡，在穩(wěn)定性研究中將隨機(jī)性結(jié)構(gòu)面的影響反映在參數(shù)弱化上。

3.1左坡模式分析

根據(jù)表1優(yōu)勢(shì)結(jié)構(gòu)面統(tǒng)計(jì)資料與邊坡開挖面的設(shè)計(jì)邊坡傾向與傾角，采用立體投影法分析邊坡的破壞模式，本文采用平均坡角(53°)進(jìn)行分析，邊坡產(chǎn)狀如圖4(a)的線8所示。由圖4(a)分析可知，左幫邊坡發(fā)生平面滑動(dòng)破壞的可能性較小，但可能會(huì)沿J1、J5結(jié)構(gòu)面交線或J2、J4結(jié)構(gòu)面交線產(chǎn)生楔形體破壞。

3.2右坡模式分析

分析右坡破壞模式，平均坡角取55°，邊坡產(chǎn)狀如圖4(b)線8所示。由圖4(b)分析可知，右坡產(chǎn)生平面滑動(dòng)J6的可能性較大，沿J1、J6或J2、J6或J4、J6結(jié)構(gòu)面交線產(chǎn)生楔形體破壞的可能性也較大。

圖4 基于赤平投影的斜坡破壞模式分析

4 基于立體投影的裂隙巖體邊坡穩(wěn)定性研究

4.1模型的建立

為了進(jìn)一步驗(yàn)證上述邊坡的破壞模式，并對(duì)邊坡的穩(wěn)定性進(jìn)行定量分析，采用3DEC離散元軟件建立三維模型對(duì)邊坡的變形破壞進(jìn)行模擬。根據(jù)研究區(qū)的地質(zhì)條件，模型考慮了強(qiáng)風(fēng)化、弱風(fēng)化上帶、弱風(fēng)化下帶(水上)弱風(fēng)化下帶(水下)。模型主要考慮引起邊坡可能破壞模式的確定性結(jié)構(gòu)面，其他隨機(jī)性結(jié)構(gòu)面反映在巖體計(jì)算參數(shù)的弱化上。分別建立了考慮左坡J15(2楔形體破壞)、J24(4楔形)，右坡J16(1楔形)、J26(3楔形)、J46(5楔形)和J6(1平面破壞)破壞模式的3DEC模型。模型以最危險(xiǎn)的方式考慮，結(jié)構(gòu)面無限延伸且剪出口位于坡角處。模型計(jì)算范圍長(zhǎng) 445 m，寬 80 m，最高點(diǎn)高 399.4 m。計(jì)算采用莫爾- 庫(kù)侖屈服條件的彈塑性模型[24]。所建3DEC模型如圖5示(以1楔形為例)。

圖5 3DEC數(shù)值模型及結(jié)構(gòu)面示意圖(1楔形)

4.2計(jì)算結(jié)果分析

本文從穩(wěn)定性系數(shù)和位移兩方面分析了溢洪道左右邊坡在自然、暴雨、地震條件下的穩(wěn)定性。為了較精細(xì)研究邊坡穩(wěn)定性，以兩坡臺(tái)階處的點(diǎn)作為監(jiān)測(cè)點(diǎn)，如圖1中所示，具體坐標(biāo)為：對(duì)左坡2、4楔形體破壞取1(164.86，24，357.78)、2(181.01，45，337.99)、3(190.11，39.11，319.36) ，對(duì)右坡1、3、5楔形體破壞取4(343.53，40，390.75)、5(329.4，40，381.87)、6(312.9，40，360.47)、7(297.1，40，341.22)、8(285.4，40，317.95)，對(duì)1平面破壞取9(343.53，0，390.75)、10(329.4，0，381.87)、11(312.9，0，360.47)、12(297.1，0，341.22)、13(285.4，0，317.95)。

(1)穩(wěn)定性系數(shù)分析

利用3DEC軟件根據(jù)強(qiáng)度折減原理，直到邊坡恰好達(dá)到臨界破壞狀態(tài)，將此時(shí)的折減系數(shù)作為邊坡穩(wěn)定性系數(shù)。表4給出了溢洪道左右邊坡在6種破壞模式下不同工況的穩(wěn)定性系數(shù)，可以看出穩(wěn)定性系數(shù)按自然、暴雨、地震工況遞減，且地震條件對(duì)穩(wěn)定性系數(shù)的減少更為明顯，說明該處受地震影響較大。計(jì)算結(jié)果表明：①左坡(由2、4楔形體破壞控制)在三種工況下的穩(wěn)定性系數(shù)均大于1.2，處于穩(wěn)定狀態(tài)。②右坡(由1、3、5楔形體和1平面滑動(dòng)破壞控制)在自然工況下，3楔形體破壞模式的穩(wěn)定性系數(shù)大于1.2，處于穩(wěn)定狀態(tài)；1、5楔形體破壞的穩(wěn)定性系數(shù)在1～1.2，處于基本穩(wěn)定狀態(tài)，邊坡有發(fā)生滑動(dòng)的可能；1平面滑動(dòng)的穩(wěn)定性系數(shù)小于1，發(fā)生滑動(dòng)的可能性較大。③在暴雨工況下，3楔形體破壞模式的穩(wěn)定性系數(shù)大于1.2，處于穩(wěn)定狀態(tài)；1楔形體穩(wěn)定性系數(shù)在1～1.2，處于基本穩(wěn)定狀態(tài)；5楔形體破壞、1平面滑動(dòng)的穩(wěn)定性系數(shù)小于1，發(fā)生滑動(dòng)的可能性較大。④在地震工況下，3楔形體破壞模式的穩(wěn)定性系數(shù)在1～1.2，處于基本穩(wěn)定狀態(tài)，1、5楔形體破壞和1平面滑動(dòng)的穩(wěn)定性系數(shù)小于1，發(fā)生滑動(dòng)的可能性較大。

不同破壞模式下不同工況的穩(wěn)定性系數(shù) 表4

(2)位移分析

用3DEC軟件可較真實(shí)、動(dòng)態(tài)地模擬邊坡塊體的變形特征及破壞過程，如圖6為沿1楔形體破壞的右坡在地震工況下的x位移云圖，從圖中可清晰看出結(jié)構(gòu)面切割出的楔形體發(fā)生了整體下滑，并且坡角處位移最大，達(dá)到 15 m。圖7、8、9為各破壞模式下邊坡在三種工況下的監(jiān)測(cè)點(diǎn)的x位移(均以滑離原坡的方向?yàn)檎??？煽闯龈鞴r下不同破壞模式的位移變化趨勢(shì)基本相同，2、3、4楔形體破壞模式位移基本不變，在幾厘米左右，1、5楔形體及1平面滑動(dòng)破壞模式的位移變化較為明顯，并且都是在坡角處變化最為明顯。暴雨工況下各破壞模式的位移較自然工況有所增加，但均為幾十厘米，而地震條件對(duì)位移影響非常明顯，變形大的在 15 m～30 m左右，變化趨勢(shì)和穩(wěn)定性系數(shù)分析較吻合。

圖6 沿1楔形體破壞的邊坡在地震工況下的X位移圖

圖7 各破壞模式下邊坡在自然工況下的監(jiān)測(cè)點(diǎn)的x位移

圖8 各破壞模式下邊坡在暴雨工況下監(jiān)測(cè)點(diǎn)的x位移

圖9 各破壞模式下邊坡在自然工況下的監(jiān)測(cè)點(diǎn)的x位移

5 結(jié) 論

(1)本文采用立體投影和3DEC的方法，對(duì)某水電站溢洪道裂隙巖體邊坡的穩(wěn)定性進(jìn)行了研究，該方法可以全面地反映平面滑動(dòng)與楔形體滑動(dòng)對(duì)邊坡穩(wěn)定性的影響。

(2)本文從穩(wěn)定性系數(shù)與位移兩個(gè)方面對(duì)溢洪道左、右?guī)瓦吰略谧匀?、暴雨與地震工況下穩(wěn)定性進(jìn)行了研究，研究結(jié)果表明：①左幫邊坡穩(wěn)定性較右?guī)秃?；②暴雨、地震?duì)邊坡的穩(wěn)定性有削弱作用，地震的影響更為明顯。對(duì)各破壞模式而言，暴雨使穩(wěn)定性系數(shù)減少了0.01～0.2，地震使之減少了0.24～0.86，暴雨使位移增加了幾十厘米，地震使位移增加了幾十厘米到幾十米。其中地震使1楔形體破壞、5楔形體破壞和1平面滑動(dòng)破壞的位移達(dá)到 15 m～30 m。③對(duì)于右?guī)瓦吰?，三種工況下穩(wěn)定性系數(shù)為3楔形體破壞>1楔形體破壞>5楔形體破壞>1平面滑動(dòng)破壞，與位移分析相類似。

(3)本文在研究過程中為了保守起見，假設(shè)結(jié)構(gòu)面為無限延伸且通過坡腳的確定性結(jié)構(gòu)面。

致謝：本文相關(guān)資料是由中國(guó)水電顧問集團(tuán)貴陽勘測(cè)設(shè)計(jì)研究院提供，在此表示衷心感謝！

[1] 付華民. 白龍江苗家壩水電站茅埡里變形體穩(wěn)定性評(píng)價(jià)及其治理措施研究[D]．成都：成都理工大學(xué)碩士學(xué)位論文，2010.

[2] 夏元友，李梅. 邊坡穩(wěn)定性評(píng)價(jià)方法研究及發(fā)展趨勢(shì)[J]．巖石力學(xué)與工程學(xué)報(bào)，2002，21(7)：1087～1091.

[3] 張倬元，王士天，王蘭生. 工程地質(zhì)分析原理(第二版)[M]. 北京：地質(zhì)出版社，1994.

[4] 祁生文，伍法權(quán)，丁振明等. 從工程地質(zhì)類比的角度看錦屏一級(jí)水電站左岸深部裂縫的形成[J]. 巖石力學(xué)與工程學(xué)報(bào)，2004，23(8).

[5] Sarma，S.K..Stability analysis of embankment sand slopes[J]. Geoteeh. Engng. ASCE，105(12)：392～398，1979.

[6] 朱大勇，錢七虎，姜弘道等. 極限平衡條分法的統(tǒng)一計(jì)算格式[C]. 全國(guó)巖土與工程學(xué)術(shù)大會(huì)論文集(上冊(cè)) ，2003.

[7] 鄭穎人，趙尚毅，張魯渝. 用有限元強(qiáng)度折減法進(jìn)行邊坡穩(wěn)定分析[J]. 中國(guó)工程科學(xué)，2002，4(1)：57～61.

[8] 劉春玲，祁生文，童立強(qiáng)等. 利用FLAC3D分析某邊坡地震穩(wěn)定性[J]. 巖石力學(xué)與工程學(xué)報(bào)，2004，23(16).

[9] 胡修文，唐輝明，劉佑榮. 三峽庫(kù)區(qū)趙樹嶺滑坡穩(wěn)定性物理模擬試驗(yàn)研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2005，24(12)：2089～2095.

[10] 李邵軍，Knappett，J. A.，馮夏庭. 庫(kù)水位升降條件下邊坡失穩(wěn)離心模型試驗(yàn)研究[J]. 巖石力學(xué)與工程學(xué)報(bào)，2008，27(8)：1586～1593.

[11] 寧宇，徐衛(wèi)亞，鄭文棠. 應(yīng)用離散元強(qiáng)度折減對(duì)復(fù)雜邊坡進(jìn)行穩(wěn)定性分析[J]. 巖土力學(xué)，2007(S1)：569～574.

[12] 邢萬波，周鐘. 錦屏一級(jí)水電站左岸壩肩邊坡的3DEC變形和穩(wěn)定性分析與認(rèn)識(shí)[J]. 水電站設(shè)計(jì)，2010，26(1).

[13] Itasca Consulting Group，Inc. 3 Dimensional Distinct Element Code User′s Guide[R]. 2003.

[14] Lmos，J.V.. A distinct element model for dynamic analysis ofjointed rock with application to dam foundation and fault motion. Minneapolis：University of Minnesota，1987.

[15] Cundall，P.A.. A computer model for simulating progressive. Large scale movements in blocky rock systems[J]. Proceedings of the International Symposium Rock Fracture.ISRM，Nancy，1971(1).

[16] Cundall，P.A.，Strack，O.D.L.. A discrete numerical model for granular assembles[J]. Geotechnique，1979，29(1)：47～65.

[17] Itasca Consulting Group，Inc，Itasca. 3DEC，theory and background，three-dimensional distinct element code，version 3.0，user’s manual[R]. USA：Itasca Consulting Group，Inc.

[18] 李世海，高波，燕琳. 三峽永久船閘高邊坡開挖三維離散元數(shù)值模擬[J]. 巖土力學(xué)，2002，23(3)：272～277.

[19] 張沖，侯艷麗，金峰等. 拱壩―壩肩三維可變形離散元整體穩(wěn)定分析[J]. 巖石力學(xué)與工程學(xué)報(bào)，2006，25(6)：1226～1 232.

[20] 孟國(guó)濤，徐衛(wèi)亞，鄭文棠等. 巖質(zhì)高邊坡開挖變形的三維離散單元法分析[J]. 河海大學(xué)學(xué)報(bào)，2007，35(4)：393～397.

[21] 鄭文棠，徐衛(wèi)亞，童富果等. 復(fù)雜邊坡三維地質(zhì)可視化和數(shù)值模型構(gòu)建[J]. 巖石力學(xué)與工程學(xué)報(bào)，2007，26(8)：1634～1644.

[22] 巨能攀，趙建軍，黃潤(rùn)秋等. 基于3DEC 的邊坡塊體穩(wěn)定性分析[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2009，28(6)：925～928.

[23] 朱永生，朱煥春，石安池等. 基于離散單元法的白鶴灘水電站復(fù)雜塊體穩(wěn)定性分析[J]. 巖石力學(xué)與工程學(xué)報(bào)，2011，30(10)：2069～2075.

[24] Brideau，M.，et al. Three-dimensional slope stability analysis of South Peak，Crowsnest Pass，Alberta，Canada [J]. Landslides，2011，8(2)：139～158.

[25] Firpoa，G.，Salvinia，R.，et al. Use of Digital Terrestrial Photogrammetry in rocky slope stability analysis by Distinct Elements Numerical Methods [J]. International Journal of Rock Mechanics and Mining Sciences，48(7)：1045～1054.

[26] Brideau，M.，et al. Evaluating Kinematic Controls on Planar Translational Slope Failure Mechanisms Using Three-Dimensional Distinct Element Modeling [J]. Geotechnical and Geological Engineering，2012，30(4)：991～1011.

[27] Kalenchuk，K. S.，Hutchinson，D. J.，Diederichs，M. S.. Geomechanical interpretation of the Downie Slide considering field data and three-dimensional numerical modeling[J]. Landslides，2013，10(6)：737～756.

[28] Curtaz，M.，et al. Terrestrial Photogrammetry and Numerical Modelling for the Stability Analysis of Rock Slopes in High Mountain Areas：Aiguilles Marbrées case [J]. Rock Mechanics and Rock Engineering，2014，47(2)：605～620.

[29] Kulatilake，P.H.S.W.，Ucpirti，H.，Wang，S，et al. Use of the distinct element method to perform stress analysis in rock with non-persistent joints and to study the effect of joint geometry parameters on the strength and deformability of rock masses[J]. Rock Mechanics and Rock Engineering，1992，25(4)：253～274.

[30] Wang，S.，Kulatilake，P.H.S.W.. Linking between joint geometry models and a distinct element method in three dimensions to perform stress analyses in rock masses containing finite size joints[J]. Japanese Society of Soil Mechanics and Foundation Engineering，1993，33(4)：88～98.

[31] 吳瓊. 復(fù)雜節(jié)理巖體力學(xué)參數(shù)尺寸效應(yīng)及工程應(yīng)用研究[D]. 武漢：中國(guó)地質(zhì)大學(xué)(武漢)博士學(xué)位論文，2012.

TheStudyfortheJointedRockMassSlopeStabilityofAHydropowerStationBasedontheStereographicProjectionand3DECMethod

Zhao Yuan1，Wang Liangqing2，Zhou Peng3

(1.Wuhan Geotechnical Engineering and Surveying Co.，Ltd.，Wuhan 430022，China； 2.Faculty of Engineering，China University of Geosciences，Wuhan 430074，China； 3.China Communications Construction Company Second Harbour Engineering Co.，Ltd National Enterprise Technology Center，Wuhan 430040，China)

The study for slope stability is the hot issue in geotechnical engineering field at home and abroad，and the traditional method for evaluating the jointed rock mass slope stability has many limits. The spillway excavation slope of a hydropower station in Tibet was taken as a case study，and the study method for slope stability of jointed rock mass based on stereographic projection and 3DEC method，with taking deterministic and random discontinuities into account，was put forward in this paper. The preferred discontinuities were determined on the basis of engineering geological research，the in-situ tests of mechanical parameters of discontinuities were carried out，and the possible sliding modes of slope with the designed slope angle were analyzed by using stereographic projection method. The 3D discontinuous geologic and mathematical mechanical model was built to analyze displacement and stability in nature，heavy storm and earthquake conditions. The achievements provide a new efficient way for jointed rock mass slope stability study，and give guidance for similar engineering.

jointed rock mass；failure mode；slope stability；stereographic projection；discontinuous numerical simulation

1672-8262(2017)05-153-07

P642

A

2017—07—14

趙淵(1989—)，女，碩士，主要從事巖體工程性質(zhì)與穩(wěn)定性方面的研究。

王亮清(1972—)，男，教授，博士生導(dǎo)師，主要從事巖土體穩(wěn)定性與地質(zhì)災(zāi)害防治的教學(xué)與科學(xué)工作。

國(guó)家自然科學(xué)基金資助(41372310)；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助(CUG090104)。