集合变换卡尔曼滤波局地化(LETKF)方案的发展及应用研究
发布时间:2020-06-04 04:50
【摘要】:集合变换卡尔曼滤波(ETKF)是一种有效的集合预报初始扰动方案,得到广泛应用。但是,有限的集合样本、相同的集合成员设置以及预报模式误差等可能会使ETKF方案中两个距离较远的状态变量产生较高的虚假相关,从而影响集合扰动的质量。为了有效解决远距离虚假相关的问题,将局地化思想引入ETKF方案。本文针对GRAPES区域集合预报系统(GRAPES REPS),研究发展了 ETKF局地化方案(LETKF),确定了 LETKF方案的局地化尺度等关键参数;通过暴雨个例和连续批量试验,从初始扰动的相关分布、能量结构、演变特征以及集合预报综合检验等多方面,分析了 LETKF方案消除远距离虚假扰动的有效性。结果表明,GRAPES区域集合预报中,LETKF局地化能够有效消除远集合扰动的距离虚假相关,进一步改善了集合扰动的质量,更加合理地捕捉到快速增长的分析误差的物理结构,较准确地再现数值模式预报误差的线性与非线性传播和演变特征。主要结论如下:(1)研究发展了 GRAPES REPS的LETKF局地化初始扰动方案。新构建的LETKF局地化方案能够有效消除区域集合扰动的远距离虚假扰动信息,提高集合扰动的整体质量。(2)在当前模式水平分辨率和15个集合成员的条件下,LETKF局地化半径为700 km的结果,对于消除集合扰动的远距离虚假、改进集合扰动质量的效果相对最好。集合扰动的相关分布、能量结构、演变特征证明了局地化方案的有效性及相应扰动结构的合理性。另外,引入低层大气信息合理增大预报的集合离散度,进一步优化了 LETKF局地化方案的性能。(3)基于个例分析和批量试验,结果表明LETKF局地化方案对改进集合预报预报质量,尤其是小雨、中雨、暴雨量级的预报具有较好效果。相比于GRAPES REPS,局地化方案的区域集合预报整体质量具有较明显的优势,尤其是对温度场的改善更为显著。不同物理量、不同层次的集合扰动、离散度的发展与传播特征并不一致。通常,对流层中高层质量较好,低层相对较差。
【图文】:
可以有效解决极地离散格点计算问题。在实际计算离散差分过程中,计算非静力模式的逡逑垂直运动方程容易受到垂直分层设置的影响,尤其是上、下边界条件对其影响更大,基逡逑于此,GRAPES模式采用了邋Charney—Philips跳层非均勾分层设计(图2.2)。逡逑图2.1邋Arakawa-C格点分布设置逡逑Charney-Phillips邋垂直层次逦Lore"K逡逑z,z,,0,w,w逦k邋=邋N逦z,z,w,w逡逑Z,Z,p,P,U,V逦k邋=邋N-Q.5逦z,z,6,p,P,u,v逡逑z,z,6>,vhw逦k^N-\.0逦z,z,w,w逡逑上边界逡逑z,z,p,P,u,\逦k邋=邋N邋-邋\邋.5逦z,z,6,p,P,u,v逡逑?逦?逦?逡逑?逦?逦?逡逑?逦?逦?逡逑z,z,6,w,w逦A:邋=邋2_0逦z,z,w,w逡逑z
确定的模式区域,理想情况是对该区域内每个格点进行局地化求解致计算量大,实际操作难以实现。为减小计算量提高计算效率,将分为N个小方块,每个方块的中心点选定为“局地化中心”(如试验设定合适的局地化半径,利用该局地半径内的观测(图中小圆分布)计算该中心点对应的的变换矩阵。该中心点所在局地区域内同的变换矩阵更新扰动,即由ETKF线性扰动变换关系Xa邋=邋XfTn:逡逑Xa=XfTLn逦(3.4)逡逑h为高度。局地化算子L是包含局地化方案中模式水平空间内局地化半径及其随高度的变化等参数。为了简化问题,这里将随高度hL简化为不随高度变化的定常值。逡逑
【学位授予单位】:南京信息工程大学
【学位级别】:硕士
【学位授予年份】:2018
【分类号】:P456.7
本文编号:2695917
【图文】:
可以有效解决极地离散格点计算问题。在实际计算离散差分过程中,计算非静力模式的逡逑垂直运动方程容易受到垂直分层设置的影响,尤其是上、下边界条件对其影响更大,基逡逑于此,GRAPES模式采用了邋Charney—Philips跳层非均勾分层设计(图2.2)。逡逑图2.1邋Arakawa-C格点分布设置逡逑Charney-Phillips邋垂直层次逦Lore"K逡逑z,z,,0,w,w逦k邋=邋N逦z,z,w,w逡逑Z,Z,p,P,U,V逦k邋=邋N-Q.5逦z,z,6,p,P,u,v逡逑z,z,6>,vhw逦k^N-\.0逦z,z,w,w逡逑上边界逡逑z,z,p,P,u,\逦k邋=邋N邋-邋\邋.5逦z,z,6,p,P,u,v逡逑?逦?逦?逡逑?逦?逦?逡逑?逦?逦?逡逑z,z,6,w,w逦A:邋=邋2_0逦z,z,w,w逡逑z
确定的模式区域,理想情况是对该区域内每个格点进行局地化求解致计算量大,实际操作难以实现。为减小计算量提高计算效率,将分为N个小方块,每个方块的中心点选定为“局地化中心”(如试验设定合适的局地化半径,利用该局地半径内的观测(图中小圆分布)计算该中心点对应的的变换矩阵。该中心点所在局地区域内同的变换矩阵更新扰动,即由ETKF线性扰动变换关系Xa邋=邋XfTn:逡逑Xa=XfTLn逦(3.4)逡逑h为高度。局地化算子L是包含局地化方案中模式水平空间内局地化半径及其随高度的变化等参数。为了简化问题,这里将随高度hL简化为不随高度变化的定常值。逡逑
【学位授予单位】:南京信息工程大学
【学位级别】:硕士
【学位授予年份】:2018
【分类号】:P456.7
【参考文献】
相关期刊论文 前10条
1 智协飞;孙晶;周文友;;2009年夏季西太平洋台风的集合预报和多模式集成预报试验[J];大气科学学报;2015年05期
2 马旭林;时洋;和杰;计燕霞;WANG Yong;;基于卡尔曼滤波递减平均算法的集合预报综合偏差订正[J];气象学报;2015年05期
3 马旭林;于月明;姜胜;李炬;;基于集合卡尔曼变换的目标观测敏感区识别系统优化及影响试验[J];大气科学学报;2014年06期
4 杜钧;李俊;;集合预报方法在暴雨研究和预报中的应用[J];气象科技进展;2014年05期
5 张涵斌;陈静;智协飞;李应林;孙云;;GRAPES区域集合预报系统应用研究[J];气象;2014年09期
6 智协飞;彭婷;李刚;王佳;王晶;;多模式集成的概率天气预报和气候预测研究进展[J];大气科学学报;2014年02期
7 龙柯吉;陈静;马旭林;纪永明;;基于集合卡尔曼变换的区域集合预报初步研究[J];成都信息工程学院学报;2011年01期
8 杜钧;陈静;;单一值预报向概率预报转变的基础:谈谈集合预报及其带来的变革[J];气象;2010年11期
9 马旭林;庄照荣;薛纪善;陆维松;;GRAPES非静力数值预报模式的三维变分资料同化系统的发展[J];气象学报;2009年01期
10 马旭林;薛纪善;陆维松;;GRAPES全球集合预报的集合卡尔曼变换初始扰动方案初步研究[J];气象学报;2008年04期
本文编号:2695917
本文链接:https://www.wllwen.com/projectlw/qxxlw/2695917.html