(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 CN 111366931 A(43)申请公布日 2020.07.03
(21)申请号 202010426361.8(22)申请日 2020.05.19
(71)申请人 范思睿
地址 610072 四川省成都市青羊区光华村
街20号南楼(72)发明人 范思睿 谭学 王维佳 (74)专利代理机构 成都智弘知识产权代理有限
公司 51275
代理人 丁亮 陈春(51)Int.Cl.
G01S 13/95(2006.01)
权利要求书2页 说明书6页 附图4页
(54)发明名称
一种基于云雷达的云降水精细化反演方法(57)摘要
本申请公开了一种基于云雷达的云降水精细化反演方法。这种方法通过提取云雷达数据，首先对云雷达数据进行质量控制，然后利用自主设计的融化层检测算法自动识别云降水融化层高度，然后采用模糊化、规则判断、相态分类检查、最大集成法集成、退模糊等步骤对云降水内部水凝物粒子进行分类识别，精细化反演水凝物粒子相态和分布，最后输出云降水精细化反演结果。实现云降水内部水凝物粒子精细化反演以及融化层高度自动识别。仅依靠云雷达数据就能够反演云降水中水凝物粒子相态及分布，提高了云降水中水凝物粒子识别和反演率，同时在普通的模糊逻辑法上增加了粒子相态检查，提高了水凝物粒子相态识别的准确性。
CN 111366931 ACN 111366931 A
权 利 要 求 书
1/2页
1.一种基于云雷达的云降水精细化反演方法，其特征在于，所述方法包括：提取云雷达数据，云雷达数据包括反射率因子、径向速度、速度谱宽和线性退偏振比；对云雷达数据进行质量控制；提取质控后的云雷达数据，将雷达参量进行距离库单元匹配；融化层检测算法自动识别云降水融化层的顶边界和底边界；
将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度；提取质控后的云雷达数据和融化层高度作为输入参量；采用模糊化、规则判断、相态分类检查、最大集成法集成和退模糊化等步骤对云降水内部水凝物粒子进行分类识别；
将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴和混合相粒子6种类型，并输出云降水精细化反演结果。
2.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述对云雷达数据进行质量控制包括：分别对所述反射率因子、径向速度、速度谱宽和线性退偏振比进行缺测数据补缺、减少噪点和对异常值进行订正。
3.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述将雷达参量进行距离库单元匹配包括：将提取到的反射率因子、径向速度、速度谱宽及线性退极化比等参量进行每个距离库单元匹配，确定每个距离库中参量一一对应。
4.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述融化层检测算法自动识别云降水融化层的顶边界和底边界包括：自动识别融化层的顶边界和底边界，作为基于模糊逻辑算法的云降水精细化反演方法的一个输入参量；融化层检测算法(MLDA)根据径向数据选取回波强度5dBz至30dBz之间，且线性退偏振比在-20dB至-10dB之间的回波点，初步识别为融化层回波点；对每个径向重复该程序，由此产生融化层点群，然后对融化层点群作3*3滤波处理，剔除孤立的融化层回波点，以顶部以下10％融化层点群位置确定为融化层顶边界MLT，以底部以上10％融化层点群位置确定为融化层底边界MLB。
5.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度包括：当云雷达没有出现亮带现象，将邻近时刻的融化层高度作为该时刻融化层高度，若存在融化层高度的邻近时刻与当前时刻间隔超过预设时间，则提取外部环境温度0℃的高度作为辅助输入参量，判别云降水的融化层高度。
6.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述模糊化即提取毫米波云雷达反射率因子、径向速度、速度谱宽、线性退偏振比、融化层高度边界信息(顶边界(MLT)和底边界(MLB))作为输入参量矩阵，通过隶属函数将输入参量矩阵转化为云降水粒子相态矩阵，选取不对称梯形函数对这5个参量进行退模糊化处理。
7.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述规则判断即利用公式对6种粒子进行分类识别，其公式为：
其中，Si为各类粒子集成值，Ti_j代表隶属度函数，i为输出粒子类型，j为输入参量类型，Wi_j表示第i种输出粒子类型，第j种参量对判断结果的权重因子，使Si最大的第i类即为识
2
CN 111366931 A
权 利 要 求 书
2/2页
别的水凝物粒子类型。
8.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述相态分类检查即检查粒子相态与融化层位置的关系，通过融化层的顶边界和底边界对粒子相态分类结果进行检查，检查主要包括融化层底边界以下不允许出现冰晶和雪晶，融化层顶边界以上不允许出现雨滴、毛毛雨。
9.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述最大集成法集成即采用最大值集成方法进行集成，把最大值结果输出。
10.根据权利要求1所述的基于云雷达的云降水精细化反演方法，其特征在于，所述退模糊化即选取最大的那类相态结果作为最终反演输出的水凝物粒子类型。
3
CN 111366931 A
说 明 书
一种基于云雷达的云降水精细化反演方法
1/6页
技术领域
[0001]本申请涉及大气科学技术领域，具体涉及一种基于云雷达的云降水精细化反演方法。
背景技术
[0002]云降水系统有其基本的一般性特征、规律，但同时又具有很强的时空变率，分析云降水内部粒子变化是研究云系结构及其演变特征的基本途径。目前，对云降水中内部结构和粒子相态、大小的观测主要采用飞机观测，飞机观测可以提供云降水内部含水量、粒子浓度、粒子谱及粒子相态等变化。然而，由于飞机观测受到气象条件和空域等因素的制约，无法获得降水过程的连续性观测资料。随着地面探测设备的发展，具有偏振性质的雷达可以测量云降水的反射率因子、径向速度、速度谱宽、线性退偏振比，结合云降水的物理特征，反演云降水粒子相态和大小变化。由于云雷达是地基观测设备，不受时间、空域、天气等因素制约，可以提供连续性观测数据。
[0003]融化层对云降水粒子相态识别和反演具有重要作用，是云降水中水凝物粒子从冰相到液相转变的关键高度。在融化层之上，粒子多为冰相；在融化层之下，粒子多为液态。水凝物粒子经过融化层会导致雷达反射率跃增，表现出亮带现象，通过自动识别亮带可以推算融化层高度和厚度，然后通过一定的规则判定对粒子相态进行识别和反演。[0004]目前，云降水中粒子相态识别和反演算法已有一些研究，主要采用的识别和反演算法有决策树识别方法、统计判决识别法、神经网络或者模糊逻辑识别法等。由于云降水中水凝物粒子特性较为复杂，不同水凝物粒子对应的雷达偏振参量信息并不是绝对排斥，而是存在某种程度的重叠，所以基于“刚性”边界条件和布尔逻辑的决策树识别法不太适合云降水中水凝物粒子分类。而决策树识别法要求所用数据不允许存在测量误差，在实际中难以达到。统计判决识别法是另一种可以解决水凝物分类问题的方法，但是基于不同水凝物的统计模型却很难建立。模糊逻辑识别法相对简便，对外界条件要求更低。发明内容
[0005]本申请为了解决上述技术问题，提出了如下技术方案：
[0006]本申请实施例提供了一种基于云雷达的云降水精细化反演方法，所述方法包括：[0007]提取云雷达数据，云雷达数据包括反射率因子、径向速度、速度谱宽和线性退偏振比；
[0008]对云雷达数据进行质量控制；[0009]提取质控后的云雷达数据，将雷达参量进行距离库单元匹配；
[0010]融化层检测算法(MLDA)自动识别云降水融化层的顶边界和底边界；
[0011]将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度；[0012]提取质控后的云雷达数据和融化层高度作为输入参量；[0013]采用模糊化、规则判断、相态分类检查、最大集成法集成和退模糊化等步骤对
4
CN 111366931 A
说 明 书
2/6页
云降水内部水凝物粒子进行分类识别；
[0014]将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴和混合相粒子6种类型，并输出云降水精细化反演结果。[0015]可选地，所述对云雷达数据进行质量控制包括：分别对所述反射率因子、径向速度、速度谱宽和线性退偏振比进行缺测数据补缺、减少噪点和对异常值进行订正。[0016]可选地，所述将雷达参量进行距离库单元匹配包括：将提取到的反射率因子、径向速度、速度谱宽及线性退极化比等参量进行每个距离库单元匹配，确定每个距离库中参量一一对应。
[0017]可选地，所述融化层检测算法(MLDA)自动识别云降水融化层的顶部和底部包括：自动识别融化层的顶部和底部，作为基于模糊逻辑算法的云降水精细化反演方法的一个输入参量；融化层检测算法根据径向数据选取回波强度5dBz至30dBz之间，且线性退偏振比在-20dB至-10dB之间的回波点，初步识别为融化层回波点；对每个径向重复该程序，产生融化层点群，然后对融化层点群作3*3滤波处理，剔除孤立的融化层回波点，以顶部以下10％融化层点群位置确定为融化层顶部MLT，以底部以上10％融化层点群位置确定为融化层底部MLB。
[0018]可选地，所述将识别出的云降水中融化层高度存储并更新为最新时刻融化层高度包括：当云雷达没有出现亮带现象，将邻近时刻的融化层高度作为该时刻融化层高度，若存在融化层高度的邻近时刻与当前时刻间隔超过预设时间(本实施例中预设时间为12小时)，则提取外部环境温度0℃高度作为辅助输入参量，判别融化层高度。[0019]可选地，所述模糊化即提取毫米波云雷达反射率因子、径向速度、速度谱宽、线性退偏振比、融化层高度边界信息(顶边界(MLT)和底边界(MLB))作为输入参量矩阵，通过隶属函数将输入参量矩阵转化为云降水粒子相态矩阵，选取不对称梯形函数对这5个参量进行退模糊化处理。
[0020]梯形函数T(x)表达式如下：
[0021]
Ti_j表示输入参量对各相态的隶属程度，越接近1表示相态为真实相态的可能性越
高，其中i代表第i个输出粒子相态类型，分别对应冰晶、雪晶、云滴、毛毛雨、雨滴、混合相粒子等，j代表第j个输入参量，分别对应反射率因子、径向速度、速度谱宽、线性退偏振比，X1、X2、X3、X4是隶属度函数的系数。[0023]可选地，所述规则判断即利用公式对6种粒子进行分类识别，其公式为：
[0024][0025]
[0022]
其中，Si为各类粒子集成值，Ti_j代表隶属度函数，i为输出粒子类型，j为输入参量类型，Wi_j表示第i种输出粒子类型，第j种参量对判断结果的权重因子，使Si最大的第i类即
5
CN 111366931 A
说 明 书
3/6页
为识别的水凝物粒子类型。[0026]可选地，所述相态分类检查即检查粒子相态类别与融化层位置的关系，通过融化层的边界信息对粒子相态分类结果进行检查，检查主要包括融化层底边界以下不允许出现冰晶和雪晶，融化层顶边界以上不允许出现雨滴、毛毛雨，具体见相态分类检查类别集。[0027](1)当H>MLT时，可以取的粒子类别包括：冰晶、雪晶；[0028](2)当MLBMLT时代表粒子处于低于0℃的环境中，一般为冰相，粒子相态可能会冰晶和雪晶；当MLB[0032]可选地，所述退模糊化即选取最大的那类相态结果作为最终反演输出的水凝物粒子类型。[0033]本申请的有益效果是：实现云降水内部水凝物精细化反演以及融化层高度自动识别，解决了现有观测手段无法获得连续性高精度云降水内部水凝物粒子分布的问题。可以根据时间浮动判别云降水融化层位置变化。仅依靠云雷达数据就能够反演云降水中水凝物粒子类型、相态变化以及分布位置，提高了云降水中水凝物粒子识别和反演率，同时在普通的模糊逻辑识别法上增加了粒子相态检查，提高了云降水中水凝物粒子相态识别准确性。
附图说明
[0034]图1为本申请实施例提供的一种基于云雷达的云降水精细化反演方法的流程示意图；
[0035]图2为本申请实施例提供的梯形隶属函数示意图；
[0036]图3为本申请实施例提供的基于云雷达的云降水内部水凝物粒子相态分类的示意图；
[0037]图4为本申请实施例提供的云雷达回波图中反射率因子(a)、径向速度(b)、速度谱宽(c)和线性退极化比(d)示意图；
[0038]图5为本申请实施例提供的利用融化层检测算法(MLDA)自动识别融化层边界信息的散点示意图；
[0039]图6为本申请实施例提供的探空资料温度对数压力示意图；
[0040]图7为本申请实施例提供的基于云雷达的云降水精细化反演结果。具体实施方式
[0041]下面结合附图与具体实施方式对本方案进行阐述。
6
CN 111366931 A[0042]
说 明 书
4/6页
图1为本申请实施例提供的一种基于云雷达的云降水精细化反演方法的流程示意
图，参见图1，所述方法包括：[0043]S101，提取云雷达数据，所述云雷达数据包括反射率因子、径向速度、速度谱宽和线性退偏振比。[0044]S102，对所述云雷达数据进行质量控制。[0045]雷达数据一般存在不同程度的质量问题，主要做以下几方面的数据质量控制工作：分别对所述反射率因子、径向速度、速度谱宽和线性退偏振比进行缺测数据补缺、减少噪点和对异常值进行订正。[0046]S103，提取质量控制后的云雷达数据，将雷达参量进行距离库单元匹配。[0047]将提取到的反射率因子、径向速度、速度谱宽及线性退极化比等参量进行每个距离库单元匹配，确定每个距离库中参量一一对应。[0048]S104，融化层检测算法(MLDA)自动识别云降水的融化层顶部和底部。
[0049]本申请自主开发的融化层检测算法(MLDA)自动识别融化层的顶部和底部，作为云降水精细化反演方法的一个输入参量。融化层检测算法(MLDA)根据径向数据选取回波强度5dBz至30dBz之间，且线性退偏振比在-20dB至-10dB之间的回波点，初步识别为融化层回波点。对每个径向重复该程序，产生融化层点群，然后对融化层点群作3*3滤波处理，剔除孤立的融化层回波点，以顶部以下10％融化层点群位置确定为融化层的顶部，以底部以上10％融化层点群位置确定为融化层底部。[0050]S105，将识别出的云降水中融化层高度信息存储并更新为最新时刻融化层高度。[0051]将上一步识别出云降水融化层高度存储并更新为最新时刻融化层高度，当云雷达没有出现亮带现象，就无法识别出云降水融化层，就将邻近时刻的融化层高度作为当下时刻的融化层高度，若有融化层高度的邻近时刻与当下时刻间隔超过12小时，就提取外部环境温度0℃高度(比如从探空资料、数值预报产品提取0℃高度)作为辅助输入参量，判别云降水融化层高度。[0052]S106，提取质控后的云雷达数据和融化层高度作为输入参量。[0053]S107，采用模糊化、规则判断、相态分类检查、最大集成法集成和退模糊化等步骤对云降水内部水凝物粒子进行分类识别。[0054]在识别云降水融化层高度的前提下，利用反射率因子、径向速度、速度谱宽、线性退偏振比、融化层高度边界(顶边界(MLT)和底边界(MLB))作为输入参量，采用模糊化、规则判断、相态分类检查、最大集成法集成、退模糊等步骤对云降水内部粒子相态进行分类识别，将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴、混合相粒子6种类型。[0055]进一步解释模糊化即提取云雷达反射率因子、径向速度、速度谱宽、线性退偏振比、融化层高度边界(顶边界(MLT)和底边界(MLB))作为输入参量矩阵，通过隶属函数将输入参量矩阵转化为云降水粒子相态矩阵，本实施例选取不对称梯形函数对这5个参量进行退模糊化处理。其梯形函数T(x)示意图如图2所示，表达式如下：
7
CN 111366931 A
说 明 书
5/6页
[0056]
其中，Ti_j表示输入参量对各相态的隶属程度，越接近1表示相态为真实相态的可
能性越高，其中i代表第i个输出粒子相态类型，分别对应冰晶、雪晶、云滴、毛毛雨、雨滴、混合相粒子等，j代表第j个输入参量，分别对应反射率因子、径向速度、速度谱宽、线性退偏振比，X1、X2、X3、X4是隶属度函数的系数。[0058]进一步解释规则判断，规则判断即利用下式对6种粒子进行分类识别：
[0059]
[0057]
其中，Si为各类粒子集成值，Tij代表隶属度函数，i为输出粒子类型，j为输入参量类型，Wi_j表示第i种输出粒子类型,第j种参量对判断结果的权重因子，使Si最大的第i类即为识别的水凝物粒子类型。
[0061]进一步解释相态分类检查，相态分类检查即检查粒子相态类别与融化层位置的关系，通过云降水融化层的顶边界和底边界对粒子相态分类结果进行检查，相态分类检查主要包括融化层底部以下不允许出现冰晶和雪晶，融化层顶部以上不允许出现雨滴和毛毛雨，具体见相态分类检查类别集。[0062](1)当H>MLT时，可以取的粒子类别包括：冰晶、雪晶；[0063](2)当MLBMLT时代表粒子处于低于0℃的环境中，一般为冰相，粒子相态可能会冰晶和雪晶；当MLB[0066]进一步解释最大集成法集成。最大集成法集成即采用最大值集成方法进行集成，即最大集成法仅把最大值结果输出。[0067]Smax＝Max(Si)[0068]进一步解释退模糊化。退模糊化即选取最大的那类相态结果作为最终反演输出的水凝物粒子类型。[0069]S108，将云降水中水凝物粒子分为冰晶、雪晶、云滴、毛毛雨、雨滴和混合相粒子6种类型，并输出云降水精细化反演结果。
[0070]下面选用具体的实施例对上述方法进行说明：[0071]选用某市云雷达垂直观测资料，具体步骤如下：[0072]提取云雷达资料，如图4所示，分别为某时间的反射率因子(a)、径向速度(b)、速度谱宽(c)、线性退偏振比(d)四个参量。
8
[0060]
CN 111366931 A[0073]
说 明 书
6/6页
图5为使用本申请自主设计的融化层检测算法(MLDA)对数据中的融化层顶边界
(MLT)和融化层底边界(MLB)进行动态识别后的结果，并与实际环境温度值(探空资料)图6进行对比，本申请自主设计识别出的云降水融化层高度和实际环境温度0℃层高度基本一致。[0074]提取到融化层高度顶边界(MLT)和融化层底边界(MLB)后，利用模糊逻辑算法对云降水粒子相态进行识别，并输出识别结果为云降水精细化反演产品，如图7所示。[0075]由上述实施例可知，本实施例提供的一种基于云雷达的云降水精细化反演方法，实现云降水内部水凝物精细化反演以及融化层高度自动识别，解决了现有观测手段无法获得连续性高精度云降水中水凝物粒子分布的问题。可以根据时间浮动判别云降水融化层位置变化。仅依靠云雷达数据就能够反演云降水中水凝物粒子类型、相态变化以及分布位置，提高了云降水粒子识别和反演率，同时在普通的模糊逻辑识别法上增加了粒子相态检查，提高了云降水内部水凝物粒子相态识别准确性。[0076]当然，上述说明也并不仅限于上述举例，本申请未经描述的技术特征可以通过或采用现有技术实现，在此不再赘述；以上实施例及附图仅用于说明本申请的技术方案并非是对本申请的，如来替代，本申请仅结合并参照优选的实施方式进行了详细说明，本领域的普通技术人员应当理解，本技术领域的普通技术人员在本申请的实质范围内所做出的变化、改型、添加或替换都不脱离本申请的宗旨，也应属于本申请的权利要求保护范围。
9
CN 111366931 A
说 明 书 附 图
1/4页
图1
10
CN 111366931 A
说 明 书 附 图
2/4页
图2
图3
11
CN 111366931 A
说 明 书 附 图
3/4页
图4
图5
12
CN 111366931 A
说 明 书 附 图
4/4页
图6
图7
13
