一种基于云雷达的云降水精细化反演方法技术

本申请公开了一种基于云雷达的云降水精细化反演方法。这种方法通过提取云雷达数据，首先对云雷达数据进行质量控制，然后利用自主设计的融化层检测算法自动识别云降水融化层高度，然后采用模糊化、规则判断、相态分类限制检查、最大集成法集成、退模糊等步骤对云降水内部水凝物粒子进行分类识别，精细化反演水凝物粒子相态和分布，最后输出云降水精细化反演结果。实现云降水内部水凝物粒子精细化反演以及融化层高度自动识别。仅依靠云雷达数据就能够反演云降水中水凝物粒子相态及分布，提高了云降水中水凝物粒子识别和反演率，同时在普通的模糊逻辑法上增加了粒子相态限制检查，提高了水凝物粒子相态识别的准确性。

A refined inversion method of cloud precipitation based on Cloud Radar

【技术实现步骤摘要】
一种基于云雷达的云降水精细化反演方法
本申请涉及大气科学
，具体涉及一种基于云雷达的云降水精细化反演方法。
技术介绍
云降水系统有其基本的一般性特征、规律，但同时又具有很强的时空变率，分析云降水内部粒子变化是研究云系结构及其演变特征的基本途径。目前，对云降水中内部结构和粒子相态、大小的观测主要采用飞机观测，飞机观测可以提供云降水内部含水量、粒子浓度、粒子谱及粒子相态等变化。然而，由于飞机观测受到气象条件和空域等因素的制约，无法获得降水过程的连续性观测资料。随着地面探测设备的发展，具有偏振性质的雷达可以测量云降水的反射率因子、径向速度、速度谱宽、线性退偏振比，结合云降水的物理特征，反演云降水粒子相态和大小变化。由于云雷达是地基观测设备，不受时间、空域、天气等因素制约，可以提供连续性观测数据。融化层对云降水粒子相态识别和反演具有重要作用，是云降水中水凝物粒子从冰相到液相转变的关键高度。在融化层之上，粒子多为冰相；在融化层之下，粒子多为液态。水凝物粒子经过融化层会导致雷达反射率跃增，表现出亮带现象，通过自动识别亮带可以推算融化层高度和厚度，然后通过一定的规则判定对粒子相态进行识别和反演。目前，云降水中粒子相态识别和反演算法已有一些研究，主要采用的识别和反演算法有决策树识别方法、统计判决识别法、神经网络或者模糊逻辑识别法等。由于云降水中水凝物粒子特性较为复杂，不同水凝物粒子对应的雷达偏振参量信息并不是绝对排斥，而是存在某种程度的重叠，所以基于“刚性”边界条件和布尔逻辑的决策树识别法不太适合云降水中水凝物粒子分类。而决策树识别法要求所用数据不允许存在测量误差，在实际中难以达到。统计判决识别法是另一种可以解决水凝物分类问题的方法，但是基于不同水凝物的统计模型却很难建立。模糊逻辑识别法相对简便，对外界条件要求更低。
技术实现思路
本申请为了解决上述技术问题，提出了如下技术方案：本申请实施例提供了一种基于云雷达的云降水精细化反演方法，所述方法包括：提取云雷达数据，云雷达数据包括反射率因子、径向速度、速度谱宽和线性退偏振比；对云雷达数据进行质量控制；提取质控后的云雷达数据，将雷达参量进行距离库单元匹配；融化层检测算法(MLDA)自动识别云降水融化层的顶边界和底边界；将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度；提取质控后的云雷达数据和融化层高度作为输入参量；采用模糊化、规则判断、相态分类限制检查、最大集成法集成和退模糊化等步骤对云降水内部水凝物粒子进行分类识别；将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴和混合相粒子6种类型，并输出云降水精细化反演结果。可选地，所述对云雷达数据进行质量控制包括：分别对所述反射率因子、径向速度、速度谱宽和线性退偏振比进行缺测数据补缺、减少噪点和对异常值进行订正。可选地，所述将雷达参量进行距离库单元匹配包括：将提取到的反射率因子、径向速度、速度谱宽及线性退极化比等参量进行每个距离库单元匹配，确定每个距离库中参量一一对应。可选地，所述融化层检测算法(MLDA)自动识别云降水融化层的顶部和底部包括：自动识别融化层的顶部和底部，作为基于模糊逻辑算法的云降水精细化反演方法的一个输入参量；融化层检测算法根据径向数据选取回波强度5dBz至30dBz之间，且线性退偏振比在-20dB至-10dB之间的回波点，初步识别为融化层回波点；对每个径向重复该程序，产生融化层点群，然后对融化层点群作3*3滤波处理，剔除孤立的融化层回波点，以顶部以下10％融化层点群位置确定为融化层顶部MLT，以底部以上10％融化层点群位置确定为融化层底部MLB。可选地，所述将识别出的云降水中融化层高度存储并更新为最新时刻融化层高度包括：当云雷达没有出现亮带现象，将邻近时刻的融化层高度作为该时刻融化层高度，若存在融化层高度的邻近时刻与当前时刻间隔超过预设时间(本实施例中预设时间为12小时)，则提取外部环境温度0℃高度作为辅助输入参量，判别融化层高度。可选地，所述模糊化即提取毫米波云雷达反射率因子、径向速度、速度谱宽、线性退偏振比、融化层高度边界信息(顶边界(MLT)和底边界(MLB))作为输入参量矩阵，通过隶属函数将输入参量矩阵转化为云降水粒子相态矩阵，选取不对称梯形函数对这5个参量进行退模糊化处理。梯形函数T(x)表达式如下：Ti_j表示输入参量对各相态的隶属程度，越接近1表示相态为真实相态的可能性越高，其中i代表第i个输出粒子相态类型，分别对应冰晶、雪晶、云滴、毛毛雨、雨滴、混合相粒子等，j代表第j个输入参量，分别对应反射率因子、径向速度、速度谱宽、线性退偏振比，X1、X2、X3、X4是隶属度函数的系数。可选地，所述规则判断即利用公式对6种粒子进行分类识别，其公式为：其中，Si为各类粒子集成值，Ti_j代表隶属度函数，i为输出粒子类型，j为输入参量类型，Wi_j表示第i种输出粒子类型，第j种参量对判断结果的权重因子，使Si最大的第i类即为识别的水凝物粒子类型。可选地，所述相态分类限制检查即检查粒子相态类别与融化层位置的关系，通过融化层的边界信息对粒子相态分类结果进行限制检查，限制检查主要包括融化层底边界以下不允许出现冰晶和雪晶，融化层顶边界以上不允许出现雨滴、毛毛雨，具体见相态分类限制检查类别集。(1)当H>MLT时，可以取的粒子类别包括：冰晶、雪晶；(2)当MLB<H<MLT时，可以取的粒子类别包括：混合相粒子；(3)当H<MLB时，可以取的粒子类别包括：雨滴、云滴、毛毛雨。其中，MLT为融化层顶高度，MLB为融化层底高度，H为云降水粒子高度，用识别出的融化层顶高度MLT来近似代表环境温度0℃，当H>MLT时代表粒子处于低于0℃的环境中，一般为冰相，粒子相态可能会冰晶和雪晶；当MLB<H<MLT时代表粒子处于0℃附近，一般既有冰相也有液态，统称为混合相粒子；当H<MLB时代表粒子处于高于0℃的环境中，一般为液态，粒子可能为雨滴、云滴、毛毛雨。可选地，所述最大集成法集成即采用最大值集成方法进行集成，把最大值结果输出。可选地，所述退模糊化即选取最大的那类相态结果作为最终反演输出的水凝物粒子类型。本申请的有益效果是：实现云降水内部水凝物精细化反演以及融化层高度自动识别，解决了现有观测手段无法获得连续性高精度云降水内部水凝物粒子分布的问题。可以根据时间浮动判别云降水融化层位置变化。仅依靠云雷达数据就能够反演云降水中水凝物粒子类型、相态变化以及分布位置，提高了云降水中水凝物粒子识别和反演率，同时在普通的模糊逻辑识别法上增加了粒子相态限制检查，提高了云降水中水凝物粒子相态识别准确性。附图说明图1为本申请实施例提供的一种基于云雷达的云降水精细化反演方法的流程示意图；图2为本申请实施例提供的梯形隶属函数示意图；图3为本本文档来自技高网...

【技术保护点】
1.一种基于云雷达的云降水精细化反演方法，其特征在于，所述方法包括：/n提取云雷达数据，云雷达数据包括反射率因子、径向速度、速度谱宽和线性退偏振比；/n对云雷达数据进行质量控制；/n提取质控后的云雷达数据，将雷达参量进行距离库单元匹配；/n融化层检测算法自动识别云降水融化层的顶边界和底边界；/n将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度；/n提取质控后的云雷达数据和融化层高度作为输入参量；/n采用模糊化、规则判断、相态分类限制检查、最大集成法集成和退模糊化等步骤对云降水内部水凝物粒子进行分类识别；/n将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴和混合相粒子6种类型，并输出云降水精细化反演结果。/n

【技术特征摘要】
1.一种基于云雷达的云降水精细化反演方法，其特征在于，所述方法包括：
提取云雷达数据，云雷达数据包括反射率因子、径向速度、速度谱宽和线性退偏振比；
对云雷达数据进行质量控制；
提取质控后的云雷达数据，将雷达参量进行距离库单元匹配；
融化层检测算法自动识别云降水融化层的顶边界和底边界；
将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度；
提取质控后的云雷达数据和融化层高度作为输入参量；
采用模糊化、规则判断、相态分类限制检查、最大集成法集成和退模糊化等步骤对云降水内部水凝物粒子进行分类识别；
将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴和混合相粒子6种类型，并输出云降水精细化反演结果。


2.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述对云雷达数据进行质量控制包括：分别对所述反射率因子、径向速度、速度谱宽和线性退偏振比进行缺测数据补缺、减少噪点和对异常值进行订正。


3.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述将雷达参量进行距离库单元匹配包括：将提取到的反射率因子、径向速度、速度谱宽及线性退极化比等参量进行每个距离库单元匹配，确定每个距离库中参量一一对应。


4.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述融化层检测算法自动识别云降水融化层的顶边界和底边界包括：自动识别融化层的顶边界和底边界，作为基于模糊逻辑算法的云降水精细化反演方法的一个输入参量；融化层检测算法(MLDA)根据径向数据选取回波强度5dBz至30dBz之间，且线性退偏振比在-20dB至-10dB之间的回波点，初步识别为融化层回波点；对每个径向重复该程序，由此产生融化层点群，然后对融化层点群作3*3滤波处理，剔除孤立的融化层回波点，以顶部以下10％融化层点群位置确定为融化层顶边界MLT，以底部以上10％融化层点群位置确定为...

【专利技术属性】
技术研发人员：范思睿，谭学，王维佳，
申请(专利权)人：范思睿，
类型：发明
国别省市：四川;51