長江上游三峽水庫入庫泥沙沙峰傳播時(shí)間研究

姜利玲 董炳江 湯成友 董溢

摘要：

為進(jìn)一步掌握長江上游三峽水庫入庫泥沙沙峰傳播規(guī)律、提高入庫泥沙沙峰預(yù)報(bào)精度，制定更加科學(xué)合理的三峽水庫沙峰調(diào)度方案，基于三峽水庫蓄水以來長江上游主要控制性水文站水沙實(shí)測資料，分析了長江上游水沙組成，并對(duì)長江上游各主要控制性水文站沙峰傳播時(shí)間進(jìn)行了研究。結(jié)果表明：2003年以來，長江上游來沙明顯減少，向家壩站、橫江站、高場站及富順站沙峰到朱沱站的平均傳播時(shí)間分別為38.7，32.0，30.9 h和20.8 h，武勝站、羅渡溪站和小河壩站沙峰到北碚站平均傳播時(shí)間分別為16.9，19.0，20.0 h，朱沱站、北碚站到寸灘站分別為20.8，7.7 h;當(dāng)向家壩至朱沱河段流量小于10 000 m3/s時(shí)，沙峰傳播時(shí)間為48.0 h，當(dāng)流量為10 000～20 000 m3/s時(shí)，河段沙峰傳播時(shí)間為32.7～46.3 h，當(dāng)流量為20 000～30 000 m3/s時(shí)，河段沙峰傳播時(shí)間為27.0～32.7 h，當(dāng)流量大于30 000 m3/s時(shí)，河段沙峰傳播時(shí)間少于27.0 h。

關(guān) 鍵 詞：

沙峰輸移; 傳播時(shí)間; 頻率分布; 數(shù)值模擬;? 三峽水庫;? 長江上游

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

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

DOI：10.16232/j.cnki.1001-4179.2021.07.002

0 引 言

三峽水庫175 m試驗(yàn)性蓄水以來，通過利用庫區(qū)洪峰沙峰異步傳播特性和入庫泥沙集中的特點(diǎn)，結(jié)合水情與泥沙預(yù)報(bào)成果，2012，2013，2018年成功實(shí)施了汛期沙峰排沙調(diào)度，即在洪峰到來時(shí)進(jìn)行攔洪削峰，減輕下游防洪壓力;當(dāng)沙峰臨近大壩時(shí)水庫加大泄量、泄洪排沙，取得了較好的排沙效果。在水庫沙峰排沙調(diào)度過程中，有效預(yù)測沙峰到達(dá)水庫和壩前的時(shí)間，對(duì)于更好地掌握沙峰調(diào)度時(shí)機(jī)、制訂更科學(xué)合理的沙峰排沙調(diào)度方案具有重要的技術(shù)支撐作用[1-2]。

長期以來水文預(yù)報(bào)都是人們關(guān)注的熱點(diǎn)問題，許多學(xué)者在泥沙輸移問題方面做了大量研究探索，取得了較多成果。張根廣等[3]提出應(yīng)多考慮采用隨機(jī)現(xiàn)象的統(tǒng)計(jì)學(xué)方法來描述泥沙起動(dòng)狀態(tài)。李慶金[4]等將回歸分析用于黃河中游高村水文站的水沙預(yù)報(bào)中。凌燕[5]、李楠[6]分別通過建立不平衡輸沙漠型和不確定性預(yù)報(bào)系統(tǒng)模型對(duì)黃河流域潼關(guān)站含沙量預(yù)報(bào)工作進(jìn)行探索。董炳江等[7-8]對(duì)沙峰與洪峰在三峽庫區(qū)傳播特性開展了分析研究;黃仁勇[9-11]建立了三峽水庫下游沙峰傳播時(shí)間公式，并發(fā)明了一種汛期出庫沙峰傳播時(shí)間的快速預(yù)測方法;王世平等[12]對(duì)三峽水庫入庫站含沙量預(yù)報(bào)進(jìn)行了初步探索。目前關(guān)于長江上游入庫沙峰傳播規(guī)律方面的研究仍較少，隨著長江上游干支流水庫群的數(shù)量不斷增加，受來沙環(huán)境的變化與水庫調(diào)度方式的不同因素影響，三峽水庫入庫沙峰傳播表現(xiàn)出更為復(fù)雜的規(guī)律。

本文基于三峽水庫建成以來2003～2018年長江上游主要控制站含沙量實(shí)測數(shù)據(jù)，采用資料分析和數(shù)值模擬等方法，系統(tǒng)地研究了長江上游三峽水庫入庫泥沙沙峰的傳播時(shí)間。

1 長江上游水沙來源

長江上游徑流的主要來源為金沙江、橫江、岷江、沱江、嘉陵江和烏江等河流，而懸移質(zhì)泥沙主要來源于金沙江和嘉陵江。1990年以來，長江上游徑流量變化不大，但受水利工程攔沙、降雨條件變化、水土保持減沙措施和河道采砂等因素的綜合影響，長江輸沙量明顯減少。21世紀(jì)之后，三峽水庫上游來沙減小的趨勢仍然持續(xù)[13]。

近年來觀測資料表明：長江上游來沙量大幅減少，三峽水庫入庫沙量維持較低水平。在三峽水庫的論證與初步設(shè)計(jì)階段，提出的“由于長江上游地區(qū)實(shí)施的水土保持工作以及水庫的運(yùn)行，三峽水庫入庫沙量將呈減少趨勢”是符合實(shí)際的。根據(jù)實(shí)測資料分析，2003～2018年，三峽水庫的主要控制站——朱沱、北碚、武隆站多年平均徑流量、輸沙量之和分別是：3 645億m3和1.54億t，較1990年以前分別減小6%和68%，較1991～2002年分別減少2%和56%。長江上游主要水文站徑流量和輸沙量多年變化情況比較結(jié)果如圖1～2所示。

受上游水電站攔沙等影響，三峽水庫入庫泥沙來源發(fā)生了明顯改變[14]。三峽水庫建成以來，2003～2018年金沙江屏山（向家壩）站年均沙量占寸灘站沙量的比例為62.6%（1956～1990年占比為53.4%、1991～2002年占比為83.4%）;岷江高場站、橫江橫江站以及沱江富順站所對(duì)應(yīng)的來沙量比重大小分別為：17.0%、4.3%和3.8%（1956～1990年分別占比為11.4%、3.0%和2.5%，1991～2002年分別占比為10.2%、4.1%和1.1%）;嘉陵江北碚站年均沙量占寸灘站沙量的比例為19.8%（1956～1990年占比為29.1%、1991～2002年占比為11.0%）。特別是2013年以來，受溪洛渡、向家壩等水電站蓄水?dāng)r沙影響，2013～2018年金沙江來沙量占寸灘沙量的比重減少為2.4%，而橫江站、高場站、富順站和北碚站來沙量占寸灘沙量比重分別增大至10.6%、22.5%、15.7%和38.5%。

2 長江上游沙峰傳播時(shí)間分析

長江干流朱沱、寸灘站，分別距三峽大壩759，606 km;嘉陵江、烏江分別在寸灘站上游約6 km、清溪場站上游約11 km匯入長江，其控制站北碚、武隆站分別距匯合口56，70 km（見圖3）。近年來，烏江來沙量較小，本文暫不對(duì)武隆站沙峰傳播時(shí)間多加分析。

本文主要采用長江上游干支流控制站向家壩站（屏山站）、橫江站、高場站、富順站、朱沱站、羅渡溪站、武勝站、小河壩站、北碚站和寸灘站洪水含沙量實(shí)測資料對(duì)沙峰傳播時(shí)間進(jìn)行分析。先利用SPSS軟件對(duì)長江上游向家壩站、橫江站、高場站及富順站到朱沱站實(shí)測沙峰傳播時(shí)間，嘉陵江流域武勝站、羅渡溪站及小河壩站到北碚站實(shí)測沙峰傳播時(shí)間，朱沱站、北碚站到三峽水庫入庫控制站寸灘站實(shí)測沙峰傳播時(shí)間進(jìn)行頻次分析;再統(tǒng)計(jì)各河段的沙峰傳播時(shí)間在不同流量級(jí)洪峰下的變化。

分析選取的沙峰和洪峰時(shí)間范圍為2003年1月至2018年12月。另外，溪洛渡和向家壩這2座特大型水庫已分別于2013年和2012年開始蓄水運(yùn)用，向家壩站作為向家壩水電站出庫控制站，受水電站蓄水?dāng)r沙影響，該站來沙量明顯減少，基于該原因，本文向家壩站分析數(shù)據(jù)為2003～2013年洪水含沙量實(shí)測數(shù)據(jù)，其余各控制站分析數(shù)據(jù)均為2003～2018年洪水含沙量實(shí)測數(shù)據(jù)。

本文以朱沱、北碚站到寸灘站沙峰傳播時(shí)間分析為例，先在上游朱沱、北碚站實(shí)測水沙過程中挑選較為明顯的沙峰過程，選取的沙峰應(yīng)同時(shí)能在下游寸灘水文站觀察到明顯的洪峰沙峰過程，即寸灘站出現(xiàn)1.0 kg/m3沙峰過程。最終從實(shí)測資料中選取了106組符合條件的沙峰輸移資料，其中從朱沱站至寸灘站沙峰輸移過程79組，北碚站至寸灘站沙峰輸移過程27組。朱沱站至寸灘站實(shí)測沙峰傳播時(shí)間變化范圍為8.0～37.0 h，朱沱站沙峰變化范圍為1.24～8.38 kg/m3，朱沱站洪峰變化范圍為5 240～56 500 m3/s，朱沱站傳播至寸灘站沙峰變化范圍為1.03～6.41 kg/m3;北碚站至寸灘站實(shí)測沙峰傳播時(shí)間變化范圍為1.0～24.0 h，北碚站沙峰變化范圍為1.48～14.60 kg/m3，北碚站洪峰變化范圍為7 770～35 700 m3/s，傳播至寸灘站沙峰變化范圍為1.06～5.32 kg/m3。該系列沙峰輸移資料的沙峰含沙量、沙峰傳播時(shí)間、洪峰流量變化范圍較廣，能夠體現(xiàn)不同的沙峰狀況，代表性較強(qiáng)。

2.1 向家壩-朱沱區(qū)間沙峰傳播時(shí)間分析

向家壩站位于金沙江干流上，位于電站壩址下游2.0 km處，距離朱沱站268.8 km;橫江橫江站、岷江高場站及沱江富順站距離朱沱站分別為279，262，217 km。實(shí)測資料分析表明：向家壩站、橫江站、高場站、富順站至朱沱站沙峰傳播時(shí)間頻率直方圖如圖4～7所示，其平均時(shí)間分別為38.7，32.0，30.9，20.8 h。

2003～2013年，向家壩站一共出現(xiàn)了94次沙峰過程，其中僅金沙江單江漲水來沙的沙峰過程有54次，還有40次沙峰過程為：金沙江漲水來沙的同時(shí)，橫江、岷江及沱江其中1～3條江也出現(xiàn)漲水來沙，即為多江同時(shí)漲水時(shí)出現(xiàn)的沙峰過程。2003～2018年橫江站、高場站及富順站分別出現(xiàn)了50，31，12次沙峰過程，其中橫江站、高場站所在河流單江漲水來沙的沙峰過程分別為8，5次，沱江富順站沒有出現(xiàn)單江漲水來沙的沙峰過程。

單江漲水來沙情況下，向家壩站、橫江站、高場站至朱沱站沙峰傳播平均時(shí)間分別為38.6，28.9，27.4 h;多江同時(shí)漲水來沙的情況下，向家壩站、橫江站、高場站、富順站至朱沱站沙峰傳播平均時(shí)間分別為38.8，32.6，31.6 h和20.8 h。

朱沱站不同流量級(jí)洪峰下向家壩站、橫江站、高場站、富順站至朱沱站沙峰傳播平均時(shí)間分析統(tǒng)計(jì)如表1所列。由表1可見：干流向家壩至朱沱沙峰傳播平均時(shí)間隨著洪峰流量增大而有所縮短，主要是因?yàn)閭鞑ニ俣入S流量增大而加快，故沙峰傳播平均時(shí)間會(huì)有所縮短;支流橫江站、高場站、富順站至朱沱站沙峰傳播平均時(shí)間受洪峰流量影響不大，朱沱站不同流量級(jí)洪峰下橫江站、高場站至朱沱站沙峰傳播平均時(shí)間基本在27.0～37.0 h，富順站至朱沱站沙峰傳播平均時(shí)間基本在16.0～25.0 h。

2.2 嘉陵江下游武勝-北碚區(qū)間沙峰傳播時(shí)間分析

嘉陵江干流武勝站、渠江羅渡溪站和涪江小河壩站距離北碚站分別約為125.1，139.4 km和142 km。根據(jù)實(shí)測數(shù)據(jù)分析，各站至北碚站沙峰傳播時(shí)間頻率直方圖如圖8～10所示，其平均時(shí)間分別為16.9，19.0，20.0 h。

2003～2018年，武勝站、羅渡溪站及小河壩站分別出現(xiàn)了20，26，20次沙峰過程，各站所在河流單江漲水來沙的沙峰過程分別為7，16，16次。單江漲水來沙情況下，武勝站、羅渡溪站及小河壩站至北碚站沙峰傳播平均時(shí)間分別為14.7，18.8 h和17.6 h;多江同時(shí)漲水來沙的情況下，各站至北碚站沙峰傳播平均時(shí)間分別為18.1，19.3 h和29.6 h。由此可見，單江漲水來沙時(shí)，沙峰傳播平均時(shí)間比多江漲水來沙沙峰傳播平均時(shí)間要短，尤其是小河壩站，單江漲水來沙比多江漲水來沙情況下沙峰傳播平均時(shí)間要短12.0 h。

北碚站不同流量級(jí)洪峰下武勝站、羅渡溪站及小河壩站至北碚站沙峰傳播平均時(shí)間統(tǒng)計(jì)如表2所列。北碚站不同流量級(jí)洪峰下單江漲水來沙比多江漲水來沙情況下沙峰傳播平均時(shí)間要短0.1～16.2 h?？傮w上，同一種來沙情況下，嘉陵江流域武勝站、羅渡溪站及小河壩站至北碚站沙峰傳播平均時(shí)間受北碚站洪峰流量影響相對(duì)較小，單江漲水來沙情況下各控制站至北碚站沙峰傳播平均時(shí)間基本在12.0～24.0 h;多江同時(shí)漲水來沙情況下各控制站至北碚站沙峰傳播平均時(shí)間基本在12.0～36.0 h。

2.3 朱沱-寸灘區(qū)間沙峰傳播時(shí)間分析

長江上游干流朱沱站和嘉陵江北碚站距離寸灘站分別約為150 km和62 km。實(shí)測資料分析表明，兩站至寸灘站沙峰傳播平均時(shí)間分別為20.8，7.7 h。兩站至寸灘站的沙峰傳播時(shí)間實(shí)測值分析的頻率直方圖如圖11～12所示。

2003～2018年，朱沱站和北碚站分別出現(xiàn)了79次和27次沙峰過程，長江單獨(dú)漲水來沙的沙峰過程為73次，嘉陵江單獨(dú)漲水來沙的沙峰過程為21次，長江、嘉陵江兩江同時(shí)漲水來沙的沙峰過程為6次。單江漲水來沙情況下，朱沱站和北碚站至寸灘站沙峰傳播平均時(shí)間分別為20.6 h和7.9 h;兩江同時(shí)漲水來沙的情況下，至寸灘站沙峰傳播平均時(shí)間分別為23.2 h和7.0 h，2種情況下沙峰傳播平均時(shí)間相差最大為2.6 h。

不同流量級(jí)洪峰下朱沱站及北碚站至寸灘站沙峰傳播平均時(shí)間分析統(tǒng)計(jì)如表3所列，長江朱沱站及嘉陵江北碚站至寸灘站沙峰傳播平均時(shí)間受洪峰流量影響相對(duì)較小，朱沱站至寸灘站沙峰傳播平均時(shí)間基本為12.0～24.0 h，北碚站至寸灘站沙峰傳播平均時(shí)間基本為12.0 h以內(nèi)。

3 數(shù)值模擬

在已有研究成果和充分認(rèn)識(shí)長江上游干流泥沙沖淤特性的基礎(chǔ)上，基于水力學(xué)及河流動(dòng)力學(xué)的基本原理，按照長江上游干支流水庫群的具體建設(shè)情況、來水來沙環(huán)境改變狀況以及水庫調(diào)度方式優(yōu)化等，建立長江上游水庫群聯(lián)合調(diào)度和新水沙條件下的向家壩下游到朱沱站河道一維水沙數(shù)學(xué)模型分析研究長江上游干流沙峰傳播時(shí)間。

一維水沙數(shù)學(xué)模型理論基礎(chǔ)及數(shù)值離散格式均比較成熟[16]，基本控制方程包括水流連續(xù)方程、水流運(yùn)動(dòng)方程、河床變形方程、水流挾沙力公式及泥沙連續(xù)方程，通常采用Preissmann四點(diǎn)偏心隱格式進(jìn)行離散，這種格式的數(shù)值穩(wěn)定性及質(zhì)量守恒性較好，國內(nèi)外應(yīng)用較為普遍。

3.1 向家壩下游到朱沱站河道水力特性

向家壩至朱沱河道長約267 km，該河段未受三峽水庫回水影響，基本屬于天然河道特性，同時(shí)從2003～2018年朱沱站的水位流量關(guān)系（見圖13）可知，朱沱站水位流量關(guān)系比較單一。

從本質(zhì)上分析，洪水是按照波動(dòng)條件傳遞的波動(dòng)過程，而河流泥沙輸移則是按質(zhì)點(diǎn)速度傳遞的對(duì)流過程。洪峰是以波的形式傳播的，而沙峰輸移則與斷面平均流速有關(guān)，二者的傳播速度是不同的。從向家壩至朱沱河段水流速度與重力波速度統(tǒng)計(jì)來看，河道基本屬于水體自由流動(dòng)區(qū)。當(dāng)河段流量為10 000 m3/s時(shí)，重力波速基本為1.5～2.3倍平均流速，當(dāng)河段流量為40 000 m3/s時(shí)，重力波速與平均流速的差距縮小，基本為1.0～1.6倍平均流速（見表4）。

3.2 不同流量下向家壩-朱沱站沙峰傳播時(shí)間

不同流量與沙峰傳播時(shí)間關(guān)系如圖14所示，可以看出，兩者乘冪關(guān)系吻合度較高，河段內(nèi)沙峰傳播時(shí)間基本可用以下公式表示：

式中：T為沙峰傳播時(shí)間，h;Q為河段流量，m3/s。

從向家壩至朱沱河道不同河段的沙峰傳播時(shí)間來看，各不同河段的沙峰傳播時(shí)間基本可按照長度比例進(jìn)行分配。

統(tǒng)計(jì)不同流量級(jí)下向家壩至朱沱河段沙峰傳播時(shí)間列于表5，向家壩至朱沱河段當(dāng)流量小于10 000 m3/s，沙峰傳播時(shí)間不少于46.3 h，當(dāng)流量位于10 000～20 000 m3/s區(qū)間時(shí)，河段沙峰傳播時(shí)間為32.7～46.3 h，當(dāng)流量位于20 000～30 000 m3/s區(qū)間時(shí)，河段沙峰傳播時(shí)間為27.0～32.7 h，當(dāng)流量大于30 000 m3/s時(shí)，河段沙峰傳播時(shí)間少于27.0 h。朱沱站不同流量級(jí)洪峰下向家壩至朱沱河段沙峰傳播時(shí)間計(jì)算成果與實(shí)測值統(tǒng)計(jì)分析成果基本一致，隨著洪峰流量增

4 結(jié) 論

（1） 近年來，三峽水庫入庫年徑流量變化不大，但入庫泥沙地區(qū)組成發(fā)生明顯變化。溪洛渡、向家壩等水電站投運(yùn)后，金沙江來沙量占寸灘站沙量的比重大幅減小，其他支流則有所增大。

（2） 2003～2018年，長江上游干流向家壩站、橫江橫江站、岷江高場站、沱江富順站至朱沱站的沙峰傳播平均時(shí)間分別為38.7，32.0，30.9，20.8h。干支流不同來水對(duì)沙峰傳播時(shí)間有一定影響，如在干支流單江漲水來沙時(shí)，向家壩站、橫江站、高場站至朱沱站沙峰傳播平均時(shí)間分別為38.6，28.9，27.4 h;干支流同時(shí)漲水來沙時(shí)，沙峰傳播平均時(shí)間則分別為38.8，32.6，31.6 h。

根據(jù)數(shù)值模擬計(jì)算成果，不同流量級(jí)下向家壩至朱沱河段沙峰傳播時(shí)間：當(dāng)流量小于10 000 m3/s時(shí)，沙峰傳播時(shí)間不少于46.3 h;流量10 000～20 000 m3/s時(shí)，傳播時(shí)間為32.7～46.3 h;當(dāng)流量為20 000～30 000 m3/s對(duì)應(yīng)傳播時(shí)間為27.0～32.7 h;當(dāng)流量大于30 000 m3/s時(shí)，沙峰傳播時(shí)間少于27.0 h。根據(jù)實(shí)測資料分析成果，當(dāng)流量小于10 000 m3/s時(shí)，沙峰傳播平均時(shí)間為57.0 h;當(dāng)流量為10 000～20 000 m3/s時(shí)，傳播平均時(shí)間為41.1 h;當(dāng)流量為20 000～30 000 m3/s時(shí)，對(duì)應(yīng)傳播平均時(shí)間為34.9 h;當(dāng)流量大于30 000 m3/s時(shí)，沙峰傳播平均時(shí)間為23.7 h。向家壩至朱沱沙峰傳播時(shí)間隨流量的增大而有所減小，不同流量級(jí)洪峰下向家壩至朱沱河段沙峰傳播時(shí)間計(jì)算成果與實(shí)測值統(tǒng)計(jì)分析成果基本一致。

（3） 嘉陵江干流武勝站、渠江羅渡溪站、涪江小河壩站至北碚站沙峰傳播平均時(shí)間分別為16.9，19.0，20.0 h。其中，嘉陵江干流、支流單江漲水來沙情況下，各控制站至北碚站沙峰傳播平均時(shí)間為12.0～24.0 h;同時(shí)漲水來沙情況下，沙峰傳播平均時(shí)間則為12.0～36.0 h。

（4） 朱沱站、北碚站至寸灘站沙峰傳播平均時(shí)間分別為20.8，7.7 h。當(dāng)長江、嘉陵江單江漲水來沙時(shí)，朱沱站和北碚站至寸灘站沙峰傳播平均時(shí)間分別為20.6 h和7.9 h;兩江同時(shí)漲水來沙的情況下，至寸灘站沙峰傳播平均時(shí)間分別為23.2 h和7.0 h，2種情況下沙峰傳播平均時(shí)間相差最大為2.6 h。

由于水沙特性以及河道邊界條件等的綜合影響，沙峰沿程表現(xiàn)出的傳播規(guī)律相對(duì)較為復(fù)雜，并且不是一成不變的。另外，沙峰的輸移不僅與流量大小、河段長度和泥沙粒徑大小有關(guān)，降雨強(qiáng)度以及電站調(diào)度方式等因素也對(duì)沙峰的輸移有重要影響，綜合考慮這些因素的影響，有待進(jìn)一步深入研究。本文分析可視為一種積極的探索和嘗試，研究成果還有待于今后進(jìn)一步的深化和完善。

參考文獻(xiàn)：

[1] 陳桂亞，董炳江，姜利玲，等.2018長江2號(hào)洪水期間三峽水庫沙峰排沙調(diào)度[J].人民長江，2018，49（19）：6-10.

[2] 張地繼，董炳江，楊霞，等.三峽水庫庫區(qū)沙峰輸移特性研究[J].人民長江，2018，49（2）：23-28.

[3] 張根廣，周雙，王愉樂，等.泥沙輸移問題研究進(jìn)展與展望[J].水利與建筑工程學(xué)報(bào)，2019，17（4）：8-15.

[4] 李慶金，孟憲靜，劉麗，等.回歸分析在高村水文站水沙預(yù)報(bào)中的應(yīng)用[J].人民黃河，2013，35（5）：14-17.

[5] 凌燕.基于泥沙運(yùn)動(dòng)理論的潼關(guān)站含沙量過程預(yù)報(bào)方法[D].西安：西安理工大學(xué)，2009.

[6] 李楠.含沙量過程不確定性預(yù)報(bào)系統(tǒng)模型研究：以潼關(guān)站為例[D].西安：西安理工大學(xué)，2009.

[7] 董炳江，喬偉，許全喜.三峽水庫汛期沙峰排沙調(diào)度研究與初步實(shí)踐[J].人民長江，2014，45（3）：7-11.

[8] 董炳江，陳顯維，許全喜.三峽水庫沙峰調(diào)度試驗(yàn)研究與思考[J].人民長江，2014，45（19）：1-5.

[9] 黃仁勇.長江上游梯級(jí)水庫泥沙輸移與泥沙調(diào)度研究[D].武漢：武漢大學(xué)，2016.

[10] 黃仁勇，王敏，藺秋生，等.一種三峽水庫汛期沙峰傳播時(shí)間的快速預(yù)測方法：中國，201610804641.1[P].2017-02-15.

[11] 黃仁勇，崔占峰，葛華，等.一種大型河道型水庫汛期出庫沙峰含沙量的快速預(yù)測方法：中國，201610804942.4[P]2018-01-09.

[12] 王世平，王渺林，許全喜，等.三峽水庫入庫站含沙量預(yù)報(bào)方法初探與試預(yù)報(bào)[J].水利水電快報(bào)，2015，36（5）：11-14.

[13] 黃仁勇，談廣鳴，范北林，等.長江上游梯級(jí)水庫聯(lián)合調(diào)度泥沙數(shù)學(xué)模型研究[J].水力發(fā)電學(xué)報(bào)，2012，31（6）：143-148.

[14] 金興平，許全喜.長江上游水庫群聯(lián)合調(diào)度中的泥沙問題[J].人民長江，2018，49（3）：1-8.

[15] 長江水利委員會(huì)水文局.2018年度三峽水庫進(jìn)出庫水沙特性、水庫淤積及壩下游河道沖刷分析[R].武漢：長江水利委員會(huì)水文局，2019.

[16] 謝鑒衡.河流模擬[M].北京：水利電力出版社，1990.

（編輯：江 文）

Study on sand peak propagation time of input sediment of Three Gorges

Reservoir in Upper Changjiang River

JIANG Liling1，DONG Bingjiang2，TANG Chengyou1，DONG Yi1

（1.Upper Changjiang River Bureau of Hydrological and Water Resources Survey，Hydrology Bureau of Changjiang Water Resources Commission，Chongqing 400014，China; 2.Bureau of Hydrology，Changjiang Water Resources Commission，Wuhan 430010，China）

Abstract：

To further grasp the law of sand peak propagation and improve the prediction accuracy of sand peak in the upper reaches of the Changjiang River，and formulate a more scientific and reasonable dispatch scheme of sand peak in the Three Gorges Reservoir，we analyzed the composition of water and sediment in the upper reaches of the Changjiang River based on the measured data of water and sediment of the main hydrological stations in the upper reaches of the Changjiang River since the impoundment of the Three Gorges Reservoir.In addition，the propagation time of sand peak at the main hydrological stations in the upper reaches of the Changjiang River was preliminarily studied.The results show that since 2003 the sediment in the upper reaches of the Changjiang River has significantly reduced，and the average propagation time of sand peak from Xiangjiaba station，Hengjiang station，Gaochang station and Fushun station to Zhutuo station is 38.7 h，32.0 h，30.9 h and 20.8 h，respectively.The propagation time from Wusheng station，Luoduxi station and Xiaoheba station to Beibeistation is 16.9 h，19.0 h and 20.0 h，respectively.The propagation time from Zhutuo station and Beibei station to Cuntan station is 20.8 h and 7.7 h，respectively.When the discharge of the reach from Xiangjiaba to Zhutuo station is less than 10 000 m3/s，the average propagation time of sand peak is 48.0 h.When the discharge is in the range of 10 000～20 000 m3/s，the average propagation time of sand peak is 32.7～46.3 h.When the discharge is in the range of 20 000～30 000 m3/s，the average propagation time of sand peak is 27.0～32.7 h.When the discharge is more than 30 000 m3/s，the average propagation time of sand peak is less than 27.0 h.

Key words：

sand peak transport;propagation time;frequency distribution;numerical simulation;Three Gorges Reservoir;upper reaches of changjiang River