劉易莊，蔣昌波，2，鄧 斌，2，王 剛

(1．長沙理工大學水利工程學院，湖南長沙 410004;2．水沙科學與水災害防治湖南省重點實驗室，湖南長沙 410004;3．中交第四航務工程勘察設計院有限公司，廣東廣州 510230)

丁壩能夠引導水流，防止河岸發(fā)生沖刷而淤高河灘，是一種廣泛使用的航道整治建筑物。水流在繞過丁壩后，水流流場及壓力場都發(fā)生顯著變化［1］，在壩頭附近形成“馬蹄形”漩渦［2］以及一系列的尾流漩渦［3］，水流呈現(xiàn)較強的三維特性。許多學者［1，4-7］對單丁壩的來流特性、回流區(qū)特性、局部沖刷、沖刷坑內水流結構以及單丁壩周邊的流速、雷諾應力分布及紊流特性等方面進行了深入研究。然而，在實際應用中，大多采用丁壩群進行航道整治，而雙丁壩是最簡單的丁壩群，可以將雙丁壩視為研究丁壩群的單元，因此，研究雙丁壩間水流特性對理解壩間泥沙輸移和能量交換、選取合理的壩間距有著重要意義。徐曉東［8］采用聲學多普勒流速儀(ADV)對雙丁壩附近水流結構進行了研究，但ADV作為一種接觸式單點測量技術，必然會對測點周邊水流流場結構造成較大干擾。粒子圖像測速(particle image velocimetry，PIV)是一種瞬態(tài)、多點、無接觸式的流體測速技術，可提供豐富的流場空間結構，因此能較好地彌補ADV測量技術的不足。齊鄂榮等［8-9］采用PIV研究了二維突起物繞流結構，主要分析了漩渦結構運動及發(fā)展過程、流場結構的時均特點;張冠卿等［11-12］采用PIV研究了單丁壩壩后水流特征，對壩后時均流速及渦量場進行了分析;Weitbrecht［13］采用PIV研究了摻混區(qū)水流特性。

以往關于PIV對丁壩的研究一般都集中在單丁壩或丁壩群的流速分布、紊動特性及漩渦結構等方面，本文采用PIV對淹沒雙丁壩間水流流場進行測量，主要分析丁壩間距與丁壩長度的比值(d/B)對壩間平面回流區(qū)流速分布、零速度線以及渦量分布的影響。

1 試驗方法

1．1 試驗條件與試驗裝置

圖1 水槽試驗裝置及PIV系統(tǒng) (單位:m)

圖2 丁壩示意圖

試驗在長沙理工大學水利實驗中心的變坡水槽(試驗中采用平坡)中進行，水槽長16 m，寬0.4 m，高0.5 m。水槽采用水泵供水，在水槽尾部設置調節(jié)水位的格柵尾門，試驗裝置見圖1。丁壩采用透明有機玻璃制成，長 B=5 cm，高 h=18 cm，厚t=1 cm，如圖2所示。流場測試采用美國TSI公司的PIV系統(tǒng)，其中光源為雙脈沖激光，采用強激光(強度約為120mJ)以保證能均勻照亮所測流場，由同步器控制觸發(fā)，采集頻率為1.04 Hz;PIV系統(tǒng)通過CCD(2048×2048像素)記錄圖像數(shù)據(jù)并傳輸?shù)絀NSIGHT3G軟件進行分析，圖像采集間隔為100μs。試驗測量了3個不同相對水深(z/h=1/6、z/h=1/2、z/h=5/6，z軸以水槽底為零點，沿水槽壁往上為正)的水平面流場，測量區(qū)域大小為200 mm×200 mm，試驗測量區(qū)域及坐標系如圖1(a)所示。

1．2 試驗工況

試驗流量Q=57.6 m3/h，水深H=0.20 m，流速us=0.20 m/s，在3種不同壩間距d(d=0.075 m、d=0.15 m、d=∞(表示單個丁壩))的情況下分別測量z=0.03 m、z=0.09 m和z=0.15 m 3個水平面流場。試驗過程中在水槽首部和尾部分別布置超聲波水位計測量試驗水深。

2 試驗結果及分析

圖3 不同相對水深水平面流場(d/B=3)

2．1 水平面流場

為探討淹沒雙丁壩不同相對水深位置的水平面(xOy平面)流場情況，選取壩間距d=0.15 m(d/B=3)、相對水深 z/h=1/6、z/h=1/2、z/h=5/6 的 3 個平面(圖3)進行流場分析。由于上游丁壩的阻擋，壩后形成相對靜水區(qū)，壩頭水流由于丁壩的束窄作用導致流速增大，水流繞過丁壩后向無丁壩一側偏轉，與壩間相對靜水區(qū)形成剪切流，使壩間出現(xiàn)漩渦。從圖3(a)可以看出在壩間近底區(qū)水流較平順，形成了穩(wěn)定規(guī)則的漩渦;隨著相對水深的增大漩渦尺度變小(圖3(b));近水面處壩間水流流速分布忽大忽小、忽左忽右，流態(tài)較為混亂，也幾乎看不到漩渦 (圖3(c))，這是由于水流受到壩頂?shù)挠绊?，離壩頂越近水流流態(tài)越不平穩(wěn)，也越不易形成漩渦。由于在z/h=1/6處存在穩(wěn)定規(guī)則的漩渦，為了更好地闡述雙丁壩間回流區(qū)水流結構特性及渦量分布規(guī)律，選擇z/h=1/6平面作為分析對象平面。

圖4為z/h=1/6平面不同壩間距的水平面流場。壩頭發(fā)生水流分離形成剪切流，使壩后形成漩渦。由圖4可以看出，當d/B=1.5時，漩渦中心在壩頭線之上，隨d/B的增大，漩渦中心向壩頭線移動，當d/B=∞時，漩渦中心幾乎與壩頭線齊平。這說明下游丁壩的存在限制了上游丁壩壩后漩渦的自由發(fā)展，而使得壩間漩渦受到下游丁壩的“擠壓”，并有向丁壩外側溢出的趨勢;隨著壩間距的增大，下游丁壩的“擠壓”作用也越弱，因此漩渦中心向壩頭線移動。

2．2 流速

圖4 不同壩間距水平面流場(z/h=1/6)

圖5為不同壩間距條件下雙丁壩附近x方向流速分布云圖(箭頭線為主流區(qū)與回流區(qū)分界線，虛線橢圓區(qū)域為負流區(qū))，可以看出由于上游丁壩的阻擋，水流立即向無丁壩側偏轉，在壩間形成回流。將過上游丁壩壩頭的流線(箭頭線)作為回流區(qū)的分界線［14］，當 d/B=1.5時，回流區(qū)最大寬度為1.6B;當d/B=1.5時，回流區(qū)最大寬度為1.8B;當d/B=∞時，回流區(qū)最大寬度為2B。這說明隨著d/B的增大，回流區(qū)寬度也隨之增大。此外，壩間回流區(qū)出現(xiàn)負向流速，當d/B=1.5時，僅在下游丁壩頭出現(xiàn)明顯負向流速，隨著d/B的增大負流區(qū)向壩間移動，且負流區(qū)面積(虛線橢圓區(qū)域)和負向流速大小也隨之增大。

圖5 壩后x方向流速云圖 (單位:m/s)

圖6 壩后特征斷面x方向流速分布

圖6 為 x=0、0.84B、1.68B、2.46B、3.24B 處 x方向流速分布(根據(jù)壩間距d=0.15 m平均選取5個斷面，豎線表示特征斷面，圓點與豎線之間的距離表示流速大小，橢圓區(qū)為壩頭逆流)?？梢钥闯鲈谙掠螇晤^處出現(xiàn)明顯逆流(橢圓圈內)，其最大負向流速分別為0.075 m/s(d/B=1.5)和0.068 m/s(d/B=3)，分別占主流流速us的37.5%和34.0%;這是由于下游丁壩仍處在上游丁壩的影響范圍內，減小了上游丁壩壩后回流長度，為了維持壩間流體連續(xù)性［1］，因此壩頭出現(xiàn)逆流;并且d/B越大其負向流速越小，這是由于距離越大，下游丁壩對上游丁壩壩后回流區(qū)影響越小。

壩間回流區(qū)內存在一條零速度線［4］，圖7給出了不同壩間距情況下零速度線位置的沿程變化曲線。由圖7可以看出，零速度線沿程變化可簡化為線性變化，但其斜率各不相同，具體擬合公式如下:

由式(1)可以看出單丁壩壩后零速度線斜率為負，這與陳稚聰?shù)龋?］研究的結論一致;而雙丁壩間零速度線斜率為正，且壩間距越大零速度線斜率越小。

圖7 零速度線沿程變化

為了更加細化分析d/B對壩間回流區(qū)流速的分布影響，結合Weitbrecht［13］的研究方法分別選取x=0.25d、x=0.50d、x=0.75d 3 條測線 x方向流速進行分析，壩間特征斷面流速分布如圖8(圖中橫線表示壩頭線，豎線表示特征斷面，十字點和圓點與豎線之間的距離表示流速大小)所示。回流區(qū)流速沿著丁壩方向先增大后減小(圖8)，且d/B越大，壩間流速也越大。當d/B=1.5時，最大流速點靠近壩頭線區(qū)域，而當d/B=3時，最大流速點則靠近邊壁(圖5(b))。這是由于當d/B=1.5時下游丁壩頭出現(xiàn)較強的回流，壩間漩渦較小，x=0.75d測線穿過的下游丁壩壩后回流流線靠近壩頭，因此在靠近壩頭線處出現(xiàn)較大負向流速;而當d/B=3時，壩間漩渦較大，x=0.75d測線穿過的下游丁壩壩后回流流線靠近邊壁，因此在邊壁處出現(xiàn)最大流速。

圖8 壩間特征斷面流速分布

在壩間回流區(qū)內，兩種壩間距壩間流速相差不大，在靠近下游丁壩(x/d=0.75測線)才出現(xiàn)明顯差異。但在摻混區(qū)，由圖8可以看出，d/B越大，摻混區(qū)流速梯度越小，且進入主流區(qū)的位置離壩頭線越遠(d/B=1.5在y/H=0.375處進入主流區(qū)，而d/B=3在y/H=0.5附近進入主流區(qū))。

2．3 渦量

渦動結構對丁壩附近沖淤以及摻混區(qū)的物質交換有著重要影響，因此有必要對雙丁壩周邊渦量結構進行分析。

垂直方向上渦量的計算公式［15］為

式中:ωz為垂向渦量;u、v分別為 x、y方向流速。

圖9為z/h=1/6時不同壩間距的平面垂向渦量。由圖9可見，水流繞過上游丁壩壩頭后分離，產生渦量值較大(負向)而范圍較小的漩渦，沿著水流方向渦量強度逐漸減小，漩渦呈“喇叭”式逐漸增大，漩渦范圍增加到一定值后，漩渦便分散成許多小塊。同時從圖9還可以看出，d/B越大漩渦越早分散;特別當d/B=1.5時漩渦一直延續(xù)到下游丁壩壩后才分散，從而加速來流水體與壩間相對靜水區(qū)的摻混。以往研究認為正渦量區(qū)只在貼近上下游丁壩壁、水槽壁以及壩間小回流區(qū)出現(xiàn)［15］，但此次研究發(fā)現(xiàn)由于下游丁壩壩頭出現(xiàn)較大逆流，渦量在下游丁壩壩頭附近出現(xiàn)較大范圍的正渦量區(qū)且渦量值最大可達0.02s－1(d/B=1.5)和0.01s－1(d/B=3)。

3 結論

a．水流繞過上游丁壩壩頭后邊界層即發(fā)生分離，壩間形成明顯回流區(qū)，d/B越大回流現(xiàn)象越明顯，漩渦中心越靠近壩頭線。

圖9 壩后渦量云圖 (單位:s－1)

b．壩間回流區(qū)寬度、負向流速區(qū)面積及負向流速大小隨著壩間距的增大而增大。下游丁壩壩頭均出現(xiàn)逆流，當 d/B=1.5時，負向流速最大可達0.075 m/s，占主流流速的37.5%;而在d/B=3時，最大負向流速為0.068 m/s，占主流流速的34%，說明d/B越大逆流流速越小，下游丁壩壩頭附近水體交換強度越弱。

c．雙丁壩間零速度線沿程基本呈線性關系，且其斜率為正并隨d/B的增大而減小。此外，d/B越大壩間回流區(qū)最大流速點越靠近邊壁，摻混區(qū)流速梯度越小。

d．上游丁壩壩頭附近出現(xiàn)較大負渦量值，負渦量區(qū)沿水流方向呈“喇叭”式擴展，而下游丁壩壩頭附近則出現(xiàn)較大正渦量值。

［1］陸永軍，周耀庭．丁壩繞流機理及其下游流場的研究［J］．河海大學學報，1990，18(1):28-36．(LU Yongjun，ZHOU Yaoting．Study on flow mechanism and velocity field near groin-like structure［J］．Journal of Hohai University，1990，18(1):28-36．(in Chinese))

［2］KAMIL H M AL I．Simulation of flow around piers［J］．Journal of Hydraulic Research，2002，40(2):161-174．

［3］陳小莉，馬吉明．受漩渦作用的水下塊石的起動流速［J］．清華大學學報:自然科學版，2005，45(3):315-318．(CHEN Xaoli，MA Jiming．Critical velocity for initiation motion of rocks under water due to vortex motion［J］．JournalofTsinghua University:Science and Technology，2005，45(3):315-318．(in Chinese))

［4］陳稚聰，黑鵬飛，丁翔．丁壩回流分區(qū)機理及回流尺度流量試驗研究［J］．水科學進展，2008，19(5):613-617．(CHEN Zhicong，HEI Pengfei，DING Xiang．Division and flow scale investigation of circulation zone around spur dike［J］．Advances in Water Science，2008，19(5):613-617．(in Chinese))

［5］曹艷敏，張華慶，蔣昌波，等．丁壩沖刷坑及下游回流區(qū)流場和紊動特性試驗研究［J］．水動力學研究與進展:A 輯，2008，23(5):560-570．(CAO Yanmin，ZHANG Huaqing，JIANG Changbo，et al．The experimental research of the flow field and turbulence characteristics in the scour and backflow region around a groin ［J］．Chinese Journal of Hydrodynamics，2008，23(5):560-570．(in Chinese))

［6］DUANJ G．Mean flow and turbulence around a laboratory spur dike ［J］．Journal of Hydraulic Engineering，2009，135:803-811．

［7］CHEN Liping，JIANG Juncheng．Experiments and numerical simulations on transport of dissolved pollutants around spur dike［J］．Water Science and Engineering，2010，3(3):341-353．

［8］徐曉東．非淹沒正交雙體丁壩的水流特性及作用尺度研究［D］．杭州:浙江大學，2013．

［9］齊鄂榮，黃明海，李煒，等．應用PIV進行二維后向臺階流流動特性的研究二維后向臺階流的旋渦結構的研究［J］．水動力學研究與進展:A 輯，2004，19(4):525-532．(QIErong，HUANG Minghai，LIWei，etal Investigation of vortex structure of the 2D backward facing step flow via PIV ［J］．Chinese Journal of Hydrodynamics，2004，19(4):525-532．(in Chinese))

［10］齊鄂榮，邱蘭，盧煒娟，等．二維矩形突起物繞流流動結構的實驗研究［J］．水動力學研究與進展:A輯，2006，21(3):388-394．(QI Erong，QIU Lan，LU Weijuan，et al．Experimental investigation of structure of flow over 2D rectangular obstacle ［J］． Chinese Journal of Hydrodynamics，2006，21(3):388-394．(in Chinese))

［11］張冠卿，林鵬智．明渠流淹沒丁壩周圍流場特性的實驗［J］．水利水電科技進展，2012，32(增刊 2):31-33．(ZHANG Guanqin，LIN Pengzhi．Flow characteristics around open channel submerged spur-dike［J］．Advances in Science and Technology of Water Resource，2012，32(Sup2):31-33．(in Chinese))

［12］賈國珍，佘俊華，王濤，等．淹沒丁壩壩后水流結構PIV試驗研究［J］．水電能源科學，2012，30(3):84-86．(JIA Guozhen， SHE Junhua， WANG Tao， etal． PIV measurement of flow structure behind submerged spur dike［J］．Water Resources and Power，2012，30(3):84-86．(in Chinese))

［13］WEITBRECHT V．Influence of dead-water zones on the dispersive mass transport in rivers［D］．Karlsruhe:Karlsruhe Institute of Technology，2004

［14］樂培九，李旺生，楊細根．丁壩回流長度［J］．水道港口，1999(2):3-9．(YUE Peijiu，LI Wangsheng，YANG Xigen．Length of whirlpool behind dikes［J］．Journal of Water and Harbor，1999(2):3-9．(in Chinese))

［15］白靜，方紅衛(wèi)，何國建．非淹沒丁壩繞流的三維大渦模擬研究［J］．力學學報，2013，45(2):151-157．(BAI Jing，F(xiàn)ANG Hongwei，HE Guojian．Study of nonsubmerged groin turbulence flow in a shallow open channel by LES［J］．Chinese Journal of Theoretical and Applied Mechanics，2013，45(2):151-157．(in Chinese))