基于體長(zhǎng)標(biāo)準(zhǔn)的大長(zhǎng)山人工魚礁區(qū)許氏平鮋資源評(píng)估

薛國進(jìn) 尹增強(qiáng) 程前 魏藝鳴 李九奇

摘 要：為了評(píng)價(jià)大長(zhǎng)山人工魚礁區(qū)許氏平鮋（Sebastes schlegelii）資源利用狀態(tài)，為當(dāng)?shù)卦S氏平鮋資源的合理開發(fā)與管理提供科學(xué)依據(jù)，采用基于體長(zhǎng)的貝葉斯評(píng)估方法（LBB）對(duì)2011-2015年大長(zhǎng)山人工魚礁區(qū)許氏平鮋資源進(jìn)行了評(píng)估。研究結(jié)果表明大長(zhǎng)山人工魚礁區(qū)許氏平鮋最優(yōu)體長(zhǎng)Lopt為31 cm，最佳開捕體長(zhǎng)Lc-opt為28 cm，平均M/K=1.5、F/K=4.82、Z/K=6.34、F/M=3.24，許氏平鮋最優(yōu)現(xiàn)存資源量與原始資源量之比B/B0（F=M，Lc=Lc-opt）=0.366、實(shí)際年平均B/B0= 0.038 6，許氏平鮋最優(yōu)的單位補(bǔ)充量漁獲量Y/R'（F=M，Lc=Lc-opt）=0.045 9 kg，實(shí)際年平均Y/R'=0.003 44 kg。

關(guān)鍵詞：大長(zhǎng)山;許氏平鮋（Sebastes schlegelii）;體長(zhǎng)分布;資源評(píng)估;LBB

許氏平鮋（Sebastes schlegelii）又名黑鮶，是近海冷水性底層魚類，以小魚、甲殼類、頭足類等為食，常棲息于巖礁底質(zhì)的海區(qū)，是我國北方人工魚礁海域的主要增殖對(duì)象[1]。作為我國北方沿海重要的經(jīng)濟(jì)魚種，其資源狀況直接影響當(dāng)?shù)貪O民的生計(jì)問題，由于近年高強(qiáng)度的過度捕撈使其資源逐漸衰退，漁民收入逐年下降，針對(duì)這些問題相關(guān)單位進(jìn)行了一系列許氏平鮋的研究，但主要集中在魚類生物學(xué)以及魚類增養(yǎng)殖方面[2]，對(duì)許氏平鮋種群數(shù)量動(dòng)態(tài)、資源利用狀態(tài)的研究相對(duì)較少[3]。而我國由于缺乏經(jīng)典漁業(yè)資源評(píng)估方法所必需的各類數(shù)據(jù)[4]，如精準(zhǔn)的歷年漁業(yè)統(tǒng)計(jì)數(shù)據(jù)、歷年資源調(diào)查數(shù)據(jù)、歷年漁獲物各齡魚年齡與體長(zhǎng)組成等，致使經(jīng)典漁業(yè)資源評(píng)估研究難以進(jìn)行，造成我國漁業(yè)資源管理存在高成本、無效率和無效性等現(xiàn)象。為了克服漁業(yè)資源評(píng)估數(shù)據(jù)缺乏的問題，漁業(yè)發(fā)達(dá)國家的科學(xué)家開發(fā)了在數(shù)據(jù)有限情況下的漁業(yè)資源評(píng)估方法，如僅需歷年漁獲量數(shù)據(jù)的cmsy方法[5]、僅需歷年漁獲物體長(zhǎng)數(shù)據(jù)的LBB[6]方法等，這些方法開辟了在有限漁業(yè)數(shù)據(jù)情況下進(jìn)行漁業(yè)資源評(píng)估的途徑，對(duì)像我國這種嚴(yán)重缺乏經(jīng)典漁業(yè)資源評(píng)估方法所需的必要數(shù)據(jù)的國家尤為重要。本文基于2011-2015年大長(zhǎng)山人工魚礁區(qū)許氏平鮋體長(zhǎng)生物學(xué)測(cè)定數(shù)據(jù)采用LBB方法對(duì)該資源利用狀態(tài)進(jìn)行評(píng)價(jià)，獲得漁業(yè)資源種群變動(dòng)的關(guān)鍵參數(shù)，以期為該資源的合理開發(fā)與管理提供科學(xué)依據(jù)。

1 材料與方法

1.1 調(diào)查時(shí)間及站位設(shè)置

國家海洋科技支撐計(jì)劃項(xiàng)目組于2011-2015年連續(xù)五年對(duì)大長(zhǎng)山人工魚礁區(qū)進(jìn)行許氏平鮋資源調(diào)查，調(diào)查站位如圖1，調(diào)查方法參照《海洋漁業(yè)資源調(diào)查規(guī)范》[7]，并對(duì)許氏平鮋進(jìn)行體長(zhǎng)測(cè)定，統(tǒng)計(jì)結(jié)果如表1所示。

2 研究結(jié)果

2.1 2011-2015年大長(zhǎng)山許氏平鮋各重要模型參數(shù)估計(jì)值

利用前面模型對(duì)2011-2015年漁獲物體長(zhǎng)分布頻率進(jìn)行擬合，所有模型計(jì)算與繪圖均采用R語言編程完成[16]，模型計(jì)算的迭代次數(shù)為20 000次，估計(jì)結(jié)果如表2所示：

由表2可知，大長(zhǎng)山許氏平鮋的漸近體長(zhǎng)L∞為46.2 cm;未開發(fā)狀態(tài)下其資源生物量最大時(shí)對(duì)應(yīng)的魚類體長(zhǎng)Lopt為31 cm;最優(yōu)開捕體長(zhǎng)Lc-opt為28 cm;自然死亡系數(shù)M與魚類生長(zhǎng)參數(shù)K之比M/K為1.5;捕撈死亡系數(shù)F與魚類生長(zhǎng)參數(shù)K之比F/K為4.82;總死亡系數(shù)Z與魚類生長(zhǎng)參數(shù)K之比Z/K為6.34;捕撈死亡系數(shù)F與自然死亡系數(shù)M之比F/M為3.24;當(dāng)捕撈死亡系數(shù)F等于自然死亡系數(shù)M且開捕體長(zhǎng)Lc等于Lc-opt時(shí)，大長(zhǎng)山許氏平鮋的資源生物量B與原始資源量B0之比Bmsy/B0為0.366;2015年評(píng)估的實(shí)際資源量與原始資源量B0之比B2015/B0為0.038 6;當(dāng)捕撈死亡系數(shù)F等于自然死亡系數(shù)M且開捕體長(zhǎng)Lc等于Lc-opt時(shí)，單位補(bǔ)充量漁獲量最優(yōu)值MSY/R為0.045 9 kg;而2015年實(shí)際平均單位補(bǔ)充量漁獲量Y2015/R為0.003 44 kg。

2.2 大長(zhǎng)山許氏平鮋各重要模型參數(shù)的時(shí)間序列分布及資源利用評(píng)價(jià)圖

歷年漁獲物體長(zhǎng)分布與最優(yōu)體長(zhǎng)比較與模型計(jì)算得到的各重要參數(shù)時(shí)間序列分布情況與繪圖表示，并分別列于圖2和圖3。

圖2顯示了大長(zhǎng)山許氏平鮋2011-2015年漁獲物體長(zhǎng)分布與最優(yōu)體長(zhǎng)Lopt比較;圖3給出了歷年漁獲物平均體長(zhǎng)Lmean與最優(yōu)體長(zhǎng)Lopt和最優(yōu)開捕體長(zhǎng)Lc-opt對(duì)比、捕撈壓力指標(biāo)F/M與最優(yōu)值（F=M）對(duì)比和資源消耗成度指標(biāo)B/B0與最優(yōu)值Bmsy/B0對(duì)比的變化趨勢(shì)，其中0.5 Bmsy為FAO規(guī)定的資源量不可逾越的警戒線，越過資源量警戒線后補(bǔ)充量將線性下降。

3 討論與分析

在有限漁業(yè)數(shù)據(jù)情況下（僅歷年漁獲物體長(zhǎng)測(cè)定數(shù)據(jù)），LBB方法通過理論計(jì)算體長(zhǎng)頻率對(duì)實(shí)測(cè)體長(zhǎng)頻率的Bayes擬合，主要提供了對(duì)最優(yōu)體長(zhǎng)Lopt、最優(yōu)開捕體長(zhǎng)Lc-opt、捕撈壓力指標(biāo)F/M、資源消耗程度指標(biāo)B/B0、最優(yōu)相對(duì)資源量指標(biāo)Bmsy/B0、單位補(bǔ)充量產(chǎn)量指標(biāo)Y2015/R和最優(yōu)單位補(bǔ)充量產(chǎn)量指標(biāo)MSY/R等進(jìn)行評(píng)估（表2）?？梢詮囊陨系脑u(píng)估指標(biāo)對(duì)大連大長(zhǎng)山許氏平鮋資源利用狀態(tài)進(jìn)行評(píng)價(jià)，首先從2011-2015年的實(shí)際開捕體長(zhǎng)Lc與其最優(yōu)開捕體長(zhǎng)Lc-opt及2011-2015年漁獲物平均體長(zhǎng)與其最優(yōu)體長(zhǎng)Lopt的對(duì)比均可以看出，各年實(shí)際開捕體長(zhǎng)和漁獲物平均體長(zhǎng)均遠(yuǎn)遠(yuǎn)小于這兩個(gè)最優(yōu)體長(zhǎng)標(biāo)準(zhǔn)，特別是2012年實(shí)際開捕體長(zhǎng)不及最優(yōu)開捕體長(zhǎng)的20%（圖2和圖3），在魚類高速生長(zhǎng)期開始捕撈，直接導(dǎo)致漁獲物小型化和漁業(yè)低效率運(yùn)行;從捕撈壓力指標(biāo)F/M來看，2011-2015年F/M平均值達(dá)到3.24（最優(yōu)值F/M=1）（表2和圖3），僅2012和2015年接近最優(yōu)值（F=M）（圖3），捕撈強(qiáng)度過大，同樣導(dǎo)致該漁業(yè)低效率;從資源消耗程度指標(biāo)B/B0對(duì)照最優(yōu)相對(duì)資源量指標(biāo)Bmsy/B0 來看（圖3），僅2012年資源消耗程度指標(biāo)B/B0達(dá)到最優(yōu)值，其它4年該指標(biāo)均處于警戒線之下，資源處于極度危險(xiǎn)狀態(tài)，較為樂觀的是2015年有逐漸恢復(fù)的趨勢(shì);最后從單位補(bǔ)充量漁獲量Y2015/R與其最優(yōu)值MSY/R對(duì)比來看，2015年實(shí)際的Y2015/R僅為最優(yōu)值MSY/R的7.5%（表2），當(dāng)前資源利用策略已造成資源的嚴(yán)重浪費(fèi)。

總之，無論從以上哪個(gè)指標(biāo)來評(píng)價(jià)，大連大長(zhǎng)山許氏平鮋資源利用均處于嚴(yán)重過度捕撈狀態(tài)，特別是從資源消耗程度指標(biāo)B/B0來看，其資源利用情況已處于警戒線之下，資源處于極度危險(xiǎn)狀態(tài)，需及時(shí)進(jìn)行人為干預(yù)，建議專門設(shè)立幼魚保護(hù)區(qū)和禁漁期，開捕體長(zhǎng)嚴(yán)格控制在28 cm以上，捕撈死亡系數(shù)嚴(yán)格控制在0.2以內(nèi)。

參考文獻(xiàn)：

[1] 唐衍力，白懷宇，盛化香，等. 海州灣前三島魚礁區(qū)許氏平鮋的分布及YPUE與近礁距離的關(guān)系[J]. 中國海洋大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，46（11）：151-157.

[2] 李蘭敏， 溫海深， 張亞晨. 卵胎生許氏平鮋個(gè)體生殖力與生物學(xué)指標(biāo)的相關(guān)分析[J]. 現(xiàn)代農(nóng)業(yè)科技， 2015（22）：267-269，279.

[3] 尹增強(qiáng)，許傳才，陳勇，等. 獐子島海域許氏平鮋生長(zhǎng)特征與資源評(píng)估[J].農(nóng)業(yè)經(jīng)濟(jì)與科技，2015，26（10）：71-72

[4] 耿喆，朱江峰，夏萌，等. 數(shù)據(jù)缺乏條件下的漁業(yè)資源評(píng)估方法研究進(jìn)展[J]. 海洋湖沼通報(bào)，2018，40（05）：132-139.

[5] Froese R， Demirel N， Coro G ， et al. Estimating fisheries reference points from catch and resilience[J]. Fish and Fisheries， 2016，18（3）：506–526.

[6] Froese R， Winker H， Coro G， et al. A new approach for estimating stock status from length frequency data[J]. ICES Journal of Marine Science，2018，75（6）： 2004-2015.

[7] 金顯仕，李顯森，趙憲勇，等.海洋漁業(yè)資源調(diào)查規(guī)范：SC/T 9403-2012[S].北京：中國農(nóng)業(yè)出版社，2013：12，25.

[8] Plummer M. A program for analysis of Bayesian graphical models using Gibbs sampling[D]. In Proceedings of the 3rd International Workshop on Distributed Statistical Computing， 2003，20-22.

[9] Von Bertalanffy. A quantitative theory of organic growth（inquiries on growth laws. II）[J]. Human Biology， 1938，10（2）：181-213.

[10] MacCall， Alec D. Quantitative fish dynamics[J]. Journal of the American Statistical Association， 2001，96（454）：781.

[11] Ntyniemi S H P， Whitlock R E， Per L TA， et al. General state-space population dynamics model for Bayesian stock assessment[J]. ICES Journal of Marine Science： Journal du Conseil， 2015，72（8）：2209-2222.

[12] Froese R， Demirel N， Coro G， et al. Estimating fisheries reference points from catch and resilience[J]. Fish and Fisheries， 2017，18（3）：506-526.

[13] Holt S. The evaluation of fisheries resources by the dynamic analysis stocks， and notes on the time factors involved[J]. International Commission for the Northwest Atlantic Fisheries， Special Publication，1958，1： 77–95.

[14] Beverton R J H， Holt S J. Manual of methods for fish stock assessment： Part 2 - tables of yield functions[J]. FAO Fisheries Technical Paper，1966，38：10.

[15] Gentleman R， Ihaka R. A language and environment for statistical computing[J]. Computing， 2011，1：12-21.

Abstract：In order to evaluate the utilization state of Sebastes schlegelii resources in the artificial reef area of Dachangshan and provide scientific basis for the reasonable development and management of Sebastes schlegelii resources，this paper adopts the bayesian evaluation method （LBB） based on body length to evaluate Sebastes schlegelii resources in the artificial reef area of Dachangshan from 2011 to 2015.Research results show that the Dachangshan artificial reef area for the Sebastes schlegelii optimal length Lopt to 31 cm， the best length Lc-opt for 28 cm，an average of M/K=1.5，F(xiàn)/K=4.82，Z/K=6.34，F(xiàn)/M=3.24，the Sebastes schlegelii optimal existing resources and the ratio of the amount of original B/B0（F=M，Lc=Lc-opt）=0.366，the actual average annual B/B0= 0.038 6，the optimal unit for the Sebastes schlegelii transshipment catches Y/R'（F=M，Lc=Lc-opt） = 0.045 9 kg. Actual annual average Y/R'=0.003 44 kg.

Key words：Dachangshan; Sebastes schlegelii;body length distribution;fish stock assessment;LBB

（收稿日期：2019-11-23）