輸水渡槽水位異常波動原因分析與改善措施研究

屈志剛 李政鵬

摘要：南水北調中線工程自2020年開始正式啟動大流量輸水，在此期間部分輸水渡槽進出口流態紊亂，槽內水位異常波動，過流能力降低。以澧河渡槽為研究對象，通過高精度三維水動力學數值模型模擬復演了澧河渡槽現狀流態異常現象，探究了渡槽水位異常波動的原因，即由于渡槽出口中墩體型寬平導致水流急劇繞流形成“卡門渦街”現象，墩尾水流出現周期性擺動，并引起兩槽流量的周期性變化，反過來又加重了墩尾水流的擺動，兩者相互激勵形成共振，放大了渡槽水位的異常波動。結合工程實際提出了各種可行的流態優化工程措施，并通過數模分析對比，得出了最適合澧河渡槽的方案，即：在渡槽進、出中墩上下游修建三角形導流墩導流，從工程性價比、結構穩定性等方面綜合考慮，推薦采用“進口長10 m+出口長30 m”的導流墩尺寸方案。相關經驗可供類似調水工程渡槽流態改善借鑒。

關 鍵 詞：輸水渡槽; 水位異常波動; 卡門渦街; 數值模擬; 流態優化; 南水北調中線工程

中圖法分類號： TV131.2

文獻標志碼： A

DOI：10.16232/j.cnki.1001-4179.2022.04.030

0 引 言

南水北調中線工程（以下簡稱中線工程）自2014年12月12日建成通水以來，經歷了汛期和冰期輸水的考驗，在設計流量下整體運行狀況一直良好，成為沿線24座城市供水的生命線，通過實施生態補水，成為助力中國生態文明建設的重要力量。

為加大對河南省、河北省、天津市等地37條河流生態補水，按照相關管理部門工作安排，2020年4月29日起，中線工程首次以設計加大流量進行輸水，以此全面檢驗中線工程大流量運行狀態和最大輸水能力。

中線工程加大流量輸水以來，部分輸水渡槽進出口流態紊亂，槽內水位異常波動，過流能力降低，甚至影響到工程運行安全，其中以澧河渡槽最為典型。通過現場調研，各級領導和專家一致認為有必要對其進行詳細研究，找出渡槽流態異常的原因，并有針對性地提出流態優化工程措施，形成一套集流態研究和施工實施的完整解決方案，為南水北調工程其他輸水建筑物流態優化提供參考和借鑒。

目前國內外許多學者圍繞輸水建筑物流態問題開展了大量的研究，提出了許多新模型、新方法，也取得了豐富的研究成果。例如王才歡等[1]以南水北調十二里河渡槽為研究對象，采用水工模型試驗揭露了大型輸水渡槽水流異常波動的原因，但由于項目周期和客觀條件所限，未對流態優化工程措施進行更深入的、定量的研究。為此，本文通過構建高精度二、三維水動力學數值分析模型，以大流量輸水期間實測水情為邊界條件，模擬復演了澧河渡槽異常流態現象，分析澧河渡槽流態異常原因，并重點對各種導流墩體型和尺寸進行了詳細的、定量的研究，提出了合理的導流墩體型和尺寸。

1 工程概況

南水北調中線工程是一項跨地區、跨流域的特大型調水工程，全長1 432 km，采用高程差自流輸水，設計多年平均調水量95億 m3。工程輸水形式以明渠為主，沿線與鐵路、公路、河流采用立體交叉，其中布置輸水渡槽27座。

為進一步發揮工程社會經濟效益，2020年4月29日起，南水北調中線工程正式啟動大流量輸水，陶岔入渠流量從350 m3/s提升至420 m3/s。大流量輸水期間，部分輸水渡槽進出口流態紊亂，槽內水位異常波動，過流能力降低，迫切需要采取有效措施改善流態，消除槽內水位異常波動，提高渡槽過流能力。

澧河渡槽全長860 m，主體建筑物包括退水閘、進口漸變段、進口節制閘、進口過渡段、渡槽段、出口過渡段、出口檢修閘、出口漸變段、出口明渠段等，澧河渡槽工程布置見圖1。

渡槽設計流量320 m3/s，加大流量380 m3/s，進口設計水位134.598 m，出口設計水位134.118 m，水頭0.480 m。

2 現狀流態問題

大流量輸水期間澧河渡槽流態紊亂表現為：

（1） 進口閘室段2孔水位交替異常波動，即波動存在相反的相位差，一孔出現波峰時另一孔剛好為波谷，最大波動幅度在0.6 m左右，波動周期8 s，該波動引起進口漸變段和上游一定范圍內渠道水位波動。

（2） 槽內水位大幅異常波動，通過持續的現場觀測發現，渡槽流量小于300 m3/s時，槽內水位波動幅度在0.1 m以內;當流量達到320 m3/s時，槽內水位波動幅度迅速增大，達到0.2～0.4 m;當流量超過330 m3/s時，槽內水體開始拍打橫梁，此時槽內水位波動幅度達到0.5～0.7 m。

（3） 在不同流量下澧河渡槽出口均存在尾流擺動，呈現典型的”卡門渦街”現象，只是隨著流量的加大，尾流擺動的范圍和頻率有所提高，擺動周期一般在6～10 s。

澧河渡槽不同流量下水情觀測結果見表1。

3 現狀流態異常原因分析

為還原澧河渡槽大流量下的各種流態異常現象，分析渡槽內水面異常波動產生的原因，本文通過構建三維數模，以典型時間（2020年6月11日）現場觀測水情為邊界條件，對澧河渡槽當天的流態現象進行數值模擬復演。

3.1 模型原理

三維數模基于流體力學基本方程組構建，采用兩方程RNG k-ε紊流模型[2-4]，其控制方程為

連續性方程：

ρt+ρuixi=0（1）

動量方程：

ρuit+xiρuiuj=-pxi+xj[（μ+μt）（uixj+ujxi）]（2）

k方程：

（ρk）t+ρkuixi=xjμ+μtσkkxi+Gk-ρε（3）

ε方程：

（ρε）t+ρεuixi=xjμ+μtσεkxj+C1εεkG-C2ερε2k（4）3E8F2EEA-38EE-4E7D-9CEE-5A4D1C0E46D2

其中，G為由平均速度梯度引起的紊動能產生項，即：

G=μt（uixj+ujxi）uixj（5）

ρ和μ是體積分數的函數，其計算公式為

ρ=αwρw+（1-αw）ρa

μ=αwμw+（1-αw）μa（6）

式中：αw為水的體積分數，ρw和ρa分別為水和空氣密度，μw和μa分別為水和空氣的分子黏性系數，其他常見量物理含義參見文獻[3]。

3.2 模型構建

選取澧河渡槽進口上游200 m至出口下游200 m作為建模范圍，總長1 136 m。

以一般部位0.4 m、重點部位加密至0.2 m進行網格剖分，以2020年6月11日當日實測流量350.34 m3/s和下游實測水位135.12 m作為上、下游邊界條件。澧河渡槽三維數模網格剖分見圖2。

3.3 模擬結果

經模擬，完整還原了2020年6月11日澧河渡槽各種流態現象。

（1） 進口節制閘前水位呈周期性波動，最大波幅約0.5 m，波動周期9 s，左右兩槽水面波動存在相反相位差，即一槽出現波峰時另一槽剛好為波谷[1]，節制閘前水位波動見圖3。

（2） 出口左右兩槽水流以中墩為中心，呈周期性擺動，擺動周期約9 s，擺動還伴隨有方向相反的局部渦流向下游傳遞，即呈現流體力學中典型的“卡門渦街”現象[5-6]，渡槽出口“卡門渦街”現象見圖4～5。

（3） 渡槽內出現超常的水面大波動現象，左右兩槽波動相位差異呈現隨機性，并間歇性出現橫梁阻水和局部漫槽現象。

3.4 槽內水面異常波動原因分析

通過三維數模成功復演現狀澧河渡槽大流量輸水流態后，再結合各種流態改善工程措施研究，共同揭露了現狀流態異常產生的根源。

由于渡槽出口中墩墩尾體型寬平，出槽水流在此急劇繞流[7]，水流與墩壁脫離，出現“卡門渦街”現象，導致尾渦周期性脫落搖擺，形成2條周期性交錯的擠壓逆行波與擴散順行波，并不斷向上游推進（見圖6）。對于上游節制閘來說，相當于在閘門下游出現周期性高低變化的水位頂托，過閘水頭時大時小，從而使得過閘流量也時大時小，這種流量周期性變化反過來又將引起左右2股出槽水流更加不均，加重了墩尾水流的擺動[8]，即節制閘過閘流量的交替變化與墩尾水流擺動是相互激勵的過程[4]，共同造成渡槽水位的劇烈波動，形成相互激勵的共振現象[9]。

此外，通過現狀流態觀測和數值模擬發現，渡槽進出口漸變段拐角位置均存在小范圍的回流區，這是由于進口漸變段收縮角和出口漸變段擴散角較大導致，其中出口回流區范圍遠大于進口。通過現狀和工程措施后的流態分析發現，進出口回流區并未對渡槽水面異常波動等整體水流紊亂產生影響，這點從流態優化措施比選結果可以得到印證，即采用導流墩就可消除渡槽水位異常波動，而進一步采用內外側組合導流墩消除拐角處的回流后，流態并未得到進一步改善。

4 流態改善工程措施研究

4.1 擬定方案

前文通過對現狀流態的數值模擬復演，找出了渡槽流態異常產生的根源，即渡槽出口“卡門渦街”現象直接造成了渡槽水位的異常波動。為消除該流態現象，借鑒其他工程經驗并結合該工程實際，提出了以下各種可行的導流墩方案，擬定了導流墩體型和初步尺寸（見表2），并在相同參數和邊界條件下進行水動力學數值模擬。從中選出既能改善渡槽流態且結構安全、經濟實用的導流墩體型，最后再對選定的導流墩體型進行詳細的尺寸研究，以確定最為合理的尺寸。

4.2 方案比選

4.2.1 體型方案比選

采用與現狀模擬相同的精度和方法構建三維數模，對不同導流墩體型方案在相同參數和邊界條件下（進口設計流量320 m3/s、出口設計水位134.11 m）進行水動力學數值模擬。結果表明：除雙曲扭面導流墩外，其他體型導流墩均可完全消除出口“卡門渦街”現象，消除渡槽內水面異常波動，但在水流平順度、流速變化梯度、水頭損失、流場等方面略有差別，整體流態改善程度略有不同，不同體型導流墩主要模擬成果見表3，主要體型出口水位縱剖面見圖7，出口流速分布見圖8。

通過從不同角度對各導流墩體型對比分析，主要結論如下：

（1） 除雙曲扭面導流墩外，其他體型導流墩均可消除出口“卡門渦街”和槽內水面異常波動。

（2） 從縱橫向水位、流速變化比較渡槽出口水流平順度為：內外側組合>等高扭面>流線型>三角形>實用型>三角+直墻型>中墩順延>雙曲扭面>現狀。

（3） 等高扭面水頭損失最小，三角形、流線型、實用型、中墩順延、內外側組合水頭損失均非常接近，相差不超過0.005 m，三角+直墻組合導流墩略大，雙曲扭面導流墩最大。

（4） 內外側組合導流墩可消除出口外側拐角位置的回流現象;中墩順延導流墩可使回流區范圍明顯減小;而等高扭面導流墩由于中心出流最集中，回流區范圍反而最大;雙曲扭面回流區范圍同現狀一樣隨尾流擺動呈不確定性;其他體型導流墩回流區范圍則相差不大。

綜上，從消除渡槽水面異常波動并最大限度改善流態考慮，由好到壞的選擇為：內外側組合>等高扭面>流線型>三角形>實用型>三角+直墻型>中墩順延>雙曲扭面，但從投資、施工難度等方面綜合考慮，仍推薦采用三角形導流墩體型。

4.2.2 尺寸方案比選

采用相同的三維數模，對不同長度的三角形導流墩在相同參數和邊界條件下進行水動力學數值模擬，以水流平順度、水頭損失、槽內水位波動等作為評價指標，以確定合理的導流墩長度。根據澧河渡槽工程布置，進口導流墩長度在5～20 m范圍內選取，出口長度在10～40 m范圍內選取。3E8F2EEA-38EE-4E7D-9CEE-5A4D1C0E46D2

（1） 進口導流墩尺寸比較。

進口導流墩主要對比墩頭墩尾水位落差、水頭損失（進口局部段）兩項指標，對比結果見表4。

（2） 出口導流墩尺寸比較。

出口導流墩主要對比出口流態、槽內水位波動、水流平順度和總水頭損失幾項指標。

流態對比結果顯示：出口安裝10 m長導流墩后， “卡門渦街”現象仍未消失，但與現狀相比尾流擺動幅度和頻率有所降低;安裝20 m長導流墩后 “卡門渦街”現象基本消失，但仍存在微小擺動;設置30 m和40 m導流墩后尾流擺動現象完全消失。隨著導流墩的加長，左右側流速分布更加對稱，出口不同尺寸三角形導流墩流速分布見圖9。

通過選取槽內典型特征點（渡槽中點），對比出口設不同長度導流墩時該特征點的水位波動，結果顯示10 m長導流墩由于未完全消除尾流擺動，使得槽內水位波動仍達到0.2 m左右，而20，30，40 m長的導流墩則均可使槽內水位波動降低到0.1 m以內，彼此并無明顯差別，出口設置不同尺寸三角形導流墩后槽內水位波動見圖10。

從水流平順度來看，隨著出口導流墩的加長，渡槽出口水躍幅度逐步降低，高低速水流銜接更加平順;從水頭損失方面比較，10 m長導流墩由于未能徹底解決尾流擺動問題，水頭損失最大，導流墩長度從20 m逐步增加到40 m，可使水頭損失進一步減小，但減小幅度非常有限，出口設不同長度三角形導流墩后水躍和水頭損失對比見表5。

通過從不同角度對進、出口導流墩長度進行對比分析，得出主要結論如下：

① 隨著進出口導流墩的加長，水流均變得更加平順;

② 進口導流墩長度在15 m時，水頭損失最小，出口導流墩長度在40 m時，水頭損失最小;

③ 出口采用20 m長導流墩后，“卡門渦街”現象基本消失，但尾流仍存在較小幅度的輕微擺動，采用更長的30 m和40 m長導流墩后，尾流擺動現象徹底消失，可完全解決現狀流態問題;

④ 出口采用20 m及以上長度導流墩后，槽內已無水位異常波動，20 m及以上尺寸導流墩槽內水位波動幅度無明顯差異;

⑤ 隨著出口導流墩的加長，左右側流速分布更加對稱。

綜上，僅從流態改善角度考慮，推薦采用“進口長15 m+出口長40 m”的導流墩方案。但該工程綜合考慮性價比、結構穩定性等其他方面后最終仍采用了“進口長10 m+出口長30 m”的導流墩方案。

5 結 論

（1） 渡槽出口中墩墩尾水流急劇繞流出現“卡門渦街”現象是造成槽內水位異常波動的根源，而由其造成的過閘流量周期性變化又加劇了這種現象，形成相互激勵的共振過程，放大了槽內水位的波動。

（2） 三角形導流墩可完全消除出口“卡門渦街”現象，解決槽內水面異常波動的主要問題，流態雖然略差于使用內外側組合、等高扭面、流線型等體型導流墩，但從工程量、施工難度等方面綜合考慮，仍推薦采用該方案。

（3） 從水流平順度、水頭損失等各方面比較，“進口長15 m+出口長40 m”的導流墩尺寸方案流態最優，但從工程性價比、結構穩定性等其他方面綜合考慮，最終仍推薦采用“進口長10 m+出口長30 m”的導流墩方案。

參考文獻：

[1] 王才歡，王偉，侯冬梅，等.大型輸水渡槽水流異常波動成因分析與對策[J].長江科學院院報，2021，38（2）：46-52.

[2] 張洛，后小霞，楊具瑞.邊寬尾墩體型對邊墻區域水流水力特性的影響研究[J].水力發電學報，2015，34（1）：85-92.

[3] 陳景仁.湍流模型及有限分析法[M].上海：上海交通大學出版社，1989.

[4] 蘇銘德，康欽軍.亞臨界雷諾數下圓柱繞流的大渦模擬[J].力學學報，1999，31（1）：100-105.

[5] 王振東.漫話卡門渦街及其應用[J].力學與實踐，2006，28（1）：88-90.

[6] 董雙嶺，吳頌平.關于卡門渦街形狀穩定性的一點分析[J].水動力學研究與進展A輯，2009，24（3）：326-331.

[7] 黃鈺期，鄧見，任安祿.黏性非定常圓柱繞流的升阻力研究[J].浙江大學學報（工學版），2003，37（5）：596-601.

[8] 張波，何國毅，班學.流場中柔性體擺動尾跡模態分析[J].南昌航空大學學報（自然科學版），2014，28（2）：18-23.

[9] 葛秀坤，趙慶賢，邵輝，等.繞流對稱性破缺誘發柔細結構振動研究：以Tacoma大橋風振致毀事故為例[J].防災減災工程學報，2011，31（4）：430-435.

（編輯：胡旭東）

Causal analysis on water level abnormal fluctuation in aqueduct and improvement measures：case

of Lihe aqueduct in Middle Route of South to North Water Transfer Project

QU Zhigang，LI Zhengpeng

（Henan Water & Power Engineering Consulting CO.，Ltd.，Zhengzhou 450008，China）

Abstract：

The Middle Route of South to North Water Diversion Project officially started large-flow water diversion in 2020.During this period，it was found that the flow pattern at inlets and exits of some conveyance aqueduct were disordered，the water level in the aqueduct fluctuated abnormally，and the diversion capacity was reduced.Taking Lihe aqueduct as the research object，high precision three-dimensional hydrodynamic numerical model was used to simulate the current flow abnormal phenomenon in Lihe aqueduct.The middle pier at aqueduct exit was relatively flat，which caused water flow sharply around and the Carmen vortex street phenomenon appeared.Pier tail flow swang periodicly causing periodic change of the two streams flow，which in turn increased the swing of the pier tail flow.The two flows motivated each other to form resonance，and enlarged the abnormal fluctuation of the water level.Combined with the actual engineering，different kinds of feasible flow regime optimization engineering measures were put forward.After numerical simualtion analysis，the most suitable engineering measures for the Lihe aqueduct were obtained，namely constructing triangular diversion piers at the upstream and downstream of the middle pier at the aqueduct entrance and exit.Considering the engineering cost performance and structural stability，it is recommended to adopt the size scheme of inlet length 10 m+outlet length 30 m.The research results can provide a reference for the flow regime improvement in the aqueduct of similar water diversion projects.

Key words：

conveyance aqueduct;abnormal fluctuation of water level;carmen vortex street;numerical simulation;flow pattern optimization;Middle Route of South to North Water Transfer Project3E8F2EEA-38EE-4E7D-9CEE-5A4D1C0E46D2