閆 科

(新疆伊犁河水利水電投資開發(fā)(集團(tuán))有限公司，新疆 烏魯木齊 830000)

引水渠是水利工程中的重要渠道構(gòu)筑物，起著輸送水源、引水發(fā)電、農(nóng)業(yè)灌溉等作用。在設(shè)計(jì)引水渠的斷面形狀時(shí)需要考慮地形因素、地質(zhì)條件等，不過(guò)最常見的仍是梯形斷面。針對(duì)梯形引水渠道施工與運(yùn)營(yíng)的問(wèn)題，眾多科研工作者通過(guò)理論解析法、有限元法、室內(nèi)試驗(yàn)法、現(xiàn)場(chǎng)監(jiān)測(cè)法等不同手段，展開了一系列的研究。周蔚蔚等[1]以某梯形引水渠為研究對(duì)象，通過(guò)COMSOL數(shù)值軟件建立梯形渠道模型，系統(tǒng)探討了不同斷面尺寸對(duì)渠道凍深、凍脹變形及應(yīng)力等影響，明確了合理梯形渠道斷面尺寸以提升渠道的抗凍脹性。肖旻等[2]在Winkler地基梁理論的基礎(chǔ)上，同時(shí)考慮了凍土與渠道接觸面作用，建立了梯形渠道彈性凍土地基梁撓曲線方程。通過(guò)假設(shè)板間連續(xù)性和引入實(shí)際邊界條件求解了該方向。以實(shí)際梯形渠道為例，分析了渠道不同截面的凍脹位移、表面應(yīng)力以及凍結(jié)李。王功[3]首先提出對(duì)邊坡系數(shù)與渠道斷面尺寸相關(guān)關(guān)系的研究中存在的問(wèn)題，對(duì)此，基于導(dǎo)數(shù)求導(dǎo)手段，求解獲得了梯形渠道水力最佳斷面函數(shù)。對(duì)函數(shù)進(jìn)一步探究，準(zhǔn)確分析了邊坡系數(shù)與渠道寬深比的關(guān)系，明確了不同條件下的水力最佳斷面。劉裕[4]以喀什噶爾河灌區(qū)阿瓦提干渠為研究對(duì)象，基于渠道襯砌結(jié)構(gòu)凍脹破壞機(jī)制，建立了其力學(xué)模型，求解獲得梯形渠道的凍脹力和凍脹位移。進(jìn)一步明確了所建立的力學(xué)模型與土體參數(shù)無(wú)直接關(guān)系。歐祖賢等[5]推導(dǎo)出渠道壁面加糙下岸邊流速系數(shù)方程，并對(duì)截面尺寸、水流流速等因素的影響規(guī)律進(jìn)行系統(tǒng)探討。研究表明，岸邊的流速系對(duì)渠道流量有顯著影響。宋禹德等[6]建立梯形渠道分流量方法，分析了主渠分水口水流形態(tài)。并通過(guò)分流試驗(yàn)分析了分水口水頭損失規(guī)律和水面變化規(guī)律?；谀程菪我汕?，通過(guò)有限元軟件建立三維數(shù)值模型，對(duì)比了4種彎曲度(15°、45°、60°和90°)下的水面線分布、水流流速等，總結(jié)了渠道彎曲度對(duì)水流流動(dòng)特性的影響規(guī)律。

1 工程概況

本文所研究的某引水干渠，全長(zhǎng)1287m，設(shè)計(jì)流量為每秒2.22m3，是該地區(qū)引水干渠中重要的一條引水渠道。渠道橫截面呈“倒梯形”，渠深2m，渠底寬2.9m，上口寬4.2m，渠底襯砌厚度12cm、渠坡襯砌厚度為12～10cm，混凝土標(biāo)號(hào)為C25，抗凍標(biāo)號(hào)為F100，抗?jié)B標(biāo)號(hào)為W6；渠底、渠坡墊層置換厚度取值均為50cm；渠道每3.3m設(shè)一道橫向伸縮縫，伸縮縫縫寬1.5cm，采用聚乙烯閉孔泡沫板夾縫，縫內(nèi)表層(過(guò)水面)灌注厚2.5cm聚氨酯密封膠防滲。

2 三維數(shù)值模型

參考該引水干渠的實(shí)際尺寸，建立了截面為“倒梯形”的彎曲干渠，梯形上寬4.2m、底寬2.9m、高度2m，彎曲處曲率半徑為8m，渠道直線段均為20m，三維數(shù)值模型如圖1所示。

圖1 三維數(shù)值模型圖

為研究彎曲度的影響，分別建立了4種彎曲度α工況：15°、45°、60°和90°。用四面體對(duì)模型進(jìn)行網(wǎng)格劃分，15°工況下共劃分了8512個(gè)網(wǎng)格、48215個(gè)節(jié)點(diǎn)；45°工況下共劃分了9621個(gè)網(wǎng)格、53251個(gè)節(jié)點(diǎn)；60°工況下共劃分了10628個(gè)網(wǎng)格、60012個(gè)節(jié)點(diǎn)；90°工況下共劃分了10723個(gè)網(wǎng)格、60158個(gè)節(jié)點(diǎn)。進(jìn)口處水流流速設(shè)為1.71m/s，水體自由流動(dòng)。

3 計(jì)算結(jié)果分析

根據(jù)計(jì)算結(jié)果可知，水流在進(jìn)入直線段已處于穩(wěn)定狀態(tài)，不同工況下水體流入彎道的條件相同。凹岸側(cè)水流流速低于凸岸側(cè)水流流速。

水面線分布隨觀測(cè)點(diǎn)距凹岸的距離的變化曲線如圖2所示。

圖2 水面線分布隨觀測(cè)點(diǎn)距凹岸的距離的變化曲線

對(duì)橫坐標(biāo)距離進(jìn)行無(wú)量綱化處理(以水面線總寬度對(duì)其進(jìn)行無(wú)量綱化)。每種彎曲度工況下分別給出6個(gè)特征斷面結(jié)果。每種彎曲度工況下不同特征斷面的水深結(jié)果見表1。從圖中可以看出，由于干渠彎道的存在，水流產(chǎn)生了壅水現(xiàn)象，即水面由凸岸向凹岸逐漸抬高。對(duì)比4個(gè)彎曲度工況可知，不同彎曲度工況下特征斷面水面線變化規(guī)律存在顯著差別。隨著渠道彎曲度減小，壅水現(xiàn)象發(fā)生位置逐漸向下游移動(dòng)。且最低水深出現(xiàn)的位置始終為在彎道出口處，而最大超高水深發(fā)生的位置向彎道出口處靠近。從表中可以看出，4中彎曲度工況下，最大超高水深基本都發(fā)生在凹岸1.6m水位的位置，可見，最大超高水深受渠道彎曲度的影響有限。

表1 不同特征斷面的水深結(jié)果

彎道中心位置的垂線流速隨水深的變化曲線如圖3所示。

圖3 彎道中心位置的垂線流速隨水深的變化曲線

縱坐標(biāo)水深進(jìn)行無(wú)量綱化處理(以直線段水深平均值對(duì)其進(jìn)行無(wú)量綱化)。每種彎曲度工況下分別給出6個(gè)特征斷面的變化曲線。從圖中可以看出，不同工況下的垂線流速沿水深先緩慢增加，隨后保持穩(wěn)定，最后在靠近渠道底部處迅速減小。具體的，水流流速最大值基本都發(fā)生在z/H為0.7附近；z/H處于0.7～1.0范圍內(nèi)，水流流速基本緩慢增大；z/H處于0.4～0.7范圍內(nèi)，水流流速基本都處于穩(wěn)定狀態(tài)；z/H處于0～0.2范圍內(nèi)，水流流速基本都顯著減小。

當(dāng)渠道彎曲度為15°時(shí)，水流發(fā)生的壅水現(xiàn)象不太明顯，水流流速隨著特征斷面角度增大而逐漸增大。隨著渠道彎曲度的增大，水流發(fā)生的壅水現(xiàn)象逐漸明顯，水流流速隨著特征斷面角度的變化無(wú)明顯規(guī)律。水流最大流速隨渠道彎曲度的增大而逐漸減小，彎道前半段的水流流速隨渠道彎曲度的增大而持續(xù)減小，彎道后半段的水流流速隨渠道彎曲度的增大而持續(xù)增大。除渠道彎曲度為15°工況的流速峰值發(fā)生在彎道進(jìn)口斷面，其余彎曲度工況的流速峰值發(fā)生在彎道出口斷面。這是因?yàn)楫?dāng)渠道彎曲度增大時(shí)，彎道內(nèi)水流受受壁面作用增大，水流受阻水位抬升現(xiàn)象愈加明顯，并且壅水現(xiàn)象發(fā)生位置逐漸向上游移動(dòng)，因而水流整體流速處于減小狀態(tài)，彎道入口處水流流速也減小。水流的壅水現(xiàn)象逐漸向上游移動(dòng)導(dǎo)致彎道后半段發(fā)生落水現(xiàn)象，因而彎道后半段水流流速變大，并再?gòu)澋莱隹谖恢眠_(dá)到峰值。

不同觀測(cè)斷面0.5m水深的水流流速隨觀測(cè)點(diǎn)距凹岸的距離的變化曲線如圖4所示。

圖4 水流流速隨觀測(cè)點(diǎn)距凹岸的距離的變化曲線

橫坐標(biāo)距離同樣進(jìn)行無(wú)量綱化處理。每種彎曲度工況下分別給出6個(gè)特征斷面的變化曲線。從圖中可以看出，不同彎曲度工況下主流都有遠(yuǎn)離凹岸向凸岸移動(dòng)的趨勢(shì)，即凹岸側(cè)水流流速小于凸岸側(cè)水流流速。但兩側(cè)水流流速在彎道入口位置接近。當(dāng)水流流進(jìn)彎道范圍內(nèi)，兩側(cè)水流流速在壁面作用下差距逐漸變大，主流靠近凸岸一側(cè)。水流在流過(guò)彎道頂后，凸岸一側(cè)的流速繼續(xù)增大并達(dá)到峰值。這是因?yàn)樗髁鬟M(jìn)彎道范圍內(nèi)，凹岸一側(cè)水流產(chǎn)生壅水現(xiàn)象導(dǎo)致水位抬升，凸岸一側(cè)水位降低，水位抬升區(qū)域內(nèi)的水流流速低于水位降低區(qū)域內(nèi)的水流流速，而水位抬升區(qū)域內(nèi)的水流擠壓主流致其在彎道入口處產(chǎn)生向凸岸一側(cè)移動(dòng)的趨勢(shì)，在彎道出口處，主流受凹岸區(qū)域內(nèi)抬升水位的擠壓流速進(jìn)一步加大。

4 數(shù)值與監(jiān)測(cè)結(jié)果對(duì)比

選取引水干渠某彎曲度約為45°位置的水流監(jiān)測(cè)結(jié)果與數(shù)值模擬彎曲度為45°工況的計(jì)算結(jié)果進(jìn)行對(duì)比。7個(gè)觀測(cè)點(diǎn)10cm水深處的水流流速的對(duì)比曲線、彎道入口和出口2處對(duì)比結(jié)果如圖5所示。

圖5 監(jiān)測(cè)值與計(jì)算值對(duì)比曲線

從圖中可以看出，彎道范圍內(nèi)水流流速的監(jiān)測(cè)值與計(jì)算值吻合較好，很好的證明了本文所建立模型的正確性。同時(shí)實(shí)際的水流流速略小于計(jì)算的水流流速，這可能是因?yàn)橐汕嬖谏倭可匙?，阻礙水體留點(diǎn)，導(dǎo)致實(shí)際水流流速較低。

5 結(jié)論

(1)由于干渠彎道的存在，水流產(chǎn)生了壅水現(xiàn)象，即水面由凸岸向凹岸逐漸抬高。且最低水深出現(xiàn)的位置始終為在彎道出口處，而最大超高水深發(fā)生的位置向彎道出口處靠近。

(2)垂線流速沿水深均先緩慢增加，隨后保持穩(wěn)定，最后在靠近渠道底部處迅速減小。小彎曲度干渠的流速峰值發(fā)生在彎道進(jìn)口斷面，大彎曲度干渠的流速峰值發(fā)生在彎道出口斷面。

(3)不同彎曲度工況下主流都有遠(yuǎn)離凹岸向凸岸移動(dòng)的趨勢(shì)。彎道范圍內(nèi)兩側(cè)水流流速在壁面作用下差距逐漸變大，主流靠近凸岸一側(cè)。水流在流過(guò)彎道頂后，凸岸一側(cè)的流速繼續(xù)增大并達(dá)到峰值。

(4)由于模型中采用了近似于渠道空流的狀態(tài)，因此模型僅適用于剛開閘通水階段。