问题描述

背景

全球海洋温度会影响某些海洋生物物种的栖息地质量。当温度变化太大而无法继续繁衍生息时,这些物种会迁移到更适合它们现在和未来生活和繁殖成功的其他栖息地。该联盟希望更好地了解如果全球海洋温度升高,苏格兰鲱鱼和鲭鱼可能从苏格兰附近的当前栖息地迁移的相关问题。这两种鱼种代表了对苏格兰渔业的重大经济贡献。鲱鱼和鲭鱼种群位置的变化可能会使苏格兰的小型渔业公司在经济上不切实际,这些公司使用没有船上冷藏的渔船来捕捞新鲜鱼并将其运送到苏格兰渔港的市场。

问题

  • 构建一个数学模型,以确定未来 50 年这两种鱼类最可能的位置。
  • 根据海水温度变化的速度,使用您的模型来预测最佳情况、最坏情况和最可能经过的时间,如果小型渔业公司继续在其当前位置运营,这些种群将距离太远,小型渔业公司无法收获。(Based upon how rapidly the ocean water temperature change occurs, use your model to predict best case, worst case, and most likely elapsed time(s) until these populations will be too far away for small fishing companies to harvest if the small fishing companies continue to operate out of their current locations.)
  • 根据您的预测分析,这些小型渔业公司是否应该改变他们的运营方式?
  • 使用您的模型来解决如果一部分渔业进入另一个国家的领海(海),您的提案会受到怎样的影响。

论文思路

在这篇论文中,作者将这些问题细化,通过一步步的分析原因,给出了每个问题的建模思路。 在建模之前,论文首先收集数据给出了对应海域的平面图,将平面图划分为一个个点域,将原先连续的海域离散化,方便计算机处理与模拟。 对于两种鱼类的位置,论文考虑到题目中影响鱼类迁移的因素为海洋温度,通过收集海域过去的温度变化,得到预测未来海域温度的模型。 之后,从每一条鱼出发,考虑它们在面对海洋不同区域的温度时,选择的迁移方向与速度,并在这个过程中加入随机因素,模拟除温度之外的影响个体迁移方向的因素。 此外,在建模时,个体的迁移变化扩展到某一个区域鱼类密度的变化,通过未来不同时间不同区域的鱼类密度,来体现迁移的结果。

之后,综合考虑温度变化以及鱼类在未来的每一个时间点的迁移趋势,得到问题一的鱼类位置预测模型 对于问题二,论文建立了小型渔业公司的营收模型,考虑捕鱼船的捕鱼范围、单次航程、捕鱼区域的鱼类密度,可以得到一次捕鱼鱼类收获,之后更具市场价格与花费,得到营收。 根据该营收模型,通过10000次的未来鱼类位置情况模拟,计算出小型渔业公司在最好、最坏、最可能情况下的未来经营情况。

根据模型结果,则可以给出适当的建议。在问题四中,限制问题三中计算方式,可重新得到结果并给出建议。

所用模型

ARIMA

论文中通过ARIMA模型,使用过去60年的海洋温度数据来预测未来50年的海洋温度数据。 通过给模型来探索复杂数据的预测情况。

迁移变化模型

对每个点进行动力学和运动学分析,就动力学而言,鱼类运动有温度和速度的影响;就运动学而言,鱼类运动有速度和位置的影响。据此,针对每一个点,可以有变化公式:

可视化如下:

最终,遍历每一个点及每一个要预测的时间点,可有如下算法:

得到未来时间点不同位置的鱼类密度

小型渔业公司收获评估模型

论文中,在建立该模型时,主要考虑了影响收获及营收的因素,考虑捕鱼过程中的支出与收获,将支出分为两类:

  • 固定支出:主要包括每次航行固定的支出,如人员费用等
  • 可变支出:主要包括每次航行中可以变化的费用,如不同航程的费用等

对它们首先定性分析,考虑不同因素的作用效果,将最终的支出与收获相结合,得到净收入的计算模型:

计算机模拟得到未来小型渔业公司的营收情况

将之前的模型综合,编写代码模拟计算10000次的可能迁移结果。 算法如下: