一种集合卡尔曼滤波静态局地化方案的自适应补偿方法
【技术领域】
[0001] 本发明属于实测海洋环境数据同化领域，具体涉及一种集合卡尔曼滤波静态局地 化方案的自适应补偿方法。
【背景技术】
[0002] 集合卡尔曼滤波和四维变分是目前国际上公认的两类先进的数据同化方法，它们 有各自的优缺点。集合卡尔曼滤波相对于变分方法的最大优点是它通过集合抽样来模拟模 式状态变量的先验概率密度分布函数，根据集合样本计算出的背景误差协方差矩阵带有模 式的动力信息，因此是流依赖的。由于计算机硬件资源的限制，目前只能采用较少的集合样 本（1〇 2量级）。而对实际的海洋数值模式来说，状态变量的维数为10 7,因此较少的集合样 本就会带来显著的抽样误差。
[0003] 本发明主要是对集合卡尔曼滤波中传统的静态局地化方案进行改进，提高集合卡 尔曼滤波的同化效果。根据前述，有限的集合样本会给集合卡尔曼滤波数据同化引入显著 的抽样误差，比如模式状态先验方差的低估以及模式状态与远距离观测之间的伪相关性 等。解决这个问题的一个有效途径是在同化过程中引入局地化方法，即一个观测仅允许影 响周围一定范围内的模式格点，且距离越远影响越小。传统的集合卡尔曼滤波局地化方案 是一种静态方法，即预先给定一个观测影响半径（即局地化因子），在整个同化过程中均采 用这个局地化因子。大量的文献记载以及试验证明，集合卡尔曼滤波的同化效果对局地化 因子的选取非常敏感。这种静态的局地化方案存在明显的局限性，因为对于实际的三维海 洋数值模式来说，要找到一个最优的局地化因子是非常耗时的，几乎是不可能的。因此，针 对静态局地化方法存在的弊端，引入一种补偿方法，从而在保证同化精度的同时降低同化 效果对局地化因子的敏感性是非常有意义和有实用价值的。
【发明内容】
[0004] 本发明的目的在于提供一种提高精度的集合卡尔曼滤波静态局地化方案的自适 应补偿方法。
[0005] 本发明的目的是这样实现的：
[0006] (1)对大气海洋环境实测数据进行预处理
[0007] 剔出实测数据中的无效数据和异常数据，然后将数据转化为同化所需的格式；
[0008] (2)对每个观测资料，根据实时的观测系统，利用不同观测要素的观测误差、观测 数目以及显著性水平，计算后续需要使用的阈值；
[0009] (3)顺序地同化所有的观测资料：
[0010] 利用静态局地化因子以及标准的集合卡尔曼滤波方法，对每个观测资料进行同 化，对模式状态的先验集合成员进行更新，得到集合卡尔曼滤波的后验集合成员；
[0011] (4)利用步骤（3)中的集合卡尔曼滤波的后验集合成员，计算其集合平均和集合 扰动；
[0012] (5)计算观测余量并更新集合成员：
[0013] 将步骤（4)得到的后验集合平均投影到观测位置，从而得到后验观测值，从原始 观测值中减去后验观测值，得到观测余量，并计算原始观测值和后验观测值之间的均方根 误差RMSE，如果RMSE大于步骤（2)得到的阈值，则：
[0014] (5. 1)提取观测余量中的多尺度信息：
[0015] 利用多重网格分析技术（MGA)提取观测余量中的多尺度信息，从而得到MGA的分 析场；
[0016] (5. 2)对后验集合平均进行更新：
[0017] 将MGA得到的分析场叠加到步骤四得到集合卡尔曼滤波的后验集合平均上，得到 最终的集合平均的分析场；
[0018] (5. 3)将最终的集合平均分析场叠加到步骤（4)得到的集合扰动上，得到最终的 集合成员的分析场。
[0019] 所述的步骤（2)包括：
[0020] 根据特定的观测系统，按照下式计算自适应补偿机制所需的阈值（0 ):
[0022] 其中，a为显著性水平（通常取为0.01)，K为观测数，r为观测误差标准差， /"(幻为自由度为K的卡方分布的（1-a)上侧分位数。
[0023] 所述的步骤（3)包括：
[0024] 对每个观测资料y°，重复进行以下步骤：
[0025] (3. 1)计算y°的先验观测集合成员
[0026] 对模式状态向量x的第i个先验集合成员，利用线性插值算子H将其投影到/的 地理位置上，从而得到y°的第i个先验集合成员；
[0027] (3. 2)计算y°的观测增量
[0028] 根据第（3. 1)步计算的y°的先验集合、观测值（y°)以及观测误差标准差（r)，按 照下式计算y°的观测增量Ayi:
[0030] 其中，表示y°的第i个先验集合成员，<表示y°的先验标准差，F表示y°的先 验集合平均；
[0031] (3. 3)对状态变量进行更新
[0032]基于一个静态局地化函数Pj，y，根据计算出的y°的观测增量，按照下式对状态变 量进行更新：
[0034]其中j表示第j个状态变量，COV〗.,.表示观测与&的先验误差协方差；pU的计算 公式为：
[0036] 其中a代表局地化因子，b表示观测y°与状态变量x」之间的物理距离。
[0037] 本发明的有益效果在于：
[0038] (1)对集合卡尔曼滤波中的静态局地化方法进行了改进，有效地提取了观测信息 中静态局地化方法未能提取的多尺度信息，显著提高了集合卡尔曼滤波的同化精度。
[0039] (2)降低了集合卡尔曼滤波对静态局地化因子的强依赖性。
【附图说明】
[0040] 图1 一种集合卡尔曼滤波静态局地化方案的自适应补偿方法的执行流程图。
【具体实施方式】
[0041] 下面结合附图对本发明做进一步描述。
[0042] 本发明提供的是一种集合卡尔曼滤波静态局地化方案的自适应补偿机制。包括 对大气海洋环境实测数据进行预处理；根据实时的观测系统，利用不同观测要素的观测误 差、观测数目以及显著性水平，计算后续需要使用的阈值；利用集合卡尔曼滤波方法顺序地 同化所有的观测资料；计算集合卡尔曼滤波的后验集合平均和集合扰动；计算观测余量并 更新集合成员。本发明对集合卡尔曼滤波中的静态局地化方法进行了改进，有效地提取了 观测信息中静态局地化方法未能提取的多尺度信息，显著提高了集合卡尔曼滤波的同化精 度。降低了集合卡尔曼滤波对静态局地化因子的强依赖性。本发明对实现海洋环境状态场 的实时构建以及预测具有重要的实用价值。
[0043] 本发明针对集合卡尔曼滤波中的静态局地化方案提出一种自适应补偿方法。与传 统的静态局地化技术相比，本发明的显著：该方法能够自适应地提取观测资料中静态局地 化方案未能提取的观测信息。对每个分析步，具体实现过程为：首先，基于一个静态局地化 因子，利用观测资料对模式的先验集合成员进行调整；其次，将模式状态的后验集合平均值 投影到观测位置上，得到观测的后验值；第三，从原始观测值中减去上述后验观测值，得到 观测余量，并计算原始观测与后验观测之间的均方根误差；第四，当上述均方根误差大于某 个阈值时，则利用一种多重网格技术（MGA)提取观测余量中的多尺度信息，从而得到相应 的分析场；第五，当上述均方根误差大于某个阈值时，将MGA的分析场叠加到第一步得到后 验集合平均上得到最终的集合卡尔曼滤波的集合平均分析场。本专利提出的方法可以大大 地改进集合卡尔曼滤波中的传统的静态局地化方案的同化效果。
[0044] 每个分析步的同化过程具体包括以下几个步骤：
[0045] 步骤一：对大气海洋环境实测数据进行预处理
[0046] 为了将实测环境数据同化到背景场数据库中，首先必须剔出实测数据中的无效数 据和异常数据，然后将数据转化为同化所需的格式。
[0047] 对每个观测资料，重复进行以下步骤：
[0048] 步骤二：根据实时的观测系统，利用不同观测要素的观测误差、观测数目以及显著 性水平，计算后续需要使用的阈值。
[0049] 步骤三：顺序地同化所有的观测资料