徐為兵, 朱克強, 張大朋, 楊然哲

(寧波大學(xué) 海運學(xué)院, 浙江 寧波 315211)

過度捕撈和近海海洋環(huán)境污染導(dǎo)致近岸海洋漁業(yè)資源枯竭。發(fā)展創(chuàng)新型漁業(yè), 大力調(diào)整漁業(yè)產(chǎn)業(yè)結(jié)構(gòu), 轉(zhuǎn)變漁業(yè)經(jīng)濟的增長方式是我國漁業(yè)的政策導(dǎo)向。深水抗風浪網(wǎng)箱是全球沿海國家開發(fā)海洋經(jīng)濟的重要戰(zhàn)略裝備, 也是我國新時期現(xiàn)代農(nóng)業(yè)向海洋領(lǐng)域拓展的重要戰(zhàn)略需要。我國深水網(wǎng)箱養(yǎng)殖事業(yè)起步比較晚, 與挪威、美國、日本等國在深水養(yǎng)殖領(lǐng)域存在相當差距。目前綜合各方面因素來看, 結(jié)構(gòu)形式較為簡單、管理較為方便而且價格較低廉的重力式網(wǎng)箱更適合我國的國情, 全國絕大多數(shù)采用的都是此類網(wǎng)箱。張雅芝[1-2]等從水產(chǎn)養(yǎng)殖的角度分析了我國海水網(wǎng)箱養(yǎng)殖長久以來存在的各種問題, 如網(wǎng)箱的布局、管理和魚類病害等, 這些問題制約著我國網(wǎng)箱養(yǎng)殖的發(fā)展。但是如果能在認識這些問題的基礎(chǔ)上, 加強對海水養(yǎng)殖容量的研究, 開發(fā)優(yōu)良品種,保護種質(zhì)資源并加強養(yǎng)殖水環(huán)境保護和病蟲防護,就能實現(xiàn)網(wǎng)箱養(yǎng)殖的可持續(xù)健康發(fā)展； 徐君義等[3]根據(jù)實際試驗海區(qū)網(wǎng)箱養(yǎng)殖的成果, 針對網(wǎng)箱的材料、魚種的選擇與運輸, 總結(jié)出了實際養(yǎng)殖過程中發(fā)現(xiàn)的一些經(jīng)驗； 王福表[4-5]則專門分析了網(wǎng)箱養(yǎng)殖水域污染源的成因和特點, 并根據(jù)實際考察給出了水污染治理的一些相關(guān)對策。同時, 作者針對網(wǎng)箱的有機質(zhì)淤積帶來的危害, 研制了一種在沉降物未到達水底前就能被收集清除的新型網(wǎng)箱。雖然目前浮式重力式網(wǎng)箱應(yīng)用較廣, 但它在強水流作用下, 網(wǎng)衣水平漂移嚴重, 網(wǎng)箱容積損失率高, 而加載保持網(wǎng)型的同時, 會使網(wǎng)衣承受的張力增加, 容易造成網(wǎng)衣撕破。因此, 在流速較大海域中, 該網(wǎng)箱的使用受到限制[6]?？梢? 對先進網(wǎng)箱的了解和研究是必不可少的, 比較典型的為美國式鋼質(zhì)升降式大型深水網(wǎng)箱, 國內(nèi)也稱之為飛碟型可潛式大型深水網(wǎng)箱。國內(nèi)外有關(guān)圓形重力式網(wǎng)箱數(shù)值模擬和模型實驗的研究成果比較豐富, 而目前為止國內(nèi)對于碟形網(wǎng)箱只有較為詳盡的結(jié)構(gòu)介紹和初步的水動力估算[6-10], 孫滿昌等[11], 李玉成等[12-13]對碟形網(wǎng)箱也僅僅是進行了模型實驗方面的研究。本文基于凝集參數(shù)法, 運用 OrcaFlex建立碟形網(wǎng)箱模型,通過數(shù)值模擬研究波流聯(lián)合作用下碟形網(wǎng)箱的纜繩受力特性和網(wǎng)箱結(jié)構(gòu)的運動響應(yīng), 以期深入了解碟形網(wǎng)箱的水動力特性, 為碟形網(wǎng)箱及其錨泊系統(tǒng)的研究提供參考。

1 數(shù)值模型

本文研究的對象為沉降式碟形網(wǎng)箱。該網(wǎng)箱為一個自張緊結(jié)構(gòu), 主要由一個八邊形浮環(huán)和一個中心立柱組成。網(wǎng)箱浮環(huán)由8個鋼管分段用法蘭連接,每個分段分別密封并加壓。立柱和浮環(huán)的材料均為鍍鋅鋼, 用16根輻繩連接。中心立柱直徑為0.92 m,長9.14 m, 浮環(huán)的公稱直徑為15 m, 單根浮環(huán)管徑為27.3 mm, 輻繩材質(zhì)為聚乙烯, 直徑13 mm。網(wǎng)箱用重達3 t的混凝土塊壓載, 通過懸掛索懸掛在中心立柱下面, 使之保證一定的穩(wěn)性。中心立柱內(nèi)有一個長達 3 m的體積可變的浮艙, 配合一定長度的懸掛索能使網(wǎng)箱下潛到特定的深度。網(wǎng)箱及錨泊系統(tǒng)的具體尺寸及浮標與網(wǎng)箱連接細節(jié)見參考文獻[14-15]。本文利用基于凝集參數(shù)法的專業(yè)水動力分析軟件OrcaFlex建立真實的網(wǎng)箱模型, 如圖1所示。有關(guān)凝集質(zhì)量法的理論可見朱克強等[16-18]的研究, 在此不再詳細介紹。由于本文研究的主要內(nèi)容是網(wǎng)箱整體和錨泊系統(tǒng)的水動力響應(yīng), 因此對網(wǎng)箱各種構(gòu)件之間復(fù)雜的連接作了必要的簡化, 如繩索之間的連接全部簡化為鉸接。同時為了節(jié)省計算時間, 本文沒有真實模擬網(wǎng)衣受力, 并選用工程上常用的 Airy波理論和均勻流進行各種組合來對網(wǎng)箱養(yǎng)殖系統(tǒng)施加不同的環(huán)境載荷, 通過時域模擬后得到了不同載荷下的網(wǎng)箱水動力結(jié)果, 結(jié)合水動力分析結(jié)果給出了深海網(wǎng)箱在布錨過程中應(yīng)注意的事項, 對實際工程實踐有一定的指導(dǎo)作用。

建模主要的思想是將整個網(wǎng)箱結(jié)構(gòu)分割成受外力作用的若干個單元, 然后將外力平均分配到每個對應(yīng)的節(jié)點, 形成一系列的運動方程。為了初步規(guī)劃和設(shè)計, 我們應(yīng)用修改后的 Morison方程(公式(1)),來計算作用在籠狀結(jié)構(gòu)小部件上的波浪力。

其中,m表示水質(zhì)點的質(zhì)量,ρ是海水的密度,Km是附加質(zhì)量系數(shù),CD是曳力系數(shù),CM是慣性系數(shù),A是投影面積, ?是網(wǎng)箱系統(tǒng)小部件排開水的體積,并且VR=V-R˙, 是水粒子速度V相對于結(jié)構(gòu)速度R˙。公式(1)右邊的第一項通常認為是曳力, 而第二項稱為慣性力。所有的這些力是由一個受到海洋環(huán)境負載的移動物體激發(fā)的。不過, 當處理一個特定的單元時, 其他外部力, 如浮力、物重、和錨繩或網(wǎng)線伸長時產(chǎn)生的張力, 計算時都應(yīng)包括在內(nèi)。因此, 由牛頓第二定律可得, 公式(1)可以表示為：

其中FD是阻力,FI的慣性力,FB浮力,FW是重力,FT是張力。

圖1 網(wǎng)箱模型示意圖Fig. 1 The model of the sea station

本文在建模過程中做了相應(yīng)的處理, 如網(wǎng)箱中作為加強筋的輻繩之間并沒有網(wǎng)衣, 因此網(wǎng)箱整體受到的水動力有所減少, 但本文探究的是網(wǎng)箱系統(tǒng)整體的運動規(guī)律, 對網(wǎng)衣細節(jié)的模擬有待后續(xù)進一步的研究。八邊形浮環(huán)之間通過6D浮標相連, 為了確保浮環(huán)作為一個整體參與水動力計算, 浮環(huán)之間的端部的強度設(shè)為無限大。浮標在此只是起到連接作用, 并不參與水動力計算。軟件中對細長型結(jié)構(gòu)均采用凝集質(zhì)量模型, 該方法不僅計算速度快, 還具有較高的計算精度, 已經(jīng)得到廣泛認可, 因此本文網(wǎng)箱系統(tǒng)模型中計算得到的錨泊張力具有一定的合理性。

2 算例分析

在深海區(qū)域, 常常會遇到臺風等惡劣天氣, 波浪和海流會同時出現(xiàn), 錨泊系統(tǒng)將限制網(wǎng)箱的運動,保證網(wǎng)箱結(jié)構(gòu)和魚群的安全。在波流聯(lián)合作用下, 網(wǎng)箱錨繩的受力將在一定程度上反映網(wǎng)箱整體結(jié)構(gòu)的受力。一旦錨泊系統(tǒng)失效, 不僅會造成魚群逃逸和死亡, 還會對網(wǎng)箱結(jié)構(gòu)造成嚴重破壞, 給養(yǎng)殖者帶來巨大的經(jīng)濟損失, 因此研究錨泊系統(tǒng)的安全性是很有必要的。

本文研究在水深為55 m的深海中, 各種波流組合對網(wǎng)箱系統(tǒng)的作用, 包括3個波高(2、4、6 m)和3個周期(6、9、12 s)以及 4 種流速(0～2 m/s)[19]。需要注意的是, 軟件中默認的波流前進方向均為 180°(圖中所示海流流向為0°)。圖2顯示的是在純波作用下,迎流面錨繩上的最大張力與波浪周期和波高H之間的關(guān)系。由圖可知, 周期相同時, 純波作用下的錨繩張力均與波高大致呈線性關(guān)系； 波高一定時, 錨繩的張力隨周期的變化并不明顯, 且周期越大, 錨繩張力的幾乎保持不變。

圖2 錨繩張力隨周期的變化Fig. 2 Changes of mooring line tension in various periods

圖3顯示了單個網(wǎng)箱系統(tǒng)在各種波流聯(lián)合作用下時迎流面錨繩張力的最大值。由圖3a～圖3c可知,波浪周期較長, 流速相同時, 相同波高作用下的網(wǎng)箱錨繩張力的最大值幾乎相同； 當波浪周期一定時(圖 3d～圖 3f), 波高越大, 網(wǎng)箱系統(tǒng)的錨繩張力最大值也越大； 且隨著流速的均勻增加, 錨繩張力最大值并不是線性增加, 而是呈現(xiàn)出增幅越來越快的趨勢。波浪周期越小, 流速越大, 水粒子的水平方向的合速度就越大, 由公式(1)可知, 網(wǎng)箱系統(tǒng)受到的水動力就越大, 波浪的波高越大, 水粒子的垂直方向上的速度就越大, 水動力也同樣會增大。

綜合對比以上各圖可以發(fā)現(xiàn), 在本文的研究中,波流參數(shù)對網(wǎng)箱錨繩最大張力影響的重要程度依次為： 流速V, 波高H, 周期T。

對于沉降式碟形網(wǎng)箱來說, 遇到臺風等惡劣天氣時迅速下沉到海底是保證網(wǎng)箱結(jié)構(gòu)安全完整的重要手段。因此研究極限海況下網(wǎng)箱下沉時的受力情況也很有必要。深海中的波浪和海流前進的方向并不總是相同, 為了研究斜入射波和海流對錨繩張力的影響, 我們假設(shè)均勻流流只沿著x軸正方向運動,但波浪入射的角度不同。

圖3 碟形網(wǎng)箱系統(tǒng)在不同波流狀態(tài)下的錨繩最大張力Fig. 3 The maximum mooring line tension induced by various wave-current conditions

圖4是在均勻流速為1.5 m/s時, 迎流面錨繩最大張力與入射波(H=6 m,T=9 s)前進角度的關(guān)系。從圖中可以看出, 錨繩的最大張力隨著入射波角度的變化呈周期性變化； 波流同向時要比波流逆向時的錨繩張力要大, 但最大張力發(fā)生在波流夾角為 135°左右的情況下, 與前人研究的結(jié)果相符； 而且從圖中可以看出, 網(wǎng)箱下沉能明顯地減少錨繩的受力,在本例中網(wǎng)箱下沉能使錨繩受力減少35%～43%。當波流夾角為 135°時, 在波浪前進的方向上 1#纜繩承受網(wǎng)箱整體受到的波浪力最大, 因此表現(xiàn)為此時 1#纜繩上的錨泊張力最大。本例計算得到的結(jié)果可以指導(dǎo)養(yǎng)殖戶在特定海洋環(huán)境下合理布置網(wǎng)箱的錨泊位置, 避免出現(xiàn)最大錨繩張力情況, 間接提高深海網(wǎng)箱系統(tǒng)的安全性, 減少不必要的損失。

圖4 錨繩張力隨入射波角度的變化Fig. 4 The maximum mooring line tension induced by various incident angles of waves

如圖 5所示, 取浮環(huán)兩側(cè)中點 A、B, 假設(shè)浮環(huán)在波流聯(lián)合作用下發(fā)生的轉(zhuǎn)動和變形幅度很小, 則浮環(huán)整體的傾角可以通過A、B兩點的連線與水平線之間的夾角θ來表示。本文計算了不同波流作用情況下網(wǎng)箱漂浮時的浮環(huán)最大傾角, 如圖6所示。

圖5 浮環(huán)傾角計算示意圖Fig. 5 The inclination calculation diagram of floating ring

圖6 碟形網(wǎng)箱系統(tǒng)在不同波流狀態(tài)下的浮環(huán)最大傾角Fig. 6 The maximum inclination induced by various wave-currents

由圖 6a～圖 6c可知, 波高一定時, 浮環(huán)的傾角與流速大致呈線性關(guān)系。周期較小時(6 s), 浮環(huán)傾角會明顯較大； 相同流速下, 周期較大(9 s和12 s)、波高較小時(2 m和4 m)浮環(huán)傾角的差值很小, 但這個差值會隨著流速的增加而逐漸變大。由圖6d圖6f可知, 周期一定時, 浮環(huán)傾角和流速呈線性關(guān)系, 隨著波高的增加而增大, 但增大的幅度隨著流速和周期的增加而減小。

3 結(jié)語

純波作用下, 波浪周期相同時, 錨繩張力隨波高線性增加； 波高一定時, 錨繩張力隨周期的變化并不明顯, 且隨著周期越來越大, 錨繩張力趨于不變。

波流聯(lián)合作用下, 波浪周期較長時, 錨繩的最大張力幾乎相同； 當波浪周期一定時, 波高越大, 錨繩最大張力也越大, 且流速越大, 錨繩最大張力隨波高增加的幅度越小。波流參數(shù)對網(wǎng)箱錨繩最大張力影響的重要程度依次為： 流速V、波高H、周期T。

海流方向一定時, 錨繩最大張力隨波浪的入射角呈周期性變化； 波流同向要比波流逆向時的錨繩張力要大, 且波流夾角為 135°時, 迎流面錨繩張力最大； 網(wǎng)箱下沉能明顯減小網(wǎng)箱運動幅度和錨繩受力。

各種工況下, 浮環(huán)傾角大小在 25°以內(nèi), 其值隨流速線性增加, 隨著波高的增加而增大, 但增大的幅度隨著流速和周期的增加而減小。

[1] 張雅芝. 我國海水魚類網(wǎng)箱養(yǎng)殖現(xiàn)狀及其發(fā)展前景[J].海洋科學(xué), 1995, 5： 21-24.

[2] 張雅芝, 蘇永全. 論我國海水魚類網(wǎng)箱養(yǎng)殖的可持續(xù)發(fā)展[J]. 海洋科學(xué), 2001, 7： 52-56.

[3] 徐君義, 金偉良. 淺?？癸L浪網(wǎng)箱養(yǎng)魚技術(shù)研究[J].海洋科學(xué), 2001, 1： 12-14.

[4] 王福表. 網(wǎng)箱養(yǎng)殖水污染及其治理對策[J]. 海洋科學(xué), 2002, 7： 24-26.

[5] 王福表. 水產(chǎn)養(yǎng)殖清淤網(wǎng)箱的研制和應(yīng)用[J]. 海洋科學(xué), 2002, 7： 12-14.

[6] 袁軍亭, 周應(yīng)祺. 深水網(wǎng)箱的分類及性能[J]. 上海水產(chǎn)大學(xué)學(xué)報, 2006, 03： 3350-3358.

[7] 郭建平, 吳常文. 美國式鋼質(zhì)升降式大型深水網(wǎng)箱結(jié)構(gòu)原理的研究探討[J]. 漁業(yè)現(xiàn)代化, 2004, 1： 28-31.

[8] 章守宇, 劉洪生. 飛碟型網(wǎng)箱的水動力學(xué)數(shù)值計算法[J].水產(chǎn)學(xué)報, 2002, 6： 519-527.

[9] 夏泰淳, 張健. 深水雙錐型網(wǎng)箱的阻力估算[J]. 上海水產(chǎn)大學(xué)學(xué)報, 2003, 1： 87-91.

[10] 劉健, 夏泰淳, 許柳雄. 深水雙錐型網(wǎng)箱的水動力估算[J].浙江海洋學(xué)院學(xué)報(自然科學(xué)版), 2007, 3： 257-265.

[11] 孫滿昌, 張健, 錢衛(wèi)國.飛碟型網(wǎng)箱水動力模型試驗與理論計算比較[J].上海水產(chǎn)大學(xué)學(xué)報, 2003, 4： 319-323.

[12] 李玉成, 宋芳, 董國海, 等. 碟形網(wǎng)箱水動力特性的研究[J].海洋工程, 2004, 4： 19-25.

[13] 李玉成, 桂福坤, 宋芳, 等. 重力式與碟形網(wǎng)箱的錨繩受力特性比較[J].海洋工程, 2005, 4： 19-24.

[14] Fredriksson D W, Muller E, Baldwin K, et al. Open ocean aquaculture engineering： system design and physical modeling[J]. Marine Technology Society Journal, 2000, 34(1)： 41-52.

[15] Tsukrov I I, Ozbay M, Swift M R, et al. Open ocean aquaculture engineering： numerical modeling[J]. Marine Technology Society Journal, 2000, 34(1)： 29-40.

[16] ZHU K Q, CAI Y, YU C L, et al. Nonlinear Hydrodynamic Response of Marine Cable-Body System Undergoing Random Dynamic Excitation[J]. Journal of Hydrodynamics. 2009, 21(6)： 851-855.

[17] 朱克強, 鄭道昌, 張玉松, 等. 隨機波浪下遙控潛水器強非線性耦合空間運動與纜索安全性研究[J]. 海洋工程, 2010, 2： 95-99.

[18] 朱克強, 朱海洋, 于春令, 等. 深海遙控潛水器多體系統(tǒng)非線性耦合動力特性模擬[J]. 海洋工程, 2008, 1： 83-87.

[19] Huang C C, Tang H J, Liu J Y. Effects of waves and currents on gravity-type cages in the open sea[J].Aquacultural Engineering, 2008, 38(2)： 105-116.