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

4 期 水 产 学 报 50 卷

表现，已在识别更稳健、更符合管理目标的管理 8 月 (后文将 2014 和 2015 年共 3 次调查统称为
措施中应用 ，其核心优势在于系统考察关键不 2015 年调查) 和 2024 年 5 月、8 月、12 月在渤海
[7]
确定性如何影响管理结果，从而降低“参数设定偏 开展的底拖网调查 (图 1)。所用网具的网口高
[8]
差”导致的决策风险 。因此，针对渤海口虾蛄这 6.0 m、宽 22.6 m，主网网目 63.0 mm，囊网网目
一不确定性来源复杂的资源 ，采用 MSE框架对 20.0 mm。每个站位的平均拖网速度为 3 kn，时间
[6]
管理方案进行稳健性检验，具有明确的方法学价 为 1 h。为保证口虾蛄渔业生物学研究的一致性及
值与现实需求。R 语言 DLMtool 程序包的开发， 结果的可比性 [19, 23] ，每个站位随机选取 25 尾个体，
使得原本高门槛、定制化、难复现的 MSE，变成 测定体长 (cm) 和体重 (g)；样本量不足 25 尾时全
了更标准化、可复制、可快速迭代的一套工作流 部测量。2015 年和 2024 年分别累计测量口虾蛄
[9]
程，尤其适合数据有限的 (data-limited) 渔业场景 ， 样品 1 346 尾和 1 719 尾。
在我国近海海域也有较多的应用该方法的案例 [10-13] 。
N
自然死亡系数 (coefficient of natural mortality, 40.5° N
M) 是渔业资源评估与管理决策中最关键却也是最
难可靠估计的生活史参数之一，M 是连接生活史 39.5°
特征与种群更新和补充的关键参数，直接影响种 渤海
Bohai Sea
群生产力判断和参考点估计等。然而在多数渔业 38.5° 黄海
Yellow Sea
种群中，M 难以通过直接观测获得，常依赖经验 调查站位
公式估算 [14-18] 。研究表明，不同的评估方法之间 37.5° survey stations
可能存在显著差异，且在数据有限时，单一方法
118° 119° 120° 121° 122° E
给出的 M 可能带来不可忽视的不确定性，从而影 图 1 渤海底拖网调查站位分布图
响渔业评估与管理建议 [14-15] 。 Fig. 1 Survey stations of bottom trawl in the Bohai Sea
围绕口虾蛄的生长、死亡与资源利用状况，
国内已累积了一定研究。在黄、渤海海域，我国 本研究获得了中国水产科学研究院黄海水产
学者已报道口虾蛄的生长参数、死亡系数、50% 研究所动物实验伦理审查委员会批准 (审批号：
性成熟体长 (L m50 ) 等 [19-22] 。此外，刘逸文等 [23] 利 YSFRI-2026009)，实验过程中严格遵守中国水产
用不同方法估算了口虾蛄 M 的取值范围。同时， 科学研究院黄海水产研究所伦理规范，并按照该
也有评估报道口虾蛄资源状态，提示渤海部分海 所动物实验伦理审查委员会制定的规章制度执行。
[24]
域存在口虾蛄过度捕捞现象 。
1.2 参数估算
然而，针对渤海口虾蛄的现有研究中，关于
不同 M 估算值对管理决策影响的研究较少，使得 口虾蛄体长-体重关系采用幂函数模型描述：
b
管理建议在不确定性条件下的风险评估受到限制。 W=a·L 。式中，W 为体重 (g)，L 为体长 (cm)，a
基于此，本研究目的：①评估渤海口虾蛄的适宜 为肥满度，b 为异速生长系数 [25] 。生长参数由季
管理措施；②分析不同 M 估算方法对管理措施稳 节性 von Bertalanffy 生长方程 (soVBGF) 估算，其
健性的影响。进而识别在不同 M 设置情景下表现 公式：
{ }
更为稳健的管理方案，旨在为渤海口虾蛄资源可 −K(t−t 0 )− CK sin[2π(t−t s )]
L t = L inf 1−e 2π
持续开发和科学管理提供可量化、可解释且面向
不确定性的决策支持工具，也为数据有限渔业在 式中，L 为 t t 龄时的体长，L in f 为渐近体长，K 为
关键生活史参数存在不确定性时的 MSE 应用提供 生长系数，t 为理论初始年龄 (本研究取 0) ；C
[19]
0
参考。 为季节振幅参数，t 为季节振荡的起始点 [26-27] 。
s
以 1 cm 间隔构建体长频率数据 (LFD)，再利用 R
1 材料方法
语言 TropFishR 和 fishboot 包，通过基于 Bootstrap
重 抽 样 的 遗 传 算 法 ELEFAN (Genetic Algorithm-
1.1 样品采集
based electronic length frequency analysis, ELEFAN
数据来源于 2014 年 10 月以及 2015 年 5 月、 GA) 估算 L 、K 等参数及其不确定性范围 [28-29] 。
inf
https://www.china-fishery.cn 中国水产学会主办 sponsored by China Society of Fisheries
2