王東旭 陳國寶 湯勇 李斌 王志超

摘要 [目的]對大亞灣南部海域漁業(yè)資源的生物資源密度和總資源量進行水聲學(xué)評估。[方法]2015年利用Simrad EY60聲學(xué)系統(tǒng)進行了春季、夏季、秋季和冬季4個航次的聲學(xué)調(diào)查，并以傳統(tǒng)拖網(wǎng)方式同步采集生物學(xué)資料，對大亞灣南部海域的生物資源密度和總資源量進行了評估。[結(jié)果]2015年4月春季調(diào)查航次的平均生物資源密度為7.38×103 kg/n mile2，總資源量為172.8 t，優(yōu)勢種為二長棘鯛（平均生物資源密度為2.37×103 kg/n mile2）；2015年8月夏季調(diào)查航次的平均生物資源密度為8.58×103 kg/n mile2，總資源量為200.9 t，優(yōu)勢種為黃鰭馬面鲀（平均生物資源密度為3.29×103 kg/n mile2）；2015年10月秋季調(diào)查航次平均生物資源密度為6.52×103 kg/n mile2，總資源量為152.6 t，優(yōu)勢種為短吻鲾（平均生物資源密度為2.59×103 kg/n mile2）；2015年12月冬季調(diào)查航次的平均生物資源密度為6.36×103 kg/n mile2，總資源量為148.9 t，優(yōu)勢種為皮氏叫姑魚（平均生物資源密度為0.82×103 kg/n mile2）。[結(jié)論]該研究結(jié)果可為大亞灣南部海域海業(yè)資源的養(yǎng)護和管理提供科學(xué)依據(jù)。

關(guān)鍵詞 水聲學(xué)評估；平均生物資源密度；總資源量；大亞灣

中圖分類號 S931 文獻標(biāo)識碼 A 文章編號 0517-6611（2017）06-0095-04

Underwater Acoustics Assessment on Fishery Resources in the Southern Daya Bay

WANG Dong-xu1，2，CHEN Guo-bao2，3*，TANG Yong1 et al （1.School of Marine Technology and Environment，Dalian Ocean University，Dalian，Liaoning 116023； 2.South China Sea Fisheries Research Institute，Chinese Academy of Fishery Sciences，Guangzhou，Guangdong 510300； 3.Scientific Observing and Experimental Station of South China Sea Fishery Resources and Environment，Ministry of Agriculture，Guangzhou，Guangdong 510300）

Abstract [Objective] To make underwater acoustics assessment on biomass density and total biomass of fishery resources in the southern Daya Bay.[Method] Using Simrad EY60 acoustics system，four times of acoustics survey was made in four seasons of 2015，and biological data were synchronously sampled by traditional trawling method.The biomass density and total biomass in southern Daya Bay were assessed.[Result] In spring survey，the average biomass density in June of 2015 was 7.38×103 kg/n mile2，its total biomass was 172.8 t，and the dominant species was Paerargyrops edita（average biomass density of 2.37×103 kg/n mile2）.In summer survey，the average biomass density in August of 2015 was 8.58×103 kg/n mile2，the total biomass was 200.9 t，the dominant species was Navodon xanthopterus（average biomass density of 3.29×103 kg/n mile2）.In autumn survey，the average biomass density in October of 2015 was 6.52×103 kg/n mile2，total biomass was 152.6 t，the dominant species was Leiognathus brevirostris（average biomass density of 2.59×103 kg/n mile2）.In winter survey，the average biomass density in December of 2015 was 6.36×103 kg/n mile2，the total biomass was 148.9 t，the dominant species was Johnius belangerii（average biomass density of 0.82×103 kg/n mile2）.[Conclusion] The research results can provide scientific basis for the maintenance and management of fishery resources in the southern Daya Bay.

Key words Underwater acoustics assessment；Average biomass density；Total biomass；Daya Bay

大亞灣是我國南海北部的重要海灣，位于廣東省東南部珠江河口東側(cè)，海灣岸線長達92 km，曲折多變，被惠州平海半島、惠陽南部沿海以及深圳大鵬半島三面環(huán)繞。海灣面積約600 km2，灣內(nèi)大小島嶼有50多個[1]，其優(yōu)越的自然環(huán)境很適合海洋生物的棲息繁衍，因而具有豐富的漁業(yè)資源。近年來，由于人類活動的影響使灣內(nèi)生態(tài)環(huán)境惡化由貧營養(yǎng)狀態(tài)發(fā)展成中營養(yǎng)狀態(tài)，生物群落組成明顯小型化，生物多樣性降低，再加上漁業(yè)的過度捕撈使灣內(nèi)漁業(yè)生物資源衰退嚴重[2-3]。為恢復(fù)大亞灣的生物資源量，2006年廣東省在大亞灣投放魚苗20多萬尾，2007年深圳市在大亞灣內(nèi)進行海洋牧場建設(shè)，投放體積9.51 km3的人工魚礁[4]，并對大亞灣進行漁業(yè)資源量的動態(tài)監(jiān)測。

傳統(tǒng)的漁業(yè)資源評估是利用拖網(wǎng)漁獲物數(shù)據(jù)來估算海域資源量[5-7]，而漁業(yè)資源聲學(xué)評估是一種新的評估方法。漁業(yè)資源聲學(xué)評估技術(shù)與其他評估方法相比具有效率高、數(shù)據(jù)采集簡便、調(diào)查范圍廣泛和不影響被調(diào)查海域的生態(tài)環(huán)境等優(yōu)點，因此近幾年得到廣泛應(yīng)用。國外在漁業(yè)資源聲學(xué)評估方面已有較為成熟的研究[8-10]，國內(nèi)于1984年引進，近年來圍繞聲學(xué)資源評估的研究日益增多，內(nèi)陸淡水方面開展了東湖、黃河小浪底和水庫方面的聲學(xué)資源評估[11-14]；海洋方面，聲學(xué)資源評估的研究大多集中在南海和北太平洋[15-20]。針對海灣水域漁業(yè)資源的聲學(xué)評估，主要有對拓林灣南澳島海洋牧場漁業(yè)資源和大亞灣楊梅坑人工魚礁區(qū)漁業(yè)資源的聲學(xué)評估等[21-22]。筆者于2015年4、8、10和12月利用Simrad EY60聲學(xué)系統(tǒng)對大亞灣南部海域進行了4次走航式漁業(yè)資源調(diào)查，結(jié)合生物學(xué)采樣數(shù)據(jù)利用聲學(xué)評估技術(shù)對大亞灣南部海域漁業(yè)資源現(xiàn)存狀況進行了評估，旨在為其漁業(yè)資源的養(yǎng)護與管理提供科學(xué)依據(jù)。

1 材料與方法

1.1 調(diào)查海域

研究調(diào)查海域如圖1所示，調(diào)查范圍為114°34.0′～114°44.0′ E，22°32.0′～22°38.0′ N，主要包括大亞灣楊梅坑人工魚礁區(qū)和大亞灣中央列島區(qū)2個區(qū)域。通常聲學(xué)調(diào)查航線依據(jù)調(diào)查區(qū)域的特征，設(shè)計成平行式或“之”字式[23]，由于該調(diào)查區(qū)域海域環(huán)境復(fù)雜，島嶼眾多，因而實際調(diào)查航線并不規(guī)則。依據(jù)海域?qū)嶋H情況，調(diào)查航線設(shè)計成環(huán)島形，基本覆蓋設(shè)計調(diào)查海域，主要圍繞大辣甲、小辣甲和牛頭洲等灣內(nèi)島嶼。設(shè)置D1～D8共8個生物學(xué)數(shù)據(jù)采集站位，這8個站位的生物學(xué)數(shù)據(jù)能夠代表調(diào)查區(qū)域的魚類信息，為后續(xù)聲學(xué)數(shù)據(jù)的處理與計算提供依據(jù)。

1.2 數(shù)據(jù)采集

調(diào)查數(shù)據(jù)采集包含聲學(xué)數(shù)據(jù)采集和生物學(xué)數(shù)據(jù)采集2個部分。聲學(xué)數(shù)據(jù)是使用Simrad EY60分裂波束科學(xué)魚探儀進行采集，采集方式為走航式聲學(xué)數(shù)據(jù)采集。聲學(xué)系統(tǒng)參數(shù)需要根據(jù)調(diào)查海域的實際校準結(jié)果進行設(shè)定，在進行聲學(xué)調(diào)查前對該聲學(xué)系統(tǒng)進行校準，校準地點為調(diào)查海域內(nèi)海況較穩(wěn)定的近岸海域。研究用的魚探儀頻率為200 kHz，換能器型號ES200-7C，序列號587，發(fā)射功率150 W，脈沖寬度0.128 ms，增益27.25 dB，波束類型為分裂波束，雙向波束寬度為20.70 dB，縱向波束角為6.69°，橫向波束角為6.74°。

聲學(xué)數(shù)據(jù)是在走航過程中錄入，通過專業(yè)軟件ER60保存在電腦中，外接GPS直接連接電腦，其位置數(shù)據(jù)和聲學(xué)數(shù)據(jù)同步保存。換能器置于導(dǎo)流罩內(nèi)通過連接鋼管固定于船只右舷中間位置，導(dǎo)流罩前后連接細鋼絲繩，以減少行船過程中導(dǎo)流罩前后晃動。導(dǎo)流罩需置于水深1 m以下，以減少船舶航行過程中水表層產(chǎn)生的氣泡和水流的噪聲影響，依據(jù)實際情況調(diào)節(jié)入水深度，避免船只左右搖擺及換能器露出水面形成數(shù)據(jù)空白。

生物學(xué)數(shù)據(jù)的采集是依據(jù)調(diào)查方案預(yù)先設(shè)定的站位，利用租用漁船在相應(yīng)位置進行拖網(wǎng)。拖網(wǎng)時間為0.5 h，計時從下網(wǎng)結(jié)束至開始起網(wǎng)。漁船船號為“粵惠灣漁16009”，發(fā)動機功率135 kW，船長15.5 m，船寬5.6 m，排水量50 t。采樣網(wǎng)具為拖網(wǎng)，網(wǎng)長4.8 m，網(wǎng)口寬2.6 m，網(wǎng)口高0.4 m，網(wǎng)目3 cm×3 cm。對大部分漁獲物進行現(xiàn)場測量，獲得生物學(xué)數(shù)據(jù)；部分漁獲物現(xiàn)場不能完成測量，冰凍后帶回實驗室進行處理。

1.3 數(shù)據(jù)分析和處理

聲學(xué)數(shù)據(jù)是利用Echoview軟件進行分析和處理，數(shù)據(jù)處理前要對全部數(shù)據(jù)進行仔細檢查，對照聲學(xué)記錄表選取有效聲學(xué)數(shù)據(jù)。在數(shù)據(jù)處理過程中，需要去除的噪聲主要包括船舶產(chǎn)生的機械噪聲、海表噪聲、海底噪聲和浮游生物噪聲[24]。該研究的聲學(xué)數(shù)據(jù)中還存在一些停船作業(yè)時段（環(huán)境調(diào)查采樣），查找聲學(xué)記錄表對應(yīng)時間和經(jīng)緯度，手動去除停船時段數(shù)據(jù)。根據(jù)實際情況，設(shè)置積分輸出上限和積分輸出下限，即海表線和海底線。2015年4月春季調(diào)查海表線設(shè)為1.5 m，海底線設(shè)為-0.2 m； 2015年8月夏季調(diào)查海表線設(shè)為1.0 m，海底線設(shè)為-0.2 m； 2015年10月秋季調(diào)查海表線設(shè)為1.0 m，海底線設(shè)為 -0.5 m； 2015年12月冬季調(diào)查海表線設(shè)為1.5 m，海底線設(shè)為-0.5 m。積分閾值設(shè)為-70 dB，以屏蔽弱散射體浮游生物等產(chǎn)生的回波信號，以0.25 n mile為航程基本輸出單元，輸出積分值NASC，即平均后向散射截面（Mean backscattering cross-section）。

根據(jù)站位信息采集的生物樣本按照《海洋漁業(yè)資源調(diào)查規(guī)范》進行現(xiàn)場生物學(xué)測量，不能進行現(xiàn)場測量的部分漁獲樣本冰凍后帶回實驗室進行處理。需要測量的生物學(xué)數(shù)據(jù)包括尾數(shù)（個）、體長（mm）、體重（g）。

生物學(xué)數(shù)據(jù)是輔助聲學(xué)數(shù)據(jù)來計算評估資源量，依據(jù)拖網(wǎng)所獲生物學(xué)數(shù)據(jù)對評估魚類進行積分值分配。

聲學(xué)數(shù)據(jù)處理結(jié)果以0.25 n mile航程為計算單元，包括輸出水層（Layer）、平均后向散射強度（Sv_mean）、平均后向散射截面積分值（NASC）、數(shù)據(jù)輸出水層的平均高度和深度（Height_mean、Depth_mean）、每個輸出單元開始“呯”和結(jié)束“呯”（Ping_E、Ping_E）、每個輸出單元開始時間節(jié)點（Time_M）和每個輸出單元開始節(jié)點的經(jīng)度與緯度（Lon_M、Lat_M）等。

參照《多種類海洋漁業(yè)資源聲學(xué)評估技術(shù)和方法》分別對4個航次調(diào)查數(shù)據(jù)進行計算，評估資源量。首先確定所評估種類目標(biāo)強度參數(shù)值，根據(jù)相關(guān)資料確定目標(biāo)強度和魚類體長存在一定的函數(shù)關(guān)系，目標(biāo)強度-體長關(guān)系公式[25]如下：TS=20logl+b20 （1）式中，TS表示魚類的目標(biāo)強度（dB），l表示魚體體長（cm），b20表示評估種類的目標(biāo)強度參數(shù)。其中，l可根據(jù)生物學(xué)數(shù)據(jù)確定，b20值通過查閱相關(guān)文獻[14，17，22]確定。評估計算過程中所包含種類的b20值如表1所示。

在Excel中輸出NASC積分總值，對積分值分配時參考生物學(xué)數(shù)據(jù)，積分值分配公式如下：SA，i=SAniiki=1nii （2）式中，SA，i表示第i種魚的積分值，SA表示積分總值，ni表示第i種魚的尾數(shù)，i表示第i種魚的平均聲學(xué)散射截面，k表示參與積分值分配魚類種類總數(shù)。

平均聲學(xué)散射截面i[23]的計算公式為：i=4π（）2·10b2010 （3）式中，表示魚的平均體長（mm），可通過生物學(xué)數(shù)據(jù)計算得出；b20為目標(biāo)強度相關(guān)參數(shù)值，通過查閱資料確定。

生物資源密度是以0.25 nmile航程為基本單元，參照相近站位的生物學(xué)拖網(wǎng)數(shù)據(jù)進行計算，具體計算公式為：ρi=iSA，i1 000 i（4）式中， ρi表示評估魚種中第i種魚的生物資源密度（kg/n mile2），i表示第i種魚的平均體重（g），SA，i表示第i種魚分配到的積分值（m2/n mile2），i表示第i種魚的平均聲學(xué)散射截面（m2）。

依據(jù)每種魚類的平均生物資源密度，計算大亞灣內(nèi)評估區(qū)域的總資源量，其計算公式為：B=A （5）式中，B表示評估區(qū)域內(nèi)某一種魚類的總資源量（kg），表示評估魚種中某一種魚類的生物資源密度（kg/n mile2），A為評估區(qū)域面積。

2 結(jié)果與分析

2.1 生物學(xué)數(shù)據(jù)分析

漁業(yè)資源量評估主要是依據(jù)聲學(xué)數(shù)據(jù)進行計算，生物學(xué)數(shù)據(jù)只是輔助聲學(xué)數(shù)據(jù)進行資源評估計算。聲學(xué)數(shù)據(jù)無法判斷水體中的生物種類以及每種所占總生物量的比例系數(shù)。通過生物學(xué)采樣數(shù)據(jù)可以確定評估種類，以及每一種占總評估種類的比例系數(shù)、平均體長和平均體重。該研究中4個航次調(diào)查都進行了生物學(xué)采樣，采樣生物種類數(shù)量的變化情況如圖2所示。

2015年4月春季調(diào)查總共捕獲魚類和頭足類40種，除去鰕虎魚類、鲆鰈魚類等貼底生活的魚類不參與積分值分配，共25種魚類和3種頭足類進行資源評估計算。漁獲物種類數(shù)量大于10尾的有10種，包括二長棘鯛、黃斑藍子魚、短吻鲾、四線天竺鯛、竹筴魚、日本紅娘魚、鹿斑鲾、六指馬鲅、少鱗鱚和花鰶。其中，二長棘鯛共捕獲1 062尾，在參與評估漁獲物中占71%，具有明顯的種群優(yōu)勢。黃斑藍子魚和短吻鲾在參與評估漁獲物中所占比例分別為9%和6%，其余種類所占比例均小于5%。

2015年8月夏季調(diào)查共捕獲64種魚類和5種頭足類，其中魚類46種和頭足類5種參與積分值分配進行資源評估計算。由于8月禁漁期剛剛結(jié)束，與春季航次相比漁獲種類和數(shù)量都有明顯增長。漁獲物數(shù)量大于10尾的種類有20種，黃鰭馬面鲀、四線天竺鯛、短吻鲾?shù)牟东@量均大于100尾。其中，黃鰭馬面鲀?yōu)樽顑?yōu)勢種，捕獲量為2 490尾，占比71%。春季航次的優(yōu)勢種二長棘鯛只捕獲到51尾，占比1.4%。

2015年10月秋季調(diào)查捕獲魚類39種和頭足類3種，參與積分值分配進行資源評估計算的魚類有23種和頭足類3種，漁獲總量與夏季航次相比有所下降，優(yōu)勢種與前2個航次相比也發(fā)生變化。優(yōu)勢種為短吻鲾，捕獲1 124尾所占比例為65%。夏季航次優(yōu)勢種黃鰭馬面鲀只捕獲36尾占比2.1%，未捕獲到春季航次優(yōu)勢種二長棘鯛。該航次中細條天竺魚和中線天竺鯛捕獲量也較大，分別為195尾（占比11%）和221尾（占比12%）。此外，金線魚、斑鰶和短蛸的捕獲數(shù)量均大于10尾。

2015年12月冬季調(diào)查漁獲物捕獲種類和數(shù)量與秋季相比都有所下降，共捕獲33種魚類和2種頭足類，參與積分值分配計算的有魚類23種和頭足類2種。未出現(xiàn)捕獲尾數(shù)超過100尾的種，其中頭足類杜氏槍烏賊捕獲48尾占比19.5%為捕獲量最大的種，黑邊天竺魚捕獲31尾（占比12.6%），在所捕獲魚類中占比最大。前3個航次中漁獲物中均出現(xiàn)具有明顯種群優(yōu)勢的種類，分別為二長棘鯛、黃鰭馬面鲀和短吻鲾，該航次未出現(xiàn)捕獲尾數(shù)占明顯優(yōu)勢的種群。捕獲尾數(shù)達到10尾的有少鱗鱚、白姑魚、金線魚、黑邊天竺魚、勒氏短須石首魚、雙帶天竺鯛、細條天竺魚和杜氏槍烏賊。

2.2 聲學(xué)數(shù)據(jù)分析

此次調(diào)查進行了春季、夏季、秋季和冬季4個航次的聲學(xué)數(shù)據(jù)采集和生物學(xué)拖網(wǎng)取樣，將聲學(xué)數(shù)據(jù)結(jié)合生物學(xué)數(shù)據(jù)，評估了大亞灣灣內(nèi)生物資源量。4個航次中參與評估的總種類有78種魚類和6種頭足類，其中春季航次出現(xiàn)25種魚類和3種頭足類評估生物種，平均生物資源密度為7.38×103 kg/n mile2；夏季航次由于禁漁期剛結(jié)束，與春季航次相比魚的種類明顯增多，共出現(xiàn)53種魚類和6種頭足類評估生物種，平均生物資源密度為8.58×103 kg/n mile2；秋季航次相比夏季航次生物種類和平均資源密度都有所下降，共有24種魚類和3種頭足類評估生物種，平均生物資源密度為6.52×103 kg/n mile2；最后冬季航次中評估生物種有23種魚類和2種頭足類，平均生物資源密度為6.36×103 kg/n mile2。4個航次平均資源密度為7.20×103 kg/n mile2左右，變化趨勢如圖3所示。

調(diào)查海域面積（A）約為23.4 n mile2 ，依據(jù)公式（5）計算大亞灣內(nèi)總資源量，結(jié)果如表2所示。

在4個航次調(diào)查評估結(jié)果中每個航次的評估種類和種類平均生物資源密度都存在一些變化。春季航次單個魚種平均生物資源密度較高的有二長棘鯛（2.37×103 kg/n mile2）、花鰶（1.20×103 kg/n mile2）和短吻鲾（0.68×103 kg/n mile2）。夏季航次評估結(jié)果中黃鰭馬面鲀（3.29×103 kg/n mile2）、海鰻（0.91×103 kg/n mile2）、多齒蛇鯔（0.78×103 kg/n mile2）和短吻鲾（0.63×103 kg/n mile2）平均生物資源密度與其他種類相比較高，但春季航次中平均資源密度較高的種類中短吻鲾?shù)钠骄镔Y源密度（0.63×103 kg/n mile2）變化不大，而花鰶（0.17×103 kg/n mile2）和二長棘鯛（0.37×103 kg/n mile2）波動較大。在秋季航次平均資源密度較高的種類與前2個航次相比也有所變化，其中短吻鲾（2.59×103 kg/n mile2）、中線天竺鯛（0.65×103 kg/n mile2）、杜氏槍烏賊（0.56×103 kg/n mile2）和斑鰶（0.56×103 kg/n mile2）較高。冬季評估結(jié)果中皮氏叫姑魚（0.82×103 kg/n mile2）、少鱗鱚（0.72×103 kg/n mile2）和金線魚（0.64×103 kg/n mile2）的平均生物資源密度較高。每個航次前6種評估種類的資源密度和資源量如表3所示。

3 結(jié)論與討論

此次調(diào)查研究中估算了大亞灣南部海域的生物資源現(xiàn)存狀況，4個季節(jié)的調(diào)查結(jié)果闡述了灣內(nèi)生物資源密度和資源量的季節(jié)性變化情況。大亞灣內(nèi)4個季節(jié)的優(yōu)勢種分別為二長棘鯛、黃鰭馬面鲀、短吻鲾和皮氏叫姑魚，平均生物資源密度為7.20×103 kg/n mile2，與李娜娜等[21]采用聲學(xué)方法評估大亞灣礁區(qū)外的平均生物資源密度（10.0×103 kg/n mile2）相比有所降低，可能與此次調(diào)查海域范圍不同有關(guān)。此次調(diào)查范圍基本涵蓋整個南部海灣而其評估范圍僅為礁區(qū)附近，由于受人工魚礁影響資源密度相比整個海灣會偏大一些。

單個魚種評估結(jié)果證明在大亞灣內(nèi)各個季節(jié)的優(yōu)勢魚種并不是固定不變的，同一魚種在一年中不同季節(jié)的平均資源密度也是不相同的。各個季節(jié)的魚類種數(shù)在評估結(jié)果中也存在明顯差異。在評估種類較少的季節(jié)并不一定灣內(nèi)生物種數(shù)就少，可能是某一種類在某一季節(jié)的資源密度較低，拖網(wǎng)采樣未捕獲到。在評估過程中對于個別魚種的評估結(jié)果不能當(dāng)作大亞灣內(nèi)這一魚種的生物資源密度，如春季的黃斑藍子魚（7.8×103 kg/n mile2）、冬季的中國鯧（5.9×103 kg/n mile2）和黑鯛（4.1×103 kg/n mile2）等與實際情況相差太大。造成這一結(jié)果的原因可能是拖網(wǎng)捕獲的樣本不能代表灣內(nèi)整體，捕獲尾數(shù)較少、個體偏大或單網(wǎng)尾數(shù)過多等情況都不能當(dāng)作灣內(nèi)實際生物分布情況。

聲學(xué)評估方法在數(shù)據(jù)分析和處理過程中雖然經(jīng)過了噪聲剔除和積分分配等步驟，但研究結(jié)果中還可能存在一些不可避免的誤差。調(diào)查采樣船只并非專業(yè)聲學(xué)調(diào)查船，租用船只在聲學(xué)儀器安裝方面還有很多欠缺。調(diào)查時儀器的穩(wěn)定性、吃水深度、換能器角度和所用電壓等方面都會影響采樣數(shù)據(jù)的準確性。參與積分分配的生物一般為非貼地類的游泳魚類，但也有一些情況中非貼地類游泳魚類只參與積分分配而計算結(jié)果不算作灣內(nèi)實際生物現(xiàn)存狀況，雖在拖網(wǎng)捕獲物中出現(xiàn)但捕獲率很低并不能代表灣內(nèi)整體情況的魚類。

參考文獻

[1] 陸忠康.《大亞灣環(huán)境與資源》出版[J].海洋科學(xué)，1990，14（6）：38.

[2] 王友紹，王肇鼎，黃良民.近20年來大亞灣生態(tài)環(huán)境的變化及其發(fā)展趨勢[J].熱帶海洋學(xué)報，2004，23（5）：85-95.

[3] 王肇鼎，練健生，胡建興，等.大亞灣生態(tài)環(huán)境的退化現(xiàn)狀與特征[J].生態(tài)科學(xué)，2003，22（4）：313-320.

[4] 陳丕茂，袁華榮，賈曉平，等.大亞灣楊梅坑人工魚礁區(qū)漁業(yè)資源變動初步研究[J].南方水產(chǎn)科學(xué)，2013，9（5）：100-108.

[5] 李純厚，賈曉平，齊占會，等.大亞灣海洋牧場生物資源養(yǎng)護效果評價[C]//十一五農(nóng)業(yè)環(huán)境研究回顧與展望：第四屆全國農(nóng)業(yè)環(huán)境科學(xué)學(xué)術(shù)研討會論文集.呼和浩特：[出版者不詳]，2011.

[6] 李娜娜，董麗娜，李永振，等.大亞灣海域魚類分類多樣性研究[J].水產(chǎn)學(xué)報，2011，35（6）：863-870.

[7] MISUND O A.Underwater acoustics in marine fisheries and fisheries research[J].Reviews in fish biology and fisheries，1997，7（1）：1-34.

[8] STANLEY D R，WILSON C A.Spatial variation in fish density at three petroleum platforms as measured with dual-beam hydroacoustics[J].Gulf of mexico science，1998，16（1）：73-82.

[9] FABI G，SALA A.An assessment of biomass and diel activity of fish at an artificial reef（Adriatic sea） using a stationary hydroacoustic technique[J].ICES Journal of Marine Science，2002，59：411-420.

[10] HWANG D J，PARK J S，LEE Y W.Estimation of fish school abundance by using an echo sounder in an artificial reef area[J].J Kor Fish Soc，2004，37（3）：249-254.

[11] 譚細暢，常劍波，陶江平，等.三峽庫首魚類分布格局的水聲學(xué)探測評估[J].生態(tài)科學(xué)，2008，27（5）：329-334.

[12] 陶江平，陳永柏，喬曄，等.三峽水庫成庫期間魚類空間分布的水聲學(xué)研究[J].水生態(tài)學(xué)雜志，2008，1（1）：25-33.

[13] 譚細暢，夏立啟，立川賢一，等.東湖放養(yǎng)魚類時空分布的水聲學(xué)研究[J].水生生物學(xué)報，2002，26（6）：585-590.

[14] 張俊，陳國寶，陳作志，等.南沙南部陸架海域漁業(yè)資源聲學(xué)評估[J].南方水產(chǎn)科學(xué)，2015，11（5）：1-10.

[15] 張俊，陳作志，陳國寶，等.南海鳶烏賊水聲學(xué)測量和評估相關(guān)技術(shù)研究[J].南方水產(chǎn)科學(xué)，2014，10（6）：1-11.

[16] 陳國寶，李永振，趙憲勇，等.南海5類重要經(jīng)濟魚類資源聲學(xué)評估[J].海洋學(xué)報，2006，28（2）：128-134.

[17] 李永振，陳國寶，趙憲勇，等.南海北部海域小型非經(jīng)濟魚類資源聲學(xué)評估[J].中國海洋大學(xué)學(xué)報（自然科學(xué)版），2005，35（2）：206-212.

[18] 陳國寶，李永振，趙憲勇，等.南海北部海域重要經(jīng)濟魚類資源聲學(xué)評估[J].中國水產(chǎn)科學(xué)，2005，12（4）：445-451.

[19] 李永振，陳國寶，孫典榮，等.南海區(qū)多魚種聲學(xué)評估工作程序[J].湛江海洋大學(xué)學(xué)報，2003，23（1）：35-40.

[20] 唐啟升，金顯仕，李富國，等.白令海阿留申海盆區(qū)狹鱈當(dāng)年生幼魚數(shù)量分布的調(diào)查研究[J].中國水產(chǎn)科學(xué)，1994，1（1）：37-47.

[21] 李娜娜，陳國寶，于杰，等.大亞灣楊梅坑人工魚礁水域生物資源量聲學(xué)評估[J].水產(chǎn)學(xué)報，2011，35（11）：1640-1649.

[22] 張俊，陳丕茂，房立晨，等.南海柘林灣-南澳島海洋牧場漁業(yè)資源本底聲學(xué)評估[J].水產(chǎn)學(xué)報，2015，39（8）：1187-1198.

[23] 趙憲勇，陳毓楨，李顯森，等.多種類海洋漁業(yè)資源聲學(xué)評估技術(shù)與方法[C]∥我國專屬經(jīng)濟區(qū)和大陸架勘測研究論文集.北京：海洋出版社，2002：341-353.

[24] 屈泰春，黃洪亮，湯勇，等.漁業(yè)聲學(xué)數(shù)據(jù)后處理中噪聲剔除的研究進展[J].漁業(yè)信息與戰(zhàn)略，2013，28（3）：208-213.

[25] FOOTE K G.Fish target strengths for use in echo-integrator surveys[J].Journal of the acoustical society of America，1987，82（3）：981-987.