一种天气雷达强度资料质量测试方法
【技术领域】
[0001] 本发明涉及气象雷达数据处理技术领域，尤其是一种天气雷达强度资料质量测试 方法。
【背景技术】
[0002] 雷达作为一种主动遥感手段可得到具有一定精度的、大范围高时空分辨率的实时 降水信息，应用雷达进行降雨监测和面雨量计算，可以提高洪水预报的精度和时效性，在洪 灾监测预报中有很好的应用前景，因此雷达在水文水资源学中的应用研究日益受到关注。
[0003] 多年来许多气象和水文学者在水文水资源研究中应用雷达降雨信息做了大量的 研究工作，这些研究为提高水文预报精度提供了许多有借鉴意义的建设性思路，也取得了 大量的积极性的成果。这些工作主要集中在利用测雨雷达回波资料进行降水估算方法研 究，如卡尔曼滤波校准法、最优插值法和变分校准法等，以及将雷达资料与雨量站资料、卫 星云图资料等相结合，构建出不同的水文河流预报系统，为流域洪水预报中采用雷达测雨 资料提供了前瞻性的研究上，但对于雷达资料本身的质量测试却少有研究，而雷达资料质 量测试应该是雷达资料应用前最为重要的一环，直接关系到最终雨量估测结果的合理性。
[0004] 影响雷达资料质量的因子有很多，其中非气象回波如地物杂波，海浪杂波和生物 杂波的干扰，和降水粒子对电磁波的衰减是两个最为主要的因子。这些因子会使雷达测量 值与其真实值之间产生偏差，并直接影响最终的降水估测产品。尽管对雷达资料质量存在 的各种问题已研究出相应的识别及订正算法，但这些算法仍存在一定的局限性：即把重心 放在数据订正而非数据测试上。如果资料质量太差，就算订正方法再先进也很难得到正确 的结果，如果据此结果来计算雨量，最终的结果就会产生很大的误差，甚至会出现严重的错 误。
【发明内容】
[0005] 本发明的目的在于提供一种能够改善雷达测雨的可信度、提高暴洪监测的能力的 天气雷达强度资料质量测试方法。
[0006] 为实现上述目的，本发明采用了以下技术方案：一种天气雷达强度资料质量测试 方法，该方法包括下列顺序的步骤：
[0007] (1)将波束展宽因子对应雷达波束展宽质量指数Fra^将波束阻挡因子对应雷达 波束阻挡质量指数F shleld;将电磁波衰减因子对应雷达电磁波衰减质量指数Fatj#垂直廓 线不均匀因子对应雷达垂直廓线不均匀质量指数F胃；
[0008] (2)将雷达波束展宽质量指数Frange、雷达波束阻挡质量指数Fshleld、雷达电磁波衰 减质量指数F att、雷达垂直廓线不均匀质量指数F胃按按相应的权重系数进行加权求和得到 雷达平均质量指数Fz:
[0009]
[0010] 其中，Cz= W range+Wshield+Watt+Wvpr ' Wrange、WsJ1 ieid、Watt、WVpr 分别为雷达波束展宽质量指 数Frange、雷达波束阻挡质量指数Fshleld、雷达电磁波衰减质量指数F att、雷达垂直廓线不均匀 质量指数F胃的权重系数；当雷达波束被地物完全遮挡，即Fshleld= 0,或雷达电磁波受降水 粒子影响而存在严重衰减，即Fatt= 0时，其F z= 0。
[0011] 所述雷达波束展宽质量指数的计算公式如下：
[0012]
[0013] 其中，是雷达最大探测距离，r为某一点离雷达的距离，二者单位都是Km。对 于雷达电磁波受到部分阻挡时，其雷达波束阻挡质量指数F shleld由下式表示：
[0014] Fshleld= 0· 5tanh (41η ⑵[Θ。- θ α] / Θ )+〇· 5 (3)
[0015] 上式中，Θ。和θ α分别为波束主瓣轴的仰角和地物相对雷达的仰角，单位均 为°，?为3dB波束宽度，单位为° ;在完全阻挡时，Fshleld=O;而在完全无阻挡时，则Fshleld =1〇
[0016] 在双程衰减下，雷达电磁波衰减质量指数Fatt由下式表示：
[0017]
[0018] 其中：1(_和K _分别为衰减量的最大和最小阈值，单位为dB ;信号衰减越严重， Fatt越接近0,反之则越接近1。
[0019] 当雷达有效照射波束低于零度层高度时，雷达探测到的回波强度Zh主要来自液 态降水粒子，其雷达垂直廓线不均匀质量指数F胃=1 ;而零度层以上存在雪、冰雹、霰等粒 子，其雷达垂直廓线不均匀质量指数F胃=0。
[0020] 所述 Wrange = 1，W shleld = 1，W att = 1，W vpr = 1 〇
[0021] 所述 rmax=60。
[0022] 所述 Knax = 6，K ηιη = I 〇
[0023] 由上述技术方案可知，本发明的优点如下：第一，本发明将导致天气雷达强度资料 质量不确定性的主要因子一一波束展宽、波束遮挡、地物杂波、电磁波衰减等因素以距离库 为单元按〇~1的取值范围量化，处理成雷达资料质量指数，并将各雷达质量指数以指定的 权重系数加权求和得到相应的平均雷达质量指数，从而直观地了解雷达资料质量受各因子 的影响程度及其整体情况，以此对雷达资料质量进行定量测试；第二，本发明充分考虑了影 响雷达强度资料质量的主要因子，方法合理，得到的测试结果直观，方法适用性强。
【附图说明】
[0024] 图1为本发明的方法流程图；
[0025] 图2、3均为本发明的算法中各参数的取值范围示意图。
【具体实施方式】
[0026] 如图1所示，一种天气雷达强度资料质量测试方法，该方法包括下列顺序的步骤：
[0027] (1)将波束展宽因子对应雷达波束展宽质量指数Fra^将波束阻挡因子对应雷达 波束阻挡质量指数F shleld;将电磁波衰减因子对应雷达电磁波衰减质量指数Fatj#垂直廓 线不均匀因子对应雷达垂直廓线不均匀质量指数F胃；
[0028] (2)将雷达波束展宽质量指数Frange、雷达波束阻挡质量指数Fshleld、雷达电磁波衰 减质量指数F att、雷达垂直廓线不均匀质量指数F胃按按相应的权重系数进行加权求和得到 雷达平均质量指数Fz:
[0029]
[0030] 其中，Cz= W range+Wshield+Watt+Wvpr ' Wrange、WsJ1 ieid、Watt、WVpr 分别为雷达波束展宽质量指 数Frange、雷达波束阻挡质量指数Fshleld、雷达电磁波衰减质量指数F att、雷达垂直廓线不均匀 质量指数F胃的权重系数；当雷达波束被地物完全遮挡，即Fshleld= 0,或雷达电磁波受降水 粒子影响而存在严重衰减，即Fatt= 0时，其Fz= 0。所述Wrangf3= 1，Wshif3ld= 1