多普勒天气雷达的对流风暴自动识别研.pdfVIP

多普勒天气雷达的对流风暴自动识别研究* 张晶晶 （九江市气象台，江西，九江 332000 ） 摘要 鉴于雷达软件版权的问题以及需要专门的显示软件，使得用户对产品的使用受到了一定的限制，设计 了基于多普勒天气雷达的对流风暴单体识别程序。利用 2011 年 6 月 14 日九江一次暴雨过程的多普勒天气 雷达体扫资料对该程序进行试算，并与中美合资敏视达公司的CINRAD-SA雷达Build 10.8软件系统生成的 风暴识别产品进行比较分析，结果表明，该程序对风暴单体具有较好的识别效果，起到较好的预警作用。 关键词: 多普勒天气雷达 核抽取 对流风暴识别算法 程序模块化 0 引言 各种灾害性强对流天气过程常伴随狂风、暴雨、强雷电甚至冰雹、龙卷等恶劣天气现象， 往往给国民经济和人民生命财产安全造成重大危害。多普勒天气雷达全天候的探测能力和反 映大气层云雨生消演变过程的动态过程的能力，大大加强了对中小尺度天气系统的探测和预 警能力，并为开展短时灾害性天气系统及洪涝灾害的监测和预警，数值预报，定时、定点的 精确天气预报奠定了基础。 多普勒天气雷达作为探测强对流风暴天气的有效工具，自二十世纪六、七十年代以来， 其探测结果为强风暴系统结构和演变的研究提供了重要依据，大大促进了关于风暴的算法研 究。美国 WSR －88D （Weather Surveillance Radar – 1988 Doppler ）的雷达软件与后来的升级 的数据处理方法，为风暴识别、跟踪和预报及强对流天气产品的业务化方案提供了成功的范 例。八十年代，美国 WSR-88D 风暴系列算法在全美 NEXRAD 业务系统中广泛投入使用， 取得较好的效果。但是，用户在使用过程中也发现当风暴单体数目较多，移动不规则，交叉、 合并时存在诸多不足,根据用户的反馈意见，Mansur[1]，Arthur Witt and J.T. Johnson[2]对该算 法做了很大的改进，我们称之为改进的强中心跟踪法。1998 年，J. T. Johnson, Pamela L. [3] 等 人提出利用七个反射率因子识别阈值来替代此前唯一的一个阈值，同时采取特征核抽取技 术，并对空间相距较近的多个风暴单体进行合并或者删除处理，该算法应用于 WSR －88D 的 Build9.0 风暴系列算法中，由此算法得到的数据被广泛用于其它产品（HI 、SS、SRM、 SRR、WER 、TVS 等）中。                                                                *作者简介：张晶晶 （1986-），女，汉，主要从事短期天气预报业务与研究，Email ：zhangjingjing217@163.com    我国的 CINRAD-SA、SB 和其他雷达的软件包目前都采用了美国的 NEXTRAD 10.8 版 本的风暴识别算法[4] 。但由于敏视达公司的气象雷达软件版权的问题以及需要专门的显示软 件，使得用户对产品的使用受到了一定的限制[5] 。鉴于此，本文使用VC++6.0[6]开发工具， 为多普勒天气雷达设计了风暴单体的识别软件，使用多阈值、核抽取以及相近单体处理等多 项新技术，以改善其算法在业务上所表现出的不足。 1 对流风暴的识别算法 风暴单体识别算法主要参考美国 WSR-88D 的 Build10.8 的 SCIT 算法[1-3]思路，把“风 暴”看作是一个三维立体结构，在识别时首先在一维径向上搜索反射率因子大于一定阈值的 连续区域，称之为风暴段(Segment)，然后依据同一层次内相邻方位间风暴段的相关性，将 满足一定条件的风暴段合成为一个二维区域，称之为风暴分量(Component)，最后根据风暴 的垂直相关性将不同层次上的风暴段合成为一个具有三维结构的风暴单体。整个过程大致可 以分为以下几部分：风暴段的搜索、风暴分量的合成及风暴单体的组成。 1.1 风暴段的搜索 算法一维部分的目的是在径向反射率因子中识别“风暴段”，段是指沿径向排列的、反 射率因子大于或等于特定阈值的一组相邻距离库。识别中主要引入三个阈值：反射率因子阈 值(REFLECTIVITY1-7)，低于反射率差值阈值(DROPOUT REF DIFF，default=5 dBZ)，低于 反射率值