呂 敏

（遼寧省丹東水文局，遼寧丹東118000）

?

水庫(kù)泥沙淤積分析及庫(kù)容測(cè)量

呂 敏

（遼寧省丹東水文局，遼寧丹東118000）

摘要：基于國(guó)內(nèi)外研究進(jìn)展，結(jié)合流體力學(xué)、水動(dòng)力學(xué)、泥沙淤積理論，對(duì)水庫(kù)泥沙淤積情況進(jìn)行分析。并研究了泥沙淤積厚度與回水抬升之間的關(guān)系，采用平均斷面法對(duì)水庫(kù)庫(kù)容進(jìn)行計(jì)算，希望為今后水庫(kù)淤積研究提供理論參考。

關(guān)鍵詞：水庫(kù)；淤積；庫(kù)容；回水抬升

水庫(kù)是我國(guó)抗洪抗旱過(guò)程中的重要水利設(shè)施，現(xiàn)有供水量的35%都由水庫(kù)提供［1］。庫(kù)容-水位關(guān)系曲線是水庫(kù)安全調(diào)度重要依據(jù)，庫(kù)容計(jì)算的準(zhǔn)確度直接關(guān)系到水庫(kù)和發(fā)電機(jī)組的安全運(yùn)行。研究表明，各類水庫(kù)病害中，水庫(kù)泥沙淤積最為嚴(yán)重，特別是水庫(kù)尾部的淤積，其導(dǎo)致庫(kù)容減小，降低了水庫(kù)的可利用率［2～3］。大量的泥沙淤積不僅會(huì)引起回水抬高，還會(huì)導(dǎo)致內(nèi)澇、水位增加等災(zāi)害，因此研究水庫(kù)泥沙淤積問(wèn)題具有重要意義。

雖然在水庫(kù)竣工時(shí)對(duì)水庫(kù)容量進(jìn)行了測(cè)量，但隨著水庫(kù)底部泥沙淤積，水庫(kù)底部的地形不斷發(fā)生變化，初始測(cè)量數(shù)據(jù)已經(jīng)失效，需要對(duì)水庫(kù)庫(kù)容進(jìn)行重新測(cè)量，以保證水庫(kù)的安全運(yùn)行。

二維泥沙輸運(yùn)方程表明，大粒徑砂礫在水庫(kù)尾部沉積，二維模型可以準(zhǔn)確的模擬泥沙沿水深、河寬的運(yùn)動(dòng)規(guī)律［4］。基于國(guó)內(nèi)外研究進(jìn)展，結(jié)合流體力學(xué)、水動(dòng)力學(xué)、泥沙淤積理論，對(duì)水庫(kù)泥沙淤積情況進(jìn)行分析。并研究了泥沙淤積厚度與回水抬升之間的關(guān)系，采用平均斷面法對(duì)水庫(kù)庫(kù)容進(jìn)行計(jì)算，計(jì)算結(jié)果滿足規(guī)范要求。

1　泥沙淤積數(shù)學(xué)模型

泥沙運(yùn)動(dòng)求解主要利用水沙運(yùn)動(dòng)規(guī)律，在此基礎(chǔ)上應(yīng)用泥沙輸運(yùn)方程，對(duì)泥沙淤積與回水問(wèn)題進(jìn)行理論分析。依據(jù)SaintVenant定理，水庫(kù)中的泥沙運(yùn)動(dòng)應(yīng)該采用一維泥沙輸運(yùn)方程進(jìn)行求解［5］。泥沙流動(dòng)首先滿足連續(xù)性方程：

式中：A—水庫(kù)進(jìn)口河流截面積，m2；Q—入庫(kù)流量，m3/s；x—流動(dòng)距離，m；t—流動(dòng)時(shí)間，s。

其次，泥沙運(yùn)動(dòng)還滿足動(dòng)量方程：

其中：Z—水庫(kù)水位，m；g—重力加速度，m3/s2；n—河流的粗糙系數(shù)；R—河流的水力半徑，m。

水庫(kù)底部變形方程：

其中：γ0—容重，kg/m3；Gsb—水庫(kù)進(jìn)口輸沙率，kg/s；αsk—恢復(fù)系數(shù)。

計(jì)算水庫(kù)淤積量前，首先要確定各參量的取值，淤積參數(shù)主要包括：粗糙率、沉降速度、挾沙系數(shù)、容重等。泥沙淤積速率：

取γs=2646kg/m3；γ=1335kg/m3；Q=0.02 m3/s；尾水水深H=0.44m；含沙量S =1kg/m3，運(yùn)用有限差分法對(duì)水庫(kù)淤積問(wèn)題進(jìn)行求解［6-7］。7h后的水庫(kù)水位情況見圖1。

由圖1可知：泥沙經(jīng)過(guò)7h的沉降后，在水平方向15～45m處已經(jīng)出現(xiàn)了明顯的泥沙淤積，30m出的淤積量最大。通過(guò)對(duì)比計(jì)算值和試驗(yàn)值，發(fā)現(xiàn)該模型計(jì)算結(jié)果基本精確，只是在沉積最大位置存在一定偏差。圖2給出了7小時(shí)后的泥沙累積淤積量。

圖1　7h后的水庫(kù)水位情況

圖2　7h后的泥沙累積淤積量

由圖2可知：沿x軸方向，累積淤積量先快速上升，隨后變成緩慢上升，拐點(diǎn)位于x=35m處。當(dāng)x＞40m后，泥沙淤積量基本飽和。從試驗(yàn)值和計(jì)算值的對(duì)比情況看，兩組數(shù)據(jù)間基本吻合，可見本文建立的水庫(kù)淤積數(shù)學(xué)模型精確度較高。只是在x＞40m后兩組數(shù)據(jù)有一定差異，測(cè)量值明顯小于理論計(jì)算值，說(shuō)明該計(jì)算模型只能用在淤積生長(zhǎng)區(qū)。

2　庫(kù)容淤積測(cè)量

由泥沙淤積數(shù)值計(jì)算可以得出，影響水庫(kù)淤積的因素多種多樣，因此需要對(duì)庫(kù)容進(jìn)行測(cè)量，以驗(yàn)證計(jì)算值的準(zhǔn)確度［6～7］。從研究情況看，泥沙淤積與回水末端隨時(shí)間的變化曲線十分類似，但是淤積起始點(diǎn)更靠近水庫(kù)大壩。研究表明，水庫(kù)泥沙淤積量取決于Jc/J0和h/H。圖3水庫(kù)三維地形分布情況。

文中的水庫(kù)庫(kù)容測(cè)量采用平均斷面法，其計(jì)算式為：

圖3　水庫(kù)三維地形分布情況

其中：Vhi，i+1—兩個(gè)相鄰斷面間水位為h時(shí)的庫(kù)容量；Shi—斷面在水位h以下的面積；Di，i+1—兩個(gè)斷面間的平均距離；N—斷面總數(shù)。

利用上述方法對(duì)某水庫(kù)進(jìn)行庫(kù)容測(cè)定，測(cè)量了不同高程下的總庫(kù)容變化情況，并將測(cè)量數(shù)據(jù)與2004年的數(shù)據(jù)進(jìn)行對(duì)比，計(jì)算出庫(kù)容減少量和庫(kù)容損失，具體計(jì)算結(jié)果見表1。

表1　水庫(kù)庫(kù)容變化情況

該水庫(kù)的設(shè)計(jì)蓄水位為330m，原設(shè)計(jì)庫(kù)容為25.8億m3，2004年測(cè)量得到的庫(kù)容為23.12 億m3，2014年測(cè)量得到庫(kù)容為21.84億m3。水庫(kù)運(yùn)行10年間，庫(kù)容一共減小了1.28億m3，損失率高達(dá)5.55%。當(dāng)水位為330m時(shí)，水庫(kù)水深約為35m，水面寬度以300m計(jì)算，則由平面定位誤差導(dǎo)致的庫(kù)容量計(jì)算誤差在0.35%左右。

該次測(cè)量的總庫(kù)容包括測(cè)試區(qū)域和估算區(qū)域，對(duì)其進(jìn)行準(zhǔn)確度計(jì)算，發(fā)現(xiàn)斷面距離的計(jì)算誤差為2.8%，庫(kù)容量的計(jì)算誤差為9.8%，由于等高距為25m，因此許多復(fù)雜的地形區(qū)域在圖3中不能清楚的表示，目前仍不能評(píng)估由此帶來(lái)的計(jì)算誤差。

3　結(jié)語(yǔ)

雖然在水庫(kù)竣工時(shí)對(duì)水庫(kù)容量進(jìn)行了測(cè)量，但隨著水庫(kù)底部泥沙淤積，水庫(kù)底部的地形不斷發(fā)生變化，初始測(cè)量數(shù)據(jù)已經(jīng)失效，需要對(duì)水庫(kù)庫(kù)容進(jìn)行重新測(cè)量，以保證水庫(kù)的安全運(yùn)行?；趪?guó)內(nèi)外研究進(jìn)展，結(jié)合流體力學(xué)、水動(dòng)力學(xué)、泥沙淤積理論，對(duì)水庫(kù)泥沙淤積情況進(jìn)行分析。并研究了泥沙淤積厚度與回水抬升之間的關(guān)系，采用平均斷面法對(duì)水庫(kù)庫(kù)容進(jìn)行計(jì)算，計(jì)算結(jié)果滿足規(guī)范要求。

（1）泥沙經(jīng)過(guò)7h的沉降后，在水平方向15～ 45m處已經(jīng)出現(xiàn)了明顯的泥沙淤積，30m出的淤積量最大，計(jì)算模型在沉積最大位置存在一定偏差。

（2）沿x軸方向，累積淤積量先快速上升，隨后變成緩慢上升，拐點(diǎn)位于x=35m處。試驗(yàn)值和計(jì)算值基本吻合，是在x＞40m后測(cè)量值明顯小于理論計(jì)算值，說(shuō)明該計(jì)算模型只能用在淤積生長(zhǎng)區(qū)。

（3）水庫(kù)運(yùn)行10年間，庫(kù)容一共減小了1.28 億m3，損失率高達(dá)5.55%。許多復(fù)雜的地形區(qū)域不能清楚的表示，由此帶來(lái)的計(jì)算誤差還有待進(jìn)一步研究。

參考文獻(xiàn)

［1］聶秋月，聶慶林，張謙厚，等.魯北地區(qū)河流泥沙淤積分析及對(duì)策探討［J］.泥沙研究，2011（01）：76-80.

［2］孫顯忠，孔令來(lái)，張松濤.金溝河流域渠首樞紐泥沙淤積原因分析及治理［J］.水利技術(shù)監(jiān)督，2011（01）：64-66.

［3］曹小武.深圳市東江水源工程取水口泥沙淤積分析及治理措施［J］.中國(guó)農(nóng)村水利水電，2011（09）：15-17.

［4］肖敬光，劉心玲.吉林省中部城市引松供水工程豐滿取水口泥沙分析［J］.水利規(guī)劃與設(shè)計(jì)，2014（04）：22-25.

［5］丁曼，于德萬(wàn).團(tuán)山子水庫(kù)泥沙淤積分析計(jì)算及防治措施［J］.水利規(guī)劃與設(shè)計(jì)，2014（04）：34-37.

［6］于立，王鑫，南向東.3S測(cè)量技術(shù)在鬧德海水庫(kù)泥沙測(cè)量中的應(yīng)用［J］.黑龍江水利科技，2014（03）：103-104.

［7］何文權(quán).金溝河渠道泥沙磨損修補(bǔ)效益分析［J］.水利技術(shù)監(jiān)督，2013（03）：29-31.

中圖分類號(hào)：TV697

文獻(xiàn)標(biāo)識(shí)碼：B

文章編號(hào)：1008-1305（2016）01-0084-03

DOI：10.3969 /j.issn.1008-1305.2016.01.029

收稿日期：2015-06-05

作者簡(jiǎn)介：呂 敏（1963年—），女，工程師。