一种面向资料同化的天气雷达径向速度质量控制方法技术

本发明专利技术公开了一种面向资料同化的天气雷达径向速度质量控制方法。该方法基于双权重算法，在参数列表中设置了科研和业务的选项，通过自主设定质量控制阈值的方式满足科研工作中操作直观、执行高效的需求，通过递归质量控制的方式实现业务工作中效果稳定、无需人工干扰的需求。该方法合理剔除观测误差较大的资料，去除模式分辨率无法分辨的小尺度变化资料及难以模拟好的资料，保证实际观测值减去背景场模拟值的概率密度函数接近高斯分布，提高了分析场的精度，有利于后续的雷达资料同化。该方法充分考虑到数值天气预报中资料同化的实际需求，并结合数值天气预报的常规流程，直接集成到数值天气预报模式中，具有简便易行、执行快速、维护方便等优点。

A data assimilation oriented radial velocity quality control method for weather radar

The invention discloses a weather radar radial velocity quality control method for data assimilation. Based on the double weight algorithm, this method sets the options of scientific research and business in the parameter list. By setting the quality control threshold independently, it can meet the needs of intuitive operation and efficient execution in scientific research. By recursive quality control, it can achieve the needs of stable effect in business work without human interference. This method can reasonably eliminate the data with large observation error, the small-scale change data that cannot be resolved by the model resolution and the data that cannot be simulated well, ensure that the probability density function of the actual observation value minus the background field simulation value is close to the Gaussian distribution, improve the accuracy of the analysis field, and facilitate the subsequent radar data assimilation. This method takes full account of the actual needs of data assimilation in numerical weather forecast, and integrates directly into the numerical weather forecast model with the conventional process of numerical weather forecast. It has the advantages of simple operation, fast execution and convenient maintenance.

【技术实现步骤摘要】
一种面向资料同化的天气雷达径向速度质量控制方法
本专利技术涉及气象雷达数据质量控制，特别是涉及一种面向资料同化的天气雷达径向速度质量控制方法。
技术介绍
多普勒天气雷达的灵敏度高，软硬件稳定可靠，具备全天候无人值守自动观测的能力，是监测中小尺度灾害性天气系统的重要设备。多普勒天气雷达可以提供高时空分辨率的径向速度、反射率因子和谱宽三种基数据。其中，径向速度可以准确提供天气系统的风场细节信息，长期以来在灾害性天气预警方面发挥了重要作用，是短时临近天气预报不可或缺的资料。同时，随着我国新一代天气雷达网建设的不断完善和数值天气预报技术的快速发展，多普勒天气雷达径向速度资料同化得到了越来越多的关注，对于数值模式初始场的改进作用也越来越显著。然而，不论是天气雷达径向速度的直接诊断分析，还是后续分析产品的制作，或是天气雷达径向速度资料同化，都需要对径向速度原始观测进行质量控制。目前，国内外学者针对天气雷达径向速度的质量控制研究主要集中于速度退模糊。近几十年来，国内外学者在天气雷达径向速度退模糊方面取得了丰硕的研究成果，然而从资料同化角度来对天气雷达径向速度进行质量控制的研究却相对稀少。目前，主流的资料同化系统都是使用自带的观测资料质量控制方法-阈值检查法，可以把|yo-yb|＞kσ(y是观测资料类型，yo是观测场，yb是背景场，k是常数，σ是标准差)的观测资料剔除掉。然而，在计算观测资料样本的平均值和标准差时，具有较大偏差的离群值会对其产生显著影响，最终会为离群资料的识别带来障碍。将这种质量控制方法应用到雷达径向速度资料中，会出现无法有效识别离群资料或者异常剔除正常资料的问题。诸多研究表明，采用双权重算法可以弱化这种负面影响，该方法主要是利用O-B(O代表yo，B代表yb)的双权重平均值和双权重标准差对离群资料进行有效识别。不过，双权重算法目前仅在GPS资料、地面观测资料、卫星臭氧资料的质量控制及同化中有初步应用，国内外尚无将其应用于雷达径向速度资料质量控制及同化的先例，难点在于如何解决雷达资料和模式背景场资料的时空匹配问题、如何合理协调科研需求和业务需求、如何与常规流程进行有机整合以适应高频次的雷达资料同化。
技术实现思路
鉴于现有技术存在的问题和缺陷，本专利技术充分考虑到数值天气预报中资料同化的实际需求，结合数值天气预报的常规流程，提供一种面向资料同化的天气雷达径向速度质量控制方法。该方法可以直接集成到数值天气预报模式中，具有简便易行、执行快速、维护方便等优点。该方法利用双权重算法并根据其特性，在参数列表中设置了科研和业务的选项，通过自主设定质量控制阈值的方式满足科研工作中操作直观、执行高效的需求，通过递归质量控制的方式实现业务工作中效果稳定、无需人工干扰的需求。本专利技术采取的技术方案是：一种面向资料同化的天气雷达径向速度质量控制方法，其特征在于，所述方法有以下步骤：(1)、利用数值天气预报前处理模块WPS对再分析资料进行处理，生成时空分辨率满足要求的met文件。(2)、将met文件输入数值天气预报初始化模块real.exe，生成背景场wrfinput文件和边界条件wrfbdy文件。(3)、利用雷达径向速度同化的观测算子将背景场从模式空间向观测空间投影，直接输出投影结果，计算距离每个有效的雷达观测值最近的多个投影值的平均值，用于平滑投影后资料中的小尺度扰动，将平均值作为与雷达观测值对应的模式背景值，保证投影后的背景场和雷达资料的时空分辨率一致。(4)、计算每一个雷达观测值与对应的模式背景值的差值。(5)、生成并导出差值数据集。(6)、计算差值数据集中每个资料点的权重函数wi：式(1)中：yo代表雷达观测值，yb代表模式背景值，M是(yo-yb)i的中位数，MAD是|(yo-yb)i-M|的中位数，c是常数。(7)、计算差值数据集中每个资料点的双权重平均值(8)、计算差值数据集中每个资料点的双权重标准差σbw(yo-yb)：(9)、基于公式(2)和公式(3)，计算差值数据集中每个资料点的Zi值：式(4)中，Zqc是识别离群资料的临界值，即质量控制阈值，根据实际需求对质量控制阈值做赋值，随后剔除数据集中Zi≥Zqc的资料。(10)、读取参数列表文件namelist中的option值，option＝1代表启用科研需求流程，option＝2代表启用业务需求流程；当option＝1时，执行以下步骤：(a1)、用户在参数列表文件中根据科研需求自主设定质量控制阈值Zqc；(b1)、删除差值数据集中Zi≥Zqc的资料。当option＝2时，执行以下步骤：(a2)、用户根据业务需求启用质量控制阈值Zqc的默认值；(b2)删除差值数据集中Zi≥Zqc的资料；(c2)计算差值数据集中剩余资料的偏态S和峰度K：式(5)中，S表示偏态值；xi表示差值数据集中第i个数据；表示差值数据集的平均值；n是差值数据集的数据量；式(6)中，K表示峰度值，xi表示差值数据集中第i个数据；表示差值数据集的平均值；n是差值数据集的数据量。(d2)、如果偏态S和峰度K满足|S|≤0.01且|K-3|≤0.1的条件，则执行步骤(11)；如果偏态S和峰度K不满足|S|≤0.01且|K-3|≤0.1的条件，则执行步骤(e2)。(e2)重新赋值质量控制阈值Zqc＝Zqc-E，其中E是质量控制阈值Zqc的递减步长，随后执行步骤(b2)。(11)、将经过质量控制的雷达资料输入数值天气预报资料同化模块。(12)、输出雷达资料同化结果。(13)、根据参数列表文件namelist中所设置的同化参数，判断同化是否结束，如果同化已经结束，则执行后续步骤；如果同化没有结束，则进行雷达资料的循环同化，跳转实现对雷达资料的循环质量控制，直至同化结束，随后执行后续步骤。(14)、将雷达资料同化结果输入数值天气预报主模块。(15)、对数值天气预报主模块生成的模式预报结果进行后处理。在所述步骤(3)中，计算距离每个有效的雷达观测值最近的上、下、左、右、前、后六个投影值的平均值。在所述步骤(10)中的(a2)步骤中、用户根据业务需求启用质量控制阈值Zqc的默认值设置为3。在所述步骤(10)中的(e2)步骤中、重新赋值质量控制阈值Zqc＝Zqc-E，其中质量控制阈值Zqc的递减步长E设置为0.02。本专利技术所产生的有益效果是：(1)本专利技术的天气雷达径向速度资料质量控制方案既可以合理剔除观测误差较大的资料，又可以去除模式分辨率无法分辨的小尺度变化资料或者观测算子难以模拟好的资料，同时还可以保证实际观测值减去背景场模拟值(O-B)的概率密度函数接近高斯分布，提高了分析场的精度，有利于后续的雷达资料同化。(2)本专利技术的天气雷达径向速度资料质量控制方案充分考虑到数值天气预报中资料同化的实际需求，并结合数值天气本文档来自技高网...

【技术保护点】
1.一种面向资料同化的天气雷达径向速度质量控制方法，其特征在于，所述方法有以下步骤：/n(1)、利用数值天气预报前处理模块WPS对再分析资料进行处理，生成时空分辨率满足要求的met文件；/n(2)、将met文件输入数值天气预报初始化模块real.exe，生成背景场wrfinput文件和边界条件wrfbdy文件；/n(3)、利用雷达径向速度同化的观测算子将背景场从模式空间向观测空间投影，直接输出投影结果，计算距离每个有效的雷达观测值最近的多个投影值的平均值，用于平滑投影后资料中的小尺度扰动，将平均值作为与雷达观测值对应的模式背景值，保证投影后的背景场和雷达资料的时空分辨率一致；/n(4)、计算每一个雷达观测值与对应的模式背景值的差值；/n(5)、生成并导出差值数据集；/n(6)、计算差值数据集中每个资料点的权重函数w

【技术特征摘要】
1.一种面向资料同化的天气雷达径向速度质量控制方法，其特征在于，所述方法有以下步骤：
(1)、利用数值天气预报前处理模块WPS对再分析资料进行处理，生成时空分辨率满足要求的met文件；
(2)、将met文件输入数值天气预报初始化模块real.exe，生成背景场wrfinput文件和边界条件wrfbdy文件；
(3)、利用雷达径向速度同化的观测算子将背景场从模式空间向观测空间投影，直接输出投影结果，计算距离每个有效的雷达观测值最近的多个投影值的平均值，用于平滑投影后资料中的小尺度扰动，将平均值作为与雷达观测值对应的模式背景值，保证投影后的背景场和雷达资料的时空分辨率一致；
(4)、计算每一个雷达观测值与对应的模式背景值的差值；
(5)、生成并导出差值数据集；
(6)、计算差值数据集中每个资料点的权重函数wi：



式(1)中：yo代表雷达观测值，yb代表模式背景值，M是(yo-yb)i的中位数，MAD是|(yo-yb)i-M|的中位数，c是常数；
(7)、计算差值数据集中每个资料点的双权重平均值



(8)、计算差值数据集中每个资料点的双权重标准差σbw(yo-yb)：



(9)、基于公式(2)和公式(3)，计算差值数据集中每个资料点的Zi值：



式(4)中，Zqc是识别离群资料的临界值，即质量控制阈值，根据实际需求对质量控制阈值做赋值，随后剔除数据集中Zi≥Zqc的资料；
(10)、读取参数列表文件namelist中的option值，option＝1代表启用科研需求流程，option＝2代表启用业务需求流程；
当option＝1时，执行以下步骤：
(a1)、用户在参数列表文件中根据科研需求自主设定质量控制阈值Zqc；
(b1)、删除差值数据集中Zi≥Zqc的资料；
当option＝2时，执行以下步骤：
(a2)、用户根据业务需求启用质量控制阈值Zqc的默认值；
(b2)删...

【专利技术属性】
技术研发人员：刘寅，赵虹，周红根，唐飞，
申请(专利权)人：江苏省气象探测中心，
类型：发明
国别省市：江苏;32