一种基于概率模型的飞机结冰强度指数预测方法
【技术领域】
[0001] 本发明属于飞机在结冰气象环境下的结冰强度指数的预测方法，具体涉及一种基 于概率模型的飞机结冰强度指数预测方法，通过计算得到结冰强度指数，为飞行员在结冰 气象环境下选择合适有效的驾驶操作提供指导，以提高飞机飞行的安全性。
【背景技术】
[0002] 当飞机在含过冷水滴的云层中飞行时，过冷水滴会在飞机部件迎风表面迅速冻结 并累积成冰，严重威胁到飞机的飞行安全。轻度结冰会降低飞机的飞行性能，导致飞机升力 下降、阻力增加，进而造成飞行姿态控制困难;严重结冰能导致飞机在小迎角下出现失速， 甚至造成机毁人亡的惨剧。据资料统计，飞机在飞行中因结冰问题而导致空难事故的概率 超过15%，近年来也多次发生因飞机结冰而造成的重大空难事故。如2009年6月法国A330客 机在大西洋上空飞行时遇到恶劣天气发生结冰引起飞机坠毁，造成228人遇难。因此，飞机 结冰问题是飞机设计重点关注的问题之一。
[0003] 目前针对飞机结冰强度的预测主要形成了统计预报方程、假霜点温度经验公式、_ 8D方法、改进的-8D方法以及积冰强度估算方法等预测方法。比如美国Dan i e 1 Corne 11等人 利用液态水含量、温度及云滴的中位数体积等参数建立积冰严重性指数来评价结冰强度。 美国国家大气科学中心考虑飞行高度上大气温度和露点温度等因素建立了RA0B积冰预报 方案。国内学者用统计预报方程和假霜点温度经验公式对某型飞机某次严重积冰事件进行 了分析，以及开发了飞机积冰预报应用软件等。
[0004] 但需要指出的是，飞机在含过冷水滴云层中飞行时，飞行速度、高度以及液态水含 量、气压等众多影响飞机结冰的参数均在一定范围内波动。同时，传感器在这些参数的测量 上也存在测量误差。这些原因致使飞机结冰预测模型中的参数具有客观的随机不确定性， 而模型中某些关键参数的微小波动可能会导致模型结果输出的巨大变动。在目前所有的结 冰预测模型和分析方法中，均把这些参数处理为确定值，忽略了由于参数的波动以及测量 误差对预测模型结果输出的影响。
【发明内容】
[0005] 要解决的技术问题
[0006] 为了避免现有技术的不足之处，本发明提出一种基于概率模型的飞机结冰强度指 数预测方法，克服现有飞机结冰强度预测技术中未考虑参数自身波动以及参数测量误差的 不足，以期在飞机结冰强度的预测上达到更高的精度，为飞行员在结冰气象环境下选择合 适有效的驾驶操作提供指导，以提高飞机飞行的安全性。
[0007] 技术方案
[0008] -种基于概率模型的飞机结冰强度指数预测方法，其特征在于步骤如下：
[0009] 步骤1、建立飞机结冰因子的预报方程：
[0010] j = 4.9335 + 0.002016r+0.0081-4.4358 x 1+ 0.283f
[0011] 其中：V为飞行速度、单位km · h'H为飞行高度、单位km，P〇和P分别为0°C和Η高度 上的气压、单位hPa，Td为飞行高度上的露点温度，Τ为飞行高度的气温，a、b为参数，具体为a = 7.5，b = 237.3;
[0012] 步骤2、建立飞机结冰事件的判别公式：
[0013]
[0014]当&>0反映飞机不结冰，&<0反映飞机会出现结冰；
[0015] ③将
[0016] 步骤3、建立&<0的联合概率密度函数hx(x):
[0017] 将判别公式中的参数￥、？、？〇、了<1、1'和!1视为随机变量，以11^2^ 334^5和16进行表 征;随机变量X对应的的概率密度函数为f ( X )，以向量X来代替变量(XI，X2，X3，X4，X5，X6 )，则 变量的联合概率密度函数hx (X)为Μ-V) = ?./: (-V,); ζ=1
[0018] 步骤4:对联合概率密度函数在&<0的变量空间积分得到结冰概率4叫<D)
[0019] & 构 < 〇) (功办
[0020] 其中：1[ ·]是指示函数，满y!
[0021] 步骤5:采用Monte Carlo法对积分公式进行抽样模拟以得到飞机结冰概率值 4积<〇)。
[0022] 有益效果
[0023] 本发明提出的一种基于概率模型的飞机结冰强度指数预测方法，采用均匀分布、 正态分布来处理飞行速度和高度、飞行高度的气压、露点温度和气温等参数的测量以及分 析中的不确定性，并采用Monte Carlo的抽样方法来进行模拟计算，最终获得结冰强度指数 的概率值，为飞行员提供合适的飞行情况，指导飞行员合理驾驶。
[0024]本发明克服现有飞机结冰强度预测技术中未考虑参数自身波动以及参数测量误 差的不足，改变在结冰预测模型中，将￥、？、?()3(1、1'、!1等参数均为确定值的传统做法，充分考 虑乂、？、？()3 (1、1'、!1等参数的随机不确定性，以及随机不确定性的传递也导致事件61具有客观 的随机不确定性，说明飞机结冰现象(&<0)以概率出现，所以采用概率论的角度来分析事 件(〇:<0)，以达到合理地评估飞机在恶劣飞行气象下出现结冰的可能性程度，使得在飞机 结冰强度的预测上达到更高的精度，为飞行员在结冰气象环境下选择合适有效的驾驶操作 提供指导，以提高飞机飞行的安全性。
【附图说明】
[0025]图l:Monte Carlo法估计概率模型下结冰强度指数的流程图
【具体实施方式】
[0026] 现结合实施例、附图对本发明作进一步描述：
[0027] 飞机结冰因子y的预报方程如（1)式所示，可用来评估飞机在恶劣飞行气象环境下 是否结冰。
[0028] y=4.9335+0.002016V+0.00818(P-P〇)-4.4358R+0.2839H (1)
[0029] 其中V为飞行速度(km · IT1)，H为飞行高度(km)，Po、P分别为0°C和Η高度上的气压 (hPa)，R是Η高度上的相对湿度（％ )，其可以通过式(2)来求解。
[0030] Μ = lQtlb^i-r^l'+T^h+T> (2)
[0031] 式中，Td为飞行高度上的露点温度，T为飞行高度的气温，a、b为参数，具体为a = 7.5,b = 237.3〇
[0032] 研究表明：当y-2.98 2 0，飞机无结冰;y-2.98<0，飞机有结冰。因此可以建立飞机 结冰事件的判别公式，如(3)式所示。
[0033]
[0034] 在(3)式中，事件(GiX))反映飞机不结冰，而事件(0:<0)说明飞机会出现结冰的 现象。
[0035] 在传统的结冰预测模型中，V、P、Po、Td、T、H等参