Page 78 - 《水产学报》2026年第3期

3 期 水 产 学 报 50 卷

系统产生影响 。 14 14
[6]
渔获量 catch
贝叶斯状态空间剩余产量模型 (JABBA) 所需 12 CPUE 12
要的仅仅为渔获量和单位捕捞努力量 (CPUE) 数 10 10
据，该模型还可以拟合多种 CPUE 的时间序列、
估 计 模 型 的 过 程 方 差 、 估 计 单 个 或 者 成 组 的 年鱼获量/(×10 4 t) catch 8 8 CPUE/[t/(船·a)]

CPUE 的时间序列的观测方差，还可以通过设置 6 6
[7]
B MS Y /K及其形状参数 m 确定鱼类生产函数 ，因 4 4
[8]
此得到了广泛的应用。符武月等 利用 JABBA 模 2 2
型 对 南 极 犬 牙 鱼 中 的 小 鳞 犬 牙 鱼 (Dissostichus
0 0
eleginoides) 和莫氏犬牙鱼 (D. mawsoni) 两个经济 1996 2000 2005 2010 2015 2020
种的资源动态进行了评估，该方法利用历年的犬 年份
year
牙鱼渔获量和标准化后的捕捞努力量数据，通过
模型对其种群状态参数进行了评估。田志盼等 [9] 图 1 1995—2020 年西北太平洋柔鱼年渔获量和
利用此方法对日本鲭 (Scomber japonicus) 资源进行 CPUE 的年间变化
了评估。赵蓬蓬等 [10] 对印度洋大眼金枪鱼 (Thun- Fig. 1 Annual catch and CPUE for O. bartramii in
nus obesus) 进行了评估。 the Northwest Pacific from 1995 to 2020
[ ( ) m−1 ]
许多学者利用不同模型对西北太平洋柔鱼进 SP = r B 1− B (1)
行了评估，Wang 等 [11] 利用环境依赖型剩余产量 m−1 K
模型 (EDSP) 对西北太平洋柔鱼进行了评估，结果 式中，SP 为剩余产量，r 为种群内禀增长率，m
显示，柔鱼渔业的生物量和发展处于良好状态。 为形状参数，B 为生物的资源量，K 为环境容纳
隋芯等 [12] 通过建立 4 种剩余产量模型，探索厄尔 量。当 m=2 时，函数为 Schaefer 形式，当 m 趋近
尼诺事件影响下西北太平洋柔鱼的资源状况变化 于 1 时为 Fox 形式，m 为其他值时为 Pella-Tom-
趋势，其研究表明传统剩余产量模型中柔鱼资源 linson 形 式 。 本 研 究 中 ， 采 用 “m=2”的 Schaefer
没有出现过度捕捞，资源状态良好。柔鱼生命周 模型对西北太平洋柔鱼进行评估。
期较短 ，一生只产卵 1 次，产卵后就会死去 ， 相关生物学参考点之间的关系：
[13]
[14]
所以其资源状态极易变动，特别近几年来，由于 rK
MSY = (2)
对柔鱼捕捞强度的增大以及海洋环境的变化，柔 4
鱼资源出现了明显的下降，因此亟需更新评估， K
B MSY = (3)
进而制定合理的管理措施。本研究以西北太平洋 2
柔鱼渔获量数据和 CPUE 数据为基础，采用 JABBA r
F MSY = (4)
模型对其资源状况进行评估，并且使用最新的西 2
北太平洋柔鱼相关数据，可以为其提供更准确的 JABBA 中，过程方程定义如下：
资源状况信息，为其管理提供科学依据。 
 φe ,y = 1
η y

 [ ]

 ( )
 r m−1 C y−1

1 材料与方法 P y =  P y−1 + P y−1 1− P y−1 −

 m−1 K


 ηy
 e ,y = 2,3,4,··· ,n 时
1.1 数据来源 (5)
西北太平洋柔鱼渔业数据来自上海海洋大学 式中， P y 为 y 年 B 和 K 的比值， φ e 为初始资源消
2
2
中国远洋渔业数据中心，时间为 1995—2020 年， 耗率， η y 为过程误差项，且 η y ~N (0， σ )， σ 为
η
η
空间范围为 38°~46°N，150°~170°E，渔业数据包 过程方差，在本研究中设其服从逆伽马分布 inver-
括年渔获量和 CPUE (图 1)。 Gamma(4, 0.01)， C y−1 为 y−1 年的渔获量。
1.2 剩余产量模型 观测方程为：
JABBA 中定义的剩余产量函数形式： I i,y = q i B i e ε y,i (6)
https://www.china-fishery.cn 中国水产学会主办 sponsored by China Society of Fisheries
2