黃 鋒，于孫相，謝秉忠，王周望，廖言豐，馬希磊，王冬明，查文華

（1.福建建工基礎(chǔ)設(shè)施建設(shè)集團(tuán)有限責(zé)任公司，福建 福州 350001；2.福建建工集團(tuán)有限責(zé)任公司，福建 福州 350001；3.華設(shè)設(shè)計(jì)集團(tuán)股份有限公司，江蘇 南京 210000；4.東華理工大學(xué)土木與建筑工程學(xué)院，江西 南昌 330013）

隨著我國(guó)基礎(chǔ)建設(shè)的不斷推進(jìn)，越來(lái)越多的山區(qū)建設(shè)項(xiàng)目日益增多。進(jìn)而人工開挖導(dǎo)致不同破碎程度的山體陡坡頻頻出現(xiàn)，加之降雨、地震和人工擾動(dòng)等外部因素的影響，此類人工削坡揭露的結(jié)構(gòu)面更容易導(dǎo)致滑坡、崩塌和剝落等地質(zhì)災(zāi)害的發(fā)生，嚴(yán)重威脅著公路的安全生產(chǎn)建設(shè)。因此對(duì)高陡破碎巖質(zhì)邊坡進(jìn)行有針對(duì)性穩(wěn)定性預(yù)測(cè)分析非常重要（樊赟赟等，2022）。陳慶峰等（2021）建立滑移-彎曲式傾倒變形力學(xué)模型，采用Rhino 6.0—FLAC3D6.0高精度耦合建模方法對(duì)斷層破碎帶區(qū)域進(jìn)行數(shù)值模擬穩(wěn)定性分析，得出隨著采場(chǎng)開拓逐漸加深，邊坡高度亦不斷增加，采場(chǎng)壩幫受斷層破碎帶的影響，區(qū)域局部失穩(wěn)變?yōu)檎w失穩(wěn)；劉永波等（2016）針對(duì)長(zhǎng)河壩水電站進(jìn)水口邊坡巖體破碎的地質(zhì)構(gòu)造，通過(guò)錨固鉆孔發(fā)現(xiàn)局部邊坡受卸荷及地質(zhì)構(gòu)造影響，坡體產(chǎn)生滑移拉裂破壞，表部形成直徑0.5～5.0 m大小不等的空腔，淺表部巖體具架空結(jié)構(gòu)，松動(dòng)巖帶厚5.0～10.0 m等特征；和大釗等（2018）通過(guò)研究斷層破碎巖土體幾何參數(shù)和力學(xué)參數(shù)對(duì)邊坡破壞模式及穩(wěn)定性的影響，得出斷層的幾何參數(shù)對(duì)邊坡的破壞模式和穩(wěn)定性有較大的影響，并且影響程度大于其力學(xué)參數(shù)的影響；巫輔宇等（2022）分別分析降雨類型、停雨間隔時(shí)間和降雨強(qiáng)度對(duì)邊坡孔隙水壓力的影響，結(jié)果表明增強(qiáng)型降雨對(duì)邊坡穩(wěn)定性系數(shù)影響最大，在總降雨量相同時(shí)，停雨間隔時(shí)間縮減對(duì)邊坡穩(wěn)定性往往不利；劉楊等（2022）為探究降雨對(duì)含緩傾軟弱夾層礦山高邊坡穩(wěn)定性的影響，運(yùn)用強(qiáng)度折減法計(jì)算得出降雨會(huì)對(duì)降雨滲流通道產(chǎn)生阻礙，進(jìn)而孔壓曲線發(fā)生畸變，邊坡豎直截面上的孔壓增長(zhǎng)率自上而下呈減小趨勢(shì)，因此在相同時(shí)間段內(nèi)浸潤(rùn)深度與坡面孔壓的增量逐步減少。

但是，上述的大多數(shù)研究成果主要考慮邊坡坡面上個(gè)別優(yōu)勢(shì)節(jié)理，無(wú)法考慮整體邊坡坡面破碎程度。本文在已有研究成果基礎(chǔ)上，以現(xiàn)場(chǎng)結(jié)構(gòu)面實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬和理論分析方法相結(jié)合，研究降雨入滲對(duì)破碎巖質(zhì)高邊坡的影響，此方法克服了傳統(tǒng)邊坡穩(wěn)定性分析只考慮優(yōu)勢(shì)節(jié)理組方法搜索效率低、工作量大、準(zhǔn)確性差等諸多問(wèn)題，為高效而精準(zhǔn)地進(jìn)行高陡破碎巖質(zhì)邊坡穩(wěn)定性分析提供了技術(shù)參考。

1 研究區(qū)概況

1.1 地形地貌與地質(zhì)構(gòu)造

研究區(qū)位于福建省福州市永泰縣境內(nèi)，穿越的地貌單元有河流階地、山間谷地、丘陵坡地等。山體陡峭地形起伏大，山地自然坡度25～35°，局部陡壁坡度可達(dá)65°。區(qū)內(nèi)邊坡最高達(dá)30 m，為高邊坡，邊坡鉆孔資料揭示，該邊坡上部為坡積粉質(zhì)黏土，下部為碎塊狀強(qiáng)風(fēng)化凝灰?guī)r和中風(fēng)化凝灰?guī)r。碎塊狀強(qiáng)風(fēng)化凝灰?guī)r有淺黃、灰黃色，巖石風(fēng)化呈碎石塊夾土狀，碎塊手可折斷，巖土等級(jí)為Ⅳ級(jí)軟巖；中風(fēng)化凝灰?guī)r有淺黃、淺灰色，厚層狀，裂隙較發(fā)育，巖土等級(jí)為Ⅴ級(jí)次堅(jiān)巖。坡面巖體呈碎裂狀或薄層狀，結(jié)構(gòu)面很發(fā)育，如圖1所示。

圖1 碎裂狀或薄層狀結(jié)構(gòu)面Fig.1 The rock mass on the slope is in cataclastic structure

1.2 結(jié)構(gòu)面統(tǒng)計(jì)

對(duì)邊坡結(jié)構(gòu)面進(jìn)行統(tǒng)計(jì)是分析邊坡破碎程度的基礎(chǔ)工作（何云勇等，2022）。開挖后邊坡坡面結(jié)構(gòu)面露頭明顯，節(jié)理面呈張開狀態(tài)且張開距離具有差異性，據(jù)現(xiàn)場(chǎng)測(cè)量結(jié)構(gòu)面張開距離為4～15 mm，且?guī)r體在早期構(gòu)造應(yīng)力、巖體自身重力、降雨入滲等因素作用及邊坡開挖卸荷作用下，巖體張開加快，使巖塊、巖屑、巖粉夾泥質(zhì)等充填物大量存在于裂隙巖塊中（王建明等，2022）。因此，此類裂隙塊體對(duì)整體巖體性質(zhì)有著實(shí)質(zhì)性的影響，控制著整體邊坡的穩(wěn)定性。

根據(jù)野外測(cè)量數(shù)據(jù)，對(duì)36組結(jié)構(gòu)面的數(shù)據(jù)進(jìn)行了詳細(xì)的統(tǒng)計(jì)分析（蔣水華等，2021）。結(jié)構(gòu)面的張開度、延續(xù)性、結(jié)合度和間距是描述結(jié)構(gòu)面的必要條件，本次統(tǒng)計(jì)主要在實(shí)測(cè)統(tǒng)計(jì)基礎(chǔ)上，針對(duì)上述36組優(yōu)勢(shì)結(jié)構(gòu)面采用K均值聚類法統(tǒng)計(jì)，分別得到了坡面巖體張開度分級(jí)圖2-a、坡面巖體跡長(zhǎng)延續(xù)性統(tǒng)計(jì)圖2-b、坡面巖體結(jié)合度統(tǒng)計(jì)圖2-c以及坡面巖體間距分級(jí)圖2-d。

由圖2-a可以得到局部張開的占6.83%，張開的占10.27%，中等寬度的占67.25%，很寬的占15.65%。因此，中等寬度—很寬裂隙的占多數(shù)。由圖2-b可以得到跡長(zhǎng)延續(xù)性很低占15.34%，低延續(xù)性占25.56%，中等延續(xù)性占43.27%，高延續(xù)性占17.83%，可見低—中等延續(xù)性裂隙占優(yōu)。由圖2-c可以得到結(jié)合度好占5.12%，結(jié)合度一般占14.36%，結(jié)合度差占60.68%，結(jié)合度很差19.84%，由此可知結(jié)合度差裂隙占優(yōu)勢(shì)。由圖2-d可以得到中等間距占76.83%，寬間距占23.17%，故中等間距裂隙的數(shù)量占優(yōu)。

圖2 破碎巖體結(jié)構(gòu)描述統(tǒng)計(jì)圖Fig.2 Description and statistics of broken rock mass structure

2 廣義Ho e k-Bro w n強(qiáng)度準(zhǔn)則

Hoek-Brown強(qiáng)度準(zhǔn)則（李良權(quán)等，2018）最初只是一個(gè)推導(dǎo)公式，理論逐漸完善后E.Hoek指出其不足，并進(jìn)行改良與修正，從1936年的混凝土強(qiáng)度理論推導(dǎo)而來(lái)，廣義的Hoek-Brown強(qiáng)度準(zhǔn)則表示為

式中：σ1、σ3分別為最大和最小有效主應(yīng)力，且壓為正；σci為完整巖塊的單軸抗壓強(qiáng)度，可由室內(nèi)試驗(yàn)確定；mb、s和a與巖體GSI指標(biāo)和巖體開挖擾動(dòng)指標(biāo)D有關(guān)，可按下式進(jìn)行取值：

式中：mi為與完整巖塊的礦物成分和結(jié)晶結(jié)構(gòu)有關(guān)的參數(shù)，反映巖石的軟硬程度，取值范圍0.001～25.0；mb為與巖體的節(jié)理化程度以及開挖擾動(dòng)有關(guān)的參數(shù)；GSI為巖體質(zhì)量指標(biāo)，根據(jù)現(xiàn)場(chǎng)觀察巖石露頭及鉆孔巖心的巖石特性進(jìn)行取值，有插值法和間接法2種便捷取值方法，取值為0～100；a和s為模型參數(shù)，a范圍取值為0.5～0.65，當(dāng)巖體為完整巖體時(shí)，a=0.5和s=1；D為現(xiàn)場(chǎng)開挖擾動(dòng)指標(biāo)（占紹祥等，2020），開挖爆破的邊坡深部巖體將會(huì)產(chǎn)生松動(dòng)，使巖體強(qiáng)度弱化與損傷，范圍取值為0～1，對(duì)于開挖未擾動(dòng)巖體，D=0。

以廣義的Hoek-Brown強(qiáng)度準(zhǔn)則為理論基礎(chǔ)，對(duì)高陡破碎巖質(zhì)邊坡進(jìn)行交叉穩(wěn)定性分析，得出邊坡在不同工況下的安全系數(shù)、應(yīng)力應(yīng)變分布規(guī)律以及位移變化情況等多個(gè)方面對(duì)破碎巖質(zhì)邊坡進(jìn)行綜合分析，確定該強(qiáng)度準(zhǔn)則對(duì)此類高陡破碎巖質(zhì)邊坡穩(wěn)定性分析計(jì)算是切實(shí)可行的。

3 邊坡穩(wěn)定性分析

3.1 計(jì)算模型與參數(shù)

本文以廣義Hoek-Brown模型（鄧小釗等，2014）為本構(gòu)模型，利用廣義的Hoek-Brown強(qiáng)度準(zhǔn)則判斷節(jié)理巖體的強(qiáng)度和變形，至少需要估計(jì)巖體3種性質(zhì)：1）完整巖石的單軸抗壓強(qiáng)度σci；2）完整巖石的Hoek-Brown常數(shù)m；3）巖體的地質(zhì)強(qiáng)度指標(biāo)GSI。

依據(jù)公式（1）、公式（2）以及研究區(qū)邊坡現(xiàn)場(chǎng)勘探工程地質(zhì)條件，完整巖石Hoek-Brown常數(shù)m（彭俊等，2013）按與片理面垂直的方向試驗(yàn)所得，通過(guò)結(jié)構(gòu)條件以及表面質(zhì)量得到地質(zhì)強(qiáng)度指標(biāo)GSI（盧書強(qiáng)等，2009），巖石的單軸抗壓強(qiáng)度σci（彭坤等，2017）通過(guò)單軸抗壓實(shí)驗(yàn)儀測(cè)得，數(shù)值模擬采用的巖土體變形參數(shù)見表1。

表1 巖土體變形參數(shù)表Tab.1 Deformation parameters of rock and soil mass

3.2 設(shè)置工況

計(jì)算工況分別為：工況Ⅰ，考慮坡面破碎旱季邊坡穩(wěn)定性，路基開挖后的自重（天然狀態(tài)）；工況Ⅱ，考慮坡面破碎雨季邊坡穩(wěn)定性，路基開挖后的自重+暴雨（飽和狀態(tài)）。工程邊坡剖面見圖3。

圖3 工程邊坡剖面示意圖Fig.3 Schematic diagram of project slope profile

3.3 邊坡穩(wěn)定性計(jì)算分析

了解降雨對(duì)于此類特殊邊坡穩(wěn)定和變形的影響，采用有限元分析軟件Midas GTS NX（黃宇，2015；張瑩，2017）對(duì)破碎邊坡降雨過(guò)程進(jìn)行了仿真分析，其中該軟件中的Hoek-Brown模型可以充分考慮坡面巖體破碎性。計(jì)算模型采用2D模型，坡面巖體本構(gòu)模型采用廣義Hoek-Brown模型，計(jì)算域內(nèi)模擬了地質(zhì)剖面圖中所有的巖層與破碎巖石，屬性為2D單元平面應(yīng)變。計(jì)算域兩側(cè)及底部約束為自動(dòng)-區(qū)域約束。巖體的初始應(yīng)力場(chǎng)按自重應(yīng)力場(chǎng)考慮，采用強(qiáng)度折減法（SRM）進(jìn)行施工階段的邊坡穩(wěn)定性分析，計(jì)算網(wǎng)格模型見圖4。

圖4 邊坡網(wǎng)格劃分Fig.4 Grid division of slope

3.3.1 位移特征分析

圖5-a、圖5-b為破碎巖質(zhì)邊坡的水平方向位移云圖，圖中數(shù)值為正，代表邊坡向臨空面方向進(jìn)行位移，由圖例及位移云圖分布直觀展現(xiàn)出邊坡受力變形時(shí)內(nèi)部位移狀態(tài)。從圖5可知，無(wú)論是工況Ⅰ還是工況Ⅱ邊坡，水平位移主要集中在坡腳。工況Ⅰ最大位移量為4.144 mm，邊坡整體穩(wěn)定；工況Ⅱ最大位移量為27.47 mm，之所以出現(xiàn)明顯差異，是因?yàn)槠扑閹r石在降雨沖刷作用下產(chǎn)生明顯滑動(dòng)，說(shuō)明破碎坡面處于欠穩(wěn)定狀態(tài)。

圖5 水平方向位移云圖Fig.5 Horizontal displacement cloud diagram

圖6-a、圖6-b為破碎巖質(zhì)邊坡的豎直位移圖，圖中負(fù)號(hào)代表沉降，由圖例及豎直位移云圖分布直觀展現(xiàn)出邊坡受力變形時(shí)內(nèi)部位移狀態(tài)。從圖6可知，無(wú)論是工況Ⅰ還是工況Ⅱ邊坡，豎直方向位移主要集中在坡腳。工況Ⅰ最大位移量為20.15 mm，邊坡整體穩(wěn)定；工況Ⅱ最大位移量為28.33 mm，出現(xiàn)明顯差異的原因，是由于破碎巖石在降雨沖刷作用下產(chǎn)生明顯沉降，說(shuō)明破碎坡面處于欠穩(wěn)定狀態(tài)。

圖6 豎直方向位移云圖Fig.6 Vertical displacement cloud diagram

3.3.2 最大剪應(yīng)力特征分析

圖7為最大剪應(yīng)力云圖。邊坡在工況Ⅰ、工況Ⅱ條件下應(yīng)力無(wú)明顯區(qū)別，自然工況下最大剪應(yīng)力為400.8 kN·m-2，降雨工況下最大剪應(yīng)力為406.7 kN·m-2。從圖7中能夠看出最大剪應(yīng)力在坡頂較小，隨埋深的增加，最大剪應(yīng)變?cè)隽棵黠@，應(yīng)力向邊坡內(nèi)部呈線性增長(zhǎng)，且最大剪應(yīng)力由上而下分布均勻，未出現(xiàn)應(yīng)力集中現(xiàn)象，表明邊坡體垂直方向主要受自身重力影響。

圖7 最大剪應(yīng)力云圖Fig.7 Cloud diagram of maximum shear stress

3.3.3 等效塑性應(yīng)變特征分析

圖8為等效塑性應(yīng)變?cè)茍D。其展現(xiàn)出邊坡潛在滑動(dòng)帶位置，工況Ⅰ下的巖質(zhì)臺(tái)階坡面位有少量塑性區(qū)，無(wú)明顯滑動(dòng)帶，安全系數(shù)為1.82，邊坡處于穩(wěn)定狀態(tài)；工況Ⅱ下的巖質(zhì)臺(tái)階貫通區(qū)明顯，坡體內(nèi)形成圓弧狀滑動(dòng)帶，邊坡穩(wěn)定性系數(shù)為1.19，處于不穩(wěn)定狀態(tài)—欠穩(wěn)定狀態(tài)。因此，該段邊坡開挖后在暴雨或連續(xù)降雨的情況下可能發(fā)生工程滑坡。

圖8 等效塑性應(yīng)變?cè)茍DFig.8 Cloud Chart of Equivalent Plastic Strain

4 結(jié)論及建議

1）研究區(qū)邊坡巖體結(jié)構(gòu)主要為破碎鑲嵌狀結(jié)構(gòu)，風(fēng)化程度較高，受降雨沖刷侵蝕嚴(yán)重，易產(chǎn)生小型的崩塌、掉塊。

2）在降雨工況作用下，其穩(wěn)定性顯著降低，邊坡整體的潛在破壞面形成，內(nèi)部發(fā)生以剪切為主的破壞。邊坡坡面沿弱層位移較大，邊坡中部移動(dòng)不明顯，邊坡水平位移表現(xiàn)為沿邊坡坡腳與弱層順層剪出，邊坡豎向位移主要表現(xiàn)為坡頂沉降，坡腳因受坡頂及坡面巖體擠壓、拉裂作用，大范圍隆起，在碎石土弱層下發(fā)生較大移動(dòng)量區(qū)域較小。

3）對(duì)于破碎巖質(zhì)高邊坡，在進(jìn)行穩(wěn)定性分析計(jì)算時(shí)，往往忽略邊坡坡面的不連續(xù)性，通過(guò)結(jié)構(gòu)面統(tǒng)計(jì)分析及廣義Hoek-Brown模型進(jìn)行穩(wěn)定性計(jì)算具有一定實(shí)際參考價(jià)值。

4）對(duì)破碎巖質(zhì)邊坡進(jìn)行治理時(shí)，建議著重注意邊坡坡腳的加固，利用截排水溝、設(shè)急流槽、平臺(tái)側(cè)溝、路塹邊溝等排水設(shè)施組成綜合地表排水系統(tǒng)，以消除降雨對(duì)破碎高邊坡穩(wěn)定性的不利影響。