李珍 江穎 安靜泊

摘 要：2018—2019年汛期，黃河來水量大，小浪底水庫(kù)連續(xù)兩年實(shí)施了低水位、長(zhǎng)歷時(shí)、大流量、高含沙的排沙運(yùn)用，取得了顯著的排沙效果?；谛±说姿畮?kù)實(shí)際排沙運(yùn)用情況，分析水庫(kù)排沙影響因素，采用回歸分析方法建立水庫(kù)排沙比計(jì)算公式，結(jié)果表明庫(kù)區(qū)蓄水量、壩前最小水深、入庫(kù)流量、水平回水長(zhǎng)度等與水庫(kù)排沙比關(guān)系較為密切，相關(guān)系數(shù)均在0.80以上，經(jīng)實(shí)測(cè)資料驗(yàn)證，排沙比計(jì)算公式可較好地模擬水庫(kù)排沙情況。進(jìn)一步結(jié)合壩區(qū)地形觀測(cè)資料，分析近壩區(qū)泥沙沖淤變化。2014年以來，受泄流底孔排沙影響，距壩4.5 km范圍內(nèi)河床比降穩(wěn)定在4.3‰～5.1‰之間;進(jìn)水塔前泥沙淤積最低高程在176～181 m之間變化，當(dāng)前左岸泥沙淤積高程逐步抬升，坡度逐漸變陡，尚未形成穩(wěn)定邊坡，而右岸受山體控制，坡度大，泥沙淤積量少。

關(guān)鍵詞：排沙比;敞泄排沙;庫(kù)壩區(qū)淤積;小浪底水庫(kù)

中圖分類號(hào)：TV697.1;TV882.1

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

doi：10.3969/j.issn.1000-1379.2021.09.006

引用格式：李珍，江穎，安靜泊.2018—2019年汛期小浪底水庫(kù)排沙運(yùn)用分析[J].人民黃河，2021，43（9）：32-37.

Analysis of the Sediment Discharging Effect of Xiaolangdi Reservoir in 2018-2019 Flood Seasons

LI Zhen， JIANG Ying， AN Jingbo

（Xiaolangdi Multi-Purpose Dam Project Management Center， Ministry of Water Resources， Zhengzhou 450000， China）

Abstract： In the flood seasons of 2018 and 2019，the inflow of Xiaolangdi Reservoir was higher than before and it was the first time that the operation of low reservoir level， long duration，large discharge， and high sediment concentration was applied for two consecutive years， which effectively scoured the sediment in the reservoir area.? This paper analyzed the influencing factor of sediment discharge and established the calculation formula of the reservoir sediment discharge ratio by the regression analysis method， based on the actual sediment discharge operation of the reservoir. The results point out that pondage， lowest water depth and inflow are most closely related to sediment discharging， with a high correlation above 0.80. The calculation formula for the sediment discharge ratio verified by the measured data can simulate the sediment discharge of the reservoir well. Based on the observation data of the topography of the reservoir dam area， the sediment scour and silting change in the dam area were analyzed， and it was proposed that the river bed gradient within 4.5 km of the reservoir area had been stabilized between 4.3‰ and 5.1‰ since 2014， under the influence by the discharge bottom hole. The lowest elevation of sediment deposition in front of the intake tower varies between 176 m and 181 m. On the right bank in front of the dam， the sediment deposition elevation gradually rises and the slope becomes steeper. Currently， a stable slope has not been formed. As the left bank in front of the dam is close to the mountain with a high slope and there is little sediment accumulation.

Key words： suspended sediment discharge ratio; open drain sediment; sediment deposition in reservoir dam area; Xiaolangdi Reservoir

小浪底水庫(kù)位于黃河中游最后一個(gè)峽谷河段，壩址距下游花園口水文站128 km，處于控制黃河下游來水來沙的關(guān)鍵部位，壩址以上控制流域面積為69.4萬km2，占黃河流域總面積的92.3%，是解決黃河下游防洪減淤等問題不可替代的關(guān)鍵性工程。工程開發(fā)任務(wù)以防洪（包括防凌）減淤為主，兼顧供水、灌溉、發(fā)電，除害興利，綜合利用[1]。小浪底水庫(kù)自1999年10月下閘蓄水以來，取得了巨大的社會(huì)、經(jīng)濟(jì)和生態(tài)效益。截至2020年4月，小浪底水庫(kù)已累計(jì)淤積泥沙32.86億m3，處于攔沙后期第一階段。水庫(kù)運(yùn)用以蓄水?dāng)r沙和調(diào)水調(diào)沙運(yùn)用為主，如何科學(xué)調(diào)控水沙，減緩庫(kù)區(qū)泥沙淤積，長(zhǎng)期發(fā)揮以防洪減淤為主的綜合利用效益是當(dāng)前小浪底水庫(kù)運(yùn)用面臨的重大技術(shù)問題。

2018—2019年汛期，黃河上中游地區(qū)降水偏多，干流洪水洪量相對(duì)較大、持續(xù)時(shí)間較長(zhǎng)。黃河防汛抗旱總指揮部抓住有利時(shí)機(jī)，積極應(yīng)對(duì)，實(shí)施小浪底水庫(kù)降低水位排沙運(yùn)用，充分利用洪水入庫(kù)過程沖刷庫(kù)區(qū)淤積的泥沙，實(shí)現(xiàn)了庫(kù)容恢復(fù)和河道多排沙的目標(biāo)。2018—2019年汛期降低水位排沙運(yùn)用期間，小浪底水庫(kù)排沙量分別為3.63億、4.67億t，排沙比分別達(dá)到237%、346%，無論排沙量還是排沙比均高于以往歷次排沙運(yùn)用。

筆者通過研究2018—2019年排沙運(yùn)用情況，結(jié)合以往水庫(kù)排沙資料，分析小浪底水庫(kù)排沙影響因素，并以此為基礎(chǔ)建立排沙比計(jì)算公式，同時(shí)研究水庫(kù)排沙對(duì)壩區(qū)漏斗淤積形態(tài)的影響，以期為今后小浪底水庫(kù)調(diào)度運(yùn)用提供參考。

1 2018—2019年汛期小浪底水庫(kù)排沙運(yùn)用情況

1.1 2018年汛期

2018年汛期，黃河流域降水量較常年偏多24%，強(qiáng)降水主要集中于上中游的龍劉區(qū)間、劉蘭區(qū)間、蘭托區(qū)間和山陜區(qū)間北部。受降水影響，黃河上游共出現(xiàn)3次編號(hào)洪水[2]，中游涇渭河、山陜區(qū)間及黃河下游大汶河也相繼出現(xiàn)了明顯的洪水過程。

為應(yīng)對(duì)渭河洪水和上游來水，在洪水來臨前，7月3日小浪底水庫(kù)開始加大下泄流量，騰庫(kù)迎洪，降低庫(kù)水位排沙減淤。7月4日5時(shí)30分，水庫(kù)異重流開始排沙出庫(kù)，7月13日庫(kù)水位最低降至211.77 m（相應(yīng)蓄水量2.0億m3），低于淤積三角洲高程10.59 m，出庫(kù)含沙量明顯增大，最大出庫(kù)含沙量達(dá)到369 kg/m3?？紤]到進(jìn)入下游的高含沙水流可能導(dǎo)致下游洪峰增大，高含沙水流出庫(kù)期間適當(dāng)減小了出庫(kù)流量。隨著入庫(kù)洪水流量減小，小浪底水庫(kù)7月22日后逐級(jí)減小下泄流量，7月27日減小至1 000 m3/s，水庫(kù)開始逐步回蓄，降水排沙階段結(jié)束。8月20日前按不超過汛限水位230 m運(yùn)用，8月21日起向后汛期汛限水位過渡。

7月4—27日小浪底水庫(kù)降水位排沙運(yùn)用期間，水庫(kù)入庫(kù)沙量為1.54億t，排沙比為265.9%，平均出庫(kù)含沙量為59.3 kg/m3。整個(gè)汛期（7—10月，下同）小浪底水庫(kù)入庫(kù)水量為241.5億m3，入庫(kù)沙量為4.84億t，出庫(kù)水量為221.9億m3，出庫(kù)沙量為4.66億t，水庫(kù)總體呈現(xiàn)微淤，淤積量0.18億t。2018年汛期小浪底水庫(kù)進(jìn)出庫(kù)水沙及庫(kù)水位變化過程見圖1。

1.2 2019年汛期

2019年汛期，黃河流域暴雨過程頻繁，強(qiáng)降水落區(qū)重疊度高，雨區(qū)主要集中在黃河源區(qū)、涇渭洛河、山陜區(qū)間。黃河源區(qū)降水持續(xù)時(shí)間長(zhǎng)，降水總量大，山陜區(qū)間降水主要發(fā)生在7月下旬至8月上旬，涇渭洛河降水主要發(fā)生在6月下旬至8月初。受降水影響，汛期流域洪水過程多、歷時(shí)長(zhǎng)、發(fā)生時(shí)間早，上、中、下游河道洪水洪量較大，編號(hào)洪水出現(xiàn)4次，其中上游3次、中游1次[3]。

為應(yīng)對(duì)渭河洪水和上游來水，小浪底水庫(kù)汛前合理調(diào)控出庫(kù)流量，于6月底降至汛限水位235 m。為增強(qiáng)水庫(kù)排沙效果，水庫(kù)7月1日后繼續(xù)大流量下泄，實(shí)施降水位排沙運(yùn)用，下泄控制流量最大為4 200 m3/s。7月20日庫(kù)水位降至210 m附近后運(yùn)行至8月1日，小浪底水庫(kù)按進(jìn)出庫(kù)平衡運(yùn)用，最低庫(kù)水位為7月26日的209.12 m。8月1日20時(shí)，水庫(kù)開始按流量1 400 m3/s下泄，庫(kù)水位逐步升高，至8月5日庫(kù)水位漲至220 m左右。為應(yīng)對(duì)中游洪水及九號(hào)臺(tái)風(fēng)“利奇馬”可能帶來的影響，8月5日水庫(kù)再次實(shí)施降水位運(yùn)用，至8月9日庫(kù)水位降至213 m附近，而后根據(jù)來水情況，逐步減小下泄流量，8月12日下泄流量減小至500 m3/s，8月21日庫(kù)水位回升至233.96 m，開始向后汛期汛限水位過渡。9月中下旬，為應(yīng)對(duì)第3號(hào)洪水，水庫(kù)又進(jìn)行了防洪運(yùn)用。

7月14日—8月12日水庫(kù)降水位排沙運(yùn)用期間，水庫(kù)入庫(kù)沙量為1.23億t，排沙比為325.6%，平均出庫(kù)含沙量為42.9 kg/m3。整個(gè)汛期小浪底水庫(kù)入庫(kù)水量為220.9億m3，入庫(kù)沙量為2.74億t，出庫(kù)水量為220.0億m3，出庫(kù)沙量為5.45億t，水庫(kù)沖刷2.71億t。2019年汛期小浪底水庫(kù)進(jìn)出庫(kù)水沙及庫(kù)水位變化過程見圖2。

2 水庫(kù)排沙影響因素分析

2.1 水庫(kù)排沙影響因素相關(guān)性分析

水庫(kù)排沙與入庫(kù)水沙、庫(kù)區(qū)地形及水庫(kù)運(yùn)用方式等密切相關(guān)。一般用出庫(kù)沙量與入庫(kù)沙量的比值即排沙比η表示水庫(kù)排沙效果，影響水庫(kù)排沙比的因素眾多，如入庫(kù)流量、出庫(kù)流量、庫(kù)區(qū)蓄水量、壩前運(yùn)用水位、前期淤積量等。以往專家開展了大量的研究，如涂啟華等[4]、張啟舜等[5]基于實(shí)測(cè)資料，得到水庫(kù)排沙比與水庫(kù)蓄水量成反比，與入庫(kù)流量成正比;焦恩澤等[6]綜合考慮進(jìn)出庫(kù)流量、入庫(kù)含沙量、水面坡降、庫(kù)容等影響因素，提出同時(shí)適用于水庫(kù)壅水排沙和明流排沙的排沙比公式;黃仁勇等[7]利用逐步回歸法研究表明V/Q（水庫(kù)蓄水量與出庫(kù)流量之比）為出庫(kù)排沙比的主要影響因素;張帥等[8]分析了小浪底水庫(kù)汛期排沙比與入庫(kù)平均含沙量、進(jìn)出庫(kù)平均流量、平均庫(kù)容等影響因素之間的相關(guān)程度;王婷等[9]通過分析水庫(kù)降水沖刷效果影響因素，認(rèn)為適時(shí)降水沖刷可以合理處理長(zhǎng)期保持有效庫(kù)容和減少黃河下游河道淤積的關(guān)系。

統(tǒng)計(jì)2004年以來小浪底水庫(kù)調(diào)水調(diào)沙期間降低水位排沙時(shí)水庫(kù)運(yùn)用指標(biāo)和水沙特征值，見表1;采用最小二乘法分析各影響因素與排沙比的相關(guān)系數(shù)，見圖3?？梢钥闯?，水庫(kù)排沙期庫(kù)區(qū)有效蓄水量、壩前平均水深、壩前最小水深、入庫(kù)流量、水平回水長(zhǎng)度、前期淤積量等因素對(duì)水庫(kù)排沙影響較大，其中有效蓄水量、壩前最小水深、水平回水長(zhǎng)度、入庫(kù)流量與水庫(kù)排沙比的相關(guān)系數(shù)均在0.80以上。

當(dāng)前小浪底水庫(kù)主要以異重流形式排沙，水庫(kù)回水末端斷面流量、含沙量過程是異重流形成并持續(xù)運(yùn)行的動(dòng)力條件，在一定的庫(kù)區(qū)河床邊界條件下，水庫(kù)運(yùn)用水位越低，回水末端以上庫(kù)區(qū)自然河段就越長(zhǎng)，調(diào)水調(diào)沙期三門峽水庫(kù)下泄大流量“清水”沖刷恢復(fù)的含沙量就越高，形成的水庫(kù)異重流挾沙能力強(qiáng)且能量大，加之水庫(kù)蓄水較小，異重流向壩前運(yùn)動(dòng)過程中能量消耗不大，排沙出庫(kù)的含沙量相應(yīng)就越高。若后續(xù)洪水流量大、持續(xù)時(shí)間長(zhǎng)，則前期形成的異重流得到源源不斷的動(dòng)力補(bǔ)充，水庫(kù)排沙效果就更明顯。

2018年、2019年小浪底水庫(kù)排沙效果明顯，主要是因?yàn)楫愔亓髋派称陂g小浪底水庫(kù)運(yùn)用水位較低，平均水位分別為218.30、215.99 m（相應(yīng)壩前平均水深分別為43.30、40.99 m），最低水位為211.80、209.00 m（相應(yīng)壩前最小水深分別為36.80、34.00 m），加之后期洪水時(shí)間長(zhǎng)，異重流排沙的后續(xù)動(dòng)力充足，因此水庫(kù)排沙比、排沙量均居歷次前列。2012年、2013年調(diào)水調(diào)沙期間小浪底水庫(kù)運(yùn)用水位較低，水庫(kù)排沙比也相對(duì)較大，但入庫(kù)水量較小，異重流后續(xù)動(dòng)力不足，水庫(kù)排沙量相對(duì)較小。

2.2 水庫(kù)排沙比計(jì)算公式

選取與小浪底水庫(kù)排沙密切相關(guān)的庫(kù)區(qū)有效蓄水量、入庫(kù)流量等因素，建立水庫(kù)排沙比計(jì)算公式，對(duì)排沙比影響較大的壩前最小水深、水平回水長(zhǎng)度與庫(kù)區(qū)有效蓄水量直接或間接相關(guān)，對(duì)排沙比的影響可通過庫(kù)區(qū)有效蓄水量表達(dá)。采用回歸分析方法建立水庫(kù)排沙比的計(jì)算公式如下：

η=1.808（Q入V有效）0.737（1）

式中：η為排沙比，%;V有效為水庫(kù)排沙期有效蓄水量，億m3;Q入為入庫(kù)流量，m3/s。

利用歷年排沙資料對(duì)式（1）進(jìn)行驗(yàn)證，排沙比公式計(jì)算值與實(shí)測(cè)值對(duì)比見圖4，計(jì)算值與實(shí)測(cè)值的擬合度達(dá)0.88，故可以認(rèn)為本排沙比公式可較好地模擬水庫(kù)排沙情況。

根據(jù)水庫(kù)排沙比計(jì)算公式可以看出，要取得較好的排沙效果，小浪底水庫(kù)排沙期間要提前預(yù)泄，降低庫(kù)水位，在洪水入庫(kù)前保持較小的蓄水量，使異重流更易潛入并較好地排泄出庫(kù)。當(dāng)然，水庫(kù)異重流的排沙效果還與庫(kù)區(qū)地形（如比降、庫(kù)區(qū)斷面形態(tài)）、入庫(kù)懸移質(zhì)泥沙組成、排沙期孔洞組合調(diào)度及泄流量等因素有關(guān)，本次針對(duì)小浪底水庫(kù)提出的水庫(kù)排沙比計(jì)算公式可為近期水庫(kù)排沙調(diào)度方案制定提供參考。今后，隨著水庫(kù)泥沙淤積，河道斷面形態(tài)發(fā)生明顯變化時(shí)，水庫(kù)排沙特性也可能發(fā)生相應(yīng)變化。

3 壩區(qū)漏斗形態(tài)變化分析

隨著小浪底水庫(kù)庫(kù)區(qū)泥沙淤積量增加，淤積三角洲頂點(diǎn)逐步向壩前推進(jìn)，淤積面不斷抬高，樞紐泄水孔洞泥沙淤堵的風(fēng)險(xiǎn)越來越大。水庫(kù)排沙比可以反映庫(kù)區(qū)沖淤狀態(tài)，但無法表達(dá)近壩區(qū)沖淤變化過程，有必要進(jìn)一步分析壩區(qū)泥沙沖淤變化。小浪底水庫(kù)距壩4.5 km范圍內(nèi)設(shè)置了31個(gè)漏斗測(cè)驗(yàn)斷面，距壩210 m以內(nèi)斷面間距為50 m，距壩210～710 m范圍內(nèi)斷面間距為100 m，距壩710 m以上的漏斗區(qū)范圍內(nèi)斷面間距為100～200 m。

3.1 縱剖面形態(tài)變化

小浪底水庫(kù)壩區(qū)縱剖面形態(tài)變化見表2、圖5、圖6。隨著小浪底庫(kù)區(qū)泥沙淤積量增加，淤積三角洲頂點(diǎn)逐步向壩前推進(jìn)，淤積三角洲前坡段（距壩約10 km范圍內(nèi)）比降變大，2019年10月比降達(dá)到4.7‰。壩區(qū)4.5 km范圍內(nèi)河床比降也逐步變大，受泄流底孔排沙影響，2014年以來比降基本穩(wěn)定在4.3‰～5.1‰。

2018年、2019年水庫(kù)近壩段（距壩15 km范圍內(nèi)）縱剖面變化見圖7和圖8。受汛期降低水位排沙影響，2018年汛期水庫(kù)淤積三角洲頂點(diǎn)向壩前推進(jìn)3.7 km，距壩4.5 km范圍內(nèi)水庫(kù)縱剖面變化不大，4.5 km以上至三角洲頂點(diǎn)河段河床發(fā)生明顯淤積。2019年汛期水庫(kù)淤積三角洲頂點(diǎn)位置未發(fā)生變化，距壩4.5 km范圍內(nèi)河床發(fā)生明顯沖刷，沖刷深度2～3 m，4.5 km以上至三角洲頂點(diǎn)河段河床發(fā)生淤積。

3.2 橫斷面形態(tài)變化

HHLD01斷面（距壩60 m）形態(tài)變化見圖9，進(jìn)水塔前斷面形態(tài)變化見圖10。可以看出，受水庫(kù)泄流排沙影響，進(jìn)水塔前泥沙淤積最低高程在176～181 m之間變化。進(jìn)水塔左岸壩前泥沙逐步淤積抬升，泥沙淤積高程已達(dá)205 m，隨著壩前淤積厚度增加，坡度逐漸變陡，尚未形成穩(wěn)定的邊坡。進(jìn)水塔右岸受山體控制，坡度大，泥沙淤積量少。

2018—2019年HHLD01斷面形態(tài)變化見圖11和圖12。2018年汛期水庫(kù)降低水位排沙運(yùn)用，壩前含沙量較高，泥沙在壩前明顯落淤，平均淤積厚度達(dá)到10 m，進(jìn)水塔前泥沙淤積最低高程為178.1～180.2 m。2019年汛期壩前泥沙平均淤積厚度為5 m左右，進(jìn)水塔前泥沙淤積最低高程為176.1～176.9 m。

4 結(jié) 論

（1）2018—2019年汛期，黃河上中游地區(qū)降水偏多，干流河道洪水洪量相對(duì)較大、持續(xù)時(shí)間較長(zhǎng)，小浪底水庫(kù)實(shí)施了降低水位排沙運(yùn)用，充分利用入庫(kù)洪水過程沖刷庫(kù)區(qū)淤積的泥沙，取得了較好的排沙效果，有效減緩了庫(kù)區(qū)泥沙淤積。2018—2019年汛期降水位排沙運(yùn)用期間，小浪底水庫(kù)排沙量分別為3.63億t、4.67億t，排沙比分別達(dá)237%、346%。

（2）分析水庫(kù)排沙期庫(kù)區(qū)有效蓄水量、壩前平均水深、壩前最小水深、入庫(kù)流量、水平回水距離、前期淤積量等因素對(duì)水庫(kù)排沙的影響，表明庫(kù)區(qū)有效蓄水量、壩前最小水深、水平回水長(zhǎng)度、入庫(kù)流量與水庫(kù)排沙比關(guān)系較為密切，相關(guān)系數(shù)均在0.80以上。采用回歸分析方法建立了排沙比計(jì)算公式，可為近期水庫(kù)排沙調(diào)度方案制定提供參考。

（3）分析庫(kù)區(qū)壩前泥沙淤積形態(tài)的變化，結(jié)果表明，隨著小浪底庫(kù)區(qū)泥沙淤積量增加，淤積三角洲頂點(diǎn)逐步向壩前推進(jìn)，淤積三角洲前坡段逐年變陡，2019年10月比降達(dá)到4.7‰。距壩4.5 km范圍內(nèi)河床比降也逐步變大，受泄流底孔排沙影響，2014年以來比降穩(wěn)定在4.3‰～5.1‰。進(jìn)水塔前泥沙淤積最低高程在176～181 m之間變化，進(jìn)水塔左岸壩前泥沙逐步淤積抬升，坡度逐漸變陡，目前尚未形成穩(wěn)定的邊坡。進(jìn)水塔右岸受山體控制，坡度大，泥沙淤積量少。

參考文獻(xiàn)：

[1] 李曉宇，李焯，郭銀.小浪底水庫(kù)2018年主汛期防洪預(yù)泄排沙效果分析[J].人民黃河，2019，41（12）：13-15，19.

[2] 魏軍，任偉，楊會(huì)穎，等.2018年汛期黃河水沙調(diào)度實(shí)踐[J].人民黃河，2019，41（5）：1-4，8.

[3] 水利部黃河水利委員會(huì).2019年汛期黃河洪水技術(shù)總結(jié)[R].鄭州：水利部黃河水利委員會(huì)，2019：12.

[4] 涂啟華，楊賚斐.泥沙設(shè)計(jì)手冊(cè)[M].北京：中國(guó)水利水電出版社，2006：133-155.

[5] 張啟舜，張振秋.水庫(kù)淤積形態(tài)及其過程的計(jì)算[J].泥沙研究，1982，7（1）：1-12.

[6] 焦恩澤，林斌文.水庫(kù)淤積的簡(jiǎn)化估算方法[J].人民黃河，1982，4（1）：9-15.

[7] 黃仁勇，談廣鳴，范北林.三峽水庫(kù)蓄水運(yùn)用后汛期洪水排沙比初步研究[J].水力發(fā)電學(xué)報(bào)，2013，32（5）：129-133，152.

[8] 張帥，夏軍強(qiáng)，李濤.小浪底水庫(kù)汛期排沙比研究[J].人民黃河，2018，40（1）：7-11.

[9] 王婷，馬迎平，張俊華，等.小浪底水庫(kù)降水沖刷效果影響因素試驗(yàn)研究[J].人民黃河，2014，36（8）：4-6.

【責(zé)任編輯 張 帥】