CN112927562B - 一种基于指定高度限制的航空器飞行高度剖面计算方法 - Google Patents
一种基于指定高度限制的航空器飞行高度剖面计算方法 Download PDFInfo
- Publication number
- CN112927562B CN112927562B CN202110116380.5A CN202110116380A CN112927562B CN 112927562 B CN112927562 B CN 112927562B CN 202110116380 A CN202110116380 A CN 202110116380A CN 112927562 B CN112927562 B CN 112927562B
- Authority
- CN
- China
- Prior art keywords
- height
- altitude
- aircraft
- flight
- sector
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/003—Flight plan management
- G08G5/0039—Modification of a flight plan
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供了一种基于指定高度限制的航空器飞行高度剖面计算方法，获取航班飞行计划、空域扇区结构和高度限制；分析限制点的高度与其前一航段飞行高度的关系，确定是否存在高低扇区转变；若不存在，依次计算当前高度飞行到指定高度所需的最小水平距离、高度限制点到其所属扇区进入边界的距离，在计算高度剖面时通过两个距离的关系决定改变高度和平飞的先后顺序；若存在，基于上述处理逻辑，考虑高度限制点与目的机场的距离，当两者距离小于设定参数，在计算高度剖面过程中优先下降到指定高度，再保持平飞。本发明方法充分考虑了空管运行的规则和实际情况，能够提高预战术阶段飞行航迹预测经过扇区的准确性，为精准的扇区流量预测提供技术支持。
Description
技术领域
本发明属于空中交通管理领域，特别涉及一种基于指定高度限制的航空器飞行高度剖面计算方法。
背景技术
随着航空运输量的持续增加，空中交通精细化运行与管理的需求愈发迫切。在交通流密集空域内仍然采用基于飞行计划结合基于经验的间隔调配逐渐不能满足空管运行的需要。为此国际民航界提出将基于四维航迹运行作为未来空中交通运行的新理念，它以航空器飞行全生命周期的四维航迹为中心，空管部门、航空公司、航空器、机场之间共享、协商和管理动态航迹，实现飞行与管制之间的协同决策。这也对空管***的4D航迹预测能力提出了更高的要求。
精确的4D飞行航迹是空中交通管理***诸多决策支持工具得以高效使用的基础，可用于空中交通流量预测与管理的各个阶段以及进行飞行冲突探测、飞机进港排序、辅助管制指挥等，提高空域资源利用率。4D飞行航迹包括：水平飞行航迹、高度剖面和速度剖面，出于计算方便，通常将这三个剖面独立处理，并最终进行整合。目前***中常用的4D航迹预测方法是采用基于航班飞行计划和航空器性能模型相结合的方式进行计算的。在该方法中，高度剖面通常相对简单，即航空器爬升到巡航高度，然后平飞到下降顶点，再下降到目的机场。巡航高度使用飞行计划中的高度或者历史经验高度，并考虑空管发布的运行限制、管制移交协议等高度限制，导致预测的高度剖面与实际有较大偏差，进而会影响速度剖面和飞行时间的计算。特别是流量预测过程中，由于高空扇区和低空扇区可能分属于不同的管制单位，高度剖面的偏差将导致扇区流量预测的不准确，影响流量管理策略的制定，导致整体运行效率的降低。目前，尚缺少一种充分考虑空管运行实际和高度限制的高度剖面准确计算方法，以支撑精确4D航迹的预测。
发明内容
发明目的：本发明要解决的技术问题是针对空管***向以4D航迹为核心的精细化运行发展需求，从准确预测航空器4D飞行航迹的角度出发，研究和发明一种基于指定高度限制的航空器飞行高度剖面计算方法。
为了解决上述技术问题，本发明公开了一种基于指定高度限制的航空器飞行高度剖面计算方法，包括如下步骤：
步骤1：获取飞行计划、空域扇区结构、进离场程序、航路点的高度限制等基础信息；
步骤2：根据离场程序计算起飞机场到第一个航路点的高度剖面，并更新航迹起始点为第一个航路点；
步骤3：从航迹起始点开始向后找到下一个高度限制点，计算航迹起始点到该高度限制点前一个扇区进入边界的高度剖面，并更新航迹起始点；
步骤4：通过航迹起始点的飞行高度、限制点指定高度和空域扇区结构之间的关系，判断航班是否由高空扇区向低空扇区转换；
步骤5：根据步骤4的判断结果，确定平飞和爬升或下降过程的先后顺序：
步骤5-1：当不存在由高空扇区向低空扇区转换情况时，根据所需爬升或下降距离和扇区边界的关系，确定平飞和爬升或下降过程的先后顺序；
步骤5-2：当存在由高空扇区向低空扇区转换情况时，通过评估高度限制点与目的机场的距离来判断是否需要优先执行下降过程；当两者之间的距离小于等于设定的参数LT时，则先下降到指定高度，再保持平飞；否则，根据所需下降距离和扇区边界的关系，确定平飞和下降过程的先后顺序；通过高度限制点与目的机场的距离来判断是否需要优先执行下降过程，使得高度剖面的计算与运行实际更加贴合。
步骤6：根据步骤5确定的高度剖面变化过程，计算航迹起始点到高度限制点的高度剖面，并更新航迹起始点；
步骤7：判断航迹起始点是否为最后一个航路点，如果是则执行步骤8；如果不是执行步骤3；
步骤8：根据进场程序计算航迹起始点到目的机场的高度剖面。
在一种实现方式中，所述步骤1中，飞行计划中包含了计算高度剖面所需的起飞机场、航路点、降落机场和巡航高度信息，还包含了航空器机型；空域扇区结构用于确定高空扇区和低空扇区，以及航路点不同高度所对应的扇区类型；进离场程序用于计算离场和进场阶段的高度剖面；航路点的高度限制来源于移交协议、流量管理或管制提出的限制；飞行计划中的第一个航路点和最后一个航路点分别对应于离场点和进场点，航路点的高度限制根据空管运行规则，在6000米以下。
在一种实现方式中，所述步骤2中，根据起飞机场的离场程序计算确定航空器从跑道到第一个航路点的飞行距离，爬升过程的计算分为加速爬升、低等表爬升和高等表爬升3个子过程，高度剖面计算公式为
式中：h为高度，t为时间，m为航空器质量，VTas为航空器真空速，g为重力加速度，T为航空器推力，D为航空器阻力；
为能量分配系数，表示按照选定的速度进行爬升时，用于爬升的推力与用于加速的推力的比值；通过航空器机型能够从欧控发布的BADA数据库获得航空器在不同阶段和高度对应的航空器推力T、航空器阻力D和航空器真空速VTas；爬升过程采用逐步迭代计算，飞机的重量在每个时间间隔
内保持不变，确定每个时间间隔
内的航空器真空速VTas，根据上述公式计算出各个时间间隔
内的爬升率，依次完成起飞机场到第一个航路点爬升过程的高度剖面计算；如果在第一个航路点之前达到指定的高度限制，则保持平飞状态至第一个航路点。
在一种实现方式中，所述步骤3中，从航迹起始点开始向后找到下一个高度限制点，航迹起始点到该高度限制点前一个扇区进入边界的高度剖面计算公式和步骤2中相同，处理逻辑如下：当已完成计算的高度剖面未曾达到巡航高度，即航班尚未进入巡航状态时，如果该高度限制点与起飞机场水平距离小于200公里，则先爬升至该限制高度，然后保持平飞；否则先爬升到巡航高度，再保持平飞；如果航班已经进入巡航阶段，则保持平飞状态至该高度限制点前一个扇区的进入边界。考虑航空器是否已经进入巡航阶段这个因素，确保高度剖面计算时使得航空器尽可能利用巡航高度飞行，与运行实际更加符合。
在一种实现方式中，所述步骤4中，当航迹起始点处于高空扇区，并且限制点的指定高度低于航迹起始点的飞行高度时，则根据空域扇区结构判断该指定高度是否处于低空扇区空域内，以此来判定航空器是否将由高空扇区飞向低空扇区。
在一种实现方式中，所述步骤5-1包括：
5-1a：计算当前高度爬升或下降到指定高度所需的最小水平距离Lmin，考虑到航空器需要提前满足高度限制，
Lmin＝Lf+20km
其中Lf为当前高度爬升或下降到指定高度的水平距离，根据步骤2中高度剖面计算公式计算获得当前高度爬升或下降到指定高度的爬升率或下降率，再根据当前高度爬升或下降到指定高度的高度差计算出当前高度爬升或下降到指定高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算得到Lf；
5-1b：计算高度限制点到其所属扇区进入边界的距离L，即高度限制点到边界线的距离；
5-1c：当L≥Lmin时，说明能够在最后一个扇区完成高度改变过程，则航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降或爬升过程到指定高度，然后平飞至高度限制点；
5-1d：当L＜Lmin时，令Lmin＝min{Lf+L，Lc}，航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降或爬升过程到指定高度，然后平飞至高度限制点；其中Lc为航迹起始点与高度限制点之间的水平距离。
在一种实现方式中，所述步骤5-2中参数LT设置为200公里；当高度限制点与目的机场的水平距离大于200公里时，高度剖面变化过程包括：
5-2a：计算当前高度下降到指定高度所需的最小水平距离Lmin，考虑到航空器需要提前满足高度限制，
Lmin＝Lf+20km
其中Lf为当前高度下降到指定高度的水平距离，根据步骤2中高度剖面计算公式计算获得当前高度下降到指定高度的下降率，再根据当前高度下降到指定高度的高度差计算出当前高度下降到指定高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算得到Lf；
5-2b：计算高度限制点到其所属扇区进入边界的距离L，即高度限制点到边界线的距离；
5-2c：当L≥Lmin时，说明能够在最后一个扇区完成高度改变过程，则航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降过程到指定高度，然后平飞至高度限制点；
5-2d：当L＜Lmin时，令Lmin＝min{Lf+L，Lc}，航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降过程到指定高度，然后平飞至高度限制点；其中Lc为航迹起始点与高度限制点之间的水平距离。
步骤3、步骤4和步骤5中，在飞行状态确定和高度剖面计算过程中，考虑了空域扇区边界这一重要运行因素。在确定平飞、下降或爬升过程先后顺序时，保持在跨扇区边界进行管制移交时的平飞状态，增加高度剖面的合理性和精细化程度。
在一种实现方式中，所述步骤6中根据步骤5确定的高度剖面变化过程，计算航迹起始点到高度限制点的高度剖面计算公式和步骤2中相同。
在一种实现方式中，所述步骤8中进场阶段的高度剖面计算公式和步骤2中相同，计算并确定从当前高度下降到目的机场高度所需的水平距离Ld，航空器从当前高度保持平飞至距离目的机场Ld的位置，然后开始下降到目的机场；根据步骤2中高度剖面计算公式计算获得从当前高度下降到目的机场高度的下降率，再根据当前高度下降到机场高度的高度差计算出当前高度下降到机场高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算获得Ld。
本发明具有如下技术效果：
1、本发明为预战术阶段的航空器高度剖面计算提供一种快速实现方法；在飞行状态确定和高度剖面计算过程中，考虑了空域扇区边界这一重要运行因素。在确定平飞、下降或爬升过程先后顺序时，保持在跨扇区边界进行管制移交时的平飞状态，增加高度剖面的合理性和精细化程度。
2、本发明为四维航迹的精确预测提供了技术支持，当从航迹起始点开始向后找到下一个高度限制点，计算航迹起始点到该高度限制点前一个扇区进入边界的高度剖面时考虑航空器是否已经进入巡航阶段这个因素，确保高度剖面计算时使得航空器尽可能利用巡航高度飞行，与运行实际更加符合；当航班存在由高空扇区向低空扇区转换情况时，通过高度限制点与目的机场的距离来判断是否需要优先执行下降过程，使得高度剖面的计算与运行实际更加贴合。
附图说明
下面结合附图和具体实施方式对本发明做更进一步的具体说明，本发明的上述和/或其他方面的优点将会变得更加清楚。
图1为本发明的方法流程图；
图2为不存在由高空扇区向低空扇区转换情况的高度剖面；
图3为由高空扇区向低空扇区转换情况的高度剖面。
具体实施方式
下面将结合附图，对本发明的实施例进行描述。
如图1所示，本发明公开了一种基于指定高度限制的航空器飞行高度剖面计算方法，包括如下步骤：
步骤1：获取飞行计划、空域扇区结构、进离场程序、航路点的高度限制等基础信息；飞行计划中包含了计算高度剖面所需的起飞机场、航路点、降落机场和巡航高度信息，还包含了航空器机型；空域扇区结构用于确定高空扇区和低空扇区，以及航路点不同高度所对应的扇区类型；进离场程序用于计算离场和进场阶段的高度剖面；航路点的高度限制来源于移交协议、流量管理或管制提出的限制等，飞行计划中的第一个航路点和最后一个航路点分别对应于离场点和进场点，根据空管运行规则会有指定的高度限制，通常在6000米以下；
步骤2：根据离场程序计算起飞机场到第一个航路点的高度剖面，并更新航迹起始点为第一个航路点；根据起飞机场的离场程序计算确定航空器从跑道到第一个航路点的飞行距离，爬升过程的计算分为加速爬升、低等表爬升段和高等表爬升3个子过程，计算公式为
式中：h为高度，t为时间，m为航空器质量；VTas为航空器真空速；g为重力加速度，T为航空器推力，D为航空器阻力；
为能量分配系数，表示按照选定的速度进行爬升时，用于爬升的推力与用于加速的推力的比值。欧控发布的Base OfAircraft Data(BADA)数据库中包含了常用机型的性能参数，包括在不同阶段和高度对应的推力、阻力和速度，通过航空器机型能够从BADA数据库获得航空器在不同阶段和高度对应的航空器推力T、航空器阻力D和航空器真空速VTas。爬升过程采用逐步迭代计算的思路，将在每个时间间隔
内假设飞机的重量保持不变，确定每个时间间隔
内的VTas，再根据公式(1)计算出各个时间间隔
内的爬升率，依次完成起飞机场到第一个航路点爬升过程的高度剖面计算，时间间隔
通常为8秒。如果在第一个航路点之前达到指定的高度限制，则保持平飞状态至第一个航路点。
步骤3：从航迹起始点开始向后找到下一个高度限制点，并计算航迹起始点到该高度限制点前一个扇区进入边界的高度剖面，并更新航迹起始点。本步骤高度剖面计算过程仍然利用步骤2的计算公式，但处理逻辑如下：1)当已完成计算的高度剖面未曾达到巡航高度，即航班尚未进入巡航状态时，如果该高度限制点与起飞机场水平距离小于200公里，则先爬升至该限制高度，然后保持平飞；否则先爬升到巡航高度，再保持平飞；2)如果航班已经进入巡航阶段，则保持平飞状态至该高度限制点前一个扇区的进入边界。
步骤4：通过航迹起始点的飞行高度、限制点指定高度和空域扇区结构之间的关系，判断航班是否由高空扇区向低空扇区转换。当航迹起始点处于高空扇区，并且限制点的指定高度低于航迹起始点的飞行高度时，则需要根据空域扇区结构判断该指定高度是否处于低空扇区空域内，以此来判定航空器是否将由高空扇区飞向低空扇区；高空扇区和低空扇区之间的边界在每个管制区内是不同的，但通常为7800米或6600米。
步骤5：根据步骤4的判断结果，确定平飞和爬升或下降过程的先后顺序：
步骤5-1：当不存在由高空扇区向低空扇区转换情况时，如图2所示，根据所需爬升或下降距离和扇区边界的关系，确定平飞和爬升或下降过程的先后顺序。具体计算过程如下：
5-1a：计算当前高度爬升或下降到指定高度所需的最小水平距离Lmin，考虑到航空器通常需要提前满足高度限制，
Lmin＝Lf+20km
其中Lf为当前高度爬升或下降到指定高度的距离，可根据公式(1)计算得到当前高度爬升或下降到指定高度的爬升率或下降率，再根据当前高度爬升或下降到指定高度的高度差计算出当前高度爬升或下降到指定高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算得到Lf；
5-1b：计算高度限制点到其所属扇区进入边界的距离L；
5-1c：当L≥Lmin时，说明可以在最后一个扇区完成高度改变过程，则航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降或爬升过程到指定高度，然后平飞至高度限制点；
5-1d：当L＜Lmin时，令Lmin＝min{Lf+L，Lc}，航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降或爬升过程到指定高度，然后平飞至高度限制点；其中Lc为航迹起始点与高度限制点之间的水平距离。
步骤5-2：当存在由高空扇区向低空扇区转换情况时，如图3所示，通过评估高度限制点与目的机场的距离来判断是否需要优先执行下降过程；当高度限制点与目的机场的水平距离小于等于200公里(设定的参数LT)时，航空器的高度剖面变化过程为：先下降到指定高度，然后保持平飞至高度限制点；当高度限制点与目的机场的水平距离大于200公里时，高度剖面变化过程的确定与步骤5-1一致，包括：
5-2a：计算当前高度下降到指定高度所需的最小水平距离Lmin，考虑到航空器需要提前满足高度限制，
Lmin＝Lf+20km
其中Lf为当前高度下降到指定高度的水平距离，根据步骤2中高度剖面计算公式计算获得当前高度下降到指定高度的下降率，再根据当前高度下降到指定高度的高度差计算出当前高度下降到指定高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算得到Lf；
5-2b：计算高度限制点到其所属扇区进入边界的距离L；
5-2c：当L≥Lmin时，说明能够在最后一个扇区完成高度改变过程，则航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降过程到指定高度，然后平飞至高度限制点；
5-2d：当L＜Lmin时，令Lmin＝min{Lf+L，Lc}，航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降过程到指定高度，然后平飞至高度限制点；其中Lc为航迹起始点与高度限制点之间的水平距离。
步骤6：根据步骤5确定的高度剖面变化过程，计算航迹起始点到高度限制点的高度剖面，并更新航迹起始点。计算航迹起始点到高度限制点的高度剖面计算公式和步骤2中相同。
步骤7：判断航迹起始点是否为最后一个航路点，如果是则执行步骤8；如果不是执行步骤3。
步骤8：根据进场程序计算航迹起始点到目的机场的高度剖面。进场阶段的高度剖面与离场高度剖面是基本对称的，仍然采用公式(1)进行计算。首先计算并确定从当前高度下降到目的机场所需的水平距离Ld，航空器从当前高度保持平飞至距离目的机场Ld的位置，然后开始下降到目的机场。根据步骤2中高度剖面计算公式计算获得从当前高度下降到目的机场高度的下降率，再根据当前高度下降到机场高度的高度差计算出当前高度下降到机场高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算获得Ld。
本发明提供了一种基于指定高度限制的航空器飞行高度剖面计算方法，具体实现该技术方案的方法和途径很多，以上所述仅是本发明的具体实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有技术加以实现。
Claims (7)
1.一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，包括如下步骤：
步骤1：获取飞行计划、空域扇区结构、进离场程序和航路点的高度限制基础信息；
步骤2：根据离场程序计算起飞机场到第一个航路点的高度剖面，并更新航迹起始点为第一个航路点；
步骤3：从航迹起始点开始向后找到下一个高度限制点，计算航迹起始点到该高度限制点前一个扇区进入边界的高度剖面，并更新航迹起始点；
步骤4：通过航迹起始点的飞行高度、限制点指定高度和空域扇区结构之间的关系，判断航班是否由高空扇区向低空扇区转换；
步骤5：根据步骤4的判断结果，确定平飞和爬升过程的先后顺序，或者平飞和下降过程的先后顺序：
步骤5-1：当不存在由高空扇区向低空扇区转换情况时，根据所需爬升或下降距离和扇区边界的关系，若需爬升，确定平飞和爬升过程的先后顺序，若需下降，确定平飞和下降过程的先后顺序；
步骤5-2：当存在由高空扇区向低空扇区转换情况时，通过评估高度限制点与目的机场的距离来判断是否需要优先执行下降过程；当两者之间的距离小于等于设定的参数LT时，则先下降到指定高度，再保持平飞；否则，根据所需下降距离和扇区边界的关系，确定平飞和下降过程的先后顺序；
步骤6：根据步骤5确定的高度剖面变化过程，计算航迹起始点到高度限制点的高度剖面，并更新航迹起始点；
步骤7：判断航迹起始点是否为最后一个航路点，如果是则执行步骤8；如果不是执行步骤3；
步骤8：根据进场程序计算航迹起始点到目的机场的高度剖面；
所述步骤5-1包括：
5-1a：计算当前高度爬升或下降到指定高度所需的最小水平距离Lmin，考虑到航空器需要提前满足高度限制，
Lmin＝Lf+20km
其中Lf为当前高度爬升或下降到指定高度的水平距离，根据步骤2中高度剖面计算公式计算获得当前高度爬升或下降到指定高度的爬升率或下降率，再根据当前高度爬升或下降到指定高度的高度差计算出当前高度爬升或下降到指定高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算得到Lf；
5-1b：计算高度限制点到其所属扇区进入边界的距离L；
5-1c：当L≥Lmin时，说明能够在最后一个扇区完成高度改变过程，则航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降或爬升过程到指定高度，然后平飞至高度限制点；
5-1d：当L＜Lmin时，令Lmin＝min{Lf+L,Lc}，航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降或爬升过程到指定高度，然后平飞至高度限制点；其中Lc为航迹起始点与高度限制点之间的水平距离；
所述步骤5-2中参数LT设置为200公里；当高度限制点与目的机场的水平距离大于200公里时，高度剖面变化过程包括：
5-2a：计算当前高度下降到指定高度所需的最小水平距离Lmin，考虑到航空器需要提前满足高度限制，
Lmin＝Lf+20km
其中Lf为当前高度下降到指定高度的水平距离，根据步骤2中高度剖面计算公式计算获得当前高度下降到指定高度的下降率，再根据当前高度下降到指定高度的高度差计算出当前高度下降到指定高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算得到Lf；
5-2b：计算高度限制点到其所属扇区进入边界的距离L；
5-2c：当L≥Lmin时，说明能够在最后一个扇区完成高度改变过程，则航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降过程到指定高度，然后平飞至高度限制点；
5-2d：当L＜Lmin时，令Lmin＝min{Lf+L,Lc}，航空器的高度剖面变化过程为：先保持平飞至距离高度限制点之前Lmin的位置，然后执行下降过程到指定高度，然后平飞至高度限制点；其中Lc为航迹起始点与高度限制点之间的水平距离。
2.根据权利要求1所述的一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，所述步骤1中，飞行计划中包含了计算高度剖面所需的起飞机场、航路点、降落机场和巡航高度信息，还包含了航空器机型；空域扇区结构用于确定高空扇区和低空扇区，以及航路点不同高度所对应的扇区类型；进离场程序用于计算离场和进场阶段的高度剖面；航路点的高度限制来源于移交协议、流量管理或管制提出的限制；飞行计划中的第一个航路点和最后一个航路点分别对应于离场点和进场点，航路点的高度限制根据空管运行规则，在6000米以下。
3.根据权利要求2所述的一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，所述步骤2中，根据起飞机场的离场程序计算确定航空器从跑道到第一个航路点的飞行距离，爬升过程的计算分为加速爬升、低等表爬升和高等表爬升3个子过程，高度剖面计算公式为
式中：h为高度，t为时间，m为航空器质量，VTas为航空器真空速，g为重力加速度，T为航空器推力，D为航空器阻力；
为能量分配系数，表示按照选定的速度进行爬升时，用于爬升的推力与用于加速的推力的比值；通过航空器机型能够从欧控发布的BADA数据库获得航空器在不同阶段和高度对应的航空器推力T、航空器阻力D和航空器真空速VTas；爬升过程采用逐步迭代计算，飞机的重量在每个时间间隔
内保持不变，确定每个时间间隔
内的航空器真空速VTas，根据上述公式计算出各个时间间隔
内的爬升率，依次完成起飞机场到第一个航路点爬升过程的高度剖面计算；如果在第一个航路点之前达到指定的高度限制，则保持平飞状态至第一个航路点。
4.根据权利要求3所述的一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，所述步骤3中，从航迹起始点开始向后找到下一个高度限制点，航迹起始点到该高度限制点前一个扇区进入边界的高度剖面计算公式和步骤2中相同，处理逻辑如下：当已完成计算的高度剖面未曾达到巡航高度，即航班尚未进入巡航状态时，如果该高度限制点与起飞机场水平距离小于200公里，则先爬升至该限制高度，然后保持平飞；否则先爬升到巡航高度，再保持平飞；如果航班已经进入巡航阶段，则保持平飞状态至该高度限制点前一个扇区的进入边界。
5.根据权利要求4所述的一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，所述步骤4中，当航迹起始点处于高空扇区，并且限制点的指定高度低于航迹起始点的飞行高度时，则根据空域扇区结构判断该指定高度是否处于低空扇区空域内，以此来判定航空器是否将由高空扇区飞向低空扇区。
6.根据权利要求5所述的一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，所述步骤6中根据步骤5确定的高度剖面变化过程，计算航迹起始点到高度限制点的高度剖面计算公式和步骤2中相同。
7.根据权利要求6所述的一种基于指定高度限制的航空器飞行高度剖面计算方法，其特征在于，所述步骤8中进场阶段的高度剖面计算公式和步骤2中相同，计算并确定从当前高度下降到目的机场高度所需的水平距离Ld，航空器从当前高度保持平飞至距离目的机场Ld的位置，然后开始下降到目的机场；根据步骤2中高度剖面计算公式计算获得从当前高度下降到目的机场高度的下降率，再根据当前高度下降到机场高度的高度差计算出当前高度下降到机场高度所需要的时间
由时间
和此阶段航空器真空速VTas和风速Vwind计算获得Ld：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110116380.5A CN112927562B (zh) | 2021-01-28 | 2021-01-28 | 一种基于指定高度限制的航空器飞行高度剖面计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110116380.5A CN112927562B (zh) | 2021-01-28 | 2021-01-28 | 一种基于指定高度限制的航空器飞行高度剖面计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112927562A CN112927562A (zh) | 2021-06-08 |
| CN112927562B true CN112927562B (zh) | 2022-03-04 |
Family
ID=76167946
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110116380.5A Active CN112927562B (zh) | 2021-01-28 | 2021-01-28 | 一种基于指定高度限制的航空器飞行高度剖面计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112927562B (zh) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114047774B (zh) * | 2021-11-11 | 2023-11-14 | 中国电子科技集团公司第二十八研究所 | 一种四旋翼无人机在低空分层空域运行仿真的验证方法 |
| CN114333426A (zh) * | 2021-11-26 | 2022-04-12 | 中国航空无线电电子研究所 | 一种基于多航路点高度窗约束的进场飞行参数优化方法 |
| CN114912688B (zh) * | 2022-05-18 | 2024-04-30 | 中国电子科技集团公司第二十八研究所 | 一种基于爬升段能量分享因子估算的航空器轨迹预测方法 |
| CN115206135B (zh) * | 2022-06-16 | 2024-02-13 | 中国电子科技集团公司第二十八研究所 | 一种不确定爬升率的航空器指令高度规划方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6134500A (en) * | 1999-06-03 | 2000-10-17 | United Air Lines, Inc. | System and method for generating optimal flight plans for airline operations control |
| FR2935495A1 (fr) * | 2008-08-29 | 2010-03-05 | Thales Sa | Systeme et procede de calcul des predictions de vol par sections verticales |
| CN101692315A (zh) * | 2009-09-25 | 2010-04-07 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
| CN101739845A (zh) * | 2009-12-18 | 2010-06-16 | 中国航空无线电电子研究所 | 基于航空数据链信息的民机航空电子验证***及其方法 |
| CN103093649A (zh) * | 2011-10-31 | 2013-05-08 | 通用电气公司 | 用于推断飞行器参数的方法和*** |
| CN104240541A (zh) * | 2014-09-09 | 2014-12-24 | 中国电子科技集团公司第二十八研究所 | 一种4d航迹生成方法 |
| CN105701552A (zh) * | 2014-10-11 | 2016-06-22 | 中国国际航空股份有限公司 | 一种飞行航路垂直剖面的计算方法 |
| CN108693890A (zh) * | 2017-04-04 | 2018-10-23 | 泰勒斯公司 | 飞行器下降阶段管理 |
| CN108766035A (zh) * | 2018-06-12 | 2018-11-06 | 云南农业大学 | 一种点密度引导下的无人机地形匹配飞行控制*** |
| CN110930770A (zh) * | 2019-11-06 | 2020-03-27 | 南京莱斯信息技术股份有限公司 | 一种基于管制意图和飞机性能模型的四维航迹预测方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2983594B1 (fr) * | 2011-12-02 | 2014-09-26 | Thales Sa | Procede de gestion d'un plan de vol vertical |
| CN109191925B (zh) * | 2018-10-17 | 2020-07-14 | 中国电子科技集团公司第二十八研究所 | 一种面向四维航迹运行的多空域轨迹规划与协商方法 |
| CN111444174B (zh) * | 2020-02-25 | 2024-06-14 | 民航数据通信有限责任公司 | 一种飞行程序航迹数据检查和控制的方法及装置 |
-
2021
- 2021-01-28 CN CN202110116380.5A patent/CN112927562B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6134500A (en) * | 1999-06-03 | 2000-10-17 | United Air Lines, Inc. | System and method for generating optimal flight plans for airline operations control |
| FR2935495A1 (fr) * | 2008-08-29 | 2010-03-05 | Thales Sa | Systeme et procede de calcul des predictions de vol par sections verticales |
| CN101692315A (zh) * | 2009-09-25 | 2010-04-07 | 民航总局空管局技术中心 | 基于实时雷达数据的高精度4d飞机航迹分析方法 |
| CN101739845A (zh) * | 2009-12-18 | 2010-06-16 | 中国航空无线电电子研究所 | 基于航空数据链信息的民机航空电子验证***及其方法 |
| CN103093649A (zh) * | 2011-10-31 | 2013-05-08 | 通用电气公司 | 用于推断飞行器参数的方法和*** |
| CN104240541A (zh) * | 2014-09-09 | 2014-12-24 | 中国电子科技集团公司第二十八研究所 | 一种4d航迹生成方法 |
| CN105701552A (zh) * | 2014-10-11 | 2016-06-22 | 中国国际航空股份有限公司 | 一种飞行航路垂直剖面的计算方法 |
| CN108693890A (zh) * | 2017-04-04 | 2018-10-23 | 泰勒斯公司 | 飞行器下降阶段管理 |
| CN108766035A (zh) * | 2018-06-12 | 2018-11-06 | 云南农业大学 | 一种点密度引导下的无人机地形匹配飞行控制*** |
| CN110930770A (zh) * | 2019-11-06 | 2020-03-27 | 南京莱斯信息技术股份有限公司 | 一种基于管制意图和飞机性能模型的四维航迹预测方法 |
Non-Patent Citations (3)
| Title |
|---|
| THALES自动化***飞行轨迹计算与投递原理浅析;陈玉平;《数字技术与应用》;20170115;全文 * |
| 航空器四维航迹实时预测与动态生成研究;刘杰;《中国优秀硕士学位论文全文数据库》;20200215;全文 * |
| 面向航空器战略航迹规划的风场信息插值方法;汤新民等;《指挥信息***与技术》;20200116;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112927562A (zh) | 2021-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112927562B (zh) | 一种基于指定高度限制的航空器飞行高度剖面计算方法 | |
| CN110930770B (zh) | 一种基于管制意图和飞机性能模型的四维航迹预测方法 | |
| EP3598262B1 (en) | Method and system for determining a climb profile | |
| US8645009B2 (en) | Method for flying an aircraft along a flight path | |
| CN110059863B (zh) | 一种基于所需到达时间的航空器四维航迹优化方法 | |
| EP2631890A2 (en) | Methods for in-flight adjusting of a flight plan | |
| CN105701552B (zh) | 一种飞行航路垂直剖面的确定方法 | |
| US11094206B2 (en) | Vertical flightpath optimization | |
| CN114283624B (zh) | 一种基于空地航迹信息共享的飞机位置预测方法 | |
| CN112669651B (zh) | 一种基于飞行动态情报中eet值修正预测过点时间的方法 | |
| CN109191921B (zh) | 一种基于4d轨迹预测的空管监视数据仿真方法 | |
| EP3598261B1 (en) | Method and system for determining a descent profile | |
| Dancila et al. | Vertical flight profile optimization for a cruise segment with RTA constraints | |
| Gong et al. | Using flight manual data to derive aero-propulsive models for predicting aircraft trajectories | |
| CN113470438B (zh) | 一种基于逻辑时序推演仿真的无冲突航班轨迹生成方法 | |
| CN112037583B (zh) | 一种航空器高度层变更程序优化引导方法 | |
| CN111142561A (zh) | 一种飞行管理***的垂直高度导引方法 | |
| CN112882488B (zh) | 航空器4d轨迹预测方法及装置 | |
| CN117406767A (zh) | 一种基于飞行计划的航空器飞行轨迹快速生成方法 | |
| De Jong et al. | Time and energy management during descent: Human vs automated response | |
| CN114333426A (zh) | 一种基于多航路点高度窗约束的进场飞行参数优化方法 | |
| CN117351783A (zh) | 一种基于机载下传数据的空管***轨迹预测修正方法 | |
| CN116466750A (zh) | 一种基于组合模型的飞机垂直轨迹预测与优化方法 | |
| CN115206135A (zh) | 一种不确定爬升率的航空器指令高度规划方法 | |
| CN111123977A (zh) | 一种飞行管理***的时间导引方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 210000 No.1, Lingshan South Road, Qixia District, Nanjing City, Jiangsu Province Applicant after: THE 28TH RESEARCH INSTITUTE OF CHINA ELECTRONICS TECHNOLOGY Group Corp. Address before: 210007 1 East Street, alfalfa garden, Qinhuai District, Nanjing, Jiangsu. Applicant before: THE 28TH RESEARCH INSTITUTE OF CHINA ELECTRONICS TECHNOLOGY Group Corp. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |