修正误差扰动的集合Kalman滤波方法在近地表土壤湿度同化中的应用:
模型建立与验证
田向军 谢正辉 师春香
中国科学院大气物理研究所
tianxj@mail.iap.ac.cn
ECCE Summer School for Advanced Study in Climate and Earth Environment
July 30-August 12, 2006, Beijing, China
提纲
引言
集合Kalman滤波算法的实现
简单的土壤水模型
理想验证试验
正在进行的相关工作___
一,引言
土壤湿度在陆气相互作用的重要性
集合Kalman滤波方法在数据同化中的优劣势
二,集合Kalman滤波
集合Kalman滤波是由Evensen(1994b)所开发的一种数据同化方法.其详尽的说明请参阅其相关的文章
集合协方差矩阵P
我们首先定义下面的集合状态矩阵
定义
定义集合扰动矩阵
由此,集合协方差矩阵 便可以定义为
观测扰动
给定一个给定的观测向量 (默认的最优观测). 是观测的个数. 定义 个具有扰动的观测向量如下
分析方程
利用集合斜方差矩阵可以将分析方程写成下面的形式
集合增量(innovation)向量定义为
分析方程可以变为
引入矩阵 ,同时定义矩阵
集合Kalman滤波的一种平方根算法
Evensen(2005)提出了一种新的平方根分析方法. 该方法能够避免观测引入的观测扰动对分析结果的影响,另外该算法在计算分析的结构的时候不需要引入任何额外的附加假设与近似,例如不需要假设观测扰动与集合状态标量扰动的无关性,也不需要对观测误差斜方差矩阵求逆,其算法简便有效.
这种新的算法是从传统的Kalman滤波分析方程的斜方差矩阵的更新出发,进行矩阵变化,得到集合状态变量扰动的更新
算法执行
1._计算出C,然后对C进行特征值分解
2.更新状态变量的集合平均值
按下列顺序进行矩阵乘法
3.计算矩阵
4.进行SVD分解
5.求解状态变量的分析集合扰动
6. 然后再加上集合分析平均值 ,便得到状态变量的分析扰动.
初始集合的生成
如果初始样本的选择适当将有助于准确的表征状态初始量的误差特征.但是,如果选择的样本偏差不是特别大,随着时间的积分一般不会对同化结果产生影响.初始的样本一般采取这样的产生方式:首先对初始状态变量给出一个较为准确的估计值,然后再在这个估计值的附近进行扰动.产生扰动的方式有很多种.最简单的是产生随机数(对标量模型),随机曲线(对1维模型),或者是随机场(对2维或者是以上维数的模型),具体模型的产生方法参加Envensen的文章.
集合的积分
对于一个离散非线性动力模型
三,用于验证模型的简单的土壤水模型
___________四, 理想验证试验
简单验证试验
___________
对上面简单的土壤水模型, 模型参数,初边值,时间步长以及空间步长取下面的值:
利用该模型进行积分,积分1800步作为模式的模拟值. 理想的观测值产生方式如下: 在上面的模型积分中,在每一个时间步随机产生一个均值为0,标准差为0.01的随机扰动加到 上, 也就是在积分的每一步都对土壤水分含量的边界值产生一个随机扰动,由此同样积分1800步,产生一套数据, 把这套数据作为观测数据. 由于实际观测中,我们往往只能得到陆地近地表的土壤湿度观测值. 因此,我们只选取每一步所产生的土壤湿度廓线靠近地面的上面四层土壤的土壤湿度作为每一步的观测数据.
一般通用的集合Kalman滤波方法的应用
利用上面描述的模型,以及上面构造的理想观测数据进行同化试验,其中模型的随机误差取为均值为0,标准差为0.03的随机向量.下面给出的是模拟的结果.
第一层土壤观测模拟与同化的结果
第二层土壤观测模拟与同化的结果
第三层土壤观测模拟与同化的结果
第四层土壤观测模拟与同化的结果
第五层土壤观测模拟与同化的结果
第六层土壤观测模拟与同化的结果
第三十层土壤观测模拟与同化的结果
一般的集合Kalman滤波在第600步的同化
与观测土壤湿度廓线的对比
到底是什么原因造成这么大的误差呢 到底是什么原因使得同化后的土壤湿度廓线如此的奇怪
这是一个问题!
一些猜测
问题的关键在于我们只有上面4层的土壤湿度观测资料,而在生成集合对模式进行扰动的时候,整个土壤30层(用于本文的这个简单的土壤水模型)都给予了0.03标准差的扰动,而没有观测数据进来,使得土壤廓线的曲线比较怪异,在海洋模型中对于这些缺乏数据的点,如果是一个平面上一般采用插值的方法,而我们这里显然不可以进行直接的插值,当然有人采用根据某些特定的经验关系给出下面土壤湿度的"观测值". 而在这里既然没有观测我们就充分相信模式,同时利用上面的观测信息去修正下面的廓线.
改进模型误差扰动的集合Kalman滤波方法的应用
一般的集合集合Kalman滤波同化失败的原因发现,在该试验中(实际情况中也多是如此),只有土壤表层的观测值,也就是说只有在土壤表层有观测信息,对于土壤其他层次的则没有观测信息.而对于一般的集合Kalman滤波而言,模型的模式误差一般取为均值为0,方差为 的随机误差向量(其维数与状态标量向量的维数一致),也就是说对状态向量的每一个变量都进行随机扰动,在我们的这个数值验证中亦是如此,因此对于深层次的土壤湿度而言,较大随机误差的扰动同时缺乏必要的观测信息,从而使得深层次的同化结果非常糟糕.既然缺乏充分的土壤深层次观测信息,
我们只能充分相信模式的模拟结果,同时土壤表层的观测信息也必然通过模式的时间积分去影响深层次的土壤水分含量.鉴于此,我们改变通常使用的模式随机误差扰动的产生方式,只在有观测信息的区域,比如前m层(m是有观测信息的层数),产生随机误差扰动,而对于m+1到n层认为误差为0.下面是改变了模式随机误差扰动后的同化结果分析.
下面是改进误差扰动方式的
模拟的结果
第一层土壤观测模拟与同化的结果
第二层土壤观测模拟与同化的结果
第三层土壤观测模拟与同化的结果
第四层土壤观测模拟与同化的结果
第三十层土壤观测模拟与同化的结果
第五层土壤观测模拟与同化的结果
第六层土壤观测模拟与同化的结果
正在进行的相关工作
建立了基于集合Kalman滤波方法的陆面土壤湿度的同化方案
发现通常的模式误差扰动方案扭曲了整个土壤湿度廓线
改进误差扰动方案后同化结果较为合理
发展了不等距空间网格的数值算法
正在发展一种新的扩展Kalman滤波方法
正在发展一种与CLM格点匹配且能量水量相容的地下水模型
谢谢,请给予指正