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

(浙江海洋大學(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)點(diǎn)，目前國內(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è)計了近海中上層浮魚礁，實(shí)際投放于象山港海洋牧場，結(jié)果表明浮魚礁能較快地吸引魚類。在實(shí)際應(yīng)用中，浮魚礁通常投放在≥10 m水深海域[8]。受到天氣和潮差影響，海水波高和水深變化會導(dǎo)致錨泊系統(tǒng)和礁體位移產(chǎn)生變化[11]。浮魚礁不是單一地投放，在投放人工浮魚礁時，要注意與其他用途的海洋區(qū)域有足夠的距離，避免由于礁體的移動而影響海域的其他功能[12]。因此，在確保浮魚礁系統(tǒng)安全的前提下，當(dāng)波高和水深發(fā)生變化時，利用數(shù)值模擬的方法，研究框架式浮魚礁的水動力特性，計算錨繩受力大小和礁體位移能為浮魚礁實(shí)際投放間隔距離提供理論參考，計算錨泊系統(tǒng)受力情況能對錨繩材料選擇爭取利益最大化提供幫助。

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

人工魚礁主要有箱形礁體、框架形礁體、三角形礁體、梯形礁體和異形礁體等[13]。為便于比較，定義ε和η為優(yōu)化參數(shù)，ε為礁體空方體積與礁體體積之比，η為礁體表面積與礁體體積之比，對于增值附著類海洋生物礁體，其η值要高些，若是以誘集魚類為目的，其ε值應(yīng)高些[8]。為保證浮魚礁具有良好的集魚效果，設(shè)計時通常要遵循的基本原則：較大的礁體表面積、良好的透空性、充分的透水性以及切實(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]。在框架浮體底端中點(diǎn)處系縛錨繩，錨繩錨碇點(diǎn)用重塊沉到水底。

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

圖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 框架式浮魚礁模型建立

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

2.1 框架浮體受力與運(yùn)動

2.1.1 框架浮體受力

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

(1)

(2)

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

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

vR=v水-va

(3)

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

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

(4)

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

因?yàn)樵谡w坐標(biāo)系下只考慮X、Z兩個方向的受力，所以直接通過向量運(yùn)算建立浮體的受力公式，如式(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ì)量點(diǎn)所包含的受力進(jìn)行累加，并將其分配到集中質(zhì)量點(diǎn)上，然后利用牛頓第二定律建立質(zhì)點(diǎn)運(yùn)動方程，簡化后的質(zhì)點(diǎn)運(yùn)動方程[15-17]為：

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

(8)

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

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

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

(9)

式中：va—相對于固定坐標(biāo)系的絕對速度，m/s；aa—相對于固定坐標(biāo)系的絕對加速度，m/s2；ve、ae—牽連速度、牽連加速度，由坐標(biāo)相對運(yùn)動導(dǎo)致，m/s，m/s2；vr—剛體中點(diǎn)相對坐標(biāo)系的相對速度，m/s；vk—科氏速度，對于剛體而言，其值為0；ar—剛體中點(diǎn)相對坐標(biāo)系的相對加速度，m/s2，對于剛性浮架而言，相對速度和相對加速度為0；ak—科氏加速度，m/s2，是牽連運(yùn)動和相對運(yùn)動之間相互影響導(dǎo)致的，對于剛體而言同樣為0。綜上分析，剛體在固定坐標(biāo)系(絕對坐標(biāo)系)中的運(yùn)動速度和加速度實(shí)際上就是牽連速度和牽連加速度，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)原點(diǎn)，實(shí)際上常取動坐標(biāo)和固定坐標(biāo)質(zhì)心重合點(diǎn)。本文因?yàn)橐紤]到錨繩和框架浮體之間的相互作用，剛體繞質(zhì)心轉(zhuǎn)動的角速度為ω，這一轉(zhuǎn)動角速度實(shí)際上即為動坐標(biāo)相對于固定坐標(biāo)的牽連角速度[18]，即ωe=ω，由于轉(zhuǎn)動存在，剛體中任意點(diǎn)的牽連速度都不一樣，剛體中任意一點(diǎn)的牽連角速度為：

ve2=ωe×r

(11)

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

va=ve1+ve2

(12)

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

2.2 錨繩受力與運(yùn)動

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

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

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

(13)

式中：M—集中質(zhì)量點(diǎn)的質(zhì)量，kg;a—各集中質(zhì)量點(diǎn)的加速度矢量，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個點(diǎn)為一組參數(shù)進(jìn)行數(shù)據(jù)分析，其中上點(diǎn)為上構(gòu)件中點(diǎn)坐標(biāo)，中點(diǎn)(是一個虛擬的點(diǎn))為框架浮體質(zhì)心坐標(biāo)，下點(diǎn)為下構(gòu)件中點(diǎn)坐標(biāo)。當(dāng)只有波高變化時，初始狀態(tài)錨繩長5 m，下端系縛在錨碇點(diǎn)，上端系縛在框架下點(diǎn)，此時框架上點(diǎn)距離水面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個點(diǎn)的運(yùn)動軌跡如圖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)角隨波高變化

定義沿波浪傳播方向?yà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ì)點(diǎn)速度很大，導(dǎo)致水質(zhì)點(diǎn)相對浮體的運(yùn)動速度較大，浮體受波浪力隨之增加，從而產(chǎn)生劇烈運(yùn)動。圖4中，X軸左右偏移量基本相等、正負(fù)偏轉(zhuǎn)角度基本相等，說明在規(guī)則波浪作用下，框架式浮魚礁的左右搖擺運(yùn)動幅度基本保持一致。

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

波高3 m，水深分別為8、9、10、11、12 m時，框架浮體上中下3個點(diǎn)的運(yùn)動軌跡以及最大轉(zhuǎn)角隨水深而變化，此時浮體最上端距離水面分別為2、3、4、5、6 m(圖5)。定義框架浮體最上端(圖1中E’H’中點(diǎn)位置)至靜水面的距離與水深的比值為框架浮體的相對入水深度。圖5表明，當(dāng)相對入水深度為1/4時，框架浮體水平偏移最大，隨著水深增加；從1/4至2/5變化時，浮體最大水平偏移迅速減??；從2/5至1/2變化時，浮體在水平方向幾乎沒有運(yùn)動，這是由于隨著水深的增加，水質(zhì)點(diǎn)運(yùn)動速度迅速減小，導(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)角度影響十分顯著，而實(shí)際海洋中由于潮差的存在，相對入水深度會發(fā)生變化，因此在實(shí)際設(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)在水底錨碇點(diǎn)處。因此，在實(shí)際安裝時要充分考慮到當(dāng)極端天氣出現(xiàn)時，錨碇點(diǎ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)在水底錨碇點(diǎn)處。綜上所述，水深的變化對錨繩最大拉力的影響十分顯著，隨著水深的增加，浮魚礁浮體附近的水質(zhì)點(diǎn)速度迅速減小，當(dāng)水深超過11 m時，水質(zhì)點(diǎn)運(yùn)動速度接近于0，此時浮魚礁浮體受波浪力幾乎為0 N。實(shí)際在海洋牧場布置浮魚礁時，必須依據(jù)當(dāng)?shù)睾Ｓ虻某蔽蛔兓?，按照低潮位時計算錨繩的最大受力，評估浮魚礁的安全性能。

圖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ā)生在水底錨碇點(diǎn)處。波高一定，當(dāng)相對入水深度小于2/5時，浮體最大水平偏移、最大偏轉(zhuǎn)角度和錨繩最大拉力隨水深增加而急劇減小；當(dāng)超過2/5時，則隨水深增加變化緩慢；當(dāng)達(dá)到1/2時，前兩項幾乎沒有變化，錨繩最大拉力和靜水時錨繩拉力幾乎相等。在實(shí)際海洋牧場建設(shè)中，投放浮魚礁應(yīng)考慮目標(biāo)海域潮差和潮位變化導(dǎo)致相對入水深度變化帶來的影響。當(dāng)相對入水深度小于2/5時，其變化對錨繩的拉力和浮體的運(yùn)動都有十分劇烈的影響。水深一定時，波高的變化也會對錨繩受力和浮體運(yùn)動產(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ù)的試驗(yàn)研究[J].海洋學(xué)報,2006,28(5):145-151.

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