[go: up one dir, main page]

CN110596787A - 一种基于x波段全固态双偏振雨量雷达的降水估测方法 - Google Patents

一种基于x波段全固态双偏振雨量雷达的降水估测方法 Download PDF

Info

Publication number
CN110596787A
CN110596787A CN201910808851.1A CN201910808851A CN110596787A CN 110596787 A CN110596787 A CN 110596787A CN 201910808851 A CN201910808851 A CN 201910808851A CN 110596787 A CN110596787 A CN 110596787A
Authority
CN
China
Prior art keywords
rainfall
radial
precipitation
radar
deblurring
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.)
Pending
Application number
CN201910808851.1A
Other languages
English (en)
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.)
Chengdu Jinjiang Electronic System Engineering Co Ltd
Original Assignee
Chengdu Jinjiang Electronic System Engineering Co Ltd
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 Chengdu Jinjiang Electronic System Engineering Co Ltd filed Critical Chengdu Jinjiang Electronic System Engineering Co Ltd
Priority to CN201910808851.1A priority Critical patent/CN110596787A/zh
Publication of CN110596787A publication Critical patent/CN110596787A/zh
Pending legal-status Critical Current

Links

Classifications

    • 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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • 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)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于X波段全固态双偏振雨量雷达的降水估测方法，主要包括以下步骤：1）在5分钟内获取X波段全固态双偏振雨量雷达对于60km范围的降水探测体积扫描基数据；2）提取雷达数据中的反射率因子dBz、径向速度V、差分传播相移φDP、差分传播相移率KDP、差分反射率因子ZDR、零延迟相关系数ρHV，进行回波强度衰减订正、径向速度退模糊、径向干扰回波处理；3）通过相态识别算法区分小雨、中雨、大雨、暴雨、固态降水、地物杂波、其它杂波等；4）根据不同的降水过程，通过自适应的降水估测算法，输出半径60km的面雨量格点数据。本方案解决了雷达波束部分遮挡（遮挡率<20%）和雨衰造成的基数据偏差问题，并提升降水估测精度。

Description

一种基于X波段全固态双偏振雨量雷达的降水估测方法
技术领域
本发明涉及雷达领域，具体涉及一种基于X波段全固态双偏振雨量雷达的降水估测方法。
背景技术
X波段全固态双偏振雨量雷达作为气象探测新型主动遥感设备，具有实时估测其探测范围内各点的降水回波强度以及一定区域内降水范围分布的能力，因此可以及时获得空间连续性较好、范围较大的定量降水资料。同时，因具有偏振体制，探测云、雨目标的差分反射率、差分传播相移率、相关系数等单偏振雷达不具有的参数，能较准确识别水凝结物。X波段雨量雷达在探测局地天气过程时，降水信息的空间分辨率和连续性均好于S和C波段天气雷达。因此，对于监测区域内的中小河流洪水、城市内涝、地质灾害等预警预报均有重大意义。
X波段雨量雷达在测雨过程时，影响降水估测精度的因素主要包括：1)降雨对电磁波存在着一定的衰减效应，雷达波长越小，其衰减越严重。X波段雨量雷达衰减最大，C波段次之，S波段雷达最小。因此，X波段雨量雷达需要进行衰减订正，对提高区域降水量估测精度具有较大的意义。2)地物等杂波的干扰；3)降水估测算法对不同类型降水的适应性。
现有技术的缺点在于：
(1)雷达波束部分遮挡(遮挡率<20％)和雨衰造成的基数据偏差；
(2)难以明确区分小雨、中雨、大雨、暴雨、固态降水(含冰雹和降雪)等，使得将水估算不准确，降水估测精度较低。
发明内容
本发明的目的在于克服现有技术的不足，提供一种基于X波段全固态双偏振雨量雷达的降水估测方法，解决现有技术中存在的雷达波束部分遮挡(遮挡率<20％)和雨衰造成的基数据偏差问题，并提升降水估测精度。
本发明的目的是通过以下技术方案来实现的：
一种基于X波段全固态双偏振雨量雷达的降水估测方法，包括以下步骤：
S1：在时间T内获取X波段全固态双偏振雨量雷达对于半径R范围的降水探测体积扫描VOL基数据；
S2：提取雷达数据中的反射率因子dBz、径向速度V、差分传播相移φDP、差分传播相移率KDP、差分反射率因子ZDR、零延迟相关系数ρHV，进行回波强度衰减订正、径向速度退模糊、径向干扰回波处理；
S3：通过相态识别算法区分小雨、中雨、大雨、暴雨、固态降水、地物杂波；
S4：根据不同的降水过程，通过自适应的降水估测算法，输出半径R的面雨量格点数据。
作为本发明的进一步改进，所述步骤S1中在时间T内完成9个仰角的扫描。
作为本发明的进一步改进，所述时间T取值为3-8min，所述半径R取值为40km-80km。
作为本发明的进一步改进，所述步骤S2中回波强度衰减订正采用用Kdp法，基于雷达观测资料中的差传播相移KDP与衰减率AH有着良好的线性关系得到:
AH＝a1KDP (1)
订正后的ZH可以表示为
其中：BW为库长，单位:km，nr为雷达探测距离库数；
公式式(2)中对KDP进行处理:将KDP＜q1及KDP＞q2的值认为是数据异常，取零，其中q1＝0.1deg/km，q2＝3.0deg/km。
作为本发明的进一步改进，所述径向速度退模糊根据模糊的速度值和最大不模糊速度Vmax，给出一系列可能的速度值，其公式如下:
V＝Vfirst+2·n·Vmax (3)
其中称Vfirst为速度的初猜值，Vmax为最大不模糊速度，n为整数，V为可能的速度值；
通过公式(3)找到最弱的风场区来决定一系列初始参考径向和距离库，然后从这些径向和距离库出发，检查相邻库的连续性并对模糊速度进行订正。
作为本发明的进一步改进，所述模糊速度进行订正的步骤包括：
S01：寻找初始参考径向；
S02：利用初始参考径向对相邻的两根径向进行径向速度退模糊处理，得到三根好的初始径向；
S03：用得到的三根好径向同时作为参考径向来进行第一轮退模糊，算法从紧挨着三根初始径向的那根径向出发，从顺时针和逆时针两个方向同时逐根径向进行退模糊；不同于其它算法，该算法用得到的三根好径向同时作为参考径向来进行第一轮退模糊。算法从紧挨着三根初始径向的那根径向出发，从顺时针和逆时针两个方向同时逐根径向进行退模糊。这也是该算法与其它算法的不同之一，这样能更好地处理切变区域周围的径向速度模糊，同时也限制了180°范围内任何可能的错误扩散。
S04：第二轮径向逐根退模糊以及误差检查，第一轮径向退速度模糊之后，对在第一轮中没有被处理的所有点进行第二轮退模糊，第二轮退模糊除了在寻找参考速度时放宽了搜寻半径以外，其它步骤及方法与第一轮退模糊相同。
作为本发明的进一步改进，所述放宽了搜寻半径是指扩展到七根径向。
作为本发明的进一步改进，所述径向干扰回波处理包括：
S11：根据某个方向上的切向回波段长度值的分布特征，可准确定位径向干扰回波；
S12：确定其具体的方位后将该方位上的虚假数据抠除和赋零值，然后利用插值法，结合相邻径向上有效的雨强数据对赋零值区域进行插值填补。
作为本发明的进一步改进，所述插值法的公式为：
其中，M为径向干扰的方位宽度其分辨率为1°；Rjn为识别出的径向干扰第n个径向，n取值为(0～M)，j为第j个距离库格点的降水信息；K为插值有效径向的宽度；N为径向干扰降水带中每个径向所在的具体方位。
作为本发明的进一步改进，所述步骤S4包括如下子步骤：
S41：当降水强度为小雨时，利用Z-R关系法估测降水可以达到较好的精度，
其中，zh为水平偏振反射率，单位为mm6/m3；
系数a、b取值为a＝300,b＝1.4，或a＝200,b＝1.6，分别对应对流云和层状云降水；
这里所指的小雨，是指本领域内通用划分的降雨等级，即降雨量为10mm/h以下。
S42：当降水强度＜20mm/h时，采用R(ZH)进行降水估测优于利用R(ZH，ZDR)、R(KDP)方法；当降水强度＞20mm/h时，采用R(KDP)方法进行降水估测；
本发明的有益效果是：
1)采用双偏振参量Kdp完成对dBz和Zdr的校准，解决雷达波束部分遮挡(遮挡率<20％)和雨衰造成的基数据偏差；
2)降水估测之前，采用模糊逻辑算法对雷达探测范围各区域进行相态识别，区分小雨、中雨、大雨、暴雨、固态降水(含冰雹和降雪)等，为下一步的降水估测算法提供输入参考条件；
3)针对不同类型、不同强度的降水采用自适应的降水估测算法，进一步提升降水估测精度。
附图说明
图1为本发明流程图；
图2为体积扫描模式VCP21示意图；
图3为原始的dBz分布和经过进行干扰处理后的dBz分布对比示意图；
图4为本发明对比方式示意图；
图5为实施例1雷达累积雨量与雨量站雨量对比图；
图6为实施例2雷达累积雨量与雨量站雨量对比图；
图7为实施例3雷达累积雨量与雨量站雨量对比图。
具体实施方式
下面结合具体实施例进一步详细描述本发明的技术方案，但本发明的保护范围不局限于以下所述。
如图1所示，一种基于X波段全固态双偏振雨量雷达的降水估测方法，包括以下步骤：
S1：在时间T内获取X波段全固态双偏振雨量雷达对于半径R范围的降水探测体积扫描VOL基数据；
S2：提取雷达数据中的反射率因子dBz、径向速度V、差分传播相移φDP、差分传播相移率KDP、差分反射率因子ZDR、零延迟相关系数ρHV，进行回波强度衰减订正、径向速度退模糊、径向干扰回波处理；
S3：通过相态识别算法区分小雨、中雨、大雨、暴雨、固态降水、地物杂波以及其他杂波；其中固态降水包括雹、大雹、雨夹雹、霰、干雪、湿雪、冰晶等，利用模糊逻辑方法中的Beta隶属函数对降水粒子相态进行识别，是目前常用的粒子相态识别方法(Fuzzy logicHydrometeor Classification，简称FHC)，能对小雨滴、大雨滴、小雹、大雹、雨夹雹、霰、干雪、湿雪、冰晶等进行详细的识别和分类，本方案不再累述。
S4：根据不同的降水过程，通过自适应的降水估测算法，输出半径R的面雨量格点数据。
作为一种优选实施例，步骤S1中在时间T内完成9个仰角的扫描，扫描方式参考图2所示的中国气象局业务天气雷达VCP21扫描模式，时间T取值为3-8min，所述半径R取值为40km-80km，在本实施例中T取值为5min，半径R取值为60km。
作为一种优选实施例，所述步骤S2中回波强度衰减订正采用用Kdp法，基于雷达观测资料中的差传播相移KDP与衰减率AH有着良好的线性关系得到:
AH＝a1KDP (1)
订正后的ZH可以表示为
其中：BW为库长，单位:km，nr为雷达探测距离库数；
公式式(2)中对KDP进行处理:将KDP＜q1及KDP＞q2的值认为是数据异常，取零，其中q1＝0.1deg/km，q2＝3.0deg/km。
作为一种优选实施例，所述径向速度退模糊根据模糊的速度值和最大不模糊速度Vmax，给出一系列可能的速度值，其公式如下:
V＝Vfirst+2·n·Vmax (3)
其中称Vfirst为速度的初猜值，Vmax为最大不模糊速度，n为整数，V为可能的速度值；
通过公式(3)找到最弱的风场区来决定一系列初始参考径向和距离库，然后从这些径向和距离库出发，检查相邻库的连续性并对模糊速度进行订正。
首先通过找到最弱的风场区来决定一系列初始参考径向和距离库，然后从这些径向和距离库出发，检查相邻库的连续性并对模糊速度进行订正。
该算法主要基于风场连续性原则，算法简单独立能自动化运行，计算性能较高。
作为本发明的进一步改进，所述模糊速度进行订正的步骤包括：
S01：寻找初始参考径向；
S02：利用初始参考径向对相邻的两根径向进行径向速度退模糊处理，得到三根好的初始径向；
S03：用得到的三根好径向同时作为参考径向来进行第一轮退模糊，算法从紧挨着三根初始径向的那根径向出发，从顺时针和逆时针两个方向同时逐根径向进行退模糊；
S04：第二轮径向逐根退模糊以及误差检查，第一轮径向退速度模糊之后，对在第一轮中没有被处理的所有点进行第二轮退模糊，第二轮退模糊除了在寻找参考速度时放宽了搜寻半径以外(扩展到七根径向)，其它步骤及方法与第一轮退模糊相同。
雨量雷达系统生成区域降水量产品中，某些时刻存在条幅状虚假降水带，经过分析为径向干扰回波所生成。这种虚假回波带污染了区域内的降水量分布，给有径向干扰的区域造成了降雨的虚警。
本方案提供的径向干扰回波识别方法是在SCIT(Storm Cell Identificationand Trackin风暴单体识别与跟踪)方法划分风暴段和风暴块的基础上，通过改变相关条件和参数实现。在识别径向干扰回波时，改变了划分回波段方向，使用沿雷达切向方向划分切向回波段。切向回波段包含了回波在某个距离上的起、止方位角，并能反映回波在该距离上的切向宽度。针对干扰回波的分析，可发现其在某个方向上距离库宽度较小且基本不变。因此对应方向上所有的切向回波段长度值基本相同且都较小。因此，根据某个方向上的切向回波段长度值的分布特征，可准确定位径向干扰回波。
作为一种优选实施例，所述径向干扰回波处理包括：
S11：根据某个方向上的切向回波段长度值的分布特征，可准确定位径向干扰回波；
S12：确定其具体的方位后将该方位上的虚假数据抠除和赋零值，然后利用插值法，结合相邻径向上有效的雨强数据对赋零值区域进行插值填补。
所述插值法的公式为：
其中，M为径向干扰的方位宽度其分辨率为1°；Rjn为识别出的径向干扰第n个径向，n取值为(0～M)，j为第j个距离库格点的降水信息；K为插值有效径向的宽度；N为径向干扰降水带中每个径向所在的具体方位，如图3所示即给出了利用本方法处理后原始的dBz分布和经过进行干扰处理后的dBz分布对比示意图。
作为一种优选实施例，所述步骤S4包括如下子步骤：
S41：当降水强度为小雨时，利用Z-R关系法估测降水可以达到较好的精度，
其中，zh为水平偏振反射率，单位为mm6/m3；
系数a、b取值为a＝300,b＝1.4，或a＝200,b＝1.6，分别对应对流云和层状云降水；
S42：当降水强度＜20mm/h时，采用R(ZH)进行降水估测优于利用R(ZH，ZDR)、R(KDP)方法；当降水强度＞20mm/h时，采用R(KDP)方法进行降水估测；
在具体实施过程中，将雨量雷达测量的累积雨量值与对应时间段对应区域内雨量站测量值进行对比。在对比过程中，考虑到风速等因素的影响，不直接采用单点雨量对比方法，采用雨量站与雨量雷达对应雨量站点周围9个格点的均值进行比较的方法，以确保对比方法的准确性。如图4所示,，以九宫格中心点代表雨量雷达探测范围内某雨量站位置，周围八个点代表雨量站周围的雨量雷达测量格点，圆圈代表雨量雷达探测范围。
选取比率(RATIO)和平均相对误差(ARE)2个参数对雷达降水估测度进行评估。
其中:Ri为雷达雨量，Gi为雨量站雨量，n为雷达探测区域内雨量站数量。
具体应用实例1
降水过程时间为：2015年06月02日20:00–2015年06月04日00:00，地点：某地雨量雷达。该过程为混合型降水过程，持续时间长，共52个小时，降水过程一直伴随着层状云降水和对流云降水，大部分雨量站点24h雨量超过100mm，属于大暴雨降水过程。该过程雨量雷达累积雨量相与雨量站雨量整体接近，两者相关性较好，相对误差较小，其评估结果见表1，评估结果如表1所示，雷达累积雨量与雨量站雨量对比见图5所示。
表1 南昌评估结果
具体应用实例2
降水过程时间为：2016年05月25日09:00-2016年05月27日00:00。地点：某地雨量雷达。该过程为稳定性层状云降水过程，持续时间较长，降水空间结构均匀，大部分雨量站点24h雨量均在50mm以下。该过程雨量雷达累积雨量与雨量站雨量整体接近，两者相关性较好，相对误差较小，评估结果如表2所示，雷达累积雨量与雨量站雨量对比见图6所示。
表2 评估结果
具体应用实例3
降水过程时间为：2015年07月9日9：00–2015年07月10日00：00，地点：某地雨量雷达。该过程为稳定性层状云降水过程，降水空间结构均匀，大部分雨量站点24h雨量均在50mm以下。该过程雨量雷达累积雨量与雨量站雨量整体接近，两者相关性较好，相对误差较小，评估结果如表3所示，雷达累积雨量与雨量站雨量对比见图7所示。
表3 评估结果
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。

Claims (10)

1.一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，包括以下步骤：
S1：在时间T内获取X波段全固态双偏振雨量雷达对于半径R范围的降水探测体积扫描VOL基数据；
S2：提取雷达数据中的反射率因子dBz、径向速度V、差分传播相移φDP、差分传播相移率KDP、差分反射率因子ZDR、零延迟相关系数ρHV，进行回波强度衰减订正、径向速度退模糊、径向干扰回波处理；
S3：通过相态识别算法区分小雨、中雨、大雨、暴雨、固态降水、地物杂波；
S4：根据不同的降水过程，通过自适应的降水估测算法，输出半径R的面雨量格点数据。
2.根据权利要求1所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述步骤S1中在时间T内完成9个仰角的扫描。
3.根据权利要求2所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述时间T取值为3-8min，所述半径R取值为40km-80km。
4.根据权利要求1所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述步骤S2中回波强度衰减订正采用用Kdp法，基于雷达观测资料中的差传播相移KDP与衰减率AH有着良好的线性关系得到:
AH＝a1KDP (1)
订正后的ZH可以表示为
其中：BW为库长，单位:km，nr为雷达探测距离库数；
公式式(2)中对KDP进行处理:将KDP＜q1及KDP＞q2的值认为是数据异常，取零，其中q1＝0.1deg/km，q2＝3.0deg/km。
5.根据权利要求4所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述径向速度退模糊根据模糊的速度值和最大不模糊速度Vmax，给出一系列可能的速度值，其公式如下:
V＝Vfirst+2·n·Vmax (3)
其中称Vfirst为速度的初猜值，Vmax为最大不模糊速度，n为整数，V为可能的速度值；
通过公式(3)找到最弱的风场区来决定一系列初始参考径向和距离库，然后从这些径向和距离库出发，检查相邻库的连续性并对模糊速度进行订正。
6.根据权利要求5所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述模糊速度进行订正的步骤包括：
S01：寻找初始参考径向；
S02：利用初始参考径向对相邻的两根径向进行径向速度退模糊处理，得到三根好的初始径向；
S03：用得到的三根好径向同时作为参考径向来进行第一轮退模糊，算法从紧挨着三根初始径向的那根径向出发，从顺时针和逆时针两个方向同时逐根径向进行退模糊；
S04：第二轮径向逐根退模糊以及误差检查，第一轮径向退速度模糊之后，对在第一轮中没有被处理的所有点进行第二轮退模糊，第二轮退模糊除了在寻找参考速度时放宽了搜寻半径以外，其它步骤及方法与第一轮退模糊相同。
7.根据权利要求6所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述放宽了搜寻半径是指扩展到七根径向。
8.根据权利要求7所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述径向干扰回波处理包括：
S11：根据某个方向上的切向回波段长度值的分布特征，可准确定位径向干扰回波；
S12：确定其具体的方位后将该方位上的虚假数据抠除和赋零值，然后利用插值法，结合相邻径向上有效的雨强数据对赋零值区域进行插值填补。
9.根据权利要求8所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述插值法的公式为：
其中，M为径向干扰的方位宽度其分辨率为1°；Rjn为识别出的径向干扰第n个径向，n取值为(0～M)，j为第j个距离库格点的降水信息；K为插值有效径向的宽度；N为径向干扰降水带中每个径向所在的具体方位。
10.根据权利要求1-9中任一项所述的一种基于X波段全固态双偏振雨量雷达的降水估测方法，其特征在于，所述步骤S4包括如下子步骤：
S41：当降水强度为小雨时，利用Z-R关系法估测降水可以达到较好的精度，
其中，zh为水平偏振反射率，单位为mm6/m3；
系数a、b取值为a＝300,b＝1.4，或a＝200,b＝1.6，分别对应对流云和层状云降水；
S42：当降水强度＜20mm/h时，采用R(ZH)进行降水估测优于利用R(ZH，ZDR)、R(KDP)方法；当降水强度＞20mm/h时，采用R(KDP)方法进行降水估测；
CN201910808851.1A 2019-08-29 2019-08-29 一种基于x波段全固态双偏振雨量雷达的降水估测方法 Pending CN110596787A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910808851.1A CN110596787A (zh) 2019-08-29 2019-08-29 一种基于x波段全固态双偏振雨量雷达的降水估测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910808851.1A CN110596787A (zh) 2019-08-29 2019-08-29 一种基于x波段全固态双偏振雨量雷达的降水估测方法

Publications (1)

Publication Number Publication Date
CN110596787A true CN110596787A (zh) 2019-12-20

Family

ID=68856267

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910808851.1A Pending CN110596787A (zh) 2019-08-29 2019-08-29 一种基于x波段全固态双偏振雨量雷达的降水估测方法

Country Status (1)

Country Link
CN (1) CN110596787A (zh)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257884A (zh) * 2020-02-18 2020-06-09 河海大学 链路长短可调的微波雨衰规律人工降雨实验方法
CN111323782A (zh) * 2020-04-16 2020-06-23 内蒙古工业大学 一种基于模糊逻辑的雷达同步监测方法
CN111856511A (zh) * 2020-07-07 2020-10-30 北方激光研究院有限公司 改变扫描夹角的dbs风场扫描方法
CN113238230A (zh) * 2021-04-12 2021-08-10 国网河南省电力公司电力科学研究院 一种针对电网生产的夏季强对流引起的大风预警方法
CN113341420A (zh) * 2021-06-07 2021-09-03 中国气象局上海台风研究所(上海市气象科学研究所) 适用不同气候态降雨的双偏振气象雷达雨滴谱反演方法
CN113391314A (zh) * 2021-06-09 2021-09-14 航天新气象科技有限公司 基于多频段气象雷达的天气探测方法、系统及存储介质
CN113466856A (zh) * 2021-08-04 2021-10-01 广州市气象台 基于x波段双偏振相控阵雷达的林火早期识别与预警方法
CN113740934A (zh) * 2021-08-18 2021-12-03 浙江省大气探测技术保障中心 一种基于s波段双偏振天气雷达降水估测方法
CN113933809A (zh) * 2021-09-30 2022-01-14 中山大学 一种基于Kmeans聚类的降水粒子识别方法及装置
EP3964854A1 (en) * 2020-09-03 2022-03-09 Vaisala Inc. Rfi mitigation in weather radar data
CN114236550A (zh) * 2021-12-21 2022-03-25 雷象科技（北京）有限公司 单波段雷达观测平台及其正交综合退径向速度模糊方法
CN114280572A (zh) * 2022-03-04 2022-04-05 浙江省气象台 去除信号干扰杂波的单雷达回波质控方法、系统及终端
CN114488160A (zh) * 2022-04-02 2022-05-13 南京师范大学 顾及三维风场影响的雷达降雨估算误差校正方法
CN114675278A (zh) * 2022-03-31 2022-06-28 山东省气象科学研究所(山东省海洋气象科学研究所、山东省气象局培训中心) 利用s波段双偏振雷达自动识别融化层的方法
CN115061105A (zh) * 2022-06-30 2022-09-16 广东纳睿雷达科技股份有限公司 一种双偏振雷达快速标定方法、装置及存储介质
CN116381692A (zh) * 2023-04-18 2023-07-04 南京市气象台 一种基于x波段双偏振雷达的降水相态识别qpe算法
CN118519118A (zh) * 2024-07-19 2024-08-20 北京城市气象研究院 基于瓦片分区z和kdp联合定量降水估测方法
CN118534466A (zh) * 2024-05-27 2024-08-23 北京理工大学珠海学院 一种基于人工智能的测雨雷达系统及方法
CN119001656A (zh) * 2024-10-23 2024-11-22 成都远望探测技术有限公司 一种多波段天气雷达多类型自适应径向干扰滤除方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316930A (zh) * 2014-10-29 2015-01-28 安徽四创电子股份有限公司 一种基于双线偏振多普勒天气雷达探测的雨强估算方法
CN107064885A (zh) * 2017-02-21 2017-08-18 水利部南京水利水文自动化研究所 X波段雨量雷达径向干扰识别与消除方法
CN107526067A (zh) * 2017-03-07 2017-12-29 中国气象局武汉暴雨研究所 全自动多普勒天气雷达径向速度退模糊算法
US20180017676A1 (en) * 2015-02-02 2018-01-18 Furuno Electric Co., Ltd. Detection apparatus and detecting method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104316930A (zh) * 2014-10-29 2015-01-28 安徽四创电子股份有限公司 一种基于双线偏振多普勒天气雷达探测的雨强估算方法
US20180017676A1 (en) * 2015-02-02 2018-01-18 Furuno Electric Co., Ltd. Detection apparatus and detecting method
CN107064885A (zh) * 2017-02-21 2017-08-18 水利部南京水利水文自动化研究所 X波段雨量雷达径向干扰识别与消除方法
CN107526067A (zh) * 2017-03-07 2017-12-29 中国气象局武汉暴雨研究所 全自动多普勒天气雷达径向速度退模糊算法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
冯亮等: "X 波段双偏振雷达水凝物粒子相态识别应用研究", 《气候与环境研究》 *
吴欢: "X波段双线偏振雷达资料的衰减订正和地物消除以及降水物类型判别研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 *
梅垚等: "青藏高原对流云的偏振雷达观测研究", 《气象学报》 *

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257884B (zh) * 2020-02-18 2020-10-16 河海大学 链路长短可调的微波雨衰规律人工降雨实验方法
CN111257884A (zh) * 2020-02-18 2020-06-09 河海大学 链路长短可调的微波雨衰规律人工降雨实验方法
CN111323782A (zh) * 2020-04-16 2020-06-23 内蒙古工业大学 一种基于模糊逻辑的雷达同步监测方法
CN111323782B (zh) * 2020-04-16 2022-05-27 内蒙古工业大学 一种基于模糊逻辑的雷达同步监测方法
CN111856511A (zh) * 2020-07-07 2020-10-30 北方激光研究院有限公司 改变扫描夹角的dbs风场扫描方法
US11474199B2 (en) 2020-09-03 2022-10-18 Vaisala, Inc. RFI mitigation in weather radar data
EP3964854A1 (en) * 2020-09-03 2022-03-09 Vaisala Inc. Rfi mitigation in weather radar data
CN113238230A (zh) * 2021-04-12 2021-08-10 国网河南省电力公司电力科学研究院 一种针对电网生产的夏季强对流引起的大风预警方法
CN113238230B (zh) * 2021-04-12 2023-07-14 国网河南省电力公司电力科学研究院 一种针对电网生产的夏季强对流引起的大风预警方法
CN113341420B (zh) * 2021-06-07 2023-07-18 中国气象局上海台风研究所(上海市气象科学研究所) 适用不同气候态降雨的双偏振气象雷达雨滴谱反演方法
CN113341420A (zh) * 2021-06-07 2021-09-03 中国气象局上海台风研究所(上海市气象科学研究所) 适用不同气候态降雨的双偏振气象雷达雨滴谱反演方法
CN113391314A (zh) * 2021-06-09 2021-09-14 航天新气象科技有限公司 基于多频段气象雷达的天气探测方法、系统及存储介质
CN113391314B (zh) * 2021-06-09 2024-08-13 航天新气象科技有限公司 基于多频段气象雷达的天气探测方法、系统及存储介质
CN113466856A (zh) * 2021-08-04 2021-10-01 广州市气象台 基于x波段双偏振相控阵雷达的林火早期识别与预警方法
CN113740934A (zh) * 2021-08-18 2021-12-03 浙江省大气探测技术保障中心 一种基于s波段双偏振天气雷达降水估测方法
CN113933809B (zh) * 2021-09-30 2022-07-29 中山大学 一种基于Kmeans聚类的降水粒子识别方法及装置
CN113933809A (zh) * 2021-09-30 2022-01-14 中山大学 一种基于Kmeans聚类的降水粒子识别方法及装置
CN114236550A (zh) * 2021-12-21 2022-03-25 雷象科技（北京）有限公司 单波段雷达观测平台及其正交综合退径向速度模糊方法
CN114236550B (zh) * 2021-12-21 2024-05-10 浙江宜通华盛科技有限公司 单波段雷达观测平台及其正交综合退径向速度模糊方法
CN114280572B (zh) * 2022-03-04 2022-08-09 浙江省气象台 去除信号干扰杂波的单雷达回波质控方法、系统及终端
CN114280572A (zh) * 2022-03-04 2022-04-05 浙江省气象台 去除信号干扰杂波的单雷达回波质控方法、系统及终端
CN114675278B (zh) * 2022-03-31 2025-03-18 山东省气象科学研究所(山东省海洋气象科学研究所、山东省气象局培训中心) 利用s波段双偏振雷达自动识别融化层的方法
CN114675278A (zh) * 2022-03-31 2022-06-28 山东省气象科学研究所(山东省海洋气象科学研究所、山东省气象局培训中心) 利用s波段双偏振雷达自动识别融化层的方法
CN114488160A (zh) * 2022-04-02 2022-05-13 南京师范大学 顾及三维风场影响的雷达降雨估算误差校正方法
CN115061105A (zh) * 2022-06-30 2022-09-16 广东纳睿雷达科技股份有限公司 一种双偏振雷达快速标定方法、装置及存储介质
CN115061105B (zh) * 2022-06-30 2023-07-14 广东纳睿雷达科技股份有限公司 一种双偏振雷达快速标定方法、装置及存储介质
CN116381692B (zh) * 2023-04-18 2024-02-23 南京市气象台 一种基于x波段双偏振雷达的降水相态识别qpe算法
CN116381692A (zh) * 2023-04-18 2023-07-04 南京市气象台 一种基于x波段双偏振雷达的降水相态识别qpe算法
CN118534466A (zh) * 2024-05-27 2024-08-23 北京理工大学珠海学院 一种基于人工智能的测雨雷达系统及方法
CN118534466B (zh) * 2024-05-27 2024-11-29 北京理工大学珠海学院 一种基于人工智能的测雨雷达系统及方法
CN118519118A (zh) * 2024-07-19 2024-08-20 北京城市气象研究院 基于瓦片分区z和kdp联合定量降水估测方法
CN119001656A (zh) * 2024-10-23 2024-11-22 成都远望探测技术有限公司 一种多波段天气雷达多类型自适应径向干扰滤除方法
CN119001656B (zh) * 2024-10-23 2024-12-20 成都远望探测技术有限公司 一种多波段天气雷达多类型自适应径向干扰滤除方法

Similar Documents

Publication Publication Date Title
CN110596787A (zh) 一种基于x波段全固态双偏振雨量雷达的降水估测方法
US9348015B2 (en) Integrated rainfall estimation method using X-band dual-polarimetric radar measurement data
Wolfensberger et al. Detection and characterization of the melting layer based on polarimetric radar scans
Lim et al. Precipitation classification and quantification using X-band dual-polarization weather radar: Application in the Hydrometeorology Testbed
Rico‐Ramirez et al. A high‐resolution radar experiment on the island of Jersey
Kalogiros et al. Evaluation of a new polarimetric algorithm for rain-path attenuation correction of X-band radar observations against disdrometer
Gorgucci et al. The influence of antenna radome on weather radar calibration and its real-time assessment
Nanding et al. Precipitation measurement with weather radars
Friedrich et al. Influence of ground clutter contamination on polarimetric radar parameters
Huuskonen et al. Radar performance monitoring using the angular width of the solar image
Picciotti et al. Coupling X-band dual-polarized mini-radars and hydro-meteorological forecast models: the HYDRORAD project
Ruzanski et al. Nowcasting rainfall fields derived from specific differential phase
Bharadwaj et al. Signal processing system for the CASA Integrated Project I radars
CN118778002A (zh) 一种双偏振天气雷达地物识别方法
Chen et al. Jensen–Shannon distance-based filter and unsupervised evaluation metrics for polarimetric weather radar processing
Capozzi et al. X-band weather radar monitoring of precipitation fields in Naples urban areas: data quality, comparison and analysis
Wang et al. Inversion and assessment of swell waveheights from HF radar spectra in the Iroise Sea
Nicol et al. Techniques for improving ground clutter identification
Pejcic et al. Polarimetric radar-based methods for evaluation of hydrometeor mixtures in numerical weather prediction models
Burgemeister et al. Reanalysis of multi-year high-resolution X-band weather radar observations in Hamburg
CN107748357B (zh) 一种fir滤波器的雷达晴天图设计方法
Kabeche et al. Quantitative precipitation estimation (QPE) in the French Alps with a dense network of polarimetric X-band radars
Iervolino et al. Roughness parameters estimation of sea surface from SAR images
CN116299475B (zh) 一种c波段双偏振天气雷达低仰角的观测方法
CN114442103B (zh) 双极化天气雷达差分传播相移率估计方法、系统、设备

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20191220

RJ01 Rejection of invention patent application after publication