Page 106 - 《水产学报》2026年第04期

4 期 冯 波，等：基于渔获量调查数据的南海北部金线鱼捕捞限额评估 50 卷
C t
CPUE t = (1) B MSY =0.5k (6)
MP t D t
模型的先验参数来自历史文献，详见表 2。
式中，C 表示该作业方式与功率段 t 年的渔获量；
t
时间序列结尾年份的种群状况判断按 Froese 等 [16]
MP 表示该作业方式与功率段 t 年的主机功率投
t
给出的定义 (表 3)。金线鱼的种群参数估计、渔获
入数；D 表示该作业方式与功率段在 t 年的作业
t
量拟合、CPUE 拟合的残差分析、生物学参考点
天数。
估计、种群状态分析和 Kobe 图概率计算用 Froese
不同作业方式与功率段的渔获量占比是指不
等 [17] 提供的 CMSY++16.R 代码在 R 4.4.3 环境下
同作业方式与功率段在 2008—2023 年的累计渔获
完成运算分析。
量除以所有作业方式与功率段在 2008—2023 年总

累计渔获量的比值。16 个不同作业方式与功率段
表 2 CMY++模型输入参数设定
的金线鱼渔获量相对于总体渔获量的占比见表 1。
Tab. 2 Input settings for CMSY++ model
其中 3 个渔获量占比较高的 C、H、N 的 CPUE 数
参数 意义 赋值 文献
据将分别代表拖网、刺网和钓具，投入 CMSY++ parameter meaning value reference
模型进行种群状态深化分析。 resilience 物种恢复力 中等 [7-8]
1.2 模型方法 r.low/r.hi 物种恢复力的范围 0.60/1.37 [7]
stb.low/stb.hi 起始年份的B/k范围 0.2/0.6 [9,15]
CMSY++模型方法是一种创新的贝叶斯状态
endb.low/endb.hi 结尾年份的B/k范围 0.2/0.6 [8]
空间方法 ，是对 CMSY 方法 [13] 的升级增强，而
[12]

CMSY 方 法 又 是 对 Catch-MSY 方 法 [14] 的 改 进 。

CMSY++模型方法由 2 个分析部分组成，这两个
表 3 种群开发利用状况评价
部 分 都 是 基 于 改 进 的 Schaefer 剩 余 产 量 模 型 ：
Tab. 3 Assessment of stock status
① CMSY 部分，处理只有渔获量的时间序列数据；
种群状况
② BSM 部分，需要额外的资源丰度指数。使用渔 B/B MSY F/F MSY stock status
获量时间序列和资源丰度指数数据来估算生产力。 ≥1 ＜1 健康
CMSY 部分随机从初始的 r 和 k 范围内选择一对 r- 0.5~1 ＜1 种群恢复过程中
k 组合，初始生物量也从第 1 年的 B/k 范围中选择， <0.5 ＜1 超出种群安全阈值
然后使用公式 (1) 计算后续几年的生物量，确保结 0.5~1 ＞1 过度捕捞
尾年份的 B/k 接近设定范围。使用马尔可夫链蒙 0.2~0.5 ＞1 严重过度捕捞
特卡洛 (MCMC) 方法反复迭代运算筛选出合适的 ＜0.2 ＞1 资源枯竭
r-k 组合。
( )
B t
B t+1 = B t +r 1− B t −C t (2) 2 结果
k
式中，B 表示 t 年时的生物量，r 表示内禀增长率；
t
k 表示环境承载量；C 表示 t 年时的捕捞产量。 2.1 捕捞生产统计
t
BSM 部分分析时，引入了 CPUE 指数，式 (1) 根据渔港抽样统计资料，南海北部的金线鱼
变： 被拖网、刺网、钓具、杂项渔具捕捞 (图 1)，拖网
渔获量占比为 3.32%~17.92%，平均为 9.89%；刺
( )
CPUE t
CPUE t+1 = CPUE t +r 1− CPUE t −qC t
qk 网渔获量占比为 51.85%~86.45%，平均为 71.87%；
(3) 钓具渔获量占比为 7.18%~28.64%，平均为 17.51%；
式中，q 为可捕系数。 杂 项 渔 具 渔 获 量 占 比 为 0.01%~3.25%， 平 均 为
CMSY++模型输出重要的种群参数计算式： 0.73%。2008—2023 年南海北部的金线鱼渔获量
总体呈一涨一落的波动特征，最低渔获量 6.16 万 t
MSY=rk/4 (4)
出 现 在 2020 年 ， 最 高 渔 获 量 27.09 万 t 出 现 在
F MSY =0.5r (5) 2016 年。

中国水产学会主办 sponsored by China Society of Fisheries https://www.china-fishery.cn
3