三峽水庫場次洪水排沙規(guī)律及影響因素

高宇 任實 王海 呂超楠 趙汗青 張成瀟

摘要：水庫泥沙淤積導(dǎo)致庫容損失，影響水庫綜合效益發(fā)揮，研究水庫排沙規(guī)律及影響因素對于改善水庫淤積狀況具有重要意義。根據(jù)2003—2021年朱沱、北碚、武隆、寸灘等水文站實測數(shù)據(jù)，合理劃分三峽水庫場次洪水過程，分析其排沙規(guī)律，并重點研究壩前水位、入出庫輸沙率、洪水滯留時間等因素對排沙的影響。結(jié)果表明：三峽水庫場次洪水過程的劃分采用統(tǒng)一的定量標(biāo)準(zhǔn)并考慮洪峰沙峰異步傳播的影響更加合理;水庫排沙主要集中在汛期場次洪水，各年汛期場次洪水排沙總量平均占全年的67.0%;場次洪水排沙比與壩前水位成反比、與入出庫輸沙率成正比，但相關(guān)性均小于與其成反比的洪水滯留時間，且當(dāng)入庫平均流量大于30 000 m3/s或壩前平均水位低于150 m時，排沙比與各因素之間的相關(guān)性均有所提升。研究成果可為優(yōu)化三峽水庫運行方式、減少水庫泥沙淤積等提供重要參考。

關(guān)鍵詞：排沙比;場次洪水;洪峰沙峰異步;洪水滯留時間;三峽水庫

中圖分類號：TV147??文獻(xiàn)標(biāo)志碼：A??文章編號：1001-6791（2024）01-0112-11

隨著人類對水資源利用需求的不斷加大，越來越多的水庫在全球范圍內(nèi)建成運用，發(fā)揮著防洪、發(fā)電、灌溉、航運等功能［1］。然而，水庫蓄水導(dǎo)致壩前水位抬高、流速減緩，水流中挾帶的泥沙易于落淤，造成水庫面臨泥沙淤塞、庫容損失等問題［2-4］，嚴(yán)重時還會影響水庫發(fā)電、航運等效益的發(fā)揮［5-6］。因此，如何減少水庫淤積成為國內(nèi)外學(xué)者廣泛關(guān)注的問題，水庫調(diào)度排沙是減少水庫淤積、延長水庫使用壽命的重要手段，研究水庫排沙規(guī)律及影響因素具有重要意義。

國外一般采用攔沙率來反映水庫攔沙效果，如Sisinggih等［7］、Jothiprakash等［8］分別探討了Sengguruh水庫、Gobindsagar水庫的攔沙率變化;而國內(nèi)一般采用排沙比來反映水庫排沙效果，如朱玲玲等［9］、張倩倩等［10］分別探討了溪洛渡水庫、三峽水庫的排沙比變化。雖然攔沙率與排沙比概念不同，但均能表征水庫的淤積程度，對于同一座水庫，攔沙率與排沙比之和為1，且不同研究普遍認(rèn)為攔沙率或排沙比主要與進(jìn)出庫流量、水庫庫容及壩前水位等相關(guān)。三峽工程作為治理和開發(fā)長江的關(guān)鍵性骨干工程，泥沙問題一直以來都深受關(guān)注［11-12］。三峽水庫蓄水運行以來，汛期入庫沙量占全年的比例達(dá)到了90%以上［13-14］，且汛期來沙又進(jìn)一步集中于短歷時的洪水過程［15-16］，故諸多學(xué)者選擇以汛期的短歷時洪水過程為研究對象來深入探究三峽水庫的排沙規(guī)律。王玉璇等［17］選取了三峽水庫2003—2018年25場次洪水過程，董占地等［18］統(tǒng)計了三峽水庫2004—2012年汛期月度排沙比，黃仁勇等［19］選取了三峽水庫2003—2011年34場次洪水過程，普遍認(rèn)為三峽水庫排沙主要與進(jìn)出庫水沙條件、壩前水位、河道形態(tài)等因素有關(guān)，且各因素之間存在交互效應(yīng)，共同影響水庫排沙。然而，上述研究對于三峽水庫汛期場次洪水過程的選取多由人為主觀判斷，經(jīng)驗性較強(qiáng)，導(dǎo)致不同研究選取的場次洪水過程存在較大差異［17，19］，且三峽水庫蓄水后，庫區(qū)內(nèi)場次洪水過程洪峰沙峰異步傳播的特性顯著［13，20］，已有研究缺乏對這一因素的考慮，故當(dāng)前對三峽水庫汛期場次洪水排沙規(guī)律的認(rèn)識仍有待于深入探討，亟需采用統(tǒng)一的方法來劃分場次洪水過程。

本文根據(jù)實測水文泥沙數(shù)據(jù)，考慮場次洪水洪峰沙峰異步傳播的影響，采用定量標(biāo)準(zhǔn)劃分三峽水庫汛期場次洪水過程，重點分析場次洪水的水庫排沙特性，探討場次洪水排沙比與進(jìn)出庫水沙條件、水庫運行條件等之間的關(guān)系。研究成果可為三峽水庫進(jìn)一步增加調(diào)度靈活性、改善庫區(qū)泥沙淤積狀況、優(yōu)化水庫調(diào)度運行方式等提供重要參考，對于保障三峽水庫綜合效益發(fā)揮具有重要意義。

1 研究區(qū)域及資料

三峽大壩位于長江上游與中游分界處的宜昌市，水庫跨湖北、重慶兩省市，控制流域面積約100萬km2。三峽水庫水文站點分布如圖1所示，江津至涪陵河段為變動回水區(qū)，長約173 km，涪陵至大壩河段為常年回水區(qū)，長約487 km。北碚水文站、武隆水文站則分別是庫區(qū)主要入?yún)R支流嘉陵江、烏江的控制站。

本文收集了三峽水庫蓄水運行以來2003—2021年朱沱、北碚、武隆、寸灘水文站的實測日均流量、懸移質(zhì)含沙量數(shù)據(jù)來反映三峽水庫的入庫水沙過程;采用2003—2021年茅坪水位站的實測日均水位數(shù)據(jù)反映壩前水位變化過程;采用2003—2021年黃陵廟水文站的實測日均流量、懸移質(zhì)含沙量數(shù)據(jù)反映出庫水沙過程，各站點分布見圖1。

2 三峽水庫蓄水以來場次洪水劃分

三峽水庫蓄水運行以來，入庫泥沙在汛期場次洪水過程集中輸移的特征顯著，且三峽庫區(qū)洪峰沙峰異步傳播的特性更加明顯。為了準(zhǔn)確反映出三峽水庫場次洪水排沙規(guī)律，本文采用統(tǒng)一的定量標(biāo)準(zhǔn)來劃分蓄水以來的場次洪水過程?；谝延邢嚓P(guān)研究方法［21-23］，采用超定量取樣法篩選洪水峰值，對峰值獨立性進(jìn)行判別，并對場次洪水的起止時間進(jìn)行判定，該方法能夠準(zhǔn)確地選取場次洪水過程，并有效減少洪水特征分析中的不確定性。在此基礎(chǔ)上考慮洪峰沙峰異步傳播的影響，進(jìn)一步確定場次洪水的入出庫水沙過程。此外，考慮到長江上游干支流水文站洪水疊加存在時間不同步的問題，僅選取長江干流代表水文站更能反映三峽水庫的入庫洪水過程［15-16，24］，且烏江武隆站的入庫徑流占比不大［25］，故本文主要選取寸灘站的流量過程來反映三峽水庫汛期場次洪水的入庫洪水過程，入庫沙量則采用朱沱、北碚、武隆3站輸沙量之和（下同）。三峽水庫場次洪水劃分的步驟如下。

2.1 場次洪水峰值的確定

場次洪水峰值的確定采用超定量取樣法，即從實測日均流量序列中提取超過某一閾值的樣本作為洪峰流量。閾值取所有樣本的某一百分位數(shù)Lk，即k%的樣本小于或等于閾值Lk，對于場次洪水洪峰流量閾值取k=90。由于洪峰閾值的確定會很大程度影響洪峰樣本的選擇，本文在此基礎(chǔ)上進(jìn)一步采用經(jīng)驗值來篩選洪峰，主要考慮寸灘站洪峰流量超過30 000 m3/s的水沙過程［16］。因此，場次洪水峰值的確定除需滿足寸灘站流量大于閾值Lk外，還需滿足寸灘站流量大于30 000 m3/s。

2.2 峰值獨立性的判別

由于樣本數(shù)據(jù)的波動，存在多個連續(xù)峰值同屬于一個獨立峰值過程的情況。因此，需要對場次洪水峰值進(jìn)行獨立性判別，以篩選得到各個獨立的峰值。峰值獨立性判別的標(biāo)準(zhǔn)包括2個峰值之間的時間間隔限制與2個峰值之間的最小值限制，如圖2及式（1）—式（2）所示。

Δtw>Δt（1）

Qmin<0.75min（Qp，Q′p）（2）

式中：Δtw為2個峰值流量的間隔時間，d;Δt為標(biāo)準(zhǔn)間隔時間，d;Qmin為2個峰值流量區(qū)間的最小流量，m3/s;Qp和Q′p分別為2個連續(xù)洪峰的峰值流量，m3/s。根據(jù)寸灘站多年來的實測洪水過程，選定Δt為7 d。

2.3 場次洪水起止時間的確定

場次洪水過程的起止時間由洪峰時間向外延伸得到，延伸標(biāo)準(zhǔn)包括起止時間到洪峰時間的間隔限制與起止時間值相對洪峰值的大小限制。對于洪峰過程的選取如圖2所示，場次洪水起止時間需滿足下式。

Qb≤0.5Qp（3）

Qe≤0.5Qp（4）

TQp-TQb≤1.5Δt（5）

TQe-TQp≤1.5Δt（6）

式中：Qb、Qe和Qp分別為場次洪水起始流量、結(jié)束流量和洪峰流量，m3/s;TQb、TQe和TQp分別為對應(yīng)的場次洪水起始時間、結(jié)束時間和峰值時間。

時間間隔一般要求不超過一定天數(shù)，即偏向于靠近峰值時間;而起止時間值要求不超過峰值某一比例，即偏向于遠(yuǎn)離峰值時間。當(dāng)峰值大小限制在時間間隔限制以外時，取時間間隔限制所在時間點;當(dāng)峰值大小限制在時間間隔限制以內(nèi)時，取峰值大小限制所在時間點;當(dāng)2個獨立峰值過程的起止時間范圍部分重疊時，以2個洪峰之間的最小值作為分割點。

2.4 洪峰沙峰異步傳播的影響

三峽水庫蓄水運行后，庫區(qū)水位抬高導(dǎo)致洪峰傳播時間減少，而沙峰輸移時間則由于水流速度減緩而增加，導(dǎo)致沙峰從寸灘站輸移至壩址的時間滯后于洪峰，且不同洪峰流量及壩前水位下的滯后時間也不相同［13，20］。張地繼等［26］采用數(shù)學(xué)模型對從寸灘到壩址的洪峰沙峰傳播時間進(jìn)行了模擬，本文在此基礎(chǔ)上計算各場次洪水的洪峰沙峰傳播時間，并對場次洪水的出庫水沙過程進(jìn)行延長，其中，為保證場次洪水排沙比計算的合理性，在統(tǒng)計出庫沙量時對起止時間均進(jìn)行相應(yīng)順延。2003—2021年三峽水庫共劃分得到55場次入庫洪水，考慮洪峰沙峰異步傳播影響前后的場次洪水入、出庫沙量年內(nèi)占比關(guān)系如圖3所示?？紤]洪峰沙峰異步傳播的影響后，三峽水庫場次洪水入、出庫沙量年內(nèi)占比的相關(guān)性有所提高，相關(guān)系數(shù)由0.671增加至0.731，且均大于7日高洪水期［15］入、出庫沙量年內(nèi)占比的相關(guān)系數(shù)0.157。對于三峽水庫場次洪水過程，雖然受到水庫調(diào)度等影響，出庫沙量往往會小于入庫沙量，但兩者之間仍存在一定的相關(guān)性，即當(dāng)入庫沙量較大時，出庫沙量也會較大。由此可見，考慮三峽庫區(qū)洪峰沙峰異步傳播的影響后，劃分得到的場次洪水過程在反映三峽水庫場次洪水的排沙規(guī)律上相對略優(yōu)。

3 三峽水庫場次洪水排沙規(guī)律

三峽水庫汛期場次洪水排沙情況統(tǒng)計如表1所示，各場次洪水排沙比變化波動較大，變化范圍為0.2%～84.0%?？紤]到三峽水庫在圍堰發(fā)電期控制庫水位為135～139 m，2006年10月進(jìn)入初期運行期，2008年汛后開始175 m試驗性蓄水，金沙江下游向家壩水電站于2012年蓄水運行，導(dǎo)致入庫沙量大幅減少［25］，對三峽水庫場次洪水排沙均產(chǎn)生一定影響，故分別統(tǒng)計了2003—2006年、2007—2012年、2013—2021年三峽水庫汛期場次洪水的排沙比均值，依次為41.7%、17.1%、18.7%。三峽水庫圍堰發(fā)電期汛期場次洪水平均排沙比達(dá)到了40%以上，而后隨著水庫分期蓄水進(jìn)程的逐步推進(jìn)，汛期入庫水沙條件及水庫運行條件等發(fā)生變化，場次洪水的水庫排沙能力有所降低，排沙比均值降至了20%以下，且2012年以后并未隨著入庫沙量的大幅減小而出現(xiàn)明顯變化。

進(jìn)一步統(tǒng)計了各年汛期場次洪水期間累計排沙的情況，如圖4所示，2003年、2004年、2005年汛期場次洪水的累計排沙比均達(dá)到了40%以上，2007—2021年汛期場次洪水累計排沙比均在30%以下。其中，各年汛期場次洪水的出庫沙量占全年出庫沙量的比例平均達(dá)到了67.0%，2018年最大達(dá)到了94.7%，表明三峽水庫排沙主要集中在汛期場次洪水過程。在金沙江下游梯級水庫蓄水?dāng)r沙導(dǎo)致三峽入庫沙量大幅減少的情況下，雖然2013年、2018年、2020年由于來水較大，出現(xiàn)了入庫沙量明顯增大的情況，汛期場次洪水來沙量遠(yuǎn)高于2013—2021年的均值7 611.1萬t，但是在2013年、2018年、2020年汛期開展沙峰排沙調(diào)度的作用下，汛期場次洪水的累計排沙比分別達(dá)到了28.4%、28.1%、24.5%，遠(yuǎn)高于2013—2021年各年汛期場次洪水累計排沙比的均值18.4%。由此表明，三峽水庫汛期沙峰排沙調(diào)度對于增大水庫排沙比具有一定效果，有效降低了水庫泥沙淤積的風(fēng)險，有利于水庫的長期使用。

4 三峽水庫場次洪水排沙影響因素

4.1 壩前水位與場次洪水排沙比的關(guān)系

壩前水位變化會改變庫區(qū)的水流流態(tài)，進(jìn)而改變水庫泥沙運動條件，影響水庫排沙。三峽水庫蓄水運行以來，依次經(jīng)歷了圍堰發(fā)電期、初期運行期、175 m蓄水運行期等不同階段，汛期場次洪水壩前平均水位逐步抬高，排沙比也相應(yīng)發(fā)生變化。圖5給出了三峽水庫汛期場次洪水排沙比與壩前平均水位的變化關(guān)系，隨著壩前平均水位的抬高，三峽水庫排沙比總體呈現(xiàn)出減小的趨勢。三峽水庫圍堰發(fā)電期，在場次洪水壩前平均水位相近的情況下，排沙比變化波動范圍仍較大，說明排沙比并非僅受壩前水位的影響，而是受到多因素的綜合作用。2003—2006年共有10場次洪水過程，其中，4場次洪水過程排沙比在30%以下，均值為22.6%，平均入庫流量僅為22 900 m3/s，其余6場次洪水過程排沙比在30%以上，均值為54.5%，平均入庫流量為30 300 m3/s，說明入庫流量等因素對排沙比也有較大影響。此外，三峽水庫在圍堰發(fā)電期基本不承擔(dān)防洪任務(wù)，由于區(qū)間支流入?yún)R，出庫流量普遍大于入庫流量，有利于水庫排沙，這也是圍堰發(fā)電期水庫排沙比較大的原因。三峽水庫進(jìn)入初期運行期和175 m蓄水運行期后，汛期承擔(dān)了防洪任務(wù)，2007—2021年場次洪水期間壩前水位均值達(dá)到了153.72 m，較圍堰發(fā)電期抬升了約18 m，2006年以后場次洪水排沙比均值僅為17.9%。

為了進(jìn)一步研究三峽水庫壩前水位對場次洪水排沙比的影響，圖6給出了不同入庫平均流量（寸灘站，Qin）下壩前平均水位與排沙比的關(guān)系。按照30 000 m3/s劃分入庫流量級后，Qin>30 000 m3/s的場次洪水壩前平均水位與排沙比的相關(guān)性大幅提高，相關(guān)系數(shù)達(dá)到了0.780;而Qin<30 000 m3/s的場次洪水則略有下降。綜合來看，三峽水庫汛期場次洪水排沙比總體隨著壩前水位的抬升而呈下降趨勢，且當(dāng)Qin>30 000 m3/s時，壩前水位與排沙比的相關(guān)性更高。

4.2 入、出庫輸沙率與場次洪水排沙比的關(guān)系

水庫排沙比是反映出庫沙量與入庫沙量比值的參數(shù)，沙量的大小取決于流量以及含沙量的大小，對單一場次洪水而言，應(yīng)綜合分析流量與含沙量對排沙比的影響，故分別統(tǒng)計了三峽水庫汛期場次洪水入庫平均輸沙率（寸灘站）、出庫平均輸沙率與排沙比的關(guān)系，如圖7所示。三峽水庫汛期場次洪水的入、出庫平均輸沙率與排沙比之間均呈現(xiàn)出一定的正相關(guān)關(guān)系，且出庫平均輸沙率與排沙比的相關(guān)性顯著大于入庫平均輸沙率，說明三峽水庫汛期場次洪水入、出庫輸沙率的增加均有利于水庫排沙，且出庫輸沙率的增加對于提高水庫排沙能力的作用更強(qiáng)。分析認(rèn)為主要是寸灘站距壩址較遠(yuǎn)，且區(qū)間存在支流入?yún)R等因素影響，導(dǎo)致入庫輸沙率對水庫排沙的影響較出庫輸沙率弱。

考慮到在不同壩前水位下入、出庫輸沙率對三峽水庫排沙比的影響會有所區(qū)別，且三峽水庫在圍堰發(fā)電期、初期運行期的汛期調(diào)洪能力有限，壩前水位不高，故本文以壩前平均水位150 m為界來劃分汛期場次洪水過程，探究不同運行水位條件下入、出庫輸沙率對水庫排沙比的影響，如圖8所示。當(dāng)壩前平均水位L<150 m時，三峽水庫汛期場次洪水入、出庫平均輸沙率與排沙比的相關(guān)性均有所提高，而當(dāng)L>150 m時，入、出庫輸沙率與排沙比的相關(guān)性均有所減小。圍堰發(fā)電期和初期運行期，三峽水庫汛期基本不進(jìn)行攔蓄或攔蓄量很小，壩前水位偏低，此時入、出庫輸沙率對水庫排沙比的影響較強(qiáng)，相關(guān)性有所提高;175 m蓄水運行后，三峽水庫會對一定量級的洪水進(jìn)行調(diào)蓄，壩前水位有所抬升，且入庫流量越大，壩前水位抬升越高，此時入、出庫輸沙率對水庫排沙比的影響受到削弱，相關(guān)性減小。

壩前水位與入、出庫輸沙率對三峽水庫場次洪水排沙比均有一定的趨勢性影響，且在劃分不同流量級或不同壩前水位后，二者與排沙比的相關(guān)性均出現(xiàn)了增大的情況，說明壩前水位與入、出庫輸沙率之間相互影響。由于水庫排沙比與壩前水位成反比，而與出庫輸沙率成正比，根據(jù)庫容關(guān)系曲線可得到水庫庫容（V），故本文引入了洪水滯留時間（V/Qout，其中Qout為出庫流量）來研究壩前水位與出庫流量對場次洪水排沙比的影響，如圖9所示。洪水滯留時間與場次洪水排沙比的相關(guān)性較單獨考慮壩前水位或入、出庫輸沙率時均有顯著提高，相關(guān)系數(shù)達(dá)到了0.810，說明二者的綜合作用對場次洪水排沙比具有重要影響。因此，今后在汛期發(fā)生洪水期間，可根據(jù)庫區(qū)排沙需求，結(jié)合三峽水庫壩前水位來調(diào)節(jié)出庫流量，以此來控制洪水滯留時間大小，為三峽水庫的排沙優(yōu)化調(diào)度提供重要參考。當(dāng)洪水滯留時間越小時，洪水沙峰過程越快排出水庫，場次洪水排沙比就越大;隨著洪水滯留時間增加，場次洪水排沙比逐漸降低，且降低的速率逐步趨緩。

基于前文劃分不同流量或水位的分析，本節(jié)討論不同入庫平均流量或不同壩前平均水位下洪水滯留時間與場次洪水排沙比的關(guān)系，如圖10所示。當(dāng)Qin>30 000 m3/s或L<150 m時，對應(yīng)的洪水滯留時間較小，相關(guān)性有所提高，表明當(dāng)洪水滯留時間較小時，場次洪水排沙比受洪水滯留時間的影響更大。此外，場次洪水L在150 m上下時，洪水滯留時間存在著相同的變化范圍，主要區(qū)間為7～8 d。但是，在洪水滯留時間相同時，L>150 m的場次洪水排沙比更大，這主要是由于其出庫平均流量較大，導(dǎo)致在L較高的情況下，場次洪水排沙比仍高于L<150 m的情況，同時也印證了三峽水庫場次洪水排沙比受到壩前水位與入、出庫輸沙率等多因素的綜合影響。

5 結(jié)? 論

本文選取了2003—2021年三峽水庫蓄水以來的入、出庫水文泥沙數(shù)據(jù)為研究對象，綜合考慮庫區(qū)洪峰沙峰異步傳播對水庫排沙的影響，分析了水庫汛期55場次洪水的排沙規(guī)律，并重點探討了壩前水位、入出庫輸沙率、洪水滯留時間等對場次洪水排沙的影響。主要結(jié)論如下：

（1） 三峽水庫場次洪水劃分在采用統(tǒng)一的定量標(biāo)準(zhǔn)基礎(chǔ)上，進(jìn)一步考慮洪峰沙峰異步傳播的影響，能夠更好地反映水庫汛期場次洪水的排沙規(guī)律。受到多重因素的綜合影響，2003—2021年場次洪水排沙比變化復(fù)雜，變化范圍為0.2%～84.0%。

（2） 三峽水庫排沙主要集中在汛期場次洪水過程，各年汛期場次洪水期間排沙總量平均占全年排沙的67.0%。三峽水庫圍堰發(fā)電期，場次洪水排沙比均值在40%以上，而進(jìn)入初期運行期及175 m蓄水運行期后，排沙比均值下降至20%以下。當(dāng)三峽水庫開展汛期沙峰排沙調(diào)度時，場次洪水排沙效果提升顯著。

（3） 三峽水庫場次洪水排沙比分別與壩前水位成反比、與入出庫輸沙率成正比關(guān)系，且當(dāng)入庫平均流量大于30 000 m3/s時，壩前水位與排沙比的相關(guān)性得到明顯提高;當(dāng)壩前平均水位小于150 m時，入出庫輸沙率與排沙比的相關(guān)性也有明顯提高。

（4） 采用洪水滯留時間反映壩前水位與出庫流量對排沙比的綜合影響，場次洪水排沙比與其成反比關(guān)系，相關(guān)性也更好，且當(dāng)洪水滯留時間較小時，兩者之間的相關(guān)性進(jìn)一步提升。

參考文獻(xiàn)：

［1］POFF N L，SCHMIDT J C.How dams can go with the flow［J］.Science，2016，353（6304）：1099-1100.

［2］鄧安軍，陳建國，胡海華，等.我國水庫淤損情勢分析［J］.水利學(xué)報，2022，53（3）：325-332.（DENG A J，CHEN J G，HU H H，et al.Analysis of reservoir siltation in China［J］.Journal of Hydraulic Engineering，2022，53（3）：325-332.（in Chinese））

［3］WHITE W R.Evacuation of sediments from reservoirs［M］.London：Thomas Telford，2001.

［4］袁晶，許全喜.金沙江流域水庫攔沙效應(yīng)［J］.水科學(xué)進(jìn)展，2018，29（4）：482-491.（YUAN J，XU Q X.Sediment trapping effect by reservoirs in the Jinsha River basin［J］.Advances in Water Science，2018，29（4）：482-491.（in Chinese））

［5］SCHLEISS A J，F(xiàn)RANCA M J，JUEZ C，et al.Reservoir sedimentation［J］.Journal of Hydraulic Research，2016，54（6）：595-614.

［6］夏軍強(qiáng)，張賢梓依，王增輝，等.黃河三門峽水庫水沙電耦合模型研究及應(yīng)用［J］.水科學(xué)進(jìn)展，2023，34（2）：238-249.（XIA J Q，ZHANG X Z Y，WANG Z H，et al.Study and application of a coupled modelling of flow-sediment transport and hydropower generation in the Sanmenxia Reservoir［J］.Advances in Water Science，2023，34（2）：238-249.（in Chinese））

［7］SISINGGIH D，WAHYUNI S，MARUF M F，et al.Evaluation of sustainable sediment management in Sengguruh Reservoir considering age and trapping efficiency［J］.International Journal of Geomate，2023，25（107）：115-122.

［8］JOTHIPRAKASH V，GARG V.Re-look to conventional techniques for trapping efficiency estimation of a reservoir［J］.International Journal of Sediment Research，2008，23（1）：76-84.

［9］朱玲玲，許全喜，董炳江，等.金沙江下游溪洛渡水庫排沙效果及影響因素［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））

［10］張倩倩，童思陳，龍孝平，等.三峽水庫排沙比與水沙因子響應(yīng)關(guān)系［J］.水電能源科學(xué)，2023，41（8）：101-104，112.（ZHANG Q Q，TONG S C，LONG X P，et al.Response relationship between sediment discharge ratio with water flow and sediment factors in Three Gorges Reservoir［J］.Water Resources and Power，2023，41（8）：101-104，112.（in Chinese））

［11］胡春宏.三峽水庫175 m試驗性蓄水十年泥沙沖淤變化分析［J］.水利水電技術(shù)，2019，50（8）：18-26.（HU C H.Analysis on sediment scouring and silting variation of Three Gorges Reservoir since 175 m trial impoundment for past ten years［J］.Water Resources and Hydropower Engineering，2019，50（8）：18-26.（in Chinese））

［12］鄭守仁.三峽工程175米試驗性蓄水運行期的泥沙觀測分析［J］.長江技術(shù)經(jīng)濟(jì)，2019，3（3）：1-7.（ZHENG S R.Analysis of sediment observation data in the trial impounding period of the Three Gorges Project to pool level of El.175-meter［J］.Technology and Economy of Changjiang，2019，3（3）：1-7.（in Chinese））

［13］張為，李昕，任金秋，等.梯級水庫蓄水對三峽水庫洪峰沙峰異步特性的影響［J］.水科學(xué)進(jìn)展，2020，31（4）：481-490.（ZHANG W，LI X，REN J Q，et al.Influence of the impoundment of cascade reservoirs on the asynchronies of flood peak and sediment peak in the Three Gorges Reservoir［J］.Advances in Water Science，2020，31（4）：481-490.（in Chinese））

［14］單敏爾，李志晶，周銀軍，等.三峽水庫入庫水沙變化規(guī)律及驅(qū)動因素分析［J］.泥沙研究，2022，47（2）：29-35.（SHAN M E，LI Z J，ZHOU Y J，et al.Variations and driving factors of water flow and sediment transport in the Three Gorges Reservoir［J］.Journal of Sediment Research，2022，47（2）：29-35.（in Chinese））

［15］高宇，任實，王海，等.金沙江下游梯級水庫運行后三峽水庫高洪水期入庫水沙特性［J］.湖泊科學(xué)，2023，35（2）：662-672.（GAO Y，REN S，WANG H，et al.Characteristics of incoming runoff and sediment load in the Three Gorges Reservoir during high flood period after the operation of the cascade reservoirs in the Lower Jinsha River［J］.Journal of Lake Sciences，2023，35（2）：662-672.（in Chinese））

［16］李思璇，楊成剛，董炳江，等.長江上游高洪水期泥沙輸移特性［J］.長江科學(xué)院院報，2021，38（12）：6-11.（LI S X，YANG C G，DONG B J，et al.Sediment transport characteristics during high floods in the upper reaches of the Yangtze River［J］.Journal of Yangtze River Scientific Research Institute，2021，38（12）：6-11.（in Chinese））

［17］王玉璇，金中武，陳鵬.三峽水庫汛期排沙效率研究［J］.人民長江，2022，53（5）：14-19.（WANG Y X，JIN Z W，CHEN P.Study on sediment transport efficiency of Three Gorges Reservoir during flood season［J］.Yangtze River，2022，53（5）：14-19.（in Chinese））

［18］董占地，胡海華，吉祖穩(wěn)，等.三峽水庫排沙比對來水來沙的響應(yīng)［J］.泥沙研究，2017，42（6）：16-21.（DONG Z D，HU H H，JI Z W，et al.Response of sediment excluding ratio to incoming flow and sediment in the Three Gorges Reservoir［J］.Journal of Sediment Research，2017，42（6）：16-21.（in Chinese））

［19］黃仁勇，談廣鳴，范北林.三峽水庫蓄水運用后汛期洪水排沙比初步研究［J］.水力發(fā)電學(xué)報，2013，32（5）：129-133，152.（HUANG R Y，TAN G M，F(xiàn)AN B L.Sediment delivery ratio of Three Gorges Reservoir in flood events after its impoundment［J］.Journal of Hydroelectric Engineering，2013，32（5）：129-133，152.（in Chinese））

［20］張幫穩(wěn)，吳保生，章若茵.三峽庫區(qū)汛期洪峰和沙峰異步運動特性的三維數(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））

［21］ZHANG Q，ZHANG L P，SHE D X，et al.Automatic procedure for selecting flood events and identifying flood characteristics from daily streamflow data［J］.Environmental Modelling & Software，2021，145：105180.

［22］DUROCHER M，MOSTOFI ZADEH S，BURN D H，et al.Comparison of automatic procedures for selecting flood peaks over threshold based on goodness-of-fit tests［J］.Hydrological Processes，2018，32（18）：2874-2887.

［23］VITTAL H，SINGH J，KUMAR P，et al.A framework for multivariate data-based at-site flood frequency analysis：essentiality of the conjugal application of parametric and nonparametric approaches［J］.Journal of Hydrology，2015，525：658-675.

［24］王玉璇，金中武，陳鵬，等.三峽水庫汛期水沙異步特性及沙峰輸移時間研究［J］.長江科學(xué)院院報，2022，39（9）：24-29，37.（WANG Y X，JIN Z W，CHEN P，et al.Water-sediment asynchrony and transport time of sand peak in flood season of Three Gorges Reservoir［J］.Journal of Yangtze River Scientific Research Institute，2022，39（9）：24-29，37.（in Chinese））

［25］郭率，周曼，胡挺，等.近年來三峽水庫水沙通量變化過程及原因分析［J］.中國農(nóng)村水利水電，2022（8）：35-40，49.（GUO S，ZHOU M，HU T，et al.An analysis of the variation process and causes of water and sediment to the Three Gorges Reservoir in recent years［J］.China Rural Water and Hydropower，2022（8）：35-40，49.（in Chinese））

［26］張地繼，董炳江，楊霞，等.三峽水庫庫區(qū)沙峰輸移特性研究［J］.人民長江，2018，49（2）：23-28，68.（ZHANG D J，DONG B J，YANG X，et al.Research on sediment peak movement characteristics in Three Gorges Reservoir［J］.Yangtze River，2018，49（2）：23-28，68.（in Chinese））

Sediment delivery regularity and its influencing factors during flood events in the Three Gorges Reservoir

Abstract：Reservoir sedimentation leads to a loss of storage capacity，thereby affecting the ability of a reservoir to fulfil its designated purpose.Accordingly，it is important to investigate the regularity of sediment delivery to constrain reservoir sedimentation.On the basis of data recorded at Zhutuo，Beibei，Wulong，Cuntan，and other hydrological stations between 2003 and 2021，individual flood events in the Three Gorges Reservoir （TGR） were identified.The regularity of sediment delivery during each flood event was then analyzed.The impacts of factors such as the water level in front of the dam，the rate of sediment transport into and out of the reservoir，and flood detention time on sediment delivery during flood events were investigated.The results indicate that it is more reasonable to classify the flood events in the TGR by considering the asynchronous propagation of flood and sediment peaks on the basis of a unified quantitative standard.Sediment delivery in the TGR occurred mainly during flood events，accounting for 67.0% of the annual sediment delivery volume.The sediment delivery ratio （SDR） during the flood events was inversely proportional to the water level in front of the dam and proportional to the sediment transport rate into and out of the reservoir.However，the correlations were weaker than that between the SDR and flood detention time，whereby the SDR was inversely proportional to the flood detention time.The correlations between the SDR and other factors were strongest during flood events when the average inflow was greater than 30 000 m3/s or the average water level in front of the dam was less than 150 m.This study provides an important reference for optimizing the operation of the TGR and reducing sediment deposition.

Key words：sediment delivery ratio;flood events;asynchronies of flood and sediment peaks;flood detention time;Three Gorges Reservoir