于定勇， 王逢雨， 鐘延超， 楊遠航

(中國海洋大學工程學院，山東 青島 266100 )

社會發(fā)展速度快，漁業(yè)資源需求量也隨之增大,近年來一些不合理的漁業(yè)發(fā)展模式,破壞了海洋生物的棲息地，導致海洋漁業(yè)資源急劇減少?！渡綎|省人工魚礁建設規(guī)劃(2014—2020年)》中提到,山東在2020年前在沿海區(qū)域將投入建設9大人工魚礁帶，40個人工魚礁群。

人工魚礁是指人們在水中經過科學選點而設置的構造物,是海洋牧場的重要組成部分,它可以為近海幼魚等水生生物提供避難所，以及促進上下層海水循環(huán)流動。放置人工魚礁后，可改變附近區(qū)域的流場效應，在海流、波浪的作用下水體上下混合和產生渦流，形成的上升流可引起海底營養(yǎng)物質的交互循環(huán)，能改善人工魚礁區(qū)周圍的海域生態(tài)環(huán)境，促進浮游生物的繁殖能力，使礁體附近形成較理想的索餌區(qū)，方便水生生物覓食，其暴露的表面可供海藻附著和生長，還可以避免使用非法漁具進行捕撈，保護漁業(yè)資源。研究人工魚礁的流場效應對增殖漁業(yè)資源具有非常重要的意義。

Falc?o等[1]對營養(yǎng)鹽通量和有機顆粒物在人工魚礁區(qū)和非人工魚礁區(qū)的差異進行了對比，分析了人工魚礁對水中營養(yǎng)鹽通量和有機物顆粒的影響。

Woo等[2]利用數(shù)值模擬方法對人工魚礁的水動力特性進行了研究，給出了包括透空立方體在內的24種不同形狀人工魚礁阻力系數(shù)值，結論表明阻力系數(shù)受礁體迎流角度的影響較明顯，與初始流速沒有太大關系。

Liu和Su[3]通過數(shù)值模擬研究了無底方型雙礁體連線垂直于來流方向和豎直堆疊布放時礁體的流場效應，結果表明雙礁體連線垂直于來流方向布設時流場效應較好。

唐衍力[4]對圓形開口魚礁進行了水槽實驗，礁體分有蓋和無蓋兩種，邊長為0.15 m、開口比為0.22，實驗得到了礁體在不同來流速度、不同迎流角度下的阻力系數(shù)。

劉彥[5]通過物理模型實驗和數(shù)值模擬方法研究了在不同來流速度下雙方體組合魚礁的流場效應，得到了組合礁體較優(yōu)的投放形式。

關長濤等[6]采用數(shù)值模擬方法分析了不同雷諾數(shù)條件下，擺放方式和布設間距對三圓管型人工魚礁流場效應的影響，得到了上升流和背渦流的規(guī)模和強度隨雷諾數(shù)、布設間距的變化情況。

總結前人工作成果可以了解到，對于各種擺放形式下的組合人工魚礁體，國內外學者主要從不同礁體形狀和海流速度對礁體水動力特性的影響進行了研究，尚缺乏對于開口雙方型礁體流場效應、阻力系數(shù)隨布設間距變化情況的科研工作，已有研究[7]表明上述流場效應在一定程度上受礁體布設間距變化的影響。為了探究礁體阻力系數(shù)和上升流特性參數(shù)隨布設間距的變化情況，本文通過數(shù)值模擬和物模實驗研究了雙方型人工魚礁體的水動力特性，為礁區(qū)布局方案的確定提供參考。

1 物模實驗

1.1 模型制作

圖1 礁體模型Fig.1 Reef model

1.2 實驗裝置及儀器

實驗在中國海洋大學水動力循環(huán)水槽中進行。水槽尺寸為5.0 m×1.2 m×1.2 m，其底部及兩側面均由透明玻璃構成，造流裝置位于水槽一端，通過調節(jié)其電機頻率改變流速。水槽中水深為1 m時可達到的水流速度為0～1.50 m/s，穩(wěn)定流速為0.15～0.80 m/s。本次實驗水深設置為0.95 m，當實際情況下海流速度為0.80 m/s時，根據重力相似準則，水槽進口流速設置約為0.18 m/s。實驗采用ADV聲學多普勒流速儀測量測點流速，采樣頻率為200 Hz，采用六分力儀測量礁體在水流中的受力，采樣頻率為100 Hz。

1.3 實驗方法

選擇開口比為0.3的方形開口雙礁體，實驗開始之前在礁體模型表面打孔，利用螺母將模型固定在直徑為0.01 m的螺桿上并垂直懸掛在水中，前方礁體模型迎流面距離水槽入口距離為2 m，通過移動后方礁體，改變兩礁體模型間的距離實現(xiàn)雙礁體布設間距的變化。模型底面與水槽邊壁相接近但不接觸，水槽邊壁作為模擬海底。假設礁體模型底部不受縫隙水流的作用，當礁體表面和水槽自由液面距離較遠時，不考慮自由液面的影響[9]。螺桿另一端與六分力儀相連接(見圖2、3)。通過信號調節(jié)器和A/D轉換儀將測量的數(shù)據信號傳輸至數(shù)據采集和處理系統(tǒng)。

圖2 模型的安裝示意圖Fig.2 Installation diagram of reef model

圖3 實驗所用六分力儀Fig.3 Six-dimentional force transducer

每次實驗開始之前，用流速儀測量距離第一個礁體模型迎流面前方0.6 m處的流速，使其穩(wěn)定流速接近0.178 m/s，然后移動流速儀測量礁體前后測點的流速，當待測點的流速測量值波動幅度較小時，采集30 s的流速數(shù)據后取平均值作為該測點的流速測量值。不同布設間距下測點位置如圖4所示。

圖4 不同布設間距下雙礁體模型測點位置示意圖(y=1.5 m)Fig.4 Sketch map of measurement points position of double reefs with different disposal spaces

2 數(shù)學模型

2.1 控制方程

假設人工魚礁體附近是不可壓縮的黏性流體進行的湍流運動，溫度幾乎不變，將能量方程忽略。

連續(xù)方程：

(1)

動量方程：

(2)

式中：ui(i=1，2，3)分別為x、y、z方向的雷諾平均速度；ρ為流體密度；p為壓強；v為運動黏性系數(shù)；fi為體積力。

2.2 湍流模型

本文采用RNGκ-ε兩方程模型計算黏性流體運動。此模型可以對湍流結構較小、分布較均勻的湍流流動[7]進行有效的模擬，對于本文人工魚礁體流場效應的研究比較適用。

湍動能κ方程：

(3)

湍流耗散率ε方程：

(4)

其中：

C1ε=1.42;C2ε=1.68；

Prandtl數(shù)：ακ=αε=1.39；

熱膨脹系數(shù)：β=0.012；

本文采用有限體積法離散控制方程，采用二階迎風格式對各方程進行空間離散，采用SIMPLEC算法進行壓力-速度耦合，采用標準差分格式進行壓力項處理，劃分非結構化網格，計算殘差值取10-5，通過數(shù)值模擬方法對方型礁體的水動力特性進行研究，分析礁體布設間距的變化對其上升流特性參數(shù)、阻力系數(shù)的影響。

2.3 魚礁結構及模擬計算區(qū)域

選取邊長L為3 m、開口比為0.3、沿水流方向布設間距分別為0.5L、1.0L、2.0L、4.0L的方形開口雙礁體，流速取為0.8 m/s。礁體結構如圖5所示，仿真計算區(qū)域見圖6。

圖5 開口比為0.3的方型人工魚礁體結構Fig.5 Sketch of cubic artificial reef with 0.3 opening ratio

圖6 布設間距為2.0L的雙礁體流場計算域示意圖Fig.6 Computational domain sketch of double reefs with disposal space of 2.0L

設定如下初邊界條件：

(1) 計算域入口設置來流速度為0.8 m/s的速度入口邊界條件(Velocity inlet)，給出邊界上湍動能κ和湍動耗散率ε，設置了邊界上各方向的速度矢量分量。

(2)計算域出口邊界設置為自由出流邊界條件(Outflow)。

(3)計算域的兩個側面設置為對稱邊界(Symmetry)。

(4)計算域頂面設置為具有與入口水流相同速度，剪切力為零的可移動壁面，礁體表面和計算域底面設置為無滑移壁面(Wall)。

2.4 模型可靠性驗證

2.4.1 網格收斂性驗證 網格收斂性驗證為了減小數(shù)值模擬過程中由于網格尺度所產生的誤差，以布設間距為2.0L的雙方型礁體舉例，采用不同的網格尺寸對計算域進行了劃分，以前面礁體的阻力系數(shù)作為變量進行了網格收斂性驗證，結果見表1。

由表1可以看出當布設間距為2.0L的雙礁體計算域最大網格尺寸為0.500 m時，收斂性較好，此時網格尺寸對礁體阻力系數(shù)影響較小，因此本文采用0.500 m作為數(shù)值模擬時的最大網格尺寸，礁體表面網格高度第一層邊界層設為0.002 m，增長率為1.1，共設置10層。

表1 不同網格尺寸模擬結果Table 1 Simulation results with different gird sizes

基于上述的數(shù)值模型，模擬研究了當布設間距分別為0.5L、1.0L、2.0L、4.0L時，開口比為0.3的雙礁體上升流特性參數(shù)和阻力系數(shù)的變化情況。

2.4.2 實驗驗證 實驗驗證礁體阻力系數(shù)是表征人工魚礁體穩(wěn)定性的重要參數(shù)。本文數(shù)模時僅考慮來流速度不變時礁體的受力情況，通過下式求得阻力系數(shù)：

(5)

式中:F為礁體沿水流方向受力(N)；ρ為海水密度(kg/m3)；A為礁體迎流面積(m2)；u為水流速度(m/s)。

將Fluent軟件模擬得到的測點對應位置處的流速和前方礁體阻力系數(shù)值與物模實驗結果進行對比，j見圖7和表2。

從圖7中可以看出在第一個礁體迎流面前方測點流速模擬值與實驗值相對誤差較小，在兩礁體間及第二個礁體背流面后方測點流速相對誤差較大，測點流速模擬值和實驗值最大相對誤差分別為10.22%和10.76%，相對誤差較大的測點出現(xiàn)在兩礁體之間和第二個礁體后方，水流流經礁體后在此處流場變得復雜，流速變化劇烈。進行物模實驗測量流速時，偶然誤差不可避免地會出現(xiàn)，同時由于此處流速較小，較小的流速變化可能帶來較大的相對誤差，除去此點，其他測點流速模擬值和實驗值相差較小，二者吻合較好。

圖7 不同布設間距雙礁體測點流速模擬值與實驗值比較Fig.7 Comparison of the simulated and measured velocities of the measuring stations of the double reefs with different disposal spaces

表2給出了前方礁體阻力系數(shù)數(shù)值模擬值和實驗值的相對關系，由表2可以得到在布設間距為2.0L時，前方礁體阻力系數(shù)模擬值與實驗值之間相對誤差較大，約為1.64%。整體來說前方礁體阻力系數(shù)模擬值與實驗值結果吻合較好。

表2 前方礁體阻力系數(shù)模擬值與實驗值比較Table 2 Comparison of the simulated and measured Cd of the front reef

由圖7和表2可以看出本文通過數(shù)值模擬得到的雙礁體測點流速和礁體受力與物模實驗結果較為吻合，表明本文模擬計算所采用的方法是可行的，數(shù)值模擬結果是可信的。

3 水動力特性研究

3.1 流場形態(tài)

通過數(shù)值模擬得到了開口比為0.3的方形開口雙礁體在布設間距為0.5L、1.0L、2.0L、4.0L時的水動力特性。本文將上升流區(qū)定義為水流z方向上速度分量與來流速度之比大于或等于5%的區(qū)域，這是由黃遠東等[10]提出的；緩流區(qū)定義為水流X方向上速度分量與來流速度之比絕對值小于40%的區(qū)域。

由圖8可以看出后方礁體產生的上升流區(qū)主要集中于迎流面上方尖角處。隨著布設間距的增加，后方礁體產生的上升流區(qū)的規(guī)模增大，但整體上升流區(qū)體積逐漸減小。兩礁體之間沒有上升流區(qū)的分布。表3中Hmax/H表示上升流最大高度/礁高，Wmax/W表示上升流水平跨度/礁寬，Vmax/V表示上升流體積/礁體體積。由表3可得，上升流特性參數(shù)值隨雙礁體布設間距的增加而逐漸減小，當布設間距為0.5L時，上升流特性參考值最大。

圖8 雙礁體縱向中軸斷面上升流區(qū)分布圖Fig.8 Velocity contour of upwelling zone on center longitudinal section of double reefs

圖9表示了開口比為0.3的雙礁體在來流速度為0.8 m/s的條件下z=1.5 m截面上的速度矢量分布。由圖9和表4可以看出當布設間距為0.5L時，緩流區(qū)最大影響范圍與后方礁體背流面距離約為4.8倍礁高，當布設間距增加至2.0L時，緩流區(qū)最大影響范圍與后方礁體背流面距離約為6.3倍礁高，當布設間距繼續(xù)增加至4.0L時，由于前方礁體對后方礁體的遮蔽效應減弱，緩流區(qū)最大影響范圍與后方礁體背流面距離減小至5.6倍礁高。模擬工況下，隨著布設間距的增加，緩流區(qū)體積逐漸增大。

表3 上升流特性參數(shù)隨布設間距的變化Table 3 Upwelling current characteristic parameters changing with different disposal spaces

圖9 不同布設間距雙礁體橫向中軸斷面速度云圖Fig.9 Velocity contour on transverse section of double reefs with different disposal spaces

表4 緩流區(qū)體積隨布設間距的變化Table 4 Volume of subcritical flow area changing with different disposal spaces

3.2 阻力系數(shù)

由圖10可以看出當兩礁體布設間距為0.5L～2.0L時，前方礁體阻力系數(shù)值變化較小，當間距為4.0L時，前方礁體阻力系數(shù)值較大，約為1.30。隨著布設間距的增加，前方礁體對后方礁體的遮蔽效應減弱，后方礁體受到的水流作用力逐漸增大，阻力系數(shù)值逐漸增大，當兩礁體布設間距增加至4.0L時，后方礁體阻力系數(shù)值增大至0.47。

圖10 雙礁體阻力系數(shù)與布設間距的關系Fig.10 Relation between the drag coefficient and disposal space

兩礁體沿水流方向布設時，后方礁體受前方礁體的遮流阻力影響程度與布設間距有關，當兩礁體布設間距與礁高的比值D/L較小時，后方礁體處于前方礁體后紊亂流場區(qū)域范圍內，受前方礁體遮流影響較嚴重，甚至可能處于前方礁體后方負壓區(qū)，受吸力而出現(xiàn)阻力負值；隨著D/L值增大，后方礁體遠離前方礁體復雜的流場區(qū)域，受前方礁體的遮流影響變得微弱，處于正壓區(qū)，阻力較大；當兩礁體間距與礁高比D/L增大到一定程度后，后方礁體將處于前方礁體流場恢復區(qū)段，遮流影響微弱到可忽略不計。

(6)

圖11 雙礁體kz與布設間距的關系Fig.11 Relation between kz and disposal space

圖11給出了礁體遮流阻力影響系數(shù)kz隨其布設間距比的變化關系。通過擬合得到雙礁體kz值與布設間距比的關系式為：

kz=0.003 2(D/L)2+0.039(D/L)+0.11。

(7)

為了提高擬合公式的可靠性，增加了3種工況，布設間距分別為3.0L、5.0L、6.0L。當布設間距比為1.0L時，kz模擬計算值為0.174，擬合公式計算值為0.152，相對誤差較大，約為12.64%。整體來說kz模擬計算值和擬合公式計算值吻合較好。由圖11可以看出kz值隨著布設間距的增加而增大，在布設間距較大時，kz值較大，說明此時兩礁體之間的流場受到礁體結構的干擾程度較小。

4 結論

本文通過數(shù)值模擬計算得到的礁體周圍測點流速及受力值與物模實驗實測值吻合較好，說明本文所采用的計算方法是可行的。因此，本文利用數(shù)值模擬研究了不同布設間距下方型人工魚礁體的水動力場，通過分析礁體布設間距對其阻力系數(shù)和上升流特性參數(shù)的影響，總結為以下三點：

(1)人工魚礁體的水動力特性隨布設間距的變化而產生一定的差異。在模擬工況下，隨著布設間距的增大，雙礁體產生的上升流體積逐漸減小，緩流區(qū)體積及遮流阻力影響系數(shù)kz均逐漸增大。

(2)當雙礁體布設間距在0.5L～6.0L之間時，遮流阻力影響系數(shù)kz與布設間距比的關系可用下式表述：

kz=0.003 2(D/L)2+0.039(D/L)+0.11。

(3)當雙礁體沿水流方向布設間距為0.5L時，其上升流特性參數(shù)值最大。當雙礁體沿水流方向布設間距為2.0L時，其緩流區(qū)最大影響范圍最廣。由此可見，流場效應較好的布設間距為0.5L～2.0L。研究結果可為礁區(qū)布局方案的確定提供參考。