框架式浮魚礁波浪水動力特性數(shù)值模擬研究

張 健， 潘 昀， 馮德軍， 桂福坤

(浙江海洋大學(xué)，國家海洋設(shè)施養(yǎng)殖工程技術(shù)研究中心,浙江 舟山 316002)

由于過度捕撈和海洋污染，漁場遭到破壞，保護(hù)和修復(fù)海洋漁場勢在必行。人工魚礁是用來改善海洋生態(tài)環(huán)境、建設(shè)海洋牧場的人工設(shè)施，能為海洋生物提供繁殖、生長、棲息的場所，同時復(fù)雜多變的空間更能為魚類提供躲避區(qū)[1-4]。人工魚礁投放后，周圍海域的流、光、音、底質(zhì)等非生物環(huán)境會發(fā)生變化。這種變化又會引起生物環(huán)境的變化，增大生物量，從而形成良好的漁場或增養(yǎng)殖場，集魚效果也非常明顯[5]。由于沉式魚礁具有投放簡單、成本低、材料來源廣、集魚效果顯著等優(yōu)點，目前國內(nèi)應(yīng)用很廣。然而，沉式魚礁在黏土、淤泥以及沙質(zhì)海區(qū)可能會導(dǎo)致礁體整體下沉，過多的淤泥還會覆蓋甚至掩埋著生在礁體上的生物，降低含氧量，阻礙光線，影響魚礁功能[6]。由于長江攜帶大量懸浮泥沙，使舟山海域附近海底大量淤泥堆積，因此，沉式魚礁并不適合舟山海域，而浮魚礁則能很好地解決這個問題。

浮魚礁是繼沉式魚礁之后的一種海洋增養(yǎng)殖設(shè)施，主要是為誘集中上層魚類而設(shè)計。浮魚礁的作用范圍和機(jī)理與沉式魚礁不同，主要是由錨塊、錨繩和浮體組成，有些浮魚礁還會配有浮球[7]。浮魚礁根據(jù)作用區(qū)域不同，可分為表層、中層、底層和深海等。表層和中層浮魚礁分別作用在水深10 m和20～50 m附近海域。日本在浮魚礁的研究應(yīng)用上處于世界領(lǐng)先地位，擁有適用于大水深和全水層的浮魚礁，且已投入應(yīng)用，并取得了很好的效果。如：礁體位置為海面的表層浮魚礁，礁體位置為海面以下 20 m的中層浮魚礁，整體高度90 m左右、應(yīng)用于大水深海域的深海浮魚礁[8]。國內(nèi)對浮魚礁的研究較少，余求妹等[9]設(shè)計探討了浮繩式網(wǎng)箱人工浮魚礁，張麗珍等[10]設(shè)計了近海中上層浮魚礁，實際投放于象山港海洋牧場，結(jié)果表明浮魚礁能較快地吸引魚類。在實際應(yīng)用中，浮魚礁通常投放在≥10 m水深海域[8]。受到天氣和潮差影響，海水波高和水深變化會導(dǎo)致錨泊系統(tǒng)和礁體位移產(chǎn)生變化[11]。浮魚礁不是單一地投放，在投放人工浮魚礁時，要注意與其他用途的海洋區(qū)域有足夠的距離，避免由于礁體的移動而影響海域的其他功能[12]。因此，在確保浮魚礁系統(tǒng)安全的前提下，當(dāng)波高和水深發(fā)生變化時，利用數(shù)值模擬的方法，研究框架式浮魚礁的水動力特性，計算錨繩受力大小和礁體位移能為浮魚礁實際投放間隔距離提供理論參考，計算錨泊系統(tǒng)受力情況能對錨繩材料選擇爭取利益最大化提供幫助。

1 框架式浮魚礁結(jié)構(gòu)參數(shù)

人工魚礁主要有箱形礁體、框架形礁體、三角形礁體、梯形礁體和異形礁體等[13]。為便于比較，定義ε和η為優(yōu)化參數(shù)，ε為礁體空方體積與礁體體積之比，η為礁體表面積與礁體體積之比，對于增值附著類海洋生物礁體，其η值要高些，若是以誘集魚類為目的，其ε值應(yīng)高些[8]。為保證浮魚礁具有良好的集魚效果，設(shè)計時通常要遵循的基本原則：較大的礁體表面積、良好的透空性、充分的透水性以及切實的可行性。本文選擇框架形礁體為研究對象。魚礁由框架浮體(圖1中12條高密度聚乙烯(HDPE)管道棱邊熱熔焊接而成的空心框架)和錨繩兩部分組成，具體尺寸見表1?？蚣芨◇w尺寸為0.4 m×0.4 m×1 m(長×寬×高)，管道直徑0.05 m，浮體密度400 kg/m3。錨繩原長5 m，直徑0.01 m，錨繩密度950 kg/m3。采用耐磨耐腐蝕的聚乙烯(PE)材質(zhì)[14]。在框架浮體底端中點處系縛錨繩，錨繩錨碇點用重塊沉到水底。

假設(shè)波浪作用下框架浮體始終不出水，并且波浪傳播方向與框架浮體六面中的ABCD面平行，那么框架浮體的運動僅在波浪傳播方向發(fā)生，即可視為二維剛體運動(圖1)。那么，框架浮體運動可以簡化為EFGH平面上的二維平動和轉(zhuǎn)動，進(jìn)而研究對象分離為框架浮體E’F’G’H’和錨繩之間的運動。

圖1 框架式浮魚礁的基本構(gòu)成及坐標(biāo)

項目材料尺寸/m直徑/m密度/(kg/m3)拖曳力系數(shù)CD慣性力系數(shù)CM錨繩彈性系數(shù)C1錨繩彈性系數(shù)C2框架浮體HDPE0.4×0.4×1.00.054000.81.2——錨繩PE5.00.019500.61.2345.3×1061.0121

2 框架式浮魚礁模型建立

波浪作用下框架浮體及其錨繩運動數(shù)值模擬方法主要采用有限單元法計算框架浮體的受力和繞質(zhì)心的力矩，進(jìn)而離散時間確定框架浮體的平動和轉(zhuǎn)動，采用有限集中質(zhì)量點法計算柔性錨繩的受力和形變后的拉力，同樣離散時間確定錨繩集中質(zhì)量點的運動速度和空間坐標(biāo)。由于Fortran軟件具有強大的計算能力，因此先用其進(jìn)行數(shù)值模擬，計算框架浮體位移和受力情況，然后用Matlab軟件對數(shù)據(jù)結(jié)果進(jìn)行可視化處理。

2.1 框架浮體受力與運動

2.1.1 框架浮體受力

數(shù)模分析采用有限單元法[15-16]，根據(jù)已有文獻(xiàn)將AB管道劃分為10份即可滿足模擬精度。由于BC管道只有AB管道的0.4倍，BC管道就被劃分為4個單元，因此，框架浮體E’F’G’H’被劃分為10×4個有限單元(高×長)?？蚣芨◇w單位受力主要包括重力、浮力、波浪力等，在系縛點(圖1中P點)處計算錨繩拉力。采用線性波理論，根據(jù)波浪要素的彌散關(guān)系，求解波浪場的水質(zhì)點速度(式1)和水質(zhì)點加速度(式2)，再使用莫里森方程計算框架浮體受到的波浪力。

(1)

(2)

式中：VX—X軸水質(zhì)點速度，m/s；VZ—Z軸水質(zhì)點速度，m/s；aX—X軸水質(zhì)點加速度，m/s2；aZ—Z軸水質(zhì)點加速度，m/s2；H—波高，m；k—波數(shù)；x—水質(zhì)點橫坐標(biāo)，m；z—水質(zhì)點縱坐標(biāo)，m；d—水深，m；ω—圓頻率，rad/s；T—周期，s；t—步長，s。

框架浮體四邊均為圓柱形桿件，因此具有圓柱的水動力性質(zhì)，其水動力系數(shù)與水質(zhì)點的相對運動速度方向有關(guān)[15-17]。在計算每個浮體單元波浪力時應(yīng)考慮浮體與波浪入射方向的夾角。本文建立整體坐標(biāo)系(X,Y,Z)，同時在浮體單元處建立局部坐標(biāo)系(δ,η,τ)，τ方向為浮體單元方向，在計算τ方向時，會出現(xiàn)4個浮體單元方向，本文定義向上和向右為正方向。δ軸在τ方向和水質(zhì)點相對速度VR組成的平面內(nèi)與τ軸垂直，η軸與τ和δ軸方向組成的平面相垂直。

vR=v水-va

(3)

式中：vR—水質(zhì)點與浮體單元的相對速度，m/s；v水—水質(zhì)點速度，m/s；va—浮體單元速度，m/s。

在整體坐標(biāo)下，局部坐標(biāo)系的δ,η,τ軸單位矢量可以通過已知的水質(zhì)點相對速度VR與浮體單元方向向量τ叉乘得到，為：

(4)

式中：eδ,eη和eτ分別為局部坐標(biāo)系(δ,η,τ)3個軸的單位矢量坐標(biāo)[15]。

因為在整體坐標(biāo)系下只考慮X、Z兩個方向的受力，所以直接通過向量運算建立浮體的受力公式，如式(5)所示。

(5)

式中：FD浮體、FI浮體—浮體單元受到的速度力和慣性力，N；FDη、FDτ—速度力FD浮體在局部坐標(biāo)系上的分量；FIη、FIτ—慣性力FI浮體在局部坐標(biāo)系上的分量。展開為式(6)和式(7)。

(6)

(7)

式中：ρw—水體密度，kg/m3；CD—局部坐標(biāo)系的速度力項系數(shù)；Aη、Aτ—浮體單元在局部坐標(biāo)系η和τ方向上的投影面積，m2；ψ—浮體單元體積，m3；CM—局部坐標(biāo)系上的慣性力項系數(shù)。

在整體坐標(biāo)系下將浮體各集中質(zhì)量點所包含的受力進(jìn)行累加，并將其分配到集中質(zhì)量點上，然后利用牛頓第二定律建立質(zhì)點運動方程，簡化后的質(zhì)點運動方程[15-17]為：

Ma=T+FD浮體+FI浮體+W+B

(8)

式中：M—集中質(zhì)量點的質(zhì)量，kg；a—各集中質(zhì)量點的加速度矢量，m/s2；T—集中質(zhì)量點所受到的張力矢量，N；W—質(zhì)點的重力矢量，N；B—質(zhì)點的浮力矢量，N。

2.1.2 框架浮體的平動與轉(zhuǎn)動

本文中浮架的運動可以看作剛體運動，剛體運動一般包括平動和轉(zhuǎn)動，這也說明牽連速度和牽連加速度包含平動和轉(zhuǎn)動兩部分。因此，根據(jù)剛體運動的速度和加速度合成原理，剛體任意點的運動速度(加速度)分別由牽連速度(加速度)、相對速度(加速度)和科氏速度(加速度)矢量合成，即：

(9)

式中：va—相對于固定坐標(biāo)系的絕對速度，m/s；aa—相對于固定坐標(biāo)系的絕對加速度，m/s2；ve、ae—牽連速度、牽連加速度，由坐標(biāo)相對運動導(dǎo)致，m/s，m/s2；vr—剛體中點相對坐標(biāo)系的相對速度，m/s；vk—科氏速度，對于剛體而言，其值為0；ar—剛體中點相對坐標(biāo)系的相對加速度，m/s2，對于剛性浮架而言，相對速度和相對加速度為0；ak—科氏加速度，m/s2，是牽連運動和相對運動之間相互影響導(dǎo)致的，對于剛體而言同樣為0。綜上分析，剛體在固定坐標(biāo)系(絕對坐標(biāo)系)中的運動速度和加速度實際上就是牽連速度和牽連加速度，va=ve，aa=ae。

(1)計算平動。單考慮平動時，剛體無轉(zhuǎn)動位移，計算公式為：

(10)

式中：r0—平動距離,m;ve1—平動部分牽連速度，m/s；ae1—平動部分牽連加速度，m/s2。

(2)計算轉(zhuǎn)動。在分析剛體的轉(zhuǎn)動時，當(dāng)動坐標(biāo)和固定坐標(biāo)靜止不重合時，剛體轉(zhuǎn)動是對應(yīng)固定坐標(biāo)原點，實際上常取動坐標(biāo)和固定坐標(biāo)質(zhì)心重合點。本文因為要考慮到錨繩和框架浮體之間的相互作用，剛體繞質(zhì)心轉(zhuǎn)動的角速度為ω，這一轉(zhuǎn)動角速度實際上即為動坐標(biāo)相對于固定坐標(biāo)的牽連角速度[18]，即ωe=ω，由于轉(zhuǎn)動存在，剛體中任意點的牽連速度都不一樣，剛體中任意一點的牽連角速度為：

ve2=ωe×r

(11)

式中：r—對應(yīng)點原點的矢徑。至此可以得出整個剛體運動的絕對速度。

va=ve1+ve2

(12)

式中：ve2—剛體中的任意一點的牽連角速度，rad/s。

2.2 錨繩受力與運動

錨繩屬于典型柔性結(jié)構(gòu)物，也可以采用集中質(zhì)量法進(jìn)行模擬。假定錨繩由有限的無質(zhì)量彈簧連接的集中質(zhì)量點構(gòu)成，通過計算集中質(zhì)量點在動力邊界條件作用下的偏移來獲得錨繩的形狀[15-17]。模型的集中質(zhì)量點設(shè)于錨繩兩端，每個集中質(zhì)量點包含錨繩上下2個錨繩(圖2)。

圖2 錨繩節(jié)點示意圖

在淺海條件下，錨繩集中質(zhì)量點受力、計算波浪力方法與框架計算方法相同。在整體坐標(biāo)系下將各集中質(zhì)量點所包含錨繩的受力進(jìn)行累加，并將其分配到集中質(zhì)量點上，然后利用牛頓第二定律建立質(zhì)點運動方程。錨繩模型簡化后的質(zhì)點運動方程[19-20]為：

(13)

式中：M—集中質(zhì)量點的質(zhì)量，kg;a—各集中質(zhì)量點的加速度矢量，m/s2;L0—錨繩原始長度，m；L—變形后的長度，m；L0—錨繩原始長度,m;C1、C2為框架浮體的材料彈性系數(shù);FD、FI—分別是錨繩單元受到的速度力和慣性力,N;ε—伸縮比；dm—錨繩直徑，m。

3 結(jié)果與分析

研究了波高為1.5、2、2.5和3 m，水深為8、9、10、11和12 m，周期為2 s工況下，浮架水平方向最大偏移、最大偏轉(zhuǎn)角度和錨繩最大受力。定義框架浮體E’F’管道為左構(gòu)件，F(xiàn)’G’管道為下構(gòu)件，G’H’管道為右構(gòu)件，H’E’管道為上構(gòu)件。為便于分析，本文在整體坐標(biāo)系下，選取框架浮體上、中、下3個點為一組參數(shù)進(jìn)行數(shù)據(jù)分析，其中上點為上構(gòu)件中點坐標(biāo)，中點(是一個虛擬的點)為框架浮體質(zhì)心坐標(biāo)，下點為下構(gòu)件中點坐標(biāo)。當(dāng)只有波高變化時，初始狀態(tài)錨繩長5 m，下端系縛在錨碇點，上端系縛在框架下點，此時框架上點距離水面2 m，水深8 m，波高變化不影響初始狀態(tài)；當(dāng)只有水深變化時，框架和錨繩初始狀態(tài)同上，水深每增加1 m，系統(tǒng)初始狀態(tài)縱坐標(biāo)向下增加1 m。

3.1 波高對浮體最大偏移和最大轉(zhuǎn)角的影響

在水深8 m，波高分別為1.5、2、2.5、3 m的工況下，框架浮體上中下3個點的運動軌跡如圖3中a、b、c、d所示。可以看出，當(dāng)波高3 m時浮體水平最大偏移出現(xiàn)在±0.3 m附近，當(dāng)波高1.5 m時浮體水平最大偏移出現(xiàn)在±0.14 m附近，波高每增加0.5 m，浮體最大偏移相應(yīng)增加約0.05 m，波高3 m的最大偏移量約為波高1.5 m時最大偏移量的2倍。相同工況下，正負(fù)方向最大轉(zhuǎn)角隨波高的變化如圖4所示。

圖3 水平最大偏移隨波高變化情況

圖4 正負(fù)方向最大轉(zhuǎn)角隨波高變化

定義沿波浪傳播方向為正方向，與波浪傳播方向相反為負(fù)方向。由圖4可知，隨著波高增加，浮體正負(fù)方向偏轉(zhuǎn)角度增加，波高每增加0.5 m，正向最大轉(zhuǎn)角增加1.2°，負(fù)向最大轉(zhuǎn)角增加1.4°。

波高的變化對浮體正負(fù)向最大偏轉(zhuǎn)角度的影響基本一致。當(dāng)波高增大后，浮體最大偏移和正負(fù)方向最大偏轉(zhuǎn)角度均增大，原因是波浪的水質(zhì)點速度很大，導(dǎo)致水質(zhì)點相對浮體的運動速度較大，浮體受波浪力隨之增加，從而產(chǎn)生劇烈運動。圖4中，X軸左右偏移量基本相等、正負(fù)偏轉(zhuǎn)角度基本相等，說明在規(guī)則波浪作用下，框架式浮魚礁的左右搖擺運動幅度基本保持一致。

3.2 水深對浮體最大偏移和最大轉(zhuǎn)角的影響

波高3 m，水深分別為8、9、10、11、12 m時，框架浮體上中下3個點的運動軌跡以及最大轉(zhuǎn)角隨水深而變化，此時浮體最上端距離水面分別為2、3、4、5、6 m(圖5)。定義框架浮體最上端(圖1中E’H’中點位置)至靜水面的距離與水深的比值為框架浮體的相對入水深度。圖5表明，當(dāng)相對入水深度為1/4時，框架浮體水平偏移最大，隨著水深增加；從1/4至2/5變化時，浮體最大水平偏移迅速減??；從2/5至1/2變化時，浮體在水平方向幾乎沒有運動，這是由于隨著水深的增加，水質(zhì)點運動速度迅速減小，導(dǎo)致浮體受波浪力也迅速減小。

圖5 最大偏移量隨水深變化情況

相應(yīng)地，圖6所示，當(dāng)相對入水深度從1/4至2/5變化時，浮體正負(fù)方向最大偏轉(zhuǎn)角度迅速減小，當(dāng)浮體從距離水面2/5到1/2變化時，浮體正向和負(fù)向的最大偏轉(zhuǎn)角度變化趨于平緩；當(dāng)?shù)竭_(dá)1/2時，浮體正負(fù)方向最大偏轉(zhuǎn)角度幾乎沒有變化。考慮到水深變化對浮體的偏移和偏轉(zhuǎn)角度影響十分顯著，而實際海洋中由于潮差的存在，相對入水深度會發(fā)生變化，因此在實際設(shè)計投放浮魚礁時，應(yīng)充分考慮漲落潮對浮魚礁系統(tǒng)的相對入水深度、最大偏移量和浮體最大偏轉(zhuǎn)角度的影響。綜上，浮體上端至水面的距離與水深的比值為0.5時，浮魚礁偏移量和正負(fù)偏轉(zhuǎn)角度約等于0，即浮魚礁保持靜止。

圖6 正負(fù)方向最大轉(zhuǎn)角隨水深變化

3.3 錨繩最大受力隨波高變化

在水深8 m，波高分別為1.5、2、2.5、3 m的工況下，錨繩的最大拉力變化如圖7所示。

圖7 錨繩最大拉力隨波高變化

錨繩最大拉力隨著波高的增加而增大，波高每增加0.5 m，錨繩最大受力約增加40 N，近似呈線性變化。計算發(fā)現(xiàn)錨繩最大拉力位置出現(xiàn)在水底錨碇點處。因此，在實際安裝時要充分考慮到當(dāng)極端天氣出現(xiàn)時，錨碇點的極限抗撥能力，保證整個浮魚礁系統(tǒng)的安全性。

3.4 錨繩最大受力隨水深變化

圖8中的兩條線分別代表當(dāng)波高3 m時錨繩的最大拉力和波高0 m時的錨繩拉力。從圖中可以看出，錨繩的最大拉力隨著水深的增加而減小，當(dāng)水深從8 m到10 m時，錨繩最大拉力呈對數(shù)函數(shù)減小，從10 m開始，隨著水深增加，錨繩最大拉力變化緩慢，當(dāng)水深11 m時，錨繩的最大拉力幾乎和靜水時的拉力相等，并在其附近小幅度震蕩，震蕩的原因與錨繩的彈性應(yīng)變系數(shù)有關(guān)。當(dāng)水深超過11 m時，錨繩的最大拉力隨著水深的增加幾乎沒有變化。同樣，錨繩最大拉力位置出現(xiàn)在水底錨碇點處。綜上所述，水深的變化對錨繩最大拉力的影響十分顯著，隨著水深的增加，浮魚礁浮體附近的水質(zhì)點速度迅速減小，當(dāng)水深超過11 m時，水質(zhì)點運動速度接近于0，此時浮魚礁浮體受波浪力幾乎為0 N。實際在海洋牧場布置浮魚礁時，必須依據(jù)當(dāng)?shù)睾Ｓ虻某蔽蛔兓凑盏统蔽粫r計算錨繩的最大受力，評估浮魚礁的安全性能。

圖8 錨繩最大拉力隨水深變化

4 結(jié)論

通過數(shù)值模擬方法研究了框架式浮魚礁在波浪條件下的水動力特性。當(dāng)水深一定時，波高每增加0.5 m，浮體最大水平偏移約增加0.05 m，浮體正負(fù)向最大偏轉(zhuǎn)轉(zhuǎn)角增加1.2°，負(fù)向最大偏轉(zhuǎn)轉(zhuǎn)角增加1.4 °，錨繩最大受力增加40 N，錨繩最大張力發(fā)生在水底錨碇點處。波高一定，當(dāng)相對入水深度小于2/5時，浮體最大水平偏移、最大偏轉(zhuǎn)角度和錨繩最大拉力隨水深增加而急劇減小；當(dāng)超過2/5時，則隨水深增加變化緩慢；當(dāng)達(dá)到1/2時，前兩項幾乎沒有變化，錨繩最大拉力和靜水時錨繩拉力幾乎相等。在實際海洋牧場建設(shè)中，投放浮魚礁應(yīng)考慮目標(biāo)海域潮差和潮位變化導(dǎo)致相對入水深度變化帶來的影響。當(dāng)相對入水深度小于2/5時，其變化對錨繩的拉力和浮體的運動都有十分劇烈的影響。水深一定時，波高的變化也會對錨繩受力和浮體運動產(chǎn)生影響，因此，在設(shè)計浮魚礁時，應(yīng)將低潮位時的水深及設(shè)計波高作為設(shè)計動力，從而選取合適的錨繩和浮體的投放位置，保證浮魚礁投放的有效性和安全性。

□

參考文獻(xiàn)

[1] 周艷波,蔡文貴,陳海剛,等.人工魚礁生態(tài)誘集技術(shù)的機(jī)理及研究進(jìn)展[J].海洋漁業(yè),2010,32(2):225-230.

[2] SHERMAN R L,GILLIAM D S,SPIELER R E. Artificial reef design:void space,complexity,and attractants [J]. ICES Journal of Marine Science,2002,59(S):196-200.

[3] 林承剛,汝少國,楊紅生,等.刺參對人工礁體設(shè)計關(guān)鍵指標(biāo)的選擇性[J].海洋科學(xué),2012,36(3):13-21.

[4] PEDRO G,MARCOS F,LUCIANO D S,et al. Impact of artificial patchy reef design on the ichthyofauna community of seasonally influenced shores at Southeastern Brazil[J]. Aquatic Ecology,2015,49(3):343-355．

[5] 陳勇,于長清,張國勝,等.人工魚礁的環(huán)境功能與集魚效果[J].大連水產(chǎn)學(xué)院學(xué)報,2002(1):64-69.

[6] 趙海濤,張亦飛,郝春玲,等.人工魚礁的投放區(qū)選址和礁體設(shè)計[J].海洋學(xué)研究,2006(4):69-76.

[7] 劉惠飛.日本人工魚礁研究開發(fā)的最新動向[J].漁業(yè)現(xiàn)代化,2002,29(1):25-27.

[8] 夏章英,盧伙勝,馮波,等.人工魚礁工程學(xué)[M].北京:海洋出版社,2011:80-81.

[9] 余求妹,馬家志,安玉,等.浮繩式網(wǎng)箱人工浮魚礁的設(shè)計優(yōu)勢及問題的探討[J].安徽農(nóng)業(yè)科學(xué),2013,41(19):8194-8195,8217.

[10]張麗珍,王江濤,胡慶松,等.近海中上層柔性浮魚礁設(shè)計與應(yīng)用[J].上海海洋大學(xué)學(xué)報,2016,25(4):613-619.

[11]余建星,韓彤,李林奇,等.海洋資料浮標(biāo)水動力特性研究現(xiàn)狀分析[J].海洋湖沼通報,2017(4):61-67.

[12]李文濤,張秀梅.關(guān)于人工魚礁礁址選擇的探討[J].現(xiàn)代漁業(yè)信息,2003(5):3-6.

[13]張劍誠,金海,王吉橋.人工漁礁建設(shè)研究現(xiàn)狀[J].水產(chǎn)科學(xué),2004(11):27-30.

[14]王彩云. HDPE高密度聚乙烯管道的應(yīng)用研究[J].企業(yè)技術(shù)開發(fā),2012,31(14):159-160.

[15]趙云鵬.深水重力式網(wǎng)箱水動力特性數(shù)值模擬研究[D].大連:大連理工大學(xué),2007.

[16]趙云鵬,李玉成,董國海,等.水流作用下重力式網(wǎng)箱網(wǎng)衣張力分布[J].漁業(yè)現(xiàn)代化,2008,35(6):58-63.

[17]桂福坤.深水重力式網(wǎng)箱水動力學(xué)特性研究[D].大連:大連理工大學(xué),2006.

[18]陳小芳.鲆鰈類方形網(wǎng)箱水動力特性數(shù)值模擬[D].大連:大連理工大學(xué),2012.

[19]陳天華,孟昂,桂福坤.波浪高度及方向?qū)吨絿W(wǎng)養(yǎng)殖系統(tǒng)網(wǎng)片水力特性的影響[J].農(nóng)業(yè)工程學(xué)報,2017,33(2):245-251.

[20]李玉成,桂福坤.平面無結(jié)節(jié)網(wǎng)衣水阻力系數(shù)的試驗研究[J].海洋學(xué)報,2006,28(5):145-151.

[21]楊正東,朱建榮,王彪,等.長江河口潮位站潮汐特征分析[J].華東師范大學(xué)學(xué)報(自然科學(xué)版),2012(3):111-119.