基于模糊逻辑的新一代天气雷达径向干扰回波识别算法

奥鹏网院作业

基于模糊逻辑的新一代天气雷达径向干扰回波识别算法*
文 浩1 张乐坚1 梁海河1 张 扬2

1.中国气象局气象探测中心，北京，100081

2.中国气象科学研究院灾害天气国家重点实验室，北京，100081

摘 要 在雷达信号处理器异常、信号受到外部电磁干扰时，雷达会观测到沿径向分布的条幅状干扰回波，这类非气象回波须在雷达估测降水等定量应用之前消除。为适用体扫数据及正在业务推行的单径向流数据，利用2017 年7—8 月收集的雷达单站PPI 数据，针对径向干扰回波，提出一种基于模糊逻辑的识别方法，该方法只提取径向上的特征参量，建立隶属函数和判定阈值，进行识别。通过对算法的评估和业务应用，结果表明：该算法对径向干扰回波和降水回波的识别和消除效果较好，与原业务算法相比，CSI 评分提高0.232，能明显提高径向干扰回波的准确识别，略微降低对降水回波的误消除；算法在业务上能够稳定运行，对于定量化应用（降水估测等）具有正贡献。

关键词 天气雷达，径向干扰回波，模糊逻辑

1 引言
中国自1999 年开展新一代天气雷达（CINRAD）布网以来，已完成200 余部多普勒天气雷达（S、C 波段）的建设和运行，在台风（高拴柱等，2007）、暴雨（孙继松等，2015）和强对流（郑媛媛等，2004）等灾害性天气的监测、预警和预报中发挥了不可替代的作用。同时，天气雷达资料已初步应用于数值预报业务（Bao，et al，2017），能明显改善6 h 降水强度和落区预报，有效提高模式分辨率（盛春岩等，2006）。为了提升天气雷达的应用效益，满足中国经济社会发展对强对流天气监测服务提出的更高要求，提高天气雷达数据质量至关重要。《2017 年全国天气雷达数据质量评估报告》表明，径向干扰回波是出现频次较高的非气象回波，占数据质量问题总数的3.25%以上，对数据质量影响严重。

径向干扰回波是一种沿径向方向分布的条幅状回波，能够延展到最大探测距离（杜言霞等，2017a）。这类回波会污染天气雷达在受影响方位扇区的气象信息，将其误认为大气观测数据（Cho，2017），可以由外部远距离单频点电磁波干扰（周红根等，2008）和内部数据处理异常引起。

过去十年，气象学界特别关注无线设备的射频干扰对天气雷达的影响（Saltikoff，et al，2016）。在信号处理端，基于I/Q（In-phase and Quadrature）数据，进行了大量研究工作，如使用单脉冲异常检测滤波器检测（FAA，1995），并用相邻脉冲的插值代替标记脉冲（Vaisala，2014）；在多普勒频率中，干扰表现为白噪声（Joe，et al，2005），可通过检测频谱噪声的异常增量并去除多余的功率（Cho，2010，2011）。中国I/Q 数据基本未保存，利用基数据的处理方法有：孤立回波消除算法和径向干扰识别算法（杨川，2012），区域滤波和检查回波径向、方位上连续性的识别方法（江源，2013），边缘识别法和径向识别法（郭春辉等，2014）等。杜言霞等（2017a，2017b）对数据进行退距离订正和衰减订正后，利用特征阈值判定识别径向干扰回波，并比较了滤波法、插值法及功率法的效果。

基于基数据的识别方法大多只做了个例研究，且一般针对较窄的径向干扰回波做分析，未进行定量分析以及开展业务化工作。根据“新一代天气雷达建设业务软件系统（Radar Operational Software Engineering，ROSE）”中径向干扰回波识别算法的评估结果，ROSE 算法能有效识别窄径向（方位角连续不超过5°）的干扰回波，对于大面积径向（方位角连续超过5°）的干扰回波无法识别（文浩等，2016）。中国气象局气象探测中心现有的业务化径向干扰回波识别算法（AO）对于此类回波识别效果也较差。同时，中国气象局已开展天气雷达单径向流数据传输工作，为了形成快速质量控制及拼图业务，满足业务应用需求，提升识别效果，文中仅提取回波径向上的特征参数，建立基于模糊逻辑的识别算法（AFL），并使用数据统计和典型个例对比分析AO、AFL 两种算法的识别效果。

2 数据与特征分析
2.1 特征分析
径向干扰由雷达信号处理器异常和信号受到外部电磁干扰引起。一般情况，所有仰角都可能出现径向干扰回波，反射率因子值一般为5—60 dBz（图1a—c），特别地，当位于雷达频率范围的太阳辐射波被雷达接收到时（NOAA，2006；Tang，et al，2014），雷达某个仰角上的1—2 个径向上会出现较弱的径向干扰回波，反射率因子值一般为5—20 dBz（图1d）。

这些情况下，雷达接收的信号是由正常观测返回的信号和灌入信号叠加而成（江源，2013）。根据反射率因子Z（单位：dBz）、回波功率与雷达最小可测功率之比B（单位：dB）的关系（俞小鼎等，2006），

pagenumber_ebook=119,pagenumber_book=117
式中，R（单位：km）为某观测点到雷达站的距离，A 是与雷达性能有关的常数。

即使接收的是稳定信号，B 转换成Z 时经过了20lgR 的距离订正（同一部雷达常数A 相同），使得Z 的大小随距离的增大有一定的增强。因此，去除距离订正后，径向干扰回波的B 在径向上由始至终差异较小，而由于降水粒子存在的位置、强度不同，降水回波的B 在径向上随距离应没有明显规律。

分析图1a 的355°方位角径向干扰回波和214°方位角降水回波在径向上Z、B 值随距离的变化，显示结果与理论分析一致（图2），即去除距离订正后，径向干扰回波的B 在整个径向上差异不大（前50 个距离库混合有超折射回波波动较大），降水回波的B 在径向上随距离没有明显规律。

pagenumber_ebook=120,pagenumber_book=118
图1 径向干扰回波示例（a、b.2017 年8 月7 日13 时48 分濮阳SA 雷达PPI，c.2017 年8 月18 日21 时28 分湖州SA 雷达PPI，d.2017 年8 月31 日21 时33 分台州SA 雷达PPI）
Fig.1 Examples of radial interference echo（a，b.PPI of Puyang SA radar at 13:48 UTC 7 August 2017，c.PPI of Huzhou SA radar at 21:28 UTC 18 August 2017，d.PPI of Taizhou SA radar at 21:33 UTC 31 August 2017）

由于干扰源类型的不同，径向干扰回波形态不同。在数据收集过程中发现，按照径向干扰回波局部水平纹理的变化可以将其分为两类，第1 类径向干扰回波局部水平纹理平滑，由发射连续波的设备造成（图1a）；第2 类径向干扰回波局部水平纹理粗糙，由发射脉冲波的设备造成（图1c）。且这两类回波一般能延展到雷达的最大探测距离。

为有效应用天气雷达流数据，对其快速质量控制，形成高质量、高时间分辨率的产品，采用回波径向特征进行分析，即文中不采用相邻仰角、相邻方位角的特征参数对径向干扰回波进行识别。

2.2 数据收集
pagenumber_ebook=121,pagenumber_book=119
图2 图1a 中355°和214°径向Z（a）与B（b）随距离变化的对比
Fig.2 Comparison of Z（a）and B（b）with distance in 355°，214° radial direction of Fig.1a

算法的建立需要一个含径向干扰回波、层云降水回波和对流云降水回波的数据库，以统计分析3 种回波的不同特征参数的差异。人为判定体扫资料径向干扰回波可以从回波的形态、径向和方位的延展性进行判断；对流云降水回波为不规则的块状结构，结构紧密，反射率因子中心一般在35 dBz 以上；层云降水回波的范围较大，连绵成片，强度较均匀，反射率因子通常为15—35 dBz。

根据中国气象局气象探测中心《每日综合气象观测快报》的人工记录，收集2017 年7—8 月42 个站共475 个径向干扰PPI 数据，420 个层云降水回波PPI 数据，420 个对流云降水回波PPI 数据。按照文件排序，将奇数文件作为统计数据，偶数文件作为验证数据，建立包含统计特征参数概率分布的统计数据库和验证算法识别率的验证数据库（表1），数据库以一个方位-距离库为一个单位。其中，RI-I、RIII、CC 和SC 分别表示第1 类径向干扰回波，第2 类径向干扰回波，对流云降水回波和层云降水回波。

表1 数据库
Table 1 Data set

pagenumber_ebook=121,pagenumber_book=119
3 方法建立
AO 基于功率法开展，对于RREF（单位：%）超过80%以上的径向，计算方位角上某一距离库未经距离订正的Bj（单位：dB）（由于同部雷达，雷达常数A 为定值，不影响B 的变化趋势，计算时忽略不计），再计算该距离库所在径向上最后10%的有效距离库的B 的平均（单位：dB），最后计算DB（单位：dB），若DB小于事先设定的阈值，则判定此距离库为径向干扰回波。RREF、Bj、pagenumber_ebook=121,pagenumber_book=119、DB的定义如式（2）—（6）

pagenumber_ebook=121,pagenumber_book=119
式中，Zj（单位：dBz）为方位角上某距离库的反射率因子，Val（单位：dBz）为有效探测值，Rj（单位：km）为方位角上某距离库与雷达的距离，NR为反射率因子的距离库数。

AO 利用阈值和单个特征参量对回波进行识别，但研究（Berenguer，et al，2006；刘黎平等，2007）表明，多个参数的有效结合能更好区分降水和非降水回波，因此，AFL 在AO 基础上利用模糊逻辑方法开展，提取能够反映径向干扰回波特征的物理参数，根据参数的概率分布建立相应的隶属函数以及径向干扰回波0—1 的取值判据，然后对判据值加权累加，当某点的判据值超过事先给定的阈值时，该点就被判定为径向干扰回波并剔除。

3.1 特征参数
根据2.1 节的特征分析，从反射率因子中提取4 个能够反映径向干扰回波和降水回波差异的特征参数：DB，表示当前径向前后距离库回波功率的一致性；RREF，表示反射率因子在当前径向上的延展性；TDBZ（单位：dB2），表示径向上局部反射率因子回波纹理的一致性（Kessinger，et al，2003）；SPIN，表示局部范围内相邻反射率因子符号变化的情况（Steiner，et al，2002；Rennie，et al，2015）。其 中，TDBZ和SPIN的定义如式（7）—（9）

pagenumber_ebook=122,pagenumber_book=120
式中，Zthresh为距离库之间反射率因子变化的阈值。

3.2 概率分布与隶属函数
pagenumber_ebook=122,pagenumber_book=120
图3 径向干扰回波和降水回波的各特征参数（a—d）概率分布
Fig.3 Probability distribution of characteristic parameters（a-d）for radial interference echo and precipitation echo

基于统计数据库，统计径向干扰回波、对流云降水回波和层云降水回波的特征参数的概率分布（图3）。分析径向干扰回波和降水回波的各特征参数差异：RI-I 和RI-II 的RREF较大，主要分布在85%—100%，该特征值与CC、SC 雷达的差异较大；RI-I 和RIII 的DB集中于0 值附近，而CC、SC 雷达的DB分布较为广泛；RI-I 的TDBZ较CC、SC 雷达的变化较小，而RI-II 的TDBZ较CC、SC 雷达的变化较大；RI-I 的SPIN较小，主要集中在0—4，RI-II 的SPIN在大值区分布较多，CC 和SC 雷达的SPIN基本分布于0—6。尽管这些特征物理量有一定的能力将径向干扰回波与降水回波相区分，但特征分布在很多区域是重叠的，所以利用单一的物理量来区分是不够有效的，需要利用多个特征参量联合识别。

概率分布函数是量化某个距离库反射率因子值是某种类型（径向干扰或降水）的置信度。因此，隶属函数的形状应该类似于概率分布。通常将隶属函数定义为简单的曲线，最常见的是三角形、梯形、分段线性以及高斯函数（Mendel，1995）。根据分布结果以及满足业务运行效率，文中使用分段线性隶属函数来表现概率分布的形状。由此，根据特征参数RREF、DB、TDBZ和SPIN的概率分布，确定各参数的隶属函数（图4）。为了在业务中实现自动化应用，设定TDBZ＞3 dB2 的距离库使用RI-II 的隶属函数，否则使用RI-I 的隶属函数。

pagenumber_ebook=123,pagenumber_book=121
图4 径向干扰回波的各特征参数（a—d）隶属函数
Fig.4 Membership functions of characteristic parameters（a-d）for radial interference echo

3.3 特征参数权重
权重的选取是利用对统计数据库识别结果的临界成功指数（CSI）评分确定的，定义如式（10）和（11）

pagenumber_ebook=123,pagenumber_book=121
式中，NALL表示参与统计的总样本数，NHIT表示准确识别的样本总数，POD 表示识别准确率，PODRI和PODPC分别表示径向干扰回波和降水回波的识别准确率。

在权重选取的过程中，使用以下先验约束条件来减少分析组合的数量：（1）特征参数RREF、DB权重占比大于0；（2）特征参数TDBZ和SPIN在RI-I 识别时权重占比分别为1 和0，在RI-II 识别时权重占比各为1，这是根据径向干扰和降水回波在参数TDBZ和SPIN的概率分布重叠区域较多决定的；（3）每个特征参数的权重不能超过权重总和的40%，因此，各参数权重的可能组合如表2 所示。

表2 RI-I 和RI-II 特征参数权重占比组合及CSI 评分
Table 2 RI-I and RI-II characteristic parameter weight proportions and CSI scores

pagenumber_ebook=124,pagenumber_book=122
根据CSI 评分结果，序号2 和序号8 的CSI 评分最高，但序号2 相较于序号8 的权重占比会造成更多降水回波误识别。在尽可能保留降水且提升径向干扰回波识别准确率的原则上，文中选择序号8，即RI-I 的特征参数RREF、DB、TDBZ和SPIN的权重占比分别为2、2、1、0；RI-II 的特征参数RREF、DB、TDBZ和SPIN的权重占比分别为1、2、1、1。

4 效果评估
为了验证AFL 对于径向干扰回波和降水回波的识别效果，开展以下工作：（1）使用验证数据库对总体识别准确率进行评估，（2）使用个例资料对识别效果进行补充说明，（3）说明算法在业务上的应用效果。

4.1 总体评估
识别效果的总体评估指标包括，算法对径向干扰回波或降水回波的识别准确率（式（10））及整体CSI 评分（式（11）），其中PODRI-I、PODRI-II、PODCC和PODSC分别表示RI-I、RI-II、CC 和SC 的识别准确率。

根据表3 的统计结果，AFL 相对于原业务算法AO，对于RI-I 和RI-II 的识别准确率都有很大提升，对于CC 的识别率略有提升（0.052），对于SC 的识别率没有改变，整体CSI 评分提高0.232。由于AFL 的特征参数的联合应用会放大3.1 节中选取的径向干扰的特征，导致类似径向干扰回波特征的降水回波的误识增多，同时也减少与径向干扰回波相区别的降水回波的误识，因此在降水回波识别准确率较高的情况下，虽有提升但不明显，4.2 节中将会结合个例做进一步分析。

表3 AO 和AFL 的识别效果统计
Table 3 Statistical analysis of recognition effect of AO and AFL

pagenumber_ebook=124,pagenumber_book=122
根据表4 的统计结果，AFL 对于SA、SB、CB雷达的综合识别效果都较好，CSI 评分0.9 以上。其中，对于降水回波的识别准确率超过0.99，CB 雷达对对流云的识别准确率甚至为1，说明在大样本量下，AFL 对于降水回波的误识情况是较少的。对于CC 雷达，由于探测范围较小（150 km），层云降水回波很容易达到探测范围，且回波纹理平滑，由此计算的特征参数与径向干扰回波十分相似，造成误识较多，识别准确率仅为0.804，CSI 评分0.746。

表4 AFL 对不同类型雷达资料的识别效果统计
Table 4 Statistics of recognition effect of AFL for different type radar data

pagenumber_ebook=124,pagenumber_book=122
4.2 个例评估
使用AO 和AFL 对个例观测回波进行识别消除，分析对比算法的识别效果以及适用性。

pagenumber_ebook=125,pagenumber_book=123
图5 AO 和AFL 对图1a、1b 的质量控制效果对比（a、b.图1a、1b 经AO 识别后的反射率因子，c、d.图1a、1b 经AFL 识别后的反射率因子）
Fig.5 Comparison of quality control effect between AO and AFL for Fig.1a and 1b（a，b.reflectivity after identification with AO for Fig.1a and 1b，respectively；c，d.reflectivity after identification with AFL for Fig.1a and 1b，respectively）

图5 给出AO 和AFL 对图1a、1b 观测回波的识别对比结果。两类算法均能有效识别0.5°和1.5°高低两个仰角的绝大部分窄径向和大面积径向的RI-I，识别面积基本相同，但0.5°仰角的RI-I 经AO 识别后会存在较少的孤立回波，以及在355°方位角附近遗留较强的回波（白色方框内），这是因为AO 仅使用单一特征参量判定，一旦特征参量不满足阈值条件时，则不能识别，会遗留较多的孤立点。

再分析图1c，雷达东北方存在强对流降水回波，西北方和东南方分别存在较窄的和大面积的RI-II，径向干扰回波的水平纹理十分粗糙。经过识别后（图6），AFL 较AO 能更好地识别窄径向和大面积径向的RI-II，同时能更好地保留强降水回波。由于AO 仅针对RI-I 的特征建立，故对于RI-II 识别效果较差。

文中还检验了AFL 对其他个例不同仰角的RII 和RI-II 的识别效果，均能较好识别（图略），因此，AFL 对于RI-I 和RI-II，包括高低仰角的窄径向、大面积径向的干扰回波，都能有效地识别和消除。

pagenumber_ebook=126,pagenumber_book=124
图6 AO（a）和AFL（b）对图1c 的质量控制效果对比
Fig.6 Comparison of quality control effect between AO（a）and AFL（b）for Fig.1c

质量控制算法的识别效果需综合考量其对于非气象回波和气象回波的识别情况。由于径向干扰回波一般能够延展到最大探测距离，当有效的距离库数越多，AFL 采用的特征参数RREF的值越大；DB也仅在雷达较远探测距离存在有效值才可以计算；对于SA、SB 和CB 雷达，一般的降水回波很少能贯穿整个方位角，故该算法一般不会对降水回波造成误识，统计结果也能说明此情况。对于图7 所示的这类长条状的降水回波，降水回波在方位角192°—200°基本已到达最远探测距离（RREF值较大），径向前后反射率因子变化不大（DB较小），特征参数满足条件被误识，但相对于AO 误识有减少。对于图8 所示的层云降水回波，由于CC 雷达最大探测距离仅为150 km，降水回波很容易达到，同时层云降水回波水平纹理平滑，在此情况下AFL 中的特征参数值与径向干扰回波高度相似，造成严重的误识别，误识的面积比AO 多。

因此，对于SA、SB 和CB 雷达，当降水回波存在较远探测距离且回波水平纹理变化较小时，AFL 可能对其造成误识别；对于CC 雷达探测的层云降水回波，AFL 存在严重误识别的情况。

4.3 业务应用效果
中国气象局气象探测中心2018 年6 月发布“综合气象观测数据质量控制系统”，该系统为国家级观测数据质量控制实时业务系统，其中涵盖对天气雷达数据的获取监视、质量控制、诊断勘误以及评估改进功能，值班员对系统数据质量无间断监控。AFL 也被应用到该系统中，以用于处理SA、SB 和CB 雷达产生的径向干扰回波，对于CC 雷达等探测范围较小的雷达未使用。为了便于值班员的监控，系统网页诊断勘误部分显示了未经质量控制和经过质量控制的全国天气雷达组合反射率拼图，若算法对于径向干扰回波识别消除失效，值班员会手动进行勘误并记录。

中国气象局气象探测中心业务系统中天气雷达拼图的时间分辨率为10 min，自2018 年7 月至2019 年5 月共生成48240 时次拼图数据，根据值班员的反馈以及记录，AFL 使用近1 年，人工质量控制441 个时次的径向干扰数据（包括对CC 雷达等未使用AFL 数据的人工勘误），占总数的0.91%。表明算法具有稳定的运行能力及良好的质量控制效果，能够识别并去除绝大多数的径向干扰回波。

pagenumber_ebook=127,pagenumber_book=125
图7 2017 年 7 月 15 日 12 时 11 分潍坊SA 雷达 PPI(a.原始反射率因子，b.AO 识别后反射率因子，c.AFL 识别后反射率因子)
Fig.7 PPI of Weifang SA radar at 12:11 UTC 15 July 2017(a.reflectivity，b.reflectivity after identification with AO，c.reflectivity after identification with AFL)

pagenumber_ebook=127,pagenumber_book=125
图8 2017 年 7 月 2 日01 时 36 分张掖CC 雷达 PPI(a.原始反射率因子，b.AO 识别后反射率因子，c.AFL 识别后反射率因子)
Fig.8 PPI of Zhangye CC radar at 01:36 UTC 2 July 2017(a.reflectivity，b.reflectivity after identification with AO，c.reflectivity after identification with AFL)

由于二次产品的质量高低均是建立在雷达基数据的基础上，AFL 能提升数据质量，对于业务上降水估测等定量应用也能产生一定正作用。2019 年6 月21—22 日，华中及华南的部分地区有大到暴雨，并伴有雷暴大风等强对流天气。图9 为6 月21 日03—06 时华东局部地区的3 h 降水估测，降水估测是在形成混合扫描反射率（肖艳姣等，2008）的基础上，利用Z-R 关系将反射率因子转换成降水率累计得到的。使用AFL 前，降水估测相应位置有微弱降水（图9a），使用AFL后条辐状的弱降水消失（图9b）。

pagenumber_ebook=128,pagenumber_book=126
图9 2019 年6 月21 日03—06 时3 h 局部降水估测AFL 前（a）后（b）拼图
Fig.9 3 h local precipitation estimation mosaic before（a）and after（b）identification with AFL from 03:00 UTC to 06:00 UTC 21 June 21 2019

5 结论
针对业务中频繁出现的径向干扰回波，根据收集的数据，提取特征参数，形成一种基于模糊逻辑的识别算法，并对算法开展效果评估，结果表明：

（1）RREF、DB、TDBZ和SPIN是4 个能够反应径向干扰回波和降水回波差异的特征参数，RI-I 和SC、CC 雷达的SPIN概率分布重合较多，因此RI-I 选择RREF、DB、TDBZ作为算法输入，RI-II 选择RREF、DB、TDBZ和SPIN作为算法输入。这些特征仅在径向方向上提取，能够与业务上天气雷达单径向流数据传输相结合，实现快速质量控制。

（2）AFL 相对于AO，整体的CSI 评分提高0.232，对于径向干扰回波的识别准确率大幅度提高（包括大面积径向干扰回波），对于降水回波的识别准确率略微提高。

（3）AFL 对于SA、SB 和CB 雷达的综合识别效果较好，CSI 评分超过0.9，当降水回波存在较远探测距离且回波水平纹理变化较小时，算法可能对其造成误识别；对于CC 雷达的识别准确率稍差，CSI 评分0.746，对于层云降水回波，存在严重误识别的情况。

（4）将AFL 应用于SA、SB 和CB 雷达资料实时业务质量控制，结果表明算法具有较高的稳定性，且对于降水估测等定量应用具有正贡献。

AFL 能够有效识别区分径向干扰回波，提高资料的可用性，但对于探测距离较短的雷达资料误识较多，严重时层云降水回波基本被识别剔除，在此情况下的定量降水估测会严重低估，因此，还需要寻找更有效的方法来识别此类雷达的径向干扰回波。由于各类非气象回波的识别用到了许多相同的特征参数，可以基于模糊逻辑建立一整套非气象回波的识别方法，提高业务运行效率及识别准确率。

参考文献

杜言霞，刘光普，邱光敏.2017a.天气雷达径向干扰回波去除方法的比较.干旱气象，35（3）：502-506. Du Y X，Liu G P，Qiu G M.2017a.Comparison of radial interference echo removal methods for weather radar.J Arid Meteor，35（3）：502-506（in Chinese）

杜言霞，陈斯智，陈州川等.2017b.天气雷达远距离单频点电磁径向干扰回波的特征分析及识别.气象研究与应用，38（4）：65-68. Du Y X，Chen S Z，Chen Z C，et al.2017b.Analysis and recognition of the characteristics of the remote single frequency electromagnetic interference radial echo of the weather radar.J Meteor Res Appl，38（4）：65-68（in Chinese）

高拴柱，矫梅燕，毕宝贵等.2007.相关法雷达反演风场在台风监测和分析中的应用.气象，33（1）：38-44.Gao S Z，Jiao M Y，Bi B G，et al.2007.Wind field retrieve from radar echo with correlation method and its application to typhoon monitoring and analysis.Meteor Mon，33（1）：38-44（in Chinese）

郭春辉，王旭，袁微.2014.多普勒天气雷达径向干扰回波的识别与消除.气象水文海洋仪器，31（2）：24-29，32. Guo C H，Wang X，Yuan W.2014.Identification and removal of radial interference echo of Doppler weather radar.Meteor，Hydrol Mar Instrum，31（2）：24-29，32（in Chinese）

江源.2013.天气雷达观测资料质量控制方法研究及其应用[D].南京：南京信息工程大学，35-37.Jiang Y.2013.Meteorological radar data quality control study and application[D].Nanjing：Nanjing University of Information Science and Technology，35-37（in Chinese）

刘黎平，吴林林，杨引明.2007.基于模糊逻辑的分步式超折射地物回波识别方法的建立和效果分析.气象学报，65（2）：252-260. Liu L P，Wu L L，Yang Y M.2007.Development of fuzzy-logical two-step ground clutter detection algorithm.Acta Meteor Sinica，65（2）：252-260（in Chinese）

盛春岩，薛德强，雷霆等.2006.雷达资料同化与提高模式水平分辨率对短时预报影响的数值对比试验.气象学报，64（3）：293-307. Sheng C Y，Xue D Q，Lei T，et al.2006.Comparative experiments between effects of Doppler radar data assimilation and inceasing horizontal resolution on short-range prediction.Acta Meteor Sinica，64（3）：293-307（in Chinese）

孙继松，雷蕾，于波等.2015.近10 年北京地区极端暴雨事件的基本特征.气象学报，73（4）：609-623. Sun J S，Lei L，Yu B，et al.2015.The fundamental features of the extreme severe rain events in the recent 10 years in the Beijing area.Acta Meteor Sinica，73（4）：609-623（in Chinese）

文浩，刘黎平，张持岸等.2016.新一代天气雷达地物回波及电磁干扰质控算法业务试运行评估.气象科学，36（6）：789-799. Wen H，Liu L P，Zhang C A，et al.2016.Operational evaluation of radar data quality control for ground clutter and electromagnetic interference.Scientia Meteor Sinica，36（6）：789-799（in Chinese）

肖艳姣，刘黎平，杨洪平.2008.基于天气雷达网三维拼图的混合反射率因子生成技术.气象学报，66（3）：470-473. Xiao Y J，Liu L P，Yang H P.2008.Technique for generating hybrid reflectivity field based on 3-D mosaicked reflectivity of weather radar network.Acta Meteor Sinica，66（3）：470-473（in Chinese）

杨川.2012.新一代天气雷达资料质量分析和处理方法研究[D].成都：成都信息工程学院，26-30.Yang C.2012.Study of method for radar data quality control with new generation weather radar[D].Chengdu：Chengdu University of Information Technology，26-30（in Chinese）

俞小鼎，姚秀萍，熊廷南等.2006.多普勒天气雷达原理与业务应用.北京：气象出版社，30-32.Yu X D，Yao X P，Xiong T N，et al.2006.Doppler Weather Radar Principle and Operational Application.Beijing：China Meteorological Press，30-32（in Chinese）

郑媛媛，俞小鼎，方翀等.2004.一次典型超级单体风暴的多普勒天气雷达观测分析.气象学报，62（3）：317-328. Zheng Y Y，Yu X D，Fang C，et al.2004.Analysis of a strong classic supercell storm with Doppler weather radar data.Acta Meteor Sinica，62（3）：317-328（in Chinese）

周红根，柴秀梅，胡帆等.2008.新一代天气雷达回波异常情况分析.气象，34（6）：112-115. Zhou H G，Chai X M，Hu F，et al.2008.Analysis of abnormal echo of new generation weather radar.Meteor Mon，34（6）：112-115（in Chinese）

Bao X W，Wu D，Lei X T，et al.2017.Improving the extreme rainfall forecast of typhoon Morakot（2009）by assimilating radar data from Taiwan Island and mainland China.J Meteor Res，31（4）：747-766

Berenguer M，Sempere-Torres D，Corral C，et al.2006.A fuzzy logic technique for identifying nonprecipitating echoes in radar scans.J Atmos Oceanic Technol，23（9）：1157-1180

Cho J Y N.2010.Signal Processing Algorithms for the Terminal Doppler Weather Radar：Build 2.Washington，DC：MIT Lincoln Laboratory，79pp

Cho J Y N.2011.Analysis of 5-GHz U-NII Device Signals Received by the PSF TDWR.Lexington，MA：MIT Lincoln Laboratory，28pp

Cho J Y N.2017.A new radio frequency interference filter for weather radars.J Atmos Oceanic Technol，34（7）：1393-1406

FAA.1995.Specification：Terminal Doppler Weather Radar with Enhancements.Federal Aviation Administration Doc.FAAE-2806c，142pp

Joe P，Scott J，Sydor J，et al.2005.Radio local area network（RLAN）and Cband weather radar interference studies∥Proceedings of the 32nd AMS Radar Conference on Radar Meteorology.Albuquerque，New Mexico：AMS

Kessinger C，Ellis S，van Andel J.2003.The radar echo classifier：A fuzzy logic algorithm for the WSR-88D∥Proceedings of the 3rd Conference on Artificial Intelligence Applications to the Environmental Science.Long Beach：Amer Meteor Soc

Mendel J M.1995.Fuzzy logic systems for engineering：A tutorial.Proc IEEE，83（3）：345-377

NOAA.2006.Federal Meteorological Handbook NO.11：Doppler radar meteorological observations part D，WSR-88D unit description and operational applications.Washington：Office of the Federal Coordinator for Meteorological Services and Supporting Research，6-7

Rennie S J，Curtis M，Peter J，et al.2015.Bayesian echo classification for Australian single-polarization weather radar with application to assimilation of radial velocity observations.J Atmos Oceanic Technol，32（7）：1341-1355

Saltikoff E，Cho J Y N，Tristant P，et al.2016.The threat to weather radars by wireless technology.Bull Amer Meteor Soc，97（7）：1159-1167

Steiner M，Smith J A.2002.Use of three-dimensional reflectivity structure for automated detection and removal of nonprecipitating echoes in radar data.J Atmos Oceanic Technol，19（5）：673-686

Tang L，Zhang J，Langston C，et al.2014.A physically based precipitationnonprecipitation radar echo classifier using polarimetric and environmental data in a real-time national system.Wea Forecasting，29（5）：1106-1119

Vaisala.2014.User's Manual：RVP900 Digital Receiver and Signal Processor.Vaisala：Vaisala Oyj，513pp

WEN Hao1 ZHANG Lejian1 LIANG Haihe1 ZHANG Yang2

1.Meteorological Observation Centre of China Meteorological Administration，Beijing 100081，China
2.State Key Laboratory of Severe Weather，Chinese Academy of Meteorological Sciences，Beijing 100081，China

Abstract When the radar signal processor is abnormal,the signal of radar is subject to external electromagnetic interference,the radar observation will have strip-shaped interference echoes along the radial direction.Such non-meteorological echoes must be identified and eliminated before quantitative applications such as radar precipitation estimation.In this study,a new algorithm for radial interference echo identification is presented based on fuzzy logic.To adapt to the current volumetric data and the upcoming single radial stream data for operational application,the characteristic parameters only on the radial direction are extracted,and the membership functions are established,so that a value that measures the probability of each bin being affected by radial interference echo is get.The algorithm is evaluated and implemented in operational system,and the results show that it can distinguish radial interference echoes from precipitation echoes very well.Compared with the original operational algorithm,the CSI score is improved by 0.232,which can significantly improve the accurate identification of radial interference echoes and slightly reduce the false elimination of precipitation echoes.It has robust operation ability and can make positive contribution to quantitative applications (such as precipitation estimation).

Key words Weather radar，Radial interference echo，Fuzzy logic