[发明专利]X波段双偏振天气雷达定量测量降水方法及系统有效
申请号: | 201910906092.2 | 申请日: | 2019-09-24 |
公开(公告)号: | CN110488393B | 公开(公告)日: | 2020-05-15 |
发明(设计)人: | 冯亮;肖辉;孙跃;李宗飞;罗丽 | 申请(专利权)人: | 中国科学院大气物理研究所 |
主分类号: | G01W1/14 | 分类号: | G01W1/14;G01S13/95;G01S7/41 |
代理公司: | 北京路浩知识产权代理有限公司 11002 | 代理人: | 张睿 |
地址: | 100029 北京*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 波段 偏振 天气 雷达 定量 测量 降水 方法 系统 | ||
1.一种X波段双偏振天气雷达定量测量降水方法,其特征在于,包括:
利用X波段双偏振雷达获取探测范围内的雷达数据,并利用激光雨滴谱仪获取所述探测范围内任一检测点的实测雨滴谱;
对所述雷达数据进行衰减订正,获取所述X波段双偏振天气雷达在所述探测范围内的雷达反射率因子和衰减率;
基于所述实测雨滴谱,获取所述检测点的实测降雨率和雷达反射率因子真值,并通过T矩阵散射模拟,获取所述检测点的衰减率真值;通过拟合所述实测降雨率与所述衰减率真值,获取所述探测范围内的降水函数;
若所述雷达反射率因子大于第一设定阈值,则直接利用所述衰减率以及所述降水函数反演降雨率;
若所述雷达反射率因子小于或等于第一设定阈值,则先通过拟合所述雷达反射率因子真值与所述衰减率真值,获取所述探测范围内的衰减率函数;再基于所述衰减率函数以及所述降水函数反演降雨率;
所述降水函数为:R=aAHb,所述衰减率函数为:AH=cZHd,其中R为降雨率,AH为衰减率,ZH为雷达反射率因子,a、b、c和d均为拟合系数。
2.根据权利要求1所述的X波段双偏振天气雷达定量测量降水方法,其特征在于,在所述对所述雷达数据进行衰减订正之前,还包括:
通过反演信噪比对所述雷达数据进行相关系数校准,若所述信噪比大于第二设定阈值时,依次对所述雷达数据进行地物杂波滤除、数据短缺点补齐、掉线径向数据补缺、差分传播相移退折叠以及二维中值滤波处理。
3.根据权利要求1所述的X波段双偏振天气雷达定量测量降水方法,其特征在于,所述利用激光雨滴谱仪获取所述探测范围内任一检测点的实测雨滴谱,其计算方法为:
其中,N(D)为实测雨滴谱的计算函数,Δt为采样时间间隔,i为雨滴直径分档标签,j为雨滴标签,mi为Δt时间内雨滴直径属于第i档的雨滴数量,D为雨滴等效体积直径,ΔD为雨滴直径分辨率,Aj为雨滴有效测量面积,vj为雨滴下落末速度。
4.根据权利要求1所述的X波段双偏振天气雷达定量测量降水方法,其特征在于,基于所述实测雨滴谱,获取所述检测点的实测降雨率,所述实测降雨率的计算公式为:
其中,Δt为采样时间间隔,i为雨滴数目标签,n为Δt时间内采样的雨滴总数目,Vi为雨滴的体积,Ai为雨滴有效测量面积。
5.根据权利要求1所述的X波段双偏振天气雷达定量测量降水方法,其特征在于,基于所述实测雨滴谱,获取所述检测点的雷达反射率因子真值,所述雷达反射率因子真值的计算公式为:
ZH=10log10(z)
其中,d表示微分,z为中间变量;
所述通过T矩阵散射模拟,获取所述检测点的衰减率真值,所述衰减率真值的计算公式为:
其中,AH为单程衰减率真值,Im为取积分后的虚部,fH(D)为水平前向散射振幅矩阵,λ为雷达波的波长;D为雨滴等效体积直径。
6.根据权利要求1所述的X波段双偏振天气雷达定量测量降水方法,其特征在于,获取所述探测范围内的降水函数,包括:
若所述雷达反射率因子大于第一设定阈值,则基于Matlab中的Curve Fitting Tool对所述实测降雨率与所述衰减率真值进行拟合,以确定所述拟合系数a和b,并获取所述降水函数;
若所述雷达反射率因子小于或等于第一设定阈值,则基于Matlab中的Curve FittingTool对所述反射率因子真值与所述衰减率真值进行拟合,获取衰减率函数AH=cZHd中的所述拟合系数c和d;并根据获取的衰减率和所述降水函数反演降雨率。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国科学院大气物理研究所,未经中国科学院大气物理研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910906092.2/1.html,转载请声明来源钻瓜专利网。