东海北部三疣梭子蟹资源评估与管理决策分析.pdf

第5 3卷 第1 1期 2 0 2 3年1 1月中 国 海 洋 大 学 学 报P E R I O D I C A L O F O C E A N U N I V E R S I T Y O F C H I N A5 3(1 1):0 5 5 0 6 4N o v.,2 0 2 3东海北部三疣梭子蟹资源评估与管理决策分析刘 雅,王 晶,栗小东,王迎宾(浙江海洋大学水产学院,浙江 舟山 3 1 6 0 2 2)摘 要:三疣梭子蟹(P o r t u n u s t r i t u b e r c u l a t u s)是中国东海北部海域重要的经济蟹类,对三疣梭子蟹资源开展评估是其可持续利用的基础。本文运用基于贝叶斯方法的S c h a e f e r剩余产量模型,对东海北部海域的三疣梭子蟹资源进行评估,确定了当前三疣梭子蟹资源的开发利用状态,估算了在不同收获率水平下未来1 0年三疣梭子蟹的生物量和年总可捕捞量,分析了管理策略实施后三疣梭子蟹资源崩溃的风险。研究显示,在基准方案和敏感性分析方案下模型参数预测值以及生物学参考点估计值比较相近。使用两种方案估算得到的最大可持续产量(M S Y)均约为2 41 04 t,M S Y对应的生物量BM S Y均约为2 5 2 1 04 t。在基准方案下,M S Y对应的捕捞死亡系数FM S Y为0.0 9 6;在敏感性分析方案下,M S Y对应的捕捞死系数FM S Y为0.0 9 7。在2 0 0 12 0 2 0年间,捕捞死亡系数均低于FM S Y,且生物量基本在BM S Y水平之上,表明近年来三疣梭子蟹未出现过度捕捞的情况。决策分析和风险分析表明,为使三疣梭子蟹资源可持续利用,将收获率设定在0.0 4左右是较合理的捕捞方案。关键词:三疣梭子蟹;贝叶斯方法;剩余产量模型;管理策略;东海北部;渔业资源中图法分类号:S 9 3 2.5 文献标志码:A 文章编号:1 6 7 2-5 1 7 4(2 0 2 3)1 1-0 5 5-1 0D O I:1 0.1 6 4 4 1/j.c n k i.h d x b.2 0 2 2 0 1 2 6引用格式:刘雅,王晶,栗小东,等.东海北部三疣梭子蟹资源评估与管理决策分析J.中国海洋大学学报(自然科学版),2 0 2 3,5 3(1 1):5 5-6 4.L i u Y a,W a n g J i n g,L i X i a o d o n g,e t a l.S t o c k a s s e s s m e n t a n d m a n a g e m e n t d e c i s i o n a n a l y s i s o f P o r t u n u s t r i t u b e r c u l a t u s i n-h a b i t i n g n o r t h e r n E a s t C h i n a S e aJ.P e r i o d i c a l o f O c e a n U n i v e r s i t y o f C h i n a,2 0 2 3,5 3(1 1):5 5-6 4.基金项目:国家重点研究发展计划项目(2 0 1 9 Y F D 0 9 0 1 3 0 4);浙江省基础公益计划项目(L G N 2 1 C 1 9 0 0 0 9);舟山市科技计划项目(2 0 2 2 C 4 1 0 0 3)资助S u p p o r t e d b y t h e N a t i o n a l K e y R e s e a r c h a n d D e v e l o p m e n t P r o g r a m o f C h i n a(2 0 1 9 Y F D 0 9 0 1 3 0 4);t h e B a s i c P u b l i c W e l f a r e P l a n o f Z h e j i a n g(L G N 2 1 C 1 9 0 0 0 9);t h e S c i e n c e a n d T e c h n o l o g y P l a n P r o j e c t o f Z h o u s h a n(2 0 2 2 C 4 1 0 0 3)收稿日期:2 0 2 2-0 3-0 2;修订日期:2 0 2 2-0 5-1 9作者简介:刘 雅(1 9 9 5),女,硕士生。E-m a i l:3 3 1 7 7 4 5 4 3 q q.c o m 通信作者:E-m a i l:y i n g b i n w a n g 1 2 6.c o m 三疣梭子蟹(P o r t u n u s t r i t u b e r c u l a t u s)隶属甲壳纲十足目梭子蟹科梭子蟹属,为广温广盐性种类,是中国最重要的海洋经济蟹类之一1。主要分布于日本、朝鲜、菲律宾、马来群岛、红海及中国的渤海、黄海、东海和南海。在东海三疣梭子蟹主要分布在大沙渔场、长江口渔场和舟山渔场2,其渔获量约占全国产量的5 0%3。东海北部海域是东海最大的三疣梭子蟹产区,其渔获量占东海总渔获量的4 0%以上4。国外学者对三疣梭子蟹的相关研究较少,A r i y a m a H5基于1 9 8 42 0 0 8年日本大阪湾区域三疣梭子蟹的渔获量数据,探究了过于频繁的渔业开发和偶发的人类活动,以及气候事件对三疣梭子蟹补充量造成的影响。我国有关三疣梭子蟹的研究主要集中于养殖领域,包括三疣梭子蟹的人工育苗6、池塘养殖7、常见病害防控8、配合饲料开发9-1 0以及优良品种选育1 1等。对野生三疣梭子蟹资源的相关研究相对有限,主要集中在资源分布、增殖容量和资源评估方面。如徐炳庆等1 2通过轮廓系数法对莱州湾三疣梭子蟹资源的时空分布特征进行了研究;栗小东等1 3基于机器学习方法,对东海北部海域三疣梭子蟹的资源分布进行了预测;林群等1 4构建了黄河口邻近海域的E c o p a t h模型,初步评估了三疣梭子蟹的增殖容量;俞存根4对东海中北部重要的经济蟹类及其现存资源量和最大可持续产量(M a x i-m u m s u s t a i n a b l e y i e l d,M S Y)进 行 了 分 析;W a n g等1 5基于传统S c h a e f e r模型对三疣梭子蟹的最大可持续产量及资源开发情况开展了评估。以上在对三疣梭子蟹进行资源评估的研究中并未考虑模型中相关参数的不确定性。在资源评估过程中,如果不考虑不确定性,将对渔业资源评估结论、渔业生产和生产管理造成负面影响。因此,实现渔业科学管理和合理开发的前提是对资源量的正确评估1 6。本研究依据2 0 0 12 0 2 0年东海北部海域三疣梭中 国 海 洋 大 学 学 报2 0 2 3 年子蟹的渔获量数据和单位捕捞努力量渔获量(C a t c h p e r u n i t e f f o r t,C P U E)数据,运用S c h a e f e r剩余产量模型,将参数基准方案(参数的先验分布均为均匀分布)与敏感性分析方案(参数与基准方案相同,但采用不同先验分布)进行对比,对东海北部海域三疣梭子蟹的资源进行评估,并设定收获率作为三疣梭子蟹资源管理策略且分析管理策略可能存在的风险。研究结果有助于掌握当前三疣梭子蟹的渔业资源变动规律,并基于目前三疣梭子蟹资源状况对未来资源的动态变化趋势进行预估。1 材料和方法1.1 数据来源东海北部海域三疣梭子蟹渔获量数据来源于 中国渔业统计年鉴1 7,时间段为2 0 0 12 0 2 0年(见图1)。除了资源量以外,其他一些因素也会导致C P U E的变动,因此需要对C P U E进行标准化,以剔除这些因素的影响。本研究使用了近2 0年的渔业生产数据,其中包含了3种三疣梭子蟹生产船型,同时考虑主要环境因子变化的影响,因此在使用广义线性模型(G e n e r a l i z e d l i n e a r m o d e l,G L M)对C P U E数据进行标准化时1 8,考虑了年份、渔船类型和平均水温3个因素的影响。具体模型如下:l m(l n C P U E f a c t o r(Y e a r)+f a c t o r(E f f o r t)+f a c t o r(V e s s e l)+f a c t o r(T e m p),d a t a=R D a t a)。(1)式中:l m是R语言中广义线性模型的函数命令;l n C P U E代表对C P U E取对数;Y e a r代表年份;E f f o r t代表捕捞努力量;V e s s e l代表3种不同的船型,分别是刺网、蟹笼、单拖网;T e m p代表底层温度。图1 2 0 0 12 0 2 0年东海北部海域三疣梭子蟹年渔获量和标准化C P U E数据F i g.1 S t a t i s t i c s o f c a t c h a n d s t a n d a r d i z e d C P U E o f P o r t u n u s t r i t u b e r c u l a t u s t h e N o r t h e r n E a s t C h i n a S e a f r o m 2 0 0 1 t o 2 0 2 01.2 剩余产量模型S c h a e f e r非平衡剩余产量模型的离散时间形式如下:Bt=Bt-1+r Bt-1(1-Bt-1K)-Ct-1;(2)It=q Btet-22,tN(0,2)。(3)式中:Bt为t年的生物量;r为内禀自然增长率;K为环境容纳量;Ct-1为t-1年的渔获量;q为可捕系数。其中r、K、q为模型的参数。It为三疣梭子蟹资源丰度的相对指数,t为误差项。在进行资源评估的过程中,一般会假设剩余产量模型中初始年的生物量B0为K,通过减少模型估算的参数数量来达到降低模型评估难度的目的1 9。但这种假设成立的前提条件是能够确定渔业开始时前一年的资源量2 0。由于中国对三疣梭子蟹的捕捞活动开始较早,因此将初始一年的资源量视为K并不合理。本研究的B0值参考W a n g等1 5的研究结果。1.3 似然函数假设C P U E正比于生物量,且观测误差服从对数正态分布,则似然函数的表达式为:L(I|)=2 0 2 0t=2 0 0 11It2 e x p(l n(It)-l n(q Bt)222 。(4)式中:L和l n均表示自然对数;I表示单位捕捞努力量渔获量;表示随机变量;表示连乘符号;q表示可捕系数;Bt表示t年的生物量;表示标准差。651 1期刘 雅,等:东海北部三疣梭子蟹资源评估与管理决策分析由于本研究采用的渔获量和C P U E数据的时间序列为2 0年,这对标准差()的估算来说存在一定的不稳定性,因此在本研究将其视为恒定值(=0.2)1 9。1.4 模型参数先验分布的设定在贝叶斯理论中,参数的后验分布会受到先验分布设定的影响2 1。因此,要对先验分布进行敏感性检验来验证样本数据提供的信息是否充足。具体方法是比较各参数在不同先验分布下得到的后验分布的差异大小,来判断样本数据提供的信息是否充足。若结果差异较大,则表明数据没有提供足够的信息1 8。因为较难搜集到三疣梭子蟹相关参数的先验信息,为了结果的合理性,故采用均匀分布作为模型参数先验概率分布的基准方案。此方案K的基准先验分布假设服从均匀分布(5 01 04 t,5 5 01 04 t),下边界设置为5 0 1 04 t,此值大于2 0 0 12 0 2 0年数据中的最大渔获量。为了尽量减少上边界值对后验概率的影响,上边界值被设定为5 5 0 1 04 t2 2。对于可捕系数q,获取其有信息的先验分布具有较大难度,在研究中一般将这种无信息的先验分布设置为均匀分布或对数均匀分布。为足以覆盖内禀增长率r的合理范围,r的上、下边界区间一般包括与评估种类相似的物种r的估计值2 3。本研究为了验证样本数据是否提供了足够的信息,用敏感性分析方案来替代基准方案中的参数先验分布,比较两种情况下得到的参数后验分布差异大小。根据三疣梭子蟹的相关研究1 5,参数r、K、q基准方案和敏感性分析方案的先验分布设定见表1。表1 剩余产量模型参数r、K、q的先验概率分布T a b l e 1 T h e p r i o r p r o b a b i l i t y d i s t r i b u t i o n o f t h e p a r a m e t e r s r,K,a n d q i n t h e s u r p l u s p r o d u c t i o n m o d e l方案S c e n a r i o内禀增长率rI n t r i n s i c r a t e o f n a t u r a l i n c r e a s e r承载力KC a r r y i n g c a p a c i t y K/1 04 t可捕系数qF i s h i n g c o e f f i c i e n t q基准方案B a s e s c e n a r i oU(0.0 1,2)U(5 0,5 5 0)U(0.0 0 0 1,0.0 5)敏感性分析方案S e n s i t i v i t y a n a l y s i s s c e n a r i oN(1,0.52)U(5 0,5 5 0)l n(q)U(-2 0,0)1.5 模型后验分布的计算贝叶斯方法不仅考虑了评估模型的不确定性,而且考虑了参数的不确定性,通过将经验信息与似然函数结合在一起,得到模型参数的后验分布,所以可以得到比传统分析方法更可靠的结果。因此,贝叶斯方法被认为是目前国际渔业资源研究领域较有代表性的方法,且适用于缺乏某些特定数据的情况下研究渔业资源动态2 4。本研究中模型参数r、K、q的后验概率分布是通过MC MC(M a r k o v c h a i n m o n t e c a r l o a l g o r i t h m)算法测算的。此算法共迭代5 0 0 0 0次,将前5 0 0 0次的结果舍弃,后每隔4 0次储存剩下的4 5 0 0 0次的运算结果。采用此算法计算时,各参数的初始值设定分别为:r=0.1 5,K=3 7 7 0 6.4,q=0.0 0 51 5。1.6 生物学参考点估算 生物学参考点是从生物学角度设定的指标,包括目标参考点(T a r g e t r e f e r e n c e p o i n t s,T R P)和限制参考点(L i m i t r e f e r e n c e p o i n t,L R P)。目标参考点是渔业管理为了持续获得某一渔获量所设定的指标,限制参考点指的是在渔业管理中为了避免出现种群的生物量过低导致种群不能被可持续利用的情况而设置的指标2 5。本研究涉及到的生物参考点有FM S Y、BM S Y、F0.1和M S Y。FM S Y和BM S Y分别指渔业达到M S Y水平时对应的捕捞死亡系数和生物量,F0.1表示平衡渔获量和捕捞死亡系数关系曲线最大斜率的1 0%对应的捕捞死亡系数。通过计算可得:FM S Y=r2,(5)F0.1=0.4 5r,(6)M S Y=r K4,(7)BM S Y=K2。(8)在本研究中,捕捞死亡系数的目标参考点Ft a r设定为F0.1,限制参考点Fl i m设定为FM S Y;BM S Y设定为资源量的目标参考点Bt a r,BM S Y/4设定为限制参考点Bl i m。可以通过以上参考点来判断当前的渔业资源状况。当捕捞死亡系数大于Fl i m时,表明该资源正在遭受过度捕捞(O v e r f i s h i n g),相反则表明该资源未遭受过度捕捞;当资源量小于Bl i m时,表明该资源已经遭到破坏,处于过度捕捞状态(O v e r f i s h e d),相反则表明该资源未处于过度捕捞状态。1.7 确定备选管理策略本文采用不同的收获率(H a r v e s t r a t e,即捕捞死亡率,规定了渔业年捕捞渔获量在资源量中的固定占比)来作为三疣梭子蟹资源管理策略。收获率分别设定为0.0 1、0.0 2、0.0 3、0.0 4、0.0 5、0.0 6、0.0 7和0.0 8。因75中 国 海 洋 大 学 学 报2 0 2 3 年此,未来t年的渔获量的计算公式如下:Ct=hiBte,N(0,12)。(9)式中:Ct为t年的渔获量;hi为收获率;Bt为t年的生物量;为误差项。1.8 管理策略实施效果模拟资源评估过程的复杂性会导致评估结果存在较多不确定性,因此通过建立指标来对不同的管理策略进行风险评估是十分必要的。本研究假设实施管理策略的时间段为2 0 2 12 0 3 0年。建立的指标有:(1)管理结束时的生物量B2 0 3 0,即2 0 3 0年资源量的期望值。(2)管理策略实施后2 0 3 0年渔获量的期望值C2 0 3 0。(3)2 0 3 0年三疣梭子蟹生物量与BM S Y之比的期望值(B2 0 3 0/BM S Y)。(4)2 0 3 0年三疣梭子蟹生物量的衰减率(D e p l e-t i o n),即实施某一管理策略后,2 0 3 0年三疣梭子蟹生物量与K的比例(B2 0 3 0/K)。(5)2 0 3 0年三疣梭子蟹生物量大于2 0 2 0年三疣梭子蟹生物量的概率P(B2 0 3 0B2 0 2 0),表示管理措施实施后资源较管理措施实施前资源恢复的概率。(6)2 0 3 0年三疣梭子蟹生物量大于BM S Y的概率P(B2 0 3 0BM S Y),表示资源恢复到健康水平的概率。(7)2 0 3 0年三疣梭子蟹生物量小于BM S Y/4的概率P(B2 0 3 0BM S Y,基准方案和敏感性分析方案的概率期望值分别为0.5 3和0.5 6。对于B2 0 3 0B2 0 2 0的概率期望值均大于0.6,B2 0 3 0BM S Y的概率期望值均大于0.8,且B2 0 3 0B2 0 2 0)P(B2 0 3 0BM S Y)P(B2 0 3 0B2 0 2 0):2 0 3 0年生物量的期望值大于2 0 2 0年生物量的期望值的概率 T h e p r o b a b i l i t y t h a t t h e e x p e c t e d v a l u e o f b i o m a s s i n 2 0 3 0 i s g r e a t e r t h a n t h e e x p e c t e d v a l u e o f b i o m a s s i n 2 0 2 0;P(B2 0 3 0BM S Y):2 0 3 0年生物量的期望值大于渔业达到M S Y水平时生物量的期望值的概率 T h e p r o b a-b i l i t y t h a t t h e e x p e c t e d v a l u e o f b i o m a s s i n 2 0 3 0 i s g r e a t e r t h a n t h e e x p e c t e d v a l u e o f b i o m a s s w h e n t h e f i s h e r y r e a c h e s t h e m a x i m u m s u s t a i n a b l e y i e l d l e v e l;P(B2 0 3 0BM S Y/4):2 0 3 0年生物量的期望值小于渔业达到M S Y/4水平时生物量的期望值的概率 T h e p r o b a b i l i t y t h a t t h e e x p e c t e d v a l u e o f b i o-m a s s i n 2 0 3 0 i s l e s s t h a n t h e e x p e c t e d v a l u e o f b i o m a s s w h e n t h e f i s h e r y r e a c h e s t h e 1/4 l e v e l o f t h e m a x i m u m s u s t a i n a b l e y i e l d。16中 国 海 洋 大 学 学 报2 0 2 3 年3.2 当前资源状况评价及未来管理策略分析 当前东海北部海域三疣梭子蟹资源状况相对乐观,在近年来并未遭受到过度捕捞且三疣梭子蟹的资源状况良好。本研究得出的捕捞死亡系数F较小,可能是由于三疣梭子蟹适应力强、存活率较高且增殖放流数量大,使得三疣梭子蟹资源量维持在较高水平,因此估算出的捕捞死亡系数F较小。本研究得出的三疣梭子蟹的M S Y约为2 4 1 04 t,这一结果不同于W a n g等1 5研究结果。W a n g等1 5研究显示,假设三疣梭子蟹的捕捞强度和放流数量保持在2 0 1 5年的水平,三疣梭子蟹的M S Y约为1 4.2 1 04 t。造成结果不同的原因可能是采用了不同的评估方法或所用数据的时间跨度不同。W a n g等1 5使用的是传统的非平衡剩余产量模型,评估所使用的数据是2 0 0 12 0 1 5年间的渔业统计数据。当采用相同时间跨度的数据(2 0 0 12 0 1 5年)时,使用贝叶斯方法得到的M S Y为1 91 04 t,与W a n g等1 51 4 1 04 t的研究结果差异明显缩小。具体分析可得,2 0 0 12 0 1 5年的年产量均值约为1 0.41 04 t,而采用2 0 0 12 0 2 0年的渔业统计数据算出的年产量均值约为1 2.11 04 t。2 0 1 5年的年产量约为1 8.5 1 04 t,与本文评估得到的M S Y值更接近。对于W a n g等1 5未使用的时间段即2 0 1 62 0 2 0年,东海北部海域三疣梭子蟹的年产量均值约为1 7.01 04 t,明显大于2 0 0 12 0 1 5年的均值,且2 0 0 12 0 1 5年C P U E均值也偏小。对数据趋势分析可得,2 0 1 32 0 2 0年的C P U E均值明显大于2 0 0 12 0 1 2年,且2 0 1 3年之后三疣梭子蟹年产量相较于之前明显增加,2 0 1 5年之后的产量虽略有下降,但总体仍保持较高水平。因为本文采用的数据时间跨度为2 0年,故后期研究海域三疣梭子蟹统计数据的变化可能会对模型产生较大影响,导致评估结果(如M S Y)数据变大。由此可见,采用不同时间段的数据做出的评估结果可能会产生差异,影响主要来自渔获量和捕捞努力量的年际变化趋势。因此,全面且准确的长时间序列数据对于评估M S Y及相关的生物学参考点至关重要。渔业资源管理的主要目标是通过合理的调控与科学的管理,从渔业资源中获得最佳的可持续利益。这就要求在对渔业管理策略进行选择时,应从多角度考虑,要在保证渔业资源可持续开发利用状态与获得高渔获量之间进行全面权衡。综合两种方案下的管理策略和风险指标分析结果可知,若要获取最大产量,应将收获率控制在0.0 6。但风险分析表明,此时的资源状况并不乐观。若将收获率设定为0.0 6开展捕捞活动,三疣梭子蟹资源会出现过度捕捞现象,甚至有可能会出现资源崩溃的情况。所以,将收获率设置为0.0 6对资源的可持续发展来说存在很大风险。从资源保护的角度来看,将三疣梭子蟹资源的收获率控制在0.0 3时,资源量会维持在较高的水平且不会出现资源崩溃的状况。从经济收益的角度来看,将收获率控制在0.0 5时,基准方案和敏感性分析方案中2 0 3 0年最大渔获量期望值与结果中最大值非常接近,且出现资源崩溃的可能性较低。但在大多数渔业已经出现衰退趋势的现状下,我们应采取使资源状况更安全的策略作为最佳管理策略。对利益及其存在的风险进行权衡后得出,0.0 4的收获率在未来1 0年的管理期内是较合理且保守的管理策略。3.3 研究的不足和展望本研究分基准和敏感性分析两种方案,使用了与贝叶斯方法相结合的S c h a e f e r剩余产量模型,对东海北部海域三疣梭子蟹的资源状况进行了评估,并估算了三疣梭子蟹资源在2 0 3 0年恢复到健康水平的概率和崩溃的概率。研究结果可为三疣梭子蟹资源养护和可持续利用提供有利支持,为制定科学合理的渔业管理策略提供参考。但是,本研究仍然存在一些不足和需改善之处。第一,模型参数的不确定性,本研究假设了C P U E数据的标准差为0.2,该值的假设对评估结果造成的影响需要进一步研究与分析。第二,本研究采用的生产数据时间段为2 0 0 12 0 2 0年,因此模型参数的结果只能较客观地反映这一时期的渔业状态。由于三疣梭子蟹属于易受环境影响的短生命周期物种,且渔业统计年鉴中对该种的统计也存在误差,因此,基于2 0 0 12 0 2 0年间估算出的模型参数对其未来的渔业资源变化进行的预测会存在一定的不确定性。第三,三疣梭子蟹是东海北部海域重要的增殖放流种类之一,每年幼蟹放流数量达到几千万尾,而本研究并未将该影响纳入模型评估当中,这无疑会对资源评估结果产生一定的影响。今后的研究中,我们将继续积累更多的统计和调查数据,进一步完善评估模型,在评估过程中充分考虑环境因子和增殖放流等因素的影响,并尝试使用生态系统模型对三疣梭子蟹的M S Y开展评估,为合理开发、有效保护三疣梭子蟹资源提供科学参考。参考文献:1 H e J,G a o Y,W a n g W,e t a l.L i m b a u t o t o m y p a t t e r n s i n t h e j u-v e n i l e s w i mm i n g c r a b(P o r t u n u s t r i t u b e r c u l a t u s)i n e a r t h p o n d sJ.A q u a c u l t u r e,2 0 1 6,4 6 3:1 8 9-1 9 2.2 宋海棠,俞存根,薛利建,等.东海经济虾蟹类M.北京:海洋出版社,2 0 0 6:8 3-8 5.S o n g H T,Y u C G,X u e L J,e t a l.T h e E c o n o m i c a l C r a b s a n d S h r i m p s i n E a s t C h i n a S e aM.B e i j i n g:C h i n a O c e a n P r e s s,2 0 0 6:8 3-8 5.3 L i u S,S u n J,H u r t a d o L A.G e n e t i c d i f f e r e n t i a t i o n o f P o r t u n u s t r i t u b e r c u l a t u s,t h e w o r l d s l a r g e s t c r a b f i s h e r y,a m o n g i t s t h r e e 261 1期刘 雅,等:东海北部三疣梭子蟹资源评估与管理决策分析m a i n f i s h i n g a r e a sJ.F i s h e r i e s R e s e a r c h,2 0 1 3,1 4 8:3 8-4 6.4 俞存根,宋海棠,姚光展.东海大陆架海域蟹类资源量的评估J.水产学报,2 0 0 4,2 8(1):4 1-4 6.Y u C G,S o n g H T,Y a o G Z.A s s e s s m e n t o f t h e c r a b s t o c k b i o-m a s s i n t h e c o n t i n e n t a l s h e l f w a t e r s o f t h e E a s t C h i n a S e aJ.J o u r n a l o f F i s h e r i e s o f C h i n a,2 0 0 4,2 8(1):4 1-4 6.5 A r i y a m a H,S e c o r D H.E f f e c t o f e n v i r o n m e n t a l f a c t o r s,e s p e c i a l l y h y p o x i a a n d t y p h o o n s,o n r e c r u i t m e n t o f t h e g a z a m i c r a b P o r t u n u s t r i t u b e r c u l a t u s i n O s a k a B a y,J a p a nJ.F i s h e r i e s S c i e n c e,2 0 1 0,7 6(2):3 1 5-3 2 4.6 廖永岩,肖展鹏,袁耀阳.三疣梭子蟹幼体和幼蟹的温度适应性J.水生生物学报,2 0 0 8,3 2(4):5 3 4-5 4 3.L i a o Y Y,X i a o Z P,Y u a n Y Y.T e m p e r a t u r e t o l e r a n c e o f l a r v a a n d j u v e n i l e o f P o r t u n u s t r i t u b e r c u l a t u sJ.A c t a H y d r o b i o l o g i c a S i n i c a,2 0 0 8,3 2(4):5 3 4-5 4 3.7 钱丽君,张德民,徐小红.应用D G G E分析三疣梭子蟹养殖塘底泥细菌的多样性J.水产学报,2 0 0 7,3 1(2):2 0 4-2 1 0.Q i a n L J,Z h a n g D M,X u X H,B a c t e r i a l d i v e r s i t y o f P o r t u n u s t r i t u b e r c u l a t u s p o n d s e d i m e n t e s t i m a t e d b y d e n a t u r i n g g r a d i e n t g e l e l e c t r o p h o r e s i sJ.J o u r n a l o f F i s h e r i e s o f C h i n a,2 0 0 7,3 1(2):2 0 4-2 1 0.8 许文军,施慧,徐汉祥,等.养殖梭子蟹血卵涡鞭虫感染的初步研究J.水生生物学报,2 0 0 7,3 1(5):6 3 7-6 4 0.X u W J,S h i H,X u H X,e t a l.P r e l i m i n a r y s t u d y o n t h e H e m a t o-d i n i u m i n f e c t i o n i n c u l t u r e d P o r t u n u s t r i t u b e r c u l a t u sJ.A c t a H y d r o b i o l o g i c a S i n i c a,2 0 0 7,3 1(5):6 3 7-6 4 0.9 张稳,谢奉军,金敏,等.饲料中n-3高不饱和脂肪酸含量对三疣梭子蟹幼蟹生长性能及脂肪酸组成的影响J.动物营养学报,2 0 1 4,2 6(5):1 2 5 4-1 2 6 4.Z h a n g W,X i e F J,J i n M,e t a l.E f f e c t s o f d i e t a r y n-3 h i g h l y u n-s a t u r a t e d f a t t y a c i d c o n t e n t o n g r o w t h p e r f o r m a n c e a n d f a t t y a c i d c o m p o s i t i o n o f j u v e n i l e s w i mm i n g c r a b(P o r t u n u s t r i t u b e r c u l a t u s)J.C h i n e s e J o u r n a l o f A n i m a l N u t r i t i o n,2 0 1 4,2 6(5):1 2 5 4-1 2 6 4.1 0 J i n M,Z h o u Q C,W a n g M Q,e t a l.D i e t a r y a r g i n i n e r e q u i r e-m e n t o f j u v e n i l e s w i mm i n g c r a b,P o r t u n u s t r i t u b e r c u l a t u sJ.A q u a c u l t u r e N u t r i t i o n,2 0 1 6,2 2(6):1 1 7 4-1 1 8 4.1 1 G a o B Q,L i u P,L i J,e t a l.E f f e c t o f i n b r e e d i n g o n g r o w t h a n d g e n e t i c d i v e r s i t y o f P o r t u n u s t r i t u b e r c u l a t u s b a s e d o n t h e f u l l-s i b-l i n g i n b r e e d i n g f a m i l i e sJ.A q u a c u l t u r e I n t e r n a t i o n a l,2 0 1 5,2 3:1 4 0 1-1 4 1 0.1 2 徐炳庆,陈玮,王田田,等.莱州湾“伏休”结束前三疣梭子蟹的资源状况及其分布特征J.水产学报,2 0 2 1,4 5(4):5 4 3-5 5 1.X u B Q,C h e n W,W a n g T T,e t a l.R e s o u r c e s t a t u s a n d d i s t r i-b u t i o n c h a r a c t e r i s t i c s o f P o r t u n u s t r i t u b e r c u l a t u s b e f o r e t h e e n d o f s u mm e r f i s h i n g m o r a t o r i u m i n L a i z h o u B a yJ.J o u r n a l o f F i s h e r i e s o f C h i n a,2 0 2 1,4 5(4):5 4 3-5 5 1.1 3 栗小东,王晶,杨春蕙,等.基于两种机器学习方法分析东海北部海域三疣梭子蟹(P o r t u n u s t r i t u b e r c u l a t u s)时空分布J.海洋与湖沼,2 0 2 1,5 2(5):1 2 8 4-1 2 9 2.L i X D,W a n g J,Y a n g C H,e t a l.S p a t i o t e m p o r a l d i s t r i b u t i o n o f P o r t u n u s t r i t u b e r c u l a t u s i n t h e N o r t h e r n E a s t C h i n a S e a b a s e d o n t w o m a c h i n e l e a r n i n g m e t h o d sJ.O c e a n o l