小浪底水庫水沙調控對淤積形態影響的數值模擬

假冬冬 江恩慧 王遠見 邵學軍

摘 要：水庫淤積形態是影響水沙調節效率的一項關鍵因素。為優化水庫調度方式，采用考慮細顆粒淤積物流動特性的水庫淤積形態數值模型，開展了小浪底水庫淤積形態對水沙調控響應的模擬分析工作。研究結果表明：三角洲形態及頂點位置隨著水庫的運行調控而發生變化，三角洲頂點附近頂坡段的沖淤調整和水庫運行低水位與三角洲頂點高程之間存在較明顯的關聯性，當水庫低水位低于三角洲頂點高程時三角洲頂坡段出現沖刷，當水庫低水位高于三角洲頂點高程時三角洲頂坡段出現淤積;淤積形態為同等淤積量的錐體時，庫區上段受河道邊界影響有沖有淤，中下段庫區淤積明顯，且淤積量較三角洲淤積形態的大;考慮人工清淤措施時，清淤量與水庫淤積總量相比占比非常小，因此淤積形態總體變化與不考慮人工清淤時基本類似，僅在清淤疏浚部位及附近局部河段有一定變化。

關鍵詞：淤積形態;水沙調控;數值模擬;小浪底水庫

中圖分類號：TV856;TV882.1

文獻標志碼：A

doi：10.3969/j.issn.1000-1379.2022.02.007

引用格式：假冬冬，江恩慧，王遠見，等.小浪底水庫水沙調控對淤積形態影響的數值模擬[J].人民黃河，2022，44（2）：32-35，44.

Abstract： Sedimentation pattern is one of the critical factors， which has great impact on the efficiency of water and sediment regulation in reservoirs. In order to improve the reservoir operation， the responses of sedimentation patterns of the Xiaolangdi Reservoir to water and sediment regulation were simulated by a mathematical model considering the effect of fine-grained sediment deposits movement. The simulation results indicate that the delta sedimentation pattern is adjusted during the process of reservoir operation. The variation pattern of the top of the delta depends on the relationship between the lowest water level of reservoir operation and the top elevation of the delta. Erosion occurs when the lowest water level is lower than the elevation of delta top， otherwise deposition will be occurred. Compared with the delta deposition morphology， the sedimentation volume of cone deposition morphology is larger with the same method of water and sediment regulation. The variation of sedimentation pattern considering dredging in the vicinity of the top of the delta is very small， except for the location of dredging.

Key words： sedimentation pattern;water and sediment regulation;numerical simulation;Xiaolangdi Reservoir

基于水庫樞紐工程的水沙關系調節，是保障黃河長久安瀾的重要手段，而水庫淤積形態則是影響水沙調節效率的一項關鍵因素。水庫淤積形態是水庫泥沙運動的結果[1]，與入庫水沙條件、水庫自身特點及其運行調度方式等因素相關。與此同時，水庫淤積形態又會影響庫區水沙輸移過程以及庫容分布，影響水庫設計功能的正常發揮。深入研究水沙調控對水庫淤積形態的影響過程與規律，不僅可以豐富水庫泥沙研究內容，也可以為水庫優化調度提供科技支撐，具有重要的理論價值和實際意義。

關于水庫淤積形態的不同類型與判別方法，我國學者提出了水庫三角洲淤積形態的具體判別方法，還在理論上給出了錐體淤積剖面近似于直線的結論，并分別提出了三角洲、錐體及帶狀3種典型淤積形態的判別方法[1-2]。小浪底水庫是黃河干流水沙調控的重要控制性工程，針對其水庫泥沙問題，許多學者開展過大量研究工作。在異重流潛入條件判別方面，李書霞等[3]通過改進描述異重流運動的動量方程，給出了潛入條件判別的新計算公式，并采用多組實測資料對該判別條件進行率定與驗證。在水庫降水庫區溯源沖刷方面，李濤等[4]利用小浪底水庫的實體模型，對溯源沖刷過程開展了多個組次的試驗，分析了侵蝕基準面、水量等因素與水庫溯源沖刷量之間的關系。張俊華等[5]通過建立小浪底水庫準二維數學模型，對水庫運用初期的淤積過程進行了模擬分析。在水庫排沙比方面，張帥等[6]基于小浪底水庫實測資料，研究了汛期排沙比與進出庫平均流量、入庫平均含沙量、平均庫容等因素之間的關系，并建立了擬合精度較高的排沙比計算公式。隨著小浪底水庫的持續運行以及實測資料的不斷積累，孫東坡等[7]分析了入庫水沙過程以及庫區泥沙淤積的變化特點。王婷等[8]對小浪底水庫運用至2016年的庫區淤積分布以及粒徑特征等進行了統計分析。張俊華等[9]從干支流倒灌機制、異重流潛入與運動規律、水庫高效輸沙調度原則、異重流高效排沙試驗等方面，總結分析了小浪底水庫高效輸沙理論與調控的關鍵技術。隨著小浪底水庫泥沙研究的深入進行，細顆粒淤積物自身的流動特性因對淤積形態存在較明顯的影響而開始受到關注[10]。水庫淤積形態是影響水沙調控效率的一項重要因素，其對水沙調控，尤其是對泥沙動態調控（如人工清淤等）的響應規律如何，還有待進一步深入研究。為此，筆者采用考慮細顆粒淤積物流動特性的水庫淤積形態數值模型[10]，針對三角洲和錐體兩種典型淤積形態，開展了小浪底水庫淤積形態對典型水沙調控響應的模擬分析，以期為水庫淤積形態演化機理和水庫優化調度的深入研究提供科技支撐。

1 小浪底水庫淤積形態數值模型簡介

1.1 小浪底水庫概況

小浪底水利樞紐工程是黃河水沙調控的關鍵控制性工程，在黃河治理與保護中具有重要戰略地位。樞紐工程位于黃河中游最后一個峽谷段的出口，控制流域面積約69.4萬km2，占黃河流域面積的92.3%，上距三門峽水利樞紐約130 km。小浪底庫區平面形態狹長彎曲，入匯支流較多，如圖1所示。水庫設計正常蓄水位275 m，原始庫容127.5億m3，長期有效庫容51億m3。該水庫主體工程于1994年9月開工，1999年10月開始下閘蓄水，2000年5月正式運行。小浪底水庫運行以后，實測資料表明庫區泥沙淤積顯著，2015年4月與1999年9月相比，水庫累計淤積的泥沙為30.49億m3（斷面法計算），約為設計攔沙庫容的42.1%[6]。

1.2 考慮細顆粒淤積物流動的淤積形態數值模型簡介

小浪底水庫近壩段淤積的泥沙粒徑較小，中值粒徑在0.008～0.015 mm之間[7，10]，此類淤積物未密實前具有明顯的流動性，是影響庫區淤積形態的重要因素之一。因此，對于小浪底水庫淤積形態的準確模擬而言，一方面要考慮傳統的水沙輸移過程，另一方面還要模擬分析已落淤的細顆粒淤積物自身的流動過程[10]。本研究所采用的小浪底水庫數值模型，模擬黃河干流河道長約123 km，即上游HH56斷面到小浪底水庫大壩。模型基于傳統泥沙輸移的三維動力學模型以及細顆粒淤積物運動模式，可較好地模擬水庫細顆粒泥沙淤積形態，通過2010年淤積形態沖淤變化的驗證表明，實測結果與計算結果吻合較好，可用于模擬淤積形態對水沙調控的響應過程，模型情況詳見文獻[10]。

2 水沙調控對小浪底水庫淤積形態的影響

本研究所指的水庫淤積形態為縱向淤積形態，即水庫沿縱向淤積的剖面形態。水庫的淤積形態外形比較復雜，其中三角洲淤積形態和錐體淤積形態是較為常見的兩種類型。水庫泥沙淤積形態與水沙調控方式（運行方式）密切相關，為此，本節采用數學模型計算的研究方法，分析三角洲淤積形態和錐體淤積形態對泥沙動態調控的響應。

小浪底水庫實測淤積形態為三角洲淤積形態，錐體淤積形態設置則以實際水庫泥沙淤積量為基礎，在同等淤積量且庫尾比降與三角洲淤積形態相當的前提下概化出錐體淤積形態。

2.1 計算水沙條件

不同水沙調控方式對水庫淤積形態的影響有所不同。結合2010年和2013年汛期小浪底水庫水沙調控過程，計算和分析淤積形態對水沙調控的響應。2010年汛期入庫平均流量為1 127 m3/s、入庫平均含沙量為29.27 kg/m3，2013年汛期入庫平均流量為1 640 m3/s、入庫平均含沙量為22.65 kg/m3，對應兩個年度的汛期壩前水位調控過程如圖2所示。

2.2 三角洲淤積形態響應分析

2010年水庫調度運行方式的三角洲淤積形態變化如圖3所示，由圖3可見，水庫的調度運行使三角洲淤積形態及其頂點位置發生變化。汛前（4月），三角洲頂點位置在距壩約24 km處，高程約為219.6 m。汛期水庫壩前水位降低運行（最低水位約211.6 m），三角洲洲面出現較明顯的沖刷，壩前段與前坡段則出現明顯的淤積。汛后（10月），三角洲頂點向下游推進到距壩約19.6 km的位置，三角洲頂點高程較汛前高程降低約4.0 m，為215.6 m。

2013年水庫調度運行方式的三角洲淤積形態變化如圖4所示，由圖4可見，隨著水庫的調度運行，三角洲淤積形態及其頂點位置發生變化。汛前（4月），三角洲頂點位置在距壩約14 km處，高程約為210.0 m。汛期水庫壩前水位降低運行（最低水位約212.8 m），水庫上游河段產生沖刷，下游河段產生淤積。汛后（10月），三角洲頂點位置沒有發生明顯變化，仍在距壩14 km處附近，三角洲頂點高程約為214.0 m，較汛前高程淤高約4.0 m。

水沙調度過程分析表明：小浪底水庫三角洲淤積形態在汛期水沙調控階段存在周期性的沖淤過程。水沙調控初始階段，壩前水位降低，庫區上游河段開始產生沖刷，下游河段產生淤積。當水位繼續下降時，沖刷河段逐漸向下游發展，三角洲前坡段向前淤積推進。三角洲頂點附近頂坡段的沖淤調整和水庫運行低水位與三角洲頂點高程之間存在較明顯的關聯性，當水庫低水位低于三角洲頂點高程時三角洲頂坡段出現沖刷（如2010年汛前三角洲頂點高程約為219.6 m，汛期壩前最低水位約為211.6 m），當水庫低水位高于三角洲頂點高程時三角洲頂坡段出現淤積（如2013年汛前三角洲頂點高程約為210.0 m，汛期壩前最低水位約為212.8 m）。

此外，從2010年典型斷面沖淤變化（見圖5）來看，近壩段模擬的斷面淤積形態均較為平坦，主要原因是近壩段淤積泥沙顆粒較細，此類淤積物在密實前流動性較強，平衡坡降很小，數值模擬過程中對細顆粒淤積物自身的失穩流動特性進行了考慮，這一模擬結果與實際觀測結果基本一致。對于淤積物粒徑相對較大的庫區上段，數值模擬過程中不考慮淤積物自身的流動特性，斷面坡度則略大（如斷面HH44），實際情況與此基本一致。

除了上述傳統水沙調控外，泥沙動態調控還包括人工擾動和輔助清淤等相關措施。本節計算對人工清淤因素進行了考慮，清淤部位不同，對淤積形態的影響也不同，本次計算方案假定在水庫三角洲頂點對應位置進行汛前集中清淤，清淤范圍為三角洲頂點及上游2 km的范圍，清淤厚度為2 m，寬度平均約1.5 km，清淤泥沙約600萬m3。

各工況計算結果基本類似，以2010年三角洲淤積形態為例，分析其淤積形態對泥沙動態調控的響應，模擬結果如圖6所示。由圖6可見，由于人工清淤量與實際水庫淤積量相比占比較小，因此淤積形態總體變化情況與不考慮人工清淤時基本類似，未發生明顯變化，僅在清淤疏浚河段及附近有一定變化，具體表現為：從高程看，以清淤段為中心，河床向壩前和庫尾方向高程略有降低，但影響范圍較小，基本為清淤長度的3倍左右;從局部河床比降變化看，在影響范圍內，清淤上段比降略有變陡，下段坡降變緩。

2.3 錐體淤積形態響應分析

錐體淤積形態設置以2010年4月實際水庫泥沙淤積量為基礎，在同等淤積量且庫尾比降與三角洲淤積形態相當的前提下概化出錐體淤積形態，如圖7所示。

2010年水庫調度運行方式的錐體淤積形態沖淤變化如圖8所示。由于錐體淤積形態的水流流速沿程減小較明顯，因此庫區上段受河道邊界影響有沖有淤，總體沖淤變化較小，中下段庫區淤積明顯，尤其是下段河道較寬處，最大淤積厚度可達10 m，壩前因汛期水位降低運行而略有沖刷。水庫排沙比由三角洲淤積形態時的35%左右降低至錐體淤積形態時的27%左右，因此錐體淤積形態的水庫淤積總量增加約8%。

本節計算對人工清淤因素進行了考慮，人工清淤方案與三角洲淤積形態時一致，即假定在水庫三角洲頂點對應位置進行汛前集中清淤，清淤泥沙同樣約為600萬m3。錐體淤積形態對泥沙動態調控的響應模擬結果如圖9所示。由圖9可見，由于人工清淤量與實際水庫淤積量相比占比非常小，因此淤積形態總體變化情況與不考慮人工清淤時基本類似，未發生明顯變化，僅在清淤疏浚部位及附近局部河段有一定變化，具體表現為：從高程看，以清淤段為中心，河床向壩前和庫尾方向高程略有降低，但影響范圍較小，基本為清淤長度的3倍左右;從局部河床比降變化看，在影響范圍內，清淤上段比降略有變陡，下段坡降變緩。這一變化特征與三角洲淤積形態基本一致。

需要指出的是，本文計算中假定在水庫三角洲頂點對應位置進行汛前集中清淤，而水庫清淤部位不同，對淤積形態的影響可能也會不同，本文僅對此問題進行了初步探討，還需結合水庫泥沙清淤利用的可能性等，開展進一步研究。此外，筆者僅在2010年、2013年汛期調度條件下進行了研究，水庫泥沙運動機理復雜，水庫淤積形態對水沙調控的響應規律還有待進一步研究。

3 結 論

采用考慮細顆粒淤積物流動特性的水庫淤積形態數值模型，針對三角洲和錐體兩種典型淤積形態，開展了小浪底水庫淤積形態對水沙調控響應的模擬分析。

（1）三角洲淤積形態及頂點位置隨著水庫的運行調控而發生變化。隨著汛期壩前水位降低運行，庫區上游河段開始產生沖刷，下游河段產生淤積;當水位繼續下降時，沖刷河段逐漸向下游發展，三角洲前坡段向前淤積推進。三角洲頂點附近頂坡段的沖淤調整和水庫運行低水位與三角洲頂點高程之間存在較明顯的關聯性，當水庫低水位低于三角洲頂點高程時三角洲頂坡段出現沖刷，當水庫低水位高于三角洲頂點高程時三角洲頂坡段出現淤積。

（2）淤積形態為錐體時，水流流速沿程減小較明顯，庫區上段受河道邊界影響有沖有淤，中下段淤積明顯，尤其是下段河道較寬處淤積厚度較大。淤積形態為錐體時淤積量較三角洲淤積形態增加約8%。

（3）考慮人工清淤措施時，由于清淤量與水庫淤積總量相比占比非常小，因此淤積形態總體變化情況與不考慮人工清淤時基本類似，僅在清淤疏浚部位及附近局部河段有一定變化，具體表現為：從高程看，以清淤段為中心，河床向壩前和庫尾方向高程略有降低，但影響范圍較小，基本為清淤長度的3倍左右;從局部河床比降變化看，在影響范圍內，清淤上段比降略有變陡，下段坡降變緩。

參考文獻：

[1] 韓其為.水庫淤積[M].北京：科學出版社，2003：1-30.

[2] 韓其為，楊小慶.我國水庫泥沙淤積研究綜述[J].中國水利水電科學研究院學報，2003，1（3）：5-14.

[3] 李書霞，夏軍強，張俊華，等.水庫渾水異重流潛入點判別條件[J].水科學進展，2012，23（3）：363-368.

[4] 李濤，張俊華，夏軍強，等.小浪底水庫溯源沖刷效率評估試驗[J].水科學進展，2016，27（5）：716-725.

[5] 張俊華，王艷平，張紅武.黃河小浪底水庫運用初期庫區淤積過程數值模擬研究[J].水利學報，2002，33（7）：110-115.

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

[7] 孫東坡，吳默溪.黃河小浪底水庫運用以來的泥沙淤積特征分析[J].華北水利水電大學學報（自然科學版），2018，39（2）：74-79.

[8] 王婷，王遠見，曲少軍，等.小浪底水庫運用以來庫區泥沙淤積分析[J].人民黃河，2018，40（12）：1-3，20.

[9] 張俊華，馬懷寶，夏軍強，等.小浪底水庫異重流高效輸沙理論與調控[J].水利學報，2018，49（1）：62-71.

[10] 假冬冬，王遠見，江恩惠，等.小浪底水庫淤積形態數值模擬[J].水科學進展，2020，31（2）：240-248.

【責任編輯 張 帥】