一种集合卡尔曼滤波静态局地化方案的自适应补偿方法
一种集合卡尔曼滤波静态局地化方案的自适应补偿方法
【技术领域】
[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] 步骤三：顺序地同化所有的观测资料 [0050] 利用一个既定的静态局地化因子以及标准的集合卡尔曼滤波方法，对每个观测资 料进行同化，从而对模式状态的先验集合成员进行更新，得到集合卡尔曼滤波的后验集合 成员。
[0051] 步骤四：利用步骤三中的集合卡尔曼滤波的后验集合成员，计算其集合平均和集 合扰动。
[0052] 步骤五：计算观测余量并更新集合成员
[0053] 将步骤四得到的后验集合平均投影到观测位置，从而得到后验观测值。从原始观 测值中减去后验观测值，得到观测余量，并计算原始观测值和后验观测值之间的均方根误 差（RMSE)。如果RMSE大于步骤二得到的阈值，则：
[0054] 步骤5. 1 :提取观测余量中的多尺度信息
[0055] 利用多重网格分析技术（MGA)提取观测余量中的多尺度信息，从而得到MGA的分 析场。
[0056] 步骤5. 2 :对后验集合平均进行更新
[0057] 将MGA得到的分析场叠加到步骤四得到集合卡尔曼滤波的后验集合平均上，得到 最终的集合平均的分析场。
[0058] 步骤5. 3 :将最终的集合平均分析场叠加到步骤四得到的集合扰动上，得到最终 的集合成员的分析场。
[0059] 本发明提出一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，具体包括以 下几个步骤：
[0060] 步骤一：对实测海洋环境数据进行预处理和质量控制
[0061] 为了尽可能的去除观测数据中包含的错误和不规范数据，需要首先对实测海洋环 境数据进行预处理和质量控制。本专利将通过区域检验、重复深度检验、深度颠倒检验、温 度数据范围检验、温度与盐度梯度检验、密度稳定性检验等六个步骤完成数据的预处理和 质量控制。
[0062] 步骤二：计算自适应补偿方法所需的阈值
[0063] 根据特定的观测系统，按照下式计算自适应补偿方法所需的阈值（0 ):
[0065] 其中，a为显著性水平（通常取为0.01)，K为观测数，r为观测误差标准差， #?(尺)为自由度为K的卡方分布的（1-a)上侧分位数。
[0066] 步骤三：对每个观测资料y°，重复进行以下步骤：
[0067] 步骤3. 1 :计算y°的先验观测集合成员
[0068] 对模式状态向量x的第i个先验集合成员，利用线性插值算子H将其投影到/的 地理位置上，从而得到y°的第i个先验集合成员。
[0069] 步骤3. 2 :计算7°的观测增量
[0070] 根据第3. 1步计算的y°的先验集合、观测值（y°)以及观测误差标准差（r)，按照 下式计算y°的观测增量Ayi:
[0072] 其中，J，f表示y°的第i个先验集合成员，<表示y°的先验标准差，，表示/的先 验集合平均。
[0073] 步骤3. 3 :对状态变量进行更新
[0074]基于一个静态局地化函数Pj，y，根据计算出的y°的观测增量，按照下式对状态变 量进行更新：
[0076]其中j表示第j个状态变量，表示观测与&的先验误差协方差；p&的计算 公式为：
[0078]其中a代表局地化因子，b表示观测y°与状态变量x』之间的物理距离。
[0079] 步骤四：计算后验观测值
[0080] 计算步骤三得到的模式状态的后验集合平均乙并利用观测算子H将其投影到观 测位置，得到后验观测值IE。
[0081] 步骤五：计算观测余量
[0082] 从原始观测向量/中减去后验观测值IK，从而得到观测余量
。此外， 计算原始观测与后验观测之间的均方根误差RMSE。当RMSE大于0时，则顺序进行下列操 作：
[0083]步骤5. 1 :利用一种多重网格分析技术（MGA)，提取观测余量中的多尺度信息， 得到MGA的分析场；
[0084] 步骤5. 2 :将MGA的分析场叠加到步骤三得到的集合平均的分析场上，得到集合平 均最终的分析场；
[0085] 步骤5. 3 :将步骤三得到的集合扰动的分析场叠加到步骤5. 2得到的集合平均的 最终分析场上，得到最终的集合分析场。
【主权项】
1. 一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，其特征是： (1) 对大气海洋环境实测数据进行预处理 剔出实测数据中的无效数据和异常数据，然后将数据转化为同化所需的格式； (2) 对每个观测资料，根据实时的观测系统，利用不同观测要素的观测误差、观测数目 以及显著性水平，计算后续需要使用的阈值； (3) 顺序地同化所有的观测资料： 利用静态局地化因子以及标准的集合卡尔曼滤波方法，对每个观测资料进行同化，对 模式状态的先验集合成员进行更新，得到集合卡尔曼滤波的后验集合成员； (4) 利用步骤（3)中的集合卡尔曼滤波的后验集合成员，计算其集合平均和集合扰动； (5) 计算观测余量并更新集合成员： 将步骤（4)得到的后验集合平均投影到观测位置，从而得到后验观测值，从原始观测 值中减去后验观测值，得到观测余量，并计算原始观测值和后验观测值之间的均方根误差 RMSE，如果RMSE大于步骤（2)得到的阈值，则： (5. 1)提取观测余量中的多尺度信息： 利用多重网格分析技术（MGA)提取观测余量中的多尺度信息，从而得到MGA的分析 场； (5. 2)对后验集合平均进行更新： 将MGA得到的分析场叠加到步骤四得到集合卡尔曼滤波的后验集合平均上，得到最终 的集合平均的分析场； (5.3)将最终的集合平均分析场叠加到步骤（4)得到的集合扰动上，得到最终的集合 成员的分析场。2. 根据权利要求1所述的一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，其 特征在于，所述的步骤（2)包括： 根据特定的观测系统，按照下式计算自适应补偿机制所需的阈值（Θ):其中，α为显著性水平（通常取为〇. 01)，K为观测数，r为观测误差标准差，；&"〇：) 为自由度为K的卡方分布的（l-α )上侧分位数。3. 根据权利要求1所述的一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，其 特征在于，所述的步骤（3)包括： 对每个观测资料y°，重复进行以下步骤： (3. 1)计算y°的先验观测集合成员 对模式状态向量X的第i个先验集合成员，利用线性插值算子H将其投影到y°的地理 位置上，从而得到y°的第i个先验集合成员； (3. 2)计算y°的观测增量 根据第（3. 1)步计算的/的先验集合、观测值（y°)以及观测误差标准差（r)，按照下 式计算y°的观测增量Δ y i:其中，表示y°的第i个先验集合成员，< 表示y°的先验标准差，，表示y°的先验 集合平均； (3. 3)对状态变量进行更新 基于一个静态局地化函数P U，根据计算出的y°的观测增量，按照下式对状态变量进 行更新：其中j表示第j个状态变量，covL表示观测与&的先验误差协方差；P u的计算公式 为：其中a代表局地化因子，b表示观测y°与状态变量X」之间的物理距离。
【专利摘要】本发明属于实测海洋环境数据同化领域，具体涉及一种集合卡尔曼滤波静态局地化方案的自适应补偿方法。本发明包括：对大气海洋环境实测数据进行预处理；对每个观测资料，根据实时的观测系统，利用不同观测要素的观测误差、观测数目以及显著性水平，计算后续需要使用的阈值；顺序地同化所有的观测资料；计算其集合平均和集合扰动；计算观测余量并更新集合成员。本发明对集合卡尔曼滤波中的静态局地化方法进行了改进，有效地提取了观测信息中静态局地化方法未能提取的多尺度信息，显著提高了集合卡尔曼滤波的同化精度。
【IPC分类】G06F19/00
【公开号】CN104899448
【申请号】CN201510312199
【发明人】赵玉新, 吴新荣, 刘厂, 付红丽, 刘利强, 王喜冬, 李刚, 张连新, 张晓爽, 张振兴
【申请人】哈尔滨工程大学, 国家海洋信息中心
【公开日】2015年9月9日
【申请日】2015年6月9日
【技术领域】
[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] 步骤三：顺序地同化所有的观测资料 [0050] 利用一个既定的静态局地化因子以及标准的集合卡尔曼滤波方法，对每个观测资 料进行同化，从而对模式状态的先验集合成员进行更新，得到集合卡尔曼滤波的后验集合 成员。
[0051] 步骤四：利用步骤三中的集合卡尔曼滤波的后验集合成员，计算其集合平均和集 合扰动。
[0052] 步骤五：计算观测余量并更新集合成员
[0053] 将步骤四得到的后验集合平均投影到观测位置，从而得到后验观测值。从原始观 测值中减去后验观测值，得到观测余量，并计算原始观测值和后验观测值之间的均方根误 差（RMSE)。如果RMSE大于步骤二得到的阈值，则：
[0054] 步骤5. 1 :提取观测余量中的多尺度信息
[0055] 利用多重网格分析技术（MGA)提取观测余量中的多尺度信息，从而得到MGA的分 析场。
[0056] 步骤5. 2 :对后验集合平均进行更新
[0057] 将MGA得到的分析场叠加到步骤四得到集合卡尔曼滤波的后验集合平均上，得到 最终的集合平均的分析场。
[0058] 步骤5. 3 :将最终的集合平均分析场叠加到步骤四得到的集合扰动上，得到最终 的集合成员的分析场。
[0059] 本发明提出一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，具体包括以 下几个步骤：
[0060] 步骤一：对实测海洋环境数据进行预处理和质量控制
[0061] 为了尽可能的去除观测数据中包含的错误和不规范数据，需要首先对实测海洋环 境数据进行预处理和质量控制。本专利将通过区域检验、重复深度检验、深度颠倒检验、温 度数据范围检验、温度与盐度梯度检验、密度稳定性检验等六个步骤完成数据的预处理和 质量控制。
[0062] 步骤二：计算自适应补偿方法所需的阈值
[0063] 根据特定的观测系统，按照下式计算自适应补偿方法所需的阈值（0 ):
[0065] 其中，a为显著性水平（通常取为0.01)，K为观测数，r为观测误差标准差， #?(尺)为自由度为K的卡方分布的（1-a)上侧分位数。
[0066] 步骤三：对每个观测资料y°，重复进行以下步骤：
[0067] 步骤3. 1 :计算y°的先验观测集合成员
[0068] 对模式状态向量x的第i个先验集合成员，利用线性插值算子H将其投影到/的 地理位置上，从而得到y°的第i个先验集合成员。
[0069] 步骤3. 2 :计算7°的观测增量
[0070] 根据第3. 1步计算的y°的先验集合、观测值（y°)以及观测误差标准差（r)，按照 下式计算y°的观测增量Ayi:
[0072] 其中，J，f表示y°的第i个先验集合成员，<表示y°的先验标准差，，表示/的先 验集合平均。
[0073] 步骤3. 3 :对状态变量进行更新
[0074]基于一个静态局地化函数Pj，y，根据计算出的y°的观测增量，按照下式对状态变 量进行更新：
[0076]其中j表示第j个状态变量，表示观测与&的先验误差协方差；p&的计算 公式为：
[0078]其中a代表局地化因子，b表示观测y°与状态变量x』之间的物理距离。
[0079] 步骤四：计算后验观测值
[0080] 计算步骤三得到的模式状态的后验集合平均乙并利用观测算子H将其投影到观 测位置，得到后验观测值IE。
[0081] 步骤五：计算观测余量
[0082] 从原始观测向量/中减去后验观测值IK，从而得到观测余量
。此外， 计算原始观测与后验观测之间的均方根误差RMSE。当RMSE大于0时，则顺序进行下列操 作：
[0083]步骤5. 1 :利用一种多重网格分析技术（MGA)，提取观测余量中的多尺度信息， 得到MGA的分析场；
[0084] 步骤5. 2 :将MGA的分析场叠加到步骤三得到的集合平均的分析场上，得到集合平 均最终的分析场；
[0085] 步骤5. 3 :将步骤三得到的集合扰动的分析场叠加到步骤5. 2得到的集合平均的 最终分析场上，得到最终的集合分析场。
【主权项】
1. 一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，其特征是： (1) 对大气海洋环境实测数据进行预处理 剔出实测数据中的无效数据和异常数据，然后将数据转化为同化所需的格式； (2) 对每个观测资料，根据实时的观测系统，利用不同观测要素的观测误差、观测数目 以及显著性水平，计算后续需要使用的阈值； (3) 顺序地同化所有的观测资料： 利用静态局地化因子以及标准的集合卡尔曼滤波方法，对每个观测资料进行同化，对 模式状态的先验集合成员进行更新，得到集合卡尔曼滤波的后验集合成员； (4) 利用步骤（3)中的集合卡尔曼滤波的后验集合成员，计算其集合平均和集合扰动； (5) 计算观测余量并更新集合成员： 将步骤（4)得到的后验集合平均投影到观测位置，从而得到后验观测值，从原始观测 值中减去后验观测值，得到观测余量，并计算原始观测值和后验观测值之间的均方根误差 RMSE，如果RMSE大于步骤（2)得到的阈值，则： (5. 1)提取观测余量中的多尺度信息： 利用多重网格分析技术（MGA)提取观测余量中的多尺度信息，从而得到MGA的分析 场； (5. 2)对后验集合平均进行更新： 将MGA得到的分析场叠加到步骤四得到集合卡尔曼滤波的后验集合平均上，得到最终 的集合平均的分析场； (5.3)将最终的集合平均分析场叠加到步骤（4)得到的集合扰动上，得到最终的集合 成员的分析场。2. 根据权利要求1所述的一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，其 特征在于，所述的步骤（2)包括： 根据特定的观测系统，按照下式计算自适应补偿机制所需的阈值（Θ):其中，α为显著性水平（通常取为〇. 01)，K为观测数，r为观测误差标准差，；&"〇：) 为自由度为K的卡方分布的（l-α )上侧分位数。3. 根据权利要求1所述的一种集合卡尔曼滤波静态局地化方案的自适应补偿方法，其 特征在于，所述的步骤（3)包括： 对每个观测资料y°，重复进行以下步骤： (3. 1)计算y°的先验观测集合成员 对模式状态向量X的第i个先验集合成员，利用线性插值算子H将其投影到y°的地理 位置上，从而得到y°的第i个先验集合成员； (3. 2)计算y°的观测增量 根据第（3. 1)步计算的/的先验集合、观测值（y°)以及观测误差标准差（r)，按照下 式计算y°的观测增量Δ y i:其中，表示y°的第i个先验集合成员，< 表示y°的先验标准差，，表示y°的先验 集合平均； (3. 3)对状态变量进行更新 基于一个静态局地化函数P U，根据计算出的y°的观测增量，按照下式对状态变量进 行更新：其中j表示第j个状态变量，covL表示观测与&的先验误差协方差；P u的计算公式 为：其中a代表局地化因子，b表示观测y°与状态变量X」之间的物理距离。
【专利摘要】本发明属于实测海洋环境数据同化领域，具体涉及一种集合卡尔曼滤波静态局地化方案的自适应补偿方法。本发明包括：对大气海洋环境实测数据进行预处理；对每个观测资料，根据实时的观测系统，利用不同观测要素的观测误差、观测数目以及显著性水平，计算后续需要使用的阈值；顺序地同化所有的观测资料；计算其集合平均和集合扰动；计算观测余量并更新集合成员。本发明对集合卡尔曼滤波中的静态局地化方法进行了改进，有效地提取了观测信息中静态局地化方法未能提取的多尺度信息，显著提高了集合卡尔曼滤波的同化精度。
【IPC分类】G06F19/00
【公开号】CN104899448
【申请号】CN201510312199
【发明人】赵玉新, 吴新荣, 刘厂, 付红丽, 刘利强, 王喜冬, 李刚, 张连新, 张晓爽, 张振兴
【申请人】哈尔滨工程大学, 国家海洋信息中心
【公开日】2015年9月9日
【申请日】2015年6月9日
最新回复(0)