一种基于海洋物理生态耦合模型的参数优化方法与流程

本发明涉及海洋数值模式领域，具体涉及一种基于海洋物理生态耦合模型的生态参数优化方法。

背景技术

我国近海生态环境持续恶化，生态环境安全和海洋经济可持续发展受到严重威胁。一方面，我们亟需研究构建可以预测预报近海高分辨率生态环境要素的海洋物理生态耦合模型，提升海洋防灾减灾能力，保障我国海洋环境安全；另一方面，海洋数值模式一直也是研究海洋现象及机制的重要手段之一，构建能很好再现海洋生态环境的海洋生物地球化学模型是海洋研究中亟待解决的重要问题。海洋是一个自然因素占优的非线性系统，对其的模拟受控于诸如对海洋系统的描述，数值计算方法，参数取值以及输入数据质量等多个方面，然而，随着计算机计算能力的提高、海洋研究的深入以及观测技术的发展，海洋生态模型的参数取值问题日益突出。主要问题表现在以下几个方面：

1.无法通过现代海洋学仪器直接测量生态参数，尽管某些参数(如浮游植物生长率、浮游动物排泄率等)可以通过实验室模拟进行测量，但是由于背景环境不同，实验室模拟不能完全代表真实的海洋生态系统情况，只能给出参数的大致参考范围。

2.即使是同一个参数，模型的复杂程度不同，参数的取值也不同，目前的模型研究中的参数值大多来自文献(他人模型参数)，只能保证参数值只是在数量级上是正确的。

3.模型中的参数具有耦合性，无法精确确定所有参数。

4.模型中许多参数(如浮游植物生长率)不是常数，在空间大尺度上，由于温度、光照等环境因素的不同，参数在不同海区的取值也不尽相同。

5.海洋生态模型应用经验方程将多样性的物种分为可用参数表示的种群，生态模型是非线性的，且含有许多参数，诸多未知量是以参数的形式表现出来的，其取值与模型的描述方法有关。

虽然参数取值受制于其本身的物理含义以及观测局限，但是我们仍能通过对现有的参数参考取值的优化选择获得尽可能理想的结果，提高对海洋的模拟能力。所以有目的针对海洋生态模型进行参数优化十分重要。目前海洋参数优化的方法主要有同化和线性优化两种，但两者都有较高的局限性，前者十分复杂，对海洋模型、计算机计算能力、模型调试者的数学基础以及可用于同化的数据要求较高，后者则因其线性化地处理非线性化生态参数的方法从原理上即具有一定局限性。所以亟需一种简单高效非线性的生态参数优化方法帮助海洋数值模式更好更快地模拟海洋真实情况。

技术实现要素：

有鉴于此，本发明提供了一种基于海洋物理生态耦合模型的生态参数优化方法，目的在于通过一种简单、高效、非线性、半自动化的参数优化方法提高对海洋生态环境的模拟能力。技术方案如下：

一种基于海洋物理生态耦合模型的生态参数优化方法，包括以下步骤：

(1)选取适合研究区域的海洋物理生态耦合模型，通过调试所选取的海洋物理生态耦合模型，获得稳定且较为真实的物理场；获得生态模块中生态参数取值范围，初步调试模型至能再现海洋生态要素尤其是营养盐的空间分布及季节变化趋势；

(2)通过敏感性实验分析所选生态模块中生态参数的敏感度，选择最为敏感的几个参数参与后续参数优化工作以提高后续参数优化效率，方法如下：将每个生态参数依次增大某个确定的百分比，通过计算叶绿素或初级生产力的变化率以及每次增大的参数的变化率之间的比值，对此比值的绝对值进行排序，判断各个参数的敏感程度，随后即基于此结果选取最为敏感的几个参数；

(3)基于步骤(2)中获得的最为敏感的几个参数，结合步骤1所获得的这些敏感参数的取值范围，对这几个参数的取值范围分别进行一定精度的细分采样取值，随后，基于每个参数的采样值，通过排列组合的方法设计多个测试方案，分别应用于步骤1中经过初步调试的海洋物理生态耦合模型，并获得对应测试方案的叶绿素或初级生产力模拟结果；

(4)基于上一步模型运行结果，通过与卫星数据或观测数据进行对比并构造代价函数，代价函数计算公式如下：

其中，m为选择的样本点数量，主要由卫星或观测点数量决定，hⁱ为i点卫星或观测值，yⁱ为模型在对应点的模拟值，当两者的值越接近，模型模拟结果越好，参数优化这一过程可以用减少代价函数这一定量过程描述；基于测试方案的参数取值以及每个测试方案的代价函数值，找到敏感参数取值与代价函数之间的一个对应的多元关系，基于梯度下降方法寻优，得到优化后的一组参数。

本发明的基于海洋物理生态耦合模型的生态参数优化方法，避免了复杂的、对数据及计算机能力要求较高的数值模式同化方法，且通过选择敏感参数减少了参数优化过程中的测试方案数量，提高了效率，此外，还通过敏感参数与构造的代价函数之间的多元关系，交叉机器学习中梯度下降方法思路，非线性、自动寻找最优参数组合，既降低了主观判断的干扰以及人工工作量，也解决了参数线性优化的局限。这是机器学习与海洋数值模式的一次结合与碰撞。

附图说明

图1为一种基于海洋物理生态耦合模型的参数优化方法流程图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面对本发明实施方式作进一步地详细描述。

本发明实施例介绍了一种基于海洋物理生态耦合模型的生态参数优化方法，目的在于通过一种简单高效半自动化的参数优化方法提高对海洋生态环境的模拟能力。具体实施步骤如下：

1.选取并构建一个海洋物理生态耦合模型，可以自行选择适合研究区域的物理模型耦合不同复杂程度的生态模块。首先完成初始场、强迫场、开边界等条件的准备，随后即通过调试物理模型(与生态参数无关)获得稳定且较为真实的物理场(能基本再现海区温度、盐度、潮汐及潮流等要素)，最后通过阅读文献获得生态模块中生态参数取值的大致范围，初步调试模型至能基本再现海洋生态要素尤其是营养盐的空间分布及季节变化趋势(趋势即可，不要求量值一致)。

2.通过敏感性实验分析所选生态模块中生态参数的敏感度，选择最为敏感的几个参数参与后续参数优化工作以提高后续参数优化效率。其中生态参数的敏感度与海洋生态环境中最为基础的要素叶绿素或初级生产力的模拟情况有关，具体计算公式如下：

sc,x＝[(cx-cx+)/cx]/[(x-x+)/x]

其中分子为叶绿素或初级生产力的变化率，分母为生态参数的变化率。

具体操作思路如下：将每个生态参数依次增大50％(可为其他值)，通过计算叶绿素或初级生产力的变化率以及每次增大的参数的变化率之间的比值，对该值的绝对值进行排序，判断各个参数的敏感程度。随后即基于该结果选取最敏感参数。

3.基于步骤2中获得的敏感参数，结合步骤1中从文献获得的这些敏感参数的取值范围，对这几个参数的取值范围分别进行一定精度的细分采样取值。随后，基于每个参数的采样值，通过排列组合的方法设计多个测试方案(case)，分别应用于步骤1中构建的海洋理生态耦合模型并获得对应测试方案的叶绿素或初级生产力模拟结果。

4.基于上一步模型运行结果，通过与卫星数据或观测数据进行对比并构造代价函数，代价函数计算公式如下：

其中，m为选择的样本点数量，因卫星或观测点数据所在区域有限，故该数量主要由卫星或观测点数量决定。h为i点卫星或观测值，y为模型在对应点的模拟值。可以看到，当h与y的值越接近，模型模拟结果越好，也即说明该代价函数可以很好地评价模型模拟结果。因此，参数优化这一过程可以用减少代价函数这一定量过程描述。

随后，基于测试方案的参数取值以及每个测试方案的代价函数值，我们可以找到敏感参数取值与代价函数之间的一个对应的多元关系，该关系是对参数与模拟结果过程的一个简化的描述。基于梯度下降方法的思路，设置一定的步长，我们可以通过该关系进行自动寻优。具体思路即，每个步长内，求取各个参数方向的梯度，梯度下降最快的方向即为下一步参数优化方向，也即该参数往此方向上的调整可以尽可能的减小代价函数。进过一定时间的计算，我们即可得到优化后的一组参数。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。