一种太赫兹与毫米波云雷达的数据联合反演方法技术

本发明专利技术公开了一种太赫兹与毫米波云雷达联合反演云物理参数的方法，建立太赫兹雷达反射率因子、毫米波雷达反射率因子与云衰减的关系；建立太赫兹频段和毫米波频段雷达反射率因子与云物理参量的关系，基于贝叶斯理论建立概率公式，迭代运算云物理参量，将物理模型值与观测值的误差作为迭代收敛的判据，得到最优解；根据反演得到的云物理参数，迭代计算云粒子的后向散射截面和衰减截面，更新雷达反射率因子和云衰减的函数关系，得到新的反演值，将云衰减系数在迭代运算中的变化作为收敛条件，得到云物理参量的双重最优解，为气候预测、天气预报、大气科学研究提供数据支撑，推动太赫兹云雷达和毫米波云雷达在气象观测系统中的有效应用。有效应用。有效应用。

【技术实现步骤摘要】
一种太赫兹与毫米波云雷达的数据联合反演方法


[0001]本专利技术涉及雷达领域，具体涉及一种太赫兹与毫米波云雷达的数据联合反演方法。

技术介绍

[0002]云是地球气候系统的重要组成部分，对地球-大气系统的辐射收支平衡和水汽循环有着重要调节作用，云粒子大小、形状、数浓度、冰水含量、粒子谱等微观特征的精度直接影响其散射辐射特性描述的准确性。联合国政府间气候变化专门委员会(IPCC)第三次报告指出，气候模式中大量的不确定性来源于云物理信息处理的不足，以及缺乏实际观测资料来进行云参数化方案验证，第五次报告指出，在气候模式中云量化仍然是一项挑战，模式对云辐射效应的估算结果存在很大的不确定性。因此，开展云物理信息的观测和反演研究，对于云物理特征分析、云辐射效应、天气预报、气候预测模型等研究具有重要意义。
[0003]毫米波云雷达(Ka和W波段)是目前的主要云遥感手段之一，可以实现层云、积云等非降水云的连续监测，但对粒径较小的云粒子穿透性过强，回波强度较弱，特别是对含水量少的薄冰云探测效果不佳，因而需要波长更接近微小粒子尺度的探测系统。THz波是指频率在0.1-10THz范围内的电磁波，与毫米波云雷达相比，该波长更接近冰云粒子尺寸，在理论上具有更高的探测能力和灵敏度，且适用于星载平台以及高纬度地区的高精度云遥感研究。利用太赫兹云雷达和毫米波云雷达观测数据开展联合反演研究，有助于获取高精度的云物理特征参数。
[0004]毫米波云雷达的反演研究多为单频反演研究，一般认为雷达反射率因子与粒子有效半径、云水含量之间存在经验关系，为提高反演精度，在经验关系中引入大气温度和观测站点的地理信息。2006年发射的CloudSat是第一个搭载94GHz云廓线雷达的卫星，其业务化反演算法基于最优估计理论开发，适用于全球云微物理参量的反演。
[0005]高频电磁波信号在传播路径上会受到衰减，使得实测雷达回波强度低于真实强度，造成了云参数的反演误差。在使用太赫兹与毫米波雷达的数据进行联合处理时，要考虑衰减订正，现有衰减订正方法多基于毫米波频段经验公式进行。
[0006]综上所述，现有云物理参数的反演和其中涉及的衰减订正方法多适用于毫米波频段，太赫兹与毫米波的联合反演局限在理论分析和探讨方面，有待深入研究。

技术实现思路

[0007]本专利技术的目的在于提供一种太赫兹与毫米波云雷达的数据联合反演方法，利用常见云遥感的毫米波频点与太赫兹频点的组合，研究利用双频雷达反射率因子差来反演冰云粒子半径和冰水含量的方法，对反演精度进行模拟计算和比对分析。毫米波(特别是W波段)与太赫兹波组合对微小冰晶粒子的探测最为敏感，反演结果精度最高。
[0008]为了达到上述目的，本专利技术采用的技术方案如下：
[0009]一种太赫兹与毫米波云雷达的数据联合反演方法，通过太赫兹云雷达和毫米波云
雷达进行联合反演，得到云参数的最优解。
[0010]可选地，所述方法具体包括：
[0011]S1、计算得到太赫兹雷达反射率因子与云衰减系数的关系，以及毫米波雷达反射率因子与云衰减系数的关系，得到α
THz
、β
THz
和α
W
、β
W
的初始值；
[0012]S2、分别对所述太赫兹雷达实测反射率因子和所述毫米波雷达实测反射率因子进行衰减订正，得到衰减订正后的每个距离库的雷达反射率因子和
[0013]S3、联合衰减订正后的太赫兹与毫米波雷达反射率因子，分别建立太赫兹频段和毫米波频段雷达回波与云物理参量的关系，由于雷达反射率因子无法得到唯一的云物理参量，待计算量多于已知量，基于贝叶斯理论建立概率公式，将云物理特征统计结果作为待反演云参数的初值，计算得到最可能的云物理参量。在数学计算上，表现为目标函数的最小化求解。
[0014]可选地，所述步骤S1具体包括：
[0015]S101、计算太赫兹频段的云衰减系数和毫米波频段的云衰减系数计算太赫兹雷达实际反射率因子和毫米波雷达实际反射率因子
[0016]S102、根据所述太赫兹频段的云衰减系数和所述太赫兹雷达实测反射率因子计算得到太赫兹频段的回归系数α
THz
、β
THz
，根据所述毫米波频段的云衰减系数和所述毫米波雷达实测反射率因子计算得到毫米波频段的回归系数α
W
、β
W
，得到α
THz
、β
THz
和α
W
、β
W
的初始值。
[0017]可选地，所述太赫兹频段的云衰减系数和所述毫米波频段的云衰减系数均通过如下公式进行计算：
[0018][0019][0020]其中，k
c
是云衰减系数；Q
ext
(D)是直径为D的单个云滴的衰减截面，等于云滴的吸收截面和散射截面之和；N(D)是粒子谱分布函数，根据云物理特征的统计分析结果，采用对数正态分布；N0是粒子数密度；D为粒子直径；D
g
为几何平均直径；w是分布宽度参数，为无量纲变量；ln代表自然对数。
[0021]可选地，所述太赫兹雷达实测反射率因子和所述毫米波雷达实测反射率因子均通过如下公式进行计算：
[0022][0023]其中，Z
e
是雷达实测反射率因子，Z
e
是雷达实测反射率因子，Q
bk
表示粒子直径D的球形粒子的后向散射截面；K由粒子的负折射指数计算得到。
k
c
(t)，当算法收敛时，差值逐渐减小，设置差值小于等于0.001dB/km为收敛判据，满足判据时，停止运算，当前得到的k
c
和x为最终解，至此完成联合计算过程。
[0041]与现有技术相比，本专利技术至少具有以下优点之一：
[0042]利用常见云遥感的毫米波频点与太赫兹频点的组合，研究利用双频雷达反射率因子差来反演冰云粒子半径和冰水含量的方法，对反演精度进行模拟计算和比对分析。毫米波(特别是W波段)与太赫兹波组合对微小冰晶粒子的探测最为敏感，反演结果精度最高。
附图说明
[0043]图1为本专利技术一实施例中太赫兹云雷达和毫米波云雷达联合反演云参数的计算流程图。
具体实施方式
[0044]以下结合附图1具体实施方式对本专利技术作进一步详细说明。根据下面说明，本专利技术的优点和特征将更清楚。需要说明的是，附图采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本专利技术实施方式的目的。为了使本专利技术的目的、特征和优点能够更加明显易懂，请参阅附图。须知，本说明书所附图式所绘示的结构、比例、大小等，均仅用以配合说明书所揭示的内容，以供熟悉此技术的人士了解与阅读，并非用以限定本专利技术实施的限定条件，故不具技术上的实质意义，任何结构的修饰、比例关系的改变或大小的调整，在不影响本专利技术所能产生的功效及所能达成的目的下，均应仍落在本专利技术所揭示的
技术实现思路
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种太赫兹与毫米波云雷达的数据联合反演方法，其特征在于，利用太赫兹云雷达和毫米波云雷达的观测数据进行联合反演，基于双重优化得到云物理参数的最优解。2.如权利要求2所述的太赫兹与毫米波云雷达的数据联合反演方法，其特征在于，所述方法具体包括：S1、计算得到太赫兹雷达反射率因子与云衰减系数的关系，以及毫米波雷达反射率因子与云衰减系数的关系，得到α
THz
、β
THz
和α
W
、β
W
的初始值；S2、分别对所述太赫兹雷达实测反射率因子和所述毫米波雷达实测反射率因子进行衰减订正，得到衰减订正后的每个距离库的雷达反射率因子和S3、联合衰减订正后的太赫兹与毫米波雷达反射率因子，基于贝叶斯理论建立概率公式，迭代计算云物理参量，将雷达反射率因子的物理模型值与观测值的误差作为迭代收敛的判据，直到得到待反演云参量的最优解；S4、重复步骤1-3，根据每次反演得到的云物理参数x(t)，迭代计算后向散射截面和衰减截面，结合粒子谱分布，得到新的雷达反射率因子Z
e
(t+1)和云衰减系数k
c
(t+1)的函数关系，其中，t表示衰减修正次数，再将新的实际雷达反射率因子Z
e
(t+1)代入反演计算，得到x(t+1)，将云衰减系数在反演过程中的变化作为收敛条件，得到双重优化的云物理参量，至此完成联合反演。3.如权利要求2所述的太赫兹与毫米波云雷达的数据联合反演方法，其特征在于，与单频雷达反射率因子作为输入相比，根据太赫兹频段和毫米波频段雷达回波与云物理参量的前向模型关系进行最优估计反演，两个不同频段的雷达反射率因子能够更好的约束待反演的云参数。具体计算公式如下：...

【专利技术属性】
技术研发人员：丁霞，王海涛，付朝伟，倪亮，
申请(专利权)人：上海无线电设备研究所，
类型：发明
国别省市：