降雨條件下采場邊坡變形機理及NPR錨索控制效果

摘" 要：邊坡存在的軟弱夾層會導致邊坡巖體性質變差，并且在不良地質情況時會導致坡體沿著軟弱結構面產生滑移。以中國某露天礦區邊坡為研究背景，通過現場地質調查，室內物理模型試驗以及數值模擬的方法對采場邊坡失穩機理及NPR 錨索控制效果進行研究。結果表明：邊坡后緣處以張裂縫和剪切裂縫為主，主要為L01和L08這2條主要裂縫，L01裂縫發育深度約為22～26 m，L08裂縫位于邊坡坡頂與后緣交界處，裂縫全長233 m，存在邊坡失穩的風險；降雨造成地表滲流，軟弱夾層含水率為5%、10%、15%、20%時，邊坡的滑動為厘米級，最大水平位移分別為 14.9，18.9，21.9，18.5 cm，最大豎向位移分別為35.8，51.7，101.8，291.0 cm；滲流對坡面巖體影響明顯，邊坡產生滑動的風險較大；NPR錨索對滑動邊坡具有很強的補償力學效應，預緊力平均在300 kN左右，比傳統錨索更具有抵抗圍巖滑動的優勢。NPR錨索對滑坡的補償力學行為可以為其他邊坡的治理提供科學依據。關鍵詞：邊坡；滲流；應力；應變；數值模擬；物理模型中圖分類號：TD 121

文獻標志碼：A

文章編號：1672-9315（2024）06-1127-13

DOI：10.13800/j.cnki.xakjdxxb.2024.0611開放科學（資源服務）標識碼（OSID）：

收稿日期：

2024-08-

15

基金項目：

國家自然科學基金項目（41941018，U2468219）；山西省基礎研究計劃（202303021212371）

第一作者：王豐年，男，山西大同人，博士，工程師，E-mail：wangfn_bj@163.com

通信作者：

陶志剛，男，河北邯鄲人，教授，博士生導師，E-mail：taozhigang@cumtb.edu.cn

Deformation mechanism of stope slope and control

effect of NPR anchor cable under rainfall conditions

WANG Fengnian1，2，HE Yanpeng3，TAO Zhigang4

（1.Transportation Industry Key Laboratory of Highway Construction and Maintenance Technology in Loess Areas，

Shanxi Provincial Transportation Technology Research and Development Co.，Ltd.，Taiyuan 030032，China；

2.School of Earth Sciences and Engineering，Hohai University，Nanjing 210098，China；

3.College of Energy Science and Engineering，Xi’an University of Science and Technology，Xi’an 710054，China；

4.State Key Laboratory for Tunnel Engineering，

China University of Mining and Technology（Beijing），Beijing 100083，China）Abstract：

The existence of weak interlayers in slope can lead to the deterioration of rock mass properties，and the slope can slip along the weak structural plane in bad geological conditions.Taking the slope of a certain open-pit mining area as the research background，the mechanism of slope instability and the control effect of NPR anchor cable were studied by in-situ geological investigation，indoor physical model test and numerical simulation.The results show that：The rear edge of the slope is dominated by tension cracks and shear cracks，mainly L01 and L08;The development depth of the L01 crack is about 22～26 m，and the L08 crack is located at the junction of the slope top and back edge，with a total length of 233 m，which has the risk of slope instability;When rainfall causes surface seepage and the moisture content of weak interlayers is 5%，10%，15%，and 20%，the sliding of the slope is centimeter level，with maximum horizontal displacements of 14.9，18.9，21.9，and 18.5 cm，and maximum vertical displacements of 35.8，51.7，101.8，and 291.0 cm，respectively;Seepage has obvious influence on the rock mass on the slope，and the risk of slope sliding is greater.NPR anchor cable has a strong compensatory mechanical effect on sliding slope，and the average preload is about 300 kN，which has more advantages than the traditional anchor cable to resist the slide of surrounding rock.The compensatory mechanical behavior of NPR anchor cable to landslide can provide scientific basis for the treatment of other slopes.

Key words：slope；seepage；stress；strain；numerical simulation；physical model

0" 引" 言

每年中國因邊坡而導致的自然災害數不勝數，造成了較多人員及財產損失，因此對危險邊坡進行防治與治理工作顯得尤為重要。抗滑樁、錨桿加固［1-2］等在邊坡領域得到了廣泛的應用，其中錨桿支護由于周期短，施工快捷［3］，因此得到了近幾年得到快速發展［4］。但由于傳統支護通常以被動支護為主，只能短期對邊坡失穩巖體進行控制，因此需要一種具有超常力學性能的材料結構對含軟弱結構面的降雨邊坡進行控制。

目前國內外學者對失穩邊坡的研究已經取得了一系列成果。數值模擬方面：楊志斌利用FLAC3D與蠕變模型相結合的研究方法對高邊坡的穩定性進行研究［5］；高丙麗通過利用Ansys對邊坡穩定性進行研究，并結合強度折減法得到邊坡失穩破壞的模式［6］；孫康等采用FLAC3D數值模擬方法對高填方邊坡進行研究，研究表明此類邊坡表層及淺層破壞為主要破壞模式，應采取適當防護［7］；

姬雨菲等通過研究發現在地震波頻率與邊坡土體頻率相近時，會使得邊坡出現共振現象，從而導致邊坡穩定性下降［8］；孫東等采用有限元模擬軟件對強風化邊坡進行研究，探究含水率變化對此類邊坡穩定性影響［9］；楊成等采用數值模擬方法探究順層邊坡在發生失穩破壞時的臨界長度，并通過與現場調查相比對，證明數值模擬方法的可行性［10］；馬玉福等基于離散元模擬方法提出一種三維計算方法，此方法對含裂隙邊坡的計算具有較大優勢，簡化了模型前處理時網格剖分的難度［11］；吳書晨等采用數值模擬研究對提出的邊坡穩定性計算模型進行驗證分析，結果表明提出的模型提高了計算效率且精度較高［12］；高亮等通過數值模擬對大廣高速邊坡的破壞模式以及防護方法進行研究［13］；范家瑋等通過數值模擬與現場調查、試驗等方法對黃土邊坡進行探究，研究總結了邊坡破壞模式并對防護措施的效果進行探究，結果表明防護體系可以有效控制邊坡變形［14］；高振宇等采用3DEC研究方法對邊坡穩定進行研究，探究凍融次數的變化對邊坡穩定性的影響，研究發現邊坡凍融次數增加會增加邊坡失穩風險［15］；李屹等通過采用有限元模擬方法對邊坡在爆破荷載下的變形情況進行探究，研究可推算出在礦區開采中安全爆破次數，從而避免災害的發生［16］；徐佳偉等通過離心試驗與數值模擬相結合的研究方法探究在降雨條件與地震條件下邊坡的穩定性，結果表明在地震后降雨會導致邊坡快速失穩破壞［17］；黃樂源等通過對數值模擬方法對不同工況下邊坡穩定性進行研究，結果表明在暴雨條件對邊坡穩定性影響較大，其變形量相比正常條件下增加約30%，應變最大增加60%［18］；陳景松等采用Plaxis2D數值模擬方法對邊坡在不同降雨條件下的穩定性進行探究，研究發現邊坡在均勻降雨條件下穩定性最差［19］；吳旭敏采用GeoStudio數值模擬軟件對邊坡在降雨條件下的穩定性進行研究，研究發現在短時間內的強降雨條件下邊坡穩定性會大幅下降，這對邊坡安全有極大影響［20］；張謙等通過數值模擬方法探究弱層邊坡穩定性的影響因素，結果表明弱層的傾角以及位置對邊坡穩定性具有影響，在傾角不變時，弱層靠近邊坡肩部則邊坡越容易失穩［21］；王建明等通過數值模擬方法探究不同高度的邊坡在卸載后的破壞模式［22］；張恒通等

采用有限元數值模擬方法對順層邊坡破壞模式進行研究，并針對邊坡的破壞模式提出防治措施［23］；李宗鴻等通過UDEC數值模擬軟件對陡傾邊坡在開挖條件下的失穩模式進行研究，并提出了該邊坡失穩模式的分析方法［24］；LI等開展了室內縮尺模型試驗，認為錨索可以有效控制巖體出現的剪切變形和剪切擴展［25］。

以上學者對邊坡的變形控制研究做出了貢獻，但大多只是從邊坡的失穩機理、破壞模式、加固方法等單一角度進行研究，并未將失穩機理、加固方法、監測預警形成一體化來進行分析，需要繼續深入研究。利用現場調查、赤平極射投影法等分析了邊坡的失穩風險，再用數值模擬的方法進一步分析了在降雨滲流加持效果下的邊坡變形機理，最后利用相似物理模型和現場試驗的方法提出了NPR錨索對失穩邊坡的控制效果。通過對具有補償力學效應的NPR錨索進行研究，為其他采場邊坡的治理提供科學依據。

1" 采場邊坡工程概況

1.1" 軟弱結構面分布

根據歷史鉆孔勘察資料可知，邊坡軟弱結構面多為白云巖和頁巖這2種軟硬程度不一致的巖石互相接觸產生，軟硬巖石的接觸處由于力學性質的差異，容易在外力（如重力、降雨等）作用下產生應力集中，演變成潛在滑動面，如圖1所示。

1.2" 研究區節理裂縫統計

黃峪口白云巖礦邊坡現場工程地質調查中選取測線法，選取了地質條件較好，地層較為平整的地段作為測線布置點。以邊坡走向為基礎，長度設置在10～50 m不等，盡量使得測線分布較為均勻。共布設4條測線，測線位置布置如圖2所示。

1.2.1" 平極射投影法節理裂隙統計

現場工程地質調查結果結合已有地質資料報告，對黃峪口采場邊坡出露巖體的節理裂隙產狀進行了大量整體性統計，將統計結果記錄在隨身攜帶的表格中，后續需要對本次調查結果進行數據整理，按10°為一組將本次勘察記錄結果進行分組，將每組的數量、走向、傾角等數據分別整理完成后采用DIPS 5.0軟件對節理和裂隙的分布情況進行了可視化分析，繪制出每一處測線上的節理裂隙走向玫瑰花圖及節理分布密度統計圖。

采場邊坡前緣測線F處節理裂隙統計，如圖3所示。邊坡前緣節理的走向多為東北到西南走向，N45°E方向上該組節理數目最多，為8條。與邊坡走向對比后發現，邊坡下緣節理裂隙走向與邊坡走向基本保持垂直，則稱該種節理為傾向節理，傾角大小在多數為0°或45°。同時在邊坡下緣處也有少部分西北到東南走向的節理，N20°W方向上節理數目最多為4條。其走向與邊坡走向基本平行或者兩者角度基本保持在較小角度范圍內，則稱該種節理為走向節理或斜向節理，且傾角大小在多數為0°，有少部分傾角可達60°。

采場邊坡后緣測線E處節理裂隙統計，如圖4所示。邊坡后緣節理裂隙的走向多為西北到東南走向，N80°W方向上該組節理數目最多，為5條。其次為N20°W方向，為4條。在邊坡上緣處也有少部分東北到西南走向的節理，N20°E方向上節理數目最多為3條。邊坡后緣節理裂隙傾角大部分為0°，少數為60°。走向與邊坡呈較小角度，稱這種節理為走向節理或斜節理。

采場邊坡前緣測線G處節理裂隙統計，如圖5所示。采場邊坡前緣在N45°E的節理裂隙數量最多，為8條。根據節理裂隙分布密度統計圖可知，該方向上裂隙傾角為90°，因此可以初步判斷采場邊坡前緣由于上部巖體質量過大，導致前緣處于極度受壓狀態，因而出現豎向的節理裂隙。同時可以發現其他走向的節理裂隙數量較為平均。

采場邊坡后緣測線H處節理裂隙統計，如圖6所示。采場邊坡后緣沿東北到西南方向走向的節理發育平均，且每組最多為2條。沿西北到東南方向的節理數量更少，只有一組統計個數為2條，走向為N70°W且分組后每組最大節理數為2，初步判斷采場邊坡后緣有破碎的風險。

1.2.2" 現場調查與分析

采場邊坡在2008年前為一采石場，主要開采白云巖作為基礎設施中的材料來使用，由于不規則的開采，使得邊坡受到的破壞嚴重，造成裂縫發育，為滑坡的出現提供了一定的致災條件。裂縫空間位置分布如圖7所示。

邊坡后緣可見大型發育裂縫共10條，分別編

號為L01，L02，…，L10，L01與L08裂縫構成了邊坡后緣的主要裂縫。L01裂縫長276 m，寬0.2至1.2 m，

發育深度約24 m。L01裂縫伴生產生了L02，

L03，L04，L05，L06及L10這6條裂縫，平均裂縫深度為1" m。裂縫發育特征見表1。

L08裂縫位于邊坡坡頂與后緣交界處，裂縫全長233 m，整體呈一條弧線。裂縫局部可見張開，張開寬度0.2～0.4 m，地面垮塌、破碎及崩積物寬度1～2 m，全段下錯，下錯高度0.5～2.5 m，平均下錯高度1 m左右。根據探槽結果，L08裂縫發育深度約為

24 m，而表層碎石土的厚度約為1～4 m，其下為基巖強風化層厚度約1.8～6.7 m，L08裂縫的發育深度也穿過基巖強風化層。L08裂縫伴生的裂縫為L07和L09，這2條裂縫均分布在L08裂縫東側（坡頂一側），長度2～3 m，寬度0.1～0.2 m，可見深0.2～0.8 m，裂縫L07發育深度0.8 m。裂縫照片如圖8所示。

邊坡后緣的L01和L08將坡體分為3級，L01裂縫上方為未擾動區域，2條裂縫之間坡面發生下錯現象。而L08裂縫下方形成陡坡，植被相對較少；坡體前緣由人工開挖形成，巖體破碎，可見表層巖體崩落。

2" 采場邊坡穩定性影響因素數值模擬

2.1" 模型的建立

邊坡是具有一定天然坡度的空間結構體，因此當發生降雨且強度較大時，雨水會沿坡面向下流動轉化為坡面徑流流出［13］。在該條件下雨水入滲滯后時間很短，基本入滲現象會隨著降雨的停止而停止。

一般性強降雨的發生過程如圖9所示。在降雨剛發生不久，由于地表巖土體處于非飽和狀態，孔隙中基本沒有水分的存在，因此水很容易進入到邊坡巖土體內部；隨著降雨的進一步持續，在坡角處地表巖土體首先出現飽和現象。隨著降雨時間的進一步增加，邊坡坡面表層巖土體逐漸進入飽和狀態，且表面徑流現象更加明顯。雨停以后，由于邊坡前緣坡角附近水流有一定的滯后現象，因此邊坡表面巖土體的暫態飽和區域的消失也是有一定的滯后現象。經過以上的變化，降落在邊坡坡面上的水的一部分會轉化為地表徑流，另一部分會沿坡體表面入滲到坡體內部，補給到內部非飽和巖土體中去，因此可能會造成浸潤線的變化。當降雨完成后，邊坡表面和內部一定深度下巖土體所占有的水分又會在太陽直射條件下以水蒸氣的形式離開邊坡。

坡體內部的滲流主要由巖土體的物理力學參數決定，其中主要影響參數有滲透系數、孔隙率等。滲透系數直接影響了邊坡在降雨條件下能進入多少水分，主要與降雨強度進行對比。若降雨強度遠遠大于滲透系數則沒有多少水分能進入到坡體內部，若降雨強度遠遠小于滲透系數則坡體內部變化幅度不大，因此根據工程類比和歷史地質資料，需要對滲透系數進行一定的優化處理。同時滲流也會影響邊坡的安全穩定，影響巖土體中有效應力的變化，從而使得邊坡穩定性發生變化，因此對參數的選擇需要再三驗證。通過現場勘察資料及室內試驗結果數據對模型參數進行選取，表2所示為選取的適當參數。根據調查推測的最危險滑動面的位置，依據此進行建模，模型320 m×20 m×166 m（長×寬×高），并對該模型進行單元網格的劃分，共劃分64 053個節點、56 768個單元，將模型導入FLAC3D中進行模擬分析。計算模型如圖10所示。

2.2" 不同含水率的軟弱夾層邊坡位移變化

不同含水率的軟弱夾層x方向位移變化，如圖11所示。隨著軟弱夾層含水率的升高，軟弱夾層處位移逐漸增大，4種工況下水平最大位移分別為14.9，18.9，21.9，18.5 cm，均發生在軟弱夾層區域。邊坡滑體滑動方向與滑床位移方向相反，可體現出邊坡沿最危險滑動面發生滑動。可見軟弱夾層含水率越大，即在自然狀態下，因為降雨、灌溉、地下水滲流等因素導致邊坡軟弱夾層巖土體含水率增加，力學性能劣化，引起邊坡水平方向位移增大，整體穩定性降低。將4種工況邊坡水平方向位移進行對比分析，軟弱夾層的存在是邊坡內部發生水平位移的誘因，當軟弱夾層性能劣化至無法承擔上部巖體荷載時，夾層上部分巖體將在水平方向上產生位移，邊坡發生滑坡。

不同含水率的軟弱夾層z方向位移變化，如圖12所示。含軟弱夾層的邊坡豎向位移以軟弱夾層為分割層，軟弱夾層以下為滑床，豎直方向位移基本為厘米級移動。軟弱夾層以上滑體部分豎直位移方向向下，含軟弱夾層邊坡5種工況下最大豎向位移分別為35.8，51.7，101.8，291.0 cm。與水平位移規律相似，含軟弱夾層的邊坡，夾層上部巖體發生滑動，且滑動位移隨著軟弱夾層含水率的升高而增加。

2.3

不同含水率的軟弱夾層邊坡最大剪應力增量變化

不同含水率的軟弱夾層最大剪應力增量云圖，如圖13所示。軟弱夾層邊坡內最大剪應力集中在分界線上方，并沿分界線分布，即該分界線處為潛在滑動面。含軟弱夾層的各邊坡最大剪應力沿軟弱夾層面分布，隨含水率的升高最大剪應力依次為 3.01×10-2，4.47×10-2，8.29×10-2，4.17×10-2 MPa。含軟弱夾層邊坡內剪應力增量主要集中在軟弱夾層面上，而無軟弱夾層較大剪應力增量分布相對分散，即分界面上部巖體自重由分界面剪切力與巖體內部自身剪切力共同承擔，含軟弱夾層邊坡沿軟弱夾層滑移，故在軟弱夾層處出現集中剪應力。

經過位移和應力分析可知看出，當含水率為20％的軟弱夾層邊坡相比較含水率為15％時最大值都有所下降，造成這種現象的原因為巖體的強度和含水率的關系不是線性增加的，不同條件下巖體對水的敏感度不同，所以有時會造成部分非線性增長情況。

3" NPR錨索控制效果

邊坡內部滑體相對于基巖相對運動的核心是作用在滑帶上的牛頓力。因此可以利用牛頓第2定律，即質量為M的物體若想由靜止狀態變為運動狀態，必須在其上作用一個合外力。

但是下滑力作為天然力學系統的內部力無法直接被現有技術裝備測量出來，但是若將天然力學系統變為人為力學系統，則該下滑力可以由換算的方法間接的測量出來。

因此何滿潮院士團隊提出了一種新的思路辦法，即將通過人為測得的力學系統與自然的相結合，通過對人為系統進行測量，從而計算出下滑力大小。其作用效果和機理如圖14所示［25］。

為了間接測量出下滑力的具體大小，需要建立合適的力學模型。在此何滿潮院士引入了“干擾力”的這一概念，采用貫穿滑體的方式，把擁有負泊松比效應的材料（NPR錨索）穿過潛在滑動面錨固在下部基巖上，并給予適當大小的張拉預應力。這時外部加入的干擾力就會與天然力學系統內部的力一起組成一套復雜力學系統。當滑帶面上的下滑力發生改變就會引起NPR錨索上的預應力的改變，通過牛頓力監測設備測量出干擾力的變化情況，就可以通過公式推導計算出滑體所受的下滑力的大小。

下滑力的計算公式為

F=k1P+k2

（1）

式中" k1=cos（α+φ）+sin（α+φ）tanφ，k2=Gcosαtanφ+cl；P為干擾力，kN；φ為NPR錨索與水平面夾角，（°）；φ為邊坡內各巖土層內摩擦角加權平均值；α為滑帶與水平面夾角，（°）；c為各巖土層黏聚力加權平均值；l為滑帶面長度，m。

3.1" NPR物理模型試驗

物理模型試驗采用隧道工程災變防控與智能建養全國重點實驗室自主開發的深部硐室物理模型試驗系統，如圖15（a）所示。模型尺寸按實際比例1∶50縮放搭建，采用含水率為15％的軟弱夾層降雨滲流后的巖體強度進行等比例配比。采用3D打印的NPR錨索和PR錨索水平打入巖體

中，模型整體開挖步數分為6步，應變采用應變片

監測且分為6層（每層5個監測點），如圖15（b）所示。

模型的強度也按實際工況的巖體強度等比例搭配。首先將重晶石粉、細河沙、石英砂、石膏粉及水按比例混合，攪拌均勻后倒入亞克力模具中，模具尺寸為40 cm×40 cm×1 cm（長×寬×高），待單元板成型后，拆掉模具。然后，將單元板置于室外進行短暫的晾曬，使板子中的水分快速蒸發，必須注意單元板不能長時間晾曬，否則很容易開裂。最后在室內經過約10 d時間的自然風干。待到單元板材干燥完全后，將板材運送至實驗室內，依據監測方案對部分石膏板黏貼應變片，在搭建過程中依據所需模型的尺寸，使用石膏板分層堆積，對于模型中不平整部分使用手鋸切割的方式對部分板材進行二次加工，使得滿足試驗的需要。最終得到的模型尺寸如圖15（b）所示。為了更好的符合實際工況，兩邊加載壓力為0.1 MPa。

3.1.1" 應變監測分析

模型應變監測曲線，如圖16所示。應變值從上到下呈現不斷減小的趨勢，且其變化幅度也較

小。0～150 min每層的監測點變化率大致相同，

150 min后各測點的應變值增長速率開始不同。在開挖完成后第1層的最大應變值約為1 102 με，第6層的最大應變值約為367 με。

由A至E點監測數據對比可知，每層靠近坡面位置的應變片受到邊坡開挖影響較大，而邊坡內部離坡面較遠位置受開挖影響較小。從第1層到第6層應變監測層中，與開挖區相距最近監測點的最大應變值分別為1 102，1 083，902，814，521和 367 με。因為第1層應變片沿著邊坡布置，所以開挖擾動作用下各個應變值相對其他層較大。隨著邊坡的繼續開挖，邊坡失穩區域繼續擴大，各層應變值繼續增加，直到試驗結束。

3.1.2" 錨索軸力監測分析

開挖階段前期，PR錨索和NPR錨索監測曲線上升幅值較為緩慢。開挖階段中期，兩者曲線呈現快速的增長的趨勢。隨后，PR錨索的軸力繼續增加至100 N，而NPR錨索由于恒阻滑移特性，發揮了本身的補償力學效應，軸力曲線保持在60 N左右。NPR錨索的補償力學行為可以對圍巖進行有效的控制，減少邊坡發生失穩的風險。

3.2" NPR錨索滑坡現場控制效果

深部滑坡牛頓力遠程監測裝備預警系統（NF-RMWS）主要由2大部分組成，即現場數據收集系統和室內數據處理系統。如圖18所示。其中現場監測系統包括：高精度牛頓力監測傳感器設備、數據發射傳輸設備模塊北斗一號、恒阻大變形預應力監測錨索、太陽能供電系統和地表位移、孔隙水壓力監測系統。

室內數據監測處理系統由數據信號接收器、服務器、信號處理器及相應軟件組成。現場監測到的預應力大小變化經過室內系統的實時處理變為下滑力的數據。之后相關軟件將下滑力數據繪制出監測曲線，供人們進行實時查詢。基于牛頓力監測原理及大量現場監測數據，可以將其分為穩定、不穩定、裂縫及滑移模式，如圖19所示。

黃峪口白云巖礦邊坡動態監測數據，如圖20所示。黃峪口白云巖礦邊坡由于得到NPR錨索的加固穩定性得到明顯提升。監測點1牛頓力在300 kN左右，總體趨勢保持穩定。監測點2牛頓力從剛開始的200 kN到第91的290 kN，隨后保持穩定。監測點3牛頓力由300 kN增長到375 kN。監測點4牛頓力由剛開始的300 kN增長到343 kN。

由上述結論可知，監測點1滿足穩定模式。監測點2和監測點4基本滿足裂縫發展模式，且為小裂縫。監測點3滿足不穩定模式。綜上所述，NPR錨索滿足邊坡加固的同時也實現了對延慶黃峪口白云巖礦邊坡的穩定性控制和遠程實時監測預警的要求。

4" 結" 論

1）邊坡后緣處分布的都是張裂縫和剪切裂縫，主要為L01和L08這2條主要裂縫，使整個坡體上緣下錯成為3級，其余裂縫都為這2條裂縫的伴生裂縫。邊坡破壞模式主要以滑移-拉裂破壞為主。

2）降雨對邊坡的影響主要是巖土體飽和孔壓和非飽和孔壓共同作用的結果，巖土體越接近于飽和其中的基質吸力就越低。雨水過大會導致邊坡前緣處基巖軟化，無法起到支撐邊坡的作用，因此會造成滑坡的發生。長時間的強降雨還可能會影響地下水位線，使其存在一定的波動，降雨造成的滲流會對邊坡坡體的穩定性會造成不利影響。

3）NPR錨索物理模型試驗和現場監測數據表明，NPR錨索具有較高的預緊力（300 kN左右），較普通錨索而言，NPR 錨索能隨著邊坡巖體大變形而發生軸向拉伸，通過恒阻體與恒阻套管之間的摩擦作用吸收緩慢變形產生的能量，并能保持恒定阻力和較大變形，為滑坡的監測預警提供一定的參考依據。

參考文獻（References）：

［1］" YIN Q，JING H W，MA G W，et al.Investigating the roles of included angle and loading condition on the critical hydraulic gradient of real rock fracture networks［J］.Rock Mechanics and Rock Engineering，2018，51（4）：1-11.［2］

LI Z，LIU S G，REN W T，et al.Multiscale laboratory study and numerical analysis of water-weakening effect on shale［J］.Advances in Materials Science and Engineering，2020，7：1-14.［3］

WANG F N，GUO Z B，QIAO X B，et al.Large deformation mechanism of thin-layered carbonaceous slate and energy coupling support technology of NPR anchor cable in Minxian tunnel：A case study［J］.Tunnelling and Underground Space Technology，2021，117：104151.［4］

TAO Z G，ZHU C，HE M C，et al.A physical modeling-based study on the control mechanisms of negative poisson’s ratio anchor cable on the stratified toppling deformation of anti-inclined slopes［J］.International Journal of Rock Mechanics and Mining Sciences，2021，138：104632.［5］

楊志斌.基于蠕變試驗下高邊坡穩定性模擬研究［J］.西安科技大學學報，2023，43（5）：988-999.

YANG Zhibin.Simulation study on stability of high slope based on creep test［J］.Journal of Xi ’an University of Science and Technology，2023，43（5）：988-999.［6］

高丙麗.水電工程左壩肩邊坡穩定性強度折減有限元分析［J］.西安科技大學學報，2014，34（5）：509-516.

GAO Bingli.Strength reduction finite element analysis on slope stability of a hydropower project left abutment［J］.Journal of Xi ’an University of Science and Technology，2014，34（5）：509-516.［7］

孫康，張浪，蔣根，等.三峽庫區巫山縣老林場高填方邊坡穩定性分析［J］.三峽大學學報（自然科學版），2024，46（5）：33-40.

SUN Kang，ZHANG Lang，JIANG Gen，et al.Stability analysis of high fill slope of Laolinchang in Wushan county，Three Gorges reservoir area［J］.Journal of China Three Gorges University（Natural Sciences），2024，46（5）：33-40.［8］

姬羽菲，邵帥，邵生俊，等.非飽和土邊坡三維地震穩定性分析［D］.西安：西安理工大學，2024.

JI Yufei，SHAO Shuai，SHAO Shengjun，et al.Three-dimensional seismic stability analysis of unsaturated soil slopes［J］.Xi’an：Xi’an University of Technology，2024.［9］

孫東，周春梅.不同含水量及滲透特性下強風化砂巖土邊坡穩定性分析［J］.安全與環境工程，2024，31（3）：178-187.

SUN Dong，ZHOU Chunmei.Stability analysis of strongly weathered sandstone soil slopes with different water contents and permeability characteristics［J］.Safety and Environmental Engineering，2024，31（3）：178-187.［10］

楊成，譙乾峰，雷松.順層巖質邊坡失穩破壞臨界長度的數值模擬分析［J］.西北地質，2019，52（4）：250-262.

YANG Cheng，QIAO Qianfeng，LEI Song.Numerical simulation analysis of critical length of instable failure for bedding rock slope［J］.Northwest Geology，2019，52（4）：250-262.［11］

馬玉福，劉建設，龔文平，等.含復雜裂隙網絡巖質邊坡三維穩定性計算方法［J］.水電能源科學，2024，42（6）：93-97.

MA Yufu，LIU Jianshe，GONG Wenping，et al.Three-dimensional stability calculation method of rock slope with complex fracture network［J］.Water Resources and Power，2024，42（6）：93-97.［12］

吳書晨，侯超群，孫志彬.初始含水率非均勻分布的無限邊坡穩定性計算模型研究［J］.合肥工業大學學報（自然科學版），2024，47（5）：690-695.

WU Shuchen，HOU Chaoqun，SUN Zhibin.Research on infinite slope stability calculation model with non-uniform distribution of initial water content［J］.Journal of Hefei University of Technology（Natural Science），2024，47（5）：690-695.［13］高亮，張騰伍，饒法強，等.公路高邊坡災變機理及滑動力響應規律分析［J］.科學技術與工程，2024，24（11）：4741-4749.

GAO Liang，ZHANG Tengwu，RAO Faqiang，et al.Analysis of disaster mechanism and sliding force response law of highway high slope［J］.Science Technology and Engineering，2024，24（11）：4741-4749.［14］

范家瑋，張玉芳，周文皎.黃土鐵路路堤滑坡破壞機制與防治技術研究［J］.防災減災工程學報，2024，44（4）：772-783.

FAN Jiawei，ZHANG Yufang，ZHOU Wenjiao.Research on the failure mechanism and prevention technology of loess railway embankment landslide［J］.Journal of Disaster Prevention and Mitigation Engineering，2024，44（4）：772-783.［15］

高振宇，苗浩東，常遠.凍融循環下地下開采對既有邊坡穩定性影響數值分析［J］.金屬礦山，2024（4）：221-227.

GAO Zhenyu，MIAO Haodong，CHANG Yuan.Numericalanalysis of the influence of underground mining on the stability of existing slopes under freeze-thaw cycle［J］.Metal Mine，2024（4）：221-227.［16］

李屹，徐繼業，冷建勇，等.循環爆破振動下邊坡變形累積效應研究［J］.礦業研究與開發，2024，44（3）：49-56.

LI Yi，XU Jiye，LENG Jianyong，et al.Study on cumulative effect of slope deformation under cyclic blasting vibration［J］.Mining Research and Development，2024，44（3）：49-56.［17］

徐佳偉，蔣瑋.地震及降雨作用下路堤邊坡變形破壞特性研究［J］.中國公路學報，2024，37（6）：76-86.

XU Jiawei，JIANG Wei.Deformation and failure characteristics of embankment subjected to earthquake and rainfall［J］.China Journal of Highway and Transportation，2024，37（6）：76-86.［18］

黃樂源，田城航，仉文崗，等.降雨作用下三維礦山邊坡的力學狀態響應［J］.土木與環境工程學報（中英文），

2024，46（6）：70-79.

HUANG Leyuan，TIAN Chenghang，ZHANG Wengang，et al.Mechanical response of 3D mining slopes under the effect of rainfall［J］.Journal of Civil and Environmental Engineering，2024，46（6）：70-79.

［19］陳景松，周春梅，吝曼卿，等.不同降雨類型下混合花崗巖路塹邊坡穩定性分析［J］.安全與環境工程，2024，31（1）：150-162.

CHEN Jingsong，ZHOU Chunmei，LIN Manqing，et al.Stability analysis of mixed granite cutting slope under different rainfall types［J］.Safety and Environmental Engineering，2024，31（1）：150-162.［20］

吳旭敏，楊溢，葉志程.巖土體邊坡在降雨入滲作用下的穩定性及影響因素分析［J］.水資源與水工程學報，2022，33（5）：189-199.

WU Xumin，YANG Yi，YE Zhicheng.Stability of rock-soil slopes under rainfall infiltration and its influencing factors［J］.Journal of Water Resources and Water Engineering，2022，33（5）：189-199.［21］

張謙，李貴平.含弱層順傾巖質邊坡穩定性影響因素分析［J］.中國安全生產科學技術，2022，18（S1）：79-85.

ZHANG Qian，LI Guiping.Analysis of the influencing factors on stability of inclined rock slopes containing weak layers［J］.Journal of Safety Science and Technology，2022，18（S1）：79-85.［22］王建明，岳星彤，李天龍，等.露天礦邊坡巖體裂隙測量與卸荷破壞的數值模擬研究［J］.中國安全生產科學技術，2022，18（8）：98-105.

WANG Jianming，YUE Xingtong，LI Tianlong，et al.Numerical simulation study on fracture measurement and unloading failure of slope rock mass in open-pit mine［J］.Journal of Safety Science and Technology，2022，18（8）：98-105.［23］張恒通，朱冬春，王鵬，等.順層巖質邊坡變形破壞機理分析［J］.公路，2021，66（3）：67-71.

ZHANG Hengtong，ZHU Dongchun，WANG Peng，et al.Analysis of deformation and failure mechanism of bedding rock slope［J］.Highway，2021，66（3）：67-71.［24］

李宗鴻，袁維，王偉，等.陡傾順層邊坡地下開挖失穩模式分析［J］.科學技術與工程，2020，20（5）：1988-1992.

LI Zonghong，YUAN Wei，WANG Wei，et al.Failure mode analysis of underground excavation of steep inclined slope［J］.Science Technology and Engineering，2020，20（5）：1988-1992.［25］

LI Y，ZHENG J，WANG F N，et al.Experimental study on the movement law of overlying rock non-pillar coal overhead mining［J］.Journal of Mountain Science，2023，20（6）：1759-1773.

（責任編輯：劉潔）