多普勒天气雷达原理与应用6雷达探测算法3
合集下载
天气雷达的基本工作原理和参数-168页文档资料
多普勒天气雷达除常规天气雷达功能 之外，还可利用降水回波频率与发射频率 之间变化的信息来测定降水粒子的径向速 度，并通过此推断风速分布，垂直气流速 度，大气湍流，降水粒子谱分布，降水中 特别是强对流降水中风场结构特征。
常规天气雷达仅能提供反射率 因子资料。多普勒天气雷达将提供 两种附加的基本资料，径向速度和 速度谱宽，它们将增强对强风暴的 探测能力，也能改进对中尺度和天 气尺度系统的预报。
体扫模式 (VCP:Volume Cover Pattern） 扫描方式确定一次体积扫中使用多少个仰角，
而具体是哪些仰角则由体扫模式来规定。WSR-88D 可有20个不同的VCP，目前只定义了其中的4个： VCP11 -- VCP11（scan strategy #1，version 1） 规定5分钟内对14个具体仰角的扫描方式。 VCP21 -- VCP21（scan strategy #2，version 1） 规定6分钟内对9个具体仰角的扫描方式。 VCP31 --- VCP31 （scan strategy #3，version 1）规定10分钟内对5个具体仰角的扫描方式。 VCP32 --- VCP32（scan strategy #3，version 2）确定的10分钟完成的5个具体仰角与VCP31相同。 不同之处在于VCP31使用长雷达脉冲而VCP32使用 短脉冲。 WSR-98D未定义VCP32。
自相干多普勒天气雷达结构框图
全相干多普勒天气雷达结构框图
fo 发射脉冲的载频 fd 多普勒频率
发射频率 Vs 多普勒频移
发射频率 多普勒频移
中国新一代天气雷达系统简介
• 1、雷达数据采集系统（RDA） • 2、雷达产品生成子系统（RPG） • 3、主用户处理器子系统（PUP）
常规天气雷达仅能提供反射率 因子资料。多普勒天气雷达将提供 两种附加的基本资料，径向速度和 速度谱宽，它们将增强对强风暴的 探测能力，也能改进对中尺度和天 气尺度系统的预报。
体扫模式 (VCP:Volume Cover Pattern） 扫描方式确定一次体积扫中使用多少个仰角，
而具体是哪些仰角则由体扫模式来规定。WSR-88D 可有20个不同的VCP，目前只定义了其中的4个： VCP11 -- VCP11（scan strategy #1，version 1） 规定5分钟内对14个具体仰角的扫描方式。 VCP21 -- VCP21（scan strategy #2，version 1） 规定6分钟内对9个具体仰角的扫描方式。 VCP31 --- VCP31 （scan strategy #3，version 1）规定10分钟内对5个具体仰角的扫描方式。 VCP32 --- VCP32（scan strategy #3，version 2）确定的10分钟完成的5个具体仰角与VCP31相同。 不同之处在于VCP31使用长雷达脉冲而VCP32使用 短脉冲。 WSR-98D未定义VCP32。
自相干多普勒天气雷达结构框图
全相干多普勒天气雷达结构框图
fo 发射脉冲的载频 fd 多普勒频率
发射频率 Vs 多普勒频移
发射频率 多普勒频移
中国新一代天气雷达系统简介
• 1、雷达数据采集系统（RDA） • 2、雷达产品生成子系统（RPG） • 3、主用户处理器子系统（PUP）
多普勒天气雷达应用研究
多普勒天气雷达应用研究【摘要】本文简要介绍了多普勒天气雷达的相关原理，并结合典型个例，应用多普勒天气雷达回波资料及常规天气资料，阐述了多普勒图像在气象保障中的应用，简要论述复杂海岸地形产生的气流在对流降水过程中的影响。
【关键词】多普勒天气雷达；闪电强度；地形1 多普勒天气雷达的相关原理1.1 PPI显示方式雷达图像的PPI显示，是指雷达天线在一系列固定仰角上扫描360。
进行取样，并经过对目标物的数据进行分析、处理而得出的结果。
在每个仰角上，沿雷达波束向外径向距离增加，离地高度也增加。
因此，当环境风场只随高度变化时，雷达扫描一周便能揭示出从地面到雷达显示范围边缘高度上所有风的信息。
1.2 零值线的意义在分析多普勒图像时.关键是要寻找到零值线，然后围绕零线进行大气流场的分析。
零值线一种情况表明此处的风向与雷达探测的径向是垂直的；另一种情况是该处真实风速为零（也可能是速度极小或处于静止状态）。
当所有高度上的风速都一样，风向从地面上（雷达站）的南风均匀地改变到显示边缘高度上的正南风。
环境风场平面图：风速固定，在地面为南风（图像中心），均匀地经西南风变为图像边缘处的西风。
在显示区的外缘，当雷达指向正北和正南时，多普勒速度值为零，这意味着在相应高度不是正西风就是正东风。
由于在显示区的西部边缘多普勒速度值是正的（朝向雷达的分量），东部边缘的多普勒速度值是负的（离开雷达的分量），那么很明显.在雷达图像显示区边缘高度上风向是由西向东的。
1.3 典型流场的多普勒模式掌握典型流场的多普勒模式，对于分析复杂天气系统的流场结构有着重要的意义，下面主要介绍基本气流模式、暖切变流场模式和冷锋（冷切变）流场模式。
1.3.1 基本气流模式基本气流（水平面上风向风速一致、风速随高度先增加后减小的西南气流）的方向是从趋近（正值）中心吹向远离（负值）中心，并和零值线所在的向径方向垂直，这就是基本气流径向速度分布模式。
其它方向的基本气流的趋近区、远离区和零值线也随之而变，但图形一样。
2DU简介双偏振多普勒天气雷达原理与应用 PPT课件
（一）几种测雨方程的形式
1、常规天气雷达
（21）
2、双线偏振雷达
（22）
3、双线偏振多普勒雷达
（二）模拟雨滴谱资料的结果（参阅有关文献）
1、测雨方程形式与前面相同，但要确定方程中的一些系数。 2、不同测雨方程求得的雨强IDP、IDR、IDK，与雨滴谱确定的真值 雨强R之间的关系，见下面给出的图11 ―图14 。
（三）双程差分传播相位变化值φDP
1、φDP的含义：设水平及垂直偏振波通过相同长度 的一个降水区（可包含非球形粒子组成），散射 回天线处的相位分别为φHH及φVV，则定义： φDP=φHH -φVV=δ+ ɸdp
其中：
（1）δ为反射相位差。 对于同一粒子δ值是固定的。 对于S波段雷达，满足瑞利散射的粒子，δ是个小 量。但冰雹时δ值较大。
1、。ZDR定义式为 ：
ZDR=10log（ZHH/ZVV）=10log ZHH-10logZVV （8） 在信号处理器RVP8中用给出的求ZDR的公式为：
ZDR=10loS gHH=
S VV
2、ZDR值的范围：一般为-0.5—— +6.0dB （1）一般雨滴呈扁旋转椭球，ZHH＞ZVV，故常为ZDR＞0 。 （2）大雨滴时 ，呈更扁的椭球形，故ZDR值可达3~5dB。 （3）冰雹 ，由于翻转作用总体效果接近球形，ZDR值在零附近 ，可以是小的负值或小的正值。
里
即KDP是双程传播相位变化值φDP随距离的变化程度。 2、若 (rm) 与 (rn) 不是相邻两库的距离，而是相隔较远的两个库之
间的距离，则KDP代表该降水段上的平均值。 3、KDP值的大小：一般KDP＜1°/Km，但含有冰核的大雨滴，KDP
（五）双线偏振雷达的退极化因子LDR
多普勒天气雷达
工作原理
多普勒雷达是世界上最先进的天气监测设备，并且已经在很多国家得到深入应用，因此，下面我们就多普勒 雷达的工作原理进行深入分析和研究，以便能够使人们对其工作原理有着更为清楚的认识。
1.1通过气象目标对雷达电磁波的散射和吸收
粒子能够对电磁波进行吸收和散射，这也是粒子对电磁波的两大基本形式。雷达探测大气的基础是由气象目 标对雷达电磁波的吸收和散射所得。如果电磁波的波束在大气传播途中遇到包括云滴、雨滴以及其他悬浮粒子和 空气分子，作为入射的电磁波波束中的有一部分会因为上述的粒子反射到不同地方，这类现象称之为散射。一部 分散射的电磁波波束会被粒子吸收，最终按照雷达的方向返回被雷达天线接收，多普勒天气雷达能够通过接收到 的电磁波束中自带的振幅和位相等数据，得出气象目标的平均速度以及发射率因子和速度谱宽等基本数据，进而 推断并计算出相对应的气象情况和其他内部结构特征。
重要意义
多普勒雷达是世界上最先进的雷达系统，有“超级千里眼”之称。相较于传统天气雷达，多普勒雷达能够监 测到位于垂直地面8-12公里的高空中的对流云层的生成和变化，判断云的移动速度，其产品信息达72种，天气预 报的精确度比以前将会有较大提高。1991至1997年，美国在全国及海外布的165台NEXRDA被称为天气雷达系统的 典范，是世界上最先进的和最精确的天气雷达系统。它所采用的多普勒信号处理技术和自动产生灾害性天气警报 的能力无与伦比。NEXRAD可以自动形成和显示丰富多彩的天气产品，极大地提高了对超级单体、湖泊效应雪、成 层雪、雷暴、降水、风切变、下击暴流、龙卷、锋面、湍流、冰雹等重大灾害性天气的监测和预报能力。对强雷 暴的侦察率是96%，对龙卷的发现率是83%，对龙卷警告的平均预警时间是18分钟，而在未建NEXRDA络之前，美国 国家上述参数的平均值分别是60%，40%和2分钟。从中可以预料CINRDA将从根本上增强探测强雷暴的能力，能较 早地探测到晴空下威胁航行的大气湍流和发生灾害性洪水的可能，并为水资源的管理决策提供极有价值的信息。 新一代天气雷达系统建设是我国20世纪末21世纪初的一项气象现代化工程，计划在全国建成S频段和C频段雷达 156部，该系统建成后，我国的气象现代化水平会上一个新的台阶。
多普勒雷达探测原理
多普勒雷达探测原理8.1.1 多普勒效应多普勒效应是奥地利物理学家J.Doppler 1842年⾸先从运动着的发声源中发现的现象，定义为"当接收者或接收器与能量源处于相对运动状态时，能量到达接收者（器）时频率的变化"。
⼀个例⼦是：当⼀辆紧急的⽕车（汽车）鸣着喇叭以相当⾼的速度向着你驶来时，声⾳的⾳调（频率）由于波的压缩（较短波长）⽽增加。
当⽕车（汽车）远离你⽽去时，这声⾳的⾳调（频率）由于波的膨胀（较长波长）⽽减低。
多普勒频率(多普勒频移)：对于⼀个运动的⽬标，向着雷达运动或远离雷达运动所产⽣的频移量是相同的，但符号不同：①如果⽬标移向雷达频移为正；②如果⽬标远离雷达频移为负。
8.1.2 径向速度径向速度简单地定义为⽬标运动平⾏于雷达径向的分量。
它是⽬标运动沿雷达径向的分量，既可以向着雷达，也可以离开雷达。
需要注意：①径向速度总是⼩于或等于实际⽬标速度；②由WSR-88D测量的速度只是⽬标向着或离开雷达的运动；③当⽬标运动垂直于雷达径向或静⽌时径向速度为零。
⽬标的实际速度与WSR-88D描述的径向速度间的关系能⽤数学⽅法描述成径向速度⽅程│Vr│=│V│•cosβ其中Vr为径向速度，V为实际速度，β为实际速度V与雷达径向之间最⼩的夹⾓。
8.1.3 多普勒天⽓雷达测速由于多普勒频移(Hz)相对发射频率(MHz)很⼩，故多普勒天⽓雷达通常不是直接测量多普勒频移，⽽是通过测量相继返回的脉冲对之间的位相差来确定⽬标物的径向速度，这种脉冲位相的变化可以⽐较容易并且⽐较准确的测量。
这种测速技术叫做"脉冲对处理"。
脉冲对处理 Pulse-Pair Method要使多普勒雷达能够提取⽬标的多普勒运动信息，必须知道每个发射波的初相位，这样就可以⽐较相继返回信号的位相。
如果每个发射波的初位相不知道，那么将⽆法知道相继返回的两个脉冲间的相移，也就⽆法对⽬标物沿雷达径向做出估计。
雷达气象学之第三章(多普勒天气雷达探测原理和方法)
2、脉冲对处理法（PPP）
在一定假设条件下对每一个距离库内的连 续两个取样值作成对处理．从而获得平均 多普勒频率和频谱宽度。此法优点在于能 实时处理．并且有一定精度，但它不能得 到频率谱。
3、相干记忆滤波器（CMF）处理法
此法只需要一个线路，在不设置距离库的 情况下同时对雷达探测范围内各个距离上 作粗略的谱分析，并能用如PSI(平面切变 线是其)等直接显示出来。但它精度不高；
垂 直 风 廓 线
补充风符号
1.风向杆 表示风的 来向。 2.风羽每 条代表风 速4米/秒， 半条代表2 米/秒，三 角旗代表 20米/秒。
谱 宽
反 射 率
三、影响速度谱宽的气象因子
• 多普勒速度谱宽表征着有效照射体内不同 大小的多普勒速度偏离其平均值的程度， 实际上它是由散射粒子具有不同的径向速 度所引起的。对气象目标物而言，影响速 度谱宽的主要因子有四个：
• 显然，雷达有效照射体中粒子直径的差别 越大，由此造成的多普勒速度谱越宽。
• 因此速度的谱宽实际上也取决于降水粒子 的谱分布。
• 当雷达水平探测时，粒子的下落末速度在 雷达波轴上的径向分量为零，所以它对多 普勒速度谱宽没有任何影响。
• 而当雷达垂直指向探测时，粒子下落末速 度即为径向速度，故由此造成的谱曾宽作 用最大。
• 在实际工作中需要了解的是有效照射体内
平均的多普勒速度和速度谱宽度，根据以
上关系式，并注意到 f 2v 关系式，则平均
多普勒速度
v
,和速度谱方差
2 v
分别为：
v 1 v v dv
Pr
2 v
1 Pr
vv
2
v dv
径向速度谱密度、平均径向速度、径向速度 谱宽三者的关系示意图
雷达气象学原理多普勒天气雷达
雷达气象学原理 多普勒天气雷达
多普勒天气雷达除此之外，还可利用 降水回波频率与发射频率之间变化的信 息来测定降水粒子的径向速度，并通过 此推断风速分布，垂直气流速度，大气 湍流，降水粒子谱分布，降水中特别是 强对流降水中风场结构特征。
以前，用常规天气雷达进行的天气预报 仅仅使用反射率因子资料。多普勒天气雷达 将提供两种附加的基本资料，径向速度和速 度谱宽，它们将增强对强风暴的探测能力， 也能改进对中尺度和天气尺度系统的预报。
多普勒频率与径向速度的关系
假设多普勒雷达发射脉冲的工作频率为f0，目标与雷达的距
离为r，则雷达波发往目标到返回天线所经过的距离为2r。这 个距离用波长来度量，相当 个波长；用弧度来衡量相当于 个弧度。若所发射的电磁波在天线处的位相为 ，那么电磁波 被散射回到天线时的相位应是
位相的时间变化率
由于目标物的径向运动引起 的雷达回波信号的频率变化，它 就是多普频移或多普勒频率。
多普勒雷达是通过直接测量多普勒 频率来得到径向速度的吗?
4.2 多普勒雷达径向速度探测方法
Pulse-Pair Method 脉冲对方法
取两个连续的脉冲然 后测量接收脉冲的相位， 这种脉冲对位相变化可以 比较容易并且比较准确地 测量
DΦ/dt 实际上就是角 速度 = w = 2πfd
假定当第一个脉冲遇到目标物时，该目标物距雷达的距离为r，则该目标物 产生的回波到达雷达时的位相为：
2、平均多普勒频移及频谱宽度
3、平均多普勒速度和速度谱宽度
注意：脉冲对方法并没有从回波信号中提取频谱或功率谱，从而 不能按以上公式计算和，而是直接对回波信号作简便计算求得。
(8.43)
影响速度谱宽的气象因子
谱宽表征着有效照射体内不同大小的多普勒速度偏离其平 均值的程度。谱宽可以用做速度估计质量控制的工具：当谱宽 增加，速度估计的可靠性就减小。对气象目标物而言，影响谱 宽的主要因子有四个：
多普勒天气雷达除此之外，还可利用 降水回波频率与发射频率之间变化的信 息来测定降水粒子的径向速度，并通过 此推断风速分布，垂直气流速度，大气 湍流，降水粒子谱分布，降水中特别是 强对流降水中风场结构特征。
以前，用常规天气雷达进行的天气预报 仅仅使用反射率因子资料。多普勒天气雷达 将提供两种附加的基本资料，径向速度和速 度谱宽，它们将增强对强风暴的探测能力， 也能改进对中尺度和天气尺度系统的预报。
多普勒频率与径向速度的关系
假设多普勒雷达发射脉冲的工作频率为f0，目标与雷达的距
离为r，则雷达波发往目标到返回天线所经过的距离为2r。这 个距离用波长来度量，相当 个波长；用弧度来衡量相当于 个弧度。若所发射的电磁波在天线处的位相为 ，那么电磁波 被散射回到天线时的相位应是
位相的时间变化率
由于目标物的径向运动引起 的雷达回波信号的频率变化，它 就是多普频移或多普勒频率。
多普勒雷达是通过直接测量多普勒 频率来得到径向速度的吗?
4.2 多普勒雷达径向速度探测方法
Pulse-Pair Method 脉冲对方法
取两个连续的脉冲然 后测量接收脉冲的相位， 这种脉冲对位相变化可以 比较容易并且比较准确地 测量
DΦ/dt 实际上就是角 速度 = w = 2πfd
假定当第一个脉冲遇到目标物时，该目标物距雷达的距离为r，则该目标物 产生的回波到达雷达时的位相为：
2、平均多普勒频移及频谱宽度
3、平均多普勒速度和速度谱宽度
注意：脉冲对方法并没有从回波信号中提取频谱或功率谱，从而 不能按以上公式计算和，而是直接对回波信号作简便计算求得。
(8.43)
影响速度谱宽的气象因子
谱宽表征着有效照射体内不同大小的多普勒速度偏离其平 均值的程度。谱宽可以用做速度估计质量控制的工具：当谱宽 增加，速度估计的可靠性就减小。对气象目标物而言，影响谱 宽的主要因子有四个：
多普勒天气雷达原理与应用6-雷达探测算法(1)
中气旋（M）与龙卷涡旋特 征（TVS）产品和算法
步骤1 1D中气旋处理：识别1D 型矢量并测试它们是 否满足切变（THS, TLS）和角动量（THM，TLM） 阈值。 步骤2 2D中气旋处理：1）将1D型矢量组成2D特征， 组成2D特征的最小型矢量的个数为可调参数TPV；2） 确定2D特征的对称性；3）计算2D特征的属性；4） 去掉在高度TFM以上的2D特征。 步骤3 3D中气旋处理：1）将2D特征组成3D特征；2） 丢弃不能组成3D特征的非对称的2D特征；3）不能组 成3D特征的2D特征归类为“非相关切变”；4）包含 少于2个对称2D特征的3D特征归类为“3D相关切变”； 5）包含2个或更多对称2D特征的三维特征归类为“中 气旋”。
表6-7 新的冰雹探测算法HDA的评分 日期 11/02/1992 17/02/1992 25/03/1992 19/04/1992 28/04/1992 28/05/1992 02/06/1992 12/06/1992 09/06/1992 01/09/1989 总体 WT （Jｍ-1s-1） 20 26 63 66 74 97 100 120 126 134 H 16 13 30 16 94 5 3 0 0 40 217 M 1 10 9 12 39 0 3 0 0 20 94 FA 33 11 18 21 32 10 6 5 0 71 207 POD（%） 94 57 77 57 71 100 50 67 70 FAR（%） 67 46 38 59 25 67 67 100 64 49 CSI（%） 32 38 53 31 57 33 25 0 31 42
6
0.084 z
W (Z )
W (Z ) 0
对于 Z L Z Z U 对于 Z Z U
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
“区域”尺寸的选择不宜太大也不宜太小，太大会导致回 波移动向量的分辨率太粗，“区域”太小则包含的数据点 太少，不足以产生稳定的相关系数。发现对于1km ×1km 的分辨率，m取值在3-7之间比较合适。
将平面直角 坐标内的二 维坐标排列 成一维，然 后计算相关 系数：
R [(
k
Z1 (k )
跟踪和外推算法
雷暴或降水的临近预报系统的基础是跟踪和外 推。主要分为两种类型：
• 单体质心跟踪和外推： 将雷暴或降水单元视为三维 单体加以识别、跟踪和外推。典型的例子有WSR-88D 和WDSS中的风暴单体识别与跟踪、以及TITAN等， 下面我们会对TITAN重点进行介绍；
• 区域跟踪和外推：对反射率因子超过某一阈值的二 维区域进行跟踪和外推。典型的例子有TREC等，我 们下面给以重点介绍。
主要思路：对于较短的时间间隔，1）一个雷暴倾向于 沿着一条直线运动；2）风暴的增长和衰减遵循线性趋 势；3）会出现对上述线性行为的随机偏差。
根据演变历史加权线性拟合外推：预报量包括以反射 率因子为权重的雷暴质心、体积、和投影椭圆的参数。
最近，改用多边形（原来为椭圆）对雷暴的水平投影 的面积和形状进行表达。
雷暴分裂：对于在t1时 刻所有的风暴，预报它 们在t2时刻的椭圆投影 的位置、取向和尺寸， 然后根据t2时刻识别的 所有雷暴判断哪些是新 的路径，即没有历史的 雷暴。如果这些没有历 史的雷暴的质心位于某 一个t1时刻雷暴的预报 的在t2时刻的投影椭圆 的区域范围内，则确认 发生了雷暴分裂。
雷暴移动的预报
TITAN
Thunderstorm Identification, Tracking, Analysis, and Nowcasting 雷暴识别、追踪、分析和临近预报 一个基于雷达观测的雷暴临近预报系统
最早版本完成于1986年，1990年代中期得到改进和完善
坐标系
算法采用三维直角坐标系统，由雷达体扫反射率因 子数据得到三维直角坐标系中的反射率因子数据。
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实N Z1 2 ) (
Z1(k) Z2 (k)
k
k
Z2 2 (k) N Z2 2 )]
k
k
其中Z1和Z2是分别是相继两个体扫t1和t2时刻的反射率因子， N是一个“区域”内数据点的数量（N=m2）。
为了避免搜索所有的数据以寻找匹配的“区域”，确 定从初始区域找相关的搜索半径为r=VΔt，其中V为 可能的最大移动速度， Δt是两个体扫之间的时间间隔。
雷暴定义：反射率因子TZ ≥ 35dBZ；体积TV ≥ 50km3
雷暴识别
考虑二维格点，先识 别x方向，再识别y方 向， 再拓展到三维 识别z方向（之上或 者之下）
阴影区：反射率≥ TZ
雷暴特征分析
反射率因子权重 质心（雷暴中心） 体积 雷暴投影到水平 面上的面积大小和形 状（最佳适应形状是 多边形和椭圆）
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
TREC
TREC 是Tracking of Radar Echo with Correlations 的缩 写，既“利用相关跟踪雷达回波”。
1978年由Rinehart和Garvey提出这项技术，用来反演雷 达回波区的气流流场。
该项技术利用交叉相关方法跟踪雷达某一个仰角扫描构 成的锥面上某一个二维回波型。即初始的算法是在由某 一仰角扫描构成的2维圆锥面上进行回波的跟踪。后来， 将该技术应用于直角坐标情况，考虑在某一等高面上的 二维直角坐标系中进行回波跟踪。首先需要将雷达体扫 资料内插到某一等高面（如2.5km）上的直角坐标系中。
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
TITAN分析和追踪预报的个例
2004 年 7 月 29 日 飑 线 的 60 分 钟 预报和实况检验
2004年“7.10”暴雨的30分钟预报 和实况检验
京津地区
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
雷暴追踪
假设T1和T2是相邻的两个雷达体扫资料时间
追踪思路： 1. 宁短不长（考虑
到体扫间隔为56分钟） 2. 特征相似（尺寸 和形状等） 3. 设置雷暴移动速 度上限
将雷暴路径的确定作为一个最优化问题来处理
雷暴合并与分裂的处理
雷暴合并：t1时刻的 雷暴多余t2时刻雷暴 数，或者有雷暴消失， 或者有雷暴合并。如 果t1时刻的2个以上雷 暴质心的预报位置在 t2时刻识别的某个雷 暴范围内，可以判断 雷暴合并。
在直角坐标系中的TREC称为CTREC。
CTREC中的跟踪方法
CTREC使用一定时间间隔内的雷达资料，将反射率因子 场分成若干个大小相当的“区域”，每个“区域”包含 m×m个像素。将这些在上一时刻的“区域”分别与下一 时刻的各个“区域”作空间交叉相关，以找出此刻与上一 个时刻的特定区域相关系数最大的“区域”，从而来确定 整个回波的移动矢量，实现回波的跟踪。
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
将平面直角 坐标内的二 维坐标排列 成一维，然 后计算相关 系数：
R [(
k
Z1 (k )
跟踪和外推算法
雷暴或降水的临近预报系统的基础是跟踪和外 推。主要分为两种类型：
• 单体质心跟踪和外推： 将雷暴或降水单元视为三维 单体加以识别、跟踪和外推。典型的例子有WSR-88D 和WDSS中的风暴单体识别与跟踪、以及TITAN等， 下面我们会对TITAN重点进行介绍；
• 区域跟踪和外推：对反射率因子超过某一阈值的二 维区域进行跟踪和外推。典型的例子有TREC等，我 们下面给以重点介绍。
主要思路：对于较短的时间间隔，1）一个雷暴倾向于 沿着一条直线运动；2）风暴的增长和衰减遵循线性趋 势；3）会出现对上述线性行为的随机偏差。
根据演变历史加权线性拟合外推：预报量包括以反射 率因子为权重的雷暴质心、体积、和投影椭圆的参数。
最近，改用多边形（原来为椭圆）对雷暴的水平投影 的面积和形状进行表达。
雷暴分裂：对于在t1时 刻所有的风暴，预报它 们在t2时刻的椭圆投影 的位置、取向和尺寸， 然后根据t2时刻识别的 所有雷暴判断哪些是新 的路径，即没有历史的 雷暴。如果这些没有历 史的雷暴的质心位于某 一个t1时刻雷暴的预报 的在t2时刻的投影椭圆 的区域范围内，则确认 发生了雷暴分裂。
雷暴移动的预报
TITAN
Thunderstorm Identification, Tracking, Analysis, and Nowcasting 雷暴识别、追踪、分析和临近预报 一个基于雷达观测的雷暴临近预报系统
最早版本完成于1986年，1990年代中期得到改进和完善
坐标系
算法采用三维直角坐标系统，由雷达体扫反射率因 子数据得到三维直角坐标系中的反射率因子数据。
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实N Z1 2 ) (
Z1(k) Z2 (k)
k
k
Z2 2 (k) N Z2 2 )]
k
k
其中Z1和Z2是分别是相继两个体扫t1和t2时刻的反射率因子， N是一个“区域”内数据点的数量（N=m2）。
为了避免搜索所有的数据以寻找匹配的“区域”，确 定从初始区域找相关的搜索半径为r=VΔt，其中V为 可能的最大移动速度， Δt是两个体扫之间的时间间隔。
雷暴定义：反射率因子TZ ≥ 35dBZ；体积TV ≥ 50km3
雷暴识别
考虑二维格点，先识 别x方向，再识别y方 向， 再拓展到三维 识别z方向（之上或 者之下）
阴影区：反射率≥ TZ
雷暴特征分析
反射率因子权重 质心（雷暴中心） 体积 雷暴投影到水平 面上的面积大小和形 状（最佳适应形状是 多边形和椭圆）
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
TREC
TREC 是Tracking of Radar Echo with Correlations 的缩 写，既“利用相关跟踪雷达回波”。
1978年由Rinehart和Garvey提出这项技术，用来反演雷 达回波区的气流流场。
该项技术利用交叉相关方法跟踪雷达某一个仰角扫描构 成的锥面上某一个二维回波型。即初始的算法是在由某 一仰角扫描构成的2维圆锥面上进行回波的跟踪。后来， 将该技术应用于直角坐标情况，考虑在某一等高面上的 二维直角坐标系中进行回波跟踪。首先需要将雷达体扫 资料内插到某一等高面（如2.5km）上的直角坐标系中。
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
TITAN分析和追踪预报的个例
2004 年 7 月 29 日 飑 线 的 60 分 钟 预报和实况检验
2004年“7.10”暴雨的30分钟预报 和实况检验
京津地区
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
安徽2002年8月24日飑线 TITAN 3 0分钟预报与实况
雷暴追踪
假设T1和T2是相邻的两个雷达体扫资料时间
追踪思路： 1. 宁短不长（考虑
到体扫间隔为56分钟） 2. 特征相似（尺寸 和形状等） 3. 设置雷暴移动速 度上限
将雷暴路径的确定作为一个最优化问题来处理
雷暴合并与分裂的处理
雷暴合并：t1时刻的 雷暴多余t2时刻雷暴 数，或者有雷暴消失， 或者有雷暴合并。如 果t1时刻的2个以上雷 暴质心的预报位置在 t2时刻识别的某个雷 暴范围内，可以判断 雷暴合并。
在直角坐标系中的TREC称为CTREC。
CTREC中的跟踪方法
CTREC使用一定时间间隔内的雷达资料，将反射率因子 场分成若干个大小相当的“区域”，每个“区域”包含 m×m个像素。将这些在上一时刻的“区域”分别与下一 时刻的各个“区域”作空间交叉相关，以找出此刻与上一 个时刻的特定区域相关系数最大的“区域”，从而来确定 整个回波的移动矢量，实现回波的跟踪。
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况
安徽2004年7月7日强对流 TITAN 3 0分钟预报与实况