CN113640803B - 一种基于回波强度和回波顶高外推的短时定量降水预报方法 - Google Patents

一种基于回波强度和回波顶高外推的短时定量降水预报方法 Download PDF

Info

Publication number
CN113640803B
CN113640803B CN202111023555.4A CN202111023555A CN113640803B CN 113640803 B CN113640803 B CN 113640803B CN 202111023555 A CN202111023555 A CN 202111023555A CN 113640803 B CN113640803 B CN 113640803B
Authority
CN
China
Prior art keywords
radar
echo
vector
data
height
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
CN202111023555.4A
Other languages
English (en)
Other versions
CN113640803A (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.)
Jiangxi Normal University
Original Assignee
Jiangxi Normal University
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 Jiangxi Normal University filed Critical Jiangxi Normal University
Priority to CN202111023555.4A priority Critical patent/CN113640803B/zh
Publication of CN113640803A publication Critical patent/CN113640803A/zh
Application granted granted Critical
Publication of CN113640803B publication Critical patent/CN113640803B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/958—Theoretical aspects
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414—Discriminating targets with respect to background clutter
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/415—Identification of targets based on measurements of movement associated with the target
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418—Theoretical aspects
    • 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)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及大气科学领域，具体是一种基于回波强度和回波顶高外推的短时定量降水预报方法，用于基于雷达的短时定量降水预报。为了达到上述目的，本发明所采用的技术方案是：雷达基数据处理，剔除雷达回波观测中的非气象回波，并将数据从极坐标格式转换至笛卡尔直角坐标格式；雷达回波外推，对雷达回波强度和雷达回波顶高同时进行外推，获得两者的预报值；雷达降水反演，基于实时雷达回波强度、雷达回波顶高和降水观测资料，拟合出不同回波顶高的Z‑R关系表达式；雷达降水预报，基于雷达回波强度和回波顶高的预报值，利用实时拟合出的不同回波顶高的Z‑R关系表达式，对外推预报的回波强度进行降水反演，获得短时定量降水预报。

Description

一种基于回波强度和回波顶高外推的短时定量降水预报方法
技术领域
本发明涉及大气科学领域，尤其是一种基于回波强度和回波顶高外推的短时定量降水预报方法。
背景技术
多普勒天气雷达资料有很高的时空分辨率，在中尺度气象业务和研究中具有其他观测资料无可比拟的优势。由于数值模式存在spin-up问题(即在非平衡初值或扰动条件下，模式需要进行调整而达到平衡态的过程)，短时(2小时内)降水的预报效果不佳。目前，短时降水的业务预报主要依赖于雷达回波的外推预报。
目前，基于雷达回波外推的短时降水预报方法主要都是先采用不同技术(如交叉相关跟踪法TREC、基于连续方程约束的COTREC、通过Barnes空间滤波处理的BTREC、光流法等)获得雷达回波的移动矢量。再在假定回波强度短时间内无明显变化的条件下，对雷达回波进行外推，获得未来2小时内不同时刻的回波预报值。最后再利用降水强度R与雷达回波强度Z的经验公式Z＝aRb(即Z-R关系)反算(反演)出未来2小时内不同时刻的雷达降水量，即预报的降水。
可见，基于雷达回波外推的短时降水预报精度，不仅受雷达回波外推的精度影响，还要受外推回波降水反演的精度影响。以往的雷达降水预报方法多集中在雷达回波外推精度的改进上，如TREC、COTREC、BTREC和光流法等，对外推回波的降水反演则多采用固定Z-R关系(即Z＝300R1.4)算法。然而，由于Z-R关系中的系数a和b主要由降水雨滴普的特征决定，且随天气系统、降水类型、地理位置、季节等因素的变化而变化。因此，即使在雷达回波外推预报精度较高的情况下，用固定Z-R关系也很难准确地获得不同地区不同季节的雷达降水预报值。近年来，有学者也基于Z-R关系中的系数a和b随季节、降水类型、地理位置等变化的特性，发展一些精度更高的Z-R关系算法，如动态Z-R关系算法、基于回波顶高分级的动态Z-R关系算法等，但这些算法不仅需要用最优法进行实时拟合，还需要额外的预报参数(如回波顶高等)，因此，很少被用于基于雷达的短时降水预报上。
发明内容
(一)解决的技术问题
本发明的目的是克服现有技术的不足，本发明提供了一种基于回波强度和回波顶高外推的短时定量降水预报方法，通过对雷达回波强度(可反映云团中水凝物的浓度)和雷达回波顶高(可反映云团内部的上升运动)的共同外推，获得两者的预报值，解决基于回波高分级动态Z-R关系算法在雷达降水预报上的应用问题，有利于提升短时降水的定量预报精度。
(二)技术方案
本发明的技术方案：一种基于回波强度和回波顶高外推的短时定量降水预报方法，用于基于雷达的短时定量降水预报。为了达到上述目的，本发明所采用的技术方案是：雷达基数据处理，剔除雷达回波观测中的非气象回波，并将数据从极坐标格式转换至笛卡尔直角坐标格式；雷达回波外推，对雷达回波强度和雷达回波顶高同时进行外推，获得两者的预报值；雷达降水反演，基于实时雷达回波强度、雷达回波顶高和降水观测资料，拟合出不同回波顶高的Z-R关系表达式；雷达降水预报，基于雷达回波强度和回波顶高的预报值，利用实时拟合出的不同回波顶高的Z-R关系表达式，对外推预报的回波强度进行降水反演，获得短时定量降水预报。
具体包括如下步骤：
S1：雷达基数据处理：对雷达基数据进行解码、质量控制和坐标转换，获得笛卡尔坐标的雷达回波资料。
步骤S2：雷达回波外推：基于网格化的雷达组合反射率数据，通过交叉相关跟踪(TREC)法获得雷达回波移动矢量(即TREC矢量)，再用Barnes空间滤波更重矢量中的噪声和不协调，获得BTREC矢量，最后再用后向外推方案对雷达回组合反射率和雷达回波顶高进行外推，获得未来2小时内雷达回波强度和回波顶高的预报值。
步骤S3：雷达降水反演(基于回波顶高分级的Z-R关系)：将经过质控处理和网格化处理后的雷达组合反射率和雷达回波顶高数据，用双线性插值法插值至降水观测站处。根据回波顶高，按1km的间距将雷达组合反射率和观测降水分成不同类，如回波顶高在0-1km、1-2km、2-3km、3-4km、4-5km、5-6km、6-7km、7-8km、8-9km、9-10km、10-11km、11-12km、12-13km、13-14km、14-15km、＞15km范围内的降水与组合反射率，再用最优法拟合出不同类的最优Z-R关系。最优法的公式为
Figure BDA0003240691920000031
其中：n为观测降水的样本数，Gi为第i个观测站的观测降水量，
Figure BDA0003240691920000032
为第i个观测站处的雷达反演降水量，Zi为第i个观测站的雷达组合反射率，aj＝1,2，……，1200，bk＝100,101，……，300，不断调整aj和bk，使CTF值最小，此时aj和bk对应的Z-R关系即为最优Z-R关系。
步骤S4：雷达降水预报。基于外推的未来2小时内不同时次的雷达组合反射率和回波顶高数据，利用当前时次(动态)拟合的不同回波顶高的Z-R关系式，计算出未来2小时内不同时刻不同地点的定量降水预报值。
本发明技术方案中，所述S1具体内容包括：
S11：雷达基数据解码：根据雷达的型号(如CINRADSA/SB/SC/CB/CC/CD等)及其基数据的存储格式，对雷达基数据进行解码，获得不同仰角的极坐标格式的雷达回波资料。
S12：雷达回波质量控制：先利用回波(反射率)垂直梯度RGDZ＝W(Z-Zup)≥20dBZ剔除雷达回波观测中的异常地物杂波，其中：权重系数W在雷达斜距为0-40km时等于1，在40-200km间线性递减至0，Z为当前仰角的雷达回波(dBZ)，Zup为上一层仰角的雷达回波(dBZ)。再利用模糊逻辑算法Px＝N/25≤0.75剔除孤立非气象回波，其中N为方位径向平面(即雷达在某仰角上的扫描平面)上x点周围5×5窗口中雷达有效观测的数量。
S13：组合反射率和回波顶高的计算：基于质量控制处理后的不同仰角的极坐标格式的雷达反射率资料，计算出组合反射率(不同仰角层的最大值)和雷达回波顶高(即，≥18dBZ反射率因子的最高高度)，雷达回波顶高(相对雷达站)的计算公式为h＝Rsinθ+R2cos2θ/(2Re)，其中R为雷达探测的斜距，Re＝8500km为标准大气下的等效地球半径，θ为雷达探测的≥18dBZ反射率因子的最高仰角。
S14：雷达资料的坐标转换：先以雷达测站为中心，构建一个600×600个格点、水平分辨率为1km×1km的笛卡尔直角坐标网格(如图2虚线网格所示)，再根据雷达波束水平距离计算公式
Figure BDA0003240691920000041
其中hr为雷达站的海拔高度，计算出雷达斜距在水平面上的投影，即图2中的L。最后，再通过x＝L·sinβ和y＝L·cosβ将极坐标系(L，β)的雷达资料转至笛卡尔直角坐标系，其中β为雷达观测的方位角。
本发明技术方案中，所述S2具体内容包括：
S21：TREC矢量的计算。将t时刻网格化(格点数600×600，水平分辨率1km×1km)的雷达组合反射率数据分成40×40个没有重叠的小区域，其中每个小区域的格点数为15×15，并将每个小区域(初始小区域)的雷达反射率数据按先x方向系后y方向的方式转换成一个一维序列Z1。将Z1序列与t+Δt(默认为12分钟)时刻的雷达组合反射率数据中任意15×15个格点组成的Z2序列进行空间相关分析，其相关系数计算公式为
Figure BDA0003240691920000051
并将相关系数最大的Z2序列对应的小区域作为初始小区域的目标区域，从初始小区域中心至目标区域中心的矢量便是TREC矢量，反映雷达回波的移动情况，如图3所示。
S22：BTREC矢量的计算。由于回波小尺度的变化及跟踪失败等原因，TREC矢量中会有一些噪声或不协调的矢量，显著地影响着外推回波的精度。将TREC矢量在x方向和y方向的分量分别作为变量F(x,y)输入公式
Figure BDA0003240691920000052
和
Figure BDA0003240691920000053
获得BTREC矢量(剔除了噪声和更正了不协调)的分量F1(x,y)。其中，x和y是雷达回波在笛卡尔直角坐标系中的位置，高斯参数
Figure BDA0003240691920000054
和
Figure BDA0003240691920000055
M为距离(x,y)点100km范围内TREC矢量的格点数，rk为点(x,y)和点(xk,yk)之间的距离，滤波参数G＝0.35和C＝300。
S23：雷达回波外推。先利用双线性插值法将15km×15km水平分辨率的BTREC矢量插值至雷达组合反射率和雷达回波顶高的网格(1km×1km)上，再用雷达回波后向外推法进行雷达组合反射率和回波顶高的外推，获得未来2小时内任意时刻(如逐6分钟)的雷达组合反射率和雷达回波顶高的外推预报值。雷达回波后向外推法为Zi,j(t+Δt)＝Zi-Δi,j-Δj(t)，其中Z为雷达组合反射率或回波顶高，i和j分别为某一格点在x和y方向的序号，Δt为外推预报时间(如6分钟，半小时或2小时等)，Δi＝int(u·Δt+0.5)为预报时间内该格点在x方向上的移动距离，Δj＝int(v·Δt+0.5)格为预报时间内该格点在y方向上的移动距离，u和v分别为BTREC矢量在x和y方向上的分量。
(三)有益效果
本发明的优点在于：基于雷达回波强度和回波顶高的共同外推，实现了在雷达降水预报中引入了能反映云团内部上升运动的回波顶高，有利于提升短时降水定量预报的精度。
附图说明
图1是本发明方法流程图；
图2是雷达极坐标资料网格化示意图；
图3是TREC矢量计算示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
请参阅图1-图3，一种基于回波强度和回波顶高外推的短时定量降水预报方法，如图1所示，包括以下步骤：
S1：雷达基数据处理。对雷达基数据进行解码、质量控制和坐标转换，获得笛卡尔坐标的雷达回波资料，具体内容包括：
S11：雷达基数据解码。根据雷达的型号(如CINRAD SA/SB/SC/CB/CC/CD等)及其基数据的存储格式，对雷达基数据进行解码，获得不同仰角的极坐标格式的雷达回波资料。
S12：雷达回波质量控制。先利用回波(反射率)垂直梯度RGDZ＝W(Z-Zup)≥20dBZ剔除雷达回波观测中的异常地物杂波，其中权重系数W在雷达斜距为0-40km时等于1，在40-200km间线性递减至0，Z为当前仰角的雷达回波(dBZ)，Zup为上一层仰角的雷达回波(dBZ)。再利用模糊逻辑算法Px＝N/25≤0.75剔除孤立非气象回波，其中N为方位径向平面(即雷达在某仰角上的扫描平面)上x点周围5×5窗口中雷达有效观测的数量。
S13：组合反射率和回波顶高的计算。基于质量控制处理后的不同仰角的极坐标格式的雷达反射率资料，计算出组合反射率(不同仰角层的最大值)和雷达回波顶高(即，≥18dBZ反射率因子的最高高度)，雷达回波顶高(相对雷达站)的计算公式为h＝Rsinθ+R2cos2θ/(2Re)，其中R为雷达探测的斜距，Re＝8500km为标准大气下的等效地球半径，θ为雷达探测的≥18dBZ反射率因子的最高仰角。
S14：雷达资料的坐标转换。先以雷达测站为中心，构建一个600×600个格点、水平分辨率为1km×1km的笛卡尔直角坐标网格(如图2虚线网格所示，其中虚线网格为笛卡尔坐标系；椭圆粗实线为雷达极坐标；x、y分别为雷达极坐标观测在笛卡尔坐标的投影；A为雷达观测点；L为雷达斜距在平面上的投影；β为雷达观测的方位角。)，再根据雷达波束水平距离计算公式
Figure BDA0003240691920000071
其中hr为雷达站的海拔高度，计算出雷达斜距在水平面上的投影，即图2中的L。最后，再通过x＝L·sinβ和y＝L·cosβ将极坐标系(L，β)的雷达资料转至笛卡尔直角坐标系，其中β为雷达观测的方位角。
S2：雷达回波外推。基于网格化的雷达组合反射率数据，通过交叉相关跟踪(TREC)法获得雷达回波移动矢量(即TREC矢量)，再用Barnes空间滤波更重矢量中的噪声和不协调，获得BTREC矢量，最后再用后向外推方案对雷达回组合反射率和雷达回波顶高进行外推，获得未来2小时内雷达回波强度和回波顶高的预报值，具体内容包括：
S21：TREC矢量的计算。将t时刻网格化(格点数600×600，水平分辨率1km×1km)的雷达组合反射率数据分成40×40个没有重叠的小区域，其中每个小区域的格点数为15×15，并将每个小区域(初始小区域)的雷达反射率数据按先x方向系后y方向的方式转换成一个一维序列Z1。将Z1序列与t+Δt(默认为12分钟)时刻的雷达组合反射率数据中任意15×15个格点组成的Z2序列进行空间相关分析，其相关系数计算公式为
Figure BDA0003240691920000081
并将相关系数最大的Z2序列对应的小区域作为初始小区域的目标区域，从初始小区域中心至目标区域中心的矢量便是TREC矢量，反映雷达回波的移动情况，如图3所示。
S22：BTREC矢量的计算。由于回波小尺度的变化及跟踪失败等原因，TREC矢量中会有一些噪声或不协调的矢量，显著地影响着外推回波的精度。将TREC矢量在x方向和y方向的分量分别作为变量F(x,y)输入公式
Figure BDA0003240691920000082
和
Figure BDA0003240691920000083
获得BTREC矢量(剔除了噪声和更正了不协调)的分量F1(x,y)。其中，x和y是雷达回波在笛卡尔直角坐标系中的位置，高斯参数
Figure BDA0003240691920000084
和
Figure BDA0003240691920000085
分，M为距离(x,y)点100km范围内TREC矢量的格点数，rk为点(x,y)和点(xk,yk)之间的距离，滤波参数G＝0.35和C＝300。
S23：雷达回波外推。先利用双线性插值法将15km×15km水平分辨率的BTREC矢量插值至雷达组合反射率和雷达回波顶高的网格(1km×1km)上，再用雷达回波后向外推法进行雷达组合反射率和回波顶高的外推，获得未来2小时内任意时刻(如逐6分钟)的雷达组合反射率和雷达回波顶高的外推预报值。雷达回波后向外推法为Zi,j(t+Δt)＝Zi-Δi,j-Δj(t)，其中Z为雷达组合反射率或回波顶高，i和j分别为某一格点在x和y方向的序号，Δt为外推预报时间(如6分钟，半小时或2小时等)，Δi＝int(u·Δt+0.5)为预报时间内该格点在x方向上的移动距离，Δj＝int(v·Δt+0.5)格为预报时间内该格点在y方向上的移动距离，u和v分别为BTREC矢量在x和y方向上的分量。
S3：雷达降水反演(基于回波顶高分级的Z-R关系)。将经过质控处理和网格化处理后的雷达组合反射率和雷达回波顶高数据，用双线性插值法插值至降水观测站处。根据回波顶高，按1km的间距将雷达组合反射率和观测降水分成不同类，如回波顶高在0-1km、1-2km、2-3km、3-4km、4-5km、5-6km、6-7km、7-8km、8-9km、9-10km、10-11km、11-12km、12-13km、13-14km、14-15km、＞15km范围内的降水与组合反射率，再用最优法拟合出不同类的最优Z-R关系。最优法的公式为
Figure BDA0003240691920000091
其中n为观测降水的样本数，Gi为第i个观测站的观测降水量，
Figure BDA0003240691920000092
为第i个观测站处的雷达反演降水量，Zi为第i个观测站的雷达组合反射率，aj＝1,2，……，1200，bk＝100,101，……，300，不断调整aj和bk，使CTF值最小，此时aj和bk对应的Z-R关系即为最优Z-R关系。
S4：雷达降水预报。基于外推的未来2小时内不同时次的雷达组合反射率和回波顶高数据，利用前一个时次(动态)拟合的不同回波顶高的Z-R关系式，计算出未来2小时内不同时刻不同地点的定量降水预报值。
尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的保护范围由所附权利要求及其等同物限定。

Claims (2)

1.一种基于回波强度和回波顶高外推的短时定量降水预报方法，其特征在于，包括以下步骤：
S1.雷达基数据处理，对雷达基数据进行解码、质量控制和坐标转换，获得笛卡尔坐标的雷达回波资料；
S2.雷达回波外推，基于网格化的雷达组合反射率数据，通过交叉相关跟踪TREC法获得雷达回波移动矢量TREC矢量，再用Barnes空间滤波更重矢量中的噪声和不协调，获得BTREC矢量，最后再用后向外推方案对雷达组合反射率和雷达回波顶高进行外推，获得未来2小时内雷达回波强度和回波顶高的预报值，所述S2的具体内容包括：
S21：TREC矢量的计算，将t时刻网格化的雷达组合反射率数据分成40×40个没有重叠的小区域，其中网格化的格点数为600×600，水平分辨率为1km×1km，每个小区域的格点数为15×15，并将每个小区域的雷达反射率数据按先x方向后y方向的方式转换成一个一维序列Z1，将Z1序列与t+Δt时刻的雷达组合反射率数据中任意15×15个格点组成的Z2序列进行空间相关分析，其相关系数计算公式为
Figure FDA0003699403080000011
并将相关系数最大的Z2序列对应的小区域作为初始小区域的目标区域，从初始小区域中心至目标区域中心的矢量便是TREC矢量，反映雷达回波的移动情况；
S22：BTREC矢量的计算，由于回波小尺度的变化及跟踪失败原因，TREC矢量中会有一些噪声或不协调的矢量，显著地影响着外推回波的精度，将TREC矢量在x方向和y方向的分量分别作为变量F(x,y)输入公式
Figure FDA0003699403080000012
获得BTREC矢量的分量F1(x,y)，其中：x和y是雷达回波在笛卡尔直角坐标系中的位置，
Figure FDA0003699403080000021
和
Figure FDA0003699403080000022
为高斯参数，M为距离(x,y)点100km范围内TREC矢量的有效格点数，rk为点(x,y)和点(xk,yk)之间的距离，滤波参数G＝0.35和C＝300；
S23：雷达回波外推，先利用双线性插值法将15km×15km水平分辨率的BTREC矢量插值至雷达组合反射率和雷达回波顶高的网格1km×1km上，再用雷达回波后向外推法进行雷达组合反射率和回波顶高的外推，获得未来2小时内任意时刻的雷达组合反射率和雷达回波顶高的外推预报值，雷达回波后向外推法为Zi,j(t+Δt)＝Zi-Δi,j-Δj(t)，其中：Z为雷达组合反射率或回波顶高，i和j分别为某一回波点在x和y方向的序号，Δt为外推预报时间，Δi＝int(u·Δt+0.5)为预报时间内回波点在x方向上的移动距离，Δj＝int(v·Δt+0.5)格为预报时间内该回波点在y方向上的移动距离，u和v分别为BTREC矢量在x和y方向上的分量；
S3.雷达降水反演基于回波顶高分级的Z-R关系，将经过质控处理和网格化处理后的雷达组合反射率和雷达回波顶高数据，用双线性插值法插值至降水观测站处，基于回波顶高，按1km的间距将雷达组合反射率和观测降水分成不同类，再用最优法拟合出不同类的最优Z-R关系；
S4：雷达降水预报，基于外推的未来2小时内不同时次的雷达组合反射率和回波顶高数据，利用当前时次动态拟合的不同回波顶高的Z-R关系式，计算出未来2小时内不同时刻不同地点的定量降水预报值。
2.根据权利要求1所述的一种基于回波强度和回波顶高外推的短时定量降水预报方法，其特征在于：所述S1的具体内容包括：
S11.雷达基数据解码，根据雷达的型号及其基数据的存储格式，对雷达基数据进行解码，获得不同仰角的极坐标格式的雷达回波资料；
S12.雷达回波质量控制，先利用回波反射率垂直梯度RGDZ＝W(Z-Zup)≥20dBZ剔除雷达回波观测中的异常地物杂波，其中权重系数W在雷达斜距为0-40km时等于1，在40-200km间线性递减至0，Z为当前仰角的雷达回波强度，Zup为上一层仰角的雷达回波强度，再利用模糊逻辑算法Px＝N/25≤0.75剔除孤立非气象回波，其中N为方位径向平面，即雷达在某仰角上的扫描平面上x点周围5×5窗口中雷达有效观测的数量；
S13.组合反射率和回波顶高的计算，基于质量控制处理后的不同仰角的极坐标格式的雷达反射率资料，计算出组合反射率不同仰角层的最大值和雷达回波顶高，即≥18dBZ反射率因子的最高高度，雷达回波顶高相对雷达站的计算公式为h＝R sinθ+R2cos2θ/(2Re)，其中R为雷达探测的斜距，Re＝8500km为标准大气下的等效地球半径，θ为雷达探测的≥18dBZ反射率因子的最高仰角；
S14.雷达资料的坐标转换，先以雷达测站为中心，构建一个600×600个格点、水平分辨率为1km×1km的笛卡尔直角坐标网格，再根据雷达波束水平距离计算公式
Figure FDA0003699403080000031
其中hr为雷达站的海拔高度，计算出雷达斜距在水平面上的投影，最后再通过x＝L·sinβ和y＝L·cosβ将极坐标系(L，β)的雷达资料转至笛卡尔直角坐标系，其中β为雷达观测的方位角。
CN202111023555.4A 2021-09-01 2021-09-01 一种基于回波强度和回波顶高外推的短时定量降水预报方法 Active CN113640803B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111023555.4A CN113640803B (zh) 2021-09-01 2021-09-01 一种基于回波强度和回波顶高外推的短时定量降水预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111023555.4A CN113640803B (zh) 2021-09-01 2021-09-01 一种基于回波强度和回波顶高外推的短时定量降水预报方法

Publications (2)

Publication Number Publication Date
CN113640803A CN113640803A (zh) 2021-11-12
CN113640803B true CN113640803B (zh) 2022-07-22

Family

ID=78424907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111023555.4A Active CN113640803B (zh) 2021-09-01 2021-09-01 一种基于回波强度和回波顶高外推的短时定量降水预报方法

Country Status (1)

Country Link
CN (1) CN113640803B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL2030904B1 (en) * 2022-02-11 2023-01-06 Meteorological Sciences Inst Of Jiangxi Province Radar echo extrapolation method based on cross tracking method and spatial filtering
CN114488352B (zh) * 2022-02-16 2022-11-15 长沙祥云瑞风信息技术有限公司 一种无缝隙精细化降水预报方法及系统
CN115792847B (zh) * 2022-11-08 2023-07-18 江西师范大学 一种基于神经网络和回波垂直信息的定量降水估测方法
CN117452369B (zh) * 2023-12-25 2024-04-05 江西师范大学 一种用于短时灾害性天气监测的回波顶高计算优化方法
CN118169780B (zh) * 2024-03-27 2024-09-24 中国气象局武汉暴雨研究所 针对不同降水回波结构的对流性短时强降水分型预警方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111897030A (zh) * 2020-07-17 2020-11-06 国网电力科学研究院有限公司 一种雷暴预警系统及方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5681575A (en) * 1992-05-19 1997-10-28 Westaim Technologies Inc. Anti-microbial coating for medical devices
FR2881436B1 (fr) * 2005-02-03 2007-04-27 Commissariat Energie Atomique Procede de determination de la diversite des lymphocytes t dans un echantillon biologique
CN104977584B (zh) * 2015-06-29 2017-08-15 深圳市气象台 一种对流天气的临近预报方法及系统
CN109975809B (zh) * 2019-04-26 2021-05-04 兰州大学 一种同化雷达及闪电资料的方法
CN110967695A (zh) * 2019-10-28 2020-04-07 兰州大方电子有限责任公司 一种基于深度学习的雷达回波外推短临预报方法
CN111650572B (zh) * 2020-07-23 2023-02-24 浪潮云信息技术股份公司 一种减小短时降水估测偏差的方法及系统
CN112116785B (zh) * 2020-08-21 2022-08-19 中国瑞林工程技术股份有限公司 一种基于强降雨气象预报下尾矿库灾害预警方法及装置
CN112363168B (zh) * 2021-01-13 2021-03-26 南京满星数据科技有限公司 一种基于雷达外推和模式预报的同化融合方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111897030A (zh) * 2020-07-17 2020-11-06 国网电力科学研究院有限公司 一种雷暴预警系统及方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An Algorithm to Enhance Nowcast of Rainfall Brought by Tropical Cyclones through Separation of Motions;W.C. Woo;《Tropical Cyclone Research and Review》;20140531;第3卷(第2期);111-121 *
河北承德地区冰雹天气雷达预警指标分析;张晓辉等;《黑龙江农业科学》;20180710(第07期);25-28 *

Also Published As

Publication number Publication date
CN113640803A (zh) 2021-11-12

Similar Documents

Publication Publication Date Title
CN113640803B (zh) 一种基于回波强度和回波顶高外推的短时定量降水预报方法
Heinselman et al. High-temporal-resolution capabilities of the national weather radar testbed phased-array radar
CN109100717B (zh) 一种多源微波遥感海面风场数据融合方法及其装置
CN108693534A (zh) Nriet x 波段雷达协同组网分析方法
Holleman et al. Analysis and correction of dual PRF velocity data
CN110879428B (zh) 一种复杂地形条件下相同波段雷达之间重叠区域组网测雨方法
Janssen et al. The shape of Doppler spectra from precipitation
CN112731564B (zh) 一种基于多普勒天气雷达数据的雷电智能预报方法
CN110161506A (zh) 一种基于多源气象观测资料的分类型冰雹落区识别方法
CN113900103A (zh) 突发性强对流灾害天气自适应快速识别预警改进算法
CN113311416B (zh) 山地小流域雷达定量降水估计方法
Torres et al. Adaptive-weather-surveillance and multifunction capabilities of the national weather radar testbed phased array radar
CN115236615B (zh) 基于t矩阵法的机载极化气象雷达降水粒子回波仿真方法
CN116953707B (zh) 一种潮位监测雷达装置及监测方法和系统
CN114898206A (zh) 一种短时强降水预报方法及计算机设备和存储介质
CN114296050B (zh) 基于激光雷达云图探测的光伏电站短期发电功率预测方法
CN111366904B (zh) 一种基于二维频谱相干合并的cfar检测方法
Kollias et al. Agile weather observations using a dual-polarization X-band phased array radar
CN112859083A (zh) 一种面向资料同化的风廓线雷达风场资料质量控制方法
CN114636981A (zh) 基于雷达回波的在线深度学习台风中心定位系统
CN113936166A (zh) 一种基于多普勒天气雷达数据的冰雹回波识别方法及系统
Yamada Characteristics of wind fields derived from the multiple-Doppler synthesis and continuity adjustment technique (MUSCAT)
Yang et al. Can CINRAD Radar with VCP-21 Mode Capture the Accumulated Rainfall Pattern and Intensity of Fast-moving Storms?
Williams et al. The ultra high resolution QuikSCAT product
Lamer et al. Evaluation of gridded Scanning ARM Cloud Radar reflectivity observations and vertical Doppler velocity retrievals.

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