基于类比的离线集合卡尔曼滤波古气候数据同化系统、方法及电子设备

技术领域

本发明涉及一种古气候场重建方法，是一种将基于类比与集合卡尔曼滤波的同化系统、方法及电子设备。

背景技术

理解古气候现象能够提升对未来气候变化的预报技巧，古气候场的重建使用古气候代用指标用于评估气候状态。多数古气候场重建方法运用多元线性回归得到气候场变量，但是受限制于代用指标和气候变量间的非线性关系。

相应地，古气候数据同化是一种结合气候模式动力系统限制和代用观测信息，找到气候系统状态的最优估计。其中，传统的基于集合的数据同化技术使用前一个时刻通过观测更新的分析场集合作为集合模式的初始条件向前积分得到下一个循环同化的先验集合，被称为‘在线’循环集合同化。但是由于模式的可预报性短于观测的时间尺度，同时集合模式积分需要耗费巨大的计算量。一种‘离线’集合数据同化方法被提出，通过已存在的气候模拟中构造出背景场集合。‘离线’与‘在线’被系统地比较过，‘在线’方法并未发现比‘离线’方法有显著的提升。

因此离线的基于集合的古气候数据同化方法被广泛应用于古气候场重建中，先验集合成员常常是从已存在的长时间气候模拟中随机抓取，同时为了使重建的气候场不被每次抓取的随机性影响，相同的抓取集合在每一次同化中被重复使用。另一种基于集合的古气候数据同化方法是使用退化的粒子滤波，从集合模式的模拟中选择最优的模拟，但是即使是具备较大的集合数目，被选择的粒子也会受到维度诅咒的影响。

发明内容

本发明针对现有技术的不足，提供了一种使用基于类比与集合卡尔曼滤波的同化方法(以下简称为AOEnKF)。AOEnKF先在一组长时间气候模拟中选择先验集合成员，选择标准是集合场与观测之间的均方根误差(RMSE)最小或者相关系数最大。得到的先验集合成员因此捕获了‘流依赖’的背景场误差协方差矩阵，再将气候代用资料同化进入先验集合成员中，即可得到后验集合成员。通过本发明所述的AOEnKF得到的先验集合成员，相比于‘在线’循环同化来说，不需要集合模式积分即可得到‘流依赖’的背景场误差协方差矩阵，计算量更小；相比于传统的‘离线’基于集合同化方法的静态的先验集合成员来说，能够更好地将代用资料信息传播给气候场。由此可知，本发明所述的同化方法能结合类比方法与集合卡尔曼滤波的优势，更好地进行古气候场的重建。

为实现上述技术目的，本发明将采取如下技术方案：

一种基于类比的离线集合卡尔曼滤波古气候数据同化方法，先根据特定标准，利用观测y在一组控制实验模拟的状态变量样本x中筛选先验集合成员，再根据集合平方根滤波的更新公式，将得到的先验集合成员通过卡尔曼增益矩阵更新后验集合平均

步骤一、给定一组控制实验模拟的状态变量样本x及观测y

步骤二、在同化观测前根据不同标准选择先验集合成员

在同化观测前，先将状态变量样本x插值到观测y中，然后依据状态变量样本x的每一个状态变量样本x

步骤三、用得到的先验集合成员同化观测

使用平方根滤波同化误差协方差矩阵为R的观测y，先更新后验集合平均

优选地，当从状态变量样本x中选择先验集合成员的标准为状态变量样本x 的每一个状态变量样本x

其中：观测y表示p×1的观测变量；y

同时，N个先验集合成员的各个样本按照对应的均方根误差RMSE升序排序。

优选地，当从状态变量样本x中选择先验集合成员的标准为状态变量样本x 的每一个状态变量样本x

其中：观测y表示p×1的观测变量；y′表示观测y减去其空间平均得到的观测空间异常场；y′

同时，N个先验集合成员的各个样本按照对应的相关系数降序进行排序。

后验集合平均

后验集合扰动

式中：

X

本发明的另一个技术目的是提供一种基于类比的离线集合卡尔曼滤波古气候数据同化系统，包括：

获取模块，用于获取需要同化的误差协方差矩阵为R的观测y；

同化模块，基于集合卡尔曼滤波的同化框架而搭建；

所述同化模块在同化观测前，先将状态变量样本x插值到观测y中，然后依据状态变量样本x的每一个状态变量样本x

所述同化模块在同化时刻，使用平方根滤波同化误差协方差矩阵为R的观测 y，先更新后验集合平均

优选地，所述同化模块从状态变量样本x中选择先验集合成员的标准为状态变量样本x的每一个状态变量样本x

其中：观测y表示p×1的观测变量；y′表示观测y减去空间平均得到的观测空间异常场；y′

优选地，所述同化模块从状态变量样本x中选择出的N个先验集合成员的各个样本，按照每一个状态变量样本x

其中：观测y表示p×1的观测变量；y

优选地，所述同化模块在同化时刻的后验集合平均

式中：

X

本发明的第三个技术目的是提供一种电子设备，包括：至少一个处理器、至少一个存储器、通信接口和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于该电子设备和其它电子设备的通信设备之间的信息传输；所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述的方法。

根据上述的技术方案，相对于现有技术，本发明具有如下的有益效果：

本发明基于类比的方法，通过一定的标准(相对于观测的RMSE或者相关系数)选择集合先验成员，代替传统的随机抓取的静态集合先验成员来更新集合均值和集合扰动。相对于‘在线’循环同化的方法减少了集合模式向前积分的巨大计算量。与通常的‘离线’重复利用随机抓取的先验集合的同化方法相比，本发明能够构造更精确的先验集合平均，同时捕获‘流依赖’的背景场误差协方差信息，能够进一步减小同化后的误差。

附图说明

图1为本发明的流程示意图。

图2为‘在线’循环集合卡尔曼滤波(CEnKF)、‘离线’集合卡尔曼滤波 (OEnKF)、基于RMSE标准类比的离线集合卡尔曼滤波(AOEnKF_E)与基于相关系数标准类比的离线集合卡尔曼滤波(AOEnKF_C)方法RMSE随不同局地化值变化的序列：圆圈表示OEnKF与其他试验的RMSE显著不同(通过了99％的置信检验)，叉号表示AOEnKF_E与AOEnKF_C的RMSE的显著不同(通过了99％的置信检验)。

图3是本发明所述电子设备的结构示意图。

具体实施方式

为了更好的理解本发明的技术方案，下面结合附图对本发明实施例进行详细描述。

应当明确，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。

在本发明实施例中使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明实施例和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。

应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，甲和/或乙，可以表示：单独存在甲，同时存在甲和乙，单独存在乙这三种情况。另外，本文中字符“/”，一般表示前后关联对象是一种“或”的关系。

实施例1

如图1所示，本发明所述的基于类比的离线集合卡尔曼滤波古气候数据同化方法，先利用观测在一组控制实验模拟的状态变量样本中根据不同的标准(模式与观测间RMSE最小或者相关系数最大)筛选先验集合成员，再用得到的先验集合成员应用集合平方根滤波进行集合平均与集合扰动的更新。应用本同化方法可以进一步改进分析场。具体步骤如下：

步骤一、给定一组控制实验模拟的状态变量及观测

1.1.给定一组控制实验模拟的状态变量样本

x

1.2.给定需要同化的观测

给定误差协方差矩阵为R的观测y，一般假设观测误差互不相关，R是一个对角矩阵，对角线元素为观测变量的误差方差。

步骤二、在同化观测前根据不同标准选择先验集合成员

在同化观测前，先将状态变量插值到观测空间，然后基于最小均方根误差 (RMSE)的标准或者基于状态变量与观测之间的相关系数最大的原则构建先验集合成员

2.1.基于最小均方根误差标准选择先验集合成员(简称为AOEnKF_E)

AOEnKF_E是选择先验集合成员

通过升序排序RMSE,最初的N个样本被选为此次同化时间里的先验集合成员

2.2.基于最大相关系数标准选择先验集合成员(简称为AOEnKF_C)

AOEnKF_C选择集合先验成员

其中撇号表示通过减去空间平均得到的空间异常场。通过将相关系数进行降序排序，最初的N个样本被选为此次同化时间里的先验集合成员

2.3.选择合适的局地化方案

选择合适的局地化函数来局地化观测的影响，减少观测与状态变量之间的虚假相关，常用的局地化函数如Gaspari和Cohn(GC；Gaspari和Cohn 1999) 函数，由单个特征尺度参数确定。局地化函数在同化时如步骤(3.1)所示应用到卡尔曼增益矩阵上。

步骤三、用得到的先验集合成员同化观测

使用平方根滤波同化误差协方差矩阵为R的观测y，先更新后验集合成员的集合平均

3.1.更新集合平均

后验集合成员的集合平均

其中：ρ是局地化矩阵，

X

3.2.更新集合扰动

第i(i＝1,…,N)个后验集合成员的扰动更新公式为：

其中：

3.3.得到后验集合成员

后验集合平均与后验集合扰动相加得到N个后验集合成员。

3.4.评估同化实验结果

当同化实验结束后得到均方根误差的时间序列进行平均作为标准评估同化实验的结果。

实施例2

本发明基于类比的离线集合卡尔曼滤波古气候数据同化方法同化观测，以Lorenz(2005)模型为例，以单尺度模式II无模式误差测试了本发明的表现，并与使用传统静态抓取，‘离线’集合卡尔曼滤波(OEnKF，Hakim 2016)方法的误差结果进行了比较。敏感试验结果表明，本发明在不同的集合大小、局地化尺度、观测误差、观测密度都优于传统的同化方法。

步骤一、给定一组控制实验模拟的状态变量及观测

L05模型单尺度模式II只包含一个大尺度慢过程变量。设X为慢过程变量，单尺度模式II可写为：

下标n代表格点序号，K为常数，F为强迫项。平流项写为：

其中∑'为一个特殊的求和算子，除首项与末项除以2外，与一般求和算子相同。 K为偶数时J＝K/2，

总格点数N＝960，常数K选择32。强迫项F为16。真值与控制实验的模式参数相同，表示控制实验中无模式误差存在。

1.1.给定一组控制实验模拟的状态变量样本

x

1.2.给定需要同化的观测

给定误差协方差矩阵为R的观测y，一般假设观测误差互不相关，R是一个对角矩阵，对角线元素为观测变量的误差方差。通过向一段时间真值的平均(～120h)添加服从正态分布N(0,R)的随机扰动来制造观测。默认的观测误差方差大小为4.0。默认的观测网络为每8个格点一个观测(共120个观测格点)。每6000个积分时步的平均(～120h)产生一个观测。

步骤二、在同化观测前根据不同标准选择先验集合成员

2.1.基于最小均方根误差标准选择先验集合成员

AOEnKF_E是选择先验集合成员

通过升序排序RMSE,最初的N个样本被选为此次同化时间里的先验集合成员

2.2.基于最大相关系数标准选择先验集合成员

AOEnKF_C选择集合先验成员

其中撇号表示通过减去空间平均得到的空间异常场。通过将相关系数进行降序排序，最初的N个样本被选为此次同化时间里的先验集合成员

2.3.选择合适的局地化方案

Gaspari和Cohn(GC；Gaspari和Cohn 1999)函数被用作局地化函数，由单个特征尺度参数确定。局地化函数在同化时如步骤(3.1)所示应用到卡尔曼增益矩阵上。GC函数的特征尺度需要选择合适的值以获得最好的结果。

步骤三、用得到的先验集合成员同化观测

已有通过AOEnKF_E或者AOEnKF_C选择得到了N个集合先验成员

3.1.更新集合平均

后验集合平均

其中ρ是局地化矩阵，

3.2.更新集合扰动

第i(i＝1,…,N)个集合成员的扰动更新公式为

3.3.得到后验集合成员

后验集合平均与后验集合扰动相加得到N个后验集合成员。

3.4.评估同化实验结果

当同化实验结束后得到均方根误差的时间序列进行平均作为标准评估同化实验的结果。每种数据同化方法进行了1000月的同化实验，图2显示了各同化方法的RMSE1000个月时间序列平均随不同局地化值的变化。CEnKF与OEnKF有着几乎相同的误差，这证明了单尺度L05模式对于初始条件几乎没有任和记忆性，模式的可预报性相比于观测的时间尺度短。而对于AOEnKFs，AOEnKF_E与 AOEnKF_C具有类似的表现，RMSE相比于OEnKF显著减小。这证明了在先验集合成员中包含的‘流依赖’的背景场误差协方差信息能够提升同化的效果。

实施例3

基于上述的古气候数据同化方法，本发明提供一种基于类比的离线集合卡尔曼滤波古气候数据同化系统，包括：获取模块，用于获取需要同化的误差协方差矩阵为R的观测y；同化模块，基于集合卡尔曼滤波的同化框架而搭建；所述同化模块在同化观测前，先将状态变量样本x插值到观测y中，然后依据状态变量样本x的每一个状态变量样本x

同化模块能够实现本发明古气候数据同化方法实施例中能够实现的各个过程，以及达到相同的有益效果，为避免重复，在此不再赘述。

实施例四

基于上述的古气候数据同化方法，本发明提供一种电子设备，包括：至少一个处理器、至少一个存储器、通信接口和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于该电子设备和其它电子设备的通信设备之间的信息传输；所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述的古气候数据同化方法。

具体地，处理器用于获取需要同化的误差协方差矩阵为R的观测y；

处理器内具有基于集合卡尔曼滤波的同化框架而搭建的同化模块；同化模块基于特定标准，利用所获取的观测y在一组控制实验模拟的状态变量样本x中筛选先验集合成员，再根据集合平方根滤波的更新公式，将得到的先验集合成员通过卡尔曼增益矩阵更新后验集合平均

本实施例中，处理器内的同化模块，同样能够实现本发明古气候数据同化方法实施例中能够实现的各个过程，以及达到相同的有益效果，为避免重复，在此不再赘述。

上面结合附图对本申请的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本申请的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，均属于本发明的保护之内。