中国毛虾(Acetes chinensis)(简称“毛虾”), 为浮游性小型虾类, 分布于我国河口港湾附近水深不超过20 m的泥沙底质浅海, 主要产地为江苏、山东、浙江、福建等沿海海域, 是我国重要的海洋经济渔业资源。我国毛虾捕捞产量自20世纪50年代以来持续上升, 2006年达到最高72万t, 此后开始急剧下降, 2018年降至42.5万t (农业农村部渔业渔政管理局等, 2006, 2019)。2020年, 我国首次实行毛虾限额捕捞, 于江苏试点。

船位监控系统(vessel monitoring system, VMS)是一种集全球卫星定位、电子地图、电子海图、计算机网络通讯和数据库技术的一体化应用系统(季民等, 2005), 记录了渔船船位、航速、航向等实时动态, 已被广泛应用于海洋渔业领域。目前, 我国共计6万余艘近海机动捕捞渔船安装了北斗船载终端, 将为我国海洋渔业管理提供极具应用前景的大数据平台; 尤其现阶段我国渔捞日志的真实性、及时性、有效性无法完全得到保障的情况下(陈静娜等, 2018; 唐议等, 2020), VMS数据无疑提供了一个全新的技术途径。

自20世纪以来, 国内外学者针对拖网、刺网、张网、灯光罩网等代表性作业渔船的VMS数据, 利用神经网络、空间聚类等方法挖掘并开展船位点状态(Lee et al, 2010; Vermard et al, 2010; 张胜茂等, 2015a, b; Wang et al, 2015; Zhang et al, 2016)、航程(Wang et al, 2016)、航次(张胜茂等, 2018; 汤先峰等, 2020)以及网次(张荣瀚等, 2019)等分析, 这为我国全面实施限额捕捞奠定了基础。管控捕捞努力量等要素是限额捕捞执行的关键, 而迄今为止, 尚未涉及VMS数据在限额捕捞中的直接应用研究。本研究基于中国毛虾限额捕捞试点, 对VMS数据进行解译、提取分析, 以期为针对性地解决我国当前单品种限额捕捞试点难点积累经验。

1 数据与方法 1.1 毛虾限额捕捞许可

根据《农业农村部关于2020年伏季休渔期间特殊经济品种专项捕捞许可和捕捞辅助船配套服务安排的通告》, 毛虾限额捕捞试点海域为江苏连云港海域; 限额捕捞时间为2020年6月15日—7月15日; 毛虾捕捞渔船全长不小于24 m、主机功率不小于110 kW; 专用网具为多锚混合张网(以下简称“毛虾网”), 属张网类, 其规格为“每片网长60 m、网高18 m, 单船最多携带25片网具”。

1.2 数据基础

北斗VMS船位点数据来源于北斗民用分理服务商, 由北斗渔船终端定时发送, 时间分辨率一般为3 min, 空间分辨率约10 m, 主要包含定位时间、经纬度、航速以及航向信息等。本研究以基于北斗渔船监控系统所获取的62艘特许捕捞毛虾网船于2020年6月15日—7月15日间VMS数据作为研究对象, 剔除毛虾网船出港前和返港后航速为零值的数据后, 共421700条有效船位点数据作为分析使用。

1.3 渔船状态划分

毛虾网一般在白天从事捕捞作业, 渔船每个航次可分为出港、捕捞作业、停航休息、转场探寻以及返港等5个阶段, 其中捕捞作业阶段再可细分为抛锚布网、等待渔获以及起网收渔获等3种作业状态(表 1)。以下是对5个阶段的详细叙述: 毛虾作业渔船通常以大于10.0 kn航速到达作业渔场。随后渔船以航速3.0—6.0 kn进行抛锚布网, 由于毛虾网的抛锚和布网是高度关联的连续行为, 故划分为同一状态, 持续时间约0.5—1.0 h; 抛锚布网完毕, 渔船处于等待渔获状态, 航速为0.0—0.6 kn, 等待渔获时长约1.0—3.0 h; 然后进入起网收渔获阶段, 航速为0.6—2.0 kn, 持续时间约为2.0—3.0 h。在捕捞区域出现毛虾渔获量无法维持捕捞成本及利润的情况, 渔船选择转场探寻或者结束捕捞返港。本文随机选取某特许捕捞渔船(以下简称“渔船111”)的10天VMS数据, 对表 1划分的渔船状态对应航速进行验证, 图 1以渔船111在6月22日0:00—24:00的航速变化情况为例进行说明。根据毛虾网船作业特点并结合航行日志划分了14段(图 1), 即: 处于抛锚布网状态的S_P段、等待渔获状态的S_D段、起网收渔获的S_S段、航行状态的S_H段以及停航休息状态的S_T段, 其各状态渔船对应的航速情况与表 1相符。

1.4 基于缓冲区、叠加分析的航次提取

由于北斗船位数据具有高精度的时空分辨率, 可通过分析港口海岸线与船位数据的空间关系进行渔船航次提取(张胜茂等, 2016b)。本研究使用ArcGIS 10.2软件以毛虾限额捕捞指定渔获物定点上岸渔港“燕尾港”及附近海岸线为基准对外缓冲1 nmile, 使缓冲线处在毛虾限额捕捞区域边界线与港口海岸线之间, 其与渔船轨迹线叠加提取的交点, 即为渔船的从渔港出发或返回渔港的位置, 从而进行渔船航次提取。

1.5 基于DBSCAN聚类算法的网次提取

DBSCAN (Density-based Spatial Clustering of Application with Noise)算法是最早由Martin等提出的一种基于密度的空间数据聚类算法(Ester et al, 1996)。该算法通过ε邻域和邻域密度阈值MinPts两个参数分析样本点分布的紧密程度, 进而得到聚类数(伍育红, 2015; 原作辉等, 2020)。其中, ε邻域为样本集中与检索点的距离不大于ε值的区域; MinPts表示在检索点ε邻域内的样本最小数量(罗丹等, 2017)。由于同一艘毛虾网船在不同时间段可能在同一区域连续多次抛锚布网进行捕捞, 仅仅基于空间密度进行聚类, 会导致提取网次比实际网次少, 因此使用DBSCAN算法提取网次时需考虑一般情况毛虾网1网次捕捞作业抛锚布网所耗时长, 即时间阈值T。

以毛虾网船1个完整的作业网次(包括抛锚布网、等待渔获、起网收渔获)来设定MinPts、ε邻域, 和时间阈值。渔船抵达作业渔场后首先进行抛锚布网, 北斗时间分辨率为3 min, 检索点3 min前后应各包含一个抛锚布网点, 即MinPts为3; 依据抛锚布网状态下的船速最大阈值3 m/s与时间分辨率的乘积540 m作为ε邻域的值; 渔船1个网次持续时长约为60 min, 即相邻网次抛锚布网点集间的时间间隔在60 min以上, 时间最小阈值为60 min。因此, 本研究设定ε=540 m、MinPts=3、时间最小阈值T=60 min进行网次提取。

1.6 捕捞时长、捕捞努力量的提取与产量分析

VMS船位点数据包含定位时间、航速及航向等信息, 由于毛虾网船不同作业阶段船位点的航速及状态持续时长之间显著的差异, 基于式(1)运用数据库查询并设定相应状态阈值, 结合毛虾捕捞仅于白天作业的属性及提取得到网次数据等特征参数进行综合判断, 提取渔船捕捞阶段中各状态的船位点。

式中, V_min和V_max为渔船处于某状态对应的速度最小阈值和最大阈值; T_min和T_max是渔船处于该状态的持续时间最小阈值和最大阈值。

结合表 1, 首先设置航速阈值为1.5 m/s ≤v < 3 m/s, 时间阈值为30 min≤t < 60 min, 提取抛锚布网状态的点集; 随后设置航速阈值为0 m/s≤v < 0.3 m/s, 时间阈值为60 min≤t < 180 min, 提取等待渔获状态的点集; 最后设置航速阈值0.3 m/s≤v < 1 m/s, 时间阈值为120 min≤t < 180 min, 提取起网收渔获状态的点集。以毛虾网船各网次的第一个处于抛锚布网状态的船位点数据作为捕捞起始点, 并以该点的定位时间作为该网次的捕捞开始时间; 以各网次的第一个处于起网收渔获状态的船位点数据作为捕捞结束点, 并以该点的定位时间作为该网次的捕捞结束时间。毛虾网船各网次的捕捞时长计算方法如式(2)所示。

式中, E_k和B_k分别为毛虾网第k网次的捕捞结束时间和捕捞开始时间, T_k为第k网次的捕捞时长(h)。

一段时间内某艘毛虾网船捕捞作业分为多个网次, 每个网次的捕捞时长不同。因此将该船每个网次的捕捞时长累加, 得到的累计捕捞时间与网口迎流总面积的乘积(m²·h)即是该船这期间内在该海域的捕捞努力量, 如式(3)所示。

式中, T_k为该毛虾网船第k网次捕捞时长(h); S_k为该船第k网次网口迎流总面积(m²); L_k为第k网次有效布网距离(m); H为毛虾网具高度(m); E为该渔船总捕捞努力量(m²·h)。

核密度分析主要用于计算在指定邻近区域范围内点、线要素测量值的单位密度。本文基于计算得到的各网次捕捞产量数据, 运用核函数(式4)进行分析, 能直观反映出各网次捕捞产量在捕捞区域的分布情况。

式中, r为捕捞产量核分析查找半径, Scale为栅格中心点到网位点的距离与查找半径的比例。

2 结果 2.1 渔船航次及网次提取

将缓冲线与渔船轨迹线叠加提取到交点, 按时间顺序排列后即每相邻2个交点对应渔船1个航次, 得到所有毛虾网船在2020年6月15日—7月15日期间共239个航次, 与渔船捕捞日志的航行记录完全吻合。运用DBSCAN算法进行密度聚类, 得到所有渔船在该段时间内捕捞作业共1942网次, 平均每天捕捞作业64.7网次。

2.2 渔船网位及捕捞时长提取

以渔船111为例, 利用DBSCAN算法提取出网次数据后, 可依据该渔船处于抛锚布网、等待渔获以及起网收渔获这3种不同状态下航速以及持续时间的差异提取相应状态(图 2), 可见处于毛虾捕捞作业状态的船位点呈现直线分布形态。

经数据库查询提取所有渔船捕捞作业各个状态点集后, 采用平均中心算法对处于抛锚布网状态的点集经纬度进行计算, 处理得到的坐标位置为该网次的网位。毛虾网船各网次捕捞时长为该网次第一个起网收渔获点定位时间与第一个抛锚布网定位时间的差值, 结果表明, 82.4%网次捕捞时长为1.5—3.5 h (图 3)。

2.3 捕捞努力量提取

通过式(2)、式(3)计算得到62艘毛虾网船总捕捞时长为4413.73 h, 总捕捞努力量为108106343 m²·h。不同渔船不同网次捕捞努力量数值上差异明显, 最小值为2246 m²·h, 最大值为173283 m²·h, 平均捕捞努力量为55667 m²·h。

2.4 毛虾捕捞产量计算

鉴于限额捕捞制度对渔船功率、网具数量、网目尺寸都有严格规定, 这里渔船性能参数不作为影响渔获量的指标。因此, 将毛虾限额捕捞入渔船只的渔获量水平分为高、中、低3个层次, 分别从3个水平中选取1艘船作为样本船。本文依据3艘样本渔船的总产量数据及捕捞努力量数据(表 2), 计算出单位捕捞努力量渔获量, 并与62艘船总捕捞努力量相乘, 最终得到62艘毛虾网船捕捞总产量计算值约为2328 t, 比其实际总捕捞产量2067 t高出12.6%。通过对所有网次的毛虾捕捞产量进行核密度分析, 得到其产量空间分布如图 4所示, 毛虾资源在空间上呈现多核心分布模式, 主要集中在A、B、C区域附近。

3 讨论

毛虾网船属张网类作业, 本研究毛虾网船各阶段及状态的划分, 是依据渔船作业习惯使用人工经验进行划分。毛虾网船1个完整的航次包含渔港出发、捕捞作业、停航休息、转场探寻以及返回渔港等5个阶段, 包含了快速航行、抛锚布网、等待渔获、起网收渔获以及停航休息等5种船位点状态。然而不同网具类型的船位点状态识别时参考的要素不同, 我国学者张胜茂等(2016a)在识别拖网船位点状态时主要参考航速及航向信息; 裴凯洋等(2021)依据航速、航向、向量偏转角及距离变化信息判断张网船位点状态。本文毛虾网在网具结构上及作业方式上与其他张网类(如帆张网、近岸小型定置张网等)存在较大差异, 是毛虾捕捞专属网具。由于毛虾网船航行捕捞作业中各个阶段状态之间航速、作业时间以及持续时长存在明显的差异, 因此主要依据以上3个要素进行多参数考量从而结合相关算法语句进行船位点状态的识别。加之, 对毛虾网船的海上实际作业情况验证发现渔船在低航速状态下航向数据均会随海流发生较大变化, 因此航向数据不能作为毛虾网船状态识别的参考要素。渔船在海域中的航速、航向均会受到海流情况、天气情况等其他因素的影响, 因此快速航行过程中会出现的部分零速及低速的特殊船位点, 需进行剔除。而基于渔船航行轨迹与渔港对外设置缓冲线叠加提取航次方法的准确性取决于卫星获取渔船船位数据时空分辨率的高低。提取不同捕捞方式渔船航次所需设定的沿岸缓冲线距离不同。张胜茂等(2016b)提取拖网航次时对中国海岸线向外设置3 nmile缓冲线, 使其包含附近的岛屿及港口, 同时确保该缓冲线处于机动渔船底拖网禁渔区线之内。本研究采用高精度北斗船位数据提取渔船航次能为毛虾渔业管理监控提供精准的数据参考, 以指定的渔获物上岸渔港“燕尾港”及附近海岸线向外设置1 nmile缓冲线距离, 使其处在限额捕捞边界线与“燕尾港”附近海岸线之间。因此, 不同捕捞方式渔船航次的提取, 港口附近海岸线向外缓冲距离的设定需确保缓冲后的线在该区域沿岸港口与可捕渔区线之间。

目前学者在量化不同类型渔船的捕捞努力量时, 采用的计算方式不同, 如拖网渔船采用累计捕捞时间与渔船功率的乘积作为捕捞努力量(张胜茂等, 2014); 张网渔船采用固定网口数对应的网具迎流面积与捕捞时长的乘积作为捕捞努力量(裴凯洋等, 2021)。实际上, 张网渔船在进行捕捞作业过程中, 不同网次布网的网具数量实际上会随着渔民对区域资源量预估、潮汐状况及时间的变化而变化, 那么采用固定网口数或网长的网具迎流面积计算捕捞努力量会产生较大偏差。本研究在以往张网VMS数据提取分析上改进优化, 精确到各网次实际布网的网具数量及有效布网网长: 通过计算毛虾网各网次最后一个布网点与第一个布网点的实际布网距离; 再以布置全部网具所产生的距离及其对应有效布网网长之间的相关关系进行转换, 即可得到该网次的有效布网网长; 最后依据式3计算得到该网次的捕捞努力量。此基于北斗船位数据提取各网次实际有效布网网长方法可应用于张网类有效布网网长及捕捞努力量统计纠偏。

近年来, 随着海洋生态文明建设的持续推进, 渔业限额捕捞制度逐渐提上日程。限额捕捞制度中总许可渔获量(TAC)的设定依据主要为现有资源量及渔业生产监测情况。因此, 获取准确的渔业资源数据与生产监测数据是实施限额捕捞制度的关键。然而目前我国渔捞日志、捕捞努力量及渔获量数据主要依据人工记录, 导致生产监测数据存在延时及准确性差等不足。渔船监控系统快速发展, 加之北斗VMS数据所具备的实时性强、全天候、自动化获取等优势, 使得VMS数据将有效服务于渔业精细化管理, 成为限额捕捞制度实施的重要支撑。目前, 我国已在多个省份实施海洋渔业单品种限额捕捞试点, 如梭子蟹、丁香鱼、海蜇、毛虾等。其中, 毛虾兼具经济价值高、汛期相对固定、海上运输时间短、渔获定点上岸易控制、加工链稳定、产品输出渠道单一等鲜明特点, 可作为限额捕捞的成功模式加以推行。本文基于北斗VMS数据采用核密度分析法提取到毛虾网船捕捞产量空间分布情况(图 4), 可见渔船捕捞毛虾资源空间分布主要呈现多核心分布模式。在不同捕捞网次作业时, 渔民以获取最大毛虾渔获量为原则, 通过探鱼器和测深仪获取毛虾分布数据, 确定下一网次位置。主要继而可探究毛虾资源分布与海水温度、盐度及潮流之间的关联性及输送机制, 也可通过叠加全国所有毛虾网船捕捞努力量及产量时空分布数据及环境调查数据, 建立毛虾资源数据分析系统, 掌握毛虾渔场时空变化规律和资源量变化, 为后续毛虾单品种捕捞限额设定提供更为精准的评估, 并为全国其他省份毛虾限额捕捞试点实施提供参考。

本文利用样本船捕捞产量推算出毛虾单位捕捞努力量渔获量, 计算62艘渔船的捕捞产量值比实际产量值高12.6%, 偏高的原因可能与作业渔船渔捞日志上报的毛虾渔获量数据偏低有关。目前我国在限额捕捞制度实施监管方面主要依据与渔获物定点上岸、观察员上船、渔民记录渔捞日志等配套制度执行; 然而, 我国渔船数量众多, 仅靠观察员上船及渔民记录上报捕捞情况进行监管, 造成观察覆盖不到位、渔民上报不及时甚至作假等现象频出, 难以达到监管效果。那么, 加强渔业捕捞监管系统的建设, 充分利用北斗船位监控系统平台, 借助VMS大数据定位作用监测限额捕捞许可渔船的航行轨迹是否在许可作业海域、渔船作业方式和投网数量是否符合要求、渔获物是否定点上岸, 并与后期产品溯源进行无缝链接, 将成为限额捕捞制度持续顺利推进的重要保障。

4 结论

本研究以中国毛虾限额捕捞管理为例, 提出基于北斗船位数据提取航次、网次、网位及捕捞时长信息用于计算海洋捕捞努力量的方法, 计算了限额捕捞期间毛虾渔获总量, 并与生产上报产量进行对比, 形成了VMS数据在限额捕捞试点的直接应用。该方法和思路一定程度上解决了我国渔船渔捞日志人工记录存在的效率低下及真实度欠缺等问题, 为后期制定单品种限额捕捞总量和捕捞配额提供了科学依据。