CN115656947B - 一种双偏振天气雷达zdr偏差实时校正方法和装置 - Google Patents

一种双偏振天气雷达zdr偏差实时校正方法和装置 Download PDF

Info

Publication number
CN115656947B
CN115656947B CN202211609704.XA CN202211609704A CN115656947B CN 115656947 B CN115656947 B CN 115656947B CN 202211609704 A CN202211609704 A CN 202211609704A CN 115656947 B CN115656947 B CN 115656947B
Authority
CN
China
Prior art keywords
radar
dual
deviation
temperature
zdr
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.)
Active
Application number
CN202211609704.XA
Other languages
English (en)
Other versions
CN115656947A (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.)
Nanjing Xinda Meteorological Science And Technology Co ltd
Original Assignee
Nanjing Xinda Meteorological Science And Technology 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 Nanjing Xinda Meteorological Science And Technology Co ltd filed Critical Nanjing Xinda Meteorological Science And Technology Co ltd
Priority to CN202211609704.XA priority Critical patent/CN115656947B/zh
Publication of CN115656947A publication Critical patent/CN115656947A/zh
Application granted granted Critical
Publication of CN115656947B publication Critical patent/CN115656947B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种双偏振天气雷达ZDR偏差实时校正方法和装置。该方法包括从双偏振天气雷达的信号处理机实时采集雷达观测数据，雷达观测数据包括差分反射率因子；实时采集双偏振天气雷达的机内温度和机外温度；根据雷达观测数据、机内温度和机外温度计算偏差，根据所述偏差生成修正ZDR偏差指令，并将修正ZDR偏差指令发送至双偏振天气雷达，以使双偏振天气雷达根据所述雷达观测数据中的差分反射率因子和偏差计算校正后的差分反射率因子。本发明可实时对ZDR进行校正，在后续使用过程中，可通过对雷达观测数据进行分析来进行筛选，以不断补充样本来更新，保证校正的准确性和代表性；提高双偏振雷达的数据质量，具有广阔的应用前景。

Description

一种双偏振天气雷达ZDR偏差实时校正方法和装置
技术领域
本发明涉及天气雷达技术领域，具体涉及一种双偏振天气雷达ZDR偏差实时校正方法和装置。
背景技术
双偏振天气雷达交替或同时发射水平和垂直偏振波，获取云雨粒子的大小、形状、数量、相态等观测信息，其探测数据被广泛应用于云物理降水研究、灾害性强天气预警、数值天气预报等多个领域。差分反射率因子（ZDR）是双偏振天气雷达的一个重要观测量，它的准确性直接决定了双偏振天气雷达的应用效果。但影响ZDR数据精度的因素很多，其中主要的是水平与垂直的双通道（天线、波导、双接收机等）差异。由于受硬件加工精度、元器件参数漂移、热噪声起伏等影响，双偏振天气雷达的水平垂直通道参数不能完全一致，所以在观测过程中需要定期校正ZDR偏差，也称为ZDR标定。
发明专利CN201310202199.1描述了采用双通道全路径收发标校技术，减少双通道收发系统误差，即修正ZDR硬件偏差。这种技术的缺点是无法实时标定天线、馈线（波导），同时存在标校信号误差的问题。此外，基于自然源信号的标定方法被广泛应用于天气雷达，但都需要一定的自然条件。太阳法在雷达天线刚好指向太阳时根据观测值校正ZDR偏差；小雨法在小雨条件下校正ZDR偏差；垂直观测法在雷达正上方有降雨或降雪时校正ZDR偏差；雨滴谱法在雨滴谱仪正上方有降雨时才可实施。
由于双偏振天气雷达水平和垂直通道的热膨胀的差异，ZDR偏差是随环境温度不断变化的，最佳方案是实时校正或标定。上述方法均需要一定的前提条件，不能实时校正。
发明内容
本发明的目的是针对现有技术存在的不足，提供一种双偏振天气雷达ZDR偏差实时校正方法和装置。
为实现上述目的，在第一方面，本发明提供了一种双偏振天气雷达ZDR偏差实时校正方法，包括：
从双偏振天气雷达的信号处理机实时采集雷达观测数据，所述雷达观测数据包括 差分反射率因子
Figure DEST_PATH_IMAGE001
；
实时采集双偏振天气雷达的机内温度
Figure 619876DEST_PATH_IMAGE002
和机外温度
Figure DEST_PATH_IMAGE003
；
根据所述雷达观测数据、机内温度
Figure 689244DEST_PATH_IMAGE002
和机外温度
Figure 384536DEST_PATH_IMAGE003
计算偏差
Figure 764571DEST_PATH_IMAGE004
， 具体如下：
Figure DEST_PATH_IMAGE005
其中，
Figure 583054DEST_PATH_IMAGE006
分别为根据样本数据拟合出的系数；
根据所述偏差
Figure DEST_PATH_IMAGE007
生成修正ZDR偏差指令，并将所述修正ZDR偏差指令发送至 双偏振天气雷达，以使所述双偏振天气雷达根据所述雷达观测数据中的差分反射率因子
Figure 342981DEST_PATH_IMAGE008
和偏差
Figure 576385DEST_PATH_IMAGE004
计算校正后的差分反射率因子
Figure DEST_PATH_IMAGE009
，具体如下：
Figure 794614DEST_PATH_IMAGE010
。
进一步的，所述雷达观测数据还包括仰角、高度、反射率因子Z和相关系数CC；
所述方法还包括：
分别判断所述仰角是否大于仰角阈值
Figure DEST_PATH_IMAGE011
、所述反射率因子Z是否大于反射 率因子阈值
Figure 151096DEST_PATH_IMAGE012
、所述相关系数CC是否大于相关系数阈值
Figure 674350DEST_PATH_IMAGE013
、所述高度是否小 于高度阈值
Figure 711445DEST_PATH_IMAGE014
；
若所述仰角大于仰角阈值
Figure DEST_PATH_IMAGE015
，且所述反射率因子Z大于反射率因子阈值
Figure 790040DEST_PATH_IMAGE012
，且所述相关系数CC大于相关系数阈值
Figure 519967DEST_PATH_IMAGE016
，且所述高度小于高度阈值
Figure DEST_PATH_IMAGE017
，则将当前雷达观测数据中的差分反射率因子
Figure 23193DEST_PATH_IMAGE018
作为当前的偏差
Figure DEST_PATH_IMAGE019
， 并将其与当前采集的机内温度
Figure 847668DEST_PATH_IMAGE020
和机外温度
Figure 265923DEST_PATH_IMAGE003
作为最新样本进行存储，并通过所 述最新样本重新拟合更新
Figure DEST_PATH_IMAGE021
。
进一步的，所述仰角阈值
Figure 150440DEST_PATH_IMAGE022
为双偏振天气雷达的天线扫描范围的最大值。
进一步的，所述反射率因子阈值
Figure DEST_PATH_IMAGE023
根据双偏振天气雷达的接收机的灵敏度 进行调整，且其大小与接收机的灵敏度呈反比。
进一步的，所述高度阈值
Figure 898820DEST_PATH_IMAGE024
为本地当前的气温零度层的平均高度减0.5 km。
在第二方面，本发明提供了一种双偏振天气雷达ZDR偏差实时校正装置，包括：
网络通讯模块，用以从双偏振天气雷达的信号处理机实时采集雷达观测数据，所 述雷达观测数据包括差分反射率因子
Figure DEST_PATH_IMAGE025
；
机内温度传感器和机外温度传感器，分别用以采集双偏振天气雷达的机内温度
Figure 995827DEST_PATH_IMAGE026
和机外温度
Figure 528308DEST_PATH_IMAGE027
主控制板，用以根据所述雷达观测数据、机内温度
Figure 340319DEST_PATH_IMAGE028
和机外温度
Figure 591040DEST_PATH_IMAGE029
计算 偏差
Figure 39208DEST_PATH_IMAGE030
，具体如下：
Figure 935183DEST_PATH_IMAGE031
其中，
Figure 912235DEST_PATH_IMAGE032
分别为根据样本数据拟合出的系数；
所述主控制板还用以根据所述偏差
Figure 119095DEST_PATH_IMAGE033
生成修正ZDR偏差指令，并控制网络 通讯模块将所述修正ZDR偏差指令发送至双偏振天气雷达，以使所述双偏振天气雷达根据 所述雷达观测数据中的差分反射率因子
Figure 105374DEST_PATH_IMAGE034
和偏差
Figure 883886DEST_PATH_IMAGE033
计算校正后的差分反射率 因子
Figure 31840DEST_PATH_IMAGE035
，具体如下：
Figure 257154DEST_PATH_IMAGE010
。
进一步的，所述雷达观测数据还包括仰角、高度、反射率因子Z和相关系数CC；
所述主控制板还用以分别判断所述仰角是否大于仰角阈值
Figure 30812DEST_PATH_IMAGE036
、所述反射率 因子Z是否大于反射率因子阈值
Figure 409970DEST_PATH_IMAGE037
、所述相关系数CC是否大于相关系数阈值
Figure 259984DEST_PATH_IMAGE038
、所述高度是否小于高度阈值
Figure 972594DEST_PATH_IMAGE039
；
若所述仰角大于仰角阈值
Figure 769517DEST_PATH_IMAGE040
，且所述反射率因子Z大于反射率因子阈值
Figure 725884DEST_PATH_IMAGE037
，且所述相关系数CC大于相关系数阈值
Figure 481219DEST_PATH_IMAGE041
，且所述高度小于高度阈值
Figure 399234DEST_PATH_IMAGE042
，则将当前雷达观测数据中的差分反射率因子
Figure 525148DEST_PATH_IMAGE008
作为当前的偏差
Figure 126899DEST_PATH_IMAGE019
， 并将其与当前采集的机内温度
Figure 53136DEST_PATH_IMAGE043
和机外温度
Figure 740338DEST_PATH_IMAGE003
作为最新样本，并利用最新样本重 新拟合更新
Figure 416082DEST_PATH_IMAGE044
；
还包括存储器，用以对所述最新样本进行存储。
进一步的，所述仰角阈值
Figure 606761DEST_PATH_IMAGE040
为双偏振天气雷达的天线扫描范围的最大值。
进一步的，所述反射率因子阈值
Figure 969478DEST_PATH_IMAGE037
根据双偏振天气雷达的接收机的灵敏度 进行调整，且其大小与接收机的灵敏度呈反比。
进一步的，所述高度阈值
Figure 143976DEST_PATH_IMAGE045
为本地当前的气温零度层的平均高度减0.5 km。
有益效果：本发明通过采集机内温度和机外温度，并通过初始的样本系数来获得初始的系数取值，可实时对ZDR进行校正；在后续使用过程中，可通过对雷达观测数据进行分析来进行筛选，以不断补充样本来更新，保证校正的准确性和代表性；可以直接应用于双偏振雷达系统，解决ZDR偏差实时校正的问题，提高双偏振雷达的数据质量，具有广阔的应用前景。
附图说明
图1是一种双偏振天气雷达ZDR偏差实时校正方法的流程示意图；
图2是一种双偏振天气雷达ZDR偏差实时校正装置的结构示意图。
具体实施方式
下面结合附图和具体实施例，进一步阐明本发明，本实施例在以本发明技术方案为前提下进行实施，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围。
如图1所示，本发明实施例提供了一种双偏振天气雷达ZDR偏差实时校正方法，包括：
从双偏振天气雷达的信号处理机实时采集雷达观测数据，雷达观测数据包括差分 反射率因子
Figure 635130DEST_PATH_IMAGE018
。此外，雷达观测数据一般还包括方位角、仰角、时间、库长、库数、高度等 参数信息，以及反射率因子Z、差分反射率因子ZDR、相关系数CC等数据信息。
实时采集双偏振天气雷达的机内温度
Figure 680315DEST_PATH_IMAGE046
和机外温度
Figure 666463DEST_PATH_IMAGE003
。机内温度
Figure 328258DEST_PATH_IMAGE020
可以是雷达机内靠近雷达接收机的位置的温度，机外温度
Figure 611384DEST_PATH_IMAGE003
可以是雷达天线底座外侧 附近的温度。
根据雷达观测数据、机内温度
Figure 776655DEST_PATH_IMAGE047
和机外温度
Figure 215595DEST_PATH_IMAGE003
计算偏差
Figure 833527DEST_PATH_IMAGE048
，具体 如下：
Figure 111887DEST_PATH_IMAGE049
其中，
Figure 397244DEST_PATH_IMAGE050
分别为根据样本数据拟合出的系数。在投入使用前，一般需要 不少于100个初始样本数据来拟合得出初始的
Figure 7086DEST_PATH_IMAGE050
，在投入使用后即可根据
Figure 101862DEST_PATH_IMAGE021
以及机内温度
Figure 455352DEST_PATH_IMAGE051
和机外温度
Figure 595215DEST_PATH_IMAGE003
计算得出偏差
Figure 641537DEST_PATH_IMAGE052
。
根据偏差
Figure 251639DEST_PATH_IMAGE033
生成修正ZDR偏差指令，并将修正ZDR偏差指令发送至双偏振天 气雷达，以使双偏振天气雷达根据所述雷达观测数据中的差分反射率因子
Figure 408820DEST_PATH_IMAGE053
和偏差
Figure 449195DEST_PATH_IMAGE033
计算校正后的差分反射率因子
Figure 406699DEST_PATH_IMAGE009
，具体如下：
Figure 752098DEST_PATH_IMAGE054
。
ZDR偏差与温度关系需要不断补充样本来更新，以保证其准确性和代表性，还可对雷达观测数据进行分析，分析的目的是筛选出一组可用于更新ZDR偏差与温度关系的样本数据。具体分析过程如下：
分别判断仰角是否大于仰角阈值
Figure 447391DEST_PATH_IMAGE055
、反射率因子Z是否大于反射率因子阈值
Figure 765109DEST_PATH_IMAGE037
、相关系数CC是否大于相关系数阈值
Figure 147374DEST_PATH_IMAGE056
、高度是否小于高度阈值
Figure 963758DEST_PATH_IMAGE057
。
若仰角大于仰角阈值
Figure 462741DEST_PATH_IMAGE058
，且反射率因子Z大于反射率因子阈值
Figure 166124DEST_PATH_IMAGE037
，且相 关系数CC大于相关系数阈值
Figure 731010DEST_PATH_IMAGE041
，且高度小于高度阈值
Figure 519843DEST_PATH_IMAGE059
，则将当前雷达观测 数据中的差分反射率因子
Figure 556938DEST_PATH_IMAGE034
作为当前的偏差
Figure 645986DEST_PATH_IMAGE033
，并将其与当前采集的机内温度
Figure 862333DEST_PATH_IMAGE060
和机外温度
Figure 653309DEST_PATH_IMAGE061
作为最新样本进行存储，并通过最新样本重新拟合更新
Figure 494095DEST_PATH_IMAGE062
。
上述仰角阈值
Figure 177929DEST_PATH_IMAGE058
可设置为双偏振天气雷达的天线扫描范围的最大值，也可 设置为默认值85°。反射率因子阈值
Figure 344337DEST_PATH_IMAGE037
可根据双偏振天气雷达的接收机的灵敏度进行 调整，且其大小与接收机的灵敏度呈反比，即接收机灵敏度越高，弱信号探测能力越好，则 阈值越小。相关系数阈值
Figure 842183DEST_PATH_IMAGE063
可直接设置默认值为0.95。高度阈值
Figure 486660DEST_PATH_IMAGE064
可根据本 地当前月份的气温零度层平均高度来调整，其中，气温零度层平均高度可由气候资料统计 而出，高度阈值
Figure 278861DEST_PATH_IMAGE065
可设置为本地当前的气温零度层的平均高度减0.5 km。
如图2所示，基于以上实施例，本领域技术人员可以轻易理解，本发明还提供了一种双偏振天气雷达ZDR偏差实时校正装置，包括网络通讯模块1、机内温度传感器2、机外温度传感器3和主控制板4。
其中，网络通讯模块1用以从双偏振天气雷达的信号处理机实时采集雷达观测数 据，雷达观测数据包括差分反射率因子
Figure 350592DEST_PATH_IMAGE001
，此外，雷达观测数据一般还包括方位角、仰 角、时间、库长、库数、高度等参数信息，以及反射率因子Z、差分反射率因子ZDR、相关系数CC 等数据信息。网络通讯模块1可基于TCP/IP协议与双偏振天气雷达的信号处理机通讯。
机内温度传感器2用以采集双偏振天气雷达的机内温度
Figure 601313DEST_PATH_IMAGE066
，机内温度传感器2 可安装于雷达机内靠近雷达接收机的位置。机外温度传感器3用以采集双偏振天气雷达的 机外温度
Figure 518322DEST_PATH_IMAGE067
，机外温度传感器3可安装在安装于机外的雷达天线底座上。作为优选实施 例，机内温度传感器2和机外温度传感器3的测量精度≤0.1K，输出频率≥1Hz，可基于RS- 232 TTL端口与主控制板4通讯，为主控制板4提供高频率的温度测量信息。
主控制板4用以根据雷达观测数据、机内温度
Figure 160437DEST_PATH_IMAGE002
和机外温度
Figure 403069DEST_PATH_IMAGE003
计算偏差
Figure 406666DEST_PATH_IMAGE048
，具体如下：
Figure 392945DEST_PATH_IMAGE068
其中，
Figure 652018DEST_PATH_IMAGE050
分别为根据多组样本数据拟合出的系数；在投入使用前，一般 需要不少于100个初始样本数据来拟合得出初始的
Figure 799971DEST_PATH_IMAGE050
，在投入使用后即可根据
Figure 759706DEST_PATH_IMAGE069
以及机内温度
Figure 273644DEST_PATH_IMAGE070
和机外温度
Figure 182695DEST_PATH_IMAGE003
计算得出偏差
Figure 750817DEST_PATH_IMAGE019
。
主控制板4还用以根据偏差
Figure 729006DEST_PATH_IMAGE052
生成修正ZDR偏差指令，并控制网络通讯模块 1将所述修正ZDR偏差指令发送至双偏振天气雷达，以使双偏振天气雷达根据所述雷达观测 数据中的差分反射率因子
Figure 785650DEST_PATH_IMAGE018
和偏差
Figure 532895DEST_PATH_IMAGE052
计算校正后的差分反射率因子
Figure 288230DEST_PATH_IMAGE071
， 具体如下：
Figure 956977DEST_PATH_IMAGE010
。
ZDR偏差与温度关系需要不断补充样本来更新，以保证其准确性和代表性，还可对 雷达观测数据进行分析，分析的目的是筛选出一组可用于更新ZDR偏差与温度关系的样本 数据。主控制板还用以分别判断仰角是否大于仰角阈值
Figure 829031DEST_PATH_IMAGE040
、反射率因子Z是否大于反 射率因子阈值
Figure 165203DEST_PATH_IMAGE023
、相关系数CC是否大于相关系数阈值
Figure 622598DEST_PATH_IMAGE072
、高度是否小于高度 阈值
Figure 778642DEST_PATH_IMAGE042
。
若仰角大于仰角阈值
Figure 449793DEST_PATH_IMAGE040
，且反射率因子Z大于反射率因子阈值
Figure 906051DEST_PATH_IMAGE023
，且相 关系数CC大于相关系数阈值
Figure 268768DEST_PATH_IMAGE016
，且高度小于高度阈值
Figure 912108DEST_PATH_IMAGE042
，则将当前雷达观 测数据中的差分反射率因子
Figure 125963DEST_PATH_IMAGE073
作为当前的偏差
Figure 171148DEST_PATH_IMAGE052
，并将其与当前采集的机内温 度
Figure 704767DEST_PATH_IMAGE046
和机外温度
Figure 100982DEST_PATH_IMAGE067
作为最新样本，并利用最新样本重新拟合更新
Figure 106810DEST_PATH_IMAGE021
；
本发明实施例还包括存储器5，存储器5用来对最新样本进行存储。
上述仰角阈值
Figure 255769DEST_PATH_IMAGE074
可设置为双偏振天气雷达的天线扫描范围的最大值，也可 设置为默认值85°。反射率因子阈值
Figure 429130DEST_PATH_IMAGE037
可根据双偏振天气雷达的接收机的灵敏度进行 调整，且其大小与接收机的灵敏度呈反比，即接收机灵敏度越高，弱信号探测能力越好，则 阈值越小。相关系数阈值
Figure 318501DEST_PATH_IMAGE063
可直接设置默认值为0.95。高度阈值
Figure 133879DEST_PATH_IMAGE064
可根据本 地当前月份的气温零度层平均高度来调整，其中，气温零度层平均高度可由气候资料统计 而出，高度阈值
Figure 419236DEST_PATH_IMAGE065
可设置为本地当前的气温零度层的平均高度减0.5 km。
以上仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，其它未具体描述的部分，属于现有技术或公知常识。在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种双偏振天气雷达ZDR偏差实时校正方法，其特征在于，包括：
从双偏振天气雷达的信号处理机实时采集雷达观测数据，所述雷达观测数据包括差分反射率因子
Figure QLYQS_1
；
实时采集双偏振天气雷达的机内温度
Figure QLYQS_2
和机外温度
Figure QLYQS_3
，所述机内温度
Figure QLYQS_4
为雷达机内雷达接收机的位置的温度，所述机外温度
Figure QLYQS_5
为雷达天线底座外侧的温度；
根据所述雷达观测数据、机内温度
Figure QLYQS_6
和机外温度
Figure QLYQS_7
计算偏差
Figure QLYQS_8
，具体如下：
Figure QLYQS_9
其中，
Figure QLYQS_10
分别为根据样本数据拟合出的系数；
根据所述偏差
Figure QLYQS_11
生成修正ZDR偏差指令，并将所述修正ZDR偏差指令发送至双偏振天气雷达，以使所述双偏振天气雷达根据所述雷达观测数据中的差分反射率因子
Figure QLYQS_12
和偏差
Figure QLYQS_13
计算校正后的差分反射率因子
Figure QLYQS_14
，具体如下：
Figure QLYQS_15
。
2.根据权利要求1所述的一种双偏振天气雷达ZDR偏差实时校正方法，其特征在于，所述雷达观测数据还包括仰角、高度、反射率因子Z和相关系数CC；
所述方法还包括：
分别判断所述仰角是否大于仰角阈值
Figure QLYQS_16
、所述反射率因子Z是否大于反射率因子阈值
Figure QLYQS_17
、所述相关系数CC是否大于相关系数阈值
Figure QLYQS_18
、所述高度是否小于高度阈值
Figure QLYQS_19
；
若所述仰角大于仰角阈值
Figure QLYQS_21
，且所述反射率因子Z大于反射率因子阈值
Figure QLYQS_22
，且所述相关系数CC大于相关系数阈值
Figure QLYQS_24
，且所述高度小于高度阈值
Figure QLYQS_25
，则将当前雷达观测数据中的差分反射率因子
Figure QLYQS_26
作为当前的偏差
Figure QLYQS_27
，并将其与当前采集的机内温度
Figure QLYQS_28
和机外温度
Figure QLYQS_20
作为最新样本进行存储，并通过所述最新样本重新拟合更新
Figure QLYQS_23
。
3.根据权利要求2所述的一种双偏振天气雷达ZDR偏差实时校正方法，其特征在于，所述仰角阈值
Figure QLYQS_29
为双偏振天气雷达的天线扫描范围的最大值。
4.根据权利要求2所述的一种双偏振天气雷达ZDR偏差实时校正方法，其特征在于，所述反射率因子阈值
Figure QLYQS_30
根据双偏振天气雷达的接收机的灵敏度进行调整，且其大小与接收机的灵敏度呈反比。
5.根据权利要求2所述的一种双偏振天气雷达ZDR偏差实时校正方法，其特征在于，所述高度阈值
Figure QLYQS_31
为本地当前的气温零度层的平均高度减0.5 km。
6.一种双偏振天气雷达ZDR偏差实时校正装置，其特征在于，包括：
网络通讯模块，用以从双偏振天气雷达的信号处理机实时采集雷达观测数据，所述雷达观测数据包括差分反射率因子
Figure QLYQS_32
；
机内温度传感器和机外温度传感器，分别用以采集双偏振天气雷达的机内温度
Figure QLYQS_33
和机外温度
Figure QLYQS_34
，所述机内温度
Figure QLYQS_35
为雷达机内雷达接收机的位置的温度，所述机外温度
Figure QLYQS_36
为雷达天线底座外侧的温度；
主控制板，用以根据所述雷达观测数据、机内温度
Figure QLYQS_37
和机外温度
Figure QLYQS_38
计算偏差
Figure QLYQS_39
，具体如下：
Figure QLYQS_40
其中，
Figure QLYQS_41
分别为根据样本数据拟合出的系数；
所述主控制板还用以根据所述偏差
Figure QLYQS_42
生成修正ZDR偏差指令，并控制网络通讯模块将所述修正ZDR偏差指令发送至双偏振天气雷达，以使所述双偏振天气雷达根据所述雷达观测数据中的差分反射率因子
Figure QLYQS_43
和偏差
Figure QLYQS_44
计算校正后的差分反射率因子
Figure QLYQS_45
，具体如下：
Figure QLYQS_46
。
7.根据权利要求6所述的一种双偏振天气雷达ZDR偏差实时校正装置，其特征在于，所述雷达观测数据还包括仰角、高度、反射率因子Z和相关系数CC；
所述主控制板还用以分别判断所述仰角是否大于仰角阈值
Figure QLYQS_47
、所述反射率因子Z是否大于反射率因子阈值
Figure QLYQS_48
、所述相关系数CC是否大于相关系数阈值
Figure QLYQS_49
、所述高度是否小于高度阈值
Figure QLYQS_50
；
若所述仰角大于仰角阈值
Figure QLYQS_51
，且所述反射率因子Z大于反射率因子阈值
Figure QLYQS_53
，且所述相关系数CC大于相关系数阈值
Figure QLYQS_55
，且所述高度小于高度阈值
Figure QLYQS_56
，则将当前雷达观测数据中的差分反射率因子
Figure QLYQS_57
作为当前的偏差
Figure QLYQS_58
，并将其与当前采集的机内温度
Figure QLYQS_59
和机外温度
Figure QLYQS_52
作为最新样本，并利用最新样本重新拟合更新
Figure QLYQS_54
；
还包括存储器，用以对所述最新样本进行存储。
8.根据权利要求7所述的一种双偏振天气雷达ZDR偏差实时校正装置，其特征在于，所述仰角阈值
Figure QLYQS_60
为双偏振天气雷达的天线扫描范围的最大值。
9.根据权利要求7所述的一种双偏振天气雷达ZDR偏差实时校正装置，其特征在于，所述反射率因子阈值
Figure QLYQS_61
根据双偏振天气雷达的接收机的灵敏度进行调整，且其大小与接收机的灵敏度呈反比。
10.根据权利要求7所述的一种双偏振天气雷达ZDR偏差实时校正装置，其特征在于，所述高度阈值
Figure QLYQS_62
为本地当前的气温零度层的平均高度减0.5 km。
CN202211609704.XA 2022-12-15 2022-12-15 一种双偏振天气雷达zdr偏差实时校正方法和装置 Active CN115656947B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211609704.XA CN115656947B (zh) 2022-12-15 2022-12-15 一种双偏振天气雷达zdr偏差实时校正方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211609704.XA CN115656947B (zh) 2022-12-15 2022-12-15 一种双偏振天气雷达zdr偏差实时校正方法和装置

Publications (2)

Publication Number Publication Date
CN115656947A CN115656947A (zh) 2023-01-31
CN115656947B true CN115656947B (zh) 2023-03-14

Family

ID=85023293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211609704.XA Active CN115656947B (zh) 2022-12-15 2022-12-15 一种双偏振天气雷达zdr偏差实时校正方法和装置

Country Status (1)

Country Link
CN (1) CN115656947B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105137407B (zh) * 2015-10-08 2017-06-16 南京信息工程大学 一种双偏振天气雷达zdr在线标定方法及装置
CN205015473U (zh) * 2015-10-08 2016-02-03 南京信息工程大学 一种双偏振天气雷达zdr在线标定装置
CN113740934B (zh) * 2021-08-18 2022-12-09 浙江省大气探测技术保障中心 一种基于s波段双偏振天气雷达降水估测方法
CN114509734B (zh) * 2022-04-19 2022-07-19 南京气象科技创新研究院 一种基于雨滴谱的双偏振天气雷达数据质量实时评估方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Kenji Nakamura等.Measurement Method for Specific Attenuation in the Melting Layer Using a Dual Ka-Band Radar System.2018,第3511-3519页. *
李斐斐 ; 邵楠 ; 虞海锋 ; 步志超 ; 郭泽勇 ; .双偏振天气雷达双通道一致性分析与校准.2020,(03),第215-219页. *
赵世颖 ; 李柏 ; 陈晓辉 ; 王振会 ; .基于交叉-平行法的双偏振雷达差分反射率硬件定标.2015,(05),第5-12页. *
邵楠 ; 步志超 ; 李柏 ; 陈玉宝 ; .清远连州CINRAD/SA-D型业务试验双偏振天气雷达技术评估.2018,(09),第77-82页. *

Also Published As

Publication number Publication date
CN115656947A (zh) 2023-01-31

Similar Documents

Publication Publication Date Title
US7353690B2 (en) Atmospheric refractivity profiling apparatus and methods
Le Bouar et al. Validation of the rain profiling algorithm “ZPHI” from the C-band polarimetric weather radar in Darwin
CN106405558B (zh) 基于双频天线的地基气象微波辐射计系统及测量方法
CN109596594A (zh) 基于拉曼-米散射激光雷达的气溶胶消光系数反演方法
Gourley et al. Absolute calibration of radar reflectivity using redundancy of the polarization observations and implied constraints on drop shapes
Bringi et al. Rainfall estimation from C-band polarimetric radar in Okinawa, Japan: Comparisons with 2D-video disdrometer and 400 MHz wind profiler
CN103293116B (zh) 一种微脉冲差分吸收激光雷达水汽时空分布自动连续探测装置
CN108562882B (zh) 一种星载sar影像几何交叉定标方法和系统
CN105137407A (zh) 一种双偏振天气雷达zdr在线标定方法及装置
CN108827878B (zh) 一种地表气压的被动微波遥感探测方法
US20210247292A1 (en) Network assisted particulate matter sensor
FR3096145A1 (fr) Procédé et dispositif de mesure de précipitations
CN112363129A (zh) 一种天气雷达差分反射率因子参量定标方法
CN103278197A (zh) 一种基于车载系统的作物生长监测装置与方法
CN107576963A (zh) 基于粒子滤波的双偏振雷达差分传播相移的估计方法
CN114509734A (zh) 一种基于雨滴谱的双偏振天气雷达数据质量实时评估方法
CN109633573A (zh) 一种对大幅宽范围内湿对流层延时误差的校正方法
CN108876172B (zh) 一种基于改进型modis植被供水指数的地表土壤含水量评估方法
CN111679293A (zh) 一种激光雷达质量控制方法
CN115656947B (zh) 一种双偏振天气雷达zdr偏差实时校正方法和装置
CN107340519B (zh) 一种气象雷达降雨探测分析系统及方法
CN110736704B (zh) 土壤水与蒸发比耦合模拟及相互转化的方法
CN115375036A (zh) 遥感、光能利用率模型与气象融合的作物成熟期预测方法
CN209086157U (zh) 一种植物水分检测系统
CN113484918B (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
GR01 Patent grant
GR01 Patent grant