老木孔庫區(qū)河段動床物理模型試驗研究

倪志輝,易 靜,張緒進,杜宗偉 ,王云莉

(1.重慶交通大學 水利水運工程教育部重點實驗室國家內河航道整治工程技術研究中心，重慶 400074；2. 重慶交通大學 西南水運工程科學研究所，重慶 400016)

0 引 言

由于我國有許多河流是含沙量高、輸沙量大的多泥沙河流，庫區(qū)的泥沙淤積問題異常嚴重且相當普遍，河工動床模型試驗成為解決江河治理及大型水利工程泥沙問題的主要手段，并大有日益受到重視的勢頭[1-2]。國內在河工動床模型試驗方面的研究起步較晚，但近年來眾多學者對此展開了大量研究并獲得了寶貴經驗。唐峰，等[3]采用實體動床模型試驗研究了三峽工程運用初期不同時期荊江重點險河段的沖淤變化過程；余明輝，等[4]通過建立動床實體模型，從三維空間角度詳細研究了取、排水河段水流泥沙運動及河道沖淤演變趨勢、取水口附近底沙運動情況等問題；陳孟榮，等[5]采用動床物理模型試驗，對黃金埠電廠補給水取水口局部區(qū)域的水流、泥沙運動及地形沖淤規(guī)律進行了研究；張勝利，等[6]通過動床物理模型試驗，探討了引灤入津工程中黎河輸水運用后河道演變規(guī)律、河道抵御洪水能力以及河道治理措施等內容。此外，還有其他學者通過動床物理模型試驗研究河床沖淤規(guī)律、泥沙淤積[7-10]。

老木孔航電樞紐工程是岷江樂山—宜賓段航電梯級開發(fā)方案中的第一級。樞紐成庫后，庫區(qū)河床邊界條件顯著改變，另外岷江、大渡河若干大型水庫相繼建成運行，進入老木孔庫區(qū)的水沙條件也相應變化。由于大壩的壅水作用，庫區(qū)流速減緩，上游來沙在庫區(qū)淤積是必然的，其淤積速度、淤積分布、淤積量與庫區(qū)航道水流條件、河道行洪能力、大橋防沖及建港環(huán)境等關系極為密切，對工程的成敗有重大影響，因此需要通過動床試驗研究庫區(qū)泥沙淤積及其影響問題。

1 河段的水沙特性

岷江是長江上游的一級支流，干流全長735 km，流域面積135 811 km2。岷江流域為降水補給河流，洪水主要由暴雨形成，暴雨最早出現(xiàn)在4 月，主要集中在6—9月，其中特大暴雨主要出現(xiàn)在6月下旬—8月中旬，洪水陡漲陡落。據五通橋站實測洪水資料統(tǒng)計，多年平均流量2 510 m3/s。年最大洪峰流量發(fā)生在6—9月，其中7—8月份最多，占全年的90.6%，6、9月出現(xiàn)的頻次大致相當。年最大洪峰最早出現(xiàn)在6月27日(1987年，洪峰流量16 900 m3/s)，最晚出現(xiàn)在9月5日(1982年，洪峰流量22 300 m3/s)，歷年實測年最大洪峰流量的最大值為35 300 m3/s(1955年7月14日)，最小值為9 640 m3/s(1969年8月11日)。

根據五通橋水文站實測及插補延長后的1953—2005年懸移質泥沙資料統(tǒng)計，并考慮大渡河下游的龔嘴及銅街子水電站的水庫攔沙作用后，五通橋水文站多年平均懸移質年輸沙量為5 240萬t。根據干支流的入庫懸移質年輸沙量統(tǒng)計，老木孔航電樞紐工程干支流的懸移質年輸沙量之和為5 722萬t。

老木孔航電工程的入庫推移質泥沙將主要來自大渡河的支流青衣江流域。岷江上游干、支流無推移質輸沙率實測資料，老木孔航電工程推移質輸沙量采用推移質輸沙率公式計算。采用五通橋水文站河段水力要素，用分組輸沙率公式計算流量與推移質輸沙率關系，推算岷江下游老木孔河段多年平均推移質年輸沙量為27萬t。

老木孔航電工程壩址位于樂山市五通橋區(qū)，岷江與大渡河匯口下游約12.3 km處的岷江干流河段。樂山至宜賓段，天然落差97 m，河道平均比降0.59‰，水面寬400 ～1 000 m左右。沿河兩岸分布有較寬的臺地及漫灘，河谷呈箱形，兩岸為低山起伏的丘陵地形。庫區(qū)河段洲灘發(fā)育、汊浩縱橫、河床寬淺、水流分散、流態(tài)紊亂，工程所在河段水沙運動十分復雜。本模型試驗河段范圍是岷江干流府河大橋—老木孔河段(長約14 km)，包括支流大渡河大件碼頭～匯流河口(長約3 km)，總長約17 km，如圖1。

圖1 模型布置Fig.1 Model layout

2 模型設計與相似性驗證

2.1 模型設計

模型設計除確保水流條件相似外，還須確保泥沙運動、河床沖淤的相似。

合適的模型沙是確保模型沖淤相似的重要條件。模型采用平面比尺λL=200和垂直比尺λh=100，并采用比重γs=1.33 t/m3的榮昌精煤作為模型沙，各相似比尺要求均能得到較好的滿足，詳見表1。

表1模型比尺匯總
Table.1Summaryofsimilarityscalesformodel

(續(xù)表1)

2.2 模型驗證

2.2.1 水面線驗證

河床糙率相似是實現(xiàn)河道水面線相似的重要條件之一，為滿足模型相似的要求，保證模型與原型阻力相似，須進行水面線驗證試驗。為此，采用實測的三級流量(Q=14 030，8 500，574 m3/s)的瞬時水面線資料進行驗證。經實測資料計算分析，原型河床糙率為0.026～0.043，相對應的模型河床糙率為0.017～0.028。同時模型采用小卵石梅花型加糙方式來滿足河床阻力相似要求。試驗結果表明：模型水位與原型水位相近，誤差多在±0.1 m以內，達到了水面線相似的要求。由于瞬時水面線實際觀測時間不一致，故此處只選取部分河段洪、中、枯三級流量同一時刻瞬時水面線驗證情況，如圖2。

圖2 水面線驗證Fig.2 Verification of water surface profile

2.2.2 流速流向驗證

根據四川省交通廳交通勘察設計研究院2010-8-22T14：00在庫區(qū)河段實測岷江流量8 500 m3/s(其中，岷江流量為Q=5 000 m3/s，大渡河流量為3 500 m3/s)時的流速流向資料，在模型上同樣位置投放浮標，并觀察浮標的運動軌跡，通過與原型比較，模型與原型浮標的運動軌跡相近，運動方向吻合，表明模型表水流流向與原型基本一致。

2.2.3 河床沖淤驗證

在前述水流條件相似的前提下，還須通過動床驗證試驗，檢驗模型河床的沖淤數(shù)量、沖淤部位及沖淤過程的相似性。根據大渡河、岷江的水沙情況，并結合沙、卵石推移質運動特點對驗證時段的水沙過程分別進行概化，其中大渡河劃分為26個時段；岷江劃分為22個時段。

模型試驗模擬了老江壩左汊驗證河段2008-8—2009-12兩次實測河道地形沖淤變化情況，試驗結果表明：模型與原型沖刷量和淤積量的相對誤差分別為14.8%和9.6%，模型沖淤部位、斷面形態(tài)和高程與原型實測基本相似。因此符合相關規(guī)程的要求，驗證河段沖淤量對照情況如表2。

表2 驗證河段沖淤量對照

綜上所述，模型的水面線、流速流向、流速分布及河床沖淤變化等均與原型基本相似，本模型較好的復演了老木孔河段的水沙運動與河床演變情況，可以用來預測老木孔樞紐庫區(qū)泥沙淤積情況。

3 試驗方案

筆者采用1 500 m堤距保留老江壩方案進行泥沙淤積動床試驗，該方案兩岸防洪堤均起于庫尾附近，至下游與老木孔樞紐閘壩岸坡段以及船閘上引航道右邊墻銜接。堤線總長17.73 km，兩岸防洪堤間距為1 500 m，保留老江壩。為不影響庫區(qū)行洪能力，在老江壩河段適當加大河底開挖深度，其堤線布置、線形及河底疏浚高程具體布置見圖3。根據國家防洪規(guī)范規(guī)定，老木孔樞紐庫區(qū)兩岸防護堤，防洪標準按20 a一遇洪水設計，設計流量Q=35 600 m3/s。

圖3 方案布置Fig.3 Layout of program

4 試驗成果及分析

4.1 庫區(qū)泥沙沖淤量分析

老木孔樞紐建成運行后，庫區(qū)河段水位升高，流速減小，比降趨緩，水流輸沙能力降低，泥沙發(fā)生累積性淤積。庫區(qū)泥沙淤積量變化見表3。老木孔樞紐庫區(qū)左右兩岸防護堤防洪標準按20 a一遇洪水設計，因此庫區(qū)泥沙淤積試驗只分析20 a內的淤積情況。

表3庫區(qū)泥沙淤積量變化

Table3Siltationamountchangeatthereservoirarea

由表3可見，水庫運用10 a和20 a全庫區(qū)泥沙累計淤積量分別為750×104和1 233×104m3。后15 a平均年淤沙量為48.2×104m3，其年均淤積量僅為前5 a平均淤積量的47.2%。試驗過程中，全程監(jiān)測了出庫沙量的變化情況，至水庫運用20 a時，出庫沙量為入庫沙量的78%左右。表明庫區(qū)河床沖淤尚未處于平衡狀態(tài)。老木孔水庫運用20 a，庫區(qū)河段累計淤沙量1 233×104m3，泥沙淤積主要分布在老江壩以下河段(常年回水區(qū))，而老江壩以上河段(變動回水區(qū))淤積不多。

相對而言老木孔庫區(qū)泥沙淤積量不大，分析其原因在于：①岷江上游干、支流建庫攔截了大量泥沙，使老木孔庫區(qū)的入庫懸沙總量大量減少，且粒徑變細，推移質泥沙主要來源于青衣江，數(shù)量較少；②左右兩岸防洪堤縮窄了庫區(qū)河道的寬度，使庫區(qū)水流仍能保持一定的流速，從而具有一定的輸沙能力，使泥沙不易落淤；③流量大于10 000 m3/s時，水庫處于敞泄沖沙狀態(tài)，庫區(qū)流速、比降較大，部分淤沙被沖走，尤其是主槽沖沙效果顯著，使其基本保持微淤狀態(tài)。

4.2 庫區(qū)泥沙淤積分布分析

成庫后庫區(qū)河道普遍淤積，但淤積分布極不均衡。限于篇幅，列舉了CS50和CS65斷面沖淤變化情況，如圖4。

圖4 庫區(qū)斷面沖淤變化Fig.4 Variance of cross-section erosion at the area of reservoir

淤沙主要集中在邊灘開挖區(qū)、岸邊緩流帶、回流區(qū)、彎道凸岸及河面放寬處等。而河道主槽淤積卻較少。重點淤積部位在老江壩尾至壩址河段左右兩岸邊灘開挖區(qū)和涌斯江口附近(包括疏浚后的江心灘)。右岸邊灘淤積厚度普遍達到2 m以上，對老江壩碼頭一期工程、船閘上引航道及樂山港規(guī)劃岸線的利用帶來了不利影響。涌斯江口附近淤積厚度也達到3 m左右，增大了生態(tài)電站的維護清淤量。其次，彎道凸岸(杜家場)及岷江河口段等(緩流區(qū))處也有一定淤積。

針對老木孔庫區(qū)泥沙淤積分布不均，分析其原因在于：與水庫運行方式、庫區(qū)河床形態(tài)及流場分布等密切相關。兩岸防洪堤修建后，庫區(qū)河道規(guī)順，邊灘降低(因開挖)，中、小流量時河床普遍淤積，但當入庫流量大于10 000 m3/s電站停機敞泄沖沙時，由于灘、槽高差較大(達5 m左右)，水流大部分集中到主河槽內，造成主槽流速加大，邊灘流速減小的態(tài)勢。結果是主槽淤沙大部分被帶走，而邊灘淤沙幾乎未動，即主槽河床呈現(xiàn)沖淤交替的微淤狀態(tài)，而邊灘則成累積性淤積的態(tài)勢。

表4重點部位泥沙淤積情況

Table4Siltationconditioninkeyparts

4.3 泥沙淤積對樂山大佛、樂山城的影響

根據國家防洪規(guī)范規(guī)定，樂山城的城市防洪標準為50 a一遇，設計洪水流量Q=43 200 m3/s(P=2%)。模型實測了岷江遭遇各級頻率洪水時庫區(qū)左岸的沿程水面線，限于篇幅，僅列舉了2級流量庫區(qū)河段水面線，見圖5。

圖5 建庫前后庫區(qū)河段水面線 Fig.5 Water surface profile of reservoir area before and after the reservoir construction

據物理模型試驗成果分析，水庫運行20 a末，在樂山大佛處，當岷江遭遇20 a一遇洪水時，淤積后的洪水位為363.86 m，較同頻率下的天然水位低0.03 m；當岷江遭遇50 a一遇洪水時淤積后的洪水位為364.84 m，較同頻率下的天然水位低0.01 m；當岷江遭遇100 a一遇洪水時淤積后的洪水位為365.76 m，與同頻率下的天然水位齊平。在樂山城(肖公咀)，當岷江遭遇20 a一遇洪水時，淤積后的洪水位為364.66 m，較同頻率下的天然水位低0.02 m；當岷江遭遇50 a一遇洪水時，樂山城(肖公嘴)水位為365.83 m，與該處天然水位一致；當岷江遭遇100 a一遇洪水時淤積后的洪水位與同頻率下的天然水位一致。

由分析可知，老木孔樞紐按擬定的水庫調度方式運行20 a，庫區(qū)泥沙淤積后各級頻率洪水的洪水位仍未超過建庫前同頻率下的天然洪水位，表明庫區(qū)泥沙淤積對樂山大佛及樂山城區(qū)的影響不大。

5 結 論

1)通過建立老木孔樞紐庫區(qū)河段泥沙沖淤模型，利用實測資料分別從水位、流速變化和泥沙淤積分布等方面，驗證了模型的相似性符合要求，可以用來預估庫區(qū)河道范圍內相應變化情況。

2)樞紐建成運行后，庫區(qū)河段水位有所升高，流速有所減少，比降將逐漸趨于平緩，水流輸沙能力也會隨之降低，泥沙沿河道發(fā)生累積性淤積，但淤積分布不均勻。泥沙淤積主要集中在岸邊緩流帶、河道邊灘開挖區(qū)、回流區(qū)、彎道凸岸、河面放寬處等，而主槽河床呈現(xiàn)沖淤交替的微淤狀態(tài)，而邊灘則成累積性淤積的態(tài)勢。

3)老木孔樞紐按擬定的水庫調度方式運行20 a，庫區(qū)泥沙淤積后各級頻率洪水的洪水位仍未超過建庫前同頻率下的天然洪水位，表明庫區(qū)泥沙淤積對樂山大佛及樂山城區(qū)的影響不大。

[1] 謝鑒衡.河流泥沙工程學:下冊[M].北京:水利出版社,1981:212-223.

Xie Jianheng.Sediment Engineering:Volume Ⅱ[M].Beijing:Hydraulic Press,1981:212-223.

[2] 張紅武,馮順新.河工動床模型存在問題及其解決途徑[J].水科學進展,2001,12(3):418-423.

Zhang Hongwu,Feng Shunxin.Problems and their solutions in physical modeling of river process in alluvial streams [J].Advances in Water Science,2001,12(3):418-423.

[3] 唐峰,李發(fā)政,姚仕明,等.三峽工程運用初期荊江河段調關彎道沖淤演變試驗研究[J].長江科學院院報,2011,28(8):77-86.

Tang Feng,Li Fazheng,Yao Shiming,et al.Physical modeling of morphological changes in Tiaoguan bend reach of Jingjiang river at the initial operation stage of Three Gorges Project [J].Journal of Yangtze River Scientific Research Institute,2011,28(8):77-86.

[4] 余明輝,吳衛(wèi)民,劉曉芳.襄樊電廠二期工程取排水河段動床物理模型試驗研究[J].泥沙研究,2007,51(1):76-81.

Yu Minghui,Wu Weimin,Liu Xiaofang.Experimental investigation on the channel bed variation in the water intake reach of Xiangfan power plant [J].Journal of Sediment Research,2007,51(1):76-81.

[5] 陳孟榮,彭楊.黃金埠電廠一期工程取水河段泥沙模型試驗[J].武漢大學學報:工學版,2010,43(3):330-339.

Chen Mengrong,Peng Yang.Sediment model tests on water intake reach of Huangjinbu power plant [J].Engineering Journal of Wuhan University,2010,43(3):330-339.

[6] 張勝利,王喜英,張炳錄,等. 黎河整治動床物理模型試驗研究[J].泥沙研究,1995,39(3):93-100.

Zhang Shengli,Wang Xiying,Zhang Binglu,et al.A movable bed model of Lihe River regulation [J].Journal of Sediment Research,1995,39(3):93-100.

[7] Chien N,Wan Z.Mechanics of Sediment Transport [M].Reston:American Society of Civil Engineers,1999.

[8] Seminara G.Fluvial sedimentary patterns [J].Annual Review of Fluid Mechanics,2010,42(1):43-66.

[9] Craig B B,Steven R A,Drew C B.Criteria for the formation of sediment plugs in alluvial rivers [J].Journal of Hydraulic Engineering,2011,137(5):569-576.

[10] 郭華,張緒進,杜宗偉,等.老木孔水庫成庫后樂山港港區(qū)河段泥沙沖淤演變研究[J].重慶交通大學學報:自然科學版,2012,31(5):1031-1035.

Guo Hua,Zhang Xujin,Du Zongwei,et al.Sediment scour and silting evolution research at the area of Leshan Portin Old Wood Hole Reservoir [J].Journal of Chongqing Jiaotong University:Natural Science,2012,31(5):1031-1035.