庫岸有限厚度粉質(zhì)黏土沖蝕演化規(guī)律及預(yù)測

陳勇 張書石 CHANDave

摘要：土體沖蝕流失是塌岸發(fā)展和庫岸再造的重要誘發(fā)因素。為準(zhǔn)確預(yù)測水流作用下岸坡土體流失進(jìn)程，自主研制了能夠?qū)崟r(shí)精確量測土樣質(zhì)量變化的新型沖蝕試驗(yàn)裝置，開展三峽庫岸粉質(zhì)黏土沖刷起動(dòng)和沖蝕進(jìn)程試驗(yàn)，分析有限厚度土體沖蝕演化規(guī)律和影響機(jī)制，并建立預(yù)測模型。結(jié)果表明：有限厚度土體的總沖蝕量隨時(shí)間增加但增速逐漸趨緩，沖蝕速率逐漸衰減，且沖蝕速率及其衰減程度均與孔隙比或流速成正比;將反映密實(shí)度的孔隙比和反映剪切力的水流速度嵌入預(yù)測模型能較好地?cái)M合粉質(zhì)黏土沖蝕演化進(jìn)程曲線。建立的有限厚度土體沖蝕進(jìn)程預(yù)測模型，為水流沖刷作用下岸坡形態(tài)變化趨勢預(yù)測研究提供參考。

關(guān)鍵詞：沖蝕試驗(yàn);有限厚度;沖蝕比;沖蝕速率;預(yù)測模型

中圖分類號(hào)：TV122

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

文章編號(hào)：1001-6791（2023）06-0938-10

收稿日期：2023-03-17;網(wǎng)絡(luò)出版日期：2023-09-27

網(wǎng)絡(luò)出版地址：https：∥link.cnki.net/urlid/32.1309.P.20230926.1758.004

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目（42272317;U21A2031）

作者簡介：陳勇（1980—），男，湖北南漳人，教授，博士，主要從事水土耦合作用機(jī)制及庫岸邊坡穩(wěn)定性分析等研究。E-mail：cyonger@126.com

水流沖刷作用下土體運(yùn)移問題，對(duì)庫岸再造、河道遷移、岸坡穩(wěn)定影響巨大，不僅導(dǎo)致水土流失、加劇生態(tài)惡化、降低水庫庫容，甚至造成財(cái)產(chǎn)損失和人員傷亡［1-3］。如何有效評(píng)價(jià)和預(yù)測土體沖蝕發(fā)展規(guī)律，是江、河、湖、庫等水體岸坡地質(zhì)災(zāi)害防治和生態(tài)環(huán)境保護(hù)中始終面臨的重要挑戰(zhàn)。土體顆粒的起動(dòng)和沖刷特性是預(yù)測沖蝕量以及岸坡或河灘形態(tài)變化的重要指標(biāo)［4-5］，國內(nèi)外學(xué)者多采用開放水槽試驗(yàn)或活塞式裝置研究海岸或河灘沉積物遷移規(guī)律和抗蝕性能，以水流濁度變化為監(jiān)測指標(biāo)，認(rèn)為土體黏粒含量、初始含水率、密度、黏聚力、表面粗糙度、鈉吸附速率等對(duì)起動(dòng)流速（臨界剪切力）和沖蝕速率影響顯著［6-10］。然而，在相同流速下相同土體的沖蝕速率隨時(shí)間發(fā)展的演化規(guī)律，也是預(yù)測沖蝕進(jìn)程的重要條件，尚有待進(jìn)一步研究。

在試驗(yàn)研究基礎(chǔ)上，描述土體沖蝕速率的數(shù)學(xué)模型日漸豐富。經(jīng)典沖蝕模型包括Wilson［11］提出的基于顆粒受力分析的剝離模型或稱威爾遜模型，以及Hanson等［12］提出的基于臨界剪應(yīng)力和可蝕性系數(shù)的超剪應(yīng)力模型，為研究不同土體在不同流速下的抗蝕性提供了總體框架。為考慮有一定埋深的黏土沉積物在重力固結(jié)后抗蝕性能更高，有學(xué)者先后將土體干密度或土樣所處的深度嵌入表面侵蝕速率預(yù)測模型中［13-15］，但尚不能預(yù)測土層厚度有限時(shí)的侵蝕發(fā)展規(guī)律。事實(shí)上，岸坡并非無限深度的均勻土體或沉積物，常見為上土下巖結(jié)構(gòu)，其沖蝕進(jìn)程隨著侵蝕發(fā)展和深度增加應(yīng)是逐漸衰減的，若按無限深度土體等速?zèng)_蝕預(yù)測庫岸侵蝕量或后退進(jìn)程，結(jié)果勢必偏保守且誤差較大。因此，分析有限厚度土體的沖蝕發(fā)展規(guī)律、建立更貼合實(shí)際的預(yù)測模型，對(duì)岸坡土體侵蝕進(jìn)程研究更有意義。

鑒于此，本文基于研制的新型沖蝕試驗(yàn)裝置，通過沖刷起動(dòng)試驗(yàn)和沖蝕進(jìn)程試驗(yàn)揭示有限厚度土體的沖蝕發(fā)展規(guī)律及其內(nèi)在原因，分析孔隙比和水流速度對(duì)粉質(zhì)黏土沖蝕演化規(guī)律的影響機(jī)制，并采用指數(shù)模型和雙曲線模型刻畫沖蝕曲線，建立同時(shí)考慮土體起動(dòng)流速、孔隙比、水流速度的預(yù)測模型，為確定有限厚度土質(zhì)岸坡的沖蝕演化趨勢提供支撐。

1 試驗(yàn)概況與土樣特性

1.1 新型試驗(yàn)裝置與方法

本研究采用自主研制的新型開放式有機(jī)玻璃矩形循環(huán)水槽試驗(yàn)裝置開展土體沖蝕試驗(yàn)，如圖1。該裝置包括上部水槽、下部水槽、帶電子天平的承重臺(tái)、通過多根傳力桿與承重臺(tái)相連的試樣盒、變頻水泵、非均勻多孔穩(wěn)流板、海綿擋板、流速儀，以及實(shí)時(shí)采集電子天平數(shù)據(jù)的軟件。水槽斷面為30×30 cm，結(jié)合便攜式流速儀測量水流速度。裝置結(jié)構(gòu)簡單、使用方便、能實(shí)時(shí)自動(dòng)獲取試樣盒內(nèi)土體質(zhì)量變化，在土樣上下游分別設(shè)置非均勻孔徑的多孔穩(wěn)流板確保流速在深度上的均勻性，并結(jié)合質(zhì)量變化和人工觀測精確獲得不同土體的起動(dòng)流速，可建立土樣顆粒沖蝕量隨時(shí)間的發(fā)展規(guī)律，從而確定土體沖蝕速率的演化特征。

1.2 試驗(yàn)土樣及試樣制備

試驗(yàn)土樣取自三峽庫岸消落帶內(nèi)，現(xiàn)場采用蠟紙密封，土樣呈紅色。測定土樣相對(duì)密度為2.739，天然含水率為20.54%，天然密度為1.62 g/cm3，孔隙比（孔隙體積與土顆粒體積之比）為1.04，孔隙率為51%。結(jié)合篩分法和密度計(jì)法對(duì)試驗(yàn)土樣進(jìn)行顆粒級(jí)配分析，顆粒級(jí)配曲線見圖2，其中粒徑1 mm以下含量為99.74%，小于0.075 mm的細(xì)粒含量為79.46%，不均勻系數(shù)（Cu）為15.6，曲率系數(shù)（Cc）為2.98，級(jí)配良好。粒徑分布特征值d60、d50、d30、d10對(duì)應(yīng)粒徑分別為0.060 8、0.053 3、0.026 6、0.003 9 mm。采用液塑限聯(lián)合法測定土樣液限為37.64、塑限為23.73，按照《巖土工程勘察規(guī)范：GB 50021—2001》將土樣定名為粉質(zhì)黏土。

根據(jù)設(shè)置的干密度，在試樣盒內(nèi)分層壓實(shí)制備重塑土試樣，保證試樣尺寸相同，如圖3，抽真空飽和24 h后取出，拆除試樣上表面和兩端擋板，并浸入靜止水中，試樣上表面和前后端與水接觸，而兩側(cè)面和底面不受水流沖刷作用。試樣上游臨水端面面積為20 cm2，僅占水流斷面（約600 cm2）的1/30，忽略土樣端面對(duì)水流流速及局部水深的影響，可較好地模擬現(xiàn)實(shí)岸坡的沖刷環(huán)境。

1.3 試驗(yàn)方案

根據(jù)三峽庫區(qū)深泓處垂線水流流速分布觀測資料，發(fā)現(xiàn)斷面水流平均流速集中于0.3～1.0 m/s［16］，試驗(yàn)中可借助水泵檔位和閥門開啟度控制水流穩(wěn)定。首先針對(duì)5組不同干密度的庫岸粉質(zhì)黏土開展沖刷起動(dòng)試驗(yàn)，從零開始逐步增大水流速度直至土樣質(zhì)量開始變化，確定其臨界起動(dòng)流速（v0）;然后針對(duì)3種不同干密度試樣開展沖蝕進(jìn)程試驗(yàn)，保持水流速度穩(wěn)定，測定了水流沖刷的3個(gè)穩(wěn)定流速（v），實(shí)時(shí)記錄試樣水下質(zhì)量變化，獲得沖蝕量隨時(shí)間的發(fā)展規(guī)律。試驗(yàn)方案見表1。

2 試驗(yàn)結(jié)果與分析

2.1 粉質(zhì)黏土起動(dòng)流速測定

砂土沖刷以單個(gè)顆粒起動(dòng)為主，而黏土沖刷主要以微團(tuán)形式起動(dòng)［17］。本研究通過試樣水下質(zhì)量變化結(jié)合人工觀測確定粉質(zhì)黏土的起動(dòng)流速和孔隙比（e）的關(guān)系，見表1，可以看出，起動(dòng)流速隨著土樣孔隙比的減小而增大，其內(nèi)在機(jī)理是孔隙比與土顆粒凈距和分子吸引力直接相關(guān)，孔隙比越小，黏聚力越大，粒團(tuán)抗分解能力越強(qiáng)，越難剝離起動(dòng)［18］。

根據(jù)試驗(yàn)結(jié)果對(duì)起動(dòng)流速和孔隙比相關(guān)性進(jìn)行擬合，發(fā)現(xiàn)粉質(zhì)黏土起動(dòng)流速與孔隙比呈負(fù)指數(shù)遞減關(guān)系，表達(dá)式為

v0=a0exp（b0e）（1）

本研究中土樣試驗(yàn)參數(shù)a0=0.525 m/s，b0=-2.564，決定系數(shù)R2=0.984 1。許琳娟等［19］對(duì)黏性土起動(dòng)流速影響因素試驗(yàn)分析中發(fā)現(xiàn)土體起動(dòng)流速與孔隙比也呈負(fù)指數(shù)遞減關(guān)系。

2.2 孔隙比和流速對(duì)沖蝕比的影響規(guī)律

土體的沖蝕發(fā)展規(guī)律與其內(nèi)部孔隙結(jié)構(gòu)和外部沖刷力密切相關(guān)，為研究庫岸土體的沖蝕演化進(jìn)程，本研究引入沖蝕比的概念：

ε=ΔMM（2）

式中：ε為沖蝕比;M和ΔM分別為土樣土顆粒總質(zhì)量和水流沖蝕導(dǎo)致的流失土顆粒質(zhì)量，g。

考慮到土樣底部與金屬底板的摩擦?xí)绊憶_蝕的演化進(jìn)程，所以每個(gè)試樣質(zhì)量沖蝕比達(dá)到95%（土樣剩余厚度1 mm等于土樣最大粒徑）即可認(rèn)為試驗(yàn)結(jié)束。不同流速下不同孔隙比土樣的沖蝕比隨時(shí)間發(fā)展曲線如圖4（a）所示，從各曲線的對(duì)比來看，當(dāng)水流速度大于起動(dòng)流速后，流速越低、孔隙比越小時(shí)，沖蝕結(jié)束所需時(shí)間呈幾何級(jí)數(shù)增加;從曲線的發(fā)展規(guī)律來看，有限厚度岸坡土體的沖蝕過程并非勻速發(fā)展，而是呈現(xiàn)增長速率逐漸變緩的規(guī)律，尤其在沖蝕比達(dá)到0.8～0.9之后，曲線呈現(xiàn)明顯的轉(zhuǎn)折，其本質(zhì)原因是由于土樣水平方向沖蝕深度不一、局部土樣已沖蝕見底、受水流作用的面積減小，且水槽底板對(duì)土顆粒的黏結(jié)力導(dǎo)致顆粒起動(dòng)和沖蝕難度急劇增大;沖蝕發(fā)展前期近乎直線，這也是以往研究中假設(shè)相同土體在相同流速下具有相同沖蝕速率的主要原因。

2.3 孔隙比和流速對(duì)沖蝕速率的影響規(guī)律

為反映土樣沖蝕比的增長速率，引入沖蝕速率的概念，如式（3）所示：

Et=ΔεΔt（3）

式中：Et為沖蝕速率，h-1;Δε為設(shè)定間隔時(shí)間Δt內(nèi)的沖蝕比差值，Δt越短，Et越精確。計(jì)算確定不同土樣在前1 h內(nèi)每間隔2～5 min的沖蝕速率變化，分析孔隙比和水流速度對(duì)Et的影響規(guī)律，如圖4（b）。

由圖4（b）可知：① 沖蝕前期階段，雖在圖4（a）中為近似直線，然而如果數(shù)據(jù)收集間隔足夠緊密時(shí)，會(huì)發(fā)現(xiàn)在短間隔內(nèi)的沖蝕速率也是逐漸衰減的，且流速越快、孔隙比越大土樣的沖蝕速率衰減越明顯。主要原因是飽和過程及起動(dòng)流速測定階段的長時(shí)間浸泡軟化使土樣表面相對(duì)松散導(dǎo)致沖蝕速率在試驗(yàn)初始階段最大，而隨著沖蝕發(fā)生試樣表面粗糙度增大不利于顆粒運(yùn)移從而造成沖蝕速率逐漸降低。② 沖蝕速率在下降過程中常有波動(dòng)而非單調(diào)降低，主要是由于黏土團(tuán)粒的遷移具有突然性和非均勻性，其與粗粒土的沖蝕發(fā)展規(guī)律呈現(xiàn)不同特征［20］。③ 相同流速下，沖蝕速率隨著孔隙比的增大而增加，亦即黏性土抗沖蝕能力隨著孔隙比的增大而減小?？紫侗仍叫。瑔挝惑w積土顆粒越多，對(duì)應(yīng)單位體積土顆粒排列越緊密，顆粒間的摩擦強(qiáng)度和黏土的黏結(jié)強(qiáng)度也就越大，抗沖蝕能力越強(qiáng)，這與已有研究［10，13，21］得出的增大土體干密度可提高臨界侵蝕應(yīng)力和抗蝕能力的結(jié)論是一致的。④ 相同孔隙比土體在各時(shí)刻沖蝕速率均隨著水流流速的增大而增加，主要是由于水流的剪切力與流速呈正相關(guān)關(guān)系。

進(jìn)一步，提取圖4（b）中部分?jǐn)?shù)據(jù)點(diǎn)作對(duì)比分析。圖5（a）對(duì)應(yīng)不同流速下孔隙比為1.1土樣在前1 h內(nèi)不同時(shí)段的平均沖蝕速率，發(fā)現(xiàn)在低流速作用下，土樣平均沖蝕速率與間隔長度的關(guān)系不明顯，即沖蝕呈勻速緩慢發(fā)展趨勢；高流速下隨著間隔時(shí)段長度的增加其平均沖蝕速率明顯減小。圖5（b）對(duì)應(yīng)流速為0.93 m/s作用下不同孔隙比土樣不同時(shí)段的平均沖蝕速率，可知較低孔隙比土樣的平均沖蝕速率在不同時(shí)段長度內(nèi)基本相同，即沖蝕呈勻速緩慢發(fā)展趨勢；而高孔隙比土樣的平均沖蝕速率在初始時(shí)刻較大，隨著間隔時(shí)段增長而顯著減小。

本研究考慮到現(xiàn)實(shí)岸坡土體并非是無限埋深的，試驗(yàn)中設(shè)置被沖蝕土樣為有限厚度2 cm，分析沖蝕的演化規(guī)律。如若增加或減少土樣厚度，沖蝕曲線仍會(huì)表現(xiàn)類似的減緩趨勢和轉(zhuǎn)折現(xiàn)象，但由于底部巖層黏結(jié)力的影響范圍有限，發(fā)生轉(zhuǎn)折時(shí)的沖蝕比可能會(huì)增大或減小。

3 粉質(zhì)黏土沖蝕預(yù)測模型

根據(jù)沖蝕進(jìn)程曲線演化規(guī)律，選取指數(shù)模型和雙曲線模型對(duì)實(shí)測數(shù)據(jù)擬合對(duì)比，并建立考慮水流速度和土樣孔隙比影響的有限厚度粉質(zhì)黏土的沖蝕預(yù)測模型。

3.1 沖蝕進(jìn)程擬合模型的建立

（1） 指數(shù)模型。實(shí)測的沖蝕比隨時(shí)間發(fā)展規(guī)律符合指數(shù)模型形式，借鑒Sterpi［22］內(nèi)侵蝕量的預(yù)測模型并修正后建立如下式：

ε=ΔMM=1-exp-tt0mnvv0（4）

式中：t為沖蝕時(shí)間，h;t0=1 h，為單位時(shí)間;m、n為量綱一參數(shù)。

對(duì)式（4）求導(dǎo)，得到?jīng)_蝕速率為

Et=dεdt=mtm-1tm0nvv0exp-tt0mnvv0（5）

（2） 雙曲線模型。de Oliveira等［23］認(rèn)為沖蝕比隨時(shí)間的發(fā)展也可用雙曲線模型近似表達(dá)：

ε=ta+bt（6）

雙曲線模型的切線模量代表任意時(shí)刻的沖蝕速率，表示為

Et=dεdt=a（a+bt）2（7）

在試驗(yàn)的起始點(diǎn)t→0，Et=E0，則E0=1/a，表明a代表該試驗(yàn)中初始沖蝕速率（E0）的倒數(shù)。

式（6）中，如果t→∞，則極限沖蝕比εult=1b。由此可知，b為雙曲線的漸近線所對(duì)應(yīng)極限沖蝕比的倒數(shù)，h-1，而本研究的有限厚度土體最終將全部沖蝕，即εult=1，故有b=1。

3.2 2種模型擬合結(jié)果對(duì)比

根據(jù)式（4）—式（7），采用2種模型分別對(duì)沖蝕進(jìn)程試驗(yàn)結(jié)果進(jìn)行擬合，僅以流速0.538 m/s為代表，對(duì)比如圖6，擬合參數(shù)見表2。由圖表可知2種模型的擬合度均較高，都能反映沖蝕比和沖蝕速率隨時(shí)間發(fā)展的基本趨勢，且指數(shù)模型的擬合效果相對(duì)更好。

進(jìn)一步，分析雙曲線模型擬合結(jié)果，發(fā)現(xiàn)沖蝕發(fā)生前期擬合結(jié)果普遍高于試驗(yàn)值，而沖蝕后期擬合結(jié)果低于試驗(yàn)值。參數(shù)a隨著孔隙比和流速的增大而減小，反映了初始沖蝕速率隨二者基本呈線性增大規(guī)律，這與de Oliveira等［23］試驗(yàn)成果趨勢基本一致。但由單一參數(shù)a同時(shí)反映2個(gè)因素的影響略顯簡單，且該模型沒有考慮臨界起動(dòng)流速的影響，物理意義不明確。

指數(shù)模型考慮了實(shí)際水流速度和土樣起動(dòng)流速的影響，可用于擬合不同土體的沖蝕演化規(guī)律。對(duì)試驗(yàn)曲線擬合時(shí)發(fā)現(xiàn)參數(shù)m值接近于常數(shù)1，為減少模型參數(shù)，設(shè)定m=1時(shí)依然得到較高的擬合度。而且，由表2可知，n值對(duì)流速變化的響應(yīng)并不明顯，僅是一個(gè)與孔隙比密切相關(guān)的變量。

根據(jù)表2，分析指數(shù)模型參數(shù)n隨孔隙比的變化規(guī)律，發(fā)現(xiàn)二者也符合指數(shù)型關(guān)系，可擬合得到參數(shù)n的預(yù)測公式為

n=0.001 3exp（3.665 7e）（8）

3.3 有限厚度土體沖蝕發(fā)展規(guī)律預(yù)測

將式（1）、式（8）帶入式（4），可得同時(shí)考慮孔隙比、起動(dòng)流速和水流流速的粉質(zhì)黏土沖蝕發(fā)展預(yù)測模型：

ε=1-exp-ηtt0vexp（δe）（9）

用參數(shù)η和δ替換原擬合參數(shù)a0、b0、m、n，得到η=0.002 5、δ=6.229 7，二者可視為土性參數(shù)。

應(yīng)用該模型，分別針對(duì)不同流速下或者不同孔隙比土體開展沖蝕發(fā)展進(jìn)程預(yù)測，得到?jīng)_蝕比及沖蝕速率隨水流沖刷時(shí)間的演化規(guī)律?？芍孩?沖蝕比和沖蝕速率預(yù)測曲線與試驗(yàn)曲線的發(fā)展趨勢相符，不同條件下土樣沖蝕比均隨時(shí)間經(jīng)歷了快速增加—緩慢增加—基本穩(wěn)定3個(gè)階段的發(fā)展趨勢，對(duì)應(yīng)的沖蝕速率均隨時(shí)間呈指數(shù)衰減并最終趨近于0。② 土體孔隙比和水流速度是影響沖蝕發(fā)展進(jìn)程的主要因素;相同流速下，土體孔隙比等間距減小時(shí)，完全沖蝕所需時(shí)間呈級(jí)數(shù)增大，而初始沖蝕速率呈幾何級(jí)數(shù)減小;相同孔隙比下，水流流速越大，導(dǎo)致土體初始沖蝕速率越大、沖蝕速率衰減更快、完全沖蝕所需時(shí)間更短。以孔隙比為0.96土樣為例，預(yù)測結(jié)果如圖7所示。

本研究中設(shè)置土樣為均質(zhì)的，如果現(xiàn)場勘察能夠得到不同深度上土體孔隙比與其埋深的關(guān)系，可在預(yù)測模型式（9）中考慮孔隙比的動(dòng)態(tài)變化，進(jìn)而預(yù)測岸坡的沖蝕發(fā)展規(guī)律。

4 結(jié)? 論

針對(duì)三峽庫岸有限厚度的粉質(zhì)黏土，開展了不同孔隙比土樣的沖刷起動(dòng)試驗(yàn)和不同流速下的沖蝕進(jìn)程試驗(yàn)，基于沖蝕發(fā)展規(guī)律建立了有限厚度土體沖蝕演化預(yù)測模型，得到以下結(jié)論：

（1） 粉質(zhì)黏土微團(tuán)臨界起動(dòng)流速與土樣孔隙比呈負(fù)指數(shù)遞減關(guān)系，有限厚度粉質(zhì)黏土完全沖蝕的耗時(shí)與土樣孔隙比或水流速度均呈反比。

（2） 受土體厚度有限的影響，沖蝕比達(dá)到0.8～0.9后隨時(shí)間發(fā)展的曲線發(fā)生明顯轉(zhuǎn)折，本質(zhì)原因是由于局部土樣已沖蝕見底、受水流作用的面積減小，且底板（類似于現(xiàn)實(shí)中的下部基巖）黏結(jié)力作用，導(dǎo)致顆粒起動(dòng)和沖蝕難度急劇增大。

（3） 沖蝕比小于0.8時(shí)，沖蝕速率隨時(shí)間發(fā)展呈負(fù)指數(shù)衰減規(guī)律，且初始階段的速率是沖蝕中期的3～5倍，其內(nèi)在原因是由于試樣密實(shí)程度隨深度增大、表面粗糙度隨顆粒沖蝕逐漸增加，粒團(tuán)分解和運(yùn)移難度增大;且孔隙比越大或流速越大時(shí)，初始沖蝕速率越大，后期衰減趨勢越明顯。

（4） 建立的指數(shù)預(yù)測模型考慮了土體臨界起動(dòng)流速、孔隙比、水流速度對(duì)沖蝕演化進(jìn)程的影響，物理意義明確，預(yù)測效果較好，僅有2個(gè)與土性相關(guān)的參數(shù)，可描述有限厚度土體的沖蝕發(fā)展規(guī)律，為水流沖蝕作用下岸坡形態(tài)變化趨勢預(yù)測研究提供了參考。

參考文獻(xiàn)：

［1］假冬冬，楊俊，陳長英，等.中國北方季節(jié)性冰凍河流岸灘崩塌數(shù)值模擬：以松花江為例［J］.水科學(xué)進(jìn)展，2021，32（5）：717-726.（JIA D D，YANG J，CHEN C Y，et al.Numerical simulation of bank erosion in a typical seasonally frozen river-case study of the Songhua River，China［J］.Advances in Water Science，2021，32（5）：717-726.（in Chinese））

［2］熊海濱，孫昭華，李明，等.清水沖刷條件下長江中游沙卵石河段局部卵石淤積成因［J］.水科學(xué)進(jìn)展，2020，31（4）：524-534.（XIONG H B，SUN Z H，LI M，et al.Causes of local gravel deposition in gravel-sand transition of the Middle Yangtze River under clear water scour［J］.Advances in Water Science，2020，31（4）：524-534.（in Chinese））

［3］KIMIAGHALAM N，GOHARROKHI M，CLARK S P，et al.A comprehensive fluvial geomorphology study of riverbank erosion on the Red River in Winnipeg，Manitoba，Canada［J］.Journal of Hydrology，2015，529：1488-1498.

［4］高曉靜，王秋生，李勇，等.重塑粗粒土-黏性土混合物沖刷特性研究［J］.水利學(xué)報(bào)，2021，52（3）：323-332.（GAO X J，WANG Q S，LI Y，et al.Study on the erodibility of reconstituted cohesive and non-cohesive soil mixtures［J］.Journal of Hydraulic Engineering，2021，52（3）：323-332.（in Chinese））

［5］鐘啟明，陳小康，梅勝堯，等.滑坡堰塞湖潰決風(fēng)險(xiǎn)與過程研究進(jìn)展［J］.水科學(xué)進(jìn)展，2022，33（4）：659-670.（ZHONG Q M，CHEN X K，MEI S Y，et al.A state of the art review on the failure risk and process of the landslide-induced dammed lake［J］.Advances in Water Science，2022，33（4）：659-670.（in Chinese））

［6］龔政，葛冉，馮騫，等.泥沙顆粒間黏結(jié)力作用及其對(duì)泥沙起動(dòng)影響研究進(jìn)展［J］.水科學(xué)進(jìn)展，2021，32（5）：801-812.（GONG Z，GE R，F(xiàn)ENG Q，et al.Cohesive forces between sediment particles and its impact on incipient motion of sediment：a review［J］.Advances in Water Science，2021，32（5）：801-812.（in Chinese））

［7］HE N，YANG L Q，XU L J，et al.Experimental study on the scouring rate of cohesive soil in the Lower Yellow River［J］.Coatings，2022，12（3）：418.

［8］BARMAN K，ROY S，DAS V K，et al.Effect of clay fraction on turbulence characteristics of flow near an eroded bank［J］.Journal of Hydrology，2019，571：87-102.

［9］ZHANG M X，YU G L，la ROVERE A，et al.Erodibility of fluidized cohesive sediments in unidirectional open flows［J］.Ocean Engineering，2017，130：523-530.

［10］KIMIAGHALAM N，CLARK S P，AHMARI H.An experimental study on the effects of physical，mechanical，and electrochemical properties of natural cohesive soils on critical shear stress and erosion rate［J］.International Journal of Sediment Research，2016，31（1）：1-15.

［11］WILSON B N.Development of a fundamentally based detachment model［J］.Transactions of the ASAE，1993，36（4）：1105-1114.

［12］HANSON G J，SIMON A.Erodibility of cohesive streambeds in the loess area of the Midwestern USA［J］.Hydrological Processes，2001，15（1）：23-38.

［13］JEPSEN R，ROBERTS J，LICK W.Effects of bulk density on sediment erosion rates［J］.Water，Air，and Soil Pollution，1997，99（1）：21-31.

［14］JACOBS W，le HIR P，van KESTEREN W，et al.Erosion threshold of sand-mud mixtures［J］.Continental Shelf Research，2011，31（10）：S14-S25.

［15］SANFORD L P，MAA J P Y.A unified erosion formulation for fine sediments［J］.Marine Geology，2001，179（1/2）：9-23.

［16］周銀軍，金中武，徐濤.三峽水庫庫區(qū)上段水流特性分析［J］.水力發(fā)電學(xué)報(bào)，2015，34（6）：131-137.（ZHOU Y J，JIN Z W，XU T.Flow structure in upper part of Three Gorges Reservoir［J］.Journal of Hydroelectric Engineering，2015，34（6）：131-137.（in Chinese））

［17］SHARMA I，THOMAS A.Experimental evaluation of internal erosion using groove erosion test［J］.International Journal of Geomechanics，2022，22（6）：04022056.

［18］李相輝，張慶松，張霄，等.一種均勻黏性土黏聚力計(jì)算方法［J］.工程科學(xué)與技術(shù)，2019，51（1）：137-143.（LI X H，ZHANG Q S，ZHANG X，et al.A method for calculating the cohesive force of uniform cohesion soil［J］.Advanced Engineering Sciences，2019，51（1）：137-143.（in Chinese））

［19］許琳娟，趙萬杰，李軍華，等.粘性土起動(dòng)流速影響因素試驗(yàn)研究［J］.水電能源科學(xué)，2020，38（3）：121-124，128.（XU L J，ZHAO W J，LI J H，et al.Experimental study on influencing factors of incipient velocity of cohesive soil［J］.Water Resources and Power，2020，38（3）：121-124，128.（in Chinese））

［20］HAM S M，KWON T H，CHANG I，et al.Ultrasonic P-wave reflection monitoring of soil erosion for erosion function apparatus［J］.Geotechnical Testing Journal，2016，39（2）：20150040.

［21］陳寶，田昌春，郭家興，等.地下水對(duì)壓實(shí)高廟子膨潤土沖蝕作用研究［J］.巖土力學(xué)，2016，37（11）：3224-3230.（CHEN B，TIAN C C，GUO J X，et al.Erosion of compacted Gaomiaozi bentonite by groundwater flow［J］.Rock and Soil Mechanics，2016，37（11）：3224-3230.（in Chinese））

［22］STERPI D.Effects of the erosion and transport of fine particles due to seepage flow［J］.International Journal of Geomechanics，2003，3（1）：111-122.

［23］de OLIVEIRA V N，de F N GITIRANA G Jr，dos ANJOS MASCARENHA M M，et al.An enhanced flume testing procedure for the study of rill erosion［J］.Water，2021，13（21）：2956.

Erosion evolution and prediction of silty clay with finite thickness

on reservoir bank

The study is financially supported by the National Natural Science Foundation of China （No.42272317;No.U21A2031）.

CHEN Yong1，ZHANG Shushi1，CHAN Dave1，2

（1. Key Laboratory of Geological Hazards on Three Gorges Reservoir Area of Ministry of Education，Three Gorges University，

Yichang 443002，China;

2. Department of Civil and Environmental Engineering，University of Alberta，

Edmonton，AB T6G 2W2，Canada）

Abstract：Soil erosion and particles loss are important factors for the development of collapse on river bank and changes in river morphology.In order to understand and predict the process of soil particles loss on river bank due to water flow，a new soil erosion test apparatus is developed that can accurately measure the quality change in soil samples in real time.Erosional tests are carried out to detect the initiation and development of erosion on a silty clay with finite thickness obtained from the bank of the Three Gorges Reservoir.A model is developed to capture the evolution of erosion and the mechanism of the erosion process.The results show that the total amount of erosion on finite thickness soils increases with time，but the erosion rate gradually decreases with time and approaches a limiting value.The rate of erosion is proportional to the void ratio，or density，of the soil and the flow velocity.The void ratio and flow velocity are incorporated in the model，reflecting the erodibility of soil and the shear forces applied on the soil′s surface，which captures the evolutional process of erosion of the silty clay.The model calculates the erosion on the finite thickness soils which can be used as a reference in determining the changes in reservoir bank morphology under the action of water flow.

Key words：soil erosion test;finite thickness;erosion ratio;erosion rate;prediction model