李天朦，趙 濤，王世興

（新疆農(nóng)業(yè)大學(xué)，新疆 烏魯木齊830052）

渠首運行的主要任務(wù)就是泄洪、引水排沙。目前對渠首工程運行方式的研究大多數(shù)為物理模型試驗方法[1－3]，對閘門調(diào)度的數(shù)值模擬研究主要為輸水渠道的水力響應(yīng)過程[4－8]，其中比較典型的數(shù)學(xué)模型有以圣·維南方程組作為渠道運行控制方程，采用Preismann四點差分隱格式對方程組進行離散化處理，采用追趕法進行迭代求解，并對其閘門附近水位下降速率、波動的傳播時間、水流波動幅度大小和影響范圍進行分析。對渠首工程閘門調(diào)度及運行方式的數(shù)值模擬研究很少[9－12]，主要應(yīng)用 VOF模型，結(jié)合k－ε紊流模型和分段線形方式構(gòu)造界面的方法，模擬大型水工建筑物內(nèi)的水流流動。

本文以新疆某引水渠首為研究對象，該渠首工程為典型的攔河閘式引水渠首[13]，其特點是可根據(jù)河道來水、來沙規(guī)律主動靈活的控制渠首運行。利用泄洪沖沙閘及進水閘的部分或全部開啟調(diào)節(jié)水位和流量，壅水沉沙、泄洪排沙，使取水口始終保持良好的引水條件，達到較優(yōu)的泄洪、引水效果。此前已對該工程進行了物理模擬試驗，獲得了大量實測數(shù)據(jù)，因此，本文在物理模型的基礎(chǔ)之上對該渠首工程的閘門調(diào)度及運行方式進行數(shù)值模擬分析，并以物理試驗結(jié)果為基礎(chǔ)，對數(shù)值模擬的分析結(jié)果進行驗證。

1 模型介紹

物理模型由上游鋪蓋段、閘室段和下游連接段三部分組成:實際工程中上游鋪蓋段（0－020.0～0＋000.0）總長20m；閘室段（0＋000.0～0＋012.5）分為攔河閘、左右岸進水閘，攔河閘包括5孔泄洪閘和2孔沖沙閘，閘室為寬頂堰式，采用弧形閘門；下游段（0＋012.5～0＋208.0）由泄水陡坡、防沖槽以及下游護坡組成，閘室段后接泄水陡坡（長度為25 m），陡坡后設(shè)置防沖槽，防沖槽后為左、右岸的河堤護坡段。該渠首工程的平面布置圖見圖1。

圖1 渠首工程平面布置圖

2 數(shù)學(xué)模型

2.1 VOF方法

目前追蹤自由流體表面的方法主要有標高函數(shù)法、標記網(wǎng)格法和體積率法[14－15]?？紤]計算內(nèi)存、計算時長及準確性等各方面因素，本文采用VOF法模擬自由表面水流。VOF方法是通過體積比函數(shù)F來構(gòu)造和追蹤自由面。若F＝1，則單元全部被指定相流體充滿；若F＝0，則單元沒有指定相流體；當0＜F＜1，則單元為相間界面單元。定義函數(shù)f（x，y，t）:在f＝1時，在（x，y）該相有流體質(zhì)點；在f＝0時，在（x，y）該相沒有流體質(zhì)點。守恒形式的傳輸方程表示為:

2.2 控制方程

本文以連續(xù)方程、動量方程、k方程以及ε方程為基本控制方程組，初步選用k－ε和RNGk－ε兩種紊流計算模型，通過模型驗證確定以k－ε為計算模型，用VOF方法追蹤自由流體表面，采用有限體積法離散控制方程組，并用點隱式高斯－塞德爾迭代方法求解代數(shù)方程組。壓力－速度耦合求解選用PISO算法，動量、紊動能k以及紊動耗散率ε選用一階迎風(fēng)格式，計算精度設(shè)置為單精度。

連續(xù)方程:

動量方程:

k方程:

ε方程:

式中:ρ為體積分數(shù)平均密度；μ為分子粘性系數(shù)；G為由平均速度梯度引起的紊動能產(chǎn)生項；μt為紊流粘性系數(shù)；p為壓力；u為速度；i＝1，2，3，｛xi＝x，y，z｝，｛ui＝u，v，w｝；常 數(shù) 值C1ε＝ 1.44，C2ε＝1.92，δk＝1.0，δε＝1.3。

3 網(wǎng)格劃分與邊界條件

3.1 網(wǎng)格劃分

網(wǎng)格的劃分應(yīng)保證計算精度的同時考慮網(wǎng)格過密對計算時長的不利影響，并且盡量與流線保持一致。網(wǎng)格的結(jié)構(gòu)類型有兩種:結(jié)構(gòu)化網(wǎng)格及非結(jié)構(gòu)化網(wǎng)格。本文采用非結(jié)構(gòu)化網(wǎng)格，并對水面以下的部分進行精細模擬，網(wǎng)格總數(shù)約25萬個。

3.2 進口邊界條件

水流為非恒定流，水流方向與X軸正向一致，Y軸表示水深，Z軸表示寬度。進口邊界條件設(shè)定為:U＝u，V＝0，W＝0，進口斷面的平均流速u由進口流量及過水斷面確定。進口邊界過水斷面以下設(shè)置為速度進口，過水斷面以上設(shè)置為壓力進口。

3.3 出口邊界條件

由于模型出口是非全斷面出流，故結(jié)合工況實際情況將出口斷面分為上下兩個斷面，上斷面設(shè)置為壓力出口，下斷面設(shè)置為自由出流，出口邊界則設(shè)置為壓力出口。

3.4 壁面邊界條件

采用無滑移條件對壁面進行處理:U＝0，V＝0，W＝0，并選用壁面函數(shù)法修正粘性底層內(nèi)的流速。

4 模型驗證

為了比較k－ε和RNGk－ε兩種紊流計算模型在模擬該引水渠首運行方面的優(yōu)劣性，本文在保證兩岸引水時，對設(shè)計流量517m3/s采用RNGk－ε計算模型，對校核流量694m3/s采用k－ε計算模型，并與對應(yīng)的實測數(shù)據(jù)作比較。

4.1 進、出口流量

設(shè)計流量和校核流量下進、出口流量計算結(jié)果與試驗結(jié)果的對比見表1。

表1 進、出口流量計算結(jié)果與試驗結(jié)果對較 單位:m3/s

結(jié)果表明，設(shè)計流量和校核流量下進、出口流量均接近試驗數(shù)據(jù)，并且在設(shè)計流量下，調(diào)整閘門開度可使兩岸引水量滿足渠道設(shè)計流量；校核流量下，在閘門全開的情況下兩岸引水量略低于渠道設(shè)計流量。計算結(jié)果與實驗結(jié)果相比，設(shè)計流量下的最大誤差為7.6%，最小誤差為0.06%，進、出口流量差0.19%；校核流量下的最大誤差為6.6%，最小誤差為0.03%，進、出口流量差0.05%。

4.2 左、右岸水邊線

設(shè)計流量下左、右岸的水邊線計算結(jié)果與試驗結(jié)果的對比，校核流量下左、右岸進水閘全關(guān)閉，兩岸的水邊線計算結(jié)果與試驗結(jié)果的對比見圖2（設(shè)計上校核下）。圖2中12～24分別代表樁號0－120.0～0＋208.0。

圖2 設(shè)計、校核流量下左、右岸的水邊線計算結(jié)果與試驗結(jié)果比較

結(jié)果表明，校核流量下左、右岸的水邊線擬合情況比設(shè)計流量要好，兩者均低于堤頂高程，并且上游整治段（0－120.0～0－080.0）右岸水位高于左岸水位，下游整治段（0＋42.0～0＋128.0）左岸水位高于右岸水位，水位沿程下降，自0＋128.0向下游，水位基本一致。

進、出口流量及左、右岸的水邊線的結(jié)果對比表明，k－ε和RNGk－ε計算模型均能達到收斂條件:（1）殘差曲線在10－3以下，呈平緩曲線；（2）進、出口流量差與進口流量的比值在0.5%以下；（3）斷面特征參數(shù)值隨迭代的進行不再發(fā)生明顯變化。因此k－ε和RNGk－ε計算模型均能夠較好地模擬該渠首工程，但從左、右岸水邊線和時間兩方面考慮，kε比RNGk－ε計算模型少消耗10%～15%的CPU時間，并且水邊線擬合情況k－ε優(yōu)于RNGk－ε計算模型，因此選擇k－ε計算模型對該引水渠首整治段流態(tài)進行分析和典型斷面流速的計算。

4.3 整治段流態(tài)

從圖3和圖4（圖4曲線為跡線）均可以看出，校核流量下左、右岸進水閘全關(guān)閉，河道全斷面過流，由于河道主流軸線與樞紐縱軸線呈一定夾角，水流呈S型進入閘前整治段，受右岸頂托和上游河道內(nèi)橋墩的影響，主流被導(dǎo)向左岸，水流相互疊加，形成菱形波。

圖3 試驗觀察的校核流量整治段流態(tài)

4.4 校核流量下閘室及護坦段流速

校核流量694m3/s下左、右岸進水閘全關(guān)閉，閘前、閘室及護坦段流速計算結(jié)果與試驗結(jié)果的對比見表2，底部流速（距河床0.2m）、斷面平均流速（2/3水深處）、表面流速（距水面0.2m）。

圖4 模擬分析的校核流量整治段流態(tài)

表2 校核流量下閘室及護坦段流速計算結(jié)果與試驗結(jié)果比較 單位:m/s

結(jié)果表明，計算結(jié)果與試驗結(jié)果均可以反映出表面流速最大，并且沿水深方向向下遞減，誤差在20%左右，原因有:（1）在急流情況下，水流波動較大，水位落差顯著；（2）儀器測量時，測桿晃動造成測量誤差；（3）流向偏角、流速脈動對流速測量的影響。

5 渠首運行方式的優(yōu)化

渠首運行方式首先應(yīng)考慮來流量和引水需求量；其次在閘門開啟時，過閘水流應(yīng)保持平穩(wěn)，避免發(fā)生不良流態(tài):如折沖水流、集中水流、回流、漩渦等，并且兼顧下游水位。若閘門開度過大，下泄水流呈遠趨水躍對消能防沖及護坦不利，若閘門開度過小，下泄水流呈淹沒水躍使閘門承受反向沖擊荷載。

當上游河道來流量小于等于253m3/s時，左、右岸進水閘和泄洪沖沙閘全部開啟，左、右岸進水閘前水位低于擋沙坎坎頂高程，不能引水；當上游河道來流量為354m3/s時，左、右岸進水閘和泄洪沖沙閘全部開啟，右岸閘前水位低于擋沙坎坎頂高程，不能引水，左岸閘前水位略高于擋沙坎坎頂高程，左岸引水渠道流量為2.5m3/s，小于設(shè)計引水流量。

因此該渠首工程運行方式原則應(yīng)包含以下三點:（1）在來流量小于等于253m3/s時，應(yīng)通過控制沖沙閘和進水閘的開度滿足渠道引水量；（2）在來流量大于253m3/s時，應(yīng)通過控制沖沙閘、泄洪閘及進水閘的開度滿足渠道引水量；（3）在設(shè)計流量和校核流量下，應(yīng)通過控制沖沙閘和泄洪閘的開度滿足泄洪量。

通過模擬典型流量，得出保證兩岸引水時的具體閘門運行方式見表3。

表3 典型來流量下閘門的運行方式 單位:m

本模型為定床試驗，考慮到發(fā)生設(shè)計洪水和校核洪水時，水流挾帶大量泥沙，必然影響兩岸進水渠道的正常運行，故建議在河道來流量達到設(shè)計洪水及校核洪水時，兩岸進水閘應(yīng)全部關(guān)閉，泄洪、沖沙閘應(yīng)全部開啟進行泄洪排沙。

6 結(jié) 論

（1）利用k－ε模型和VOF法，對新疆某渠首工程閘門調(diào)度與運行方式進行了三維紊流數(shù)值模擬研究，得到沿程水邊線及典型斷面流速。研究表明數(shù)值模擬與物理試驗得出的結(jié)果吻合情況較好，因此驗證了該模型能夠很好的模擬渠道非恒定流。

（2）數(shù)值模擬結(jié)果表明:在洪水期整治段內(nèi)最高水位均低于整治段堤頂高程。優(yōu)化渠首運行方式的方法主要是根據(jù)不同時期的來水量，控制左、右進水閘及泄洪沖沙閘的開啟高度:在來流量小于等于五年一遇洪水時，控制沖沙閘和進水閘滿足渠道引水量；在來流量大于五年一遇洪水時，控制全部閘門滿足渠道引水量；在設(shè)計流量和校核流量下，泄洪沖沙閘全開滿足泄洪量。優(yōu)化后的閘門運行方式，在正常引水期均能滿足渠道設(shè)計引水流量和加大流量要求，在洪水期滿足樞紐泄洪能力要求。

[1]趙經(jīng)華，侍克斌，顧佳俊.呼圖壁河青年渠首工程整體水工模型試驗研究[J].人民長江，2009，40（13）:43－45.

[2]吉永軍，張小帥.葉爾羌河中游渠首樞紐運行管理相關(guān)問題探討[J].水利與建筑工程學(xué)報，2012，10（6）:187－190.

[3]陳成林，馬 亮，趙經(jīng)華，等.引水排沙和泄空排沙在臺蘭河引水渠首中的應(yīng)用試驗研究[J].水資源與水工程學(xué)報，2011，22（6）:100－102.

[4]呂宏興，宋松柏，馬孝義，等.灌溉渠道閘門調(diào)控過程中的非恒定流研究[J].農(nóng)業(yè)工程學(xué)報，2002，18（6）:18－22.

[5]方神光，吳保生，傅旭東.南水北調(diào)中線干渠閘門調(diào)度運行方式探討[J].水力發(fā)電學(xué)報，2008，27（5）:93－97.

[6]張 成，張尚弘.閘門調(diào)控下大型輸水系統(tǒng)的水力響應(yīng)特征研究[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報，2011，19（增刊）:98－108.

[7]張 成，鄭立松.閘門調(diào)控下的水動力學(xué)模型建立與應(yīng)用[J].應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報，2012，20（增刊）:159－168.

[8]丁志良，談廣鳴，陳 立，等.輸水渠道中閘門調(diào)節(jié)速度與水面線變化研究[J].南水北調(diào)與水利科技，2005，3（6）:46－50.

[9]余金星，陳 虹.攔河閘引水口三維流場數(shù)值模擬[J].水利與建筑工程學(xué)報，2007，5（1）:25－29.

[10]李 玲，陳永燦，李永虹.三維VOF模型及其在溢洪道水流計算中的應(yīng)用[J].水力發(fā)電學(xué)報，2007，26（2）:83－87.

[11]張光碧，鄧 軍，張法星，等.VOF模型在有支流匯入的河道復(fù)雜流場中的應(yīng)用[J].河海大學(xué)學(xué)報:自然科學(xué)版，2007，35（5）:592－595.

[12]李志勤，李 嘉，易文敏，等.VOF模型在黃河潘家臺灘整治工程中的應(yīng)用[J].四川大學(xué)學(xué)報:自然科學(xué)版，2004，36（6）:18－23.

[13]邱秀云.水力學(xué)[M].新疆:新疆電子出版社，2004:76.

[14]王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社，2004:123－125.

[15]江 帆，黃 鵬.Fluent高級應(yīng)用與實例分析[M].北京:清華大學(xué)出版社，2008:140－145.