[发明专利]基于太平洋震荡指数的日本冬生群鱿鱼资源丰度预测方法

权利要求书

1.一种基于太平洋震荡指数的日本冬生群鱿鱼资源丰度预测方法，其特征在于，其包括以下步骤：

S1、获取前N年日本冬生群鱿鱼分布的北海道的太平洋一侧海域的每月的太平洋震荡指数PDO值；

S2、利用时间序列分析方法，对日本冬生群鱿鱼资源丰度CPUE与前N年每月的PDO值进行相关性分析，获取在统计上相关的P＜0.05的月PDO值，该些月PDO值作为影响日本冬生群鱿鱼资源丰度的气候因子；

S3、利用多元线性方程建立多个日本冬生群鱿鱼资源丰度预测模型并计算统计上P值，其公式为：

CPUE＝a+b1*x1+b2*x2+b3*x3+……+bn*xn

式中，CPUE为日本中小型鱿钓船日产量，作为日本冬生群鱿鱼资源丰度衡量指标，a为常数，b1、b2、b3、……、bn为方程的系数；x1、x2、x3、……、xn为影响资源丰度的月PDO值；

S4、在上述多个日本冬生群资源丰度预测模型中，选择统计上P值最小的模型作为最优模型。

2.如权利要求1所述的日本冬生群鱿鱼资源丰度预测方法，其特征在于，在步骤S2中，根据日本冬生群鱿鱼资源丰度CPUE与前2年各月的PDO值的相关性分析，日本冬生群鱿鱼资源丰度CPUE与前2年的10-12月PDO值相关性显著，且呈现负相关，其相关系数分别为-0.4506(P＜0.05)、-0.4985(P＜0.05)、-0.5878(P＜0.01)；

根据日本冬生群鱿鱼资源丰度CPUE与前1年各月PDO值的相关性分析，日本冬生群鱿鱼资源丰度CPUE与前1年的1-4月PDO值相关性显著，且呈现负相关，其相关系数分别为-0.4665(P＜0.05)、-0.4365(P＜0.05)、-0.4295(P＜0.05)、-0.5072(P＜0.01)；

根据日本冬生群鱿鱼资源丰度CPUE与同年的各月PDO值的相关性分析认为，日本冬生群鱿鱼资源丰度CPUE与同年的1-4月PDO值相关性显著，且呈现负相关，其相关系数分别为-0.4746(P＜0.05)、-0.4837(P＜0.05)、-0.5458(P＜0.01)、-0.5570(P＜0.01)。

3.如权利要求2所述的日本冬生群鱿鱼资源丰度预测方法，其特征在于，在步骤S3中，

1)预测模型之一

以前2年的11-12月PDO值(PDO_t-2，12、PDO_t-2，11)、同年的2-3月PDO值(PDO_t，2，PDO_t，3)和上一年度3月份PDO值(PDO_t-1，3)作为预测因子，建立日本冬生群鱿鱼资源丰度预测模型为：

CPUE＝1.2048+0.0330*PDO_t-2，12-0.1811*PDO_t-2，11-0.2260*PDO_t，₂+0.0749*PDO_t，3-0.0196*PDO_t-1，3

其F值为4.5183，P＝0.0069＜0.01；

2)预测模型之二

以前2年的11-12月PDO值(PDO_t-2，12、PDO_t-2，11)、同年的2-3月PDO值(PDO_t，2，PDO_t，3)作为预测因子，建立日本冬生群鱿鱼资源丰度预测模型为：

CPUE＝1.1968+0.0273*PDO_t-2，12-0.1865*PDO_t-2，11-0.2290*PDO_t，2+0.0740*PDO_t，3

其F值为5.9135，P＝0.0026＜0.01；

3)预测模型之三

以同年的2-3月PDO值(PDO_t，2，PDO_t，3)和前一度3月份PDO值(PDO_t-1，3)作为预测因子，建立日本冬生群鱿鱼资源丰度预测模型为：

CPUE＝1.3257-0.1120*PDO_t，2-0.0461*PDO_t，3-0.1307*PDO_t-1，3

其F值为5.1699，P＝0.0078＜0.01；

4)预测模型之四

以同一年的2-3月PDO值(PDO_t，2，PDO_t，3)和上一年度3月份PDO值(PDO_t-1，3)作为预测因子，建立日本冬生群鱿鱼资源丰度预测模型为：

CPUE＝1.3093-0.0792*PDO_t，2-0.1223*PDO_t，3

其F值为5.1233，P＝0.0149＜0.05；

5)预测模型之五

以前2年的11-12月PDO值(PDO_t-2，12、PDO_t-2，11)作为预测因子，建立日本海冬生群鱿鱼资源丰度预测模型为：

CPUE＝1.1582+0.0383*PDO_t-2，12-0.2125*PDO_t-2，11

其F值为5.8894，P＝0.0089＜0.01。

4.如权利要求3所述的日本冬生群鱿鱼资源丰度预测方法，其特征在于，在步骤S4中，选择前2年的10-12月PDO值、前1年1-4月PDO值、同年1-4月PDO值作为气候预测因子，日本海冬生群鱿鱼资源丰度预测模型为：CPUE＝1.1968+0.0273*PDO_t-2，12-0.1865*PDO_t-2，11-0.2290*PDO_t，₂+0.0740*PDO_t，3，PDO_t-2，12、PDO_t-2，11分别是前2年的12、11月PDO值，PDO_t，2、PDO_t，3分别是同年的2-3月PDO值。