本发明涉及海洋渔业和海洋科学技术领域，具体而言，本发明涉及一种基于栖息地指数模型的秋刀鱼中心渔场预测方法。
背景技术：
渔情预报是基于鱼类与海洋环境要素间关系及规律，对未来一定时期和水域内水产资源状况进行预报，主要包括渔期、中心渔场、资源丰度及可能达到的渔获量等，是渔业生产的重要环节。
秋刀鱼属冷水性高度洄游中上层鱼类，广泛分布于西北太平洋亚热带和温带的亲潮、黑潮水域。秋刀鱼一年中春、秋、冬季均可产卵，生命周期短，资源恢复力强，年间资源量波动显著，有学者认为该现象并非是由过度捕捞造成，而与生境条件影响下的亲体繁殖力、卵及幼鱼成活率有直接联系。目前，秋刀鱼捕捞国家和地区主要有中国、日本、韩国、俄罗斯、中国台湾省以及瓦努阿图，年总渔获量在35-60万吨。鉴于秋刀鱼的经济、文化、生态等价值，现已被北太平洋渔业委员会列为优先管理鱼种之一。现有研究表明秋刀鱼资源丰度和渔场分布易受栖息地环境条件影响，随着气候变化的逐渐加剧，海洋-气候的综合作用对该鱼种渔场预报提出了更高要求。
渔情及渔获量变动往往与海洋环境密切相关，海表面温度(seasurfacetemperature，简称sst)、海表面水平梯度(seasurfacetemperaturegradient，简称sstg)、海表面高度(seasurfaceheight，简称ssh)、叶绿素浓度(seasurfacechlorophyllaconcentration，简称chl-a)等都会影响鱼类分布。
在诸多环境要素中，水温变化一致被认为对鱼类产卵、仔鱼发育和成活率等有重要意义，同时还影响鱼类饵料生物的消盛，这些都会直接或间接地影响到鱼类资源量的分布、洄游移动和空间集群等。一般情况下，只有在适宜的温度条件下才能形成优良渔场，因此依据鱼类的特征温度推测渔场位置准确性更高。影响秋刀鱼渔场形成的主要因子是水温，主要反映在渔场和渔期变化上。
秋刀鱼渔场主要分布在黑潮、亲潮及两者交汇海域，来源于北赤道流的高温、高盐黑潮水系与起源于白令海的低温、低盐亲潮水系交汇于北海道东部海域，其混合区域构成了亚极海洋锋面(40°n附近)，宽度约2-4个纬距，该锋面在160°e以西海域较为明显，以东海域不明显。
亚极海洋锋面南部锋(一般在36°-37°n)和北部锋(一般在42°-43°n)间的区域形成混合区，锋面北侧亚极环流的持续性气旋分散作用致使冷水上扬，使得底部营养盐上升，加上锋面附近不同水体携运营养盐，致使浮游生物大量繁殖，对秋刀鱼渔场形成提供了基础保障。此外，秋刀鱼渔场与暖涡分布有着密切关系，其变动会影响秋刀鱼洄游路线及渔场形成，而海面高度和涡旋流关系密切，所以可以通过海面高度探究秋刀鱼渔场形成机制。
国内外学者对秋刀鱼渔情预报进行过相关研究。ito在2013年通过基于生态系统的生物能模型研究了气候变暖对秋刀鱼资源影响情况，认为水温升高会延迟秋刀鱼向南迁徙时期，并因浮游生物含量减小而影响幼鱼生长和资源密度状况。tseng则根据2006-2010年卫星遥感影像分析了西北太平洋秋刀鱼栖息地环境因子，表明适合秋刀鱼栖息的sst为14°～16℃、chl-a为0.4～0.6mg/m3、初级生产力浓度为600～800mg/cm3。在探究秋刀鱼资源空间变动时，广义线性模型解释率之35.7％较垂直广义生产模型解释率之20.5％更高。
chen在2011年结合2006-2008年台湾秋刀鱼生产数据和卫星遥感温度数据，研究北太平洋秋刀鱼最适温度范围，并预测秋刀鱼潜在的栖息地分布。结果表明，秋刀鱼偏好温度范围为12°～18.5℃，并随月份不断变化；此外，秋刀鱼栖息地范围随水温增加会向北偏移，栖息地最南端可由40.24°n移动到46.15°n。tian等在2004通过分析秋刀鱼与水温、厄尔尼诺-南方涛动等海洋气候因子间关系，认为大尺度的气候和海洋环境变动可强烈影响秋刀鱼的年间资源量和捕捞产量。张孝民根据2013年7～12月秋刀鱼生产调查数据和海洋环境数据，利用海洋渔业地理信息系统和数理统计探究秋刀鱼渔场在时间和空间上的分布情况。韩士鑫报道了通过卫星红外观测黑潮暖流环和其周围冷水变化对秋刀鱼渔场形成的短期预报情况。tseng等认为sst前锋是影响秋刀鱼分布和提高渔场预报准确度的基础。tian等指出，秋刀鱼的幼体生长和存活率很大程度上受海洋条件的影响。
渔情预报是渔场学研究的重点，准确的预报模型可有效缩短寻找渔场时间，降低捕捞成本、提高渔获产量，并可指导秋刀鱼生产企业进行合理安排布局。目前，相关国家和地区逐渐加大秋刀鱼资源开发力度，更需要选择合理的预测模型，以有效提高渔船生产效率，为政府和企业秋刀鱼年度生产规划提供参考依据。
现有常见预测模型主要包括：范例推理、贝叶斯方法、决策树模型、神经网络模型、线性回归模型、等权重的栖息地指数模型、最大熵模型等。由于海洋环境的复杂性、区域性和易变性，致使渔场分布会有较大的时空波动，所以准确地对中心渔场进行预测存在诸多困难，但却至关重要。现有预测模型在处理某些特殊问题时会有一定的局限性，致使目前仍缺乏统一的渔情速报预测方法，以上各常见模型的不足之处主要包括：
范例推理：应用该模型前，需要事先建立大量渔场范例，以供预报推理时选择类似的范例，对范例数量的要求较高，必须有足够多的样本；此外，该模型所提取的范例主要是对历史数据的总结，很难应用于预测新渔场；
贝叶斯方法：需要将环境数据进行人为分级处理，数据尺度大或结构复杂时工作量大，且只限用于处理单个环境因素；某些情况下无法适用，如捕捞数据统计精度有限，会因与先验概率周期不一致使得对渔场预报精度大大降低；
决策树模型：当数据量过大时可能会出现过度拟合现象，且顶层的错误会向下传播，从而影响下面所有层的分裂情况；
神经网络模型：对知识的表达、存储和推理计算都是隐式的，无法解释本身的推理依据和推理过程，存在收敛速度慢、局部极值的缺点；
线性回归模型：渔场与海洋环境因子间关系一般具有随机性和模糊性，有时模型无法满足对数据的假设，或很难建立相关性很高的回归方程，现实条件下环境因子与鱼类分布大多是呈非线性关系，某些线性拟合模型可能无法满足；
等权重的栖息地指数模型：多种环境因子对鱼类分布的影响效果并不是完全一致的，不同的环境因子应赋予不同的权重，但以往栖息地指数模型中环境因子的权重主要由专家知识赋值或默认为相等权重，无法反映出各环境因子对栖息地质量的贡献率差异；
最大熵模型：根据鱼类当前的分布情况去预测其未来栖息地，更适合洄游能力弱且群体结构单一的物种，对某些洄游能力强或种群结构复杂群体的应用性减小。
技术实现要素：
为了寻找更为有效的渔情速报预测的实现方案，本发明实施例提供了一种基于栖息地指数模型的秋刀鱼中心渔场预测方法，该基于栖息地指数模型的秋刀鱼中心渔场预测方法包括如下步骤：基于栖息地指数模型的秋刀鱼中心渔场预测方法，其特征在于，所述基于栖息地指数模型的秋刀鱼中心渔场预测方法包括如下步骤：
获取历史的n个月份秋刀鱼的生产数据和与所述生产数据对应的环境因子，所述生产数据包括作业日期、经纬度信息、渔获量、捕捞努力量；所述环境因子包括sst、sstg和ssh；所述n为大于0的自然数；
利用所述环境因子并基于产量密度分别建立各个月份的适应性指数与sst、sstg、ssh间关系，其中，产量密度方程如下：
其中，y表示sst、sstg、ssh在不同月份m和预置组距c(k)下的适应性指数
x为sst、sstg、ssh在各自组距c(k)内的平均值
a、b、c为产量密度方程中的参数，m为n个月份中的各个月份，k为1、2、3；
用非线性最小二乘法拟合所述环境因子与所述适应性指数的关系并估计产量密度方程中参数a、b、c；
基于brt计算的sst、sstg、ssh相对所述适应性指数的贡献率确定sst、sstg、ssh在栖息地指数中的权重值wsst，m、wsstg，m、wssh，m，并得到栖息地指数如下：
优选地，所述n为6，所述m为6、7、…、11；所述获取历史的n个月份秋刀鱼的生产数据和与所述生产数据对应的环境因子为获取2013—2015年6-11月份秋刀鱼的生产数据和与所述生产数据对应的环境因子。
优选地，所述利用所述环境因子并基于产量密度分别建立各个月份的适应性指数与sst、sstg、ssh间关系包括如下步骤：
从预设的sst、sstg、ssh的组距范围中基于预置步长按照拟合度r2确定sst、sstg、ssh的最优组距；
设定捕捞努力量最高渔区的适应性指数为1，计算其它渔区的适应性指数，计算方式如下：
其中，si为适应性指数，effortn为某渔区n的捕捞努力量；effortmax为所有渔区中最高的捕捞努力量；
利用产量密度方程分别建立各月份的适应性指数与sst、sstg、ssh间关系，以得到sst、sstg、ssh的最适值。
优选地于，所述sst、sstg、ssh的最优组距分别为：sst的最优组距为0.3℃，sstg的最优组距为0.005℃/km，ssh的最优组距为0.02m。
优选地，所述基于brt计算的sst、sstg、ssh相对所述适应性指数的贡献率确定sst、sstg、ssh在栖息地指数中的权重值wsst，m、wsstg，m、wssh，m为在r编程环境中使用gbm函数确定。
优选地，所述权重值wsstg，m为40％，所述权重值wsst，m和所述权重值wssh，m均为30％。
优选地，所述得到栖息地指数之后包括如下步骤：
基于栖息地指数和实时获取的环境因子预测秋刀鱼栖息地分布。
与现有技术相比，本发明实施例一种基于栖息地指数模型的秋刀鱼中心渔场预测方法具有如下有益效果：
本发明实施例一种基于栖息地指数模型的秋刀鱼中心渔场预测方法采用产量密度方程建立适应性指数，结合brt方法合理选取各月份的环境因子权重，分别计算了不同作业月份秋刀鱼适宜栖息地各环境因子范围，预测了秋刀鱼的栖息地指数模型分布，为生产作业提供技术支持，该方法可提高秋刀鱼中心渔场预测的准确性，可行性较高。
本发明附加的方面和优点将在下面的描述中部分给出，这些将从下面的描述中变得明显，或通过本发明的实践了解到。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能解释为对本发明的限制。
本发明实施例一种基于栖息地指数模型的秋刀鱼中心渔场预测方法，其包括如下步骤s101～s107：
步骤s101：获取历史的n个月份秋刀鱼的生产数据和与所述生产数据对应的环境因子，所述生产数据包括作业日期、经纬度信息、渔获量、捕捞努力量；所述环境因子包括sst、sstg和ssh；所述n为大于0的自然数，其中，n优选为6。
示例地，获取历史的n个月份秋刀鱼的生产数据和与生产数据对应的环境因子为获取2013—2015年6-11月份秋刀鱼的生产数据和与生产数据对应的环境因子。
在一些实施方式中，sst、ssh直接由卫星遥感下载获取，sst空间分辨率为0.1°x0.1°，ssh空间分辨率为0.25°x0.25°，两者时间分辨率均为天。sstg则通过sst、ssh计算得到，具体如下：
假设某位置处的sst值为ssti，j，其相邻四个网格位置处sst值依次为ssti+1，j，ssti-1，j，ssti，j+1，ssti，j-1，则sstg计算方式如下：
其中，i、j为正整数，分别表示网格数据行号和列号；△x表示第j–1列与第j+1列间在经度方向距离，△y表示第i–1行与第i+1行间在纬度方向距离，单位均为km；sstgi，j表示当前位置的sstg值，单位为℃/km。
步骤s103：利用所述环境因子并基于产量密度分别建立各个月份的适应性指数与sst、sstg、ssh间关系，其中，产量密度方程如下：
其中，y表示sst、sstg、ssh在不同月份m和预置组距c(k)下的适应性指数
值得注意的是，当获取的是2013—2015年6-11月份秋刀鱼的生产数据和与生产数据对应的环境因子时，m取值为6、7、…、11。
在一些实施方式中，定义cpue为每艘船每天捕捞产量，设定0.5°×0.5°为一个渔区。根据渔船作业日期和经纬度信息，查找该日期内环境数据网格上最近位置经纬度所对应的数据值，对渔业数据和环境数据进行匹配，按月份分别计算各渔区的捕捞努力量及cpue、sst、sstg及ssh四要素的平均值。
具体地，利用环境因子并基于产量密度分别建立各个月份的适应性指数与sst、sstg、ssh间关系包括如下步骤：
从预设的sst、sstg、ssh的组距范围中基于预置步长按照拟合度r2确定sst、sstg、ssh的最优组距；
设定捕捞努力量最高渔区的适应性指数为1，计算其它渔区的适应性指数，计算方式如下：
其中，si为适应性指数，effortn为某渔区n的捕捞努力量；effortmax为所有渔区中最高的捕捞努力量；
利用产量密度方程分别建立各月份的适应性指数与sst、sstg、ssh间关系，以得到sst、sstg、ssh的最适值。
值得注意的时，sst、sstg、ssh的最适值通过对产量密度方程求导得到：
当y'＝0时，x＝-b/2c，根据函数性质可以计算每个环境因子的最适值。
在一些实施方式中，预设的sst、sstg、ssh的组距范围以及预置步长如下：sst的组距范围为0.1℃-0.3℃，步长为0.1℃；sstg的组距范围为0.002℃/km-0.005℃/km，步长为0.001℃/km；ssh的组距范围为0.02m-0.04m，步长为0.01m。
优选地，sst、sstg、ssh的最优组距分别为：sst的最优组距为0.3℃，sstg的最优组距为0.005℃/km，ssh的最优组距为0.02m。
步骤s105：用非线性最小二乘法拟合所述环境因子与所述适应性指数的关系并估计产量密度方程中参数a、b、c。
步骤s107：基于brt计算的sst、sstg、ssh相对所述适应性指数的贡献率确定sst、sstg、ssh在栖息地指数中的权重值wsst，m、wsstg，m、wssh，m，并得到栖息地指数如下：
增强回归树(boostingregressiontrees，brt)通过随机选择与自学习的方式产生多重回归树，提高了模型稳定性和预测精度，可用于处理非线性、相关和相互作用的变量关系，能解释生物与环境间的复杂关系，且不受多重共线性、缺失预测变量和异常值的影响，每个环境因子对鱼类分布的相对重要性可通过其方差贡献率来确定，以此进行合理权重赋值。
在一些实施方式中，基于brt计算的sst、sstg、ssh相对所述适应性指数的贡献率确定sst、sstg、ssh在栖息地指数中的权重值wsst，m、wsstg，m、wssh，m为在r编程环境中使用gbm函数确定。
在一些实施方式中，权重值wsst，m所占百分比范围为31.99％-34.87％，权重值wsstg，m占比范围为34.71％-43.79％，权重值wssh，m范围为24.22％-31.37％。
优选地，权重值wsstg，m为40％，权重值wsst，m和权重值wssh，m均为30％。
在一些实施方式中，得到栖息地指数之后包括如下步骤：
基于栖息地指数和实时获取的环境因子预测秋刀鱼栖息地分布。
与现有技术相比，本发明实施例一种基于栖息地指数模型的秋刀鱼中心渔场预测方法具有如下有益效果：
本发明实施例一种基于栖息地指数模型的秋刀鱼中心渔场预测方法采用产量密度方程建立适应性指数，结合brt方法合理选取各月份的环境因子权重，分别计算了不同作业月份秋刀鱼适宜栖息地各环境因子范围，预测了秋刀鱼的栖息地指数模型分布，为生产作业提供技术支持，该方法可提高秋刀鱼中心渔场预测的准确性，可行性较高。
以上所述仅是本发明的部分实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。