Page 137 - 《水产学报》2026年第2期

2 期 水 产 学 报 50 卷

涠洲岛和斜阳岛附近渔场 (B D 5 层)、涠洲岛西南渔 min) 归一化处理，并使用 Kriging 插值法预测未调
场 (C D 5 层) 共计 3 层 (方案 D5)；⑥按 30′×30′将海 查区域的资源密度。本实验按照 5′×5′划分最小调
域划分为 12 个调查渔区，在每个渔区内部设计若 查单元，去除不能调查的区域，共计产生 341 个
干站位 (方案 D6)。 潜在站位 (图 2)。
本研究拟投入模拟的站位数分别为 9、12、 C
D = (1)
15 个。站位数分配原则按随机分配 (D1)、比例分 aq
[25]
配 (D2~D5)、配额分配 (D6) 进行 ，各方案站位
D norm = D− D min (2)
数设计见表 1。 D max − D min

2
表 1 站位设计方案和站位数分配情况 式 (1) 和 (2) 中，D 为资源密度 (kg/km )，C 为平
均每小时拖网渔获量 (kg/h)，a 为网具每小时实际
Tab. 1 Sample size allocation for station designs
2
拖网面积 (km /h)，q 为可捕系数，取值范围为
站位设计 潜在调查 拟定采样
station 层 站位数 比例/% 站位数 size 0~1，表示渔业资源总量与可能捕捞量的比率
stratum percentage
design sample size 9 12 15 (本研究取 0.5)，D m 表示归一化后的资源密度，
D1 341 100.00 9 12 15 nor
D mi n 和 D ma x 为资源密度最小值和最大值。
D2 A D2 129 37.83 3 5 6
1.3 模拟流程
141 41.35 4 5 6
B D2
71 20.82 2 2 3 本研究假设历次渔业资源调查能够反映北部
C D2
D3 A D3 45 13.20 1 2 2 湾东北部重要产卵场 4 种生物类群资源密度情况，
98 28.74 3 3 4 并基于上述假设从站位设计和抽样努力角度开展
B D3
198 58.06 5 7 9
C D3 模拟优化研究。模拟流程：第 1 步，根据调查航
D4 A D4 143 41.94 4 5 6 次计算每种类群的资源密度，对于未调查区域使
85 24.93 2 3 4
B D4 用 Kriging 插值法进行预测，由此得到每种类群单
113 33.14 3 4 5
C D4 个航次的资源密度空间分布，取该航次每种类群
D5 A D5 79 23.17 2 3 3 的平均资源密度作为“真值”；第 2 步，拟投入调
156 45.75 4 5 7
B D5 查的站位数分别设置为 9、12、15 个，每种情景
106 31.09 3 4 5
C D5 均使用 D1~D6 共 6 种站位设计进行 1 000 次的模
D6 A 23 6.74 1 1 拟调查，每模拟 1 次计算 1 次平均资源密度，使
B 36 10.56 1 2 用相对估计误差 (REE) 和相对偏差 (RB) 进行准确
C 31 9.09 1 1 度和精确度估算；第 3 步，将第 2 步依次遍历 6
D 21 6.16 1 1 个航次和所有类群，并且重复抽样 1 000 次，每
E 14 4.11 1 1 重复 1 次计算 1 次 REE 和 RB，绘制 REE 和 RB
F 36 10.56 1 2 箱线图；第 4 步，根据 REE 和 RB 的分布情况，
G 36 10.56 1 2 分析不同站位设计对鱼类、虾类、蟹类和头足类
H 17 4.99 1 1 资源密度估值的影响，并使用设计效应 (deff) 分
I 22 6.45 1 1
析 D2~D6 相对于 D1 的调查效率。
J 36 10.56 1 2
数据处理、模拟抽样和图形绘制使用 R 统
K 36 10.56 1 2
计软件 (version 4.4.1) 和 ArcGIS Pro 软件 (version
L 33 9.68 1 1 [27-29]
3.1.4) 完成。REE、RB 和 deff 的计算公式 ：
注：D1~D6表示方案D1~方案D6。 √
∑ ) 2
Notes: D1-D6 denotes design 1 - design 6. R ( N estimated − N true /R
r=1 ir M
REE = (3)
以历次调查数据为基础，使用扫海面积法估 N true
M
2 [26]
算单个航次每种生物类群的资源密度 (kg/km ) 。 ∑ R
N estimated − N true
此外，为减少不同生物类群资源密度数值过大所 r=1 ir M
RB = true (4)
造成的影响，将资源密度进行最大最小值 (max- N M
https://www.china-fishery.cn 中国水产学会主办 sponsored by China Society of Fisheries
4