本发明涉及头足类渔场预报技术领域，特别是涉及一种基于不同权重的海洋环境因子的摩洛哥头足类渔场预报方法。
背景技术：
摩洛哥临北大西洋与地中海，是地中海与大西洋之间相互沟通联系的桥梁。海洋渔业是摩洛哥主要外汇来源，在其国内经济发展中处于关键位置。在摩洛哥的渔业资源中，头足类具有最高经济效益的。在摩洛哥，距离我国渔船第一次开展底拖网作业，已经过去20多年了，但底拖网主要捕捞对象一直以来还都是章鱼、鱿鱼这些头足类。国内的众多专家们曾对摩洛哥的渔业发展、渔具渔法以及水域的中上层鱼类资源等各方面开展了许多的调查与讨论。但是，国内与摩洛哥的头足类渔场分布有关的文献，这方面的研究仍是不多。
技术实现要素：
本发明针对现有技术存在的问题和不足，提供一种基于不同权重的海洋环境因子的摩洛哥头足类渔场预报方法。
本发明是通过下述技术方案来解决上述技术问题的：
本发明提供一种基于不同权重海洋环境因子的摩洛哥头足类渔场预报方法，其特点在于，其包括以下步骤：
s1、获取2012-2015年摩洛哥渔场的渔获生产统计数据，该渔获生产统计数据包括作业时间、经度、纬度、作业水深、作业网次以及作业产量；
s2、获取该渔获生产统计数据对应的海洋环境数据，该海洋环境数据包括海表温度sst、海面高度距平均值ssha和叶绿素浓度chl-a，该海洋环境数据以月为时间分辨率、0.5°×0.5°为空间分辨率，统计时间为2012-2015年的1-3月和11-12月；
s3、以作业网次、作业产量所占比率以及每个区间内的平均网次产量作为中心渔场的指标，研究其与sst、ssha、chl-a以及水深之间的关系；
s4、建立不同环境因子的适应性指数si，采用专家赋值法赋予适应性指数si的值，设定作业网次最高的海域为中心渔场的分布概率最高的海域，则适应性指数si赋值为1；当没有作业网次则适应性指数si赋值为0；当作业网次高于均值的海域，则适应性指数si赋值为0.5；当作业网次低于均值的海域，则适应性指数si赋值为0.1；
s5、利用公式hsi＝xsst*isi_sst+xssha*isi_ssha+xchl-a*isi-chl-a+xdepth*isi_depth计算5种不同权重方案下的栖息地综合指数his；
式中：isi_sst表示基于海表温度的适应性指数；isi_ssha表示基于海面高度距平均值的适应性指数；isi-chl-a表示基于叶绿素浓度的适应性指数；isi_depth表示基于水深的适应性指数；xsst、xssha、xchl-a、xdepth分别表示海表温度、海面高度距平均值、叶绿素浓度、水深所占权重的值；
其中，权重方案1：xsst为0.25，xssha为0.25，xchl-a为0.25，xdepth为0.25；
权重方案2：xsst为0，xssha为0.9，xchl-a为0，xdepth为0.1；
权重方案3：xsst为0.1，xssha为0.1，xchl-a为0，xdepth为0.8；
权重方案4：xsst为0.9，xssha为0.1，xchl-a为0，xdepth为0；
权重方案5：xsst为0.4，xssha为0.4，xchl-a为0.1，xdepth为0.1；
分别分析5种不同权重方案中hsi值与设定阈值的大小，从而比较得出最佳权重方案，来用于预测摩洛哥头足类渔场的中心渔场。
较佳地，步骤s3包括：
以sst相隔1℃为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场sst的范围；
以ssha相隔10cm为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场ssha的范围；
以chl-a的含量在0.01～1.0、1.0～2.0、2.0～5.0、5.0～20.0、20.0～50.0mg/m3为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场chl-a的范围；
以水深相隔10m为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场水深的范围。
在符合本领域常识的基础上，上述各优选条件，可任意组合，即得本发明各较佳实例。
本发明的积极进步效果在于：
本发明通过不同权重对摩洛哥头足类栖息地模型的影响，获得影响头足类栖息地分布的主要环境因子，为后续中心渔场预测提供基础。
附图说明
图1为本发明较佳实施例的摩洛哥头足类渔场预报方法的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本实施例提供一种基于不同权重海洋环境因子的摩洛哥头足类渔场预报方法，其包括以下步骤：
步骤101、获取2012-2015年摩洛哥渔场的渔获生产统计数据，该渔获生产统计数据包括作业时间、经度、纬度、作业水深、作业网次以及作业产量。
摩洛哥渔场的渔获生产统计数据来源于上海蒂尔远洋渔业公司。
步骤102、获取该渔获生产统计数据对应的海洋环境数据，该海洋环境数据包括海表温度sst、海面高度距平均值ssha和叶绿素浓度chl-a，该海洋环境数据以月为时间分辨率、0.5°×0.5°为空间分辨率，统计时间为2012-2015年的1-3月和11-12月。
步骤103、以作业网次、作业产量所占比率以及每个区间内的平均网次产量作为中心渔场的指标，研究其与sst、ssha、chl-a以及水深之间的关系。
1)以sst相隔1℃为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场sst的范围；
2)以ssha相隔10cm为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场ssha的范围；
3)以chl-a的含量在0.01～1.0、1.0～2.0、2.0～5.0、5.0～20.0、20.0～50.0mg/m3为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场chl-a的范围；
4)以水深相隔10m为间距，计算出每个间距内的作业网次、作业产量比率以及平均网次产量，并依此得出最适宜的中心渔场水深的范围。
步骤104、建立不同环境因子的适应性指数si，采用专家赋值法赋予适应性指数si的值，设定作业网次最高的海域为中心渔场的分布概率最高的海域，则适应性指数si赋值为1；当没有作业网次则适应性指数si赋值为0；当作业网次高于均值的海域，则适应性指数si赋值为0.5；当作业网次低于均值的海域，则适应性指数si赋值为0.1。
一般渔船在进行作业对渔场进行判别时，都是根据船长的经验以及探鱼仪的影像，所以作业网次可被看作是发现鱼类的指标，用来表示为栖息地的适应性指数。
表1栖息地适应性指数的确定标准
步骤105、利用公式hsi＝xsst*isi_sst+xssha*isi_ssha+xchl-a*isi-chl-a+xdepth*isi_depth计算5种不同权重方案下的栖息地综合指数his。
计算相关海洋环境因子不同权重下，栖息地综合指数(habitatsuitabilityindex，hsi)，hsi在0到1之间变化。一般认为hsi大于0.6的区域为中心渔场分布的海域。
式中：isi_sst、isi_ssha、isi-chl-a、isi_depth分别是各自基于表温、海面高度距平均值、叶绿素浓度、水深的适应性指数。xsst、xssha、xchl-a、xdepth分别是表温、海面高度距平均值、叶绿素浓度、水深所占权重的值。总共设有不同权重的5种方案，各方案如下表2。
表2基于中心渔场相关的不同环境因子设定权重值
分别分析5种不同权重方案中hsi值与设定阈值的大小，从而比较得出最佳权重方案，来用于预测摩洛哥头足类渔场的中心渔场。
利用2012年-2015年1-3月、11-12月统计的数据对不同权重方案进行比较，hsi值可分为0-0.2、0.2-0.4、0.4-0.6、0.6-0.8、0.8-1.0，在此基础上分别进行统计分析5种不同权重方案中hsi值>0.6以及hsi值<0.4的大小关系，从而比较得出其中最佳方案，来用于预测中心渔场。
根据上述方法，下面基于具体的统计数据进行分析：
1、生产情况分析
1)渔场分布与表温sst的关系
分析结果表明，头足类渔场分布与表温关系密切，不同月份有不同的适宜sst范围。1月份作业主要分布在sst范围为16～19℃海域，高平均网次产量的适宜sst范围为16～17℃和18～19℃，平均产量为130～153kg。2月份作业主要分布在sst范围为15～19℃海域，高平均网次产量的适宜sst范围为16～17℃和18～19℃，平均产量为122～147kg。3月份作业主要分布在sst范围为15～17℃海域，高平均网次产量的适宜sst范围为15～16℃，平均产量为89.16kg。11月份作业主要分布在sst范围为18～23℃海域，高平均网次产量的适宜sst范围为19～23℃，平均产量为162～185kg。12月份作业主要分布在sst范围为16～21℃海域，高平均网次产量的适宜sst范围为20～21℃，平均产量为457kg。
2)渔场分布与ssha的关系
分析结果表明，头足类渔场分布与海面高度距平值关系密切，不同月份有不同的适宜ssha范围。1月份作业主要分布在ssha范围为-60～-20cm海域，高平均网次产量的适宜ssha范围为-60～-30cm，平均产量为124～143kg。2月份作业主要分布在ssha范围为-60～-30cm海域，高平均网次产量的适宜ssha范围为-60～-40cm，平均产量为123～137kg。3月份作业主要分布在ssha范围为-60～-30cm海域，高平均网次产量的适宜ssha范围为-60～-40cm，平均产量为96～101kg。11月份作业主要分布在ssha范围为-50～10cm海域，高平均网次产量的适宜ssha范围为-40～-30和-10～0cm，平均产量为189～209kg。12月份作业主要分布在ssha范围为-50～10cm海域，高平均网次产量的适宜ssha范围为-50～-40cm，平均产量为558.69kg。
3)渔场分布与叶绿素浓度的关系
分析结果表明，头足类渔场分布与叶绿素浓度关系密切，不同月份有不同的适宜水深范围。1月份作业主要分布在chl-a范围为0.01～50mg/m3海域，高平均网次产量的适宜chl-a范围为1.0～5.0mg/m3，平均产量为96～127kg。2月份作业主要分布在chl-a范围为0.01～20mg/m3海域，高平均网次产量的适宜chl-a范围为1.0～20.0mg/m3，平均产量为119～128kg。3月份作业主要分布在chl-a范围为0.01～50mg/m3海域，高平均网次产量的适宜chl-a范围为1.0～2.0和5.0～50mg/m3，平均产量为99～110kg。11月份作业主要分布在chl-a范围为0.01～20mg/m3海域，高平均网次产量的适宜chl-a范围为0.01～5.0mg/m3，平均产量为169～176kg。12月份作业主要分布在chl-a范围为0.01～50mg/m3海域，高平均网次产量的适宜chl-a范围为2.0～5.0mg/m3，平均产量为256.24kg。
4)渔场分布与水深的关系
分析结果表明，头足类渔场分布与作业水深关系密切，不同月份有不同的适宜水深范围。1月份作业主要分布在水深范围为20～90m海域，高平均网次产量的适宜水深范围为20～40m，平均产量为131～140kg。1月份作业主要分布在水深范围为20～100m海域，高平均网次产量的适宜水深范围为20～50m和60～70m，平均产量为117～141kg。3月份作业主要分布在水深范围为20～80m海域，高平均网次产量的适宜水深范围为70～80m，平均产量为169kg。11月份作业主要分布在水深范围为30～80m海域，高平均网次产量的适宜水深范围为30～40m，平均产量为246.49kg。12月份作业主要分布在水深范围为20～80m海域，高平均网次产量的适宜水深范围为20～50m，平均产量为217～283kg。
2、适应性指数(si)建立
依据表1，分别建立各月份间基于sst、ssha、chl-a和海底水深的适应性指数(表3)。由表3可知，1月si最高的sst、ssha、chl-a和水深分别为17～18℃、-50～-40cm、2.0～5.0mg/m3、30～40m；2月si最高的sst、ssha、chl-a和水深分别为16～17℃、-50～-40cm、2.0～5.0mg/m3、30～40m；3月si最高的sst、ssha、chl-a和水深分别为16～17℃、-50～-40cm、2.0～5.0mg/m3、20～30m；11月si最高的sst、ssha、chl-a和水深分别为19～20℃、0～10cm、0.01～1.0mg/m3、60～70m；12月si最高的sst、ssha、chl-a和水深分别为18～19℃、-40～-30cm、2.0～5.0mg/m3、60～70m。其最适的sst、ssha、chl-a和水深各月有所不同。
表3各月基于sst、ssha、chl-a和海底水深的适应性指数
3、基于栖息地综合指数(hsi)相关因子的权重方案比较
通过与中心渔场相关的不同环境因子设定的权重值的栖息地指数(表3)，按不同hsi来汇总2012-2015年间1-3月、11-12月间的作业网次、产量所占比率和平均网次产量，得出5种方案的平均值(表4)。
由表4可知，5种方案中，方案3的网次、产量比重都是最小的，分别为42.97％和38.53％，平均网次产量相较于其他方案，只有130.17kg，因此可知，方案3权重的设定是最差的。方案2与方案4所得数值相近且都低于方案1和方案5(表4)，因此可知，这两种方案也有所不及。方案1与方案5中hsi大于0.6上的网次、产量比重比较接近，方案1分别为59.69％和60.2％，方案5分别为58.38％和60.96％，但通过表5可以比较发现，平均网次产量以及hsi大于0.8的网次和产量比重相较于方案1更好，因此方案5权重的设定是最佳的。
表45种方案中各月作业网次、产量比率以及平均网次产量的平均值
根据2012-2015年间上海某远洋渔业公司的生产统计数据，结合海表温度(sst)、海面高度距平值(ssha)、叶绿素质量浓度(chl-a)以及水深资料，研究不同权重下摩洛哥头足类栖息地模型。
研究认为，摩洛哥头足类栖息地分布与表温、海面高度距平值、水深等环境因子关系密切，每月适宜的环境因子存在差异；作业渔场分布海域中sst范围是15～23℃，ssha范围是-60～10cm，叶绿素浓度是0～50mg/m3，水深范围是20～100m，而其中最合适的sst是16～18和19～20℃，最合适的ssha是-50～-30cm，最合适的叶绿素含量是1.0～5.0mg/m3，最合适水深是30～40和60～70m。模型分析认为，方案5的权重时最适，其sst、ssha、chl-a和水深的权重因子分别为0.4、0.4、0.1和0.1，说明在栖息地指数模型中，sst和ssha影响最大，水深次之，叶绿素再次之。
虽然以上描述了本发明的具体实施方式，但是本领域的技术人员应当理解，这些仅是举例说明，本发明的保护范围是由所附权利要求书限定的。本领域的技术人员在不背离本发明的原理和实质的前提下，可以对这些实施方式做出多种变更或修改，但这些变更和修改均落入本发明的保护范围。