均質透水丁壩二維水流模擬研究

游立新，許 慧，魏 帥，陳梓偲

(1.中交第二航務工程勘察設計院有限公司，武漢 430060；2.水利部交通運輸部國家能源局南京水利科學研究院，南京 210029；3.山東省淮河流域水利管理局規劃設計院，濟南 250100；4.河海大學，南京210098)

丁壩作為常見的航道整治建筑物之一，其作用主要是整治淺灘、穩定航槽、確保航深。根據透水性能的不同，可分為實體和透水丁壩。實體丁壩壩頭處水流分離嚴重，壩頭沖深較大，壩后水流復雜，而透水丁壩具有一定的透水率，可減少壩頭水流分離程度和大漩渦的產生，減弱壩頭局部沖刷深度[1-3]。此外，透水丁壩具有一定孔隙率，可為魚類和水生生物提供庇護和棲息環境，有利于生態環境多樣化發展，對環境保護起到積極作用。因此，透水丁壩是一種極具發展前景的整治建筑物。常見的透水丁壩有三類：(1)由樁、木頭構成的樁柱式丁壩；(2)由鉛絲籠或土工織物布包裹粗砂材料構成的丁壩；(3)在實體丁壩(不透水丁壩)中布設過水管。

1 背景

自1970年黃河中下游河道治理開始試用透水丁壩以來，目前全國各大流域的河道治理及主要航道整治工程中均出現有透水丁壩工程。透水丁壩附近水流流動是一個三維問題，表現為壩頭集中繞流、下潛流及壩頭后方豎軸漩渦這三者在一定弗勞德數范圍內相互作用而形成的復雜渦街結構。在實際工程中，往往更關心建壩后大范圍水位、流場變化及影響范圍，故采用二維水流數學模式已經可以滿足實際需要。許多學者模擬了實體丁壩附近二維水流運動[4-11]，而對于透水丁壩，專門的研究則相對較少[12-13]。本文基于正交曲線貼體坐標系下平面二維水流數學模型方程，通過引入通度系數的概念考慮丁壩的透水率，對考慮通度系數后的方程進行離散和求解，并采用水槽試驗資料對模型進行了驗證。

2 二維數模中透水丁壩的處理

2.1 控制方程組

正交曲線貼體坐標系下，水流運動控制方程式可表示為

水流連續方程

(1)

ξ方向動量方程

(2)

η方向動量方程

(3)

(4)

(5)

(6)

2.2 通度系數概念的引入

相比于實體丁壩，透水丁壩壩體具有一定的透水率。透水率是指壩體透空體積與丁壩體積之比，反映了透水丁壩結構阻水特性，在數值模擬時透水率需作為重要影響參數加以考慮。受戴會超[14]和金忠清[15]用通度系數的概念處理叉管不規則邊界的啟發，引入通度系數的概念處理透水丁壩。通度系數是指流體可通過區域在整個區域中所占的比例，在計算區域離散化以后，網格單元會有不同的通度系數。通度系數又可分為體通度系數和面通度系數兩種[14-15]。體通度系數表示一個單元內流體所占據的體積分數，用PV表示；面通度系數則表示一個單元的一個側面上可供流體通過的面積分數。對于三維問題，一個控制體有6個面通度系數，用Pe、Pw、Pn、Ps、Pt、Pb表示；對于二維問題，一個控制體有4個面通度系數，分別用Pe、Pw、Pn、Ps表示。

由于丁壩的寬度跟計算區域長度相比較小，完全可以將丁壩概化為占有一個網格單元寬度的一列控制體。各水力要素布置如圖1所示。

圖1 水力要素網格布置圖Fig.1 Layout of hydraulic elements grid 注：○表示水位Z節點位置，→表示u節點位置，↑表示v節點位置，□表示水深H節點位置。

丁壩透水以后，有效阻水面積減少，阻水作用降低。假設丁壩在某一方向的透水率為φ，則丁壩有效阻水面積=丁壩投影斷面面積×φ。考慮沿水流方向和垂直水流方向透水率一致的情況下，則

(7)

式中：DG表示壩高，(Z-HD)P、(Z-HD)w仍代表當地網格中心控制點和網格線上水位與地形差值，其余方向類似定義。上式對淹沒透水丁壩也是適用的，丁壩壩頂以上部分完全通水，以Pw為例，通度系數隨水位Z的抬高而增大。

2.3 考慮通度系數后的方程離散和求解

將通度系數視作空間分布的變量，在交錯網格中，體通度系數PV與其他變量布置在主控制體中，面通度系數Pe、Pw、Pn、Ps與速度分量一起布置在主控制體的相應側面上(圖2)。

圖2 通度系數和其他變量在 交錯網格中布置Fig.2 Layout of block up coefficient and other variables in staggered grids

水流運動方程組(1)、(2)、(3)可寫成如下通用微分方程

(8)

式中：Φ表示通用變量， Γ為廣義擴散系數，S為源項。

根據守恒定律，推導單元離散方程時應在相應的體積分上乘以該單元的體通度系數PV；在相應的面積分上乘以相應交界面的面通度系數Pe、Pw、Pn、Ps，由此離散方程變為

(9)

式中：

(10)

(11)

(12)

(13)

(14)

(15)

各控制面上對流、擴散強度及派克里特數Fe、De、Pe、Fw、Dw、Pw、Fn、Dn、Pn、Fs、Ds、Ps表達式見文獻[10]。

采用水位校正法SIMPLER程式求解耦合水位-速度場，求解代數方程組時仍采用交替方向隱式迭代法，反復迭代直至流場結果滿足精度，通過控制水位校正方程中質量源項的最大值來確定收斂情況。

3 模型驗證

3.1 水槽布置

采用王振[16]2012年鉛絲籠透水丁壩定床水槽試驗結果對采用通度系數法的平面二維數學模型進行驗證。試驗在總長為20 m，寬0.5 m，高為0.44 m的矩形玻璃水槽中進行，試驗段選取在水槽中間6～12 m范圍內。經率定水槽綜合糙率n取0.009 9，進出口水位流量條件由三角堰和尾門控制，其中進口流量設置為0.015 m3/s，出口水深為0.065 m。沿水流方向在水槽左側距進口9 m處布置正挑透水丁壩模型，丁壩尺寸長L=0.15 m，寬B=0.05 m，壩高H0=0.10 m。鉛絲籠丁壩沿水流方向和垂直水流方向透水率一致，試驗測得透水率均取20%。丁壩平面位置及測量斷面布置如圖3所示，每個斷面布設5～11根垂線。

圖3 丁壩平面布置和測量斷面布置Fig.3 Plane layout of spur dike and survey section

3.2 流場和水位對比

圖4 透水丁壩附近平面流場圖Fig.4 Flow field near the permeable spur dike

透水丁壩附近區域流場如圖4所示。由圖可見，透水丁壩作用后的水流流態與實體丁壩類似，壩頭形成分離流，壩后形成回旋流。在壩軸斷面，由于丁壩具有一定的透水性，少量的水通過壩體間的空隙流向下游，但流速較小；在壩頭，由于丁壩的束窄作用使壩軸斷面單寬流量增大，主流流速增大；在丁壩軸線下游側，由于丁壩的影響，在壩后形成回流區，主流由于慣性作用進一步收縮，流速繼續增大，到最小收縮斷面處達到最大，隨后主流擴散，在一定距離后重新趨于均勻。

丁壩上下游各個斷面的流速(u/u0)、水深(h/h0)驗證見圖5和圖6，其中u、h為測點流速、水深；u0、h0為無壩時斷面平均流速、水深；坐標原點設置在透水丁壩壩根處，x軸方向為水流運動方向，y軸方向為垂直水流方向，并指向河內。由圖可見，計算值和實測值符合較好，模型能夠反映透水丁壩作用影響下流速和水位分布。在丁壩上游，越靠近丁壩，丁壩側水面逐漸壅高，流速減小，丁壩對岸側流速則相應增加，丁壩的阻水作用逐漸增強；當水流到達丁壩位置時，丁壩束窄作用最強，比降和流速迅速增大，丁壩處水面出現明顯降落；水流流過丁壩后，主流流速仍有繼續增大趨勢，壩后出現與水流流向相反的流速，形成回流，隨著斷面向下游移動，主流逐漸擴散到無丁壩影響的斷面水流狀態。通過驗證，說明采用通度系數可較好模擬沿水流方向與垂直水流方向透水率一致的丁壩附近水流運動。

圖5 透水丁壩斷面流速計算與實測對比Fig.5 Comparison between calculation and measurement of flow velocity in permeable spur dike section

圖6 透水丁壩縱向水面線計算與實測對比Fig.6 Comparison between calculation and measurement of longitudinal water surface profile of permeable spur dike

4 不同透水率對丁壩附近水流影響

圖7 不同透水率時丁壩附近流場圖Fig.7 Flow field near spur dike with different permeability

不同透水率的丁壩對壩后回流的影響是不同的，其保護范圍也不一樣。為研究不同透水率對丁壩附近流速分布影響，在前述3.1節水槽和丁壩尺度基礎上，以不透水丁壩(透水率為0%)作為對照，模擬5組不同透水率(0%、20%、40%、60%、80%)工況下丁壩附近流場。

當透水率為0時，丁壩上游流速有較明顯的減小，下游壩后有明顯的回流區；隨著透水率逐漸增加，丁壩阻水作用減弱，相應地透過壩體的水流增加，與壩后水流相互碰撞，導致壩后回流區域逐漸縮小，流速增大；當透水率增至80%時，丁壩上游斷面流速分布均勻，壩后回流現象基本消失。

圖8 不同透水率時丁壩附近斷面 流速對比圖Fig.8 Comparison of cross-sectional flow velocity near spur dike with different permeability

由丁壩軸線斷面流速分布可看出(圖8)，不同透水率時流速比(其中U0為無壩時流速，橫坐標為距丁壩側壁面的相對位置)均大于1，原因是透水率增加后，一部分流量通過壩體間隙穿過，壩頭流速比下降。可見，丁壩透水后能有效降低壩頭流速，緩解壩頭沖刷現象。

由壩下游典型斷面流速分布可看出(圖8)，丁壩透水率越低，回流區長度越長，回流區域越大；不同透水率時丁壩對岸側流速比差別比較明顯，流速比為1.7左右；透水率較大時流速比維持在1.1左右，說明隨透水率增加，丁壩挑流作用逐漸減弱，過壩后的水流逐漸趨于均勻。

5 結語

(1)基于正交曲線貼體坐標系下平面二維水流數學模型，引入通度系數的概念考慮丁壩的透水率，并對考慮通度系數后的方程進行詳細的離散推導和求解。

(2)采用已有的水槽試驗資料，對模型進行了初步驗證，模擬透水丁壩作用后的流場和水位場，結果表明計算和實測符合較好。丁壩透水率越大，阻水緩流作用越弱，丁壩束窄區域主流區流速降低，壩軸及下游各斷面流速變化幅度減小；壩后回流范圍縮小，流速與主流區流速趨于均勻。

(3)受水槽試驗限制，目前僅模擬了各方向透水率相同的丁壩附近水流運動，下一步將考慮對各部位透水率不同的丁壩水流進行模擬對比。