一种基于概率分布的航段运行时间计算方法
【专利摘要】本发明涉及一种基于概率分布的航段运行时间计算方法，步骤包括通过分析空域实际运行情况和机型性能数据，对城市对航段运行过程中运行阶段进行划分和对机型进行分类；对不完整的航班历史运行记录数据和不正常航班的历史运行记录数据清洗；按照运行阶段设定初步样本空间，其中运行阶段中滑出/滑入段再按航站楼位置设定样本空间，运行阶段中飞行段再按机型的巡航速度以及所属航空公司设定样本空间，并在初步样本空间的基础上分析各个运行阶段的运行时间概率分布情况，确定各个运行阶段的有效样本范围；按航站楼计算滑入/滑出阶段运行时间；计算飞行段运行时间；综合滑出、飞行，滑入三个阶段的运行时间得到航段运行综合时间。
【专利说明】一种基于概率分布的航段运行时间计算方法
【技术领域】
[0001] 本发明公开了一种基于概率分布的航段运行时间计算方法，涉及城市对航段运行 时间的度量研究领域。
【背景技术】
[0002] 航段标准运行时间是飞行计划编排过程中一项重要的参考依据。预先飞行计划编 排人员使用航班标准运行时间审查航空公司飞行计划编排中航段衔接是否合理，飞行计划 编排是否科学，机场运营是否高效以及安全性是否可靠等。
[0003] 随着我国国民经济的快速发展，日益增长的航空运输需求与当前预先制定的飞行 计划的航段标准运行时间矛盾越加显著。目前，航段运行时间都是由各航空公司在每年两 次航班换季前的航班时刻协调过程中提出。各航空公司所提的航段运行时间都是从自身 出发，不同航段、不同机型、不同的巡航方式等都会造成航段运行时间的不一致。由于预先 飞行计划管理部门缺乏技术手段和数据统计资料，确定班机航段标准运行时间缺乏数据保 证，无法均衡各个航空公司的飞行计划编排需求，导致机场时隙资源的浪费以及航班延误 的增加。
[0004] 总体而言，目前关于航段运行时间计算方法研究现状如下：
[0005] (1)在依据历史运行数据，基于其概率分布的航段运行时间的研究尚处于空白，对 预先航班计划的科学管理缺乏支持；
[0006] (2)航段运行时间是由各公司依据各自不同情况给出建议，其标准与出发点不一 致，造成其局限性；
[0007] (3)某些航段运行时间计算，以航空器性能及航路航线设计，进行仿真计算，评估 航段运行时间，为理想的目标时间，未能反映机场、空域当前的保障服务能力，因而对实际 运行的指导性降低；
[0008] (4)航段运行的滑入段与滑出段的时间采用固定值计算，无法体现实际运行情况。
【发明内容】
[0009] 本发明所要解决的技术问题是：针对现有技术的缺陷，进行相应的设计和优化，提 供一种基于概率分布的航段运行时间计算方法，通过对航段运行阶段的拆分，设定基于概 率分布的样本空间，并根据层次单排序法对各个航空公司权重进行确定，从而计算出综合 的航段运行时间，为空域运行于流量管理的战略与预战术阶段提供数据依据。
[0010] 本发明为解决上述技术问题采用以下技术方案：
[0011] 一种基于概率分布的航段运行时间计算方法，具体步骤如下：
[0012] 步骤一、根据空域实际运行情况，将城市对航段的运行过程进行划分；对同一条运 行航段上的不同机型，依据各机型性能数据进行分类；
[0013] 步骤二、对航班的历史运行记录进行数据清洗，数据清洗的内容包括对不完整数 据的清洗以及对不正常航班数据的清洗；
[0014] 步骤三、对航段运行各阶段基于概率分布的样本空间进行设定，具体包括：
[0015] (301)航段运行各阶段样本空间初步设定：将运行阶段中滑出/滑入段按照航站 楼位置设定样本空间；将运行阶段中飞行段按照机型的巡航速度以及所属航空公司设定样 本空间；
[0016] (302)基于概率分布对步骤（301)中初步设定的样本空间进行有效样本范围的确 定，去除样本中的极端情况和低概率情况，依据航班运行历史数据，分别分析城市对航段滑 出、飞行、滑入三个阶段的运行时间概率分布情况，所述时间概率分布情况符合高斯曲线分 布特征；
[0017] 步骤四、计算滑入/滑出段运行时间；
[0018] 按起飞机场与目的地机场的航站楼区分，基于步骤三中确定的滑出与滑入的有效 样本空间，直接求取滑出与滑入时间的平均值；
[0019] 步骤五、计算飞行段运行时间；
[0020] 步骤六、综合度量航段运行时间；
[0021] 根据步骤一中对航段运行过程阶段的划分，分阶段运行时间共同构成了航段综合 运行时间，即：综合飞行时间=滑出段运行时间+飞行段运行时间+滑入段运行时间。
[0022] 作为本发明的进一步优选方案，所述步骤一中：
[0023] 城市对航段的运行过程分为三个阶段，时间由撤轮挡开始到加轮挡结束，依次包 括滑出段、飞行段和滑入段，其中飞行段具体包括起飞、爬升、巡航、下降、着陆；
[0024] 所述依据各机型性能数据进行分类，以经济巡航速度为性能参考，具体包括： [0· 8M, 0· 89M]、[0· 7M, 0· 79M]、[0· 6M, 0· 69M]以及其他类型，其中 Μ 表示音速。
[0025] 作为本发明的进一步优选方案，所述步骤二中：
[0026] 不完整数据包括：实际的撤轮挡时间、起飞时间、落地时间、加轮挡时间缺失的历 史运行记录数据；
[0027] 不正常航班数据包括：
[0028] (201)航班计划、领航计划报以及实际运行这三者中，起飞机场不一致或是降落地 机场不一致的航班数据；
[0029] (202)撤轮挡时间、起飞时间、落地时间、加轮挡时间时序顺序不合理以及数值异 常的航班数据。
[0030] 作为本发明的进一步优选方案，所述步骤（302)具体包括：
[0031] 基于对航段运行各阶段运行时间高斯分布的特征，确定有效样本范围，所述有效 样本范围包括置信区间和置信度；
[0032] 所述分析各个阶段的时间概率分布情况具体为：根据所设定的样本空间，画出样 本的时间概率分布图，其中，横轴为飞行时间，纵轴为样本数据占总数据的比例，依据横轴 最短距离包容的原则，按照置信度选取总数据量中的对应数据作为样本的置信区间。
[0033] (501)基于步骤三中确定的飞行段的有效样本空间，直接求取航段所含有的各个 航空公司的飞行段运行时间平均值；
[0034] (502)计算各个航空公司在航段上的权重；
[0035] (503)根据（502)中得出的各个航空公司在航段上的权重，计算飞行段运行时间。
[0036] 作为本发明的进一步优选方案，在所述步骤（502)中，依据该航段上各航空公司 历史飞行流量，同时结合该航段上各航空公司飞行时间偏差、延误情况的历史数据，采用层 次单排序法计算各个航空公司在该航段上的权重，具体步骤包括：
[0037] (5021)建立单排序的层次结构模型；
[0038] (5022)构造判断矩阵；
[0039] (5023)层次单排序及一致性检验。
[0040] 本发明采用以上技术方案与现有技术相比，具有以下技术效果：采用本发明所公 开的方法制定的航段时间，考虑了地面滑行和飞行阶段的运行特点，避免了传统单一航段 运行时间所带来的局限性与公平性，可以客观的反映在当前空域环境和空管保障水平下城 市对航班运行的实际情况。将此航段时间应用于飞行计划编排过程中，可作为战略与预战 术阶段的一项重要的参考依据，使飞行计划编排中的航段衔接更加合理，减少机场时隙资 源浪费。同时，有助于解决各航空公司在航班时刻协调过程中存在的差异性，平衡各个航空 公司的飞行计划编排需求。
【专利附图】
【附图说明】
[0041] 图1为本发明的核心流程图；
[0042] 图2为航段运行各阶段基于概率分布的样本空间设定的流程图；
[0043] 图3为广州至北京飞行时间的高斯分布图；
[0044] 图4为基于概率分布设定有效样本范围的示例图；
[0045] 图5为使用层次单排序法计算飞行段时间的流程图；
[0046] 图6为使用层次单排序法建立的层次结构模型。
【具体实施方式】
[0047] 下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
[0048] 图1为本发明技术方案的核心流程图，所述基于概率分布特性的航段运行时间计 算方法具体实施步骤如下：
[0049] 步骤一：城市对航段运行过程阶段划分及对机型分类。
[0050] (101)城市对航段是由具有一定商业载量的航空器在两个城市间从事定期运输服 务而形成的航空运输线，根据空域实际运行情况，对城市对航段运行阶段（以下简称运行 阶段）进行定义，运行阶段指从撤轮挡开始，到加轮挡结束，依次包括滑出、起飞、爬升、巡 航、下降、着陆、滑入七个阶段，其中起飞、爬升、巡航、下降、着陆定义为飞行段，即运行阶段 可分三个阶段：滑出段、飞行段、滑入段；
[0051] (102)机型速度分类：同一条运行航段上，不同机型对航段运行时间的影响主 要由其速度所决定，而不同机型的速度的差别主要体现在航路巡航阶段，因此依据各机 型性能数据，按照机型飞行正常巡航马赫数对机型进行分类，依次分为：[0. 8M，0. 89M]、 [0· 7M, 0· 79M]、[0· 6M, 0· 69M]、及其他类型。
[0052] 步骤二：航班的历史运行记录数据清洗。
[0053] (201)对不完整数据的清洗：对实际的撤轮挡时间、起飞时间、落地时间、加轮挡 时间缺失的历史运行记录数据进行去除；
[0054] (202)对不正常航班数据的清洗：
[0055] 对航班计划、FPL报（领航计划报）以及实际运行三者中的起飞机场不一致或是 降落地机场不一致的航班数据应进行去除。
[0056] 对撤轮挡时间、起飞时间、落地时间、加轮挡时间时序顺序不合理以及数值异常等 的航班数据应进行去除。
[0057] 步骤三：航段运行各阶段基于概率分布的样本空间设定，其流程图如图2所示。
[0058] (301)航段运行各阶段样本空间初步设定：即按照运行阶段设定样本空间，其中 运行阶段中滑出/滑入段再按航站楼位置设定样本空间，运行阶段中飞行段再按机型的巡 航速度以及所属航空公司设定样本空间。
[0059] 其中飞机的滑出段、飞行段、滑入段的运行时间计算规则如下：
[0060] 滑出段：
[0061]
【权利要求】
1. 一种基于概率分布的航段运行时间计算方法，其特征在于，具体步骤如下： 步骤一、根据空域实际运行情况，将城市对航段的运行过程进行划分；对同一条运行航 段上的不同机型，依据各机型性能数据进行分类； 步骤二、对航班的历史运行记录进行数据清洗，数据清洗的内容包括对不完整数据的 清洗以及对不正常航班数据的清洗； 步骤三、对航段运行各阶段基于概率分布的样本空间进行设定，具体包括： (301) 航段运行各阶段样本空间初步设定：将运行阶段中滑出/滑入段按照航站楼位 置设定样本空间；将运行阶段中飞行段按照机型的巡航速度以及所属航空公司设定样本空 间； (302) 基于概率分布对步骤（301)中初步设定的样本空间进行有效样本范围的确定， 去除样本中的极端情况和低概率情况，依据航班运行历史数据，分别分析城市对航段滑出、 飞行、滑入三个阶段的运行时间概率分布情况； 步骤四、计算滑入/滑出段运行时间； 按起飞机场与目的地机场的航站楼区分，基于步骤三中确定的滑出与滑入的有效样本 空间，直接求取滑出与滑入时间的平均值； 步骤五、计算飞行段运行时间； 步骤六、综合度量航段运行时间； 根据步骤一中对航段运行过程阶段的划分，将各阶段运行时间之和作为航段综合运行 时间，即：综合飞行时间=滑出段运行时间+飞行段运行时间+滑入段运行时间。
2. 如权利要求1所述的一种基于概率分布的航段运行时间计算方法，其特征在于，所 述步骤一中： 城市对航段的运行过程分为三个阶段，时间从撤轮挡开始到加轮挡结束，依次包括滑 出段、飞行段和滑入段，其中飞行段具体包括起飞、爬升、巡航、下降、着陆； 所述依据各机型性能数据进行分类，以经济巡航速度为性能参考，具体包括： [0· 8M, 0· 89M]、[0· 7M, 0· 79M]、[0· 6M, 0· 69M]，其中，Μ 表示音速。
3. 如权利要求1所述的一种基于概率分布的航段运行时间计算方法，其特征在于：所 述步骤二中： 不完整数据包括：实际的撤轮挡时间、起飞时间、落地时间、加轮挡时间缺失的历史运 行记录数据； 不正常航班数据包括： (201) 航班计划、领航计划报以及实际运行这三者中，起飞机场不一致或是降落地机场 不一致的航班数据； (202) 撤轮挡时间、起飞时间、落地时间、加轮挡时间时序顺序不合理以及数值异常的 航班数据。
4. 如权利要求1所述的一种基于概率分布的航段运行时间计算方法，其特征在于，所 述步骤（302)具体包括： 基于对航段运行各阶段运行时间高斯分布的特征，确定有效样本范围，所述有效样本 范围包括置信区间和置信度； 所述分析各个阶段的时间概率分布情况具体为：根据所设定的样本空间，画出样本的 时间概率分布图，其中，横轴为飞行时间，纵轴为样本数据占总数据的比例，依据横轴最短 距离包容的原则，按照置信度选取总数据量中的对应数据作为样本的置信区间。
5. 如权利要求1所述的一种基于概率分布的航段运行时间计算方法，其特征在于，所 述步骤五具体包括： (501) 基于步骤三中确定的飞行段的有效样本范围，直接求取航段所含有的各个航空 公司的飞行段运行时间平均值； (502) 计算各个航空公司在航段上的权重； (503) 根据（502)中得出的各个航空公司在航段上的权重，计算飞行段运行时间。
6. 如权利要求5所述的一种基于概率分布的航段运行时间计算方法，其特征在于：在 所述步骤（502)中，依据该航段上各航空公司历史飞行流量，同时结合该航段上各航空公 司飞行时间偏差、延误情况的历史数据，采用层次单排序法计算各个航空公司在该航段上 的权重，具体步骤包括： (5021) 建立单排序的层次结构模型； (5022) 构造判断矩阵； (5023) 层次单排序及一致性检验。
【文档编号】G06Q10/06GK104156805SQ201410342987
【公开日】2014年11月19日 申请日期:2014年7月17日 优先权日:2014年7月17日
【发明者】袁立罡, 谢华, 胡明华, 王鹏鹏, 王兵, 李 杰, 丛玮 申请人:南京航空航天大学