CN113740934A - 一种基于s波段双偏振天气雷达降水估测方法 - Google Patents

一种基于s波段双偏振天气雷达降水估测方法 Download PDF

Info

Publication number
CN113740934A
CN113740934A CN202110949908.7A CN202110949908A CN113740934A CN 113740934 A CN113740934 A CN 113740934A CN 202110949908 A CN202110949908 A CN 202110949908A CN 113740934 A CN113740934 A CN 113740934A
Authority
CN
China
Prior art keywords
precipitation
radar
rainfall
data
polarization
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110949908.7A
Other languages
English (en)
Other versions
CN113740934B (zh
Inventor
汪章维
陈昊
王晗
卢兴来
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Atmospheric Measurement Technology Security Center
Original Assignee
Zhejiang Atmospheric Measurement Technology Security Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang Atmospheric Measurement Technology Security Center filed Critical Zhejiang Atmospheric Measurement Technology Security Center
Priority to CN202110949908.7A priority Critical patent/CN113740934B/zh
Publication of CN113740934A publication Critical patent/CN113740934A/zh
Application granted granted Critical
Publication of CN113740934B publication Critical patent/CN113740934B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/14—Rainfall or precipitation gauges
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Atmospheric Sciences (AREA)
  • Environmental Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Ecology (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Hydrology & Water Resources (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于S波段双偏振天气雷达降水估测方法，涉及降水估测领域。目前基于双偏振雷达的降水估测方法，主要建立基于偏振量的不同组合与降水率之间的关系方法，使用小时雨量计进行订正，时间尺度较大，影响最后的订正效果。本方法使用雷达有效覆盖范围内的实时分钟雨量计数据，并且在进行液态降水估测时采用双阈值，选择不同的降水估测算法计算，对双偏振雷达的逐体扫降水率产品进行实时校准订正。其订正时间尺度更小，订正后的逐体扫雨强产品的精度更高，使用订正后的雨强产品累计得到的小时降水量产品更准确，且长期稳定性更好。

Description

一种基于S波段双偏振天气雷达降水估测方法
技术领域
本发明涉及降水估测领域，尤其涉及一种基于S波段双偏振天气雷达降水估测方法。
背景技术
天气雷达是大气探测、云降水物理研究的重要手段，被广泛应用于气象监测预警和气象科学研究。天气雷达不仅可以为中小尺度天气过程提供准确的观测资料，还能够通过降水产品反映空间降水的不均匀性。雷达定量降水估测业务是现代天气业务体系的重要组成部分，也是定量降水预报、强降水短临预警的基础，对山洪地质灾害和中小河流洪水精细化预报起重要支撑作用，因此提高雷达降水估测精度具有十分重要的意义。
随着双偏振多普勒天气雷达的快速发展及在云降水领域的广泛应用，双偏振雷达的多偏振量极大提高了降水估测及水凝物分类识别的能力。国内外学者对于S波段双偏振天气雷达在降水估测方面的算法做了大量研究，目前基于双偏振雷达的降水估测方法，主要建立基于偏振量(ZH、ZDR、KDP)的不同组合与降水率之间的关系方法，即选择组合关系法，如JPOLE算法，CSU-ICE算法，CSU-HIDRO算法等。这些主流方法的思路大致分为三类，第一类是基于反射率与降水率之间的简单关系(R(ZH))计算出降水率强度来选择不同的降水估测关系，即采用R(ZH)关系对降水强度进行简单分型，然后将降水强度划分为几个不同的强度范围，如小雨(R<2.5mm/h)，中雨(2.5mm/h<R<8.0mm/h)，大雨(8mm/h<R<16mm/h)，暴雨(R>16mm/h)等降水类型，最后，根据粗算的降水率R所对应的同降水强度类型，采用不同的降水率与偏振量的组合关系进行二次精确计算。第二类是根据偏振量(ZH、ZDRv、KDP)阈值来选择不同参数组合计算降水强度。第三类算法将水凝物的相态识别加入到降水估测中，首先使用模糊逻辑法对水凝物进行分类(如：地物，冰晶，小雨，大雨，大雨滴等)，再根据ZH、ZDR和KDP的阈值来选择降水估测的关系式。
现有S波段双偏振天气雷达的降水估测方法主要有以下几个缺点：第一：使用提前划分好的降水强度类型选择不同的偏振量与降水率之间的关系组合的做法太过于简单化，估测精度会大大受到粗略计算的降水强度的影响，如冰晶粒子粗略计算的降水率会偏大，最终精确计算得到的降水率较高，但实际却没有液态降水。第二，降水率与偏振量之间的组合方法中大量使用了偏振量差分反射率因子(ZDR)。但是差分反射率因子由于存在标定不及时将出现较大的偏差，在弱降水时对降水估测结果产生较大的影响，使得计算结果不准确。第三，目前双偏振雷达降水估测方法由于受到观测参量标定等影响，其估测精度的稳定性较难保证，难以提供长期稳定、准确的降水场估测资料。第四：目前也存在使用雨量计对雷达降水估测进行校准订正的方法，但是基本都是使用小时雨量计进行订正，时间尺度较大，影响最后的订正效果。
发明内容
本发明要解决的技术问题和提出的技术任务是对现有技术方案进行完善与改进，提供一种基于S波段双偏振天气雷达降水估测方法，以提升S波段双偏振天气雷达对降水强度的估计准确度为目的。为此，本发明采取以下技术方案。
一种基于S波段双偏振天气雷达降水估测方法，包括以下步骤：
1)获取S波段双偏振雷达体扫基数据；
2)对S波段双偏振雷达体扫基数据进行质量控制；
3)使用基于模糊逻辑的双偏振雷达定量降水估测算法计算雷达降水估测数据；
4)获取S波段雷达覆盖范围内自动站分钟级雨量计数据；
5)采用分钟级雨量计数据对雷达降水估测数据进行逐体扫实时订正，计算逐体扫实时校准订正后的降水估测数据，得到校准订正后雷达逐体扫瞬时降水强度产品QPE；
6)通过逐体扫实时校准订正后的降水估测数据，累计后得到双偏振雷达小时降水量估测数据，即双偏振雷达小时降水量产品OHP。
使用雷达有效覆盖范围内的实时分钟级雨量计数据，对双偏振雷达的逐体扫降水率产品进行实时校准订正，校准订正的数据时间分辨率更高，其订正时间尺度更小，从本质上改善小时降水量的系统性偏差问题，订正后的逐体扫降水强度的精度更高，使用订正后的降水强度累计得到的小时降水量更准确，且长期稳定性更好，极大的改善降水估测的精度，提升了S波段双偏振天气雷达对降水强度的估计准确度。
而且，相比单偏振基于ZR关系进行降水率计算的方法，双偏振量能提供更多的降水微物理信息，并且差分相移率KDP受到雨区衰减和波束遮挡的影响较小，与强降水的线性关系较好。因此，基于双偏振量的降水反演能够获得更准确的降水估计结果。
作为优选技术手段：步骤2)中，数据质量控制包括：
201)反射率ZH使用中值滤波器以及沿着径向方向上采用5个距离库的滑动平均值的方法来消除异常值和减少随机波动；
202)根据相关系数ρhv小于0.85的阈值来剔除包括地物杂波、湍流散射、异常杂波的非气象回波；
203)对差分相移率KDP进行质量控制，首先根据地物杂波的差分相位均值来标定系统的初始相位，然后平滑滤波，最后通过线性编程的方法来计算KDP，线性规划的方法可以确保KDP的非负性，同时符合其物理含义。
S波段双偏振雷达数据质量控制主要考虑以下三个方面对数据带来的影响：(1)非气象回波(地物杂波、湍流散射、异常干扰)；(2)差分相移率(KDP)受到前向散射和后向散射所造成的噪点对降水估测的影响；(3)反射率的噪点对降水估测的影响，通过数据的质量控制，可有效降低以上问题可能对降水估测结果带来影响的数据情况。
作为优选技术手段：步骤4)中，自动站分钟级雨量计数据采用雷达覆盖150km范围内的实时分钟级雨量计数据。数据覆盖范围可靠性高，实时校准订正，其订正时间尺度小。
作为优选技术手段：步骤3)中，计算雷达降水估测数据时，以模糊逻辑水凝物分类为基础，将降水类型大致分为液态和非液态，对其中的反射率ZH和差分相移率KDP进行判定，然后根据KDP和ZH的阈值选择对应的降水率计算公式R(ZH)＝0.019ZH 0.761和R(KDP)＝44.84KDP 0.763进行计算。基于模糊逻辑的水凝物分类方法，根据其结果只对其中液态降水类型进行定量降水估测，并且在进行液态降水估测时采用双阈值法，即只对其中的反射率ZH和差分相移率KDP进行判定，选择不同的降水估测算法，可有效实现降水估测；使用模糊逻辑水汽分类且使用较为简单的结果(固态和液态)，并且在降水率计算的过程中主要使用R(ZH)和R(KDP)，减少了数据质量对降水率计算的影响。
作为优选技术手段：步骤5)中，采用分钟级雨量计数据对雷达降水估测数据进行逐体扫实时订正包括以下步骤：
501)寻找到雨量计在体扫时间内(雷达基数据体扫结束减去体扫开始时间)的累计雨量并转换为其对应的小时降水强度Rg，同时根据雨量计的经纬度信息寻找其垂直方向上降水强度对应距离库的估测雨量Rr；
502)使用降水强度Rg和估测雨量Rr二者的比值Rg/Rr作为订正因子，将订正因子通过基于薄板样条TPS算法插值到整个雷达降水估测范围内的每个格点，即形成了雷达降水强度中每个格点所对应的订正因子；
503)通过将每个格点中的订正因子乘以对应的降水强度的对应格点，计算逐体扫的校准订正实时降水估测数据，得到校准订正后雷达逐体扫瞬时降水强度产品QPE。从而有效实现降水强度的订正。
有益效果：使用雷达有效覆盖范围内的实时分钟级雨量计数据，对双偏振雷达的逐体扫降水率产品进行实时校准订正，校准订正的数据时间分辨率更高，其订正时间尺度更小，使用模糊逻辑水汽分类且使用较为简单的结果(固态和液态)，并且在降水率计算的过程中主要使用R(ZH)和R(KDP)，减少了数据质量对降水率计算的影响，从本质上改善小时降水量的系统性偏差问题，订正后的逐体扫降水强度的精度更高，使用订正后的降水强度累计得到的小时降水量更准确，且长期稳定性更好，极大的改善降水估测的精度，提升了S波段双偏振天气雷达对降水强度的估计准确度，更符合业务系统的需求。
附图说明
图1是本发明方法过程示意图。
图2是本发明计算逐体扫实时降水强度过程示意图。
图3是本发明实例中选取雷达站和地面站雨量计分布图。
具体实施方式
以下结合说明书附图对本发明的技术方案做进一步的详细说明。
如图1所示，一种基于S波段双偏振天气雷达降水估测方法，其过程包括以下步骤：
S1:获取S波段双偏振雷达体扫基数据。
S2:对S波段双偏振雷达体扫基数据进行质量控制。
数据质量控制包括：
S201:反射率ZH使用中值滤波器以及沿着径向方向上采用5个距离库的滑动平均值的方法来消除异常值和减少随机波动。
S202:根据相关系数(ρhv)小于0.85的阈值来剔除包括地物杂波、湍流散射、异常杂波的非气象回波。
S203:对差分相移率(KDP)进行质量控制，首先根据地物杂波的差分相位均值来标定系统的初始相位，然后平滑滤波，最后通过线性编程的方法来计算KDP，线性规划的方法可以确保KDP的非负性，同时符合其物理含义。
S波段双偏振雷达数据质量控制主要考虑以下三个方面对数据带来的影响：(1)非气象回波(地物杂波、湍流散射、异常干扰)；(2)差分相移率(KDP)受到前向散射和后向散射所造成的噪点对降水估测的影响；(3)反射率的噪点对降水估测的影响，通过数据的质量控制，可有效降低以上问题可能对降水估测结果带来影响的数据情况。
S3:使用基于模糊逻辑的双偏振雷达定量降水估测算法计算雷达降水估测数据；
为了有效实现降水估测，如图2所示，计算雷达降水估测数据时，以模糊逻辑水凝物分类为基础，将降水类型大致分为液态和非液态，对其中的反射率(ZH)和差分相移率(KDP)进行判定，然后根据KDP和ZH的阈值选择不同的降水率计算公式R(ZH)＝0.019ZH 0.761和R(KDP)＝44.84KDP 0.763进行计算。基于模糊逻辑的水凝物分类方法，根据其结果只对其中液态降水类型进行定量降水估测，并且在进行液态降水估测时采用双阈值法，即只对其中的反射率(ZH)和差分相移率(KDP)进行判定，选择不同的降水估测算法，可有效实现降水估测；使用模糊逻辑水汽分类且使用较为简单的结果(固态和液态)，并且在降水率计算的过程中主要使用R(ZH)和R(KDP)，减少了数据质量对降水率计算的影响。
S4:获取S波段雷达覆盖范围内自动站分钟级雨量计数据。
自动站分钟级雨量计数据采用雷达覆盖150km范围内的实时分钟级雨量计数据。数据覆盖范围可靠性高，实时校准订正，其订正时间尺度小。
S5:采用分钟级雨量计数据对雷达降水估测数据进行逐体扫实时订正，计算逐体扫实时校准订正后的降水估测数据，得到校准订正后雷达逐体扫瞬时降水强度产品QPE。
其包括以下步骤：
S501:寻找到雨量计在体扫时间内(雷达基数据体扫结束减去体扫开始时间)的累计雨量并转换为其对应的小时降水强度Rg，同时根据雨量计的经纬度信息寻找其垂直方向上降水强度对应距离库的估测雨量Rr；
S502:使用降水强度Rg和估测雨量Rr二者的比值Rg/Rr作为订正因子，将订正因子通过基于薄板样条TPS算法插值到整个雷达降水估测范围内的每个格点，即形成了雷达降水强度中每个格点所对应的订正因子；
S503:通过将每个格点中的订正因子乘以对应的降水强度的对应格点，计算逐体扫的校准订正实时降水估测数据，得到校准订正后雷达逐体扫瞬时降水强度产品QPE。有效实现降水强度的订正。
基于薄板样条TPS算法是一种数值插值方法，其原理为：在空间区域R2内分布n个已知点集，Pi，i＝1,2,…,n.其坐标表示为：(xi,yi,z(xi,yi))，如果z(xi,yi)具有二次连续导数，那么一个通过所有给定的点的“最小弯曲”光滑表面。弯曲能量可以表示为：
Figure BDA0003218285610000091
当该弯曲能量为最小时，所得到的插值函数就是TPS插值函数。此时，插值函数为：
Figure BDA0003218285610000092
所给出的点集为3D点坐标，分别用Pi和vi表示，其中Pi表示已给出的n个点中的第i个，U(r)为的核函数：
U(r)＝r2 log(r2) (3)
|Pi-(x，y)|表示目标点(x，y)点到已知点Pi的距离。
(2)式中各项的矩阵表达式如下(4)、(5)、(6)、(7)和(8)所示，最终式(9)为式(2)的矩阵表达式，其中，rij表示点Pi到点Pj的距离：
Figure BDA0003218285610000093
Figure BDA0003218285610000094
Figure BDA0003218285610000095
Figure BDA0003218285610000096
Figure BDA0003218285610000101
Figure BDA0003218285610000103
Figure BDA0003218285610000102
至此，根据表达式(10)解出的a和w构建插值函数f(x，y)。弯曲能量也可以由矩阵表示：
If＝wTKw (11)
现在就可以根据a和w，对任意点(x，y)进行v插值。
S6：通过逐体扫实时校准订正后的降水估测数据，累计得到双偏振雷达小时降水量估测。
为了保证本方法的有效性，还可包括步骤：
定时对小时降水量预测数据进行精度评估，若精度低于设定阈值时，对订正时参数进行修正；精度评估采用多个评估指标，包括：相对误差(RE)、均方根误差(RMSE)和皮尔逊相关系数(CC)；其中：相对误差(RE)能够较好的反映数据的可信度，RE越小，可信度越高；RMSE表示均方根误差，均方根误差能够很好地反映雷达估测值和雨量计实测值之间的离散程度，其反映了估测值偏离实测值的程度，其值越小离散程度越低，表明数据越集中，算法的稳定性越好；CC表征了估测值样本与实测值样本之间的相关性，相关系数越接近1，表明估测值与实测值之间的相关性越好；
三个指标具体定义如下：
相对误差(RE)：
Figure BDA0003218285610000111
均方根误差(RMSE)：
Figure BDA0003218285610000112
皮尔逊相关系数(CC)：
Figure BDA0003218285610000113
式中n为雷达-雨量计匹配数据点个数，Rgage,i表示第i组雨量计小时降水量，Rradar,i表示第i组雨量计对应雷达估测小时降水量，
Figure BDA0003218285610000114
表示n个雨量计实测小时降水量均值，
Figure BDA0003218285610000115
表示同n个雨量计对应雷达小时降水量估测均值。
以下实例内容为进一步验证本方法的优势。
一、获取数据：本技术方案选取杭州地区2020年汛期内多次大范围降水过程计算验证，以9月17日一次强降水过程数据为代表进行分析。该数据主要包括杭州S波段双偏振天气雷达观测资料和地面自动站雨量计资料。
地面站资料选取了在本次降水过程中雷达覆盖范围内约1000个有效降水的雨量计数据参与算法订正。其中，选取50多个雨量计用于评估，其不参与订正过程。所选取的雷达和雨量计分布情用来评估分析况如图3中圈内所示，降水信息统计表如表1所示。
表1降水过程信息统计
Figure BDA0003218285610000121
二、数据处理：
通过实时QPE(quantitative precipitation estimation)校准订正实时降水强度(简称QPE-ADJUST)，通过自动获取雷达覆盖150km范围内的分钟级雨量计资料，采用薄板样条(Thin Plate Splines，TPS)插值法对天气雷达的降水场资料进行校准，将雷达单个体扫时间内的定量降水估测产品QPE与该体扫时段内的分钟级雨量计资料进行降水强度比对、订正，修正雷达单位体扫时间内的降水强度。最后，通过时间累积完成1小时降水、3小时降水等降水估测产品。
三、评估分析：
采用多种评估指标，分别为相对误差(RE)、均方根误差(RMSE)和皮尔逊相关系数(CC)，分别对不同降水估测算法进行验证分析。
分别使用不同的方法：ZR关系法(基于单偏振ZR关系的降水估测方法)、CSU-HIDRO法(基于偏振量的雷达降水估测方法)、OHP-ADJUST法(基于小时雨量计订正的双偏振降水估测法)和本申请的方法：QPE-ADJUST(基于分钟雨量计资料订正的双偏振降水估测法)对9月17号的雷达资料进行处理，估测方法得到2020年9月17日5点至6点的小时降水量估测(OHP)结果。表2给出了本次降水过程的各方法的误差统计数据。其结果表明，本发明的方法明显优于其他方法，这说明雨量计的订正起到了正贡献。本技术方案的雨量计订正采用逐体扫订正和逐小时订正，其降水估计精度最好。
表2不同降水估计算法误差统计结果
Figure BDA0003218285610000131
因此，本方法能够在保证降水估计准确性的同时，还能够提高实时降水强度估计的准确性。
在本技术方案中，定时对小时降水量预测数据进行精度评估时，设定阈值为30％，由于其满足要求，故不对订正参数进行修正。
以上图1-3所示的一种基于S波段双偏振天气雷达降水估测方法是本发明的具体实施例，已经体现出本发明实质性特点和进步，可根据实际的使用需要，在本发明的启示下，对其进行形状、结构等方面的等同修改，均在本方案的保护范围之列。

Claims (5)

1.一种基于S波段双偏振天气雷达降水估测方法，其特征在于包括以下步骤：
1)获取S波段双偏振雷达体扫基数据；
2)对S波段双偏振雷达体扫基数据进行质量控制；
3)使用基于模糊逻辑的双偏振雷达定量降水估测算法计算雷达降水估测数据；
4)获取S波段雷达覆盖范围内自动站分钟级雨量计数据；
5)采用分钟级雨量计数据对雷达降水估测数据进行逐体扫实时订正，计算逐体扫实时校准订正后的降水估测数据，得到校准订正后雷达逐体扫瞬时降水强度产品QPE；
6)通过逐体扫实时校准订正后的降水估测数据，累计后得到双偏振雷达小时降水量估测数据，即双偏振雷达小时降水量产品OHP。
2.根据权利要求1所述的一种基于S波段双偏振天气雷达降水估测方法，其特征在于：步骤2)中，数据质量控制包括：
201)反射率ZH使用中值滤波器以及沿着径向方向上采用5个距离库的滑动平均值的方法来消除异常值和减少随机波动；
202)根据相关系数ρhv小于0.85的阈值来剔除包括地物杂波、湍流散射、异常杂波的非气象回波；
203)对差分相移率KDP进行质量控制，首先根据地物杂波的差分相位均值来标定系统的初始相位，然后平滑滤波，最后通过线性编程的方法来计算KDP，线性规划的方法可以确保KDP的非负性，同时符合其物理含义。
3.根据权利要求1所述的一种基于S波段双偏振天气雷达降水估测方法，其特征在于：步骤4)中，自动站分钟级雨量计数据采用雷达覆盖150km范围内的实时分钟级雨量计数据。
4.根据权利要求1所述的一种基于S波段双偏振天气雷达降水估测方法，其特征在于：步骤3)中，计算雷达降水估测数据时，以模糊逻辑水凝物分类为基础，将降水类型大致分为液态和非液态，对其中的反射率ZH和差分相移率KDP进行判定，然后根据KDP和ZH的阈值选择对应的降水率计算公式R(ZH)＝0.019ZH 0.761和R(KDP)＝44.84KDP 0.763进行计算。
5.根据权利要求1所述的一种基于S波段双偏振天气雷达降水估测方法，其特征在于：步骤5)中，采用分钟级雨量计数据对雷达降水估测数据进行逐体扫实时订正包括以下步骤：
501)寻找到雨量计在体扫时间内的累计雨量并转换为其对应的小时降水强度Rg，同时根据雨量计的经纬度信息寻找其垂直方向上降水强度对应距离库的估测雨量Rr；
502)使用降水强度Rg和估测雨量Rr二者的比值Rg/Rr作为订正因子，将订正因子通过基于薄板样条TPS算法插值到整个雷达降水估测范围内的每个格点，即形成了雷达降水强度中每个格点所对应的订正因子；
503)通过将每个格点中的订正因子乘以对应的降水强度的对应格点，计算逐体扫的校准订正实时降水估测数据，得到校准订正后雷达逐体扫瞬时降水强度产品QPE。
CN202110949908.7A 2021-08-18 2021-08-18 一种基于s波段双偏振天气雷达降水估测方法 Active CN113740934B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110949908.7A CN113740934B (zh) 2021-08-18 2021-08-18 一种基于s波段双偏振天气雷达降水估测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110949908.7A CN113740934B (zh) 2021-08-18 2021-08-18 一种基于s波段双偏振天气雷达降水估测方法

Publications (2)

Publication Number Publication Date
CN113740934A true CN113740934A (zh) 2021-12-03
CN113740934B CN113740934B (zh) 2022-12-09

Family

ID=78731632

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110949908.7A Active CN113740934B (zh) 2021-08-18 2021-08-18 一种基于s波段双偏振天气雷达降水估测方法

Country Status (1)

Country Link
CN (1) CN113740934B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115656947A (zh) * 2022-12-15 2023-01-31 南京信大气象科学技术研究院有限公司 一种双偏振天气雷达zdr偏差实时校正方法和装置
CN116381692A (zh) * 2023-04-18 2023-07-04 南京市气象台 一种基于x波段双偏振雷达的降水相态识别qpe算法
CN116303368B (zh) * 2023-04-24 2023-07-21 中国人民解放军国防科技大学 一种双偏振雷达体扫数据插值方法、装置、设备及介质

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002350560A (ja) * 2001-05-30 2002-12-04 Foundation Of River & Basin Integrated Communications Japan レーダ雨量補正・配信システム
FR2827391A1 (fr) * 2001-07-11 2003-01-17 Centre Nat Rech Scient Technique pour l'estimation de pluie avec radar lorsque l'attenuation est negligeable
KR101431707B1 (ko) * 2013-11-26 2014-09-22 한국건설기술연구원 엑스밴드 이중편파 레이더 관측자료를 이용한 통합형 강우량 산정 방법
CN105548984A (zh) * 2015-12-10 2016-05-04 中国人民解放军理工大学 一种基于模糊逻辑的双偏振多普勒天气雷达地物杂波抑制方法
CN106950614A (zh) * 2017-02-28 2017-07-14 中船重工鹏力（南京）大气海洋信息系统有限公司 一种区域自动气象站小时雨量数据质量控制方法
CN107015229A (zh) * 2017-05-22 2017-08-04 新疆维吾尔自治区人工影响天气办公室 基于双偏振天气雷达的人工防雹作业指挥方法
CN108680920A (zh) * 2018-04-28 2018-10-19 深圳市雅码科技有限公司 一种基于双偏振雷达的灾害天气识别预警系统及方法
CN110018448A (zh) * 2019-03-30 2019-07-16 中国人民解放军陆军工程大学 基于任意取向微波链路的双偏振雷达衰减订正方法
CN110082842A (zh) * 2019-05-24 2019-08-02 北京敏视达雷达有限公司 一种降水估测方法及装置
CN110488393A (zh) * 2019-09-24 2019-11-22 中国科学院大气物理研究所 X波段双偏振天气雷达定量测量降水方法及系统
CN110596787A (zh) * 2019-08-29 2019-12-20 成都锦江电子系统工程有限公司 一种基于x波段全固态双偏振雨量雷达的降水估测方法
CN110691988A (zh) * 2017-01-09 2020-01-14 诺维美特公司 用于估计降水率、特别是降雨率或降雪率的方法
KR102101436B1 (ko) * 2019-09-25 2020-04-16 대한민국 이중편파 기상레이더 기반의 우박 및 강우 영역별 강우 추정 장치 및 그 방법
CN111738339A (zh) * 2020-06-23 2020-10-02 中国民航大学 一种针对双偏振气象雷达低分辨数据的降水粒子分类方法
CN112526641A (zh) * 2020-12-10 2021-03-19 重庆市气象台 一种实时标识雨量观测值质量的方法、系统及设备
JP2021092407A (ja) * 2019-12-09 2021-06-17 コリア インスティテュート オフ コンストラクション テクノロジー 超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法，及び推定装置

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002350560A (ja) * 2001-05-30 2002-12-04 Foundation Of River & Basin Integrated Communications Japan レーダ雨量補正・配信システム
FR2827391A1 (fr) * 2001-07-11 2003-01-17 Centre Nat Rech Scient Technique pour l'estimation de pluie avec radar lorsque l'attenuation est negligeable
KR101431707B1 (ko) * 2013-11-26 2014-09-22 한국건설기술연구원 엑스밴드 이중편파 레이더 관측자료를 이용한 통합형 강우량 산정 방법
CN105548984A (zh) * 2015-12-10 2016-05-04 中国人民解放军理工大学 一种基于模糊逻辑的双偏振多普勒天气雷达地物杂波抑制方法
CN110691988A (zh) * 2017-01-09 2020-01-14 诺维美特公司 用于估计降水率、特别是降雨率或降雪率的方法
CN106950614A (zh) * 2017-02-28 2017-07-14 中船重工鹏力（南京）大气海洋信息系统有限公司 一种区域自动气象站小时雨量数据质量控制方法
CN107015229A (zh) * 2017-05-22 2017-08-04 新疆维吾尔自治区人工影响天气办公室 基于双偏振天气雷达的人工防雹作业指挥方法
CN108680920A (zh) * 2018-04-28 2018-10-19 深圳市雅码科技有限公司 一种基于双偏振雷达的灾害天气识别预警系统及方法
CN110018448A (zh) * 2019-03-30 2019-07-16 中国人民解放军陆军工程大学 基于任意取向微波链路的双偏振雷达衰减订正方法
CN110082842A (zh) * 2019-05-24 2019-08-02 北京敏视达雷达有限公司 一种降水估测方法及装置
CN110596787A (zh) * 2019-08-29 2019-12-20 成都锦江电子系统工程有限公司 一种基于x波段全固态双偏振雨量雷达的降水估测方法
CN110488393A (zh) * 2019-09-24 2019-11-22 中国科学院大气物理研究所 X波段双偏振天气雷达定量测量降水方法及系统
KR102101436B1 (ko) * 2019-09-25 2020-04-16 대한민국 이중편파 기상레이더 기반의 우박 및 강우 영역별 강우 추정 장치 및 그 방법
JP2021092407A (ja) * 2019-12-09 2021-06-17 コリア インスティテュート オフ コンストラクション テクノロジー 超短距離二重偏波レーダの多重高度観測資料を用いた降雨強度推定方法，及び推定装置
CN111738339A (zh) * 2020-06-23 2020-10-02 中国民航大学 一种针对双偏振气象雷达低分辨数据的降水粒子分类方法
CN112526641A (zh) * 2020-12-10 2021-03-19 重庆市气象台 一种实时标识雨量观测值质量的方法、系统及设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
井高飞、罗丽等: "双偏振雷达在地面自动雨量观测资料质量控制中的应用研究", <<气象>> *
宋文婷等: "基于S 波段双偏振雷达资料的降水粒子类型识别算法及应用", <<大学气象学报>> *
陈超等: "ＣＩＮ ＲＡ Ｄ－Ｓ Ａ 双偏振雷达资料在降水", <<气象>> *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115656947A (zh) * 2022-12-15 2023-01-31 南京信大气象科学技术研究院有限公司 一种双偏振天气雷达zdr偏差实时校正方法和装置
CN116381692A (zh) * 2023-04-18 2023-07-04 南京市气象台 一种基于x波段双偏振雷达的降水相态识别qpe算法
CN116381692B (zh) * 2023-04-18 2024-02-23 南京市气象台 一种基于x波段双偏振雷达的降水相态识别qpe算法
CN116303368B (zh) * 2023-04-24 2023-07-21 中国人民解放军国防科技大学 一种双偏振雷达体扫数据插值方法、装置、设备及介质

Also Published As

Publication number Publication date
CN113740934B (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
CN113740934B (zh) 一种基于s波段双偏振天气雷达降水估测方法
Goudenhoofdt et al. Evaluation of radar-gauge merging methods for quantitative precipitation estimates
CN111308581B (zh) 一种基于时空局部模型的雷达-雨量计联合降水估测方法
Berndt et al. Geostatistical merging of rain gauge and radar data for high temporal resolutions and various station density scenarios
CN110346844A (zh) Nriet基于云分类和机器学习的定量降水估测方法
Chumchean et al. Correcting of real-time radar rainfall bias using a Kalman filtering approach
Kirstetter et al. Toward a framework for systematic error modeling of spaceborne precipitation radar with NOAA/NSSL ground radar–based National Mosaic QPE
Alfieri et al. Time-dependent ZR relationships for estimating rainfall fields from radar measurements
CN110082842B (zh) 一种降水估测方法及装置
Chumchean et al. Radar rainfall error variance and its impact on radar rainfall calibration
Goudenhoofdt et al. Generation and verification of rainfall estimates from 10-yr volumetric weather radar measurements
CN109001701B (zh) 基于双偏振参数特征库匹配的雷达动态定量估测降水方法
CN110263838B (zh) 一种多传感器降水估计融合方法
CN113341389B (zh) 一种基于最大似然估计的微波散射计风场反演方法及系统
Kwon et al. Inter-comparison of radar rainfall rate using Constant Altitude Plan Position Indicator and hybrid surface rainfall maps
Pulkkinen et al. Probabilistic radar-gauge merging by multivariate spatiotemporal techniques
Wardah et al. Radar rainfall estimates comparison with kriging interpolation of gauged rain
CN116047524A (zh) 一种双偏振天气雷达数据质量实时评估方法及系统
Woodley Optimizing the measurement of convective rainfall in Florida
Kimpara et al. Evaluation of the radar-based quantitative precipitation estimation composite in Viet Nam
Yousefi Estimation of bias-corrected high-resolution radar precipitation maps using the radar and rain gauge network over Turkey
CN117805826B (zh) 基于mim网络和雷达拼图的分钟降水估测方法及系统
CN118011358B (zh) 一种双波段雷达雨滴谱反演方法及装置
CN114594444B (zh) 一种双偏振雷达降水率估测方法及系统
Smith Improving the Accuracy of Bias Corrected, Gridded Precipitation for Drought Monitoring

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Chen Hao

Inventor after: Wang Zhangwei

Inventor after: Lu Xinglai

Inventor after: Wang Han

Inventor before: Wang Zhangwei

Inventor before: Chen Hao

Inventor before: Wang Han

Inventor before: Lu Xinglai

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant