机载气象雷达高空冰晶冰水含量预测方法技术

本发明专利技术提供一种机载气象雷达高空冰晶冰水含量预测方法，包括以下步骤：步骤一、建立高空冰晶三维网格化数据库；步骤二、将高空冰晶数据按照温度进行分层，并将每个温度层数据按反射率因子进行分区；步骤三、选择训练样本并形成训练样本数据集；步骤四、建立反射率因子与冰水含量之间的映射关系模型，实现特征参量矩阵和标签矩阵的多元回归重构；步骤五、利用最小二乘回归的多因变量分析估计回归系数；步骤六、对所有的温度层和反射率区间重复步骤三至步骤五，并构建回归系数集合；步骤七、根据回归系数集合对机载气象雷达回波数据的冰水含量进行预测，当预测值大于门限阈值时，雷达进行高空冰晶显示告警。行高空冰晶显示告警。行高空冰晶显示告警。

【技术实现步骤摘要】
机载气象雷达高空冰晶冰水含量预测方法


[0001]本专利技术涉及雷达
，具体涉及一种机载气象雷达高空冰晶冰水含量预测方法。

技术介绍

[0002]高空冰晶通常是飞机在高空巡航阶段飞行时经常遭遇的一种气象危害。当飞机长时间在冰晶高聚集区空域飞行时，飞机机翼结冰以后会引起阻力增加，升力下降，操控性能变差；冰晶进入发动机后融化成液态水并再结冰后慢慢增长，从而造成发动机损坏或者停车；冰晶还会进入大气数据传感器，形成堵塞从而造成相关传感器数据异常，带来飞机失控。
[0003]根据国际航空运输协会(International Air Transport Association，IATA)统计分析，2011
‑
2015年，42％的飞行失控由气象威胁引起，其中40％是由于结冰引起的；根据美国国家航空航天局(National Aeronautics and Space Administration，NASA)统计，过去20年来，航空工业记录过超过200次由于冰晶引起的涡扇发动机高空停车事故。FAA和EASA已在相关民航法规中将高空冰晶列为重点关注的气象危害。欧盟委员会于2012年启动了一项“机载HAIC研究计划”，联合空客、Rockwell Collins、Honeywell、Thales等二十余家机构共同开展HAIC科学研究；RTCA SC
‑
230小组正在制定冰晶检测的技术性能标准；机载气象雷达的两大主要厂商Rockwell Collins和Honeywell也正在开展HAIC探测技术研究。
[0004]高空冰晶一般在环境温度为
‑
5℃～
‑
55℃、海拔高度为9000～40000英尺的高空出现，通常与对流云层有关，例如孤立的积雪云层，雷暴，或甚至热带气旋。这种对流天气的强上升气流可将液态水提升数千英尺到大气中，从而形成高空冰晶。根据冰晶形成的机理，冰晶可以是不同的尺寸和形状，从简单的球形到复杂的雪花状结构。对于机载气象雷达而言，高空冰晶都具有较低的反射率因子值，很难通过反射率因子判断高空冰晶的存在或其浓度。
[0005]现有技术对高空冰晶的存在可以由特定空间点处的冰水含量(IWC)来指示，利用模型估计IWC实现冰晶预测。但是该方法需要建立冰晶粒子的分布模型，设置合理的关键参数，才能得到IWC与Z的良好的对应关系，同时该方法对IWC的估计比较单一，无法适应不同粒子特性下的IWC的估计，在不同条件下，估计误差增大。

技术实现思路

[0006]有鉴于此，本说明书实施例提供一种机载气象雷达高空冰晶冰水含量预测方法，为高空巡航阶段的飞机飞行安全提供保障。
[0007]本说明书实施例提供以下技术方案：一种机载气象雷达高空冰晶冰水含量预测方法，包括以下步骤：
[0008]步骤一、建立高空冰晶三维网格化数据库；
[0009]步骤二、将高空冰晶数据按照温度进行分层，并将每个温度层数据按反射率因子
进行分区；
[0010]步骤三、选择训练样本并形成训练样本数据集；
[0011]步骤四、建立反射率因子与冰水含量之间的映射关系模型，实现特征参量矩阵和标签矩阵的多元回归重构；
[0012]步骤五、利用最小二乘回归的多因变量分析估计回归系数；
[0013]步骤六、对所有的温度层和反射率区间重复步骤三至步骤五，并构建回归系数集合；
[0014]步骤七、根据回归系数集合对机载气象雷达回波数据的冰水含量进行预测，当预测值大于门限阈值时，雷达进行高空冰晶显示告警。
[0015]进一步地，步骤二包括：
[0016]步骤2.1、将高空冰晶的三维网格化数据，按照温度进行分层，确定温度分层的间隔或者范围，确保各个温度层内粒子分布的一致性；
[0017]步骤2.2、将温度分层的数据，再根据反射率值进行分区，使各反射率区间的反射率与冰水含量的关系呈线性关系。
[0018]进一步地，步骤三包括：
[0019]取出反射率区间对应的数据；
[0020]统计反射率区间的冰水含量数据中值；
[0021]在冰水含量中值基础上进行插值计算，形成样本选择所需的双门限；
[0022]取出反射率区间双门限内L个与冰水含量中值最近邻的数据作为有效训练样本，其中L为训练样本最大个数。
[0023]进一步地，步骤四包括：
[0024]步骤4.1、将训练样本中的反射率、高度或者温度参数作为特征参量矩阵X，将冰水含量作为标签矩阵Y；
[0025]步骤4.2、建立多元线性回归函数Y＝XW+F，其中，W为标签函数与特征参量矩阵的最终的回归系数，每个特征参量对应一个系数权值，F为残差系数。
[0026]进一步地，步骤五具体为：
[0027]步骤5.1、利用最小二乘回归的多因变量分析估计回归系数W和残差系数F，并将回归系数W和残差系数F按照温度层和反射率区间进行保存；
[0028]步骤5.2、将反射率区间未参与训练的样本作为测试样本，利用回归系数W对数据进行测试，对预测出的冰水含量与真实冰水含量进行对比，将预测误差引入残差中，实现回归系数的优化，并更新反射率区间的回归系数。
[0029]进一步地，步骤七具体为：
[0030]步骤7.1、对各距离门回波数据进行处理得到反射率；
[0031]步骤7.2、根据载机高度、波束宽度和波束指向角参数，计算各距离门数据对应高度；
[0032]步骤7.3、根据载机上装载的温度传感器获取载机周围环境温度，并利用温度
‑
高度变化曲线，计算该距离门数据对应的温度；
[0033]步骤7.4、根据各距离门计算得到的温度和反射率值，提取系数集合中相应区间的回归系数；
[0034]步骤7.5、基于多特征参量的冰水含量值估计，根据公式IWC＝XW+F，计算各距离门数据对应的冰水含量值。
[0035]与现有技术相比，本说明书实施例采用的上述至少一个技术方案能够达到的有益效果至少包括：本专利技术提升了IWC的估计精度及稳健性，适用于不同条件下的高空冰晶危险性评估，并可在雷达显示画面中对具有冰晶威胁的区域进行告警显示，可有效降低高空巡航阶段飞机结冰的风险，提升了飞行安全性。
附图说明
[0036]为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
[0037]图1是本专利技术实施例的流程示意图；
[0038]图2是本专利技术实施例对某个区间数据的最小二乘回归参数估计流程图。
具体实施方式
[0039]下面结合附图对本申请实施例进行详细描述。
[0040]需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种机载气象雷达高空冰晶冰水含量预测方法，其特征在于，包括以下步骤：步骤一、建立高空冰晶三维网格化数据库；步骤二、将高空冰晶数据按照温度进行分层，并将每个温度层数据按反射率因子进行分区；步骤三、选择训练样本并形成训练样本数据集；步骤四、建立反射率因子与冰水含量之间的映射关系模型，实现特征参量矩阵和标签矩阵的多元回归重构；步骤五、利用最小二乘回归的多因变量分析估计回归系数；步骤六、对所有的温度层和反射率区间重复步骤三至步骤五，并构建回归系数集合；步骤七、根据回归系数集合对机载气象雷达回波数据的冰水含量进行预测，当预测值大于门限阈值时，雷达进行高空冰晶显示告警。2.根据权利要求1所述的机载气象雷达高空冰晶冰水含量预测方法，其特征在于，所述步骤二包括：步骤2.1、将高空冰晶的三维网格化数据，按照温度进行分层，确定温度分层的间隔或者范围，确保各个温度层内粒子分布的一致性；步骤2.2、将温度分层的数据，再根据反射率值进行分区，使各反射率区间的反射率与冰水含量的关系呈线性关系。3.根据权利要求2所述的机载气象雷达高空冰晶冰水含量预测方法，其特征在于，所述步骤三包括：取出所述反射率区间对应的数据；统计所述反射率区间的冰水含量数据中值；在冰水含量中值基础上进行插值计算，形成样本选择所需的双门限；取出所述反射率区间双门限内L个与冰水含量中值最近邻的数据作为有效训练样本，其中L为训练样本最大个数。4.根据权利要求3所述的机载气象雷达高空冰晶冰水含量预测方法，其...

【专利技术属性】
技术研发人员：陈娟，孟凡旺，蒋银丰，
申请(专利权)人：中国航空工业集团公司雷华电子技术研究所，
类型：发明
国别省市：