一种基于集合资料同化技术的初始扰动方法技术

本发明专利技术涉及一种基于集合资料同化技术的初始扰动方法，包括如下步骤：1）将扰动信息构造WRFDA 3D-Var的三个集合分量，通过迭代的方法极小化代价函数来获取统计最优的真实大气状态的估计值X；2）设定不相关的全球集合预报初始条件，形成新的代价函数，3）在WRF3DVAR同化系统中加入新观测；4）从ECMWF的全球尺度集合预报中提取多个包含各个气象要素的垂直廓线；5）为每个集合成员都运行所述新观测的WRF3DVAR同化系统，同化观测集合和大尺度模式场的集合，生成进行集合预报的最终集合成员。有益效果为：采用集合变分同化方法得到的集合成员能够将风暴尺度和大尺度的扰动信息有机地结合起来，代价函数以模式为动力约束，使得集合扰动具有物理和动力的协调性，并且使初始扰动和侧边界扰动相互协调一致。

【技术实现步骤摘要】

本专利技术涉及气象预报，尤其涉及。
技术介绍
过去的一系列研宄表明，风暴尺度集合预报是可行和有效的，但是目前各风暴尺 度集合预报系统的扰动方法还有待完善。如美国风暴分析和预报中心（CAPS)的风暴尺度 集合预报系统，通过提高成员数和分辨率能够提高预报评分，但是其集合扰动方案是固定 的，即各集合成员的初值扰动采用固定的扰动振幅，并且与特定的物理方案匹配；部分集合 成员的侧边界条件由NAM模式短期预报提供，成员间相同的侧边界条件将在一定程度上限 制离散度的发展；另一方面，该系统未考虑侧边界扰动与初值扰动间扰动尺度的不匹配问 题。上述典型问题在各风暴尺度集合预报系统中均有所体现，初值扰动、物理（参数）扰 动、侧边界扰动与初值扰动间的有机结合问题均对于风暴尺度集合预报的成功与否起到关 键作用。 由于初始误差随时间的演变在斜压不稳定和对流不稳定中具有显著地差异，中期 集合扰动方法构造的初始扰动无法在对流系统中快速的增长，也即初始扰动结构与风暴系 统的发展不相适应，进而导致集合成员离散度低于合理的水平；另一方面，风暴尺度的初值 扰动和侧边界扰动的尺度不相适应也将限制初始扰动和集合离散度的增长。多物理过程和 多模式集合有助于改善离散度问题，但是模式的微物理方案中的降水预报变量和微物理过 程的阈值条件判断是不同的，对降水强度和落区范围预报均有显著地影响，针对实时发生 的不同类型风暴系统，采用固定不变的集合配置组合性地选取不同模式微物理方案是盲目 的。
技术实现思路
本专利技术目的在于克服以上现有技术之不足，提供一种基于集合资料同化技术的初 始扰动方法，具体有以下技术方案实现： 所述基于集合资料同化技术的初始扰动方法，基于WRF3DVAR同化系统，包括如下 步骤： 1)将扰动信息构造WRFDA3D_Var的三个集合分量：扰动的观测集合、来自于风暴 尺度集合预报短期预报的背景条件集合以及以观测形式被同化的全球集合预报初始条件 集合，通过迭代的方法极小化代价函数来获取统计最优的真实大气状态的估计值X，X的代 价函数J如式⑴所示， 其中，J表示代价函数，B表示背景误差协方差，R表示被观测误差协方差，y°表示 观测，H(X)表示模式对y°的预报； 2)设定不相关的全球集合预报初始条件，则新的代价函数将表示为： 其中，xls表示准观测矩阵，为xls的经过去相关处理的矩阵，下标i表示第i个 集合成员，H是将Xi插值到X;格点上的线性算子，R#为单位矩阵。 3)在WRF3DVAR同化系统中加入新观测，如式⑶所示，新观测的梯度如式⑷所 示， 4)从ECMWF的全球尺度集合预报中提取多个包含各个气象要素的垂直廓线，并将 所述垂直廓线按WRF3DVAR所用的prebuff格式进行输出； 5)为每个集合成员都运行一次步骤3)中所述新观测的WRF3DVAR同化系统，同化 观测集合和大尺度模式场的集合，生成进行集合预报的最终集合成员。 所述基于集合资料同化技术的初始扰动方法的进一步设计在于，所述新观测中包 括& 3的I/O代码、代价函数的代码、以及开关、检验、和错误判断的功能性的程序代码。 所述基于集合资料同化技术的初始扰动方法的进一步设计在于，所述步骤3)中 乂 13的扰动误差协方差矩阵R通过使用特征值分解方式实现R的逆的计算。 所述基于集合资料同化技术的初始扰动方法的进一步设计在于，所述R的特征值 分解如式（5)、式（6)所示， r= (xlsvTr_1/2)r(xlsvTr_1/2)T (5) R1/2=X1sVt (6) 所述基于集合资料同化技术的初始扰动方法的进一步设计在于，所述步骤1)中 的扰动信息通过网上下载获取，包括变量U、V、T、Q，其中U、V为水平风分量，T为温度，Q为 水汽混合比。 所述基于集合资料同化技术的初始扰动方法的进一步设计在于，所述步骤2)中 包括对X1=做去相关处理操作，如式（7)所示，(7) 所述X〗按WRF3DVAR所用的prebuff格式进行输出，通过WRF3DVAR直接按常规观 测读入X〗数据。 本专利技术的优点如下： 本专利技术构造集合资料同化的三个集合分量：扰动的观测集合；来自于风暴尺度集 合预报短期预报的背景条件集合；以观测形式被同化的全球集合预报初始条件集合。其中， 为准确地得到大尺度信息，将在三维变分资料同化的代价函数中引入新的观测J。项，用于 同化全球集合预报的初始条件。采用集合变分同化方法得到的集合成员能够将风暴尺度和 大尺度的扰动信息有机地结合起来，同时，代价函数以模式为动力约束，使得集合扰动具有 物理和动力的协调性，并且使初始扰动和侧边界扰动相互协调一致。【具体实施方式】 下面对本专利技术方案进行详细说明。 该方法包括如下步骤：1)将扰动信息构造WRFDA3D-Var的三个集合分量：扰动的观测集合、来自于风暴 尺度集合预报短期预报的背景条件集合以及以观测形式被同化的全球集合预报初始条件 集合，通过迭代的方法极小化代价函数来获取统计最优的真实大气状态的估计值X，X的代 价函数J如式⑴所示， 其中，J表示代价函数，B表示背景误差协方差，R表示被观测误差协方差，y°表示 观测，H(X)表示模式对y°的预报； 2)设定不相关的全球集合预报初始条件，则新的代价函数将表示为： 其中，Xls表示准观测矩阵，为Xls的经过去相关处理的矩阵，下标i表示第i个 集合成员，H是将\插值到X；；格点上的线性算子，R#为单位矩阵。 3)在WRF3DVAR同化系统中加入新观测，如式（3)所示，新观测的梯度如式⑷所 示， (4)〇 4)从ECMWF的全球尺度集合预报中提取多个包含各个气象要素的垂直廓线，并将 所述垂直廓线按WRF3DVAR所用的prebuff格式进行输出； 5)为每个集合成员都运行一次步骤3)中所述新观测的WRF3DVAR同化系统，同化 观测集合和大尺度模式场的集合，生成进行集合预报的最终集合成员。 以下给出一具体实施例： (I. 1)首先进入http://nomads.ncdc.noaa.gov/data.php?name=access，查 找GENSONCDC的002格式的集合预报数据，水平分辨率1°。点击FTP或HTTP进入http: // nomads,ncdc.noaa.gov/data/gens/进行下载。网站内数据按年月日分目录存放。一共21 个集合成员，每6小时一次。进行业务预报时，仅需下载最新一个时次的21个成员的数据。 将预报数据划分成WRFDA3D-Var的三个集合分量：扰动的观测集合、来自于风暴尺度集合当前第1页1 2 本文档来自技高网...

【技术保护点】
一种基于集合资料同化技术的初始扰动方法，基于WRF3DVAR同化系统，其特征在于具体包括如下步骤：1)将扰动信息构造WRFDA 3D‑Var的三个集合分量：扰动的观测集合、来自于风暴尺度集合预报短期预报的背景条件集合以及以观测形式被同化的全球集合预报初始条件集合，通过迭代的方法极小化代价函数来获取统计最优的真实大气状态的估计值X，X的代价函数J如式(1)所示，J=XTB-1X+[H~(x)-yo]TR~-1[H~(x)-yo]---(1)]]>其中，J表示代价函数，B表示背景误差协方差，R表示被观测误差协方差，yO表示观测，H(X)表示模式对yO的预报；2)设定不相关的全球集合预报初始条件，则新的代价函数将表示为：J=XTB-1X+[H~(X)-yo]TR~-1[H~(X)-yo]+(R-1/2HXi-Xls*)TR*-1(R-1/2HXi-Xls*)---(2)]]>其中，Xls表示准观测矩阵，为Xls的经过去相关处理的矩阵，下标i表示第i个集合成员，H是将Xi插值到格点上的线性算子，R*为单位矩阵。3)在WRF3DVAR同化系统中加入新观测，如式(3)所示，新观测的梯度如式(4)所示，Jc=(R-1/2HXi-Xls*)TR*-1(R-1/2HXi-Xls*)---(3)]]>▿Jc=(R-1/2H)T(R-1/2HXi-Xls*)---(4).]]>4)从ECMWF的全球尺度集合预报中提取多个包含各个气象要素的垂直廓线，并将所述垂直廓线按WRF3DVAR所用的prebuff格式进行输出；5)为每个集合成员都运行一次步骤3)中所述新观测的WRF3DVAR同化系统，同化观测集合和大尺度模式场的集合，生成进行集合预报的最终集合成员。...

【技术特征摘要】

【专利技术属性】
技术研发人员：闵锦忠，王世璋，王勇，
申请(专利权)人：南京信息工程大学，
类型：发明
国别省市：江苏;32