鄭洪春 閆國強 孫位 馮冰洋 唐瑤

收稿日期：2023-05-15；接受日期：2023-08-29

基金項目：國家自然科學(xué)基金項目（42077234）

作者簡介：鄭洪春，男，高級工程師，博士，主要從事庫區(qū)地質(zhì)災(zāi)害及移民管理研究。E-mail：zheng_hongchun@ctg.com.cn

通信作者：閆國強，男，工程師，博士，主要從事庫區(qū)地質(zhì)災(zāi)害防治與岸坡劣化機理研究。E-mail：yan_guoqiang@ctg.com.cn

Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

文章編號：1001-4179（2024） 06-0159-08

引用本文：鄭洪春，閆國強，孫位，等.

考慮劣化效應(yīng)的順層灰?guī)r岸坡滑移-拉裂研究

［J］.人民長江，2024，55（6）：159-166.

摘要：三峽庫區(qū)在高程145～175 m之間的周期性水位調(diào)控，促使庫區(qū)巫峽段灰?guī)r岸坡劣化損傷加劇，穩(wěn)定性衰減。通過引入劣化深度L0以及巖層滑移厚度H建立了滑移-拉裂計算模型，研究巫峽段順層灰?guī)r岸坡受侵蝕劣化脫空時的滑移-拉裂破壞。研究結(jié)果表明：當(dāng)L0不變時，隨著H增加，岸坡易于失穩(wěn)且失穩(wěn)距離x明顯增大；隨著L0變大，同一H巖層面滑移失穩(wěn)距離x明顯變大，可增加約3倍有余，巖層面參數(shù)的劣化衰減使順層岸坡更易發(fā)生深厚層、大規(guī)模的失穩(wěn)破壞；當(dāng)L0較小時，易沿坡肩后部“滑移-拉裂”破壞，隨著L0逐漸變大，破壞模式變?yōu)槠录缜昂蟛客瑫r發(fā)生滑移-拉裂失穩(wěn)。針對巖層面易受庫水損傷的順層岸坡，治理方案建議為：坡肩及后緣滑移塊錨索（桿）加固+坡肩前緣巖層出露面“掛網(wǎng)+噴漿+隨機錨桿”組合方案，以達到“固頭”以及隔絕“劣化損傷”通道的目的。研究成果可為庫區(qū)內(nèi)順層灰?guī)r岸坡巖體劣化穩(wěn)定性分析以及工程治理提供一定技術(shù)參考。

關(guān)? 鍵? 詞：順層灰?guī)r岸坡； 穩(wěn)定性分析； 巖體劣化； 劣化深度； 三峽庫區(qū)

中圖法分類號： P642.22

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

DOI：10.16232/j.cnki.1001-4179.2024.06.022

0? 引 言

由于三峽水庫庫水調(diào)度，形成了高程145～175 m之間的水位消落帶［1］。在庫水位升降循環(huán)中，變幅帶巖體處于一種浸泡-風(fēng)干交替狀態(tài)，這種浸泡-風(fēng)干交替作用對巖體結(jié)構(gòu)來說是一種“劣化損傷”作用［2-4］，促使巖溶岸坡消落帶巖體節(jié)理裂隙發(fā)育。而節(jié)理裂隙與結(jié)構(gòu)面處往往是巖溶以及巖體劣化最活躍場所［5-7］，進一步加速了碳酸鹽岸坡變形破壞演化與巖體劣化進程［8］。在三峽庫區(qū)巫峽段廣泛分布有順層灰?guī)r岸坡，順層岸坡本身作為極易誘發(fā)大型致命滑坡的結(jié)構(gòu)類型［9］，同時疊加巖層面（結(jié)構(gòu)面）劣化效應(yīng)，極易發(fā)生沿層面的滑移失穩(wěn)破壞［10-11］。

中國首例順層巖質(zhì)水庫滑坡為1961年湖南省柘溪塘巖光滑坡，方量140萬m3，滑速達25 m/s，涌浪21 m，死亡40余人［12］，經(jīng)濟損失巨大。1963年震驚世界的瓦伊昂滑坡即為典型的順層巖質(zhì)滑坡，方量高達 2.6億m3，滑速達20～30 m/s，雖然災(zāi)前顯示出明顯滑動征兆，布置了專業(yè)監(jiān)測儀器，但仍未能避免悲劇的發(fā)生［13］。不少學(xué)者［14-17］對瓦伊昂滑坡的啟動及快速滑移機制進行了探索，認(rèn)為滑帶受庫水作用軟化、力學(xué)參數(shù)急劇衰減是導(dǎo)致岸坡失穩(wěn)并快速滑移的主要原因［18-20］。2003年7月13日三峽水庫初蓄水1月后發(fā)生千將坪順層巖質(zhì)滑坡，方量約1 500萬m3，造成14人死亡，10人失蹤，經(jīng)濟損失慘重。不少學(xué)者［21-23］對千將坪滑坡進行研究總結(jié)后發(fā)現(xiàn)，千將坪滑坡的失穩(wěn)啟動機制與瓦伊昂滑坡頗為相似。順層巖質(zhì)邊坡變形失穩(wěn)的主控界面是巖層層面［24-25］，其最典型的兩種破壞模式分別為：滑移-拉裂與滑移-彎曲［8，26］。傳統(tǒng)的順層岸坡穩(wěn)定性評價較少考慮巖層面（結(jié)構(gòu)面）的劣化損傷特性，在邊界約束條件不變的情況下，穩(wěn)定系數(shù)為一定值，這與巫峽段順層灰?guī)r岸坡實際不符。隨著庫水周期的干濕循環(huán)、冷熱交替、物理化學(xué)溶蝕等，即使是相同的邊界約束，由于岸坡巖層面自身強度參數(shù)、各力學(xué)指標(biāo)的衰減劣化，岸坡的穩(wěn)定性呈下降趨勢。因此，對于巖體劣化特征明顯的巫峽段順層灰?guī)r岸坡的穩(wěn)定性評價，亟待引入一種考慮巖體劣化效應(yīng)的順層岸坡概化模型及穩(wěn)定性計算方法。

本文結(jié)合三峽庫區(qū)巫峽段順層灰?guī)r岸坡劣化特征，建立考慮巖體劣化效應(yīng)的順層滑移-拉裂破壞模型，引入巖層面劣化深度L0，探討了不同劣化深度L0以及巖層厚度H對順層灰?guī)r岸坡的穩(wěn)定性影響，以及對應(yīng)的極限滑移失穩(wěn)距離x，最后針對該類型劣化岸坡提出了相應(yīng)的防治建議。

1? 巫峽段順層岸坡宏觀劣化特征

2008～2021年間對三峽庫區(qū)巫峽段巖溶岸坡進行了長期跟蹤調(diào)查，調(diào)查區(qū)兩岸總長43.5 km，順層岸坡長為20.68 km，占比48%，其中巫峽右岸順層岸坡占右岸長度比高達62%，左岸順層岸坡占比相

對較小，為33%。特別是青石村至抱龍河口段兩岸

順層灰?guī)r岸坡劣化現(xiàn)象極為顯著，此處兩岸劣化帶順層岸坡長達9.29 km，占巫峽順層岸坡總長45%（表1、圖1）。

岸坡巖體在長時間物理化學(xué)作用下，巖體強度逐漸劣化，會促使已有節(jié)理裂隙進一步加寬變深。庫水沿著新生擴展裂隙進行物理化學(xué)侵蝕，又進一步加劇了巖體劣化。巫峽區(qū)順層岸坡底部受侵蝕劣化脫空，形成良好臨空面，上部巖層出露，存在潛在的沿層面滑移-拉裂破壞風(fēng)險（圖2（a））。部分岸坡由于30 m消落帶長期的水巖劣化作用，導(dǎo)致岸坡趾部軟化，在后部上覆山體壓力作用下，形成了滑移-彎曲隆起［26］（圖2（b））。根據(jù)前期調(diào)查研究［6-7］，巖體劣化主要沿著結(jié)構(gòu)面進行，且由表層向岸坡深部劣化作用趨弱，存在著一個臨界劣化深度L0。順層巖質(zhì)岸坡之間的巖層面以及分布的泥化夾層本身就是易滑結(jié)構(gòu)，疊加周期循環(huán)的“水-巖-結(jié)構(gòu)面”劣化作用，岸坡穩(wěn)定狀態(tài)大大削弱，極易形成大規(guī)?；破茐摹Ｏ抻谄?，本文主要針對考慮巖體劣化效應(yīng)的滑移-拉裂模型進行研究。

2? 考慮巖體劣化效應(yīng)的順層滑移-拉裂破壞模型

根據(jù)巫峽段順層灰?guī)r岸坡宏觀劣化現(xiàn)象（圖2（a））概化順層滑移-拉裂破壞模型如圖3所示。其中β為坡角，α為巖層傾角，各巖層近似平行，后部拉裂面BiC近似與巖層面垂直。任意選取第i層巖面與坡面交界處為原點Ai（0，0），順巖層方向作x軸，以垂直x軸向上作y軸，建立直角坐標(biāo)系。頂部坡肩位置D0已知為（xD0，yD0）。據(jù)已有研究成果［2，6-7，25］，巖體劣化的強弱不僅與巖性以及所處的工況環(huán)境有關(guān)，同時與巖體埋深以及與天然界面的距離有關(guān)，巖體埋深越深，距外部劣化作用面越遠，其劣化效應(yīng)越弱［27］。這里為了便于研究，采用近似線性相關(guān)的空間劣化函數(shù)式，即任意第i層巖面中黏聚力Ci、摩擦系數(shù)tanΦi、后緣拉裂面抗拉強度σti距外界劣化作用原點Ai（0，0）距離x越遠，受到劣化作用越弱，強度越高，直至一個臨界劣化深度L0時，劣化作用可忽略不計，此時的巖體物理參數(shù)可采用原位新鮮巖體的強度參數(shù)值。由圖3建立巖體參數(shù)空間劣化函數(shù)式如下：tanΦi=kΦix+tanΦi，0≤x≤L0tanΦi0，x＞L0（1）

Ci=kCix+Ci，0≤x≤L0Ci0，x＞L0（2）

σti=kσtx+σti，0≤x≤L0σti0，x＞L0（3）

式中：kΦi，kCi，kσt為第i層巖面各物理參數(shù)對應(yīng)空間線性劣化系數(shù)；Ci，tanΦi，σti為第i層巖面經(jīng)過庫水多次循環(huán)劣化作用后的強度數(shù)值；tanΦi0，Ci0，σti0為原巖層面強度參數(shù)；x值為沿所在x軸距原點Ai（0，0）點的距離。

對比巫峽劣化岸坡與順層滑移-拉裂概化模型（圖3）不難發(fā)現(xiàn)，依據(jù)極限滑移距離x所在位置可劃分為兩種破壞類型：當(dāng)0≤x≤xD0沿著坡肩D0前緣AiD0界面滑移破壞（圖3（a））；當(dāng)x≥xD0沿著坡肩D0后緣D0C界面滑移破壞（圖3（b））。下文對此兩種破壞模式分別求解其滑移面深度x。

2.1? 沿坡肩D0前緣滑移拉裂破壞

沿坡面走向取單位長度，對三角塊體AiBiC進行分析，對其下滑力FSi與阻滑力FRi分別進行求解。當(dāng)FSi與FRi相等時即為臨界破壞點，令FRi≤FSi即可解得失穩(wěn)滑移距離x的范圍。

由圖3（a）可知：AiBi=x，BiC=xtan（β-α），WAiBiC=12γAiBiBiC；拉裂面BiC上抗力有Fσt =BiCσti；滑移面AiBi上抗力有Fpt =WAiBiCcosαtanΦi+xCi；滑移面AiBi上下滑力有FSi=WAiBiCsinα；則對于三角塊體AiBiC而言，總阻滑力FRi = Fσt + Fpt；將相關(guān)參數(shù)代入下滑力FSi與阻滑力FRi有：FSi =WAiBiCsinα

=12γAiBiBiCsinα=12γx2tan（β-α）sinα（4）

FRi=Fσt + Fpt

=BiCσti+WAiBiCcosαtanΦi+xCi（5）

式中：γ為第i層巖層面上的加權(quán)平均重度，γ=i1γihi/i1hi；tanΦi與Ci為滑移面AiBi上加權(quán)平均強度參數(shù)，且由式（1）～（3）知，當(dāng)0≤x≤L0與x＞L0時，兩參數(shù)的加權(quán)平均值不一致。

當(dāng)0≤x≤L0時，為

tanΦi=tanΦi+kΦix+tanΦi2=tanΦi+kΦi2xCi=Ci+kCix+Ci2=Ci+kΦi2x（6）

當(dāng)x＞L0時，為

tanΦi=tanΦi0+tanΦiL02x+tanΦi0x-L0xCi=Ci0+CiL02x+Ci0x-L0x（7）

將上述約束條件分別代入到式（5），且令FRi ≤FSi，則化簡可得：當(dāng)0≤x≤L0時，為

14γtanβ-αcosαkΦix2+

12tanβ-α2kσt+γcosαtanΦi-γsinα+12kCix+

12tanβ-ασti+Ci≤0（8）

當(dāng)x＞L0時，為

12γtanβ-αcosαtanΦi0-sinαx2+14tanβ-α4σti0+γcosαtanΦi-tanΦi0L0+Ci0x+Ci-Ci02L0≤0（9）

由方程式（8）、（9）可以看出，滑移距離x值與坡角β、巖層傾角α、臨界劣化深度L0以及物理力學(xué)參數(shù)密切相關(guān)。

2.2? 沿坡肩D0后緣滑移拉裂破壞

沿坡面走向取單位長度，對梯形塊體AiBiCD0進行分析，對其下滑力FSi與阻滑力FRi分別進行求解。當(dāng)FSi與FRi相等時即為臨界破壞點，令FRi ≤FSi即可解得失穩(wěn)滑移距離x的范圍。

由圖3（b）可知：AiBi=x，BiC=i1hi=H，WAiBiCD0=12γHxD0+γHx-xD0；拉裂面BiC上抗力有Fσt =BiCσti；滑移面AiBi上抗力有Fpt=WAiBiCD0×cosαtanΦi+xCi；滑移面AiBi上下滑力有FSi=WAiBiCD0×sinα；則對于梯形塊體AiBiCD0而言，總阻滑力FRi = Fσt + Fpt；將相關(guān)參數(shù)代入下滑力FSi與阻滑力FRi有：FSi=WAiBiCD0sinα

=12γHxD0+γHx-xD0sinα（10）

FRi=Fσt+ Fpt

=BiCσti+WAiBiCD0cosαtanΦi+xCi（11）

與2.1節(jié)類似，由式（1）知當(dāng)0≤x≤L0與x＞L0時兩參數(shù)的加權(quán)平均值不一致：當(dāng)0≤x≤L0時，為

tanΦi=tanΦi+kΦix+tanΦi2=tanΦi+kΦi2xCi=Ci+kCix+Ci2=Ci+kΦi2x（12）

當(dāng)x＞L0時，為

tanΦi=tanΦi0+tanΦiL02x+tanΦi0x-L0xCi=Ci0+CiL02x+Ci0x-L0x（13）

將上述約束條件分別代入到式（7），且令FRi ≤ FSi，則化簡可得：當(dāng)0≤x≤L0時，為

12γHcosαkΦi+12kCix2+Hkσt+Ci+γHtanΦicosα-14kΦixD0cosα-sinαx+12γHxD0sinα-cosαtanΦi+Hσti≤0（14）

當(dāng)x＞L0時，為

γHtanΦi0cosα+Ci0-γHsinαx2+? 12γHcosαtanΦi-tanΦi0L0-cosαxD0tanΦi0+? xD0sinα+Hσti0+12Ci-Ci0L0x+? 14γHxD0tanΦi0-tanΦiL0cosα≤0（15）

由方程式（14）、（15）可以看出滑移距離x值與坡角β、巖層傾角α、臨界劣化深度L0、滑移層面與坡肩D0的垂直距離H以及物理力學(xué)參數(shù)密切相關(guān)。

3? 考慮巖體劣化效應(yīng)的順層滑移-拉裂實例

巫峽段某順層灰?guī)r岸坡計算模型如圖4所示。其中有β=75°、α=45°、γ=28 kN/m3，根據(jù)鉆孔資料可知劣化深度L0大約在6～16 m不等［6-7］，且岸坡坡肩上部巖層偏薄、向下偏厚（圖4）。根據(jù)前期室內(nèi)以及相關(guān)原位試驗，測得巖層層間軟弱結(jié)構(gòu)面力學(xué)參數(shù)劣化前有：tanΦi0=tan37°，Ci0=28.6 kPa，σti0=28 kPa。劣化后力學(xué)參數(shù)為：tanΦi=tan15°，Ci=10 kPa，σti=0 kPa。為了探討劣化深度L0以及距離坡肩不同位置H巖層面極限失穩(wěn)距離x以及破壞形式，獲得規(guī)律性的認(rèn)知，這里采取了相應(yīng)簡化約定：巖層假設(shè)為無限遠、各個巖層面的力學(xué)參數(shù)取值相同。通過Matlab編程進行了多個劣化深度L0以及滑移層厚H的探討，計算結(jié)果見表2。

由計算結(jié)果（表2）不難看出，考慮劣化后的順層岸坡穩(wěn)定性與失穩(wěn)距離x受劣化深度L0、滑移厚度H關(guān)系密切。當(dāng)劣化深度L0不變時，隨著巖層滑移厚度H的增加極限失穩(wěn)距離x也逐漸增大。特別當(dāng)H=6 m時，3種不同劣化深度L0的岸坡呈現(xiàn)出失穩(wěn)距離x無窮大的效果（由于巖層不可能無限長，此處僅是概念模型上的無限遠）。取劣化深度L0=16 m進行分析，當(dāng)H=1 m時，巖層不會滑移；當(dāng)H=2 m時，巖層沿坡肩后部滑移，且滑移范圍3.46≤x≤13.883 2；當(dāng)H=4 m時，巖層坡肩前后都可以發(fā)生失穩(wěn)滑移，坡肩前部滑移范圍4.550 1≤x≤6.93，坡肩后部滑移范圍6.93≤x≤25.962 3；當(dāng)H=6 m時，巖層坡肩前后都可以發(fā)生失穩(wěn)滑移，坡肩前部滑移范圍4.550 1≤x≤10.39，坡肩后部滑移范圍為x≥16至無限遠（圖5）。綜上可見，當(dāng)考慮巖體劣化后，隨著滑移層厚H變大，滑移失穩(wěn)極限距離和范圍也隨之增大，且變形破壞模式發(fā)生了改變：由開始的僅在坡肩后緣滑移-拉裂破壞，轉(zhuǎn)換成坡肩前后可以同時發(fā)生滑移-拉裂破壞（圖6）。當(dāng)巖層滑移厚度H一定時，隨著劣化深度L0的增加，極限滑移失穩(wěn)范圍x也越大。取巖層滑移厚度H=4 m進行橫向?qū)Ρ龋▓D7），發(fā)現(xiàn)隨著劣化深度L0的增加，巖層極限滑移失穩(wěn)范圍x明顯變大，由最初的6.93～8.00 m滑移范圍（L0=8 m）擴展為6.93～25.96 m（L0=16 m），滑移失穩(wěn)距離增加了3.25倍。綜上所述，巖層面參數(shù)的劣化衰減，使得順層岸坡更容易發(fā)生深厚層、大規(guī)模的失穩(wěn)破壞。

4? 討 論

傳統(tǒng)的順層岸坡穩(wěn)定性評價沒有考慮巖體特別是巖層面的劣化特性，在邊界條件不變的情況下，其穩(wěn)定系數(shù)為一定值。但是在引入劣化深度L0以及巖層滑移厚度H后，即使是相同的邊界約束條件下，順層灰?guī)r岸坡的穩(wěn)定狀態(tài)也是呈現(xiàn)下降趨勢，與現(xiàn)場實際調(diào)查相吻合。為定量化研究這一規(guī)律，采用Matlab編程對H=6 m，x≥xD0，且x≥L0時的滑移失穩(wěn)界限進行求解（圖8）。由圖8可以看出，隨著劣化深度L0的延展，失穩(wěn)界限x呈指數(shù)快速下降，進一步說明失穩(wěn)滑移界限對于劣化深度L0有極高的敏感性。特別是當(dāng)L0=12.12 m時，滑移失穩(wěn)界限x=12.12 m逐步向坡肩前緣演化，可以看成對于該工況（H=6 m）時，劣化深度L0=12.12 m岸坡會發(fā)生滑坡坡肩前后緣的貫通性破壞，不妨定義該劣化深度L0為“貫通破壞臨界深度”。

值得注意的是，本文對巖層面的力學(xué)參數(shù)劣化采用了近似線性的空間劣化函數(shù)式（1）進行概化，與實際情況略有出入，但可大體反映出巖體劣化的主要特征。文獻［7］指出巖體劣化主要沿岸坡結(jié)構(gòu)面進行不均勻劣化，但整體呈現(xiàn)出由表及里劣化趨弱的狀態(tài)（圖9）。文獻［6］進一步總結(jié)出兩種典型的岸坡空間劣化類型，特別指出當(dāng)巖溶岸坡巖體由大型結(jié)構(gòu)面控制時，相對結(jié)構(gòu)面劣化而言，巖石的劣化可忽略不計［28-29］，僅需對結(jié)構(gòu)面進行力學(xué)性質(zhì)函數(shù)動態(tài)變化（圖10）。本文研究的灰?guī)r岸坡即屬于典型受控于軟弱巖層面（結(jié)構(gòu)面）、泥化夾層的順層岸坡，其破壞失穩(wěn)模式可認(rèn)為主要沿巖層面滑移破壞。特別是當(dāng)巖層面受水巖劣化作用后，強度參數(shù)大幅度下降，沿層面破壞趨勢加劇。本文計算結(jié)果也同時表明，當(dāng)巖層面力學(xué)參數(shù)劣化后，極易發(fā)生沿坡肩后部的“厚層—巨厚層”滑移拉裂破壞，規(guī)模大，危險性較高。對于巖層面劣化損傷的順層岸坡治理主要從兩方面統(tǒng)籌考慮：坡肩及后部進行錨索（桿）加固，將潛在“滑移塊”錨固于下伏穩(wěn)定巖石之上；對坡肩前部的巖層出露面進行“掛網(wǎng)+噴漿+隨機錨桿”支護；一方面對破碎巖體進行支護加固，更重要的是阻隔巖層面與庫水作用的“劣化損傷通道”，阻止或減緩劣化深度L0的進一步發(fā)展（圖11）。

本文通過求解順層灰?guī)r岸坡極限滑移距離x對劣化順層岸坡滑移-拉裂破壞模式進行了詮釋。當(dāng)滑移面距離x、后緣拉裂邊界BiC條件較為明晰時，也可通過總阻滑力FRi與下滑力FSi式子之比求解對應(yīng)的穩(wěn)定系數(shù)，即：Fos=FRi/FSi。需要注意的是，當(dāng)后緣拉裂邊界BiC受節(jié)理切割貫通張開時，應(yīng)根據(jù)實際考慮后緣拉裂縫BiC是否充水以及充水高度的不利影響?？梢悦黠@看到，巖層面物理力學(xué)參數(shù)隨著水巖劣化損傷作用下降、劣化深度L0逐步變大，在坡體幾何參數(shù)、邊界約束條件不變的情況下，阻滑力FRi（即式（3）、（7））明顯呈快速下降趨勢，對應(yīng)穩(wěn)定系數(shù)Fos也隨之下降。這與現(xiàn)場調(diào)查的順層灰?guī)r岸坡由于劣化損傷，穩(wěn)定性呈衰減趨勢相吻合，表明考慮巖體劣化效應(yīng)的順層滑移-拉裂破壞模型的合理性。

5? 結(jié)論與建議

本文結(jié)合巫峽段順層灰?guī)r岸坡宏觀劣化特征，建立了考慮巖體劣化效應(yīng)的順層滑移-拉裂破壞模型，主要結(jié)論如下：（1） 巫峽青石至抱龍段順層灰?guī)r岸坡受庫水長期的物理化學(xué)作用，宏觀劣化現(xiàn)象明顯，岸坡底部受侵蝕劣化脫空，巖層出露、臨空面良好，易發(fā)生滑移-拉裂破壞。

（2） 模型中引入劣化深度L0以及巖層滑移厚度H，當(dāng)劣化深度L0不變時，隨著滑移厚度H增加，岸坡滑移失穩(wěn)距離x明顯增大。隨著劣化深度L0逐漸發(fā)展，同一H巖層面滑移失穩(wěn)距離x明顯變大。表明巖層面參數(shù)的劣化衰減，使得順層岸坡更容易發(fā)生深厚層、大規(guī)模的失穩(wěn)破壞。

（3） 順層劣化岸坡極易發(fā)生沿坡肩后部的滑移-拉裂破壞，當(dāng)劣化深度L0逐漸發(fā)展至“貫通破壞臨界深度”時，破壞模式轉(zhuǎn)變?yōu)槠录缜昂蟛客瑫r發(fā)生滑移-拉裂失穩(wěn)。對于此類巖層面易受庫水損傷的順層岸坡，治理方案建議采用：坡肩及后緣進行滑移塊錨索（桿）加固+坡肩前緣巖層出露面“掛網(wǎng)+噴漿+隨機錨桿”組合方案，達到“固頭”以及隔絕“劣化損傷”通道的目的。

（4） 在邊界約束不變的前提下，本模型可較為真實反映順層灰?guī)r岸坡隨劣化深度L0發(fā)展穩(wěn)定系數(shù)下降的趨勢，與現(xiàn)場調(diào)查一致。后續(xù)工作中建議加強順層灰?guī)r岸坡消落帶處巖層面原位損傷探測工作，建立巖體劣化時空演化解析解，為后續(xù)巖體劣化評估預(yù)測及穩(wěn)定性分析奠定基礎(chǔ)。

參考文獻：［1］? 陳輝，趙云發(fā)，梁志明.金沙江下游-三峽梯級水庫聯(lián)合調(diào)度技術(shù)及實踐［J］.人民長江，2022，53（11）：203-210.

［2］? 鄧華鋒，周美玲，李建林，等.水-巖作用下紅層軟巖力學(xué)特性劣化規(guī)律研究［J］.巖石力學(xué)與工程學(xué)報，2016，35（增2）：3481-3491.

［3］? 汪標(biāo)，易慶林，鄧茂林，等.降雨與庫水作用下譚家河滑坡變形響應(yīng)規(guī)律分析［J］.人民長江，2023，54（4）：141-149.

［4］? HALE P A.A laboratory investigation of the effects of cyclic heating and cooling，wetting and drying，and freezing and thawing on the compressive strength of selected sandstones［J］.Environmental and Engineering Geoscience，2003，9（2）：117-130.

［5］? DEREK F，PAUL W K.Karst hydrogeology and geomorphology［M］.Chichester：John Wiley & Sons Ltd，2007.

［6］? 黃波林，殷躍平，張枝華，等.三峽工程庫區(qū)巖溶岸坡消落帶巖體劣化特征研究［J］.巖石力學(xué)與工程學(xué)報，2019，38（9）：1786-1796.

［7］? 殷躍平，閆國強，黃波林，等.三峽水庫消落帶斜坡巖體劣化過程地質(zhì)強度指標(biāo)研究［J］.水利學(xué)報，2020，51（8）：1-14.

［8］? 張倬元，王蘭生，王士天，等.工程地質(zhì)分析原理［M］.北京：地質(zhì)出版社，2009.

［9］? HUNGR O，LEROUEIL S，PICARELLI L.The varnes classification of landslide types，an update［J］.Landslides，2014，11：167-194.

［10］閆國強，殷躍平，黃波林，等.三峽庫區(qū)順層灰?guī)r岸坡劣化-潰屈災(zāi)變機制研究［J］.巖土力學(xué)，2022，43（9）：2568-2580.

［11］秦盼盼，黃波林，秦臻，等.考慮巖體劣化的三峽庫區(qū)順層巖質(zhì)邊坡破壞分析［J］.人民長江，2022，53（9）：147-154.

［12］杜伯輝.柘溪水庫塘巖光滑坡：中國典型滑坡［M］.北京：科學(xué)出版社，1998.

［13］MULLER L.The Vajiont slide［J］.Engineering Geology，1987，24：513-523.

［14］SEMENZA E，GHIROTTI M.History of 1963 Vaiont slide：the importance of the geological factors to recognise the ancient landslide［J］.Bulletin of Engineering Geology，2000（59）：87-97.

［15］VOIGHT B，F(xiàn)AUST C.Frictional heat and strength loss in some rapid slides［J］.Geotechnique，1992，42（4）：641-643.

［16］MOREGENSTERN N.Stability charts for earth slopes during rapid drawdown［J］.Geotechnique，1963，13（2）：121-131.

［17］PARONUZZI P，RIGO E，BOLLA A.Influence of filling-drawdown cycles of the Vajont Reservoir on Mt.Toc slope stability［J］.Geomorphology，2013，191：75-93.

［18］MüLLER L.The Vaiont catastrophe-a personal review［J］.Engineering Geology，1987，24：423-444.

［19］HENDRON A J，PATTON F D. The Vaiont Slide，a geotechnical analysis based on new geologic observations of the failure surface［R］.Vicksburg，MS：Waterways Experiment Station，1985.

［20］肖詩榮，魏瑞琦，李瑩，等.意大利瓦依昂滑坡研究綜述［J］.人民長江，2023，54（4）：130-140.

［21］文寶萍，申鍵，譚建明.水在千將坪滑坡中的作用機理［J］.水文地質(zhì)工程地質(zhì)，2008，35（3）：12-17.

［22］WANG F W，ZHANG Y M，HUO Z T，et al.The July 14，2003，Qianjiangping landslide，Three Gorges Reservoir，China［J］.Landslides，2004，1（2）：157-162.

［23］李曉，張年學(xué)，廖秋林，等.庫水位漲落與降雨聯(lián)合作用下滑坡地下水動力場分析［J］.巖石力學(xué)與工程學(xué)報，2004，24（2）：3714-3720.

［24］殷躍平.三峽庫區(qū)邊坡結(jié)構(gòu)及失穩(wěn)模式研究［J］.工程地質(zhì)學(xué)報，2005，13（2）：145-154.

［25］馮君.順層巖質(zhì)邊坡開挖穩(wěn)定性及其支護措施研究［D］.成都：西南交通大學(xué)，2002.

［26］閆國強，黃波林，王勛，等.基于巖體劣化順層灰?guī)r岸坡滑移-彎曲失穩(wěn)機理和評價［J］.工程地質(zhì)學(xué)報，2021，29（3）：668-679.

［27］張景昱，宛良朋，潘洪月，等.考慮水-巖作用特點的典型岸坡長期穩(wěn)定性分析［J］.巖土工程學(xué)報，2017，39（10）：1851-1858.

［28］閆國強，黃波林，代貞偉，等.三峽庫區(qū)巫峽段典型巖體劣化特征研究［J］.水文地質(zhì)工程地質(zhì)，2020，47（4）：62-72.

［29］YAN G Q，HUANG B L，QIN Z，et al.Rock mass deterioration model of bank slope based on high-precision 3D multiperiod point clouds in the Three Gorges Reservoir，China［J］.Quarterly Journal of Engineering Geology and Hydrogeology，2022，55（1）：qjegh2020-100.

（編輯：鄭 毅）

Study on slip-tensile failure of consequent bedding limestone bank slope considering deterioration effect

ZHENG Hongchun，YAN Guoqiang，SUN Wei，F(xiàn)ENG Bingyang，TANG Yao

（China Three Gorges Corporation，Chengdu 610041，China）

Abstract：

The periodic water level regulation in the Three Gorges reservoir area between 145～175 m promotes the deterioration damage and stability attenuation of the limestone bank slope in Wuxia section of Three Gorges Reservoir.In this paper，a slip-tensile calculation model was established by introducing deterioration depth L0 and slip-thickness H to study the slip-tensile risk caused by erosion and deterioration of bedding limestone bank slope in the Wuxia section.The results showed that when L0 was unchanged，the bank slope was more prone to instability and the instability distance x significantly increases with the increase of H.As L0 increased，the slip instability distance x of the same H rock layer obviously increased，which can be increased by more than three times.The deterioration attenuation of the parameters of the rock strata made the bedding bank slope more prone to deep layer and large-scale instability failure.When L0 was small，it was easy to slide and crack along the back of the slope shoulder.As L0 gradually increased，the failure mode was transformed to sliding tensile instability that occurred simultaneously at the front and rear of slope shoulder.It was suggested to adopt the combined scheme of sliding block anchor cable （rod） reinforcement for the slope shoulder and the back edge+"net hanging+shotcrete+random bolt" support for the rock formation at the front edge of the slope shoulder to achieve the purpose of "fixing head" and plugging "deterioration damage" channel.The research results provide a certain technical reference for the stability analysis of rock mass deterioration and engineering treatment of bedding limestone bank slopes in reservoir areas.

Key words：

consequent bedding limestone slope;

stability analysis;

rock mass deterioration;

deterioration depth; Three Gorges Reservoir area