陳霄航 連玉喜 黃 耿 項(xiàng) 濤 張?zhí)昧?劉家壽 葉少文 李鐘杰
(1. 中國科學(xué)院水生生物研究所淡水生態(tài)與生物技術(shù)國家重點(diǎn)實(shí)驗(yàn)室, 武漢 430072; 2. 中國科學(xué)院大學(xué), 北京 100049;3. 安慶師范大學(xué)水生生物保護(hù)與水生態(tài)修復(fù)安徽省高校工程技術(shù)研究中心, 安慶 246133)
水聲學(xué)是近代聲學(xué)的一個(gè)重要分支, 隨著水聲學(xué)儀器的更新?lián)Q代和計(jì)算機(jī)技術(shù)的飛速發(fā)展, 漁業(yè)資源聲學(xué)評(píng)估法廣泛應(yīng)用于魚類資源調(diào)查及評(píng)估[1]。與傳統(tǒng)方法相比, 水聲學(xué)調(diào)查魚類資源具有快速高效、覆蓋范圍廣、不破壞生境等優(yōu)點(diǎn), 在此基礎(chǔ)上能夠得到豐富的魚類行為、分布、密度等信息。在漁業(yè)資源的水聲學(xué)評(píng)估中, 魚體目標(biāo)強(qiáng)度(Target Strength, TS)是衡量魚體對(duì)聲波反射能力的一個(gè)物理量, 它從回聲強(qiáng)度的角度描述目標(biāo)的聲學(xué)特性,是將回聲積分值轉(zhuǎn)換為魚體大小和魚類資源量估值的關(guān)鍵參數(shù)[2]。用回聲積分法計(jì)算魚類密度時(shí),一般是將單個(gè)魚體的平均目標(biāo)強(qiáng)度代表成魚類的總目標(biāo)強(qiáng)度, 通過計(jì)算單位水體內(nèi)魚類平均體積散射強(qiáng)度(Volume Backscattering Strength, SV)與單個(gè)魚體目標(biāo)強(qiáng)度的比值來估算魚類密度[3]。因此, 魚類目標(biāo)強(qiáng)度的精確測(cè)定對(duì)于漁業(yè)資源聲學(xué)評(píng)估的可靠性影響重大[4]。
影響魚類目標(biāo)強(qiáng)度的主要因素包括魚體形狀、是否有鰾、鰾的大小、魚體相對(duì)入射聲波的姿態(tài)傾角, 以及探測(cè)聲波的頻率等[5—7]。實(shí)驗(yàn)測(cè)定法是研究魚類目標(biāo)強(qiáng)度的一個(gè)重要途徑, 使用聲學(xué)實(shí)驗(yàn)裝置對(duì)受控魚類進(jìn)行實(shí)驗(yàn)測(cè)量, 或利用船載聲學(xué)儀器對(duì)自然狀態(tài)下的魚類進(jìn)行現(xiàn)場(chǎng)測(cè)定并通過船后的拖網(wǎng)取樣比對(duì)[8]。采用實(shí)驗(yàn)法測(cè)定魚體目標(biāo)強(qiáng)度通常受到場(chǎng)地、可測(cè)量頻率以及測(cè)試對(duì)象差異等諸多條件的限制, 對(duì)一般性規(guī)律的把握較為困難, 但具有針對(duì)性強(qiáng)、精度高的優(yōu)點(diǎn), 在特定種類的漁業(yè)資源評(píng)估中具有顯著的應(yīng)用價(jià)值。
由于魚體目標(biāo)強(qiáng)度測(cè)定對(duì)實(shí)驗(yàn)條件要求較高,繩系實(shí)驗(yàn)法早在70年代便有記載[9], 多使用死魚或麻醉魚。迄今為止國際上僅對(duì)幾十種魚類的目標(biāo)強(qiáng)度進(jìn)行過研究, 遠(yuǎn)遠(yuǎn)不能滿足漁業(yè)資源聲學(xué)評(píng)估的需求。我國自1984年將水聲學(xué)方法引入漁業(yè)資源評(píng)估以來, 關(guān)于魚類目標(biāo)強(qiáng)度測(cè)量研究的報(bào)道較少[7,10—12]。目前國內(nèi)仍主要使用國際上已有的目標(biāo)強(qiáng)度與體長關(guān)系式, 考慮到魚類種類組成和聲學(xué)特性的差異, 會(huì)導(dǎo)致資源量估算結(jié)果出現(xiàn)一定偏差[13]。通過實(shí)驗(yàn)方法測(cè)定在我國廣泛分布和具有重要經(jīng)濟(jì)價(jià)值魚類的目標(biāo)強(qiáng)度, 建立它們的水聲目標(biāo)強(qiáng)度與個(gè)體大小關(guān)系式, 將有助于提高魚類種群結(jié)構(gòu)和資源量評(píng)估的準(zhǔn)確性和可靠性。
鯉(Cyprinus carpio)和鳙(Aristichthys nobilis)是廣泛分布的鯉科經(jīng)濟(jì)魚類, 也是我國養(yǎng)殖的主要品種, 據(jù)2016年中國漁業(yè)統(tǒng)計(jì)年鑒統(tǒng)計(jì), 鯉、鳙約占全國淡水養(yǎng)殖主要魚類產(chǎn)量22%[14]。鯉科魚類更是占我國淡水魚類總數(shù)目的一半以上[15]。鯉、鳙和大多數(shù)鯉科魚類有著相似的紡錘體體形, 從棲息空間上看, 鯉生活于水體的下層和底層, 鳙生活于水體的中上層, 因而開展這兩種魚類的水聲目標(biāo)強(qiáng)度測(cè)定研究具有一定的典型性和實(shí)用性。有鑒于此, 本文采用繩系控制實(shí)驗(yàn)法和Simrad EY60型分裂波束回聲探測(cè)儀, 對(duì)實(shí)驗(yàn)魚鯉和鳙的目標(biāo)強(qiáng)度進(jìn)行單體測(cè)定, 建立目標(biāo)強(qiáng)度與體長的關(guān)系式并對(duì)實(shí)際應(yīng)用效果進(jìn)行評(píng)價(jià), 旨在探討魚類目標(biāo)強(qiáng)度測(cè)量的實(shí)驗(yàn)方法, 完善我國以鯉科魚類為主的內(nèi)陸水聲學(xué)評(píng)估技術(shù), 為漁業(yè)資源的可持續(xù)利用和保護(hù)提供科學(xué)依據(jù)。
魚體目標(biāo)強(qiáng)度測(cè)定實(shí)驗(yàn)地點(diǎn)選在湖北省宜昌市隔河巖水庫一處僻靜的庫灣(N 30.404579°, E 110.980000°), 該水域開闊, 平均水深約40 m, 水面平靜, 無航運(yùn)和漁業(yè)活動(dòng)干擾, 水中氣泡較少。實(shí)驗(yàn)平臺(tái)離岸1 km, 所在位置的水深47 m。經(jīng)探測(cè),半徑60 m內(nèi)水下無強(qiáng)聲波反射體, 實(shí)驗(yàn)水域近似滿足聲學(xué)自由場(chǎng)條件。
實(shí)驗(yàn)裝置結(jié)構(gòu)如圖 1所示。實(shí)驗(yàn)平臺(tái)搭建在網(wǎng)箱上(5 m × 5 m × 10 m, 網(wǎng)目2a = 1 cm), 用錨固定網(wǎng)箱, 提高其穩(wěn)定性。測(cè)量實(shí)驗(yàn)在網(wǎng)箱內(nèi)進(jìn)行, 可避免其他魚類闖入干擾。通過不銹鋼輔助支臂, 將Simrad EY60型分裂波束回聲探測(cè)儀的換能器(頻率120 kHz)固定在網(wǎng)箱中央位置、水面以下0.5 m處。用直徑0.6 mm尼龍絲(釣魚線)將活體實(shí)驗(yàn)魚用長約2 cm魚線一端系住下頜, 另一端連在換能器下方9 m處的塑料環(huán)上, 任其圍繞塑料環(huán)游動(dòng)。塑料環(huán)下方0.5 m處系上5 kg鉛錘, 保證整個(gè)裝置垂直。

圖 1 魚類目標(biāo)強(qiáng)度測(cè)定實(shí)驗(yàn)平臺(tái)示意圖Fig. 1 Experimental platform for measuring fish target strength
材料魚鯉和鳙均就地取自隔河巖水庫, 實(shí)驗(yàn)前在網(wǎng)箱中暫養(yǎng)。根據(jù)體長范圍, 實(shí)驗(yàn)時(shí)按照4個(gè)體長組選取健康無損傷個(gè)體, 用于目標(biāo)強(qiáng)度的單體測(cè)定。實(shí)驗(yàn)魚鯉26尾, 全長19.9—29.6 cm, 標(biāo)準(zhǔn)長16.6—24.8 cm, 體重174.5—460.9 g, 相應(yīng)4個(gè)體長組的體征參數(shù)如表 1所示。實(shí)驗(yàn)魚鳙26尾, 全長49.2—74.2 cm, 標(biāo)準(zhǔn)長40.5—63.2 cm, 體重1200—4700 g,相應(yīng)4個(gè)體長組的體征參數(shù)如表 2所示。

表 1 實(shí)驗(yàn)魚鯉的規(guī)格大小(平均值±標(biāo)準(zhǔn)差)Tab. 1 Body size of experimental C. carpio (Mean ± SD)
水聲學(xué)測(cè)量水聲學(xué)測(cè)量實(shí)驗(yàn)時(shí)間為2017年5月1—20日, 期間平均水溫25℃、電導(dǎo)率223.2 μS/cm、溶解氧8.14 mg/L、pH 8.38、流速小于0.1 m/s。在實(shí)驗(yàn)平臺(tái)上, 使用Simrad EY60型分裂波束回聲探測(cè)儀(頻率120 kHz)測(cè)量魚體目標(biāo)強(qiáng)度測(cè)量。換能器-3 dB波束為7°×7°, 通過不銹鋼支臂固定在網(wǎng)箱中央位置, 入水深度0.5 m, 垂直向下。在實(shí)際測(cè)量前, 采用Simrad公司提供的直徑為23 mm標(biāo)準(zhǔn)銅球, 按標(biāo)準(zhǔn)方法[16]對(duì)回聲探測(cè)儀參數(shù)設(shè)置進(jìn)行校正。測(cè)量時(shí)功率設(shè)定為50 w, 脈沖發(fā)射頻率10 ping/s, 脈沖寬度64ms。筆記本電腦與回聲探測(cè)儀連接, 實(shí)時(shí)顯示和存儲(chǔ)數(shù)據(jù)。
將每種實(shí)驗(yàn)魚按照體長組分類, 從小到大逐尾逐次用尼龍線固定連接在換能器下方的塑料環(huán)上(圖 1), 使得魚體圍繞塑料環(huán)保持正常姿態(tài)游動(dòng), 發(fā)射脈沖并使魚體使其處于換能器波束中央, 每尾魚記錄5min。
數(shù)據(jù)分析使用Echoview軟件(version 5.2)對(duì)采集得到的水聲學(xué)數(shù)據(jù)進(jìn)行分析。根據(jù)魚體目標(biāo)信號(hào)所處水深, 界定回聲映像圖中的分析區(qū)域,見圖中方框區(qū)域內(nèi), 該區(qū)域以外不予分析(圖 2)。對(duì)分析區(qū)域進(jìn)行除噪之后, 用單體回波探測(cè)(Single Echo Detection, SED)方法分析單個(gè)魚體信號(hào), 主要參數(shù)設(shè)置如表 3所示。

表 2 實(shí)驗(yàn)魚鳙的規(guī)格大小(平均值±標(biāo)準(zhǔn)差)Tab. 2 Body size of experimental A. nobilis (Mean ± SD)

圖 2 實(shí)驗(yàn)魚體目標(biāo)強(qiáng)度信號(hào)的回聲映像圖Fig. 2 Target strength echogram of experimental fishes
通過Echoview得到的是單體魚信號(hào)TS值(單位為dB), 故在計(jì)算單體魚平均TS和各體長組魚的平均TS時(shí), 需要先將每個(gè)回波數(shù)據(jù)的TS轉(zhuǎn)換為反向散射截面σ(單位為m2), 對(duì)聲學(xué)截面進(jìn)行平均后, 再將其轉(zhuǎn)換為TS。相關(guān)計(jì)算公式如下:

TS與魚體全長(TL, 單位為cm)之間回歸關(guān)系依據(jù)公式(4)進(jìn)行擬合[17]:

為了評(píng)價(jià)所得魚類目標(biāo)強(qiáng)度(TS)與全長(TL)關(guān)系式的實(shí)際應(yīng)用效果, 選擇一個(gè)面積較小且有較詳細(xì)魚產(chǎn)量記錄的封閉水域(鹽龍湖)作為試驗(yàn)水體,使用本文TS-TL關(guān)系式和其他幾個(gè)常用TS-TL關(guān)系式, 對(duì)該湖的魚類資源水聲學(xué)調(diào)查數(shù)據(jù)進(jìn)行分析,估算全湖魚類資源量并與捕撈產(chǎn)量數(shù)據(jù)作比較, 在此基礎(chǔ)上評(píng)價(jià)不同關(guān)系式在魚類資源評(píng)估中的可靠性和適用性。
鹽龍湖位于江蘇省鹽城市龍崗鎮(zhèn)(N 33.334568°,E 120.022081°, 圖 3), 水域面積2.228 km2, 平均水深4.2 m, 無沉水植被分布。于2017年2月14日在該湖開展魚類資源水聲學(xué)調(diào)查, 將Simrad EY60(頻率120 kHz)換能器固定于船體右側(cè)、水面以下0.5 m處, 水聲學(xué)探測(cè)從鹽龍湖東南角開始, 至西北端結(jié)束, 全湖調(diào)查路線如圖 3所示, 總航程約為6 km, 水聲學(xué)探測(cè)時(shí)航速控制為5-6 km/h, 探測(cè)覆蓋率為4.0, 接近較為可信評(píng)估的覆蓋率。
圖 4所示為實(shí)驗(yàn)魚鯉和鳙的體重與體長(全長、標(biāo)準(zhǔn)長)關(guān)系擬合曲線, 符合冪函數(shù)(W= a ×Lb)分布。運(yùn)用最小二乘法得到鯉和鳙的體重與全長關(guān)系式的b值分別為3.178和3.227, 經(jīng)t檢驗(yàn)均與3無顯著差異 (P>0.05), 說明兩種魚的生長屬于等速生長。

表 3 分裂波束單體回波探測(cè)方法的參數(shù)設(shè)置Tab. 3 Parameter used in a single echo detection by split beam
圖 5為鯉和鳙目標(biāo)強(qiáng)度與全長關(guān)系擬合曲線,實(shí)驗(yàn)所測(cè)26尾鯉和26尾鳙的TS值分別為-59.35—-45.20 dB和-34.22— -17.49 dB, 回歸得到的兩種魚目標(biāo)強(qiáng)度(TS, dB)與全長(TL, cm)關(guān)系式分別為:
鯉:TS= 29.84×lgTL- 95.23 (n= 26,R2= 0.74,P< 0.01);
鳙:TS= 35.88×lgTL- 90.33 (n= 26,R2= 0.83,P< 0.01)。

圖 4 實(shí)驗(yàn)魚(a)鯉和(b)鳙的體重與體長關(guān)系Fig. 4 Relationships between body weight and length for C. carpio (a) and A. nobilis (b)

圖 3 鹽龍湖水聲學(xué)調(diào)查路線(圖中箭頭示調(diào)查走航方向)Fig. 3 Hydroacoustic sampling route in Lake Yanlong
使用本研究鳙目標(biāo)強(qiáng)度與全長關(guān)系式、Love[18]、Foote[16]和任玉芹等[12]的TS-TL關(guān)系式, 對(duì)2017年2月14日采集的鹽龍湖水聲學(xué)調(diào)查數(shù)據(jù)進(jìn)行分析, 分別得到該湖魚類資源總量的估值為89360、36630、493650和8244470 kg, 折合生物量估值為40.1、16.4、221.6和3700.4 g/m2(表 4); 另一方面, 得到2017年3月該湖“趕、攔、刺、張”聯(lián)合漁法的捕撈產(chǎn)量為73610 kg (即33.0 g/m2), 其中鳙、鰱約占70%。綜合這兩方面信息, 計(jì)算捕撈量與4個(gè)TS-TL關(guān)系式所得魚類生物量的比值(C/B)。
對(duì)于本研究TS-TL關(guān)系式, C/B值為82.4%, 較為合理; 對(duì)于Love[18]的TS-TL關(guān)系式, C/B值超過了100%, 說明生物量估值偏低; 對(duì)于Foote[16]的TSTL關(guān)系式, C/B值較低(14.9%), 說明生物量估值有所偏高; 對(duì)于任玉芹等[12]的TS-TL關(guān)系式, C/B值僅為0.9%, 說明生物量估值太高。

圖 5 實(shí)驗(yàn)魚(a)鯉和(b)鳙的目標(biāo)強(qiáng)度與全長關(guān)系Fig. 5 Relationships between target strength and total length forC. carpio (a) and A. nobilis (b)

表 4 采用4種不同的目標(biāo)強(qiáng)度與全長關(guān)系式評(píng)估鹽龍湖魚類資源量Tab. 4 Assessment of fish resource in Lake Yanlong using four different TS-TL equations
目前在我國內(nèi)陸水域魚類資源聲學(xué)評(píng)估中, 在將目標(biāo)強(qiáng)度轉(zhuǎn)換為魚體長度時(shí)多借鑒國外的經(jīng)驗(yàn)公式, 使用較多的有Love[18]建立的多種魚類全長和目標(biāo)強(qiáng)度的回歸關(guān)系式、Foote[16]建立的有鰾魚類全長和目標(biāo)強(qiáng)度的回歸關(guān)系式等, 然而, 由于魚類種類組成及其聲學(xué)特性的差異, 可能導(dǎo)致魚類資源評(píng)估的偏差[19]。本文選擇我國廣泛分布的鯉科經(jīng)濟(jì)魚類鯉和鳙作為研究對(duì)象, 具有很強(qiáng)的針對(duì)性和代表性; 在環(huán)境干擾小的網(wǎng)箱水體里開展繩系控制實(shí)驗(yàn), 現(xiàn)場(chǎng)測(cè)定單體實(shí)驗(yàn)魚游動(dòng)狀態(tài)下的目標(biāo)強(qiáng)度,建立了這2種魚的目標(biāo)強(qiáng)度與個(gè)體大小關(guān)系式, 并取得較好的實(shí)際應(yīng)用效果, 有助于提高以鯉科魚類為主體的漁業(yè)資源評(píng)估的準(zhǔn)確性和可靠性。
在具體實(shí)驗(yàn)過程中, 魚體形狀對(duì)目標(biāo)強(qiáng)度的影響在于相同長度的魚體具有不同的聲學(xué)散射截面[20], 本研究實(shí)驗(yàn)魚鯉和鳙為典型的紡錘體體形,與鰱、草魚、青魚、鯽等鯉科常見經(jīng)濟(jì)魚類的體形相近; 另一方面, 鰾是影響魚類目標(biāo)強(qiáng)度的重要結(jié)構(gòu)因素, 占據(jù)魚體散射能力的90%—95%[21], 鯉、鳙和其他鯉科魚類均為有鰾魚, 鯉主要棲息于水體下層和底層, 鳙棲息于水體的中上層, 它們?cè)诳臻g生態(tài)位上各具代表性。因此, 以鯉、鳙以及近似體形和具鰾結(jié)構(gòu)魚類為主的水域, 在尚未建立各種魚綜合TS-TL回歸關(guān)系的狀況下, 可以借鑒本研究的TS-TL關(guān)系式, 用于魚類種群結(jié)構(gòu)和資源量的聲學(xué)評(píng)估。
魚體目標(biāo)強(qiáng)度與魚體相對(duì)入射聲波的姿態(tài)傾角也存在著密切聯(lián)系, 魚體姿態(tài)傾角的差異可引起聲學(xué)散射截面的改變, 進(jìn)而引起目標(biāo)強(qiáng)度的變化。由姿態(tài)傾角引起的目標(biāo)強(qiáng)度差異可高達(dá)30 dB[22]。本研究用0.6 mm的釣魚線將實(shí)驗(yàn)魚下頜系住, 并未對(duì)魚體組織產(chǎn)生任何損傷, 圍繞水下固定的塑料環(huán)自由活動(dòng), 因此較為接近自然游動(dòng)姿態(tài), 在此過程中通過分裂波束回聲探測(cè)儀不斷錄入數(shù)據(jù), 達(dá)到減小實(shí)驗(yàn)誤差的目的。
水體背景噪音是影響魚類目標(biāo)強(qiáng)度測(cè)定的主要外界因素[1]。背景噪聲一般來源于聲學(xué)設(shè)備、調(diào)查船只、非目標(biāo)生物、環(huán)境噪聲等, 在目標(biāo)強(qiáng)度測(cè)量過程中會(huì)降低回聲信號(hào)的信噪比, 增加數(shù)據(jù)提取的難度, 當(dāng)背景噪音增大到一定程度時(shí), 魚體的反射信號(hào)就很難被檢測(cè)出來。本文實(shí)驗(yàn)選擇隔河巖水庫的一個(gè)僻靜庫灣, 水環(huán)境總體穩(wěn)定, 周邊無航運(yùn)、漁業(yè)活動(dòng)等人為干擾, 同時(shí)通過網(wǎng)箱將實(shí)驗(yàn)魚與其他魚類隔離開, 這些措施可有效地降低水體背景噪音, 提高目標(biāo)魚類聲學(xué)數(shù)據(jù)測(cè)量的精準(zhǔn)度。
回聲探測(cè)儀的聲波頻率是影響魚類目標(biāo)強(qiáng)度測(cè)定的另一重要因素, 當(dāng)外源聲波頻率很低, 接近魚鰾的固有頻率時(shí), 魚鰾會(huì)產(chǎn)生低頻共振, 引起反射強(qiáng)度的變化[23]。目前漁業(yè)資源評(píng)估中常用的換能器頻率為38—200 kHz, 本文實(shí)驗(yàn)采用的換能器頻率適中(120 kHz), 測(cè)試目標(biāo)強(qiáng)度受低頻共振的影響較小, 具有較高的穩(wěn)定性。
為了評(píng)價(jià)所獲魚類目標(biāo)強(qiáng)度與體長關(guān)系式的應(yīng)用效果, 我們選擇一個(gè)小型封閉水體(鹽龍湖)大捕撈之前的回聲探測(cè)數(shù)據(jù)作為分析實(shí)例, 估算全湖魚類資源總量, 并與聯(lián)合漁法大捕撈的產(chǎn)量進(jìn)行比較。根據(jù)捕撈經(jīng)驗(yàn), 使用聯(lián)合漁法捕撈的捕撈量約占資源量60%—70%, 結(jié)果表明通過本研究TS-TL關(guān)系式得到的魚類總產(chǎn)量和總生物量比值(82.4%)較為合理, 基本符合“趕、攔、刺、張”聯(lián)合漁法的捕撈選擇性(捕撈對(duì)象主要為中上層的鳙、鰱)和高效性(尤其是在該湖形狀簡(jiǎn)單且無沉水植被的情況下), 因而在一定程度上反映了本研究目標(biāo)強(qiáng)度測(cè)定及其與體長關(guān)系式的可靠性和合理性。
本研究存在的不足之處在于, 由于實(shí)驗(yàn)期間活魚來源的限制, 材料魚鯉和鳙的體長范圍相對(duì)較窄,今后仍需補(bǔ)充更多體長組魚類來完善實(shí)驗(yàn)方案, 在本文實(shí)驗(yàn)條件和結(jié)果的基礎(chǔ)上進(jìn)一步描述目標(biāo)強(qiáng)度隨體長的變化特征; 同時(shí)亟需開展更多常見魚類(如鰱、草魚、鯽、團(tuán)頭魴、鱖、中華鱘等)的目標(biāo)強(qiáng)度測(cè)量實(shí)驗(yàn), 確定不同形態(tài)和生態(tài)類型魚類的目標(biāo)強(qiáng)度與體長關(guān)系式, 以滿足我國內(nèi)陸水域魚類生態(tài)學(xué)研究、漁業(yè)資源可持續(xù)利用和瀕危魚類保護(hù)等領(lǐng)域?qū)μ岣咚晫W(xué)評(píng)估精準(zhǔn)度的不斷要求[24]。