[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日