雷达基数据处理与雷达回波和地闪资料叠加设计.pdfVIP

雷达基数据处理与雷达回波和地闪资料叠加设计 余海1，张廷龙1，高皴1，傅仁壮2 (1海南省防雷中心．海南海口，570203：2海南省气象服务中心，海南海口，570203) 摘要：本文主要介绍基丁MATLABGUl编程技术实现对海u新一代s波段的天气雷达(CINRAD／SA)基 数据的读取显示及雷达回波与地闪定位数据叠加。利用MATLABGUI界猫编程手段，设计制作了雷达基数 据读取显示软件，该软件能够直观、方便地查看雷达凹波的空间分布及三维视图。并且实现了雷达回波与 对应时段地闪定位数挺的叠加并且叠加结果基本一致，能较好地结合雷达回波和地闪定位资料实现对雷器 天气过程的分析。通过对比雷达回波与地闪定位数据的叠加图，本文还简要分析了“威马逊”台风引起的闪 电活动。结粜表明：地闪活动多发生在台风后方的螺旋云带，随着台风向西北方向移动，地cj_；l活动由鼯南 海面发展至海南岛昌江、白沙等地，并在7月18日19．2l点时段，海南岛的地闪活动发展最为旺盛，随后 逐渐沿东南．西北方向发展到海南岛西北部海域。当台风于19曰凌晨7点10分左右登陆广1西，海面闪电活 、动再次发展至最旺盛阶段，随后闪电活动开始减弱。 关键词：新一代S波段天气雷达，基数据读取及显示，数据叠加，威马逊，螺旋云带 1． 引言 美国在20世纪80年代初开始设计为气象业务使用的多普勒天气雷达WSR88D，称为 新一代天气雷达(NEXRAD)【I】并设计了丰富的应用处理软件‘2剖。我国在WSR-88D的基 础上设计生产了CINRADS波段和C波段的雷达。目前新一代天气雷达(CINRAD)已在 全国范围内架设并用于气象业务工作，利用编程手段实现CINRAD／SA雷达基数据的读取及 显示，提取反射率、径向速度及谱宽等基本数据，计算出云顶亮温、强回波区顶高及回波体 积等相关气象参数，能更加高效地提高科研及业务水平及效率，国内许多学者开展了这部分 工作[9-131。运用雷达数值产品开展的科学研究和气象业务服务，为研究暴天气过程和雷电预 警预报技术提供了重要的思路与方法。国外有学者乖j用雷达资料分析了云闪频次与雷达圈波 反射率的关系㈤1和地闪频次与反射率的相关性‘‘16】，国内学者对于雷暴的分析或单独利用雷 达数据或结合闪电探测技术与雷达回波‘17-21】。 本文在结合前人的基础上，利用MATLABGUI编程技术除了实现了基数据的读取、显 示等基本功能，在此基础上，将地闪定位数据叠加到对应时段的雷达回波图中，实现对雷暴 天气过程中闪电活动与雷达回波的相关性的分析，为下一步雷电预警预报系统的开发做准 备。并利用2014年7月18目的“威马逊”超级台风的雷达同波及地闪定位资料进行了叠加技 术的验证。 基会项目：国家自然科学基会项目和海南省气象局科技创新项目(HN2013MSl7)苁同资助 2．基数据处理 CINRAD／SA雷电基数据有着统一的数据存储格式，每个体扫存储为一个二进制文件。 每个体扫文件按径向数据的方式进行排列，排列顺序按照仰角由低至高排列。径向数据的长 度固定为2432个字节，程序设计过程中对每一个体扫文件按径向数据长度进行读取处理， 读取每个径向数据中记录的扫描时间、体扫或仰兔扫描开始结束标示符、反射率谱宽速度开 始位置、距离库数及相应数据，，如图1所示的径向数据格式。 图lCINRAD／SA基数据径向数据格式 RadialdataformatofCINRAD，SAbasedata Fig．1 表1．仰角、方位角及径向数据状态说明 Table．1Instruction andradialdatastate ofElevation，azimuth 利用MATLAB中的fread函数按照二进制的形式读取基数据，基数据的处理首先需要 判断径向数据状态，当状态值为3时，表示体扫开始的第一条径向数据，这时程序才可进入 下一步的仰凫逐层扫描，当径向数据状态为0时，仰角扫描开始，并按照表l中的仰角和方 位角说明计算出每一条径向数据的仰角和方位角。由于SA雷达在低层每个仰角上扫描两次， 程序中，在保存仰角数据，对同一仰