雷达气象学(10)——双偏振雷达及其变量

news2024/9/22 1:14:07

文章目录

10.1 双偏振雷达的优势

回顾本系列的开篇文章《雷达电磁波的散射》就知道,电磁波是横波,我们把电场振动所在平面称为偏振通道。常规天气雷达通过发射水平或垂直偏振通道来探测降水粒子,因此在探测过程中仅对粒子单一方向上的特性比较敏感。

在这里插入图片描述

上图第一张子图中,常规天气雷达所发射的电磁波方向为水平方向,因此雷达探测的降水粒子特性更多地与降水粒子的水平方向上的特性有关。然而,在许多情况下,降水粒子在水平方向上的特性与垂直方向上存在一定的差异。为了能更准确地探测到降水粒子的结构,在常规天气雷达的基础上新增了另一个偏振通道(如上图第二张子图所示),可用于同时探测降水粒子的两个方向上的特性。这样,相比常规天气雷达,双偏振雷达能获得更多关于降水粒子特性的变量。

双偏振雷达有两种观测模式:

  • 交替发射和接收模式:脉间交替发射水平和垂直偏振波,同时接收水平和垂直偏振波。这种模式比较常用。
  • 同时发射和接收模式:同时发射水平和垂直偏振波,同时接收水平和垂直偏振波。这种模式对硬件的要求较高。

10.2 双偏振变量——反映降水粒子特性的变量

在这里插入图片描述

双偏振雷达发射电磁波时,目标取样体积内的水凝物粒子会对电磁波进行后向散射;与此同时,电磁波在探测到目标物之前,也会在降水介质中传播,传播过程中降水粒子会对电磁波产生作用,称为传播效应

双偏振雷达中,与目标区域直接相关的变量有:

  • 反射率因子 Z H Z_H ZH:可以由雷达直接测量
  • 差分反射率 Z D R Z_{DR} ZDR:可以由雷达直接测量
  • 后向散射相位 δ h v \delta_{hv} δhv
  • 协相关系数 ρ h v \rho_{hv} ρhv(CC):可以由雷达直接测量

与传播效应相关的变量有:

  • 差分相位移 ϕ D P \phi_{DP} ϕDP:可以由雷达直接测量
  • 比差分相移 K D P K_{DP} KDP:可间接测量
  • 衰减 A H A_H AH:可间接测量
  • 差分衰减 A D P A_{DP} ADP:可间接测量
  • 退极化比 LDR

这些变量有什么用呢?如下表所示:

双偏振变量用途
差分反射率 Z D R Z_{DR} ZDR判断降水粒子大小
协相关系数 ρ h v \rho_{hv} ρhv(CC)区分降水与非降水,并识别降水相态
比差分相移 K D P K_{DP} KDP判断雨滴数目

下面介绍几个比较重要的双偏振变量。

10.3 差分反射率 Z D R Z_{DR} ZDR

10.3.1 差分反射率的定义

之前在《雷达电磁波的散射》讲过,对于球形粒子,若满足瑞利散射的条件下,反射率因子为:

Z = λ 4 π 5 ∣ K ∣ 2 ∑ i = 1 N σ Z = \frac{\lambda^4}{\pi^5 |K|^2} \sum_{i=1}^{N} \sigma Z=π5K2λ4i=1Nσ

式中, λ \lambda λ 为电磁波波长, K = m 2 − 1 m 2 + 2 K = \frac{m^2 - 1}{m^2 + 2} K=m2+2m21 m m m 为粒子的复折射系数, σ \sigma σ 为降水粒子的后向散射截面。

但是,对于椭球形粒子,反射率因子并不简单:

Z h v = λ 4 π 5 ∣ K ∣ 2 ∑ i = 1 N σ h v Z_{hv} = \frac{\lambda^4}{\pi^5 |K|^2} \sum_{i=1}^{N} \sigma_{hv} Zhv=π5K2λ4i=1Nσhv

这个 σ h v \sigma_{hv} σhv 实际上与水凝物的形状、分布、电磁参数有关,由上式可知:

  • 当水凝物的水平尺寸大于垂直尺寸时,有 σ h > σ v \sigma_h > \sigma_v σh>σv,所以有 Z h > Z v Z_h > Z_v Zh>Zv
  • 当水凝物的水平尺寸小于垂直尺寸时,有 σ h < σ v \sigma_h < \sigma_v σh<σv,所以有 Z h < Z v Z_h < Z_v Zh<Zv

因此,水平通道的反射率因子与垂直通道的反射率因子并不相等。为了反映两者之间的差距,使用水平反射率因子垂直反射率因子之比,即差分反射率 Z d r Z_{dr} Zdr 来表示:

Z d r = Z h Z v Z_{dr} = \frac{Z_h}{Z_v} Zdr=ZvZh

跟反射率因子一样,差分反射率也是使用分贝(dB)作为单位(使用分贝作单位时,下标一律使用大写字母),即:

Z D R = 10 log ⁡ 10 ( Z h Z v ) = Z H − Z V Z_{DR} = 10 \log_{10} \left( \frac{Z_h}{Z_v} \right) = Z_H - Z_V ZDR=10log10(ZvZh)=ZHZV

降水粒子的形状和取向、粒子的介电常数是影响差分反射率的两个关键因素。首先是降水粒子的形状和取向,具体来说:

在这里插入图片描述

  • Z D R > 0 d B Z_{DR} > 0 \mathrm{dB} ZDR>0dB:水平椭球形粒子;
  • Z D R = 0 d B Z_{DR} = 0 \mathrm{dB} ZDR=0dB:球形粒子;
  • Z D R < 0 d B Z_{DR} < 0 \mathrm{dB} ZDR<0dB:垂直椭球形粒子。

而当粒子的介电常数减小时,差分反射率也会随之减小。现在我们来简单介绍不同水凝物(降水粒子)的差分反射率。

10.3.2 雨滴的差分反射率

在这里插入图片描述

雨滴的形状一般呈扁球体(水平椭球形粒子),其形状与大小密切相关。如上图左边所示为雨滴轴长比与雨滴直径的关系,由此可以发现,雨滴越大,形状就越扁, Z D R Z_{DR} ZDR 越大,且大于 0dB。

在这里插入图片描述

在实际的降水观测中,大雨滴和小雨滴一般是并存的,因此 Z D R Z_{DR} ZDR 会受到大雨滴数量和小雨滴数量比例的影响,另外雨滴的数量和分布也会影响 Z D R Z_{DR} ZDR。通常而言,强降雨伴随着数量较多的大雨滴,而对于弱降水来说,小雨滴占多数。

如果我们把影响 Z D R Z_{DR} ZDR 的因素综合起来,便引出了雨滴谱的概念:单位体积中雨滴的数量(或数密度)随雨滴直径的变化关系称为雨滴谱,它是雨滴生成、下落、增长、破碎、蒸发等过程的综合结果。正是不同的雨滴谱分布导致了不同的 Z D R Z_{DR} ZDR

下面我们具体说说为什么双偏振雷达能识别不同的雨滴谱。不妨考虑两个不同的雨滴谱:

雨滴谱12
雨滴直径1mm3mm
雨滴数量7291
雨强5.588mm/h0.508mm/h

由反射率因子的计算公式 Z = ∑ i = 1 N D i 6 Z=\sum_{i=1}^{N} D_i^6 Z=i=1NDi6 可知,两者的反射率均为 29dBZ。相同的反射率意味着什么?意味着对于传统天气雷达(即单偏振的天气雷达),这两种不同的雨滴谱探测出来的结果是一样的,我们无法分辨它们的区别,尽管两者的雨强相差甚远,而双偏振雷达却能够区分两者的差别,这就是传统天气雷达所不能触及的优势了。

10.3.3 冰、雪和冰雹的差分反射率

对于冰,在下落过程中,因受到空气阻力作用,其等效形状依然是扁的,因此 Z D R > 0 d B Z_{DR} > 0 \mathrm{dB} ZDR>0dB。不过由于冰晶本身的尺寸较小,差分反射率也是比较小的。

对于干雪,在下落过程中,它的等效形状依然是扁的,但是其结构松散,可看做冰和空气所组成的粒子,因此介电常数较小, Z D R Z_{DR} ZDR 也较小,一般接近于 0dB。

对于湿雪,可以将其看做水、冰和空气三者的混合体。融化一部分的湿雪,其粒子形状接近于干雪,湿雪融化越多,介电常数就越大,且接近于水得介电常数,因此反射率因子和 Z D R Z_{DR} ZDR 都会增大,但是 Z D R > 0 d B Z_{DR} > 0 \mathrm{dB} ZDR>0dB

对于冰雹,由于在下落过程中会不断地翻滚,其形状也会慢慢靠近球形,其 Z D R Z_{DR} ZDR 会接近于 0dB。

下面给出几个例子。第一个例子是南京信息工程大学 C 波段双偏振雷达的低仰角观测:

在这里插入图片描述

这次降水的类型为层云降水,可以看到左图的反射率因子普遍较弱,主要对应的降水粒子为雪。从右边的 Z D R Z_{DR} ZDR 图可知,其值接近于 0dB,而下面有一块区域的 Z D R Z_{DR} ZDR 相对较大,可初步说明这片区域的雪处于融化状态,因为湿雪会使 Z D R Z_{DR} ZDR 增大。

第二个例子是南京信息工程大学 C 波段双偏振雷达的高仰角观测:

在这里插入图片描述

大气高空中主要为冰晶粒子,由右边的 Z D R Z_{DR} ZDR 图可知,这些冰晶粒子的形状各异,对应不同的 Z D R Z_{DR} ZDR,且最大有 3dB 左右。随着下落过程的持续,聚合作用使得冰晶形成雪,粒子形状逐渐变为球形, Z D R Z_{DR} ZDR 将变小。

第三个例子是广州 S 波段双偏振雷达的观测:

在这里插入图片描述

左上图为 PPI 反射率因子图,可看到最大回波强度有 70dBZ,而且还能看到三体散射回波 TBSS 位于回波核心区的后侧,其沿着雷达径向方向延伸,说明存在着冰雹天气。而在右上的 PPI- Z D R Z_{DR} ZDR 图可以看到,回波核心区即冰雹区的 Z D R Z_{DR} ZDR 接近 0dB,而冰雹区周围的降水区则大于 0dB。

同样的特征发生在 RHI 上,右下图为 RHI 反射率因子图,冰雹区延伸到 8km 高度;而在右下的 RHI- Z D R Z_{DR} ZDR 图中,对应区域的 Z D R Z_{DR} ZDR 接近 0dB。

上述这几个例子提醒我们,不是有了双偏振雷达就高枕无忧了,很多时候也要参考单通道雷达的探测,要两只手一起抓才能准确判断降水类型。

10.3.4 层云降水的差分反射率

在这里插入图片描述

在这个例子中,两张图均为 PPI,左边的反射率因子图和右边的 Z D R Z_{DR} ZDR 图均出现了零度层亮环。关于零度层亮带的详细介绍已于本系列的《反射率因子图分析(气象回波篇)》讲过,现在来分析为什么在 Z D R Z_{DR} ZDR 图上会出现这样的零度层亮带:

  • 降水粒子在 0˚C 层以上:较大的降水粒子大多为冰晶和雪花,所以在右图的亮环以外 Z D R Z_{DR} ZDR 较弱;
  • 降水粒子在 0℃ 层附近:此时降水粒子刚发生融化,其表面出现一层水膜。同时,冰晶和雪花在下降过程中,有强烈的碰并聚合作用,导致粒子逐渐变成球形,介电常数增大,所以 Z D R Z_{DR} ZDR 也随之增大,在图中呈现一个亮环;
  • 降水粒子在 0˚C 层往下:当冰晶和雪花继续下降至完全融化为水滴时, Z D R Z_{DR} ZDR 也随之减小。

10.3.5 对流降水的差分反射率

在这里插入图片描述

这是 2016 年 5 月 15 日广州的一次飑线过程,两张图均为 RHI。左图为反射率因子图,可发现对流高度达到 8km,强回波区突破零度层,表明可能有霰粒子或冰雹粒子。在右边 Z D R Z_{DR} ZDR 图所对应区域也有一个强 Z D R Z_{DR} ZDR 区,称为 ZDR 柱,其 Z D R Z_{DR} ZDR 值随高度增加而减小。

10.4 相关系数 ρ h v \rho_{hv} ρhv

相关系数 ρ h v \rho_{hv} ρhv 在雷达应用中也可写为 CC,即 Correlation Coefficient 的缩写,无单位,用于描述同时接收到的前后脉冲变化情况的 Z H Z_H ZH Z V Z_V ZV 之间相关程度的物理量,一般可用于识别降水和非降水回波。下图为常见水凝物的 CC:

在这里插入图片描述

  • 小雨滴的形状接近于球形,因此 CC 接近于 1;而随着雨滴的增大,其形状变为扁平,CC 也会有所降低。
  • 对于冰雹和霰,当粒子较小时,CC 一般大于 0.95;当粒子开始融化时,CC 开始降低。
  • 对于干雪,CC 较大且接近 1;对于湿雪,为融化中的雪,CC 降低,甚至可以低于 0.9。
  • 对于杂波和地物杂波,CC 明显小于 0.8,这也是为什么 CC 能识别降水和非降水回波的缘故。

10.5 比差分相移 K D P K_{DP} KDP

我们知道,电磁波在真空中以光速传播,而在电介质中的传播速度会降低,因此对于同样的传播距离,真空中的电磁波的波数要比电介质中的电磁波少,或者说两者存在着相位差。同样的道理,对于非球形粒子,水平偏振波和垂直偏振波的传播速度是不同的,所以随着传播距离的增加,两者的波数也会存在差异,或者说两者存在着相位差,这个相位差就称为差分相移 ϕ D P \phi_{DP} ϕDP,如下图所示。

在这里插入图片描述

ϕ D P \phi_{DP} ϕDP 对较大尺寸水滴非常敏感,而且粒子数量越多,它的数值变化越快,但它是累积效应的参数,即随距离的增加而线性增加,其数值在不同位置、不同时刻不能直接对比使用。于是引入了一个新的参数:比差分相移。

在单位距离上,由电介质所导致的水平通道和垂直通道的相位差,就称为比差分相移 K D P K_{DP} KDP,单位为度/km,其物理意义是单位距离上降水介质对 ϕ D P \phi_{DP} ϕDP 的贡献,其计算公式如下:

K D P = 1 2 ∂ ϕ D P ∂ r K_{DP} = \frac{1}{2} \frac{\partial \phi_{DP}}{\partial r} KDP=21rϕDP

下图为 S 波段各种水凝物粒子的 K D P K_{DP} KDP 值:

在这里插入图片描述

  • 降雨随着粒子的直径和浓度增大, K D P K_{DP} KDP 逐渐增大。
  • 对于同样的雨滴谱分布,电磁波波长越短, K D P K_{DP} KDP 越大。
  • 对于冰雹、霰、雪, K D P K_{DP} KDP 较小且接近于 0;但是当它们开始融化时,由于粒子外部包裹着一层水膜, K D P K_{DP} KDP 显著增大。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2051541.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Chrome快捷键提高效率

浏览效率提高快捷建 快速切换标签页 Ctrl 数字&#xff08;1或者2&#xff09;&#xff0c;标签页数字从左到右为顺序&#xff0c;1开始。快速切换标签页。 Ctrl1 到 Ctrl8 切换到标签栏中指定位置编号所对应的标签页 Ctrl9切换到最后一个标签页 CtrlTab 或 CtrlPgDown 切…

Elasticsearch中磁盘水位线的深度解析

❃博主首页 &#xff1a; 「码到三十五」 &#xff0c;同名公众号 :「码到三十五」&#xff0c;wx号 : 「liwu0213」 ☠博主专栏 &#xff1a; <mysql高手> <elasticsearch高手> <源码解读> <java核心> <面试攻关> ♝博主的话 &#xff1a…

Clobbering DOM attributes to bypass HTML filters

目录 寻找注入点 代码分析 payload构造 注入结果 寻找注入点 DOM破坏肯定是出现在js文件中&#xff0c;我们首先来看源码 /resources/labheader/js/labHeader.js这个源码没什么问题我们重点关注在下面两个源码上 /resources/js/loadCommentsWithHtmlJanitor.js这个源码中重…

从关键新闻和最新技术看AI行业发展(第二十九期2024.7.29-8.11) |【WeThinkIn老实人报】

写在前面 【WeThinkIn老实人报】旨在整理&挖掘AI行业的关键新闻和最新技术&#xff0c;同时Rocky会对这些关键信息进行解读&#xff0c;力求让读者们能从容跟随AI科技潮流。也欢迎大家提出宝贵的优化建议&#xff0c;一起交流学习&#x1f4aa; 欢迎大家关注Rocky的公众号&…

docker 安装mino服务,启动报错: Fatal glibc error: CPU does not support x86-64-v2

背景 docker 安装mino服务&#xff0c;启动报错&#xff1a; Fatal glibc error: CPU does not support x86-64-v2 原因 Docker 镜像中的 glibc 版本要求 CPU 支持 x86-64-v2 指令集&#xff0c;而你的硬件不支持。 解决办法 降低minio对应的镜像版本 经过验证&#xff1a;qu…

sanic + webSocket:股票实时行情推送服务实现

&#x1f49d;&#x1f49d;&#x1f49d;欢迎莅临我的博客&#xff0c;很高兴能够在这里和您见面&#xff01;希望您在这里可以感受到一份轻松愉快的氛围&#xff0c;不仅可以获得有趣的内容和知识&#xff0c;也可以畅所欲言、分享您的想法和见解。 推荐&#xff1a;「storm…

Linux中的内核编程

Linux内核是操作系统的核心组件&#xff0c;负责管理系统的资源、提供硬件抽象和执行系统调用。内核编程是一项涉及操作系统核心的高级任务&#xff0c;它允许开发人员直接与系统内核进行交互&#xff0c;实现更高效、更特定的功能。本文将深入探讨Linux中的内核编程&#xff0…

2021年上半年网络工程师考试上午真题

2021年上半年网络工程师考试上午真题 网络工程师历年真题含答案与解析 第 1 题 以下关于RISC和CISC计算机的叙述中&#xff0c;正确的是&#xff08; &#xff09;。 (A) RISC不采用流水线技术&#xff0c;CISC采用流水线技术(B) RISC使用复杂的指令&#xff0c;CISC使用简…

事件驱动架构的事件版本管理

有一种办法&#xff1a;发送会议邀请给所有团队&#xff0c;经过101次会议后&#xff0c;发布维护横幅&#xff0c;所有人同时点击发布按钮。或... 可用适配器&#xff0c;但微调。没错&#xff01;就像软件开发中90%问题一样&#xff0c;有种模式帮助你找到聪明解决方案。 1…

C Primer Plus(中文版)第13章编程练习,仅供参考

第十三章编程练习 对于文件的操作是程序开发过程中必不可少的。首先&#xff0c;来看一下第一题&#xff0c;对13.1程序进行修改&#xff0c;输入文件名&#xff0c;而不是命令行参数。完整程序代码以及运行结果如下&#xff1a; #include<stdio.h> #include<stdlib…

【数据结构篇】~单链表(附源码)

【数据结构篇】~链表 链表前言链表的实现1.头文件2.源文件 链表前言 链表是一种物理存储结构上非连续、非顺序的存储结构&#xff0c;数据元素的逻辑顺序是通过链表中的指针链接次序实现的。 1、链式机构在逻辑上是连续的&#xff0c;在物理结构上不一定连续​ 2、结点一般是从…

Java二十三种设计模式-命令模式(18/23)

命令模式&#xff1a;将请求封装为对象的策略 概要 本文全面探讨了命令模式&#xff0c;从基础概念到实现细节&#xff0c;再到使用场景、优缺点分析&#xff0c;以及与其他设计模式的比较&#xff0c;并提供了最佳实践和替代方案&#xff0c;旨在帮助读者深入理解命令模式并…

【xr-frame】微信小程序xr-frame典型案例

微信小程序xr-frame典型案例 在之前的工作中&#xff0c;我大量使用XR-Frame框架进行AR开发&#xff0c;并积累了一些案例和业务代码。其中包括2D图像识别、手部动作识别、Gltf模型加载、动态模型加载、模型动画等内容。小程序部分使用TypeScript编写&#xff0c;而XR-Frame组…

利用puppeteer将html网页生成图片

1.什么是puppeteer&#xff1f; Puppeteer是一个Node库&#xff0c;它提供了一个高级API来通过DevTools协议控制Chromium或Chrome。 可以使用Puppeteer来自动化完成浏览器的操作&#xff0c;官方给出的一些使用场景如下&#xff1a; 生成页面PDF抓取 SPA&#xff08;单页应用…

3.Windows Login Unlocker-忘记电脑密码也可以解决

想要解锁Windows系统的开机密码&#xff0c;但官网的传统方法只适合Windows本地账户&#xff0c;对微软账户或PIN码()束手无策&#xff1f;别担心&#xff0c;小编之前推荐过的「Windows Login Unlocker」软件能为您排忧解难。这款出色的工具不仅能够轻松绕过各种Windows密码&a…

C语言-写一个用矩形法求定积分的通用函数,分别求积分区间为[0,1]sinx,cosx,e的x方的定积分

一、题目要求&#xff1a; 二、思路 ①数学方面:矩形法求定积分的公式 将积分图形划分成为指定数量的矩形&#xff0c;求取各个矩形的面积&#xff0c;然后最终进行累加得到结果 1.积分区间: [num1, num2] 2.分割数量:count 每个矩形的边长:dx(num2-num1)/count 3.被积分…

智慧草莓基地管理系统--论文pf

TOC springboot359智慧草莓基地管理系统--论文pf 绪论 1.1 研究背景 当前社会各行业领域竞争压力非常大&#xff0c;随着当前时代的信息化&#xff0c;科学化发展&#xff0c;让社会各行业领域都争相使用新的信息技术&#xff0c;对行业内的各种相关数据进行科学化&#xf…

泰安鲁菜根普照店重装开业:传承与创新并举 品味舌尖上的泰安

8月17日&#xff0c;泰安鲁菜根普照店重装开业。店内装修以泰山文化为主题&#xff0c;营造出一种浓郁的地方文化特色氛围。这家以泰山文化为底蕴&#xff0c;以大汶河传统民俗食材为基础的餐饮企业&#xff0c;一直致力于发掘和传承泰山周边及汶河两岸的传统特色美食&#xff…

【漫谈C语言和嵌入式009】探索LVDT:线性可变差动变压器的工作原理与应用

引言 在现代工业和工程领域中&#xff0c;精确的位移测量是许多系统正常运行的关键。线性可变差动变压器&#xff08;LVDT, Linear Variable Differential Transformer&#xff09;是一种广泛使用的位移传感器&#xff0c;因其高精度、可靠性和耐用性在各种应用中得到了普遍认可…

iOS(OC)学习第1天-怎么设置UI

xCode 版本为 新建工程 xCode->iOS->App->输入工程名称 项目结构 storyboard文件&#xff1a;可以通过拖动方式添加UI的UI布局文件&#xff0c;理解为Android的XML布局文件 名字说明Main.storyboard首页LaunchScreen.storyboard引导页 添加布局 引导页 设置组件的属…