陳曦 喻文振 彭春林 熊新宇 何為 羅立哲

摘要：金沙江旭龍水電站徐龍石料場(chǎng)開采高度高達(dá)280 m，距大壩等重要水工建筑物較近，為分析在不同開挖步下，該料場(chǎng)高陡巖質(zhì)邊坡的變形及穩(wěn)定性特征，建立了料場(chǎng)分級(jí)開挖有限元模型。結(jié)果表明：① 徐龍料場(chǎng)開采過(guò)程中邊坡變形較大區(qū)域由原坡腳向開挖新形成的坡腳不斷變化；② 隨著開挖的進(jìn)行，徐龍石料場(chǎng)最大主應(yīng)變較大區(qū)域發(fā)生“分岔”，最不利滑動(dòng)圓弧由原邊坡坡腳逐漸移動(dòng)至開挖新形成的坡腳處；③ 不同開挖步下的徐龍石料場(chǎng)整體穩(wěn)定性滿足要求。研究成果可以為類似的石料場(chǎng)邊坡分級(jí)開挖穩(wěn)定性分析提供參考。

關(guān)鍵詞：石料場(chǎng)； 高陡巖質(zhì)邊坡； 分級(jí)開挖； 旭龍水電站

中圖法分類號(hào)：TV223

文獻(xiàn)標(biāo)志碼：A

DOI：10.15974/j.cnki.slsdkb.2024.05.012

文章編號(hào)：1006-0081（2024）05-0066-06

0 引言

料場(chǎng)是水電站大壩澆筑的重要骨料來(lái)源。為提高開采利用率、減小用地面積，水電工程料場(chǎng)邊坡通常較為高陡［1-2］。目前國(guó)內(nèi)關(guān)于水電工程料場(chǎng)的研究多關(guān)注料場(chǎng)的開采規(guī)劃，例如石清文［1］根據(jù)優(yōu)先利用回采料的原則，分析了枕頭壩二級(jí)水電站料源規(guī)劃；尹岳降等［2］以白鶴灘水電站為例，提出了基于分形理論的人工骨料比表面積計(jì)算方法；宋寅［3］采用有限元方法分析了猴子巖水電站色龍溝料場(chǎng)邊坡穩(wěn)定性；徐敬武等［4］基于大崗山水電站，研究了棱子壩人工骨料場(chǎng)儲(chǔ)量情況；沈貴基等［5］研究了爆破對(duì)錦屏一級(jí)水電站大奔流溝料場(chǎng)開采的影響；劉大顯等［6］分析了大奔流溝料場(chǎng)高邊坡錨索錨固力增大的原因。但是，現(xiàn)有的研究更注重實(shí)際工程的應(yīng)用，而較少考慮分級(jí)開挖過(guò)程中的受力與變形特征及其對(duì)料場(chǎng)邊坡穩(wěn)定性的影響［7-10］。

金沙江上游旭龍水電站徐龍石料場(chǎng)為典型花崗巖料場(chǎng)，邊坡開采高度達(dá)280 m，為典型高陡巖質(zhì)邊坡。不同于一般的巖質(zhì)邊坡，徐龍石料場(chǎng)邊坡具有一定的特殊性：① 受開采、征地范圍限制影響，料場(chǎng)較為高陡；② 料場(chǎng)邊坡整體穩(wěn)定性受開挖過(guò)程影響，在不同的開挖步，邊坡相應(yīng)的安全系數(shù)不同并隨開挖過(guò)程不斷發(fā)生變化。徐龍料場(chǎng)邊坡一旦發(fā)生失穩(wěn)破壞，極有可能造成嚴(yán)重的生命財(cái)產(chǎn)、經(jīng)濟(jì)損失。本文以徐龍石料場(chǎng)為例，建立高陡巖質(zhì)邊坡有限元模型，分析開挖過(guò)程中料場(chǎng)邊坡變形及穩(wěn)定性特征，為料場(chǎng)邊坡安全施工提供科學(xué)依據(jù)。

1 工程概況及地質(zhì)條件

徐龍石料場(chǎng)（圖1）位于金沙江旭龍水電站壩址區(qū)上游約700 m處的徐龍溝附近，料場(chǎng)巖性以花崗巖為主，料場(chǎng)開采區(qū)投影面積8.51萬(wàn)m2，開采總量549萬(wàn)m3，剝離量117.25萬(wàn)m3，花崗巖有用料開采量431.75萬(wàn)m3，主要為旭龍水電站混凝土拱壩澆筑提供骨料來(lái)源。料場(chǎng)邊坡開口線頂部高程約為2 510 m，設(shè)計(jì)終采高程為2 230 m，邊坡開采高度達(dá)280 m，遠(yuǎn)高于GB 50330-2013《建筑邊坡工程技術(shù)規(guī)范》適用的巖質(zhì)邊坡高度上限（30 m），為典型高陡巖質(zhì)邊坡。徐龍石料場(chǎng)毗鄰大壩上游2號(hào)索道橋（溝通左右岸的重要交通樞紐），下距高、低導(dǎo)流洞入口、截流圍堰約400 m。料場(chǎng)表面覆蓋層厚度約5～15 m，由于料場(chǎng)與截流備料場(chǎng)分居金沙江左右岸，一旦產(chǎn)生崩塌有可能形成高危堰塞湖，從而威脅下游圍堰、大壩等水工建筑物安全。根據(jù)SL 386-2007《水利水電工程邊坡設(shè)計(jì)規(guī)范》，料場(chǎng)邊坡破壞對(duì)周圍重要水工建筑物影響較大，邊坡級(jí)別確定為1級(jí)。

徐龍石料場(chǎng)為順江發(fā)育的單面山體，屬高中山地貌（圖2）。金沙江流向SE140°，枯季江水位2 160 m左右；高程2 170～2 190 m為江邊一級(jí)階地，寬度20～80 m；高程2 190～2 400 m為陡崖，坡度65°～80°，其上為地形坡度35°左右的斜坡，局部為高20 m左右的陡坎，直至山頂，山頂高程3 900 m左右。

料場(chǎng)內(nèi)地層由第四系沖洪積物（Qpal）和三疊紀(jì)印支期侵入的花崗巖（γ15）組成。第四系沖洪積物（Qpal）分布在山坡腳與江邊之間的緩坡上，為一級(jí)階地，其物質(zhì)成分為含碎塊石、漂石、粉砂質(zhì)黏土夾花崗巖碎屑，厚度一般為20～30 m?；◢弾r（γ15）呈灰白色或淺灰色，主要礦物為角閃石、長(zhǎng)石和石英，多為中細(xì)粒、少量中粗粒和細(xì)粒結(jié)構(gòu)，塊狀構(gòu)造，巖質(zhì)堅(jiān)硬?；◢弾r內(nèi)夾有少量灰黑色混合巖，呈條帶狀、角礫狀、腸狀等不規(guī)則狀分布。根據(jù)巖礦試驗(yàn)成果，區(qū)內(nèi)花崗巖根據(jù)巖石礦物成分又可進(jìn)一步分為蝕變石英閃長(zhǎng)巖、花崗閃長(zhǎng)巖、斜長(zhǎng)角閃巖、石英閃長(zhǎng)巖、斜長(zhǎng)花崗巖、角閃片麻巖以及黑云母花崗巖等。地表花崗巖一般呈微風(fēng)化狀，偶見(jiàn)有弱風(fēng)化帶，弱風(fēng)化帶不連續(xù)，厚度為2.0～10.0 m。

2 有限元模型

徐龍石料場(chǎng)開挖共形成13級(jí)馬道，相鄰兩級(jí)馬道高差為15 m，高程2 305 m以上單級(jí)邊坡一般坡比為1∶0.25；坡頂層卸荷帶邊坡坡比為1∶0.5；高程2 305 m以下單級(jí)邊坡坡比為1∶0.3?？紤]邊坡支護(hù)施工便利，大部分馬道寬為3 m；高程2 410 m的馬道加寬至10 m，并設(shè)置一道被動(dòng)防護(hù)網(wǎng)；高程2 305 m馬道加寬至8 m，作為料場(chǎng)后期安全巡視平臺(tái)。開挖邊坡的支護(hù)在分層開挖過(guò)程中逐層進(jìn)行，上層支護(hù)保證下層開挖安全。根據(jù)上述要求，本文建立徐龍石料場(chǎng)邊坡有限元模型如圖3（a）所示。由于開挖過(guò)程中馬道從上至下逐級(jí)形成，故將料場(chǎng)開采分為13級(jí)開挖體。由于開挖體體型較不規(guī)則，采用有限元軟件Abaqus四面體一階線性單元對(duì)邊坡進(jìn)行網(wǎng)格劃分，并對(duì)開挖體附近進(jìn)行局部網(wǎng)格加密，得到共計(jì)約12萬(wàn)個(gè)四面體單元，見(jiàn)圖3（b）。

強(qiáng)度折減法是分析邊坡穩(wěn)定性的常用手段。強(qiáng)度折減法認(rèn)為影響邊坡穩(wěn)定性的參數(shù)主要是黏聚力c和內(nèi)摩擦角φ，通過(guò)不斷調(diào)整折減系數(shù)反復(fù)計(jì)算，直至其達(dá)到臨界破壞，得到的折減系數(shù)即為穩(wěn)定性系數(shù)Fs。假設(shè)此時(shí)黏聚力和內(nèi)摩擦角分別為ccr和φcr，由于邊坡處于臨界狀態(tài)，所對(duì)應(yīng)的安全系數(shù)Kcr=1，則原始邊坡對(duì)應(yīng)的穩(wěn)定性系數(shù)Fs可由式（3）求得。

式中：φ和φ′分別為折減前后的內(nèi)摩擦角；c和c′分別為折減前后的黏聚力；Fi為i折減時(shí)步對(duì)應(yīng)的折減系數(shù)；K為邊坡安全系數(shù)；c0和φ0分別為初始黏聚力和初始內(nèi)摩擦角。

有限差分軟件FLAC3D采用強(qiáng)度折減法計(jì)算邊坡穩(wěn)定安全系數(shù)。將建立的開挖體模型和劃分的網(wǎng)格單元導(dǎo)入有限差分軟件FLAC3D，進(jìn)行邊坡穩(wěn)定性分析。錨桿采用結(jié)構(gòu)性錨索（Cable）單元進(jìn)行模擬，上一級(jí)邊坡開挖完成后施加邊坡錨桿支護(hù)，按照“一挖一支護(hù)”依次進(jìn)行。料場(chǎng)巖土體采用Mohr-Coulomb本構(gòu)模型，料場(chǎng)邊坡邊界條件為：離坡面較遠(yuǎn)的后緣側(cè)和底部垂直于邊界法向方向的位移為0。根據(jù)前期地質(zhì)勘探成果，該料場(chǎng)以單一花崗巖為主，巖土物理力學(xué)參數(shù)如表1所示。

邊坡開挖穩(wěn)定性計(jì)算前，首先進(jìn)行地應(yīng)力平衡。初始地應(yīng)力平衡結(jié)果如圖4所示。開挖過(guò)程采用生死單元進(jìn)行模擬，每開挖一級(jí)，對(duì)開挖體進(jìn)行“挖除”操作。

3 模擬結(jié)果分析

3.1 料場(chǎng)邊坡分段開挖位移演化特征

圖5為不同開挖步對(duì)應(yīng)的徐龍石料場(chǎng)邊坡分級(jí)開挖及支護(hù)位移云圖，紅色為邊坡變形較大區(qū)域，顏色越深，邊坡位移越大。前9個(gè)開挖步，邊坡位移較大區(qū)域靠近坡腳，和臨空面相交（對(duì)應(yīng)的料場(chǎng)開采高程為2 450～2 320 m）；第10個(gè)開挖步（對(duì)應(yīng)的開采高程為2 320 m）后，各開挖步邊坡位移較大區(qū)域穿過(guò)開挖形成的新的坡腳。因此，在料場(chǎng)開采過(guò)程中，邊坡變形特征逐漸發(fā)生變化，變形較大區(qū)域由原坡腳逐漸轉(zhuǎn)至開挖新形成的坡腳。

圖6為開挖區(qū)位移云圖（已扣除地應(yīng)力平衡產(chǎn)生的位移）。隨著開挖的進(jìn)行，開挖擾動(dòng)區(qū)的范圍不斷增大。前7個(gè)開挖步，開挖區(qū)下部變形擾動(dòng)范圍較上部變形擾動(dòng)區(qū)范圍大；后5個(gè)開挖步，開挖區(qū)上部土體變形擾動(dòng)范圍較下部擾動(dòng)范圍大。開挖形成的平臺(tái)處擾動(dòng)位移量最大。

3.2 料場(chǎng)邊坡分段開挖及支護(hù)最大主應(yīng)變演化特征

圖7為不同開挖步料場(chǎng)邊坡最大主應(yīng)變?cè)茍D。料場(chǎng)邊坡最大主應(yīng)變區(qū)域可以近似視為邊坡滑動(dòng)最不利圓弧滑動(dòng)位置。與圖6所反映的結(jié)果類似，前9個(gè)開挖步最大主應(yīng)變較大區(qū)域通過(guò)料場(chǎng)邊坡坡腳，這表明最不利滑動(dòng)圓弧也通過(guò)邊坡坡腳，和邊坡的典型破壞特征相吻合。隨著開挖的進(jìn)行，第10個(gè)開挖步和第11個(gè)開挖步最大主應(yīng)變較大區(qū)域發(fā)生“分岔”現(xiàn)象，其最不利滑動(dòng)圓弧分別通過(guò)原邊坡坡腳和開挖新形成的坡腳處，表明邊坡整體穩(wěn)定性由新開挖形成的坡腳控制。

3.3 料場(chǎng)邊坡分段開挖穩(wěn)定性系數(shù)

圖8為采用有限元強(qiáng)度折減法計(jì)算的第2個(gè)開挖步到第13個(gè)開挖步邊坡穩(wěn)定性系數(shù)FOS（考慮到第一個(gè)開挖步開挖量較小故暫不展示）。邊坡抗滑穩(wěn)定性系數(shù)變化趨勢(shì)見(jiàn)圖9。在料場(chǎng)開采過(guò)程中，邊坡穩(wěn)定性系數(shù)并非隨著開挖的進(jìn)行呈單調(diào)增大或減小，而是有一定波動(dòng)。前5個(gè)開挖步，邊坡開挖量相對(duì)較小，穩(wěn)定性系數(shù)在2.16左右，變化較??；第6～8個(gè)開挖步，邊坡穩(wěn)定性系數(shù)有一定增加（最大為2.19），這可能是由于料場(chǎng)前期開挖造成下滑力減小量大于抗滑力減小量；第9個(gè)開挖步到第13個(gè)開挖步，料場(chǎng)邊坡穩(wěn)定性系數(shù)由2.17逐漸降低為1.91，表明開挖引起的下滑力減小量逐漸小于抗滑力減小量。因此，應(yīng)當(dāng)重點(diǎn)對(duì)后4個(gè)開挖步進(jìn)行邊坡變形監(jiān)測(cè)?？傮w來(lái)看，開挖過(guò)程中邊坡穩(wěn)定性系數(shù)大于規(guī)范所要求的臨界值，邊坡整體安全。

4 結(jié)論

本文采用有限差分軟件FLAC3D計(jì)算了不同開挖階段金沙江上游旭龍水電站徐龍石料場(chǎng)邊坡變形及穩(wěn)定性特征，主要結(jié)論如下。

（1） 徐龍石料場(chǎng)開挖過(guò)程中，最不利滑動(dòng)圓弧由原坡腳逐漸向新生成的坡腳移動(dòng)，邊坡整體穩(wěn)定性由原坡腳控制逐漸變?yōu)橛尚麻_挖形成的坡腳控制，開挖擾動(dòng)范圍不斷增大。

（2） 徐龍石料場(chǎng)開采過(guò)程中穩(wěn)定性系數(shù)滿足規(guī)范要求，應(yīng)當(dāng)重點(diǎn)對(duì)后4個(gè)開挖步進(jìn)行邊坡變形監(jiān)測(cè)。料場(chǎng)邊坡發(fā)生整體滑移失穩(wěn)可能性較小，整體穩(wěn)定性較好，對(duì)周圍臨建水工建筑物影響也較小。

（3） 徐龍石料場(chǎng)邊坡開挖過(guò)程中邊坡變形特征不斷變化，邊坡下挖時(shí)應(yīng)加強(qiáng)監(jiān)測(cè)預(yù)警，并及時(shí)與數(shù)值分析結(jié)果進(jìn)行比較。

參考文獻(xiàn)：

［1］石清文.枕頭壩二級(jí)水電站料源規(guī)劃分析［J］.人民長(zhǎng)江，2022，53（增2）：89-91.

［2］尹岳降，朱子晗，盧文波，等.基于分形理論的人工骨料比表面積計(jì)算方法——以白鶴灘水電站為例［J］.人民長(zhǎng)江，2020，51（4）：196-201.

［3］宋寅.猴子巖水電站色龍溝料場(chǎng)邊坡穩(wěn)定性分析［J］.人民長(zhǎng)江，2015，46（增1）：145-147.

［4］徐敬武，鄧忠文，吳灌州.大崗山水電站棱子壩人工骨料場(chǎng)勘察［J］.人民長(zhǎng)江，2014，45（22）：40-42.

［5］沈貴基，彭義，孫明瑞.錦屏一級(jí)水電站大奔流溝料場(chǎng)開采爆破技術(shù)［J］.人民長(zhǎng)江，2013，44（14）：66-68，83.

［6］劉大顯，孫云志，陜碩.大奔流溝料場(chǎng)高邊坡錨索錨固力增大原因分析［J］.人民長(zhǎng)江，2013，44（14）：18-21.

［7］李勤軍，鄢雙紅.大奔流溝料場(chǎng)高邊坡支護(hù)設(shè)計(jì)研究與實(shí)踐［J］.人民長(zhǎng)江，2013，44（14）：22-25.

［8］何楊，張志強(qiáng)，向勇.李家?guī)r水庫(kù)導(dǎo)流泄洪放空洞進(jìn)水口邊坡支護(hù)方案三維有限元模擬［J］.水利水電快報(bào)，2023，44（12）：48-53.

［9］陳紅如，李堅(jiān)，李冬冬，等.滇中引水石鼓水源工程強(qiáng)風(fēng)化巖質(zhì)邊坡穩(wěn)定分析與治理［J］.水利水電快報(bào)，2023，44（11）：64-72.

［10］姚行.楊房溝水電站大壩邊坡壓力分散型預(yù)應(yīng)力錨索張拉驗(yàn)收試驗(yàn)研究［J］.水利水電快報(bào)，2023，44（增1）：47-50.

編輯：高小雲(yún)

Simulation on staged excavation of high and steep quarry slope of Xulong Hydropower Station

CHEN Xi1，YU Wenzhen1，PENG Chunlin2，XIONG Xinyu1，HE Wei1，LUO Lizhe1

（1.Changjiang Survey，Planning，Design and Research Co.，Ltd.，Wuhan 430010，China；2.China South-to-North Water Diversion Jianghan Water Network Construction and Development Co.，Ltd.，Wuhan 430040，China）

Abstract：

The mining height of Xulong quarry of Xulong Hydropower Station in the Jinsha River is up to 280 m，which is close to the dam and other important hydraulic structures.A finite element model of staged excavation of the quarry slope was established to analyze the deformation and stability characteristics of the high and steep slope under different excavation steps.The results showed that：① In the mining process of Xulong quarry，the large deformation area of the slope gradually changes from the original slope foot to the newly formed slope foot.② With the progress of excavation，bifurcation occured in the large maximum principal strain area of Xulong quarry，and the most unfavorable sliding arc gradually moved from the original slope foot to the newly formed slope foot.③ The overall stability of Xulong quarry met the requirements under different excavation steps.The research results can provide a reference for stability analysis of staged excavation of slope for similar rock quarry.

Key words：

quarry； high and steep rock slope； staged excavation； Xulong Hydropower Station

收稿日期：2023-10-13

基金項(xiàng)目：長(zhǎng)江設(shè)計(jì)集團(tuán)自主科研項(xiàng)目（CX2017Z32）

作者簡(jiǎn)介：陳曦，男，工程師，博士，主要從事水利水電工程施工組織設(shè)計(jì)及研究工作。E-mail：chenxi5@cjwsjy.com.cn

通信作者：喻文振，男，高級(jí)工程師，碩士，主要從事水利水電工程施工組織設(shè)計(jì)及研究工作。E-mail：yuwenzhen@cjwsjy.com.cn