一种多波段双偏振雷达的冰晶凇附增长估算方法技术

本发明专利技术涉及一种多波段双偏振雷达的冰晶凇附增长估算方法，利用冰晶增长的三维模型，生成不同凇附增长的冰晶参数，利用T矩阵计算所建立冰晶模型的后向散射特性，进一步利用雷达前向模拟系统仿真天顶指向的多波段偏振雷达观测，之后利用神经网络模型映射多波段偏振雷达观测和凇附质量分数的关系，实现基于多波段偏振雷达观测的凇附增长的估算。本发明专利技术不需要其他辅助参数，直接由雷达观测数据来估算，且由于估算模型建立过程中，考虑了大气温度和湿度环境因素、冰晶的微物理参数、冰晶凇附增长粒子的后向散射特性和雷达辐射传输的物理过程，并且估算过程中使用了多波段双偏振雷达观测，因此获得的冰晶凇附质量分数精度较高。因此获得的冰晶凇附质量分数精度较高。因此获得的冰晶凇附质量分数精度较高。

【技术实现步骤摘要】
一种多波段双偏振雷达的冰晶凇附增长估算方法


[0001]本专利技术涉及气象雷达探测
，尤其涉及一种多波段双偏振雷达的冰晶凇附增长估算方法。

技术介绍

[0002]冰晶生长主要有凝华、丛集和凇附三个过程，其中凇附过程主要指冰晶与云中的过冷水滴碰撞，并在下降过程时冻结；凇附是冰晶增长的重要过程之一，它对降水量和强度有重要影响，已有研究表明，超过40%的累积降雪可能是由于冰晶凇附过程造成的，并且凇附增长程度和降雪强度直接相关，这意味着凇附是冰晶、雪花生长和最终形成降水的关键过程。此外，凇附过程可影响气溶胶的间接效应，对全球大气的云
‑
气溶胶相互作用有重要影响。
[0003]在目前对地球大气的认识中，降水是一个关键的不确定因素。特别是在数值天气预报和气候模型中，降水的表示是目前最薄弱的环节之一，部分原因是在模拟控制降水颗粒形成和生长的潜在微物理过程方面存在缺陷。比如，冰相微物理增长过程、速率有较大不确定性，这对数值天气预报和气候模型中的微物理方案提出了重大挑战。目前，数值预报模型在凇附参数化存在较大困难，主要原因之一是缺乏对凇附过程的观测和定量估算，因此，如何对凇附增长过程进行估算，是现阶段需要考虑的问题。
[0004]需要说明的是，在上述
技术介绍
部分公开的信息只用于加强对本公开的背景的理解，因此可以包括不构成对本领域普通技术人员已知的现有技术的信息。

技术实现思路

[0005]本专利技术的目的在于克服现有技术的缺点，提供了一种多波段双偏振雷达的冰晶凇附增长估算方法，解决了现有技术不能对凇附增长过程进行估算的问题。
[0006]本专利技术的目的通过以下技术方案来实现：一种多波段双偏振雷达的冰晶凇附增长估算方法，所述估算方法包括：S1、获取待观测区域的大气温度和湿度空气球数据，根据大气温度、湿度和位势高度数据，结合三维云模型判断冰晶粒子存在的位置信息，并根据获取ECMWF数据得到三维模型的初始猜测值；S2、设置冰晶粒子的参数，生成不同冰晶凇附质量分数的冰晶粒子宏微观参数数据集；S3、根据雷达的频率、谱宽、极化参数和不同冰晶凇附质量分数的冰晶粒子宏微观数据，利用T矩阵方法计算不同凇附质量分数冰晶粒子的后向散射特征参数；S4、将观测得到的几何信息、雷达仪器参数、ECMWF数据和后向散射特征参数输入到雷达前向模拟系统，模拟指定雷达性能参数条件下的雷达观测数据集，并获取不同冰晶凇附质量分数和雷达观测值的时空匹配的数据集；S5、根据不同冰晶凇附质量分数和雷达模拟观测值构建神经网络估算模型，并对
不同冰晶粒子的凇附质量分数进行估算。
[0007]所述根据获取ECMWF数据得到三维模型的初始猜测值包括：获取ECMWF云量和比云冰含水量，根据探空气球的位置和观测时间对ECMWF数据进行双线性插值，将插值后的云量、比云液态含水量和比云冰含水量作为三维云模型的初始猜测值。
[0008]所述设置冰晶粒子的参数，生成不同冰晶凇附质量分数的冰晶粒子宏微观参数数据集包括：利用步骤S1得到的冰晶粒子位置和冰水含量作为初始值，设置冰晶粒子的形状和有效粒子半径参数，生成基准冰晶粒子的宏微观参数数据集；调整基准冰晶粒子的凇附质量分数，生成不同冰晶凇附质量分数的冰晶粒子宏微观参数数据集。
[0009]所述利用T矩阵方法计算不同凇附质量分数冰晶粒子的后向散射特征参数包括：将电磁波的入射和散射场描述为矢量球面波函数，并通过T矩阵将函数的展开项进行关联；将冰晶粒子的尺度参数、形状、折射率、粒子谱分布和雷达频率输入到T矩阵中计算得到冰晶粒子的后向散射特征参数。
[0010]所述时空匹配的数据集包括：冰晶凇附质量分数，模拟的多波段双线偏振雷达水平通道反射率因子Zhh，垂直通道反射率因子Zvv，水平和垂直正交偏振波Zvh ，差分反射率因子ZDR，差分传播相移率KDP，相关系数CC，反射率因子的功率谱Spectra_Zhh、Spectra_Zvv、Spectra_Zvh，观测角度，频率和发射功率。
[0011]所述根据不同冰晶凇附质量分数和雷达模拟观测值构建神经网络估算模型，并对不同冰晶粒子的凇附质量分数进行估算包括：将雷达模拟观测值作为神经网络的输入，将不同冰晶凇附质量分数作为神经网络的输出，构建包含两个隐藏层的前向三层神经网络；对神经网络进行训练和验证评估，确定隐藏层节点和最佳的权重函数，得到构建好的神经网络估算模型；将真实雷达观测的数据集输入到构建好的神经网络估算模型中进行不同冰晶粒子的凇附质量分数的估算。
[0012]所述凇附质量分数计算公式包括：，其中，ICW
ur
是冰晶粒子发生凇附之前的质量，ICW是冰晶粒子发生凇附之后的质量。
[0013]所述宏微观参数包括：冰晶粒子数浓度、冰晶形状特征、冰晶粒子尺度谱、冰晶有效粒子半径、冰晶对应的云量和云厚。
[0014]所述后向散射特征参数包括散射效率、散射相函数、单次散射反照率、散射效率因子和吸收效率因子。
[0015]本专利技术具有以下优点：一种多波段双偏振雷达的冰晶凇附增长估算方法，基于Ka，X和W波段的双偏振雷达观测数据开展，能够估算凇附质量分数反演，进而实现凇附增长的估算；不需要其他辅助参数，直接由雷达观测数据来估算，且由于估算模型建立过程中，考虑了大气温度和湿度环境因素、冰晶的微物理参数、冰晶凇附增长粒子的后向散射特性和雷达辐射传输的物理过程，并且估算过程中使用了多波段双偏振雷达观测，因此获得的冰
晶凇附质量分数精度较高。
附图说明
[0016]图1 为本专利技术的流程示意图；图2 为本专利技术的不同凇附质量分数情况下冰晶粒子凇附增长后向散射计算方法的流程图；图3 为本专利技术不同环境参数、冰晶参数和仪器参数情况下多波段双线偏振雷达模拟观测计算方法的流程图。
具体实施方式
[0017]为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此，以下结合附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的保护范围，而是仅仅表示本申请的选定实施例。基于本申请的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。下面结合附图对本专利技术做进一步的描述。
[0018]如图1所示，本专利技术涉及一种多波段双偏振雷达的冰晶凇附增长估算方法，利用冰晶增长的三维模型，生成不同凇附增长的冰晶参数，利用T矩阵计算所建立冰晶模型的后向散射特性，进一步利用雷达前向模拟系统仿真天顶指向的多波段偏振雷达观测，之后利用神经网络模型映射多波段偏振雷达观测和凇附质量分数的关系，实现基于多波段偏振雷达观测的凇附增长的估算；具体包括以下内容：步骤1、获取Ka，X和W波段双偏振雷达观测区域的大气温度和湿度探空气球数据，根据大气温度、湿度和位势本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种多波段双偏振雷达的冰晶凇附增长估算方法，其特征在于：所述估算方法包括：S1、获取待观测区域的大气温度和湿度空气球数据，根据大气温度、湿度和位势高度数据，结合三维云模型判断冰晶粒子存在的位置信息，并根据获取ECMWF数据得到三维模型的初始猜测值；S2、设置冰晶粒子的参数，生成不同冰晶凇附质量分数的冰晶粒子宏微观参数数据集；S3、根据雷达的频率、谱宽、极化参数和不同冰晶凇附质量分数的冰晶粒子宏微观数据，利用T矩阵方法计算不同凇附质量分数冰晶粒子的后向散射特征参数；S4、将观测得到的几何信息、雷达仪器参数、ECMWF数据和后向散射特征参数输入到雷达前向模拟系统，模拟指定雷达性能参数条件下的雷达观测数据集，并获取不同冰晶凇附质量分数和雷达观测值的时空匹配的数据集；S5、根据不同冰晶凇附质量分数和雷达模拟观测值构建神经网络估算模型，并对不同冰晶粒子的凇附质量分数进行估算。2.根据权利要求1所述的一种多波段双偏振雷达的冰晶凇附增长估算方法，其特征在于：所述根据获取ECMWF数据得到三维模型的初始猜测值包括：获取ECMWF云量和比云冰含水量，根据探空气球的位置和观测时间对ECMWF数据进行双线性插值，将插值后的云量、比云液态含水量和比云冰含水量作为三维云模型的初始猜测值。3.根据权利要求1所述的一种多波段双偏振雷达的冰晶凇附增长估算方法，其特征在于：所述设置冰晶粒子的参数，生成不同冰晶凇附质量分数的冰晶粒子宏微观参数数据集包括：利用步骤S1得到的冰晶粒子位置和冰水含量作为初始值，设置冰晶粒子的形状和有效粒子半径参数，生成基准冰晶粒子的宏微观参数数据集；调整基准冰晶粒子的凇附质量分数，生成不同冰晶凇附质量分数的冰晶粒子宏微观参数数据集。4.根据权利要求1所述的一种多波段双偏振雷达的冰晶凇附增长估算方法，其特征在于：所述利用T矩阵方法计算不同凇附质量分数冰晶粒子的后向散射特征参数包括：将电磁波的入射和散射场描述为矢量球面波函数，并通过T矩阵将函数的展开项进行关联；将冰晶粒子的尺度参数、形...

【专利技术属性】
技术研发人员：罗继成，
申请(专利权)人：成都远望科技有限责任公司，
类型：发明
国别省市：