鄭國誕 ，陳韜霄 ，馬 駿 ，唐子文 ，陳 剛

（1.浙江省水利河口研究院，浙江 杭州 310020；2.浙江省海洋規(guī)劃設(shè)計研究院，浙江 杭州 310020；3.浙江省河口海岸重點實驗室，浙江 杭州 310020；4.杭州市港航管理局，浙江 杭州 310014）

杭州市桐廬縣現(xiàn)狀水廠取水口位于富春江電站下游約5.4 km，現(xiàn)存在取水規(guī)模不能滿足規(guī)劃供水需求及取水保證率不足的問題，擬將桐廬水廠取水口上移0.9 km，至富春江電站下游約4.5 km，具體位置見圖1。取水口的設(shè)置將引起河床局部沖刷，由于取水口結(jié)構(gòu)復(fù)雜，且取水的時候，還存在局部流態(tài)改變，無法用公式直接計算局部沖刷值，需利用水槽模型試驗研究洪水作用下取水口附近局部沖刷深度。

圖1 取水口位置示意圖

1 概 況

1.1 河段概況

工程處于錢塘江近口段，以徑流為主，近口段的范圍為富春江大壩至聞家堰之間，總長77.0 km，1960年新安江水庫建成后，具有明顯的削減洪峰作用，使流量在年內(nèi)的區(qū)域平均分配，經(jīng)統(tǒng)計，梅汛期洪峰流量可削減20%左右，枯水期流量可增加22%左右。建庫后出現(xiàn)最多的洪水為2017年梅汛期洪水，富春江電站最大下泄流量14 900 m3/s，為1968年12月水庫蓄水以來出現(xiàn)的最大下泄流量，基本達到10 a一遇，下游之江站實測最大流量17 400 m3/s。

1.2 工程概況

取水口設(shè)計正常取水量15.75萬m3/d（1.82 m3/s），雙管進水時，每個進水口設(shè)計進水量0.91 m3/s；單管進水時，每個進水口進水量1.82 m3/s。設(shè)計最大取水量26.25萬m3/d（3.04 m3/s），雙管進水時，每個進水口設(shè)計進水量1.52 m3/s；單管進水時，每個進水口進水量3.04 m3/s。頭部采用樁架固定吸水喇叭口的形式，設(shè)計2只取水喇叭口；取水管道為2根DN1 200鋼管，直徑為1.20 m，從取水頭部穿過富春江江堤至取水泵房，取水頭部分2格，每格可獨立運行，在頂部和側(cè)面設(shè)進水口，并在進水口設(shè)置格柵，頭部結(jié)構(gòu)示意見圖2，頭部模型示意見圖3。

圖2 取水口頭部結(jié)構(gòu)示意圖

2 模型設(shè)計

由于取水口結(jié)構(gòu)相對復(fù)雜，且在工作時還存在取水工況，局部流態(tài)復(fù)雜，用一般的經(jīng)驗公式無法計算局部沖刷情況，也不適合利用變態(tài)的大比尺物理模型，宜通過小比尺水槽模型實驗進行最大沖刷深度的研究。

2.1 模型比尺

試驗在長40.00 m，寬3.70 m的水槽中進行。綜合考慮模型沙的選擇、取水口的尺寸以及試驗場地等情況，確定模型幾何比尺為50。為了能較好地模擬原型水流運動，除了滿足重力相似、水流連續(xù)相似還滿足河床泥沙起動相似以及水下休止角相似[2]。

2.2 模型沙選擇

根據(jù)地質(zhì)勘測資料顯示，工程區(qū)域可沖刷河床基本為砂礫卵石，厚8.20 ～ 8.30 m。河床中值粒徑為42.550 0 mm，懸沙中值粒徑在0.018 6 ～ 0.019 0 mm，經(jīng)結(jié)合張瑞瑾起動公式和秦榮昱計算并考慮研究水域的水深情況，在水深為10.00 ～ 20.00 m下，中值粒徑為42.550 0 mm、平均粒徑為11.110 0 mm的泥沙起動流速為1.60 ～ 1.95 m/s。按泥沙運動的起動相似模型率，試驗中模型沙起動流速應(yīng)為0.23 ～ 0.27 m/s，宜使用原型沙進行試驗。經(jīng)過對多種原型沙的綜合比較，選擇其中值粒徑約為0.300 0 mm，起動流速為 0.25 ～ 0.30 m/s（模型水深 20 ～ 40 cm）[3]。

3 試驗條件

3.1 水流條件[5]

工程位置靠近富春江電站，以徑流為主，徑流的年際變化較大且具有連續(xù)豐、枯水文年交替出現(xiàn)的情況，上游建庫以后，具有明顯的削減洪峰作用，使流量在年內(nèi)的分配區(qū)域平均，經(jīng)統(tǒng)計，梅汛期洪峰流量可削減20%左右，枯水期流量可增加22%左右，富春江電站100 a一遇洪峰流量為23 100 m3/s，下游的桐廬水文站100 a一遇水位為15.25 m，根據(jù)并行化的平面二維數(shù)學(xué)模型[4]計算可知，取水口100 a一遇水位為17.40 m，取水口處最大流速為2.75 m/s。

3.2 起沖高程

取水口斷面靠右岸江道為航道，較為穩(wěn)定，沖淤幅度不大。靠左岸江道斷面形態(tài)變化較大，整體呈現(xiàn)逐年下降趨勢。圖4為2000年以來工程斷面形態(tài)圖，從圖4可知，2004 — 2011年取水口位置河床高程逐年加深，2011年高程接近-1.30 m，2011年以后取水口位置高程穩(wěn)定，2017 —2018年略有淤積，淤積幅度在0.20 m左右。取下包絡(luò)線作為試驗的起沖高程，取水口位置高程為-1.30 m。

圖4 取水口歷年斷面形態(tài)圖

3.3 實驗工況

試驗之前先進行水流條件調(diào)試，試驗段水流條件嚴格按照相似比尺轉(zhuǎn)化為模型水流條件進行控制。水槽模型上、下游邊界采用水位控制，首先控制試驗斷面水位至設(shè)計條件，然后調(diào)整上下游水位過程，最終使試驗斷面水深和流速過程與模型所需值相符。實驗工況共5組，分別為：單孔取水3.04 m3/s、單孔取水1.82 m3/s、雙孔取水3.04 m3/s、雙孔取水1.82 m3/s以及不取水。取水口設(shè)計流量的大小，按照相似比尺換算，選用配有流量計的水泵控制確保取水過程的相似性。

4 試驗成果

實驗現(xiàn)象可見，取水口鋼管樁附近有明顯的局部沖刷坑形成，主要由于鋼管樁等阻水建筑物的建設(shè)，來流受阻后部分轉(zhuǎn)向河底，產(chǎn)生下降水流，并在近河床處形成一橫軸反向旋流，此旋流和來自上游底部集中的縱向水流結(jié)合在一起，形成繞樁靠近河底流向下游的馬蹄形旋渦，在馬蹄形旋渦的作用下，樁周圍河床上的泥沙被沖起帶向背流側(cè)，形成沖刷坑，表現(xiàn)為樁基兩側(cè)以及群樁內(nèi)部區(qū)域沖刷后高程較低，上下游區(qū)域高程略高的形態(tài)。當取水口工作時，四周流態(tài)更加復(fù)雜，繞流、環(huán)流和上升流等多種流態(tài)并存[6]，馬蹄形旋渦強度加強，則沖刷坑深度越深，形態(tài)越大。

4.1 沖刷過程分析

從各種實驗工況的沖刷過程看，雖然取水口取水流量有所不同，但是其沖刷坑的形成、發(fā)展至最終的相對平衡過程往往類似。沖刷實驗過程中對樁基附近河床高程的監(jiān)測表明，初期樁基附近河床快速下切形成沖刷坑，但隨著沖刷的發(fā)展，水深加大，流速降低，水流的挾沙能力逐漸下降，沖刷坑內(nèi)泥沙運移變緩，沖刷率得以迅速減小并漸趨穩(wěn)定，逐漸達到相對沖淤平衡，形成相對穩(wěn)定的沖刷深度，穩(wěn)定的時間為2.0 ～ 2.5 h，相當于洪峰持續(xù)作用18.0 h左右（見圖5）。

4.2 沖刷形態(tài)及沖刷深度分析

圖6為不取水情況下沖刷坑輪廓線示意圖。從圖6可知，取水口周邊受馬蹄形旋渦形成的床面沖刷坑形態(tài)類似于馬蹄形，兩側(cè)深槽向水流方向延伸，樁基底部河床較周邊明顯沖蝕，沖刷坑橫向?qū)挾燃s為2樁基距離的3 ～ 4倍，不同取水方案下的沖刷坑形態(tài)較類似，均為如圖6所示的馬蹄形，只是沖刷坑的橫向?qū)挾扰c縱向?qū)挾炔煌?/p>

圖6 不取水情況下沖刷坑輪廓線示意圖

圖7為不同取水方案下，順水流方向沿外側(cè)樁基的沖刷坑縱剖面圖。從圖7可知，縱剖面圖呈現(xiàn)一個明顯的沖刷漏斗形態(tài)，前端沖刷坑坡度為1∶1.6 ～ 1∶1.8，坡角接近此泥沙粒徑的水下休止角，形成沖刷坑最深點的最大沖刷深度工況從大到小的依次為：單孔取水3.04 m3/s、單孔取水1.82 m3/s、雙孔取水3.04 m3/s、雙孔取水1.82 m3/s以及不取水，最大沖刷深度及沖刷坑范圍見表1。從表1可知，由于取水口的取水效應(yīng)，沖刷深度加強，單孔取水3.04 m3/s時，取水口內(nèi)流速最大，沖刷深度也最大，為3.12 m，比不取水時沖刷深了將近1.00 m，且2根輸水管道均裸露，不取水時，裸露上游側(cè)的輸水管道。圖8和圖9分別為不取水以及單孔取水3.04 m3/s情況下的模型實驗圖。

表1 各工況沖刷結(jié)果綜合表

圖7 沿水流方向的沖刷坑縱剖面圖

圖8 不取水情況下試驗圖

圖9 單孔取水3.04 m3/s試驗圖

5 結(jié) 語

新建桐廬縣水廠取水口位于富春江電站下游4.50 km，屬于錢塘江流域近口段，以徑流為主，利用數(shù)學(xué)模型計算取水口位置100 a一遇流速的水深情況，并且通過河勢分析得到取水口位置最低高程作為起沖高程，并在水槽中進行沖刷實驗，得到不同工況下取水口局部最大沖刷深度以及沖刷坑范圍，單孔取水3.04 m3/s時，沖刷深度為所有工況中最大，為3.12 m，比不取水時沖刷深了將近1.00 m，且2根輸水管道均裸露，建議設(shè)計時考慮單孔取水3.04 m3/s為最不利的設(shè)計工況進行設(shè)計。