基于高分辨率高程数据统计分析新一代天气雷达组网的地...
基于高分辨率高程数据统计分析新一代天气雷达组网的地形遮挡影响*柳云雷1 李昌兴2 张乐坚2 文 浩2 王曙东3
1. 大连市气象局气象装备保障中心，大连， 116001
2. 中国气象局气象探测中心，北京， 100081
3. 中国气象局公共气象服务中心，北京， 100081
摘 要 新一代天气雷达由于受到地形限制产生波束遮挡导致波束能量衰减，从而造成雷达探测回波强度偏弱、雷达定量估测降水结果失真，因此对于雷达波束遮挡情况的统计和分析是一项重要的基础研究工作。利用SRTM （Shuttle Radar Topography Mission）数字高程数据对中国目前业务运行的212 部新一代天气雷达波束遮挡情况进行模拟计算分析。计算结果包括雷达单站遮蔽角、VCP21 模式0.5°、1.5°、2.4°、3.4°、4.3°仰角波束遮挡率、混合扫描及分区混合扫描波束遮挡率、雷达单站探测范围覆盖情况；计算并绘制全国天气雷达组网遮挡率拼图，统计全国天气雷达组网遮挡情况；利用2019 年8 月广东省11 部天气雷达基数据对比验证单站及组网遮挡计算结果。结果表明雷达组网探测面积覆盖率超过70%，整体覆盖效果较好，遮挡计算结果与实际数据对比验证结果高度一致，对雷达数据订正、降水估测等产品具有正贡献。
关键词 新一代天气雷达， 波束遮挡， 遮蔽角， 雷达组网遮挡
1 引 言
天气雷达作为探测降水系统的主要技术手段之一，其在短时临近预报和强对流天气的监测、预警中发挥着重要作用。经过20 多年的建设和发展，中国目前业务运行的新一代天气雷达（CINRAD）已达到212 部，所组成的雷达观测网基本覆盖中国沿海及降水较多的区域，国土面积覆盖率超过70%，在气象监测、预警和气象防灾、减灾工作中扮演着重要角色（许小峰，2003）。
降水估计是基于基本反射率因子与降水率的正相关关系（俞小鼎等，2006），地形遮挡导致波束能量衰减，造成雷达探测降水回波强度偏弱，是估测降水结果失真的主要原因之一（Westrick，et al，1999；Dinku，et al，2002；Gabella，2004；勾亚彬等，2014）。中国地势复杂，城市高大建筑大量建设，对雷达周边探测环境产生较大影响（范思睿，2016），导致中国天气雷达组网探测依旧存在盲区，遮挡问题仍难以解决（景号然等，2019；李柏等，2013），所以对于天气雷达探测环境遮挡情况的分析和计算受到长期关注和研究（王曙东等，2011）。
中国对雷达地形遮挡的研究从早期的人工测量，发展为地理信息系统（GIS）计算，以及目前的数字高程模型（DEM）数据分析。早期通过人工测量站址四周遮挡的经纬度、海拔高度，计算雷达单站遮挡仰角和等射束高度图（邓志等，1999）。万玉发等（2000）研发了相关业务软件，该方法测量结果准确可靠，但总体工作量较大。赵瑞金等（2003）利用地理信息系统的地形数据模拟计算雷达单站遮挡情况，通过与人工测量方法进行对比发现，地理信息系统的地形数据计算效果更优、效率大幅度提高。目前雷达遮挡相关分析主要基于数字高程模型数据和组网观测，利用地形高程数据计算波束阻挡率和组网拼图有效数据区域，提出波束阻挡的反射率强度订正方法（杨洪平等，2009）。王曙东等（2011）利用航天飞机雷达地形测绘数据对158 部雷达进行了地形遮挡和覆盖能力研究，得到青藏高原和西北区域覆盖较少，东部地区覆盖较好的统计结果，通过对云南组网雷达探测冰雹覆盖情况进行研究，结果证实受地形遮挡存在部分盲区（石宝灵等，2018）。针对雷达地形遮挡的研究也推动了订正算法和业务应用的发展，包括地形遮挡区域雷达波束部分遮挡的订正算法研究有效提高了定量估测降水精度（勾亚彬等，2017；蒋秋菲等，2017）；分析雷达单站业务体扫模式VCP21 的波束阻挡系数，统计雷达单站各仰角层遮挡比例（朱丹等，2018）；研发遮挡订正系统改善基数据可用性（郭春辉等，2018）等工作。
美国强风暴实验室利用高分辨率SRTM（Shuttle Radar Topography Mission）数据计算雷达波束遮挡（Farr，et al，2007），结合反射率垂直廓线（VPR）建立雷达质量指数（RQI）（Zhang，et al，2012），通过 NSBM（Nonstandard Blockage Mitigation）系统订正雷达波束遮挡（Tang，et al，2013），并利用组网天气雷达配合地面观测数据搭建了MRMS（ Multi-Radar Multi-Sensor system） 系 统（Qi，et al，2017；Zhang，et al，2016），可每 5 min 生成一组全国组网的雷达质量指数用于基本反射率、定量估测降水等产品的订正计算。
由于利用高分辨率高程地形数据对雷达探测环境进行遮挡计算的计算量巨大，中国目前大多相关研究都基于各省或几部雷达，同时天气雷达站的新建和迁址也对组网雷达的覆盖有效范围产生影响。文中利用SRTM 90 m 分辨率高程数据对目前组网业务运行的212 部天气雷达进行单站和组网遮挡情况分析，统计雷达覆盖情况。绘制单站遮蔽角、径向遮挡、混合扫描遮挡、分区混合扫描遮挡以及分区混合扫描对应仰角；绘制雷达组网遮挡图并对整体数据进行统计分析，给出全国雷达单站及组网遮挡情况的统计分析结果，利用2019 年8 月广东省雨量计观测数据和11 部雷达基数据对计算结果进行实际数据验证。统计分析结果对雷达基本反射率、雷达定量估测降水等产品的遮挡订正具有正贡献。
2 资料与分析方法
2.1 数据资料来源
雷达探测环境遮挡计算主要依赖于雷达覆盖区域的地形高程资料，高程资料的准确性以及雷达站经纬度、馈源海拔高度对于遮挡计算具有很大的影响。所用的SRTM 地形数据为SRTM3 公开数据，由美国太空总署和国防部国家测绘局联合测量。数据版本为国际热带农业中心利用新的插值算法得到的SRTM v4.1 地形数据。该版本数据平均误差-0.35 m，绝对误差小于16 m，高程数据误差引起的遮挡计算误差对近距离遮挡影响较为明显，当距离较远时由于波束本身发散其带来的影响有所减弱，以SA 雷达0.5°仰角为例在水平距离0.96 km之内可能引起大于10%的误差，水平距离超过10 km时由于高程数据引起的误差可小于1%，该数据是目前国际上使用的主要高程数据之一，总体可靠性较高。数据来源于中国科学院计算机网络信息中心地理空间数据云平台（www.gscloud.cn，可公开下载）。雷达探测环境遮挡计算过程中所使用的雷达站经纬度和海拔高度采用最新更新数据，并与卫星地图进行核对以确保其准确可靠。广东省2019年8 月雨量数据来自香港天文台发布的中国气象局国家气象中心数据（www.hko.gov.hk/tc/wxinfo/gdrf/gd_prfac.shtml，公开发布）。
2.2 雷达测高公式
雷达发射电磁波的传播路径计算需要考虑地球曲率与地球大气折射对其产生的影响，按照标准大气情况下，测高示意如图1 所示。
雷达遮蔽角和径向遮挡图中黑色部分为雷达探测可覆盖区域，白色部分为遮挡区域，仰角遮挡图外部曲线愈向内凹说明该方位遮挡仰角越大，同时对应方向的径向探测距离越短，从图中可以看到仰角遮挡与径向遮挡是相互匹配的，仰角遮挡越大径向探测距离在对应仰角层也会受到影响。拉萨雷达站地处西藏高原，馈源高度3.67 km，喜马拉雅山脉及高原大地形都对雷达产生很大的遮挡影响，其仰角遮挡最大方位位于303°，遮挡角度达到9.33°，其低层5 层仰角均存在不同程度遮挡。兰州等其他3 部雷达在低层仰角均存在不同程度的遮挡，随着仰角抬升遮挡情况有所改善。雷达仰角遮挡图配合径向遮挡图或等射束高度图能够较好地反映雷达周边探测环境的好坏，对于建站选址评估、后期回波数据订正具有重大效益。
pagenumber_ebook=160,pagenumber_book=707
pagenumber_ebook=160,pagenumber_book=707
图 1 基于等效地球半径雷达测高示意
Fig. 1 Altitude measurement of radar based on equivalent radius of the Earth
式中，h 为天线高度（km），δ 为仰角（°），R 为斜距（km），Rm 为等效地球半径（km），H 为波束中心轴线在斜距R 处离地面的高度（km），在标准大气情况下，等效地球半径（Rm）为真实地球半径的4/3 倍。
2.3 雷达波束遮挡比例
雷达波束遮挡比例（以下简称为遮挡率）采用SRTM 高程数据和雷达站基础信息进行匹配计算。对雷达探测范围进行采样计算，根据雷达站信息计算目标点经纬度与探测高度，对比该经纬度地形数据高度，计算该点某仰角时波束截面遮挡情况，计算时采用标准大气等效地球半径（图2），h 为天线馈源高度，Rm 为等效地球半径，AR 为探测距离、θ 为仰角值（°）、GR 为雷达测量高度，GT 为SRTM 高程数据的海拔高度，BE 为过T 点垂直于波束方向的法线交AR 于D 点。
pagenumber_ebook=160,pagenumber_book=707
pagenumber_ebook=160,pagenumber_book=707
图 2 雷达波束遮挡比例示意
Fig. 2 Blockage percentage of radar beam
式中，R、Rm、OA、GR、GT 单位均为 km，θ、α 单位为°。计算过程中采样点数量决定最终计算量的大小，通过匹配不同参数进行实验后，综合考虑计算效率以及分辨率产生的偏差影响，采用距离分辨率90 m，方位角度分辨率0.01°，波束宽度按照0.95°进行遮挡比例计算。
3 天气雷达单站遮挡分析
天气雷达单站遮挡分析同样基于SRTM 高程数据对各个雷达站点的遮蔽角、径向遮挡、混合扫描遮挡、分区混合扫描遮挡情况进行计算统计。
3.1 遮蔽角及径向遮挡情况
雷达仰角遮挡图（也称为遮蔽角图）绘制主要用于表征标准大气折射情况下，雷达扫描区域范围内某一方向上由于地形或建筑等障碍物所产生的最大遮挡角度，即当波束低于该角度时就会产生部分或全部波束被遮挡。作为建站选址的一项重要参数，可以模拟新建站点的波束遮挡情况，也可用于对在用站点的波束遮挡情况进行分析。
通过计算全国212 部业务运行天气雷达的单站仰角遮挡及径向遮挡并绘图，选取拉萨（CD）、兰州（CC）、韶关（SA）、宜春（SA）4 站 VCP21 模式下低层 5 个仰角（0.5°、1.5°、2.4°、3.4°、4.3°）进行举例，其中雷达有效最大探测距离按照C 波段150 km、SA/SB 230 km，计算结果如图3 所示。
江西省茶文化历史悠久，自古就有“唐载茶经，宋称绝品，明清入贡，中外驰名”的美誉。江西省的地形以丘陵、山地为主，属中亚热带温暖湿润季风气候，雨水充足，这种自然环境十分适合茶叶的种植和生长。
3.2 雷达各仰角波束遮挡率
雷达波束遮挡包括部分遮挡与完全遮挡，波束遮挡率可较好反映某一方位径向上的遮挡具体情况，相比于遮蔽角和径向遮挡图更精细化地反映波束实际遮挡情况（图4）。
以宜春雷达站为例，从其径向遮挡及遮蔽角（图3d）中可以看出，其北部约有28°方位范围、南部约有129°方位范围的遮挡情况。通过计算波束遮挡率（图4）可以发现0.5°仰角存在较多的波束部分遮挡，在径向遮挡中并未直观体现，随着仰角抬升北部遮挡明显减弱，2.4°仰角时仅南部有约53°方位范围存在遮挡，在4.3°仰角时仅有极少量部分遮挡情况。因此，可以通过波束遮挡率在雷达估算定量降水时对部分波束遮挡所带来的回波减弱情况进行相应回波强度订正。
3.3 混合扫描遮挡率分析
混合扫描方式在雷达降水算法中广泛应用，不考虑地形情况下，按照0—20 km 使用3.4°仰角、20—35 km 使用 2.4°仰角、35—50 km 使用 1.5°仰角、50—230 km 使用0.5°仰角数据进行反射率因子预处理。根据混合扫描特性，匹配相应雷达波束遮挡情况进行分析（图5）。
根据图中结果可知，混合扫描遮挡率基本趋近于雷达0.5°仰角的波束遮挡率，不同之处主要在于50 km 以内波束遮挡情况会由于高仰角数据的引入而有所改善，同时混合扫描方式并未考虑地形影响，所以对于探测范围50 km 以外的波束遮挡情况并未减弱。
3.4 分区混合扫描遮挡率分析
在考虑地形影响的情况下，美国亚利桑那州凤凰城WSR-88D 雷达通过引入分区混合扫描方式来减少雷达波束遮挡所带来的影响（Fulton，et al，1998）。其按照混合扫描方式进行计算的同时，对于低层4 层仰角某径向的波束遮挡率低于10%的不进行反射率因子订正，11%—29%反射率因子订正+1 dBz，30%—43%反射率因子订正+2 dBz，44%—55%反射率因子订正+3 dBz，56% —60%反射率因子订正+4 dBz，对于遮挡部分大于60%的使用上一层仰角对应值。根据此理论进行模拟计算，分区混合扫描遮挡率如图6 所示。
pagenumber_ebook=162,pagenumber_book=709
图 3 雷达遮蔽角与低层 5 个仰角的径向遮挡图（a. 拉萨，b. 兰州，c. 韶关，d. 宜春）
Fig. 3 Defilade angle and radial blockage of radar beam （a. Lhasa，b. Lanzhou，c. Shaoguan，d. Yichun）
pagenumber_ebook=163,pagenumber_book=710
pagenumber_ebook=164,pagenumber_book=711
图 5 混合扫描方式雷达波束遮挡率（a. 拉萨，b. 兰州，c. 韶关，d. 宜春）
Fig. 5 Blockage percentage of radar beam for hybrid scan mode （a. Lhasa，b. Lanzhou，c. Shaoguan，d. Yichun）
通过分区混合扫描方式，对于遮挡较为严重区域其波束遮挡率有所改善，该方式可以在低层仰角受到较大的地形波束遮挡时，利用高层仰角反射率因子导出降水率，但应当注意高层仰角的数据对于低层仰角的近地面条件的代表性存在不足，尤其是对于距离雷达较远的探测区域和层状云降水区域。该方法是对遮挡造成的低仰角数据缺失进行填补的技术处理方式之一，并不能完全解决雷达波束遮挡所带来的影响。
3.5 天气雷达单站覆盖情况统计
基于SRTM 高程数据及波束遮挡率的计算，分别统计雷达各个仰角层的探测区域面积覆盖情况，以便于直观了解雷达探测效果。
统计结果利用不同颜色代表不同覆盖面积比例，标签数量代表该覆盖面积所对应雷达站数量及在全部站点中所占比例。根据VCP21 模式低层4 层仰角的统计结果（图7）显示，在0.5°仰角（图7a）情况下，全国212 部雷达有29 部实际探测面积小于理论值的25%，占雷达总数的13.68%；31 部大于25%且小于50%，占雷达总数的14.62%；39 部大于50%且小于75%，占雷达总数的18.40%；113 部大于75%，占雷达总数的53.30%。1.5°仰角（图7b）时，遮挡情况明显改善，覆盖面积大于理论值75%的共计178 部，占雷达总数的83.96%。在2.4°（图7c）和 3.4°（图 7d）仰角情况下，分别仅有 12 部和 6 部雷达覆盖面积小于理论值的75%。
pagenumber_ebook=166,pagenumber_book=713
图 6 分区混合扫描方式雷达波束遮挡率 （a1、b1、c1、d1）和仰角（a2、b2、c2、d2）（a. 宜春，b. 兰州，c. 拉萨，d. 韶关）
Fig. 6 Blockage percentage of radar beam for subarea hybrid scan mode （a1，b1，c1，d1） and the elevation （a2，b2，c2，d2）（a. Yichun，b. Lanzhou，c. Lhasa，d. Shaoguan）
pagenumber_ebook=166,pagenumber_book=713
图 7 雷达探测覆盖面积比例统计
Fig. 7 Statistics of radar coverage area proportion
4 天气雷达组网遮挡统计分析
4.1 雷达组网遮挡情况分析
将212 部天气雷达分别匹配SRTM 高程数据进行不同仰角的波束遮挡计算，绘制全国范围的波束遮挡率拼图，选取0.5°和1.5°仰角拼图为例，如图8 所示。
图中每个黑色圆点代表一部天气雷达，在0.5°仰角时探测遮挡面积比例为29.91%，华中、华东及南部沿海地区覆盖情况较好，基本不存在盲区；华北、东北覆盖情况一般，存在部分遮挡的情况；西北、西南地区由于青藏高原等大地形的影响，遮挡较为严重，覆盖区域较差。在1.5°仰角时遮挡比例下降至9.85%，角度抬升及组网雷达之间相互补充，整体覆盖效果较好，仅在西北和西南地区存在少量遮挡，雨区基本全覆盖。雷达组网遮挡面积统计结果如表1。
随着仰角抬升波束遮挡情况势必减弱，从0.5°到1.5°效果较为明显，遮挡面积减少20 个百分点。虽然仰角抬升遮挡有所减少但应当注意探测高度和距离也相应发生变化，VCP21 模式下其他高层仰角的组网和单站的探测覆盖效果均优于低层仰角。雷达组网遮挡拼图时仅考虑同经纬度覆盖情况，对于高差未进行修订。
4.2 高差阈值对于雷达组网遮挡的影响
雷达波束遮挡在组网拼图计算时，由于雷达站海拔高度的不同，随着探测距离的增加，来自不同雷达站的同经纬度数据高度上可能也相距甚远。在实际雷达回波拼图处理过程中，由于高差所带来的回波不连续性处理也是雷达数据处理的重要环节之一。文中在制作遮挡拼图过程中，通过设定不同的高差阈值，计算不同高度差阈值时天气雷达的组网遮挡情况，选取0.5°仰角高差阈值0.5、3 和5 km为例结果如图9。
考虑高度差阈值情况下，雷达的探测覆盖效果整体有所降低。不同高差阈值雷达组网遮挡比例统计如表2。
pagenumber_ebook=167,pagenumber_book=714
图 8 0.5°（a）、1.5°（b） 仰角雷达组网遮挡拼图
Fig. 8 Composite images of radar network blockage area with different elevation（a. 0.5°，b. 1.5°）
表 1 雷达组网遮挡面积比例统计
Table 1 Statistics of radar network blockage area proportion
pagenumber_ebook=167,pagenumber_book=714
pagenumber_ebook=168,pagenumber_book=715
图 9 0.5（a）、3（b）和 5（c） km 高差阈值雷达组网遮挡率拼图
Fig. 9 Composite images of radar network blockage percentage with different threshold values of height（a. 0.5 km，b. 3 km，c. 5 km）
表 2 不同高差阈值雷达组网遮挡面积比例统计
Table 2 Statistics of radar network blockage area proportion with differentthreshold values of height
pagenumber_ebook=169,pagenumber_book=716
0.5 km 阈值时，0.5°仰角遮挡面积比例达到46.88%，在华中及华南沿海具有较多的部分遮挡情况，东部沿海遮挡略有增加但整体依旧覆盖良好。随着高差阈值的增加组网遮挡情况有所改善，3 km阈值时遮挡面积比例为37.08%，华南沿海地区遮挡明显减少，5 km 阈值时南部沿海地区基本无遮挡，华北遮挡情况也有明显的改善。可以看出雷达组网遮挡情况受高差阈值影响较大。处理雷达拼图数据时根据回波顶高、降水类型等信息选取不同高差阈值匹配雷达遮挡情况有利于减少遮挡，同时提高回波连续性和空间一致性。
5 实际雷达基数据对比验证
利用广东省11 部天气雷达2019 年8 月VCP21模式0.5°仰角基数据反演降水量并进行累计，通过降水累计验证雷达单站及组网遮挡区域与模拟算法的一致性，其中降雨率算法采用默认公式Z=300R1.4，未做地物杂波剔除及回波订正减少算法影响原始数据结果。选取韶关、汕尾、梅州为例对比单站遮挡率与降水累计分布如图10。
计算遮挡率发现韶关站仅南部有少量无遮挡区域，汕尾西北部、梅州东南部存在大量遮挡区域，通过降水累计可明显看出上述遮挡区域对应的降水累计显著偏弱。相同时段韶关、汕尾、梅州雨量计降水统计如表3。
通过对比雷达与雨量计数据，汕尾、梅县雨量计统计数据与雷达降水累计对比结果基本一致，韶关雷达降水累计由于受到遮挡影响小于雨量计降水统计。
表 3 2019 年8 月韶关、汕尾、梅州雨量计降水量统计
Table 3 Statistics of rain gauge data collected at Shaoguan，Shanwei and Meizhou in August 2019
pagenumber_ebook=169,pagenumber_book=716
雷达组网遮挡率、降水累计与雨量计数据如图11 所示，其中雨量计数据（图11c）与表3 中数据存在少量差异，源于表3 中数据为每日20 时至次日20 时累计结果，与香港天文台数据存在统计时间差异但总体差距很小不影响与雷达估测对比验证。少。两者之间具体数值存在一定差异，分析原因主要包括：雷达估测降水使用普适Z-R 关系 Z =300R1.4，不能完全真实反映降水滴谱分布情况；组网观测的探测高度不同会对降水估测产生误差；雷达基数据的质量控制算法也会影响降水估测质量。通过与雨量计数据对比，可以验证遮挡区域降水偏弱情况的存在，以及雷达组网相比于单站对遮挡区域有较大改善。
雷达单站及组网遮挡率与降水累计对比可发现，广东单站0.5°仰角受遮挡的区域通过相互填补基本不存在遮挡区，广东省雷达组网探测基本全省覆盖。对比验证所用数据未与周边省市雷达站点进行组网计算所以广东省外存在遮挡区域，广东省外探测盲区与11 部雷达组网遮挡率计算结果一致，连州西偏北区域存在遮挡盲区，25°N 以北的淡紫色弱雨区与遮挡率匹配一致，在前面全国拼图中广东省以外区域遮挡影响明显改善。雷达组网估测与雨量计数据进行对比，整体趋势基本一致，主要降水集中区域位于粤中部，东北部降水相对较
6 结 论
pagenumber_ebook=170,pagenumber_book=717
图 10 2019 年 8 月雷达单站 0.5°仰角遮挡率（a1、b1、c1）与估测降水累计（a2、b2、c2）分布对比（a. 韶关，b. 汕尾，c. 梅州）
Fig. 10 Comparison of blockage percentage at 0.5° elevation angle（a1，b1，c1） and accumulation of precipitation （a2，b2，c2）estimates between radar stations in August 2019 （a. Shaoguan，b. Shanwei，c. Meizhou）
pagenumber_ebook=171,pagenumber_book=718
图 11 广东省 2019 年 8 月雷达组网 0.5°仰角遮挡率（a）、估测降水累计（b）与雨量计数据（c）对比
Fig. 11 Comparison of blockage percentage at 0.5° elevation angle（a） and accumulation of precipitation estimates （b） of Guangdong radar network with rain gauge data （c） in August 2019
文中利用SRTM 数字高程模型90 m 分辨率高程数据对目前组网业务运行的212 部天气雷达进行单站和组网遮挡情况分析，统计VCP21 模式下低仰角层雷达覆盖情况，根据统计结果显示，低层5 层仰角（0.5°、1.5°、2.4°、3.4°、4.3°）探测范围覆盖面积分别为70.09%、90.15%、95.39%、97.60%、98.32%。对于中国东部、南部沿海易发生台风暴雨的区域能够全覆盖；对于内陆主要雨区具有较高的覆盖率；西部地区雷达布网较少，青藏高原等大地形区域覆盖效果较差，但上述地区为弱降水区，对于大范围的降水观测影响较小。通过广东省2019 年8 月的雨量计数据和11 部雷达的降水估测数据对比验证遮挡区域的计算准确性，遮挡区域计算结果与实际数据一致性较高，证明计算结果较为准确。雷达组网相比于单站受遮挡情况有较大改善，更接近雨量计数据，证实遮挡计算结果有利于为基数据遮挡订正提供数值参考依据。
综上所述，中国天气雷达组网探测环境效果较好，地形遮挡通过组网观测方式能够较好进行控制处理。通过统计单站相关数据，掌握雷达站探测环境情况，也可为新建站选址提供数据依据。雷达组网遮挡拼图可为后期多源数据融合时雷达回波强度处理及数据空间一致性提供订正依据。
雷达测高和波束遮挡计算在文中均采用标准大气模式进行计算，仅能代表历史数据平均情况下的计算结果，对于一般的地形遮挡情况误差在可接受范围内，如有更高精度需求应当考虑温度、湿度等对雷达波束传播路径的影响。由于计算量巨大，雷达探测距离波束展宽情况计算有待完善。
雷达波束遮挡计算对于高程数据的分辨率和更新及时性依赖性较高，随着城市发展新建高大建筑对于雷达探测环境也产生不利影响。因此，文中所使用的地形高程数据会造成在某些方位未能计算出相应的遮挡情况。获取更加准确的遮挡信息可通过与测绘部门合作叠加更高分辨率的地形数据进行更高精度的雷达波束遮挡计算，配合更长时间的基数据累积进行分析，以提高雷达探测回波强度计算的准确性。
参考文献
邓志，刘文文，黄美金. 1999. 雷达等射束高度图制作方法. 气象，25（11）：47-49. Deng Z，Liu W W，Huang M J. 1999. Drawing radar chart with ISO-beam height. Meteor Mon，25（11）：47-49 （in Chinese）
范思睿. 2016. 基于SRTM 数据的波束遮挡能量耗损率计算方法和效果分析. 高原山地气象研究，36（4）：79-85. Fan S R. 2016. Generation and verification of beam blockage coefficient based on SRTM data. Plateau Mountain Meteor Res，36（4）：79-85 （in Chinese）
勾亚彬，刘黎平，杨杰等. 2014. 基于雷达组网拼图的定量降水估测算法业务应用及效果评估. 气象学报，72（4）：731-748. Gou Y B，Liu L P，Yang J， et al. 2014. Operational application and evaluation of the quantitative precipitation estimates algorithm based on the multi-radar mosaic. Acta Meteor Sinica，72（4）：731-748 （in Chinese）
勾亚彬，汪章维，刘黎平等. 2017. 雷达波束部分遮挡识别应用及效果评估.高原气象，36（1）：229-240. Gou Y B，Wang Z W，Liu L P，et al. 2017.Application and evaluation of radar partial blockage identification.Plateau Meteor，36（1）：229-240 （in Chinese）
郭春辉，袁微，雷卫延. 2018. 肇庆天气雷达径向遮挡订正系统设计. 气象水文海洋仪器，35（1）：57-60. Guo C H，Yuan W，Lei W Y. 2018. Design on radial blocking correction system of Zhaoqing weather radar. Meteor，Hydrol Mar Instrum，35（1）：57-60 （in Chinese）
蒋秋菲，王旭，何建新等. 2017. 天气雷达遮挡区域识别. 成都信息工程大学学 报 ，32（5）：469-473. Jiang Q F，Wang X，He J X，et al. 2017.Identification of weather radar shielding region. J Chengdu Univ Inf Technol，32（5）：469-473 （in Chinese）
景号然，谢晓琳，郑伟等. 2019. 基于SRTM 数据的天气雷达探测环境分析研究. 气象，45（6）：871-876. Jing H R，Xie X L，Zheng W，et al. 2019.Research on weather radar environment detection based on SRTM elevation data. Meteor Mon，45（6）：871-876 （in Chinese）
李柏，古庆同，李瑞义等. 2013. 新一代天气雷达灾害性天气监测能力分析及未来发展. 气象，39（3）：265-280. Li B，Gu Q T，Li R Y，et al. 2013.Analyses on disastrous weather monitoring capability of CINRAD and future development. Meteor Mon，39（3）：265-280 （in Chinese）
石宝灵，王红艳，刘黎平. 2018. 云南多普勒天气雷达网探测冰雹的覆盖能力. 应用气象学报，29（3）：270-281. Shi B L，Wang H Y，Liu L P. 2018.Overage capacity of hail detection for Yunnan Doppler weather radar network. J Appl Meteor Sci，29（3）：270-281 （in Chinese）
王曙东，裴翀，郭志梅等. 2011. 基于SRTM 数据的中国新一代天气雷达覆盖和地形遮挡评估. 气候与环境研究，16（4）：459-468. Wang S D，Pei C，Guo Z M，et al. 2011. Evaluations on Chinese next generation radars coverage and terrain blockage based on SRTM data. Climatic Environ Res，16（4）：459-468 （in Chinese）
许小峰. 2003. 中国新一代多普勒天气雷达网的建设与技术应用. 中国工程科 学 ， 5（6）： 7-14. Xu X F. 2003. Construction， techniques and application of new generation Doppler weather radar network in China.Eng Sci，5（6）：7-14 （in Chinese）
俞小鼎，姚秀萍，熊廷南等. 2006. 多普勒天气雷达原理与业务应用. 北京：气象出版社，217-231. Yu X D，Yao X P，Xiong T N，et al. 2006. The Theory and Operational Application of Doppler Weather Radar. Beijing：China Meteorological Press，217-231 （in Chinese）
杨洪平，张沛源，程明虎等. 2009. 多普勒天气雷达组网拼图有效数据区域分析. 应用气象学报，20（1）：47-55. Yang H P，Zhang P Y，Cheng M H，et al. 2009. The valid mosaic data region of the CINRAD network. J Appl Meteor Sci，20（1）：47-55 （in Chinese）
张培昌，杜秉玉，戴铁丕. 2001. 雷达气象学. 北京：气象出版社，122-124.Zhang P C，Du B Y，Dai T P. 2001. Radar Meteorology. Beijing：China Meteorological Press，122-124 （in Chinese）
朱丹，谷军霞，师春香等. 2018. 新一代天气雷达布网设计的有效覆盖和地形遮挡分析. 气象，44（11）：1434-1444. Zhu D，Gu J X，Shi C X，et al.2018. Analysis of effective coverage and terrain blockage for next generation radar network design. Meteor Mon， 44（11）： 1434-1444 （in Chinese）
Dinku T，Anagnostou E N，Borga M. 2002. Improving radar-based estimation of rainfall over complex terrain. J Appl Meteor，41（12）：1163-1178
Farr T G，Rosen P A，Caro E，et al. 2007. The shuttle radar topography mission. Rev Geophys，45（2）：RG2004
Fulton R A，Breidenbach J P，Seo D J，et al. 1998. The WSR-88D rainfall algorithm. Wea Forecasting，13（2）：377-395
Gabella M. 2004. Improving operational measurement of precipitation using radar in mountainous terrain. Part Ⅱ：Verification and applications. IEEE Geosci Remote Sens Lett，1（2）：84-89
Qi Y C， Zhang J. 2017. A physically based two-dimensional seamless reflectivity mosaic for radar QPE in the MRMS system. J Hydrometeorol，18（5）：1327-1340
Tang L，Zhang J，Qi Y C，et al. 2013. Non-standard blockage mitigation for national radar QPE products∥Preprints，36th Conf. on Radar Meteorology. Breckenridge：Amer Meteor Soc，354
Westrick K J，Mass C F，Colle B A. 1999. The limitations of the WSR-88D Radar Network for quantitative precipitation measurement over the Coastal Western United States. Bull Amer Meteor Soc， 80（11）：2289-2298
Zhang J，Qi Y，Howard K，et al. 2012. Radar Quality Index （RQI）：A combined measure of beam blockage and VPR effects in a national network∥Proceedings，International Symposium on Weather Radar and Hydrology. Exeter：IAHS Publ.，351：388-393
Zhang J， Howard K， Langston C， et al. 2016. Multi-Radar Multi-Sensor（MRMS）quantitative precipitation estimation：Initial operating capabilities. Bull Amer Meteor Soc，97（4）：621-638
LIU Yunlei1 LI Changxing2 ZHANG Lejian2 WEN Hao2 WANG Shudong3
1. Dalian Meteorological Equipment Support Center，Dalian 116001，China
2. Meteorological Observation Centre of China Meteorological Administration，Beijing 100081，China
3. Public Meteorological Service Center of China Meteorological Administration，Beijing 100081，China
Abstract The weather radar beam energy is usually attenuation due to the terrain blockage,which will lead to lower intensity of radar reflectivity and unrealistic radar quantitative precipitation estimation (QPE). Therefore, it is an important basic work to statistics and analyze of radar beam blockage. This paper simulates and analyzes the observational environments of the 212 operational CINRAD radars using Shuttle Radar Topography Mission (SRTM) digital elevation data. The shielding angle, the Beam Blockage Percentages (BBPs) at 0.5°, 1.5°, 2.4°, 3.4° and 4.3° elevation angles in VCP21 mode, the BBPs in hybrid scan mode and subarea hybrid scan mode，and the effective coverage of observational range are first calculated for each weather radar. The terrain blockage effect on the CINRAD radar network is then analyzed statistically and the BBP mosaic image of the CINRAD radar network is mapped. The 11 CINRAD radars base data in Guangdong province in August 2019 are used to verify the analysis result of single station and radar network beam blockage. The results indicate that the effective observational coverage of the CINRAD radar network region exceeds 70%. The analysis results of BBPs are highly consistent with the actual data, and have positive contributions to radar reflectivity correction, QPE correction and other products.
Key words CINRAD，Beam blockage，Shielding angle，Radar network blockage
2020-01-02 收稿，2020-03-11 改回.
柳云雷，李昌兴，张乐坚，文浩，王曙东. 2020. 基于高分辨率高程数据统计分析新一代天气雷达组网的地形遮挡影响. 气象学报，78（4）：705-720
Liu Yunlei， Li Changxing， Zhang Lejian， Wen Hao， Wang Shudong. 2020. Statistical analysis of terrain blockage impacts on the CINRAD network based on DEM data. Acta Meteorologica Sinica， 78（4）:705-720
* 资助课题：“十三五”中国气象局新一代天气雷达探测工程项目。
作者简介：柳云雷，主要从事气象装备保障与雷达资料应用研究。E-mail：liuyunlei@tom.com
通信作者：李昌兴，主要从事多源数据融合与卫星、雷达资料应用研究。E-mail：licx@cma.gov.cn
中图法分类号 P412.25
页:
[1]