高含沙水流對亭口水庫泄洪洞泄流能力影響的試驗研究與分析探討

李 樂，張耀哲

(1.西北農林科技大學，陜西 楊凌 712100；2.咸陽市水利工作隊，陜西 咸陽 712000)

1 引言

多沙河流的泄洪建筑物經常遭遇高含沙水流，與一般的清水或低含沙量洪水相比，這一特殊工況對過流能力的影響如何，經常引起工種技術人員的思考，但很少見到相關的分析探討。它所涉及到的主要問題是高含沙水流阻力規律的研究。近年來學術界的研究已基本達成共識，即水庫河道的高含沙水流為賓漢非牛頓體。本文以亭口水庫工程為例分析和探討了高含沙水流對水庫泄洪建筑物過流能力的影響，得出了實際工程運行中高含沙水流工況對排沙洞泄流能力的影響可忽略不計的初步結論，為確定亭口水庫工程排沙洞布置方案提供了依據。

2 關于賓漢體高含沙水流阻力規律的研究

2.1 賓漢非牛頓流體

流體的流型問題較為復雜，總體而言當流體流型不能再用牛頓流體的本構方程描述時，稱之為非牛頓流體。為便于研究，一般將挾沙水流顯示為非牛頓流體運動特征統稱為高含沙水流。學者通過大量的試驗得出結論，挾沙水流中粘性細顆粒含量是影響相同濃度時挾沙水流是否為非牛頓體的重要因素。非牛頓流體一般分為賓漢體、偽塑性體、膨脹體。賓漢體模型的本構方程如下式：

(1)

賓漢體流變特性的主要規律有：①隨著濃度的增加、細顆粒含量增多，剛度系數增大。由于“剛度系數”是反映高含沙水流內在阻力的尺度，因此濃度和級配組成是影響高含沙水流阻力的一個主要因素。因此泥沙濃度、顆粒級配是影響高含沙水流阻力的主要因素。②與“剛度系數”一樣，極限切應力值的變化也與泥沙濃度、細顆粒泥沙含量有關，其值也隨含沙量和細顆粒泥沙含量的增加而增加。

2.2 賓漢體高含沙水流阻力規律問題研究

多年來，對于高含沙水流的阻力問題主要是通過室內水槽及管道試驗進行研究。但目前由于高含沙水流的模型率問題沒有得到解決，將室內小尺度水槽和管道實驗的成果仍無法有效通過比尺放大推演到原型。盡管高含沙水流阻力問題非常復雜研究成果甚微，但一些定性的研究結果對于我們討論相關工程問題仍具有指導意義。

首先，已有的共性結論是特定條件下高含沙水流具有一定的減阻作用，其中0.015 mm粒徑以下粘性顆粒的含量，對減阻的效果有直接影響，級配適宜的高含沙水流可以在天然河道中和灌溉渠道中長距離輸送說明了這一點。

其次，用系列模型試驗法(不同比尺的模型由小到大研究同一中物理現象，以觀測尺度對特定物理規律的影響)對高含沙水流的阻力變化規律進行定性研究的結論是高含沙水流的阻力比尺是小于1的，也就是說受極限切應力τb影響，在室內的小尺度模型中，水流的阻力等參數比原型反映得更大。

3 亭口水庫壩區高含沙渾水的流變學特征

3.1 亭口水庫排沙洞基本情況

亭口水庫工程地處陜西省咸陽市長武縣境內涇河一級支流黑河上，壩址距離河口2.0 km，為Ⅱ等大(2)型水庫，是彬長礦區的骨干水源工程，以工業和城鎮生活供水為主，兼顧發電、防洪等綜合效益。壩址處河床段為攔河壩，在右岸依次布置泄洪排沙洞、輸水洞、溢洪道。樞紐右岸緊靠大壩布置泄洪排沙洞，主要任務是在運行期起到泄洪、排沙的作用，在施工期兼顧導流的任務，洞室圍巖以砂巖、頁巖為主。泄洪排沙洞主要由進口放水塔、洞身段、出口明渠段、挑流鼻坎段、護坦段以及下游先鋒渠等部分組成[1]。根據亭口水庫工程調水調沙的運行方式，排沙洞常遇的高含沙水流多由異重流或渾水水庫排沙形成，具有泥沙粒徑細、粘性顆粒含量高等特點，水流進入壩區后，經過排沙洞下泄的高含沙水流中泥沙顆粒的上限粒徑約為0.018 mm，在常遇工況下排沙洞含沙量下限大于等于80 kg/m3。

3.2 亭口水庫壩區高含沙水流特性

黑河亭口水庫所屬流域是黃河中游細泥沙源區，其所在的區域室內試驗研究及野外原型觀測分析的結果表明，當水流泥沙粒徑小于0.01 mm時的含量同懸移質泥沙的含量成正向關系，細顆粒含量同懸移質含沙量的數量關系可以表示為：

S<0.01=135.31lgS-270.6

(2)

式中：S<0.01為水流中泥沙粒徑小于0.01 mm的含沙量；S為含沙水體中懸移質的含沙量，kg/m3。

按此共性的規律特征配置沙樣并進行不同濃度的流變試驗，結果表明粒徑小于0.01 mm的細顆粒泥沙含量大于50 kg/m3時，即可用賓漢體模型的本構方程描述其流型。

3.3 亭口水庫所屬流域高含沙水流的流變參數測定

根據亭口水庫所屬地域高含沙水流的沙樣特征，共配置了三種沙樣(按照不同粘性顆粒含量)，在含沙濃度不同的情況下進行流變試驗，可以估算出亭口水庫泄洪排沙洞在常遇高含沙水流的情況下流變參數的變化范圍。如表1所示，極限切應力τb項可由3組沙樣觀測到，表明了在所測試的含沙濃度范圍內的水流流型為典型的非牛頓賓漢流體，其流變特性符合賓漢體流變特性的共性規律，現有的高含沙非牛頓賓漢流體理論可以對該問題進行分析討論。

表1 亭口水庫排沙洞常遇高含沙水流流變參數

4 清水與高含沙量渾水對排沙洞泄流能力的影響分析

4.1 對比試驗觀測方法的設計

由于當前的學術水平無法將一般的水工模型試驗所得到的數據推演到原型直接應用，故本文依據小尺度模型試驗所觀測的結果和高含沙水流系列模型的試驗研究成果，對高含沙水流影響排沙洞過流能力的問題作初步分析和評估。試驗方案主要考慮了相同流量不同含沙量之間的對比，并通過觀測上游渠道的量水、在前池加沙、排沙洞出口測試濃度這種過程控制的方法進行試驗，通過對上游相同來流量的觀測，并根據泄洪排沙洞內在不同含沙量下水庫水位的壅高程度，得出過流能力受高濃度含沙量的影響程度。試驗采用的排沙洞體型及進出口邊界條件與《亭口水庫工程排沙洞水工模型試驗報告》推薦方案相同。模型尺度為原型的1/50。

4.2 試驗工況及觀測結果

試驗按下泄流量大小分為六組，每組又分為三種含沙量，觀測水庫水位壅高程度。每組設清水試驗做參照，估算由于水位壅高值導致的流量變化，繼而給出相同模型尺度下的含沙量對過流能力的相對影響程度，在此基礎上再根據系列模型試驗的初步研究結論，評估原型工程應用時高含沙水流對過流能力的影響。觀測結見表2。

表2 相同來流量及運行條件不同含沙量庫水位壅高統計

對比試驗結果表明，與清水試驗相比，在高含沙量水流情況下，壩前的水庫水位均有壅高；在不同工況下，模型內由測針讀出的水位壅高值均分布在0 cm～3 mm范圍內，由此推演至原型水位雍高值在0 cm～12 cm之間，但是在相同流量不同含沙量和不同流量間之間沒有反映出明顯的規律[2]。此次試驗采用的是上游量水，但在本次的渾水試驗中前池庫底有薄層泥沙淤積，它會所引起整體體積增大，導致水位升高，考慮到這個原因，水庫水位的實際壅高值小于試驗結果。根據表2所得到的結果，參照推薦的弧門全開時、用清水試驗方案得出的水位流量關系曲線，上述水庫水位的壅高值等同于過流能力減少了0～3%。但是，考慮到試驗觀測手段的誤差影響，尤其是考慮到高含沙水流阻力參數模型比尺小于1的共性規律，傾向于認為高含沙量對排沙洞泄流能力的影響可以忽略不計。

5 結論

1)亭口水庫所在地域的高含沙水流在所測試的含沙濃度范圍內的水流流型為典型的非牛頓賓漢流體，可用現有的高含沙非牛頓賓漢流體的相關理論進行本工程問題的分析討論。

2)根據小尺度模型試驗所得到的結果，結合高含沙水流系列模型試驗的相關研究成果，對高含沙水流影響排沙洞過流能力的問題作初步分析和評估?？紤]到試驗觀測手段的誤差影響，尤其是考慮到高含沙水流阻力參數模型比尺小于1的共性規律，傾向于認為高含沙量水流對排沙洞泄流能力的影響可以忽略不計，所得此結論可以在工程實際的運行管理中參考使用。