辽宁省海洋渔业产量结构调整的初步研究

尹增强, 张国胜, 李九奇, 冯维山

尹增强, 张国胜, 李九奇, 冯维山. 辽宁省海洋渔业产量结构调整的初步研究[J]. 南方水产科学, 2005, 1(6): 55-62.
引用本文: 尹增强, 张国胜, 李九奇, 冯维山. 辽宁省海洋渔业产量结构调整的初步研究[J]. 南方水产科学, 2005, 1(6): 55-62.
YIN Zengqiang, ZHANG Guosheng, LI Jiuqi, FENG Weishan. The preliminary study on the structure regulation of marine fishery yield in Liaoning province[J]. South China Fisheries Science, 2005, 1(6): 55-62.
Citation: YIN Zengqiang, ZHANG Guosheng, LI Jiuqi, FENG Weishan. The preliminary study on the structure regulation of marine fishery yield in Liaoning province[J]. South China Fisheries Science, 2005, 1(6): 55-62.

辽宁省海洋渔业产量结构调整的初步研究

基金项目: 

辽宁省教育厅A类项目 20104330

详细信息
    作者简介:

    尹增强(1973-), 男, 硕士, 讲师, 从事渔业资源开发与管理的研究工作。Email: zqyin@163.com

  • 中图分类号: S932

The preliminary study on the structure regulation of marine fishery yield in Liaoning province

  • 摘要:

    根据辽宁省海洋渔业统计资料,分析了辽宁省海洋渔业产量结构,运用剩余产量模型估算了最大持续产量,运用ARIMA(p, d, q)时间序列模型与GM(1, 1)模型分析预测了2004~2010年海水养殖产量、远洋渔业产量与各海区的捕捞产量。结果表明:辽宁省近海捕捞产量已经超过最大持续产量,如不采取有效措施仍将有一定程度的增长,在渤海与黄海的产量到2010年将比2003年增加23.40万t与23.79万t;远洋渔业与海水养殖业,其产量在2010年将比2003年增加8.59万t与242.40万t。文章对辽宁省海洋渔业产量结构调整进行了讨论。

    Abstract:

    Based on marine fishery statistics data in Liaoning province, an analysis of marine fishery yield structure was made, the surplus yield model was used to calculate the maximum sustainable yield, and the forecasted yield of marine culture, pelagic fishing and each domestic fishing sea over the period from 2004 to 2010 was calculated by using autoregressive integrated moving average model [ARIMA(p, d, q)] and grey model theory[GM(1, 1)]. The result indicated that offshore fishing yield in Liaoning province, which had overrun maximum sustainable yield, would increase on the certain degree over the period from 2004 to 2010 if some effective measures wasn′t be adopted. Compared with the number of 2003 in Bohai Sea and Yellow Sea, the marine fishing yield would increase by 234.0 thousand tons and 237.9 thousand tons. Furthermore, compared with the number of 2003 in mariculture and pelagic fishing, the yield would increase by 85.9 thousand tons and 2 424.0 thousand tons. And a discussion about the structure regulation of marine fishery yield was made.

  • 刀额新对虾(Metapenaeus ensis)隶属于节肢动物门、甲壳纲、十足目、对虾科、新对虾属[1-2]。新对虾属是南海海域经济虾类中最重要的属,既是捕捞种类,也是养殖品种[3]。刀额新对虾是南海海域新对虾属中重要的经济种,资源量丰富、分布广、繁殖期长[2],是南海海域虾拖网的主要捕捞对象,渔汛旺期其产量占虾拖网虾类总产量的50%以上[4-5]。由于虾拖网网囊网目尺寸小、选择性差,导致捕获的刀额新对虾幼虾比例较高,有时甚至高达90%[5],严重影响虾资源的合理利用。因此,有必要对南海海域虾拖网的网囊进行优化设计,以降低对刀额新对虾幼虾的捕捞比例,减少捕捞作业对渔业资源补充群体的不良影响。

    明确拖网网囊对主要捕捞种类的选择性是渔具管理、渔具优化和资源利用的基础[6-7]。目前,已对南海海域刀额新虾进行过食性、营养级和群体形态学等方面的专题研究[8-9];关于网囊对刀额新对虾的选择性有过少量报道[10-12],但是关于虾拖网网囊网目对刀额新对虾选择性的系统、专题研究尚未见报道。本文在整理和总结南海海域虾拖网网囊系列选择性试验数据的基础上,专题开展网囊对刀额新对虾的选择性研究,为南海海域刀额新对虾资源的合理利用和虾拖网渔具的设计优化提供科学依据。

    试验于2014—2017年分4个航次进行(表1),船名分别为“粤阳东渔12057”(船长16 m、主机功率79 kW)和“粤阳东渔12081”(船长21 m、主机功率98 kW)。试验区域为南海区近海渔场,经纬度范围为112°03'E~112°22′E、21°37'N~21°50′N,水深6~13 m,底质为沙泥。

    表  1  网囊选择性试验基本信息
    Table  1  Basic information of tested codends
    航次
    trial
    时间
    time
    网囊
    codend
    试验渔船
    vessel
    方法
    method
    网目内径/mm
    mesh opening
    网囊规格 (圆周×纵向)
    codend specification
    (circumferernce×vertical)
    有效网次
    valid haul
    12014.08D25粤阳东渔 12057YYDY12057套网法 covered21.10±0.9358×606
    12014.08D30粤阳东渔 12057YYDY12507套网法 covered26.64±0.7348×508
    22015.08S35+D18粤阳东渔 12081YYDY12081套网法 coveredS: 32.70±0.63
    D: 15.07±0.36
    S: 29×43,D: 80×428
    32016.08S25+D25粤阳东渔 12081YYDY12081套网法 coveredS: 21.80±0.42
    D: 21.80±0.42
    S: 40×60,D: 58×3011
    32016.08S30+D25粤阳东渔 12081YYDY12081套网法 coveredS: 27.14±0.45
    D: 21.80±0.42
    S: 35×50,D: 58×3012
    32016.08S35+D25粤阳东渔 12081YYDY12081套网法 coveredS: 32.42±0.18
    D: 21.80±0.42
    S: 29×43,D: 58×3010
    42017.09S35+D25粤阳东渔 12081YYDY12081裤网法 trouserS: 32.42±0.18
    D: 21.80±0.42
    S: 23×29,D: 81×2011
    42017.09S35+D30粤阳东渔 12081YYDY12081裤网法 trouserS: 32.42±0.18
    D: 27.14±0.45
    S: 23×29,D: 67×1710
    42017.09S35+D35粤阳东渔 12081YYDY12081裤网法 trouserS: 32.42±0.18
    D: 32.42±0.18
    S: 23×29,D: 58×148
    合计 total84
     注:S. 方形网目;D. 菱形网目
     Note: S. square mesh; D. diamond mesh
    下载: 导出CSV 
    | 显示表格

    试验网具为南海区最常见的双桁杆虾拖网(以下简称虾拖网),其原主尺度为10.50 m×7.15 m(2.2 m)。关于虾拖网的结构特征、渔法特点和渔获组成等信息,可参阅相应的参考文献[5,13-14]

    2014—2017年,采用套网法和裤式拖网法(简称裤网法)依次对2组菱目网囊和6组菱目与方目混合网囊(简称混目网囊)进行了84个有效网次作业试验(表1)。套网法的试验网囊拉直长度约为1.5 m;裤网法试验网囊和对照网囊的拉直长度约为1.2 m。

    试验网囊的网目尺寸分4档:18 mm、25 mm、30 mm和35 mm。其中18 mm是南海海域虾拖网普遍使用的网囊网目尺寸;25 mm为现行虾拖网管理规定中的网囊最小网目尺寸;30 mm和35 mm为试验网网囊尺寸。按照网目尺寸对网囊进行编号,例如D30表示网目尺寸为30 mm的菱目网囊;S30+D25表示方目段网目尺寸为30 mm,菱目段网目尺寸为25 mm的混目网囊。套网法试验中套网的网目尺寸为15 mm,内径为(12.55±0.28) mm,圆周286目,纵向136目。裤网法中对照网囊第一段网目尺寸为24 mm,第二段为22 mm;对照网囊第二段网目内径为(19.49±0.63) mm。试验网具和网囊的具体参数详见表1图1图2

    图  1  套网法网囊选择性试验渔具网衣展开图
    Fig. 1  Schematic diagram of trawl and tested codends with covered codend method
    图  2  裤网法网囊选择性试验渔具网衣展开图
    Fig. 2  Schematic diagram of trawl and tested codends with trouser trawl method

    为了方便操作,将试验网设置于渔船船首方向右舷第一顶网处。试验中不设置站点,也不规定渔船的拖曳速度、时间等作业参数,由船长根据实际情况控制,使试验尽量与传统虾拖网渔船作业保持一致。试验期间渔船拖速2.2~2.7 kn,网次作业时间1.5~3.5 h,平均1.75 h,网次平均拖曳距离9.01 km。

    起网后,分别称量试验网囊和套网(或对照网囊)的渔获总质量,然后进行分类和鉴定,统计每个种类的渔获数量,并进行生物学测量。当刀额新对虾渔获数量≤50尾时,全部测量其长度;当渔获数量>50尾时,随机抽取50尾进行长度测量。在统计渔获种类的体长分布时,对于取样种类按取样比例进行加权处理。为了保证试验数据的准确性,尽量避免渔获物的2次取样。

    以航次为单位统计虾拖网网囊和套网(或对照网囊)的渔获总数和总质量,分别计算各网囊刀额新对虾的数量和质量比例。将每个网次中网囊和套网(或对照网囊)的刀额新对虾体长数据进行排序,以5 mm为单位进行分组,分别统计各体长组刀额新对虾的数量。

    网囊对刀额新对虾的选择性估算按如下2个步骤进行:1)以网次为单位,使用Logistic曲线进行选择性模型拟合,估算网次选择性参数;2)将选择性模型拟合度较好的网次体长数据进行叠加,再次使用Logistic曲线进行选择性模型拟合,估算网囊对刀额新对虾的平均选择性参数,同时考虑网囊网次间差异对选择性的影响。

    选择性参数的估算是在SELECT模型[6,15-18]的框架下,通过比较试验网囊和套网(或对照网囊)的刀额新对虾渔获数量,得到估算结果。在第i网次中,用NTij表示体长为lj的刀额新对虾被试验网囊捕获的数量,用NCij表示套网(或对照网囊)中刀额新对虾(体长亦为lj)的渔获数量,则体长为lj的刀额新对虾被试验网囊捕获的概率(φij)可表示为:

    $${{\textit{φ}} _{ij}} = \frac{{{N_{Tij}}}}{{{N_{Tij}} + {N_{Cij}}}}$$ (1)

    对于套网法试验,使用logistic方程拟合试验网囊对体长为l的刀额新对虾的选择率r(l),其计算公式为:

    $$r(l) = \frac{{\exp(a + bl)}}{{1 + \exp(a + bl)}}$$ (2)

    式中ab为待估参数。

    对于裤网法试验,引入相对作业强度(p),表示1尾刀额新对虾进入试验网囊的概率,则进入对照网囊的概率为1−p。体长为l的刀额新对虾被试验网囊捕获的概率可表示为:

    $$\varphi(l) = \frac{{p\exp(a + bl)}}{{(1 - p) + \exp(a + bl)}}$$ (3)

    根据选择性参数ab,计算试验网囊对刀额新对虾的选择性指标50%选择体长(L50)和选择范围(SR),相应的计算公式为:

    $${L_{50}} = - a/b$$ (4)
    $${\rm SR} = {L_{75}} - {L_{25}} = 2\ln 3/b$$ (5)

    通过极大似然法估算选择性参数abp,其似然方程[7]为:

    $$L(a,b,p) = \sum\limits_i {\sum\limits_j {[{N_{Tij}}\ln \varphi({l_j}) + {N_{Cij}}\ln(1 - \varphi({l_j}))]} } $$ (6)

    式(6)的求解通过MS-Excel的“规划求解”功能实现[19-20]。参数abp和选择性指标L50、SR的标准差通过Delta方法估算[6,17,21]

    选择性模型拟合度的优劣可根据残差值(deviance)和自由度(degree of freedom, df )的大小关系判断。通常当残差值小于自由度(P>0.05),则认为模型拟合度为优。但当试验数据呈离散分布时,即使拟合度为优也会出现残差值远大于自由度(P<0.05)的情况。这时需要根据各体长组的实际残差值分布图作进一步判断[6,12,22]

    通过比较联合网次模型拟合的皮尔逊卡方统计量(Q)和自由度(d)的大小来判断网囊选择性是否存在网次间差异。当不存在网次间差异时,Q服从自由度为d的卡方分布(P>0.05),反之则说明存在网次间差异[6,17,21-25]。如果存在网次间差异,计算过度离散叠加估算值(replication estimate of dispersion, REP),并对联合网次估算的选择性参数标准差进行修正,即乘以(REP)1/2Q和REP的计算公式为:

    $$Q = \sum\limits_l {\sum\limits_h {\frac{{{{(n_{lt}^h - n_{l + }^h{y_l})}^2}}}{{n_{l + }^h{y_l}{{(1 - {y_l})}^{}}}}} } $$ (7)

    式中$n_{lt}^h$表示第h网次中试验网囊捕获的体长为l的刀额新对虾数量,$n_{l + }^h$表示第h网次中试验网囊与套网(或对照网囊)捕获的体长为l的刀额新对虾总数量,yl表示h网次中试验网囊捕获的体长为l的刀额新对虾的数量比例。

    $${\rm REP} = \frac{Q}{d}$$ (8)

    式中d=(H−1)×体长组数,H为联网次的总网次数。

    试验渔获质量为294.68 kg,刀额新对虾总渔获质量为100.71 kg,占渔获总质量的34.18%。各试验网囊中的总渔获质量及刀额新对虾所占比例详见表2

    表  2  试验渔获产量基本信息
    Table  2  Basic information of fish catch
    方法
    method
    网囊
    codend
    质量/kg mass质量百分比/%
    ratio of mass
    刀额新对虾
    M. ensis
    总渔获
    total catch
    套网法 coverD258.8027.9231.50
    D3013.4345.0129.83
    S35+D189.5028.5333.28
    S25+D2515.3623.3665.77
    S30+D254.9433.6114.70
    S35+D2520.1337.6053.53
    裤网法 trouserS35+D2518.4732.4756.89
    S35+D307.1842.9116.73
    S35+D352.9223.2712.53
    合计 total100.73294.6834.18
    下载: 导出CSV 
    | 显示表格

    试验中刀额新对虾的渔获总数量为14 138尾,各网囊和套网(或对照网囊)的数量、体长范围、众数体长及其比例等数据详见表3,体长分布见图3图4

    表  3  刀额新对虾的渔获数量及体长分布
    Table  3  Number of M. ensis caught by each codend
    网囊
    codend
    渔获尾数
    number of catch
    体长范围/mm
    length range
    众数体长/mm
    modal length
    网囊
    codend
    套网/对照网
    cover/control
    网囊
    codend
    套网/对照网
    cover/control
    网囊
    codend
    套网/对照网
    cover/control
    D25964253~12333~5888~93 (38.38%)
    D301 698353~11353~6383~88 (51.94%)
    S35+D181 8102863~12863~9878~83 (37.07%)73~78 (67.86%)
    S25+D252 30738458~10328~9378~83 (48.63%)43~48 (39.84%)
    S30+D254993363~10363~9883~88 (46.49%)73~78 (42.42%)
    [S35+D25]13 00811758~11828~9878~83 (41.42%)73~78 (42.74%)
    [S35+D25]21 1211 36663~13363~13383~88 (58.25%)83~93 (65.37%)
    S35+D3028327063~13373~13393~103 (55.83%)93~103 (60.00%)
    S35+D3511313263~12868~13393~98 (47.79%)93~98 (46.97%)
    合计 total11 8032 335
     注:[S35+D25]1. 套网法中的S35+D25网囊;[S35+D25]2. 裤网法中的S35+D25网囊
     Note: [S35+D25]1. the S35+D25 codend for the covered codend method; [S35+D25]2. the S35+D25 codend for the trouser trawl method
    下载: 导出CSV 
    | 显示表格
    图  3  套网法选择性试验各网囊捕捞刀额新对虾的体长分布
    Fig. 3  Length distribution of M. ensis caught by each codend with covered codend method
    图  4  裤网法选择性试验各网囊捕捞刀额新对虾的体长分布
    Fig. 4  Length distribution of M. ensis caught by each codend with trouser trawl method

    D25和D30网囊由于套网的刀额新对虾样本量过少,无法进行选择性分析;S35+D18网囊仅1个单网次数据获得有效解,模型的拟合度较好(P>0.05),刀额新对虾的L50和SR分别为51.25 mm和16.48 mm;S25+D25网囊共5个单网次数据获得有效解,网次数据拟合度均较好(P>0.05),联合网次拟合度P<0.05,但模型残差值分布正常,所以认为是由于体长数据的过度离散分布造成,联合网次估算的刀额新对虾平均L50和SR分别为60.84 mm和14.31 mm,存在网次间差异,对估算参数的标准差进行了修正;S30+D25网囊仅2个单网次数据获得有效解,单网次和联合网次拟合度较好(P>0.05),联合网次估算的刀额新对虾平均L50和SR分别为63.21 mm和12.84 mm,不存在网次间差异;S35+D25网囊共5个单网次数据获得有效解,1个单网次和联合网次拟合的P<0.05,但模型残差值分布正常,联合网次估算的刀额新对虾平均L50和SR分别为64.53 mm和9.75 mm,存在网次间差异,对估算参数的标准差进行了修正(表4)。各网囊对刀额新对虾的平均选择性曲线详见图5

    表  4  套网法选择性试验参数估算
    Table  4  Selective parameters of tested codends with covered codend method
    时间
    time
    网囊
    codend
    网次
    NH
    选择性指标
    selective index
    选择性参数
    selective parameter
    拟合度
    goodness of fit
    网次间差异
    estimate of REP
    尾数
    number
    L50/
    mm
    SDSR/
    mm
    SDaSDbSDDdfPQdPREP网囊
    codend
    套网
    cover
    2014D25c-b9642
    2014D30c-b1 6983
    2015S35+D18951.258.9416.487.20–6.834.140.130.064.29580.8320311
    2016S25+D25259.655.8421.517.41–6.092.660.100.045.91970.5513724
    453.8912.6850.4426.94–2.351.790.040.0210.60870.1612044
    653.317.7930.9611.92–3.781.980.070.036.28390.7112731
    960.861.196.051.07–22.124.210.360.0614.094140.4436764
    1060.941.705.181.21–25.846.390.420.102.288141.0027875
    c-b60.841.7414.311.90–9.341.440.150.0251.61714<0.05155.9439<0.054.001 029238
    S30+D25461.487.8911.786.01–11.477.250.190.109.811100.461064
    563.183.4713.933.23–9.972.800.160.048.82260.1823317
    c-b63.212.9712.842.66–10.822.700.170.043.65790.9320.22210.510.9634221
    S35+D25166.161.223.560.89–40.8610.770.620.157.000110.8044316
    568.581.617.761.83–19.424.920.280.079.94660.1318817
    762.433.459.002.31–15.254.700.240.0622.1664<0.0535615
    865.193.1010.642.78–13.464.110.210.0514.34370.0527918
    971.903.5514.865.32–10.634.200.150.055.30360.514012
    c-b64.531.159.751.00–14.541.720.230.0237.4416<0.0572.31710.431.021 30678
     注:NH. 网次;c-b. 联合网次;L50. 50%选择体长;SR. 选择范围;ab. 选择性参数;SD. 标准差;D. 残差值;dof和d. 自由度;Q. 皮尔逊卡方统计量;REP. 过度离散叠加估算值.  Note: NH. number of hauls; c-b. combined hauls; L50. 50% retention length; SR. selection range; a and b are selective parameters; SD. standard error; D. value of model deviance; dof and d indicate the degree of freedom; Q. Pearson chi-square statistic; REP. replication estimation of dispersion
    下载: 导出CSV 
    | 显示表格
    图  5  套网法混目网囊对刀额新对虾的平均选择性曲线
    Fig. 5  Mean selective curves of tested codends for M. ensis with covered codend method

    S35+D25网囊共8个单网次数据获得有效解,4个单网次和联合网次拟合的P<0.05,但模型残差值分布正常,联合网次估算的刀额新对虾平均L50、SR和相对作业强度(p)分别为75.43 mm、6.93 mm和0.51,存在网次间差异,对选择性参数的标准差进行了修正;S35+D30网囊共5个单网次数据获得有效解,单网次和联合网次的拟合度较好(P>0.05),联合网次估算的刀额新对虾平均L50、SR和p分别为82.38 mm、6.39 mm和0.52,不存在网次间差异;S35+D35网囊共4个网次获得有效解,所有单网次和联合网次估算的拟合度均较好(P>0.05),联合网次估算的刀额新对虾平均L50、SR和p分别为95.39 mm、20.44 mm和0.64,不存在网次间差异(表5)。各网囊对刀额新对虾的平均选择性曲线详见图6

    表  5  裤网法选择性试验参数估算
    Table  5  Selective parameters of tested codends with trouser trawl method
    时间
    time
    网囊
    codend
    网次
    NH
    选择性指标
    selective index
    选择性参数
    selective parameter
    拟合度
    goodness of fit
    网次间差异
    estimate of REP
    尾数
    number
    L50/mmSDSR/mmSDaSDbSDPSDDdfPQdPREP网囊
    codend
    对照
    control
    2017S35+D25278.163.153.744.45–45.8853.610.590.700.410.078.26290.512443
    375.353.037.174.69–23.0814.440.310.200.460.078.96860.1855101
    473.521.403.452.58–46.7734.530.640.480.410.0339.62510<0.05149244
    678.467.907.466.76–23.1219.450.290.270.770.1010.982100.362812
    772.795.326.927.07–23.1222.960.320.320.510.069.76680.285457
    875.982.388.553.66–19.527.890.260.110.540.0447.5316<0.05244288
    973.521.474.223.61–38.3232.460.520.450.510.0329.7207<0.05182193
    1185.594.7410.743.67–17.515.220.200.070.730.0825.51710<0.05116108
    c-b75.431.366.932.14–23.937.090.321.000.510.0353.55912<0.05152.0081<0.051.888521046
    S35+D30178.846.407.688.90–22.5725.180.290.330.650.0815.36180.053020
    586.794.966.926.25–27.5624.110.320.290.500.079.579110.573948
    781.183.615.014.56–35.6031.570.440.400.440.0512.484110.336290
    888.475.7011.327.00–17.179.820.190.120.600.0814.706110.206465
    992.098.499.807.41–20.6414.180.220.170.710.145.98680.652319
    c-b82.382.016.392.72–28.3211.680.340.150.520.0314.30880.0744.21500.700.88218242
    S35+D35183.597.528.729.18–21.0621.320.250.270.550.1010.743100.382324
    292.138.379.629.67–21.0319.670.230.230.650.146.20070.522422
    3110.7831.9720.4319.45–11.918.420.110.100.770.296.39790.701719
    5100.5237.5829.7831.29–7.425.580.070.080.650.3010.246110.513647
    c-b95.3910.4320.4410.51–10.254.440.110.060.640.1111.599110.3922.73340.930.67100112
    下载: 导出CSV 
    | 显示表格
    图  6  裤网法混目网囊对刀额新对虾的平均选择性曲线
    Fig. 6  Mean selective curves of tested codends for M. ensis with trouser trawl method

    目前,南海海域传统虾拖网均使用菱目网囊作业,且其网目尺寸较小(一般<20 mm)。据相关调查,当虾拖网网囊网目尺寸为18 mm时,刀额新对虾的幼虾比例>70%[5]。本研究2014年对菱目网囊D25和D30进行试验,结果表明对于刀额新对虾而言,网目尺寸增大,网囊选择性未见显著提高。由于套网的样本量太少,上述2种菱目网囊对刀额新对虾的选择性参数也无法估算。

    鉴于菱目网囊选择性差的结论,2015年研究团队首次尝试使用混目网囊(S35+D18)进行试验;该混目网囊菱目段网目尺寸与渔船所用网囊一致(18 mm),方目段网目尺寸为35 mm。结果表明,副渔获物(鱼类)的逃逸率显著提高,但目标种类的逃逸率仍然较低;刀额新对虾的L50太小(53.85 mm),远小于南海区刀额新对虾的首次性成熟体长(80 mm)[26],网囊网目尺寸需要进一步优化[11]。本研究采用的是单网次数据处理方法,得出S35+D18网囊对刀额新对虾的L50仅为51.52 mm,小于联合网次数据处理得到的值。此结论也进一步说明S35+D18网囊的网目尺寸需要进一步放大。

    2016年使用3组混目网囊(S25+D25、S30+D25和S35+D25)进行试验,网囊菱目段网目尺寸均为25 mm,方目段网目尺寸分别为25 mm、30 mm和35 mm。根据前期联合网次分析法的结果,S25+D25和S35+D25网囊对刀额新对虾的L50分别为63.93 mm和67.96 mm; S30+D25网囊的刀额新对虾体长组数据不收敛,无法取得有效解[12]。本研究得出的S25+D25、S30+D25和S35+D25网囊对刀额新对虾的平均L50分别为60.84 mm、63.21 mm和64.53 mm。使用本研究的数据处理方法能够获得S30+D25网囊对刀额新对虾的选择性。这从侧面说明了如果在数据处理时完全忽略网囊的网次间差异,简单地将所有网次的体长组数据进行无条件叠加处理可能会得到错误的结论。本试验采用“2部曲”(two-step procedure)法:先以网次为单位进行数据处理,然后根据模型的拟合结果将一些网次数据剔除,在考虑网囊网次间差异的基础上叠加网次数据,求得平均选择性。这种数据处理方法得出的结论更加准确、客观和科学,国外一些相关研究也普遍采用这种方法进行数据处理[27-29]

    2017年,笔者课题组对3组混网囊(S35+D25、S35+D30和S35+D35)进行试验。此次网目尺寸的选择是以2014—2016年试验结论为基础,将方目段网目尺寸保持为35 mm,菱目段网目尺寸分别为25 mm、30 mm和35 mm。同时,为了避免套网法的“覆盖效应”,在试验方法上首次尝试采用目前国际上普遍使用的裤网法。根据试验结果,S35+D25、S35+D30和S35+D35网囊对刀额新对虾的平均L50随菱目段网目尺寸的增大而不断增大,SR则是先减少后增大,p不断增大。根据陈丕茂[26]的研究,南海区刀额新对虾的首次性成熟体长为80 mm。假设将刀额新对虾的可捕规格定为80 mm,所有试验网囊中仅S35+D30和S35+D35网囊的L50大于该可捕标准。但是,S35+D30网囊对刀额新对虾的SR远小于S35+D35的(6.39 mm vs 20.44 mm),说明S35+D30网囊的选择性尖锐度较好。另外,考虑到网目尺寸过大会影响到渔船的生产效益,同时兼顾渔民的心理接受程度,笔者认为将S35+D30网囊作为南海区虾拖网网囊最小网目尺寸标准较为理想。

    本研究中的试验数据均在渔船上即时收集。由于网次间间隔短,测量时间有限,对于渔获数量较多的种类(>50尾),先统计该种类的总数量后,随机抽取50尾渔获进行体长测量。选择性数据处理时,以随机测量的50尾体长数据为基础,按照相应的取样比例进行加权处理,最后得出网次渔获种类的体长分布数据。但有些前期的研究则直接使用实测体长组数据进行选择性参数估算[12,15-16]。使用实测体长组数据和加权处理后的数据进行估算的结果是否存在显著差异?哪个方法计算的结果更加准确和科学?这将是未来的一个研究方向。

    本研究存在的不足之处在于4个航次试验存在一定的时空差异。虽然网囊对刀额新对虾的选择性主要由网囊结构、网目尺寸和刀额新对虾的体型特征决定,但是为了减少时空差异对选择性研究的影响,今后可尝试将所有试验网囊在同一渔场、同一渔船和同一网次进行试验。南海区虾拖网通常是一船拖曳多顶网具。比如,本研究中租用的“粤阳东渔12081”作业中拖曳12顶网,一些大功率渔船甚至最多可拖曳32顶网。利用南海区虾拖网这种独特的作业方式,可将本研究中的所有试验网囊共同进行试验。同时,在试验方法上也可以将套网法、裤网法等进行有机结合,优势互补,相互印证。这样获得的试验数据会更加全面,得出的结论能够更好地为渔具管理和渔业资源的合理利用提供科学依据。

  • 图  1   辽宁省1961~2003年海洋渔业产量结构百分比例变化

    Figure  1.   The variation of the percent of marine fishery yield structure in Liaoning province from 1961 to 2003

    图  2   辽宁省渤海捕捞产量与模型拟合产量图

    Figure  2.   Marine fishing yield, its fitted yield and its error of Liaoning province in Bohai Sea

    图  3   辽宁省黄海捕捞产量与模型拟合产量图

    Figure  3.   Marine fishing yield, its fitted yield and its error of Liaoning province in Yellow Sea

    图  4   辽宁省东海捕捞产量、模型拟合产量与误差图

    Figure  4.   Marine fishing yield, its fitted yield and its error of Liaoning province in East China Sea

    图  5   辽宁省海水养殖产量、模型拟合产量与误差图

    Figure  5.   Mariculture yield, its fitted yield and its error of Liaoning province

    表  1   辽宁省1961~2003年海洋渔业产量与捕捞渔船功率

    Table  1   Marine fishery yield in Liaoning province from 1961 to 2003 104 t, 104 kW

    年份
    year
    海水养殖产量
    mariculture yield
    远洋渔业产量
    pelagic fishing yield
    国内海洋捕捞业domestic marine fishing area
    渔船总功率
    total power
    总产量
    total yield
    渤海产量
    Bohai Sea yield
    黄海产量
    Yellow Sea yield
    东海产量
    East China Sea yield
    1961 1.60 4.06 16.97 4.30 12.67 0
    1962 1.40 4.69 17.44 4.60 12.84 0
    1963 1.10 4.49 16.56 5.38 11.18 0
    1964 0.80 4.60 17.89 5.22 12.66 0
    1965 1.20 4.61 15.45 4.47 10.98 0
    1966 1.60 4.95 15.36 3.86 11.49 0
    1967 1.80 5.10 13.76 3.86 9.90 0
    1968 1.60 5.15 14.07 4.43 9.64 0
    1969 1.80 5.30 16.64 5.61 11.03 0
    1970 2.20 5.46 19.98 5.79 14.19 0
    1971 3.60 5.97 23.47 6.39 17.09 0
    1972 4.20 7.42 26.73 6.26 20.46 0
    1973 4.60 9.13 31.46 8.60 22.86 0
    1974 4.30 11.00 33.42 9.69 23.73 0
    1975 5.90 14.11 36.84 10.74 26.11 0
    1976 7.50 16.55 36.93 9.99 26.93 0
    1977 10.30 18.16 36.00 10.87 25.13 0
    1978 10.70 19.38 35.49 11.32 24.17 0
    1979 9.80 20.51 33.43 12.00 18.20 3.23
    1980 11.30 23.73 29.95 9.58 15.47 4.91
    1981 14.10 24.02 29.44 8.63 15.69 5.12
    1982 15.20 24.33 34.51 8.71 18.75 7.05
    1983 15.90 25.88 33.67 9.38 17.82 6.47
    1984 17.60 28.23 35.18 11.02 15.57 8.59
    1985 16.70 0.08 31.60 38.43 15.03 12.67 10.74
    1986 24.50 0.41 39.25 38.32 13.89 14.05 10.39
    1987 32.10 1.99 42.00 41.78 14.45 16.65 10.68
    1988 43.20 3.45 45.25 42.11 16.62 14.72 10.78
    1989 48.10 1.84 52.06 44.87 16.69 16.09 12.10
    1990 51.40 2.37 54.31 46.94 17.57 19.47 9.90
    1991 54.90 1.73 56.56 49.70 20.76 18.75 10.19
    1992 65.50 3.56 59.03 53.56 25.26 19.79 8.50
    1993 73.80 7.12 58.76 58.82 28.00 24.06 6.76
    1994 76.50 6.74 58.55 68.04 30.66 34.82 2.56
    1995 87.30 9.31 63.32 81.84 39.74 39.27 2.83
    1996 111.28 16.36 65.83 107.71 50.85 48.92 7.94
    1997 112.31 20.31 71.66 125.36 60.63 57.03 7.71
    1998 120.80 20.35 83.31 140.26 71.56 61.51 7.19
    1999 139.16 22.12 89.14 135.57 66.97 62.63 5.97
    2000 152.12 16.85 96.73 133.34 65.02 64.53 3.79
    2001 160.89 16.00 102.67 131.45 64.14 60.67 6.64
    2002 178.20 18.80 108.61 121.42 59.25 56.04 6.14
    2003 196.20 25.80 110.94 121.86 55.28 63.42 3.16
    注:由于辽宁省在南海的产量很低且大多数年份无产量故未考虑Note: Marine fishing yield of Liaoning province from South China Sea is not taken into account because of its low and no yield for the most years.
    下载: 导出CSV

    表  2   辽宁省国内海洋渔业最大持续产量(MSY)与相应捕捞努力量(fMSY)估算值及统计检验

    Table  2   The calculated values and statistics tests of maximum sustainable yield(MSY) and the corresponding fishing effort(fMSY) in Liaoning province

    模式model MSY (104 t) fMSY (104 kW) 参数parameter 统计检验statistics test
    Schaefer 105.85 72.75 a=2.91 b=0.02 F=37.65>F0.01=7.35
    Fox 99.40 100.00 c=2.70 d=0.01 F=35.15>F0.01=7.35
    平均mean 102.63 86.38
    下载: 导出CSV

    表  3   辽宁省国内海洋渔业产量ARIMA(p, d, q)预测模型及模型检验汇总表

    Table  3   The ARIMA(p, d, q) models and model test of domestic marine fishery yield in Liaoning province

    名称name 模型检验model test
    残差序列的白噪声检验white noise test of residual series SBC AIC 相关系数(R) correlation coefficient
    渤海Bohai Sea
    ARIMA(1, 1, 1)
    白噪声 -49.70 -54.92 0.9904
    黄海Yellow Sea
    ARIMA(1, 1, 1)
    白噪声 -41.86 -47.07 0.9839
    东海East China Sea
    ARIMA(2, 0, 1)
    白噪声 -42.56 -87.92 0.9098
    海水养殖maricuture
    ARIMA(1, 1, 1)
    白噪声 -22.01 -27.22 0.9968
    注:残差序列为白噪声时才可应用ARIMA模型预测;SBC是Schwarz贝叶斯准侧,AIC是Akaike信息准则,SBC与AIC越小,相关系数(R)越接近1,模型越好
    Note: Applying ARIMA model is based on the premise that white noise test of residual series can be passed. SBC—Schwarz Bayesian criterion, AIC—Akaike information criterion; SBC and AIC more small and coefficient correlation is more close to 1, the accuracy of model prediction is more high.
    下载: 导出CSV

    表  4   辽宁省渤黄东海与远洋渔业捕捞产量GM(1, 1)预测模型及模型检验表

    Table  4   The GM(1, 1) models and model tests of marine fishing yield in Liaoning province

    名称
    name
    GM(1, 1)预测模型
    predicted model of GM(1, 1)
    模型检验model test
    后验差比(c)
    ratio of observed data variance and error variance
    小误差概率(P)
    probability of mini-error
    预测精度检验
    accuracy of model prediction
    渤海Bohai Sea $\stackrel{\Lambda}{Y}^{(1)}(t)=21.51 \mathrm{e}^{0.084(t-1961)}-17.21$ 0.39 0.90 合格
    黄海Yellow Sea $\stackrel{\Lambda}{Y}^{(1)}(t)=112.84 \mathrm{e}^{0.052(t-1961)}-100.17$ 0.50 0.90 勉强
    东海East China Sea GM(1, 1)预测模型及其n次残差序列预测模型的预测精度检验均不合格
    远洋渔业Pelagic fishing $\hat{q}^{(6)}(t)=-14.00 \mathrm{e}^{-0.27(t-1996)}+15.84$ 0.30 1.00
    注:GM(1, 1)的预测精度检验标准:当c<0.35,P>0.95时,好;当c<0.45,P>0.8时,合格;当c<0.5,P>0.7时,勉强;当c≥0.65,P≤0.7时,不合格Note: The test criterion of GM(1, 1): c<0.35,P>0.95, the result is excellent; c<0.45,P>0.8, the result is qualified; c<0.5,P>0.7, the result is unconvinced; c≥0.65,P≤0.7, the result is disqualified.
    下载: 导出CSV

    表  5   辽宁省2004~2010年海洋渔业预测产量

    Table  5   The calculated yield of marine fishery in Liaoning from 2004 to 2010 104t

    年份
    year
    海水养殖
    mariculture
    远洋渔业
    pelagic fishing
    渤海
    Bohai Sea
    黄海
    Yellow Sea
    东海
    East China Sea
    2004 220.82 20.54 56.16 67.79 2.95
    2005 246.67 21.90 58.81 71.26 3.02
    2006 277.57 23.60 62.16 74.41 3.00
    2007 310.12 25.66 65.88 77.50 3.01
    2008 348.92 28.12 69.88 80.64 3.01
    2009 389.89 31.01 74.15 83.87 3.02
    2010 438.60 34.39 78.68 87.21 3.02
    下载: 导出CSV

    表  6   辽宁省海洋渔业预测产量的增量

    Table  6   The increment of marine fishery predicted yield in Liaoning province 104t

    年份
    year
    与2003年产量相比compared with the yield in 2003 与MSY相比国内海洋捕捞产量增量
    the increment of domestic fishing predicted yield compared with MSY
    海水养殖
    mariculture
    远洋渔业
    pelagic fishing
    国内海洋捕捞业
    domestic fishing area
    2004 24.62 -5.26 5.04 24.27
    2005 50.47 -3.90 11.23 30.46
    2006 81.37 -2.20 17.71 36.94
    2007 113.92 -0.14 24.53 43.76
    2008 152.72 2.32 31.67 50.90
    2009 193.69 5.21 39.18 58.41
    2010 242.40 8.59 47.05 66.28
    下载: 导出CSV
  • [1] 陈新军, 周应祺. 中国渔业产量结构的灰色关联分析[J]. 中国渔业经济, 2001, (6): 30-33. doi: 10.3969/j.issn.1009-590X.2002.02.016
    [2] 史岩, 张小栓, 穆维松, 等. 渔业产量增长与结构变化研究[J]. 莱阳农学院学报, 2004, 21(1): 71-73. doi: 10.3969/j.issn.1674-148X.2004.01.020
    [3] 中华人民共和国农业部水产司. 中国渔业统计四十年[M]. 北京: 海洋出版社, 1991. https://xueshu.baidu.com/usercenter/paper/show?paperid=e52f3ea9f175ff3fb3eea98cfdee88d9&site=xueshu_se&hitarticle=1
    [4] 中华人民共和国农业部渔业局. 中国渔业统计汇编(1989-1993)[M]. 北京: 海洋出版社, 1996. https://xueshu.baidu.com/usercenter/paper/show?paperid=f3221ab34bbf0d5dec6c8076a1743e0f&site=xueshu_se
    [5] 农业部渔业局. 中国渔业统计年鉴2000[M]. 北京: 海洋出版社, 2001. https://www.shujuku.org/china-fishery-statistical-yearbook.html
    [6] 农业部渔业局. 中国渔业统计年鉴2001[M]. 北京: 海洋出版社, 2002. https://xueshu.baidu.com/usercenter/paper/show?paperid=889db05498e664468332c6c60ee4bcc5&site=xueshu_se
    [7] 农业部渔业局. 中国渔业统计年鉴2002[M]. 北京: 海洋出版社, 2003. https://www.shujuku.org/china-fishery-statistical-yearbook.html
    [8] 农业部渔业局. 中国渔业统计年鉴2003[M]. 北京: 海洋出版社, 2004. https://www.shujuku.org/china-fishery-statistical-yearbook.html
    [9] 詹秉义. 渔业资源评估[M]. 北京: 中国农业出版社, 1993.173-179. https://xueshu.baidu.com/usercenter/paper/show?paperid=7d780225e511d43a52f368e49751c36b&site=xueshu_se&hitarticle=1
    [10] 邓聚龙. 灰色预测与决策[M]. 武汉: 华中理工大学出版社, 1992.97-191. https://xueshu.baidu.com/usercenter/paper/show?paperid=af984f29bbfc96db20f6f034d65c77f6&site=xueshu_se&hitarticle=1
    [11] 田诤(译). 时间序列的理论与方法[M]. 北京: 高等教育出版社, 2001.1-257. https://xueshu.baidu.com/usercenter/paper/show?paperid=e6df18ac8819c3cc2cf7781734658658&site=xueshu_se
    [12] 中国海洋年鉴编纂委员会. 中国海洋年鉴2001[M]. 北京: 海洋出版社, 2001.300-310. https://xueshu.baidu.com/usercenter/paper/show?paperid=d49fa9365247d40756e645496835eb01&site=xueshu_se
    [13] 大连市史志办公室. 大连年鉴2000[M]. 上海: 浦东电子出版社, 2000.101-104.
  • 期刊类型引用(5)

    1. 付鑫洋,司烨仪,周鑫,刘莉莉. 我国拖网网囊网目选择性研究现状. 农村经济与科技. 2023(02): 73-76 . 百度学术
    2. 晏磊,杨炳忠,李亚男. 南海北部拖网方形和菱形混合网目网囊的渔获性能研究. 水产科学. 2023(04): 690-697 . 百度学术
    3. 杨炳忠,晏磊,李杰,王腾,张鹏. 南海区双船大网目拖网网囊网目对带鱼的选择性. 中国水产科学. 2023(10): 1271-1281 . 百度学术
    4. 李亚男,杨炳忠,张鹏,李杰,王腾,晏磊. 南海北部拖网对蓝圆鲹的选择性研究. 南方水产科学. 2022(03): 170-176 . 本站查看
    5. 许庆昌,李显森,孙珊,樊钢洲. 秋季莱州湾弓子网网囊网目尺寸选择性研究. 渔业科学进展. 2020(06): 19-27 . 百度学术

    其他类型引用(0)

图(5)  /  表(6)
计量
  • 文章访问数:  5193
  • HTML全文浏览量:  143
  • PDF下载量:  3265
  • 被引次数: 5
出版历程
  • 收稿日期:  2005-08-01
  • 修回日期:  2005-08-19
  • 刊出日期:  2005-12-19

目录

/

返回文章
返回