柳松濤 李廣寧 石凱 孫雙科 石小濤

摘要：為優(yōu)化魚道池室結(jié)構(gòu)，提升魚類洄游效果，以常規(guī)豎縫式魚道為原型，以齊口裂腹魚（Schizothorax prenanti）為試驗(yàn)對象，通過數(shù)值模擬和對照性過魚試驗(yàn)，比較分析了4種池室內(nèi)部增加與未增加樁柱的水流結(jié)構(gòu)。結(jié)果表明，在魚道主流中間增設(shè)樁柱結(jié)構(gòu)，可以有效減小回流區(qū)面積，減幅達(dá)49%，并可減小豎縫處的紊動能，最大減少18%；增加魚類上溯路線分布，提高了池室空間利用率。2種結(jié)構(gòu)的上溯軌跡分析發(fā)現(xiàn)，齊口裂腹魚多選擇回流區(qū)外側(cè)上溯，避免進(jìn)入回流區(qū)中心，即選擇較高的水力應(yīng)變與紊動能區(qū)域作為上溯方向，且上溯路線中偏好較低水力因子區(qū)域，上溯偏好流速為0.01～0.09 m/s，紊動能為0.001～0.008 m2/s2，水力應(yīng)變?yōu)?.2～3.0/s。體型Ⅱ魚道通過在主流位置增加樁柱結(jié)構(gòu)，形成了多股水流，有效減小回流區(qū)的占比，提高了魚道池室空間利用率，為魚類的上溯路徑提供了更多選擇。

關(guān)鍵詞：豎縫式魚道；樁柱結(jié)構(gòu)；數(shù)值模擬；水力學(xué)特性

中圖分類號：S956? ? ? ? 文獻(xiàn)標(biāo)志碼：A? ? ? ? 文章編號：1674-3075（2024）02-00148-11

受大壩和水庫修建的影響，魚類的洄游通道受到阻隔，作為緩解措施的魚道得到廣泛運(yùn)用（董哲仁，2006；陳凱麒等，2012）；其中，豎縫式魚道（Vertical Slot Fishway，VSF）因其消能效果充分，能適應(yīng)較大水位變幅受到廣泛關(guān)注（劉志雄等，2010；Mao，2018）。然而，我國豎縫式魚道實(shí)際運(yùn)行效果并不理想，主要體現(xiàn)在魚類通過時(shí)間較長且通過率較低。為此，有必要以提高魚道通過效率為目標(biāo)，對常規(guī)池室結(jié)構(gòu)開展進(jìn)一步探索研究。

國內(nèi)外對VSF結(jié)構(gòu)開展了許多研究，Rajaratnam等（1986）最早提出豎縫式魚道設(shè)計(jì)結(jié)構(gòu)，在對18種不同體型池室研究后發(fā)現(xiàn)，10：8為最佳池室長寬比，池室內(nèi)回流區(qū)流速與湍流強(qiáng)度適中，可為魚類提供休息區(qū)（徐體兵和孫雙科，2009；Rajaratnam et al，2010）。對比不同魚道坡度的流場分布發(fā)現(xiàn)，當(dāng)坡度<5%，流場呈現(xiàn)二維分布（Wu，1999）。對豎縫寬度（b）、池室寬度（B）、導(dǎo)板長度（P）等參數(shù)研究發(fā)現(xiàn)，當(dāng)0.1≤b/B≤0.2、0.2

本文提出了4種增設(shè)樁柱結(jié)構(gòu)的魚道池室體型方案，并開展數(shù)值模擬計(jì)算和過魚試驗(yàn)，從而優(yōu)選流態(tài)較好、過魚效率較高的池室體型。對優(yōu)選池室進(jìn)行齊口裂腹魚（Schizothorax prenanti）對比性過魚試驗(yàn)，分析其上溯行為規(guī)律，進(jìn)而識別魚類行為與水力因子間的響應(yīng)關(guān)系，以期為優(yōu)化魚道結(jié)構(gòu)及提升過魚效率提供理論依據(jù)。

1? ?研究方法

1.1? ?試驗(yàn)設(shè)計(jì)

1.1.1? ?模型設(shè)計(jì)? ?魚道模型以常規(guī)豎縫式魚道為原型，按重力相似準(zhǔn)則設(shè)計(jì)，模型比尺為1：4。試驗(yàn)水槽長10.2 m、寬0.6 m，底坡坡度J=1%，設(shè)置11級試驗(yàn)池室。魚道池室長L=0.75 m，寬 B=0.6 m，導(dǎo)板長P=0.15 m，豎縫寬b=0.075 m，導(dǎo)向角θ=45°，導(dǎo)板厚d和隔板厚b均為0.075 m。導(dǎo)板迎水面和隔板背水面坡度均為1：1，隔板迎水面和導(dǎo)板背水面坡度均為1：3（Li et al，2021）。設(shè)計(jì)4種池室，分別為體型I、Ⅱ、Ⅲ、Ⅳ，其差異為樁柱數(shù)量及位置不同，樁柱直徑為豎縫寬度b=0.075 m，樁柱布置及相關(guān)參數(shù)見圖1。

試驗(yàn)?zāi)Ｐ陀盟捎米匝h(huán)閉合系統(tǒng)，模型下游設(shè)回水池，試驗(yàn)時(shí)采用水泵供水，進(jìn)水池中設(shè)置消能柵，魚道池室運(yùn)行水深H0=0.3 m。模型下游末端池室為放魚池，內(nèi)部上下游各安裝一道攔魚柵，上游側(cè)攔魚柵可防止試驗(yàn)魚提前進(jìn)入試驗(yàn)段，下游側(cè)攔魚柵可防止試驗(yàn)魚進(jìn)入回水池。

1.1.2? ?試驗(yàn)方法? ?試驗(yàn)所用齊口裂腹魚采購自四川某漁場，選擇體長（10±2） cm的幼魚，該體長下其臨界游泳能力大于豎縫處最大流速0.39 m/s，確保不會因試驗(yàn)魚游泳能力不足影響試驗(yàn)結(jié)果（Cai et al，2018）。試驗(yàn)前暫養(yǎng)于8 m3的方形水池中，采用深井水24 h循環(huán)供水，并進(jìn)行曝氣處理，溶解氧濃度保持在7 mg/L以上（Branco et al，2013）。試驗(yàn)期間暫養(yǎng)池水溫（16±2）℃，pH值（7.3±2.0），光照為室內(nèi)自然光。試驗(yàn)前2 d停止喂食（Mallen-Cooper，1994）。

試驗(yàn)分為2組，設(shè)置試驗(yàn)組Ⅰ篩選出較優(yōu)的魚道結(jié)構(gòu)，分別在試驗(yàn)池室第3、5、7、9號布置體型Ⅰ～Ⅳ，根據(jù)各級池室通過時(shí)間，選擇相對較優(yōu)的池室結(jié)構(gòu)布置體型。該試驗(yàn)未設(shè)置對照試驗(yàn)，作為初步篩選較優(yōu)池室結(jié)構(gòu)方案，有利于縮小試驗(yàn)周期，避免試驗(yàn)期間水溫等環(huán)境因素變幅過大對過魚試驗(yàn)結(jié)果的影響。試驗(yàn)組Ⅱ設(shè)置對照試驗(yàn)，以篩選的較優(yōu)結(jié)構(gòu)作為試驗(yàn)組，標(biāo)準(zhǔn)魚道作為對照組，進(jìn)一步分析樁柱結(jié)構(gòu)對魚類上溯行為的影響，試驗(yàn)布置見圖2。相比于單一試驗(yàn)，對照試驗(yàn)?zāi)鼙ＷC在相同條件下進(jìn)行過魚試驗(yàn)，降低水溫、光照和聲音等無關(guān)因素對試驗(yàn)結(jié)果的影響（郭子琪等，2021；黃婕等；2021；譚紅林等，2021）。

試驗(yàn)前調(diào)整魚道運(yùn)行水深為0.3 m，測量并記錄試驗(yàn)魚道中的水溫、含氧量、pH值等相關(guān)參數(shù)，確保試驗(yàn)中循環(huán)水含氧量大于6 mg/L，pH值在（7.75±2.5）（Branco et al，2013）。每天隨機(jī)選取20尾魚作為試驗(yàn)用魚，每尾魚僅用于1次試驗(yàn)，試驗(yàn)后將魚放置于另一水池暫養(yǎng)。魚道流態(tài)穩(wěn)定后開始放魚，關(guān)閉攔魚柵，待試驗(yàn)魚適應(yīng)15 min后，撤去攔魚柵。自試驗(yàn)魚進(jìn)入11號池室時(shí)開始計(jì)時(shí)，待其游出1號池室時(shí)停止計(jì)時(shí)，記錄上溯時(shí)間和上溯視頻。若試驗(yàn)魚超過20 min未通過魚道，則記為上溯失敗。本試驗(yàn)共放魚106尾，其中順利通過11級池室67尾，作為數(shù)據(jù)分析對象，避免各體型池室通過率偏差對統(tǒng)計(jì)數(shù)據(jù)的影響。試驗(yàn)組Ⅱ中，將2組魚同時(shí)放入下游放魚池并開始試驗(yàn)，試驗(yàn)過程中分別記錄2組魚道的通過時(shí)間和上溯視頻。試驗(yàn)完成后，隨機(jī)選取另2尾魚繼續(xù)下一組對照試驗(yàn)。

1.2? ?數(shù)值模擬

1.2.1? 數(shù)學(xué)模型? ?本研究使用Fluent軟件對4種樁柱結(jié)構(gòu)進(jìn)行數(shù)值模擬。采用標(biāo)準(zhǔn)k-ε湍流模型進(jìn)行湍流計(jì)算，VOF模型追蹤魚道自由水面，結(jié)合壓力-速度耦合SIMPLE算法進(jìn)行求解計(jì)算（Li et al，2021）。計(jì)算域包括11級池室，如圖3所示。采用結(jié)構(gòu)性六面體網(wǎng)格，魚道池室網(wǎng)格尺寸0.05 m，豎縫及樁柱處網(wǎng)格尺寸0.04 m，垂向網(wǎng)格0.025 m，計(jì)算區(qū)域高度1 m。魚道進(jìn)口及出口邊界條件設(shè)置為壓力進(jìn)口，給定水深0.3 m；頂部邊界條件設(shè)置為壓力進(jìn)口，相對壓強(qiáng)設(shè)置為0；底面及邊壁設(shè)置為無滑移壁面。迭代計(jì)算時(shí)間步長為0.01 s，每步迭代最大次數(shù)設(shè)置為30（Anastasios et al，2018; Fuentes-Pérez et al，2018）。

1.2.2? ?模型驗(yàn)證? ?數(shù)值模擬結(jié)果與物理模型試驗(yàn)進(jìn)行對比，以驗(yàn)證模擬結(jié)果的可靠性。魚道流場為典型的二維分布（Wu et al，1999），為此選取0.5H0橫截面作為數(shù)據(jù)提取斷面，以斷面Ⅰ（X=15 cm）、斷面Ⅱ（X=40 cm）和斷面Ⅲ（X=60 cm）作為驗(yàn)證斷面，對體型Ⅱ的模型進(jìn)行驗(yàn)證。計(jì)算結(jié)果表明，各體型驗(yàn)證斷面數(shù)值模擬與實(shí)測結(jié)果流速分布吻合程度較高，流速誤差在0.03 m/s以內(nèi)，數(shù)值模擬結(jié)果與物理模型試驗(yàn)結(jié)果基本一致（圖4）。

2? ?結(jié)果與分析

2.1? ?池室水力特性

2.1.1? ?回流區(qū)尺度? ?按數(shù)值計(jì)算的流線圖劃分主流區(qū)及回流區(qū)，結(jié)果見圖5。將回流區(qū)的分布圖導(dǎo)入CAD，計(jì)算回流區(qū)面積a，與魚道池室面積A相比進(jìn)行無量綱化處理，得出回流區(qū)在魚道中的占比a/A。統(tǒng)計(jì)5種結(jié)構(gòu)的回流區(qū)相關(guān)數(shù)據(jù)如表1。

體型Ⅰ中，主流經(jīng)過樁柱產(chǎn)生向池室左側(cè)分支（順?biāo)鞣较蜃髠?cè)），主流范圍更加廣泛，增加可能的魚類上溯路線，同時(shí)左側(cè)大回流區(qū)被削減。體型Ⅱ中，主流由于中間樁柱作用主流范圍有所擴(kuò)大，偏向樁柱左側(cè)。體型Ⅲ、Ⅳ在魚道回流區(qū)增設(shè)樁柱，其中體型Ⅳ由于樁柱靠近池室左側(cè)，減小了左側(cè)回流區(qū)對主流的壓縮，回流區(qū)范圍更加廣泛。

相較于標(biāo)準(zhǔn)魚道，布置樁柱結(jié)構(gòu)的池室回流區(qū)總面積占比都有一定程度減小，其中體型Ⅱ減小幅度最大，達(dá)49%。由于回流區(qū)面積減小，對主流的約束作用減弱，其分布范圍更加廣泛，體型Ⅰ、Ⅱ、Ⅳ主流均有偏向左側(cè)墻壁的趨勢。4種結(jié)構(gòu)中，回流區(qū)的變化主要集中在左側(cè)，體型Ⅰ中左側(cè)大回流區(qū)被分割為4個(gè)小回流區(qū)，體型Ⅲ、Ⅳ回流區(qū)也有所增加。體型Ⅱ中大回流區(qū)縮減為位于上游隔板側(cè)的小回流區(qū)，減小為原來的26.1%，這也是其回流區(qū)占比減小明顯的主要原因。

2.1.2? ?池室流場水力參數(shù)分布? ?使用Tecplot提取各體型池室數(shù)值模擬結(jié)果，選擇Z=0.5H0平面作為數(shù)據(jù)分析斷面，分別做流速、紊動能和總水力應(yīng)變云圖，觀察不同體型池室各水力參數(shù)分布差異。

流速是影響魚類上溯的主要水力要素，也是魚道設(shè)計(jì)的重要判斷標(biāo)準(zhǔn)（Hou et al，2019）。5種體型中，高流速區(qū)主要在主流附近區(qū)域，低流速區(qū)主要位于兩側(cè)回流區(qū)。主流區(qū)流速為 0.20～0.38 m/s，回流區(qū)流速較小，僅0～0.1 m/s，其最大流速均位于豎縫靠近隔板側(cè)。除在體型Ⅰ中池室流速有明顯增加外，其他3種體型魚標(biāo)準(zhǔn)魚道流速分布差異不大，各體型魚道的流速場分布如圖6-a。

紊動能TKE反應(yīng)出水流流速脈動振幅，若在實(shí)際流場中魚類所經(jīng)過的區(qū)域水流紊動過大，則會導(dǎo)致魚類上溯體能消耗加快，進(jìn)一步影響其游泳能力（Enders et al，2003）。5種體型池室的紊動能如圖6-b所示。紊動能較高區(qū)域都位于豎縫處，回流區(qū)紊動能低于0.006 m2/s2，豎縫及部分主流區(qū)域則在0.010～0.025 m2/s2。

對比5種體型池室的紊動能分布發(fā)現(xiàn)，體型Ⅱ最大紊動能0.018 m2/s2低于標(biāo)準(zhǔn)魚道0.022 m2/s2，減小幅度為18%；體型Ⅰ豎縫處最大紊動能小于標(biāo)準(zhǔn)魚道。5種魚道中標(biāo)準(zhǔn)魚道豎縫處紊動能為0.018～0.022 m2/s2的區(qū)域面積最大，主流區(qū)紊動能在0.006～0.008 m2/s2的面積更加廣泛。體型Ⅲ、體型Ⅳ豎縫處紊動能有明顯的增加，最大紊動能為0.024 m2/s2?？傮w來說，在主流區(qū)增設(shè)樁柱的方案可以有效減小魚道內(nèi)的紊動能，回流區(qū)增設(shè)樁柱結(jié)構(gòu)會增加豎縫處的紊動能。

總水力應(yīng)變用以反映流場區(qū)域內(nèi)控制體變形情況，可以綜合反映各方向流速變化的梯度（Nestler et al，2008；齊亮等，2012）。本研究5種魚道池室中，總水力梯度相對較高的區(qū)域均位于主流區(qū)域，且各體型豎縫處總水力應(yīng)變明顯高于池室內(nèi)部；其中，體型Ⅰ豎縫處水力梯度最高，達(dá)到18/s；體型Ⅲ、Ⅳ豎縫處水力梯度分布相似，導(dǎo)板端部較小部分存在水力梯度達(dá)到18/s；體型Ⅱ與標(biāo)準(zhǔn)魚道豎縫處水力梯度相似，處于相對較低水平。各魚道池室體型總水力應(yīng)變?nèi)鐖D6-c所示。

2.2? ?過魚效果

2.2.1? ?通過時(shí)間? ?試驗(yàn)組Ⅰ過魚試驗(yàn)于9月上、中旬進(jìn)行，平均水溫15.4℃，溶解氧為7.0 mg/L以上。試驗(yàn)共放魚106尾，其中順利通過11級池室共67尾，統(tǒng)計(jì)其通過時(shí)間中位數(shù)及平均數(shù)如圖7-a所示。體型Ⅱ通過時(shí)間平均數(shù)相比其他體型時(shí)間更短，為 14 s，同時(shí)體型Ⅱ通過時(shí)間中位數(shù)與其他體型相比也有優(yōu)勢，通過時(shí)間為 9 s。為此，以齊口裂腹魚為過魚對象時(shí)推薦體型Ⅱ作為初步篩選的池室結(jié)構(gòu)。試驗(yàn)組Ⅱ過魚試驗(yàn)于10月中旬至11月中旬進(jìn)行，平均水溫15.6℃，溶解氧為7.0 mg/L以上。體型Ⅱ與標(biāo)準(zhǔn)魚道通過時(shí)間中位數(shù)與平均數(shù)差異并不大。試驗(yàn)組Ⅱ各體型魚道通過時(shí)間如圖7-b所示。

2.2.2? ?通過率? ?通過率是判斷魚道過魚效果的重要判斷指標(biāo)。試驗(yàn)組Ⅰ中未詳細(xì)統(tǒng)計(jì)各體型魚道的通過率，參考樣本為順利通過11級池室的魚類，各體型魚道通過率一致，以避免各池室布置位置差異對過魚效果的影響。試驗(yàn)組Ⅱ?qū)φ赵囼?yàn)共放魚168尾，標(biāo)準(zhǔn)魚道成功上溯121尾，通過率PrⅠ=72%；體型Ⅱ成功上溯115尾，通過率PrⅡ=68%。試驗(yàn)組Ⅱ中標(biāo)準(zhǔn)魚道中試驗(yàn)魚通過率略高于體型Ⅱ。

2.2.3? ?生存曲線? ?采用生存分析方法，分析試驗(yàn)組Ⅱ中不同池室體型對齊口裂腹魚通過效果的影響。生存分析不僅綜合考慮了通過時(shí)間與通過率指標(biāo)，還考慮到隨通過時(shí)間遞增通過率的變化。隨后利用時(shí)序檢驗(yàn)（log-rank test）的方法比較不同生存曲線分布區(qū)別是否顯著。將試驗(yàn)組數(shù)據(jù)導(dǎo)入統(tǒng)計(jì)分析軟件GraphPad Prism 7對生存曲線進(jìn)行比較。結(jié)果表明，試驗(yàn)組Ⅱ中標(biāo)準(zhǔn)魚道和體型Ⅱ過魚效果無顯著差異（P=0.361>0.05），各試驗(yàn)組生存曲線如圖7-c所示。

2.3? ?上溯行為

2.3.1? ?上溯軌跡熱點(diǎn)? ?為進(jìn)一步深入分析2種結(jié)構(gòu)對魚類上溯的影響，對其上溯的軌跡進(jìn)行統(tǒng)計(jì)分析。由攝像機(jī)記錄魚類在池室中的上溯視頻，采用LoggerPro對魚類上溯過程進(jìn)行打點(diǎn)，通過多條上溯軌跡疊加繪制熱點(diǎn)圖8。

試驗(yàn)組Ⅱ中，體型Ⅱ魚類上溯軌跡主要集中在樁柱上部，之后沿隔板側(cè)通過豎縫，上溯軌跡明顯，部分魚類沿樁柱下端邊墻貼壁上溯。標(biāo)準(zhǔn)魚道中大部分魚類選擇沿下部分邊墻上溯，在導(dǎo)板處轉(zhuǎn)向通過豎縫，部分魚類靠近池室右側(cè)沿隔板側(cè)通過豎縫，結(jié)果與Shi等（2022）的研究結(jié)論一致。其中樁柱結(jié)構(gòu)中魚道上溯軌跡分布更加豐富，而標(biāo)準(zhǔn)魚道中魚類上溯多集中于導(dǎo)板一側(cè)，位于左側(cè)大回流區(qū)少有魚類通過，樁柱結(jié)構(gòu)豐富了魚類上溯路徑的選擇，有利于提高池室空間利用率。

觀察體型Ⅱ魚類主要上溯路徑上的行為發(fā)現(xiàn)，魚類多沿隔板側(cè)通過豎縫，通過豎縫后存在明顯的調(diào)整上溯方向行為，大多數(shù)魚類調(diào)整方向朝向樁柱側(cè)，之后沿該方向移動至樁柱左側(cè)。在樁柱附近魚類又一次調(diào)整方向垂向隔板側(cè)，沿該方向行一段距離后，魚類轉(zhuǎn)向靠近隔板側(cè)，貼近隔板行進(jìn)至豎縫處后調(diào)整方向沿隔板背水面通過豎縫，完成單級池室的上溯，魚類上溯過程如圖9。

2.3.2? ?水力因子與上溯軌跡耦合? ?由上溯軌跡熱點(diǎn)圖8可以發(fā)現(xiàn)主要上溯路線差異明顯，提取魚類的主要上溯路線與水力因子、回流區(qū)分布耦合，分析魚類在上溯軌跡上差異的成因，結(jié)果如圖10。剔除軌跡線中規(guī)律不明顯的其他軌跡，得出體型Ⅱ主要軌跡線占比為52%，次要軌跡占比30%，標(biāo)準(zhǔn)魚道中主要軌跡線占比78%，次要軌跡占比12%。

魚類上溯過程中都偏向由回流區(qū)域外側(cè)通過，避免了進(jìn)入回流區(qū)中中心，同時(shí)偏好于較小回流區(qū)側(cè)通過，可能是魚類偏好于流向單一的通道上溯，回流區(qū)中流向復(fù)雜不利于魚類辨明上溯方向。從流速云圖的耦合中可以發(fā)現(xiàn)，魚類偏好避免通過豎縫下游段的高流速區(qū)域，行進(jìn)路線中流速大多處于較低水平。

兩者軌跡線的差異主要集中在魚類通過豎縫后上溯方向的選擇，標(biāo)準(zhǔn)魚道中大部分魚類通過豎縫后由于慣性作用繼續(xù)沿當(dāng)前方向行進(jìn)至邊墻，而體型Ⅱ中魚類大多選擇沿隔板側(cè)直行，之后調(diào)整方向朝向樁柱側(cè)行進(jìn)，在隔板側(cè)存在主動選擇上溯方向的行為，可能是魚類偏向高紊動能、高水力應(yīng)變區(qū)域上溯。標(biāo)準(zhǔn)魚道中，豎縫上游端突出存在較高水力應(yīng)變與高紊動能區(qū)域，其與魚類的主要上溯路線存在較高的吻合度；而體型Ⅱ中高水力應(yīng)變高紊動能區(qū)集中在隔板側(cè)，且其范圍較窄導(dǎo)向性并強(qiáng)烈，為此魚類的主要上溯路線偏向隔板側(cè)，且主要軌跡占比沒有標(biāo)準(zhǔn)魚道高。在體型Ⅱ中樁柱左側(cè)，豎縫下游段存在較高的水力應(yīng)變與紊動能區(qū)域，魚類的主要上溯軌跡中在該區(qū)域都有方向調(diào)整的傾向，表明魚類偏好于高水力應(yīng)變高紊動能區(qū)域進(jìn)行方向調(diào)整。

2.3.3? ?水力因子上溯偏好范圍? ?對魚類軌跡分布與池室水力因子分布進(jìn)行定量分析，將池室按照熱點(diǎn)圖中網(wǎng)格劃分，提取每個(gè)網(wǎng)格內(nèi)水力參數(shù)的平均值，結(jié)合魚類在該網(wǎng)格分布數(shù)量繪制直方圖11。對魚類在各網(wǎng)格中的數(shù)量分布與水力因子進(jìn)行相關(guān)性分析發(fā)現(xiàn)，無論標(biāo)準(zhǔn)魚道還是體型Ⅱ魚道，魚類數(shù)量的分布與流速、紊動能、水力應(yīng)變具有顯著的相關(guān)性，大部分P<0.01（表2）。雖然2種體型過魚軌跡存在明顯的差異，但魚類偏好的水力因子范圍存在一致性，都偏向于較低水力因子區(qū)域上溯，裂腹魚在上溯過程中偏好水力應(yīng)變0.2～3.0/s，偏好流速為0.01～0.09 m/s，偏好紊動能為0.001～0.008 m2/s2。本研究與譚均軍等（2017）研究得出草魚、鰱等幼魚偏好水力因子范圍有所不同，可能是不同種類不同體長魚類偏好水力因子范圍存在差異。與本研究團(tuán)隊(duì)在裂腹魚的其他試驗(yàn)研究中結(jié)論一致，都偏好在低流速、低紊動能區(qū)域上溯（Li et al，2021） 。

本試驗(yàn)表明，體型Ⅱ魚道與標(biāo)準(zhǔn)魚道從通過時(shí)間及通過率的角度看來，兩者效果相當(dāng)。體型Ⅱ魚道通過在主流位置增加樁柱結(jié)構(gòu)，形成了多股水流，有效減小了回流區(qū)的占比，提高了魚道池室空間利用率，為魚類的上溯路徑提供了更多選擇。鑒于本文僅開展了齊口裂腹魚的相關(guān)試驗(yàn)研究，尚未能證實(shí)體型Ⅱ魚道的過魚優(yōu)勢，后續(xù)擬進(jìn)一步開展其他魚種的過魚試驗(yàn)，以進(jìn)一步論證體型Ⅱ魚道的合理性。

3? ?結(jié)論

在豎縫式魚道池室內(nèi)部增設(shè)樁柱結(jié)構(gòu)，并采用數(shù)值模擬和放魚試驗(yàn)方法驗(yàn)證4種樁柱結(jié)構(gòu)的可靠性，評價(jià)樁柱結(jié)構(gòu)實(shí)際過魚效果，得出如下結(jié)論：

（1）體型Ⅱ魚道在主流中間處布置樁柱結(jié)構(gòu)可以有效的減小回流區(qū)的占比，有利于減小池室內(nèi)紊動能，同時(shí)相較于本文其他結(jié)構(gòu)，齊口裂腹魚的通過時(shí)間較短。

（2）體型Ⅱ魚道與標(biāo)準(zhǔn)魚道過魚效果相當(dāng)，但具有更高的空間利用率，形成了具有流速梯度的多股水流，為魚類的上溯路徑提供了更多選擇。

（3）齊口裂腹魚在上溯過程中偏好于回流區(qū)外側(cè)通過，都選擇避免豎縫下游段的高流速區(qū)域，偏好于高水力應(yīng)變高紊動能區(qū)域進(jìn)行方向調(diào)整。

（4）魚類上溯軌跡偏好與流速、紊動能、水力應(yīng)變顯著相關(guān)（P≤0.01），偏好流速0.01～0.09 m/s，偏好紊動能0.001～0.008 m2/s2，偏好水力應(yīng)變0.2～3.0/s。

參考文獻(xiàn)

陳凱麒，常仲農(nóng)，曹曉紅，等，2012. 我國魚道的建設(shè)現(xiàn)狀與展望[J]. 水利學(xué)報(bào)， 43（2）：182-188.

董哲仁，2006. 筑壩河流的生態(tài)補(bǔ)償[J]. 中國工程科學(xué)， 8（1）：5-10．

郭子琪，李廣寧，郄志紅，等，2021. 水溫對豎縫式魚道中齊口裂腹魚上溯行為影響試驗(yàn)研究[J]. 中國水利水電科學(xué)研究院學(xué)報(bào)， 19（2）：255-261.

黃婕，林晨宇，石小濤，等，2021. 基于齊口裂腹魚趨光特性的涵洞式魚道光環(huán)境優(yōu)化[J]. 生態(tài)學(xué)雜志， 40（7）：2155-2163.

劉志雄，周赤，黃明海， 2010. 魚道應(yīng)用現(xiàn)狀和研究進(jìn)展[J]. 長江科學(xué)院院報(bào)， 27（4）：28-31，35.

齊亮，楊宇，王悅，等，2012. 魚類對水動力環(huán)境變化的行為響應(yīng)特征[J]. 河海大學(xué)學(xué)報(bào)（自然科學(xué)版）， 40（4）：438-445.

譚紅林，譚均軍，石小濤，等，2021. 魚道進(jìn)口誘魚技術(shù)研究進(jìn)展[J]. 生態(tài)學(xué)雜志， 40（4）：1198-1209.

譚均軍，高柱，戴會超，等，2017. 豎縫式魚道水力特性與魚類運(yùn)動特性相關(guān)性分析[J]. 水利學(xué)報(bào)， 48（8）：924-932，944.

王猛，史德亮，陳輝，等，2015. 豎縫式魚道池室結(jié)構(gòu)變化對水力特性的影響分析[J]. 長江科學(xué)院院報(bào)， 32（1）：79-83.

徐體兵，孫雙科， 2009. 豎縫式魚道水流結(jié)構(gòu)的數(shù)值模擬[J]. 水利學(xué)報(bào)， 40（11）：1386-1391．

張國強(qiáng)，孫雙科，2012. 豎縫寬度對豎縫式魚道水流結(jié)構(gòu)的影響[J]. 水力發(fā)電學(xué)報(bào)， 31（1）：151-156．

Anastasios I S， Mitsopoulos G， Rutschmann P， et al， 2018. Verification of a 3D CFD model for vertical slot fish-passes[J]. Environmental Fluid Mechanics， 18（6）：1435-1461.

Aurélien B， Gérard P， Damien C， et al， 2019. Experimental-based methodology to improve the design of vertical slot fishways[J]. J Hydraul Eng， 145（9）：04019031.

Branco P， Santos J M， Katopodis C， et al， 2013. Effect of flow regime hydraulics on passage performance of Iberian chub （Squalius pyrenaicus Günther， 1868） in an experimental pool-and-weir fishway[J]. Hydrobiologia， 714（1）：145-154.

Cai L， Katopodis C， Johnson D， et al， 2018. Case study： targeting species and applying swimming performance data to fish lift design for the Huangdeng Dam on the upper Mekong River[J]. Ecological Engineering， 122：32-38.

Calluaud D， Pineau G， Texier A， et al， 2014. Modification of vertical slot fishway flow with a supplementary cylinder[J]. Journal of Hydraulic Research， 52（5）：614-629.

Eva C Enders， Daniel Boisclair， André G Roy， 2003. The effect of turbulence on the cost of swimming for juvenile Atlantic salmon （Salmo salar）[J]. Canadian Journal of Fisheries & Aquatic Sciences， 60（4）：1149-1160.

Fuentes-Pérez J F， Silva A T， Tuhtan J A， et al， 2018. 3D modelling of non-uniform and turbulent flow in vertical slot fishways[J]. Environmental Modelling & Software， 99：156-169.

Guiny E， Ervine D A， Armstrong J D， 2005. Hydraulic and biological aspects of fish passes for Atlantic salmon[J]. J Hydraulic Eng， 131：542-553.

Hou Y Q， Yang Z， An R， et al， 2019. Water flow and substrate preferences of Schizothorax wangchiachii （Fang， 1936）[J]. Ecological Engineering， 138（9）：1-7.

Kynard B， Horgan M， Pugh D， et al， 2008. Using juvenile sturgeons as a substitute for adults： a new way to develop fish passage for large fish[C]//American Fisheries Society Symposium.

Li G N， Sun S K， Zhang C， et al， 2019. Evaluation of flow patterns in vertical slot fishways with different slot positions based on a comparison passage experiment for juvenile grass carp[J]. Ecological Engineering， 133：148-159.

Li G N， Sun S K， Liu H T， et al， 2021. Schizothorax prenanti swimming behavior in response to different flow patterns in vertical slot fishways with different slot positions[J]. Science of the Total Environment， 754：142142.

Mallen-Cooper M， 1994. Swimming ability of adult golden perch， Macquaria ambigua （Percichthyidae）， and adult silver perch， Bidyanus bidyanus （Teraponidae）， in an experimental vertical-slot fishway[J]. Marine and Freshwater Research， 45（2）：191-198.

Mao X， 2018. Review of fishway research in China[J]. Ecological Engineering， 115：91-95.

Nestler J M ， Goodwin R A， Smith D L， et al， 2008. Optimum fish passage and guidance designs are based in the hydrogeomorphology of natural rivers[J]. River Research and Applications， 24（2）：148-168.

Rajaratnam N， Gary Van der Vinne， Katopodis C， 1986. Hydraulics of vertical slot fishways[J]. Journal of Hydraulic Engineering， 112（10）：909-927．

Rajaratnam N， Katopodis C， Solanki S， 2010. New designs for vertical slot fishways[J]. Journal of Civil Engineering， 36（1）：36-46.

Shi K， Li G N， Liu S T， et al， 2022. Experimental study on the passage behavior of juvenile Schizothorax prenanti by configuring local colors in the vertical slot fishways[J]. Science of The Total Environment， 843：156989.

Tarrade L， Texier A， David L， et al， 2008. Topologies and measurements of turbulent flow in vertical slot fishways[J]. Hydrobiologia， 609（1）：177-188.

Wang R W， David L， Larinier M， 2010. Contribution of experimental fluid mechanics to the design of vertical slot fish passes[J]. Knowledge and Management of Aquatic Ecosystems， 396（2）：202-221.

Wu S， Rajaratnam N， Katopodis C， 1999. Stricture of flow in vertical slot fishway[J]. Journal of Hydraulic Engineering， 25（4）：351-360.

（責(zé)任編輯? ?萬月華）

收稿日期：2022-07-14? ? ? 修回日期：2022-11-14

基金項(xiàng)目：國家自然科學(xué)基金（51679261）；流域水循環(huán)模擬與調(diào)控國家重點(diǎn)試驗(yàn)室自主研究課題（SKL2020TS04）。

作者簡介：柳松濤，1997年生，男，碩士研究生，研究方向?yàn)轸~道水力學(xué)。E-mail：202008150021032@ctgu.edu.cn

通信作者：李廣寧，1993年生，男，高級工程師，主要從事水力學(xué)及河流動力學(xué)研究。E-mail：lgnchina@163.com