呂澤硯，朱文斌，唐偉堯，管青龍，劉莉莉，張紋瑄，傅鑫洋

(浙江海洋大學水產(chǎn)學院，浙江舟山 316022)

人工魚礁是人工投放于海底的漁業(yè)工程設施，可以有效地吸引魚群，為魚類提供適宜棲居、繁殖及躲避敵害的場所，從而有效地增殖魚類，養(yǎng)護漁業(yè)資源，提高漁業(yè)資源修復效率，改善水域生態(tài)環(huán)境[1-5]。許多沿海國家建立專門的組織，系統(tǒng)地實施沿海漁業(yè)修復工程項目，并加大人工魚礁的科學研究和開發(fā)投入。

由于人工魚礁的投放可以改變海域局部流場，所產(chǎn)生的上升流及渦流的影響可以改變沉積物中營養(yǎng)物質在水域中的分布[6-10]。因此，在人工魚礁流場效應研究過程中，掌握其周圍流態(tài)的變化尤為關鍵。目前，關于人工魚礁的流場效應研究主要分為數(shù)值模擬與模型試驗，數(shù)值模擬主要集中在對不同結構特性和不同工況下人工魚礁所在的海域流場分布進行分析，而模型試驗主要分為水槽試驗[11-13]、PIV 模型試驗等。劉洪生等[14]和李珺等[15]借助水槽實驗，驗證了人工魚礁流場效應數(shù)值模擬研究的可行性；隨著數(shù)值模擬技術在人工魚礁流場效應研究的深入，在礁體開孔情況[16-19]、礁體形狀[20-22]、礁體布設間距[23]、投放個數(shù)與組合方式[24-26]等諸要素對人工礁體流場效應的影響研究方面，取得了一定進展。其中單個礁體的研究均表明礁體形狀的變化對繞流情況有顯著影響。

人工魚礁的設計需要同時考慮流場效應以及造礁材料的成本問題。魚礁按外觀可分為正方體形、多面體形、錐形、圓筒形、半圓形等多種形狀，對錐形魚礁而言，傾斜角是影響魚礁周圍流場的因素之一，不同的傾斜角會對上升流效應及流場變化有較大影響。當前，國內關于不同傾斜角情況下的魚礁流場效應研究未見報道。因此，本研究以四角錐臺形魚礁為例，在保持魚礁等底等高的情況下，預設5 種不同的斜面傾角，建立計算流體動力學(CFD)模型，通過有限元分析軟件(ANSYS)對其周圍流態(tài)進行仿真模擬，揭示不同傾斜角變化對人工魚礁流場效應的影響，以期為四角錐臺形人工魚礁的設計提供參考依據(jù)。

1 材料與方法

1.1 計算流體動力學模型

該模擬計算的流體范圍為長30 m、寬16 m、高16 m 的長方體區(qū)域，來流方向水平從左至右(圖1)。計算所用魚礁模型為正四角錐臺形人工魚礁，底面正方形邊長為3 m，高度為3 m，在等底等高的情況下，通過改變頂面正方形邊長以改變斜面傾角,由此得到5 個魚礁模型(表1)。表1 中傾斜角在文章中指的是礁體側面與底面所成的夾角，即斜面傾角，邊長比(α)為正四角錐臺形魚礁上表面正方形邊長(D)與下表面正方形邊長之比(L)，即α=D/L。等距改變頂面邊長得到5 個邊長比為0.2、0.4、0.6、0.8、1，對應5 種情況的傾斜角分別為68.2°、73.3°、78.7°、84.3°、90°。由于本研究只考慮傾斜角作為單一變化因素，考慮到開孔狀態(tài)對流場的干擾，故模型做無孔設計。圖2 和圖3 分別為傾斜角73.3°(即上底邊1.2 m，上下正方形邊長比為0.4)的人工魚礁的正視圖與俯視圖。魚礁底部所置地面設為Z=0 m 的X-Y 平面，坐標原點在魚礁底部中心點處。魚礁迎流面底邊離計算域入口8.5 m(礁體高度的2.8 倍)，魚礁背流面底邊離計算域出口18.5 m(礁體高度的6.2 倍)，礁體上表面離計算域頂部13 m(礁體高度的4.3 倍)，礁體中心點到計算域兩側的距離為8 m(礁體高度的2.6 倍)。

圖1 流體計算域Fig.1 Fluid computational domain

圖2 礁體俯視圖Fig.2 Vertical view of artificial reef

圖3 魚礁正視圖Fig.3 Front view of artificial reef

表1 四角錐臺形人工魚礁的結構參數(shù)表Tab.1 Structure parameters of artificial reef

1.2 控制方程

該研究為探究魚礁周圍的流場效應的時間平均流特性，假定模型中流體為粘性不可壓縮流體。水流控制方程采用定常、不可壓縮流動下的連續(xù)性方程和Navier-Stokes 方程(N-S 方程)[27]：

連續(xù)方程：

動量方程：

式中：u、v、w 分別為x、y、z 方向上的速度分量，m·s-1；v 為流體的運動學黏性系數(shù)，m2·s-1；t 為時間，s；ρ 為流體密度，kg·m-3；p 為壓強，Pa；fx、fy、fz分別為x、y、z 方向上的單位質量力，m·s-2。

1.3 計算區(qū)域的網(wǎng)格劃分

在計算區(qū)域的網(wǎng)格劃分上，應用計算軟件的自適應網(wǎng)格劃分方式，劃分成四面體非結構化網(wǎng)格。礁體附近區(qū)域最大層數(shù)設為8，增長率設為1.2，使數(shù)值模擬時可捕捉礁體的幾何細節(jié)；在其他區(qū)域網(wǎng)格最大尺寸設為0.6 m。

1.4 邊界條件

邊界條件的設置，進口條件均設置為速度進口。通常情況下，近岸海水流速為1 m·s-1以內，以杭州灣為例，潮流平均速度約為0.2 m·s-1，而在漲落潮時平均速度達0.8 m·s-1[28]。故設置4 種流速分別為0.2、0.4、0.6 以及0.8 m·s-1。出口條件設置為壓力出口。計算域兩側面以及頂面設置為對稱邊界條件。所有固體均采用無滑移條件。

1.5 數(shù)值算法

計算在ANSYS fluent 平臺上完成。壓力速度耦合用SIMPLEC 算法，壓力項處理則用standard 差分格式[35]，各方程的空間離散均采用QUICK 格式[36]，計算殘差值取10-5，計算的最大迭代步數(shù)在多次計算過程中平均500 步以后轉為平穩(wěn)，所以取500 步。

2 結果與分析

上升流會把深水區(qū)大量的海水營養(yǎng)鹽(磷酸鹽、硝酸鹽等)帶到表層，為魚類生長帶來豐富的餌料，因此，上升流顯著的海區(qū)經(jīng)常會形成著名的漁場。此外，水流流經(jīng)礁石時，在礁石后面會產(chǎn)生一個漩渦區(qū)，即背渦流，渦流速度動態(tài)穩(wěn)定。許多研究表明，魚類常常聚集在珊瑚礁周圍，因為其背流區(qū)內的不同水層更容易相互混合，其營養(yǎng)鹽以及溶解鹽也會高于正常水域，可以為魚類提供有利的生長環(huán)境[29-30]，且背渦流的影響效應與背渦區(qū)內部流態(tài)的范圍大小、漩渦強度及大小等因素相關。因此，上升流以及背渦流的整體規(guī)模、渦量強度、速度方向等參數(shù)是衡量人工魚礁效益的重要指標。

2.1 不同速度下的上升流與背渦流

文章以表1 中模型c(邊長比0.6，傾斜角78.7°)為計算模型，對不同速度下的礁體周圍流態(tài)進行了數(shù)值模擬。圖4、圖5、圖6 分別為不同流速(0.2、0.4、0.6、0.8 m·s-1)下礁體周圍水流在y=0 m 平面(側面圖)的速度矢量分布、渦量分布以及z=0.5 m 平面(俯視圖)的速度矢量分布。為對各個速度下的情況進行合理比較，將各種情況下的比色棒所代表的速度統(tǒng)一大小區(qū)間。由于4 種來流速度是逐漸增大的，因此圖4 和圖6 中的4 幅圖在主色調上有明顯差異。

從速度矢量分布來看，不同來流速度，相同傾斜角的情況下，魚礁周圍流場的變化趨勢相近，具體表現(xiàn)為流體變化的速度方向和影響范圍相似，但是流速大小有明顯差異(圖4、圖6)。另外，4 種速度下，渦量最大處均出現(xiàn)于魚礁迎流面與上表面交匯處，還可以清晰地看出渦量較大處的規(guī)模隨來流速度的增大同步增長，且其分布于后方延長處(圖5)。

圖4 4 種不同來流速度下y=0 平面上的速度矢量分布圖(α=0.6)Fig.4 Velocity vector distribution on the plane y=0 at 4 different inflow velocities (α=0.6)

圖5 4 種不同來流速度下y=0 平面上的渦量分布云圖(α=0.6)Fig.5 Vorticity distribution cloud maps on the y=0 plane at four different inflow velocities (α=0.6)

圖6 4 種不同來流速度下z=0.5 m 平面上的速度矢量分布圖(α=0.6)Fig.6 Velocity vector distribution on the plane z=0.5 m at 4 different inflow velocities (α=0.6)

2.2 不同傾斜角下的上升流、背渦流、側向流

根據(jù)礁體傾斜角的不同，迎流面的上后方會形成規(guī)模大小不等的上升流，由圖7 可知，傾斜角大的魚礁繞流的影響高度以及寬度都大于傾斜角小的繞流，且上升流規(guī)模大小隨著傾斜角的增大而增大，這是由于隨著迎流面積的增大，阻水能力增強所導致的。

從側視角度(y=0 m 平面)來看，魚礁后方的背渦流是一個順時針漩渦，且隨著魚礁傾斜角的增大，它的基本形態(tài)不變，規(guī)模略微增大，渦心則向上偏移，且水平高度始終低于魚礁頂面高度(圖7)。

圖7 5 種不同傾斜角下y=0 m 平面上的速度矢量分布圖(v=0.4 m·s-1)Fig.7 Velocity vector distribution on the plane y=0 m at 5 different inclination angles (v=0.4 m·s-1)

魚礁周圍水體不僅會形成上升流，還會形成側向流，且該平面內的速度最大處也出現(xiàn)在側向流內，側向流規(guī)模隨著傾斜角的增大而增大，如圖8 所示。在礁體迎流面兩側后方均會形成2 個方向相反的漩渦，這兩個漩渦會形成1 個大漩渦區(qū)(即背渦區(qū))，且在相同來流速度下，傾斜角越大，其背渦區(qū)內的兩個相反漩渦均更明顯。背渦區(qū)的寬度以及長度均隨傾斜角的增大而增大，且增幅較大。另外，當傾斜角越大，其后方的流態(tài)更為雜亂，這可能是由于傾斜角大的魚礁在z=0.5 m 平面上方的體積更大，導致其對與z=0.5 m平面魚礁后方的流態(tài)影響更為明顯。

圖8 5 種不同傾斜角下z=0.5 m 平面上的速度矢量分布圖(v=0.4 m·s-1)Fig.8 Velocity vector distribution in the plane z=0.5 m with 5 different inclination angles(v=0.4 m·s-1)

圖9 為5 種不同傾斜角下y=0 m 的平面上的渦量分布云圖，該圖顯示不同傾斜角的魚礁渦量最大值均出現(xiàn)在迎流面與魚礁上表面相交處，但渦量范圍差異并不明顯，且除迎流面上方與上表面相交處的一小段弧形外，其余區(qū)域的渦量幾乎為0。

圖9 5 種不同傾斜角下y=0 m 平面上的渦量分布云圖(v=0.6 m·s-1)Fig.9 Vorticity distribution cloud maps on y=0 m plane under 5 different inclination angles (v=0.6 m·s-1)

3 討論與分析

3.1 相同傾斜角下的產(chǎn)生的展向渦和流向渦與背渦流的聯(lián)系

當來流受到無孔礁體斜面的阻擋，流場有一個主要的流向時，流經(jīng)礁體上表面的水體在其上方形成三角狀渦旋區(qū)，產(chǎn)生了一個橫向的渦，該渦旋與SWEARINGEN，et al[37]提到的對流向渦的概念相符合。而發(fā)現(xiàn)因無孔礁體斜面阻擋，來流在礁體斜面形成一對附著在其上的渦旋區(qū)，渦的方向和此流向垂直的渦旋也與WILLIAMSON[34]提到的展向渦描述相符。在穩(wěn)定流中，由于魚礁的不同面間構成交線(如迎流面與頂面)導致流體分離。隨著流體連續(xù)流動，不斷有動態(tài)穩(wěn)定漩渦產(chǎn)生在分離點下游的壁面上，且不斷從壁面上脫落，形成尾跡流。而脫落的漩渦將隨流體對流發(fā)生粘性擴散，正是因為粘性擴散作用，漩渦的渦流強度在運動過程中不斷減弱，直至消亡。從礁體壁面脫落的漩渦帶有大量能量，該能量的產(chǎn)生與消耗具體表現(xiàn)是在礁體上施加一個阻力。同一傾斜角的魚礁，在不同來流速度作用下，其周圍所產(chǎn)生的流速大小、渦量大小有明顯差異，將展向渦和流向渦與背渦流流場進行對比觀察后發(fā)現(xiàn)，展向渦寬度與背渦流流場寬度近似相等，流向渦高度大于展向渦高度，流向渦渦心位置處的渦流高度與背渦流的高度近似相等，展向渦和流向渦長度近似相等，并且與背渦流的長度近似相等，與李曉磊等[31]的研究描述高度相似。而速度方向及影響規(guī)模等變化的差異不大，這與黃遠東等[20]的研究結果一致。

3.2 傾斜角對上升流影響

OGAWA[2]、黃遠東等[17]、何文榮等[22]和姜少杰等[32]將水流Z 方向速度分量與來流速度之比大于5%的水域作為上升流區(qū)域，本研究中亦據(jù)此選取上升流數(shù)據(jù)。

通過分析y=0 m 平面內的Z 方向的速度分量分布，求得上升流參數(shù)Hm/H(最大高度Hm 與礁體高度H 之比),得到上升流作用高度與礁體傾斜角之間的關系。圖10 為在4 種不同來流速度下，各礁體所形成的最大上升流高度與礁體高度之比的折線圖。根據(jù)前人的研究，得到梯形人工魚礁Hm/H 值為2.49，箱形人工魚礁Hm/H 值為2.63，姜少杰等通過對正六棱柱式人工魚礁的模擬計算得到Hm/H為2.85，而文章中正四角錐臺形魚礁Hm/H 范圍在2.25～2.8 之間，與前人的結論相近。當來流速度相同時，隨著傾斜角的增加，曲線a(傾斜角68.2°)變化到曲線b(傾斜角73.3°)的過程中，上升流作用最大高度(Hm)與礁體高度(H)之比Hm/H 增速較大，而在曲線d(傾斜角84.3°)之后，雖然Hm/H 依然在增加，但是增速明顯放緩。a、b 兩種情況在低速條件下(0.2～0.4 m·s-1),其速度增加使Hm/H 有一定的提高，但大多數(shù)情況下，隨著速度的增大，上升流最大高度與魚礁高度之比Hm/H 幾乎不變化，即來流速度的變化與上升流最大高度無明顯的相互關系。

圖10 不同傾斜角、不同來流速度下Hm/HFig.10 The ratio of maximum upwelling height to reef height under different inclination angles and different inflow velocities

3.3 傾斜角對背渦流影響

本研究中數(shù)值計算試驗的背渦流模擬結果與TENEDA[30]采用流體顯示技術，在低雷諾數(shù)時圓柱繞流產(chǎn)生的圓柱尾渦(雙渦)的結果相吻合，TENEDA的試驗結論目前仍是數(shù)值模擬工作者驗證算法準確性的標準之一。漩渦的流速分布為由外向內逐漸減小，渦心處流速趨近于零，該結論與李曉磊等[31]的研究結果一致。

對背渦流的最大渦量值進行分析，如圖11 所示，隨著來流速度的增大，最大渦量值也會隨之增大，且其呈線性關系。而來流流速相同的情況下，隨著傾斜角的增大，最大渦量值呈現(xiàn)出一種先增后減的趨勢，且當傾斜角為78.7°時，相較于其余4 種情況，最大渦量值達到峰值。通過與何文榮等[22]的金字塔型魚礁后方背渦流情況進行比較發(fā)現(xiàn)，文章中5 種傾斜角的魚礁的背流區(qū)高度均高于金字塔型魚礁的背流區(qū)，這是由于本研究中的魚礁的背流面面積均大于金字塔型魚礁，且迎流面與地面的夾角同樣大于金字塔型魚礁。

圖11 不同傾斜角、不同來流速度下在y=0 平面的最大渦量值(s-1)Fig.11 The maximum vorticity value in the y=0 plane under different inclination angles and different inflow velocities (s-1)

4 結論

文章基于ANSYS 有限元分析軟件，對不同傾斜角的四角錐臺形人工魚礁在水流中的三維流場效應進行了數(shù)值仿真。文中對不同流速情況下同一傾斜角魚礁礁體周圍流態(tài)情況進行了分析，并在同一流速下改變礁體傾斜角，就上升流和背渦流的水動力學特性進行比較。由此得出以下結論：

(1)在不同來流速度下，同種傾斜角魚礁周圍流場流速大小、渦量大小有明顯差異，但是流體變化趨勢即速度方向和影響范圍相似。

(2)通過分析y=0 m 平面以及z=0.5 m 平面上流體的速度分布，揭示傾斜角對人工魚礁的流場效應具有顯著影響：1)礁體迎流面與上表面交匯處形成流動分離點，在礁體上表面形成小漩渦區(qū)，在礁體后方形成尾渦區(qū)。2)最大上升流高度與礁體高度之比以及上升流的規(guī)模大小均隨傾斜角的增大而增大，其值與來流速度大小無明顯關聯(lián)。3)魚礁的傾斜角越大，其背流區(qū)規(guī)模越大。4)隨著傾斜角的增大，魚礁背部的尾渦基本形態(tài)不變，規(guī)模略微增大，渦心向上偏移，且水平高度始終低于魚礁頂面高度。5)傾斜角越大，魚礁后方流態(tài)越雜。

(3)通過分析y=0 平面上流體的渦量分布，得出：1)魚礁迎流面上方與上表面相交處為渦量值最大處，其他區(qū)域的渦量幾乎為0。2)來流速度相同的情況下，隨傾斜角的增大，最大渦量值呈現(xiàn)一種先增大后減小的趨勢，相較于另外4 種情況，當傾斜角為模型c 情況下(78.7°)，最大渦量值達到峰值。

本研究尚存在一定不足之處，首先是未能開展模型試驗與數(shù)值模擬結果進行對比驗證。此外，為單純研究傾斜角的影響，魚礁被設計為無孔模型，避免開孔后流體干擾，但是，實際應用中的魚礁勢必為有孔設計。以上問題在后續(xù)的研究工作中將進一步修正和完善，為人工魚礁的優(yōu)化設計提供更全面的參考依據(jù)。