數值模擬在水庫泥沙淤積管理中的應用分析

陳　偉

(新疆烏蘇市水利局，新疆 烏蘇 833000)

?

數值模擬在水庫泥沙淤積管理中的應用分析

陳偉

(新疆烏蘇市水利局，新疆 烏蘇 833000)

利用實驗方法和數值模擬法對水庫泥沙淤積規律進行了分析，得到了不同時間段各位置的河床高度變化情況。將水庫根據泥沙淤積情況分成特征明顯的三個區域，并分析了這三個區域內的特征，從而為水庫的泥沙淤積治理提供參考。文章比較了Realizable k-ε和標準k-ε兩種數學模型的模擬結果，驗證了利用數值模擬法來預測水庫泥沙淤積情況的可行性，并對兩者的精確度進行了比較。

數值模擬；水庫；泥沙淤積；Realizable k-ε

1　概　述

水庫大多存在泥沙淤積的問題，其相關研究有諸多方面的應用，如保證水庫庫容、預測水庫調水能力、判斷水庫使用壽命等。

目前水庫淤積方面的研究有很多，但由于水流與泥沙運動規律的復雜性，該問題還有很大的研究空間。泥沙淤積的主要研究方法是物理模擬法，即應用較小尺度、且簡化的物理模型來模擬水庫中的泥沙淤積規律，然而，物理模擬法要求的時間及成本很高，因此泥沙淤積的計算研究具有十分重要的實際意義。目前的計算方法主要是一維、二維的數學模型，且大多的數學模型主要是基于概率統計原理而建立，而基于物理意義的三維數值模擬則較少。數值模擬中最常用的模型是標準k-ε模型，而Realizable k-ε模型則在許多問題中被證明要優于標準。因此文章旨在分析數值模擬在水庫泥沙淤積管理中的應用，比較標準和Realizable k-ε模型在該問題中的適用性及預測性能，從而為水庫泥沙淤積管理工作提供科學指導。

2　方　法

測量數據主要在一個長50 m、寬0.5 m，深1.3 m的渠道中獲得。渠道初始河底坡度為0.002， 在渠道尾端設置后擋板，其高度為0.4 m，可保證波形水躍現象與泥沙區末端出現。渠道流量由半徑為10cm的管道來提供、并由流量計來測量，其流量為50 L/s。

利用標準的測針來測量水位及地形變化情況，其誤差范圍約為0.5%，自上游至下游共進行20個位置點的測量，而時間方面則進行四組測量，分別是初始、2 h后，8 h后、以及16h后的狀態。采用不同的泥沙進行實驗，其最下大小為0.12 mm，最大為64 mm，這是許多類似實驗的標準配置方法，泥沙的供應速率為11.45 kg/min。

傳統的水庫泥沙淤積數學模型主要是基于統計原理，其物理意義并不明確，而基于數值模擬技術的三維模擬則被證明是更好的分析工具，它基于物理定律、結果相對精確、且廣泛適用性較強。然而不同的湍流模型卻可能產生差別較大的結果，因此在應用某一模型在一具體問題之前要對該模型進行適用性及精確度的驗證。

目前應用最為廣泛的是標準k-ε模型，但近幾年Realizable k-ε模型也越來越受關注，它在湍流黏度以及能量消散率的計算方面相對標準k-ε更有優勢。在許多問題中Realizable k-ε模型被證明比標準k-ε模型更有優勢，如閆曉惠等人證明了在階梯式溢洪道流的問題分析中前者相對后者更為精確，但兩者在水庫泥沙淤積問題中的比較卻幾乎沒有被進行過，因此文章將對這2種模型在水庫泥沙淤積預測中的適用性進行分析，從而為水庫泥沙淤積管理提供技術參考。

3　結　果

3.1泥沙淤積

不同時間各測量點的河床高度數據，其中橫坐標為水平距離，單位為m，0m處為壩體，50 m處為尾部；縱坐標為河床高度，其中黑實線代表的是初始河床高度(其坡度為0.002)，藍色菱形表示測量值，紅色方塊代表Realizable k-ε模型的模擬值，綠色三角形代表標準k-ε模型的模擬值，見圖1。

圖1(a)　2h后河床高度分布情況

圖1(b)　8h后河床高度分布情況

圖1(a)為2h后河床高度分布情況。由圖可見，在壩體處其河床高度明顯上升，從0.1 m上升至大約0.4 m，其上升幅度約為4倍，說明在該位置處淤積情況嚴重。從0 m到12 m左右，河床高度大體成線性規律下降，在12 m左右回歸到初始高度，之后的位置河床高度基本保持不變。在15～40m范圍內個別位置的河床高度略低于初始高度，說明河床侵蝕情況同時在發生。可以看出，在水庫淤積的初始階段，受影響最大的是壩址處，而之后位置的影響逐漸減少，因此在實際水庫淤積管理中要加強壩址處的清淤工作，而在遠離壩體處的位置則可能發生河床侵蝕現象，因此在實際管理中要對這些位置加強保護。通過比較測量與模擬數據可以看出，3組數據有多處重合，比較接近，因此可以認為這兩種模型都可以很好地預測出該水庫淤積過程，但Realizable k-ε模型的數據更接近與實驗值，可以認為它在該階段相對于標準k-ε更為精確。

圖1(b)表示的是8h后河床高度的分布情況，其變化規律明顯與初始的2 h階段不同。其壩體處的河床高度從大約0.4～0.63 m左右，依然為淤積最高點，但6 h內只上升了20%左右，可見在淤積過程的中期，壩體處的淤積速率有所放緩。從壩體之后各位置的河床高度大體也呈線性變化，但主要分成了2個階段：在0～25 m左右范圍內，河床高度下降的速度相對較慢，即在該范圍內各點的淤積速率較大，而在25～30 m左右的距離則急速下降，說明該范圍內各點處在淤積的過渡狀態，而30 m之后各點的河床高度則很少高于初始高度，說明該范圍內沒有明顯的淤積現象，有的點低于初始河床高度，說明該位置處存在河床侵蝕現象。測量數據與模擬數據的比較證明兩種模型的預測同樣精確，但在淤積范圍內，模型傾向于高估其高度，而在侵蝕范圍內，模型傾向于低估其高度。

圖1(c)　16h后河床高度分布情況

16h后河床高度的分布情況如圖1(c)所示。可見，其壩體處的高度上升至0.88m左右，8h內升高了大約0.25m，其增速不快，同樣可以看出在0～25 m左右范圍內的各點的增速同樣有所放緩，但25～28 m范圍內其增速有所提高而且在28m左右的位置開始高度分布線開始急速下降，但各點高度都高于河床初始高度，說明淤積范圍已擴散至整個研究范圍。模型可以很好地復制實驗結果，證明兩者的精確性，Realizable k-ε模型的模擬值大多高于測量值，而標準k-ε的結果相對于測量數據則略有波動。

3.2模型比較

為更好地比較Realizable和標準k-ε模型，可以用散點圖來進行詳細地分析，如圖2所示。圖中橫坐標為測量值，縱坐標為模擬值，其中藍色菱形標志為Realizable k-ε模型結果、紅色方塊標志為標準k-ε模型結果，黑線為等值線。當一點坐落于等值線上時，說明該值完全等于測量值，相反，當一點遠離等值線是，則該點所代表的模擬值與測量值差別越大。

圖2　模擬值與測量值的比較

由圖2可見，所有的點基本都比較靠近等值線，說明兩個模型在水庫泥沙淤積問題中的預測能力都較強，但相對來說代表Realizable結果的點更接近等值線，因此其精確度較高。大多點坐落于等值線的左上方，說明大多數的模擬值略高于測量值，個別點坐落于右下方，說明這些位置處的模擬值偏低。仔細觀察圖1可以發現，在泥沙淤積范圍內，模型傾向于高估河床高度，而在侵蝕范圍內則傾向于低估河床高度。

為更好地量化兩個模型在水庫泥沙淤積預測中的差異性，計算各個模型在不同階段與測量值的均方根偏差，其結果如表1所示。當均方根偏差值越大時則代表其誤差越大，反之亦然。誤差比由Realizable k-ε模型模擬結果值除以標準k-ε模型的模擬結果而求得，當其值越小代表前者相對于后者精確度越高。

表1　Realizable k-ε和標準k-ε2種模型的均方根偏差值

由表1可見，在所有工況下，Realizable k-ε模型的均方根偏差值皆<標準k-ε模型的值，即誤差比都<1，證明Realizable k-ε模型相對于標準k-ε模型更為精確，但該誤差比隨時間而增長，在16h后Realizable k-ε模型的均方根偏差值只比標準ke的小0.007左右，可見前者依然具有優勢，但不如初始階段明顯。

4　分　析

河床高度最高的始終在壩體處，可見該處淤積情況最為嚴重。由3個階段的數據可見，水庫泥沙淤積在不同階段內其淤積速度最快的點在不斷向遠離壩體方向移動，即淤積范圍在不斷擴散。所有時間段內的淤積都可以分為3個部分：

1)第1個線性范圍可以認為是淤積范圍，在該范圍內淤積速度最快。

2)第2個線性范圍可以認為是過渡范圍，在該范圍內淤積剛開始發生。

3)而之后的部分則可以認為是非淤積范圍，在該范圍內淤積幾乎沒有發生、甚至會有小程度的河床侵蝕想象出現。

在具體的水庫泥沙淤積管理中可以首先判斷這3個范圍，再因時、因地制宜，采用最合理、高效的清淤或者保護手段。

兩種模型都十分接近測量值，可以證明模型在水庫泥沙淤積管理中的應用是可行的。但2個模型都傾向于高估淤積范圍內的河床高度、低估非淤積范圍內的河床高度，即略微高估了淤積和河床侵蝕的程度，原因可能是2個模型都沒有考慮濃度變化對湍流的影響。相對于標準k-ε模型，Realizable k-ε模型更加精確，尤其在初始階段其誤差只約為前者的37%，因此可以代替標準k-ε模型成為水庫泥沙淤積管理中的主要模型。

5　總　結

利用實驗數據與計算數據對水庫泥沙淤積規律進行了分析，證明了數值模擬在該問題中的適用性，為合理、高效地進行水庫泥沙淤積管理提供科學指導并啟發相關技術人員改進計算能力。主要可以得到一下幾點結論:

1)壩體處泥沙淤積情況最嚴重。

2)水庫可以根據泥沙淤積情況可以分為3個部分：淤積范圍、過渡范圍以及非淤積范圍，

3)淤積范圍隨時間而不斷擴散。

4)淤積區的淤積現象與非淤積區的河床侵蝕現象同時存在。

5)利用數值模擬方法進行水庫泥沙淤積分析是可行的。

6)Realizable k-ε模型模型在該問題中比標準k-ε模型更為精確，尤其在淤積的初始階段。

[1]王彥芳，馮民權，尹宏偉，等.汾河水庫泥沙淤積數值模擬[J].水資源與水工程學報,2011(1):101-104.

[2]閆曉惠，陳新，李華煜.基于標準和Realizable kEpsilon湍流模型的階梯式溢洪道流的數值模擬[J].水利科技與經濟,2015(10):29-31.

[3]朱恒峰，康慕誼，趙文武，等.水利水保措施對延河流域侵蝕、泥沙輸移和沉積的影響[J].水土保持研究,2007(14):1-4.

[4]錢雪梅.新疆克孜爾水庫泥沙淤積現狀及水庫排沙調度分析[J].水利科技與經濟.2013(02) :61-63.

[5]柴曉利.石佛寺水庫低水位蓄水帶來的泥沙淤積問題[J].東北水利水電，2012(06) :66-67.

Application and Analysis of Value Simulation in Reservoir Sediment Deposition Management

Chen Wei

(Wusu City Water Conservancy Bureau, Wusu 833000,China)

Experimental method and valuable simulation method were used to analyze the reservoir silt sedimentation law, the height change of river bed during various time intervals were obtained. Based on the situation of silt sedimentation, reservoir was divided into three areas with obvious characteristics, so as to provide references for the treatment of sediment deposition in the reservoir. By comparing the simulated results of two mathematical models including Realizable k-ε and standard k-ε,the feasibility using value simulation method to forecast reservoir sediment deposition was verified,meanwhile, both accuracies were compared.

valuable simulation; reservoir; silt sedimentation; Realizable k-ε

1007-7596(2016)06-0122-04

2016-04-16

陳偉(1975-)，男，新疆烏蘇人，工程師，從事水務工作。

TV145；TV697.22

B