基于體長標準的大長山人工魚礁區許氏平鮋資源評估

薛國進 尹增強 程前 魏藝鳴 李九奇

摘 要：為了評價大長山人工魚礁區許氏平鮋（Sebastes schlegelii）資源利用狀態，為當地許氏平鮋資源的合理開發與管理提供科學依據，采用基于體長的貝葉斯評估方法（LBB）對2011-2015年大長山人工魚礁區許氏平鮋資源進行了評估。研究結果表明大長山人工魚礁區許氏平鮋最優體長Lopt為31 cm，最佳開捕體長Lc-opt為28 cm，平均M/K=1.5、F/K=4.82、Z/K=6.34、F/M=3.24，許氏平鮋最優現存資源量與原始資源量之比B/B0（F=M，Lc=Lc-opt）=0.366、實際年平均B/B0= 0.038 6，許氏平鮋最優的單位補充量漁獲量Y/R'（F=M，Lc=Lc-opt）=0.045 9 kg，實際年平均Y/R'=0.003 44 kg。

關鍵詞：大長山;許氏平鮋（Sebastes schlegelii）;體長分布;資源評估;LBB

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

1 材料與方法

1.1 調查時間及站位設置

國家海洋科技支撐計劃項目組于2011-2015年連續五年對大長山人工魚礁區進行許氏平鮋資源調查，調查站位如圖1，調查方法參照《海洋漁業資源調查規范》[7]，并對許氏平鮋進行體長測定，統計結果如表1所示。

2 研究結果

2.1 2011-2015年大長山許氏平鮋各重要模型參數估計值

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

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

2.2 大長山許氏平鮋各重要模型參數的時間序列分布及資源利用評價圖

歷年漁獲物體長分布與最優體長比較與模型計算得到的各重要參數時間序列分布情況與繪圖表示，并分別列于圖2和圖3。

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

3 討論與分析

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

總之，無論從以上哪個指標來評價，大連大長山許氏平鮋資源利用均處于嚴重過度捕撈狀態，特別是從資源消耗程度指標B/B0來看，其資源利用情況已處于警戒線之下，資源處于極度危險狀態，需及時進行人為干預，建議專門設立幼魚保護區和禁漁期，開捕體長嚴格控制在28 cm以上，捕撈死亡系數嚴格控制在0.2以內。

參考文獻：

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

[2] 李蘭敏， 溫海深， 張亞晨. 卵胎生許氏平鮋個體生殖力與生物學指標的相關分析[J]. 現代農業科技， 2015（22）：267-269，279.

[3] 尹增強，許傳才，陳勇，等. 獐子島海域許氏平鮋生長特征與資源評估[J].農業經濟與科技，2015，26（10）：71-72

[4] 耿喆，朱江峰，夏萌，等. 數據缺乏條件下的漁業資源評估方法研究進展[J]. 海洋湖沼通報，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] 金顯仕，李顯森，趙憲勇，等.海洋漁業資源調查規范：SC/T 9403-2012[S].北京：中國農業出版社，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/K=4.82，Z/K=6.34，F/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）