一种强对流天气下飞机安全航迹规划方法与流程

技术特征：
1.一种强对流天气下飞机安全航迹规划方法，其特征在于，依次包括以下步骤：a：使用网格法建立三维空域模型，三维空域模型中的x轴、y轴和z轴分别表示经度、纬度和高度；b：将三维空域模型按照设定的长度划分为若干大小相同的立方体，将每一个立方体视为一个网格，然后对划分出的所有网格依次进行编号；c：获取气象数据并对气象数据进行垂直累积液态水含量计算，然后根据垂直累积液态水含量的计算结果由低到高将强对流天气划分为6个恶劣天气等级；d：初始化蚁群算法中搜索路径的参数，利用蚁群算法在建立的三维空域模型内寻找全局最优路径，所得到的全局最优路径即为强对流天气下飞机的安全航迹；在使用蚁群算法时，利用不同坐标位置处的恶劣天气等级计算蚂蚁的转移概率，并在每个蚁群搜索周期结束后更新网格的信息素浓度；其中，参数包括网格中蚂蚁的初始个数m、周期的最优蚂蚁序号k、网格序号s、当前蚂蚁坐标g(x,y,z)、启发式因子η、信息素挥发系数ρ、路径上对应网格的信息素浓度τ、蚂蚁禁忌表tabu以及起始点和目标点所在的网格；网格序号s和当前蚂蚁坐标g(x,y,z)的映射关系由下列公式确定：x
i
＝(s
i
‑
z
i
*n
x
*n
y
‑
1)％n
x
1y
i
＝(s
i
‑
z
i
*n
x
*n
y
‑
1)％n
y
1其中，z
i
表示第i个航迹点所在的高度，i表示当前航迹点的编号，s
i
表示第i个航迹点所在的网格序号，n
x
表示三维空域模型中在x轴方向上的网格总数，n
y
表示在y轴方向上的网格总数，x
i
表示第i个航迹点所在的经度，(s
i
‑
z
i
*n
x
*n
y
‑
1)％n
x
表示(s
i
‑
z
i
*n
x
*n
y
‑
1)对n
x
取余，其中％为取余符号，y
i
表示第i个航迹点所在的纬度。2.根据权利要求1所述的强对流天气下飞机安全航迹规划方法，其特征在于：所述的步骤b中，对划分出的每一个网格，以三维空域模型坐标原点所在的网格开始，先在x轴方向上按照经度坐标由小至大的顺序，再在y轴方向上按照纬度坐标由小至大的顺序，最后在z轴方向上按照高度坐标由小至大的顺序，依次对所有的网格进行编号，分别为定义为网格1，网格2，
……
，网格m。3.根据权利要求1所述的强对流天气下飞机安全航迹规划方法，其特征在于：所述的步骤c中，当垂直累积液态水含量数值范围为(5,20]dbz时，为第1等级；当垂直累积液态水含量数值范围为(20,35]dbz时，为第2等级；当垂直累积液态水含量数值范围为(35,45]dbz时，为第3等级；当垂直累积液态水含量数值范围为(45,50]dbz时，为第4等级；当垂直累积液态水含量数值范围为(50,65]dbz时，为第5等级；当垂直累积液态水含量数值范围为(65,70]dbz时，为第6等级。4.根据权利要求1所述的强对流天气下飞机安全航迹规划方法，其特征在于，所述的步骤d包括以下具体步骤：d1：初始化蚁群算法中搜索路径的参数；将网格中蚂蚁的初始个数设置为m，将周期的最优蚂蚁序号设置为k，k＝1，2，3，
…
，m，
网格序号分别网格1至网格m，当前蚂蚁的坐标为g(x,y,z)，x,y和z分别表示经度、纬度和高度；启发式因子为θ，信息素挥发系数为ρ，路径上对应网格的信息素浓度为τ，设置蚂蚁禁忌表tabu，设置起始点网格和目标点网格，起始点网格为飞机的初始位置所在网格，目标点网格为飞机航线目的地位置所在网格；初始的信息素浓度τ由正态分布公式生成的正态分布概率密度函数确定，以航线为中心向两边呈高斯分布：正态分布公式为：其中，f(x)为正态分布概率密度函数，ψ为平衡系数，π为圆周率，σ为标准差，exp()函数为以e为底数的指数函数，x为航线的经度坐标，μ为均值；d2：利用蚁群算法在建立的三维空域模型内寻找全局最优路径；在第一个蚁群搜索周期内：第一步：以飞机所在初始位置的经度、纬度和高度坐标确定该飞机在x轴、y轴和z轴上的位置并作为起始点；然后通过转移概率公式，分别计算起始点网格内m只蚂蚁中的每一只蚂蚁到各个网格的转移概率，再根据计算出的转移概率结果中的最大值确定该最大值所对应的q1只蚂蚁及q1只蚂蚁各自所对应的下一个航迹点，再将所确定的q1只蚂蚁所在的网格(此时为起始点所在网格)放入禁忌表tabu内并对禁忌表tabu进行更新，最后将所确定的q1只蚂蚁分别移动至各自所对应的下一个航迹点，q1为自然数且q1<m；转移概率计算公式如下：其中，表示蚂蚁k在t时刻从航迹点i到j的航段(i,j)上剩余的信息素浓度，i<j；τ
ij
(t)表示t时刻航迹点i到j的航段(i,j)上的剩余信息素的浓度，α表示表示蚂蚁在运动过程中积累的信息素对蚂蚁路径选择的相对重要性，λ
ij
(t)表示t时刻随着恶劣天气等级的变化航迹点i到j所对应航段(i,j)路径被选择的概率，u
ij
(t)表示t时刻航迹点i到j所在网格处的恶劣天气的气象因子，u
ij
(t)＝{1，2，3，4，5，6}；表示避免天气阻塞触发因素对蚂蚁路径选择的相对重要性；η
ij
(t)表示t时刻航迹点i到j的航段(i,j)上的飞行可见度的启发式函数，ε
ij
(t)表示t时刻航迹点i到j的航段(i,j)上的飞机方向变化因子，ε
ij
(t)的取值范围为(0,1)，ε
ij
(t)＝1表示飞机处于机动转弯飞行；β表示飞行方向选择对蚂蚁路径选择的相对重要性，τ
is
(t)表示t时刻航迹点i到网格s的路径上的信息素浓度，h
is
(t)表示t时刻航迹点i到网格s的高度，λ
is
(t)表示t时刻随着恶劣天气等级的变化航迹点i到网格s被选择的概率，η
is
(t)表示t时刻航迹点i到网格s的飞行可见度的启发式函数，s∈allowedroutes表示所选取的网格s在可允许路线的集合中选取，allowed
k
表示下一次蚂蚁k可以选择的网格的集合，allowed
k
＝{1，2，
…
，c}，c表示自然数；第二步：对已经移动至下一个航迹点的q1只蚂蚁，利用转移概率公式，计算每只蚂蚁到所有不在禁忌表tabu内的网格的转移概率，再根据计算出的转移概率结果中的最大值确定
该最大值所对应的q2只蚂蚁及q2只蚂蚁各自所对应的下一个航迹点，再将所确定的q2只蚂蚁当前所在的网格(此时为q2只蚂蚁各自当前所在网格)放入禁忌表tabu内并对禁忌表tabu进行更新，最后将所确定的q2只蚂蚁分别移动至所对应的下一个航迹点；q2为自然数且q2<q1；第三步：对已经移动至下一个航迹点的q2只蚂蚁，利用转移概率公式，计算每只蚂蚁到所有不在禁忌表tabu内的网格的转移概率，再根据计算出的转移概率结果中的最大值确定该最大值所对应的q3只蚂蚁及q3只蚂蚁各自所对应的下一个航迹点，再将所确定的q3只蚂蚁当前所在的网格(此时为q3只蚂蚁各自当前所在网格)放入禁忌表tabu内并对禁忌表tabu进行更新，最后将所确定的q3只蚂蚁分别移动至所对应的下一个航迹点；q3为自然数且q3<q2；重复执行第三步，并在每次重复执行后，判断是否有蚂蚁已经到达目标点网格；如果有蚂蚁已经到达目标点网格，则结束该蚁群搜索周期，将已到达目标点网格的一只或多只蚂蚁所途径的网格作为该蚁群搜索周期的最优路径，并将已到达目标点网格的唯一一只蚂蚁或多只蚂蚁中的任意一只作为该蚁群搜索周期的最优蚂蚁k；如果没有蚂蚁到达目标点网格，则继续重复执行第三步，直至有蚂蚁到达目标点网格；第四步：根据信息素浓度变化公式增强该蚁群搜索周期的最优路径上的网格即蚂蚁k所途径的网格的信息素浓度，并减弱其他路径上的网格的信息素浓度；信息素浓度变化计算公式如下所示：信息素浓度变化计算公式如下所示：信息素浓度变化计算公式如下所示：其中，表示蚂蚁k在t到t d时刻在航段(i,j)上留下的信息素浓度，d为自然数，ρ是信息素挥发系数，表示蚂蚁k在t到t d时刻在航段(i,j)上留下的信息素浓度增量，信息素浓度增量，表示蚂蚁k在t时刻在航段(i,j)上留下的信息素浓度，q是信息素常数；l
k
表示蚂蚁k在该蚁群搜索周期内经过的飞行路线的长度和恶劣天气的气象因子的和，l
k
表示蚂蚁k在该蚁群搜索周期内完成的飞行路线的长度；表示蚂蚁k在t到t d时刻沿飞行路线机动转弯的飞机方向变化因子的总和，计算公式为：表示蚂蚁k在t到t d时刻在飞行变化过程中遇到的恶劣天气的气象因子的总和，计算公式为：第五步：保留变化后的最优路径和其他路径上的网格的信息素浓度，并在新的信息素浓度下重复执行第一步至第四步，直至完成第二个蚁群搜索周期内的最优路径搜索，得到第二个蚁群搜索周期下的最优路径和最优蚂蚁；
第六步：重复执行第五步，并在每次重复执行后，判断迭代次数是否到达设定的最大迭代次数；若未达到设定的最大迭代次数，则继续执行第五步进行迭代；若达到设定的最大迭代次数，则根据每一个蚁群搜索周期的最优路径，分别求取该最优路径上所有相邻两个网格中航迹点的直线距离之和，作为该蚁群搜索周期的最优路径长度；第七步：根据得到的每一个蚁群搜索周期的最优路径长度，以数值最小的最优路径长度作为全局最优路径长度，并将该数值最小的最优路径长度所对应的最优路径作为全局最优路径。5.根据权利要求4所述的强对流天气下飞机安全航迹规划方法，其特征在于：所述的步骤d2中，使用改进后的转移概率公式，选择转移概率的最高值确定对应只蚂蚁及蚂蚁各自所对应的下一个航迹点；改进后的转移概率公式如下：其中，ω为方向选择优先级参数，0<ω<1。6.根据权利要求5所述的强对流天气下飞机安全航迹规划方法，其特征在于：所述的方向选择优先级参数ω可以通过三维神经网络训练得到；训练方法如下：a：从原始航迹数据和原始气象数据中分别读取航迹数据和雷达回波数据，并将航迹数据和雷达回波数据连接起来，生成三维神经网络训练数据集，并将三维神经网络训练数据集划分为多组训练输入数据；b：建立包含子网络a和子网络b的三维神经网络模型；其中，子网络a的结构从上至下依次为:子网络a第一循环数据特征提取层、子网络a第一最大池化层、子网络a第二循环数据特征提取层、子网络a第二最大池化层、子网络a第三循环数据特征提取层、子网络a第三最大池化层、子网络a第四循环数据特征提取层、子网络a第四最大池化层、子网络a第五循环数据特征提取层和dropout层；子网络a第一循环数据特征提取层至子网络a第五循环数据特征提取层的卷积核大小均为3，卷积步长均为1，卷积核通道数依次为64、128、256、512和1024；子网络b的结构从上至下依次为：子网络b第一单次数据特征提取层，子网络b第一数据特征平均融合层，子网络b第一循环数据特征提取层，子网络b第二数据单次特征提取层，子网络b第二数据特征平均融合层，子网络b第二循环数据特征提取层，子网络b第三单次数据特征提取层，子网络b第三数据特征平均融合层，子网络b第三循环数据特征提取层，子网络b第四单次数据特征提取层，子网络b第四数据特征平均融合层和子网络b第四循环数据特征提取层和子网络b第五单次数据特征提取层；子网络b第一单次数据特征提取层至子网络b第五单次数据特征提取层的卷积核大小均为3，卷积步长均为1，卷积核通道数依次为1024、512、256、128、64；子网络b第一循环数据特征提取层至子网络b第四循环数据特征提取层的卷积核大小均为3，卷积步长均为1，卷积核通道数依次为1024、512、256、128。c：利用步骤a中生成的三维神经网络训练数据集，对步骤b中构建的三维神经网络模型
进行训练，最终得到训练后的三维神经网络；d：按照步骤a的方法，将当前雷达回波数据和当前航迹点数据输入到训练后的三维神经网络，预测出下一时刻飞机飞行到达的航迹点所在的网格的方向选择优先级参数ω的值，从而确定下一时刻飞机选择的最优飞行方向。7.根据权利要求6所述的强对流天气下飞机安全航迹规划方法，其特征在于，所述的步骤a包括以下具体步骤：a1：首先，对多个原始航迹数据依次进行采集，采集时将航迹数据中的经度、维度和高度分别作为一个维度存入三维矩阵中；然后对多个原始多普勒雷达数据中的雷达回波数据依次进行采集，采集时将雷达回波数据中的高度、距离和方位角分别作为一个维度存入三维矩阵中；再利用concat方法将相同时间段的航迹数据和雷达回波数据连接起来，按照采集顺序分别将连接后的三维矩阵以对应的原始航迹数据的名称存储为txt文本，作为预处理后的航迹和对应的雷达回波数据；每个原始航迹数据的名称均包含航迹数据和对应的雷达回波数据的记录时间；a2：依据名称中存储的航迹数据和对应的雷达回波数据的记录时间，分别计算预处理后的航迹数据和对应的雷达回波数据中相邻两个航迹和对应的雷达回波数据的时间差，若相邻两个航迹和对应的雷达回波数据的时间差大于设定的时间差阈值t，则判定这两个相邻的航迹和对应的雷达回波数据中存在航迹数据缺失，然后按照插值填补法对缺失的航迹和对应的雷达回波数据进行填补，最终得到三维神经网络的训练数据集；a3：将步骤a2中得到的三维神经网络的训练数据集中，以相邻的n个航迹和对应的雷达回波数据作为一组输入数据，将三维神经网络的训练数据集划分为多组训练输入数据。8.根据权利要求6所述的强对流天气下飞机安全航迹规划方法，其特征在于，所述的步骤b中：子网络a和子网络b中的单次数据特征提取层，均以待进行单次特征提取的输入数据作为输入，利用数据特征提取模块对航迹和对应的雷达回波数据进行单次特征提取，然后将提取得到的经单次特征提取后的特征矩阵进行输出；数据特征提取模块从上至下依次为三维卷积操作、三维批归一化操作和relu激活函数；单次特征提取时数据特征提取模块的输入输出关系为：其中，input1
n
表示待进行单次特征提取的三维神经网络的训练数据集中的第n组训练输入数据，output1
n
表示经单次特征提取后的特征矩阵中的第n组输出数据，kernel表示卷积核大小，padding表示填充值，step表示卷积运算的步长；n为自然数；循环数据特征提取层，首先以待进行循环特征提取的输入数据作为输入，利用数据特征提取模块对多普勒雷达回波数据进行一次特征提取，然后将经一次特征提取后得到的特征矩阵进行输出；然后再以经一次特征提取后得到的特征矩阵作为输入，再次利用数据特征提取模块对多普勒雷达回波数据进行二次特征提取，最终将经循环特征提取后得到的特征矩阵进行输出；一次特征提取时数据特征提取模块的输入输出关系为：
二次特征提取时数据特征提取模块的输入输出关系为：其中，input2
n
表示待进行循环特征提取的三维神经网络的训练数据集中的第n组训练输入数据，output2
n
表示经一次特征提取后的特征矩阵中的第n组输出数据，output3
n
表示经循环特征提取后的特征矩阵中的第n组输出数据。9.根据权利要求6所述的强对流天气下飞机安全航迹规划方法，其特征在于，所述的步骤c中：在使用子网络b中的数据特征平均融合层对两个待融合的特征矩阵进行融合时：首先，对两个待融合的特征矩阵进行处理，将两个特征矩阵中的雷达回波数据中特征值小于零的特征点的特征值全部替换为零，以简化融合过程中的计算步骤，提高神经网络的运行速度。然后，若两个特征矩阵的相应位置的特征值均为零，则融合后的上采样矩阵该位置的特征值还为零；若两个特征矩阵的相应位置的特征值中，一个特征矩阵的特征值是零，另一个特征矩阵的特征值非零，则融合后的上采样矩阵的该位置的特征值采用对应的非零特征值；若两个特征矩阵的相应位置的特征值都是非零，则融合之后的上采样矩阵的特征值为两个非零特征值的平均数。10.根据权利要求6所述的强对流天气下飞机安全航迹规划方法，其特征在于，所述的步骤c中：在三维神经网络进行训练时，子网络a第一循环数据特征提取层，对所输入的数据进行循环特征提取得到特征矩阵a1，并将输出的特征矩阵a1输入至子网络a第一最大池化层中；然后由子网络a第一最大池化层对特征矩阵a1进行降维操作以减少特征矩阵中特征值为零或小于零的特征点个数，并将得到的降维后的特征矩阵a2输入至子网络a第二循环数据特征提取层；然后，按照子网络a第一循环数据特征提取层和子网络a第一最大池化层的相同操作，由子网络a逐层对上一层得到的特征矩阵进行特征提取和降维操作，最后由子网络a第五循环数据特征提取层将经过五次特征提取后的特征数据输入至dropout层，由dropout层进行正则化操作；最终，将由dropout层输出的特征矩阵a10输入至子网络b第一单次数据特征提取层；子网络b第一单次数据特征提取层对输入的特征矩阵a10进行单次特征提取，然后将输出的特征矩阵b1输入至子网络b第一数据特征平均融合层；子网络b第一数据特征平均融合层将子网络b第一单次数据特征提取层输出的特征矩阵b1和子网络a第一循环数据特征提取层输出的特征矩阵a1进行融合，并将融合后的特征矩阵b2输入至子网络b的第一循环特征提取层，子网络b的第一循环特征提取层对融合后的特征矩阵进行循环特征提取后输出特征矩阵b3，然后将特征矩阵b3输入至子网络b的第二数据单次特征提取层；再由子网络b的第二数据单次特征提取层对对输入的特征矩阵b3进行单次特征提取，然后将输出的特征
矩阵b4输入至子网络b第二数据单次特征提取层；然后，按照子网络b第一单次数据特征提取层、子网络b第一数据特征平均融合层和子网络b第一循环数据特征提取层的相同操作，逐层对输入的特征矩阵进行特征提取，最后由第五单次数据特征提取层输出最终生成的特征矩阵；按照上述方法对三维神经网络进行训练，最终得到训练后的三维神经网络。

技术总结
本发明公开了一种强对流天气下飞机安全航迹规划方法，依次包括以下步骤：A：使用网格法建立三维空域模型；B：将三维空域模型按照设定的长度划分为若干网格，然后对划分出的所有网格依次进行编号；C：获取气象数据并对气象数据进行垂直累积液态水含量计算，然后根据垂直累积液态水含量的计算结果由低到高将强对流天气划分为6个恶劣天气等级；D：初始化蚁群算法中搜索路径的参数，利用蚁群算法在建立的三维空域模型内寻找全局最优路径。本发明能够在实时动态变化的强对流天气环境下，快速规划飞机安全航迹，确保飞行安全，提高飞行效率。提高飞行效率。提高飞行效率。


技术研发人员：张磊 陈会平 贾培艳 翟凯 孟坤颖 马夏青 姜帅民 肖琦 李星辰
受保护的技术使用者：中国民用航空中南地区空中交通管理局河南分局
技术研发日：2021.07.22
技术公布日：2021/9/24