水利工程河道泥沙沖淤規律及水沙調節仿真

李東陽, 倪菲菲

（河南黃河勘測規劃設計研究院有限公司, 河南鄭州 450003）

隨著數字化、物聯網、通信技術的不斷發展, 人們對環境監控及管理的要求不斷提高, 實現精細化水利管理已成為河流治理發展的趨勢。研究河流泥沙沖淤[1-2］的發展過程及其形式, 對河口河流防洪, 提出合理的沖淤治理措施具有重要意義。

目前, 眾多學者已經針對泥沙沖淤進行了相關研究, 并取得了豐碩的成果。一方面, 有學者針對泥沙輸移模型[3］、細粒沉積模型[4］、沖淤過程[5］、調節水沙關系模型[6］等進行研究, 為多泥沙河流改善治理提供了理論支持。另一方面, 也有許多學者針對調節水沙關系進行實踐研究。這些研究的操作方法在內容和形式上都有所不同, 但都充分利用了下游河道在高流量河段的輸沙能力。同時, 調節儲存和排放的時間和體積等, 將盡可能多的泥沙通過下游河道從水庫輸送到大海。

本文以黃河下游為例, 對小浪底水庫的運行方式進行了優化, 以提高下游多泥沙河流的輸沙能力, 延長水庫的使用壽命。

1 相關理論

1.1 泥沙沖淤響應方程

多泥沙河流泥沙沖淤響應方程如下所示：

式（1）中, Wslj表示河流的泥沙流出量, Wxhw表示河流的來水量, Wsxhw表示河流的泥沙流入量。k1和k2為系數, c為一個常數。

然而, 多泥沙河流的引水往往伴隨著泥沙的分流, 因此沖淤計算也必須考慮到引水的影響。選定河流的侵蝕量或沉積量可通過以下方程計算：

式（2）中, Wslit為侵蝕量或沉積量, 其中沉積量為正值, 侵蝕量為負值。Wsin是與引水相關的泥沙分流量。Wsxhw和Wslj與式（1）相同。

1.2 輸沙平衡方程

輸沙平衡方程可描述為

式（3）中, Jp為平衡輸沙梯度；Qs為汛期平均輸沙率；d50為采用粒徑劑法測定的河道懸移質泥沙粒徑中值, 其單位為mm；n為河流糙率；B和h分別為水面寬度和水深, 且對應于汛期平均流量Q（單位為m）。一般情況下, 對于含沙河床, 有B=38.6Q0.31, h=0.081Q0.44。此外, 對于大多數多泥沙河流而言, k為90～160之間的經驗系數（黃河流段通常為140）。

令維持輸沙平衡所需的泥沙量為Qs, 則Qs可通過汛期泥沙量Ws除以汛期持續時間Δt來計算。因此, Ws可計算為

1.3 多泥沙河流系統模型

多泥沙河流系統模型（Silt-laden River System Simulation, SRSS）的數學模型以日懸移質為基礎運行, 考慮了當前水庫的水沙調節效應, 便于模擬含沙河流下游的實際輸沙過程。模擬懸移質的控制方程包括以下4個方程：水流連續性方程、水流運動方程、泥沙連續性方程（按泥沙組）、河床變形方程。各方程計算過程如下：

式（5）-（8）中, Q為水流量；A為橫截面積；Z為水位；J為能量斜率；Sk是泥沙濃度；Ad是侵蝕和沉積的橫截面積；ql是每單位流量長度的橫向流出（流入）, 其中流出為正, 流入為負；Ul是主流方向橫向流出（流入）速度的分量；qs是每單位流量長度的橫向輸沙率, 其中流出為正, 流入為負；γs是沉積物的干容重；Zb為河床高程；α為飽和度恢復系數；ω為沉降速度；S*是水流挾沙能力；X為流量長度；t為時間；k為沉積物群數。進一步可將水流挾沙能力方程描述為

式（9）中, γ為渾水的體積密度；Κ為卡門常數；ωm為混合挾沙能力的沉降速度；K和m分別為挾沙能力系數和指數, 其中K=2.5, m=0.62。

此外, 沉降速度ωm可根據河流流態確定。對于停滯狀態（沉積物直徑d≤0.1 mm）, 有

2 基于SRSS的水沙調節方法

本節將對上游水庫的運行方式進行優化, 以最大限度地減少多泥沙河流下游侵蝕, 同時最大限度地增加輸沙能力和入海輸沙量。含沙河流下游的上游水庫運行規則優化過程包含5個步驟, 具體描述如下：

步驟1收集水庫首次運行以來所記錄的所有設計參數和數據；

步驟2分析水庫蓄水、調水調沙和攔沙對河道幾何形態和輸沙效率的影響；

步驟3提出上游水庫的精細化運行模式, 以提高下游多泥沙河流的輸沙效率；

步驟4計算上游水庫精細調度模式下下游多泥沙河流的輸沙能力。然后, 將其與維持生態效益所需的泥沙閾值進行比較。如果輸沙量超過閾值, 繼續調整運行方式；否則, 執行步驟5；

步驟5將最新的優化運行模式與現有的和以前的優化運行模式進行比較, 不斷調整輸沙平衡, 確保防洪和生態健康。

3 案例仿真與分析

3.1 研究區域

研究區域為沿著黃河下游從小浪底水庫延伸到渤海。小浪底水庫控制流域面積為69.4萬平方千米, 占黃河流域總面積的92.3%。該水庫的工程開發任務側重于防洪、防冰堵和減淤, 同時考慮到供水、灌溉和發電要求。小浪底水庫總庫容為126.5億立方米, 其中攔沙能力為75.5億立方米。

3.2 泥沙平衡輸移閾值研究

3.2.1 泥沙沖淤響應方程驗證

根據小浪底1960—2015年測得的年輸沙量, 泥沙量計算的回歸分析參數如表1所示。

表1 小浪底泥沙量計算回歸分析參數Tab.1 Regression analysis parameters for Xiaolangdi sediment calculation

通過回歸分析確定挾沙能力系數和指數, 則泥沙沖淤響應方程描述為

此外, 相關系數R2為0.976, 標準差σ為1.12, 表明與測量結果具有較高的相關性。

3.2.2 多泥沙河流系統模型驗證

本研究利用1987—2015年收集的數據重新驗證了相關模型參數。經驗證, SRSS中的飽和恢復系數在侵蝕狀態下為0.02～0.04, 在沉積狀態下為0.002～0.009。值得注意的是, 對于坡度較小和水深較深的斷面, 系數值較大, 反之亦然。對于大多數截面, 系數值為0.005～0.009。通過驗證模型得到的下游沖淤計算結果與圖1所示的實測值進行了比較, 可以看出, 沉積和侵蝕狀態的計算值與實測值一致。說明該模型能可靠地模擬水庫下游河流沖淤平衡的臨界狀態。

圖1 下游沖淤計算及實測對比結果Fig.1 Comparison between the calculation value and the measured value of downstream erosion and deposition

3.3 泥沙調節方法驗證

假定小浪底水庫來沙3億噸。如圖2所示, 在現行運行模式下, 小浪底水庫預計將于2050年停止攔沙, 而優化運行模式預計將使攔沙期延長15年。如圖3所示, 在現行運行模式下, 50年內累計將沉積2.16億噸泥沙, 但在優化運行模式下, 河道將得到更好的沖刷, 可能減少5.43億噸泥沙沉積。

圖2 不同運行方式下游沖淤對比結果Fig.2 Comparison results of downstream erosion and deposition under different operation modes

圖3 不同運行方式下游輸沙對比結果Fig.3 Comparison results of downstream sediment transportation under different operation modes

4 結 論

本文基于多泥沙河流系統對黃河小浪底水庫下游沖淤過程進行了研究與分析, 提出了一種SRSS的泥沙調節方法, 從而最大限度地減少多泥沙河流下游侵蝕, 同時最大限度地增加輸沙能力和入海輸沙量。