張幫穩(wěn) 鄧安軍 王黨偉 馮勝航 史紅玲 呂瑞茹

摘要：水庫(kù)蓄水后，水深增加使得河道干支流水動(dòng)力條件發(fā)生改變，干流泥沙倒灌造成支流淤積，影響水庫(kù)有效庫(kù)容和航運(yùn)安全。以溪洛渡水庫(kù)為例，基于實(shí)測(cè)資料和三維水沙數(shù)學(xué)模型，研究庫(kù)區(qū)干流和支流的水沙輸移規(guī)律，分析支流河口段泥沙淤積特性及攔門沙形成機(jī)理。結(jié)果表明：① 庫(kù)區(qū)支流河道淤積泥沙來(lái)自支流上游和干流倒灌，其中，支流上游來(lái)沙是庫(kù)區(qū)支流河道淤積的主要來(lái)源。② 干流倒灌庫(kù)區(qū)支流河口形成對(duì)流和環(huán)流區(qū)，泥沙易于淤積，西溪河和牛欄江坡度較緩，存在形成攔門沙的風(fēng)險(xiǎn)較大;美姑河坡度較大，存在形成攔門沙的風(fēng)險(xiǎn)較小。③ 對(duì)于西溪河、牛欄江和美姑河，干流倒灌的距離隨著流量的增加而增加，高水位運(yùn)行時(shí)，水位對(duì)干流倒灌的距離影響不明顯;而對(duì)于金陽(yáng)河和西蘇角河，干流倒灌的距離隨著流量的增加基本不變，干流倒灌的距離隨著水位的增加而增加。研究成果為新水沙情勢(shì)下庫(kù)區(qū)支流攔門沙風(fēng)險(xiǎn)防控提供依據(jù)。

關(guān)鍵詞：攔門沙;干流倒灌;對(duì)流區(qū);支流;溪洛渡水庫(kù)

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

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

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

收稿日期：2023-06-07;網(wǎng)絡(luò)出版日期：2023-10-31

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

基金項(xiàng)目：國(guó)家自然科學(xué)基金資助項(xiàng)目（U2040217;52209104）

作者簡(jiǎn)介：張幫穩(wěn)（1987—），男，山東菏澤人，高級(jí)工程師，博士，主要從事水庫(kù)泥沙治理及數(shù)值模擬方面研究。

E-mail：zhangbw@iwhr.com

通信作者：鄧安軍，E-mail：denganj@iwhr.com

水庫(kù)蓄水后，水深增加使得干支流水沙動(dòng)力特性發(fā)生改變，導(dǎo)致庫(kù)區(qū)泥沙不斷淤積，不僅關(guān)系到水庫(kù)防洪、發(fā)電、航運(yùn)等綜合效益的發(fā)揮，而且對(duì)水庫(kù)及下游生態(tài)環(huán)境具有重要的影響［1-3］。水庫(kù)干支流交匯處易形成攔門沙，阻礙水庫(kù)支流與干流水沙交匯，影響水庫(kù)的總庫(kù)容，在一定程度上削弱了水庫(kù)的調(diào)節(jié)能力［4-5］。

研究者對(duì)攔門沙形成因素進(jìn)行了分析，指出異重流、來(lái)水來(lái)沙、干支流交匯段河勢(shì)、水庫(kù)運(yùn)行方式等影響著攔門沙的形成［5-8］。胡春宏等［6］指出黃河口來(lái)水來(lái)沙過(guò)程變異對(duì)黃河口攔門沙的發(fā)育演變有著重要的影響;王遠(yuǎn)見等［7］指出小浪底水庫(kù)的出庫(kù)流量和庫(kù)區(qū)淤積形態(tài)是影響攔門沙形成的2個(gè)最重要因素;王婷等［8］指出干流異重流倒灌導(dǎo)致小浪底支流泥沙淤積集中在河口附近;朱玲玲等［5］指出三峽水庫(kù)支流河口段淤積的泥沙以干流倒灌為主，也存在支流來(lái)沙淤積，隨著干支流來(lái)沙量趨于減少，支流河口淤積形成攔門沙坎的風(fēng)險(xiǎn)較小。水庫(kù)蓄水使得庫(kù)區(qū)水動(dòng)力條件減弱是庫(kù)區(qū)干支流泥沙淤積的根本原因。庫(kù)區(qū)不同干支流交匯處水流、泥沙、地形相互作用是一個(gè)非常復(fù)雜的水動(dòng)力現(xiàn)象，特別是在河漫灘、彎曲河段和沙波地形情況下會(huì)產(chǎn)生二次流和流體分離等三維水動(dòng)力特性，支流攔門沙的形成往往是多種因素共同作用的結(jié)果。隨著計(jì)算機(jī)科學(xué)技術(shù)和高性能計(jì)算機(jī)的發(fā)展，三維數(shù)值模型逐漸應(yīng)用到工程實(shí)際當(dāng)中，已有研究者采用三維水沙數(shù)值模型研究了三峽大壩附近的流場(chǎng)、懸移質(zhì)輸移及壩前庫(kù)區(qū)的泥沙淤積問題［9-11］。目前，對(duì)攔門沙形成機(jī)理多以實(shí)測(cè)資料分析為主，數(shù)值模型研究的較少，缺乏較為深入的認(rèn)識(shí)。

2012年向家壩、2013年溪洛渡水庫(kù)相繼蓄水運(yùn)行，據(jù)資料分析［12］，金沙江干支流及區(qū)間的來(lái)沙基本被這2座電站攔截，且大部分泥沙淤積在溪洛渡庫(kù)區(qū)。本文以溪洛渡水庫(kù)為例，利用三維水沙數(shù)學(xué)模型對(duì)汛期洪水傳播和泥沙輸移過(guò)程進(jìn)行數(shù)值模擬，分析庫(kù)區(qū)干支流水沙輸移特性，探討支流淤積泥沙來(lái)源及攔門沙形成的風(fēng)險(xiǎn)，給出不同水位不同流量下干流泥沙倒灌支流的距離。研究成果以期為揭示新水沙情況下庫(kù)區(qū)干支流淤積規(guī)律和攔門沙的形成風(fēng)險(xiǎn)提供依據(jù)。

1 數(shù)學(xué)模擬方法

1.1 水流控制方程

三維水動(dòng)力學(xué)模型的控制方程［13-14］采用基于Reynolds的時(shí)均N-S方程，滿足靜壓假定和Boussinesq渦黏性假定。在笛卡爾坐標(biāo)系下，不可壓縮流體N-S方程的連續(xù)性方程為

ux+vy+wz=0（1）

動(dòng)量守恒方程為：

dudt=fu-gηx-1ρ0pax-gρ0∫ηzbρxdz+zkmvuz+kmh2ux2+2uy2（2）

dvdt=fv-gηx-1ρ0pay-gρ0∫ηzbρydz+zkmvvz+kmh2vx2+2vy2（3）

式中：x和y分別為水平坐標(biāo);z為垂向坐標(biāo)，向上為正;u、v、w分別為x、y、z方向的流速;t為時(shí)間;f為柯氏力系數(shù);η為自由水面;zb為河床底高程;ρ0和ρ分別為參考密度和混合流體的密度;g為重力加速度;kmh和kmv分別為水平與垂直向渦黏性系數(shù)，其中，垂向渦黏性系數(shù)基于紊流Generic Length Scale（GLS）模型［15］進(jìn)行封閉求解，水平渦黏性系數(shù)采用常數(shù)化處理;pa為自由水面大氣壓強(qiáng)。

1.2 懸移質(zhì)泥沙輸移方程

懸移質(zhì)泥沙輸運(yùn)方程基于對(duì)流擴(kuò)散理論，將泥沙視為連續(xù)性介質(zhì)，并假設(shè)懸移質(zhì)泥沙顆粒和水體運(yùn)動(dòng)在垂直方向上存在速度差，且等于泥沙顆粒的沉速。懸移質(zhì)泥沙輸運(yùn)方程［9］為

Ct+uCx+vCy+（w+ω）Cz=zKmvCz+Kmh2Cx2+2Cy2+Qs（4）

式中：C為泥沙顆粒濃度;ω為泥沙顆粒的沉降速度;Kmv為泥沙顆粒的垂向擴(kuò)散系數(shù)，通常假定與水流紊動(dòng)黏性系數(shù)呈倍數(shù)關(guān)系，可通過(guò)紊流閉合模型求解或采用經(jīng)驗(yàn)關(guān)系估計(jì)：Ksv=Kmv/σs，σs為Schmidt數(shù)，取值通常在0.6～1.2之間;Kmh為泥沙顆粒的水平向擴(kuò)散系數(shù)，泥沙顆粒在水平上擴(kuò)散遠(yuǎn)小于垂向擴(kuò)散，通常忽略不計(jì);Qs為懸移質(zhì)和床沙的泥沙交換量。在本研究中，由于金沙江下游白鶴灘水庫(kù)建設(shè)，粗顆粒泥沙基本淤積在白鶴灘庫(kù)區(qū)，故模型中沒有考慮推移質(zhì)泥沙運(yùn)動(dòng)。

1.3 河床沖刷演變方程

在水體紊流作用下，懸移質(zhì)和床面泥沙在垂向上進(jìn)行相互交換。泥沙交換過(guò)程可看作床面附近泥沙顆粒在垂向上泥沙通量的交換［16］，根據(jù)泥沙質(zhì)量守恒，懸移質(zhì)和床面泥沙之間交換表達(dá)式為：

Qs=Ds+Es（5）

Ds=-ωCbz（6）

Es=E0（1-p）（b/cr-1）?? b>cr（7）

式中：Ds和Es分別為沉降通量和沖刷通量;Cb為數(shù)值模型中最底部一層網(wǎng)格的泥沙濃度;E0為經(jīng)驗(yàn)沖刷率系數(shù)，大小根據(jù)床面泥沙顆粒特性確定;p為床面泥沙層泥沙顆粒的體積分?jǐn)?shù);cr為泥沙顆粒臨界起動(dòng)剪切應(yīng)力;b為水流底部的剪切應(yīng)力。

懸移質(zhì)和床面泥沙交換及凈輸移導(dǎo)致河床表面地形的變化，床面地形變化的公式為

Δh=（Ds-Es）Δtρs（1-p）（8）

式中：Δh為泥沙床面高程的變化量。

2 數(shù)學(xué)模型的建立及驗(yàn)證

2.1 數(shù)值模型的邊界條件

溪洛渡庫(kù)區(qū)長(zhǎng)約190 km，庫(kù)區(qū)內(nèi)支流主要包括西溪河、牛欄江、金陽(yáng)河、美姑河和西蘇角河;水文站主要有白鶴灘、昭覺、大沙店和美姑，壩前水位站為黃角堡站，見圖1。在數(shù)值模擬研究中，由于缺少金陽(yáng)河和西蘇角河的實(shí)測(cè)數(shù)據(jù)，本文中作為閉合邊界考慮，入口邊界條件為白鶴灘、昭覺、大沙店、美姑4個(gè)水文站實(shí)測(cè)流量和含沙量。需要指出的是，支流昭覺、大沙店、美姑3站缺少2012年以后實(shí)測(cè)含沙量數(shù)據(jù)，通過(guò)對(duì)比分析歷年白鶴灘和各支流7月份流量和含沙量的大小，在本文研究中采用2018年7月白鶴灘站和2011年7月支流的流量和含沙量實(shí)測(cè)數(shù)據(jù)作為支流入流邊界條件，見圖2。出口邊界條件為黃角堡的水位，見圖3。地形采用2016年實(shí)測(cè)散點(diǎn)地形，通過(guò)網(wǎng)格進(jìn)行插值，見圖1。水平向網(wǎng)格采用三角形網(wǎng)格，網(wǎng)格大小在變動(dòng)回水區(qū)和支流為16 m，在干流常年回水區(qū)為20 m，垂向上采用LSC2（Localized Sigma Coordinates with Shaved Cells）網(wǎng)格［17］，15 m為1層，平均層數(shù)為19層。

2.2 數(shù)學(xué)模型的驗(yàn)證

模型模擬了汛期溪洛渡庫(kù)區(qū)水沙的傳播過(guò)程，模型在前期研究溪洛渡水庫(kù)洪峰和沙峰異步傳播規(guī)律過(guò)程中進(jìn)行了驗(yàn)證［18］，三維數(shù)值模擬效果見圖4。主要對(duì)水平網(wǎng)格大小、垂向網(wǎng)格大小及壩前邊界條件的設(shè)置進(jìn)行了結(jié)果的敏感性分析，這里僅給出壩前流量和含沙量的模擬值和實(shí)測(cè)值的對(duì)比結(jié)果，見圖5，詳細(xì)的模型驗(yàn)證見參考文獻(xiàn)［18］。

3 溪洛渡庫(kù)區(qū)干支流水沙輸移規(guī)律及淤積分析

3.1 干流水沙輸移規(guī)律及淤積分析

圖6給出了溪洛渡庫(kù)區(qū)干流不同時(shí)刻深泓線處流速和泥沙濃度的瞬時(shí)分布，從圖6（a）中可以看出，上游庫(kù)區(qū)水深較小，流速較大，壩前常年回水區(qū)水深較大，流速較小;從圖6（b）可以看到，上游庫(kù)區(qū)泥沙濃度較大，隨著洪水向下游庫(kù)區(qū)傳播，水深增加，流速減小，水流挾沙能力逐漸減小，泥沙逐漸落淤，壩前常年回水區(qū)泥沙濃度較小。2018年8月5—31日入庫(kù)泥沙量為2 500萬(wàn)t，出庫(kù)泥沙量為700萬(wàn)t，通過(guò)分析對(duì)坪鎮(zhèn)的流量和含沙量，其中，對(duì)坪鎮(zhèn)以上泥沙淤積量為4 300萬(wàn)t，占庫(kù)區(qū)淤積量的17.7%，可以得出在場(chǎng)次洪水傳播過(guò)程中，泥沙主要淤積在常年回水區(qū)。

3.2 西溪河水沙輸移規(guī)律及淤積分析

西溪河位于庫(kù)區(qū)變動(dòng)回水區(qū)，受水庫(kù)運(yùn)行水位影響較大。圖7給出了西溪河不同時(shí)刻深泓線處流速和含沙量的瞬時(shí)分布。從圖7（a）和圖2中可以看出，庫(kù)區(qū)干流流量較西溪河大，庫(kù)區(qū)干流流速較大，在交匯處存在庫(kù)區(qū)干流倒灌西溪河的現(xiàn)象。從圖7（b）可以看出，在7月15日之前，溪洛渡水庫(kù)運(yùn)行水位較低，隨著庫(kù)區(qū)干流水流倒灌西溪河，泥沙隨水流倒灌西溪河河口附近;7月15日之后，溪洛渡水庫(kù)運(yùn)行水位較高，西溪河來(lái)水來(lái)沙量增加，泥沙可以從源頭流向河口附近;西溪河河口仍存在庫(kù)區(qū)干流泥沙的倒灌，但是倒灌的距離有所減小。因此，對(duì)于支流西溪河的泥沙淤積，上游來(lái)沙仍然是河道淤積的主要來(lái)源;庫(kù)區(qū)干流倒灌攜帶的泥沙主要淤積在河口附近，同時(shí)西溪河坡度較小，形成攔門沙的風(fēng)險(xiǎn)較大。

3.3 牛欄江水沙輸移規(guī)律及淤積分析

牛欄江位于庫(kù)區(qū)變動(dòng)回水區(qū)和常年回水區(qū)的之間，圖8給出了牛欄江不同時(shí)刻深泓線處流速和含沙量的瞬時(shí)分布。從圖8（a）中可以看到，庫(kù)區(qū)干流流量較牛欄江大，庫(kù)區(qū)干流流速較大，在交匯處存在庫(kù)區(qū)干流倒灌牛欄江的現(xiàn)象。從圖8（b）中可以看到，在牛欄江和干流交匯處，存在泥沙倒灌現(xiàn)象;牛欄江來(lái)沙隨著水流從上游逐漸向下游輸移，輸移的距離隨著庫(kù)區(qū)運(yùn)行水位的增加而減小。因此，對(duì)于支流牛欄江的泥沙淤積，上游來(lái)沙仍然是河道淤積的主要來(lái)源;庫(kù)區(qū)干流倒灌攜帶的泥沙主要淤積在河口附近，同時(shí)牛欄江坡度較小，形成攔門沙的風(fēng)險(xiǎn)較大。

3.4 美姑河水沙輸移規(guī)律及淤積分析

美姑河位于常年回水區(qū)，受庫(kù)區(qū)運(yùn)行水位影響較小。圖9給出了美姑河不同時(shí)刻深泓線處流速和含沙量的瞬時(shí)分布。從圖9（a）中可以看到，庫(kù)區(qū)干流流速較美姑河河口處流速大，在交匯處存在庫(kù)區(qū)干流倒灌牛欄江的現(xiàn)象;美姑河支流流速隨著來(lái)流量的增加而增大。從圖9（b）中可以看到，在美姑河和干流交匯處，也存在著庫(kù)區(qū)干流泥沙倒灌美姑河的現(xiàn)象;美姑河泥沙隨著水流從上游逐漸向下游輸移，輸移的距離隨著來(lái)流量的增加而增加。因此，對(duì)于支流美姑河的泥沙淤積，上游來(lái)沙仍然是河道淤積的主要來(lái)源;庫(kù)區(qū)干流倒灌攜帶的泥沙主要淤積在河口附近，但是干流庫(kù)區(qū)倒灌水流含沙量較小，總體的淤積量不大，加之美姑河坡度較大，形成攔門沙的風(fēng)險(xiǎn)較小。

3.5 溪洛渡庫(kù)區(qū)支流攔門沙形成機(jī)理分析

圖10 給出了溪洛渡庫(kù)區(qū)干流與西溪河、牛欄江和美姑河支流交匯處瞬時(shí)三維流速場(chǎng)分布，可以看出，干流的流速明顯大于支流流速，與支流的來(lái)流方向相反，干流在支流口向支流倒灌，在支流口上端形成一個(gè)對(duì)流區(qū)，流速減小，易于泥沙落淤，同時(shí)在干流與支流交匯口形成一個(gè)環(huán)流區(qū)，如圖11給出了水庫(kù)支流攔門沙形成機(jī)理示意。

4 溪洛渡庫(kù)區(qū)干流倒灌支流距離的影響因素分析

溪洛渡水庫(kù)汛期防洪控制水位為560 m、正常運(yùn)行水位為600 m。溪洛渡水庫(kù)干流泥沙倒灌與干支流流速的相對(duì)大小有關(guān)，干流的流速與庫(kù)水位和入庫(kù)流量有關(guān)。本文通過(guò)設(shè)置不同的研究方案，分析了不同庫(kù)水位、不同入庫(kù)流量對(duì)干流泥沙倒灌支流距離的影響。表1給出了不同入庫(kù)流量、不同庫(kù)水位下干流泥沙倒灌支流的距離。

從表1的數(shù)據(jù)分析可知，對(duì)于西溪河、牛欄江和美姑河，干流倒灌的距離隨著干流流量的增加而增加，高水位運(yùn)行時(shí)，水位對(duì)干流倒灌的距離影響不明顯;對(duì)于金陽(yáng)河和西蘇角河，干流倒灌的距離隨著流量的增加基本不變，干流倒灌的距離隨著水位的增加而增加，可能與干流倒灌支流的流向和干流流向的交匯角有關(guān)，金陽(yáng)河和西蘇角河與干流的流向交匯角小于90°，順著主流的方向。

5 結(jié)? 論

本文利用三維水沙數(shù)學(xué)模型對(duì)汛期溪洛渡水庫(kù)洪水傳播和泥沙輸移過(guò)程進(jìn)行數(shù)值模擬，分析庫(kù)區(qū)干支流水沙輸移特性，主要探討支流攔門沙形成的風(fēng)險(xiǎn)及干流倒灌支流距離的影響因素，主要結(jié)論如下：

（1） 溪洛渡庫(kù)區(qū)干流流量較大，支流流量較小，干流流速較支流流速大，存在干流倒灌支流的現(xiàn)象;干流倒灌庫(kù)區(qū)支流時(shí)在河口形成對(duì)流和環(huán)流區(qū)，泥沙易于淤積，形成攔門沙。

（2） 溪洛渡庫(kù)區(qū)支流河道淤積泥沙來(lái)自支流上游來(lái)沙和干流倒灌泥沙，支流上游來(lái)沙是支流河道淤積的主要來(lái)源;干流較高濃度泥沙倒灌支流，在支流河口附近發(fā)生淤積，西溪河和牛欄江坡度較緩，有利于干流倒灌，形成攔門沙的風(fēng)險(xiǎn)較大;倒灌美姑河河口水流含沙量較小，總體的淤積量不大，形成攔門沙的風(fēng)險(xiǎn)較小。

（3） 溪洛渡庫(kù)區(qū)西溪河、牛欄江和美姑河，干流倒灌的距離隨著流量的增加而增加，高水位運(yùn)行時(shí)，水位對(duì)干流倒灌的距離影響不明顯;對(duì)于金陽(yáng)河和西蘇角河，干流倒灌支流流向和干流流向一致，干流倒灌的距離隨著流量的增加基本不變，隨著水位的增加而增加。

參考文獻(xiàn)：

［1］HUANG L，LI X，F(xiàn)ANG H W，et al.Balancing social，economic and ecological benefits of reservoir operation during the flood season：a case study of the Three Gorges Project，China［J］.Journal of Hydrology，2019，572：422-434.

［2］WU P，WANG N R，ZHU L J，et al.Spatial-temporal distribution of sediment phosphorus with sediment transport in the Three Gorges Reservoir［J］.Science of the Total Environment，2021，769：144986.

［3］TANG X Y，TONG S C，HUANG G X，et al.Characteristics of sedimentation and channel adjustment linked to the Three Gorges Reservoir［J］.International Journal of Sediment Research，2021，36（2）：177-189.

［4］王延貴，胡春宏.官?gòu)d水庫(kù)淤積特點(diǎn)及攔門沙整治措施［J］.泥沙研究，2003（6）：25-30.（WANG Y G，HU C H.Study on sedimentation and the mouth bar control in Guanting Reservoir［J］.Journal of Sediment Research，2003（6）：25-30.（in Chinese））

［5］朱玲玲，許全喜，張歐陽(yáng)，等.三峽水庫(kù)支流河口淤積及攔門沙形成風(fēng)險(xiǎn)研究［J］.中國(guó)科學(xué)：技術(shù)科學(xué)，2019，49（5）：552-564.（ZHU L L，XU Q X，ZHANG O Y，et al.Sedimentation at estuary of 66 tributaries in the Three Gorges Reservoir［J］.Scientia Sinica Technologica，2019，49（5）：552-564.（in Chinese））

［6］胡春宏，張治昊.水沙過(guò)程變異條件下黃河口攔門沙的演變響應(yīng)與調(diào)控［J］.水利學(xué)報(bào)，2006，37（5）：511-517，524.（HU C H，ZHANG Z H.Evolution and control of estuary sandbar in Yellow River responding to variation of flow-sediment condition［J］.Journal of Hydraulic Engineering，2006，37（5）：511-517，524.（in Chinese））

［7］王遠(yuǎn)見，謝蔚，王婷，等.小浪底水庫(kù)庫(kù)區(qū)支流攔門沙形成的主因分析［J］.泥沙研究，2016（6）：51-58.（WANG Y J，XIE W，WANG T，et al.Key factors influencing evolution of sand bars in tributaries of the Xiaolangdi Reservoir［J］.Journal of Sediment Research，2016（6）：51-58.（in Chinese））

［8］王婷，陳書奎，馬懷寶，等.小浪底水庫(kù)1999—2009年泥沙淤積分布特點(diǎn)［J］.泥沙研究，2011（5）：60-66.（WANG T，CHEN S K，MA H B，et al.Distribution of deposition in Xiaolangdi Reservoir［J］.Journal of Sediment Research，2011（5）：60-66.（in Chinese））

［9］張幫穩(wěn)，吳保生，章若茵.三峽庫(kù)區(qū)汛期洪峰和沙峰異步運(yùn)動(dòng)特性的三維數(shù)值模擬［J］.水科學(xué)進(jìn)展，2021，32（3）：408-417.（ZHANG B W，WU B S，ZHANG R Y.Three-dimensional numerical of the asynchronous movement characteristics of flood and sediment peaks in the Three Gorges Reservoir［J］.Advances in Water Science，2021，32（3）：408-417.（in Chinese））

［10］FANG H W，RODI W.Three-dimensional calculations of flow and suspended sediment transport in the neighborhood of the dam for the Three Gorges Project （TGP） reservoir in the Yangtze River［J］.Journal of Hydraulic Research，2003，41（4）：379-394.

［11］JIA D D，SHAO X J，ZHANG X N，et al.Sedimentation patterns of fine-grained particles in the dam area of the Three Gorges Project：3-D numerical simulation［J］.Journal of Hydraulic Engineering，2013，139（6）：669-674.

［12］朱玲玲，許全喜，董炳江，等.金沙江下游溪洛渡水庫(kù)排沙效果及影響因素［J］.水科學(xué)進(jìn)展，2021，32（4）：544-555.（ZHU L L，XU Q X，DONG B J，et al.Study on the effect and influencing factors of sand discharge of Xiluodu Reservoir in the Lower Jinsha River［J］.Advances in Water Science，2021，32（4）：544-555.（in Chinese））

［13］ZHANG Y L，BAPTISTA A M.SELFE：a semi-implicit Eulerian-Lagrangian finite-element model for cross-scale ocean circulation［J］.Ocean Modelling，2008，21（3/4）：71-96.

［14］ZHANG Y J，YE F，STANEV E V，et al.Seamless cross-scale modeling with SCHISM［J］.Ocean Modelling，2016，102：64-81.

［15］UMLAUF L，BURCHARD H.A generic length-scale equation for geophysical turbulence models［J］.Journal of Marine Research，2003，61（2）：235-265.

［16］WARNER J C，SHERWOOD C R，SIGNELL R P，et al.Development of a three-dimensional，regional，coupled wave，current，and sediment-transport model［J］.Computers & Geosciences，2008，34（10）：1284-1306.

［17］ZHANG Y J，ATELJEVICH E，YU H C，et al.A new vertical coordinate system for a 3-D unstructured-grid model［J］.Ocean Modelling，2015，85：16-31.

［18］ZHANG B W，DENG A J，DONG X Y，et al.Numerical analysis of propagation characteristics of sediment peak and its implications for the Xiluodu Reservoir［J］.Water Supply，2022，22（10）：7490-7512.

Study on formation mechanism of barrier sandbar of tributaries

in Xiluodu Reservoir

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

ZHANG Bangwen1，2，DENG Anjun1，2，WANG Dangwei1，2，F(xiàn)ENG Shenghang3，SHI Hongling1，2，LYU Ruiru1，2

（1. State Key Laboratory of Watershed Water Cycle Simulation and Regulation，China Institute of Water Resources and

Hydropower Research，Beijing 100038，China;

2. Key Laboratory of Water and Sediment Science and River Management，

Ministry of Water Resources，China Research Institute of Water Resources and Hydropower，Beijing 100038，China;

3. China Three Gorges Corporation，Wuhan 430010，China）

Abstract：After the reservoir is impounded，the increase of water depth changes the hydrodynamic conditions of the river channel，the backfilling of the incoming sediment from main stream to the tributaries resulted the sedimentation of tributary mouth，affecting the effective storage capacity of the reservoir and the safety of navigation.In this paper，taking Xiluodu Reservoir as an example，the characteristics of water and sediment transport in the main stream and tributaries was studied，and the formation mechanism of barrier sandbar at the tributary mouth was analyzed based on the measured data and the three-dimensional mathematical model.The results is indicated that：the deposited sediment of tributaries mainly comes from the upstream incoming sediment and the backfilling of the main stream sediment，and the upstream incoming sediment is the main source of the sedimentation of tributaries;the main stream backflows into tributaries near the tributary mouth，and a convective and circulation zone are formed，where the sediment is easy to be deposited.The slope of Xixi River and Niulan River is slow，and the risk of the barrier sandbar formation is large，while the slope of Meigu River is large，and the risk of the barrier sandbar formation is small;for Xixi River，Niulan River and Meigu River，the distance of main stream backfilling into the tributary increases with the increase of discharge，and the influence of water level on the distance of main stream backfilling into the tributary is not obvious when the Xiluodu Reservoir is operating at high water level.In Jinyang River and Xisujiao River，the distance of main flow backfilling into the tributary is basically unchanged with the increase of discharge，and the distance of main stream backfilling into the tributary increases with the increase of water level.The research results provide a basis for controlling the risk of barrier sandbar formation at the tributary mouth of reservoir under the new condition of water and sediment.

Key words：barrier sandbar;mainstream backfilling;convective zone;tributary;Xiluodu Reservoir