[发明专利]基于稠密雨量站定量降水估测方法

说明书

本发明提出了一种基于稠密雨量站定量降水估测方法，将一个选定的雷达测雨区域划分为很多小的多边形区域，每个区域使用其各自的Z‑R关系，当不同区域为不同降水性质的时候，避免了使用同一Z‑R关系所带来的误差；当雷达数据质量出现错误时，比如波速阻挡导致的雷达反射率因子偏小等等情况，该方法可以根据实时资料，通过对Z‑R关系的调整，对由于雷达质量造成的降水估算偏差给予一定的订正。

技术领域

本发明涉及气象雷达估测降水领域，尤其涉及一种基于稠密雨量站定量降水估测方法。

背景技术

雷达定量估测降水是短时临近预警的一个重要技术参考，结合雷达定量估算降水(QPE)，出现了具有各种不同思路和技术方法的研究成果，诸如变分法、卡尔曼滤波法以及概率配对法等。

目前，湖北省雷达定量估测降水的主要技术是以雨量计为代表的点测量和以天气雷达为代表的面测量结合起来，采用雷达和雨量计实时同步积分结合方法(RASIM)。该方法在理论探讨、降水估测技术、质量控制、误差分析与评估标准等方面与国内外其他方法相比具有先进性和实用性，并且该方法使用方便，易于推广，具有较好的正确率，已经成为我国短时临近预报系统(SWAN)中的重要功能模块，其广泛应用，效果良好。但是该技术在一个雷达测量范围内仅使用一个Z-R关系，不能准确反映降水局地性分布的变化特征。

发明内容

有鉴于此，本发明提出了一种能准确反映降水局地性分布的基于稠密雨量站定量降水估测方法。

本发明的技术方案是这样实现的：本发明提供了一种基于稠密雨量站定量降水估测方法，包括以下步骤，

S1，对于一个选定的雷达测雨区域，根据该雷达测雨区域内覆盖的雨量站的分布，利用泰森多边形法对该雷达测雨区域其进行划分，得到多个多边形区域；

S2，计算每个多边形区域的Z-R关系；

S3，根据每个多边形区域的Z-R关系，计算选定的雷达测雨区域内的降雨强度R。

在以上技术方案的基础上，优选的，所述步骤S2中，对于选定的多边形区域S，其覆盖的雨量站有i个，计算该多边形区域S内在h时间段内i个雨量站的降雨量总和Q_Sh，根据以下公式计算该多边形区域的Z-R关系，即A_Sh的值，

其中，Z为雷达反射率因子，R为降雨强度，A为Z-R关系的系数，A＝10-2000,b为Z-R关系的指数，b＝1.0-2.0；

A_Sh为选定的多边形区域S的公共系数；

Z_Sh为选定的多边形区域S的雷达总反射率因子；

Z_ih为选定的多边形区域S内的雨量站i上空在h时间段内对应的0.5°仰角的雷达反射率因子。

进一步优选的，b＝1.5。

进一步优选的，步骤S3中，根据以下公式计算选定的多边形区域S内的降雨强度R，

Z_jh为选定的多边形区域S内的网格点j上空在h时间段内对应的0.5°仰角的雷达反射率因子。

更优选的，步骤S3中，根据雷达测雨多边形区域的公共系数Ash和每个网格点上空的0.5度仰角的反射率因子，计算所有多边形区域内的降雨强度R，得到雷达测雨区域多个多边形区域的降雨强度R。

再优选的，根据雷达测雨区域多个多边形区域的降雨强度R，绘制雷达测雨区域1km分辨率的降雨强度R按多个多边形区域的分布图。